Coastal Aquifer Research Articles

Uncertainty quantification in coastal aquifers using the multilevel Monte Carlo method

AbstractWe are solving a problem of salinisation of coastal aquifers. As a test case example, we consider the Henry saltwater intrusion problem. Since porosity, permeability and recharge are unknown or only known at a few points, we model them using random fields and random variables. The Henry problem describes a two‐phase flow and is non‐linear and time‐dependent. The solution to be found is the expectation of the salt mass fraction, which is uncertain and time‐dependent. To estimate this expectation, we use the well‐known multilevel Monte Carlo (MLMC) method. The MLMC method takes just a few samples on computationally expensive (fine) meshes and more samples on cheap (coarse) meshes. Then, by building a telescoping sum, the MLMC method estimates the expected value at a much lower computational cost than the classical Monte Carlo method. The deterministic solver used here is the well‐known parallel and scalable ug4 solver.

Modelling the Spatial Distribution and Future Trends of Seawater Intrusion due to Aquaculture Activities in Coastal Aquifers of Nizampatnam, Andhra Pradesh

Coastal aquifers suffer from saltwater intrusion due to several natural and anthropogenic activities and they pose a severe problem in many coastal regions of India. Anthropogenic activities like saltwater aquaculture activities and salt pans are one of the major causes for saltwater intrusion into freshwater aquifers in coastal regions. Along with these activities, over-exploitation of fresh groundwater for domestic and agriculture purpose also leads to infiltration of saltwater into shallow freshwater aquifers. To address this effect, it is necessary to understand the spatial distribution and extent of the saltwater intrusion in coastal areas. In the present study, SEAWAT model is used to delineate the saltwater intrusion effect and predict the extent of saltwater intrusion resulting from anthropogenic activities in the Nizampatnam coastal area, Andhra Pradesh. The simulation was carried out for the period 2016 to 2018 revealing a change in the extent of the saltwater intrusion up to 2 km into the inland areas from the coast as well as in the vicinity of fishponds. Based on the hydrological conditions observed during the study period and leveraging the calibrated model, saltwater intrusion prediction was carried out for the next two decades. The results show that anthropogenic activities will lead to saltwater intrusion into freshwater aquifers up to 6 km into inland areas from coast at a rate of 200 m/year. Therefore, the sustainable management plans are essential for protecting the freshwater aquifers from the saltwater intrusion in the Nizampatnam coastal area.

Spatio-temporal variability of seawater mixing in the coastal aquifers based on hydrogeochemical fingerprinting and statistical modeling

This study discusses monitoring and characterization of spatial and temporal variability to comprehend groundwater salinization based on hydrogeochemical fingerprinting and statistical modeling in the coastal belt of Digha-Shankarpur-Tajpur-Mandarmani, West Bengal, India. An integrated study of hydrogeochemical, bulk magnetic susceptibility, multivariate statistical, and geochemical modeling methods is implemented. The major cationic and anionic concentrations in groundwater are in the order Na+ > Ca2+ > Mg2+ > K+ and Cl− > HCO3− > SO4− > NO3− > F− respectively. The major water types are dominated by (Ca2+ – HCO3−) followed (Ca2+ - Mg2+ - Cl−), (Ca2+ - Na+ - HCO3−), (Na+ - HCO3−), and (Na+ - Cl). The results showed that the groundwater quality continuously declined steadily from pre-monsoon 2020 to pre-monsoon 2022. The deterioration of groundwater is due to an interplay of multiple factors, i.e., water-rock interaction, including ion-exchange, seawater mixing, and anthropogenic actions. Furthermore, it is also found that the regions showing higher seawater mixing index and oversaturated with carbonate minerals are also areas where groundwater is unsuitable for irrigation. The findings are beneficial in assisting local communities and legislators in designing appropriate management and mitigation techniques to arrest seawater intrusion in coastal regions.

Vertical leaching of paleo-saltwater in a coastal aquifer–aquitard system of the Pearl River Delta

