Marine electromagnetic (EM) signals largely depend on three factors: flow velocity, Earth's main magnetic field, and seawater's electrical conductivity (which depends on the local temperature and salinity). Because of this, there has been recent interest in using marine EM signals to monitor and study ocean circulation. Our study utilizes voltage data from retired seafloor telecommunication cables in the Pacific Ocean to examine whether such cables could be used to monitor circulation velocity or transport on large oceanic scales. We process the cable data to isolate the seasonal and monthly variations and then evaluate the correlation between the processed data and numerical predictions of the electric field induced by an estimate of ocean circulation. We find that the correlation between cable voltage data and numerical predictions strongly depends on both the strength and coherence of the model velocities flowing across the cable, the local EM environment, as well as the length of the cable. The cable within the Kuroshio Current had good correlation between data and predictions, whereas two of the cables in the Eastern Pacific Gyre – a region with both low flow speeds and interfering velocity directions across the cable – did not have any clear correlation between data and predictions. Meanwhile, a third cable also located in the Eastern Pacific Gyre showed good correlation between data and predictions – although the cable is very long and the speeds were low, it was located in a region of coherent flow velocity across the cable. While much improvement is needed before utilizing seafloor voltage cables to study and monitor oceanic circulation across wide regions, we believe that with additional work, the answer to the question of whether or not seafloor voltage cables can be used to study large-scale circulation may eventually be yes.
Introduction
Evaluating and predicting the ocean state is crucially important for reconciling and mitigating the impact of climate change on our planet. Oceanic electromagnetic (EM) signals may be directly related to physical parameters of the ocean state, including flow velocity, temperature, and salinity. This has been known for centuries: in 1832, Michael Faraday was the first to attempt an experiment of measuring the voltage induced by the brackish water of the Thames River . His study was inconclusive, but since then, marine EM signals have been detected by both ground and satellite measurements .
Marine electromagnetic fields are produced because saline ocean water is a conducting fluid with a mean electrical conductivity of σ=3–4 S m-1. As this electrically conductive fluid passes through Earth's main magnetic field (|F|≈20–70 µT), it induces electric fields, electric currents, and secondary magnetic fields. The electric current produced by a specific oceanic flow depends on the flow's velocity, the Earth's main magnetic field, and the seawater electrical conductivity, which in turn depends on salinity and temperature. Thus, ideally, three physical oceanic parameters could be extracted from marine EM studies: velocity, salinity, and temperature. However, extracting multiple parameters would require using multiple oceanic electromagnetic signals (e.g., the signals from multiple tidal modes and perhaps also from circulation) .
In practice, velocity is the only quantity so far determinable from marine EM data. This was accomplished using a passive seafloor telecommunications cable which recorded the voltage difference between Florida and Grand Bahama Island, a distance of approximately 200 km . As the Florida Current passed over the cable, a voltage was induced, and this voltage was directly related to the depth-integrated velocity across the cable (i.e., they determined the transport volume). Since 1985, the National Oceanic and Atmospheric Administration (NOAA) has been using submarine cables to monitor the transport of the Florida Current through the Straits of Florida .
While data from seafloor voltage cables have been used to study a variety of geopotential fields , NOAA's work in the Straits of Florida is the only case of a seafloor voltage cable being reliable to determine the overlying oceanic flow. Numerical work suggests that cables spanning larger regions should still strongly correlate with the flow velocities ; however, there are many challenges in using longer cables. These challenges are largely due to the myriad of processes which may also induce marine electromagnetic fields, especially across the length of the cable, such as secular variation , variations in ionospheric tides , geomagnetic storms, or longer period ionospheric and magnetospheric signals . Additionally, because the cable voltage is produced from the electric field integrated along the entire cable length, the longer the cable is, the more challenging it is to decompose the total contribution to the cross-cable ocean transport in any particular section of the cable.
This study aims to provide a “first step” answer to the following question: can seafloor voltage cables be used to study large-scale circulation? To investigate whether it may eventually be feasible to use large-scale voltage cables for monitoring ocean flows, we evaluate the correlation between data from large-scale seafloor voltage cables and numerical predictions of the electric field induced by 3-D ocean circulation velocity fields. While this work builds on studies using seafloor voltage cables to monitor flow velocity in ∼ 100 km wide passages, this study aims to examine this application in basin-wide seafloor voltage cables.
The seafloor voltage cables used in this study. The HAW1N and HAW1S cables run parallel to each other.
CableStarting locationEnding locationLength (km)Time spanHAW1N and HAW1SPoint Arena, CA, USAHanauma Bay, HI, USA3805Apr 1990–Dec 2001HAW3San Luis Obispo, CA, USAMakaha, HI, USA3946Aug 1994–Jul 2000OKINinomiya, Honshu, JapanOkinawa, Japan1447Apr 1999–Dec 2001Data and data processing
This study used hourly data from four seafloor voltage cables (detailed in Table ): three retired AT&T cables (the HAW cables) and one cable managed by the University of Tokyo's Earthquake Research Institute (the OKI cable). The HAW1N and HAW1S cables are 3805 km long and run parallel to each other from Point Arena, California, to Hanauma Bay, Hawaii. As shown in Fig. , the parallel cables have very similar data, providing a unique and helpful situation for testing the data processing methods and for comparing the observations to numerical predictions. These three cables have been used in previous studies, including those examining geopotential variations , ionospheric phenomena , oceanic tides , and electrical conductivity of the lithosphere and the mantle .
The voltage data of HAW1N versus HAW1S are shown in a correlation scatterplot. As shown by the line of best fit (y=x-0.002), the data from the two cables match very closely.
The first step in processing the hourly data was the removal of geomagnetically noisy days (i.e., days where the geomagnetic Ap index was greater than or equal to 20; see for more on the Ap index). In this way, we reduce the contribution from magnetic field variations of magnetospheric
origin and their induced counterparts. This shrunk the amount of available data by 16.1 %–21.6 % for each cable. Further reduction of the datasets by using only night-side data is impossible for the HAW
cables, which span multiple time zones, and impractical for the OKI cable due to a significant decrease in the
dataset size and increase in variance. Next, to remove tidal signals, the 12 dominant daily tidal modes were fit to the data via least squares and then subtracted.
The following tidal periods were used: 4 (S6), 4.8 (S5), 6 (S4), 8 (S3), 11.967236 (K2), 12 (S2), 12.421 (M2), 12.6583 (N2), 23.934472 (K1), 24 (S1), 24.066 (P1), and 25.891 h (O1). Because the datasets have many gaps exceeding 24 h in length (for example, see Fig. ), bandpass filtering was not used. The data were then smoothed using cubic splines. For seasonal variations, we used 90 d knots between splines, and for monthly variations, we used 30.5 d between knots. Although the daily variations should directly relate to barotropic wind-forced processes , because of both the data's hourly time sampling and frequent data gaps, as well as challenges in producing daily numerical predictions, we chose to focus on monthly and seasonal variations. Each step of the data processing is shown in Fig. . As
the final step, the mean value is removed from all time series.
Each step of the data processing is shown here using HAW3 as an example: (a) the raw time series, (b) the time series with days of Ap >20 removed and tidal signals also removed, and (c) the smoothed time series produced by splines with 90 d knots.
A weakness of this data processing is that it does not prevent the inclusion of induced signals due to seasonal changes in ionospheric electromagnetic tidal strength. While we removed tidal signals from a least-squares fit, we applied this fit to the entire dataset and did not attempt to remove seasonal changes in ionospheric tides. Seasonally, ionospheric tides can significantly change amplitude , and the horizontal components of these tides are likely to induce signals at the ground , however, attempting to constrain seasonal changes in tidal strength is challenging. Ideally, the least-squares fit could be conducted on shorter intervals of the data, but this worsens the accuracy of the least-squares inversion. Ionospheric field models could be used, but this would also introduce unknown error quantities. Thus, we did not attempt to remove seasonal changes in tidal amplitude but remind the reader that these signals may influence the monthly and seasonal variations.
The contribution of the main field secular variation is not removed from the data, as it is included in the numerical calculations described in the next section.
Numerical predictions of the ocean circulation's electric field
We numerically predict the electromagnetic signals produced by ocean circulation using the ElmgTD time-domain numerical solver of the electromagnetic induction equation :
μ0∂B∂t+∇×1σ∇×B=μ0∇×u×F.
Here, B(r;t) is the induced magnetic field, u(r;t) is the velocity, μ0 is the magnetic permeability of vacuum, σ(r;t) is the electrical conductivity, and F(r;t) is the main geomagnetic field. The observable electric field E(r;t) is obtained from the induced magnetic field by post-processing,
E=1μ0σ(∇×B)-u×F.
The ElmgTD time-domain solver is based on spherical harmonic parameterization in lateral coordinates and uses 1-D finite elements for radial discretization. The model is fully 3-D, also incorporating the vertical stratification of the ocean electrical conductivity and of the velocities as well as accounting for the effect of variable bathymetry. Moreover, the seasonal variations in the ocean electrical conductivity and the secular variations in the main field are taken into account. The solution includes both the poloidal and toroidal components of the induced magnetic field , thereby allowing for the inductive and galvanic coupling between the oceans and the mantle as well as self-induction within the oceans. Numerically, the linear system is solved by the preconditioned iterative BiCGStab(2) scheme with massive parallelization applied across the time levels.
Monthly values of the horizontal and vertical components of ocean velocity from the data-assimilated Estimating the Circulation and Climate of the Ocean (ECCOv4r4) model were input into the ElmgTD solver to compute the electromagnetic fields that they induce from January 1997 to November 2001. Along with the monthly velocity values from ECCO, monthly values from the International Geomagnetic Reference Field (IGRF) were used for the main field, and monthly climatological data from NOAA's World Ocean Atlas (WOA) were used to describe the global seawater electrical conductivity σ. The conductivity model also includes the coastal and ocean sediments
on the seafloor with thickness distribution and conductivity values following .
Figure illustrates these inputs used for the ElmgTD numerical solver. The vertical
velocity is not shown here; although it is included in our calculations, as it represents
only a minimum additional computational burden, its effect on the induced fields is negligible.
Underlying these inputs, the electrical conductivity of the mantle follows the 1-D global profile obtained by the inversion of satellite data .
The surface velocities from ECCO are shown (a) for the zonal (U) component and (b) for the meridional (V) component. The labeled, thick black lines denote the seafloor voltage cables used in this study. A snapshot of the IGRF vertical main field, Fz, from 17 January 1997 is illustrated in panel (c), and panel (d) depicts the NOAA World Ocean Atlas seawater electrical conductivity's January climatology in the surface layer.
In the present calculations, we truncate the spherical harmonic expansion at degree 240, corresponding to approximately 0.75∘× 0.75∘ resolution. The radial parameterization within the oceans uses 50 shell layers, following the irregular discretization of the ECCO model. The seawater monthly conductivities from NOAA's WOA were interpolated to the same grid via bilinear formula in angular coordinates and weighted averaging in radial coordinate, which preserves the total conductance.
The model was run from January 1997 through to the end of November 2001. Global results were extracted from the middle of every month (e.g., 17 January, 15 February, 18 March, and 17 April 1997), but daily results were extracted along the transect of the cables' paths.
To compare numerical predictions with the processed seafloor cable observations, the electric
field was integrated along the seafloor between the endpoints of each cable. For each cable element,
the electric field component along the cable direction was calculated in the lowermost ocean
discretization layer. The linear trend was finally removed from each time series of predicted
cable voltages.
Results and discussion
Figures , , and summarize the processed voltages and their numerical
predictions from the ElmgTD ECCO-based simulation for individual cables.
Panel a in each figure shows the time series of cable voltages processed with the 90 d
knotted spline fit and the 30.5 d knotted spline fit in red and green, respectively. In the case
of the HAW1 cables, the HAW1N and HAW1S branches are distinguished by solid and dashed lines, respectively, and the blue line shows the results of the numerical predictions. A linear trend was removed from all shown time series.
Panel b in Figs. , , and shows the numerical predictions of the voltage gradient (i.e., the electric field) on the
seafloor, along the respective cables, before integration. Finally, panel c in Figs. , , and displays
the transport T⟂ of the ECCO model across each cable. Note that by transport in this
context we denote the vertically integrated velocity component perpendicular to the cable for each cable element
position and time (hence the unit of m2 s-1). Although it is not a direct input to the numerical simulations (contrary
to the velocities in individual ECCO layers), it serves as a useful proxy for discussions below.
The results for the OKI cable. (a) The smoothed time
series of cable voltages using 30.5 and 90 d knot separation are shown in red and green, respectively. The blue line
corresponds to the predictions obtained by the numerical model.
(b) The time development of the voltage gradient along the cable length from the
3-D model. In panel (c), we plot the ECCOv4r4 vertically integrated transport
across the cable in a similar way. The cable orientation follows Table , from Honshu to Okinawa.
The results for the HAW1 cables. The HAW1N and HAW1S cables are distinguished by the solid and dashed
lines in panel (a). The cable orientation is from California to Hawaii.
Otherwise, the description corresponds to Fig. .
Looking first at the common features of the results for all cables, we note, as expected from
basic geometrical considerations, a general similarity between the voltage gradient along the
cable and the water transport across the cable T⟂. We can use these to discuss the
effect of individual currents on the numerical predictions. However, while the ocean flows are certainly
the dominant term controlling the induced electric fields, the additional contributions of other effects
yield a much richer spatiotemporal structure. The main field variations in both space and time can have a
linear impact on the large-scale features, as implied by the forcing term of the EM induction
Eq. (). Moreover, the local variations in seawater conductivity, the bathymetry,
and the sediment thickness affect the electric field in a non-linear way. In particular, the
toroidal magnetic mode, which corresponds to the poloidal electric currents and stems
from the galvanic coupling between the ocean and the underlying solid Earth, can play an
important role .
The results for the HAW3 cable. The cable orientation is from California to Hawaii.
The description corresponds to Fig. .
Upon closer inspection of the OKI cable results, the importance of the Kuroshio Current stands out,
at the distance of 300–600 km from Honshu (Fig. c). It produces
the largest contribution to the predicted voltages by far (Fig. a, b).
In terms of spatial distribution, the positions along the cable where
the largest contributions to the electric field are induced do not match with
the peak positions of the cross-cable transport. This discrepancy can be attributed
to the electrically strongly heterogeneous environment caused by large bathymetry changes
in the vicinity of the Ryukyu arc.
The ECCO model suggests an increase in the transport in the last months of 2000,
which is consequently responsible for the increased voltage in the numerical model. However,
no such increase is present in the observed voltages, and this discrepancy remains an
open question. If we trust the OKI voltages, it is possible that the ECCO model is
overestimating the Kuroshio strength by the end of 2000.
In the case of HAW1N and HAW1S, the numerical model predicts significantly smaller amplitudes
of cable voltage variations when compared with the observations (Fig. ).
The California Current is
the main contributor to the total voltages, at distances up to 1000 km from the Californian coast.
The spatiotemporal distribution of the cross-cable transport and the induced voltages is
in good agreement due to an uncomplicated electrical conductivity distribution in the
deep ocean.
The ocean transports across the HAW1 cables demonstrate larger seasonal variations than in the case
of Kuroshio. However, the lack of significant contributions perpendicular to the cable as well as the
changing direction of these flows both along the cable and in time yield poor agreement
of the total integrated voltage with the observations.
The HAW3 cable, in contrast, shows good agreement between the predicted
and observed voltages (Fig. ). The numerical model is again dominated by the
California Current, which is here closer to the coast. The HAW3 cable lies a bit to the
south of the HAW1N and HAW1S cables, and it is also within the low-speed region of the Eastern Pacific Gyre. The transport across the cable in the central Pacific is more coherent,
yielding slightly stronger signals than in the case of the HAW1 cables. Again,
the pattern of the cross-cable transport is well matched with the spatiotemporal map of
the induced voltages.
In Table , we calculated two sets of correlation coefficients.
In the second column of the table, the voltages predicted by the numerical model
were correlated with the total ECCO-based water transport (in m3 s-1) across
the respective cables:
P⟂=∫startendT⟂dl.
These values are independent of the actual cable voltage measurements and can provide
an upper limit on what can be achieved by the interpretation of long-cable voltages
in terms of ocean flows. Large correlations were obtained for the HAW1 and HAW3
cable locations, whereas the integrated flow across the OKI cable was poorly correlated
with the predicted voltage. This stresses the importance of the accurate modeling
of the induced electric field in strongly heterogeneous areas.
In the third column of Table , we show the correlation coefficients between
the predicted and observed voltages using the 30.5 d knot separation
datasets. Due to the gaps present in the data, the Gaussian kernel method was applied.
It is obvious that the discrepancies between the predicted and observed voltages
are still large, and significant efforts are required both on the side of data processing
and numerical modeling to reconcile the results. The OKI cable in particular presents
an interesting case. Although the total cross-cable transport is less correlated
with the predicted voltages than in the case of both HAW1 cables, the agreement
with the observations is considerably better. This again points to the role
of local EM effects.
For individual seafloor cables, we show the correlation coefficients between the cable voltages
predicted by the numerical model and the total ECCO-derived water transport across the cable P⟂ in the second column.
The third column shows the correlation coefficients between the predicted and observed voltages for 30.5 d spline knot
separation.
On the side of numerical modeling, one could devise a comparison study between
different ocean models. Indeed, we have used our model to predict the magnetic fields
of the LSOMG (Large Scale Ocean Model for Geophysics) model in the past , and we have also attempted the calculation
of the cable voltages for the eddy-resolving GLORYS (Global Ocean reanalysis and Simulation) ocean model (not shown here). One
problem related to this approach is the volume of computational resources necessary to carry out
the calculations. As the cable voltages are sensitive to local electric fields, the usual
simplifications of the EM induction solver, based on the thin-sheet approximation, or representing
the oceans by a single layer with integrated water transports and electrical conductances,
are problematic . The single 5-year calculation of the full physical
model presented here, with 50 ocean layers and spherical harmonic truncation degree 240,
required about 105 CPU hours to complete. Semi-global or regional modeling tools with
local refinement ability are needed for more accurate numerical studies.
The qualitative comparison of the induced voltages and water transports along the cables,
as presented in this paper, could be made more exact by applying the principal component analysis methodology. When applied only to the water transports provided by different
ocean models, it could reduce the burden of calculating a detailed 3-D EM response to each model
and allow a more focused interpretation of the observed voltages. We plan to carry out such analysis
in the future.
The studies by evaluating transport in the Straits of Florida from seafloor voltage cable data had correlation values corresponding to much higher values than those of this study. As shown in Fig. 20 of , the correlation squared values ranged from 0.61 to 0.94. However, Larsen's study was fundamentally different: the seafloor voltage cable was an order of magnitude shorter than the cables considered in this study and the Gulf Stream within the Straits of Florida has large speeds as well as coherent velocities flowing perpendicularly to the cables. Therefore, overall, the study had a more ideal signal-to-noise ratio.
Conclusions
We present an evaluation of using seafloor voltage cables for monitoring circulation across oceanic basins. We compare processed seafloor voltage cable data with the numerical predictions produced using an electromagnetic induction solver, fed by flow velocity estimates from the data-assimilated ECCO model and seawater electrical conductivity climatologies from the NOAA World Ocean Atlas. We find that the correlation between cable voltage data and numerical predictions strongly depends on both the amplitude and direction of the flow velocities across the cable.
Due to the computational constraints, the calculations of the ocean-induced electric field presented here are limited to a single realization of ocean velocity estimates: the ECCO model. Therefore, beside the unmodeled or uncorrected signals in the seafloor cable voltages, a first-order source of discrepancy between the numerical prediction and observation is the inaccuracy of the flow velocity estimates. Therefore, an extended analysis and comparison of the cable voltage calculations driven by other ocean circulation models is desirable in the future as well as
the consideration of direct velocity observations .
While much improvement is needed before utilizing seafloor voltage cables to study and monitor ocean circulation across large regions, we believe that seafloor voltage cables can eventually be used to study and monitor large-scale ocean flow. The cables used in this study were installed for telecommunication purposes – there was no regard for whether these cables would be best suited to monitor ocean currents. Flow information can most reliably be extracted from seafloor voltage cable data when the flow has mostly unidirectional, perpendicular velocities across the cable. For our study, the OKI cable was in the area with the largest velocities, but because it is oriented mostly parallel to the Kuroshio Current, its correlation would likely greatly improve if it was instead perpendicular to the current's flow.
If voltage cables were strategically placed on the seafloor between Antarctica and Chile (a distance of ∼ 700 km) or between Antarctica and New Zealand (a distance of ∼ 1300 km), the correlation between data and predictions could be quite high, due to both the shorter cable length (compared with the HAW1 and HAW3 cables) and the relatively uniform and large flow velocities. Indeed, seafloor voltage cables may be a very effective method for measuring and continuously monitoring the flow of the Antarctic Circumpolar Current – which is definitely something worth investigating.
Using existing cables, the correlation between data and numerical predictions will likely also improve if the methodology is enhanced to remove induced signals from
seasonal variations in ionospheric signals.
Data availability
The data and numerical predictions discussed in this study are freely available for download at
https://geomag.colorado.edu/OCEM (Velímský et al., 2021).
Author contributions
NRS and MCN conceived the questions and methodology of this study
and wrote the initial version of the paper.
NRS and NPT worked on the processing of cable data.
MCN supervised NRS and NPT on work related to this project and provided useful feedback on improving the paper. JV carried out the numerical
modeling and also contributed to the paper, in particular to the revised version.
Competing interests
The authors declare that they have no conflict of interest.
Acknowledgements
The computational resources were provided by the Ministry of Education, Youth and Sport of the Czech Republic, from the Large Infrastructures for Research, Experimental Development and Innovations project “IT4Innovations National Supercomputing Center – LM2015070” (project ID OPEN-13-21). We thank Zoltan Szuts and the anonymous reviewer for their helpful comments.
Financial support
This research has been supported by the NASA Earth and Space Science Fellowship (grant no. 80NSSC17K0450), the CIRES Innovative Research Program 2017 and the Grant Agency of the Czech Republic (grant no. P210/17-03689S).
Review statement
This paper was edited by Erik van Sebille and reviewed by Zoltan Szuts and one anonymous referee.
ReferencesBaringer, M. O. and Larsen, J. C.: Sixteen years of Florida Current Transport at 27∘ N, Geophys. Res. Lett., 28, 3179–3182, 2001.Chave, A. D., Filloux, J. H., and Luther, D. S.: Electromagnetic induction by ocean currents: BEMPEX, Phys. Earth Planet. In., 53, 350–359, 10.1016/0031-9201(89)90021-6, 1989.
Chave, A. D., Luther, D. S., Lanzerotti, L. J., and Medford, L. V.:
Geoelectric field measurements on a planetary scale: oceanographic and
geophysical applications, Geophys. Res. Lett., 19, 1411–1414, 1992.
Cox, C. S., Filloux, J. H., and Larsen, J. C.: Electromagnetic studies of
ocean currents and electrical conductivity below the ocean-floor, in: The
Sea, Wiley, New York, ISBN: 0674-01732-3, 637–693, 1971.Denig, W. F.: Geomagnetic kp and ap Indices, available at: http://www.ngdc.noaa.gov/stp/GEOMAG/kp_ap.html (last access: 15 February 2021), 2015.Everett, M. E., Constable, S., and Constable, C. G.: Effects of near-surface conductance on global satellite induction responses,
Geophys. J. Int., 153, 277–286, 10.1046/j.1365-246X.2003.01906.x, 2003.Faraday, M.: The Bakerian Lecture, Experimental Researches in Electricity, Terrestrial Magneto-electric Induction,
Philos. T. R. Soc. Lond., 122, 163–194, 10.1098/rstl.1851.0001, 1832.Finlay, C. C., Maus, S., Beggan, C. D., Bondar, T. N., Chambodut, A., Chernova, T. A., Chulliat, A., Golovkov, V. P., Hamilton, B., Hamoudi, M., Holme, R., Hulot, G., Kuang, W., Langlais, B., Lesur, V., Lowes, F. J., Lühr, H., Macmillan, S., Mandea, M., McLean, S., Manoj, C., Menvielle, M., Michaelis, I., Olsen, N., Rauberg, J., Rother, M., Sabaka, T. J., Tangborn, A., Tøffner-Clausen, L., Thébault, E., Thomson, A. W. P., Wardinski, I., Wei, Z., and Zvereva, T. I.: International Geomagnetic Reference Field: the eleventh generation, Geophys. J. Int., 183, 1216–1230,
10.1111/j.1365-246X.2010.04804.x, 2010.
Flosadóttir, Á. H., Larsen, J. C., and Smith, J. T.: Motional
induction in North Atlantic circulation models,
J. Geophys. Res., 102, 10353–10372, 1997.Forget, G., Campin, J.-M., Heimbach, P., Hill, C. N., Ponte, R. M., and Wunsch, C.: ECCO version 4: an integrated framework for non-linear inverse modeling and global ocean state estimation, Geosci. Model Dev., 8, 3071–3104, 10.5194/gmd-8-3071-2015, 2015.
Fujii, I. and Utada, H.: On Geoelectric Potential Variations Over a Planetary
Scale, PhD thesis, The University of Tokyo, Tokyo, Japan, 81 pp., 2000.Fukumori, I., Wang, O., Fenty, I., Forget, G., Heimbach, P., and Ponte, R. M.: ECCO version 4, release 3, Tech. Rep., JPL/Caltech and NASA Physical Oceanography, 10 pp., 1721.1/110380, 2017.Grayver, A. V., Munch, F. D., Kuvshinov, A. V., Khan, A., and Sabaka, T. J.:
Joint inversion of satellite-detected tidal and magnetospheric signals
constrains electrical conductivity and water content of the upper mantle and
transition zone, Geophys. Res. Lett., 44, 6074–6081,
10.1002/2017GL073446, 2017.
Irrgang, C., Saynisch, J., and Thomas, M.: Ensemble simulations of the
magnetic field induced by global ocean circulation: estimating the
uncertainty, J. Geophys. Res.-Oceans, 121, 1866–1880, 2016a.Irrgang, C., Saynisch, J., and Thomas, M.: Impact of variable seawater conductivity on motional induction simulated with an ocean general circulation model, Ocean Sci., 12, 129–136, 10.5194/os-12-129-2016, 2016b.Irrgang, C., Saynisch, J., and Thomas, M.: Utilizing oceanic electromagnetic
induction to constrain an ocean general circulation model: A data
assimilation twin experiment, JAMES, 9, 1703–1720,
10.1002/2017MS000951, 2017.
Koyama, T.: A study on the electrical conductivity of the mantle by voltage
measurements of submarine cables, PhD thesis, University of Tokyo, Tokyo, Japan, 130 pp., 2001.
Lanzerotti, L. J., Thomson, D. J., Meloni, A., Medford, L. V., and Maclennan,
C. G.: Electromagnetic study of the Atlantic continental margin using a
section of a transatlantic cable, J. Geophys. Res., 91, 7417–7427, 1986.
Lanzerotti, L. J., Medford, L. V., Kraus, J. S., Maclennan, C. G., and
Hunsucker, R. D.: Possible measurements of small-amplitude TID's using
parallel, unpowered telecommunications cables, Geophys. Res. Lett.,
19, 253–256, 1992a.
Lanzerotti, L. J., Sayres, C. H., Medford, L. V., Kraus, J. S., and Maclennan, C. G.: Earth potential over 4000 km between Hawaii and California, Geophys. Res. Lett., 19, 1177–1180, 1992b.Lanzerotti, L. J., Medford, L. V., Maclennan, C. G., and Thomson, D. J.:
Studies of Large-Scale Earth Potentials Across Oceanic Distances,
AT&T Tech. J., 74, 73–84, 10.1002/j.1538-7305.1995.tb00185.x, 1995.
Lanzerotti, L. J., Medford, L. V., Maclennan, C. G., Kraus, J. S., Kappenman,
J., and Radasky, W.: Trans-atlantic geopotentials during the July 2000 solar
event and geomagnetic storm, Solar Physics, 204, 351–359, 2001.
Larsen, J. C.: Electric and Magnetic Fields Induced by Deep Sea Tides,
Geophys. J. Roy. Astr. S., 16, 47–70, 1968.Larsen, J. C.: Transport measurements from in-service undersea telephone
cables, IEEE J. Oceanic Eng., 16, 313–318,
10.1109/48.90893, 1991.Larsen, J. C.: Transport and heat flux of the Florida Current at 27∘ N derived from cross-stream voltages and profiling data: theory and observations, Philos. T. Roy. Soc. A, 338, 169–236,
10.1098/rsta.1992.0007, 1992.
Larsen, J. C. and Sanford, T. B.: Florida current volume transports from
voltage measurements, Science, 227, 302–304, 1985.
Malin, S. R. C.: Separation of lunar daily geomagnetic variations into parts
of ionospheric and oceanic origin,
Geophys. J. Roy. Astr. S., 21, 447–455, 1970.Manoj, C., Kuvshinov, A., Neetu, S., and Harinarayana, T.: Can undersea
voltage measurements detect tsunamis?, Earth Planets Space, 62,
353–358, 10.5047/eps.2009.10.001, 2010.Meinen, C. S., Smith, R. H., and Garcia, R. F.: Evaluating pressure gauges as a potential future replacement for electromagnetic cable observations of the
Florida Current transport at 27∘ N,
J. Oper. Oceanogr., 0, 1–11, 10.1080/1755876X.2020.1780757, 2020.Pedatella, N. M., Liu, H., and Richmond, A. D.: Atmospheric semidiurnal lunar tide climatology simulated by the Whole Atmosphere Community Climate Model, J. Geophys. Res., 117, A06327, 10.1029/2012JA017792, 2012.Rehfeld, K., Marwan, N., Heitzig, J., and Kurths, J.: Comparison of correlation analysis techniques for irregularly sampled time series, Nonlin. Processes Geophys., 18, 389–404, 10.5194/npg-18-389-2011, 2011.
Sabaka, T. J., Tyler, R. H., and Olsen, N.: Extracting ocean-generated tidal
magnetic signals from Swarm data through satellite gradiometry, Geophys. Res. Lett., 43, 3237–3245, 2016.Šachl, L., Martinec, Z., Velímský, J., Irrgang, C., Petereit, J., Saynisch, J., Einšpigel, D., and Schnepf, N. R.: Modelling of electromagnetic signatures of global ocean circulation: physical approximations and numerical issues, Earth Planets Space, 71, 58, 10.1186/s40623-019-1033-7, 2019.Sanford, T. B.: Motionally induced electric and magnetic fields in the sea,
J. Geophys. Res., 76, 3476–3492, 10.1029/JC076i015p03476, 1971.Schnepf, N. R.: Going electric: Incorporating marine electromagnetism into
ocean assimilation models, J. Adv. Model. Earth Sy., 9,
1772–1775, 10.1002/2017MS001130, 2017.Schnepf, N. R., Nair, M., Maute, A., Pedatella, N. M., Kuvshinov, A., and
Richmond, A. D.: A Comparison of Model-Based Ionospheric and Ocean Tidal
Magnetic Signals With Observatory Data, Geophys. Res. Lett., 45,
7257–7267, 10.1029/2018GL078487, 2018.
Shimizu, H., Koyama, T., and Utada, H.: An observational constraint on the
strength of the toroidal magnetic field at the CMB by time variation of
submarine cable voltages, Geophys. Res. Lett., 25, 4023–4026,
1998.
Sleijpen, G. L. G. and Fokkema, D. R.: BiCGstab(ell) for Linear Equations
involving Unsymmetric Matrices with Complex Spectrum,
Electron. T. Numer. Ana., 1, 11–32, 1993.
Spain, P. and Sanford, T. B.: Accurately monitoring the Florida Current with motionally-induced voltages, J. Mar. Res., 7, 843–870, 1987.Szuts, Z. B., Bower, A. S., Donohue, K. A., Girton, J. B., Hummon, J. M.,
Katsumata, K., Lumpkin, R., Ortner, P. B., Phillips, H. E., Rossby, H. T.,
Shay, L. K., Sun, C., and Todd, R. E.: The Scientific and Societal Uses of
Global Measurements of Subsurface Velocity, Frontiers in Marine Science, 6,
358, 10.3389/fmars.2019.00358, 2019.Tyler, R. H., Maus, S., and Lühr, H.: Satellite observations of magnetic fields due to ocean tidal flow, Science, 299, 239–241,
10.1126/science.1078074, 2003.Tyler, R. H., Boyer, T. P., Minami, T., Zweng, M. M., and Reagan, J. R.:
Electrical conductivity of the global ocean, Earth Planets Space, 69, 156,
10.1186/s40623-017-0739-7, 2017.Vanyan, L. L., Utada, H., Shimizu, H., Tanaka, Y., Palshin, N. A., Stepanov,
V., Kouznetsov, V., Medzhitov, R. D., and Nozdrina, A.: Studies on the
lithosphere and the water transport by using the Japan Sea submarine cable
(JASC): 1. Theoretical considerations, Earth Planets Space, 50, 35–42,
10.1186/BF03352084, 1998.Velímský, J.: Determination of three-dimensional distribution of
electrical conductivity in the Earth's mantle from Swarm satellite data:
Time-domain approach, Earth Planets Space, 65, 1239–1246,
10.5047/eps.2013.08.001, 2013.Velímský, J. and Martinec, Z.: Time-domain, spherical
harmonic-finite element approach to transient three-dimensional geomagnetic
induction in a spherical heterogeneous earth,
Geophys. J. Int., 161, 81–101, 10.1111/j.1365-246X.2005.02546.x, 2005.
Velímský, J., Šachl, L., and Martinec, Z.: The global
toroidal magnetic field generated in the Earth's oceans,
Earth Planet. Sc. Lett., 509, 47–54, 10.1016/j.epsl.2018.12.026, 2019.Velímský, J., Schnepf, N. R., Nair, M. C., and Thomas, N. P.: Download Ocean Circulation Electromagnetism (OCEM)
Data, available at: https://geomag.colorado.edu/OCEM,
last access: 17 February 2021.