Ground water最新文献

Many sedimentary aquifers consist of small layers of coarser and finer material. When groundwater flow in these aquifers is modeled, the hydraulic conductivity may be simulated as homogeneous but anisotropic throughout the aquifer. In practice, the anisotropy factor, the ratio of the horizontal divided by the vertical hydraulic conductivity, is often set to 10. Here, numerical experiments are conducted to determine the effective anisotropy of an aquifer consisting of 400 horizontal layers of which the homogeneous and isotropic hydraulic conductivity varies over two orders of magnitude. Groundwater flow is simulated to a partially penetrating canal and a partially penetrating well. Numerical experiments are conducted for 1000 random realizations of the 400 layers, by varying the sequence of the layers, not their conductivity. It is demonstrated that the effective anisotropy of the homogeneous model is a model parameter that depends on the flow field. For example, the effective anisotropy for flow to a partially penetrating canal differs from the effective anisotropy for flow to a partially penetrating well in an aquifer consisting of the exact same 400 layers. The effective anisotropy also depends on the sequence of the layers. The effective anisotropy values of the 1000 realizations range from roughly 5 to 50 for the considered situations. A factor of 10 represents a median value (a reasonable value to start model calibration for the conductivity variations considered here). The median is similar to the equivalent anisotropy, defined as the arithmetic mean of the hydraulic conductivities divided by the harmonic mean.

Streamflow records are biased toward large streams and rivers, yet small headwater streams are often the focus of ecological research in response to climate change. Conventional flow measurement instruments such as acoustic Doppler velocimeters (ADVs) do not perform well during low-flow conditions in small streams, truncating the development of rating curves during critical baseflow conditions dominated by groundwater inflow. We revisited an instantaneous solute tracer injection method as an alternative to ADVs based on paired measurements to compare their precision, efficiency, and feasibility within headwater streams across a range of flow conditions. We show that the precision of discharge measurements using salt dilution by slug injection and ADV methods were comparable overall, but salt dilution was more precise during the lowest flows and required less time to implement. Often, headwater streams were at or below the depth threshold where ADV measurements could even be attempted and transects were complicated by coarse bed material and cobbles. We discuss the methodological benefits and limitations of salt dilution by slug injection and conclude that the method could facilitate a proliferation of streamflow observation across headwater stream networks that are highly undersampled compared to larger streams.

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 km² 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 m³/d and Biscayne Bay (5696 m³/d), and well fields (14,066 m³/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 CO₂ (sCO₂) removes water from brine held in pumice stone at levels well above the solubility of water in sCO₂. The higher water removal results from a combination of passive emulsification of water in sCO₂ and viscous fingering of sCO₂ 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 sCO₂ as well as its movement in the subsurface. The finding that the water concentration in sCO₂ 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.

Understanding fate and transport processes for per- and poly-fluoroalkyl substances (PFAS) is critical for managing impacted sites. "PFAS Salting Out" in groundwater, defined herein, is an understudied process where PFAS in fresh groundwater mixes with saline groundwater near marine shorelines, which increases sorption of PFAS to aquifer solids. While sorption reduces PFAS mass discharge to marine surface water, the fraction that sorbs to beach sediments may be mobilized under future salinity changes. The objective of this study was to conceptually explore the potential for PFAS Salting Out in sandy beach environments and to perform a preliminary broad-scale characterization of sandy shoreline areas in the continental U.S. While no site-specific PFAS data were collected, our conceptual approach involved developing a multivariate regression model that assessed how tidal amplitude and freshwater submarine groundwater discharge affect the mixing of fresh and saline groundwater in sandy coastal aquifers. We then applied this model to 143 U.S. shoreline areas with sandy beaches (21% of total beaches in the USA), indirectly mapping potential salinity increases in shallow freshwater PFAS plumes as low (<10 ppt), medium (10-20 ppt), or high (>20 ppt) along groundwater flow paths before reaching the ocean. Higher potential salinity increases were observed in West Coast bays and the North Atlantic coastline, due to the combination of moderate to large tides and large fresh groundwater discharge rates, while lower increases occurred along the Gulf of Mexico and the southern Florida Atlantic coast. The salinity increases were used to estimate potential perfluorooctane sulfonic acid (PFOS) sorption in groundwater due to salting out processes. Low-category shorelines may see a 1- to 2.5-fold increase in sorption of PFOS, medium-category a 2.0- to 6.4-fold increase, and high-category a 3.8- to 25-fold increase in PFOS sorption. The analysis presented provides a first critical step in developing a large-scale approach to classify the PFAS Salting Out potential along shorelines and the limitations of the approach adopted highlights important areas for further research.

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.