Search

The sustainability of limited freshwater resources in coastal settings requires an understanding of the processes that affect them. This is especially relevant for freshwater lenses of oceanic islands. Yet, these processes are often obscured by dynamic oceanic water levels that change over a range of time scales. We use regression deconvolution to estimate an Oceanic Response Function (ORF) that accounts for how sea-level fluctuations affect measured groundwater levels, thus providing a clearer understanding of recharge and withdrawal processes. The method is demonstrated using sea-level and groundwater- level measurements on the island of Norderney in the North Sea (Northwest Germany). We expect that the method is suitable for any coastal groundwater system where it is important to understand processes that affect freshwater lenses or other coastal freshwater resources.

Water movement above the water table, yet below the root zone, is predominantly downward because of gravitational forces. This vertical movement is frequently affected by lower permeability geologic units that divert water horizontally. Also, vapor-driven transport of water can be important in very dry climates. Shallow soil-water is affected by plant uptake of water, and also by evaporation from the soil surface. These processes are characteristic of unsaturated media, i.e., zones where the soil materials are not saturated with water. In these cases, the surface tension of water causes adherance to soil surface, and soil-water pressures are generally negative. Each soil material has a distinctive relationship between the negative fluid pressure (or tension) and water content. The unsaturated hydraulic conductivity is lower than the hydraulic conductivity under saturated conditions because not all available pore space contributes to fluid flow. Deep soil-water movement is also affected by macropores (large pores) that can provide preferential flow paths under very wet conditions, but may not significantly contribute to flow otherwise. Keywords: deep percolation; unsaturated zone; vadose zone; unsaturated hydraulic conductivity; characteristic curves; capillary rise

Determining subsurface hydraulic and geomechanical properties crucially underpins groundwater resource investigation and management. While standard practice relies on active testing, passive approaches require less effort and cost but are underutilised. We present the new Python package named HydroGeoSines (HGS) which quantifies hydraulic and poroelastic subsurface properties using the groundwater response to natural forces (such as Earth tides and atmospheric pressure changes) embedded in standard measurements. All implemented methods are drawn from the peer-reviewed literature. The package includes basic handling of time series, such as joining and aligning records and handling gaps. HGS uses standard atmospheric and groundwater pressure records to estimate the Barometric Response Function (BRF) groundwater state of confinement, hydraulic conductivity, specific storage, barometric efficiency (BE) and porosity. If Earth tides are required, they can be calculated on-the-fly using the PyGTide package which is based on ETERNA and included. HGS allows easy compensation and correction of pressure or hydraulic heads from barometric pressure or Earth tide influences. Further, HGS includes import from and export to common data formats as well as visualisation of data and results. We demonstrate the use of HGS using example datasets from around the world. Since HGS unlocks sophisticated methods for use by anyone with Python skills, we anticipate that it will support subsurface investigations and add value to standard monitoring practice.

Brackish to saline groundwater in arid environments encourages the development and sustainability of inland freshwater lenses (IFLs). While these freshwater resources supply much-needed drinking water throughout the Arabian Peninsula and other drylands, little is understood about their sustainability. This study presents a numerical model using the SEAWAT programming code (i.e., MODFLOW and the Modular Three-Dimensional Multispecies Transport Model (MT3DMS)) to simulate IFL transient evolution. The numerical model is based on a physical laboratory model and calibrated using results from simulations conducted in a previous study of the Raudhatain IFL in northern Kuwait. Data from three previously conducted physical model simulations were evaluated against the corresponding numerical model simulations. The hydraulic conductivities in the horizontal and vertical directions were successfully optimized to minimize the objective function of the numerical model simulations. The numerical model matched observed IFL water levels at four locations through time, as well as IFL thicknesses and lengths (R2 = 0.89, 0.94, 0.85). Predicted lens degradation times corresponded to the observed lenses, which demonstrated the utility of numerical models and physical models to assess IFL geometry and position. Improved understanding of IFL dynamics provides water-resource exploration and development opportunities in drylands throughout the Arabian Peninsula and elsewhere with similar environmental settings.