Marine transgressive and regressive processes have promoted the formation of deltas and estuaries from the postglacial period, and the coastal aquifer-aquitard system can be identified based on the complicated sedimentary environments in the river delta. The groundwater system is usually composed of freshwater and paleo saltwater, and the vertical leaching of paleo saltwater from aquitard to fresh aquifer occurs due to density difference and gravity effect. In this study, radium quartet and radon were used as the tracers for groundwater due to their different half-lives. We constructed a radium reactive transport model considering the vertical leaching of paleo saltwater in a deltaic aquifer-aquitard system and then applied it in the Pearl River Delta (PRD). The activity ratio of 224Ra to 228Ra was used to characterize whether the groundwater system was controlled by α-recoil or sufficient 232Th under fully weathered condition, and results showed that the deltaic system was well weathered and the vertical leaching of paleo saltwater inhibited the α-recoil input and increased the decay of the sorbed parent nuclide in the aquifer. After incorporating the vertical leaching of saltwater in the model, the value of retardation factor in the shallow aquitard evidently increased, meaning that the leaching of salty porewater into the aquifer greatly affected the adsorption–desorption process of radium isotopes in the aquitard, and the value of retardation factor reduced to the lower limit in the basal aquifer, indicating that the desorption process was determined by the vertical leaching of saltwater here. After considering the effect of vertical leaching, the estimated water transit time in the aquifer was decreased and the groundwater discharge rate was increased, and the converse results were obtained in the aquitard. The proposed radium reactive transport model advances the understanding of the salinization process in coastal aquifers through vertical migration and facilitates the understanding of freshening process of paleo saltwater in the river delta with similar hydrogeological settings.

Geo-environmental factors controlling groundwater hydrogeochemistry, salinization and quality in El Qaa plain coastal aquifer, Southwest Sinai, Egypt

El Qaa plain is one of the promising areas for agriculture projects, owing to its soil and groundwater resources. So, 14 groundwater samples were collected, chemically analyzed, hydrochemical (ratios, models and plots) and GIS techniques were applied, to evaluate the factors governing this water chemistry and suitability for irrigation. The results indicated the expansion of the agricultural area 8 times in the last 2 decades. The water can be divided into two groups; (A) freshwater (TDS ≤ 1000) in the eastern part of the study area within the clastic weathered basement rocks part of the aquifer, and (B) saline water (TDS > 1000) in the western part within the carbonates and evaporites rich sediments and close to the Gulf of Suez. The deduced factors controlling groundwater chemistry were precipitation, evaporites dissolution, seawater intrusion and ion exchange process. The calculated sodium absorption ratio (SAR) was < 10 indicating the good quality of this water for irrigation without sodicity hazard. Water of Group (A) (electric conductivity < 2250µS/cm) can be used safely for irrigation without any impact of salinity on all types of soil and plants. Group (B) water (electric conductivity > 2250µS/cm) has salinity hazard and can be used for irrigation on well-permeable soils and salt-tolerant plants. The concentrations of the studied potentially toxic elements (PTEs) were within the permissible limit for irrigation. Expansion in agricultural activities has adverse impact on groundwater quality. So, eco-friendly fertilizers should be applied in this promising area with the application of hydrologic pumping advises to safe aquifer quality from the degradation.

Salinity origin in the coastal aquifer of the Southern Venice lowland

Groundwater salinization can be natural and anthropogenic in origin, although it often results from a combination of both, especially in low-lying coastal regions that are hydraulically controlled. This study proposes a method to assess the origin of salinity using environmental tracers in porewater, like Cl− and Br−, combined with depositional facies associations detected in sediment cores. Such integrated approach was tested in a target area south of the Venice Lagoon (Italy), where groundwater salinization is triggered by multiple mechanisms due to the complexity of the hydro-geomorphological environment. Batch tests were performed on sediment core samples from boreholes to quantify major anions and total inorganic N. Cl− and Br− porewater concentrations coupled with sedimentary facies association provided insights into the origin of groundwater salinity from a variety of sources, including past and present seawater intrusion, agricultural leaching, and evaporites. The strengths and limitations of the integrated approach are discussed to provide a pathway for improving water resource management and planning measures to prevent groundwater salinization in coastal areas.

Mitigating seawater intrusion in coastal aquifers: Novel approach with treated wastewater injection and groundwater circulation

