Miami-Dade County (MDC) has over 112,000 septic systems, some of which are at risk of compromise due to water table rise associated with sea level rise. MDC is surrounded by protected water bodies, including Biscayne Bay, with environmentally sensitive ecosystems and is underlain by highly transmissive karstic limestone. The main objective of the study is to provide first estimates of the locations and magnitudes of septic return flows to discharge endpoints. This is accomplished by leveraging MDC's county-scale surface-groundwater model using pathline analysis to estimate the transport and discharge fate of septic system flows under the complex time history of groundwater flow response to pumping, canal management, storms, and other environmental factors. The model covers an area of 4772 km2 in Southeast Florida. Outputs from the model were used to create a 30-year (2010 to 2040) simulation of the spatial-temporal pathlines from septic input locations to their termination points, allowing us to map flow paths and the spatial distribution of the septic flow discharge endpoints under the simulated conditions. Most septic return flows were discharged to surface water, primarily canals 52,830 m3/d and Biscayne Bay (5696 m3/d), and well fields (14,066 m3/d). Results allow us to identify "hotspots" to guide water quality sampling efforts and to provide recommendations for septic-to-sewer conversion areas that should provide most benefit by reducing nutrient loading to water bodies.

Numerical modeling of the recovery of moisture by injecting warm air in the unsaturated zone in a 100 m × 100 m plot of agricultural land in Kuwait, a country located in an arid environment, was conducted to provide "proof of concept" of the technique. If technically and economically feasible, it will be a potential additional source of water that could be exploited for farming activities and other uses. The COMSOL software was used to develop the model and, based on the results of the scenario runs, the effects of different hydraulic and operational parameters, including that of well spacing, on moisture recovery were assessed. In general, the results suggested that the recovery should increase with the increase in the hydraulic conductivity of the unsaturated zone, the amount of heat input, and the pressure differential between the unsaturated zone and the well head. Within the period examined (0 to 11 days), the recovery decreases with the increase in the soil moisture content, possibly due to the fall in relative permeability to moisture-rich air with the increased water contents in the pore spaces, although the effects may change over a longer period as water contents decrease with moisture recovery. The moisture recovery from the unsaturated zone through the injection of warm air appears to be a feasible proposition from this study that should be demonstrated through a pilot scale experiment in the field.

Supercritical CO2 (sCO2) removes water from brine held in pumice stone at levels well above the solubility of water in sCO2. The higher water removal results from a combination of passive emulsification of water in sCO2 and viscous fingering of sCO2 through the saturated pumice. This leads to higher levels of salt deposition than that expected from solubility considerations alone. These deposits could impact the injectivity of sCO2 as well as its movement in the subsurface. The finding that the water concentration in sCO2 is not necessarily capped at the solubility limit should influence the parametrization of injection models.

In this paper, we review the derivation of the Gauss-Levenberg-Marquardt (GLM) algorithm and its extension to ensemble parameter estimation. We explore the use of graphical methods to provide insights into how the algorithm works in practice and discuss the implications of both algorithm tuning parameters and objective function construction in performance. Some insights include understanding the control of both parameter trajectory and step size for GLM as a function of tuning parameters. Furthermore, for the iterative Ensemble Smoother (iES), we discuss the importance of noise on observations and show how iES can cope with non-unique outcomes based on objective function construction. These insights are valuable for modelers using PEST, PEST++, or similar parameter estimation tools.

Water-table maps are fundamental to hydrogeological studies and a manual, hand-drawn method is still commonly used to produce them. Despite this, the accuracy and variability of such maps have received little attention in international literature. In a unique experiment, 63 groundwater professionals drew water-table equipotential contours based on the same dataset of point measurements and were asked to infer flow directions and predict groundwater elevations at predefined locations. The root mean squared error (RMSE) for the average map calibration data was 10.5 m, which is accuracy comparable to numerical groundwater models. This study confirmed that to produce hand-drawn water-table maps, practitioners seek to not only fit the spatial data, but also to conform to their own cognitive model of hydrogeological concepts and processes. The calibration accuracy increased with experience; from a RMSE of 13.3 m for practitioners with 0-3 years of experience to a RMSE of 9.2 m for those with four or more years. Despite considerable variability in the style of the hand-drawn water-table maps, the maps were consistent in their representation of the dominant regional groundwater flow directions. There was less consensus, however, in predicting the direction of surface water-groundwater interaction for a stream reach. Hand-drawn water-table mapping remains useful and valid, especially as a starting point for hydrogeological conceptualization, yet further work is required to resolve issues around transparency, repeatability, and reproducibility.

Analytical and semi-analytical models for stream depletion with transient stream stage drawdown induced by groundwater pumping are developed to address a deficiency in existing models, namely, the use of a fixed stream stage condition at the stream-aquifer interface. Field data are presented to demonstrate that stream stage drawdown does indeed occur in response to groundwater pumping near aquifer-connected streams. A model that predicts stream depletion with transient stream drawdown is developed based on stream channel mass conservation and finite stream channel storage. The resulting models are shown to reduce to existing fixed-stage models in the limit as stream channel storage becomes infinitely large, and to the confined aquifer flow with a no-flow boundary at the streambed in the limit as stream storage becomes vanishingly small. The model is applied to field measurements of aquifer and stream drawdown, giving estimates of aquifer hydraulic parameters, streambed conductance, and a measure of stream channel storage. The results of the modeling and data analysis presented herein have implications for sustainable groundwater management.

