Temporal Variation and Spatial Structure of the Kuroshio-Induced Submesoscale Island Vortices Observed from GCOM-C and Himawari-8 Data

: Dynamics of ocean current-induced island wake has been an important issue in global oceanography. Green Island, a small island located o ﬀ southeast of Taiwan on the Kuroshio path was selected as the study area to more understand the spatial structure and temporal variation of well-organized vortices formed by the interaction between the Kuroshio and the island. Sea surface temperature (SST) and chlorophyll-a (Chl-a) concentration data derived from the Himawari-8 satellite and the second generation global imager (SGLI) of global change observation mission (GCOM-C) were used in this study. The spatial SST and Chl-a variations in designed observation lines and the cooling zone transitions on the left and right sides of the vortices were investigated using 250 m spatial resolution GCOM-C data. The Massachusetts Institute of Technology general circulation model (MITgcm) simulation conﬁrmed that the positive and negative vortices were sequentially detached from each other in a few hours. In addition, totals of 101 vortexes from July 2015 to December 2019 were calculated from the 1-h temporal resolution Himawari-8 imagery. The average vortex propagation speed was 0.95 m / s. Totals of 38 cases of two continuous vortices suggested that the average vortex shedding period is 14.8 h with 1.15 m / s of the average incoming surface current speed of Green Island, and the results agreed to the ideal Strouhal-Reynolds number ﬁtting curve relation. Combined with the satellite observation and numerical model simulation, this study demonstrates the structure of the wake area could change quickly, and the water may mix in di ﬀ erent vorticity states for each observation station.


Background
Island wakes have been studied for many years in the global oceans. Different driving forces cause this phenomenon and these wakes have different characteristics according to the island's scale and water depth [1][2][3][4]. Green Island, a small island located off southeast of Taiwan, is on the path of the Kuroshio to the East China Sea (Figure 1). Green Island is an obstacle on the "ocean highway". Kuroshio passes Green Island and causes a change in flow fields, such as the Von Kármán vortex street [5][6][7][8][9]. In the lee of the island, a high chlorophyll-a (Chl-a) concentration and low sea surface temperature (SST) wake region can be formed [8]. This recirculation area contained a cyclonic/anticyclonic vorticity pair accompanied by a density overturn and water upwelling [9]. The increased surface Chl-a concentration and decreased SST are induced by upwelling from deep-layer waters in the island wakes [5,9]. Like a supply depot for marine life, the SST fronts between the wake area and the Kuroshio could provide Previous researches have used high-resolution satellite imagery, acoustic Doppler current profilers, cruise observations, and numerical models to analyze the Green Island wakes [5,[8][9]]. An island-induced ocean vortex train (IOVT) [7] is formed leeward of Green Island by the incoming Kuroshio [8]. When an ocean vortex forms, observations from the moored acoustic Doppler current profiler (ADCP) suggest that the velocity on the western side of Green Island increases [8]. In addition, the IOVT shows variability under wind forcing [8,13]. Southerly wind helps the wake expand because it lies in the same direction as the incoming current, while the northerly wind compresses the development of the wake [8,13]. The island wakes influence the area of SST drops, which could extend 35 km downstream of Green Island, and vortex shedding could propagate as far as 80 km downstream of the island [9]. Research on underwater island wakes suggests that a depth range with a high inverted Reynolds number is supposed to have density overturns and a high Thorpe scale with a turbulent kinetic energy dissipation rate of O(10 -6 -10 -5 )W kg ⁄ , which corresponds to an eddy diffusivity of O(10 -2 -10 -1 )m 2 s ⁄ [9]. Numerical model research has suggested that the vortex street features are adapted by inertial and barotropic instabilities [14], and the shedding period of the vortex is synchronized to a tidal period [15].

Objectives
Although there are many studies on the Green Island wakes, the dynamics of these wakes remain difficult to be determined. Cruise and instrument observations can measure detailed data on the underwater wakes, but the advantages of remote sensing could be used to obtain accurate and simultaneous data for sea surface wakes. Before the next measurement, it is necessary to clarify the spatial scale and dynamic processes of the wakes on the island surface using satellite data. Only by determining the feedback of the surface seawater to the wake can a researcher accurately design the observation line and measurement frequency by cruise. Previous researches have used high-resolution satellite imagery, acoustic Doppler current profilers, cruise observations, and numerical models to analyze the Green Island wakes [5,8,9]. An island-induced ocean vortex train (IOVT) [7] is formed leeward of Green Island by the incoming Kuroshio [8]. When an ocean vortex forms, observations from the moored acoustic Doppler current profiler (ADCP) suggest that the velocity on the western side of Green Island increases [8]. In addition, the IOVT shows variability under wind forcing [8,13]. Southerly wind helps the wake expand because it lies in the same direction as the incoming current, while the northerly wind compresses the development of the wake [8,13]. The island wakes influence the area of SST drops, which could extend 35 km downstream of Green Island, and vortex shedding could propagate as far as 80 km downstream of the island [9]. Research on underwater island wakes suggests that a depth range with a high inverted Reynolds number is supposed to have density overturns and a high Thorpe scale with a turbulent kinetic energy dissipation rate of O(10 −6 -10 −5 )W/kg, which corresponds to an eddy diffusivity of O(10 −2 -10 −1 )m 2 /s [9]. Numerical model research has suggested that the vortex street features are adapted by inertial and barotropic instabilities [14], and the shedding period of the vortex is synchronized to a tidal period [15].

