Spatial and seasonal controls on dissolved organic matter composition in shallow aquifers under the rapidly developing city of Patna, India.

The distribution and composition of dissolved organic matter (DOM) affects numerous (bio)geochemical processes in environmental matrices including groundwater. This study reports the spatial and seasonal controls on the distribution of groundwater DOM under the rapidly developing city of Patna, Bihar (India). Major DOM constituents were determined from river and groundwater samples taken in both pre- and post-monsoon seasons in 2019, using excitation-emission matrix (EEM) fluorescence spectroscopy. We compared aqueous fluorescent DOM (fDOM) composition to satellite-derived land use data across the field area, testing the hypothesis that the composition of groundwater DOM, and particularly the components associated with surface-derived ingress, may be controlled, in part, by land use. In the pre-monsoon season, the prominence of tryptophan-like components likely generated from recent biological activity overwhelmed the humic-like and tyrosine-like fluorescence signals. Evidence from fluorescence data suggest groundwater in the post-monsoon season is composed of predominantly i) plant-derived matter and ii) anthropogenically influenced DOM (e.g. tryptophan-like components). Organic tracers, as well as Eh and Cl-, suggest monsoonal events mobilise surface-derived material from the unsaturated zone, causing dissolved organic carbon (DOC) of more microbial nature to infiltrate to >100 m depth. A correlation between higher protein:humic-like fluorescence and lower vegetation index (NDVI), determined from satellite-based land use data, in the post-monsoon season, indicates the ingression of wastewater-derived OM in groundwater under the urban area. Attenuated protein:humic-like fluorescence in groundwater close to the river points towards the mixing of groundwater and river water. This ingress of surface-derived OM is plausibly exacerbated by intensive groundwater pumping under these areas. Our approach to link the composition of aqueous organics with land use could easily be adapted for similar hydrogeochemical settings to determine the factors controlling groundwater DOM composition in various contexts.

Emerging organic contaminants in the River Ganga and key tributaries in the middle Gangetic Plain, India: Characterization, distribution & controls.

The presence and distribution of emerging organic contaminants (EOCs) in freshwater environments is a key issue in India and globally, particularly due to ecotoxicological and potential antimicrobial resistance concerns. Here we have investigated the composition and spatial distribution of EOCs in surface water along a ∼500 km segment of the iconic River Ganges (Ganga) and key tributaries in the middle Gangetic Plain of Northern India. Using a broad screening approach, in 11 surface water samples, we identified 51 EOCs, comprising of pharmaceuticals, agrochemicals, lifestyle and industrial chemicals. Whilst the majority of EOCs detected were a mixture of pharmaceuticals and agrochemicals, lifestyle chemicals (and particularly sucralose) occurred at the highest concentrations. Ten of the EOCs detected are priority compounds (e.g. sulfamethoxazole, diuron, atrazine, chlorpyrifos, perfluorooctane sulfonate (PFOS), perfluorobutane sulfonate, thiamethoxam, imidacloprid, clothianidin and diclofenac). In almost 50% of water samples, sulfamethoxazole concentrations exceeded predicted no-effect concentrations (PNECs) for ecological toxicity. A significant downstream reduction in EOCs was observed along the River Ganga between Varanasi (Uttar Pradesh) and Begusarai (Bihar), likely reflecting dilution effects associated with three major tributaries, all with considerably lower EOC concentrations than the main Ganga channel. Sorption and/or redox controls were observed for some compounds (e.g. clopidol), as well as a relatively high degree of mixing of EOCs within the river. We discuss the environmental relevance of the persistence of several parent compounds (notably atrazine, carbamazepine, metribuzin and fipronil) and associated transformation products. Associations between EOCs and other hydrochemical parameters including excitation emission matrix (EEM) fluorescence indicated positive, significant, and compound-specific correlations between EOCs and tryptophan-, fulvic- and humic-like fluorescence. This study expands the baseline characterization of EOCs in Indian surface water and contributes to an improved understanding of the potential sources and controls on EOC distribution in the River Ganga and other large river systems.

Environmental tracers and groundwater residence time indicators reveal controls of arsenic accumulation rates beneath a rapidly developing urban area in Patna, India.

