Occurrence and concentrations of traditional and emerging contaminants in onsite wastewater systems and water supply wells in eastern North Carolina, USA.

Onsite wastewater treatment systems (OWTSs) and private wells are commonly used in Eastern North Carolina, USA. Water from private wells is not required to be tested after the initial startup, and thus persons using these wells may experience negative health outcomes if their water is contaminated with waste-related pollutants including bacteria, nitrate or synthetic chemicals such as hexafluoropropylne oxide dimer acid and its ammonium salt (GenX). Water samples from 18 sites with OWTSs and groundwater wells were collected for nitrate, Escherichia coli (E. coli), total coliform, and GenX concentration analyses. Results showed that none of the 18 water supplies were positive for E. coli, nitrate concentrations were all below the maximum contaminant level of 10 mg L-1, and one well had 1 MPN 100 mL-1 of total coliform. However, GenX was detected in wastewater collected from all 18 septic tanks and 22% of the water supplies tested had concentrations that exceeded the health advisory levels for GenX. Water supplies with low concentrations of traditionally tested for pollutants (nitrate, E. coli) may still pose health risks due to elevated concentrations of emerging contaminants like GenX and thus more comprehensive and routine water testing is suggested for this and similar persistent compounds.

Risk factors for Cryptosporidium contamination in Minnesota public supply wells.

In a recent monitoring study of Minnesota's public supply wells, Cryptosporidium was commonly detected with 40% of the wells having at least one detection. Risk factors for Cryptosporidium occurrence in drinking water supply wells, beyond surface water influence, remain poorly understood. To address this gap, physical and chemical factors were assessed as potential predictors of Cryptosporidium occurrence in 135 public supply wells in Minnesota. Univariable analysis, regression techniques, and classification trees were used to analyze the data. Many variables were identified as significant risk factors in univariable analysis and several remained significant throughout the succeeding analysis techniques. These factors fell into general categories of well use and construction, aquifer characteristics, and connectedness to the land surface, well capture zones, and land use therein, existence of potential contaminant sources within 200-feet of the well, and variability in the chemical and isotopic parameters measured during the study. These risk categories, and the specific variables and threshold values we have identified, can help guide future research on factors influencing Cryptosporidium contamination of wells and can be used by environmental health programs to develop risk-based sampling plans and design interventions that reduce associated health risks.

Integrated management of groundwater quantity, physicochemical properties, and microbial quality in West Nile delta using a new MATLAB code and geographic information system mapping.

Groundwater is an excellent alternative to freshwater for drinking, irrigation, and developing arid regions. Agricultural, commercial, industrial, residential, and municipal activities may affect groundwater quantity and quality. Therefore, we aimed to use advanced methods/techniques to monitor the piezometric levels and collect groundwater samples to test their physicochemical and biological characteristics. Our results using software programs showed two main types of groundwater: the most prevalent was the Na-Cl type, which accounts for 94% of the groundwater samples, whereas the Mg-Cl type was found in 6% of samples only. In general, the hydraulic gradient values, ranging from medium to low, could be attributed to the slow movement of groundwater. Salinity distribution in groundwater maps varied between 238 and 1350 mg L-1. Although lower salinity values were observed in northwestern wells, higher values were recorded in southern ones. The collected seventeen water samples exhibited brackish characteristics and were subjected to microbial growth monitoring. Sample WD12 had the lowest total bacterial count (TBC) of 4.8 ± 0.9 colony forming unit (CFU mg L-1), while WD14 had the highest TBC (7.5 ± 0.5 CFU mg L-1). None of the tested water samples, however, contained pathogenic microorganisms. In conclusion, the current simulation models for groundwater drawdown of the Quaternary aquifer system predict a considerable drawdown of water levels over the next 10, 20, and 30 years with the continuous development of the region.

Do animal husbandry operations contaminate groundwater sources with antimicrobial resistance: systematic review.

