Mineralogical and Micro-structural Investigation into the Mechanical Behaviour of a Soft Calcareous Mudstone

The construction industry in Abu Dhabi is thriving and its coastline has some of the most ambitious structures in the world. Whilst the subsurface evaporitic and calcareous soft rocks of this region are of great geological interest, they are relatively poorly understood from a geotechnical engineering perspective, forcing foundation designs to be overly conservative. Understanding the stiffness of the underlying geology at small strains is of great importance for the accurate estimation of ground movements around excavations and foundations, and yet routine post-SI laboratory testing programmes tend to focus on basic rock mechanics tests such as UCS tests. These procedures are generally unsuitable for use with calcareous rocks due to their friable and moisture sensitive nature, and rarely obtain parameters representative of actual in situ behaviour. The calcareous mudstone investigated in this paper has mechanical and structural characteristics falling between those of a soil and those typical of a rock and, as such, requires a geotechnical testing approach that combines methods from both soil and rock mechanics disciplines. The mineralogical, micro-structural and mechanical characteristics of this lithology have been examined via a suite of testing techniques, including XRPD, SEM, advanced triaxial with bender elements, along with industry standard procedures. Shearing, tensile and consolidation behaviours have been explored. Examination of the micro- and macro-scale features of this material shows it to be highly structured, with strength and stiffness being controlled by inter-granular bonding of Dolomite grains, as well as by mean effective stress state and rate of strain. The presence of fibrous Palygorskite acts to reduce the degree of bonding, causing specimens rich in this clay mineral to have a more ductile mechanical behaviour.


Introduction
Many of the coastal areas of Abu Dhabi in the United Arab Emirates (UAE) are heavily built up urban conurbations and Abu Dhabi Island is home to some of the tallest buildings in the world. The young soft rocks that make up the subsurface of this part of the Middle East are well studied from a geological standpoint. However, research into their complex mechanical behaviour is less common (Cook 1999;Epps 2011). The geotechnical properties of the soft rocks encountered around the Abu Dhabi coastline are often poorly understood and their unique physical characteristics cause problems when attempting to sample and test them, making the application of many standard laboratory testing methodologies inappropriate (Wright et al. 1999;Kanji 2014). This ultimately leads to overly conservative geotechnical design solutions being used to compensate for characteristic uncertainties.
Since the end of the Palaeozoic, the climate of the UAE has remained broadly similar to that of the present day, with its solid and drift geology defined by repeated cycles of glacio-eustatic triggered marine transgressions in and around the Arabian Gulf and subsequent sub-aerial erosion and transportation (Walkden and Williams 1998). These cycles have allowed the formation of vast areas of carbonate and evaporite sediments from the Gulf's saline waters, including fossiliferous limestones and dolomites along with siltstones, 1 3 sandstones and inter-bedded gypsum, halite and anhydrite as well as the calcareous mudstones studied in this paper (Walkden and Williams 1998;Evans et al. 1973).
These mudstones were likely deposited within a low energy but chemically dynamic hyper-saline environment, and are predominantly biogenic/chemical in origin, as opposed to purely detrital. In lithologies such as these, the time gap between the formation of clasts and their bonding materials can be small, or even occur contemporaneously. Weak cementitious bridges may precipitate within voids far more rapidly than, for example, the formation of iron oxide or alumina bonds in terrestrial and freshwater sediments, introducing a degree of structure that needs to be taken into account when attempting to characterise their behaviour (Aversa and Evangelista 1998;Epps 2011;Truong et al. 2012;Fereidooni 2016). Sediments formed in hyper-saline waters are noted to be particularly sensitive to changing physico-chemical environments, which can further complicate the process of carrying out routine geotechnical laboratory tests (which typically demand full saturation), or rock testing procedures (often requiring samples to be dry (Burland and Lord 1970;Hawkins 1998)). As effective stress triaxial tests are carried out on fully saturated samples (often using deaired water), it is important that the minerals making up the samples do not dissolve over prolonged exposure to water (Jeong et al. 2007;Ciantia et al. 2015).
The young soft rocks of this study contain a percentage of evaporite minerals within their granular/cement structure that are susceptible to dissolution. Therefore, the sensitivity of these soft rocks to water ingress must be assessed and the type and quantity of the dissolvable minerals (as well as non-soluble minerals) within the soft rock structure must be understood (Liu et al. 2019). In addition, the calcareous mudstones investigated in this work are highly friable and even brittle at lower moisture contents, making the measurement of local deformations using standard strain gauges designed for application to hard rocks (or even apparatus designed for soils) particularly difficult (Lagioia and Nova 1995;Lagioia et al. 1998).
This paper presents both the methodologies and findings of a custom suite of characterisation, optical and mechanical tests carried out upon one extensively occurring soft rock lithology from the Abu Dhabi coastline. In particular, a testing programme including the use of X-ray Powder diffraction (XRPD) and scanning electron microscopy (SEM), along with standard soil and rock characterisation tests is described in the paper (1) to better highlight the degree to which dissolution may be occurring; (2) to gain insight into the mineralogy of these soft rocks and their cementing agents; (3) to assess and characterise their physical and mechanical properties, with a particular focus on the smallstrain stiffness behaviour in triaxial compression. This is of particular relevance to the design of tall buildings and underground structures in the region.

