Multi-Sensor Observations of Submesoscale Eddies in Coastal Regions

: The temporal and spatial variation in submesoscale eddies in the coastal region of Lianyungang (China) is studied over a period of nearly two years with high-resolution (0.03 ◦ , about 3 km) observations of surface currents derived from high-frequency coastal radars (HFRs). The centers and boundaries of submesoscale eddies are identiﬁed based on a vector geometry (VG) method. A color index (CI) representing MODIS ocean color patterns with a resolution of 500 m is used to compute CI gradient parameters, from which submesoscale features are extracted using a modiﬁed eddy-extraction approach. The results show that surface currents derived from HFRs and the CI-derived gradient parameters have the ability to capture submesoscale processes (SPs). The typical radius of an eddy in this region is 2–4 km. Although no signiﬁcant di ﬀ erence in eddy properties is observed between the HFR-derived current ﬁelds and CI-derived gradient parameters, the CI-derived gradient parameters show more detailed eddy structures due to a higher resolution. In general, the HFR-derived current ﬁelds capture the eddy form, evolution and dissipation. Meanwhile, the CI-derived gradient parameters show more SPs and ﬁll a gap left by the HFR-derived currents. This study shows that the HFR and CI products have the ability to detect SPs in the ocean and contribute to SP analyses. values of the eddies’ kinematic parameters. It illustrates that the average vorticity of anticyclonic eddies is larger than that of cyclonic eddies; moreover, the instability of anticyclonic eddies is also derived from the CI images is slightly smaller than those found from the HFR velocity ﬁelds. Thus, an advantage of high-resolution data is that they can resolve the eddy structures in more detail. This study is the ﬁrst to extract submesoscale signals from MODIS CI data, with a resolution of 500 m, to investigate submesoscale coastal eddies. Comparing HFR surface currents with MODIS CI data, they can be mutually veriﬁed and they are complementary. Submesoscale eddies in CI images are detected near the location, where there are indeed some anticyclonic or cyclonic eddies from recent HFR records, revealing their consistency. So HFR data signiﬁcantly prove the ability of submesoscale signals extracted from CI-derived gradient parameters, while CI data enhance the spatial resolution and the scope of the observations. Based on the VG eddy detection algorithm, submesoscale eddies near the coastal region of the Huang Sea, close to the city of Lianyungang, China have been well analyzed. This methodology is potentially useful for investigations of ecosystem distributions in coastal biological communities in other areas of the global ocean.


