Influence of Expendable Current Profiler Probe on Induced Electric Field of Ocean Currents

The expendable current profiler (XCP) is a new instrument that is internationally used to rapidly monitor ocean currents in marine environments. The most crucial part of this instrument is the XCP probe. Since the probe is of high electrical resistance, it acts almost like an insulator with respect to seawater. Placing it into the induced electric field (IEF) of seawater therefore yields a certain level of influence over the electric field.Therefore, in order to improve the accuracy of XCPmeasurements, the conditions associated with this influence can be used to guide the design of XCP probes; at the same time, these can also serve as reference points in order to provide technical support for the processing of XCP data on ocean currents. To this end, computer-based numerical simulations and laboratory-based physical simulations are used in this study. The results showed that after an XCP probe (diameter: 5 cm; length: 52 cm) was inserted into seawater, the voltage difference of ocean currents at both ends of the electric field sensor placed above the XCP probe increased by a factor of 1.85 (as compared to the case in which there is no influence from the probe).


Introduction
The United States began researching expendable temperature and velocity profilers (XTVPs) as early as the 1970s and in 1978 successfully developed the first XTVP [1,2].Between 1979 and 1980, a research team led by Sanford and Sippican company made coproduction and sea trial of several hundred XTVP probes and obtained preliminary detection results [3,4].The company subsequently changed the name of this device from "XTVP" to "expendable current profiler" (XCP) and launched production of the instrument, which became widely used in marine surveys, scientific research, and national defense [4,5].The XCP is a type of expendable profiling instrument for marine environments and can obtain the profile information of ocean currents rapidly [6].For the first time in China, we have conducted an in-depth study of various XCP technologies [7,8].After independent research and development, China's first set of XCP equipment was manufactured using precision design [7].Multiple marine tests indicated that placing the XCP probe into seawater changed the IEF of ocean currents.The XCP can be deployed through a probe launch or be manually cast from the carrying platform of ships, submarines, and aircraft.It can quickly measure ocean currents and temperature profiles while sinking and can calculate water depth based on the probe's sinking velocity [9,10].Data are then transferred to the carrying platform through wired or wireless communication modes, and the real-time data of ocean currents and temperature variation with respect to water depth is obtained after data processing.To verify the accuracy of the XCP measurements, an acoustic Doppler current profiler (ADCP) is used for comparison.The used ADCP is an OS-75K from the RDI Company.The results of seawater experiments showed that the ocean current velocities measured by XCP and ADCP are very similar [7].
We found that placing this instrument in seawater led to changes in the induced electric field (IEF) that is generated by the movement of seawater.Hence, in-depth studies on the influence of XCP probes on the IEF of ocean currents are required in order to improve detection accuracy [11,12].In this study, computer-based numerical simulations and laboratory-based physical simulations were used to examine the influence of the XCP probe on the IEF of ocean currents and to obtain the corresponding coefficients of influence.

Integral Equation Method.
The ellipsoidal probe was placed in homogeneous seawater with an electrical conductivity of  1 .The probe itself has an electrical conductivity of  2 ( ⃗ ), which is a function of .⃗  represents the radius vector.In addition, because the influence of magnetic permeability is usually minimal compared to that of electrical conductivity, it was assumed that  =  0 .
It was assumed that the electric dipole source  was located somewhere in space and that the time-harmonic factor was   .We began with Maxwell's equations for the frequency domain as follows: The response of homogeneous seawater was defined as the primary field and was represented by the subscript ""; the primary field also satisfied the following Maxwell equations: At this stage, the following equations were obtained by subtracting the formulas in (2) from the corresponding formulas in (1): where  is the practical conductivity value of the model.The conductivity values inside and outside of the probe equal  2 ( ⃗ ) and  1 , respectively.At this stage, (4) can be rewritten as follows: If the difference between the total and primary field is considered to be the secondary field (represented by the subscript ""), then (5) can be simplified as follows: where the following equation holds: This is known as the scattered current, which would exist only within the probe.Equation (3) can be similarly simplified: The secondary field can be treated as being caused by the scattered current   .Because   = −ẑ + (1/ŷ)∇(∇ ⋅ ), the secondary field in the seawater can be calculated using the following equation: where   and   are the secondary vector potential and secondary scalar under Lorentz condition, respectively: In (10), (,   ) is Green's function for the scalar quantity of total space, which was derived from the following equation: where When the probe is in half-space, an additional item must be added to (9) to reflect the influence of the interface.This additional item has been described by Hohman and Wannamaker et al. under the conditions of homogeneous and layered ground, respectively [13][14][15].This leads to the following expression for the secondary field: As described previously, the electromagnetic field is made up of two portions, the primary and secondary fields: Substituting (12), which represents the secondary field, into (13) yields the total electric field's singular Fredholm integral equation of the second kind: where   (,   ) represents dyadic Green's function.This function, which was required because the direction of the electric field at location  would be different from that of the source current at location   , can be derived using the following equation: where ∇  represents the derivative for the prime coordinate system and  is the unit dyadic.
The probe was partitioned into  cubic units, each of which had a length of Δ [16].The electric conductivity within each unit was constant.Because the scattered current   within each unit would also be constant, the integral equation (3) can be approximated as follows: where  2 and   represent the electric conductivity and field of the th unit, respectively.
When Γ(,   ) is used to represent the integral term of ( 16), the latter can be written as follows: where Γ is dyadic Green's function for a small current integration.It is different from , which is dyadic Green's function for a relatively infinitesimal current element.
Under these conditions, the electric field at the centre of the th unit can be written as follows: After transposition, it can be written as where the following holds: Here,  is a 3 × 3 unit matrix, and 0 is the zero tensor.When every  is expressed using (20), the following block matrix equation can be obtained: in the matrix [], each element is itself a 3 × 3 matrix: Equation ( 22) can then be used to solve for the value of the electric field at the centre of each unit within the probe.At this stage, the electric field at any position outside the probe can be obtained using ( 16).

