Discussion on “Experimental study on fracture toughness and tensile strength of a clay” [Engineering Geology 94 (2007) 64–75]

The paper [Wang, J.-J., Zhu, J.-G., Chiu, C.F., Zhang, H., 2007. Experimental study on fracture toughness and tensile strength of a clay. Engineering Geology 94, 64–75.] focuses on two important fracture parameters of soils: tensile strength and fracture toughness. These parameters control the behaviour of soils in a wide range of situations, from the design of a simple footing to much complicated fracture behaviour of clay liners or covers. The authors have done extensive laboratory work to determine these two parameters and their laborious and complicated experimental program needs praise. However, some of the points raised in their conclusions, based on the analysis and comparison with the data from the literature, need to be discussed.     

Quantitative textural analysis of packings of elongate crystals

The spatial distribution of grains in a solidifying igneous rock controls the physical properties of the crystal mush, and in turn is controlled by the rate of crystal growth and accumulation. A predominant non-spherical habit for igneous minerals brings into question the use of spherical particles in reference packings used for quantification of spatial distribution. Furthermore, variations of crystal clustering/ordering with length scale require spatial statistics which take into account the distribution of particles beyond nearest neighbours. Using random close packings of spherocylinders, we demonstrate the importance of aspect ratio for the aggregation index (usually known as R) and show that packings of spherical particles have more structure than packings of rods. The spatial distribution functions demonstrate that the plagioclase grains in the colonnade from the Holyoke basalt are clustered on a length scale of 0.5 mm. Understanding the controls on grain spatial distribution in igneous rocks will depend on the application of these techniques to well-understood environments.     

THERIA_G: a software program to numerically model prograde garnet growth

We present the software program THERIA_G, which allows for numerical simulation of garnet growth in a given volume of rock along any pressure–temperature–time (P–T–t) path. THERIA_G assumes thermodynamic equilibrium between the garnet rim and the rock matrix during growth and accounts for component fractionation associated with garnet formation as well as for intracrystalline diffusion within garnet. In addition, THERIA_G keeps track of changes in the equilibrium phase relations, which occur during garnet growth along the specified P–T–t trajectory. This is accomplished by the combination of two major modules: a Gibbs free energy minimization routine is used to calculate equilibrium phase relations including the volume and composition of successive garnet growth increments as P and T and the effective bulk rock composition change. With the second module intragranular multi-component diffusion is modelled for spherical garnet geometry. THERIA_G allows to simulate the formation of an entire garnet population, the nucleation and growth history of which is specified via the garnet crystal size frequency distribution. Garnet growth simulations with THERIA_G produce compositional profiles for the garnet porphyroblasts of each size class of a population and full information on equilibrium phase assemblages for any point along the specified P–T–t trajectory. The results of garnet growth simulation can be used to infer the P–T–t path of metamorphism from the chemical zoning of garnet porphyroblasts. With a hypothetical example of garnet growth in a pelitic rock we demonstrate that it is essential for the interpretation of the chemical zoning of garnet to account for the combined effects of the thermodynamic conditions of garnet growth, the nucleation history and intracrystalline diffusion. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.

Application of precise 142Nd/144Nd analysis of small samples to inclusions in diamonds (Finsch,South Africa) and Hadean Zircons (Jack Hills,Western Australia)

