An Eocene/Oligocene blueschist‐/greenschist facies P–T loop from the Cycladic Blueschist Unit on Naxos Island,Greece: Deformation‐related re‐equilibration vs. thermal relaxation

Geothermobarometric and geochronological work indicates a complete Eocene/early Oligocene blueschist/greenschist facies metamorphic cycle of the Cycladic Blueschist Unit on Naxos Island in the Aegean Sea region. Using the average pressure–temperature (P–T) method of thermocalc coupled with detailed textural work, we separate an early blueschist facies event at 576 ± 16 to 619 ± 32°C and 15.5 ± 0.5 to 16.3 ± 0.9 kbar from a subsequent greenschist facies overprint at 384 ± 30°C and 3.8 ± 1.1 kbar. Multi‐mineral Rb–Sr isochron dating yields crystallization ages for near peak‐pressure blueschist facies assemblages between 40.5 ± 1.0 and 38.3 ± 0.5 Ma. The greenschist facies overprint commonly did not result in complete resetting of age signatures. Maximum ages for the end of greenschist facies reworking, obtained from disequilibrium patterns, cluster near c. 32 Ma, with one sample showing rejuvenation at c. 27 Ma. We conclude that the high‐P rocks from south Naxos were exhumed to upper mid‐crustal levels in the late Eocene and early Oligocene at rates of 7.4 ± 4.6 km/Ma, completing a full blueschist‐/greenschist facies metamorphic cycle soon after subduction within c. 8 Ma. The greenschist facies overprint of the blueschist facies rocks from south Naxos resulted from rapid exhumation and associated deformation/fluid‐controlled metamorphic re‐equilibration, and is unrelated to the strong high‐T metamorphism associated with the Miocene formation of the Naxos migmatite dome. It follows that the Miocene thermal overprint had no impact on rock textures or Sr isotopic signatures, and that the rocks of south Naxos underwent three metamorphic events, one more than hitherto envisaged.     

Computational model coupling mode II discrete fracture propagation with continuum damage zone evolution

We propose a numerical method that couples a cohesive zone model (CZM) and a finite element‐based continuum damage mechanics (CDM) model. The CZM represents a mode II macro‐fracture, and CDM finite elements (FE) represent the damage zone of the CZM. The coupled CZM/CDM model can capture the flow of energy that takes place between the bulk material that forms the matrix and the macroscopic fracture surfaces. The CDM model, which does not account for micro‐crack interaction, is calibrated against triaxial compression tests performed on Bakken shale, so as to reproduce the stress/strain curve before the failure peak. Based on a comparison with Kachanov's micro‐mechanical model, we confirm that the critical micro‐crack density value equal to 0.3 reflects the point at which crack interaction cannot be neglected. The CZM is assigned a pure mode II cohesive law that accounts for the dependence of the shear strength and energy release rate on confining pressure. The cohesive shear strength of the CZM is calibrated by calculating the shear stress necessary to reach a CDM damage of 0.3 during a direct shear test. We find that the shear cohesive strength of the CZM depends linearly on the confining pressure. Triaxial compression tests are simulated, in which the shale sample is modeled as an FE CDM continuum that contains a predefined thin cohesive zone representing the idealized shear fracture plane. The shear energy release rate of the CZM is fitted in order to match to the post‐peak stress/strain curves obtained during experimental tests performed on Bakken shale. We find that the energy release rate depends linearly on the shear cohesive strength. We then use the calibrated shale rheology to simulate the propagation of a meter‐scale mode II fracture. Under low confining pressure, the macroscopic crack (CZM) and its damaged zone (CDM) propagate simultaneously (i.e., during the same loading increments). Under high confining pressure, the fracture propagates in slip‐friction, that is, the debonding of the cohesive zone alternates with the propagation of continuum damage. The computational method is applicable to a range of geological injection problems including hydraulic fracturing and fluid storage and should be further enhanced by the addition of mode I and mixed mode (I+II+III) propagation. Copyright © 2016 John Wiley & Sons, Ltd.     

Macropores as preferential flow paths in meander bends

Macropores are subsurface connected void spaces caused by processes such as fracture of soils, micro‐erosion, and fauna burrows. They are common near streams (e.g. hyporheic and riparian zones) and may act as preferential flow paths between surface and groundwaters, affecting hydrologic and biogeochemical processes. We tested the hydrologic function of macropores by constructing an artificial macropore within the saturated zone of a meander bend (open macropore, ‘OM’) and later filling its upstream end (partially filled macropore, ‘PFM’). For each treatment, we injected saline tracer at an upgradient monitoring well within the meander and monitored downgradient hydraulics and tracer transport. Pressure transducers in monitoring wells indicated hydraulic gradients within the meander were 32% higher perpendicular to and 6% higher parallel to the macropore for the OM than for the PFM. Additionally, hydraulic conductivities measured via falling head tests were 29 to 550 times higher along the macropore than in nearby sediment. We used electrical conductivity probes in wells and electrical resistivity imaging to track solute transport. Transport velocities through the meander were on average 9 and 21% higher (per temporal moment analysis and observed tracer peak, respectively) for the OM than for the PFM. Furthermore, temporal moments of tracer breakthrough analysis indicated downgradient longitudinal dispersion and breakthrough tracer curve tailing were on average 234% and 182% higher for the OM, respectively. This suggests the OM enabled solute transport at overall shorter timescales than the matrix but also increased tailing. Our results demonstrate the importance of macropores to meander bend hydrology and solute transport. Copyright © 2012 John Wiley & Sons, Ltd.     

