Implementation of a linear optimization water allocation algorithm into a fully integrated physical hydrology model

Connections between groundwater and surface water are well-established and remain the focus of significant ongoing research. Within the water management community, stream-aquifer interactions are widely acknowledged and often subject to regulation. Yet, the tools most commonly used by water managers often have limited physical hydrology and a resulting inability to represent feedbacks between elements of the physical system. This paper presents the development of a water allocation module (WAM) for an integrated hydrologic model, ParFlow. The module presented here is designed to maximize the satisfaction of demand subject to a system of priorities, preferences and constraints. Because ParFlow solves the groundwater and surface water systems simultaneously, the coupled management model can explicitly incorporate groundwater surface water interactions into management decisions. Proof of concept for the new module is demonstrated for a hypothetical water management scenario, utilizing moisture dependent irrigation, in a real world basin, the Little Washita.     

A novel platform to evaluate the dampening of water and solute transport in an experimental channel under unsteady flow conditions

This paper presents a novel platform to study the dampening of water and solute transport in an experimental channel under unsteady flow conditions, where literature data are scarce. We address the question about what could be the smallest size of experimental platform that is useful for research, project studies, and teaching activities and that allows to do rational experiments characterized by small space occupation, short experimental duration, high measurement precision, high quality and reproducible experimental curves, low water and energy consumption, and the possibility to test a large variety of hydrograph scenarios. Whereas large scale hydraulic laboratories have focused their studies on sediment transport, our platform deals with solute transport. The objectives of our study are (a) building a platform that allows to do rational experiments, (b) enriching the lack of experimental data concerning water and solute transport under unsteady state conditions, and (c) studying the dampening of water and solute transport. We studied solute transport in a channel with lateral gain and lateral loss under different experimental configurations, and we show how the same lateral loss flow event can lead to different lateral loss mass repartitions under different configurations. In order to characterize water and solute dampening between the input and the output of the channel, we calculate dampening ratios based on peak coordinates of time flow curves and time mass curves and that express the decrease of peak amplitude and the increase of peak occurrence time between the input and output curves. Finally, we use a solute transport model coupling the diffusive wave equation for water transfer and the advection–diffusion equation for solute transport in order to simulate the experimental data. The simulations are quite good with a Nash–Sutcliffe efficiency NSE > 0.98 for water transfer and 0.84 < NSE < 0.97 for solute transport. This platform could serve hydrological modellers because it offers a variety of measured parameters (flow, water height, and solute concentration), at a fine time step under unsteady flow conditions.     

Integrated surface–groundwater flow modeling: A free-surface overland flow boundary condition in a parallel groundwater flow model

Interactions between surface and groundwater are a key component of the hydrologic budget on the watershed scale. Models that honor these interactions are commonly based on the conductance concept that presumes a distinct interface at the land surface, separating the surface from the subsurface domain. These types of models link the subsurface and surface domains via an exchange flux that depends upon the magnitude and direction of the hydraulic gradient across the interface and a proportionality constant (a measure of the hydraulic connectivity). Because experimental evidence of such a distinct interface is often lacking in field systems, there is a need for a more general coupled modeling approach.     

A comparison between two statistical and a physically-based model in snow water equivalent mapping

Snow water equivalent (SWE) estimates at the end of the winter season have been compared for the 2002–2006 period in a 200 km² mountainous area in Switzerland, using three different models. The first model, ALPINE3D, is a physically based process-oriented model, which solves the snowpack energy and mass balance equations. The other two models, SWE-SEM and HS-SWE, are statistical algorithms interpolating snow data on a grid. While SWE-SEM interpolates local estimates of SWE, HS-SWE converts interpolated snow depth maps into maps of SWE using a regionally-calibrated conversion model. We discuss similarities and differences among the models’ results, both in terms of total volume, and spatial distribution of SWE. The comparison shows a general good agreement of the results of the three models, with a mean difference in the total volumes between the two statistical models of ∼8%, and between the physical model and the statistical ones of ∼−3% to −10%.     

Water flow and solute transport model of potentially toxic elements through unsaturated zone at regional wellfield Kosnica

The objective of this work was to build a prognostic water flow model and potentially toxic elements (lead, cadmium, zinc) transport model in the unsaturated zone. Research was conducted in the catchment area of Kosnica regional wellfield, where the unsaturated zone is characterised by Fluvisol. Lower sorption capacities were determined in the first horizons for all three potentially toxic elements. Correlation coefficient of the measured and simulated values of tracer concentration is 0.58 for the AC horizon and 0.84 for the 2C/C1 horizon. Based on calibrated water flow and transport parameters, a prognostic water flow model and potentially toxic elements (lead, cadmium, zinc) transport model in the unsaturated zone was built. In case of an accidental spill of potentially toxic elements with concentrations of 1000 mg/l, the risk of contamination of the aquifer is present. Copyright © 2016 John Wiley & Sons, Ltd.     