The Dougherty Plain in southwest Georgia is a flat, karstic, depressional-landscape dominated by irrigated and dry-land agriculture devoted to row-crops and pasture with interspersed wetlands and forests. Stormwater runoff rarely discharges into perennial rivers and streams, except during large storms that induce hydrologic connectivity between fields, wetlands, and streams (event return period is less than one per year).We report the hydrologic and water-quality effects of a 173-mm rainfall event that generated three weeks (Feb 15 to Mar 9, 2014) of continuous flows through and between three normally isolated wetlands. A suite of water-quality parameters (physical, nutrients, and pathogen indicators) was monitored daily from offsite (agricultural) and onsite (forested) sources at two sites along one flowpath and five sites along a second at the Joseph W Jones Ecologic Research Center at Ichauway.Decreasing sediment, nutrient, and pathogen concentrations were observed as water moved across the forested landscapes with embedded wetlands. Two physical parameters (specific conductance and turbidity) were strongly-to-moderately correlated (r > 0.8, 0.5, respectively) with laboratory-measured parameters (e.g., nutrients, suspended solids, pathogens), which suggest their utility for routine stormwater monitoring and prioritizing sample collection for laboratory analyses at this site. Keywords: Longleaf-pine, Isolated wetlands, Stormflow, Agricultural runoff, Water quality, Dougherty plain, Nutrients, Pathogens

Accurate assessments of groundwater resources in major aquifers across the globe are crucial for sustainable management of freshwater reservoirs. Observations from the Gravity Recovery and Climate Experiment (GRACE) satellite have become invaluable as a means to identify regions groundwater change. While there is a large body of research that focuses on downscaling coarse (1°) GRACE products, few studies have attempted to spatially downscale GRACE to produce fine resolution (5 km) maps that are more useful to resource managers. This study trained a boosted regression tree model to statistically downscale GRACE total water storage anomaly to monthly 5 km groundwater level anomaly maps in the karstic upper Floridan aquifer (UFA) using multiple hydrologic datasets. Evaluation of spatial predictions with existing groundwater wells indicated satisfactory performance (R = 0.79, NSE = 0.61). Results demonstrate that groundwater levels were stable between 2002–2016 but varied seasonally. The data also highlights areas where groundwater pumping is exacerbating UFA water-level declines. While results demonstrate the applicability of machine learning based methods for spatial downscaling of GRACE data, future studies should account for preferential flowpaths (i.e., conduits, lineaments) in karstic systems.

Study region: The Lower Flint River Basin (LFRB): a karst catchment in southwestern Georgia, USA. Study focus: Using the U.S. EPA Reach File 3 data set, we generated stream reach azimuths for all tributaries of the Lower Flint River Basin (LFRB) in southwestern Georgia, USA, then compared these results to regional bedrock jointing orientations and stream chemistry (indicating incoming groundwater discharge) in one tributary of the LFRB, Ichawaynochaway Creek. Our objective was to determine if stream bearing might be a useful predictor of increased groundwater discharge in streams of the LFRB where groundwater development has significantly impacted baseflows. New hydrological insights for the region: We identified a dominant N-S trend in 44% of reaches in tributaries of the LFRB with lesser E-W, NNW, NW and NE trends. Bedrock joints and stream reaches in Ichawaynochaway Creek (a tributary of the Flint River) shared similar azimuth trends. When we compared stream reach orientation with known locations of enhanced groundwater inputs (previously detected by Rugel and others) we found that 55% of the time reaches in Ichawaynochaway Creek with increased groundwater discharge followed NW or NNW bearings (mean N49W). Further investigation to replicate these results in other tributaries of the LFRB is warranted and may help inform management strategies which could protect both ecological and economic interests in this region. Keywords: Grondwater discharge, Karst, Water resources, Groundwater-surface water interaction

Evapotranspiration (ET) typically accounts for 60–70% of precipitation in rural basins of the Southeastern United States. Since 1930, substantial reforestation of former croplands has occurred in the Piedmont and Appalachian Highlands in this area, leading to an expected increase in ET and reduction in baseflow. This study examines relationships between basin vegetative cover, abiotic factors, and water-budget partitioning in 45 USGS-gaged rural basins in the Southeastern US. Data are for the 1982–2014 water years with watersheds having ≥40% forest cover, crystalline-rock aquifers, minimal basin water export, and no large reservoirs. Long-term annual ET is calculated using the water-budget equation (ET = P-Q), which ranges from 641 to 971 mm/yr. (median 824). Vegetative cover and other basin variables are regressed against ET to quantify the effects of vegetative and forest types. Budyko analysis is employed to compare the watersheds and to evaluate factors affecting residuals. Regression analysis indicates that ET behavior is best explained by abiotic factors (i.e., precipitation and temperature) but forest-cover type also has some effect. Evergreen forest cover is less common than deciduous or mixed forest but has a positive relationship with ET, while deciduous and total forest have negative relationships with ET. Comparison of water-balance and Budyko-estimated ET indicates that deciduous and total forest are associated with negative residuals while evergreen is not significant. These results show that: forest cover effects on basin ET are complicated; forest-cover type is important for water-yield management in this region, and abiotic basin characteristics exert stronger control than forest cover on ET.

The authors of a new book on groundwater hydraulics and hydrology describe a theoretical foundation for modeling time-periodic groundwater flow.