Modelling and experimental study of coupled exchange of p,p′‐DDE and colloids in a flume simulated stream–streambed system

Stream–subsurface exchange plays a significant role in the fate and transport of contaminants in streams. It has been modelled explicitly by considering fundamental processes such as hydraulic exchange, colloid filtration, and contaminant interactions with streambed sediments and colloids. The models have been successfully applied to simulate the transport of inorganic metals and nutrients. In this study, laboratory experiments were conducted in a recirculating flume to investigate the exchange of a hydrophobic organic contaminant, p,p′‐dichloro‐diphenyl‐dichloroethane (DDE), between a stream and a quartz sand bed. A previously developed process‐based multiphase exchange model was modified by accounting for the p,p′‐DDE kinetic adsorption to and desorption from the bed sediments/colloids and was applied to interpret the experimental results. Model input parameters were obtained by conducting independent small‐scale batch experiments. Results indicate that the immobilization of p,p′‐DDE in the quartz sand bed can occur under representative natural stream conditions. The observed p,p′‐DDE exchange was successfully simulated by the process‐based model. The model sensitivity analysis results show that the exchange of p,p′‐DDE can be sensitive to either the sediment sorption/desorption parameters or colloidal parameters depending on the experimental conditions tested. For the experimental conditions employed here, the effect of colloids on contaminant transport is expected to be minimal, and the stream–subsurface exchange of p,p′‐DDE is dominated by the interaction of p,p′‐DDE with bed sediment. The work presented here contributes to a better mechanistic understanding of the complex transport process that hydrophobic organic contaminants undergo in natural streams and to the development of reliable, predictive models for the assessment of impacted streams. Copyright © 2015 John Wiley & Sons, Ltd.     

Water uptake of riparian plants in the lower Lhasa River Basin,South Tibetan Plateau using stable water isotopes

Riparian plants can adapt their water uptake strategies based on climatic and hydrological conditions within a river basin. The response of cold-alpine riparian trees to changes in water availability is poorly understood. The Lhasa River is a representative cold-alpine river in South Tibet and an under-studied environment. Therefore, a 96 km section of the lower Lhasa River was selected for a study on the water-use patterns of riparian plants. Plant water, soil water, groundwater and river water were measured at three sites for δ¹⁸O and δ²H values during the warm-wet and cold-dry periods in 2018. Soil profiles differed in isotope values between seasons and with the distance along the river. During the cold-dry period, the upper parts of the soil profiles were significantly affected by evaporation. During the warm-wet period, the soil profile was influenced by precipitation infiltration in the upper reaches of the study area and by various water sources in the lower reaches. Calculations using the IsoSource model indicated that the mature salix and birch trees (Salix cheilophila Schneid. and Betula platyphylla Suk.) accessed water from multiple sources during the cold-dry period, whereas they sourced more than 70% of their requirement from the upper 60–80 cm of the soil profile during the warm-wet period. The model indicated that the immature rose willow tree (Tamarix ramosissima Ledeb) accessed 66% of its water from the surface soil during the cold-dry period, but used the deeper layers during the warm-wet period. The plant type was not the dominant factor driving water uptake patterns in mature plants. Our findings can contribute to strategies for the sustainable development of cold-alpine riparian ecosystems. It is recommended that reducing plantation density and collocating plants with different rooting depths would be conducive to optimal plant growth in this environment.     

Numerical and theoretical research on stress distribution in the loosening zone of the trapdoor problem

The traditional theory of soil arching effect was developed on the assumption that stress distribution in the loosening zone is uniform. However, because of the deflection of principal stress' direction, the stress distribution in the loosening zone is actually ununiform. For the evaluation of principal stress axis deflection and stress redistribution, a discrete element method numerical model of trapdoor problem is established for the simulation of soil arching effect. Based on the numerical results, an arc shape of major principal stress trajectory and uniform horizontal stress distribution at the same depth of the loosening zone are adopted. An analytical model is raised to estimate the average loosening earth pressure acting on the trapdoor and stress distribution in the loosening zone at a limit state. In addition, comparison studies are carried out between the predictions of the proposed solutions and discrete element method numerical results as well as available model test results, thereby validating the accuracy of the proposed theoretical model. Both numerical and theoretical results indicate that the vertical stress distribution in the loosening zone is obviously ununiform. The load acting in the middle of loosening zone is transferred toward two sides so that the vertical stress distribution in loosening zone is concave.     

