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.     

Mixing of event and pre‐event water in a shallow Entisol in sloping farmland based on isotopic and hydrometric measurements,SW China

Water percolation and flow processes in subsurface geologic media play an important role in determining the water source for plants and the transport of contaminants or nutrients, which is essential for water resource management and the development of measures for pollution mitigation. During June 2013, the dynamics of the rainwater, soil water, subsurface flows and groundwater in a shallow Entisol on sloping farmland were monitored using a hydrometric and isotopic approach. The results showed that effective mixing of rainwater and soil water occurred in hours. The rebound phenomenon of δD profiles in soils showed that most isotope‐depleted rainwater largely bypassed the soil matrix when the water saturation in the soil was high. Preferential‐flow, which was the dominant water movement pattern in the vadose zone, occurred through the whole soil profile, and infrequent piston‐flow was mainly found at 20–40 cm in depth. The interflow in the soil layer, composed of 75.2% rainwater, was only generated when the soil profile had been saturated. Underflow in the fractured mudrock was the dominant flow type in this hillslope, and outflow was dominated by base flow (groundwater flow) with a mean contribution of 76.7%. The generation mechanism of underflow was groundwater ridging, which was superimposed upon preferential‐flow composed mainly of rainwater. The quick mixing process of rainwater and soil water and the rapid movement of the mixture through preferential channels in the study soil, which shows a typical bimodal pore size distribution, can explain the prompt release of pre‐event water in subsurface flow. Water sources of subsurface flows at peak discharge could be affected by the antecedent soil water content, rain characteristics and antecedent groundwater levels. 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.     

A real-time inertial-aided cycle slip detection method based on ARIMA-GARCH model for inaccurate lever arm conditions

GPS Solutions - The GPS satellite transmitter antenna phase center offsets (PCOs) can be estimated in a global adjustment by constraining the ground station coordinates to the current International...     

Numerical modeling of toxic nonaqueous phase liquid removal from contaminated groundwater systems: mesh effect and discretization error estimation

Numerical modeling has now become an indispensable tool for investigating the fundamental mechanisms of toxic nonaqueous phase liquid (NAPL) removal from contaminated groundwater systems. Because the domain of a contaminated groundwater system may involve irregular shapes in geometry, it is necessary to use general quadrilateral elements, in which two neighbor sides are no longer perpendicular to each other. This can cause numerical errors on the computational simulation results due to mesh discretization effect. After the dimensionless governing equations of NAPL dissolution problems are briefly described, the propagation theory of the mesh discretization error associated with a NAPL dissolution system is first presented for a rectangular domain and then extended to a trapezoidal domain. This leads to the establishment of the finger‐amplitude growing theory that is associated with both the corner effect that takes place just at the entrance of the flow in a trapezoidal domain and the mesh discretization effect that occurs in the whole NAPL dissolution system of the trapezoidal domain. This theory can be used to make the approximate error estimation of the corresponding computational simulation results. The related theoretical analysis and numerical results have demonstrated the following: (1) both the corner effect and the mesh discretization effect can be quantitatively viewed as a kind of small perturbation, which can grow in unstable NAPL dissolution systems, so that they can have some considerable effects on the computational results of such systems; (2) the proposed finger‐amplitude growing theory associated with the corner effect at the entrance of a trapezoidal domain is useful for correctly explaining why the finger at either the top or bottom boundary grows much faster than that within the interior of the trapezoidal domain; (3) the proposed finger‐amplitude growing theory associated with the mesh discretization error in the NAPL dissolution system of a trapezoidal domain can be used for quantitatively assessing the correctness of computational simulations of NAPL dissolution front instability problems in trapezoidal domains, so that we can ensure that the computational simulation results are controlled by the physics of the NAPL dissolution system, rather than by the numerical artifacts. Copyright © 2014 John Wiley & Sons, Ltd.     

Quantifying the impact of landscape changes on hydrological variables in the alpine and cold region using hydrological model and remote sensing data

Quantifying the impact of landscape on hydrological variables is essential for the sustainable development of water resources. Understanding how landscape changes influence hydrological variables will greatly enhance the understanding of hydrological processes. Important vegetation parameters are considered in this study by using remote sensing data and VIC-CAS model to analyse the impact of landscape changes on hydrology in upper reaches of the Shule River Basin (URSLB). The results show there are differences in the runoff generation of landscape both in space and time. With increasing altitude, the runoff yields increased, with approximately 79.9% of the total runoff generated in the high mountains (4200–5900 m), and mainly consumed in the mid-low mountain region. Glacier landscape produced the largest runoff yields (24.9% of the total runoff), followed by low-coverage grassland (LG; 22.5%), alpine cold desert (AL; 19.6%), mid-coverage grassland (MG; 15.6%), bare land (12.5%), high-coverage grassland (HG; 4.5%) and shrubbery (0.4%). The relative capacity of runoff generation by landscapes, from high to low, was the glaciers, AL, LG, HG, MG, shrubbery and bare land. Furthermore, changes in landscapes cause hydrological variables changes, including evapotranspiration, runoff and baseflow. The study revealed that HG, MG, and bare land have a positive impact on evapotranspiration and a negative impact on runoff and baseflow, whereas AL and LG have a positive impact on runoff and baseflow and a negative impact on evapotranspiration. In contrast, glaciers have a positive impact on runoff. After the simulation in four vegetation scenarios, we concluded that the runoff regulation ability of grassland is greater than that of bare land. The grassland landscape is essential since it reduced the flood peak and conserved the soil and water.     

Varying water utilization of Haloxylon ammodendron plantations in a desert‐oasis ecotone

Haloxylon ammodendron is a desert shrub used extensively in China for restoring degraded dry lands. An understanding of the water source used by H. ammodendron plantations is critical achieving sustainable vegetation restoration. We measured mortality, shoot size, and rooting depth in 5‐, 10‐, 20‐, and 40‐year‐old H. ammodendron plantations. We examined stable isotopic ratios of oxygen (δ¹⁸O) in precipitation, groundwater, and soil water in different soil layers and seasons, and in plant stem water to determine water sources at different shrub ages. We found that water acquisition patterns in H. ammodendron plantations differed with plantation age and season. Thus, the main water source for 5‐year‐old shrubs was shallow soil water. Water sources of 10‐year‐old shrubs shifted depending on the soil water conditions during the season. Although their tap roots could absorb deep soil water, the plantation main water sources were from soil water, and about 50% of water originated from shallow and mid soil. This pattern might occur because main water sources in these plantations were changeable over time. The 20‐ and 40‐year‐old shrubs acquired water mainly from permanent groundwater. We conclude that the main water source of a young H. ammodendron plantation was soil water recharged by precipitation. However, when roots reached sufficient depth, water originated mainly from the deep soil water, especially in the dry season. The deeply rooted 20‐ and 40‐year‐old shrubs have the ability to exploit a deep and reliable water source. To achieve sustainability in these plantations, we recommend a reduction in the initial density of H. ammodendron in the desert‐oasis ecotone to decelerate the consumption of shallow soil water during plantation establishment.