Coupling of a simultaneous heat and water model with a distributed hydrological model and evaluation of the combined model in a cold region watershed

Snow and frozen soil prevail in cold regions worldwide, and the integration of these processes is crucial in hydrological models. In this study, a combined model was developed by fully coupling a simultaneous heat and water model with a geomorphologically based distributed hydrological model. The combined model simulates vertical and lateral water transfer as well as vertical heat fluxes and is capable of representing the effects of frozen soil and snowmelt on hydrological processes in cold regions. This model was evaluated by using in situ observations in the Binggou watershed, an experimental watershed for cold region hydrology of the Watershed Allied Telemetry Experimental Research Project. Results showed that the model was able to predict soil freezing and thawing, unfrozen soil water content, and snow depth reasonably well. The simulated hydrograph was in good agreement with the in situ observation. The Nash–Sutcliffe coefficient of daily discharge was 0.744 for the entire simulation period, 0.472 from April to June, and 0.711 from June to November. This model can improve our understanding of hydrological processes in cold regions and assess the impacts of global warming on hydrological cycles and water resources. Copyright © 2012 John Wiley & Sons, Ltd.     

Data-based modelling approach for variable density flow and solute transport simulation in a coastal aquifer

Data-based models, namely artificial neural network (ANN), support vector machine (SVM), genetic programming (GP) and extreme learning machine (ELM), were developed to approximate three-dimensional, density-dependent flow and transport processes in a coastal aquifer. A simulation model, SEAWAT, was used to generate data required for the training and testing of the data-based models. Statistical analysis of the simulation results obtained by the four models show that the data-based models could simulate the complex salt water intrusion process successfully. The selected models were also compared based on their computational ability, and the results show that the ELM is the fastest technique, taking just 0.5 s to simulate the dataset; however, the SVM is the most accurate, with a Nash-Sutcliffe efficiency (NSE) ≥ 0.95 and correlation coefficient R ≥ 0.92 for all the wells. The root mean square error (RMSE) for the SVM is also significantly less, ranging from 12.28 to 77.61 mg/L.     

Numerical method of moments for solute transport in a nonstationary flow field conditioned on hydraulic conductivity and head measurements

The unconditional stochastic studies on groundwater flow and solute transport in a nonstationary conductivity field show that the standard deviations of the hydraulic head and solute flux are very large in comparison with their mean values (Zhang et al. in Water Resour Res 36:2107–2120, 2000; Wu et al. in J Hydrol 275:208–228, 2003; Hu et al. in Adv Water Resour 26:513–531, 2003). In this study, we develop a numerical method of moments conditioning on measurements of hydraulic conductivity and head to reduce the variances of the head and the solute flux. A Lagrangian perturbation method is applied to develop the framework for solute transport in a nonstationary flow field. Since analytically derived moments equations are too complicated to solve analytically, a numerical finite difference method is implemented to obtain the solutions. Instead of using an unconditional conductivity field as an input to calculate groundwater velocity, we combine a geostatistical method and a method of moment for flow to conditionally simulate the distributions of head and velocity based on the measurements of hydraulic conductivity and head at some points. The developed theory is applied in several case studies to investigate the influences of the measurements of hydraulic conductivity and/or the hydraulic head on the variances of the predictive head and the solute flux in nonstationary flow fields. The study results show that the conditional calculation will significantly reduce the head variance. Since the hydraulic head measurement points are treated as the interior boundary (Dirichlet boundary) conditions, conditioning on both the hydraulic conductivity and the head measurements is much better than conditioning only on conductivity measurements for reduction of head variance. However, for solute flux, variance reduction by the conditional study is not so significant.     

Simulation of field‐scale water flow and bromide transport in a cracked clay soil

Single domain models may seriously underestimate leaching of nutrients and pesticides to groundwater in clay soils with shrinkage cracks. Various two‐domain models have been developed, either empirical or physically based, which take into account the effects of cracks on water flow and solute transport. This paper presents a model concept that uses the clay shrinkage characteristics to derive crack volume and crack depth under transient field conditions. The concept has been developed to simulate field average behaviour of a field with cracks, rather than flow and transport at a small plot. Water flow and solute transport are described with basic physics, which allow process and scenario analysis. The model concept is part of the more general agrohydrological model SWAP, and is applied to a field experiment on a cracked clay soil, at which water flow and bromide transport were measured during 572 days. A single domain model was not able to mimic the field‐average water flow and solute transport. Incorporation of the crack concept considerably improved the simulation of water content and bromide leaching to the groundwater. Still deviations existed between the measured and simulated bromide concentration profiles. The model did not reproduce the observed bromide retardation in the top layer and the high bromide dispersion resulting from water infiltration at various soil depths. A sensitivity analysis showed that the amounts of bromide leached were especially sensitive to the saturated hydraulic conductivity of the top layer, the solute transfer from the soil matrix to crack water flow and the mean residence time of rapid drainage. The shrinkage characteristic and the soil hydraulic properties of the clay matrix showed a low sensitivity. Copyright © 2000 John Wiley & Sons, Ltd.     

