An alternative Boussinesq equation considering the effect of hysteresis on coastal groundwater waves

This paper presents an alternative Boussinesq equation considering hysteresis effect via a third‐order derivative term. By introducing an improved moisture–pressure retention function, this equation describes, with reasonable precision, groundwater propagation in coastal aquifers subject to Dirichlet boundary condition of different oscillation frequencies. Test results confirmed that it is necessary to consider horizontal and vertical flows in unsaturated zone, because of their variable influences on hysteresis. Hysteresis in unsaturated zone can affect the water table wave number of groundwater wave motion, such as wave damping rate and phase lag. Oscillations with different periods exert different hysteresis effect on wave propagation. Truncation/shrinkage of unsaturated zones also affects the strength of hysteresis. These impacts can be reflected in the alternative Boussinesq equation by adjusting the parameter representing the variation rate of moisture associated with pressure change, as opposed to traditional computationally expensive hysteresis algorithms. The present Boussinesq equation is simple to use and can provide feasible basis for future coupling of groundwater and surface water models. Copyright © 2016 John Wiley & Sons, Ltd.     

Development of a soil moisture‐based distributed hydrologic model for determining hydrologically based critical source areas

A simple grid cell‐based distributed hydrologic model was developed to provide spatial information on hydrologic components for determining hydrologically based critical source areas. The model represents the critical process (soil moisture variation) to run‐off generation accounting for both local and global water balance. In this way, it simulates both infiltration excess run‐off and saturation excess run‐off. The model was tested by multisite and multivariable evaluation on the 50‐km² Little River Experimental Watershed I in Georgia, U.S. and 2 smaller nested subwatersheds. Water balance, hydrograph, and soil moisture were simulated and compared to observed data. For streamflow calibration, the daily Nash‐Sutcliffe coefficient was 0.78 at the watershed outlet and 0.56 and 0.75 at the 2 nested subwatersheds. For the validation period, the Nash‐Sutcliffe coefficients were 0.79 at the watershed outlet and 0.85 and 0.83 at the 2 subwatersheds. The per cent bias was less than 15% for all sites. For soil moisture, the model also predicted the rising and declining trends at 4 of the 5 measurement sites. The spatial distribution of surface run‐off simulated by the model was mainly controlled by local characteristics (precipitation, soil properties, and land cover) on dry days and by global watershed characteristics (relative position within the watershed and hydrologic connectivity) on wet days when saturation excess run‐off was simulated. The spatial details of run‐off generation and travel time along flow paths provided by the model are helpful for watershed managers to further identify critical source areas of non‐point source pollution and develop best management practices.     

Continuous measurement of whole‐tree water balance for studying urban tree transpiration

Street and garden trees in urban areas are often exposed to advection of strong vapour pressure deficit (VPD) air that can raise the whole‐tree transpiration rate (E_T), known as the oasis effect. However, urban trees tend to have small soil volume compared with natural conditions, and so they are believed to strongly regulate stomata. E_T characteristics of such urban trees have not been well understood because of a lack of reliable measurement methods. Therefore, we propose a novel weighing lysimeter method and investigate the whole‐tree water balance of an isolated container‐grown Zelkova serrata to examine (a) which biotic and abiotic factors determine E_T and (b) which spatial and temporal information is needed to predict E_T under urban conditions. Whole‐tree water balance and environmental conditions were measured from 2010 to 2012. Although leaf area substantially increased in the study period, daily E_T did not vary much. E_T increased with VPD almost linearly in 2010 but showed saturation in 2011 and 2012. Root water uptake lagged E_T by 40 min in 2012. These results suggest that the small planter box interfered with root growth and that hydraulic supply capacities did not increase sufficiently to support leaf area increase. From analysis of water balance, we believe that neglecting soil drought effects on street trees without irrigation in Japan will overestimate E_T over 4–5 sunny days at the longest. This is unlike previous studies of forest.     

Using hysteretic behaviour and hydrograph classification to identify hydrological function across the “hillslope–depression–stream” continuum in a karst catchment