Site Description and Geological Setting
The geology of the area around Abu Dhabi Island is highly variable both laterally and vertically due to the past and present dynamic sedimentary environment (Kässler 1973). Water salinity is high and average water depth low, producing a system capable of precipitating carbonates and evaporites rapidly within a short geological timescale. This area is being heavily developed with ambitious construction projects, and yet serious engineering concerns ranging from corrosive and aggressive groundwater, dissolution of evaporites/associated subsidence (often triggered by anthropogenic activities such as ground-source heating/cooling) and piling through Caprock continue to complicate design (Epps 2011). The calcareous mudstone studied was formed in a low energy, hyper-saline environment, with a modern analogy being the area's lagoonal terraces. Significant thicknesses of this material are noted across the area of study, often interbedded with gypsum. Samples tested in this paper arose from two site investigations in Abu Dhabi. The two sites are close to one another from a geological standpoint, and they are both underlain by similar superficial deposits and solid geology. The full extents of the sites are shown in Fig. 1, and the sites will be referred to as Site 1 and Site 2.

Site 1
This investigation covered a long, narrow stretch of land, and the nature of the paleogeography (comprising islands, lagoons and channels) resulted in a laterally discontinuously layered geology. Figure 1a shows a representative cross-section of Site 1. Generalised ground conditions encountered can be summarised as follows: medium dense to dense, poorly graded silty sand with occasional cemented layers, overlaying weak to very weak sandstone/calcarenite layers. Below this, a weak mudstone, inter-bedded with gypsum, is encountered. Both on-land and offshore boreholes were undertaken and produced cores of approximately 76 mm diameter. Cores were generally carried out in run lengths of 2.50 m, with total core recovery of between 90 and 95% within beds where sub-samples were taken.

Site 2
This site is located to the south of Site 1, and ground conditions are described as medium dense, silty, gravelly, fine to medium siliceous carbonate sand overlaying extremely weak to very weak sandstone/calcarenite. The deepest layer encountered is described as an extremely weak to very weak siltstone inter-bedded with gypsum and mudstone, which continues to the full depth of the exploratory holes. Figure 1b shows a representative cross-section of Site 2. It should be noted that the material logged as siltstone at Site 2 is mechanically and characteristically similar to that of the mudstone layers noted at Site 1. Indeed, particle size distribution tests and plasticity indices confirm their comparable nature and predominantly clay sized grains. For the purpose of clarity, the soft rock referred to as mudstone at Site 1 and erroneously as siltstone at Site 2 will be referred to as calcareous mudstone going forwards, as this designation better describes the material of interest within this paper.

Experimental Results
Eight core runs of the calcareous mudstone were obtained from boreholes within Site 1 along with eight core runs from Site 2, from which sub-samples were extracted for the various experimental testing as required. Whilst all calcareous mudstone cores appeared superficially similar, preliminary testing revealed that samples could be reliably divided into two distinct sub-types (namely Type A and Type B) based upon their behaviour and characteristics. For example, Type A samples failed in triaxial compression following a high peak stress with a high stiffness, whilst Type B samples failed at larger strains and lower stresses under the same conditions. On this basis, specimen names have been simplified as A1 through A9 for Type A materials and B1 through B7 for Type B materials for ease of reading, and all samples have been subjected to a broad suite of experimental testing, the results of which are presented in this paper ordered from micro to macro. Optical techniques such as XRPD and SEM are presented first, as they focus on very small-scale features, followed by the results of standard characterisation tests and ultimately the findings of advanced triaxial tests.