Introduction
Submesoscale processes (SPs) are frequently observed as eddies, filaments and fronts, where the approximate scale ranges are 0.1-10 km in the horizontal, 0.01-1 km in the vertical, and between hours and days in time over open oceans [1]. The most fundamental characteristic difference between submesoscale and mesoscale is that the Rossby (Ro) and Richardson numbers (Ri) are both close to O(1) for SPs. SPs have large vertical velocities which are typically several times larger than those associated with mesoscale processes; thus, they play a critical role in vertical transport [2,3], enhance bio-physical interactions, and provide a possible link between mixing and dissipation. Submesoscale buoyancy flux in SPs enhances seasonal restratification and the mixing of the mixed layer [4]. Submesoscale  Dong et al. [37] used satellite-measured SST images to obtain the characteristics of ocean mesoscale eddies, which extracted most mesoscale signals limited by the coarse resolution of Remote Sensing System SST production. Here, Moderate Resolution Imaging Spectroradiometer (MODIS) ocean color patterns of resolution up to 500 m are applied to extract submesoscale signals, based on a similar approach proposed by Dong et al. [37]. A color index (CI, units: non-dimensional) represents the MODIS ocean color patterns, which are derived from an empirical approach proposed by Hu [38]. The empirical approach has been proven to validly remove sun glint and clouds and increase the data coverage. The CI has a significant correlation with the MODIS band-ratio Chl-a (< 1mg m −3 ) and is a good index to represent coastal color changes. Liu et al. [39] pointed out that mesoscale anticyclonic eddies in the Gulf of Mexico Loop Current system were accompanied by low CI values which indicated clear waters. Local CI increases (decreases) when submesoscale cyclonic (anticyclonic) eddies occur, while local SST decreases (increases) when submesoscale cyclonic (anticyclonic) eddies occur. This could help us understand the process by which eddies drive the upwelling of colder, nutrient rich waters that boost plankton production, the aggregation and retention of particles, etc. Because SST and CI have opposite responses to cyclonic and anticyclonic eddies, CI-derived gradient parameters are calculated by a modified method based on an automated approach to extract submesoscale signals. Zhang et al. [40] used two thresholds to determine the candidate boundaries for eddy detection. These two thresholds were determined from 625 MODIS CI images over the Florida Straits, which may not be applicable for other regions. Our modified method focuses on the gradient information only, without the need to use thermal wind approximation, such as in Dong et al. [37]. This is because the thermal wind approximation is not applicable for submesoscale eddies. In this study, the modified method of calculating gradient parameters goes beyond originally intended limits of the applicability of Dong et al. [37]. With the availability of HFR ground-truth data, this method is evaluated when scaled down to submesoscale. The gradient parameters are used to indicate the relative changes in CI values across eddy features; they contain meridional and zonal values to show eddy types (anticyclonic or cyclonic).
In this study, the VG method is modified to extract submesoscale signals from MODIS CI fields and to identify submesoscale eddies from HFR-derived currents fields and CI-derived gradient parameters. The primary goal of this paper is to analyze the coastal surface currents measured by HFRs with high-resolution data and MODIS ocean color patterns. The paper is organized as follows: in Section 2 we describe the dataset, including the HFR data and the high-resolution MODIS ocean color patterns, the eddy detection method, and the dynamics quantities. Section 3 presents results concerning our ability to extract submesoscale signals by CI-derived gradient parameters, and the eddy properties from HFR-derived current fields and MODIS CI-derived gradient parameters. Discussion is provided in Section 4. Final summary and conclusions are presented in Section 5.

HFR-Derived Currents
Surface currents were measured from June 2017 to March 2019 for the coastal area of Huang Sea, close to the city of Lianyungang, China, using two HFRs. The temporal sampling rate is about 20 min. The final radial current maps recorded simultaneously from two HFRs were synthesized to construct total vector maps on a Cartesian grid. The composite radar system mapped the Lianyungang coastal surface currents on 34 × 21 grids covering the surface area from 119.64 • to 120.63 • E and from 34.38 • to 34.98 • N, with a spacing of 0.03 • (about 3 km) in both directions. Figure 1b shows the observational range covered by the HFRs. The depth of the study area varied from about 1 to 33 m.
The numbers of observed data are shown in Figure 2. Part of the data are missing because of terrible weather and equipment failures. Although the temporal coverage rate is not always continuous, only a small number of submesoscale eddies are observed to survive for more than one day [13,22]. Therefore, the data are sufficient to show the entire submesoscale eddy processes. continuous, only a small number of submesoscale eddies are observed to survive for more than one day [13,22]. Therefore, the data are sufficient to show the entire submesoscale eddy processes.