¹⁴⁶Sm–¹⁴²Nd and ¹⁴⁷Sm–¹⁴³Nd systematics were investigated in garnet inclusions in diamonds from Finsch (S. Africa) and Hadean zircons from Jack Hills (W. Australia) to assess the potential of these systems as recorders of early Earth evolution. The study of Finsch inclusions was conducted on a composite sample of 50 peridotitic pyropes with a Nd model age of 3.3 Ga. Analysis of the Jack Hills zircons was performed on 790 grains with ion microprobe ²⁰⁷Pb/²⁰⁶Pb spot ages from 3.95 to 4.19 Ga. Finsch pyropes yield 100 × ?¹⁴²Nd = ? 6 ± 12 ppm, ?¹⁴³Nd = ? 32.5, and ¹⁴⁷Sm/¹⁴⁴Nd = 0.1150. These results do not confirm previous claims for a 30 ppm ¹⁴²Nd excess in South African cratonic mantle. The lack of a ¹⁴²Nd anomaly in these inclusions suggests that isotopic heterogeneities created by early mantle differentiation were remixed at a very fine scale prior to isolation of the South African lithosphere. Alternatively, this result may indicate that only a fraction of the mantle experienced depletion during the first 400 Myr of its history. Analysis of the Jack Hills zircon composite yielded 100 × ?¹⁴²Nd = 8 ± 10 ppm, ?¹⁴³Nd = 45 ± 1, and ¹⁴⁷Sm/¹⁴⁴Nd = 0.5891. Back-calculation of this present-day ?¹⁴³Nd yields an unrealistic estimate for the initial ?¹⁴³Nd of ? 160 ?-units, clearly indicating post-crystallization disturbance of the ¹⁴⁷Sm–¹⁴³Nd system. Examination of ^146,147Sm–^142,143Nd data reveals that the Nd budget of the Jack Hills sample is dominated by non-radiogenic Nd, possibly contained in recrystallized zircon rims or secondary subsurface minerals. This secondary material is characterized by highly discordant U–Pb ages. Although the mass fraction of altered zircon is unlikely to exceed 5–10% of total sample, its high LREE content precludes a reliable evaluation of ¹⁴⁶Sm–¹⁴²Nd systematics in Jack Hills zircons.     

A transdisciplinary analysis of water problems in the mountainous karst areas of Morocco

Water management and engineering in the karstic High Atlas of Morocco are difficult tasks under the prevailing geological, hydrogeological, geomorphological, vegetational and climatic conditions. It is important to be able to understand and predict the characteristics and availability of water for future water planning in the region under changing climatic and agricultural conditions. An interdisciplinary analysis of problems and adequate hydrological modelling tools developed by geologists, hydrologists and biologists are necessary. The karst areas of the High Atlas Mountains are characterised by impermeable triassic basalt underlying substantial subsurface reservoirs with high potential discharge rates. The karst groundwater aquifers are extensive but largely unknown in dimension, probably with a hierarchical network of groundwater flow paths. It is estimated that approximately 70% of the surface water is directly lost to groundwater. Steep landslide- and debris flow prone slopes exist next to coarse-grained, highly porous river beds. Infrequent, high intensity rainfall or snowmelt causes a particularly high flood risk to these karst areas. In addition, agriculture and land use changes have degraded the karst areas. The most important driving forces for degradation include permanent overgrazing even during droughts and the use of firewood by a continually growing population. Large scale degradation of vegetation has occurred in the oro-mediterranean (mountainous Mediterranean) zone, between 2600 and 3400 m which coincides with the most important zone for karstic groundwater creation. The combination of high amounts of groundwater flow and rapid surface flow due to sparse vegetation has increased the problems of flood flow.     

Assessing soil erosion in Mediterranean karst landscapes of Lebanon using remote sensing and GIS

Maps showing the potential for soil erosion at 1:100,000 scale are produced in a study area within Lebanon that can be used for evaluating erosion of Mediterranean karstic terrain with two different sets of impact factors built into an erosion model. The first set of factors is: soil erodibility, morphology, land cover/use and rainfall erosivity. The second is obtained by the first adding a fifth factor, rock infiltration. High infiltration can reflect high recharge, therefore decreasing the potential of surface runoff and hence the quantity of transported materials. Infiltration is derived as a function of lithology, lineament density, karstification and drainage density, all of which can be easily extracted from satellite imagery. The influence of these factors is assessed by a weight/rate approach sharing similarities between quantitative and qualitative methods and depending on pair-wise comparison matrix.The main outcome was the production of factorial maps and erosion susceptibility maps (scale 1:100,000). Spatial and attribute comparison of erosion maps indicates that the model that includes a measure of rock infiltration better represents erosion potential. Field investigation of rills and gullies shows 87.5% precision of the model with rock infiltration. This is 17.5% greater than the precision of the model without rock infiltration. These results indicate the necessity and importance of integrating information on infiltration of rock outcrops to assess soil erosion in Mediterranean karst landscapes.     

