Three-dimensional structure of mesoscale eddies and their interaction with Kuroshio based on observations from a CPIES array

Mesoscale eddies are widespread in the oceans; however, tracking and observing them remains a challenge. Additionally, there are few studies of eddy–Kuroshio interactions based on large-scale in situ data. Data from a large array of 39 CPIESs deployed around the Luzon Strait in 2018–2019 (12 deployed east of Taiwan within a 3 × grid) were used to reveal the three-dimensional structure of mesoscale eddies and their interaction with Kuroshio. The results show a strengthening (weakening) of the tilt of the pycnocline, which has also caused a deeper (shallower) pycnocline in eastern Taiwan following an anticyclonic (cyclonic) interaction with Kuroshio, which resulted in a strengthening (weakening) of the Kuroshio. During the anticyclonic eddy impinging Kuroshio, the downward (upward) movement at the center of the anticyclonic eddy (Kuroshio) results in a positive (negative) anomaly structure for temperature, with corresponding positive (negative) anomalies for salinity above 600 m in the intermediate layer and negative (positive) anomalies below 600 m. In contrast, when the cyclonic eddies interact with the Kuroshio, the temperature and salinity anomalies change simultaneously, with the temperature showing an overall negative anomaly and the salinity showing a negative and positive anomaly above and below 600 m, respectively. The research in this study provides the basis for further research on energy exchange during eddy–Kuroshio interactions.