MODIS Ocean Color Patterns
The MODIS ocean color patterns obtained from MODIS/Aqua satellite measurements over our study area (119.64°-120.63°E, 34.38°-34.98°N) have a spatial resolution of 0.004° × 0.004° and a daily temporal sampling rate (under cloud-free conditions) for the period from June 2017 to March 2019. Its original spatial resolution is about 500 m, but a 3 × 3 median filter was applied to remove speckle noise. Details about MODIS color index (CI) can be seen in Hu [38]. Briefly, CI is defined as the reflectance magnitude at 555 nm relative to two neighboring wavelengths (469 and 645 nm). For most ocean waters (chlorophyll concentration < 1 mg m -3 ) CI is highly correlated with MODIS chlorophyll data product (Chl-a) derived from the NASA default algorithms. For more productive waters, CI can be used as an index to show relative color patterns. The CI algorithm takes a similar band-subtraction approach to the OCI algorithm, but the former uses 469, 555, 645 nm band combinations while the latter uses 443, 547, 667 nm band combinations. The former bands do not saturate in severe sun glint, thick aerosols, and thin clouds, but the latter bands do saturate. The advantage of CI over MODIS Chl-a is that, because of the three land bands used to derive CI, CI does not saturate over sun glint, thin clouds, or thick aerosols, thereby providing more spatial and temporal coverage than MODIS chlorophyll to reveal the relative ocean color patterns. As shown in Figure 3, CI patterns from three subsequent images show the evolution of eddy patterns.
The CI dataset has advantages in research activities related to SPs because of its higher resolution and coverage compared to other MODIS SST or Chl-a datasets. The area surrounded by dotted lines in Figure 3 shows decreases in CI, corresponding to anticyclonic eddies shown in Figure 4.

MODIS Ocean Color Patterns
The MODIS ocean color patterns obtained from MODIS/Aqua satellite measurements over our study area (119.64 • -120.63 • E, 34.38 • -34.98 • N) have a spatial resolution of 0.004 • × 0.004 • and a daily temporal sampling rate (under cloud-free conditions) for the period from June 2017 to March 2019. Its original spatial resolution is about 500 m, but a 3 × 3 median filter was applied to remove speckle noise. Details about MODIS color index (CI) can be seen in Hu [38]. Briefly, CI is defined as the reflectance magnitude at 555 nm relative to two neighboring wavelengths (469 and 645 nm). For most ocean waters (chlorophyll concentration < 1 mg m -3 ) CI is highly correlated with MODIS chlorophyll data product (Chl-a) derived from the NASA default algorithms. For more productive waters, CI can be used as an index to show relative color patterns. The CI algorithm takes a similar band-subtraction approach to the OCI algorithm, but the former uses 469, 555, 645 nm band combinations while the latter uses 443, 547, 667 nm band combinations. The former bands do not saturate in severe sun glint, thick aerosols, and thin clouds, but the latter bands do saturate. The advantage of CI over MODIS Chl-a is that, because of the three land bands used to derive CI, CI does not saturate over sun glint, thin clouds, or thick aerosols, thereby providing more spatial and temporal coverage than MODIS chlorophyll to reveal the relative ocean color patterns. As shown in Figure 3, CI patterns from three subsequent images show the evolution of eddy patterns.
The CI dataset has advantages in research activities related to SPs because of its higher resolution and coverage compared to other MODIS SST or Chl-a datasets. The area surrounded by dotted lines in Figure 3 shows decreases in CI, corresponding to anticyclonic eddies shown in Figure 4. Remote Sens. 2019, 11, x FOR PEER REVIEW 6 of 17 MODIS ocean color patterns do not only focus on the local CI changes in time, but may also be used to estimate divergence and convergence through the CI-derived gradient fields. In our study, we define a modified method to extract CI gradient parameters from the MODIS CI data based on the method proposed by Dong et al. [37].
A CI gradient parameter = ( , ) is defined as