Objectives
Although there are many studies on the Green Island wakes, the dynamics of these wakes remain difficult to be determined. Cruise and instrument observations can measure detailed data on the underwater wakes, but the advantages of remote sensing could be used to obtain accurate and simultaneous data for sea surface wakes. Before the next measurement, it is necessary to clarify the spatial scale and dynamic processes of the wakes on the island surface using satellite data. Only by determining the feedback of the surface seawater to the wake can a researcher accurately design the observation line and measurement frequency by cruise. The purpose of this study is to use the satellite data to obtain novel results of the Green Island wake in order to make up for the shortcomings of in-situ measurements and the inability for the model to be measured. They include (1) designing the imaginary lines to describe the Chl-a and SST response to the wake from the second generation global imager (SGLI) data on the global change observation mission (GCOM-C), (2) using the Himawari-8 SST data to calculate the wake vortex detachment propagating speed and using the maximum cross-correlation (MCC) method to estimate the incoming velocity of the wake, (3) calculating and discussing the Strouhal number-Reynolds number relationship of wake dynamics.
Understanding the spatial structure of surface wakes can help determine the biological hotspots in the wake area and improve knowledge of oceanic wakes. The evolution process of a wake can be separated into four types, including a wake occurring alone, a wake occurring with a tail stretching downstream, an S-shaped meandering wake, and a wake with a small cyclonic/anticyclonic vortex pair downstream [9]. However, the detailed spatial structures of each type of wake remain unknown. The SGLI observations with 250 m spatial resolution will improve our understanding of ocean change mechanisms through long-term monitoring. The use of imaginary design observation lines helps to understand the sea surface characteristics of the wake region at the same time but different locations.
Using high spatial resolution sun-synchronous satellite data could successfully detect the SST drop in the cold wake area but could not further track the process of vortex movement in continuous an hour interval due to the lack of a high temporal resolution [9]. To study the dynamics of island wakes, the most important factor is to solve the setting of the two major parameters, the Strouhal number (St) and the Reynolds number (Re), as well as the propagation velocity of the cold vortex (U e ). However, it is very difficult to measure these two parameters with a cruise observation in the ocean. For the Green Island area, there must be a certain number of research surveys to measure the current velocity and SST variations in front of and downstream of the island. Previous studies applied parameters such as the estimated distance between two consecutive vortices and the radius of vortices from Synthetic Aperture Radar (SAR) images to the empirical equation to obtain the incoming current velocity and the propagation velocity of the vortex [16]. However, SAR imagery is a snapshot image, which does not have a high temporal resolution or a sufficient image, making it difficult to engage in a long-term observation of wake dynamics. With advances in satellite imagery, the Himawari-8 geostationary satellite data provide a useful solution to the small cold vortex detachment process and detect the vortex trajectory in wake dynamics. Understanding the relationship between St and Re and the speed and period of cold vortex detachment can be useful for improving numerical simulations.
Recognizing the Green Island wake not only requires us to physical discuss wake dynamics but also to investigate the development of ocean environmental sustainability. The results could help marine fisheries and chemists better understand the status of this area.

Materials and Methods
The process of this study was divided into three steps. Firstly, the in-situ observations were used to present the underwater structure in the wake region. Secondly, the detailed structure of the surface wake was analyzed using the SGLI data from the GCOM project. Finally, the Himawari-8 imagery was used to analyze the vortices movement process, and the MITgcm (Massachusetts Institute of Technology General Circulation Model) numerical model was used to simulate the wake variation.
The MCC method was applied to hourly Himawari-8 SST images to calculate the incoming surface current velocity for Green Island (22.5 • -22.7 • N, 122.3 • -122.5 • E). The principle of the MCC method is to find the maximum correlation between a template image of SST patterns and the selected search image in the next time interval. In this way, we can determine how far the center point of the template image has moved in an hour and estimate the flow velocity. The MCC method needs to adjust the parameter selection settings for different research scopes [17][18][19]. In this study, the pixel resolution of 2 km was increased to 1 km by interpolation. Considering the complex flow field variations and spatial scale of the SST variation in southeast Taiwan, the four template images (7 × 7, 9 × 9, 11 × 11, and Remote Sens. 2020, 12, 883 4 of 18 13 × 13 pixels) and the sides of the search image (7 km to each direction from the center of the template image) were chosen. In this study, the acceptance criteria for the cross-correlation coefficient should be larger than 0.9. With this parameter set, the minimum and the maximum speeds in the estimate were about 0.28 m/s and 1.94 m/s, respectively. The maximum flow velocity per template image was selected as the flow velocity in this area. Finally, the average maximum speed of the four template images during each vortex shedding case was calculated.

Sea Surface Temperature and Chlorophyll-a Concentration
The GCOM-C, carrying an SGLI conducts surface and atmospheric measurements, such as measurements of clouds, aerosols, ocean color, vegetation, and snow and ice. The SGLI is an optical sensor capable of multi-channel observations at wavelengths from near-UV to thermal infrared and obtains global observation data once every 2 or 3 days. The SGLI data are provided by the Japan Aerospace Exploration Agency/National Aeronautics and Space Administration. The data are available from January 2018, with a 250 m spatial resolution. The standard product and algorithm handbooks can be found on https://suzaku.eorc.jaxa.jp/GCOM_C. Himawari-8 is the 8th of the Himawari geostationary Japanese weather satellites operated by the Japan Meteorological Agency. Himawari-8 carries the Advanced Himawari Imager with a wide spectral range and very high spatial and 10-minute temporal resolutions. The research products of SST and Chl-a used in this paper were supplied by the P-Tree System, Japan Aerospace Exploration Agency (JAXA). The standard product can be found on https://www.eorc.jaxa.jp/ptree/index.html. The SST data (2 km spatial resolution) and Chl-a (5 km spatial resolution) are available from 7 July 2015, with a 1-h temporal resolution. The Himawari-8 data have good accuracy and have been applied to calculate short-term sea surface currents [17,20].

