www.nature.

com/scientificreports

OPEN Aquifer systems extending far

offshore on the U.S. Atlantic margin
Chloe Gustafson   1, Kerry Key   1 & Rob L. Evans   2

Low-salinity submarine groundwater contained within continental shelves is a global phenomenon.

Received: 24 September 2018
Mechanisms for emplacing offshore groundwater include glacial processes that drove water into
Accepted: 20 May 2019 exposed continental shelves during sea-level low stands and active connections to onshore hydrologic
Published: xx xx xxxx systems. While low-salinity groundwater is thought to be abundant, its distribution and volume
worldwide is poorly understood due to the limited number of observations. Here we image laterally
continuous aquifers extending 90 km offshore New Jersey and Martha’s Vineyard, Massachusetts, on
the U.S. Atlantic margin using new shallow water electromagnetic geophysical methods. Our data
provide more continuous constraints on offshore groundwater than previous models and present
evidence for a connection between the modern onshore hydrologic system and offshore aquifers. We
identify clinoforms as a previously unknown structural control on the lateral extent of low-salinity
groundwater and potentially a control on where low-salinity water rises into the seafloor. Our data
suggest a continuous submarine aquifer system spans at least 350 km of the U.S. Atlantic coast and
contains about 2800 km3 of low-salinity groundwater. Our findings can be used to improve models of
past glacial, eustatic, tectonic, and geomorphic processes on continental shelves and provide insight
into shelf geochemistry, biogeochemical cycles, and the deep biosphere.

Offshore groundwater plays an important role in continental shelf systems, with implications for past and present
environments. Submarine groundwater discharge, the flow of water from continents to the ocean via onshore
outcrops and submarine channels, affects the biologic and chemical processes of the ocean, which in turn may
regulate nutrient fluxes important for global carbon cycling1. Sequestered submarine groundwater emplaced
in continental shelves during past sea-level low stands may provide information regarding past ocean salinity
changes2 and may also be used to constrain paleo-hydrologic models that give insight into past glacial extent and
sea-level changes3–6. Whether modern or ancient, low-salinity groundwater also represents a potential resource
in regions where onshore freshwater resources have diminished.
There is ample evidence of low-salinity submarine groundwater contained within the world’s continental
shelves, with direct or indirect observations occurring on all continents7. Several of these observations have taken
place on the U.S. Atlantic coast. During the 1970’s, the AMCOR drilling project revealed low-salinity groundwa-
ter up to 120 km offshore New Jersey. This dataset resulted in the development of a model of offshore low-salinity
groundwater characterized by a wedge-shaped feature, gradually increasing in salinity seaward and with depth8.
A more focused drilling survey completed by IODP Expedition 313 in 2010 showed more complex groundwater
geometry offshore New Jersey9, proving the earlier AMCOR-based model too simple. Further northeast, low
salinities seen in well data from onshore Martha’s Vineyard10 and Nantucket11 have also suggested low-salinity
submarine groundwater is present, although there has been no offshore drilling to confirm this hypothesis. While
borehole data provide direct samples, they are only spatial point measurements and do not provide a laterally
continuous characterization of the subsurface necessary for understanding offshore groundwater distributions.
Paleo-hydrologic modeling combines stratigraphy with past sea-level and ice-sheet history to better
understand offshore groundwater. Modeling studies have attempted to characterize the spatial distribution of
low-salinity groundwater offshore the U.S. Atlantic coast3–6, especially in regions with no direct offshore observa-
tions, such as Martha’s Vineyard. These models all provide different predictions of groundwater characterization
due to variations in model forcings and hydrological parameters, resulting in significant uncertainty in the nature
of the groundwater that may be present within the U.S. Atlantic margin.
Marine electromagnetic (EM) methods are geophysical methods capable of providing continuous measure-
ments required for characterizing offshore groundwater systems. EM methods are sensitive to bulk electrical

1
Lamont-Doherty Earth Observatory, Columbia University, Palisades, New York, USA. 2Department of Geology
and Geophysics, Woods Hole Oceanographic Institution, Woods Hole, Massachusetts, USA. Correspondence and
requests for materials should be addressed to C.G. (email: [email protected])

Scientific Reports | (2019) 9:8709 | https://doi.org/10.1038/s41598-019-44611-7 1

www.nature.com/scientificreports/ www.nature.com/scientificreports

