The Salt Water Intrusion Meetings, or SWIMs, are a series of meetings that focus on seawater intrusion in coastal aquifers and other salinisation processes. 2018 marks the 50th year of the SWIM and the 25th biennial meeting. The SWIM proceedings record half a century of research progress on site characterisation, geophysical and geochemical techniques, variable-density flow, modelling, and water management. The SWIM is positioning itself to remain a viable platform for discussing the coastal aquifer management challenges of the next 50 years.

Since the end of the Cambodian Civil War in 1998, the population of the Oddar Meanchey province has drastically increased despite the lack of adequate infrastructure, including basic amenities such as drinking-water supply. To improve the access to drinking water, governmental and aid agencies have focussed on drilling shallow boreholes. The use of groundwater for irrigation is also a growing concern to cope with the occasional late arrival of the rainy season or to produce food during the dry season. Since the groundwater resource in the province has not been documented, a 4-year study was undertaken (2011–2014), aiming to estimate the capability of groundwater to supply domestic needs and supplementary irrigation for rice production. Aquifer properties were estimated by combined use of hydrogeological techniques with the geophysical magnetic resonance sounding method. Groundwater storage and recharge were estimated based on new developments in the application of the geophysical method for quantifying specific yield. The median groundwater storage of the targeted sandstone aquifer is 173 mm, the recharge is diffuse and annually ranges from 10 to 70 mm, and the transmissivity is low to medium. Simulations of pumping indicate that the aquifer can easily supply 100 L of drinking water per capita daily, even considering the estimated population in 2030. However, the shallow aquifer can generally not deliver enough water to irrigate paddy fields of several hectares during a 2-month delay in the onset of the monsoon.

Coastal groundwater flow is driven by an interplay between terrestrial and marine forcings. One of the distinguishing features in these settings is the formation of a freshwater lens due to the density difference between fresh and saline groundwater. The present study uses data collected on Sable Island, Canada—a remote sand island in the northwest Atlantic Ocean—to highlight the potential of exploiting freshwater lens geometry for calibration of numerical groundwater flow models in coastal settings. Three numerical three-dimensional variable-density groundwater flow models were constructed for different segments of the island, with only one model calibrated using the freshwater–saltwater interface derived from an electromagnetic geophysical survey. The other two (uncalibrated) models with the same parameterisation as the calibrated model successfully reproduced the interpreted interface depth and location of freshwater ponds at different parts of the island. The successful numerical model calibration, based solely on the geophysically derived interface depth, is enabled by the interface acting as an amplified version of the water table, which reduces the relative impact of the interpreted depth uncertainty. Furthermore, the freshwater–saltwater interface is far more inertial than the water table, making it less sensitive to short-term forcings. Such “noise-filtering” behaviour enables the use of the freshwater–saltwater interface for calibration even in dynamic settings where selection of representative groundwater heads is challenging. The completed models provide insights into island freshwater lens behaviour and highlight the role of periodic beach inundation and wave overheight in driving short-term water-table variability, despite their limited impact on the interface depth.

A hydrogeophysical field experiment was conducted on a karst hydrosystem in the south of France to investigate groundwater transfer and storage variability at a scale of a few hundred meters. A 200-m-long N/S tunnel going through limestone provided the unique opportunity to set up measurements with original configurations inside the unsaturated zone. Three geophysical methods were used: gravimetry, electrical, and seismic. Two-dimensional electrical resistivity and seismic velocity images were retrieved by surrounding the medium with electrodes and geophones, both at the surface and inside the tunnel to improve the sensitivity in depth. This gave information about the weathering state but also about the limestone content and associated porosity characteristics, as the methods are sensitive to distinct properties with different resolution patterns. A time-lapse gravity surface-to-tunnel profile supplied information on the seasonal water mass changes and its variations along the tunnel. Besides, tracers were injected on each side of the profile from the surface and the restitution was sampled in the tunnel drip flows. A contrasting hydrological behavior was evidenced on each side of the tunnel from temporal gravity measurements and tracing tests. The analysis of the whole dataset allowed for better interpretation of the imaged structures, with different hydrological functioning. This study demonstrates the variability of the karst behavior at the scale of a few hundred meters and the benefits of a multi-method approach coupling hydrological and geophysical measurements. This kind of experiment provides fundamental understanding of systems that cannot be directly observed.