Antimicrobial resistance (AMR) is a critical global health concern. Animal husbandry operations are AMR hotspots due to heavy antibiotic use and dissemination of animal waste into the environment. In this systematic review, we examined the impact of swine, poultry, and cattle operations on AMR in groundwater. We searched PubMed, Web of Science, CAB Direct, and the North Carolina State University Agricultural and Environmental Science databases in June 2022. The search returned 2487 studies. Of the 23 eligible studies, 17 were conducted in high-income countries (primarily the USA, also Canada, Saudi Arabia, Cyprus), and 6 were conducted in a single upper-middle-income country (China). Studies investigated facilities for swine (13), poultry (4), cattle (3), and multiple types of animals (3). The sampling distance ranged from onsite to > 20 km from facilities; the majority of studies (19) sampled onsite. Most studies collected samples from monitoring wells; only 5 studies investigated private drinking water wells. AMR in groundwater was associated with animal husbandry operations in 74% (17/23) of all studies, 65% (11/17) of studies in high-income countries, and 100% (6/6) of studies in China. Contamination was mostly found in onsite wells, especially downgradient of waste lagoons, but also in offsite private wells up to 2-3 km away. Few studies reported weather data, but AMR contamination appeared to increase with rainy conditions. Future studies should sample private wells at varying distances from animal husbandry operations under different weather conditions and include low- and middle-income countries where food animal production is intensifying.

Private groundwater contamination and risk management: A comparative scoping review of similarities, drivers and challenges across two socio-economically developed regions.

Consolidation of multi-domain risk management research is essential for strategies facilitating the concerted government (educational) and population-level (behavioural) actions required to reduce microbial private groundwater contamination. However, few studies to date have synthesised this literature or sought to ascertain the causal generality and extent of supply contamination and preventive responses. In light of the Republic of Ireland (ROI) and Ontario's high reliance and research focus on private wells and consequent utility for empirical comparison, a scoping review of pertinent literature (1990-2022) from both regions was undertaken. The SPICE (Setting, Perspective, Intervention, Comparison, Evaluation) method was employed to inform literature searches, with Scopus and Web of Science selected as primary databases for article identification. The review identified 65 relevant articles (Ontario = 34, ROI = 31), with those investigating well user actions (n = 22) and groundwater quality (n = 28) the most frequent. A markedly higher pooled proportion of private supplies in the ROI exhibited microbial contamination (38.3 % vs. 4.1 %), despite interregional similarities in contamination drivers (e.g., weather, physical supply characteristics). While Ontarian well users demonstrated higher rates of historical (≥ 1) and annual well testing (90.6 % vs. 71.1 %; 39.1 % vs. 8.6 %) and higher rates of historical well treatment (42.3 % vs. 24.3 %), interregional levels of general supply knowledge were analogous (70.7 % vs. 71.0 %). Financial cost, organoleptic properties and residence on property during supply construction emerged as predictors of cognition and behaviour in both regions. Review findings suggest broad interregional similarities in drivers of supply contamination and individual-level risk mitigation, indicating that divergence in contamination rates may be attributable to policy discrepancies - particularly well testing incentivisation. The paucity of identified intervention-oriented studies further highlights the importance of renewed research and policy agendas for improved, targeted well user outreach and incentivised, convenience-based services promoting routine supply maintenance.

Novel approach to roof rainwater harvesting and aquifer recharge in an urban environment: Dry and wet infiltration wells comparison.