Diagnostic
The VG method for the eddy detection algorithm was applied to both the HFR-derived current fields and CI-derived gradient parameters to detect eddy structures. Several eddy properties were calculated after one eddy structure was identified. To understand the surface flow kinematics in the study region, various diagnostic values were calculated. The eddy kinetic energy (EKE) at each point was computed using the classical relation where and are the velocities derived with HFRs. In addition, the eddy vorticity( ), the shearing deformation rate ( ), the stretching rate ( ), the total deformation rate ( ) and the divergence ( ) were calculated using the following relations: MODIS ocean color patterns do not only focus on the local CI changes in time, but may also be used to estimate divergence and convergence through the CI-derived gradient fields. In our study, we define a modified method to extract CI gradient parameters from the MODIS CI data based on the method proposed by Dong et al. [37].
A CI gradient parameter V CI = U x , U y is defined as The meridional and zonal derivatives for each datapoint, G x and G y , are calculated from a Sobel gradient operator [41] in order to smooth the images as, follows where * denotes the 2-D convolution operation and A is a 2-D matrix of the CI data.
Since the rotation of cyclonic (anticyclonic) eddies is counterclockwise (clockwise) in the northern hemisphere, the CI values remain higher (lower) for cyclonic (anticyclonic) eddies. The parameter V CI U x , U y contains the gradient information of the CI data. Figure 5 shows examples of the gradient parameters derived from the CI data. The highest (lowest) CI values correspond to the centers of the cyclonic (anticyclonic) eddies.

Eddy Detection Algorithm
In this study, we have used the eddy detection algorithm based on the VG method proposed by Nencioli et al. [33]. This method is suitable for measuring the velocity field. It has been successfully applied to different study regions, such as the leeward side of Hawaii [42], the subtropical zonal band of the North Pacific Ocean [43], the Taiwan Strait [18], and the Kuroshio Extension Region [44]. We invite the reader to refer to Nencioli et al. [33] for more details.
There are two key steps in using the VG-based method to define an eddy structure: (a): In order to find the center of an eddy with four constraints, first, the u(v) component of the HF current velocity along the north-south (east-west) section must reverse sign across the eddy center, and its magnitude must increase away from the center, where the u and v components are the meridional and zonal components of the HF current velocity vector, respectively. Second, the center is where the velocity is at a local minimum. Finally, the velocity's change in direction must exhibit a constant sense of rotation and the directions of the two neighboring velocities must lie within the same quadrant of the two adjacent quadrants to the eddy center. The detection algorithm of eddy centers is applied to the HFR-derived, and CI-derived, velocity fields; (b): In order to obtain the eddy boundary, the outermost enclosed streamline is chosen as the eddy boundary around the center. The streamline at one position (i, j) is computed as where u and v are meridional and zonal components of the velocity. This approach is consistent with the criterion requiring that the values of the velocity components should remain increasing along the eddy boundary. In Figure 4, the blue dots and lines indicate, respectively, the detected centers and the boundaries of anticyclonic eddies. The radius is defined as the mean distance between the center and the perimeter of a given shape.

Diagnostic
The VG method for the eddy detection algorithm was applied to both the HFR-derived current fields and CI-derived gradient parameters to detect eddy structures. Several eddy properties were calculated after one eddy structure was identified. To understand the surface flow kinematics in the study region, various diagnostic values were calculated. The eddy kinetic energy (EKE) at each point was computed using the classical relation where u and v are the velocities derived with HFRs. In addition, the eddy vorticity(ζ), the shearing deformation rate (S s ), the stretching rate (S n ), the total deformation rate (S) and the divergence (ψ) were calculated using the following relations:

Features of Submesoscale Eddies (Eddy Properties) from HFR Observations
Submesoscale eddies were obtained from the HFR velocity fields, with a total of 1102 cyclonic eddies and 886 anticyclonic eddies identified, based on the VG eddy detection method. Submesoscale eddies cover most of the study region. Their number density is shown in Figure 6, which indicates that the occurrence of eddies mostly happens in the center of the study region and also in the southwestern coastal part. Moreover, the region where eddies frequently appear is off the cape at (39 • 40 N, 119 • 50 E) and the seamount in the center. In addition, a small number of eddies are detected in the study region as a result of being close to the outer edge of the HFR observation region. However, the identified eddies in the region do not show significant monthly variability due to differences in the validated data obtained each month (Figure 7).

