Groundwater recharge assessment in a rural,arid, mid-mountain basin in North-Central Chile

ABSTRACT

A basic component of any hydrogeological study is the magnitude and temporal variation of groundwater recharge. This can be difficult to assess accurately, particularly in arid and semi-arid rainfed mid-mountain zones, as is the situation in the rural, low population density zones of North-Central Chile. In this study, recharge in the Punitaqui Basin, North-Central Chile, was characterized, contrasting the results of two methods: a modified Thornthwaite-Mather (MTM) and discharge recession analysis (DRA). We found a recharge rate of between 1 and 4% of average annual precipitation. Average recharge estimated by the MTM method is consistently higher than that estimated by DRA. Also, DRA tends to smooth the recharge values, resulting in a lower inter-annual variation coefficient. Both methods identified a threshold value of total annual precipitation, above which recharge can be expected to occur, of the order of 180 mm year^?1, consistent with values reported in similar areas.     

Seasonalizing mountain system recharge in semi-arid basins-climate change impacts

Climate variability and change impact groundwater resources by altering recharge rates. In semi-arid Basin and Range systems, this impact is likely to be most pronounced in mountain system recharge (MSR), a process which constitutes a significant component of recharge in these basins. Despite its importance, the physical processes that control MSR have not been fully investigated because of limited observations and the complexity of recharge processes in mountainous catchments. As a result, empirical equations, that provide a basin-wide estimate of mean annual recharge using mean annual precipitation, are often used to estimate MSR. Here North American Regional Reanalysis data are used to develop seasonal recharge estimates using ratios of seasonal (winter vs. summer) precipitation to seasonal actual or potential evapotranspiration. These seasonal recharge estimates compared favorably to seasonal MSR estimates using the fraction of winter vs. summer recharge determined from isotopic data in the Upper San Pedro River Basin, Arizona. Development of hydrologically based seasonal ratios enhanced seasonal recharge predictions and notably allows evaluation of MSR response to changes in seasonal precipitation and temperature because of climate variability and change using Global Climate Model (GCM) climate projections. Results show that prospective variability in MSR depends on GCM precipitation predictions and on higher temperature. Lower seasonal MSR rates projected for 2050-2099 are associated with decreases in summer precipitation and increases in winter temperature. Uncertainty in seasonal MSR predictions arises from the potential evapotranspiration estimation method, the GCM downscaling technique and the exclusion of snowmelt processes.     

Impacts of model initialization on an integrated surface water–groundwater model

Integrated hydrologic models characterize catchment responses by coupling the subsurface flow with land surface processes. One of the major areas of uncertainty in such models is the specification of the initial condition and its influence on subsequent simulations. A key challenge in model initialization is that it requires spatially distributed information on model states, groundwater levels and soil moisture, even when such data are not routinely available. Here, the impact of uncertainty in initial condition was explored across a 208 km² catchment in Denmark using the ParFlow.CLM model. The initialization impact was assessed under two meteorological conditions (wet vs dry) using five depth to water table and soil moisture distributions obtained from various equilibrium states (thermal, root zone, discharge, saturated and unsaturated zone equilibrium) during the model spin‐up. Each of these equilibrium states correspond to varying computation times to achieve stability in a particular aspect of the system state. Results identified particular sensitivity in modelled recharge and stream flow to the different initializations, but reduced sensitivity in modelled energy fluxes. Analysis also suggests that to simulate a year that is wetter than the spin‐up period, an initialization based on discharge equilibrium is adequate to capture the direction and magnitude of surface water–groundwater exchanges. For a drier or hydrologically similar year to the spin‐up period, an initialization based on groundwater equilibrium is required. Variability of monthly subsurface storage changes and discharge bias at the scale of a hydrological event show that the initialization impacts do not diminish as the simulations progress, highlighting the importance of robust and accurate initialization in capturing surface water–groundwater dynamics. Copyright © 2015 John Wiley & Sons, Ltd.     

Addressing snow model uncertainty for hydrologic prediction

Streamflow forecasting methods are moving towards probabilistic approaches that quantify the uncertainty associated with the various sources of error in the forecasting process. Multi-model averaging methods which try to address modeling deficiencies by considering multiple models are gaining much popularity. We have applied the Bayesian Model Averaging method to an ensemble of twelve snow models that vary in their heat and melt algorithms, parameterization, and/or albedo estimation method. Three of the models use the temperature-based heat and melt routines of the SNOW17 snow accumulation and ablation model. Nine models use heat and melt routines that are based on a simplified energy balance approach, and are varied by using three different albedo estimation schemes. Finally, different parameter sets were identified through automatic calibration with three objective functions. All models use the snow accumulation, liquid water transport, and ground surface heat exchange processes of the SNOW17. The resulting twelve snow models were combined using Bayesian Model Averaging (BMA). The individual models, BMA predictive mean, and BMA predictive variance were evaluated for six SNOTEL sites in the western U.S. The models performed best and the BMA variance was lowest at the colder sites with high winter precipitation and little mid-winter melting. An individual snow model would often outperform the BMA predictive mean. However, observed snow water equivalent (SWE) was captured within the 95% confidence intervals of the BMA variance on average 80% of the time at all sites. Results are promising that consideration of multiple snow structures would provide useful uncertainty information for probabilistic hydrologic prediction.     

