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.     

Hydrological cycle and water balance estimates for the megadune–lake region of the Badain Jaran Desert,China

The hydrology and water balance of megadunes and lakes have been investigated in the Badain Jaran Desert of China. Field observations and analyses of sand layer water content, field capacity, secondary salt content, and grain size reveal 3 types of important natural phenomenon: (a) vegetation bands on the leeward slope of the megadunes reflect the hydrological regime within the sandy vadose zone; (b) seepage, wet sand deposits, and secondary salt deposits indicate the pattern of water movement within the sandy vadose zone; (c) zones of groundwater seeps and descending springs around the lakes reflect the influence of the local topography on the hydrological regime of the megadunes. The seepage exposed on the sloping surface of the megadunes and gravity water contained within the sand layer confirm the occurrence of preferential flow within the vadose zone of the megadunes. Alternating layers of coarse and fine sand create the conditions for the formation of preferential flows. The preferential flows promote movement of water within the sand layer water that leads to deep penetration of water within the megadunes and ultimately to the recharging of groundwater and lake water. Our results indicate that a positive water balance promotes recharge of the megadunes, which depends on the high permeability of the megadune material, the shallow depth of the surface sand layer affected by evaporation, the occurrence of rainfall events exceeding 15 mm, and the sparse vegetation cover. Water balance estimates indicate that the annual water storage of the megadunes is about 7.5 mm, accounting for only 8% of annual precipitation; however, the shallow groundwater per unit area under the megadunes receives only 3.6% of annual precipitation, but it is still able to maintain a dynamic balance of the lake water. From a water budget perspective, the annual water storage in the megadunes is sufficient to serve as a recharge source for lake water, thereby enabling the long‐term persistence of the lakes. Overall, our findings demonstrate that precipitation is a significant component of the hydrological cycle in arid deserts.     

Simulation of site-scale water fluxes in desert and natural oasis ecosystems of the arid region in Northwest China

The study of water fluxes is important to better understand hydrological cycles in arid regions. Data-driven machine learning models have been recently applied to water flux simulation. Previous studies have built site-scale simulation models of water fluxes for individual sites separately, requiring a large amount of data from each site and significant computation time. For arid areas, there is no consensus as to the optimal model and variable selection method to simulate water fluxes. Using data from seven flux observation sites in the arid region of Northwest China, this study compared the performance of random forest (RF), support vector machine (SVM), back propagation neural network (BPNN), and multiple linear regression (MLR) models in simulating water fluxes. Additionally, the study investigated inter-annual and seasonal variation in water fluxes and the dominant drivers of this variation at different sites. A universal simulation model for water flux was built using the RF approach and key variables as determined by MLR, incorporating data from all sites. Model performance of the SVM algorithm (R² = 0.25–0.90) was slightly worse than that of the RF algorithm (R² = 0.41–0.91); the BPNN algorithm performed poorly in most cases (R² = 0.15–0.88). Similarly, the MLR results were limited and unreliable (R² = 0.00–0.66). Using the universal RF model, annual water fluxes were found to be much higher than the precipitation received at each site, and natural oases showed higher fluxes than desert ecosystems. Water fluxes were highest during the growing season (May–September) and lowest during the non-growing season (October–April). Furthermore, the dominant drivers of water flux variation were various among different sites, but the normalized difference vegetation index (NDVI), soil moisture and soil temperature were important at most sites. This study provides useful insights for simulating water fluxes in desert and oasis ecosystems, understanding patterns of variation and the underlying mechanisms. Besides, these results can make a contribution as the decision-making basis to the water management in desert and oasis ecosystems.     

Changes in potential evapotranspiration and surface runoff in 1981–2010 and the driving factors in Upper Heihe River Basin in Northwest China

Changes in potential evapotranspiration and surface runoff can have profound implications for hydrological processes in arid and semiarid regions. In this study, we investigated the response of hydrological processes to climate change in Upper Heihe River Basin in Northwest China for the period from 1981 to 2010. We used agronomic, climatic and hydrological data to drive the Soil and Water Assessment Tool model for changes in potential evapotranspiration (ET₀) and surface runoff and the driving factors in the study area. The results showed that increasing autumn temperature increased snow melt, resulting in increased surface runoff, especially in September and October. The spatial distribution of annual runoff was different from that of seasonal runoff, with the highest runoff in Yeniugou River, followed by Babaohe River and then the tributaries in the northern of the basin. There was no evaporation paradox at annual and seasonal time scales, and annual ET₀ was driven mainly by wind speed. ET₀ was driven by relative humidity in spring, sunshine hour duration in autumn and both sunshine hour duration and relative humility in summer. Surface runoff was controlled by temperature in spring and winter and by precipitation in summer (flood season). Although surface runoff increased in autumn with increasing temperature, it depended on rainfall in September and on temperature in October and November. Copyright © 2016 John Wiley & Sons, Ltd.     

Influence of heterogeneity on unsaturated hydraulic properties: (1) local heterogeneity and scale effect

Spatial heterogeneity is ubiquitous in nature, which may significantly affect the soil hydraulic property curves. The models of a closed‐form functional relationship of soil hydraulic property curves (e.g. VG model or exponential model) are valid at point or local scale based on a point‐scale hydrological process, but how do scale effects of heterogeneity have an influence on the parameters of these models when the models are used in a larger scale process? This paper uses a two‐dimensional variably saturated flow and solute transport finite element model (VSAFT2) to simulate variations of pressure and moisture content in the soil flume under a constant head boundary condition. By changing different numerical simulation block sizes, a quantitative evaluation of parameter variations in the VG model, resulting from the scale effects, is presented. Results show that the parameters of soil hydraulic properties are independent of scale in homogeneous media. Parameters of α and n in homogeneous media, which are estimated by using the unsaturated hydraulic conductivity curve (UHC) or the soil water retention curve (WRC), are identical. Variations of local heterogeneities strongly affect the soil hydraulic properties, and the scale affects the results of the parameter estimations when numerical experiments are conducted. Furthermore, the discrepancy of each curve becomes considerable when moisture content becomes closer to a dry situation. Parameters estimated by UHC are totally different from the ones estimated by WRC. Copyright © 2012 John Wiley & Sons, Ltd.