Search

With the help of a physically based recharge-groundwater flow model and robust detrended fluctuation analysis (r-DFAn), the effect of local (catchment-scale) forcing on groundwater levels’ scaling behavior in a riparian aquifer in Wallingford, UK, is investigated. The local forcings investigated in this study include the rainfall’s temporal scaling behavior (which is simulated by changing rainfall’s intermittency parameter in a β-lognormal multiplicative random cascade model), the aquifer’s physical parameters (saturated hydraulic conductivity, specific yield, the empirical coefficients of the water retention curve, and the river stage’s scaling behavior). Groundwater level’s scaling behaviour was found to be most sensitive to rainfall’s fractal behaviour. Additionally, there is preliminary evidence suggesting that changes to the rainfall’s local scaling behaviour (i.e., change to the series’ scaling that induces crossovers) affects the groundwater’s and the recharge’s local scaling behaviour.

This paper details the design, development and evaluation of a 3D printed, rechargeable, Ag-AgCl electrode to measure Self-Potential (SP) in laboratory experiments. The challenge was to make a small, cheap, robust and stable electrode that could be used in a wide range of applications. The new electrodes are shown to offer comparable performance with custom-machined laboratory standards, and the inclusion of 3D printing (both Fused Filament Fabrication (FFF) and stereolithography (SLA)) makes them more versatile and less expensive than laboratory standards. The devices have been used in both low-pressure experiments using beadpacks, and high-pressure experiments using natural rock samples. Designs are included for both male and female connections to laboratory equipment.

The aim of this paper is to classify and segment roofs using vertical aerial imagery to generate three-dimensional (3D) models. Such models can be used, for example, to evaluate the rainfall runoff from properties for rainwater harvesting and in assessing solar energy and roof insulation options. Aerial orthophotos and building footprints are used to extract individual roofs and bounding boxes, which are then fed into one neural network for classification and then another for segmentation. The approach initially implements transfer learning on a pre-trained VGG16 model. The first step achieves an accuracy of 95.39% as well as a F1 score of 95%. The classified images are segmented using a fully convolutional network semantic segmentation model. The mask of the segmented roof planes is used to extract the coordinates of the roof edges and the nexus points using the Harris corner detector algorithm. The coordinates of the corners are then used to plot a 3D Level of Detail 2 (LOD2) representation of the building and the roof height is determined by calculating the maximum and minimum height of a Digital Surface Model LiDAR point cloud and known building height data. Subsequently the wireframe plot can be used to compute the roof area. This model achieved an accuracy of 80.2%, 96.1%, 96.0%, 85.1% and 91.1% for flat, hip, gable, cross-hip and mansard roofs, respectively.

Petroleum hydrocarbon contamination is a global problem which can cause long-term environmental damage and impacts water security. Natural source zone depletion (NSZD) is the natural degradation of such contaminants. Chemotaxis is an aspect of NSZD which is not fully understood, but one that grants microorganisms the ability to alter their motion in response to a chemical concentration gradient potentially enhancing petroleum NSZD mass removal rates. This study investigates the distribution of potentially chemotactic and hydrocarbon degrading microbes (CD) across the water table of a legacy petroleum hydrocarbon site near Perth, Western Australia in areas impacted by crude oil, diesel and jet fuel. Core samples were recovered and analysed for hydrocarbon contamination using gas chromatography. Predictive metagenomic profiling was undertaken to infer functionality using a combination of 16 S rRNA sequencing and PICRUSt2 analysis. Naphthalene contamination was found to significantly increase the occurrence of potential CD microbes, including members of the Comamonadaceae and Geobacteraceae families, which may enhance NSZD. Further work to explore and define this link is important for reliable estimation of biodegradation of petroleum hydrocarbon fuels. Furthermore, the outcomes suggest that the chemotactic parameter within existing NSZD models should be reviewed to accommodate CD accumulation in areas of naphthalene contamination, thereby providing a more accurate quantification of risk from petroleum impacts in subsurface environments, and the scale of risk mitigation due to NSZD.

