A Conceptual Framework for Integration Development of GSFLOW Model: Concerns and Issues Identified and Addressed for Model Development Efficiency

In Coupled Groundwater and Surface-Water Flow (GSFLOW) model, the three-dimensional finite-difference groundwater model (MODFLOW) plays a critical role of groundwater flow simulation, together with which the Precipitation-Runoff Modeling System (PRMS) simulates the surface hydrologic processes. While the model development of each individual PRMS and MODFLOW model requires tremendous time and efforts, further integration development of these two models exerts additional concerns and issues due to different simulation realm, data communication, and computation algorithms. To address these concerns and issues in GSFLOW, the present paper proposes a conceptual framework from perspectives of: Model Conceptualization, Data Linkages and Transference, Model Calibration, and Sensitivity Analysis. As a demonstration, a MODFLOW groundwater flow system was developed and coupled with the PRMS model in the Lehman Creek watershed, eastern Nevada, resulting in a smooth and efficient integration as the hydrogeologic features were well captured and represented. The proposed conceptual integration framework with techniques and concerns identified substantially improves GSFLOW model development efficiency and help better model result interpretations. This may also find applications in other integrated hydrologic modelings

Examination of the seepage face boundary condition in subsurface and coupled surface/subsurface hydrological models

A seepage face is a nonlinear dynamic boundary that strongly affects pressure head distributions, water table fluctuations, and flow patterns. Its handling in hydrological models, especially under complex conditions such as heterogeneity and coupled surface/subsurface flow, has not been extensively studied. In this paper, we compare the treatment of the seepage face as a static (Dirichlet) versus dynamic boundary condition, we assess its resolution under conditions of layered heterogeneity, we examine its interaction with a catchment outlet boundary, and we investigate the effects of surface/subsurface exchanges on seepage faces forming at the land surface. The analyses are carried out with an integrated catchment hydrological model. Numerical simulations are performed for a synthetic rectangular sloping aquifer and for an experimental hillslope from the Landscape Evolution Observatory. The results show that the static boundary condition is not always an adequate stand-in for a dynamic seepage face boundary condition, especially under conditions of high rainfall, steep slope, or heterogeneity; that hillslopes with layered heterogeneity give rise to multiple seepage faces that can be highly dynamic; that seepage face and outlet boundaries can coexist in an integrated hydrological model and both play an important role; and that seepage faces at the land surface are not always controlled by subsurface flow. The paper also presents a generalized algorithm for resolving seepage face outflow that handles heterogeneity in a simple way, is applicable to unstructured grids, and is shown experimentally to be equivalent to the treatment of atmospheric boundary conditions in subsurface flow models

2D approach for modelling self-potential anomalies. Application to synthetic and real data

The aim of this work is to present a 2-D Matlab code based on the finite element method for providing numerical modelling of both groundwater flow and self-potential signals. The distribution of the self-potential is obtained by starting with the solution of the groundwater flow, then computing the source current density, and finally calculating the electrical potential. The reliability of the algorithm is tested with synthetic case studies in order to simulate both the electric field resulting from the existence of a leak in the dam and SP signals associated with a pumping test in an unconfined aquifer. In addition, the algorithm was applied to field data for the localization of piping sinkholes. The results show that the outputs of the algorithm yielded satisfactory solutions, which are in good agreement with those of previous studies and field investigations. In details, the synthetic data and SP anomalies calculated by using the code are very close in terms of sign and magnitude, while real data tests clearly indicated that the computed SP signals were found to be consistent with the measured values

A state-dependent parameterization of saturated-unsaturated zone interaction

The relevance of groundwater as an important source of root zone moisture by means of capillary rise is increasingly being recognized. This is partly reflected in many current land surface schemes, which increasingly replace a one-way (i.e., downward) drainage of water by a two-way interaction flux between the root zone and a groundwater system. A fully physically correct implementation of this two-way saturated-unsaturated interaction flux requires transient simulations using the highly nonlinear Richards' equation, which is a computationally demanding approach. We test a classic simple approximation that computes the root zone¿groundwater interaction flux as the net effect of a downward drainage flux and an upward capillary rise flux against the Darcy equation for quasi steady state conditions. We find that for a wet root zone and/or shallow groundwater, the errors within this approximation are significant and of the same magnitude as the interaction flux itself. We present a new closed-form parameterization of the Darcy equation¿based fluxes that accounts both for root zone soil moisture and depth to the water table. Parameter values for this parameterization are listed for 11 different, widely applied soil texture descriptions. The high numerical efficiency of the proposed method makes it suitable for inclusion into demanding applications, e.g., a Monte Carlo framework, or high spatial resolution

Remote Detection of Saline Intrusion in a Coastal Aquifer Using Borehole Measurements of Self-Potential

Funded by NERC CASE studentship . Grant Number: NE/I018417/1Peer reviewedPublisher PD

Spatio-temporal stochastic resonance induces patterns in wetland vegetation dynamics

Water availability is a major environmental driver affecting riparian and wetland vegetation. The interaction between water table fluctuations and vegetation in a stochastic environment contributes to the complexity of the dynamics of these ecosystems. We investigate the possible emergence of spatial patterns induced by spatio-temporal stochastic resonance in a simple model of groundwater-dependent ecosystems. These spatio-temporal dynamics are driven by the combined effect of three components: (i) an additive white Gaussian noise, accounting for external random disturbances such as fires or fluctuations in rain water availability, (ii) a weak periodic modulation in time, describing hydrological drivers such as seasonal fluctuations of water table depth, and (iii) a spatial coupling term, which takes into account the ability of vegetation to spread and colonize other parts of the landscape. A suitable cooperation between these three terms is able to give rise to ordered structures which show spatial and temporal coherence, and are statistically steady in time.Comment: 9 pages, 7 figure

A new lattice Boltzmann model for interface reactions between immiscible fluids

In this paper, we describe a lattice Boltzmann model to simulate chemical reactions taking place at the interface between two immiscible fluids. The phase-field approach is used to identify the interface and its orientation, the concentration of reactant at the interface is then calculated iteratively to impose the correct reactive flux condition. The main advantages of the model is that interfaces are considered part of the bulk dynamics with the corrective reactive flux introduced as a source/sink term in the collision step, and, as a consequence, the model’s implementation and performance is independent of the interface geometry and orientation. Results obtained with the proposed model are compared to analytical solution for three different benchmark tests (stationary flat boundary, moving flat boundary and dissolving droplet). We find an excellent agreement between analytical and numerical solutions in all cases. Finally, we present a simulation coupling the Shan Chen multiphase model and the interface reactive model to simulate the dissolution of a collection of immiscible droplets with different sizes rising by buoyancy in a stagnant fluid