ORIGINAL RESEARCH article
Sec. Geomagnetism and Paleomagnetism
Volume 11 - 2023 | https://doi.org/10.3389/feart.2023.1194230
Impact of variable seawater conductivity on ocean wave-induced electromagnetic fields simulated with finite difference method
- 1College of Marine Geo-sciences and Key Lab of Submarine Geo-sciences and Prospecting Techniques of Ministry of Education, Ocean University of China, Qingdao, China
- 2National Engineering Research Center of Offshore Oil and Gas Exploration, Beijing, China
Electric fields generated by the motion of ocean waves through the Earth’s ambient geomagnetic fields and the induced secondary magnetic field can be observed at the seafloor and at the sea-surface, and even in the air. Most of current studies on ocean wave-induced electromagnetic fields assume that seawater conductivity is constant, and ocean waves are treated as regular waves with a fixed amplitude and frequency. However, these assumptions are inconsistent with actual ocean conditions. In this paper, we present a finite difference algorithm for simulating the ocean wave-induced electromagnetic fields with variable seawater conductivity. We investigate impacts of variable seawater conductivity on the electromagnetic fields induced by the wind waves and swell as well as mixed ocean waves, which are treated as the superposition of a number of regular waves with different frequencies and amplitudes, and analyze the characteristics of the induced electromagnetic fields.
Oceanic waves can carry charged ions dissolved in seawater, and electromagnetic fields can be induced by the dynamo interaction of ocean waves with the geomagnetic field. The ocean wave-induced electromagnetic fields can be observed at various locations, including the seafloor, sea surface, and, sometimes, even at satellite altitudes (Crews and Futterman, 1962; Cox et al., 1978; Minami, 2017). These electromagnetic fields can provide sufficient information for the inversion of the ocean wave spectrum, which is widely used in ocean engineering and seakeeping considerations in ship design (Cieutat et al., 2003; Techet, 2005; Nielsen and Dietz, 2020). However, they will seriously interfere with natural magnetotelluric (MT) signals and reduce the quality of MT data.
In the past decades, the study of ocean wave-induced electromagnetic fields focused on the plane surface gravity waves with a fixed amplitude and frequency, which are usually considered the approximation of ocean waves. Weaver (1965) and Fraser (1966) calculated and observed the magnetic fields induced by surface gravity waves in the infinite ocean. Larsen (1971) developed a general theory of ocean wave-induced electromagnetic fields, especially long and intermediate surface gravity waves. Miles et al. (1977) established an analog model of magnetic fields induced by surface gravity waves in the laboratory. Ochadlick (1989) measured the magnetic fields generated by surface gravity waves at an aircraft. Lilley and Weitemeyer (2004a) calculated the apparent aeromagnetic wavelengths of the magnetic signals of the ocean swell. Lilley et al. (2004b) also observed the magnetic fields generated by ocean waves near the sea surface. Semkin and Smagin (2012) investigated the effect of self-induction exhibited by the ocean wave-induced electromagnetic fields. Shimizu and Utada (2015) demonstrated that the ocean wave-induced electromagnetic fields were rarely affected by the conductive subseafloor media.
It is well-known that actual ocean waves are not regular waves with a fixed amplitude and frequency but rather the sum of multiple-frequency amplitudes that exhibit a particular sea state with a significant wave height and peak period (Rasool et al., 2021). Recently, a more realistic and accurate ocean wave-induced electromagnetic field model has been proposed by combining surface gravity waves with the wave spectrum, which provides the distribution of ocean wave energy empirically (Chave, 1983; Ailliot et al., 2013; Ryabkova et al., 2019), and thus, connections between ocean wave-induced electromagnetic fields and the ocean wave spectrum are established. The electromagnetic fields generated by wind waves can be solved efficiently using wind speed and significant wave height, and other parameters (Yaakobi et al., 2011). Meanwhile, many studies focus on the wind-driven electromagnetic fields at different wind speeds (Zhu and Xia, 2014; Zhan and Pan, 2019). To the best of our knowledge, the electromagnetic fields generated by both the swell and mixed ocean waves have not been discussed in the literature.
In the aforementioned studies, seawater conductivity is often assumed to be constant. However, seawater conductivity is variable in the realistic ocean (Zheng et al., 2018). Previous studies have shown that there is an approximately linear relationship between seawater conductivity and temperature, and the salinity of seawater also influences the conductivity of seawater. Both the salinity and temperature of ocean water can exhibit significant spatial and temporal variations, thereby leading to variations in sea water conductivity. The mean conductivity of seawater is approximately 3.3–4 S/m, while it may reach 6 S/m above the main thermocline. Chave and Luther (1990) and Irrgang et al. (2016) investigated the effect of vertically varying seawater conductivity on ocean current-induced electromagnetic fields. As far as we know, the impact of changes in seawater conductivity on the ocean wave-induced electromagnetic fields has not been well-investigated.
In this paper, we simulate the ocean wave-induced electromagnetic fields with variable seawater conductivity using the finite difference method. The impact of variable seawater conductivity will be investigated, especially within the main thermocline at different seasons and latitudes. Then, the wave spectra of wind waves, swell, and mixed waves are presented and the characteristics of the electromagnetic fields are analyzed.
FD simulation of surface gravity wave-induced electromagnetic fields
The electromagnetic fields induced by surface gravity waves can be obtained by solving the motional induction equation. Here, we adopt two-dimensional surface waves propagating in the x-direction, with no variation in the y-direction. The z-axis points downward, and the origin of the Cartesian coordinate system is at the averaged sea surface. We assume that seawater is incompressible and flow is irrotational. This means that there exists a velocity potential
The gradient of the velocity potential
The frequency domain Maxwell equations in the quasi-stationary state approximation can be written as
Compared to the spatial variability of ocean waves
Since the impact of the conductivity of the subseafloor material on the electromagnetic fields within the frequency range of ocean waves at and above the seabed is negligible (Larsen, 1971; Minami and Toh, 2013; Shimizu and Utada, 2015), an air–seawater–seafloor three-layer conductivity model, as shown in Figure 1, is considered here. Assuming that the conductivity of air is
For the length scale of ocean waves, the horizontal variation in seawater conductivity is negligible compared to its vertical variations (Chave and Luther, 1990; Tyler et al., 2017). Hence, Eq. 10 can be simplified as (Minami et al., 2021)
The magnetic fields in the air can be expressed as
Similarly, the magnetic fields in the seafloor layer can be expressed as
At the air–ocean interfaces
The boundary value problems (11) and (14) can be solved using the finite difference method. We consider that the one-dimensional grid and the space from the sea surface to the seafloor are spatially discretized into
Using the Taylor series expansion, the second derivative of the normal
Eq. 11 can then be approximated using the symmetric difference equation
The first derivative of the magnetic field component
The boundary conditions in Eq. 14 become
One obtains a linear system of N-1 equations, which can be written in the following matrix form:
The aforementioned equations can be solved by using a direct solver to obtain
It should be noted that the electric fields calculated in this section are in situ electric fields. The geomagnetic electro-kinetograph field
The magnetic induction intensity
Validation of FD simulation
In order to validate the accuracy of the finite difference (FD) algorithm described in the previous section, we simulated the wave-induced electromagnetic fields in the model, as shown in Figure 1. The conductivity in the air is set to be
Figure 2 shows the amplitudes of ocean wave-induced electrical and magnetic field components obtained using the FD method. For comparison, the analytic solutions calculated using the formulae used by Shimizu and Utada (2015) are also shown. The FD numerical results agree very well with the analytic solutions. The relative errors of all three components (Ey, Bx, and Bz) are below 0.5%. Figure 2shows that 1) the amplitudes of both
FIGURE 2. Amplitudes of electromagnetic fields induced by ocean waves; the blue triangles and the black solid line represent the FD results and analytic solution, respectively. Bx (A), Bz (B), Ey (C), and relative errors (D).
Impact of variable seawater conductivity on ocean wave-induced electromagnetic fields
Generally, the conductivity of seawater decreases with the depth of the ocean. It may exhibit either a linear or exponential decrease. In this section, we evaluate how seawater conductivity affects the electromagnetic field generated by ocean waves in both linear and exponential models. Additionally, we utilize empirical formulas to calculate the ocean wave-induced electromagnetic field while also factoring in the effects of the thermocline layer. In the following numerical simulation, the amplitude and frequency of the surface gravity wave are set to be 1 m and 0.1 Hz, respectively.
The effect of the depth-averaged seawater conductivity
To investigate the impacts of seawater conductivity distribution on ocean wave-induced electromagnetic fields, we conducted modeling studies for two cases. In the first case, the seawater conductivity is set to be a constant of 3 S/m, 4 S/m, and 5 S/m. In the second case, seawater conductivity linearly varies with depth z, with a depth-averaged conductivity of 3 S/m, 4 S/m, and 5 S/m (Figure 3A). The ambient geomagnetic field is
FIGURE 3. (A) Seawater conductivity profiles of six conductivity distributions; vertical profiles with different depth-averaged conductivity distributions are presented for (B) Bx, (C) Ey, and (D) Bz.
Figures 3B–D show the ocean wave-induced electromagnetic fields with different seawater conductivity distributions. One can see that 1) maintaining a constant depth-averaged seawater conductivity but changing from a uniform to a linear distribution with depth results in an increase in the induced electromagnetic fields. Here, the variations in seawater conductivity distribution led to an increase of up to 25% in the electromagnetic fields; 2) for a fixed seawater conductivity distribution (i.e., a constant gradient of linear conductivity variation), changes in the depth-averaged seawater conductivity also affect the induced electromagnetic field. It is evident that the induced electromagnetic field increases with the increase in the depth-averaged seawater conductivity.
The effect of the seawater conductivity gradients
To investigate the effects of the seawater conductivity gradients, we consider the model with the same depth-averaged conductivity but different conductivity gradients. The depth-averaged conductivity gradient
Figure 4A shows the variation in seawater conductivity for depth-averaged conductivity gradients of 0.0025, 0.006, and 0.011. As shown in Figures 4B–D, the induced electromagnetic fields increase with the increase in the depth-averaged conductivity gradient.
FIGURE 4. (A) Seawater conductivity with different depth-averaged conductivity gradients; (B) Bx, (C) Ey, and (D) Bz.
The effect of the thermocline
The thermocline is the transition layer between the warmer mixed water near the surface and the cooler deep water below. It plays an important role in marine ecology, meteorological forecasting, underwater communication, and aquaculture. In the body of the thermocline, there is a sudden temperature change. The temperature decreases rapidly from the mixed layer to the cooler deep layer.
The temperature t, salinity S, and pressure p in seawater affect the variations in seawater conductivity. According to the 1978 Practical Salinity Scale (Fofonoff and Millard, 1983), the seawater conductivity
The conductivity ratio
As a high linear correlation exists between temperature and conductivity in the thermocline, the conductivity in it may include the temperature dependence only and a constant salinity of 36 psu can be considered in Eq. 29 (Tyler et al., 2017; Zheng et al., 2018).
The thermocline largely depends on the seasons and latitudes. Figure 5 displays the variation in temperature and seawater conductivity in the thermocline for different latitudes (Figure 5A) and different seasons (Figure 5B). The temperature profiles vary at different latitudes as the surface water is warmer near the equator and colder at the poles. In low-latitude tropical regions, the sea surface is much warmer, leading to a highly pronounced thermocline, while in polar regions, the temperature is fairly constant at all depths. Seasons also have an impact on the vertical variation in seawater temperature and conductivity. In winter, the sea surface temperature is lower than that in summer, resulting in a thicker thermocline. In contrast, in summer, due to solar radiation, the thermocline is thinner and has more pronounced variations.
FIGURE 5. (A) Representative temperature and corresponding conductivity profiles for tropical, mid-latitude, and polar regions, and for (B) winter and summer.
The geomagnetic field also varies with latitudes. The intensity of the magnetic field is generally higher near the magnetic poles and lower near the equator. At the magnetic equator, where the proximity to the geographic equator is high, the magnetic field lines intersect Earth’s surface, resulting in a smaller magnetic inclination angle. However, as the latitude increases, the magnetic inclination angle gradually becomes larger (Table 1).
TABLE 1. Geomagnetic field intensity, magnetic declination, and inclination vary with latitudes (polar, mid-latitude, and tropical regions) and seasons (summer and winter).
Figure 6 shows the electromagnetic fields in the thermocline for tropical, mid-latitude, and polar regions, and for winter and summer, respectively. Figure 6 shows that 1) compared to tropical and polar regions, the mid-latitude regions have the highest amplitude values of induced electromagnetic fields (Figure 6A). This is because tropical regions have relatively high seawater conductivity but low geomagnetic field values, while mid-latitude regions have slightly lower seawater conductivity than tropical regions but significantly higher geomagnetic field values; 2) due to the relatively insignificant seasonal variations in the geomagnetic field, the values of induced electromagnetic fields are larger in summer than those in winter (Figure 6B).
FIGURE 6. (A) Induced electromagnetic fields of ocean waves in thermoclines for tropical, mid-latitude, and polar regions and (B) for winter and summer, respectively.
Numerical simulation of electromagnetic fields induced by ocean wave based on sea wave spectrums within inhomogeneous ocean
The motion of realistic ocean waves can be described as a stationary random process. The ocean waves can be simulated by the wave spectrum, which represents the statistical characteristic of ocean wave motion (Longuet-Higgins, 1962; Grainger et al., 2021). Thus, the height and velocity potential of the ocean waves can be simulated by linear superposition:
We consider the wave-induced electromagnetic field in a specific area of the South China Sea, taking into account the actual temperature, salinity, and density of seawater. The area is located at latitude 18°N and longitude 119.336°E. The water depth in this region is 294 m. The seawater conductivity can be calculated using Eqs 25–30. Figure 7 shows the distribution of temperature, salinity, and seawater conductivity in this area. The geomagnetic induction intensity in this region is 40,174.6 nT with a magnetic declination of 47° and a magnetic inclination of 1.1°.
FIGURE 7. Distribution of temperature (the blue dashed line), salinity (the blue solid line), and seawater conductivity (the greed solid line) in a certain area of the South China Sea.
Electromagnetic fields induced by wind waves
The JONSWAP spectrum is commonly used to describe the spectrum of wind waves (Hasselmann et al., 1973):
Figure 8 shows that with the increase in both the wind speed and wind fetch, the magnitude of the spectrum increases but the peak frequency decreases.
FIGURE 8. JONSWAP spectrum for different wind speeds and fetches V = 15 m/s and x = 9,000 m (the blue solid line); V = 10 m/s and x = 9,000 m (the red solid line); V = 8 m/s and x = 9,000 m (the yellow solid line); V = 8 m/s and x = 6,000 m (the purple dashed line); and V = 8 m/s and x = 3,000 m (the green dashed line).
The wind-induced electromagnetic fields at the sea surface can be obtained by combining Eqs 20, 23, 32 with the JONSWAP spectrum (34) and are depicted in Figure 9. One can observe that 1) the amplitude of the induced electric and magnetic fields increases with the increase in wind speed; 2) the magnitude of the electrical and magnetic fields increases with the expansion of the wind fetch; 3) as the wind speed increases and the wind fetch expands, the dominant peak frequency gradually shifts toward lower frequencies; and 4) the wind-induced electromagnetic fields are considerably smaller in magnitude at the seafloor (dashed lines) than those near the sea surface (solid lines).
FIGURE 9. Electromagnetic fields generated with the (A) wind speed of 8 m/s, 10 m/s, and 15 m/s (a wind fetch of 9,000 m) and the (B) wind fetch of 3,000, 6,000, and 9,000 m (a wind speed of 8 m/s).
Electromagnetic fields induced by a swell
Wind waves travel in a great circle route after being generated, and after moving out of the area of the wind fetch, the waves are called swell waves and can travel thousands of kilometers. As short-wavelength waves carry less energy and dissipate faster, swell waves often have a relatively long wavelength (Gao et al., 2022). The swell spectrum is given by (Lucas and Guedes Soares, 2015)
FIGURE 10. Swell spectrum for three different peak shape parameters γS = 2, 5, and 8, and HS = 5 m, and ωP = 0.3142 rad/s.
Figure 10 shows that with the increase in the shape peak parameter
Figure 11 shows the electromagnetic fields generated by the ocean swell at both the sea surface (A) and seafloor (B). Figure 11 shows that 1) as the shape peak parameter
FIGURE 11. Electromagnetic field spectrum for peak shape parameters γS = 2, 5, and 8 at the sea surface (A) and seafloor (B).
The electromagnetic fields induced by mixed ocean waves
In some cases, the realistic ocean waves are mixed with wind waves and swell, called mixed ocean waves. The former refers to wind waves in equilibrium with the local wind, while the latter is defined as the swell waves generated elsewhere and not significantly affected by the local wind at that time (Hwang et al., 2012; Garcia-Gabin, 2015).
In order to describe complicated sea states due to the coexistence of wind waves and ocean swell, a double-peaked spectrum is proposed by combining the wind wave system and the swell system (Guedes Soares, 1984, 1991; Rossi et al., 2021). Ochi and Hubble (1976) presented a double-peaked spectrum as a combination of two gamma distributions as follows:
Due to the difference in the effect caused by the wind waves and swell, the spectrum of the mixed ocean waves can be divided into three types: the strong wind and weak swell type, weak wind and strong swell type, and equivalent wind and swell type. The classification of the three types can be determined by the
Table 2 shows the parameters of the mixed ocean waves. Figure 12 shows that the energy in the spectrum of weak wind and strong swell type is usually concentrated at low frequencies, but the energy in the spectrum of the strong wind and weak swell type is usually concentrated at high frequencies. The energy distribution is more balanced in the spectrum of the equivalent wind and swell type.
TABLE 2. Parameters of mixed ocean waves; type 1: strong wind and weak swell type; type 2: equivalent wind and swell type; and type 3: weak wind and strong swell type.
FIGURE 12. Double-peaked spectrum of the weak wind and strong swell type (the blue line), the equivalent wind and swell type (red line), and the strong wind and weak swell type (the green line).
The electromagnetic fields induced by three types of mixed ocean waves are shown in Figure 13. From Figure 13, we can see that 1) the three-type mixed wave spectrum can be easily recognized by the x- and z-components of the induced magnetic intensity strength spectrum observed at the sea surface (Figure 13A); 2) due to the high-frequency filtering effect of the ocean, the electromagnetic field spectra at the seafloor (Figure 13B) mainly retain the low-frequency information on the wave spectrum and may not reflect the realistic spectrum of the mixed ocean waves.
FIGURE 13. Electromagnetic fields induced by three types of mixed ocean waves at the sea surface (A) and seafloor (B).
In this paper, we presented a finite difference algorithm for simulating ocean wave-induced electromagnetic fields with varying seawater conductivity. Our method allows us to simulate and investigate the influence of arbitrary variation in seawater conductivity on wave-induced electromagnetic fields. Our numerical examples show that seawater conductivity has a significant effect on wave-induced electromagnetic fields. We also investigated the effects of the thermocline on wave-induced electromagnetic responses at different latitudes and seasons. We observed the maximum amplitude of induced electromagnetic fields in the mid-latitude regions, while we observed the minimum amplitude of induced electromagnetic fields in the polar regions. This is mainly influenced by both the geomagnetic field and temperature of the thermocline at different latitudes. Furthermore, in mid-latitude regions, the thermocline conductivity exhibits seasonal variations, with a greater decrease during summer than that during winter, leading to a larger induced electromagnetic field. The latitude and seasonal variations in the geomagnetic field and seawater conductivity are crucial for the accurate simulation of ocean wave-induced electromagnetic fields.
Moreover, we simulate the electromagnetic field spectra for wind waves, ocean swell, and mixed ocean waves in the inhomogeneous ocean. We find that the energy of the wind wave-induced electromagnetic field is predominantly concentrated in the high-frequency band, while the energy of the swell wave-induced electromagnetic field is mainly concentrated in the low-frequency band. Additionally, we simulate three types of mixed wave-induced electromagnetic fields based on the combination of wind and swell strengths, and analyze their characteristics.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding author.
JG implemented the algorithms, drew the figures, and wrote the manuscript. YL led the research, supervised the study, and edited the manuscript. All authors contributed to the article and approved the submitted version.
This work was supported by the National Science Foundation of China (91958210 and U2241201).
The authors thank Ji Cai for helping in manuscript editing.
Conflict of interest
YL is a adjunct professor of the National Engineering Research Center of Offshore Oil and Gas Exploration.
The remaining author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Gao, X., Ma, X., Ma, Y., Huang, X., Zheng, Z., and Dong, G. (2022). Spectral characteristics of swell-dominated seas with in situ measurements in the coastal seas of Peru and Sri Lanka. J. Atmos. Ocean. Technol. 39, 755–770. doi:10.1175/JTECH-D-21-0143.1
Hasselmann, K., Barnett, T. P., Bouws, E., Carlson, H., Cartwright, D. E., Enke, K., et al. (1973). Hamburg, Germany: Deutches Hydrographisches Institut.Measurements of wind-wave growth and swell decay during the joint north sea wave project (JONSWAP)
Hwang, P. A., Ocampo-Torres, F. J., and García-Nava, H. (2012). Wind Sea and swell separation of 1D wave spectrum by a spectrum integration method. J. Atmos. Ocean. Technol. 29, 116–128. doi:10.1175/JTECH-D-11-00075.1
Irrgang, C., Saynisch, J., and Thomas, M. (2016). Impact of variable seawater conductivity on motional induction simulated with an ocean general circulation model. Ocean. Sci. 12, 129–136. doi:10.5194/os-12-129-2016
Lilley, F. E. M., Hitchman, A. P., Milligan, P. R., and Pedersen, T. (2004a). Sea-surface observations of the magnetic signals of ocean swells. Geophys. J. Int. 159, 565–572. doi:10.1111/j.1365-246X.2004.02420.x
Rasool, S., Muttaqi, K. M., and Sutanto, D. (2021). Modelling Ocean waves and an investigation of ocean wave spectra for the wave-to-wire model of energy harvesting. Eng. Proc. 12, 51. doi:10.3390/engproc2021012051
Rossi, G. B., Crenna, F., Berardengo, M., Piscopo, V., and Scamardella, A. (2021). Investigation on spectrum estimation methods for bimodal sea state conditions. Sensors (Basel) 21, 2995. doi:10.3390/s21092995
Ryabkova, M., Karaev, V., Guo, J., and Titchenko, Yu. (2019). A review of wave spectrum models as applied to the problem of radar probing of the sea surface. J. Geophys. Res. Oceans 124, 7104–7134. doi:10.1029/2018JC014804
Techet, A. H. (2005). 13.42 design principles for ocean vehicles. https://ocw.mit.edu/courses/2-22-design-principles-for-ocean-vehicles-13-42-spring-2005/fea036f3255fa4824279668c7ff8c6cd_r1_lti.pdf.
Yaakobi, O., Zilman, G., and Miloh, T. (2011). Detection of the electromagnetic field induced by the wake of a ship moving in a moderate sea state of finite depth. J. Eng. Math. 70, 17–27. doi:10.1007/s10665-010-9410-z
Keywords: electromagnetic fields, variable seawater conductivity, wind waves, swell, mixed ocean waves
Citation: Ge J and Li Y (2023) Impact of variable seawater conductivity on ocean wave-induced electromagnetic fields simulated with finite difference method. Front. Earth Sci. 11:1194230. doi: 10.3389/feart.2023.1194230
Received: 26 March 2023; Accepted: 19 June 2023;
Published: 24 July 2023.
Edited by:Jin Li, Hunan Normal University, China
Reviewed by:Ren Zhengyong, Central South University, China
Benyu Su, China University of Mining and Technology, China
Nian Yu, Chongqing University, China
Xiaoyue Cao, Yangtze University, China
Copyright © 2023 Ge and Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yuguo Li, firstname.lastname@example.org