3D Characterization of a Coastal Freshwater Aquifer in SE Malta (Mediterranean Sea) by Time-Domain Electromagnetics

Electromagnetic (EM) geophysical methods are well equipped to distinguish electrical resistivity contrasts between freshwater-saturated and seawater-saturated formations. Beneath the semi-arid, rapidly urbanizing island of Malta, offshore groundwater is an important potential resource but it is not known whether the regional mean sea-level aquifer (MSLA) extends offshore. To address this uncertainty, land-based alongshore and across-shore time-domain electromagnetic (TDEM) responses were acquired with the G-TEM instrument (Geonics Ltd., Mississauga, ON, Canada) and used to map the onshore structure of the aquifer. 1-D inversion results suggest that the onshore freshwater aquifer resides at 4–24 m depth, underlain by seawater-saturated formations. The freshwater aquifer thickens with distance from the coastline. We present 2D and 3D electromagnetic forward modeling based on finite-element (FE) analysis to further constrain the subsurface geometry of the onshore freshwater body. We interpret the high resistivity zones that as brackish water-saturated bodies are associated with the mean sea-level aquifer. Generally, time-domain electromagnetic (TDEM) results provide valuable onshore hydrogeological information, which can be augmented with marine and coastal transition-zone measurements to assess potential hydraulic continuity of terrestrial aquifers extending offshore

3D characterisation and quantification of an offshore freshened groundwater system in the Canterbury Bight

Although offshore freshened groundwater (OFG) systems have been documented in numerous continental margins worldwide, their geometry, controls and emplacement dynamics remain poorly constrained. Here we integrate controlled-source electromagnetic, seismic reflection and borehole data with hydrological modelling to quantitatively characterise a previously unknown OFG system near Canterbury, New Zealand. The OFG system consists of one main, and two smaller, low salinity groundwater bodies. The main body extends up to 60 km from the coast and a seawater depth of 110 m. We attribute along-shelf variability in salinity to permeability heterogeneity due to permeable conduits and normal faults, and to recharge from rivers during sea level lowstands. A meteoric origin of the OFG and active groundwater migration from onshore are inferred. However, modelling results suggest that the majority of the OFG was emplaced via topographically-driven flow during sea level lowstands in the last 300 ka. Global volumetric estimates of OFG will be significantly revised if active margins, with steep coastal topographies like the Canterbury margin, are considered

Controlled source electromagnetic measurements offshore the Maltese Islands : implications for offshore freshened groundwater in a carbonate shelf

Carbonate lithologies host considerable quantities of the Earth’s freshwater resources and partially supply a significant amount of the global population with drinkable water. Although they comprise substantial amounts of the coastlines, it is not known if these carbonate lithologies can sustain freshened groundwater offshore, and if this can help meet future water demands in coastal regions. To date, predicting volumes of freshened groundwater within marine carbonates has been challenging. Here, we integrate controlled source electromagnetic profiles with seismic reflection and core log data to derive a lithological model for the eastern carbonate margin of the Maltese Islands, one of the most water-starved countries in the world. Electrical resistivity models are used to guide lithological inference where seismic data provide limited information due to the superimposed seafloor multiple. We show that resistivity values within the Upper Coralline and Globigerina Limestone formations exceed the measured resistivity of seawater-saturated core log samples by at least a factor of four. This could be indicative of offshore freshened groundwater that occupies the pore space of the low permeability limestone along the eastern Maltese shelf. To validate this observation without further ground-truthing data, we use extensive forward modelling to show that a similar resistivity footprint can be achieved by localized interbedded low-porosity or highly cemented units. However, the spatial extent of such units across the entire eastern Maltese margin is geologically improbable. This points to the occurrence of offshore freshened groundwater that was likely emplaced during the last sea-level lowstand.peer-reviewe

Electrical resistivity anomalies offshore a carbonate coastline : evidence for freshened groundwater?

Carbonate lithologies host considerable quantities of the Earth's freshwater resources and partially supply a quarter of the global population with drinkable water. In addition, carbonates constitute substantial amounts of the global coastlines, yet it is not known if and how they can sustain freshened groundwater offshore. Here, we use controlled-source electromagnetic, seismic reflection, and core sample data to derive a lithological model for the eastern margin of the Maltese Islands and identify four distinct resistivity anomalies within the Upper Coralline and Globigerina Limestone formations. The anomalies hosted in the former are likely associated with low porosities, whereas the anomaly within the latter is indicative of pore fluid freshening. Hydrogeological modeling suggests that freshened pore fluids, emplaced during sea-level lowstands and preserved in low permeability units, are potentially still found within carbonate shelves. However, resource potential is low due to its relict nature and low permeability host environment.peer-reviewe

2D Inversion of Marine Time Domain CSEM Data: Developments and Application to an Offshore Aquifer Investigation in the Canterbury Basin, New Zealand

Time-Domain Controlled Source Electromagnetic (TD-CSEM) methods are commonly applied to detect and delineate sub-surface resistivity anomalies on land and in the marine environment. The most general form of the application utilizes a horizontal electric dipole (HED) current excitation with an alternating current step function ranging from periods of a couple 100 ms to several seconds. TD-CSEM measurements are typically conducted along a transect in a towed inline configuration. Accordingly, the data is mainly sensitive to resistivity variations located between the source and the receiver. Therefore, 2D inverse modelling is often considered sufficient to adequately represent the measured data in regions of simple geometry/bathymetry. Although recording of TD-CSEM data is quite sophisticated in terms of instrumentation, the interpretation in academia is still based on simplified assumptions of 1D resistivity models that may fail to represent the true resistivity distribution of the sub-surface. Here, we present a development of the freely available frequency-domain code MARE2DEM, which now allows the inversion of TD-CSEM data. The theoretical code developments are tested in terms of comparative 2D forward modelling studies and have been applied to an offshore groundwater investigation from the Canterbury Basin, New Zealand. 2D inversion models from almost 300 km of data collected are presented and compared to stitched 1D inversion models. Overall, the 2D inversion models offer a more complete picture of the subsurface resistivity distribution and allow to map the lateral extent of the aquifer across the basin