Petrographic Study: XRPD
Two sets of qualitative XRPD tests on powdered bulk samples were carried out to characterise both mineralogy and the potential effect of dissolution. A PANalytical X'Pert Pro Multipurpose Diffractometer (fitted with an X'Celerator and secondary monochromator) was used to identify likely mineral phases present within the calcareous mudstone. Samples were scanned in continuous mode using the X'Celerator RTMS detector. Samples displaying both Type A and Type B behaviour were selected to see if a mineralogical discrepancy might account for their mechanical differences. In addition, sub-samples of all specimens were prepared by exposing them first to flowing deionised water, to assess the effect of dissolution on a mineralogical scale.
XRPD results are shown in Fig. 2a and b with notable mineral phases highlighted and notable 2-Theta angle peaks annotated. However, it should be noted that not all mineral 2-Theta peaks appear on this image for the sake of readability. All samples show similar mineralogical finger prints, displaying comparable peak angles and peak heights between samples taken, indicating that they are qualitatively composed of similar minerals. Results reveal the presence of Albite, Anorthoclase, Quartz, Halite and the clay mineral Palygorskite in all samples. However, Dolomite is clearly the prevailing mineral type with very high intensity counts. Interestingly, no Calcite is noted in any of the samples tested. No Gypsum or Anhydrite is detected in any of the samples tested, despite the in situ proximity of the sample cores to known gypsum beds.
Dolomite is a common mineral to be found associated with sediments and rocks formed in hyper-saline water bodies. However, its method of formation is a great source of controversy within the literature, as Dolomite is not forming today at the earth's surface except within extreme environments such as hot springs (Eugster and Hardie 1978;Last 1990;Sanz-Montero et al. 2006). Although the aim of the work is not to explain the origin of the Dolomite deposits at the studied sites, it is interesting to note the variety of processes and theories that might be responsible for its presence, ranging from the replacement of Calcite and Gypsum by Dolomite to the microbially mediated precipitation of Dolomite by algae and sulphate-reducing bacteria (Pierre and Rouchy 1988;Machel 2001;Sanz-Montero et al. 2006). The minor quantities of Quartz and Feldspars present in the XRPD analysis are likely to be due to occasional grains of (b) (a) detrital debris, transported from an adjacent high energy sub-environment into the low energy sub-environment that allowed the carbonate muds to form. Halite is also present in samples which have not been exposed to flowing deionised water (Fig. 2a), likely present due to the recent precipitation of residual saline pore-water from the sample drying during preparation for testing, although the potential for crystalline Halite existing within the calcareous mudstone's structure cannot be ruled out. Halite is not seen in samples that have been exposed to flowing deionised water ( Fig. 2b), owing to its dissolvable nature, potentially creating problems during effective stress triaxial tests which require samples to be saturated, typically with deaired deionised water. One mineral of significant interest is Palygorskite, which is present in all samples. Palygorskite is a hydrous magnesium alumina (and iron) silicate, a fibrous clay mineral, the origin of which is debated at length within the literature (Haydn and Huitang 2005).
Palygorskite, sometimes known as Attapulgite, is a member of the Hormite family of clays and has many industrial uses, due to its sorptive capacity, including lubricants, adhesives, paints and pharmaceuticals (Haydn and Huitang 2005;Rhouta et al. 2013). It tends to be found in carbonate sedimentary rocks and sediments deposited in high temperature/salinity marine/lagoonal waters (Isphording 1973;Rodriguez-Navarro et al. 1998;Machel 2001;Rhouta et al. 2013) as well as forming in peri-marine and non-marine environments by detrital sedimentation under varying salinity conditions. Stagnation of the source water, along with its alkalinity, the presence of Si-and Mg-rich fluids and intense evaporation are all considered to be critical to its formation (Inglès and Anadón 1991). The modern-day lagoons around the Abu Dhabi coast fluctuate in temperature and salinity throughout the year and, whilst some are connected to the main body of the Arabian Gulf, others are isolated allowing them to become hyper-saline. It is likely that these conditions were also present in the environments that formed the soft rocks of this study and that localised high concentrations of Palygorskite may have occurred within the carbonate muds prior to their lithification. The highest peak indicating Palygorskite during XRPD analysis, like many clay minerals, occurs at a low 2-Theta angles. Within this region there is much noise associated with the limitations of the test itself, making Palygorskite difficult to quantitatively measure.
XRPD analysis is concluded as being unable to highlight any mineralogical difference between the Type A and Type B material types, meaning that any differences must be due to the relative abundances of the minerals present and not due solely to the presence or lack of specific minerals.
Quantitative XRPD tests on powdered bulk samples were, therefore, carried out to characterise the samples mineralogy. Samples were again selected from the Type A and Type B specimens to ascertain if changes in the relative abundance of minerals might account for the mechanical differences noted during preliminary triaxial tests. Semi-quantitative Rietveld refinement was used to quantify the relative abundances of minerals by examination of the peak shape and size during XRPD analysis (Rietveld 1969), the results of which are summarised in Table 1 and reveal that the Type A samples have overall higher percentage fractions of Dolomite (with an average of 88.5%) and lower percentages of Palygorskite (on average equal to 6.9%), whilst having similar amounts of Quartz and Halite to the Type B samples. The Type B samples, by contrast, have significantly lower percentages of Dolomite (average of 78.0%), but higher amounts of Palygorskite (on average equal to 16.3%). These trends correspond well to the findings of the PSD analysis where upon Type B samples have a greater percentage of clay sized grains than Type A specimens.