RIPGIS-NET: a GIS tool for riparian groundwater evapotranspiration in MODFLOW

RIPGIS-NET, an Environmental System Research Institute (ESRI's) ArcGIS 9.2/9.3 custom application, was developed to derive parameters and visualize results of spatially explicit riparian groundwater evapotranspiration (ETg), evapotranspiration from saturated zone, in groundwater flow models for ecohydrology, riparian ecosystem management, and stream restoration. Specifically RIPGIS-NET works with riparian evapotranspiration (RIP-ET), a modeling package that works with the MODFLOW groundwater flow model. RIP-ET improves ETg simulations by using a set of eco-physiologically based ETg curves for plant functional subgroups (PFSGs), and separates ground evaporation and plant transpiration processes from the water table. The RIPGIS-NET program was developed in Visual Basic 2005, .NET framework 2.0, and runs in ArcMap 9.2 and 9.3 applications. RIPGIS-NET, a pre- and post-processor for RIP-ET, incorporates spatial variability of riparian vegetation and land surface elevation into ETg estimation in MODFLOW groundwater models. RIPGIS-NET derives RIP-ET input parameters including PFSG evapotranspiration curve parameters, fractional coverage areas of each PFSG in a MODFLOW cell, and average surface elevation per riparian vegetation polygon using a digital elevation model. RIPGIS-NET also provides visualization tools for modelers to create head maps, depth to water table (DTWT) maps, and plot DTWT for a PFSG in a polygon in the Geographic Information System based on MODFLOW simulation results.     

Impact of land-surface elevation and riparian evapotranspiration seasonality on groundwater budget in MODFLOW models

Riparian groundwater evapotranspiration (ETg) constitutes a major component of the water balance especially in many arid and semi-arid environments. Although spatial and temporal variability of riparian ETg are controlled by climate, vegetation and subsurface characteristics, depth to water table (DTWT) is often considered the major controlling factor. Relationships between ETg rates and DTWT, referred to as ETg curves, are implemented in MODFLOW ETg packages (EVT, ETS1 and RIP-ET) with different functional forms. Here, the sensitivity of the groundwater budget in MODFLOW groundwater models to ETg parameters (including ETg curves, land-surface elevation and ETg seasonality) are investigated. A MODFLOW model of the hypothetical Dry Alkaline Valley in the Southwestern USA is used to show how spatial representation of riparian vegetation and digital elevation model (DEM) processing methods impact the water budget when RIPGIS-NET (a GIS-based ETg program) is used with MODFLOW’s RIP-ET package, and results are compared with the EVT and ETS1 packages. Results show considerable impact on ETg and other groundwater budget components caused by spatial representation of riparian vegetation, vegetation type, fractional coverage areas and land-surface elevation. RIPGIS-NET enhances ETg estimation in MODFLOW by incorporating vegetation and land-surface parameters, providing a tool for ecohydrology studies, riparian ecosystem management and stream restoration.     

Complexity in microbial metabolic processes in soil nitrogen modeling: a case for model averaging

Model uncertainty is rarely considered in the field of biogeochemical modeling. The standard biogeochemical modeling approach is to proceed based on one selected model with the “right” complexity level based on data availability. However, other plausible models can result in dissimilar answer to the scientific question in hand using the same set of data. Relying on a single model can lead to underestimation of uncertainty associated with the results and therefore lead to unreliable conclusions. Multi-model ensemble strategy is a means to exploit the diversity of skillful predictions from different models with multiple levels of complexity. The aim of this paper is two fold, first to explore the impact of a model’s complexity level on the accuracy of the end results and second to introduce a probabilistic multi-model strategy in the context of a process-based biogeochemical model. We developed three different versions of a biogeochemical model, TOUGHREACT-N, with various complexity levels. Each one of these models was calibrated against the observed data from a tomato field in Western Sacramento County, California, and considered two different weighting sets on the objective function. This way we created a set of six ensemble members. The Bayesian Model Averaging (BMA) approach was then used to combine these ensemble members by the likelihood that an individual model is correct given the observations. Our results demonstrated that none of the models regardless of their complexity level under both weighting schemes were capable of representing all the different processes within our study field. Later we found that it is also valuable to explore BMA to assess the structural inadequacy inherent in each model. The performance of BMA expected prediction is generally superior to the individual models included in the ensemble especially when it comes to predicting gas emissions. The BMA assessed 95% uncertainty bounds bracket 90–100% of the observations. The results clearly indicate the need to consider a multi-model ensemble strategy over a single model selection in biogeochemical modeling study.