Seismic Reﬂection Methods in O ﬀ shore Groundwater Research

: There is growing evidence that passive margin sediments in o ﬀ shore settings host large volumes of fresh and brackish water of meteoric origin in submarine sub-surface reservoirs. Marine geophysical methods, in particular seismic reﬂection data, can help characterize o ﬀ shore hydrogeological systems and yet the existing global database of industrial basin wide surveys remains untapped in this context. In this paper we highlight the importance of these data in groundwater exploration, by reviewing existing studies that apply physical stratigraphy and morpho-structural interpretation techniques to provide important information on—reservoir (aquifer) properties and architecture, permeability barriers, paleo-continental environments, sea-level changes and shift of coastal facies through time and conduits for water ﬂow. We then evaluate the scientiﬁc and applied relevance of such methodology within a holistic workﬂow for o ﬀ shore groundwater research.

The presence of low salinity water in the shallow offshore might produce an economic advantage for countries that rely heavily on desalinization as their main source of freshwater. The economic costs of this expensive process could considerably decrease by the use of extracted water with lower salt content. Therefore, there is a mounting interest in producing accurate maps of offshore groundwater resources, determine their likely responses to climate and sea-level changes and assess how much water they could sustainably supply, especially in proximity of water-stressed coastal cities [21].  [4]), with case studies shown in this paper. In red: reserves proven by observational data, in green: reserves with brackish component highlighted by pore water analysis, in blue: significant meteoric offshore groundwater with indirect evidence based on onshore paleo-groundwater [4]. In black dashed circles: case studies that use offshore geophysical data for groundwater analysis and are described in this paper. The yellow transparent overlay shows the highest density coverage of seismic reflection data, mostly from industry-funded and confidentiality protected sources. This map only provides an approximation of data distribution and it represents a sub-set of data globally acquired (source: https://www.globalseismiclibrary.com/, https://pwdatalibraries.maps.arcgis.com/).

Distribution of Offshore Groundwater
Offshore groundwater is composed of meteoric-derived water, compaction-driven fluids and thermobaric-derived fluids [22]. Shallow and proximal sediments (at 10 s to few 100 s m of seabed and burial depths) on continental shelves contain pore fluids that are dominated by water of meteoric origin, flowing in a basinward direction [23]. Deeper sediments up to few 1000 s of meters of burial depth contain pore fluids, of higher salinity, expelled by compaction in a landward and upward direction. Occasionally, the presence of fresh or low salinity fluids in the marine subsurface is linked to the expulsion of diagenetic fluids from sediments such as evaporites, clays and silica [24][25][26][27][28]. Thermobaric derived fluids are present at higher depths (see summary in Reference [22]).
Of these offshore groundwater systems, topographically driven meteoric recharge (TDM), that is, a combination of rainfall recharge and topographically driven flow, is the most important driver for the emplacement of fresh or low salinity (also defined as 'freshened') water in continental shelves and upper slopes of passive margins [29][30][31]. In these systems, the offshore reach of modern coastal aquifers into the submarine realm is constrained by a flow cell system with recharge from land. This depends on whether the onshore water table is high enough to provide the necessary driving force and whether there is hydraulic connectivity between the offshore aquifer and the onshore recharge area (see e.g., References [4,32]). Land-derived groundwater can be transported by conduits in sub-sea sediments or by pore-water flow in sediments of sufficient porosity and permeability (see e.g., Reference [33]). Groundwater may discharge directly on the seafloor through submarine groundwater discharge (SGD) [30,[34][35][36][37][38][39].
These processes alone, however, do not explain the large volumes of low-salinity groundwater that are found across continental shelves [4]. Offshore low-salinity groundwater accumulations, not necessarily connected to onshore recharge areas, can also be stored in confined paleo-aquifers. These are generally paleo-freshwater of meteoric origin, trapped in sediments during past times of lower sea level (see e.g., Reference [40]). These accumulations do not appear in global resources calculations but potentially also represent a large volume of deeply stored fresh-water.
Submarine groundwater is mostly studied along the shallower part of passive continental margins, with some notable exceptions on convergent margins [38,41]. On continental margins, submarine groundwater is usually stored in sedimentary rocks, in two main lithological groups-carbonate and siliciclastic. Both groups are covered by the case studies presented in this paper. The geological age of the sediments hosting the aquifers is variable but generally not older than Cenozoic. This could be related to data bias towards shallower sedimentary sequences, the higher potential for offshore aquifers hosting TDM and to the much higher probability of preservation/permanence of freshened water aquifers in the last few thousands of years [40,42].

Methods and Approaches to the Study of Offshore Groundwater
A combination of geological, geophysical, hydrogeological, petrophysical and hydrogeological data provide essential information for the study of offshore groundwater. These include: • Active source, multi-channel reflection seismic (hereby referred to simply as 'seismic data'); • Electromagnetics (EM, CSEM); • Seabed mapping tools (multibeam, water column backscatter, side-scan sonar and LIDAR); • Thermal remote sensing; • Sediment core and pore water geochemical analyses; • Borehole geophysics, imaging and hydrogeology (e.g., spontaneous potential, nuclear magnetic resonance, resistivity, flow-meters, etc.).
Offshore data acquisition presents logistical and technical challenges and it is generally more costly than onshore operations. For example, in clastic systems, coring is often problematic at shallow depth on continental shelves, due to the difficulty of drilling through poorly consolidated and coarse-grained sediments and possible presence of overpressured horizons [15]. Consequently, many studies of large-scale offshore groundwater circulation largely rely on modelling. However, accurate flow models primarily depend on the correct estimation of rock permeability and stratigraphic organization, which define the presence or absence of groundwater reservoirs and the modes of groundwater/seawater exchange.
Seismic reflection data, which are the focus of this paper, constitute an effective tool to map reservoir properties (especially in siliciclastic reservoirs) at a large scale, to reconstruct past sea level and environmental changes and support the analysis of the response of aquifers to these events. Thus, seismic reflection data can play an important role in the global identification of offshore freshened groundwater and potential resources, as one of the components of a workflow which includes the methods described above.
Our main aims are-to describe how seismic data can be used within this workflow, to study offshore groundwater and help identify potential low-salinity aquifers; to review the key case studies where seismic reflection data have been integrated into offshore groundwater analysis (either acquired on purpose or subsequently made available for such studies); and to discuss the limitations and future potential of such approach, both for scientific and applied purposes.

Seismic Reflection Data
Seismic data provide information on subsurface geology both in two dimensions (2D seismic), through individual vertical cross sections and in three-dimensions (3D seismic) through areal coverage of subsurface geological structures, generally down to kilometers in depth. Chirp profiling systems are used to acquire ultra-high resolutions for shallower depths [43][44][45][46][47]. In this review we present examples of both high-resolution, shallow penetration 2D datasets and lower resolution, deeper penetration 2D datasets.
Shallow, high-resolution (few centimeters to few meters) seismic data have been used, often serendipitously, to help characterize offshore groundwater systems, mostly located on passive continental margins [4,16,17]. The focus of this seismic method has been on stratigraphic architecture, facies prediction and identification of past sea level and environmental changes that might have affected paleo-aquifers. Conversely, deep penetration and lower resolution (in the order of 10 s of meters) seismic reflection data are routinely acquired and processed in industry (mostly oil and gas) and more occasionally, in academia. They have a deeper sub-surface imaging capability (several kilometers) and are used to interpret geology, structure and stratigraphy of continental margins (Figure 1), while they have been rarely applied to the analysis of offshore groundwater. In the few studies where this has been done, seismic data have been used for the interpretation of the large scale basin evolution and architecture, including the definition of lithological and stratigraphic units and structural elements, rather than direct imaging of aquifers. This global database already exists; therefore, it can potentially spur new studies in previously unexplored areas.
Seismic and sequence stratigraphic methods are broadly based on the identification and mapping of reflection geometries, terminations and seismic facies, systems tracts and sequences/sequence boundaries. Horizon-based interpretation can be combined with attribute analysis, both pre-and post-stack, which include trace attributes (e.g., amplitude, magnitude, continuity, frequency) or horizon/time attributes (e.g., dip, azimuth, curvature). Attribute analysis may support the interpretation of stratigraphic bodies, facies and discontinuities. Well tie, when available, leads to the cross-borehole physical properties of the sediments throughout the extent of the seismic data. Figure 2 shows an example of this methodology applied to offshore groundwater research. Quantitative seismic methods are instead used to understand amplitude variations, mostly relevant as a validation of hydrocarbon anomalies but also to provide additional information for reservoir characterization [68]. The most important of these techniques include post-stack amplitude analysis (bright-spot and dim-spot analysis for fluid content identification), offset-dependent amplitude analysis (AVO analysis), acoustic and elastic impedance inversion and forward seismic modeling [69].

Seismic Morpho-Structural Interpretation
On the seabed and in the subsurface, seismic data can image specific geomorphological and structural features which are relevant for the interpretation of the presence of freshened groundwater environments [47]. Examples are (paleo)channels and canyons [72][73][74] (Figure 3), sinkholes and other dissolution-related karstic features [75][76][77], pockmarks, pipes [78,79] and faults. Example of use of seismic stratigraphic and attribute analysis for lithology/facies prediction in groundwater research (modified after Reference [70]). (A) Amplitude of a seismic line across the New Jersey continental shelf (see case study in Section 3.1.1 of this paper). (B) Cosine of phase attribute with major systems tracts boundaries (including sequence boundaries and transgressive surfaces/maximum regression surfaces) labelled to highlight the correlation of key seismic sequence stratigraphic relationships. (C) Original systems tracts interpretation by Reference [70] overlain on seismic data. (D) Stochastically generated facies distribution. FSST falling stage systems tract, TST transgressive systems tract, HST highstand systems tract, LST lowstand systems tract. Horizons m1 to m4.1 correspond to interpreted sequence boundaries used to constrain the shelf model (nomenclature after [71]).
Buried channels and canyons can contribute to submarine groundwater discharge by serving as conduits for pore water flow [74], if they breach confining units and their infill is composed of permeable sediments. An example is offshore North Carolina, where based on simulation modeling and offshore geophysical data, the authors suggest that pore water of intermediate salinity is present in paleochannel sediments (Figure 3).
In karstic terrains, substratum dissolution can create sinkholes and other morphologies that favor groundwater flow. If the sinkholes are buried, seismic data can be used to recognize whether they are active or fossil and if they have been subject to multiple collapse phases ( Figure 4). The presence of poly-phase solution-subsidence⁄collapse structures associated to deep palaeo-karst conduits can be linked to repeated sea-level low-stands ( Figure 4) [76]. Eustatic variations can also favor dissolution, as changes in hydraulic gradient, related to the elevation difference, potentially lead to increase in karstification [80] (Figure 5). Limestone corrosion resulting from mixing of fresh and salt waters occurs in proximal offshore environments [81][82][83]. The presence of these features is detected on seismic data and can be used as an indication of paleo freshwater tables in limestones.
Geosciences 2020, 10, x FOR PEER REVIEW 6 of 34 Figure 3. Example of combination of high resolution seismic data (CHIRP seismic profiles, upper panels) and electromagnetic data (as apparent porosities, lower panels) illustrating the offshore extent of a paleochannel with no seabed expression and its relationship with groundwater (modified from Reference [74]). The data are from offshore North Carolina. In (A) (2.8 km offshore) the channel geometry is visualized through seismic data as concave upwards high amplitude reflections between the seabed and 0.025 s TWT. The related EM response of the channel shows an increase in porosity reflecting coarser grained material within the channel, particularly at the channel floor. (B) Transect located 1 km offshore, where the channel is located between the seabed and 0.015 s TWT, with an EM response more complex than (A), with an initial drop in porosity at the channel flanks, followed by an increase in mid-channel. The 40-m receiver shows a smooth decrease in porosity starting about 100 m either side of the channel. In the seismic sections, 6ms corresponds approximately to 5 m; resolution of data, considering frequency of 1-7 kHz and a typical shallow siliciclastic sediment velocity of 1600-1800 m/s, is in the order of few 10 s of cm.
Buried channels and canyons can contribute to submarine groundwater discharge by serving as conduits for pore water flow [74], if they breach confining units and their infill is composed of permeable sediments. An example is offshore North Carolina, where based on simulation modeling and offshore geophysical data, the authors suggest that pore water of intermediate salinity is present in paleochannel sediments (Figure 3). . Example of combination of high resolution seismic data (CHIRP seismic profiles, upper panels) and electromagnetic data (as apparent porosities, lower panels) illustrating the offshore extent of a paleochannel with no seabed expression and its relationship with groundwater (modified from Reference [74]). The data are from offshore North Carolina. In (A) (2.8 km offshore) the channel geometry is visualized through seismic data as concave upwards high amplitude reflections between the seabed and 0.025 s TWT. The related EM response of the channel shows an increase in porosity reflecting coarser grained material within the channel, particularly at the channel floor. (B) Transect located 1 km offshore, where the channel is located between the seabed and 0.015 s TWT, with an EM response more complex than (A), with an initial drop in porosity at the channel flanks, followed by an increase in mid-channel. The 40-m receiver shows a smooth decrease in porosity starting about 100 m either side of the channel. In the seismic sections, 6ms corresponds approximately to 5 m; resolution of data, considering frequency of 1-7 kHz and a typical shallow siliciclastic sediment velocity of 1600-1800 m/s, is in the order of few 10 s of cm. The spring discharge is visible in the seismic image and the associated collapse feature (modified from Reference [77]). The horizontal scale is missing in the original figure, however, as a reference, the active spring vent related depression is ca. 150 m wide at the sea floor. (B) In the Gulf of Lions, buried collapse feature on a CHIRP profile (vertical resolution 2 to 5.2 kHz, i.e., 30-50 cm), interpreted as solution-subsidence⁄collapse associated with deep paleo-karst conduits [76]. LGM = last glacial maximum.
In karstic terrains, substratum dissolution can create sinkholes and other morphologies that favor groundwater flow. If the sinkholes are buried, seismic data can be used to recognize whether they are active or fossil and if they have been subject to multiple collapse phases ( Figure 4). The presence of poly-phase solution-subsidence⁄collapse structures associated to deep palaeo-karst conduits can be linked to repeated sea-level low-stands ( Figure 4) [76]. Eustatic variations can also favor dissolution, as changes in hydraulic gradient, related to the elevation difference, potentially lead to increase in karstification [80] (Figure 5). Limestone corrosion resulting from mixing of fresh and salt waters occurs in proximal offshore environments [81][82][83]. The presence of these features is detected on seismic data and can be used as an indication of paleo freshwater tables in limestones. The spring discharge is visible in the seismic image and the associated collapse feature (modified from Reference [77]). The horizontal scale is missing in the original figure, however, as a reference, the active spring vent related depression is ca. 150 m wide at the sea floor. (B) In the Gulf of Lions, buried collapse feature on a CHIRP profile (vertical resolution 2 to 5.2 kHz, i.e., 30-50 cm), interpreted as solution-subsidence⁄collapse associated with deep paleo-karst conduits [76]. LGM = last glacial maximum.  [80]. This part of the seafloor was exposed for up to 100 ka during the last glacial cycle [84], allowing meteoric dissolution to take place. Changes in hydraulic gradients related to the elevation difference between the water table onshore and offshore terrain likely resulted in higher topographically-driven flow than today, which could have also contributed to karstification. Large non-eustatic sea level drops can also be recorded by coastal karstic springs and their linked conduits. An example is the dramatic drop in sea-level that occurred in the Mediterranean during the Messinian salinity crisis, ca. 6 Million years ago [85][86][87][88] which left seismic evidence of conduits  [80]. This part of the seafloor was exposed for up to 100 ka during the last glacial cycle [84], allowing meteoric dissolution to take place. Changes in hydraulic gradients related to the elevation difference between the water table onshore and offshore terrain likely resulted in higher topographically-driven flow than today, which could have also contributed to karstification. Large non-eustatic sea level drops can also be recorded by coastal karstic springs and their linked conduits. An example is the dramatic drop in sea-level that occurred in the Mediterranean during the Messinian salinity crisis, ca. 6 Million years ago [85][86][87][88] which left seismic evidence of conduits extending to a considerable depth below the last glacial maximum lowstand.
Pockmarks and pipes are diagnostic features for the occurrence of focused and vertical fluid flow [78]. This type of fluid flow is caused by upwards flowing gas and/or pore water [89,90]. Pockmarks are also present in lacustrine environments, such as in Lake Neuchatel (Switzerland), where they have been related to the expulsion of fresh pore waters and to active fluid flow by karstic groundwater discharge from the Jura Mountain front into the Swiss Plateau Hydrological system [91]. Finally, faults and associated structural interpretation of seismic data can indicate the presence of fractured reservoirs and potential flow conduits or lateral permeability barriers.

Case Studies
We present a series of key studies that show how seismic data have been used to aid interpretation of offshore groundwater systems, with a focus on low-salinity (fresh or freshened) aquifers. For each of these case studies, we focus on results or contributions from seismic stratigraphic methods and morpho-structural interpretation, based on the criteria and examples described in Sections 2.1 and 2.2.

North Atlantic Margin
The North Atlantic margin (location in Figure 1) is a site of primary interest in the field of submarine groundwater research. In this area, offshore groundwater bodies hosting large amounts of low salinity submarine pore water have been discovered [15,92].

New Jersey
In this case study, we show how seismic methods have been used for reservoir and seal characterization and fault conduit interpretation in a siliciclastic passive margin setting.
Geological and hydrogeological setting: The sedimentary units of the northern Atlantic passive continental margin (location in Figures 1 and 6) are composed of a wedge of Cretaceous through Pleistocene sediments. The shelf extends 200 km out up to the continental slope (e.g., References [93][94][95][96][97][98][99]. The deepest part of the sedimentary wedge consists of carbonates and minor amounts of evaporites [100]. The upper part consists of a siliciclastic prograding shelf and slope system built since the Oligocene. Beneath the coastal plain, a multilayered aquifer system saturated with low-salinity pore water is prevalently hosted in unconsolidated siliciclastic deposits of Miocene age [6,[101][102][103] ( Figure 6B). The upper unconfined aquifer locally contains brackish or salty water at the coastline and this saline contamination is due to onshore pumping and salt-water intrusion.
Offshore, beneath the shelf, the presence of pore water (fresher than seawater) is known from a series of boreholes completed during the 1970s and 1980s [92]. The modern onshore hydraulic system is thought to extend offshore, at least beneath the inner shelf where freshwater has been recovered from boreholes in both fine-grained and coarse-grained sediments [92]. Resistivity models derived from shallow water electromagnetic geophysical method support the presence beneath the inner shelf of a low salinity aquifer geometrically connected with the shore and overlain by a low-permeability clay unit acting as a confining unit ( Figure 6) [104]. The average porosities of sediments hosting low salinity waters range from 30% to 60% and permeability is estimated between 10 and 1000 mD [17]. The origin of the low salinity groundwater beneath the middle shelf of this area is still debated, with two prevailing hypotheses-first, that it is linked to present-day active dynamic connections with onshore aquifers [19] and second, that it is due to trapped meteoritic and/or sub-ice-sheet waters recharged during lowered sea-level periods [16,92,95,105].
Use of seismic data: The seismic and sequence stratigraphy of this margin is defined in detail with high-resolution chronostratigraphic control (e.g., References [71,106] (Figure 7). Seismic data interpretation shows that pore water salinity is primarily controlled by clinothem geometry, which determines the upwards extent of the low-salinity aquifer beneath the inner shelf ( Figure 6B).
Further offshore, in the middle shelf, recent IODP drilling showed that the organization of the offshore reservoirs is also multilayered, with low salinity water (<15 g/L) intervals of different thickness alternating vertically with salty intervals (salinity > 15 g/L) [15]. A direct relationship exists between the sedimentary facies of the siliciclastic system and the salinity distribution of the pore waters, with fresher water preferentially stored in the fine-grained sediments and salty water in coarse grained intervals [17]. Borehole-calibrated seismic data analysis ( Figure 6A) shows also that at these sites there is a good correlation between pore water salinity and seismic facies, with the exception of the emplacement of deep upwelling brines [73].
A package of discontinuous seismic reflections corresponds to the recent sand-dominated clinoform topsets, with locally interbedded clay layers and defines the uppermost 200 m-thick interval (ca. 0.3 s TWT). This interval hosts alternating salty and fresh pore-waters [17].
The deeper Miocene gently dipping seaward topsets are composed of continuous, high-amplitude reflections that are clay and silt dominated. These deposits preferentially host thick fresh-water intervals ( Figure 7). The upper discontinuous to steeply dipping foresets consist of unconsolidated clean sand deposits and preferentially store salty pore-waters. Finally, the lower part of the foresets is dominated by fine-grained deposits which are fresh-water bearing.
The relationship between salinity distribution and seismic facies allowed the development a 2D geometrical model of the reservoir and salinity distribution along a dip section of the shelf margin [17] (Figure 6). In summary, both beneath the inner and middle shelf, fresh water distribution is primarily controlled by clinothem geometry and associated lithology distribution ( Figure 7) [17]. However, beneath the inner shelf fresh water is stored in both fine-grained and coarse-grained sediments, whereas beneath the middle shelf, freshened groundwater is preferentially stored in deeper fine grained and low-permeability deposits (clay and silt dominated) and salty waters recovered in coarse-grained deposits (unconsolidated sands). Low permeability barriers such as fine-grained maximum flooding/transgression sediments and cemented diagenetic layers potentially play a role in preserving the emplaced fresh water [17,104,105,107]. Conversely, the faults interpreted in the deeper parts of the seismic sections and based on the results from the resistivity modelling, appear to connect the shallow freshened water aquifers with upwelling fluids from dissolved deep evaporites, facilitating their salinization. onshore and offshore boreholes with yellow and pink circles [15,17]. Location of seismic lines in the New Jersey ( Figure 7B,C) and Martha's Vineyard (Figures 8 and 9) case studies are indicated with red lines. (B) Seismic profile Oc270-529 and CH0698-19 (R/V Oceanus cruise 270 and R/V Cape Hatteras cruise 0698) combined to form a dip profile that intersects the IODP Expedition 313 holes (location in Figure 7A), allowing both detailed facies analysis and larger scale geometry analysis (AMCOR-Atlantic Margin Coring Project; TWT-two-way travel time; cdp-common depth point). Superimposed pore-water chlorinity from borehole data illustrate the clear correlation between seismic facies and salinity [17]. Seismic data vertical resolution is ca. 5 m with the best imaging between 0 and 1.3 s. (C) Resistivity models obtained from jointly inverting surface-towed CSEM and seafloor MT data for the shore-to-shelf profiles off New Jersey with seismic and well log data (modified from Reference [104]). Shaded colors show high resistivity as yellow hues and low resistivity as blue hues. White triangles show seafloor MT receiver locations and dashed white lines show the extent of the surface-towed CSEM data used in the inversions. R3 and R4 are resistive zones interpreted as low-salinity aquifers. Resistivity models overlay co-located seismic reflection images and pore water salinity data. Salinity data are plotted as white lines on linear scales. Black dashed lines represent a salinity value of 15. Seismically-imaged confining units and clinoforms that influence groundwater distribution patterns are shown as light grey dashed lines.

Martha's Vineyard (Massachusetts)
In this case study (location in Figures 1 and 6) seismic methods have been used to understand reservoir properties, aquifer onshore-offshore connectivity and paleo-coastline definition, in a siliciclastic, passive margin setting.
Geological and hydrogeological setting: The hydrogeology of the continental-shelf sedimentary units beneath Nantucket Island and offshore into Martha's Vineyard, present many analogies with the case study of New Jersey [95,96,[108][109][110][111][112] (location in Figures 1 and 6). The stratigraphic and hydrostratigraphic units span the Cretaceous to the Holocene ( Figure 7) and show significant variability in thickness downdip and along strike. On Nantucket Island, freshwater extends to depths greater than 500 m below sea level [93,113]. . Geo-hydrostratigraphic cross section for Nantucket, Massachusetts (modified from Reference [95]). The position and thickness of the Laurentide Ice Sheet at the Last Glacial Maximum at ca. 21 ka are estimated from Reference [114]. The cross section is constructed by using borehole lithologic and hydrogeologic data from References [94,115,116].
A 514-m-deep borehole (USGS 6001, from Reference [47]) (location in Figure 6A) penetrating the entire Cretaceous sedimentary package and two shallower wells (100 m), completed in the Miocene-Pleistocene sands, all had pore-water salinities of less than 1 ppt within the permeable intervals. The relatively thick clay and silt layers in well USGS 6001 exhibited higher salinity levels (30-70% seawater) and displayed a geometry consistent with ongoing vertical diffusion. The permeability measured on Figure 7. Geo-hydrostratigraphic cross section for Nantucket, Massachusetts (modified from Reference [95]). The position and thickness of the Laurentide Ice Sheet at the Last Glacial Maximum at ca. 21 ka are estimated from Reference [114]. The cross section is constructed by using borehole lithologic and hydrogeologic data from References [94,115,116].
A 514-m-deep borehole (USGS 6001, from Reference [47]) (location in Figure 6A) penetrating the entire Cretaceous sedimentary package and two shallower wells (100 m), completed in the Miocene-Pleistocene sands, all had pore-water salinities of less than 1 ppt within the permeable intervals. The relatively thick clay and silt layers in well USGS 6001 exhibited higher salinity levels (30-70% seawater) and displayed a geometry consistent with ongoing vertical diffusion. The permeability measured on Nantucket Island in low-salinity aquifers range from 10-780 mD and it allows flow of water through sediments [98,107]. The exact extent of the low-salinity groundwater offshore Massachusetts isis unknown, as no offshore drilling has confirmed the predicted models. Similarly to the New Jersey case study, seismic reflection images co-rendered on resistivity models ( Figure 9) reveal consistent architectural controls on the low-salinity groundwater distribution [104]. Continuous, well-defined stratigraphic boundaries mark the upward extent of low-salinity groundwater and the onset of Oligo-Miocene clinothem deposition is thought to limit the seaward extent of the resistive bodies that are interpreted as fresh-water bearing. The authors suggest that the

Use of seismic data:
Recent studies [117,118] use high-resolution, multichannel 2D seismic data to define the stratigraphy and constrain a three-dimensional, variable-density model that couples fluid flow, heat and solute transport for the continental shelf ( Figure 9). Clastic Oligo-Miocene units are separated from the overlying Pleistocene-Holocene units by a prominent glacial unconformity (U1 in Figure 8), interpreted as marine oxygen isotope stage (MIS) 12 [119]. U2 ( Figure 8) is a shallow sequence boundary that formed during the last sea level fall (40-30 ka B.P.) and marks the onset of siliciclastic sedimentation and glacial outwash during the late Pleistocene and Holocene [117,120]. According to the authors, freshwater distribution on this margin is strongly dependent on depositional architecture and ice sheet history. Model results indicate that the late Pleistocene ice sheet was responsible for much of the emplaced freshwater. Through seismic data analysis the authors are able to estimate the maximum extent of the late Pleistocene ice sheet to near the shelf-slope break (U2 in Figure 8).
Similarly to the New Jersey case study, seismic reflection images co-rendered on resistivity models ( Figure 9) reveal consistent architectural controls on the low-salinity groundwater distribution [104]. Continuous, well-defined stratigraphic boundaries mark the upward extent of low-salinity groundwater and the onset of Oligo-Miocene clinothem deposition is thought to limit the seaward extent of the resistive bodies that are interpreted as fresh-water bearing. The authors suggest that the resistive body R1 (freshwater aquifer) ( Figure 10) connects onshore, as shown by comparison with the well ENW-50 and is thus consistent with seaward transport of freshwater from onshore recharge through coarse sediments.
Geosciences 2020, 10, x FOR PEER REVIEW 14 of 34 Figure 9. Resistivity models obtained from jointly inverting surface-towed CSEM and seafloor MT data for the shore-to-shelf profiles of Martha's Vineyard with seismic data. Shaded colors show high resistivity as yellow hues and low resistivity as blue hues. White triangles show seafloor MT receiver locations and dashed white lines show the extent of the surface-towed CSEM data used in the inversions. Co-located seismic reflection image as well as pore water salinity data are also shown. Salinity data are plotted as white lines on linear scales. Onshore well ENW-50 salinity data are plotted in yellow. Black dashed lines represent a salinity value of 15. Seismically imaged confining units and clinoform structures that influence groundwater distribution patterns are shown as labeled light grey dashed lines (modified from Reference [104]). Location of line in Figure 6A.

Florida
In this case study, seismic methods have been used for characterization of reservoir properties, aquifer onshore-offshore connectivity, paleo-continental environments, identification of dissolutionrelated karstic features and associated fossil aquifers.
Geological and hydrogeological setting: The Floridan Plateau is a broad, flat limestone and dolomite platform separating the Atlantic Ocean from the Gulf of Mexico [121] (Figure 10). This Plateau hosts the Floridan aquifer, one of the most productive aquifer systems worldwide ( Figure 10). In northeastern Florida, water within this aquifer system can be under artesian pressure, depending on the permeability of overlying strata [122][123][124]. The aquifer can be subdivided into two moderately high Figure 9. Resistivity models obtained from jointly inverting surface-towed CSEM and seafloor MT data for the shore-to-shelf profiles of Martha's Vineyard with seismic data. Shaded colors show high resistivity as yellow hues and low resistivity as blue hues. White triangles show seafloor MT receiver locations and dashed white lines show the extent of the surface-towed CSEM data used in the inversions. Co-located seismic reflection image as well as pore water salinity data are also shown. Salinity data are plotted as white lines on linear scales. Onshore well ENW-50 salinity data are plotted in yellow. Black dashed lines represent a salinity value of 15. Seismically imaged confining units and clinoform structures that influence groundwater distribution patterns are shown as labeled light grey dashed lines (modified from Reference [104]). Location of line in Figure 6A. The submarine karst morphologies and the narrow and linear distribution of sinkholes on the two terraces may indicate the past or present location of mixing zones, where limestone dissolution and solution-collapse processes are enhanced [75]. The formation of some of the sinkholes observed in the Straits of Florida has been related to late Miocene-Pleistocene sea-level lowstands resulting in an elevated head and stimulating artesian flow of groundwater out of the aquifer [128,129].
Subaqueous formation of sinkholes in the same region is also associated with submarine discharge of fresh groundwater flowing from the nearby platform [81] (Figure 10). These sinkholes are located on the seabed at water depths greater than 600 m, which is deeper than any Neogene eustatic change [124]. Glacial sea-level lowstands during the Pleistocene would have presumably stimulated submarine artesian flow in the area and caused the sinkhole formation, with only a 10m increase in head within the onshore Floridan aquifer [81].

Mediterranean Region
The Mediterranean region ( Figure 11) offers a series of examples of coastal aquifers and their shallow offshore extension, in both karstic and clastic (carbonate and terrigenous) sediments. We show a case study in Southern France, where seismic data derived from both academia (shallow penetration, high resolution) and oil and gas industry (deeper penetration, lower resolution) have been used for reservoir architecture and characterization, pending ground-truthing from borehole sampling.  [75,81]). B and C indicate the approximate position of two sinkholes that are observed on bathymetric imaging in the Miami terrace and are currently located at 200-400 m water depth and extend 10s of kilometers into the offshore [75]. Resolution is not indicated in the cited work but it is likely to be in the order of 10s of cm to m considering the 7 kHz source of the seismic profiler.

Florida
In this case study, seismic methods have been used for characterization of reservoir properties, aquifer onshore-offshore connectivity, paleo-continental environments, identification of dissolution-related karstic features and associated fossil aquifers.
Geological and hydrogeological setting: The Floridan Plateau is a broad, flat limestone and dolomite platform separating the Atlantic Ocean from the Gulf of Mexico [121] (Figure 10). This Plateau hosts the Floridan aquifer, one of the most productive aquifer systems worldwide ( Figure 10). In northeastern Florida, water within this aquifer system can be under artesian pressure, depending on the permeability of overlying strata [122][123][124]. The aquifer can be subdivided into two moderately high permeability systems, referred to as the Upper and Lower Floridan. A series of relatively impermeable Miocene and younger beds confine the aquifer. The continental shelf is about 110-km wide [125] and a large volume of freshened groundwater under artesian pressure is known to exist beneath the Atlantic Ocean [100,126]. A geologic exploratory well drilled 43 km offshore from Jacksonville penetrated an artesian aquifer at 250 m depth below sea level. Historic submarine discharge has been observed to occur off the coast of Florida [121,127], with active springs, located a few kilometers offshore, associated fractures and sinkholes ( [122] and references therein).
Use of seismic data: Seismic data have been used in this region both to identify karstic vents on the seafloor and karstic morphology in the subsurface. A fresh water discharge vent is visible on seismic profiles, close to an ancient collapse feature ( Figure 4A) [77]. In the Straits of Florida, in the southern part of the case study area, detailed bathymetric and high-resolution seismic records reveal well-defined karstic morphology and multiple large collapse features on the Miami and Poutales Terraces, currently located at 200-400 m water depth and at 10 s to 100 km from the coastline ([75], Figure 10). These features are linked to sustained submarine discharge of artesian water [75,81] (see also Section 2.2).
The submarine karst morphologies and the narrow and linear distribution of sinkholes on the two terraces may indicate the past or present location of mixing zones, where limestone dissolution and solution-collapse processes are enhanced [75]. The formation of some of the sinkholes observed in the Straits of Florida has been related to late Miocene-Pleistocene sea-level lowstands resulting in an elevated head and stimulating artesian flow of groundwater out of the aquifer [128,129].
Subaqueous formation of sinkholes in the same region is also associated with submarine discharge of fresh groundwater flowing from the nearby platform [81] (Figure 10). These sinkholes are located on the seabed at water depths greater than 600 m, which is deeper than any Neogene eustatic change [124]. Glacial sea-level lowstands during the Pleistocene would have presumably stimulated submarine artesian flow in the area and caused the sinkhole formation, with only a 10 m increase in head within the onshore Floridan aquifer [81].

Mediterranean Region
The Mediterranean region ( Figure 11) offers a series of examples of coastal aquifers and their shallow offshore extension, in both karstic and clastic (carbonate and terrigenous) sediments. We show a case study in Southern France, where seismic data derived from both academia (shallow penetration, high resolution) and oil and gas industry (deeper penetration, lower resolution) have been used for reservoir architecture and characterization, pending ground-truthing from borehole sampling. Geological and hydrogeological setting: The majority of Mediterranean coastal porous aquifers are hosted in clastic sediments and are located in Pliocene to Quaternary deposits [9,[132][133][134][135]. These aquifers are "coastal detritic formations" following the classification for European aquifers by Reference [136] ( Figure  12). They are often multi-layered and mainly studied onshore. As littoral zones are important places of economic development, they are often linked to intensive water needs for agriculture, drinking water supply, tourism and industry. These aquifers are extremely sensitive to sea-level and climatic changes, due to the Mediterranean Sea configuration as a landlocked basin (e.g., Reference [137]). As a consequence of these factors and of over-exploitation, the quality of coastal aquifers in the region is currently heavily affected by seawater intrusion (Figure 11) [2, 138]. Seawater intrusion at shallow depth can also occur naturally as a result of offshore strata geometries ( Figure 12A).
In Southern France, the Plio-Quaternary calcareous sands are one of the most important aquifers for the freshwater supplies in the Gulf of Lions coastal areas, such as the complex multi-layered aquifer located below the Roussillon coastal plain and potentially extending into the offshore (Figure 12). This aquifer is hosted in marine and continental Pliocene and Quaternary coarse grained sediments. Each aquifer layer is characterized by different hydrogeological properties [132]. Geochemical and isotopic analyses performed on groundwater samples taken from two coastal boreholes [135] showed groundwater residence ranging from about 15 ky in the deep Pliocene aquifer to 5 ky in the intermediate Pliocene aquifer [135].
Use of seismic data: The extent of coastal aquifers into the offshore clastic Plio-Pleistocene shelfal reservoirs has been mapped using high-resolution, shallow seismic and deep, lower resolution industry data, tied to onshore boreholes [73] (Figure 12). The offshore record mostly consists of prograding Geological and hydrogeological setting: The majority of Mediterranean coastal porous aquifers are hosted in clastic sediments and are located in Pliocene to Quaternary deposits [9,[132][133][134][135]. These aquifers are "coastal detritic formations" following the classification for European aquifers by Reference [136] (Figure 12). They are often multi-layered and mainly studied onshore. As littoral zones are important places of economic development, they are often linked to intensive water needs for agriculture, drinking water supply, tourism and industry. These aquifers are extremely sensitive to sea-level and climatic changes, due to the Mediterranean Sea configuration as a landlocked basin (e.g., Reference [137]). As a consequence of these factors and of over-exploitation, the quality of coastal aquifers in the region is currently heavily affected by seawater intrusion (Figure 11) [2, 138]. Seawater intrusion at shallow depth can also occur naturally as a result of offshore strata geometries ( Figure 12A).
In Southern France, the Plio-Quaternary calcareous sands are one of the most important aquifers for the freshwater supplies in the Gulf of Lions coastal areas, such as the complex multi-layered aquifer located below the Roussillon coastal plain and potentially extending into the offshore (Figure 12). This aquifer is hosted in marine and continental Pliocene and Quaternary coarse grained sediments. Each aquifer layer is characterized by different hydrogeological properties [132]. Geochemical and isotopic analyses performed on groundwater samples taken from two coastal boreholes [135] showed groundwater residence ranging from about 15 ky in the deep Pliocene aquifer to 5 ky in the intermediate Pliocene aquifer [135].
Use of seismic data: The extent of coastal aquifers into the offshore clastic Plio-Pleistocene shelfal reservoirs has been mapped using high-resolution, shallow seismic and deep, lower resolution industry data, tied to onshore boreholes [73] (Figure 12). The offshore record mostly consists of prograding clinothems. As for the New Jersey margin case study, beneath the coastal plain, a multilayered aquifer system saturated with fresh pore water is hosted in unconsolidated siliciclastic sediments. The upper part may locally contain brackish or salty water at the coastline. The present day extent of the low-salinity groundwater beneath the shelf is currently uncalibrated. However due to the similarity in geological setting and hydrogeological conditions with the New Jersey shelf, it can be inferred that fresh water distribution might be primarily controlled by clinothem geometry and associated lithology distribution.

East Africa
This case study demonstrates the use and value of oil & gas exploration data for deep groundwater exploration purposes in a coastal setting in Tanzania. Models obtained from seismic reflection data were used to define the geometry and reservoir architecture of a deep coastal aquifer system (the Kimbiji Aquifer) near the city of Dar es Salaam [139,140]. The geological model of the Kimbiji Aquifer is verified by drilling, both onshore and in the shallow offshore. The hydrogeological model is verified by drilling and testing activity onshore and is guided offshore by established hydrogeological principles and numerical groundwater modelling techniques.
Geological and hydrogeological setting: The Kimbiji Aquifer comprises a thick sequence of paleo-deltaic sediments of mainly Miocene age [20,[139][140][141] overlying Lower Tertiary carbonates ( Figure 13). The sediments, which are associated with a rapidly subsiding basin, are layered and heterogeneous in nature and incorporate sands, gravels, clays and marls.
Seismic reflection profiles show how the stratigraphic sequence hosting the Kimbiji Aquifer onshore extends in the offshore domain, where it thickens beneath the Indian Ocean. At the coastline, the verified thickness is 1.5 km [140].  [46]. In the dip profile, the Upl (upper Pliocene) reflections are gently dipping toward the sea at an angle of approximately 0.5°. The strike profile images the details of the Upl unit, which consist of low frequency, sub-parallel continuous reflections. Upl is eroded at the top by the RES and deeply incised in the axis of paleo-valleys inherited from the last glacial maximum. The Holocene permeable sands discordantly overlie the RES. Upl outcrops locally at the sea floor. Beneath the lido, salt water saturates the upper ~30 m and is thought to result from a connection of the aquifers with the sea either directly in places where Upl outcrops offshore or indirectly via the paleo-channels that may have breached a confining unit. (B) Pliocene-Quaternary correlations between onshore and offshore (seismic profile LRM04) formations in the Roussillon Basin area (modified from Reference [131]). Plio-Quaternary stratigraphy and systems tracts: PR1-8 (PR1-3: Zanclean Gilbert deltas foresets and bottomsets, PR4-6: Zanclean topsets, foresets and bottomsets, Pr7: Piacenzian topsets, foresets and bottomsets, Pr8: Gelasian/Quaternary topsets, foresets and bottomsets). (C) Hydrogeological conceptual model of the Roussillon aquifer system without scale and location [130]. Sea water intrusion occurs in the unconfined superficial Quaternary aquifer in the offshore domain (more than 5 km from the  [46]. In the dip profile, the Upl (upper Pliocene) reflections are gently dipping toward the sea at an angle of approximately 0.5 • . The strike profile images the details of the Upl unit, which consist of low frequency, sub-parallel continuous reflections. Upl is eroded at the top by the RES and deeply incised in the axis of paleo-valleys inherited from the last glacial maximum. The Holocene permeable sands discordantly overlie the RES. Upl outcrops locally at the sea floor. Beneath the lido, salt water saturates the upper~30 m and is thought to result from a connection of the aquifers with the sea either directly in places where Upl outcrops offshore or indirectly via the paleo-channels that may have breached a confining unit. (B) Pliocene-Quaternary correlations between onshore and offshore (seismic profile LRM04) formations in the Roussillon Basin area (modified from Reference [131]). Plio-Quaternary stratigraphy and systems tracts: PR1-8 (PR1-3: Zanclean Gilbert deltas foresets and bottomsets, PR4-6: Zanclean topsets, foresets and bottomsets, Pr7: Piacenzian topsets, foresets and bottomsets, Pr8: Gelasian/Quaternary topsets, foresets and bottomsets). (C) Hydrogeological conceptual model of the Roussillon aquifer system without scale and location [130]. Sea water intrusion occurs in the unconfined superficial Quaternary aquifer in the offshore domain (more than 5 km from the coastline, based on morphobathymetric map observations). Submarine groundwater discharge could occur in the open sea, although it has not been directly observed.
The complex, layered aquifer system becomes increasingly confined with depth. Deep artesian heads were encountered in several exploration wells, both near the coast and on structurally elevated positions within the Kimbiji Aquifer. Onshore drilling campaigns in 2006-2007 and 2013-2017 demonstrated that high-yielding production wells can be developed, abstracting low-salinity groundwater of potable quality. The apparent age of groundwater that was pumped from depths of 600 m was approximately 40,000 y, which is consistent with simulated flow velocities using a regional numerical groundwater model.
Numerical groundwater modelling also indicates that large-scale abstraction of up to 200,000 m 3 /d may be sustained without inducing saline intrusion for a period of 100 years [140]. Under a hydrogeological scenario that is defined by current, measured groundwater and seawater levels, the simulated fresh-saline groundwater interface extends approximately 25 km offshore (Figure 14). Under a different scenario which accounts for the deeper seawater levels that prevailed during the Pleistocene and Holocene [142,143], a 19,000-year simulation of subsequent sea level rise simulated an interface which extends approximately 40 km offshore ( Figure 14). Thus, a significant volume of fresh groundwater is inferred to extend a significant distance offshore.   Figure 13C. The Kimbiji Aquifer is hosted in deltaic sediments of (mainly) Miocene age. A calibrated high-amplitude seismic reflector representing top Eocene carbonates served as a marker horizon (orange) to define the base of the Kimbiji Aquifer system. The stratigraphic sequence is extensively faulted, due to the development of the East African Rift system and giving rise to a sedimentary basin, which thickens towards and extends beneath the Indian Ocean. (Seismic data courtesy of F. Ruden at Ruden AS) (Figures modified from: [139,141,144]).  [140]). (B) Simulated extent of the fresh-saline groundwater interface along cross-section A-A', whereby the red "conservative scenario" is defined by current measured groundwater and sea water levels and the blue "historical scenario" accounts for lower sea levels during the Pleistocene and Holocene.

Canterbury Bight, New Zealand
In this case study, shallow penetration and high resolution seismic data have been used for reservoir characterization and identification of faults and buried valleys in a siliciclastic passive margin setting.
Geological and hydrogeological background: An offshore freshened groundwater system hosted in unconsolidated, clastic sediments has been reported in the Canterbury Bight (located off the eastern coast of the South Island of New Zealand) [31] (location in Figure 1; Figure 15). The Canterbury Bight is a 95 km wide and up to 140-150 m deep continental shelf. It comprises a 1 km thick progradational succession of shelf-slope deposits punctuated by advances of the braid plain during periods of sea level fall [145]. Eustatically-driven transgressive-regressive cycles have controlled the sedimentary history of the Canterbury Bight since the middle Miocene [146][147][148].
Use of seismic data: Seismic reflection profiles across the Canterbury Bight were divided into five different facies using multi-attribute seismic facies classification ( Figure 15C-F). The latter entailed two  [140]). (B) Simulated extent of the fresh-saline groundwater interface along cross-section A-A', whereby the red "conservative scenario" is defined by current measured groundwater and sea water levels and the blue "historical scenario" accounts for lower sea levels during the Pleistocene and Holocene.
Extending east from Dar es Salaam and the Kimbiji Aquifer region, the area between coastal mainland Tanzania and the island of Unguja (Zanzibar) is characterized by basin thicknesses up to 6000 m of Neogene/Pleistocene/Holocene sands ( Figure 13). The basin axis is oriented roughly N-S, with a water depth rarely exceeding 40 m and with a total basin thickness exceeding 10,000 m in places. Logs from hydrocarbon exploration wells in the West Zanzibar Basin indicate bulk resistivity values in sand sequences pointing to the potential presence of fresh water [141].
Use of seismic data: Hydrocarbon exploration data are not conventionally or traditionally used in water exploration activity. In the case of the Kimbiji Aquifer, both the discovery [139] and study phase of the aquifer system [140] relied on access to oil & gas exploration data. The seismic reflection data, especially, was instrumental in defining the scale and magnitude of the aquifer system and its possible connection with deeper geological formations. The interpretation of the seismic data were aided by oil & gas well completion records, which included both geological and wireline logs, as well as drill stem test data.
In the case of the Kimbiji Aquifer study, the access to 1140 km of seismic data was facilitated by the Dar es Salaam Water and Sewerage Authority under an agreement with the Tanzanian Petroleum Development Corporation. The data that could be accessed was from the early 1980s and although data and record formats were of variable nature and quality, careful review and analysis allowed older records to be 'recycled' for water exploration purposes.

Canterbury Bight, New Zealand
In this case study, shallow penetration and high resolution seismic data have been used for reservoir characterization and identification of faults and buried valleys in a siliciclastic passive margin setting.
Geological and hydrogeological background: An offshore freshened groundwater system hosted in unconsolidated, clastic sediments has been reported in the Canterbury Bight (located off the eastern coast of the South Island of New Zealand) [31] (location in Figure 1; Figure 15). The Canterbury Bight is a 95 km wide and up to 140-150 m deep continental shelf. It comprises a 1 km thick progradational succession of shelf-slope deposits punctuated by advances of the braid plain during periods of sea level fall [145]. Eustatically-driven transgressive-regressive cycles have controlled the sedimentary history of the Canterbury Bight since the middle Miocene [146][147][148].
Use of seismic data: Seismic reflection profiles across the Canterbury Bight were divided into five different facies using multi-attribute seismic facies classification ( Figure 15C-F). The latter entailed two steps: (i) sub-division based on amplitude characteristics, lateral continuity, reflector geometry and two seismic attributes-instantaneous frequency and envelope [149]; (ii) determination of a depth-travel time relationship from the sonic logs and correlation of the features in the borehole logs, recorded in the depth domain, with those in the seismic reflection data, recorded in the time domain; a synthetic seismogram was constructed in the boreholes from the sonic log and the density curve was calculated from the resistivity log using Archie's relationship [150]. Four of the five facies (clay, silt, fine sand and coarse sand) could be correlated with lithologies in the boreholes, whereas the fifth one was interpreted as gravel in view of the similarity of the seismic facies with that of gravel onshore.
The stratigraphic framework that emerges from the classification of the seismic reflection profiles is an alternation of lowstand fluvial gravels and sands, which become thicker towards the shore and highstand sands, silts and clays, which are more dominant in the deeper sections. This distribution of facies is consistent with that of sediment transport models [145].
The CSEM resistivity data, locally calibrated by pore water salinity measurements from boreholes, point to an offshore groundwater system that consists of one main and two smaller, low salinity groundwater bodies ( Figure 15G,H). The main body has a maximum thickness of at least 250 m and extends up to 60 km from the coast and a seawater depth of 110 m. The maximum groundwater volume estimated from the CSEM resistivity data for a porosity of 40% is 213 km 3 .
The offshore freshened groundwater occurs in sedimentary layers that mainly include silt and sand and occasionally gravel and clay. More importantly, a recent study [31] reports a spatial variability in groundwater salinity along the shelf. This can be attributed to specific sedimentary structures (e.g., high permeability conduits and corridors; the latter are visible as 5-10 km long isolated bodies of parallel, continuous, high amplitude reflectors occurring at multiple depths and interpreted as buried valleys) and normal faults primarily acting as barriers and potentially to recharge from rivers during sea level lowstands.  A meteoric origin of the offshore freshened groundwater and active migration from onshore are inferred based on the geochemical characteristics of the pore water in the borehole data and the recent groundwater age at the coast. However, 2D models of groundwater flow and solute transport, based on the interpretation of the seismic reflection profiles, indicate that recharge from onshore aquifers can only account for a small fraction of the offshore freshened groundwater at present. The majority of the offshore freshened groundwater was emplaced during the last three glacial cycles. Topographically-driven shore-normal flow was higher than at present during sea level lowstands, due to an increase in the hydraulic head and steep onshore gradients. Such flow was reduced during sea level highstands, when lateral differences in salinity on the shelf drove groundwater laterally or shoreward.

Gippsland Basin, Australia
In this case study, seismic data provided the framework for the characterization of aquifer/reservoir architecture, mapping of regional hydrostratigraphic units at a basinal scale and their onshore-offshore connectivity, in a mixed carbonatic-siliciclastic setting.
Geological and hydrogeological background: Evidence of aquifers connected to onshore groundwater have been found offshore southern Australia [151,152] (Figure 16, general location in Figure 1), through the use of combined seismic data, onshore and offshore wells, historical production and pressure data. Data have been acquired through the years in successive campaigns focused on oil and gas exploration and the mining industry. This case study thus offers an example of a multiple usage aquifer and the consequences of its exploitation-significant withdrawal of fluid resulted in a complex pattern of pressure decline in the aquifer, potentially associated with its usage [153].

Gippsland Basin, Australia
In this case study, seismic data provided the framework for the characterization of aquifer/reservoir architecture, mapping of regional hydrostratigraphic units at a basinal scale and their onshore-offshore connectivity, in a mixed carbonatic-siliciclastic setting.
Geological and hydrogeological background: Evidence of aquifers connected to onshore groundwater have been found offshore southern Australia [151,152] (Figure 16, general location in Figure  1), through the use of combined seismic data, onshore and offshore wells, historical production and pressure data. Data have been acquired through the years in successive campaigns focused on oil and gas The basin extends from onshore to offshore and hosts depositional sequences ranging in age from Early Cretaceous to Holocene. The main deep aquifer is hosted in the Latrobe Group sediments, of Cretaceous to Oligocene age, composed of siliciclastic sediments, ranging from alluvio-fluvial sediments at the base to marginal-marine, deltaic and increasingly marine sediments, sealed by the confining aquitard of the Lake Entrance formation (marine shales of late Oligocene age) and by the Strzelecky aquitard below (sandstones and shales of Albian age) ( Figure 16) [151,152,154]. The The basin extends from onshore to offshore and hosts depositional sequences ranging in age from Early Cretaceous to Holocene. The main deep aquifer is hosted in the Latrobe Group sediments, of Cretaceous to Oligocene age, composed of siliciclastic sediments, ranging from alluvio-fluvial sediments at the base to marginal-marine, deltaic and increasingly marine sediments, sealed by the confining aquitard of the Lake Entrance formation (marine shales of late Oligocene age) and by the Strzelecky aquitard below (sandstones and shales of Albian age) ( Figure 16) [151,152,154]. The Latrobe Group forms a major freshwater aquifer (Latrobe aquifer) in the onshore part of the basin. This sedimentary unit also contains substantial coal resources that are being mined in the onshore part of the basin [151].
Use of seismic data: Approximately two thirds of the basin are located offshore. The Latrobe group has been identified as the most important reservoir for oil and gas in both onshore and offshore parts of the basin. Seismic data have been used to reconstruct stratigraphic cross sections and define the vertical and lateral extent of the Latrobe reservoir/aquifer and the various hydrostratigraphic units in the offshore.
The salinity of the formation water in the Latrobe aquifer varies from less than 0.5 g/L in the onshore area to more than 30 g/L total dissolved solids (TDS) in the Central Deep region of the offshore Gippsland Basin. In the deep Latrobe aquifer a lower salinity gradient and wide zone of brackish formation water extends to about 20 km offshore. Based on both wells and seismic data, the freshwater hydraulic-head distributions and the architecture of the sedimentary units suggest that groundwater flows from the northern and western margins towards the central part of the basin. Freshwater hydraulic heads in the western part of the basin imply that there is a source of water in the Central Deep where the Latrobe Aquifer occurs at depths > 2500 m below sea level and subcrops beneath the seafloor.

Summary of Use and Limitations of Seismic Reflection Data in Offshore Freshened Water Research
The selected case studies highlight how seismic data has been used, indirectly, in the identification and mapping of geological features hosting aquifers in the offshore realm. Seismic data is an important tool in a workflow that includes other geophysical and geological methods (see Section 1.2), in order to understand the quality (salinity) of offshore groundwater. In this context, lower resolution, deep penetration seismic data would be applied, initially, in an exploratory and scouting phase, with higher resolution and shallow penetration data subsequently applied to define the details of shallow aquifers and their geometry in any prospective area. Therefore, existing deeper penetration industry data can be used to widen the search for offshore freshened groundwater in a global context, if implemented by additional data and techniques.
So far, in the described examples and other case studies worldwide, seismic stratigraphic and sequence stratigraphic interpretation has predominantly provided value for the analysis of ( Figure 17): • Reservoir (aquifer) properties and architecture, for example, through spatial mapping of the permeability/porosity of sedimentary bodies, analysis of onshore-offshore aquifer connectivity, identification of the open or confined nature of offshore aquifers and their active or 'fossil' setting.
In carbonate rocks, seismic data can also indicate areas of porosity and permeability enhancement due to karstification processes (Sections 3.1-3.5) • Seal (aquitard/aquiclude), through the identification of permeability barriers, such as laterally extensive fine-grained sediments or tightly cemented layers (Section 3.1.1) • Paleo-continental environments, through the detection of potential freshwater-bearing environments, such as fluvio-lacustrine systems or subaerially exposed karstic terranes (Section 3.1.3), • Paleo-coastline and change to active aquifers: absolute and relative sea-level changes and shift of coastal facies through time and therefore potentially associated coastal aquifers (e.g., References [44,67]); The onlap shifts related to sea-level change are used to document shifting of the mixing zone and saltwater intrusion into coastal plain groundwater systems in the past (Section 3. The results of the considered case studies are currently biased towards aquifers hosted in siliciclastic sedimentary units, as these settings have been better surveyed and calibrated compared to offshore karstic aquifers, whose architecture remains much less predictable. In clastic continental margins settings, clinothems exert an important structural control on the lateral and vertical extent of low-salinity groundwater (e.g., New Jersey and Martha's Vineyard case studies). Coarser grainsize sediments, with high porosity and permeability, are naturally suited for storage of buried fluids (e.g., Section 3.1.2). According to facies and sequence stratigraphic models, these sediments are most likely to be deposited during regressions and are often separated by maximum flooding surfaces. Clinoform topsets or other high porosity and permeability bodies (such as basin floor fans) are reservoirs for low salinity waters. This is the case in the Pliocene sands on the Mediterranean offshore Southern France and on the New Jersey inner shelf. However, due to their high permeability, they are also subject to density driven sea-water intrusion, as it is the case in the New Jersey middle shelf, that is, in a more distal area with limited connection to the onshore recharge.
In the case on the New Jersey inner shelf, the modern onshore hydraulic system is thought to extend offshore in both the coarse and fine grained fraction of the clinothems. Further offshore, the presence of salt water in the coarse grained fractions of the topsets has tentatively been explained by seawater penetration as a result of fast vertical density-driven flow [17]. In this case, the vertical flow would be expected to be greater than any horizontal, land-to-sea directed flow within the connected reservoirs. Diffusion is also sometimes accounted for the low salinity water hosted in mid-low permeability sediments including silts and clays. In summary, the emplacement and preservation of bodies of freshened groundwater is complex and depends on the timing of emplacement, 3D reservoir geometry, present-day connection to the sea and the presence of tight layers which may aid the preservation of freshwater aquifers.
As a general rule, the preservation of low-salinity pore waters is inversely proportional to the age/depth of sediments. Well defined stratigraphic boundaries such as unconformities, lateral and vertical facies changes and more specifically, maximum flooding and transgressive surfaces marking a shift to fine-grained facies or diagenetic boundaries, can seal the aquifer and constrain the upward extent of low-salinity aquifers. The role of low-permeability barriers (intra or extra-formational) to separate superimposed aquifers is shown in the Tanzania/Zanzibar and New Jersey case studies. This phenomenon has also been observed elsewhere, such as in Hong Kong [156][157][158].
Evaporites, which are considered a good seal for petroleum fluids, when in contact with undersaturated pore fluids, are prone to quick dissolution and produce salinization of offshore freshened groundwater [159,160]. This upwelling-driven salinization process is recorded for Figure 17.
The results of the considered case studies are currently biased towards aquifers hosted in siliciclastic sedimentary units, as these settings have been better surveyed and calibrated compared to offshore karstic aquifers, whose architecture remains much less predictable. In clastic continental margins settings, clinothems exert an important structural control on the lateral and vertical extent of low-salinity groundwater (e.g., New Jersey and Martha's Vineyard case studies). Coarser grain-size sediments, with high porosity and permeability, are naturally suited for storage of buried fluids (e.g., Section 3.1.2). According to facies and sequence stratigraphic models, these sediments are most likely to be deposited during regressions and are often separated by maximum flooding surfaces. Clinoform topsets or other high porosity and permeability bodies (such as basin floor fans) are reservoirs for low salinity waters. This is the case in the Pliocene sands on the Mediterranean offshore Southern France and on the New Jersey inner shelf. However, due to their high permeability, they are also subject to density driven sea-water intrusion, as it is the case in the New Jersey middle shelf, that is, in a more distal area with limited connection to the onshore recharge.
In the case on the New Jersey inner shelf, the modern onshore hydraulic system is thought to extend offshore in both the coarse and fine grained fraction of the clinothems. Further offshore, the presence of salt water in the coarse grained fractions of the topsets has tentatively been explained by seawater penetration as a result of fast vertical density-driven flow [17]. In this case, the vertical flow would be expected to be greater than any horizontal, land-to-sea directed flow within the connected reservoirs. Diffusion is also sometimes accounted for the low salinity water hosted in mid-low permeability sediments including silts and clays. In summary, the emplacement and preservation of bodies of freshened groundwater is complex and depends on the timing of emplacement, 3D reservoir geometry, present-day connection to the sea and the presence of tight layers which may aid the preservation of freshwater aquifers.
As a general rule, the preservation of low-salinity pore waters is inversely proportional to the age/depth of sediments. Well defined stratigraphic boundaries such as unconformities, lateral and vertical facies changes and more specifically, maximum flooding and transgressive surfaces marking a shift to fine-grained facies or diagenetic boundaries, can seal the aquifer and constrain the upward extent of low-salinity aquifers. The role of low-permeability barriers (intra or extra-formational) to separate superimposed aquifers is shown in the Tanzania/Zanzibar and New Jersey case studies. This phenomenon has also been observed elsewhere, such as in Hong Kong [156][157][158].
Evaporites, which are considered a good seal for petroleum fluids, when in contact with undersaturated pore fluids, are prone to quick dissolution and produce salinization of offshore freshened groundwater [159,160]. This upwelling-driven salinization process is recorded for example, offshore New Jersey. Evaporites (in particular, halite) have a very specific seismic expression and can be clearly identified on seismic data.
In carbonate systems, poro-permeability prediction based on seismic and sequence stratigraphy is more complex and highly variable. However, seismic data can image the presence of karstic geomorphological features, of freshwater vents associated to karstic conduits and of subaerially exposed shelves which might have hosted a freshwater aquifer in the past.
Although we have highlighted the benefits of integrating seismic analysis into offshore freshened water research, the use of seismic data in such studies is constrained by a number of factors that might limit its application to this field of study. These factors include-(1) the high economic cost of seismic data acquisition and processing, (2) the constraints posed by environmental regulations that vary across country boundaries and (3) the technical challenges for data acquisition in specific areas such as in shallow waters, for example, the transition zone and proximal part of continental shelves (few to 10 s of meters of seabed depth).

Integration with Petroleum Geology Studies
The large seismic database used in oil exploration could potentially bridge the previously described 'data gap', due to its spatial coverage and vast geographic availability [161]. Seismic data interpretation in offshore petroleum studies is mostly focused on hydrocarbon-bearing basins, normally at a deeper burial depth compared to the described freshened groundwater systems. However, the case studies described in the Gippsland basin and in Tanzania [144,151,152] provide clear examples of the how the integration of petroleum exploration data with offshore hydrogeology analysis can lead to a better understanding of onshore-offshore groundwater systems. Similarly, in Suriname [162] industry data support the characterization of lithologies acting as a seal both for petroleum fluids and for freshwater diffusion.
More in general, industry seismic and drilling data can provide further information on offshore aquifers through data on fluid pressures, on presence of permeability barriers, on the 'hydrodynamic' or 'hydrostatic' nature of aquifers and on the presence of fossil low salinity waters [22,23,25,161,[163][164][165]. An example of deep low salinity water in a karst oil reservoir is recorded in the offshore drilling of the Rospo Mare structure (Adriatic Sea, about 30 km off the coast of Italy) at about 1300 m below sea-level, where heavy oil is trapped in Cretaceous karstic limestones. Most of these karst cavities are filled with fresh or slightly brackish water [18,166,167], which is thought to be directly connected with the Maiella limestone-rich mountains, where the aquifer is recharged, some 70 km to the NW of Rospo-Mare. Interestingly, this example is located in an active tectonic setting and differs from the classic passive margin case studies (e.g., New Jersey). Therefore, the role played by relative sea-level changes, such as in back-arc or foreland basins, must be considered in relation to how the tectonic tilting or drowning might create favorable conditions for the recharge from landmasses [59,61,62].
Although there is a significant amount of seismic data already acquired on continental margins worldwide, a series of factors make the data still relatively inaccessible or of limited use to the scientific community studying offshore groundwater. These are primarily-(1) the limits imposed by confidentiality and proprietary nature of most of industry-derived data; (2) the differences in acquisition design, normally focused on deeper targets; and related to this, (3) the large heterogeneity of data quality and resolution.

Applicability and Future Avenues of Research
In this review we have described a number of case studies focused on the application of seismic methods in offshore freshened groundwater analysis. The same concepts can be applied to other areas of continental margins worldwide, where extensive literature and a wealth of case studies exist and have already defined the seismic and sequence stratigraphic setting and their structural contexts. However, the focus of these studies (e.g., reservoir characterization, paleoclimate, sea-level, paleoenvironments, structural evolution of continental margins) has not usually included offshore groundwater. Comprehensive, large-scale mapping of the sedimentary units hosting potential aquifers, based on these data, can be used to improve calculations of the volumes of fresh/brackish waters entrapped in continental margins, building on previous attempts [4,16,93,102].
Further avenues in offshore groundwater research can be created by extracting additional information from seismic data and by expanding the analysis to different areas of marine basins. Extracting additional information from seismic data can involve for example, a wider use of attribute analysis including quantitative methods. Quantitative seismic interpretation can be used to map permeability and porosity variation across a study area and provide an enhanced implementation of EM models, for example, through seismic inversion analysis. Porosity and permeability are key when modelling EM data to map offshore groundwater salinity and so far they can only be determined from borehole data, which have a very limited spatial extent. Increased accessibility of 3D seismic and well calibration would facilitate the use of these seismic attributes for offshore freshened groundwater mapping.
Proximal and shallow aquifers of freshened water, generally connected to onshore systems, are currently the most commonly observed sites for the occurrences of offshore low-salinity groundwater. A frontier area for freshened groundwater mapping is represented by deep offshore reservoirs, as the few case studies using hybrid industry/academic data have suggested. According to the literature, there are a series of factors that can lead to the presence of low-salinity water in deep and distal submarine sediments, far beyond the expected current limits. These factors include-(a) recharge during past sea-level lowstand [95], (b) entrapment of connate water in rapidly subsiding basins [139], (c) release of water by the dissociation of gas hydrates and other diagenetic reactions [79], glacier-related injection ( [20] and references therein). Seismic data analysis can bring insights in these deeper aquifers, which could potentially host significant amount of freshened groundwater; however, they would be considered as paleo-groundwater resources and their exploitation regulated by economic, environmental, sustainability and recovery considerations.
Theoretically the methods reviewed here could be also applied to the analysis of freshened groundwater in shallower setting, such as the highly populated and urbanized land-sea "transition zone." However, mapping this 'white gap' (between coast and shallow offshore area) can prove challenging because acquisition of seismic data is limited by technical and socio-economic factors. Overcoming these challenges will make seismic data in the land-sea transition zone increasingly available and provide new tools for understanding the global distribution of offshore freshened groundwater close to coastal communities.