Features of Submesoscale Eddies (Eddy Properties) from HFR Observations
Submesoscale eddies were obtained from the HFR velocity fields, with a total of 1102 cyclonic eddies and 886 anticyclonic eddies identified, based on the VG eddy detection method. Submesoscale eddies cover most of the study region. Their number density is shown in Figure 6, which indicates that the occurrence of eddies mostly happens in the center of the study region and also in the southwestern coastal part. Moreover, the region where eddies frequently appear is off the cape at (39°40′ N, 119°50′E) and the seamount in the center. In addition, a small number of eddies are detected in the study region as a result of being close to the outer edge of the HFR observation region. However, the identified eddies in the region do not show significant monthly variability due to differences in the validated data obtained each month (Figure 7).

Features of Submesoscale Eddies (Eddy Properties) from HFR Observations
Submesoscale eddies were obtained from the HFR velocity fields, with a total of 1102 cyclonic eddies and 886 anticyclonic eddies identified, based on the VG eddy detection method. Submesoscale eddies cover most of the study region. Their number density is shown in Figure 6, which indicates that the occurrence of eddies mostly happens in the center of the study region and also in the southwestern coastal part. Moreover, the region where eddies frequently appear is off the cape at (39°40′ N, 119°50′E) and the seamount in the center. In addition, a small number of eddies are detected in the study region as a result of being close to the outer edge of the HFR observation region. However, the identified eddies in the region do not show significant monthly variability due to differences in the validated data obtained each month (Figure 7).  The submesoscale characteristics have a radius of 1.3 to 8 km and an asymmetric normalized vorticity with a maximum between 0.3 and 1.5 (Figure 8). There is little difference between the size distribution of cyclonic and anticyclonic eddies. Figure 8 shows that eddies with a radius between 2 and 4 km are frequent, whereas eddies with a radius larger than 10 km are very infrequent. The mean radius values are 2.70 and 2.59 km for the anticyclonic and cyclonic eddies, respectively. Table 1 lists the mean values of the eddies' kinematic parameters. It illustrates that the average vorticity of anticyclonic eddies is larger than that of cyclonic eddies; moreover, the instability of anticyclonic eddies is also larger, according to the values of total deformation. Moreover, since surface current fields in this region are influenced by tides and topography, the vorticity of the surface current field is not large compared with western boundary region currents [45]. In addition, no significant difference is found for the radii of anticyclonic and cyclonic eddies.
The submesoscale characteristics have a radius of 1.3 to 8 km and an asymmetric normalized vorticity with a maximum between 0.3 and 1.5 (Figure 8). There is little difference between the size distribution of cyclonic and anticyclonic eddies. Figure 8 shows that eddies with a radius between 2 and 4 km are frequent, whereas eddies with a radius larger than 10 km are very infrequent. The mean radius values are 2.70 and 2.59 km for the anticyclonic and cyclonic eddies, respectively. Table 1 lists the mean values of the eddies' kinematic parameters. It illustrates that the average vorticity of anticyclonic eddies is larger than that of cyclonic eddies; moreover, the instability of anticyclonic eddies is also larger, according to the values of total deformation. Moreover, since surface current fields in this region are influenced by tides and topography, the vorticity of the surface current field is not large compared with western boundary region currents [45]. In addition, no significant difference is found for the radii of anticyclonic and cyclonic eddies.