Groundwater security is a pressing environmental and societal issue, particularly due to significantly increasing stressors on water resources, including rapid urbanization and climate change. Groundwater arsenic is a major water security and public health challenge impacting millions of people in the Gangetic Basin of India and elsewhere globally. In the rapidly developing city of Patna (Bihar) in northern India, we have studied the evolution of groundwater chemistry under the city following a three-dimensional sampling framework of multi-depth wells spanning the central urban zone in close proximity to the River Ganges (Ganga) and transition into peri-urban and rural areas outside city boundaries and further away from the river. Using inorganic geochemical tracers (including arsenic, iron, manganese, nitrate, nitrite, ammonium, sulfate, sulfide and others) and residence time indicators (CFCs and SF6), we have evaluated the dominant hydrogeochemical processes occurring and spatial patterns in redox conditions across the study area. The distribution of arsenic and other redox-sensitive parameters is spatially heterogenous, and elevated arsenic in some locations is consistent with arsenic mobilization via reductive dissolution of iron hydroxides. Residence time indicators evidence modern (<~60-70 years) groundwater and suggest important vertical and lateral flow controls across the study area, including an apparent seasonal reversal in flow regimes near the urban center. An overall arsenic accumulation rate is estimated to be ~0.003 ± 0.003 µM.yr-1 (equivalent to ~0.3 ± 0.2 µg.yr-1), based on an average of CFC-11, CFC-12 and SF6-derived models, with the highest rates of arsenic accumulation observed in shallow, near-river groundwaters also exhibiting elevated concentrations of nutrients including ammonium. Our findings have implications on groundwater management in Patna and other rapidly developing cities, including potential future increased groundwater vulnerability associated with surface-derived ingress from large-scale urban abstraction or in higher permeability zones of river-groundwater connectivity.

Emerging organic contaminants in groundwater under a rapidly developing city (Patna) in northern India dominated by high concentrations of lifestyle chemicals.

Aquatic pollution from emerging organic contaminants (EOCs) is of key environmental importance in India and globally, particularly due to concerns of antimicrobial resistance, ecotoxicity and drinking water supply vulnerability. Here, using a broad screening approach, we characterize the composition and distribution of EOCs in groundwater in the Gangetic Plain around Patna (Bihar), as an exemplar of a rapidly developing urban area in northern India. A total of 73 EOCs were detected in 51 samples, typically at ng.L-1 to low µg.L-1 concentrations, relating to medical and veterinary, agrochemical, industrial and lifestyle usage. Concentrations were often dominated by the lifestyle chemical and artificial sweetener sucralose. Seventeen identified EOCs are flagged as priority compounds by the European Commission, World Health Organisation and/or World Organisation for Animal Health: namely, herbicides diuron and atrazine; insecticides imidacloprid, thiamethoxam, clothianidin and acetamiprid; the surfactant perfluorooctane sulfonate (and related perfluorobutane sulfonate, perfluorohexane sulfonate, perfluorooctanoic acid and perfluoropentane sulfonate); and medical/veterinary compounds sulfamethoxazole, sulfanilamide, dapson, sulfathiazole, sulfamethazine and diclofenac. The spatial distribution of EOCs varies widely, with concentrations declining with depth, consistent with a strong dominant vertical flow control. Groundwater EOC concentrations in Patna were found to peak within ∼10 km distance from the River Ganges, indicating mainly urban inputs with some local pollution hotspots. A heterogeneous relationship between EOCs and population density likely reflects confounding factors including varying input types and controls (e.g. spatial, temporal), wastewater treatment infrastructure and groundwater abstraction. Strong seasonal agreement in EOC concentrations was observed. Co-existence of limited transformation products with associated parent compounds indicate active microbial degradation processes. This study characterizes key controls on the distribution of groundwater EOCs across the urban to rural transition near Patna, as a rapidly developing Indian city, and contributes to the wider understanding of the vulnerability of shallow groundwater to surface-derived contamination in similar environments.

Exposure-time based modeling of nonlinear reactive transport in porous media subject to physical and geochemical heterogeneity.

Transport of reactive solutes in groundwater is affected by physical and chemical heterogeneity of the porous medium, leading to complex spatio-temporal patterns of concentrations and reaction rates. For certain cases of bioreactive transport, it could be shown that the concentrations of reactive constituents in multi-dimensional domains are approximately aligned with isochrones, that is, lines of identical travel time, provided that the chemical properties of the matrix are uniform. We extend this concept to combined physical and chemical heterogeneity by additionally considering the time that a water parcel has been exposed to reactive materials, the so-called exposure time. We simulate bioreactive transport in a one-dimensional domain as function of time and exposure time, rather than space. Subsequently, we map the concentrations to multi-dimensional heterogeneous domains by means of the mean exposure time at each location in the multi-dimensional domain. Differences in travel and exposure time at a given location are accounted for as time difference. This approximation simplifies reactive-transport simulations significantly under conditions of steady-state flow when reactions are restricted to specific locations. It is not expected to be exact in realistic applications because the underlying assumption, such as neglecting transverse mixing altogether, may not hold. We quantify the error introduced by the approximation for the hypothetical case of a two-dimensional, binary aquifer made of highly-permeable, non-reactive and low-permeable, reactive materials releasing dissolved organic matter acting as electron donor for aerobic respiration and denitrification. The kinetically controlled reactions are catalyzed by two non-competitive bacteria populations, enabling microbial growth. Even though the initial biomass concentrations were uniform, the interplay between transport, non-uniform electron-donor supply, and bio-reactions led to distinct spatial patterns of the two types of biomass at late times. Results obtained by mapping the exposure-time based results to the two-dimensional domain are compared with simulations based on the two-dimensional, spatially explicit advection-dispersion-reaction equation. Once quasi-steady state has been reached, we find a good agreement in terms of the chemical-compound concentrations between the two approaches inside the reactive zones, whereas the exposure-time based model is not able to capture reactions occurring in the zones with zero electron-donor release. We conclude that exposure-time models provide good approximations of nonlinear bio-reactive transport when transverse mixing is not the overall controlling process and all reactions are essentially restricted to distinct reactive zones.

