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.
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.
Groundwater hydrographs contain a rich set of information on the dynamics of aquifer systems and the processes and properties that influence them. While the importance of seasonal cycles in hydrologic and environmental state variables is widely recognized there has yet to be a comprehensive analysis of the seasonal dynamics of groundwater across the United States. Here we use time series of groundwater level measurements from 997 wells from the National Groundwater Monitoring Network to identify and describe groundwater seasonal cycles in unconfined aquifers across the United States. We use functional data analysis to obtain a functional form fit for each site and apply an unsupervised clustering algorithm to identify a set of five distinct seasonal cycles regimes. Each seasonal cycle regime has a distinctive shape and distinct timing of its annual minimum and maximum water level. There are clear spatial patterns in the occurrence of each seasonal cycle regime, with the relative occurrence of each regime strongly influenced by the geologic setting (aquifer system), climate, and topography. Our findings provide a comprehensive characterization of groundwater seasonal cycles across much of the United States and present both a methodology and results useful for assessing and understanding unconfined groundwater systems.