The tectono‐metamorphic evolution of the very low‐grade hangingwall constrains two‐stage gneiss dome formation in the Montagne Noire (Southern France)

The Montagne Noire in the southernmost French Massif Central is made of an ENE‐elongated gneiss dome flanked by Palaeozoic sedimentary rocks. The tectonic evolution of the gneiss dome has generated controversy for more than half a century. As a result, a multitude of models have been proposed that invoke various tectonic regimes and exhumation mechanisms. Most of these models are based on data from the gneiss dome itself. Here, new constraints on the dome evolution are provided based on a combination of very low‐grade petrology, K–Ar geochronology, field mapping and structural analysis of the Palaeozoic western Mont Peyroux and Faugères units, which constitute part of the southern hangingwall of the dome. It is shown that southward‐directed Variscan nappe‐thrusting (D₁) and a related medium‐P metamorphism (M₁) are only preserved in the area furthest away from the gneiss dome. The regionally dominant pervasive tectono‐metamorphic event D₂/M₂ largely transposes D₁ structures, comprises a higher metamorphic thermal gradient than M₁ (transition low‐P and medium‐P metamorphic facies series) and affected the rocks between c. 309 and 300 Ma, post‐dating D₁/M₁ by more than 20 Ma. D₂‐related fabrics are refolded by D₃, which in its turn, is followed by dextral‐normal shearing along the basal shear zone of both units at c. 297 Ma. In the western Mont Peyroux and Faugères units, D₂/M₂ is largely synchronous with shearing along the southern dome margin between c. 311 and 303 Ma, facilitating the emplacement of the gneiss dome into the upper crust. D₂/M₂ also overlaps in time with granitic magmatism and migmatization in the Zone Axiale between c. 314 and 306 Ma, and a related low‐P/high‐T metamorphism at c. 308 Ma. The shearing that accompanied the exhumation of the dome therefore was synchronous with a peak in temperature expressed by migmatization and intrusion of melts within the dome, and also with the peak of metamorphism in the hangingwall. Both, the intensity of D₂ fabrics and the M₂ metamorphic grade within the hangingwall, decrease away from the gneiss dome, with grades ranging from the anchizone–epizone boundary to the diagenetic zone. The related zonation of the pre‐D₃ metamorphic field gradients paralleled the dome. These observations indicate that D₂/M₂ is controlled by the exhumation of the Zone Axiale, and suggest a coherent kinematic between the different crustal levels at some time during D₂/M₂. Based on integration of these findings with regional geological constraints, a two‐stage exhumation of the gneiss dome is proposed: during a first stage between c. 316 and 300 Ma dome emplacement into the upper crust was controlled by dextral shear zones arranged in a pull‐apart‐like geometry. The second stage from 300 Ma onwards was characterized by northeast to northward extension, with exhumation accommodated by north‐dipping detachments and hangingwall basin formation along the northeastern dome margin.     

Effects of differing coverage of moss‐dominated soil crusts on hydrological processes and implications for disturbance in the Mu Us Sandland,China

To study the effects of biological soil crusts (BSCs) on hydrological processes and their implications for disturbance in the Mu Us Sandland, the water infiltration, evaporation and soil moisture of high coverage (100% BSCs), middle coverage (40% BSCs) and low coverage (0% BSCs, bare sand) of moss‐dominated crusts were conducted in this study, respectively. The conclusions are as follows: (1) the main effects of moss‐dominated crusts in the Mu Us Sandland on the infiltration of rainwater were to reduce the infiltration depths and to retain the limited rainwater in shallow soil; (2) moss‐dominated crusts have no significant effects on daily evaporation when the volumetric water content at 4 cm depth in 100% BSCs (VWC4) was over 24.7%, on enhanced daily evaporation when the VWC4 ranged from 6.5% to 24.7% and on reduced daily evaporation when the VWC4 was less than 6.5%; and (3) decreasing the coverage of moss‐dominated crusts (from 100% to 40%) did not significantly change its effects on infiltration, evaporation and soil moisture. Our results demonstrated that for the growth and regeneration of shrubs, which were dominated by Artemisia ordosica in the Mu Us Sandland, high coverage of moss‐dominated crusts has negative effects on hydrological processes, and these negative effects could not be significantly reduced by decreasing the coverage of moss‐dominated crusts from 100% to 40%. Therefore, for the sustained and healthy development of shrub communities in the Mu Us Sandland, it is necessary to take appropriate measures for the well‐developed BSCs in the sites with high vegetation coverage in the rainy season. Copyright © 2015 John Wiley & Sons, Ltd.     