Ocean Currents
OSCAR (Ocean Surface Current Analysis Real-time) contains near-surface ocean current estimates derived using quasi-linear and steady flow momentum equations and combining geostrophic, Ekman, and Stommel shear dynamics, and a complementary term from the surface buoyancy gradient, detailed calculations can be referred to the User's Handbook [21]. These data were collected from various satellites and in situ instruments and directly estimated from the sea surface height, surface vector wind, and SST. These data were generated by the Earth Space Research (ESR) and are available from October 1992 to present, with a 1/3 degree spatial resolution and a 5-day temporal resolution. The OSCAR third-degree resolution ocean surface currents (OSCAR_L4_OC_third-deg) data which can be accessed through the https://podaac.jpl.nasa.gov/dataset/OSCAR_L4_OC_third-deg from NASA PODAAC (Physical Oceanography Distributed Active Archive Center). The OSCAR data are highly associated and accurate with the global tropical moored buoy array, and it has also been successfully used to explore the island wake around Palau in the western tropical North Pacific Ocean [1,22,23].

Numerical Model
MITgcm is a numerical model designed for the study of the atmosphere, ocean, and climate. The source code for this model can be download from http://mitgcm.org/. We used this model to solve the Navier-Stokes equations under hydrostatic and Boussinesq approximations using the finite volume method [24]. The south, north, and east directions of the model domain are open boundaries that use the Orlanski radiation condition [25]. The HYbrid Coordinate Ocean Model and the Navy Coupled Ocean Data Assimilation (HYCOM + NCODA) Global 1/12 degree analysis (GLBu0.08) was used for the initial background data and the boundary driving force in the model. In this study, the case of the summer season was based on the average of June to August from 2013 to 2018. These data can be download from https://www.hycom.org/. The study area in the model set is from 120.85 • E to 122.5 • E and from 22.25 • N to 23.5 • N. The horizontal resolution is 500 m. The model uses z-coordinates. There are 80 levels in the vertical direction. The top layer's thickness is 5 m, and the other thicknesses increase at a 5.26% increase rate for each layer. The maximum depth is 5513 m. The horizontal eddy viscosity is 40 m 2 /s. The Laplacian diffusion of heat and salt laterally is 4 m 2 /s. The vertical diffusion of temperature and salt depends on the buoyancy frequency [26].