Ability of Submesoscale Signals Extracted by CI-Derived Gradient Parameters
The CI-derived gradient parameters were calculated by the method described in Section 2. Figure 5 shows some examples of the gradient parameters at three different times. According to the location surrounded by the black dashed lines in Figure 3, the local CI is changed from large to small values. When anticyclonic eddies lasted for a short period, as shown in Figure 4, the anticyclonic eddies are found near the location of the black dashed lines in Figure 3; the size of anticyclonic eddies in the CI images is smaller than that from HFRs. In addition, a big difference is that there is one eddy that does exist in the center of HFR map and some other small eddies in this region. This illustrates that CI images could find the submesoscale signals due to physical processes, but they also contain other submesoscale information. The gradient parameters capture the submesoscale signals due to submesoscale anticyclonic eddies, and also other signals from other physical processes or biological processes. These processes mostly result from the sinking process accompanying anticyclonic eddies shown in Figure 3. The size and location of the anticyclonic eddies in CI gradient parameter fields have little difference from those in HFR-derived current fields because: (1) the Chl-a field aroused by the physical eddy field tends to be weaker than the current field and appears to be shifted by the force of the advection term; (2) the two datasets have different spatial resolutions and different properties (CI and velocity vectors). In addition, our method relies on a single image and gradient information, which is different from other studies about tracer-derived velocity products [46][47][48]. Nevertheless, some limitations do exist in both our method and the CI imagery when compared with HFR imagery. In fact, it is known in fluid mechanics that (a) submesoscale patterns could have taken a long time to emerge due to purely biological processes; (b) submesoscale patterns could be an example of "fossil turbulence" [49]. In both cases, the eddy patterns may only show up in the CI imagery as opposed to the HFR imagery. On the other hand, submesoscale eddies can exist without any color signature. Even with these known limitations, as long as eddy patterns are captured in the CI imagery, the CI-derived gradient parameters can successfully capture submesoscale signals.
Therefore, the CI-derived gradient parameters show similar eddy structures and provide even more submesoscale eddies. Although these gradient parameters do not provide realistic velocity fields, they are evidence of SPs and can also provide more detailed information about the eddy structures. This enhances the observation and analysis of submesoscale phenomena, implying that a CI (MODIS color pattern) with a resolution up to 500 m is a good tool to extract submesoscale signals from relatively larger observation areas, compared with HFR datasets.

Characteristics of Submesoscale Eddies (Eddy Properties) from CI-derived Gradient Fields
The results indicate that more submesoscale eddies can be detected with MODIS CI patterns with a resolution of 500 m than with the HFR data with a resolution of 3 km, although the former have a much lower temporal frequency (i.e., at most once per day versus every 20 min from HFR). A total of 4990 anticyclonic eddies and 4968 cyclonic eddies were identified with the abovementioned eddy detection algorithm from CI-derived gradient fields, which is more than four times those identified from the HFR velocity fields. The monthly distribution of the number of anticyclonic and cyclonic eddies is shown in Figure 9. We don't have sufficient CI data to report results during summer due to frequent clouds.
The CI-derived parameters (units: s −1 ) aren't the actual velocities, and so the vorticity of the CI parameters only focuses on its relative value. In Figure 10a, the mean radius values for submesoscale eddies are 1.85 and 1.56 km for the anticyclonic and cyclonic eddies, respectively. The mean radius values are smaller than those computed by the HFR-derived eddies. This indicates that eddy structures with a higher resolution are more suitable to detect SPs, and therefore better able to observe submesoscale signals. The strength of anticyclonic and cyclonic eddies is almost equal, as shown in Figure 10c. structures with a higher resolution are more suitable to detect SPs, and therefore better able to observe submesoscale signals. The strength of anticyclonic and cyclonic eddies is almost equal, as shown in Figure 10c.  structures with a higher resolution are more suitable to detect SPs, and therefore better able to observe submesoscale signals. The strength of anticyclonic and cyclonic eddies is almost equal, as shown in Figure 10c.