Introduction
The Kuroshio is a strong western boundary current that originates from the northward branch of the North Equatorial Current bifurcation off the east coast of the Philippines (Zhang et al., 2001;Qiu and Chen, 2010b). Kuroshio flows northward mainly along Luzon and flows into the East China Sea and the south coast of Japan via waterways east of Taiwan, while bringing heat and salinity from low equatorial latitudes to mid-latitudes, thus affecting climate and ocean-atmosphere interactions along the flow path (Nitani, 1972;Kwon et al., 2010).
Mesoscale eddies are spatially scaled eddies of 50-500 km and temporally scaled eddies of days to months in the ocean, which are mainly characterized as circulations that remain closed for longer periods of time (Chelton and Schlax, 1996;Chelton et al., 2007;Chelton et al., 2011). Depending on the direction of rotation, mesoscale eddies are divided into anticyclonic eddies and cyclonic eddies, in the ideal case, while the central water mass of an anticyclonic (cyclonic) eddy exhibits a downward (upward) vertical motion, while the situation in the ocean is very complex and is related to the life cycle of the eddy (Koszalka et al., 2009;Buongiorno Nardelli, 2013;McGillicuddy, 2016;Viúdez, 2018). The instability of currents, atmospheric forcing, and topography all play an important role in the formation of mesoscale eddies. As a carrier of oceanic material and energy transport, mesoscale eddies are commonly found in the oceans, and their kinetic energies are one order of magnitude larger than the average kinetic energies in most of the ocean, and they make important contributions to oceanic circulation, ocean energy balance, water mass distribution, thermal salinity, and nutrient transport (Chelton et al., 2011;Yang et al., 2013;Dong et al., 2014;Xu et al., 2014;Zhang et al., 2014;Zhang et al, 2016). In the subtropical countercurrent (STCC) region of the mid-latitude western Pacific, a large number of eddies are generated due to strong shear and oblique pressure instability (Qiu, 1999). These eddies propagate westward since their generation to near the western boundary of the Pacific Ocean and some eddies can intrude the South China Sea (Guo et al., 2017). The complex oceanic background circulation east of Taiwan is not only the main current area of the Kuroshio, but also one of the areas where the mesoscale eddies in the northwest Pacific are highly active ( Figure 1). It has been shown that the mesoscale eddies are one of the main factors causing the intraseasonal variability of the Kuroshio east of Taiwan, with a period of 100 days (Zhang et al., 2001;Liu et al., 2004;Qiu and Chen, 2010a;Chang et al., 2015;Jan et al., 2015;Ren et al., 2020). The existing studies on the influence of mesoscale eddies on Kuroshio mainly rely on numerical models, and the evolution of the eddy-Kuroshio interaction combined with the measured data remains less studied due to the insufficient data of in situ long-time observation arrays, and the dynamic mechanism of mesoscale eddies affecting Kuroshio still needs to be further clarified.
The study of mesoscale eddies mainly relies on satellite altimeter data, combined with eddy tracking algorithms to obtain long-term surface eddy dynamics information (Chelton et al., 2007;Chelton et al, 2011). With the large number of Argo profiling buoys deployed in the global oceans, a large amount of upper ocean temperature and salt data have been accumulated, and the combination of satellite altimeter data, in situ survey data, and reconstruct eddy methods has greatly facilitated the study of mesoscale eddies. In the northwest Pacific Ocean, we have studied the three-dimensional structure of mesoscale eddies and their heat transport and have found that the eddies can affect the depth of 1,000 m (Qiu et al., 2005;Hu et al., 2011;Nan et al., 2011;Liu et al., 2012;Yang et al., 2013;Zhang et al., 2013;Chu et al., 2014;Zhang et al., 2014a;Zhang et al., 2015;Shu et al., 2016;Dong et al., 2017;Chen et al., 2018;He et al., 2018). Zhang et al (2016) used the S-MEE project to deploy a mooring system submersible array in the northeastern part of the South China Sea and combined it with synthetic analysis to describe the three-dimensional structure of the mesoscale eddies and yield the process of eddy energy dissipation. In the latest study,  combined satellite observations and global Argo data to find a general vertically tilted structural feature of the global eddy, further contributing to our understanding of the three-dimensional structure of the mesoscale eddy (Li et al., 2022a;Li et al, 2022b). Despite the large number of studies on mesoscale eddies, tracking eddies and their dynamical processes from in situ observations is a challenge because eddies vary in both time and space, and most studies can only focus on the state of the eddy at a single moment in time.
The Current-Pressure Inverted Echo Sounder (CPIES), which integrates current and pressure sensors, has been widely used in oceanographic surveys because of its small size and easy deployment. Combined with the gravest empirical mode (GEM) method, we can obtain four-dimensional temperature, salinity, density, and current data of the whole water layer in the observation area. Based on the data of a large-scale CPIES array, this paper will construct mesoscale three-dimensional structure characteristics, reveal the isopycnic variation process of one of the eddy current interaction process mechanisms, and provide a basis for the study of mass transport and energy dissipation of mesoscale eddies. The surface geostrophic currents and CPIES array distribution around east of Taiwan. Red triangles denote CPIES stations; color shading is the bottom topography from the ETOPO1 Global Relief Model, and average surface geotropic currents (black arrows) between January 2018 and July 2019 are from AVISO data, covering the observation time of CPIES. The CPIES array was deployed from January 2018 and July 2019. The topographic data were from the following website: https://www.ngdc.noaa.gov/mgg/global/etopo1sources.html 2 Materials and methods