The existence of submarine fresh groundwater has been recorded at continental shelves worldwide. The dynamic preservation and lifetime of fresh groundwater in the offshore environment remains an open hydrogeological problem. The mechanisms and time scales of fresh groundwater preservation are examined using numerical simulations based on a geologically representative model of the New Jersey shelf, USA. Utilizing two-dimensional depth-migrated seismic and well data, a detailed hydrogeological model is built, with a vertical resolution of 10 m. The model captures the highly heterogeneous shelf environment and accounts for porosity compaction trends derived from core data. The results show transient coupled simulations of groundwater flow and heat and salt transport from the late Pleistocene until present day and projected 18,000 years into the future. They reveal freshwater preservation patterns and yield simulated borehole salinity profiles broadly consistent with field observations. The simulations show that freshwater intervals of a thickness of 200–300 m and lateral extent of tens of kilometers may have been preserved from the Last Glacial Maximum until today. It was found that approximately 30–45% of the initial freshwater volume remains preserved after 12,000 years, depending on the recharge boundary condition. The preserved volume ranges between 15 and 30% after 30,000 years. These results improve the understanding of submarine preservation of fresh groundwater through an interdisciplinary approach which integrates seismic imaging, hydrogeological modeling and high-performance numerical simulation.

In coral islands, groundwater is a crucial freshwater resource for terrestrial life, including human water supply. Response of the freshwater lens to expected climate changes and subsequent vegetation alterations is quantified for Grande Glorieuse, a low-lying coral island in the Western Indian Ocean. Distributed models of recharge, evapotranspiration and saltwater phytotoxicity are integrated into a variable-density groundwater model to simulate the evolution of groundwater salinity. Model results are assessed against field observations including groundwater and geophysical measurements. Simulations show the major control currently exerted by the vegetation with regards to the lens morphology and the high sensitivity of the lens to climate alterations, impacting both quantity and salinity. Long-term changes in mean sea level and climatic conditions (rainfall and evapotranspiration) are predicted to be responsible for an average increase in salinity approaching 140 % (+8 kg m⁻³) when combined. In low-lying areas with high vegetation density, these changes top +300 % (+10 kg m⁻³). However, due to salinity increase and its phytotoxicity, it is shown that a corollary drop in vegetation activity can buffer the alteration of fresh groundwater. This illustrates the importance of accounting for vegetation dynamics to study groundwater in coral islands.

Relict rock glaciers are complex hydrogeological systems that might act as relevant groundwater storages; therefore, the discharge behavior of these alpine landforms needs to be better understood. Hydrogeological and geophysical investigations at a relict rock glacier in the Niedere Tauern Range (Austria) reveal a slow and fast flow component that appear to be related to the heterogeneous structure of the aquifer. A numerical groundwater flow model was used to indicate the influence of important internal structures such as layering, preferential flow paths and aquifer-base topography. Discharge dynamics can be reproduced reasonably by both introducing layers of strongly different hydraulic conductivities or by a network of highly conductive channels within a low-conductivity zone. Moreover, the topography of the aquifer base influences the discharge dynamics, which can be observed particularly in simply structured aquifers. Hydraulic conductivity differences of three orders of magnitude are required to account for the observed discharge behavior: a highly conductive layer and/or channel network controlling the fast and flashy spring responses to recharge events, as opposed to less conductive sediment accumulations sustaining the long-term base flow. The results show that the hydraulic behavior of this relict rock glacier and likely that of others can be adequately represented by two aquifer components. However, the attempt to characterize the two components by inverse modeling results in ambiguity of internal structures when solely discharge data are available.