In urban environments there is a severe reduction of infiltration and groundwater recharge due to the existence of large impervious areas. During rain events, large volumes of water that could have recharged groundwater and surface water bodies are diverted into the municipal drainage system and lost from the freshwater storage. Moreover, extreme rain events impose high peak flows and large runoff volumes, which increase the risk of urban floods. Recent studies have suggested the use of rainwater harvesting for groundwater recharge, as a plausible solution for these challenges in dense urban environments. While the benefits of this approach are well understood, research on its practical, engineering, and hydrological aspects is relatively limited. The objective of the present study was to examine the use of infiltration wells for groundwater recharge with harvested rainwater collected from building rooftops under Mediterranean climate conditions. Two types of wells with similar hydraulic and technical properties were examined: a well that reaches the groundwater (wet well); and a well that discharges the harvested water into the unsaturated zone (dry well). Infiltration capacities of the wells were compared in controlled experiments conducted during summer months, and in operational recharge of harvested rainwater, during winter. Both dry and wet wells were found to be suitable for purposes of groundwater recharge with rooftop-harvested rainwater. Infiltration capacity of the wet well was about seven times greater than the infiltration capacity of the dry well. While the infiltration capacity of the wet well was constant throughout the entire length of the study (∼10 m3/h/m), the dry well infiltration capacity improved during winter (from 0.5 m3/h/m to 1.5 m3/h/m), a result of development of the dry well with time. Considering Tel-Aviv, Israel, as a case study for a dense modern city in a Mediterranean climate, it is demonstrated herein that the use of infiltration wells may reduce urban drainage by ∼40 %.

A brief note on substantial sub-daily arsenic variability in pumping drinking-water wells in New Hampshire.

Large variations in redox-related water parameters, like pH and dissolved oxygen (DO), have been documented in New Hampshire (United States) drinking-water wells over the course of a few hours under pumping conditions. These findings suggest that comparable sub-daily variability in dissolved concentrations of redox-reactive and toxic arsenic (As) also may occur, representing a potentially critical public-health data gap and a fundamental challenge for long-term As-trends monitoring. To test this hypothesis, discrete groundwater As samples were collected approximately hourly during one day in May and again in August 2019 from three New Hampshire drinking-water wells (2 public-supply, 1 private) under active pumping conditions. Collected samples were assessed by laboratory analysis (total As [AsTot], As(III), As(V)) and by field analysis (AsTot) using a novel integrated biosensor system. Laboratory analysis revealed sub-daily variability (range) in AsTot concentrations equivalent to 16 % - 36 % of that observed in the antecedent 3-year bimonthly trend monitoring. Thus, the results indicated that, along with previously demonstrated seasonality effects, the timing and duration of pumping are important considerations when assessing trends in drinking-water As exposures and concomitant risks. Results also illustrated the utility of the field sensor for monitoring and management of AsTot exposures in near-real-time.

Estimating Lithium Concentrations in Groundwater Used as Drinking Water for the Conterminous United States.

Lithium (Li) concentrations in drinking-water supplies are not regulated in the United States; however, Li is included in the 2022 U.S. Environmental Protection Agency list of unregulated contaminants for monitoring by public water systems. Li is used pharmaceutically to treat bipolar disorder, and studies have linked its occurrence in drinking water to human-health outcomes. An extreme gradient boosting model was developed to estimate geogenic Li in drinking-water supply wells throughout the conterminous United States. The model was trained using Li measurements from ∼13,500 wells and predictor variables related to its natural occurrence in groundwater. The model predicts the probability of Li in four concentration classifications, ≤4 µg/L, >4 to ≤10 µg/L, >10 to ≤30 µg/L, and >30 µg/L. Model predictions were evaluated using wells held out from model training and with new data and have an accuracy of 47-65%. Important predictor variables include average annual precipitation, well depth, and soil geochemistry. Model predictions were mapped at a spatial resolution of 1 km2 and represent well depths associated with public- and private-supply wells. This model was developed by hydrologists and public-health researchers to estimate Li exposure from drinking water and compare to national-scale human-health data for a better understanding of dose-response to low (<30 µg/L) concentrations of Li.

Simulation-Optimization Approach for Siting Injection Wells in Urban Area with Complex Hydrogeology.