Petrographic Study: SEM
SEM analysis was performed to corroborate the findings of the XRPD study and to obtain high resolution, high magnification images of the minerals present to qualitatively assess their relative abundances and observe their 3D micro-structure, with an aim to compare Type A and Type B samples. The study required the preparation of small broken fragments, approximately 1 cm 3 in size and displaying at least one fresh rough surface. To maximise image clarity and resolution, samples were gold coated.
Rounded (peloidal), fused Dolomite grains can be clearly observed, for example in sample A1 (Fig. 3), along with High-resolution SEM images successfully highlighted significant differences between Type A and Type B samples in terms of both their structure and the presence (or lack) of the inter-clast clay mineral Palygorskite. A clear example of this is shown in Fig. 4 by comparing a typical Type A sample (A2) and a Type B specimen (B2). The Type A image shows clasts of Dolomite fused tightly together with very little or no inter-clast minerals present (alpha), whereas the Type B image contains high quantities of a clay mineral with a micro fibrous morphology located both on and between clasts (beta). The fibrous clay mineral in question is believed to be Palygorskite, as suggested by the XRPD analysis and due to its superficial morphology being markedly similar to that of the Palygorskite tested by Neaman and Singer (2000) and Rhouta et al. (2013). It appears as mono-mineralic assemblages of straight fibres at random orientations within void spaces and between larger grains and it is thought to be formed via its precipitation from solution (Singer 1981). Merkl (1989) identified that Palygorskite may appear as short length (<10 μ m) and long length (>10 μ m) particles, signifying low and high magnesium content, respectively (Haydn and Huitang 2005).
At high magnifications, as shown in sample B5 (Fig. 5), several samples are noted to contain areas covered in frequent hollows, through which the underlying fibrous Palygorskite may be seen and the grains themselves are texturally rough, as opposed to amorphous. The Dolomite appears to have grown over the Palygorskite fibres (overgrowths), again likely due to processes of neomorphism (Jameson 1992). The existence of Palygorskite filled hollows supports the hypothesis that the clay mineral is stopping the Dolomite grains from fusing fully.