An innovative mode of groundwater recharge to a buried esker aquifer is considered. The current conceptual model affords a natural safeguard to underlying aquifers from the overlying muds. A hypothesis of groundwater recharge to a buried esker aquifer via preferential pathways across its overlying muds is tested here by heuristic numerical one-dimensional and two-dimensional modeling simulations. The hypothesis has been tested against two other conventionally accepted scenarios involving: (1) distal esker outcrop areas and (2) remote shallow-bedrock recharge areas. The main evidence comes from documented recharge pressure pulses in the overlying mud aquitard and in the underlying esker hydraulic-head time series for the Vars-Winchester esker aquifer in Eastern Ontario, Canada. These perturbations to the potentiometric surface are believed to be the aquifer response to recharge events. The migration rate of these pressure pulses is directly related to the hydraulic diffusivity of the formation. The measured response time and response amplitude between singular radar precipitation events and well hydrographs constituted the heuristic model calibration targets. The main evidence also includes mud-layering deformation (water escape features) which was observed in seismic surveys of the over-esker muds. These disturbed stratigraphic elements provide a realistic mechanism for migrating water to transit through the muds. The effective hydraulic conductivities of these preferential pathways in the muds were estimated to be between 2 × 10⁻⁶ and 7 × 10⁻⁶ m/s. The implications of these findings relate to the alleged natural safeguard of these overlying muds.

Three-dimensional geological and groundwater flow models of a submarine groundwater discharge (SGD) site at Hanko (Finland), in the northern Baltic Sea, have been developed to provide a geological framework and a tool for the estimation of SGD rates into the coastal sea. The dataset used consists of gravimetric, ground-penetrating radar and shallow seismic surveys, drill logs, groundwater level monitoring data, field observations, and a LiDAR digital elevation model. The geological model is constrained by the local geometry of late Pleistocene and Holocene deposits, including till, glacial coarse-grained and fine-grained sediments, post-glacial mud, and coarse-grained littoral and aeolian deposits. The coarse-grained aquifer sediments form a shallow shore platform that extends approximately 100–250 m offshore, where the unit slopes steeply seawards and becomes covered by glacial and post-glacial muds. Groundwater flow preferentially takes place in channel-fill outwash coarse-grained sediments and sand and gravel interbeds that provide conduits of higher hydraulic conductivity, and have led to the formation of pockmarks on the seafloor in areas of thin or absent mud cover. The groundwater flow model estimated the average SGD rate per square meter of the seafloor at 0.22 cm day⁻¹ in autumn 2017. The average SGD rate increased to 0.28 cm day⁻¹ as a response to an approximately 30% increase in recharge in spring 2020. Sensitivity analysis shows that recharge has a larger influence on SGD rate compared with aquifer hydraulic conductivity and the seafloor conductance. An increase in recharge in this region will cause more SGD into the Baltic Sea.

Accurate characterization of aquifers remains challenging for large-scale systems because of the spatial heterogeneity of hydraulic properties and temporal variability of hydrologic inputs. This study highlights the importance of integrating geological, hydrogeological and geophysical approaches to characterize an aquifer and evaluate groundwater productivity. Data from geological maps, drill logs, a pumping test, vertical electrical soundings (VES) and different field hydrogeological studies were combined and applied to a heavily extracted aquifer system—Lake Haramaya watershed, Ethiopia. From the geological characterization, the aquifer was found to be a single heterogeneous and anisotropic unconfined unit. Significant differences were found between the three-dimensional geological models of the aquifer developed from the drill logs and VES data; the VES data were likely affected by moisture content. The pumping-test and VES data were combined to estimate transmissivity (T; 126.5 ± 25.8 m²/day) and hydraulic conductivity (K; 4.1 ± 1.0 m/day). This combined use allowed for a reduction in uncertainty (40.1% for T and 33.3% for K) compared with values estimated from the VES data alone. The combined approach also allowed for much greater spatial coverage and a higher resolution characterization of the aquifer. The available volume of groundwater resource in the system was estimated at ~0.62 ± 0.09 km³. The groundwater extraction rate was ~30,120 m³/day, approximately double the estimated sustainable yield of the aquifer (15,720 m³/day). This showed that the current exploitation rate could exhaust groundwater resources in 27–32 years and should be reduced by 50% to ensure sustainability of the groundwater resource.