A quantitative interpretation of recent experimental results on stable carbon isotope fractionation by aerobic CH4-oxidizing bacteria

A quantitative model of recent laboratory experiments on carbon isotope fractionation by methane-oxidizing bacteria is proposed. The simulated experimental apparatus consists of a bacterial culture with a constant liquid volume, a gas headspace and a methane bubbling mechanism. The relative effects of bacterial growth and transport phenomena that do not depend on cell density are clarified. In all calculations, gas-liquid mass transfer is defined by unconstrained model parameters. Limited mass transfer from the culture into the headspace, rather than the incomplete dissolution of substrate-rich bubbles, seems to have caused an apparent decrease in the measured carbon isotope fractionation. The experimenters attributed this fractionation shift to a growing imbalance among kinetic rates as methane consumption by bacteria increases. Model predictions support this interpretation but also show that changes in carbon isotope fractionation in the course of the experiments cannot be unambiguously correlated with bacterial cell density unless gas-liquid mass transfer parameters are calibrated. Simulations of other laboratory experiments indicate that a reported change in carbon isotope fractionation could, in part at least, be the result of experimental conditions rather than the emergence of a different methane oxidation pathway postulated by the experimenters. A careful evaluation of mass transfer from the liquid culture into the gas headspace is warranted in this type of experiments since isotope fractionation factors are likely to be used in a wide variety of environmental contexts.     

Origin of the Oceanic Lithosphere

In a global examination of mid-ocean ridge basalt (MORB) glasscompositions, we find that Na₈–Fe₈–depth variationsdo not support modeling of MORBs as aggregates of melt compositionsgenerated over a large range of temperature and pressure. However,the Na₈–Fe₈ variations are consistent with the compositionalsystematics of solidus melts in the plagioclase–spinellherzolite transition in the CaO–MgO–Al₂O₃–SiO₂–Na₂O–FeO(CMASNF) system. For natural compositions, the P–T rangefor melt extraction is estimated to be 1·2–1·5GPa and 1250–1280°C. This P–T range is a closematch with the maximum P–T conditions for explosive pressure-releasevaporization of carbonate-bearing melts. It is proposed thatfracturing of the lithosphere induces explosive formation andescape of CO₂ vapor. This provides the vehicle for extractionof MORBs at a relatively uniform T and P. The upper portionof the CO₂-bearing and slightly melted seismic low-velocityzone flows toward the ridge, rises at the ridge axis to lower-lithospheredepths, melts much more extensively during this rise, and releasesMORB melts to the surface driven by explosively escaping CO₂vapor. The residue and overlying crust produced by this meltingthen migrate away from the ridge axis as new oceanic lithosphere.The entire process of oceanic lithosphere creation involvesonly the upper 140 km. When lithospheric stresses shift fractureformation to other localities, escape of CO₂ ceases, the vehiclefor transporting melt to the surface disappears, and ridgesdie. Inverse correlations of Na₈ vs Fe₈ for MORB glasses areexplained by mantle heterogeneity, and positive variations superimposedon the inverse variations are consistent with progressive extractionof melts from short, ascending melting columns. The uniformlylow temperatures of MORB extraction are not consistent withthe existence of hot plumes on or close to ocean ridges. Inthis modeling, the southern Atlantic mantle from Bouvet to about26°N is relatively homogeneous, whereas the Atlantic mantlenorth of about 26°N shows significant long-range heterogeneity.The mantle between the Charlie Gibbs and Jan Mayen fracturezones is strongly enriched in FeO/MgO, perhaps by a trappedfragment of basaltic crust. Iceland is explained as the productof this enrichment, not a hot plume. The East Pacific Rise,Galapagos Ridge, Gorda Ridge, and Juan de Fuca Ridge samplemantle that is heterogeneous over short distances. The mantlebeneath the Red Sea is enriched in FeO/MgO relative to the mantlebeneath the northern Indian Ocean.