Managed aquifer recharge has become a standard water resources management practice to promote the development of locally sustainable water supplies and combat water scarcity. However, installation of injection wells for replenishment purposes in urban areas with complex hydrogeology faces many challenges, such as limited land availability, potential impacts on municipal production wells and known subsurface contamination plumes, and complex spatially variable hydraulic connections between aquifer units. To assess the feasibility and cost-effectiveness of injecting advanced treated water (ATW) into a complex urban aquifer system, a Simulation-Optimization (SO) model was developed to automate a systematic search for the most cost-effective locations to install new wells for injecting various quantities of ATW, if feasible. The generalized workflow presented here uses an existing MODFLOW groundwater model-along with advanced optimization routines that are publicly available-to flexibly accommodate a multiobjective function, complex constraints, and specific project requirements. The model successfully placed wells for injection of 1 to 4 MGD of ATW in aquifers underlying the study area. The injection well placement was primarily constrained by avoiding excessive impact on environmental sites with underlying groundwater plumes. The largest costs were for well installation and piping to the wells from the existing ATW pipes. This workflow is readily adaptable to other sites with different complexities, decision variables, or constraints.

Automated Estimation of Aquifer Parameters from Arbitrary-Rate Pumping Tests in Python and MATLAB.

Inspired by the analysis by Mishra et al. (2012) of variable pumping rate tests using piecewise-linear reconstructions of the pumping history, this article contains a derivation of the convolutional form of pumping tests in which the pumping history may take any possible form. The solution is very similar to the classical Theis (1935) equation but uses the Green's function for a pumped aquifer given by taking the time derivative of the well function W ( u ( t ) ) . This eliminates one integration inside another and renders the convolution including the pumping history about as computationally demanding as calculating the well function alone, so that the convolution can be completed using handy mathematical software. It also allows nonlinear well losses, and because an easily-computed deterministic model exists for all data points and pumping history, an objective function may include all data, so that errors are reduced in calculating any nonlinear-well losses. In addition, data from multiple observation wells may be used simultaneously in the inversion. We provide codes in MATLAB and Python to solve for drawdown resulting from an arbitrary pumping history and compute the optimal aquifer parameters to fit the data. We find that the subtleties in parameter dependencies and constructing an appropriate objective function have a substantial effect on the interpreted parameters. Furthermore, the optimization from step-drawdown tests is typically nonunique and strongly suggests that a Bayesian inversion should be used to fully estimate the joint probability density of the parameter vector.

Assessment of wellbore integrity failure risk and hazardous zones in depleted reservoirs underground gas storage during the operation processes.

Depleted reservoirs are widely used for underground gas storage due to their advantages of low construction cost and easy development. Under the influence of complex geological conditions and frequent operations, the integrity of the wells in depleted reservoirs is prone to failure, which would potentially lead to gas leakage. In this study, by using a finite element-based computational fluid dynamics model, we have developed evaluation criteria for assessing the severity of the occurred wellbore integrity failure and the risk of the un-occurred wellbore integrity failures respectively to identify hazardous zones potentially prone to wellbore integrity failure. The study results indicate that the gas storage wellbore integrity failure is prone to occur inside the wellbore structure in the direction of the minimum ground stress near the lower boundary of the formation interlayer. The wellbore integrity failure hazardous zones are mainly concentrated at the formation interlayer boundaries. The practical guidelines and solutions derived from current research results can provide an accurate direction for monitoring and protecting work of wellbore integrity and avoid environment pollution problems caused by natural gas leakage.

Delineation of a PFOA Plume and Assessment of Data Gaps in its Conceptual Model Using PlumeSeeker™.

An accurate conceptual site model (CSM) and plume-delineation at contamination sites are pre-requisites for successful remediation and for satisfying regulators and stakeholders. PlumeSeeker™ is well-suited for assessing data gaps in CSMs by using available site data and for identifying the optimal number and locations of sampling locations to delineate contaminant plumes. It is an enhancement of a university research code for plume delineation using geostatistical and stochastic modeling integrated with the groundwater modeling software MODFLOW-SURFACT™. PlumeSeeker™ increases the overall confidence in the location of the plume boundary through a variance-reduction approach that selects existing- or new monitoring wells for sampling based on minimizing the uncertainty in plume boundary and on new field information. Applicable at sites with or without existing monitoring wells, PlumeSeeker™ is particularly powerful for optimally allocating project resources (labor, well installation, and laboratory costs) between existing wells and sampling at new locations. An application of PlumeSeeker™ at Lakehurst, the naval component of Joint Base McGuire-Dix-Lakehurst in New Jersey, demonstrates how the cost of delineating the migration pathway of a perfluorooctanoic acid (PFOA) plume can be minimized by requiring only 9 new sampling locations in addition to samples from 2 existing wells for achieving a 70% reduction in plume uncertainty. In addition, the use of available site data in three different scenarios identified CSM data-gaps in the source area and in the interaction between Manapaqua Branch and groundwater, where the observed high concentration in this area could have resulted from a combination of groundwater migration and induced infiltration.