In-Situ Observation
In-situ observations were conducted in November 2012 using R/V Ocean Researcher I (OR1) in the lee of Green Island (the datasets are available from http://www.odb.ntu.edu.tw/en/). Seven stations (A1-A7) were selected to detect the vertical layer of the island wake from the southeast coast of Taiwan across the Kuroshio and the wake area north of Green Island (see Figure 2a). The shipboard acoustic Doppler current profiler (Sb-ADCP) and Conductivity Temperature Depth (CTD) surveys were also collected during the cruise. The AquaTracka III fluorometer is used for in-situ detection of Chl-a, dissolved compounds re-emit a fraction of this energy as fluorescence at longer wavelengths when they absorb light and the intensity of fluorescence is directly proportional to the concentration. The measurement range of the 75 kHz Sb-ADCP on board R/V OR1 was 16.56-650 m, with a bin size of 8 m. It recorded the current velocity with 2 min ensembles averaged over 30 pings.
Remote Sens. 2020, 12, x FOR PEER REVIEW 5 of 18 horizontal eddy viscosity is 40 m 2 s ⁄ . The Laplacian diffusion of heat and salt laterally is 4 m 2 s ⁄ . The vertical diffusion of temperature and salt depends on the buoyancy frequency [26].

In-Situ Observation
In-situ observations were conducted in November 2012 using R/V Ocean Researcher I (OR1) in the lee of Green Island (the datasets are available from http://www.odb.ntu.edu.tw/en/). Seven stations (A1-A7) were selected to detect the vertical layer of the island wake from the southeast coast of Taiwan across the Kuroshio and the wake area north of Green Island (see Figure 2a). The shipboard acoustic Doppler current profiler (Sb-ADCP) and Conductivity Temperature Depth (CTD) surveys were also collected during the cruise. The AquaTracka III fluorometer is used for in-situ detection of Chl-a, dissolved compounds re-emit a fraction of this energy as fluorescence at longer wavelengths when they absorb light and the intensity of fluorescence is directly proportional to the concentration. The measurement range of the 75 kHz Sb-ADCP on board R/V OR1 was 16.56-650 m, with a bin size of 8 m. It recorded the current velocity with 2 min ensembles averaged over 30 pings.

Field Experiment
In this field experiment, the area between A4 and A5 was significantly affected by the island wake (Figure 2b,c), the isotherms and isohalines were uplifted from a depth of 100-120 m to the surface and the high Chl-a concentration water also upwelled from a depth of 70 m to the surface (Figure 2d-f). The wake water has lower SST, higher salinity, and higher Chl-a than the Kuroshio water. Moreover, the cruise observations could not measure the whole area of the wake at the same time, and the hydrological structure may be changed due to the dynamic nature of the island wakes. Therefore, remote sensing data with high temporal and spatial resolutions are recommended to discuss the surface structures of wakes.

Spatial Structure of Island Wake
In this section, GCOM-C SGLI data with a 250 m grid resolution are used to clarify these properties. Notably, these high-resolution cloud-free images are used to observe the spatial changes of SST and Chl-a at different locations throughout the wake region at the same time, and we designed imaginary observation lines for each case. In case 1, the classic and most common wake type with ten east-west observation lines is used to discuss the wake from the lee of Green Island to downstream (Figures 3 and 4). In case 2, an S-shaped meandering wake with three lines along the wake and one imaginary line are used to discuss the spatial-structural differences between two cold vortices ( Figure 5). Case 3 offers a good comparative example to show the change of vorticity in the wake region during the vortex transition ( Figure 6). A numerical model simulation is used to illustrate this phenomenon ( Figure 7). Case 4 presents a typical theoretical von Kármán vortex street, but such phenomena rarely occur in an SST response ( Figure 8).

Field Experiment
In this field experiment, the area between A4 and A5 was significantly affected by the island wake (Figure 2b,c), the isotherms and isohalines were uplifted from a depth of 100-120 m to the surface and the high Chl-a concentration water also upwelled from a depth of 70 m to the surface (Figure 2d-f). The wake water has lower SST, higher salinity, and higher Chl-a than the Kuroshio water. Moreover, the cruise observations could not measure the whole area of the wake at the same time, and the hydrological structure may be changed due to the dynamic nature of the island wakes. Therefore, remote sensing data with high temporal and spatial resolutions are recommended to discuss the surface structures of wakes.

Spatial Structure of Island Wake
In this section, GCOM-C SGLI data with a 250 m grid resolution are used to clarify these properties. Notably, these high-resolution cloud-free images are used to observe the spatial changes of SST and Chl-a at different locations throughout the wake region at the same time, and we designed imaginary observation lines for each case. In case 1, the classic and most common wake type with ten east-west observation lines is used to discuss the wake from the lee of Green Island to downstream (Figures 3 and 4). In case 2, an S-shaped meandering wake with three lines along the wake and one imaginary line are used to discuss the spatial-structural differences between two cold vortices ( Figure  5). Case 3 offers a good comparative example to show the change of vorticity in the wake region during the vortex transition ( Figure 6). A numerical model simulation is used to illustrate this phenomenon ( Figure 7). Case 4 presents a typical theoretical von Kármán vortex street, but such phenomena rarely occur in an SST response ( Figure 8).             The typical Green Island wake (Figure 3) presents with a cold wake extending about 25 km behind the island. We assumed that the 10 lines cross the wake (L1 to L10 in Figure 3e). If a cruise survey is performed, only one of these lines could be observed in one hour due to the speed limit of the voyage. The lowest SST and the highest Chl-a concentration in the wake region were found at the imaginary lines of L1 to L4. L1 is about 1 km behind the island, and the SST of the wake area (27.3 °C) can reach 2.1 °C lower than that of the Kuroshio water (29.4 °C) as seen in Figure 4. The Chl-a concentration reached 0.23 mg m 3 ⁄ , which was three times that of the Kuroshio water (0.07 mg m 3 ⁄ ). We can use this case to quantify the magnitude of the horizontal gradient of SST and Chl-a. From L1 to L4, the SST gradient is about 0.5 °C/km, and the Chl-a concentration is about 0.025 mgm -3 km ⁄ . Notably, the dynamic process of a wake can change rapidly, and the seawater structure of the entire wake area will likely change drastically within one hour. It can be seen here that from L7 to L10, there is a small vortex that may detach soon and propagate downstream. Around 23°-23.2° N, 121.7°-121.8° E, two small vortices were captured (Figure 3c,f). However, we are unsure whether this was caused by the shedding of the vortex train or by the coastal current. The typical Green Island wake (Figure 3) presents with a cold wake extending about 25 km behind the island. We assumed that the 10 lines cross the wake (L1 to L10 in Figure 3e). If a cruise survey is performed, only one of these lines could be observed in one hour due to the speed limit of the voyage. The lowest SST and the highest Chl-a concentration in the wake region were found at the imaginary lines of L1 to L4. L1 is about 1 km behind the island, and the SST of the wake area (27.3 • C) can reach 2.1 • C lower than that of the Kuroshio water (29.4 • C) as seen in Figure 4. The Chl-a concentration reached 0.23 mg/m 3 , which was three times that of the Kuroshio water (0.07 mg/m 3 ). We can use this case to quantify the magnitude of the horizontal gradient of SST and Chl-a. From L1 to L4, the SST gradient is about 0.5 • C/km, and the Chl-a concentration is about 0.025 mgm −3 /km. Notably, the dynamic process of a wake can change rapidly, and the seawater structure of the entire wake area will likely change drastically within one hour. It can be seen here that from L7 to L10, there is a small vortex that may detach soon and propagate downstream. Around 23 • -23.2 • N, 121.7 • -121.8 • E, two small vortices were captured (Figure 3c,f). However, we are unsure whether this was caused by the shedding of the vortex train or by the coastal current.
Unlike the low-SST wake in Figure 3, which was concentrated only in a certain range behind the island, Figure 5 shows a case of a vortex train, which the low SST and high Chl-a concentration continue to extend downstream. The imagery line L1 in Figure 5b has the same purpose as L1 in Figure 3e. In this case, a Chl-a concentration of 0.66 mg/m 3 was observed, which was eight times the Kuroshio water, and the difference in SST was 2 • C. The distance between the two vortices is about 30 km (L2 to L4 in Figure 5e). L3 and L4 lines in Figure 5b passed through two vortices. The data sequence shows the situation of a low-high-low SST and high-low-high Chl-a concentration. Since the downstream vortex is affected by seawater mixing, the SST is usually higher than that of the wake close to the island.
A classic state in wake dynamics is illustrated in Figure 6, which was not captured in previous studies on Green Island. With high-resolution satellite data, we can observe the cooling zone transitions on the left and right side of the island wake region (usually corresponding to the magnitude of the vorticity). From the data in Figure 6e, it can be seen that there are two peaks in the low-SST wake region, 28.8 • C and 28.9 • C, respectively, and the corresponding Chl-a concentrations are 0.24 mg/m 3 and 0.20 mg/m 3 . The left side of the wake region has a lower SST and a higher Chl-a concentration. From the data in Figure 6f, two peaks in the low-SST wake region can be found (28.3 • C and 27.5 • C), and the corresponding Chl-a concentrations are 0.14 mg/m 3 and 0.16 mg/m 3 . The right side of the wake region has lower SST and higher Chl-a concentration. To understand the wake evolution more clearly, we used the MITgcm model to present the variations in Rossby number (Ro) during the wake's dynamic process (Figure 7): where U and L are the velocity and length scales of the phenomenon, f is the Coriolis frequency, and ζ is the relative vorticity. The quantified vorticity using Ro helps to interpret the eddy strength.
A very small Ro is typical to general circulation. A very high Ro can be seen in rapidly rotating eddies. The simulation results show that positive and negative vorticity are sequentially detached from each other in a few hours. This clearly illustrates the results we captured from the GCOM-C images ( Figure 6). This also demonstrates the importance of being careful when making cruise measurements, because each wake can change quickly, and the water may mix in different vorticity states for each observation station. An ideal von Carmen vortex street can be seen in Figure 8a,b. Four sub-mesoscale vortices detached from the island to downstream. Vortex 3 (Figure 8c,d) and vortex 4 (Figure 8e,f) were selected for observations, and two imaginary lines that crossed the highest Chl-a position of the vortex were designed. The size of the vortex's core is about 5 km, which is close to the size of Green Island. The maximum Chl-a concentration in both vortices was 0.22 mg/m 3 . An interesting phenomenon was also seen here. In the vortex core, the Chl-a concentration did not exactly match that of the SST. As a rotating vortex, the SST may increase during the propagation process due to the mixing of the Kuroshio. The SST variation of the vortex in the model experiment is difficult to simulate because it is not known which factors alongside seawater mixing could affect the vortex when it detaches downstream. In the satellite observations, it was often found that vortex street rarely appeared, which was very different from the model simulations [9]. In the next section, the geosynchronous Himawari-8 satellite data were used to analyze and discuss the dynamic process of the vortex train.

Temporal Variation and Vortex Trajectory
The previous section presented the spatial structure of biological hotspots in the wake area, and the dynamic process of the vortex is mentioned in this section. To study the dynamics of island wakes, the most important factor is to solve the setting of the two major parameters, the Strouhal number (St) and the Reynolds number (Re), as well as the propagation velocity of the vortex (U e ). The St is a dimensionless number describing the oscillating flow mechanisms as: where T is the vortex shedding period, U 0 is the incoming current speed, and L is the characteristic length. The Re is an important dimensionless quantity in fluid mechanics used to help predict the flow patterns in different fluid flow situations as: where ν h is the horizontal eddy viscosity.
With continuous cloud-free data, the trajectory and velocity of the vortex could be tracked and calculated. A vortex has the characteristics of a closed isotherm of SST, and the lowest SST in the eddy can represent its center position [9]. In this way, the trajectory of the vortex can be tracked. An example of two consecutive vortices within 24 h of SST images is shown in Figure 9, with the temporal ranging from 21:00 (UTC) on 12 July 2016 to 20:00 (UTC) on 13 July 2016. We can see how the wake vortex detached and passed downstream. The first vortex was formed at 23:00 (UTC) on 12 July 2016 (the red star in Figure 9), and the second vortex was formed at 11:00 (UTC) on 13 July 2016 (the red dot in Figure 9 Figure 10c presents a histogram of property velocity statistics for 101 vortex cases. As mentioned in previous studies [8], the angle of the vortex distribution and the angle of the incoming velocity in front of the island are highly correlated. Therefore, based on the trajectory distribution in Figure 10a, the angle of the incoming current field in front of the island can be estimated. More than half of the vortex cases have propagation velocities between 0.8 and 1.2 m/s. The probability distribution of the trajectory is shown in Figure 10b. Most of the vortices with low SST and high Chl-a concentrations are concentrated within 30 km northeast of Green Island. The biological hotspots in the wake area mainly appear in this area of 900 km 2 and propagate downstream with the vortex stream. Based on the MCC method, the calculation results show that the average incoming flow velocity U 0 that caused the vortex shedding cases in Figure 10a   Next, we use the case of vortex shedding to calculate St, and we must determine the time interval between two continuous vortices. Therefore, a total of 38 cases were provided in this study for calculations. The results show that the average vortex shedding period is 14.8 hours. Based on the settings of L = 5500 m and ν h = 100 m 2 s ⁄ , the average St and Re are 0.114 and 64, respectively, which are consistent with the fitted curve values [27] as shown in Figure 11.  Next, we use the case of vortex shedding to calculate St, and we must determine the time interval between two continuous vortices. Therefore, a total of 38 cases were provided in this study for calculations. The results show that the average vortex shedding period is 14.8 hours. Based on the settings of L = 5500 m and ν h = 100 m 2 s ⁄ , the average St and Re are 0.114 and 64, respectively, which are consistent with the fitted curve values [27] as shown in Figure 11. Next, we use the case of vortex shedding to calculate St, and we must determine the time interval between two continuous vortices. Therefore, a total of 38 cases were provided in this study for calculations. The results show that the average vortex shedding period is 14.8 h. Based on the settings of L = 5500 m and ν h = 100 m 2 /s, the average St and Re are 0.114 and 64, respectively, which are consistent with the fitted curve values [27] as shown in Figure 11. Remote Sens. 2020, 12, x FOR PEER REVIEW 13 of 18 Figure 11. The Strouhal number (St) versus the Reynolds number (Re) diagram. The point for this study is expressed as the mean value with one standard deviation.

The Relationship Between Re and St
The relationship between Re and St plays a major role in determining island wake patterns, evolution, and vortex shedding. The fitting curve relation between Re and St for cylinder wakes were formulated as [27]: (5) In theory, the Re value depends only on the size of U 0 when focusing on the same obstacle. However, Green Island is an irregular island, and there were differences in the choice of the diameter (L) of the obstacle. According to the measurements from the National Land Surveying and Mapping Center, Taiwan, the size of Green Island is 15.09 km 2 , 4.6 km wide, 5.2 km long, and the longest side is about 6 km. Therefore, it is appropriate to choose L = 5-6 km for research on Green Island. However, some researchers have instead chosen L = 7 km, which indicates a width of 100 m for the island in the water [5,14,15]. It is worth noting that the choice of L will significantly affect the values of St and Re. In addition, the value selection of horizontal eddy viscosity (ν h ) when calculating the Re value is a difficult problem to be solved. Since this value is not easy to measure with observations, previous studies have used estimated values. In the open ocean, the ν h varies from 10 2 to 10 5 m 2 s ⁄ P [28]. For the mesoscale or sub-mesoscale ocean processes, the ν h is generally chosen as 10 3 to 10 4 m 2 s ⁄ P [16]. Next, we sort the data of various researchers and recommend what types of values are most suitable for studying Green Island. Table 1 shows the corresponding data for each point in Figure 11. In this study, L and ν h were selected as 5500 m and 100 m 2 s ⁄ , respectively. The results of 38 cases showed that the cases fit the equation of the fitted curve. With the average U 0 = 1.15 m/s, the average Re is 64, and the average St is 0. 114. We recommend that L be selected from 5-6 km. If a larger value of L is selected, the St may be too large (>0.3). A ν h value of 100 m 2 s ⁄ P is suitable, as a ν h value that is too small may cause the Re value to be too large. Although previous studies were based on observational data, there was only one sample. Take the case of this study as an example. If the ν h value was chosen to be 15 m 2 s ⁄ P [14] or even 0.2-7 m 2 s ⁄ P [15], the points of the case could be far from the fitted curve equation.

The Relationship Between Re and St
The relationship between Re and St plays a major role in determining island wake patterns, evolution, and vortex shedding. The fitting curve relation between Re and St for cylinder wakes were formulated as [27]: In theory, the Re value depends only on the size of U 0 when focusing on the same obstacle. However, Green Island is an irregular island, and there were differences in the choice of the diameter (L) of the obstacle. According to the measurements from the National Land Surveying and Mapping Center, Taiwan, the size of Green Island is 15.09 km 2 , 4.6 km wide, 5.2 km long, and the longest side is about 6 km. Therefore, it is appropriate to choose L = 5-6 km for research on Green Island. However, some researchers have instead chosen L = 7 km, which indicates a width of 100 m for the island in the water [5,14,15]. It is worth noting that the choice of L will significantly affect the values of St and Re. In addition, the value selection of horizontal eddy viscosity (ν h ) when calculating the Re value is a difficult problem to be solved. Since this value is not easy to measure with observations, previous studies have used estimated values. In the open ocean, the ν h varies from 10 2 to 10 5 m 2 /s [28]. For the mesoscale or sub-mesoscale ocean processes, the ν h is generally chosen as 10 3 to 10 4 m 2 /s [16]. Next, we sort the data of various researchers and recommend what types of values are most suitable for studying Green Island. Table 1 shows the corresponding data for each point in Figure 11. In this study, L and ν h were selected as 5500 m and 100 m 2 /s, respectively. The results of 38 cases showed that the cases fit the equation of the fitted curve. With the average U 0 = 1.15 m/s, the average Re is 64, and the average St is 0. 114. We recommend that L be selected from 5-6 km. If a larger value of L is selected, the St may be too large (>0.3). A ν h value of 100 m 2 /s is suitable, as a ν h value that is too small may cause the Re value to be too large. Although previous studies were based on observational data, there was only one sample. Take the case of this study as an example. If the ν h value was chosen to be 15 m 2 /s [14] or even 0.2-7 m 2 /s [15], the points of the case could be far from the fitted curve equation. Table 1. Corresponding values for each point in Figure 10.

L (m)
U 0 (m/s) ν h (m 2 /s) Re Next, we will discuss the cases in our sample with Re > 40 (vortex street occurrence) but St < 0.1095 (the minimum value of the fitted curve equation). According to our parameter settings, if Re is greater than 40, U 0 must be larger than 0.73 m/s. The minimum U 0 in 38 cases is 0.85 m/s. To make St > 0.1095, when U 0 = 1 m/s, T < 13.95 h, when U 0 = 1.2 m/s, T < 11.62 h, when U 0 = 1.4 m/s, T < 9.96 h. According to the field observation data [14], the vortex shedding period is 12 h. However, there are many cases in this study that showed a longer period for vortex shedding, which results in a lower St. According to our observations, a low SST vortex formed over a long period of time, but it did not detach for a long time. This verifies the results of previous studies that have used about 15 years of MODIS SST data to determine that about 87% of cases feature wakes occurring alone [9]. However, the cause of this situation is still unknown. It may be due to local wind, or it may be due to the seawater mixing at lower depths.
Another interesting feature is the vortex propagation velocity. According to numerical simulations [6], U e is between 0.3 and 0.51 m/s. Based on a near-shore observation [15], U e is 0.347 m/s. These two results indicate that U e is only 1/3 to 1/2 of U 0 . However, this study indicates that the average U e is 0.95 m/s, which is about 83% of the U 0 value. In the island vortex of the Luzon Strait, a previous study calculated that U e is about 89% of U 0 [16]. We believe that the calculations in this study are reasonable. The previous results may be greatly underestimated because the complete period of the vortex movement is not measured.

Seasonal Changes in U 0 and Chlorophyll-a Concentrations
The velocity data for OSCAR based on the satellite altimeter, temperature gradients, and wind fields, it can explore seasonal variations. In addition to the Kuroshio being possibly affected by mesoscale eddy invasions [29], the seasonal variation of velocity also significantly affects wake evolution. Figure 12a-c show the distribution of the flow field of the Kuroshio annually and in summer and winter, respectively. Figure 12d shows the velocity changes each month. Table 2 summarizes the incoming flow speeds from 2010 to 2019. The fastest incoming flow speed was in July, which was 2.4 times the slowest in November. The Kuroshio is variable in different seasons. In addition to analyzing the changes in SST when observing the development of wakes [9], it is also important to observe changes in Chl-a concentration in the wake region, which may affect the aggregation of phytoplankton and fish populations. Based on the average Chl-a concentration (0.15 mg/m 3 ) in the wake region, we analyzed the percentage of data that exceeds this value in different seasons ( Figure 13). Since the surface water is relatively cold, and the wind is relatively strong in winter, more nutritious seawater in the deeper layer was mixed into the euphotic zone near the surface of the ocean, which nourishes the growth of phytoplankton. Therefore, such a broadly high Chl-a concentration is not related to the wake development. In summer, due to the vigorous development of wakes, there is a good chance that high Chl-a concentration could be produced within 15 km of the lee of the island.  Of the 101 vortex trajectories in this study, 78 cases appeared in summer (June to August). Velocity and many missing values in winter due to clouds may be the main factors behind this phenomenon. According to theoretical experiments, when 5 < Re < 40, a fixed pair of symmetric vortices occurred, which is also the most commonly observed cold wake formations in the lee of Green Island, when 40 < Re < 200, a laminar vortex street occurred, which corresponds to the vortex shedding trains observed in this study. We used the MITgcm model to establish the SST of wake     Of the 101 vortex trajectories in this study, 78 cases appeared in summer (June to August). Velocity and many missing values in winter due to clouds may be the main factors behind this phenomenon. According to theoretical experiments, when 5 < Re < 40, a fixed pair of symmetric vortices occurred, which is also the most commonly observed cold wake formations in the lee of Green Island, when 40 < Re < 200, a laminar vortex street occurred, which corresponds to the vortex shedding trains observed in this study. We used the MITgcm model to establish the SST of wake Of the 101 vortex trajectories in this study, 78 cases appeared in summer (June to August). Velocity and many missing values in winter due to clouds may be the main factors behind this phenomenon. According to theoretical experiments, when 5 < Re < 40, a fixed pair of symmetric vortices occurred, which is also the most commonly observed cold wake formations in the lee of Green Island, when 40 < Re < 200, a laminar vortex street occurred, which corresponds to the vortex shedding trains observed in this study. We used the MITgcm model to establish the SST of wake under the same ocean conditions in summer but with different incoming speeds ( Figure 14). The three incoming speed conditions correspond to Re values of 70, 118, and 156, respectively. When Re = 70, no cold wake is generated, when Re = 118, the vortex street is generated, when Re = 156, the area of a cold wake is significantly expanded, and the SST drop is also significantly increased. The physical parameters obtained in this study can improve the establishment of numerical models of wakes. In addition to changes in the incoming speed, the ocean conditions, which include the values of the sea temperature, seawater stratified structures, Kuroshio meander [29], and local wind effects can affect the evolution of wakes. We will conduct simulation studies of different situations in a future study.
Remote Sens. 2020, 12, x FOR PEER REVIEW 16 of 18 under the same ocean conditions in summer but with different incoming speeds ( Figure 14). The three incoming speed conditions correspond to Re values of 70, 118, and 156, respectively. When Re = 70, no cold wake is generated, when Re = 118, the vortex street is generated, when Re = 156, the area of a cold wake is significantly expanded, and the SST drop is also significantly increased. The physical parameters obtained in this study can improve the establishment of numerical models of wakes. In addition to changes in the incoming speed, the ocean conditions, which include the values of the sea temperature, seawater stratified structures, Kuroshio meander [29], and local wind effects can affect the evolution of wakes. We will conduct simulation studies of different situations in a future study.

Uncertainties, Errors, and Accuracies
The data of SST and Chl-a used in this study were processed by JAXA. The SST algorithm is based on the method developed for Himawari-8 SST [30]. SGLI SSTs were validated by the comparison with buoy data with 0.8 K difference, the Himawari-8 had an average SST difference of 0.18 K with a tropical atmosphere-ocean array [20]. The Chl-a concentration algorithm is developed based on the empirical algorithms for convenient use with the other sensor products, the estimated errors are −60 to +150% with in-situ data and MODIS data (https://suzaku.eorc.jaxa.jp/GCOM_C/). In this study, GCOM-C data presented significant variations in SST and Chl-a concentrations in the wake region, and the SST drop and changes in Chl-a were similar to previous in-situ measurements [9]. In addition, this study objectively used the lowest SST in the vortex to track its movement, which was hardly affected by the accuracy of the data. All data have been well managed based on cloud masking classification and decision. In this study, we chose cloud-free data for analysis.

Conclusions
This study used 250 m spatial resolution GCOM-C data and 1-hour temporal resolution Himawari-8 imagery to analyze the surface structure and dynamic processes of the Green Island wakes. These two sets of satellite data help us more understand the characteristic of sub-mesoscale eddies in the Kuroshio region. Based on the GCOM-C data, we designed observation lines that were different from the cruise observations because every point on the imaginary lines was certain to occur at the same time. Details of the spatial structure of the wake are revealed in this study, including detailed SST and Chl-a changes, the fronts of SST and Chl-a between the wake and Kuroshio water, the distance between two consecutive vortices, the vorticity transition response to sea surface characteristics, and the different structures of SST and Chl-a in the same vortex. Based on the Himawari-8 imagery, the incoming current speed and the propagation speed of the vortex could be calculated. In total, 101 vortex cases from July 2015 to December 2019 were calculated. About 77% of the cases appeared in summer. The average vortex propagation speed was 0.95 m/s. More than half

Uncertainties, Errors, and Accuracies
The data of SST and Chl-a used in this study were processed by JAXA. The SST algorithm is based on the method developed for Himawari-8 SST [30]. SGLI SSTs were validated by the comparison with buoy data with 0.8 K difference, the Himawari-8 had an average SST difference of 0.18 K with a tropical atmosphere-ocean array [20]. The Chl-a concentration algorithm is developed based on the empirical algorithms for convenient use with the other sensor products, the estimated errors are −60 to +150% with in-situ data and MODIS data (https://suzaku.eorc.jaxa.jp/GCOM_C/). In this study, GCOM-C data presented significant variations in SST and Chl-a concentrations in the wake region, and the SST drop and changes in Chl-a were similar to previous in-situ measurements [9]. In addition, this study objectively used the lowest SST in the vortex to track its movement, which was hardly affected by the accuracy of the data. All data have been well managed based on cloud masking classification and decision. In this study, we chose cloud-free data for analysis.

Conclusions
This study used 250 m spatial resolution GCOM-C data and 1-h temporal resolution Himawari-8 imagery to analyze the surface structure and dynamic processes of the Green Island wakes. These two sets of satellite data help us more understand the characteristic of sub-mesoscale eddies in the Kuroshio region. Based on the GCOM-C data, we designed observation lines that were different from the cruise observations because every point on the imaginary lines was certain to occur at the same time. Details of the spatial structure of the wake are revealed in this study, including detailed SST and Chl-a changes, the fronts of SST and Chl-a between the wake and Kuroshio water, the distance between two consecutive vortices, the vorticity transition response to sea surface characteristics, and the different structures of SST and Chl-a in the same vortex. Based on the Himawari-8 imagery, the incoming current speed and the propagation speed of the vortex could be calculated. In total, 101 vortex cases from July 2015 to December 2019 were calculated. About 77% of the cases appeared in summer. The average vortex propagation speed was 0.95 m/s. More than half of the vortices had propagation speeds between 0.8 and 1.2 m/s. The average incoming surface current speed of Green Island was 1.15 m/s, which was calculated by the maximum cross-correlation method. In the 38 cases which have two continuous vortices, the average vortex shedding period was 14.8 h, the corresponding average Strouhal number was 0.114, and the Reynolds number was 64. In this study, the St-Re fitting curve relation was used to discuss the calculation results of the Green island wake compared to previous studies, and the results suggest that the size of Green Island is suitably selected from 5 to 6 km. Further, the 100 m 2 /s of the horizontal eddy viscosity is still suitable. FIn addition, there is a good chance that a Chl-a larger than 0.15 mg/m 3 will be produced within 15 km at the lee of the island. This study used new remote sensing data to successfully observe and analyze the dynamic processes of the sub-mesoscale vortices. In the future, Himawari-8, the SGLI of GCOM-C, and the Sentinel-3 data will help establish a wake database to be used for ocean sustainability development and assistance in the protection of fishery resources.