Development of a rain‐on‐snow detection algorithm using passive microwave radiometry

Currently observed climate warming in the Arctic has numerous consequences. Of particular relevance, the precipitation regime is modified where mixed and liquid precipitation can occur during the winter season leading to rain‐on‐snow (ROS) events. This phenomenon is responsible for ice crust formation, which has a significant impact on ecosystems (such as biological, hydrological, ecological and physical processes). The spatially and temporally sporadic nature of ROS events makes the phenomenon difficult to monitor using meteorological observations. This paper focuses on the detection of ROS events using passive microwave (PMW) data from a modified brightness temperature (T_B) gradient approach at 19 and 37 GHz. The approach presented here was developed empirically for observed ROS events with coincident ground‐based PMW measurements in Sherbrooke, Quebec, Canada. It was then tested in Nunavik, Quebec, with the Advanced Microwave Scanning Radiometer for the Earth Observing System (AMSR‐E). We obtained a detection accuracy of 57, 71 and 89% for ROS detection for three AMSR‐E grid cells with a maximum error of 7% when considering all omissions and commissions with regard to the total number of AMSR‐E passes throughout the winter period. Copyright © 2016 John Wiley & Sons, Ltd.     

The exponential decline in saturated hydraulic conductivity with depth: a novel method for exploring its effect on water flow paths and transit time distribution

The strong vertical gradient in soil and subsoil saturated hydraulic conductivity is characteristic feature of the hydrology of catchments. Despite the potential importance of these strong gradients, they have proven difficult to model using robust physically based schemes. This has hampered the testing of hypotheses about the implications of such vertical gradients for subsurface flow paths, residence times and transit time distribution. Here we present a general semi‐analytical solution for the simulation of 2D steady‐state saturated‐unsaturated flow in hillslopes with saturated hydraulic conductivity that declines exponentially with depth. The grid‐free solution satisfies mass balance exactly over the entire saturated and unsaturated zones. The new method provides continuous solutions for head, flow and velocity in both saturated and unsaturated zones without any interpolation process as is common in discrete numerical schemes. This solution efficiently generates flow pathlines and transit time distributions in hillslopes with the assumption of depth‐varying saturated hydraulic conductivity. The model outputs reveal the pronounced effect that changing the strength of the exponential decline in saturated hydraulic conductivity has on the flow pathlines, residence time and transit time distribution. This new steady‐state model may be useful to others for posing hypotheses about how different depth functions for hydraulic conductivity influence catchment hydrological response. Copyright © 2016 John Wiley & Sons, Ltd.     

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.