Continental Scale Hydrostratigraphy: Comparing Geologically Informed Data Products to Analytical Solutions.

This study synthesizes two different methods for estimating hydraulic conductivity (K) at large scales. We derive analytical approaches that estimate K and apply them to the contiguous United States. We then compare these analytical approaches to three-dimensional, national gridded K data products and three transmissivity (T) data products developed from publicly available sources. We evaluate these data products using multiple approaches: comparing their statistics qualitatively and quantitatively and with hydrologic model simulations. Some of these datasets were used as inputs for an integrated hydrologic model of the Upper Colorado River Basin and the comparison of the results with observations was used to further evaluate the K data products. Simulated average daily streamflow was compared to daily flow data from 10 USGS stream gages in the domain, and annually averaged simulated groundwater depths are compared to observations from nearly 2000 monitoring wells. We find streamflow predictions from analytically informed simulations to be similar in relative bias and Spearman's rho to the geologically informed simulations. R-squared values for groundwater depth predictions are close between the best performing analytically and geologically informed simulations at 0.68 and 0.70 respectively, with RMSE values under 10 m. We also show that the analytical approach derived by this study produces estimates of K that are similar in spatial distribution, standard deviation, mean value, and modeling performance to geologically-informed estimates. The results of this work are used to inform a follow-on study that tests additional data-driven approaches in multiple basins within the contiguous United States.

Urban Thirst and Rural Water: The Saga of the Southern Nevada Groundwater Development Project.

In 1989, the Southern Nevada Water Authority (SNWA) launched the Southern Nevada Groundwater Development Project-a bold plan to construct a series of deep wells in east-central Nevada to pump groundwater and send it to the Las Vegas region through 300 miles of pipeline. Before starting work on the project, SNWA conducted an environmental impact study and secured water rights in the valleys. Applications for additional new water rights were filed with Nevada State Engineer on the basis of uncaptured evapotranspiration. The SNWA spent decades and millions of dollars studying the hydrogeology of the region and developing computer models to demonstrate that the project would not have an unduly negative impact on the ecology or water users in the region. The project was opposed by environmental groups, native American tribes, and existing water rights holders. One of the protestants was the Cleveland Ranch in Spring Valley. Using the SNWA's own groundwater model, the ranch argued that the project would result in substantial harm to the ranch's water rights which included springs, wells, and a stream, and that the project would result in perpetual groundwater mining, which is forbidden by Nevada state policy. The Nevada State Engineer approved the project, but the decision was eventually reversed by Seventh District Court, which sided with the ranch and ruled that the project would never be sustainable and is therefore not compatible with Nevada policy. The project was formally abandoned in 2020.

Radiocarbon integrity of dissolved inorganic carbon (DIC) samples stored in plastic and glass bottles: implications for reliable groundwater age dating.

Various approaches based on the natural variations of carbon isotopes (14C and 13C) in dissolved inorganic carbon (DIC) are routinely used to study groundwater dynamics and to estimate recharge rates by deriving groundwater ages. However, differences in 14C activities in groundwater samples collected repeatedly from the same wells and discordantly young 14C groundwater ages compared to noble gases led some authors to question the validity of radiocarbon dating. Poor sampling protocols and storage effects (14C contamination) for radiocarbon analysis are a critical factor in explaining age determination discrepancies. We evaluated the impact of storage protocols on carbon isotope exchange with atmospheric carbon dioxide by comparing glass versus standard plastic field sampling bottles for various storage times before radiocarbon and 13C analyses. The 14C bias after 12 months in pre-evacuated glass vials was minimal and within analytical precision. However, storage of DIC samples in plastic sampling bottles led to marked changes in 14C and 13C contents (up to ∼15 pmC and ∼ 5 ‰, respectively, after 12 months), meaning contamination led to younger groundwater age estimations than it should have been. Protocols for sampling and storing DIC samples for radiocarbon using pre-evacuated glass bottles help avoid atmospheric 14CO2 contamination and microbial activity.