In cockpit karst landscapes, fluxes from upland areas contribute large volumes of water to low-lying depressions and stream flow. Hydrograph hysteresis and similarity between monitoring sites is important for understanding the space–time variability of hydrologic responses across the “hillslope–depression–stream” continuum. In this study, the hysteretic feature of hydrographs was assessed by characterizing the loop-like relationships between responses at upstream sites relative to subsurface discharge at the outlet of a small karst catchment. A classification of hydrograph responses based on the multi-scale smoothing Kernel -derived distance classifies the hydrograph responses on the basis of similarities between hillslope and depression sites, and those at the catchment outlet. Results demonstrate that the temporal and spatial variability of hydrograph hysteresis and similarity between hillslope flow and outlet stream flow can be explained by the local heterogeneity of depression aquifer. Large depression storage deficits emerging in the highly heterogeneous aquifer produce strong hysteresis and multiple relationships of upstream hydrographs relative to the outlet subsurface discharge. In contrast, when depression storage deficits are filled during consecutive rainfall events, depression hydrographs at the high permeability sites are almost synchronous or exhibit a monotonous function with the hydrographs at the outlet. This reduced hydrograph hysteresis enhances preferential flow paths in fractured rocks and conduits that can accelerate the hillslope flow to the outlet. Therefore, classification of hydrograph similarities between any upstream sites and the catchment outlet can help to identify the dominant hydrological functions in the heterogeneous karst catchment.     

Piecewise prediction model for watershed‐scale erosion and sediment yield of individual rainfall events on the Loess Plateau,China

Establishing a universal watershed‐scale erosion and sediment yield prediction model represents a frontier field in erosion and soil/water conservation. The research presented here was conducted on the Chabagou watershed, which is located in the first sub‐region of the hill‐gully area of the Loess Plateau, China. A back‐propagation artificial neural model for watershed‐scale erosion and sediment yield was established, with the accuracy of the model, then compared with that of multiple linear regression. The sensitivity degree of various factors to erosion and sediment yield was quantitatively analysed using the default factor test. On the basis of the sensitive factors and the fractal information dimension, the piecewise prediction model for erosion and sediment yield of individual rainfall events was established and further verified. The results revealed the back‐propagation artificial neural network model to perform better than the multiple linear regression model in terms of predicting the erosion modulus, with the former able to effectively characterize dynamic changes in sediment yield under comprehensive factor conditions. The sensitivity of runoff erosion power and runoff depth to the erosion and sediment yield associated with individual rainfall events was found to be related to the complexity of surface topography. The characteristics of such a hydrological response are thus closely related to topography. When the fractal information dimension is greater than the topographic threshold, the accuracy of prediction using runoff erosion power is higher than that of using runoff depth. In contrast, when the fractal information dimension is smaller than the topographic threshold, the accuracy of prediction using runoff depth is higher than that of using runoff erosion power. The developed piecewise prediction model for watershed‐scale erosion and sediment yield of individual rainfall events, which introduces runoff erosion power and runoff depth using the fractal information dimension as a boundary, can be considered feasible and reliable and has a high prediction accuracy. Copyright © 2013 John Wiley & Sons, Ltd.     

A new numerical approach to the hyperstatic reaction method for segmental tunnel linings

This paper proposes a numerical approach to the hyperstatic reaction method (HRM) for the analysis of segmental tunnel linings. The influence of segmental joints has been considered directly using a fixity ratio that is determined on the basis of the rotational stiffness. The parameters necessary for the calculation are presented. A specific implementation has been developed using a FEM framework. This code is able to consider the three‐dimensional (3D) effect of segment joints in successive rings on the tunnel lining behaviour. The present HRM allows one to take an arbitrary distribution of segment joints along the tunnel boundary into consideration. In addition, the rotational stiffness of segment joints has been simulated using nonlinear behaviour, as it is closer to the true behaviour of a joint than linear or bilinear behaviour. The numerical results of three hypotheses on ring interaction, which allow the 3D effect of a segmental tunnel lining to be taken into account, have been compared with data obtained from the shield‐driven tunnel of the Bologna–Florence high‐speed railway line project. The numerical results presented in the paper show that the proposed HRM can be used to effectively estimate the behaviour of a segmental tunnel lining. Copyright © 2014 John Wiley & Sons, Ltd.     

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.     

Borehole failure analysis in a sandstone under anisotropic stresses

Borehole failure under anisotropic stresses in a sandstone is analyze numerically for various borehole sizes using a nonlinear elastic–plastic constitutive model for a Cosserat continuum. Borehole failure is identified as macroscopic failure of the borehole through the development of shear bands and breakouts. The results compare well both qualitatively and quantitatively with experimental results from polyaxial tests on Red Wildmoor sandstone. They show that the hole size effect of the borehole failure strength is independent of the far‐field stress anisotropy and follows a ? power law of the hole size. A similar scale effect equation with a ? power law is proposed for the scale effect of the maximum plastic shear strain at failure. This equation can be useful for better predicting hole‐size‐dependent failure with standard codes based on classical continua. The effect of stress anisotropy on the borehole failure stress is found to be independent of the hole size. The failure stress decreases linearly to 40% as the stress anisotropy increases. However, the maximum plastic shear strain at failure is stress anisotropy independent and therefore the critical plastic shear strain for failure is only hole‐size dependent. Copyright © 2009 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.