Seawater intrusion (SWI) is a natural phenomenon that negatively impacts the potability of groundwater and is expected to worsen with rising sea levels due to climate change. Artificial recharge of freshwater is a commonly-used remediation method to mitigate SWI and improve freshwater supply security in affected coastal aquifers. However, limited freshwater availability can limit the effectiveness of this approach. This study proposes a novel mitigation measure, called Inj-GCW, which combines the injection of reclaimed water with the use of groundwater circulation wells (GCW) to enhance the effectiveness of artificial recharge in controlling SWI. GCW is a dual-screened well with isolated screens that extract and inject groundwater into the aquifer. The performance of the proposed measure was quantitatively evaluated using an illustrative simplified unconfined coastal aquifer. Based on the findings, the design parameters were estimated for a field-scale case study of the Nile Delta aquifer (NDA), a large Mediterranean coastal aquifer in Egypt. The study adopts a future scenario that considers Sea-level rise due to climate change and projected population growth by 2100. Results demonstrate that introducing of GCWs into the saltwater wedge, along with injection, effectively retreated the saltwater wedge, due to enhanced velocity, seaward fluxes, and dilution of contamination. Inj-GCW measure led to an 8.9% reduction in SWI and a 5.2% decrease in aquifer salinity compared to the expected intrusion in 2100. Furthermore, the Inj-GCW measure resulted in a 2.2% higher repulsion rate and 0.3% reduction in total salt mass compared to injection alone. The Inj-GCW measure presents a promising solution to SWI challenges in the NDA and other coastal aquifers facing similar issues. The formation of a brackish water bubble at the injection well screen of the GCW and the generated vertical groundwater circulation cells acted as a hydraulic barrier and contributed to the proposed method effectiveness.

Efficient representation of transient tidal overheight in a coastal groundwater flow model using a phase-averaged tidal boundary condition

Ocean tides cause water table overheight near the seaward boundary of coastal aquifers which can have a large impact on groundwater flows and saltwater intrusion. Despite this, regional-scale coastal groundwater flow models often neglect the effects of ocean tides or alternatively assume that the influences of ocean tides and sea level are constant over time. These simplifications are often required because simulation of the phase-resolved tidal signal at a coastal boundary is computationally demanding. Our objective was to derive a phase-averaged tidal boundary condition for application along gently sloping coastlines that enables simulation of real, complex tidal signals combined with sea-level changes due to meteorological effects. The boundary condition extends an existing analytical solution of tidal overheight by including an empirical correction function for conditions where the assumptions of the existing solution are violated. The correction function was developed by conducting a parametric study on an idealized two-dimensional cross-sectional coastal aquifer model and considering the effects of parameters including horizontal hydraulic conductivity, vertical anisotropy, specific yield, aquifer thickness, and beach slope. The performance of the new phase-averaged tidal boundary condition was assessed using a three-dimensional groundwater flow model of the island of Spiekeroog, Northwest Germany, which comprises complex coastal morphology with non-planar beach slopes. Simulations applying the new boundary condition were compared to simulations that adopted a phase-resolved tidal boundary condition and to observed groundwater levels. Accounting for time-varying changes in the tidal signal and sea level was found to be critical to simulate observed data and to adequately reproduce the transient tidal overheight simulated by the phase-resolved model. The performance of the new phase-averaged boundary condition in the regional-scale model varied depending on how the coastal morphology was represented in the boundary condition with local morphological features increasingly important for lower intertidal vertical infiltration capacity (low isotropic hydraulic conductivity or high vertical anisotropy) or specific yield. In conclusion, the new boundary condition presented overcomes current limitations in simulating transient tidal overheight in regional-scale groundwater models.

An adaptive multi-fidelity optimization framework based on co-Kriging surrogate models and stochastic sampling with application to coastal aquifer management

Surrogate modelling has been used successfully to alleviate the computational burden that results from high-fidelity numerical models of seawater intrusion in simulation-optimization routines. Nevertheless, little attention has been given to multi-fidelity modelling methods to address cases where only limited runs with computationally expensive seawater intrusion models are considered affordable imposing a limiting factor for single-fidelity surrogate-based optimization as well. In this work, a new adaptive multi-fidelity optimization framework is proposed based on co-Kriging surrogate models considering two model fidelities of seawater intrusion. The methodology is tailored to the needs of solving pumping optimization problems with computationally expensive constraint functions and utilizes only small high-fidelity training datasets. Results from both hypothetical and real-world optimization problems demonstrate the efficiency and practicality of the proposed framework to provide a steep improvement of the objective function while it outperforms a comprehensive single-fidelity surrogate-based optimization method. The method can also be used to locate optimal solutions in the region of the global optimum when larger high-fidelity training datasets are available.