An easily implemented agro-hydrological procedure with dynamic root simulation for water transfer in the crop–soil system: Validation and application

Models for water transfer in the crop–soil system are key components of agro-hydrological models for irrigation, fertilizer and pesticide practices. Many of the hydrological models for water transfer in the crop–soil system are either too approximate due to oversimplified algorithms or employ complex numerical schemes. In this paper we developed a simple and sufficiently accurate algorithm which can be easily adopted in agro-hydrological models for the simulation of water dynamics. We used a dual crop coefficient approach proposed by the FAO for estimating potential evaporation and transpiration, and a dynamic model for calculating relative root length distribution on a daily basis. In a small time step of 0.001 d, we implemented algorithms separately for actual evaporation, root water uptake and soil water content redistribution by decoupling these processes. The Richards equation describing soil water movement was solved using an integration strategy over the soil layers instead of complex numerical schemes. This drastically simplified the procedures of modeling soil water and led to much shorter computer codes. The validity of the proposed model was tested against data from field experiments on two contrasting soils cropped with wheat. Good agreement was achieved between measurement and simulation of soil water content in various depths collected at intervals during crop growth. This indicates that the model is satisfactory in simulating water transfer in the crop–soil system, and therefore can reliably be adopted in agro-hydrological models. Finally we demonstrated how the developed model could be used to study the effect of changes in the environment such as lowering the groundwater table caused by the construction of a motorway on crop transpiration.     

Evaluating the source and seasonality of submarine groundwater discharge using a radon-222 pore water transport model

Pore water radon (²²²Rn) distributions from Indian River Lagoon, Florida, are characterized by three zones: a lower zone where pore water ²²²Rn and sediment-bound radium (²²⁶Ra) are in equilibrium and concentration gradients are vertical; a middle zone where ²²²Rn is in excess of sediment-bound ²²⁶Ra and concentration gradients are concave-downward; and an upper zone where ²²²Rn concentration gradients are nearly vertical. These ²²²Rn data are simulated in a one-dimensional numerical model including advection, diffusion, and non-local exchange to estimate magnitudes of submarine groundwater discharge components (fresh or marine). The numerical model estimates three parameters, fresh groundwater seepage velocity, irrigation intensity, and irrigation attenuation, using two Monte Carlo (MC) simulations that (1) ensure the minimization algorithm converges on a global minimum of the merit function and the parameter estimates are consistent within this global minimum, and (2) provide 90% confidence intervals on the parameter estimates using the measured ²²²Rn activity variance. Model estimates of seepage velocities and discharge agree with previous estimates obtained from numerical groundwater flow models and seepage meter measurements and show the fresh water component decreases offshore and varies seasonally by a factor of nine or less. Comparison between the discharge estimates and precipitation patterns suggests a mean residence time in unsaturated and saturated zones on the order of 5 to 7 months. Irrigation rates generally decrease offshore for all sampling periods. The mean irrigation rate is approximately three times greater than the mean seepage velocity although the ranges of irrigation rates and seepage velocities are the same. Possible mechanisms for irrigation include density-driven convection, wave pumping, and bio-irrigation. Simulation of both advection and irrigation allows the separation of submarine groundwater discharge into fresh groundwater and (re)circulated lagoon water.     

Application of a Bayesian model to infer the contribution of coalbed natural gas produced water to the Powder River,Wyoming and Montana

The Powder River Basin (PRB) of Wyoming and Montana contains significant coal and coal bed natural gas (CBNG) resources. CBNG extraction requires the production of large volumes of water, much of which is discharged into existing drainages. Compared to surface waters, the CBNG produced water is high in sodium relative to calcium and magnesium, elevating the sodium adsorption ratio (SAR). To mitigate the possible impact this produced water may have on the quality of surface water used for irrigation, the State of Montana passed water anti‐degradation legislation, which could affect CBNG production in Wyoming. In this study, we sought to determine the proportion of CBNG produced water discharged to tributaries that reaches the Powder River by implementing a four end‐member mixing model within a Bayesian statistical framework. The model accounts for the ⁸⁷Sr/⁸⁶Sr, δ¹³C_DIC, [Sr] and [DIC] of CBNG produced water and surface water interacting with the three primary lithologies exposed in the PRB. The model estimates the relative contribution of the end members to the river water, while incorporating uncertainty associated with measurement and process error. Model results confirm that both of the tributaries associated with high CBNG activity are mostly composed of CBNG produced water (70–100%). The model indicates that up to 50% of the Powder River is composed of CBNG produced water downstream from the CBNG tributaries, decreasing with distance by dilution from non‐CBNG impacted tributaries from the point sources to ~10–20% at the Montana border. This amount of CBNG produced water does not significantly affect the SAR or electrical conductivity of the Powder River in Montana. Copyright © 2013 John Wiley & Sons, Ltd.