CPIES data
An inverted echo sounder (IES) equipment with a current sensor (C) and a pressure sensor (P) is called CPIES, which consists of a fishing float, a glass float, an Aanderaa current meter, a 50-m communication cable, a PIES body, and a heavy anchor. The working mode is mainly to measure the sound propagation time from the sea surface to the sea bottom. Since its inception in the 1970s, IES has been widely used in marine research with the help of technological developments and the progress of data processing methods (Meinen et al., 2003;Donohue et al., 2008;. This has relied heavily on the ground-transfer empirical modal GEM method proposed by  for data inversion methods. This method establishes an empirical relationship between sound propagation time (t) and temperature, salinity, and streamflow functions by integrating all historical measured temperature and salinity profiler for the study area, creating a Lagrangian matrix of scatter points and projecting the existing temperature and salinity data onto this two-dimensional space. Also, an observation array consisting of multiple CPIES can be based on the GEM method to obtain an absolute current profile. A reference layer needs to be chosen to create GEM, and this reference layer requires a depth that can basically cover the depth of baroclinic influence. Considering the weak baroclinic below 1,000 m east of Taiwan, this study takes 1,000 m as the reference layer. As a comparison, the 1,000 m selected by Anders et al. (2017) in this region is used as the reference layer.
From December 2017 to July 2018, a total of 40 CPIES were deployed around the Luzon Strait, and 39 sites were successfully recovered in July 2019. Among them, 12 of the CPIES were deployed in the area east of Taiwan with a 3 × 4 size array. All CPIES observations were then preprocessed, which included being despiked, tidally corrected, and dedrifted, following Kennelly et al. (2007). All the final raw data were low-pass filtered with a 72-h window and a fourth-order Butterworth filter. For further information on the data processing and inversion process, see Zhao et al. (2022) and Zheng et al. (2022).

Satellite altimeter data
An Archiving, Validation, and Interpretation of Satellite Data in Oceanography (AVISO) altimetry data (http://marine.copernicus.eu/ services-portfolio/access-to-products/) set was used in this study to track the mesoscale eddies. Sea surface height anomaly data are used in the manuscript, the data spatial resolution is 1/4°× 1/4°, and the data set extended from June 2018 to July 2019, which provided satisfactory overlap with the CPIES array data.

Historical data
In order to establish the relationship between the GEM array inversion temperature, salinity, and current, as well as the acoustic wave propagation time (t) and the pycnocline, it is necessary to draw on as many existing historical data profiles as possible. Historical hydrographic data from the World Ocean Database (WOD) and Argo were selected from the locations shown in Figure 2, with longitudes of 121.3-127°E and latitudes of 20-23.8°N. The locations selected were larger than the CPIES deployment sites. The total number of profiles is 3,155, of which 1,367 profiles are less than 1,500 m, 1,660 profiles are between 1,500 and 2,000 m, and 128 profiles are greater than 2,000 m.

Relationship between sea level anomaly, t, nd pycnocline depth
Changes in the pycnocline depth have been shown in many studies to be one of the mechanisms of interaction between mesoscale eddies and Kuroshio (Ren et al., 2020). The time series of t measurement from CPIES station C04 and sea level anomaly (SLA) are shown in Figure 3, and the result reveals a highly negative correlation between SLA with a coefficient of -0.9. This means that the increase in sea surface height caused by the anticyclonic eddies (cyclonic eddies) can reduce (increase) the sound speed round trip propagation time.
To quantify the variation of the pycnocline at the interaction between eddies and Kuroshio, the measured propagation (t) time series needs to be converted to the pycnocline variation, and here, we use the depth of the maximum density gradient to represent the pycnocline variation. Pycnocline depth is defined in this paper by the maximum density gradient (dr/dz) below the seasonal thermocline (Tsai et al., 2015;Jan et al., 2017), calculated for each density profile located in the regions in the red point area from the historical Argo data ( Figure 2). Thus, it is necessary to establish the relationship between the propagation time t and the pycnocline depth. As mentioned before, in order to invert the temperature, salinity, and other data, the measurement time t by PIES needs to be converted into the reference (t) which is chosen to be 1,000 m. According to this, a function of the pycnocline depth (d pyn ) and t 1,000 are established based on historical CTD profile data. The result is plotted against the corresponding t 1,000 in each profile shown in Figure 4, the mean pycnocline is s=25.9 kg m −3 and the depth is approximately 300 m at eastern Taiwan. The t 1,000 and pycnocline depth exhibit a linear relationship with a correlation coefficient of 0.8, which means that pycnocline depth increases (decreases) with the increase (decrease) in t 1,000 . Therefore, we can obtain a simple corresponding function of the pycnocline depth and t 1,000 by linear fitting, shown by the black line in Figure 3B. Thus, we can invert the time series pycnocline depth at each time through t 1,000 inverted from t measured by CPIES. Combining the relationship between t 1,000 and SLA and pycnocline depth obtained from Figure 3, we can conclude that the positive (negative) SLA caused by the anticyclonic eddies (cyclonic eddies), which corresponds to a smaller (larger) t 1,000 is able to cause a deeper (shallower) pycnocline depth at the same time.

Interaction between mesoscale eddy and Kuroshio
During the observation period, multiple eddies propagating east of Taiwan are successfully captured by our CPIES array shown in Figure 5, and based on the relationship between the eddies and the pycnocline depth obtained in the previous section, we have the opportunity to directly study the interaction between mesoscale eddies and Kuroshio. At the anticyclonic eddy moment, t 1,000 forms a closed region of low values corresponding to the eddy. The minimum value of t 1,000 ithe center of the eddy is 1.316 s, and t 1,000 ll be larger the closer to the edge of the anticyclonic eddy. Additionally, it can be found that the distribution of t 1,000 shows obvious east-west differences, and t 1,000 in the region near the western Kuroshio is larger than that in the region controlled by the anticyclonic eddy, and the maximum value of t 1000 at the Kuroshio location is close to 1.320 s. The depth of the pycnocline at the center of the anticyclonic eddy is calculated from the linear relationship between t 1,000 and d pyn and is determined to be approximately 450 m, while the pycnocline depth at the Kuroshio side is approximately 410 m, while at the cyclonic eddy impinging with Kuroshio, t 1,000 forms a closed region of high values corresponding to the eddy shown in Figure 4A, with a maximum t 1,000 value of 1.322 s at the center of the eddy and decreasing toward the edge of the eddy. The distribution of t 1,000 shows a small eastwest difference, with a t 1,000 value of 1.321 s near the Kuroshio side, but is still smaller than the interior of the eddy. According to the relationship, the depth of the pycnocline in the center of the cyclonic eddy and Kuroshio side is approximately 365 m and 370 m, respectively.
Thus, the pycnocline depth at the anticyclonic eddy is significantly greater than at the cyclonic eddy, and the pycnocline is nearly parallel at the cyclonic eddy moment. Based on the possible variation mechanism proposed by Jan et al. (2015) and our calculations, the interaction between the eddy and the Kuroshio can be revealed by the schematic shown in Figure 6. Under normal conditions, the pycnocline of the Kuroshio is tilted to the lower right (white line in Figure 6) and the Kuroshio flows northward. While in the influence of the anticyclonic eddy ( Figure 6A), the overall pycnocline becomes significantly deeper and the tilt of the pycnocline increases significantly, exhibiting a "see-saw" effect, thus resulting in an increase in the Kuroshio velocity, a wider Kuroshio width, and a greater Kuroshio thickness, even up to 800-1,000 m, while at the cyclonic eddy influence ( Figure 6B), the pycnocline rises upward and reduces the tilt of the pycnocline, similar to a see-saw reverting to an intermediate state, corresponding to a decrease in the velocity and width of Kuroshio, as well as reducing the thickness of Kuroshio. The relationship between t SLA, and pycnocline; the SLA is the sea level anomaly from the AVISO product; d pyn is the pycnocline depth. (A) The measurement of t nd SLA at the CPIES site C02; the blue line and the orange line are the t and SLA, respectively. (B) The corresponding relationship between t 1000 and the pycnocline depth (d pyn ) calculated from historical CTD profiles in eastern Taiwan (red dots); the black solid line is the fitted curve.

Discussion
To study the changes of the hydrological background east of Taiwan under the interaction of eddies and Kuroshio, the average temperature and salinity characteristics at the study site are shown in Figure 7. It can be clearly seen that the temperature vertical gradient (red curve) shows the existence of two thermoclines, the seasonal thermoclines above 200 m and the permanent thermoclines within

The variation of anticyclonic eddies' interaction with Kuroshio
An anticyclonic eddy is captured by the CPIES array shown in Figure 8, where the temperature and salinity anomalies of the eddies are defined as the measured values at each point minus the average of each layer, and the center of each layer of the eddy is identified by the minimum velocity point, which is located at 123.5°E, 22.35°N. The left front of the anticyclonic eddy is adjacent to the Kuroshio, and the radius of this eddy is estimated to be approximately 70 km. The maximum velocity of the eddy at the surface layer is approximately 0.7 m/s, and the eddy extends to a depth of 1,000 m. Combined with Figure 4A, it can be seen that the Kuroshio velocity west of 122.5°E is 0.7 m/s. Ren et al. (2020)  when interacting with anticyclonic eddies, probably due to the fact that the CPIES array is located on the eastern side of the Kuroshio main axis and thus the velocities are relatively small. Taking a latitudinal section along the center of the eddy at 22.35°N shown in Figure 9, the temperature and salinity anomalies and velocity clearly show that the variation is completely different between the left and right sides of 122.5°E. The isotherm, isohaline, and pycnocline of the anticyclonic eddy are concave, while those in the Kuroshio side are convex. The largest positive temperature anomaly is located at the center of the eddy, with two positive temperature anomaly cores located at approximately 200 m and 500 m, with values of approximately 0.6°C and 1.4°C, respectively, which are basically located near the seasonal thermocline and the main thermocline. The salinity anomalies also show two cores below the subsurface layer, which are located at depths of 450 m and 850 m, with a maximum positive salinity anomaly of approximately 0.08 psu at 450 m. The characteristics of the salinity anomaly structure are related to the subsurface highsalinity water and intermediate waters. At the center of the anticyclonic eddy, downward movement of saline water from the subsurface will result in positive salinity anomalies at depths of 400-600 m. The downward transport of low-saline water in the intermediate layer will result in negative salinity anomalies below 600 m.
For the Kuroshio region west of 122.5°E, due to the see-saw like changes in the interaction between the anticyclonic eddy process, the variation of the Kuroshio is the opposite of the variation of the anticyclonic eddy. The Kuroshio side produces an overall negative temperature anomaly, and the salinity anomaly of the Kuroshio shows a negative-positive structure below 200 m. The upward movement of the low-salinity intermediate water at depths of 400-600 m results in a negative salinity anomaly at depths shallower than 600 m, with a maximum value of −0.06 psu at 450 m. Positive salinity anomalies of 0.04 psu are due to the upward movement of highly saline water in the deep sea below 700 m. Figure 9C shows a velocity profile at the moment of the anticyclonic eddies, which shows in detail the characteristics of the velocity structure distribution inside the Kuroshio and anticyclonic eddies. The velocity at the center of the anticyclonic eddies is almost 0 m/s, with the velocity gradually increasing from the center outward. The left side of the eddy is connected to the Kuroshio, and the larger velocity range in the horizontal direction extends to approximately 122.75°E, thus illustrating the increase in the amplitude of the Kuroshio caused by the anticyclonic eddies.

The variation of cyclonic eddies interaction with Kuroshio
A cyclonic eddy structure is shown in Figure 10; the center of this eddy is located at 123.75°E, 22.35°N, and the radius is estimated to be approximately 60 km from the sea surface height. The maximum velocity of the eddy is approximately 0.7 m/s, and can exceed 0.1 m/s at 1,000 m. Based on the latitudinal section of the eddy center at 22.35°N (Figure 11), the isotherm, isohaline, and pycnocline are found to be nearly horizontal, with a slight upward convexity at the center of the cyclonic eddy. The elevation of the pycnocline is evident; for example, the 26 kg m −3 density contour is located between 450 m and 500 m in a downward concave structure at the anticyclonic eddy, while at the cyclonic eddy, it is almost entirely above 400 m and is nearly horizontal. It can also be seen that the negative temperature anomaly has two cores, corresponding to two thermoclines with a maximum value of −0.8°C, located between 100 and 200 m. The cyclonic eddy-induced salinity anomaly shows a negative-positive structure below the subsurface layer. The core of the negative salinity anomaly is located at approximately 350 m with a maximum value of −0.06 psu, while the core of the positive salinity is located at approximately 800 m with a maximum value of 0.04 psu. The salinity anomalies can also be explained as follows: the upward movement of eddy center causes the low saline water between 400-600m to upward lift, resulting in negative salinity anomalies above the intermediate layer, while the salinity anomalies below 600m are positive due to the upward lift of the deeper high saline water.
For the cyclonic eddy impinging on the Kuroshio, the reduction in the tilt of the pycnocline has already been analyzed and the corresponding velocity distribution ( Figure 11) clearly depicts this reduction in velocity. The velocity of the Kuroshio decreases to approximately 0.2 m/s and the amplitude of the Kuroshio decreases significantly. This is in contrast to the velocity of the eastern side of the cyclonic eddy, which is approximately 0.4 m/s. When a cyclonic eddy interacts with Kuroshio, it weakens the tilt of the pycnocline, and if the cyclonic eddy is strong enough, it can even induce the depth of the pycnocline on the Kuroshio side to be lower than that on the eddy side, resulting in a reversal of the Kuroshio direction.

Conclusion
Using the extra-large CPIES array deployed in the area from July 2018 to July 2019, the relationship between propagation time t 1,000 , SLA, and pycnocline variability was established in conjunction with historical data. This revealed that pycnocline variability is one of the mechanisms of the eddy-current interaction process. The results show that a positive (negative) SLA corresponds to a smaller (larger) t 1,000 propagation time and thus to a deeper (shallower) depth of the pycnocline. During the impact between anticyclonic eddy and Kuroshio, the downward movement of the water at the inner The average temperature and salinity characteristics in the study area; DT is the temperature vertical gradient.
center of the eddy depresses the pycnocline, while the upwelling of the water at the center of the Kuroshio generates a convex pycnocline. The interacting process increases the isopycnic tilt, thus increasing the Kuroshio velocity. The influence of the cyclonic eddy on the Kuroshio reduces the tilt of the pycnocline, resulting in a reduction in the Kuroshio velocity or amplitude.
Based on the temperature, salinity, and current data from the CPIES inversion, the three-dimensional structure of the mesoscale eddies and Kuroshio was constructed for the first time from the eddy interaction processes captured by the large observation array, showing that the influence depth of the mesoscale eddies and Kuroshio can exceed 1,000 m. During the interaction between anticyclonic eddies and A B

FIGURE 8
Three-dimensional structure of the anticyclonic eddy. (A) Temperature anomaly and (B) salinity anomaly; the vector is the geostrophic current, and the color shading represents the temperature and salinity anomaly.

FIGURE 9
The temperature, salinity, and velocity along the 22.3°N section through the center of the anticyclonic eddy from (A, C). The color shading in (A-C) indicates temperature anomalies, salinity anomalies, and meridional velocity, respectively. The solid black lines in (A-C) indicate isothermal, isohaline, and isopycnic lines, respectively.
Kuroshio, the temperature in the eddy center shows a positive anomaly, and the salinity shows a positive-negative structure. The corresponding Kuroshio variability is opposite to that of the anticyclonic eddy and shows negative temperature anomalies; negative salinity anomalies are seen above 600 m and positive salinity anomalies are observed below 600 m due to the uplift of the deeper high-saline water.
When a cyclonic eddy impinges the Kuroshio, pycnocline variability is almost horizontally distributed. The structures of the temperature and salinity anomalies within the cyclonic eddy and Kuroshio are similar, with the temperature showing an overall negative anomaly and the salinity anomaly showing a negativepositive structure below the subsurface layer. Both temperature and A B

FIGURE 10
Three-dimensional structure of the cyclonic eddy. (A) Temperature anomaly and (B) salinity anomaly; the vector is the geostrophic current, and the color shading represents the temperature and salinity anomaly.

FIGURE 11
The temperature, salinity, and velocity along the 22.3°N section through the center of the cyclonic eddy from (A-C). The color shading in (A-C) indicates temperature anomalies, salinity anomalies, and meridional velocity, respectively. The solid black lines in (A-C) indicate isothermal, isohaline, and isopycnic lines, respectively. salinity changes are caused by the vertical movement of the water, such that the low-saline intermediate water above 600 m produces a negative salinity anomaly, while the upwelling of the deeper highsalinity water below 600 m causes a positive salinity anomaly. These results also demonstrate the feasibility of applying CPIES to the tracking and observing of mesoscale eddies, which is important for the study of heat and salt transport and for understanding how mesoscale eddies influence energy exchange.

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 authors.

Author contributions
QR conceived the study. FY managed the research. QR, YL and RW performed the data analysis. FN, JW, XD, ZC, CZ, RZ, HZ and X-HZ carried out raw data analysis and were in charge of the necessary resources. QR led the manuscript writing and produced the figures, with relevant input from all authors. All authors contributed to the article and approved the submitted version.

Funding
This research is supported by the National Key Research and Development Program of China (No. 2022YFC3104104), the National Natural Science Foundation of China (No. 42206032), and the Natural Science Foundation of Shandong Province (No. ZR2022QD045). All data and code used in the current study are available from the corresponding author upon reasonable request.