Inverse modelling using PS‐InSAR data for improved land subsidence simulation in Las Vegas Valley,Nevada

Las Vegas Valley has had a long history of groundwater development and subsequent surface deformation. InSAR interferograms have revealed detailed and complex spatial patterns of subsidence in the Las Vegas Valley area that do not coincide with major pumping regions. This research represents the first effort to use high spatial and temporal resolution subsidence observations from InSAR and hydraulic head data to inversely calibrate transmissivities (T), elastic and inelastic skeletal storage coefficients (S_ke and S_kv) of the developed‐zone aquifer and conductance (CR) of the basin‐fill faults for the entire Las Vegas basin. The results indicate that the subsidence observations from InSAR are extremely beneficial for accurately quantifying hydraulic parameters, and the model calibration results are far more accurate than when using only groundwater levels as observations, and just a limited number of subsidence observations. The discrepancy between distributions of pumping and greatest levels of subsidence is found to be attributed to spatial variations in clay thickness. The Eglington fault separates thicker interbeds to the northwest from thinner interbeds to the southeast and the fault may act as a groundwater‐flow barrier and/or subsidence boundary, although the influence of the groundwater barrier to this area is found to be insignificant. Copyright © 2016 John Wiley & Sons, Ltd.     

Linking tracers,water age and conceptual models to identify dominant runoff processes in a sparsely monitored humid tropical catchment

Assessing catchment runoff response remains a key research frontier because of limitations in current observational techniques to fully characterize water source areas and transit times in diverse geographical environments. Here, we report a study that combines empirical data with modelling to identify dominant runoff processes in a sparsely monitored humid tropical catchment. The analysis integrated isotope tracers into conceptual rainfall–runoff models of varying complexity (from 5 to 11 calibrated parameters) that are able to simulate discharge and tracer concentrations and track the evolving age of stream water exiting the catchment. The model structures can be seen as competing hypotheses of catchment functioning and were simultaneously calibrated against uncertain streamflow gaugings and a 2‐year daily isotope rainfall–runoff record. Comparison of the models was facilitated using global parameter sensitivity analysis and the resulting effect on calibration. We show that a variety of tested model structures reproduced water and tracer dynamics in stream, but the simpler models failed to adequately reproduce both. The resulting water age distributions of the tested models varied significantly with little similarity between the stream water age and stored water age distributions. The sensitivity analysis revealed that only some of the more complex models (from eight parameters) could be better constrained to infer more plausible water age distributions and catchment storage estimates. These models indicated that the age of water stored in the catchment is generally older compared with the age of water fluxes, with evapotranspiration age being younger compared with streamflow. However, the water age distributions followed a similar temporal behaviour dominated by climatic seasonality. Stream water ages increased during the dry season (greater than 1 year) and decreased with increased streamflow (a few weeks old) during the wet season. We further show that the ratios of the streamwater age to stored water age distribution and the water age distribution of actual evapotranspiration to the stored water age distribution from constrained models could potentially serve as useful hydrological indicators of catchment functioning. Copyright © 2016 John Wiley & Sons, Ltd.     

A two‐dimensional coupled hydromechanical discontinuum model for simulating rock hydraulic fracturing

Within the framework of our discontinuous deformation analysis for rock failure algorithm, this paper presents a two‐dimensional coupled hydromechanical discontinuum model for simulating the rock hydraulic fracturing process. In the proposed approach, based on the generated joint network, the calculation of fluid mechanics is performed first to obtain the seepage pressure near the tips of existing cracks, and then the fluid pressure is treated as linearly distributed loads on corresponding block boundaries. The contribution of the hydraulic pressure to the initiation/propagation of the cracks is considered by adding the components of these blocks into the force matrix of the global equilibrium equation. Finally, failure criteria are applied at the crack tips to determine the occurrence of cracking events. Several verification examples are simulated, and the results show that this newly proposed numerical model can simulate the hydraulic fracturing process correctly and effectively. Although the numerical and experimental verifications focus on one unique preexisting crack, because of the capability of discontinuous deformation analysis in simulating block‐like structures, the proposed approach is capable of modeling rock hydraulic fracturing processes. Copyright © 2014 John Wiley & Sons, Ltd.     

Dissolved inorganic carbon evolution in neutral discharge from mine tailings piles