Theoretical Model.
A theoretical model for the XCP probe (Figure 1(a)) was first established for conducting numerical simulations [17].In the figure, AB and A  B  refer to the electric dipoles.The used current was 1000 A, and the used frequencies were 0, 0.0001, and 1 Hz.The XCP probe was placed in an infinite amount of seawater with a resistivity of 0.33 Ω⋅m.A Cartesian coordinate system was then established by assuming that seawater flow only occurs at the sea surface.The origin is located at the sea surface, while the -axis points in the direction of the ocean currents, and there are no ocean currents along the -axis.The and -axes were both located at the sea surface, while the -axis pointed upward and was perpendicular to the sea surface.The -, -, and -axes conform to the right-hand rule.The conductivity of the probe was assumed to be 0 S/m, while the observation plane was located at  = −1000 cm.The length of the probe's minor axis was, respectively, at 1, 2, 4, 5, 6, 10, and 20 cm, while that of its major axis along the -direction was 52 cm.Mesh decomposition was carried out for the probe model [18], resulting in 40 × 40 × 52 grids (Figure 1(b)).

Probe's Influence on Electric Field Distribution.
The integral equation method [13,14] was used to carry out theoretical calculations of the probe model's electromagnetic response.Figures 2(a) and 2(b) show the distribution characteristics of the   component of the disturbance field caused by the XCP probe along the and -axis-label on the observation plane [19].Along the direction of the -axis, the electrical current field was affected by the high-resistance probe, causing a repulsion of the current or electric field.As can be seen from the figures, the electric field decreases as the observation points moved nearer to the probe along the -axis, which exhibits negative anomalous characteristics.When the repulsed current line passes through the left and right sides of the XCP probe, the current density correspondingly increases.This results in positive anomalous characteristics.The influence of the XCP probe was minor and would eventually disappear when the observation points were located far enough from the probe.Based on the boundary conditions that current density was continuous in the normal direction at the outer and inner sides of the XCP probe, we can obtain the following equations: The terms  1 and  2 are the scattered currents inside and outside of the XCP probe, respectively;  1 and  2 are the electric field intensities inside and outside of the XCP probe, respectively; and  1 and  2 are the densities inside and outside of the XCP probe, respectively.Because the resistivity of the XCP probe was higher than that of seawater, the electric field within the probe was much greater than the external anomalous electric field.
The distribution characteristics of the   component's anomalous field and current-induced electric field, known as background field ratios, are shown in Figures 2(c  There are positive anomalous characteristics on the left and right sides of the probe, where the ratio of the anomalous field to the background field is 7.
The simulation results below indicated that the amplitude of the electric field on the observation plane was significantly influenced by the probe.The distribution was symmetrical along the -axis, whereas the ratio of the anomalous field to the background field varied with respect to the location of the observation points.In this situation, in-depth computational analyses must be carried out in order to determine whether placing the probe in seawater would influence the measured voltage.