Water mining from the deep critical zone by apple trees growing on loess

There have been significant recent advances in understanding the ecohydrology of deep soil. However, the links between root development and water usage in the deep critical zone remains poorly understood. To clarify the interaction between water use and root development in deep soil, we investigated soil water and root profiles beyond maximum rooting depth in five apple orchards planted on farmland with stand ages of 8, 11, 15, 18, and 22 years in a subhumid region on the Chinese Loess Plateau. Apple trees rooted progressively deeper for water with increasing stand age and reached 23.2 ± 0.8 m for the 22‐year‐old trees. Soil water deficit in deep soil increased with tree age and was 1,530 ± 43 mm for a stand age of 22 years. Measured root deepening rate was far great than the reported pore water velocity, which demonstrated that trees are mining resident old water. The deficits are not replenished during the life‐span of the orchard, showing a one‐way mining of the critical zone water. The one‐way root water mining may have changed the fine root profile from an exponential pattern in the 8‐year‐old orchard to a relative uniform distribution in older orchards. Our findings enhance our understanding of water‐root interaction in deep soil and reveal the unintended consequences of critical zone dewatering during the lifespan of apple trees.     

Root‐zone moisture replenishment in a native vegetated catchment under Mediterranean climate

The root‐zone moisture replenishment mechanisms are key unknowns required to understand soil hydrological processes and water sources used by plants. Temporal patterns of root‐zone moisture replenishment reflect wetting events that contribute to plant growth and survival and to catchment water yield. In this study, stable oxygen and hydrogen isotopes of twigs and throughfall were continuously monitored to characterize the seasonal variations of the root‐zone moisture replenishment in a native vegetated catchment under Mediterranean climate in South Australia. The two studied hillslopes (the north‐facing slope [NFS] and the south‐facing slope [SFS]) had different environmental conditions with opposite aspects. The twig and throughfall samples were collected every ~20 days over 1 year on both hillslopes. The root‐zone moisture replenishment, defined as percentage of newly replenished root‐zone moisture as a complement to antecedent moisture for plant use, calculated by an isotope balance model, was about zero (±25% for the NFS and ± 15% for the SFS) at the end of the wet season (October), increased to almost 100% (±26% for the NFS and ± 29% for the SFS) after the dry season (April and May), then decreased close to zero (±24% for the NFS and ± 28% for the SFS) in the middle of the following wet season (August). This seasonal pattern of root‐zone moisture replenishment suggests that the very first rainfall events of the wet season were significant for soil moisture replenishment and supported the plants over wet and subsequent dry seasons, and that NFS completed replenishment over a longer time than SFS in the wet season and depleted the root zone moisture quicker in the dry season. The stable oxygen isotope composition of the intraevent samples and twigs further confirms that rain water in the late wet season contributed little to root‐zone moisture. This study highlights the significant role of the very first rain events in the early wet season for ecosystem and provides insights to understanding ecohydrological separation, catchment water yield, and vegetation response to climate changes.     

New P–T constraints on the Tamayen glaucophane‐bearing rocks,eastern Taiwan: Perple_X modelling results and geodynamic implications

New pseudosection modelling was applied to better constrain the P–T conditions and evolution of glaucophane‐bearing rocks in the Tamayen block of the Yuli belt, recognized as the world's youngest known blueschist complex. Based on the predominant clinoamphibole, textural relationships, and mineral compositions, these glaucophane‐bearing high‐P rocks can be divided into four types. We focused on the three containing garnet. The chief phase assemblages are (in decreasing mode): amphibole + quartz + epidote + garnet + chlorite + rutile/titanite (Type‐I), phengite + amphibole + quartz + garnet + chlorite + epidote + titanite + biotite + magnetite (Type‐II), and amphibole + quartz + albite + epidote + garnet + rutile + hematite + titanite (Type‐III). Amphibole exhibits compositional zoning from core to rim as follows: glaucophane → pargasitic amphibole → actinolite (Type‐I), barroisite → Mg‐katophorite/taramite → Fe‐glaucophane (Type‐II), glaucophane → winchite (Type‐III). Using petrographic data, mineral compositions and Perple_X modelling (pseudosections and superimposed isopleths), peak P–T conditions were determined as 13 ± 1 kbar and 550 ± 40 °C for Type‐I, 10.5 ± 0.5 kbar and 560 ± 30 °C for Type‐II (thermal peak) and 11 ± 1 kbar and 530 ± 30 °C for Type‐III. The calculations yield higher pressures and temperatures than previously thought; the difference is ~1–6 kbar and 50–200 °C. The three rock types record similar P–T retrograde paths with clockwise trajectories; all rocks followed trajectories with substantial pressure decrease under near‐isothermal conditions (Type‐I and Type‐III), with the probable exception of Type‐II where decompression followed colder geotherms. The P–T paths suggest a tectonic environment in which the rocks were exhumed from maximum depths of ~45 km within a subduction channel along a relative cold geothermal gradient of ~11–14 °C km^?1.