Comparing Darcy’s Law and the Brinkman Equation for Numerical Simulations of Saltwater Intrusion

Saltwater intrusion into coastal aquifers presents a significant global challenge to fresh groundwater resources. Numerical modelling represents a valuable tool to study this phenomenon. Darcy’s Law is widely applied to groundwater studies and is extended into the Brinkman Equation to account for kinetic dissipations due to viscous shear. However, their comparative performance and accuracy in density-driven flows remain unclear. To determine the circumstances where the Brinkman Equation is required, numerical simulations with both models were implemented in hypothetical coastal aquifer scenarios. The results revealed that the largest discrepancies between the two models occur inside the dispersion zone during the break-through period, with concentration differences of up to 2.5%. The mixing of freshwater and saltwater induces rapid density and velocity variations. Brinkman’s viscous term moderates the rate of change and decreases the intrusion length by up to 6.1 m in a 180 m intrusion case. Furthermore, higher permeability and a lower recharge rate both strengthen the viscous effects in most sandy coastal aquifers. The Brinkman Equation excels at capturing intricate flow patterns with large variations. Therefore, it is necessary to be employed for studies on freshwater–saltwater interfaces and other similar conditions including groundwater–surface water interfaces, non-isothermal flows, and complex geological conditions.

Study of Salinity in Market Garden Production on the Togolese Coast

Agriculture accounts for over 40% of Togo's GDP, and employs almost 65% of its working population. Market gardening, one of its sub-sectors, remains an important source of income and food self-sufficiency for rural and peri-urban populations. Among the peri-urban production zones, the coastline, by virtue of its geographical position, is sensitive to the salinisation of water and soil, as well as to sea spray. Salinity has a negative impact on market garden production. This study was therefore carried out to investigate the salinity of water and soil used in market garden production on the Togolese coast. To do this, measurements were taken of the electrical conductivity (EC) and pH of irrigation water and soil sampled in the Ablogamé and Agodéké districts, two production areas on the coastal strip. In each of these areas, water and soil samples were taken at three sites, with three samples per site, giving a total of six sites for analysis. At Ablogamé, the values for irrigation water ranged from 1113.67±3.06 µS/cm, 1971.67±9.24 µS/cm, to 1926.33±6.51 µS/cm and for Agodéké, they ranged from 756.67 ± 1.53, 1016.00 ± 1.00 to 1178.33 ± 0.58 µS/cm. This shows that, according to CIRAD and WHO standards, the water analysed is classified as "salty", with a variability of EC within the same production zone, the exception being one of the Agodeke sites, which has soft water according to WHO standards. The FAO, for its part, classifies all the waters as having slight to moderate restrictions for irrigation. Furthermore, the EC values obtained are moderately and slightly high. In addition, those of the Agodéké sites, located at around 83 to 108 m from the sea, in the erosive zone, are generally lower than those of the Ablogamé sites located in the accumulation zone, at around 288 m to 406 m from the sea. The pH of the water is slightly alkaline. As for the soil, the EC values are below the norms at all the sites studied, thus demonstrating the non-saline nature of the market-garden production areas on the Togolese coastal strip. The soil pH is neutral to slightly alkaline. There is therefore a dynamism in the Togolese coastal aquifers that requires further studies to improve the performance of vegetable production on this coast.

Estimating the Maximum Safe Pumping Rate in Sloping Unconfined Coastal Aquifers