Calculating the Probe's Influence on Measured Voltages.
Electric field distributions obtained from forward modelling were required for the analysis of voltage changes.These were then used to calculate the electric potential difference between the two electrodes [20].The electric potential difference was actually obtained by integrating the electric field's intensity vectors measured along a particular path between one electrode and the other.(The two electrodes were located 5 cm apart at both ends of the probe's minor axis and constituted two observation points on a straight line.) The solid circular line represents the cross-section passing through the centre of the instrument, whereas the dots represent the outer side of the instrument.The outer diameter of the instrument was 5 cm, and the major axis of the ellipsoid was 52 cm.
When the electric field vectors  ⇀   and  ⇀   were separately applied to (25) below, the electric potential differences caused by the background and anomalous fields between the two electrodes were obtained as ΔΦ  and ΔΦ  , respectively, after integration.ΔΦ  corresponds to the electric potential difference between the two electrodes when the probe was not inserted, whereas ΔΦ  + ΔΦ  corresponds to the electric potential difference between the electrodes after insertion of the probe.This can be used to calculate the increase of  between the original electric potential difference of the two electrodes and that after insertion of the probe, as shown in (26).
By changing the angle between the two electrodes and the external electric field, we were able to observe the way in which  varied for various azimuths.The calculation results when the diameter of the probe's minor axis was 5 cm are shown in Figure 3.It can be seen that the electric potential difference caused by the background and anomalous fields exhibited similar patterns as the azimuth varied.The calculated  values indicated that measurements of the electric field were significantly influenced after insertion of the probe.The average of the various  values was calculated, which yielded a value of 1.85 in this example.

Variations in Diameter of Probe's Minor Axis and the
Influence on .The calculated values of  are shown in Table 1.It can be seen that changes in the value of  were bigger when the minor axis was 5 cm long; the values of  were relatively lower for lengths of 20 and 52 cm.It turns that all of the conditions are in the range of the errors permitted.

Variations in Signal Frequency of Electric Dipole Source
and the Influence on .Numerical simulations of the electric field were performed using the model in which the length of the probe's minor axis was 5 cm; the signal frequencies of the electric dipole source were then varied.The frequencies used during the simulations were 0, 0.0001, and 1 Hz.Data concerning the electric field obtained from the simulations were then used to calculate the corresponding changes in the value of  as a function of the azimuth; the results are shown in Figure 4.The values of  were relatively large (average: 1.85) when the frequency of the transmission signal was 0 Hz.These decreased slightly (average: 1.84) at the frequency of 0.0001 Hz but were relatively small (average: 1.77) at the frequency of 1 Hz.These results indicate that the voltage increase rate exhibited a decreasing trend as the frequency increased.

Physical Simulations of Probe's Influence on Electric Field of Ocean Currents.
In order to further understand the trends discussed above and to confirm the probe's influence on IEF measurements, physical simulations were carried out.

Test Environment.
The electric ionization currents of industry create too much interference when a regular water tank is used to measure electric fields.Hence, physical simulations for this study were conducted indoors using a large plastic container.In order to simulate a genuine seawater environment, the conductivity of seawater was adjusted to be 3.3 S/m.The main used instruments and equipment consisted of a signal recovery 7265 DSP lock-in amplifier, which operated over a frequency range of 1 mHz to 250 kHz; model 7265 offers full-scale voltage sensitivities down to 2 nV and current sensitivities to 2 fA.We also used Matrix MPS-3003L-3 (voltage display precision: three and a half A/D conversion digital display, ±0.5% + 2 words; current display precision: three and a half A/D conversion digital display, ±1% + 2 words).Furthermore, we used Tektronix TDS 2002 (vertical resolution: 8 bits; vertical sensitivity: 2 MV to 5 V/div; DC vertical precision: ±3%), Agilent 34420A (display resolution: 7 1/2 ; sensitivity: 100 pV/nΩ), Hewlett Packard 33120A (accuracy at 1 kHz: ±1% of specified output), Victor VC9801A + (DC voltage: ±(0.5% + 3); AC voltage: ±(0.8% + 5); DC current: ±(0.8% + 10); AC current: ±(1.0%+ 15)).A personal computer (PC), two copper plates that supplied electricity, two small Ag|AgCl nonpolarizable electrodes, and a solid highresistance cylinder (outer diameter: 5 cm) were also used.The underwater measuring environment dictates that the signal of the ocean currents' electric field must undergo a transmission process from a liquid to a solid medium.Normal nonpolarizable electrodes create electrochemical noise when the electrodes come into mutual contact, which is extremely unfavourable when observing the weak IEF signals within ocean currents [21].For this reason, the first step in carrying out physical simulations was to search for an electrode material with a small yet steady polarization potential when placed in a marine environment.Many past experiments have shown that when silver and silver chloride in powder form are mixed according to a specific formula and then made into electrodes using metallurgical processes, the latter exhibit good electrochemical properties when placed in seawater [22].This could be explained through the conductive mechanism of Ag|AgCl itself.First, compared to other electrode materials, it is easier to refine Ag to its pure state under laboratory conditions, thereby eliminating any "battery effect" [23] caused by impurities (one of the noise sources).Pure Ag also has better electrochemical stability in an environment with a generally constant temperature.(The temperature of seawater can be almost constant across an extremely short period of time and within a specific area.)Second, Clions are the material carriers of seawater conductivity.After AgCl comes into contact with seawater, the same chemical composition is also the main carrier that conducts electricity at the contact surface between the solid and liquid phases [24].