Using travel times to simulate multi-dimensional bioreactive transport in time-periodic flows.

In travel-time models, the spatially explicit description of reactive transport is replaced by associating reactive-species concentrations with the travel time or groundwater age at all locations. These models have been shown adequate for reactive transport in river-bank filtration under steady-state flow conditions. Dynamic hydrological conditions, however, can lead to fluctuations of infiltration velocities, putting the validity of travel-time models into question. In transient flow, the local travel-time distributions change with time. We show that a modified version of travel-time based reactive transport models is valid if only the magnitude of the velocity fluctuates, whereas its spatial orientation remains constant. We simulate nonlinear, one-dimensional, bioreactive transport involving oxygen, nitrate, dissolved organic carbon, aerobic and denitrifying bacteria, considering periodic fluctuations of velocity. These fluctuations make the bioreactive system pulsate: The aerobic zone decreases at times of low velocity and increases at those of high velocity. For the case of diurnal fluctuations, the biomass concentrations cannot follow the hydrological fluctuations and a transition zone containing both aerobic and obligatory denitrifying bacteria is established, whereas a clear separation of the two types of bacteria prevails in the case of seasonal velocity fluctuations. We map the 1-D results to a heterogeneous, two-dimensional domain by means of the mean groundwater age for steady-state flow in both domains. The mapped results are compared to simulation results of spatially explicit, two-dimensional, advective-dispersive-bioreactive transport subject to the same relative fluctuations of velocity as in the one-dimensional model. The agreement between the mapped 1-D and the explicit 2-D results is excellent. We conclude that travel-time models of nonlinear bioreactive transport are adequate in systems of time-periodic flow if the flow direction does not change.

On the validity of travel-time based nonlinear bioreactive transport models in steady-state flow.

Travel-time based models simplify the description of reactive transport by replacing the spatial coordinates with the groundwater travel time, posing a quasi one-dimensional (1-D) problem and potentially rendering the determination of multidimensional parameter fields unnecessary. While the approach is exact for strictly advective transport in steady-state flow if the reactive properties of the porous medium are uniform, its validity is unclear when local-scale mixing affects the reactive behavior. We compare a two-dimensional (2-D), spatially explicit, bioreactive, advective-dispersive transport model, considered as "virtual truth", with three 1-D travel-time based models which differ in the conceptualization of longitudinal dispersion: (i) neglecting dispersive mixing altogether, (ii) introducing a local-scale longitudinal dispersivity constant in time and space, and (iii) using an effective longitudinal dispersivity that increases linearly with distance. The reactive system considers biodegradation of dissolved organic carbon, which is introduced into a hydraulically heterogeneous domain together with oxygen and nitrate. Aerobic and denitrifying bacteria use the energy of the microbial transformations for growth. We analyze six scenarios differing in the variance of log-hydraulic conductivity and in the inflow boundary conditions (constant versus time-varying concentration). The concentrations of the 1-D models are mapped to the 2-D domain by means of the kinematic (for case i), and mean groundwater age (for cases ii & iii), respectively. The comparison between concentrations of the "virtual truth" and the 1-D approaches indicates extremely good agreement when using an effective, linearly increasing longitudinal dispersivity in the majority of the scenarios, while the other two 1-D approaches reproduce at least the concentration tendencies well. At late times, all 1-D models give valid approximations of two-dimensional transport. We conclude that the conceptualization of nonlinear bioreactive transport in complex multidimensional domains by quasi 1-D travel-time models is valid for steady-state flow fields if the reactants are introduced over a wide cross-section, flow is at quasi steady state, and dispersive mixing is adequately parametrized.