AbstractUnconfined coastal aquifers with a sloping aquifer bed are ubiquitous. However, most analytical solutions previously developed to estimate pumping‐induced seawater intrusion considered a horizontal aquifer bed. In this study, we first developed a steady‐state analytical solution to quantify the maximum safe pumping rate in sloping unconfined coastal aquifers with a fixed‐flux inland boundary condition, using the single potential approach. The analytical solution corrected by an empirical factor is validated against the numerical simulation, which reproduces the corresponding numerical results well. It is found that coastal aquifers with a higher aquifer bed toward inland (i.e., a positive sloping angle) yield a higher maximum safe pumping rate than that of coastal aquifers with a horizontal aquifer base, since the elevated freshwater head prevents seawater intrusion and enhances groundwater pumping. Specifically, for the aquifer bases with the sloping angle of 0.01 and −0.01, the maximum safe pumping rate is increased by 42.6% and decreased by 48.4%, respectively, in comparison to the case of a horizontal aquifer base, suggesting that neglecting the aquifer base slop can result in a significant error in estimating the maximum safe pumping rate. Moreover, the sensitivity analysis reveals that the maximum safe pumping rate increases with a lower hydraulic conductivity and a larger dispersivity. Our results provide valuable insights into the effect of the aquifer base slope on the maximum safe pumping in coastal aquifers, and analytical solutions developed can serve as a powerful tool for quick assessment of pumping‐induced seawater intrusion in sloping unconfined coastal aquifers.

Decreasing groundwater temperature relieves seawater intrusion in coastal aquifers

Fresh submarine groundwater discharge (FSGD), as an important natural and renewable water resource, occurs along much of the global coast but finally gets lost into the sea. Exploiting FSGD, however, could induce seawater intrusion in coastal aquifers. Here, we show that decreasing groundwater temperature (by extracting fresh groundwater from and recharging cooled and more viscous water back into coastal aquifers) can relieve seawater intrusion and lead to safe exploitation of FSGD, based on laboratory experiments and numerical simulations. The performance of groundwater cooling can be approximately quantified by an analytical solution which shows its dependence on the FSGD extraction ratio and temperature contrast. The seawater retreat ratio (ratio of the seawater retreat distance to the length of seawater intrusion prior to water extraction) can reach 47.4% with groundwater cooling and without fresh groundwater exploitation. Furthermore, without incurring further seawater intrusion, groundwater cooling can lead to 35.7% FSGD exploitation.

MODELING COASTAL WATER TABLE FLUCTUATIONS USING PFLOTRAN

Coastal aquifers are highly dynamic groundwater systems. Sea level rise will cause a rise in coastal groundwater tables resulting in increased risk of shallow or emergent groundwater (Befus et al., 2020). Marine water level fluctuations cause the beach groundwater table to oscillate over a relatively large range. Understanding these oscillations is crucial, as shallow (i.e., high) water tables may impact subsurface infrastructure, mobilize sediment, and increase liquefaction risks. Although the impacts of tides and wave setup on coastal water tables have been studied (e.g., Nielsen, 1990; Housego et al, 2021), the cumulative impacts of wave runup, partially saturated flow, complex beach topography, and dual tidal forcing for bay-backed regions have not been explored. This work numerically models beach water table fluctuations which are compared to in-situ swash and beach groundwater observations at Cardiff State Beach in Encinitas, CA.

Integrated clustering analysis for delineating seawater intrusion and heavy metals in Arabian Gulf Coastal groundwater of Saudi Arabia

The intrusion of seawater (SWI) into coastal aquifers is a major concern worldwide, affecting the quantity and quality of groundwater resources. The region of Saudi Arabia that lies along the eastern coast has been affected by SWI, making it crucial to accurately identify and monitor the affected areas. This investigation aimed to map the degree of seawater intrusion in a complex aquifer system in the study area using an integrated clustering analysis approach. The study collected 41 groundwater samples from wells penetrating multi-layered aquifers, and the samples were analyzed for physicochemical properties and major ions. Clustering analysis methods, including Hierarchical Clustering Analysis (double-clustering) (HCA-DC), K-mean (KMC), and fuzzy k-mean clustering (FKM), were employed to evaluate the spatial distribution and association of the groundwater properties. The results revealed that the analyzed GW samples were divided into four clusters with varying degrees of SWI. Clusters A, B, C, and D contained GW samples with very low (fsea of 1.9%), high (fsea of 14.9%), intermediate (fsea of 7.9%), and low (fsea of 5.2%) degrees of SWI, respectively. FKM clustering exhibited superior performance with a silhouette score of 0.83. Additionally, the study found a direct correlation between the degree of SWI and increased concentrations of boron, strontium, and iron, demonstrating SWI's impact on heavy metal levels. Notably, the boron concentration in cluster B, which endured high SWI, exceeded WHO guidelines. The study demonstrates the value of clustering analysis for accurately monitoring SWI and associated heavy metals. The findings can guide policies to mitigate SWI impacts and benefit groundwater-dependent communities. Further research can help develop effective strategies to mitigate SWI effects on groundwater quality and availability.