Characterisation Tests: Atterberg Limits
Atterberg limits were carried out according to BS1377-2 (1990 primarily to distinguish between Type A and Type B specimens. Liquid and plastic limits show the calcareous mudstones encountered in both sites to be highly comparable (see Table 2). Type A and Type B samples form tight clusters of data points, irrespective of depth or location, with few outliers. Type A samples have a mean plastic limit of 23.8%, liquid limit of 40.6% and a mean plasticity index of 16.8%, designating them as being of intermediate to low plasticity. Type B samples have higher values of plastic and liquid limit, equal on average to 27.7, 52.4%, respectively, and therefore, a higher average plasticity index of 24.7%, classifying them as an intermediate to high plasticity material. This increase in plasticity is interpreted to be a result of higher concentrations of Palygorskite in Type B samples, supporting the findings of the SEM analysis. The classification of the soils according to the Casagrande plasticity chart is shown in Fig. 6.

Characterisation Tests: Particle Size Distribution, Dispersivity and Slake Durability Tests
PSD tests were carried according to BS1377-2 (1990) on samples from both sites, with those from Site 1 reported in Fig. 7a and those from Site 2 reported in Fig. 7b. Samples have between 97.3 and 100% fine particles (less than 63μm in diameter) from both sites, with clay fraction between 15 and 40%, depending upon specific sample. Type B samples have, on average, 9.1% more grains passing the 0.002 mm particle size and less coarse silt sized particle than those of the Type A samples. PSD tests corroborate the findings of the SEM and Atterberg analysis in that Type B samples have a significantly higher fraction of grains smaller than 24 m than those of Type A samples. Dispersivity (Crumb) tests performed following BS1377-5 (1990 indicate the material to be highly sensitive to the presence of water in unconfined conditions, breaking down almost entirely within a short time period and producing colloids. This feature is confirmed via slake durability testing (Franklin and Chandra 1972;Brown 1986;ASTM 1990;ISRM 1979) which shows both very low first (32-38%) and second cycle (1-5%) durability indices (ID) of both the Type A and Type B materials from both sites, with no significant differences between materials being noted (Results summary shown in Table 3). Whilst this occurs for samples under zero confining pressure, it has not been observed on a macro-scale during preliminary triaxial testing under confining stress, despite samples being fully saturated.

Mechanical Study
Two sets of effective stress tests, CDTs (Consolidated Drained Triaxial), were undertaken on samples of both Type A and Type B calcareous mudstone as part of this study, the first examining the effect of changing confining stress and the second examining the effect of changing strain rate. CDT tests were conducted to investigate the long term behaviour of the soft rocks to support future designs of deep foundations at the two selected sites. The CDT tests were carried out according to the BS1377-8 (1990 albeit with slower strain rates to allow multiple data points to be obtained within the small strain and very small strain range.

Effect of Confining Stress
Suite of eight consolidated drained triaxial (CDT) tests were carried out upon Type A and Type B specimens. Samples were isotropically consolidated to initial mean effective stresses ( p ′ 0 ) of 235, 470, 940 and 1300 kPa prior to shearing at a set deformation rate equal to 0.001% axial strain per minute under drained conditions. In addition, Bender Elements (BE) were employed for the measurement of small-strain shear stiffness ( G 0 ) at varying confining stresses. Table 4 summarises the initial properties of the samples along with the conditions under which they were consolidated and sheared.
Type A samples comprise A1, A2, A3 and A4 and Type B samples comprise B1, B2, B3 and B4. Figure 8 shows the stress-strain (a) and volumetric behaviour (b) observed during the CDT tests performed with increasing mean confining pressure. All tests on Type A samples display peak behaviour, with visual inspection of samples post failure revealing well defined shear planes, occurring at shear strains no greater than 1.5%. Type A samples also have a significantly higher maximum deviator stress than the corresponding Type B specimens for any given initial stress state. Type B samples stop displaying peak behaviour between the 470 and 940 kPa p ′ 0 tests and, instead, show some degree of strain hardening, ductile behaviour. Post-failure, these samples also show shear planes, although they have a more pronounced barrelling at their centre. This likely indicates the degradation of cementitious bonds as isotropic mean effective consolidation states are increased (beyond 470 kPa). Failure occurs in Type B samples at much greater shear strains of between 2.2 and  Volumetric behaviour observed during the 235, 470 and 940 kPa Type A tests is seen to shift from being initially contractant to dilatant as shear strains increase. Only the sample A4 is noted to be purely contractant throughout the test, indicating that the high mean effective confining stresses of 1300 kPa may have caused a degradation of the material bonding, allowing a ductile behaviour. The volumetric response of Type B samples also indicates a shift from brittle to ductile behaviour as the higher mean effective stress states are reached, however this transition is seen to occur at significantly lower values of p ′ 0 (as this behaviour is apparent from the 940 kPa test onwards), corresponding to   Figure 9 presents the triaxial results in terms of stress path and stress ratio q∕p � vs deviatoric strain. The failure points for both Type A and Type B are indicated by the crosses in Fig. 9a. It can be seen that both types samples show similar failure strength. Type A samples display higher shear stiffness at lower shear strains than Type B samples tested at the same mean effective stress. The stiffness of Type A samples degrades  rapidly with the onset of shear strains until at around 0.04% where a reasonably steady plateau is noted in all but sample A4, creating a characteristic S shaped curve (Fig. 10). The stiffness of sample A4 degrades steadily from G max indicating an isotropic de-bonding caused by the high effective confining pressure, corresponding to entirely contractant volumetric behaviour. A further drop in shear stiffness is noted at between 0.6 and 1.2% shear strains, corresponding to the point immediately before peak deviator stresses are met. It is also at this point that samples consolidated to 235, 470 and 940 kPa shift from contractant to dilatant volumetric straining. Upon entering large strains, all Type A samples' stiffness converge and have a similar residual value.
Type B samples do not show the stiffness plateau noted in the Type A samples, even at low mean effective stresses. Their stiffness instead reduces at a relatively constant rate from their G max without displaying the degradation pattern. Type B stiffnesses also converge at large strains. G 0 was assessed with Bender element (BE) tests performed on Type A samples. Several techniques were considered for the interpretation of these results, concluding that the 1st bump method (e.g. Viggiani and Atkinson 1995) provides an upper boundary of G 0 values, whilst the peak to peak (e.g. Jovicic and Coop 1996) and the start to start techniques give lower boundary estimates (e.g. Fonseca et al. 2009). The measured G 0 values are shown in Fig. 10 alongside Type A stiffness

Effect of Strain Rate
All specimens used in this study were isotropically consolidated to a mean effective stress of 940 kPa, as this stress state represents the approximate change from brittle to ductile behaviour noted in the previous sections. A repeat test at 0.001%/min strain rate was carried out on sample B5 (Table 5), to allow for cross-comparison with the test performed on sample B3 discussed in the previous section. A CDT test was carried out at a rate of 0.0001% per minute on sample B6 and another test was undertaken at a rate of 0.01% per minute on sample B7 to assess the effect of strain rate. Stress-strain and volumetric plots are reported in Fig. 11. The comparison between the tests carried out at the same strain rate (i.e. samples B3 and B5) highlights the reliability and repeatability of the testing undertaken. The two plots display a very similar behaviour, despite the samples have been retrieved from different sites and the difference in their depth of origin is more than 24 m. A markedly higher maximum deviatoric stress is observed in the faster 0.01%/ min test along with a brittle strength reduction, whilst the 0.0001%/min test shows a significantly less stiff response. It should be noted that the 0.0001%/min test had to be stopped prematurely due to time constraints. Whilst all volumetric responses are overall contractant (as discussed for Type B samples in the previous section), the results indicate a gradual change in behaviour as shearing rate increases, with the slowest test (sample B6) having a continuously contractant response and the fastest test (sample B7) beginning to shift towards a dilatant behaviour after achieving its maximum volumetric strain of 0.7% at approximately 1.5% shear strain.
Small-strain stiffness analysis of the CDT results, reported in Fig. 12, show a similar G max (of proximately 400-500 MPa) for all three rates of strain, although the strain at which stiffness begins to degrade is progressively shifted to the right as shearing rate increases. It can also be noted that the degradation curve of sample B5 is in good agreement with that of sample B3, tested at the same strain rate and initial mean effective stress, thus confirming the repeatability of the adopted testing procedure.

Discussion and Conclusions
The calcareous mudstones studied in this paper can be considered soft rocks, a material type that displays attributes of both soils and rocks in terms of their mechanical and structural characteristics. All samples looked superficially similar in appearance and bulk properties. The specimens disintegrated rapidly when exposed to water in unconfined conditions during durability and dispersivity tests, but they retained their shape during triaxial testing under confining stress, despite being fully saturated.
From a micro-scale point of view, XRPD analysis showed samples to consist of predominantly clay to silt sized Dolomite grains with minor constituents of Quartz and Feldspar (appearing as rare sand sized grains), Halite (often lining Fig. 12 Secant shear stiffness degradation with changing strain rate of Type B samples grain surfaces) and Palygorskite (a clay mineral occurring within void spaces). Relative abundances of minerals varied significantly, with samples displaying high quantities of Dolomite (88-89%) and low percentages of Palygorskite (4-9%), i.e. Type A, having distinctly different mechanical characteristics to samples with higher Palygorskite (10-23%) and lower Dolomite (73-84%) quantities (Type B). Moreover, SEM tests revealed the micro-structure of the calcareous mudstones on a granular scale and showed Type A samples to display tightly bonded mosaics of Dolomite grains with small amounts of fibrous Palygorskite within voids. It is believed that continuous dissolution and re-precipitation at the contact points between Dolomite grains has led to their interlocked appearance, probably contributing to their higher strength and stiffness than Type B samples, as seen by Mazzullo (1992) and Jameson (1992). On the contrary, the presence of a greater quantity of Palygorskite in Type B samples is thought to inhibit the processes of dissolution and precipitation, preventing the strong mosaics of Dolomite grains to form.
Particle size distribution and Atterberg limits confirmed the interpretation of the XRPD and SEM analysis, indicating that Type A samples were characterised by lower void ratios, they were typically less plastic and had a smaller percentage of material passing the 2 m grain size fraction than their Type B counterparts, which instead showed higher plasticity and a greater percentage of fines.
In terms of macro-scale mechanical behaviour, the soft calcareous mudstones were observed to have strength and stiffness that are highly dependent upon mean effective stress (similar to the behaviour of hard soils/soft rocks noted by Heymann (1998) and soft sandstones investigated by Kate and Gokhale (1998)) and strain rate (comparable to structured soils (e.g. Leroueil 1998), hard soils (e.g. Silvestri and d'Onofrio 1998) and fine grained soils (e.g. Vardanega and Bolton 2011)). However, significant differences emerged as the triaxial experiments were undertaken, resulting in the identification of two distinct patterns: Type A samples showed a brittle and dilatant response during shearing, while Type B a more ductile and contractant behaviour. Considering the evidences obtained from the XRPD and SEM testing, it is believed that the presence (or absence) of large quantities of inter-clast Palygorskite within the specimens is responsible for the differences observed in their strength and stiffness response. The decrease in material strength with increased Palygorskite content is thought to be related to earlier diagenetic processes having been inhibited, thus resulting in weaker bonded Dolomite grains. In fact, both Type A and B samples shared a similar effective angle of internal friction at failure and post failure, but Type B specimens showed a lower effective cohesion due to a poorer degree of bonding. This is also reflected on the stiffness response, with the better bonded Type A samples being stiffer than their Type B counterparts at any given mean effective stress. Therefore, Palygorskite relative abundances seem to be a key controlling factor of the overall behaviour of the investigated rocks during shearing.
Concluding, soft calcareous mudstones should be thought of as highly structured/bonded/cemented soils, similar to materials such as Vallericca Clay (Amorosi and Rampello 2007), Pappadai Clay (Cotecchia and Chandler 1997;Leroueil 1998) or even London Clay (Gasparre et al. 2007), in terms of the degree of additional strength and stiffness due to the presence of micro-structure. Deviations in behaviour from structured soils occur in terms of the higher shear stiffness observed at lower mean effective stress, likely caused by the presence of developed bonding/cementation. In terms of practical advice, testing methodologies adopted for these soft rocks should largely follow the framework developed for cemented soils rather than that of rock mechanics tests. In addition, recommendations for further research include non-destructive testing methods, such as Needle Penetration Test, X-ray Computed Tomography Scanning and Infrared Thermography Techniques, which are being increasingly used for the determination of the mechanical properties of soft rocks.
Within this perspective, future work will be focused on the validation and calibration of existing advanced soil constitutive models accounting for structure degradation against the experimental data presented in this paper.