We measured the concentrations of dissolved inorganic carbon (DIC) and major ions and the stable carbon isotope ratios of DIC (δ¹³C_DIC) in two creeks discharging from carbonate‐rich sulphide‐containing mine tailings piles. Our aim was to assess downstream carbon evolution of the tailings discharge as it interacted with the atmosphere. The discharge had pH of 6.5–8.1 and was saturated with respect to carbonates. Over the reach of one creek, the DIC concentrations decreased by 1.1 mmol C/l and δ¹³C_DIC increased by ~4.0‰ 200 m from the seep source. The decrease in the DIC concentrations was concomitant with decreases in the partial pressure of CO_2(aq) because of the loss of excess CO_2(aq) from the discharge. The corresponding enrichment in the δ¹³C_DIC is because of kinetic isotope fractionation accompanying the loss of CO_2(g). Over the reach of the other creek, there was no significant decrease in the DIC concentrations or notable changes in the δ¹³C_DIC. The insignificant change in the DIC concentrations and the δ¹³C_DIC is because the first water sample was collected 160 m away from the discharge seep, not accessible during this research. In this case, most of the excess CO_2(aq) was lost before our first sampling station. Our results indicate that neutral discharges from tailings piles quickly lose excess CO_2(aq) to the atmosphere and the DIC becomes enrich in ¹³C. We suggest that a significant amount of carbon cycling in neutral discharges from tailings piles occur close to the locations where the discharge seeps to the surface. Copyright © 2015 John Wiley & Sons, Ltd.     

Landform‐oriented flow‐routing algorithm for the dual‐structure loess terrain based on digital elevation models

The loess landform in the Loess Plateau of China is with typical dual structure, namely, the upper smooth positive terrain and the lower cliffy negative terrain (P–N terrain for short). Obvious differences in their morphological feature, geomorphological mechanism, and hydrological process could be found in the both areas. Based on the differences, a flow‐routing algorithm that separately addresses the dual‐structure terrain would be necessary to encompass this spatial variation in their hydrological behaviour. This paper proposes a mixed flow‐routing algorithm to address aforementioned problems. First, the loess landform surface is divided into P–N terrains based on digital elevation models. Then, specific catchment area is calculated with the new algorithm to simulate the water flows in both positive and negative terrain areas. The mixed algorithm consists of the multiple flow‐routing algorithm (multiple‐flow direction) for positive areas and the D8 algorithm for negative areas, respectively. The approach is validated in two typical geomorphologic areas with low hills and dense gullies in the northern Shaanxi Loess Plateau. Four indices are used to examine the results, which show that the new algorithm is more suitable for loess terrain in simulating the spatial distribution of water accumulation, as well as in modeling the flow characteristics of the true surface by considering the morphological structures of the terrain. Copyright © 2013 John Wiley & Sons, Ltd.     

Long‐term streamflow trends in Hawai'i and implications for native stream fauna

Climate change has fundamentally altered the water cycle in tropical islands, which is a critical driver of freshwater ecosystems. To examine how changes in streamflow regime have impacted habitat quality for native migratory aquatic species, we present a 50‐year (1967–2016) analysis of hydrologic records in 23 unregulated streams across the five largest Hawaiian Islands. For each stream, flow was separated into direct run‐off and baseflow and high‐ and low‐flow statistics (i.e., Q10 and Q90) with ecologically important hydrologic indices (e.g., frequency of flooding and low flow duration) derived. Using Mann–Kendall tests with a running trend analysis, we determined the persistence of streamflow trends through time. We analysed native stream fauna from ~400 sites, sampled from 1992 to 2007, to assess species richness among islands and streams. Declines in streamflow metrics indicated a general drying across the islands. In particular, significant declines in low flow conditions (baseflows), were experienced in 57% of streams, compared with a significant decline in storm flow conditions for 22% of streams. The running trend analysis indicated that many of the significant downward trends were not persistent through time but were only significant if recent decades (1987–2016) were included, with an average decline in baseflow and run‐off of 10.90% and 8.28% per decade, respectively. Streams that supported higher native species diversity were associated with moderate discharge and baseflow index, short duration of low flows, and negligible downward trends in flow. A significant decline in dry season flows (May–October) has led to an increase in the number of no‐flow days in drier areas, indicating that more streams may become intermittent, which has important implications for mauka to makai (mountain to ocean) hydrological connectivity and management of Hawai'i's native migratory freshwater fauna.     

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.