Self‐potential (SP) measurements can be used to characterize and monitor, in real‐time, fluid movement and behavior in the subsurface. The electrochemical exclusion‐diffusion (EED) potential, one component of SP, arises when concentration gradients exist in porous media. Such concentration gradients are of concern in coastal and contaminated aquifers and oil and gas reservoirs. It is essential that estimates of EED potential are made prior to conducting SP investigations in complex environments with heterogeneous geology and salinity contrasts, such as the UK Chalk coastal aquifer. Here we report repeatable laboratory estimates of the EED potential of chalk and marls using natural groundwater (GW), seawater (SW), deionized (DI) water, and 5 M NaCl. In all cases, the EED potential of chalk was positive (using a GW/SW concentration gradient the EED potential was ca. 14 to 22 mV), with an increased deviation from the diffusion limit at the higher salinity contrast. Despite the relatively small pore size of chalk (ca. 1 μm), it is dominated by the diffusion potential and has a low exclusion efficiency, even at large salinity contrasts. Marl samples have a higher exclusion efficiency which is of sufficient magnitude to reverse the polarity of the EED potential (using a GW/SW concentration gradient the EED potential was ca. −7 to −12 mV) with respect to the chalk samples. Despite the complexity of the natural samples used, the method produced repeatable results. We also show that first order estimates of the exclusion efficiency can be made using SP logs, supporting the parameterization of the model reported in Graham et al. (2018, https://doi.org/10.1029/2018WR022972), and that derived values for marls are consistent with the laboratory experiments, while values derived for hardgrounds based on field data indicate a similarly high exclusion efficiency. While this method shows promise in the absence of laboratory measurements, more rigorous estimates should be made where possible and can be conducted following the experimental methodology reported here.

Reliable estimates of the sustainable yield of supply boreholes are critical to ensure that groundwater resources are managed sustainably. Sustainable yields are dependent on the pumped groundwater level in a borehole, how this relates to vertical aquifer heterogeneity, and features of the borehole itself. This paper presents a 3D radial flow model (SPIDERR), based on the Darcy-Forchheimer equation, for simulating the groundwater level response in supply boreholes in unconfined, heterogeneous aquifers. The model provides a tool for investigating the causes of non-linear behaviour in abstraction boreholes, which can have a significant impact on sustainable yields. This is demonstrated by simulating a variable-rate pumping test in a Chalk abstraction borehole. The application suggests the non-linear response to pumping is due to a combination of factors: a reduction in well storage with depth due to changes in the borehole diameter, a reduction in hydraulic conductivity with depth, and non-Darcian flow.

Due to the COVID-19 pandemic, citizens of the United Kingdom were required to stay at home for many months in 2020. In the weeks before and months following lockdown, including when it was not being enforced, citizens were advised to stay at home where possible. As a result, in a megacity such as London, where long-distance commuting is common, spatial and temporal changes to patterns of water demand are inevitable. This, in turn, may change where people’s waste is treated and ultimately impact the in-river quality of effluent receiving waters. To assess large scale impacts, such as COVID-19, at the city scale, an integrated modelling approach that captures everything between households and rivers is needed. A framework to achieve this is presented in this study and used to explore changes in water use and the associated impacts on wastewater treatment and in-river quality as a result of government and societal responses to COVID-19. Our modelling results revealed significant changes to household water consumption under a range of impact scenarios, however, significant impacts on pollutant concentrations in household wastewater were only seen in central London. Pollutant concentrations in rivers were shown to be most sensitive in the tributaries of the River Thames, highlighting the vulnerability of smaller rivers and the important role that they play in diluting pollution. Ammonia and phosphates were found to be the pollutants that rivers were most sensitive to because their main source in urban rivers is domestic wastewater that was significantly altered during the imposed mobility restrictions. A model evaluation showed that we can accurately validate individual model components (i.e., water demand generator). A model evaluation showed that we can accurately validate individual model components (i.e., water demand generator) and emphasised need for continuous water quality measurements. Ultimatly, the work provides a basis for further developments of water systems integration approaches to project changes under never-before seen scenarios.

Salinity is a pervasive problem in the coastal low-lying area of the deltas including Bangladesh located in one of the largest delta, Ganges-Brahmaputra-Meghna delta. This delta is susceptible to episodic cyclones since it is nearly every 3 years hit by tropical cyclones in the early monsoonal season (April to June) or the early dry season (October to November). These successive cyclones associated with low-lying reclaimed lands that trigger extensive flooding and result in excess salinity in soil and surface water, which have led to low agricultural productivity. Salinity in drinking water causes negative effects on human health such as cardiovascular disease. A fully coupled surface-subsurface model was used to investigate the impact of the episodic cyclonic surges on the drinking pond and groundwater salinities in the coastal reclaimed lands of Dacope Upazila in southwest Bangladesh. We considered 5 scenarios: a cyclone hit the land in the monsoon season with remediation (clean-up the pond at (1) 7 days, (2) 3 months), a cyclone hit the land in the dry season with remediation (clean-up the pond at (3) 7 days (4) 9 months) and (5) the recurrent intervals of cyclones hit the land every 8 years. The hydrological parameters were calibrated from the fieldwork at DAB site in using in situ field observations. The results show that the episodic cyclones caused inevitable salinity to near-surface groundwater, and in pond water because of post-event delayed emptying of ponds and reversal of hydraulic head gradient. However, rapid remediation after a surge event may help avoid serious salinity in drinking water. The result of scenario 5 indicates that near-surface groundwater salinity progressively developed and move downward over time. The episodic surge events might be one of the reasons that cause shallow groundwater salinity in coastal Bangladesh. This study improves our understanding of salinization processes and how to manage drinking water ponds after a storm surge induced flooding in deltaic coastal settings.

Saltwater intrusion (SWI) in coastal aquifers poses a significant hazard to freshwater security for many of the world’s population centers. SWI is challenging to monitor and model due to the physical complexity of real aquifers. Self-Potential (SP) has been an important method for monitoring the subsurface for many years. Previous studies have suggested that borehole measurements of SP could be used to identify saline interface movement and provide advance warning of imminent saline breakthrough at an abstraction borehole. SP produced during SWI comprises the combined effects of electro-kinetic potential, arising from transport of excess charge in response to water potential (head) gradients, and exclusion-diffusion potential, arising from transport of excess charge in response to ion (salt) concentration gradients. SP can have advantages over other geophysical methods, such as electrical resistivity tomography and borehole fluid electrical conductivity measurements, because the effect of moving saltwater fronts can be determined using a relatively small number of localized probes.We quantitatively investigate the relationship between SP and SWI using experimental and numerical modelling with the aim of reproducing experimentally measured SP response via simulation. Building on well-established methods, a novel laboratory setup has been developed to optically monitor SWI in a thin homogenous aquifer while simultaneously recording SP data at multiple probe points. A Matlab solver is used to calculate SP data from simulated hydrodynamic SWI data computed by the fixed-grid finite element software SUTRA. Similarly, finite element SWI simulations using adaptive meshing are carried out using the IC-FERST software, which directly computes hydrodynamic and SP solutions. We compare these numerical results with experimental data and show similarity in SP signal trends as functions of brine movement near probe locations. We conclude with a discussion of the merits of SP modelling and its suitability for interpreting SP signals for monitoring and characterization of saltwater intrusion in coastal aquifers.

A new multi-scale groundwater modelling methodology is presented to simulate pumped water levels in abstraction boreholes within regional groundwater models, providing a robust tool for assessing the sustainable yield of supply boreholes and improving our understanding of groundwater availability during drought. A 3D borehole-scale model, which solves the Darcy-Forchheimer equation in cylindrical co-ordinates to simulate both linear and non-linear radial flow to a borehole in a heterogeneous aquifer, is embedded within a Cartesian grid, using a hybrid radial-Cartesian finite difference method. The local-scale model is coupled to a regional groundwater model, ZOOMQ3D, using the OpenMI model linkage software, providing a flexible and efficient tool for assessing the behaviour of a groundwater source within its regional hydrogeological context during historic droughts and under climate change. The advantages of the new method are demonstrated through application to a Chalk supply borehole in the UK.

In order to enable community groups and other interested parties to evaluate the effects of flood management, water conservation and other hydrological issues, better localised mapping is required. Although some maps are publicly available many are behind paywalls, especially those with three dimensional features. In this study London is used as a test case to evaluate, machine learning and rules-based approaches with opensource maps and LiDAR data to create more accurate representations (LOD2) of small-scale areas. Machine learning is particularly well suited to the recognition of local repetitive features like building roofs and trees, while roads can be identified and mapped best using a faster rules-based approach. In order to create a useful LOD2 representation, a user interface, processing rules manipulation and assumption editor have all been incorporated. Features like randomly assigning sub terrain features (basements) - using Monte-Carlo methods - and artificial sewage representation enable the user to grow these models from opensource data into useful model inputs. This project is aimed at local scale hydrological modelling, rainfall runoff analysis and other local planning applications. The goal is to provide turn-key data processing for small scale modelling, which should help advance the installation of SuDs and other water management solutions, as well as having broader uses. The method is designed to enable fast and accurate representations of small-scale features (1 hectare to 1km2), with larger scale applications planned for future work. This work forms part of the CAMELLIA project (Community Water Management for a Liveable London) and aims to provide useful tools for local scale modeller and possibly the larger scale industry/scientific user.

The Permo-Triassic Sandstone aquifers of the Eden Valley, Cumbria UK, are a key water resource for public water supply in NW England as well as local agriculture and industries. Permo-Triassic Sandstone aquifers are characterised as having large storage capacities and moderate transmissivities, however, in the Eden Valley these characteristics vary greatly on a range of scales i.e. granulation seams (deformation bands) that are millimetres thick but have been shown to extend for hundreds of metres on analogous sandstones; silicified layers that are several metres thick and extending 10s to 100s of metres laterally; and lithological variation and faulting have been shown to juxtapose hydrogeological units with different hydraulic properties. Complex heterogeneous superficial deposits overlay 75% of the Permo-Triassic Sandstone aquifers and comprise glacial till, glacio-fluvial outwash deposits, river terrace deposits and alluvium. The lateral and vertical continuity of these superficial deposits is highly uncertain. The complex geological and superficial deposits in the Eden Valley impose a control on flow processes and impact sub-surface runoff. Specifically, lenses of high conductivity sands and gravels within low conductivity clay till deposits coupled with the presence of low conductivity strata at ground level suggests that indirect recharge is an important sub-surface runoff component. Therefore, the magnitude and location of recharge to the Permo-Triassic Sandstone aquifers is highly uncertain. Published recharge estimates rely on baseflow separation techniques and thus do not distinguish between indirect and direct recharge. This highlights the uncertainty regarding the sub-surface flow processes active in the Eden Valley. A methodology for characterizing the surface water – groundwater interaction spatially and temporally in an ungauged upland sub-catchment is presented. A non-invasive approach has been implemented to investigate the relationship between the surface water and groundwater systems in the Eden Valley. This involved the design and installation of low-cost ultrasonic sensors that measure stream stage. The sensors have been installed at key locations within sub-catchments that incorporate limestone pavements, geological contacts and along fault trends in the headwaters of the Eden Valley. Flow gauging has been conducted along the reach of these streams to investigate the spatial variation in discharge. Data from the low-cost sensors and flow gauging have been used to estimate the magnitude of volumetric water exchange between the surface water and groundwater systems, as well as characterise this relationship spatially and temporally. The thickness and composition of the superficial deposits along these stream reaches will be investigated via passive seismic survey. The superficial investigation and the volumetric water balance will be used to estimate indirect recharge in the upper Eden catchment. The results of which will be compared to localised recharge estimates calculated from groundwater level timeseries. This comparison will indicate the importance of indirect recharge within sub-surface runoff processes. This ongoing research is a vital step in quantifying the relationship between the surface water and groundwater systems in a complex upland catchment. A knowledge of the active sub-surface runoff processes highlighted are key for reliably assessing the long-term security of groundwater resources in the Eden Valley.

Leaching of nitrogen-containing compounds (e.g., ammonia, nitrate) from pit latrines and seepage tanks into groundwater may pose health risks, given that groundwater is a significant source for drinking water in many low-income countries. In this study, three communities within Kathmandu, Nepal (Manohara, Kupondole, and Lokanthali) were visited to investigate the impact of pit latrines on groundwater quality, with a focus on understanding the fate of nitrogen-containing compounds specifically. Well water samples were analyzed over two seasons (wet and dry) for their nitrogen content, dissolved oxygen (DO), chemical oxidation demand (COD), and oxidation-reduction potential (ORP), and samples collected from within the nearby pits were also analyzed to determine nitrogen content and COD. Hand dug wells were found to be more likely receptors of contamination than tube wells, as expected, with inter-well variations related to the relative redox conditions in the wells. Increased pit-emptying frequency was related to lower levels of nitrogen in the latrines and in the nearest wells, suggesting this may be an effective strategy for reducing the risks of groundwater contamination in such settings, all else being equal.

Two years of self‐potential (SP) measurements were made in a monitoring borehole in the coastal UK Chalk aquifer. The borehole SP data showed a persistent gradient with depth, and temporal variations with a tidal power spectrum consistent with ocean tides. No gradient with depth was observed at a second coastal monitoring borehole ca. 1 km further inland, and no gradient or tidal power spectrum were observed at an inland site ca. 80 km from the coast. Numerical modeling suggests that the SP gradient recorded in the coastal monitoring borehole is dominated by the exclusion‐diffusion potential, which arises from the concentration gradient across a saline front in close proximity to, but not intersecting, the base of the borehole. No such saline front is present at the two other monitoring sites. Modeling further suggests that the ocean tidal SP response in the borehole, measured prior to breakthrough of saline water, is dominated by the exclusion‐diffusion potential across the saline front, and that the SP fluctuations are due to the tidal movement of the remote front. The electrokinetic potential, caused by changes in hydraulic head across the tide, is one order of magnitude too small to explain the observed SP data. The results suggest that in coastal aquifers, the exclusion‐diffusion potential plays a dominant role in borehole SP when a saline front is nearby. The SP gradient with depth indicates the close proximity of the saline front to the borehole and changes in SP at the borehole reflect changes in the location of the saline front. Thus, SP monitoring can be used to facilitate more proactive management of abstraction and saline intrusion in coastal aquifers.

Crop simulation models are valuable tools for quantifying crop yield response to water, and for devising strategies to improve agricultural water management. However, applicability of the majority of crop models is limited greatly by a failure to provide open-access to model source code. In this study, we present an open-source version of the FAO AquaCrop model, which simulates efficiently water-limited crop production across diverse environmental and agronomic conditions. Our model, called AquaCrop-OpenSource (AquaCrop-OS), can be run in multiple programming languages and operating systems. Support for parallel execution reduces significantly simulation times when applying the model in large geospatial frameworks, for long-run policy analysis, or for uncertainty assessment. Furthermore, AquaCrop-OS is compliant with the Open Modelling Interface standard facilitating linkage to other disciplinary models, for example to guide integrated water resources planning., This research was funded in part by USDA NIFA Award number 2012-67003-23227 and by a PhD studentship from the Grantham Institute at Imperial College London. Further support was provided by the Daugherty Water for Food Global Institute at the University of Nebraska.

Worldwide, there is growing recognition of the need to reduce agricultural groundwater use in response to rapid rates of aquifer depletion. To date, however, few studies have evaluated how benefits of conservation vary along an aquifer's depletion pathway. To address this question, we develop an integrated modeling framework that couples an agro-economic model of farmers' field-level irrigation decision-making with a borehole-scale groundwater flow model. Unique to this framework is the explicit consideration of the dynamic reductions in well yields that occur as an aquifer is depleted, and how these changes in intraseasonal groundwater supply affect farmers' ability to manage production risks caused by climate variability and, in particular, drought. For an illustrative case study in the High Plains region of the United States, we apply our model to analyze the value of groundwater conservation activities for different initial aquifer conditions. Our results demonstrate that there is a range of initial conditions for which reducing pumping will have long-term economic benefits for farmers by slowing reductions in well yields and prolonging the usable lifetime of an aquifer for high-value irrigated agriculture. In contrast, restrictions on pumping that are applied too early or too late will provide limited welfare benefits. We suggest, therefore, that there are ‘windows of opportunity’ to implement groundwater conservation, which will depend on complex feedbacks between local hydrology, climate, crop growth, and economics. This article is protected by copyright. All rights reserved.

Pervasive salinity in soil and water is affecting agricultural yield and the health of millions of delta dwellers in Asia. This is also being exacerbated by climate change through increases in sea level and tropical storm surges. One consequence of this has been a widespread introduction of salt water shrimp farming. Here, we show, using field data and modeling, how changes in climate and land use are likely to result in increased salinization of shallow groundwater in SE Asian mega-deltas. We also explore possible adaptation options. We find that possible future increase of episodic inundation events, combined with salt water shrimp farming, will cause rapid salinization of groundwater in the region making it less suitable for drinking water and irrigation. However, modified land use and water management practices can mitigate the impacts on groundwater, as well as the overlying soil, from future salinization. The study therefore provides guidance for adaptation planning to reduce future salinization in Asian deltas.

Arsenic is a carcinogen known for its acute toxicity to organisms. Geothermal waters are commonly high in arsenic, as shown at the Bjarnarflag Power Plant, Iceland (∼224 μg/kg of solvent). Development of geothermal energy requires adequate disposal of arsenic-rich waters into groundwater/geothermal systems. The outcome of arsenic transport models that assess the effect of geothermal effluent on the environment and ecosystems may be influenced by the sensitivity of hydraulic parameters. However, previous such studies in Iceland do not consider the sensitivity of hydraulic parameters and thereby the interpretations remain unreliable. Here we used the Lake Mývatn basaltic aquifer system as a case study to identify the sensitive hydraulic parameters and assess their role in arsenic transport. We develop a one-dimensional reactive transport model (PHREEQC ver. 2.), using geochemical data from Bjarnarflag, Iceland.In our model, arsenite (H3AsO3) was predicted to be the dominant species of inorganic arsenic in both groundwater and geothermal water. Dilution reduced arsenic concentration below ∼5 μg/kg. Adsorption reduced the residual contamination below ∼0.4 μg/kg at 250 m along transect. Based on our modelling, we found volumetric input to be the most sensitive parameter in the model. In addition, the adsorption strength of basaltic glass was such that the physical hydrogeological parameters, namely: groundwater velocity and longitudinal dispersivity had little influence on the concentration profile. Keywords: Geothermal, Groundwater, Hydrogeochemistry, Arsenic, Modelling

Land surface models (LSMs) are prospective starting points to develop a global hyper-resolution model of the terrestrial water, energy, and biogeochemical cycles. However, there are some fundamental limitations of LSMs related to how meaningfully hydrological fluxes and stores are represented. A diagnostic approach to model evaluation and improvement is taken here that exploits hydrological expert knowledge to detect LSM inadequacies through consideration of the major behavioural functions of a hydrological system: overall water balance, vertical water redistribution in the unsaturated zone, temporal water redistribution, and spatial water redistribution over the catchment's groundwater and surface-water systems. Three types of information are utilized to improve the model's hydrology: (a) observations, (b) information about expected response from regionalized data, and (c) information from an independent physics-based model. The study considers the JULES (Joint UK Land Environmental Simulator) LSM applied to a deep-groundwater chalk catchment in the UK. The diagnosed hydrological limitations and the proposed ways to address them are indicative of the challenges faced while transitioning to a global high resolution model of the water cycle.

Monitoring of self-potentials (SPs) in the Chalk of England has shown that a consistent electrical potential gradient exists within a coastal groundwater borehole previously affected by seawater intrusion (SI) and that this gradient is absent in boreholes further inland. Furthermore, a small but characteristic reduction in this gradient was observed several days prior to SI occurring. We present results from a combined hydrodynamic and electrodynamic model, which matches the observed phenomena for the first time and sheds light on the source mechanisms for the spatial and temporal distribution of SP. The model predictions are highly sensitive to the relative contribution of electrochemical exclusion and diffusion potentials, the exclusion efficiency, in different rock strata. Geoelectric heterogeneity, largely due to marls and hardgrounds with a relatively high exclusion efficiency, was the key factor in controlling the magnitude of the modeled SP gradient ahead of the saline front and its evolution prior to breakthrough. The model results suggest that, where sufficient geoelectric heterogeneity exists, borehole SP may be used as an early warning mechanism for SI.

Salinity, both in soil and water, is a ubiquitous problem in coastal Bangladesh, particularly in the southwest. Salinity varies at the local scale (5–10 m), but the relative roles of land use and surface geology on salinity variation in near-surface (2) were explored in the southwest coastal region of Bangladesh. The sediments in the near-surface at both sites are fine and hydrometer analyses of cored samples indicate the dominance of silty clay (70%) along with very fine sand (5%), sandy clay (15%) and clay (10%) materials. Salinity variation in near-surface groundwater tends to follow land use rather than surface geology at the scale of our investigations. The study provides evidence of the influence of land use on the near-surface salinity variation and indicates the importance of land-use planning in salinity management in coastal areas.

Summary Previous research has found that irrigation water demand is relatively insensitive to water price, suggesting that increased pumping costs due to declining groundwater levels will have limited effects on agricultural water management practices. However, non-linear changes in well yields as aquifer saturated thickness is reduced may have large impacts on irrigated production that are currently neglected in projections of the long-term sustainability of groundwater-fed irrigation. We conduct empirical analysis of observation data and numerical simulations for case studies in Nebraska, USA, to compare the impacts of changes in well yield and groundwater depth on agricultural production. Our findings suggest that declining well pumping capacities reduce irrigated production areas and profits significantly, whereas increased pumping costs reduce profits but have minimal impacts on the intensity of groundwater-fed irrigation. We suggest, therefore, that management of the dynamic relationship between well yield and saturated thickness should be a core component of policies designed to enhance long-term food security and support adaptation to climate change.

Pit latrines are one of the most common forms of onsite sanitation facilities in many developing countries. These latrines are suitable as a means of isolating human waste, however, conditions within pits often lead to nitrification of the contained waste. In areas with a near-surface aquifer, the potential for nitrate pollution arising from pit latrines cannot be ignored. In this study, site visits were made to three densely populated, peri-urban areas near three West African cities (Dakar, Abidjan, Abomey-Calavi) to gather relevant information about the latrines in use and the soil and groundwater underneath the sites. Modelling was then conducted to demonstrate the potential for nitrate pollution of the groundwater from the latrines in such settings. The depth from the bottom of the pits to the water table was considered as 5, 10 or 30 m, to represent the range of aquifer depths at the study sites. Nitrate half-lives ranging from 500 to 1500 days were considered, and time scales from 6 months to several years were modelled. The results highlighted the high likelihood of nitrate pollution of groundwater reaching levels exceeding the World Health Organization guideline value for nitrate in drinking water of 50 mg/L after as short a period as two years for the aquifer situated 5 m below the pits, when considering moderate to long nitrate half-lives in the subsurface. Careful siting of latrines away from high water table areas, more frequent pit emptying, or switching to urine diversion toilets may be effective solutions to reduce nitrate passage from pit latrines into groundwater, although these solutions may not always be applicable, because of social, technical and economic constraints. The study highlights the need for more reliable data on the typical nitrate concentrations in pit latrines and the nitrate half-life in different subsurface conditions.

Estimating groundwater recharge rates is important for water resource management studies. Modeling approaches to forecast groundwater recharge typically require observed historic data to assist calibration. It is generally not possible to observe groundwater recharge rates directly. Therefore, in the past, much effort has been invested to record soil moisture content (SMC) data, which can be used in a water balance calculation to estimate groundwater recharge. In this context, SMC data is measured at different depths and then typically integrated with respect to depth to obtain a single set of aggregated SMC values, which are used as an estimate of the total water stored within a given soil profile. This article seeks to investigate the value of such aggregated SMC data for conditioning groundwater recharge models in this respect. A simple modeling approach is adopted, which utilizes an emulation of Richards’ equation in conjunction with a soil texture pedotransfer function. The only unknown parameters are soil texture. Monte Carlo simulation is performed for four different SMC monitoring sites. The model is used to estimate both aggregated SMC and groundwater recharge. The impact of conditioning the model to the aggregated SMC data is then explored in terms of its ability to reduce the uncertainty associated with recharge estimation. Whilst uncertainty in soil texture can lead to significant uncertainty in groundwater recharge estimation, it is found that aggregated SMC is virtually insensitive to soil texture.

In order to determine objectively the fractal behaviour of a time series, and to facilitate potential future attempts to assess model performance by incorporating fractal behaviour, a multi-order robust detrended fluctuation analysis (r-DFAn) procedure is developed herein. The r-DFAn procedure allows for robust and automated quantification of mono-fractal behaviour. The fractal behaviour is quantified with three parts: a global scaling exponent, crossovers, and local scaling exponents. The robustness of the r-DFAn procedure is established by the systematic use of robust regression, piecewise linear regression, Analysis of Covariance (ANCOVA) and Multiple Comparison Procedure to determine statistically significant scaling exponents and optimum crossover locations. The MATLAB code implementing the r-DFAn procedure has also been open sourced to enable reproducible results. r-DFAn will be illustrated on a synthetic signal after which is used to analyse high-resolution hydrologic data; although the r-DFAn procedure is not limited to hydrological or geophysical time series. The hydrological data are 4 year-long datasets (January 2012 to January 2016) of 1-min groundwater level, river stage, groundwater and river temperature, and 15-min precipitation and air temperature, at Wallingford, UK. The datasets are analysed in both time and fractal domains. The study area is a shallow riparian aquifer in hydraulic connection to River Thames, which traverses the site. The unusually high resolution datasets, along with the responsive nature of the aquifer, enable detailed examination of the various data and their interconnections in both time- and fractal-domains.

Integrated hydro-economic models have been widely applied to water management problems in regions of intensive groundwater-fed irrigation. However, policy interpretations may be limited as most existing models do not explicitly consider two important aspects of observed irrigation decision making, namely the limits on instantaneous irrigation rates imposed by well yield and the intraseasonal structure of irrigation planning. We develop a new modeling approach for determining irrigation demand that is based on observed farmer behavior and captures the impacts on production and water use of both well yield and climate. Through a case study of irrigated corn production in the Texas High Plains region of the United States we predict optimal irrigation strategies under variable levels of groundwater supply, and assess the limits of existing models for predicting land and groundwater use decisions by farmers. Our results show that irrigation behavior exhibits complex nonlinear responses to changes in groundwater availability. Declining well yields induce large reductions in the optimal size of irrigated area and irrigation use as constraints on instantaneous application rates limit the ability to maintain sufficient soil moisture to avoid negative impacts on crop yield. We demonstrate that this important behavioral response to limited groundwater availability is not captured by existing modeling approaches, which therefore may be unreliable predictors of irrigation demand, agricultural profitability, and resilience to climate change and aquifer depletion.

A still neglected consequence of climate change is increasing salinity levels in coastal areas. This is due to a combination of factors including sea-level rise and tropical cyclones transporting saltwater upstream. The salination of drinking water is a problem that could affect millions of people living in coastal areas globally. Excess dietary intake of salt is strongly associated with high blood pressure and gastric cancer. Saltwater intrusion can also have a detrimental effect on crops and agriculture. Measurements of salinity in coastal areas are urgently needed to estimate the extent of the problem and the correlation of saline waters with health outcomes. In situ measurements of salinity cannot capture the complicated spatial and temporal changes that take place in surface water bodies. In this article, we will discuss an effective alternative approach to mapping the salinity distribution of surface water using satellites.

BACKGROUND: Millions of coastal inhabitants in Southeast Asia have been experiencing increasing sodium concentrations in their drinking-water sources, likely partially due to climate change. High (dietary) sodium intake has convincingly been proven to increase risk of hypertension; it remains unknown, however, whether consumption of sodium in drinking water could have similar effects on health. OBJECTIVES: We present the results of a cohort study in which we assessed the effects of drinking-water sodium (DWS) on blood pressure (BP) in coastal populations in Bangladesh. METHODS: DWS, BP, and information on personal, lifestyle, and environmental factors were collected from 581 participants. We used generalized linear latent and mixed methods to model the effects of DWS on BP and assessed the associations between changes in DWS and BP when participants experienced changing sodium levels in water, switched from "conventional" ponds or tube wells to alternatives [managed aquifer recharge (MAR) and rainwater harvesting] that aimed to reduce sodium levels, or experienced a combination of these changes. RESULTS: DWS concentrations were highly associated with BP after adjustments for confounding factors. Furthermore, for each 100 mg/L reduction in sodium in drinking water, systolic/diastolic BP was lower on average by 0.95/0.57 mmHg, and odds of hypertension were lower by 14%. However, MAR did not consistently lower sodium levels. CONCLUSIONS: DWS is an important source of daily sodium intake in salinity-affected areas and is a risk factor for hypertension. Considering the likely increasing trend in coastal salinity, prompt action is required. Because MAR showed variable effects, alternative technologies for providing reliable, safe, low-sodium fresh water should be developed alongside improvements in MAR and evaluated in "real-life" salinity-affected settings. https://doi.org/10.1289/EHP659.

Background – Millions of people living in coastal areas in South-east Asia have been experiencing an increasing number of land-falling cyclones and subsequent inundations: these inundations cause e...

We report the first measured values of the streaming potential coupling coefficient in chalk samples saturated with natural groundwater, and preliminary field measurements of the spontaneous potential (SP), at both ambient and pumped conditions, at a test site in the Berkshire Chalk aquifer in the southern UK. The ultimate aim of the work is to use measurements of SP, in conjunction with borehole data, to characterize groundwater flow and aquifer properties. Laboratory measurements yield a value of the streaming potential coupling coefficient of −60 ± 4 mV MPa −1 and a corresponding zeta potential of −13 ± 1 mV. A negative zeta potential contrasts with previous published open-system measurements on artificial calcite, and may reflect the presence of organic material in the natural chalk samples or HCO 3 and SO 4 ions in the groundwater. Field measurements at ambient conditions show temporal variations in SP consistent with flow processes within the aquifer, but no coherent spatial variations. Measurements during water abstraction demonstrate that voltages at the ground surface and in monitoring boreholes become more positive during pressure drawdown and more negative during pressure build-up, consistent with the negative values of streaming potential coupling coefficient and zeta potential observed in the laboratory. Moreover, the magnitude of the change in voltage is similar to that estimated using the laboratory value of the coupling coefficient. Our results suggest that measurements of SP may make a valuable contribution to characterizing groundwater flow in the UK Chalk aquifer.

Groundwater in Chalk catchments is a major resource that also helps support internationally important habitats and ecosystems. Its dual porosity and dual permeability properties, coupled with large-scale structural features (such as hard rock layers and marls), produce a highly complex hydrogeological system. Recent impacts from groundwater flooding as well as vulnerability to drought have raised questions over the ability of traditional approaches to model these aquifers. Current work on near-surface hydrological processes has highlighted the importance of the soil and weathered zone for controlling recharge rates. In addition, karst-like features, sedimentary deposits and valley bottom processes govern stream–aquifer interaction and present a challenge in their representation in any modelling system. Methods that have, and are being, developed to incorporate these features, and their use in modelling Chalk catchments, are described. These are required in order to address major challenges, such as groundwater flooding and drought impacts, both of which could become more frequent and intense as a result of climate change.

Summary There is evidence that, under certain conditions, rapid preferential recharge via the fracture network can occur in Chalk aquifers. This has potentially important implications for contaminant migration through the Chalk unsaturated zone, CUZ, and for groundwater flooding in Chalk catchments. In the case of groundwater flooding, deficiencies in modelling aquifer response have been attributed to inadequate representation of flow processes in the CUZ. In this paper we consider two complementary approaches for assessing controls on preferential recharge to Chalk aquifers: an empirical approach and a physically-based modelling approach. We show that the main controls on preferential recharge to Chalk aquifers are the characteristics of rainfall events, in terms of duration and intensity, the physical properties of the near surface, and the antecedent soil moisture in the near surface. We demonstrate a number of deficiencies when past models of the CUZ are applied to the problem of simulating preferential recharge, notably that the assumption of instantaneous equilibrium between fractures and matrix is not valid, particularly during extreme recharge events. In order to simulate preferential recharge, fractures and matrix must be modelled as separate but interacting domains. This was achieved using a dual continua model. The model was computationally demanding, but was able to reproduce observed behaviour, including apparently hysteretic soil moisture characteristic relationships in the near surface, and rapid preferential recharge fluxes in response to high intensity rainfall events.

An integrated set of different measurements has been used to study the behavior of groundwater in an observation well in a fractured rock formation, the UK Chalk, under pumped and ambient conditions. Under pumped conditions, the response of the open borehole was relatively straightforward with flow mainly concentrated along four discrete flow horizons. Furthermore, excellent correspondence was observed between the three methods of borehole flow velocity measurement: impeller flowmeter, heat-pulse flowmeter and dilution testing. Under ambient conditions, the system appeared more complicated. Specifically, in the upper half of the borehole, the impeller flowmeter exhibited substantial downward flow and the heat-pulse flowmeter exhibited almost negligible upward flow, whilst dilution testing indicated significant dilution. It was concluded that this was due to cross-flow occurring over the upper 29 m. Analysis of drawdown data, recovery data and a Drost analysis of the ambient cross-flow data yielded aquifer transmissivity estimates of 2,049, 2,928 and > 4,388 m2/day respectively. The discrepancy between the drawdown and recovery estimates was attributed to non-linear head-losses associated with turbulence and inertial effects. The difference between the pumping test and Drost results was explained by the flow during the pumping test bypassing this aforementioned 29 m region of rock.

As with most fractured rock formations, Chalk is highly heterogeneous. Therefore, meaningful estimates of model parameters must be obtained at a scale comparable with the process of concern. These are frequently obtained by calibrating an appropriate model to observed concentration–time data from radially convergent tracer tests. Arguably, an appropriate model should consider radially convergent dispersion and Fickian matrix diffusion. Such a model requires the estimation of at least four parameters. A question arises as to whether this level of model complexity is supported by the information contained within the calibration data. A four-parameter model was developed for the analysis of two Chalk aquifer radially convergent tracer tests. The parameters included an advective travel time, t a , a characteristic fracture diffusion time, t cf , a characteristic matrix block diffusion time, t cb and a Peclet number, P. Because the tracer test duration was less than 500 h, t cb was impossible to identify. Further analysis showed that a large set of correlated values of P, t a and t cf would lead to equally good model fits. To resolve this ambiguity, more and better quality data are needed at the very start of the breakthrough curve, to constrain the mechanical dispersion parameter, P.

Summary Groundwater from unconfined chalk aquifers constitutes a major water resource in the UK. The unsaturated zone in such systems plays a crucial role in the hydrological cycle, determining the timing and magnitude of recharge, and the transport and fate of nutrients. However, despite more than three decades of study, our physical understanding of this system is incomplete. In this research, state of the art instrumentation provided high temporal resolution readings of soil moisture status, rainfall and actual evaporation from two sites in the Pang and Lambourn catchments (Berkshire, UK), for a continuous two year period (2004/5). A parsimonious, physically based model for the flow of water through the chalk unsaturated zone, including a novel representation of the soil and weathered chalk layers, was developed. The parameters were identified by inverse modelling using field measurements of water content and matric potential. The model was driven by rainfall and evaporation data, and simulated fluxes throughout the profile (including the discrete matrix and fracture components), down to the water table (but not the water table response). Results showed that the model was able to reproduce closely the observed changes in soil moisture status. Recharge was predominantly through the matrix, and the recharge response was strongly attenuated with depth. However, the activation of fast recharge pathways through fractures in the chalk unsaturated zone was highly sensitive to rainfall intensity. Relatively modest increases in rainfall led to dramatically different recharge patterns, with potentially important implications for groundwater flooding. The development and migration of zero flux planes with time and depth were simulated. The simulations also provided strong evidence that, for water table depths greater than 5 m, recharge fluxes persist throughout the entire year, even during drought conditions, with important implications for the calculation of specific yield from baseflow estimates and the representation of recharge in groundwater models.

An exact solution for transient Forchheimer flow to a well does not currently exist. However, this paper presents a set of approximate solutions, which can be used as a framework for verifying future numerical models that incorporate Forchheimer flow to wells. These include: a large time approximation derived using the method of matched asymptotic expansion; a Laplace transform approximation of the well-bore response, designed to work well when there is significant well-bore storage and flow is very turbulent; and a simple heuristic function for when flow is very turbulent and the well radius can be assumed infinitesimally small. All the approximations are compared to equivalent finite-difference solutions.

The Chalk unsaturated zone is crucial in controlling the delivery of nitrate to Chalk streams and groundwater abstraction wells. In this paper, results from a dual-permeability numerical model of the Chalk unsaturated zone are used to illustrate the relative roles of matrix and fracture flow. A major challenge arises in representing the Chalk unsaturated zone within catchment-scale models for nutrient management. These have generally been based on simple conceptual reservoirs or compartments to represent soils and groundwater. A more appropriate conceptualization has recently been developed and applied to the Lambourn catchment within a catchment-scale nutrient model (INCA-Chalk). Preliminary results from this work are discussed, which clearly illustrate the decadal time scales that need to be considered in the context of nutrient management and the Water Framework Directive.

The unsaturated zone exerts a major control on the delivery of nutrients to Chalk streams, yet flow and transport processes in this complex, dual-porosity medium have remained controversial. A major challenge arises in characterising these processes, both at the detailed mechanistic level and at an appropriate level for inclusion within catchment-scale models for nutrient management. The lowland catchment research (LOCAR) programme in the UK has provided a unique set of comprehensively instrumented groundwater-dominated catchments. Of these, the Pang and Lambourn, tributaries of the Thames near Reading, have been a particular focus for research into subsurface processes and surface water–groundwater interactions. Data from LOCAR and other sources, along with a new dual permeability numerical model of the Chalk, have been used to explore the relative roles of matrix and fracture flow within the unsaturated zone and resolve conflicting hypotheses of response. From the improved understanding gained through these explorations, a parsimonious conceptualisation of the general response of flow and transport within the Chalk unsaturated zone was formulated. This paper summarises the modelling and data findings of these explorations, and describes the integration of the new simplified unsaturated zone representation with a catchment-scale model of nutrients (INCA), resulting in a new model for catchment-scale flow and transport within Chalk systems: INCA-Chalk. This model is applied to the Lambourn, and results, including hindcast and forecast simulations, are presented. These clearly illustrate the decadal time-scales that need to be considered in the context of nutrient management and the EU Water Framework Directive.