Test Contents.
The tests were carried out as the copper plates were charged and discharged and as the conditions for supplying electricity and making measurements were varied; the electric field was stabilized during these situations.Under direct current conditions, the voltage was measured before and after the high-resistance cylinder was placed in the brine.Under alternating current conditions, the voltage and supply current were measured before and after the high-resistance cylinder was placed in the brine.

Tests Using Electric
Field with Direct Current.The voltage source Matrix supplied 2 V of constant voltage to the copper plates in the simulation water tank; the copper plates were connected to the PC and Agilent 34420A.This process involved first placing the high-resistance cylinder into the water for 30 min of data collection, followed by another 30 min of data collection after removal of the cylinder.Five groups of data were measured during the tests.The measurements for Group 1, as shown in Figures 5(a  period of 2-3 min.During the measurement process, the supply current signals were also recorded at the appropriate intervals.The purpose of this was to normalize the electric currents, thereby eliminating the influence of current changes on the measurement data.The measured voltages under the two test conditions were approximately 78 and 145 mV, respectively.The maximum relative error caused by the potential difference between the nonpolarizable electrodes within 30 min did not exceed 1%.Thus, the influence of the potential difference between electrodes on the measurements was ignored during subsequent calculations and analyses.The measurements for Group 2 are shown in Figures 5(c) and 5(d).For each group of data, the average of two voltages was used to calculate the ratio of the electrode voltage with the high-resistance cylinder to that without the cylinder (i.e., the voltage increase rate ).Data for the five sets of measurements are shown in Table 2.

Tests Using Electric
Field with Alternating Current.The 7265 DSP lock-in amplifier by signal recovery was used for these tests.The alternating signals from the lock-in amplifier were loaded onto the copper plates of the water tank.The size of the signal was adjusted to control the strength of the signal received by the electric field sensor.The test data are shown in Table 3.

Conclusions
This study investigated the influence that the probes used for XCP detection have on the IEF of ocean currents.Computerbased numerical simulations were used as the basis for physical tests that were made in simulated marine environments.Based on the assumption that both simulations were carried out under similar conditions, the conclusions from the theoretical analysis are as follows.The amplitude of the

Figure 1 :
Figure 1: (a) Schematic diagram of the expendable current profiler (XCP) probe.(b) Mesh decomposition results for the probe model.

Figure 2 :
Figure 2: (a) Distribution of the   component of the anomalous field on the plane  = 1000 cm.(b) Distribution of the   component of the anomalous field along the and -axis-label on the plane  = 1000 cm.(c) Ratio of the   component's anomalous field to the background field on the plane  = 1000 cm.(d) Ratio of the   component's anomalous field to the background field along the and -axis-label on the plane  = 1000 cm.

Figure 3 :
Figure 3: Changes in electric potential difference and voltage enhancement factor as a function of the azimuth.(a) Electric potential difference caused by the background field.(b) Electric potential difference caused by the anomalous field.(c)  as a function of the azimuth.

Figure 4 :
Figure 4: Changes in voltage increase rate as a function of the azimuth for different signal frequencies.

Figure 5 :
Figure 5: (a) Measured voltage between electrodes for Group 1 after being supplied with electricity (with high-resistance cylinder).(b) Measured voltage between electrodes for Group 1 after being supplied with electricity (without high-resistance cylinder).(c) Measured voltage between electrodes for Group 2 after being supplied with electricity (with high-resistance cylinder).(d) Measured voltage between electrodes for Group 2 after being supplied with electricity (without high-resistance cylinder).

Table 1 :
Average value of  under different conditions.

Table 2 :
Test data for influence of high-resistance cylinder on electric field.