Discussion
The small sizes of the eddies detected from both HFR-derived currents and CI-derived gradient parameters are different from other coastal submesoscale eddies using HFRs [13,33]. Here, although eddies have high relative vorticities, they are weaker than the submesoscale island wake at Green Island in the Kuroshio, as simulated by Liu and Chang [45], who found that values for the eddy vorticity can reach up to 10 f (f is the Coriolis parameter), and the wake exhibits frequently detached recirculation, containing upwelling and von Kármán vortex streets. In our study region, cyclonic (anticyclonic) eddies were often accompanied by an increase (decrease) in CI (representing Chl-a and other in-water constituents). Jang et al. [16] reported similar studies about eddy properties from HFRs but didn't consider coastal responses of Chl-a because part of their study region did not have valid Geostationary Ocean Color Imager (GOCI) images. In this study, it is demonstrated that the efficient way to understand coastal SPs is via the application of HFR-derived currents and CI-derived gradient parameters.
The potential mechanism for the generation of submesoscale eddies in coastal areas can be inferred from geophysical factors, tides and wind stress. Firstly, a primary source for eddies' generation consists of geophysical factors, including coastline changes and bottom topography. Topographic features, such as islands, headlands and seamounts, generally exhibit a turbulent-drag interaction between currents and the topography. The formation and separation of bottom boundary layers affects the generation and injection of vertical vorticity. In this study area, multiple regions of eddies are mainly located along the coastline and near the seamount. The seamount changes the local vorticity and thus increases the probability of submesoscale eddies. The current fields near the coastline become unstable due to lateral boundaries, which also contribute to the generation of submesoscale eddies. Secondly, eddy generation is one possible cause of asymmetric advection by the rotating tidal flow [5]. The major axis of the M2 semidiurnal current is~29 cm/s near Lianyungang, which greatly influences the local currents. Thirdly, local eddies respond to the wind direction and the wind stress curl, thereby providing some indication of when and where eddies may be formed, evolve and ultimately dissipate [2]. Detailed discussions of these mechanisms will be presented in future papers, using numerical simulations.
In summary, although HFRs have some limitations due to their limited coastal coverage, they do meet the requirements for the spatio-temporal resolutions of SPs. Moreover, the gradient parameters derived from daily CI images have the ability to fill the gap left by HFRs, with the advantages of high space coverage. The latter method therefore shows potential for further research on SPs of global oceans.

Conclusions
In our study, CI-derived gradient parameters are calculated by a modified method from MODIS CI data. Two kinds of dataset, HFR-derived current fields and CI-derived gradient parameters, were applied to identify submesoscale eddies using the VG eddy detection algorithm, and thus to analyze the eddy properties. Both the HFR-derived velocity fields and the CI-derived gradient parameters are able to efficiently observe submesoscale eddies with data resolutions of 3 km and 500 m, respectively. The CI gradient parameters are able to capture more detailed submesoscale eddy features due to the higher spatial resolution of the CI images (500 m). However, the eddies' distribution and sizes are found to have some differences because of different sampling rates, observational variables, response delay, fossil turbulence, and possibly colorless eddies. In conclusion, the study reveals that HFR surface currents and MODIS CI are good options to investigate SPs.
In the coastal Lianyungang region, a total of 1102 cyclonic eddies and 886 anticyclonic eddies were recorded from the HFR-derived velocity fields, and a total of 4990 anticyclonic eddies and 4968 cyclonic eddies were recorded from the MODIS CI-derived gradient parameters, based on the VG eddy detection algorithm. The typical radius of eddies derived from the HFR velocity fields was 2.70 and 2.59 km for the anticyclonic and cyclonic eddies, respectively. The typical radius of eddies derived from the CI images is slightly smaller than those found from the HFR velocity fields. Thus, an advantage of high-resolution data is that they can resolve the eddy structures in more detail.
This study is the first to extract submesoscale signals from MODIS CI data, with a resolution of 500 m, to investigate submesoscale coastal eddies. Comparing HFR surface currents with MODIS CI data, they can be mutually verified and they are complementary. Submesoscale eddies in CI images are detected near the location, where there are indeed some anticyclonic or cyclonic eddies from recent HFR records, revealing their consistency. So HFR data significantly prove the ability of submesoscale signals extracted from CI-derived gradient parameters, while CI data enhance the spatial resolution and the scope of the observations. Based on the VG eddy detection algorithm, submesoscale eddies near the coastal region of the Huang Sea, close to the city of Lianyungang, China have been well analyzed. This methodology is potentially useful for investigations of ecosystem distributions in coastal biological communities in other areas of the global ocean.