We study the contribution of typically uncertain subsurface flow parameters to gravity changes that can be recorded during pumping tests in unconfined aquifers. We do so in the framework of a Global ...Sensitivity Analysis and quantify the effects of uncertainty of such parameters on the first four statistical moments of the probability distribution of gravimetric variations induced by the operation of the well. System parameters are grouped into two main categories, respectively, governing groundwater flow in the unsaturated and saturated portions of the domain. We ground our work on the three‐dimensional analytical model proposed by Mishra and Neuman (2011), which fully takes into account the richness of the physical process taking place across the unsaturated and saturated zones and storage effects in a finite radius pumping well. The relative influence of model parameter uncertainties on drawdown, moisture content, and gravity changes are quantified through (a) the Sobol' indices, derived from a classical decomposition of variance and (b) recently developed indices quantifying the relative contribution of each uncertain model parameter to the (ensemble) mean, skewness, and kurtosis of the model output. Our results document (i) the importance of the effects of the parameters governing the unsaturated flow dynamics on the mean and variance of local drawdown and gravity changes; (ii) the marked sensitivity (as expressed in terms of the statistical moments analyzed) of gravity changes to the employed water retention curve model parameter, specific yield, and storage, and (iii) the influential role of hydraulic conductivity of the unsaturated and saturated zones to the skewness and kurtosis of gravimetric variation distributions. The observed temporal dynamics of the strength of the relative contribution of system parameters to gravimetric variations suggest that gravity data have a clear potential to provide useful information for estimating the key hydraulic parameters of the system.
Key Points
Effect of hydrogeological parameter uncertainty on drawdown, water content and gravity changes during pumping tests in unconfined aquifers
The strength of the relative contribution of saturated and unsaturated zone parameters to gravimetric variations markedly varies over time
Gravimetric information are mostly sensitive to specific yield and aquifer specific storage, especially at early pumping times
Full text
Available for:
BFBNIB, FZAB, GIS, IJS, KILJ, NLZOH, NUK, OILJ, SAZU, SBCE, SBMB, UL, UM, UPUK
In recent years, wildfires in the western United States have occurred with increasing frequency and scale. Climate change scenarios in California predict prolonged periods of droughts with even ...greater potential for conditions amenable to wildfires. The Sierra Nevada Mountains provide 70% of water resources in California, yet how wildfires will impact watershed‐scale hydrology is highly uncertain. In this work, we assess the impacts of wildfires perturbations on watershed hydrodynamics using a physically based integrated hydrologic model in a high‐performance‐computing framework. A representative Californian watershed, the Cosumnes River, is used to demonstrate how postwildfire conditions impact the water and energy balance. Results from the high‐resolution model show counterintuitive feedbacks that occur following a wildfire and allow us to identify the regions most sensitive to wildfires conditions, as well as the hydrologic processes that are most affected. For example, whereas evapotranspiration generally decreases in the postfire simulations, some regions experience an increase due to changes in surface water run‐off patterns in and near burn scars. Postfire conditions also yield greater winter snowpack and subsequently greater summer run‐off as well as groundwater storage in the postfire simulations. Comparisons between dry and wet water years show that climate is the main factor controlling the timing at which some hydrologic processes occur (such as snow accumulation) whereas postwildfire changes to other metrics (such as streamflow) show seasonally dependent impacts primarily due to the timing of snowmelt, illustrative of the integrative nature of hydrologic processes across the Sierra Nevada‐Central Valley interface.
An integrated hydrologic model is used to simulate watershed hydrodynamics following land cover changes due to a wildfire. Differences between present‐day and postwildfire groundwater pressure show nonlinear increases and decreases that are not spatially limited to burn scar areas.
Full text
Available for:
BFBNIB, FZAB, GIS, IJS, KILJ, NLZOH, NUK, OILJ, SAZU, SBCE, SBMB, UL, UM, UPUK
Water resources are impacted by water‐energy balance fluxes at the land surface, most notably evapotranspiration (ET), the largest component of surface energy balance. While integrated hydrologic ...models show promise in quantifying the nonlinear dynamics at this interface, the model results are plagued by parametric uncertainties. Given the high computational demand of running multiple parameter spaces of these models, little is known about how these uncertainties propagate into land surface processes. In this work, we perform a global sensitivity analysis by computing Sobol and AMAE indices using a surrogate model constructed with a polynomial chaos expansion to assess the impacts of the subsurface physical properties on ET. We do so by modeling a semisynthetic test case. Our results show that the effects of vertical hydraulic conductivity, porosity, and the water retention curve parameter Van Genuchten α mainly control ET. However, we note that while evaporation (E) shows behavior similar to the ET with high sensitivity to the parameters controlling the flow in the unsaturated zone, transpiration (T) is very sensitive to the saturated zone parameters and groundwater flow, especially during periods without rain. Our results also show that heterogeneities in land cover impact the flow in the saturated zone. Furthermore, our work demonstrates that reduced‐order models can be developed to make integrated models more accessible for rigorous sensitivity analyses and calibration purposes.
Key Points
A global sensitivity analysis is performed to assess the connectivity between land surface and subsurface processes
Evapotranspiration is mainly controlled by the effects of vadose zone parameters while transpiration shows sensitivity to saturated zone parameters
Deeper water tables show a disconnection between land and subsurface, although not an insensitivity to hydrodynamic parameters
Full text
Available for:
BFBNIB, FZAB, GIS, IJS, KILJ, NLZOH, NUK, OILJ, SAZU, SBCE, SBMB, UL, UM, UPUK
High-Mountain Asia (HMA) exhibits one of the highest increases in vegetation greenness on Earth, subsequently influencing the exchange of water and energy between the land surface and the atmosphere. ...Given the strong interactions between the hydrosphere, the biosphere, and the cryosphere, understanding the drivers of greening in this highly complex region with significant land cover heterogeneity is essential to assess the changes in the regional water budget. Here we perform a holistic multivariate remote sensing analysis to simultaneously examine the primary components of the terrestrial water cycle from 2003 to 2020 and decipher the principal drivers of greening in HMA. We identified three drivers of greening: (1) precipitation drives greening in mid and low elevation areas covered by evergreen and mixed forests (e.g., Irrawaddy basin), (2) decreases in snow enhance greening in most of the hydrologic basins, and (3) irrigation induces greening in irrigated lands (Ganges-Brahmaputra and Indus).
Located in northwestern Niger, the Tim Mersoï Basin (TMB) is an important mining region in the scale of West Africa. Groundwater is considered the main source of fresh water in the basin, especially ...for mining activities. It, therefore, appears essential to monitor their responses to these activities. However, no study has been carried out in the Tim Mersoï Basin. This study aims to evaluate the groundwater storage changes (GWSC) of the TMB and to analyze the spatio-temporal evolution of the Tarat aquifer under the effect of mining activities in the Arlit region. For this purpose, Gravity Recovery And Climate Experiment (GRACE), Global Land Data Assimilation System (GLDAS), and in-situ data were used. The results show a variation of the GWS from 2002 to 2019 of about −0.1310 cm/year on the scale of the basin and −0.0109 cm/year in the Arlit mining area. The GWSC at the basin scale and the one at the Arlit region scale were shown to be linked with an RMSE between the two datasets of 0.79. This shows the potential of GRACE for contextualizing studies in small areas. The study also highlighted that the groundwater flow direction was highly modified; the drawdown of the Tarat water table was more than 50 m in the areas heavily impacted by mining activities, with an increasing intensity from the northwest to the southeast of Arlit.
•Analyses of gravity changes while pumping in 3D randomly heterogeneous aquifers.•Probability distributions of head and gravimetric variations during well operation.•Randomly heterogeneous domains ...consider uncertain geostatistical system descriptors.•The variance of heterogeneous fields affects mean and variance of gravity changes.•The correlation scale of the field affects extreme values of gravity changes.
We evaluate the relative importance of the uncertainty related to parameters characterizing partially saturated groundwater flow on head and gravity changes associated with pumping tests taking place in homogeneous and heterogeneous porous media. We frame our study in a Global Sensitivity Analysis setting and assess the way imperfect knowledge of such parameters influences the probability distribution (pdf) of head and gravimetric variations recorded during well operation. We rely on a set of detailed computational analyses and conceptualize uncertain model parameters as random quantities. Randomly heterogeneous domains are treated by considering main geostatistical descriptors (i.e., variance and correlation scale) of three-dimensional spatial distributions of system properties as affected by uncertainty. We quantify the effects of the latter on the resulting pdf of (ensemble) mean and variance of head and gravity changes through a numerical Monte Carlo approach. While all uncertain parameters are influential to gravity changes in both homogeneous and randomly heterogeneous scenarios, consistent with the integral nature of gravity observations, our study enables us to quantify their relative importance. Values of Ensemble mean and variance of head and gravity changes associated with randomly heterogeneous fields are generally more influenced by the variance rather than by the correlation scale of the spatially heterogeneous parameters considered. Uncertainty in the correlation scale is more influential to the shape, and hence on extreme values, of the probability distribution of these moments.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NLZOH, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UILJ, UL, UM, UPCLJ, UPUK, ZAGLJ, ZRSKP
In recent years, wildfires in the western United States have occurred with increasing frequency and scale. Climate change scenarios in California predict prolonged periods of droughts with even ...greater potential for conditions amenable to wildfires. The Sierra Nevada Mountains provide 70% of water resources in California, yet how wildfires will impact watershed-scale hydrology is highly uncertain. In this work, we assess the impacts of wildfires perturbations on watershed hydrodynamics using a physically based integrated hydrologic model in a high-performance-computing framework. A representative Californian watershed, the Cosumnes River, is used to demonstrate how postwildfire conditions impact the water and energy balance. Results from the high-resolution model show counterintuitive feedbacks that occur following a wildfire and allow us to identify the regions most sensitive to wildfires conditions, as well as the hydrologic processes that are most affected. For example, whereas evapotranspiration generally decreases in the postfire simulations, some regions experience an increase due to changes in surface water run-off patterns in and near burn scars. Postfire conditions also yield greater winter snowpack and subsequently greater summer run-off as well as groundwater storage in the postfire simulations. Finally, comparisons between dry and wet water years show that climate is the main factor controlling the timing at which some hydrologic processes occur (such as snow accumulation) whereas postwildfire changes to other metrics (such as streamflow) show seasonally dependent impacts primarily due to the timing of snowmelt, illustrative of the integrative nature of hydrologic processes across the Sierra Nevada-Central Valley interface.
Full text
Available for:
BFBNIB, FZAB, GIS, IJS, KILJ, NLZOH, NUK, OILJ, SAZU, SBCE, SBMB, UL, UM, UPUK
•Predictor-corrector (P-C) is a robust method to solve the 1D Richards equation.•Mass is conserved for the P-C method with the proposed post-allocation procedure.•Courant number can be used to ...adaptively adjust time step size.•Different time stepping strategies are compared on various test scenarios.
The predictor-corrector-type (P-C) numerical solution to the 1D Richards equation only requires one matrix inversion operation per time step, making it attractive in terms of computational cost. However, the mass conservation could be violated at the saturated-unsaturated interface. A new post-allocation procedure is designed for the P-C method, which redistributes moisture after the corrector step to achieve strict mass balance. A novel adaptive time-stepping strategy is proposed to further improve model efficiency and robustness. It adjusts time step size based on both moisture difference and the Courant number. The proposed solution method and time control strategies are tested and compared with an analytical solution, the previous P-C solution and other existing iterative solutions. The new method shows good conservation property and good agreements to the existing solutions. Compared to the iterative methods that occasionally experience convergence issues, the proposed P-C method is more robust. The new time-control strategy improves computational efficiency compared to the original P-C method, but it remains less efficient than iterative methods for most of the tested scenarios because of its explicit treatment of the corrector step.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NLZOH, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UILJ, UL, UM, UPUK, ZAGLJ, ZRSKP
Les eaux souterraines constituent une réserve d’eau potable non négligeable, leur alimentation se fait en majeure partie par les précipitations, appelée recharge des nappes. Du fait de leur grande ...importance, la compréhension du fonctionnement de ces ressources en eau est plus que jamais indispensable. Celle-ci passe par l’élaboration de modèles mathématiques. Ces outils nous offrent une meilleure appréhension et une bonne prévision des phénomènes physiques. Les systèmes hydrogéologiques sont généralement très complexes et caractérisés par des dynamiques hydriques très variables dans le temps et dans l’espace. Cette complexité a attiré l’attention de nombreux hydrogéologues et un grand nombre de modèles très sophistiqués a été développé afin de décrire ces systèmes avec précision. Cependant, la prise en compte de la recharge de ces réservoirs reste toujours un défi dans la modélisation hydrogéologique. En effet, le plus souvent, les modèles hydrogéologiques simulent l’écoulement dans la nappe tout en considérant la recharge comme une constante sur le domaine et indépendante du système. De plus, elle est souvent calculée de façon simplifiée. Or, la recharge traduisant la quantité des précipitations atteignant les nappes est une composante hydrologique complexe et variable car elle interagit avec les nappes et dépend des conditions climatiques, du couvert végétal et du transfert de l’eau dans le sol. Ce présent travail vise à intégrer cette recharge variable et complexe aux modèles hydrogéologiques. À cet effet, un modèle couplé a été développé. Une première partie de ce modèle permet de calculer la recharge des nappes en modélisant les interactions précipitations-sol et l’hydrodynamique dans le sol. Cette modélisation a été effectuée en utilisant des modèles conceptuels simples basés sur des lois empiriques (Gardénia, Nash) et des modèles physiques résolvant l’équation de Richards. La recharge ainsi calculée est intégrée à la deuxième partie du modèle simulant l’hydrodynamique dans les nappes décrite par l’équation de diffusivité. Des méthodes numériques précises et robustes ont été utilisées pour résoudre les équations du modèle mathématique : les éléments finis non conformes ont été utilisés pour résoudre l’équation de diffusivité et l’équation de Richards est résolue sous sa forme mixte par une méthode itérative en temps. En somme, ce modèle couplé permet de décrire les variations de niveaux de nappe à partir des données météorologiques connaissant les paramètres caractéristiques de cet aquifère. ...
Groundwater is the main available water resource for many countries; they are mainly replenished by water from precipitation, called groundwater recharge. Due to its great importance, management of groundwater resources is more essential than ever, and is achieved through mathematical models which offer us a better understanding of physical phenomena as well as their prediction. Hydrogeological systems are generally complex thus characterized by a highly variable dynamic over time and space. These complexities have attracted the attention of many hydrogeologists and many sophisticated models that can handle these issues and describe these systems accurately were developed. Unfortunately, modeling groundwater recharge is still a challenge in groundwater resource management. Generally, groundwater models are used to simulate aquifers flow without a good estimation of recharge and its spatial-temporal distribution. As groundwater recharge rates show spatial-temporal variability due to climatic conditions, land use, and hydrogeological heterogeneity, these methods have limitations in dealing with these characteristics. To overcome these limitations, a coupled model which simulates flow in the unsaturated zone and recharge as well as groundwater flow was developed. The flow in the unsaturated zone is solved either with resolution of Richards equation or with empirical models while the diffusivity equation governs flow in the saturated zone. Robust numerical methods were used to solve these equations: we apply non-conforming finite element to solve the diffusivity equation and we used an accurate and efficient method for solving the Richards equation. ...