Wind-modulated groundwater discharge along a microtidal Arctic coastline

Groundwater discharge transports dissolved constituents to the ocean, affecting coastal carbon budgets and water quality. However, the magnitude and mechanisms of groundwater exchange along rapidly transitioning Arctic coastlines are largely unknown due to limited observations. Here, using first-of-its-kind coastal Arctic groundwater timeseries data, we evaluate the magnitude and drivers of groundwater discharge to Alaska’s Beaufort Sea coast. Darcy flux calculations reveal temporally variable groundwater fluxes, ranging from −6.5 cm d−1 (recharge) to 14.1 cm d−1 (discharge), with fluctuations in groundwater discharge or aquifer recharge over diurnal and multiday timescales during the open-water season. The average flux during the monitoring period of 4.9 cm d−1 is in line with previous estimates, but the maximum discharge exceeds previous estimates by over an order-of-magnitude. While the diurnal fluctuations are small due to the microtidal conditions, multiday variability is large and drives sustained periods of aquifer recharge and groundwater discharge. Results show that wind-driven lagoon water level changes are the dominant mechanism of fluctuations in land–sea hydraulic head gradients and, in turn, groundwater discharge. Given the microtidal conditions, low topographic relief, and limited rainfall along the Beaufort Sea coast, we identify wind as an important forcing mechanism of coastal groundwater discharge and aquifer recharge with implications for nearshore biogeochemistry. This study provides insights into groundwater flux dynamics along this coastline over time and highlights an oft overlooked discharge and circulation mechanism with implications towards refining solute export estimates to coastal Arctic waters.

Simulation of groundwater flow and seawater intrusion in response to climate change and human activities on the coastal aquifer of Gaza Strip, Palestine

This study has investigated the potential impact of human activities and climate change on the groundwater budget, water levels, and seawater intrusion in the coastal aquifer of the Gaza Strip (State of Palestine) over the next two decades. The impact of a proposed measure to use alternative freshwater provision from desalinated seawater on the future groundwater quantity and quality was also analyzed. Following extensive analysis of available observed data, a three-dimensional groundwater flow model, coupled with variable-density saltwater flow and transport components, was utilized for the investigations. Compared with the benchmark scenario (SC0), the climate change scenario (SC1) suggests that over the next two decades, an average annual aquifer recharge of 6.3 Mm3 can be expected, while the human activities scenario (SC2) indicates that the groundwater levels will decline at a rate of 0.09 m/year with expected urban area expansion. The combined human activities and climate change scenario (SC3) indicate severe groundwater storage depletion and seawater intrusion over the next two decades. The alternative freshwater provision scenario (SC4) indicates a strongly positive response in groundwater recovery (quantity and quality) over the next decades. The findings of this study emphasize strongly that the human activity impact, rather than climate change, is the driving force of groundwater depletion and that groundwater recovery interventions will be crucial in the future and should be implemented urgently.

A Hydrogeological Conceptual Model Refines the Behavior of a Mediterranean Coastal Aquifer System: A Key to Sustainable Groundwater Management (Grombalia, NE Tunisia)

