Abstract:Knowing the distribution between fresh and saline groundwater is imperative for sustainable and integrated management of water resources in coastal areas. The airborne electromagnetic (AEM) method is increasingly used for hydrogeological mapping over large areas via bulk electrical resistivity. However, accurately and reliably mapping the fresh/saltwater interface (FSI) requires accurate knowledge about the transition zone. The objective is to quantify the uncertainty in using AEM data to inform on the depth of the FSI. The study mimics a dual-moment time-domain SkyTEM sounding recorded in the Belgian coastal plain based on borehole data. It quantifies uncertainty using a differential evolution adaptive Metropolis algorithm to sample the posterior distribution. The results indicate the importance of reliable altitude, pitch and roll logging. Gathering prior knowledge about the transition zone, for example, through borehole logs, significantly improves the estimation of the FSI. The Resolve frequency-domain system, especially in context with very shallow to shallow FSIs, is more suitable for salinity mapping than the time-domain SkyTEM used in the field survey. The depth of the FSI may be defined via various threshold values. The uncertainty of three different thresholds is studied. The FSI based on the middle of the transition zone is the most reliable, while the FSI based on the 1500 mg/L total dissolved solids threshold is the least robust.

What problem does this paper attempt to address?

### Problems the paper attempts to solve This paper aims to quantify the uncertainty of the fresh water/saline water interface (FSI) depth using airborne electromagnetic (AEM) data. Specifically, the main objectives of the study include: 1. **Quantify uncertainty**: Quantify the uncertainty of determining the FSI depth using AEM data by sampling the posterior distribution with the Differential Evolution Adaptive Metropolis algorithm (DE - MCMC). 2. **Evaluate the influence of different factors**: - **Flight altitude and its uncertainty**: Study the influence of flight altitude and its measurement error on the FSI depth estimation. - **Tilt of the transmitter - receiver frame**: Analyze the influence of the pitch and roll angles of the transmitter - receiver frame on the results. - **The role of prior knowledge**: Explore the improvement of FSI depth estimation through prior knowledge such as borehole logs. 3. **Compare different systems**: Compare the applicability of the frequency - domain electromagnetic system (FDEM) and the time - domain SkyTEM system in shallow FSI imaging. 4. **Define different thresholds of FSI**: Study the reliability and robustness of defining FSI based on the mid - point of the transition zone, the 1500 mg/L total dissolved solids (TDS) threshold, and the optimistic threshold. ### Research background In coastal areas, understanding the distribution of fresh water and salt water is crucial for the sustainable management and protection of water resources. The airborne electromagnetic method (AEM) is widely used in hydrogeological mapping because it can quickly obtain underground conductivity information over a large area. However, accurate and reliable mapping of the fresh water/saline water interface (FSI) requires precise knowledge of the transition zone. By simulating SkyTEM detection data in the Belgian coastal plain, this paper explores how to quantify the uncertainty of determining the FSI depth using AEM data and puts forward some improvement suggestions. ### Main methods 1. **Probabilistic inversion**: Adopt the Bayesian method, combine the prior distribution and the likelihood function, and generate a set of models that conform to the observed data. Use the DE - MCMC algorithm to sample the posterior distribution. 2. **Parameterization**: Based on borehole logs, parameterize the underground structure as a fresh - water lens mainly composed of sandy sediments, a saline - water lens, and a clay half - space, as well as the transition zone between them. The transition zone is described by a logistic function. 3. **FSI threshold definition**: Three different FSI definition methods are studied: based on the mid - point of the transition zone, the 1500 mg/L TDS threshold, and the optimistic threshold (0.5 S/m water conductivity or 0.125 S/m volume conductivity in sandy sediments). ### Results 1. **Non - uniqueness problem**: The study found that even with perfect knowledge of altitude, the non - uniqueness of the inverse problem still significantly affects the estimation of the FSI depth. In particular, there is a bimodal distribution in the conductivity estimation of the fresh - water lens. 2. **Flight altitude uncertainty**: The uncertainty of flight altitude has a significant impact on the FSI depth estimation. A higher flight altitude will lead to greater uncertainty and deviation. 3. **Performance of different thresholds**: The FSI definition based on the mid - point of the transition zone is the most reliable, while the FSI definition based on the 1500 mg/L TDS threshold is the least robust. 4. **Importance of prior knowledge**: Prior knowledge such as borehole logs can significantly improve the accuracy of FSI depth estimation. ### Conclusions This study provides an important reference for future water resource management in coastal areas by quantifying the uncertainty of AEM data in determining the FSI depth. The study shows that the uncertainty of flight altitude, the use of prior knowledge, and the selection of FSI thresholds have a significant impact on the results. It is recommended to consider these factors in future investigations to improve the accuracy and reliability of FSI depth estimation.

Quantitative imaging of the fresh/saltwater interface with airborne electromagnetics: examining different sources of uncertainty

Using an airborne electromagnetic method to map saltwater intrusion in the northern Salinas Valley, California

Water Flux Characterization Through Hydraulic Head and Temperature Data Assimilation: Numerical Modeling and Sandbox Experiments

Effect of Groundwater Extraction and Artificial Recharge on the Geophysical Footprints of Fresh Submarine Groundwater Discharge in the Western Belgian Coastal Area

A Likelihood-Based Triangulation Method for Uncertainties in Through-Water Depth Mapping

3D frequency-domain modeling of marine controlled source electromagnetic responses with topography using finite volume method

Uncertainty Quantification in Geostatistical Modelling of Saltwater Intrusion at a Coastal Aquifer System

Characterizing Offshore Freshened Groundwater Salinity Patterns Using Trans‐Dimensional Bayesian Inversion of Controlled Source Electromagnetic Data: A Case Study From the Canterbury Bight, New Zealand

Application of the Electromagnetic Method to the Spatial Distribution of Subsurface Saline and Fresh Water in the Coastal Mudflat Area of Jiangsu Province

Seafloor massive sulphide exploration using deep-towed controlled source electromagnetics: Navigational uncertainties

Mapping Underwater Topography for Coastal Estuary Ecosystem Monitoring from Quickbird High-Resolution Multispectral Imagery

Mapping the Hydrogeological Structure of a Small Danish Island Using Transient Electromagnetic Methods

Quantitative Estimation of Soil Salinity Using UAV-Borne Hyperspectral and Satellite Multispectral Images

Quantifying uncertainty of salt body shapes recovered from gravity data using trans-dimensional Markov chain Monte Carlo sampling

Unified Physical Mechanism of Frequency-Domain Controlled-Source Electromagnetic Exploration on Land and in Ocean

Mapping Spatial Variability of Soil Salinity in a Coastal Paddy Field Based on Electromagnetic Sensors

Modeling dissolved organic carbon in inland waters using an unmanned aerial vehicles-borne hyperspectral camera

Successive bootstrapping deep learning approach and airborne EM-borehole data fusion to understand salt water in the Mississippi River Valley Alluvial Aquifer

Integrated Airborne Geophysical Data for Mapping Saltwater in the Mississippi River Valley alluvial Aquifer, Northeast Louisiana

Improved imaging of a karst aquifer using focused source electromagnetic and differentially normalized method: a qualitative analysis

Footprint-driven stochastic approximation inversion for time-domain airborne electromagnetic data