Methane Emissions from Non-producing Oil and Gas Wells and the Potential Role of Seismic Activity: A Case Study in Northeast British Columbia, Canada.

Increasing seismic activity due to fluid injections for oil and gas production may be contributing to leakage along non-producing oil and gas wells and emitting methane, a potent greenhouse gas. However, the extent to which nearby seismicity may drive or exacerbate methane emissions and cause well integrity issues is unknown. Therefore, we analyze field evaluations at 448 non-producing oil and gas wells in Northeast British Columbia (NEBC) and geospatially analyze oil and gas well and fluid injection data alongside locations of 3515 earthquakes from 2001 to 2021 and 130 faults. Through analysis of ground and helicopter-based field evaluations of non-producing wells in NEBC, we show that methane emission rates of non-producing wells average at 8301 mg/h/well but vary by 10 orders of magnitude. We find that higher methane emission rates (milligrams of methane/h/well) are observed at wells with larger flowing pressures at the wellhead during completion (kPa) and with shorter distances (m) to earthquakes, particularly at plugged wells. These results imply that seismicity may increase the likelihood of non-producing well integrity issues and methane leakage, thereby also exacerbating groundwater contamination and environmental degradation risks.

E. coli contamination of drinking water sources in rural and urban settings: an analysis of 38 nationally representative household surveys (2014-2021).

The world is not on track to achieve universal access to safely managed water by 2030, and access is substantially lower in rural areas. This Sustainable Development Goal target and many other global indicators rely on the classification of improved water sources for monitoring access. We aimed to investigate contamination in drinking water sources, comparing improved and unimproved sources in urban and rural settings. We used data from Multiple Indicator Cluster Surveys, which tested samples from the household water source and a glass of water for Escherichia coli contamination across 38 countries. Contamination was widespread and alarmingly high in almost all countries, settings, and water sources, with substantial inequalities between and within countries. Water contamination was found in 51.7% of households at the source and 70.8% in the glass of water. Some improved sources (e.g., protected wells and rainwater) were as likely to be contaminated as unimproved sources. Some sources, like piped water, were considerably more likely to be contaminated in rural than urban areas, while no difference was observed for others. Monitoring water contamination along with further investigation in water collection, storage, and source classification is essential and must be expanded to achieve universal access to safely managed water.

A parsimonious approach to predict regions affected by sewer-borne contaminants in urban aquifers.

Leaky urban drainage networks (UDNs) exfiltrating wastewater can contaminate aquifers. Detailed knowledge on spatiotemporal distributions of water-dissolved, sewer-borne contaminants in groundwater is essential to protect urban aquifers and to optimize monitoring systems. We evaluated the effect of UDN layouts on the spreading of sewer-borne contaminants in groundwater using a parsimonious approach. Due to the UDN's long-term leakage behavior and the existence of non-degradable sewer-borne contaminants (equivalent to a conservative and constant contaminant source), we employed a concept of horizontal line sources to mimic the UDN layout. This does not require the consideration of bio-degradation processes or temporal delay and effectively bypasses the vadose zone, thus reducing computational requirements associated with a full simulation of leakages. We used a set of synthetic leakage scenarios which were generated using fractals and are based on a real-world UDN layout. We investigated the effects of typical leakage rates, varying groundwater flow directions, and UDN's layouts on the shape of the contaminant plume, disregarding the resulted concentration. Leakage rates showed minimal effects on the total covered plume area, whereas 89% of the variance of the plume's geometry is explained by both the UDN's layout (e.g., length and level of complexity) and groundwater flow direction. We demonstrated the potential of applying this approach to identify possible locations of groundwater observation wells using a real UDN layout. This straightforward and parsimonious method can serve as an initial step to strategically identify optimal monitoring systems locations within urban aquifers, and to improve sewer asset management at city scale.