The Mediterranean coastal aquifer system of the Grombalia basin (NE Tunisia) offers immense potential as a source of fresh water for agriculture, industry, and drinking water supply. Nonetheless, due to its intricate hydrogeological characteristics and the prevailing issue of groundwater salinity, comprehending its groundwater system behavior becomes crucial for the effective and sustainable management of this aquifer system. Based on the hydrogeological characterization of the Grombalia basin, a novel 3D hydrogeological conceptual model was developed to enhance the understanding of its complex aquifer system. The integration of insights from geological, hydrogeological, hydrodynamic, and hydrochemical components facilitated the construction of the hydrogeological conceptual model. Although the model’s validity faced initial uncertainties due to spatial interpolation of lithological sequences, this study’s thorough and encompassing hydrogeological investigation overcame these limitations. As a result, a more informed comprehension of the aquifer system complexities was achieved. This study reveals that the basin is underlain by an extensive, cohesive Mio–Plio–Quaternary aquifer system. The model demonstrates vertical and lateral hydrogeological continuity between the Quaternary and underlying Mio–Pliocene deposits, enabling groundwater flow and exchange between these layers. Over-abstraction of the Mio–Plio–Quaternary aquifer system has led to a significant drop in piezometric levels and raised the risk of seawater intrusion. These findings emphasize the critical necessity of taking into account the interconnections among hydrogeological units to ensure sustainable groundwater management. The developed conceptual model offers a key tool for understanding the hydrodynamic functioning of the Grombalia aquifer system with a view toward guiding future groundwater management strategies. The application of this approach in the Grombalia basin suggests its potential applicability to other regional aquifers facing comparable challenges.

Seawater–Groundwater Interaction Governs Trace Metal Zonation in a Coastal Sandy Aquifer

AbstractTrace metals in the groundwater of coastal sandy aquifers significantly influence the coastal ecosystem, yet the spatiotemporal controls of these metals remain unclear. In this paper, we comprehensively revealed the distribution patterns, key controlling factors, potential ecological risks, and fluxes to the ocean of groundwater trace metals (As, Ba, Cr, Cd, Fe, Mn, Pb, and Zn) in a coastal deep sandy aquifer. The results showed a clear zonation of trace metals in the groundwater in relation to the mixing extent between seawater and terrestrial freshwater. The freshwater zone exhibited a relatively low concentration of trace metals, whereas the freshwater‐seawater transition zone showed a substantial quantity of dissolved Fe, Mn, As, and Ba. Seawater‐groundwater interactions significantly affected the Fe, Mn, As, and Ba concentrations through redox potential and pH gradients. The tide‐driven saline water zone was vulnerable to oceanic environments and anthropogenic activities, resulting in the enrichment of trace metals such as Zn and Cd. Driven by recirculated submarine groundwater discharge (SGD), the concentrations of trace metals in the density‐driven saline circulation zone were found to be higher than those in the surrounding areas. The ecological risk index suggested that the freshwater‐seawater transition zone posted the highest ecological risks. Trace metal fluxes (i.e., Fe, Mn, and As) via SGD significantly contributed to the total input into the sea, which may have potential impacts on the coastal environments. Our study highlighted the importance of seawater‐groundwater interactions on trace element cycling in coastal sandy aquifers.

Temporal and spatial distribution patterns of upper saline plumes and seawater-groundwater exchange under tidal effect

Under tidal influence, an active exchange process occurs between seawater and groundwater in the coastal aquifer system. However, in-situ studies and research on the temporal and spatial distribution patterns of upper saline plumes (USPs) and seawater-groundwater exchange process are lacking. In this study, we selected the Shilaoren beach in Qingdao as the study area, and used pseudo-3D electrical resistivity tomography (ERT) and hydrogeological measurement methods to monitor the spatial distribution of groundwater salinity and the seawater and groundwater exchange process. We found that during the ebbing tide, USPs were divided into several smaller USPs driven by groundwater discharge. Subsequently, during the flooding tide, high-salinity seawater recharge was observed, and several smaller USPs gradually aggregated into a complete USP. In coastal aquifers, the exchange of seawater and groundwater not only occurs through the beach surface and seepage face, but also includes high-permeability underground channels. Seawater and groundwater exchange occurs more rapidly through these underground channels randomly distributed at depth, which also leads to spatial differences in the seawater and groundwater exchange process. Additionally, affected by underground channels, the groundwater discharge flux shows spatial heterogeneity. In areas with underground channels, the groundwater discharge flux significantly increases. Therefore, the non-seepage face area has a similar groundwater discharge flux to the main groundwater discharge area, namely, the seepage face. In the intertidal zone of the study area (70 × 60 m), the total groundwater discharge flux through the seepage face is 3407.6 m3, and no significant difference is observed in the groundwater discharge flux at various parts of the seepage face. In the non-seepage face area, the groundwater discharge is highly spatially heterogeneous; the total groundwater discharge flux is 3252.0 m3, and approximately 60% of this occurs through densely developed underground channels.