In groundwater modeling studies, accurate spatial and intensity identification of water sources and sinks is of critical importance. Precise construction data about wells (water sinks) are particularly difficult to obtain. The collection of well log data is expensive and laborious, and government records of historic well log data are often imprecise and incomplete with respect to the precise location or pumping rate. In many groundwater modeling studies, such as groundwater quality assessments, a precise representation of the horizontal and vertical distribution of well screens is required to accurately estimate contaminant breakthrough curves. The number of wells under consideration may be very large, for example, in the assessment of nonpoint source pollution. In this paper, we propose an imputation framework that allows for proper reconstruction of missing well data. Our approach exploits available information and tolerates data gaps and imprecisions. We demonstrate the value of this method for a subregion of the Central Valley aquifer (California, USA). We show that our framework imputes missing values that preserve statistical properties of available data and that remain consistent with the known spatial distribution of well screens and pumping rates in the three-dimensional aquifer system.

This study advances a methodology to estimate effective apertures of fractures in glacial tills based on dye tracer infiltration tests and numerical simulations. The approach uses the visible penetration depth of the dye tracer along fracture flow paths as primary information to calculate effective fracture apertures. Further data used in the calculation are the dye tracer input concentration and retardation, the duration of the tracer injection, and the hydraulic gradient applied to control the infiltrating water fluxes. The method does not require measurement of hydraulic conductivity for the fractured till and enables direct observation of flow and transport patterns within the fractures (e.g., uniform flow and dye tracer distribution, channeling due to aperture variability, and presence of biogenic macropores in fractures). The approach was successfully verified by using the estimated effective fracture aperture values in Large Undisturbed Columns (LUCs) to consistently simulate both the observed LUC effluent breakthrough of a conservative bromide tracer and the water fluxes with the hydraulic gradient applied in the experiments. Sensitivity analyses revealed that estimation of small effective fracture apertures (<10 μm) required accurate determination of the dye tracer retardation factor. By contrast, in the case of larger effective apertures (>20 μm), the sensitivity of the estimated effective fracture aperture to variations in the porous material and solute transport parameters was low compared to the dominant sensitivity to the water flow through the fractures (cubic relation between flow and aperture). The proposed approach may be extended beyond laboratory applications and assist in characterizing field-scale fracture networks.

Faults can fundamentally change a groundwater flow regime and represent a major source of uncertainty in groundwater studies. Much research has been devoted to uncertainty around their location and their barrier-conduit behavior. However, fault timing is one aspect of fault uncertainty that appears to be somewhat overlooked. Many faulted models feature consistent layer offsets, thereby presuming that block faulting has occurred recently and almost instantaneously. Additionally, barrier and/or conduit behavior is often shown to extend vertically through all layers when a fault may in fact terminate well below-ground surface. In this study, we create three plausible geological interpretations for a transect in the Perth Basin. Adjacent boreholes show stratigraphic offsets and thickening which indicate faulting; however, fault timing is unknown. Flow modeling demonstrates that the model with the most recent faulting shows profoundly different flow patterns due to aquifer juxtaposition. Additionally, multiple realizations with stochastically generated parameter sets for layer, fault core, and fault damage zone conductivity show that fault timing influences flow more than layer or fault zone conductivity. Finally, fault conduit behavior that penetrates aquitards has significant implications for transport, while fault barrier behavior has surprisingly little. This research advocates for adequate data collection where faults may cause breaches in aquitards due to layer offsets or conduit behavior in the damage zone. It also promotes the use of multiple geological models to address structural uncertainty, and highlights some of the hurdles in doing so such as computational expense and the availability of seamless geological-flow modeling workflows.

Open pit mining frequently requires regional water tables to be lowered to access ore deposits. When mines close, dewatering ceases allowing the water table to recover. In arid and semi-arid mining regions, the developing pit lakes are predominantly fed by groundwater during this recovery phase and pit lakes develop first into "terminal sinks" for the surrounding groundwater system. With time, the re-establishment of regional hydraulic gradients can cause pit lakes to develop into throughflow systems, in which pit lake water outflows into adjacent aquifers. In this study, we use numerical groundwater modeling to aid process understanding of how regional hydraulic gradients, aquifer properties, net evaporation rates, and pit geometry determine the hydraulic evolution of groundwater-fed pit lakes. We find that before the recovery of the regional water table to its new equilibrium, pit lakes frequently transition to throughflow systems. Throughflow from pit lakes to downstream aquifers can develop within two decades following cessation of dewatering even under low hydraulic gradients (e.g., 5 × 10-4) or high net evaporation rates (e.g., 2.5 m/year). Pit lakes remain terminal sinks only under suitable combinations of high evaporation rates, low hydraulic gradients, and low hydraulic conductivities. In addition, we develop an approximate analytical solution for a rapid assessment of the hydraulic status of pit lakes under steady-state conditions. Understanding whether pit lakes remain terminal sinks or transition into throughflow systems largely determines the long-term water quality of pit lakes and downstream aquifers. This knowledge is fundamental for mine closure and planning post-mining land use.