Changes in paleo-groundwater levels revealed by water wells and their relationship with climate variations in imperial Southern China.

Based on records of the bottom elevations of 511 ancient water wells from published archaeological reports, we reconstructed the paleo-groundwater levels (PGWL) in urban areas of Chengdu, Changsha, Nanjing, Suzhou, Suqian, Yancheng, Fuzhou, and Guangzhou cities in southern China. Our PGWL reconstruction shows that PGWL varied in two patterns. In the inland monsoon region (Chengdu and Changsha), there was a low PGWL in Jin (AD 266-420) and South Song (AD 1127-1279), and a high PGWL in Tang (AD 618-907) and Ming (AD 1368-1644). In the coastal region (Yancheng, Fuzhou, and Guangzhou), there was a low PGWL in Jin (AD 266-420) and Ming (AD 1368-1644) but a high PGWL in Tang (AD 618-907) and Song (AD 960-1279). Via cross-wavelet transform and wavelet transform coherence analyses, we found that monsoon and temperature significantly drove the PGWL fluctuations at the inter-centennial scale. East Asian Summer Monsoon-induced precipitation has continuously affected cities in the inland monsoon area represented by Chengdu and Changsha over the past 2,500 years. It has also intermittently affected Nanjing and Suzhou when EASM intensified. In parallel, temperature influenced the PGWL in coastal cities such as Yancheng, Fuzhou, and Guangzhou via the changes in the sea level. Also, the temperature affected the PGWL in relatively inland cities during climatic anomalies such as the Medieval Warm Period and Little Ice Age. This study demonstrates the value of archaeological records in learning how climatic factors influence the PGWL variation and its mechanism.

Assessing the impact of precipitation on hardrock aquifer system using standard precipitation index and groundwater resilience index: a case study of Purulia, West Bengal, India.

Groundwater stored in the aquifers provides water security during natural hazards, e.g. clean water access during floods and droughts. Groundwater drought, a phenomenon closely linked with rainfall (climate) variability, is less researched, especially in India. This study aims to detect precipitation and groundwater droughts and comprehend the groundwater response to long-term precipitation trends (25 years). As a case study, the drought-affected and groundwater-depleted Purulia district in West Bengal, India, which is a part of the Chotanagpur plateau, was selected. Precipitation and groundwater droughts (in aquifer types of shallow, moderate and deep) are detected using the Standard Precipitation Index (SPI) and Groundwater Resilience Index (GRI). During the 25 year study period (1996-2020), Purulia had 13 (52%) rainfall deficiency years, with an annual average rainfall of 1382 mm. SPI detected four severe droughts and the most severe occurring in 2010-2011 (1.50). GRI found that aquifermedium had a 71% [Formula: see text] conditions and are the most resilient and aquiferdeep experienced maximum extreme drought events and is the most stressed. The cross-correlation coefficients (CCCs) between rainfall and groundwater is moderate in deep, shallow, and medium aquifers, with CCCs - 0.43, - 0.59, and - 0.49, respectively. Positive CCCs are found for seasonal lags of - 3, - 4, and - 7. The study found that during the monsoon, average depth to groundwater level is 1 - 4 m and it drops to 8 - 10 m during the lean period, more than 85% of wells are vulnerable to extreme droughts (SPI > 1.5), aquifer's response to rainfall is aquifershallow > aquifermoderate > aquiferdeep, and aquifer's may be arranged as aquifermoderate > aquifershallow > aquiferdeep depending on their drought resistance. This study, with the use of statistical tools and long term data, will aid in the management of groundwater at varying depths by creating basis for understanding the groundwater response to rainfall events.