Figure 1.  Electromagnetic survey of groundwater on the U.S. Atlantic continental shelf. Surface-towed CSEM
profiles (white dashed lines) and seafloor MT stations (white triangles) were collected during surveys offshore
New Jersey and Martha’s Vineyard. Yellow solid lines show the portions of the survey profiles where the EM data
support low-salinity pore water. Circles denoting onshore and offshore wells8–11,42–45 are colored yellow where
low pore fluid salinities (<15) were observed offshore or where aquifers have been identified via geophysical
data onshore. Magenta circles denote where only high salinity (>15) pore waters have been encountered.
Wells featured in the text are labeled with a larger font size. Yellow crosses and the yellow dashed line show our
inferred spatial extent of the low-salinity aquifer system. The white dotted line denotes the terminal moraine
of the Laurentide ice-sheet at the Last Glacial Maximum30. This map was created using GMT-5.4.3 (http://gmt.
soest.hawaii.edu/projects/gmt/wiki/Download). The topography source data for this figure is the Global Multi-
Resolution Topography Data (https://www.gmrt.org).

resistivity, which is primarily controlled by porosity and pore fluid salinity12 in offshore sediments. Many authors
have called for the utilization of EM methods to better understand offshore groundwater systems 5–7 and
active-source EM studies over small spatial scales have successfully identified offshore groundwater13,14. Our
study implemented technology originally produced for hydrocarbon exploration15 in order to establish the first
large-scale marine EM survey of low-salinity groundwater. We used a recently developed surface-towed con-
trolled source EM (CSEM) system16 in combination with the passive-source magnetotelluric (MT) method to
map low-salinity groundwater systems in two locations on the US Atlantic coast: offshore Martha’s Vineyard and
New Jersey (Fig. 1). CSEM and MT methods are preferentially sensitive to resistors and conductors, respectively
(Methods), so we use both data sets to obtain a more robust resistivity model.
Offshore New Jersey, our profiles were designed to better understand the spatial extent of low salinities already
measured by the nearby AMCOR wells (6011, 6020; Fig. 1) and the co-located IODP Exp. 313 wells (M0027,
M0028, M0029; Fig. 1), while offshore Martha’s Vineyard we collected a profile to test whether similar low-salinity
anomalies are present, as suggested by nearby onshore wells (6001 and ENW-50 (Fig. 1)). At both locations we
collected continuous surface-towed CSEM data along 130 km long profiles extending from the shoreline to the
shelf edge. Eight MT receivers were deployed 10–20 km apart along each main transmission profile to constrain
the broad-scale structure. Offshore New Jersey, two additional receivers were deployed 20 km northeast of the
seaward edge of the main transmission profile and we also collected CSEM data along seven additional intersect-
ing profiles (Fig. 1).
We carried out independent and joint inversions of the MT and CSEM data using a standard regularized 2-D
inversion code that yielded resistivity models for each profile (Methods and ref.17). We use our resistivity mod-
els to infer the distributions of low- and high- salinity water in the subsurface. Here we focus on the top 800 m of
the models where the resolution is highest and where we expect low-salinity water. For the purpose of discussion

Scientific Reports | (2019) 9:8709 | https://doi.org/10.1038/s41598-019-44611-7 2

Scientific Reports |
www.nature.com/scientificreports/

(2019) 9:8709 | https://doi.org/10.1038/s41598-019-44611-7

Figure 2.  Resistivity models obtained from jointly inverting surface-towed CSEM and seafloor MT data for the shore-to-shelf profiles off of Martha’s Vineyard and New Jersey with seismic and
well log data. a-d, 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. The left column shows the resistivity models for (a) Martha’s Vineyard and (c) New Jersey overlain with black lines showing resistive
zones interpreted as low-salinity aquifers and white lines showing the deeper conductive zones interpreted as saline or briny water. The deeper Jurassic basement seen in seismic reflection data
collected off Martha’s Vineyard is shown by a dark grey line. The right column displays co-located seismic reflection images on top of the (b) Martha’s Vineyard and (d) New Jersey resistivity
models as well as pore water salinity data8–11. 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. Solid pink and purple lines on salinity plots are the depths and corresponding salinity values used in Fig. 3. Seismically-imaged confining units and clinoform structures that
www.nature.com/scientificreports

influence groundwater distribution patterns are shown as labeled light grey dashed lines.

3
www.nature.com/scientificreports/ www.nature.com/scientificreports

SW

Fresh Brackish Saline Brine

10 2

Bulk Resistivity (ohm-m)

Resistive Zone
10 1 30
%
60 po
% ro
po sit
ro y
sit
y

10 0
Conductive Zone
6011
ENW-50 M0027 M0028

10 -1 10 0 10 1 10 2
Salinity
Figure 3.  Dependence of bulk resistivity on salinity. Bulk resistivity calculated as a function of pore water
salinity and porosity using Archie’s Law is shown as thick black lines. Salinity is shown using the unitless
Practical Salinity Scale 197818. Porosity values ranging from 30–60% are representative of the range encountered
in wells M0027–M00299,20. Shaded colors denote resistivity with the same color scale as Fig. 2. Grey horizontal
bars show resistive and conductive zones identified in Fig. 2. The combination of observed resistivity and
porosity indicates that there is low-salinity pore water in the resistive zone and saline or briny pore water in the
conductive zone. Water salinity classifications are shown on the top of the horizontal axis with seawater (SW)
specifically noted at the value 35. Select low- and high- salinity ranges from co-located wells are plotted as pink
and purple bars, respectively. The salinity ranges are limited to depths where the corresponding wells intersect
the resistive and conductive zones, shown in Fig. 2b,d. These ranges show good agreement with our imaged
resistivity for both the resistive and conductive zones.

within this text, we define low- and high- salinity water as having salinity less than and greater than 15, respec-
tively, where salinity is a unitless quantity as defined by the Practical Salinity Scale18.
Offshore Martha’s Vineyard, our model (Fig. 2a,b) contains two resistive zones and one conductive zone.
The resistive and conductive zones range from 10–70 ohm-m and 0.6–1 ohm-m, respectively. Resistive zone R1,
located within the top 400 m, gently dips seaward out to 90 km offshore and is approximately 200 m in thickness.
R1 is partially underlain by conductive zone C1, which spans from 40–120 km offshore and reaches depths of at
least 800 m. The deeper near-shore resistive zone, R2, is consistent with a Jurassic basement rock interpreted in a
co-located seismic image6.
Off New Jersey, our resistivity model (Fig. 2c,d) reveals two resistive zones within the upper 400 m that are
underlain by a conductive zone. The resistive and conductive zones range from 10–110 ohm-m and 0.4–1 ohm-m,
respectively. The near-shore resistive zone, R3, shallowly dips downward from the shore out to 60 km offshore,
reaching a maximum depth of about 400 m. The mid-shelf resistive zone, R4, is imaged in the top 50–200 m and
extends from 60 km to 90 km offshore, with no obvious dip. The underlying conductive zone C2, ranges from
450–750 m depth and spans from 60 km offshore to the end of the profile. Additional resistivity images for our
shoreline parallel profiles show that the resistive and conductive zones extend in the shore-parallel directions for
the entire lengths of the crossing profiles (Extended Data Fig. 8).
We infer pore water salinity by combining measured porosities with our resistivity models and employing
standard relationships between pore fluid resistivity, salinity and porosity12,18. These relationships demonstrate
that areas of high and low resistivity should correspond to areas of low- and high- salinity pore water, respectively
(Methods). Borehole sediment analyses from IODP Exp. 313 shows that low-salinity waters exist in sediments
with average porosities ranging from 30–60%19, so we use these values in our conversion from resistivity to salin-
ity. The resistivity values of R1, R3, and R4 suggest pore water salinity values range from about 0.2–9 (Fig. 3), in
good agreement with salinity ranges observed in intersecting wells9,20,21. The deeper conductive zones in both
profiles require salinities that exceed that of seawater (Fig. 3). Our imaged resistivity variations could also be
explained by porosity variations alone; however, if the pore fluid has seawater salinity, resistive zones R1, R3 and
R4 would have only 2–10% porosity, which is significantly lower than the borehole observations. Instead, our
preferred interpretation is that resistive zones R1, R3 and R4 represent low-salinity groundwater and conductive
zones represent saline or briny groundwater in the subsurface.
We interpret our resistive zones as submarine aquifers given observed or predicted permeabilities for
both New Jersey and Martha’s Vineyard. Previous work combining IODP Exp. 313 sediment and seismic

Scientific Reports | (2019) 9:8709 | https://doi.org/10.1038/s41598-019-44611-7 4

www.nature.com/scientificreports/ www.nature.com/scientificreports

onshore rech
arge
clinoforms

faults
basement rock upwelling brines

Figure 4.  Conceptual model of offshore groundwater. Offshore low-salinity aquifers are fed by the adjacent
onshore hydrologic system and are vertically constrained by a confining unit, until clinoform structures are
present. Low-salinity water is also present above the clinoforms and may originate from the larger low-salinity
aquifer. High-salinity groundwater exists deeper and further seaward and is caused by groundwater interaction
with underlying salt deposits. Arrows denote groundwater flow paths.

profile analyses offshore New Jersey predicts permeabilities ranging from 10–1000 md for offshore New Jersey19.
Permeability measurements from onshore Nantucket (6001) range from 10–780 md11. Although permeabilities
vary, all permeability values are high enough to allow the flow of water through sediments19,22, justifying our use
of the term “aquifer” to describe our resistive zones.
Direct pore water salinity measurements offshore New Jersey correspond well with our resistivity model
(Fig. 2d). The upper and lower boundaries of R3 show excellent agreement with low-salinity measurements from
wells 60118 and M00279,20,while high-salinity measurements from wells M0027, M0028, and M00299,20 are con-
sistent with the low resisivities of C2. The correlation between direct salinity measurements and resistivity is
not as strong for R4, however, wells 60208 and M00299,20 both show a sharp increase in pore water salinity at the
upper boundary of R4. Further seaward, well 60108 is also compatible with our resistivity model, as there are no
low-salinity measurements present where we do not image high resistivity.
No well data exists offshore Martha’s Vineyard, however, we can compare our offshore resistivity model with
nearby onshore well data from Martha’s Vineyard23 (ENW-50) and Nantucket21 (6001) (Fig. 2b). Well ENW-50
shows low salinities over its entire depth range, which extends from the surface down to 228 m. This correlates
well with the location of R1 and strongly suggests that R1 connects onshore. This region of low-salinity is also
seen in 6001, but is expressed at greater depths (below 320 m) due to the thinning of sediments from Natucket to
Martha’s Vineyard10,24.
Seismic reflection images6,25 co-rendered on our models (Fig. 2b,d) reveal consistent structural controls on
the low-salinity groundwater distribution. Continuous, well-defined stratigraphic boundaries mark the upward
extent of low-salinity groundwater and the onset of clinoform structures mark the seaward extent for R1 and R3.
While it is known that groundwater systems are typically constrained vertically by a confining unit, the corre-
spondence of the low-salinity groundwater seaward terminus with the onset of clinoforms for both R1 and R3
demonstrates a previously undetermined structural control on offshore groundwater distribution. Notably, both
confining units mark a change in geology; the top bound of R1 is a known unconformity26 and the top bound of
R3 marks the transition from continuous to discontinuous seismic reflectors, indicating marine and non-marine
sediment deposition, respectively.
Seismic reflection images combined with our resistivity models and well data strongly suggest that both R1
and R3 connect onshore. Resistors R1 and R3 have high-resistivity values (and therefore low-salinity) at the
shoreward extent of both profiles, suggesting a connection to the onshore system. Offshore New Jersey, the high-
est resistivity values occur at the most landward extent of our profile, which is consistent with lateral transport of
freshwater from onshore recharge. The seismically imaged aquifer confining units also continue to the shoreward
bounds of our images (Fig. 2b and ref.19), providing further evidence for an onshore connection. Previous work
combining IODP Exp. 313 sediment and seismic profile analyses indicates R3 exists within low-permeability
sediments that continue onshore New Jersey19, where there are several known aquifers27. Moreover, onshore
analysis of pore water sediments from the Atlantic City well show low-salinity water in low-permeability seaward
dipping sediments at depths corresponding to R328. Additionally, stable hydrogen and oxygen isotope data20 from
IODP Exp. 313 wells indicate the low-salinity water within R3 is a mix of modern seawater and modern terrestrial
meteoric water, which requires an onshore connection. Similar analyses have not been done for R1, as there are
no borehole measurements offshore Martha’s Vineyard. However, as previously mentioned, onshore boreholes
ENW-50 and 6001 provide direct evidence of onshore freshwater at depths corresponding with R1, making an
onshore connection highly plausible.

Scientific Reports | (2019) 9:8709 | https://doi.org/10.1038/s41598-019-44611-7 5

www.nature.com/scientificreports/ www.nature.com/scientificreports

Resistive zone R4 does not exhibit the same structural controls seen for R1 and R3. Above R3’s seismically
inferred confining unit and clinoforms, R4 exists within less-continuous seismic reflectors, indicating the sedi-
ments hosting this groundwater were deposited in a non-marine environment, whereas the continuous seismic
reflections in R3 indicate marine deposition. Also, unlike R3, the low-salinity water imaged by R4 exists primarily
in units with predicted high permeability19. The difference in permeability and depositional environment of the
host lithologies indicates the low-salinity water imaged by R4 may have a different origin than low-salinity water
in R3. Alternatively, R3 and R4 may be connected. Similar to R3, stable hydrogen and oxygen isotope data20
from IODP Exp. 313 wells indicate the low-salinity water within R4 are a mix of modern seawater and mod-
ern terrestrial meteoric water. This supports the idea that R3 and R4 are connected, or suggests the low-salinity
water imaged in R4 originates from a different onshore source, potentially connecting with known Long Island
aquifers27. A connection between R3 and R4 could indicate that there is a relationship between clinoforms and
upwelling groundwater, as R4 occurs above the seismically imaged clinoforms and we image high-resistivity val-
ues between R3 and R4. Low-salinity water may rise from R3 to R4 via unresolved faults, as low-salinity water is
gravitationally buoyant compared to denser saline water.
For both surveys, the conductive zones appear to have no corresponding relationship with seismically imaged
stratigraphic boundaries. However, underlying faults seismically imaged offshore New Jersey allow for the migra-
tion of buried salt deposits29, indicating the conductive zones in our results arise from groundwater interacting
with these deposits.
Our results demonstrate a massive offshore aquifer system, extending 90 km offshore, spans the New Jersey
coastline up to Martha’s Vineyard, and likely spans beyond the bounds of our survey. To the southwest of our New
Jersey profile, wells 6008 and 6009 (Fig. 1) have identified groundwater with salinities less than fifteen8, estab-
lishing a significant southward extension of the low-salinity groundwater. Between our two surveys, observed
onshore aquifers exist within seaward dipping sediments on Long Island that extend offshore. Previous hydro-
logic models show the extent of low-salinity groundwater terminates near the shoreline27; however, we assert
that some aquifer water must extend several 10’s of kilometers offshore (Fig. 1) given the extent of the Laurentide
ice-sheet during the Last Glacial Maximum30 (LGM) (Fig. 1) and its relationship to low-salinity groundwater that
we have observed offshore Martha’s Vineyard. Subglacial runoff beneath and in front of the ice-sheet during the
LGM is thought to have deposited freshwater around Martha’s Vineyard, as suggested by observed and modeled
overpressures3–6,21,31, and likely emplaced low-salinity groundwater from offshore Long Island to the northeast of
our Martha’s Vineyard survey. The yellow cross pattern in Fig. 1 shows our inferred low-salinity aquifer extent on
the U.S. Atlantic margin.
The modeled size, continuity, and onshore connectedness of the low-salinity groundwater advance ideas of the
offshore groundwater systems on the U.S. Atlantic coast. Previously, borehole data combined with seismic imag-
ing were used to infer a system of multiple, localized water reservoirs within the mid-shelf offshore New Jersey19,
whereas our models reveal a more regional view of widespread water distribution across the shelf. Offshore
Martha’s Vineyard there was greater uncertainty due to an absence of data sensitive to offshore groundwater prior
to our study; previous offshore models relied solely on paleo-hydrologic modeling, resulting in various hypoth-
esized systems, ranging from isolated freshwater packets6 to diffusive tongues of salinity contours3–5. Our study
shows laterally continuous low-salinity systems with well-defined boundaries in both survey locations, providing
the first look at the groundwater system offshore Martha’s Vineyard and redefining the New Jersey system. We
present a generalized conceptual model of offshore groundwater, including structural controls, in Fig. 4.
The spatial distribution of the low-salinity groundwater shown in our study provides important constraints
for understanding continental shelf processes. Our results can be used to improve paleo-hydrologic modeling,
as the interactions between stratigraphy, sea-level changes, and ice-sheet dynamics are all imprinted in sub-
marine groundwater systems. The sequestered low-salinity groundwater we image offshore Martha’s Vineyard
may be used to improve estimates of past ocean salinity changes2 as water stored underground and in ice-sheets
impacts ocean composition. Both of our profiles may also enhance understanding of sub-seafloor ecology and
microbiology, which are influenced by exchanges between low-salinity groundwater and the overlying ocean32,33.
Furthermore, onshore connectivity of the low-salinity groundwater may represent a significant, yet widely unac-
counted for, nutrient exchange mechanism between the terrestrial and ocean environments important for global
carbon cycling1.
The low-salinity groundwater we have imaged represents a significant resource with potential for future off-
shore groundwater production. Given the 90 km offshore extent, 350 km distance between profiles, about 200 m
thickness and 45% average porosity, we estimate the aquifer system contains at least 2800 km3 of water with
salinity less than 15. If we consider the potential northeast and southwest extensions beyond our profiles, there
may be several times more groundwater underlying the northeast portion of the U.S. Atlantic continental shelf,
representing a freshwater resource that rivals the largest onshore aquifers. Utilization of such groundwater would
require desalinization and effects such as upconing should be considered prior to extraction.
We have shown continuous aquifers extend from onshore New Jersey and Martha’s Vineyard out to 90 km
offshore. Our data show good agreement with onshore and offshore borehole data, demonstrating marine elec-
tromagnetic data are useful for mapping offshore groundwater. Correlating our resistivity images with seismic
data shows aquifers connecting onshore are vertically constrained by confining units and laterally constrained
by the onset of clinoforms. Our results improve the overall knowledge of groundwater offshore the U.S. Atlantic
coast and may be used in hydrologic and oceanographic models of the region in order to better understand past,
present, and future continental shelf dynamics.

Scientific Reports | (2019) 9:8709 | https://doi.org/10.1038/s41598-019-44611-7 6

www.nature.com/scientificreports/ www.nature.com/scientificreports

Methods
The marine MT method measures the horizontal components of naturally-occurring low frequency EM fields
on the seafloor to produce an impedance response sensitive to electrical resistivity from the seafloor to the deep
crust and mantle, depending on the frequency range. The CSEM method uses a towed horizontal electrical dipole
source to broadcast EM fields to one or more receivers that measure the attenuation and phase shift of the trans-
mitted fields to generate an Earth response sensitive to the seafloor resistivity at depths of meters to a few kilom-
eters, depending on the source-receiver offset and frequency band.
During the ten day cruise, we collected MT stations at each survey location by deploying ten broadband sea-
floor EM receivers34 for up to 4 days. Surface-towed CSEM data were collected by towing a 336.9 m long dipole
antenna behind the ship and broadcasting a 88 A current using a doubly symmetric square waveform35 with a
0.25 Hz fundamental frequency. This waveform has higher signal-to-noise ratios for higher frequencies compared
to a square wave and other common waveforms. The dipole moment of the antenna - the product of the antenna
length and its transmitted current - is 2.96 × 104 Am. This large dipole moment combined with short receiver
offsets ensures the CSEM signal is stronger than MT signal34; MT signal manifests as lower power random noise
in the measured CSEM signal. Four broadband EM receivers were surface-towed behind the vessel at offsets 600,
870, 1120, and 1380 m. Each surface-towed receiver measured the inline horizontal electric field on a 2 m dipole
positioned about 0.67 m below the sea surface16. The receiver orientations and positions were measured using
recording electronic compasses and GPS logging systems. All CSEM tows off New Jersey were completed success-
fully. Off Martha’s Vineyard, the shortest offset towed receiver was fouled by fishing gear and stopped producing
reliable data around 10 km into the profile. At about 70 km along the profile (Fig. 1), the towed array became
further damaged when snagging fishing gear and a nearby buoy line and so the surface towed data collection was
terminated. Despite incurred damages, we collected 70 km of high quality CSEM data.
Extended Data Figs 1–3 show the inversion models and data from Martha’s Vineyard and Extended
Data Figs 4–6 show the inversion models and data from New Jersey. CSEM response functions for the four
surface-towed EM receivers were estimated with a robust stacking method35. We used 240 s stacking windows
corresponding to 250–300 m lateral distance between transmitter stack points. This yielded high quality ampli-
tude and phase responses for each receiver as a function of position and frequency harmonic. While the trans-
mitter waveform yielded about 8 usable frequency harmonics, here we only consider data at the two strongest
low-frequency harmonics, 0.75 and 1.75 Hz, as they contain the highest sensitivity to the depth range of interest
(Extended Data Figs 3 and 6). The transmitted CSEM signal is well above the noise floor (Extended Data Fig. 7).
The CSEM data from both survey regions exhibit high amplitudes and low phase lags on the inner to mid-shelf
sections of the profiles, which is consistent with low attenuation associated with high resistivity seafloor in these
regions. Conversely, the amplitudes significantly decrease and the phase lags increase on the outer shelf sections
of the profiles, indicating the seafloor is more conductive in these regions.
MT impedance tensors were obtained by processing two to four day recordings of the seafloor horizontal
electric and magnetic fields at each station using a robust multivariate array processing method36. This resulted
in generally high-quality MT responses at periods from 0.01 to 1000 s. Broadband seafloor MT data collected
in deep water (around 1 km) can typically only be obtained at periods of 1 s and longer due to the attenuation
of the shorter periods by the conductive seawater37, but here the water depths (20–100 m) were too shallow to
significantly attenuate the shorter period fields at 0.01 to 1 s. Noisy responses due to ocean waves in the 2–10 s
period band were omitted from the data, depending on the station. Further, data at periods longer than 10 s were
omitted from the analysis here because the depth sensitivity at periods of 10 to 1000 s greatly exceeds the depth
range of interest for studying the shallow aquifers. Analysis of the Martha’s Vineyard MT profile for deep crustal
and mantle structure using the full data bandwidth is presented in ref.38.
The resulting MT responses in the 0.01–10 s band have nearly equal off-diagonal impedance components,
suggesting the shallow structure is one-dimensional within the range of sensitivity of each instrument (Extended
Data Figs 2 and 5). Off Martha’s Vineyard, the MT responses at nearshore stations M01–M06 have characteristi-
cally high apparent resistivities that exceed 1 ohm-m and display a positive peak at periods less than 1 s, indicat-
ing a shallow resistive feature is present beneath these stations. The further offshore stations (M07–M10) have
apparent resistivities at periods below 1 s that are closer to 1 ohm-m and without the strong peak, suggesting the
absence of significant resistive structures beneath these stations. Off New Jersey, similar high apparent resistivities
below 1 s period are seen for stations N01–N05, suggesting shallow resistive structure is present beneath these
stations, while the outer-shelf stations (N06–N08) have much lower apparent resistivities.
Two-dimensional (2-D) resistivity models were obtained by inverting surface-towed CSEM and MT data
jointly and independently using MARE2DEM17, a freely-available goal-oriented adaptive finite-element inver-
sion code for 2-D electromagnetic modeling and inversion. MARE2DEM uses a standard non-linear regularized
inversion approach to find the smoothest model fitting the data. We manually removed clear outliers in both the
MT and surface-towed CSEM data as a final processing step before inverting the data. MT data uncertainties were
subject to a 10% error floor while CSEM data were subject to a 1% error floor. We modeled the 336.9 m long trans-
mitter antenna as a finite-length bipole by integrating point dipole sources along the length of the antenna using
an efficient Gauss quadrature approach. The 2 m long towed CSEM receiver dipoles and 10 m long seafloor MT
receiver dipoles were modeled using point dipole approximations. The main New Jersey profile consists of a total
of 6830 CSEM data and 360 MT data. The Martha’s Vineyard profile has a total of 3260 CSEM data and 336 MT
data. We constructed inversion parameter grids consisting of 31,004 cells for the New Jersey profile and 33,739
cells for the Martha’s Vineyard profile. Due to the extreme width to depth ratio in the inversion models and the
~20 km wide MT station spacing, both shore-to-shelf profiles were subject to a horizontal to vertical roughness
penalty ratio of 300 while the crossing profiles were subject to a ratio of 30. All resistivity inversions fit the joint
data set to a root-mean-square (RMS) misfit of 1.0.

Scientific Reports | (2019) 9:8709 | https://doi.org/10.1038/s41598-019-44611-7 7

www.nature.com/scientificreports/ www.nature.com/scientificreports

Separate and joint inversions of the CSEM and MT data sets show the relative sensitivities to the seafloor
structure (Extended Data Figs 1 and 4). The MT-only inversions show subtle resistive features in the upper 400 m
from the start of the profiles to about 90 km offshore; however, the resistivity only reaches about 10 ohm-m. This is
not surprising as MT data are more sensitive to thin conductive layers than thin resistive layers and MT inversions
are well-known to underestimate the resistivity of thin resistive layers39. CSEM data are much more sensitive to
resistive layers and thus the CSEM-only inversions find much higher resistivities. Further, they have much better
lateral resolution due to the relatively short spacing between data points. However, given the limited number of
towed receivers, their relatively short offsets from the transmitter, and the diffusive nature of low frequency EM
fields, the resistive features are quite smoothly resolved in depth. For both survey regions, the joint inversions of
both CSEM and MT data are significantly improved over either of the individual data set inversions, showing the
highest lateral and depth resolution and also the strongest resistivity contrasts. Thus, our preferred interpretations
are based on the joint inversion models.
We must also consider the possibility that the seafloor conductivity is anisotropic. While MT data are primar-
ily sensitive to horizontal resistivity, the vertical loops of current created by the CSEM transmitter make the towed
receiver CSEM data sensitive to both horizontal and vertical resistivity anisotropy. However, all of our independ-
ent and joint inversions for isotropic resistivity fit the respective data sets well, so from a data fitting point of view
anisotropy is not required. Further, we do not see any anisotropic artifacts such as alternating horizontal bands
of resistors and conductors that are characteristic of isotropic inversions of strongly anisotropic data40. Therefore,
we conclude that isotropic modeling is sufficient for these data.
The distal edge of the resistive zone 90 km offshore Martha’s Vineyard is in close proximity to the seaward
limit of the towed CSEM data and so it is worthwhile to examine whether the data require the resistor to end at
this location or if the resistor ends because of the truncation of CSEM data. We tested this by running forward
models of our preferred inversion model, but with the resistor extending all the way to the shelf edge, representing
about 4.5 km2 or 9 km2 of additional cross sectional area, so that its impact on the MT responses on the outer shelf
could be tested (Extended Data Fig. 9). We tested a range of resistivity values (5–100 ohm-m) for these exten-
sions and tabulated the resulting RMS misfit increases for the three MT stations without overlapping CSEM data
(M07-M10) in Extended Data Table 1. All resistivity thickness combinations result in an RMS increase greater
than our error floor of 10%, indicating it is unlikely fresh to brackish water exists further beyond our interpreta-
tion (Fig. 2a).
For our comparison of salinity and bulk resistivity we used a two-stage process. First bulk resistivity was
converted into pore fluid resistivity using Archie’s Law, which relates bulk resistivity and porosity to pore fluid
resistivity12. We use a cementation exponent of m = 2 for our Archie’s Law calculation, which is a typical value
for seafloor sediments comprised of sands and clays41. Porosity was assumed to be 30–60% based on the range
encountered in the boreholes off New Jersey9. The pore fluid resistivity was then converted into salinity using the
Practical Salinity Scale 197818, which is commonly used to convert standard oceanographic seawater conductivity
(the inverse of resistivity) measurements into a unitless salinity value. We use a temperature of 10 °C in converting
the pore fluid resistivity to salinity, which is reasonable based on temperature measurements from sea-bottom
expendable bathy thermographs (XBTs) during data collection. If we use m = 1.5 or m = 3 for our Archie’s calcu-
lations, the salinity values corresponding to our resistive zones slightly decrease and increase, respectively (see
Extended Data Table 2), but do not change our overall interpretation of submarine groundwater distributions.

Data Availability
Upon publication, all data that were analyzed in this study will be made available at https://emlab.ldeo.columbia.
edu/index.php/projects/freshwater/.

Code Availability
The MARE2DEM code used to invert the data is available at http://mare2dem.ucsd.edu.

References
1. Moore, W. S. The effect of submarine groundwater discharge on the ocean. Annu. Rev. Mar. Sci 2, 59–88 (2010).
2. Adkins, J. F., McIntyre, K. & Schrag, D. P. The salinity, temperature, and d18O of the glacial deep ocean. Science 298, 1769–1773
(2002).
3. Person, M. et al. Pleistocene hydrogeology of the Atlantic continental shelf, New England. Geological Society of America Bulletin 115,
1324–1343 (2003).
4. Marksamer, A. J. et al. In Subsurface Hydrology: Data Integration for Properties and Processes (eds Hyndman, D. W., Day-Lewis, F. D.
& Singha, K.) 171, 143–159 (American Geophysical Union, 2007).
5. Cohen, D. et al. Origin and Extent of Fresh Paleowaters on the Atlantic Continental Shelf, USA. Ground Water 48, 143–158 (2010).
6. Siegel, J. et al. Influence of late Pleistocene glaciations on the hydrogeology of the continental shelf offshore Massachusetts, USA.
Geochem. Geophys. Geosyst. 15, 4651–4670 (2014).
7. Post, V. E. A. et al. Offshore fresh groundwater reserves as a global phenomenon. Nature 504, 71–78 (2013).
8. Hathaway, J. C. et al. Preliminary summary of the 1976 Atlantic margin coring project of the US Geological Survey. 76–844 (U.S.
Geological Survey: U.S. Geological Survey Open-File Report, 1976).
9. Mountain, G. S. et al. Proceedings of the Integrated Ocean Drilling Program, Volume 313. Integrated Ocean Drilling Program
Management International, Inc., Tokyo, https://doi.org/10.2204/iodp.proc.313.101.2010 (2010).
10. Hall, R. E., Poppe, L. J. & Ferrebee, W. M. A Stratigraphic Test Well, Martha’s Vineyard, Massachusetts. (Washington D.C., U.S.
Geological Survey, 1980).
11. Folger, D. W., Hathaway, J. C., Christopher, R. A., Valentine, P. C. & Poag, C. W. Stratigraphic Test Well, Nantucket Island,
Massachusetts. (Arlington, VA, U.S. Geological Survey, 1978).
12. Hoefel, F. G. & Evans, R. L. Impact of Low Salinity Porewater on Seafloor Electromagnetic Data: A Means of Detecting Submarine
Groundwater Discharge? Estuarine, Coastal and Shelf Science 52, 179–189 (2001).
13. Evans, R. L. Using CSEM techniques to map the shallow section of seafloor: From the coastline to the edges of the continental slope.
Geophysics 72, WA105–WA116 (2007).

Scientific Reports | (2019) 9:8709 | https://doi.org/10.1038/s41598-019-44611-7 8

www.nature.com/scientificreports/ www.nature.com/scientificreports

14. Haroon, A., Lippert, K., Mogilatov, V. & Tezkan, B. First application of the marine differential electric dipole for groundwater
investigations: A case study from Bat Yam, Israel. Geophysics 83, B59–B76 (2018).
15. Constable, S. Review paper: Instrumentation for marine magnetotelluric and controlled source electromagnetic sounding.
Geophysical Prospecting 61, 505–532 (2013).
16. Sherman, D., Kannberg, P. & Constable, S. Surface towed electromagnetic system for mapping of subsea Arctic permafrost. Earth
and Planetary Science Letters 460, 97–104 (2017).
17. Key, K. MARE2DEM: a 2-D inversion code for controlled-source electromagnetic and magnetotelluric data. Geophysical Journal
International 207, 571–588 (2016).
18. Lewis, E. L. & Perkin, R. G. Salinity: Its definition and calculation. Journal of Geophysical Research Solid Earth 83, 466–478 (1978).
19. Lofi, J. et al. Fresh-water and salt-water distribution in passive margin sediments: Insights from Integrated Ocean Drilling Program
Expedition 313 on the New Jersey Margin. Geosphere 9, 1009–1024 (2013).
20. van Geldern, R. et al. Stable isotope geochemistry of pore waters and marine sediments from the New Jersey shelf: Methane
formation and fluid origin. Geosphere 9, 96–112 (2013).
21. Kohout, F. A. et al. Fresh ground water stored in aquifers under the continental shelf: implications from a deep test, Nantucket
Island, Massachusetts. Water Resources Bulletin 13, 373–386 (1977).
22. Cossé, R. Oil and gas field techniques. Basics of reservoir engineering. (Paris, Éditions Technip, 1993).
23. Delaney, D. F. Ground-water hydrology of Martha’s Vineyard, Massachusetts. Hydrologic Atlas 618 Plate 1, https://doi.org/10.3133/
ha618 (Arlington VA, U.S. Geological Survey,1980).
24. Kohout, F. A. et al. In The Geology of North America (eds Sheridan, R. E. & Grow, J. A.) 2, 463–480 (Boulder, CO, 1988).
25. Mountain, G. & Proust, J.-N. & Expedition 313 Science Party. The New Jersey Margin Scientific Drilling Project (IODP Expedition
313): Untangling the Record of Global and Local Sea-Level Changes. Scientific Drilling 10, 26–34 (2010).
26. Siegel, J. et al. Geophysical evidence of a late Pleistocene glaciation and paleo-ice stream on the Atlantic Continental Shelf offshore
Massachusetts, USA. Marine Geology 303, 63–74 (2012).
27. Masterson, J. P. et al. Hydrogeology and Hydrologic Conditions of the Northern Atlantic Coastal Plain Aquifer System From Long Island,
New York, to North Carolina. (Reston, VA, U.S. Geological Survey, 2015).
28. Pucci, A. A., Szabo, Z. & Owens, J. P. Variations in pore-water quality, mineralogy, and sedimentary texture of clay-silts in the lower
Miocene Kirkwood formation, Atlantic City, New Jersey. 150X, 317–341 (College Station, TX, Ocean Drilling Program, 1997).
29. Manheim, F. T. & Hall, R. E. Deep evaporitic strata off New-York and New-Jersey-Evidence from interstitial water chemistry of drill
cores. Journal of Research of the US Geologic Survey 4, 697–702 (1976).
30. Uchupi, E., Driscoll, N., Ballard, R. D. & Bolmer, S. T. Drainage of late Wisconsin glacial lakes and the morphology and late
quaternary stratigraphy of the New Jersey-southern New England continental shelf and slope. Marine Geology 172, 117–145 (2001).
31. Siegel, J., Lizarralde, D., Dugan, B. & Person, M. A. Glacially generated overpressure on the New England continental shelf:
Integration of full-waveform inversion and overpressure modeling. Journal of Geophysical Research Solid Earth 119, 3393–3409
(2014).
32. Berndt, C. Focused fluid flow in passive continental margins. Philosophical Transactions of the Royal Society A: Mathematical,
Physical and Engineering Sciences 363, 2855–2871 (2005).
33. Schippers, A. et al. Prokaryotic cells of the deep sub-seafloor biosphere identified as living bacteria. Nature 433, 858–861 (2005).
34. Constable, S. Review paper: Instrumentation for marine magnetotelluric and controlled source electromagnetic sounding.
Geophysical Prospecting 61, 505–532 (2013).
35. Myer, D., Constable, S. & Key, K. Broad-band waveforms and robust processing for marine CSEM surveys. Geophysical Journal
International 184, 689–698 (2011).
36. Egbert, G. D. Robust multiple-station magnetotelluric data processing. Geophysics Journal International 130, 475–496 (1997).
37. Key, K. Marine Electromagnetic Studies of Seafloor Resources and Tectonics. Surveys in Geophysics 33, 135–167 (2011).
38. Attias, E., Evans, R. L., Naif, S., Elsenbeck, J. & Key, K. Conductivity structure of the lithosphere-asthenosphere boundary beneath
the eastern North American margin. Geochem. Geophys. Geosyst. 18, 676–696 (2017).
39. Key, K. W., Constable, S. C. & Weiss, C. J. Mapping 3D salt using the 2D marine magnetotelluric method: Case study from Gemini
Prospect, Gulf of Mexico. Geophysics 71, B17–B27 (2006).
40. Du, Z. & Key, K. Case study: North Sea heavy oil reservoir characterization from integrated analysis of towed-streamer EM and
dual-sensor seismic data. The Leading Edge 77, 608–615 (2018).
41. Revil, A. & Cathles, L. M. III. Permeability of shaly sands. Water Resources Research 35, 651–662 (1999).
42. Proceedings of the Ocean Drilling Program, Scientific Results. 150X, (College Station, TX, Ocean Drilling Program, 1997).
43. Austin, J. A., Jr, Christie-Blick, N. & Malone, M. J. Proceedings of the Ocean Drilling Program, Initial Reports. 174A, (College Station,
TX, Ocean Drilling Program, 1998).
44. Miller, K. G. et al. Proceedings of the Ocean Drilling Program, Initial Reports. 174AX, (College Station, TX, Ocean Drilling Program,
1998).
45. Mountain, G. S., Miller, K. G. & Blum, P. Proceedings of the Ocean Drilling Program, Initial Reports. 150, (College Station, TX, Ocean
Drilling Program, 1994).

Acknowledgements
This project was supported by National Science Foundation awards 1458392 and 1459035. We thank the captain
and crew of the R.V. Marcus G. Langseth for a successful cruise and the Marine EM Lab at Scripps Institution of
Oceanography for providing the instrumentation. We also thank Chris Armerding, Marah Dahn, John DeSanto,
Jimmy Elsenbeck, Matt Folsom, Keiichi Ishizu, Jeff Pepin, Charlotte Wiman and Georgie Zelenak for participating
in the cruise. We thank Greg Mountain for providing depth migrated seismic imaging of the New Jersey shelf
and Brandon Dugan for providing depth migrated seismic imaging of the Martha’s Vineyard shelf. We thank
Michael Steckler and Robin Bell for helpful conversations.

Author Contributions
K.K. and R.E. designed the experiment and collected the data. C.G. and K.K. processed and modeled the data. All
authors contributed to writing the manuscript.

Additional Information
Supplementary information accompanies this paper at https://doi.org/10.1038/s41598-019-44611-7.
Competing Interests: The authors declare no competing interests.
Publisher’s note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and
institutional affiliations.

Scientific Reports | (2019) 9:8709 | https://doi.org/10.1038/s41598-019-44611-7 9

www.nature.com/scientificreports/ www.nature.com/scientificreports

Open Access This article is licensed under a Creative Commons Attribution 4.0 International
License, which permits use, sharing, adaptation, distribution and reproduction in any medium or
format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Cre-
ative Commons license, and indicate if changes were made. The images or other third party material in this
article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the
material. If material is not included in the article’s Creative Commons license and your intended use is not per-
mitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the
copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2019

Scientific Reports | (2019) 9:8709 | https://doi.org/10.1038/s41598-019-44611-7 10