An official website of the United States government
Abstract Modern integrated hydrologic models (IHMs) are powerful tools for investigating coupled hydrologic system dynamics. The tradeoff for this realism is a high computational burden and large numbers of parameters in each cell, few of which can be specified with a high degree of confidence. These factors combined make uncertainty quantification (UQ) a problem for IHM‐based simulations, yet without rigorous UQ, it is not clear how much confidence can be placed on conclusions made with IHMs. Previous work evaluated steady‐state cases where the permeability field was a random variable, and a logical continuation is to consider transient conditions. This work assesses the confidence of an IHM representation of a first‐order basin in central Idaho, USA, using an ensemble of 250 permeability realizations under three different recharge forcing signals. The results show that surface water is simulated with high confidence across all the permeability realizations, but the groundwater system and changes to it have lower confidence. However, uncertainty in changes to the groundwater system decrease with time since an increase in the recharge, meaning that the farther one gets from a “peak” in the flow (of any size) the more confident one can be in the response (i.e., smaller inter‐quartile range). The ensemble was also used to assess how many realizations were needed to capture expected behaviors of the ensemble and their range of variability. Unsurprisingly, groundwater requires larger ensembles than surface flows, but the size of the ensembles necessary for convergence were smaller than initially expected.
more »
« less
Impacts of Permeability Uncertainty in a Coupled Surface‐Subsurface Flow Model Under Perturbed Recharge Scenarios
Abstract Coupled simulations of surface and variably saturated subsurface flow, termed integrated hydrologic models (IHMs), can provide powerful insights into the complex dynamics of watersheds. The system of governing equations solved by an IHM is non‐linear, making them a significant computational burden and challenging to accurately parameterize. Consequently, a large fraction of the IHM studies to date have been “numerical hypothesis testing” studies, but, as parallel computing continues to improve, IHMs are approaching the point where they might also be useful as predictive tools. For this to become reality, the predictive uncertainty of such highly parameterized simulations must be considered. However, uncertainty is seldom considered in the IHM literature, likely due to the long runtimes of the complex simulations. The questions considered herein are how much uncertainty is there in an IHM for a common watershed simulation scenario, and how likely is it that any one realization of a system will give the same relative change as any other due to a perturbation in recharge? A stochastic ensemble of 250 permeability field realizations was used to show that uncertainty in a high‐mountain headwaters systems is dominated by the subsurface. Recharge perturbation scenarios echo these results, but the uncertainty ofchangesin streamflow or groundwater pressure heads were significantly smaller than the uncertainty in their base‐case values. The main finding is that IHMs do provide confident, predictive estimates ofrelativechanges in watersheds, even when uncertainty in specific simulation outputs may be high.
more »
« less
- Award ID(s):
- 2049687
- PAR ID:
- 10493856
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Water Resources Research
- Volume:
- 60
- Issue:
- 3
- ISSN:
- 0043-1397
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
In forested, seasonally dry watersheds, winter rains commonly replenish moisture deficits in the vadose zone before recharging underlying hillslope groundwater systems that sustain streamflow. However, the relative inaccessibility of the subsurface has hindered efforts to include the role of storage deficits, primarily generated by plant-water uptake, in moderating groundwater recharge. Here, we compare groundwater recharge inferred from the storage-discharge relationship with independent, distributed estimates of vadose zone storage deficits across 12 undisturbed California watersheds, thereby tracking the evolution of the deficit-recharge relationship without intensive field instrumentation. We find accrued deficits during the dry season alone insufficiently explain differences in the wet season partitioning of rainfall due to the non-monotonic behavior of the deficit during the subsequent wet season. Tracking the deficit at the storm event-scale within the wet season, however, reveals a characteristic response in groundwater to increasing rainfall not captured in the seasonal analysis, and may improve estimates of the rainfall required to generate recharge and streamflow on a per-storm basis. Our findings demonstrate the potential for existing public datasets to better capture water partitioning within the subsurface using a combined deficit-recharge approach, though our analysis is currently limited to basins with select characteristics. CODE AVAILABLE ON GITHUB: https://github.com/noah-beniteznelson/recharge_deficitmore » « less
-
Abstract Observations show increases in river discharge to the Arctic Ocean especially in winter over the last decades but the physical mechanisms driving these changes are not yet fully understood. We hypothesize that even in the absence of a precipitation increase, permafrost degradation alone can lead to increased annual river runoff. To test this hypothesis we perform 12 millennium-long simulations over an idealized hypothetical watershed (1 km 2 ) using a distributed, physically based water balance model (Water flow and Balance Simulation Model, WaSiM). The model is forced by both a hypothetical warming defined by an air temperature increase of 7.5 ∘ C over 100 years, and a corresponding cooling scenario. To assess model sensitivity we vary soil saturated hydraulic conductivity and lateral subsurface flow configuration. Under the warming scenario, changes in subsurface water transport due to ground temperature changes result in a 7%–14% increase in annual runoff accompanied by a 6%–20% decrease in evapotranspiration. The increase in runoff is most pronounced in winter. Hence, the simulations demonstrate that changes in permafrost characteristics due to climate warming and associated changes in evapotranspiration provide a plausible mechanism for the observed runoff increases in Arctic watersheds. In addition, our experiments show that when lateral subsurface moisture transport is not included, as commonly done in global-scale Earth System Models, the equilibrium water balance in response to the warming or cooling is similar to the water balance in simulations where lateral subsurface transport is included. However, the transient changes in water balance components prior to reaching equilibrium differ greatly between the two. For example, for high saturated hydraulic conductivity only when lateral subsurface transport is considered, a period of decreased runoff occurs immediately after the warming. This period is characterized by a positive change in soil moisture storage caused by the soil moisture deficit developed during prior cooling.more » « less
-
Abstract. Mountainous regions act as the water towers of the worldby producing streamflow and groundwater recharge, a function that isparticularly important in semiarid regions. Quantifying rates of mountainsystem recharge is difficult, and hydrologic models offer a method toestimate recharge over large scales. These recharge estimates are prone touncertainty from various sources including model structure and parameters.The quality of meteorological forcing datasets, particularly in mountainousregions, is a large source of uncertainty that is often neglected ingroundwater investigations. In this contribution, we quantify the impact ofuncertainty in both precipitation and air temperature forcing datasets onthe simulated groundwater recharge in the mountainous watershed of theKaweah River in California, USA. We make use of the integrated surface water–groundwater model, ParFlow.CLM, and several gridded datasets commonly usedin hydrologic studies, downscaled NLDAS-2, PRISM, Daymet, Gridmet, andTopoWx. Simulations indicate that, across all forcing datasets, mountain front recharge is an important component of the water budget in themountainous watershed, accounting for 9 %–72 % of the annual precipitation and ∼90 % of the total mountain system recharge to theadjacent Central Valley aquifer. The uncertainty in gridded air temperatureor precipitation datasets, when assessed individually, results in similarranges of uncertainty in the simulated water budget. Variations in simulatedrecharge to changes in precipitation (elasticities) and air temperature(sensitivities) are larger than 1 % change in recharge per 1 % change inprecipitation or 1 ∘C change in temperature. The total volume ofsnowmelt is the primary factor creating the high water budget sensitivity, and snowmelt volume is influenced by both precipitation and air temperatureforcings. The combined effect of uncertainty in air temperature andprecipitation on recharge is additive and results in uncertainty levels roughly equal to the sum of the individual uncertainties depending on thehydroclimatic condition of the watershed. Mountain system recharge pathwaysincluding mountain block recharge, mountain aquifer recharge, and mountainfront recharge are less sensitive to changes in air temperature than changesin precipitation. Mountain front and mountain block recharge are moresensitive to changes in precipitation than other recharge pathways. Themagnitude of uncertainty in the simulated water budget reflects theimportance of developing high-quality meteorological forcing datasets in mountainous regions.more » « less
-
Modern hydrologic models have extraordinary capabilities for representing complex process in surface-subsurface systems. These capabilities have revolutionized the way we conceptualize flow systems, but how to represent uncertainty in simulated flow systems is not as well developed. Currently, characterizing model uncertainty can be computationally expensive, in part, because the techniques are appended to the numerical methods rather than seamlessly integrated. The next generation of computers, however, presents opportunities to reformulate the modeling problem so that the uncertainty components are handled more directly within the flow system simulation. Misconceptions about quantum computing abound and they will not be a “silver bullet” for solving all complex problems, but they might be leveraged for certain kinds of highly uncertain problems, such as groundwater (GW). The point of this issue paper is that the GW community could try to revise the foundations of our models so that the governing equations being solved are tailored specifically for quantum computers. The goal moving forward should not just be to accelerate the models we have, but also to address their deficiencies. Embedding uncertainty into the models by evolving distribution functions will make predictive GW modeling more complicated, but doing so places the problem into a complexity class that is highly efficient on quantum computing hardware. Next generation GW models could put uncertainty into the problem at the very beginning of a simulation and leave it there throughout, providing a completely new way of simulating subsurface flows.more » « less
