Statistical Characteristics of Cyclonic Warm-Core Eddies and Anticyclonic Cold-Core Eddies in the North Pacific Based on Remote Sensing Data

A (an) cyclonic (anticyclonic) eddy is usually associated with a cold (warm) core caused by the eddy-induced divergence (convergence) motion. However, there are also some cyclonic (anticyclonic) eddies with warm (cold) cores in the North Pacific, named cyclonic warm-core eddies (CWEs) and anticyclonic cold-core eddies (ACEs) in this study, respectively. Their spatio-temporal characteristics and regional dependence are analyzed using the multi-satellite merged remote sensing datasets. The CWEs are mainly concentrated in the northwestern and southeastern North Pacific. However, besides these two areas, the ACEs are also concentrated in the northeastern Pacific. The annual mean number decreases year by year for both CWEs and ACEs, and the decreasing rate of the CWEs is about two times as large as that of the ACEs. Moreover, the CWEs and ACEs also exhibit a significant seasonal variation, which are intense in summer and weak in winter. Based on the statistics of dynamic characteristics in seven subregions, the Kuroshio Extension region could be considered as the most active area for the CWEs and ACEs. Two possible mechanisms for CW-ACEs generation are discussed by analyzing two cases.


Introduction
As satellite altimeter data becomes more and more abundant, it is found that the energetic mesoscale eddies are ubiquitous in the global ocean surface [1][2][3][4].Subsequently, the implementation of the Argo project greatly extends the mesoscale eddy research studies from the sea surface to the inner ocean [5][6][7][8][9].These previous studies show that the mesoscale eddies play significant roles in several aspects of the oceanic and atmospheric environment.
In most of the previous studies, the cyclonic (anticyclonic) eddy is often associated with a cold (warm) core, so the cyclonic (anticyclonic) eddy is also known as cold (warm) eddy.This result is usually obtained from the composite analysis, which assumes that the same polarity eddies (cyclonic or anticyclonic) share similar characteristics.Although the composite analysis captures the main features in the statistics, and eliminate some random and measurement errors, it could cover up some interesting characteristics as well.For example, cyclonic (anticyclonic) eddies with warm (cold) cores are not identified using this method.
The cyclonic warm-core eddy (CWE) and anticyclonic cold-core eddy (ACE) are commonly referred to as CW-ACE in this study.The definition of the CW-ACE can be found in Section 2.2.2.The CW-ACE usually occur less frequently than the "normal" eddies.Therefore, there are only a few literature studies concerning the CW-ACE, and most are based on case studies in high latitude areas [55][56][57][58][59].
An anticyclonic cold-eddy with a diameter of about 16 km was found at the Chukchi Shelf in the western Arctic Ocean [57].This ACE formed at the continental slope in September 2004, and was centered at the depth of 160 m.The nutrient concentration in this ACE core was higher than that in other waters of similar density in the deep Canada Basin.Based on a CTD and shipboard ADCP survey data from 13 to 29 May 1994, another ACE with a cold/low-salinity core and low potential vorticity (named as anticyclonic Oyashio eddy) was observed on the eastern coast of Japan [60].This ACE, which originated from the Okhotsk Sea, played a significant role in the oceanic circulation between the Okhotsk Sea and the western subarctic gyre.Combining the vertical profile of temperature and salinity with satellite altimeter data, Itoh and Yasuda [58] investigated the water mass characteristics of anticyclonic warm-and cold-core eddies in the western boundary region of the subarctic North Pacific.The results suggested that ACEs mainly appear around the Oyashio southward intrusion areas and farther north near the Kuril Islands.
Through cross-contrasting different detection results from different data sources (sea surface temperature (SST), sea surface height anomaly (SSHA), and surface drifter trajectories data), the CW-ACEs almost existed in each latitudinal band of the Kuroshio Extension (KE) region, and ACEs were significantly more than CWEs [61].In addition, CWEs mainly appeared in the south side of the KE, while ACEs mainly concentrated in the north.Besides these studies, an intrathermocline lens with strong anticyclonic circulation and cold surface temperature anomaly (−1 • C) (i.e., ACE) was reported in the Iceland Basin [62].
About the formation mechanism of the CW-ACEs, most studies mainly concentrated on ACEs and their relationship with the Okhotsk Sea.Previous studies mentioned five situations for the generation of CW-ACE.(1) ACE forms from mode water created by deep winter mixing in the Iceland Basin region [21].(2) ACE forms south of the Bussol' Strait and intensifies with the supply of Okhotsk Sea's cold and fresh water [55,63].(3) The "normal" anticyclonic eddy generated in the Kuroshio region moves northeastward until the Bussol' Strait.With the supply of Okhotsk Sea water, it then amplifies and changes to ACE [21].(4) The CW-ACE form from the baroclinic instability of the shelf-edge current on the Chukchi Sea continental slope area [64][65][66].
Due to the limited observational data in the past, the previous literature was almost based on the case analysis, and the understanding of the CW-ACEs was restricted.In order to fill that gap, the systematic statistical analysis in this study is executed in the North Pacific (99 The major oceanic circulations and annual mean eddy kinetic energy (EKE) in this area are schematically shown in Figure 1.The Equatorial Countercurrent (ECC) is seen flowing eastward in the south of the North Pacific.To the north of the ECC is the North Equatorial Current (NEC) flowing from east to west.When the NEC reaches east of the Philippine Islands, it bifurcates into two limbs.The southward branch is Mindanao Current (not shown in Figure 1) and the northward branch forms the famous western boundary current-Kuroshio.The Kuroshio flows roughly northeast, past the south of Japan.At about 35 • N, the Kuroshio Current turns to the east and re-enters the North Pacific, forming the Kuroshio Extension [67].To the north of the KE is the southwestward flowing Oyashio Current (OC).The OC splits into two branches at about 42 • N; one turns east to join the Subarctic Current, and the other branch continues flowing southward and joins the KE.To the south of the KE is the Subtropical Countercurrent (STCC) flowing from west to east.At the eastern boundary, the eastward-flowing North Pacific Current also splits into two currents, the northward Alaska Current (AC) and the southward California Current (CC).The northernmost current in the North Pacific is the Alaska Stream, flowing southwestward along the Aleutian Islands.Besides, there are three energetic wind jets [68], also shown in Figure 1: the Panama Gulf (point A), Papagayo Gulf (point B), and Tehuantepec Gulf (point C).
In addition to the major currents, the spatial variation of the EKE is also significant in Figure 1.As pointed out in the previous literature [53,61], the KE region is the area with the highest EKE in the North Pacific.Another high EKE activity area is located along the STCC region, extending from the east of the Luzon Strait all the way to the Hawaiian Islands.On the contrary, the lowest EKE appears at the northeast region of the North Pacific.Based on the color shown in Figure 1, the value of the EKE in the KE region is on the order of O (10 3 cm 2 s −2 ) and on the order of O (10 1 cm 2 s −2 ) in the northeast region of the North Pacific.In addition to the major currents, the spatial variation of the EKE is also significant in Figure 1.As pointed out in the previous literature [53,61], the KE region is the area with the highest EKE in the North Pacific.Another high EKE activity area is located along the STCC region, extending from the east of the Luzon Strait all the way to the Hawaiian Islands.On the contrary, the lowest EKE appears at the northeast region of the North Pacific.Based on the color shown in Figure 1, the value of the EKE in the KE region is on the order of O (10 3 cm 2 s −2 ) and on the order of O (10 1 cm 2 s −2 ) in the northeast region of the North Pacific.The main goal of this study is to demonstrate the statistical characteristics of the CW-ACE in the North Pacific.Considering that different dynamic backgrounds in the subregions may lead to distinct characteristics of the CW-ACE, this paper divides the North Pacific into seven subregions (for detailed range, please refer to Appendix A) and discusses the CW-ACE dynamic characteristics in each subregion as well.The same seven subregions were used by Cheng et al. [69].This work is expected to improve the understanding of the CW-ACE phenomena in the North Pacific.The paper is organized as follows.The data and methods used in this study are introduced in Section 2. Section 3 presents the cases of the CW-ACE, along with the statistical spatial distribution, general characteristics, inter-annual and seasonal variation, and the different CW-ACEs characteristics in seven subregions of the North Pacific.Section 4 discusses the potential generation mechanism of the CW-ACE.Finally, the summary of this study is given in Section 5.

Data
Two satellite remote sensing datasets are used: a new version of SSHA data and a satellite-measured SST data.The Archiving, Validation, and Interpretation of Satellite Oceanographic (AVISO), multiple satellite-merged SSHA data (http://www.aviso.oceanobs.com/) is applied to extract and track the mesoscale eddies in the North Pacific.This dataset merges several satellite altimeter measurements to obtain a daily product, which has a spatial resolution of 1/4 • by 1/4 • .The data from 1 January 1993, to 31 December 2017, is used in this study.For more details about this new version of SSHA dataset, please refer to Pujol et al. [70].
The sea surface geostrophic velocity anomalies is derived from the geostrophic current formula: (u , v ) = g f (− ∂h ∂y , ∂h ∂x ), which assumes the balance between pressure gradient force and Coriolis force.In the above formula, u and v are the zonal and meridional components of the geostrophic velocity anomaly; h , g, and f are the sea surface height anomaly, gravitational acceleration parameter, and the Coriolis parameter, respectively.This original geostrophic velocity anomaly data is linearly interpolated onto a 1/8 • by 1/8 • grid before using the eddy detection method.
To extract the CW-ACEs, daily SST data with a horizontal resolution of 1/4 • by 1/4 • is used.The remotely-sensed SST is derived from the optimum interpolation Advanced Very High Resolution Radiometer (AVHRR-only-v2) SST analysis product (https://www.ncei.noaa.gov/thredds/catalog/OisstBase/NetCDF/AVHRR/catalog.html)provided by National Centers for Environmental Information (NCEI).The same period and interpolation as SSHA data are applied for the SST data.For more details, please refer to Reynolds et al. [71].

Eddy Detection Method
An automatic eddy detection method [72] is employed to detect and track mesoscale eddies.This method is based on the vector geometry characteristic of the mesoscale eddy velocity.Compared with the Okubo-Weiss method [73,74] and the Winding Angle method [75], it has a higher successful (lower excessive) identification rate [72], and can detect the mesoscale eddy well [15,[76][77][78][79][80][81].The detailed method is given in Appendix B.

Definition of Cyclonic Warm-Core and Anticyclonic Cold-Core Eddies
The mesoscale eddies are identified satisfying the following two criteria: one is the lifespan longer than or equal to 30 days, and the other is the mean eddy radius not less than 25 km.Besides the above two conditions, the CW-ACEs are defined satisfying two extra conditions.One condition is the temperature anomaly (T , the mean temperature difference between the eddy and background area) should be warmer (cooler) than or equal to 0.1 • C (−0.1 • C) for the cyclonic (anticyclonic) eddy.The background area is defined as an annular region surrounding the eddy, with its inner boundary being precisely the eddy boundary (solid curve in Figure 2), and the outer boundary is 1.5 times the radius of the eddy (dashed curve).The other condition is that the ratio (γ p ) of the abnormal points to all the points within the eddy should be larger than or equal to 60%, i.e., γ p ≥ 60%.The abnormal point indicates that the local temperature is warmer (cooler) than the background mean temperature for the cyclonic (anticyclonic) eddy.
From January 1993 to December 2017, 5,706,237 cyclonic eddies (snapshot) and 5,658,354 anticyclonic eddies were identified in the North Pacific, corresponding to 55,345 cyclonic and 52,578 anticyclonic eddies with lifespan.Based on the definition of the CW-ACE, there are 251,586 (4.4%) CWEs and 346,509 (6.1%) ACEs in the research period.

Cyclonic Warm-Core and Anticyclonic Cold-Core Eddy Cases
According to the definition of CW-ACE in this study, four cases (2 CWEs and 2 ACEs) are shown in Figure 2. The temperature at all the points within these CWE (ACE) cases is warmer (cooler) than the background mean temperature, i.e., γ p = 100% in all of the four cases.

Cyclonic Warm-Core and Anticyclonic Cold-Core Eddy Cases
According to the definition of CW-ACE in this study, four cases (2 CWEs and 2 ACEs) are shown in Figure 2. The temperature at all the points within these CWE (ACE) cases is warmer (cooler) than the background mean temperature, i.e.,  p = 100% in all of the four cases.Figure 2a displays day 43 of a CWE, which appeared on 1 December 1998.The eddy center is located at (167.125°W, 55.875°N), the mean radius is about 56.3 km, the amplitude (the maximum absolute value of SSHA) reaches 7.9 cm, and the maximum  T reaches 0.48 °C.This eddy survived 45 days, and turned from "normal" cyclonic eddy to CWE on day 35.In the other case, the eddy appeared on 21 August 2010 (Figure 2b).For day 267, its center is located at (112.625°W, 22.625°N), the mean radius is about 52.5 km, the amplitude is 8.4 cm, and the maximum  T reaches 0.42 °C.This eddy survived 281 days, and turned from "normal" cyclonic eddy to CWE on day 25.
Besides the CWEs, there are ACEs in the North Pacific as well.Figure 2c presents an ACE case occurring on 27 April 2001.Its center is located at (150.375°W, 15.125°N), the mean radius is about 82.5 km, the amplitude reaches 12.7 cm, and the maximum  T is  0.54 °C.This eddy survived 92 days and turned from "normal" anticyclonic eddy to ACE on day 85, and the figure shows the eddy on day 87.The second ACE case emerges on 19 June 2005 (Figure 2d).Its center is located at (155.5°E, 42.0°N), the mean radius is about 55.5 km, the amplitude reaches 14.4 cm, and the maximum  T is  0.77 °C.This eddy survived 105 days, and it turned from "normal" anticyclonic eddy to ACE on day 93. Figure 2d shows the eddy on day 99. Figure 2a displays day 43 of a CWE, which appeared on 1 December 1998.The eddy center is located at (167.125 • W, 55.875 • N), the mean radius is about 56.3 km, the amplitude (the maximum absolute value of SSHA) reaches 7.9 cm, and the maximum T reaches 0.48 • C.This eddy survived 45 days, and turned from "normal" cyclonic eddy to CWE on day 35.In the other case, the eddy appeared on 21 August 2010 (Figure 2b).For day 267, its center is located at (112.625 • W, 22.625 • N), the mean radius is about 52.5 km, the amplitude is 8.4 cm, and the maximum T reaches 0.42 • C.This eddy survived 281 days, and turned from "normal" cyclonic eddy to CWE on day 25.
Besides the CWEs, there are ACEs in the North Pacific as well.Figure 2c presents an ACE case occurring on 27 April 2001.Its center is located at (150.375 • W, 15.125 • N), the mean radius is about 82.5 km, the amplitude reaches 12.7 cm, and the maximum T is −0.54 • C.This eddy survived 92 days and turned from "normal" anticyclonic eddy to ACE on day 85, and the figure shows the eddy on day 87.The second ACE case emerges on 19 June 2005 (Figure 2d).Its center is located at (155.5 • E, 42.0 • N), the mean radius is about 55.5 km, the amplitude reaches 14.4 cm, and the maximum T is −0.77• C.This eddy survived 105 days, and it turned from "normal" anticyclonic eddy to ACE on day 93. Figure 2d shows the eddy on day 99.
From the four cases, it seems that the CW-ACE can appear in different seasons and different subregions of the North Pacific.Therefore, it is necessary to systematically study the spatio-temporal characteristics and regional dependence of the CW-ACE.In the Discussion section, their generation mechanisms are introduced.

Spatial Distribution Characteristics
In order to illustrate the spatial distribution of the CW-ACEs, Figure 3 presents the annual mean number of the ACEs (Figure 3a) and their occurrence frequency (γ b , Figure 3b) in each 5  From the four cases, it seems that the CW-ACE can appear in different seasons and different subregions of the North Pacific.Therefore, it is necessary to systematically study the spatio-temporal characteristics and regional dependence of the CW-ACE.In the Discussion section, their generation mechanisms are introduced.

Spatial Distribution Characteristics
In order to illustrate the spatial distribution of the CW-ACEs, Figure 3 presents the annual mean number of the ACEs (Figure 3a) and their occurrence frequency (  b , Figure 3b) in each 5° by 5° bin;  b is defined as the ratio of  reaches 194 (17%) per year, and the minimum number is only 5 (less than 1%) per year in the NEC region.
Similarly, the annual mean number of the ACEs and the associated occurrence frequency are shown in Figure 4.The ACEs are mainly concentrated in the northwest (north of 25°N, west of 180°, mainly in the OC and the KE region), southeast (10°N~30°N, east of 145°W, mainly in the CC region),   The previous studies suggested that the CW-ACEs could be generated from "normal" eddies [55], and thus "When do 'normal' eddies turn to CW-ACEs?" is a rather interesting question.To address this question, Figure 5 gives the proportion of generation ratio (  g ) of the CW-ACEs in the different lifespan stages and the annual mean occurrence frequency of CW-ACEs (  s ) in these stages (Figure 5b).Here,  g is defined as the ratio of   The previous studies suggested that the CW-ACEs could be generated from "normal" eddies [55], and thus "When do 'normal' eddies turn to CW-ACEs?" is a rather interesting question.To address this question, Figure 5 gives the proportion of generation ratio (γ g ) of the CW-ACEs in the different lifespan stages and the annual mean occurrence frequency of CW-ACEs (γ s ) in these stages (Figure 5b).Here, γ g is defined as the ratio of N abn,generation to N tot , where N abn,generation is the number of newly generated CW-ACEs in the target normalized lifespan stage and N tot is the total number of cyclonic From Figure 5a, the CWEs (blue curve) mainly generate in the beginning stage and the decaying stage, while fewer CWEs generate in the mature stage.This characteristic may come from the eddy instability during the beginning and decaying stages, while it is more coherent in the mature stage.Liu et al. [76] pointed out that the eddy's radius, vorticity, and shape almost do not change in the mature stage.This characteristic prevents the generation of the CW-ACEs (for more details, please refer to the Discussion section).
It is clearly seen that  s of the CWEs (Figure 5b, blue curve) is high in the beginning and decaying stages, while low in the mature stage.This indicates that some CWEs generate in the beginning stage, but turn into "normal" eddies afterwards (at least no longer meet the definition of the CWE), which reduces the proportion in the mature stage.However, some new CWEs are generated in the decaying stage, thus increasing the proportion again.
On the contrary, the proportion of the ACEs keeps growing continually from 0.1~0.2normalized lifespan stage until dissipation (Figure 5b, red curve).Therefore, once an ACE forms, it is harder for it to turn back to "normal" again.Besides, the "normal" anticyclonic eddies continue to transform into ACEs in the following stages, leading to a higher and higher proportion until the end.Such different evolutions may imply different generation mechanisms between the CWEs and ACEs.Further studies of the three-dimensional structure of CW-ACEs are therefore needed.

Cyclonic Warm-Core and Anticyclonic Cold-Core Eddy Survival Time
As in Figure 5 and the previous literature [55], the "normal" eddies can transform into CW-ACEs, and reverse back.In other words, CW-ACEs may only occupy some periods of the whole lifespan.Figure 6 gives the annual mean distribution of  T , which is defined as the ratio of  From Figure 5a, the CWEs (blue curve) mainly generate in the beginning stage and the decaying stage, while fewer CWEs generate in the mature stage.This characteristic may come from the eddy instability during the beginning and decaying stages, while it is more coherent in the mature stage.Liu et al. [76] pointed out that the eddy's radius, vorticity, and shape almost do not change in the mature stage.This characteristic prevents the generation of the CW-ACEs (for more details, please refer to the Discussion section).
It is clearly seen that γ s of the CWEs (Figure 5b, blue curve) is high in the beginning and decaying stages, while low in the mature stage.This indicates that some CWEs generate in the beginning stage, but turn into "normal" eddies afterwards (at least no longer meet the definition of the CWE), which reduces the proportion in the mature stage.However, some new CWEs are generated in the decaying stage, thus increasing the proportion again.
On the contrary, the proportion of the ACEs keeps growing continually from 0.1~0.2normalized lifespan stage until dissipation (Figure 5b, red curve).Therefore, once an ACE forms, it is harder for it to turn back to "normal" again.Besides, the "normal" anticyclonic eddies continue to transform into ACEs in the following stages, leading to a higher and higher proportion until the end.Such different evolutions may imply different generation mechanisms between the CWEs and ACEs.Further studies of the three-dimensional structure of CW-ACEs are therefore needed.

Cyclonic Warm-Core and Anticyclonic Cold-Core Eddy Survival Time
As in Figure 5 and the previous literature [55], the "normal" eddies can transform into CW-ACEs, and reverse back.In other words, CW-ACEs may only occupy some periods of the whole lifespan.Figure 6 gives the annual mean distribution of γ T , which is defined as the ratio of T abn,span to T tot , where T abn,span is the time span during which the eddy is an CWE or ACE for each eddy track, and T tot corresponds to the whole eddy lifespan, i.e., γ T = T abn,span T tot × 100%.In order to present the figure clearly, a logarithmic scale (base 10) for γ T is used along the y-axis.About (69.5 ± 0.8)% of the ACEs survives no longer than 10% of the whole lifespan, and only (1.3 ± 0.6) out of ten thousand persist longer than 90%.For the CWEs, the survival time of about (63.6 ± 1.1)% is no more than 10% of the whole lifespan, and only 2.9 ± 0.8 out of ten thousand CW-ACEs survive longer than 90%.
present the figure clearly, a logarithmic scale (base 10) for  T is used along the y-axis.About (69.5 ± 0.8)% of the ACEs survives no longer than 10% of the whole lifespan, and only (1.3 ± 0.6) out of ten thousand persist longer than 90%.For the CWEs, the survival time of about (63.6 ± 1.1)% is no more than 10% of the whole lifespan, and only 2.9 ± 0.8 out of ten thousand CW-ACEs survive longer than 90%.

Inter-Annual and Seasonal Variation
In some respects, the number of CW-ACE indirectly reflect the mesoscale eddies activity.The time series of CW-ACE number and the associated occurrence frequency (  y ) are shown in Figure 7a,b, respectively;  y is defined as the ratio of

Inter-Annual and Seasonal Variation
In some respects, the number of CW-ACE indirectly reflect the mesoscale eddies activity.The time series of CW-ACE number and the associated occurrence frequency (γ y ) are shown in Figure 7a,b, respectively; γ y is defined as the ratio of N abn,year to N tot,year , where N abn,year is the number of CWEs (ACEs) in the target year, and N tot,year is the total number of cyclonic (anticyclonic) eddies in the same year, i.e., γ y = N abn,year N tot,year × 100%.As a whole, the CW-ACE number and the associated occurrence frequency γ y show a decreasing pattern year by year.The linear decreasing rate of the CWEs (ACEs) number is −162.1 ± 26.9 (−80.3 ± 35.8) (Figure 7a, dashed line) per year.It means that the CWEs decay about twice as fast as the ACEs.The linear slope of γ y is −0.11 ± 0.03 for CWEs and −0.05 ± 0.03 for ACEs (Figure 7b, dashed line).Similarly, the decreasing rate of the CWEs occurrence frequency is about twice as fast as the ACEs.
As pointed out by Zhang et al. [26], the horizontal advection of eddy currents can lead to an asymmetric temperature anomaly dipole near the surface, whereas the vertical advection of upwelling and downwelling results in a temperature core at deeper levels in the large SST gradient area.In other words, the sea water temperature within eddies are determined by two reasons, the surrounding area sea water temperature and the eddy-induced vertical movement.Considering that the vertical distribution of the sea water temperature usually does not significantly change, this paper assume that the decreasing trend of the CW-ACEs may be caused by the SST gradient changes.Figure 7c gives the time series of the annual mean SST gradient, where a similar variation trend to the CW-ACEs decreasing pattern is observed.The linear decreasing rate of the SST gradient is (−0.076 ± 0.008) × 10 −3 • C per year.The correlation coefficients between the SST gradient curve and CWE, ACE, and CW-ACE curves are 0.69, 0.54, and 0.68, respectively, and are all statistically significant on the 90% confidence level.This means the SST gradient changes can be considered as the main reason for the CW-ACEs decreasing trend.
The number of CW-ACEs and the corresponding occurrence frequency (γ se ) show significant seasonal variation (Table 1).Here, γ se is defined as the ratio of N abn,season to N tot,season , where N abn,season is the number of the CWEs (ACEs) appearing in the target season, and N tot,season is the total number of the cyclonic (anticyclonic) eddies in the same season, i.e., γ se = N abn,season N tot,season × 100%.The CW-ACEs are active in summer (June, July, August), while inactive in winter (December, January, February).The maximum number (γ se ) of CWEs is in summer, up to 3618.6 ± 357.8 ((11.43 ± 1.11)%), and the minimum number (γ se ) is only 1896.6 ± 311.6 ((6.28 ± 0.87)%) in winter.Correspondingly, for the ACEs, the maximum and minimum number (γ se ) are 5029.4± 480.1 ((14.37 ± 1.23)%) and 2231.4 ± 275.9 ((7.39 ± 0.80)%), appearing in summer and winter, respectively.In order to show the seasonal variation clearly, the climatological monthly number, and the associated occurrence frequency (γ m ) of the CW-ACEs are displayed in Figure 8.Similarly, γ m is defined as the ratio of N abn,month to N tot,month , where N abn,month is the number of the CWEs (ACEs) appearing in the target month, and N tot,month is the total number of the cyclonic (anticyclonic) eddies in the same month, i.e., γ m = N abn,month N tot,month × 100%.The CWEs are the most active in July, with a maximum number (γ m ) up to 1355.7 ± 26.6 ((12.1 ± 0.24) %), while they are the most inactive in February, with a minimum number (γ m ) of only 603.9 ± 26.3 ((6.1 ± 0.23) %).Likewise, for the ACEs, the maximum number (γ m ) is also in July, reaching 1757.4 ± 32.4 ((14.9 ± 0.25) %), and the minimum monthly mean number (γ m ) is only 702.8 ± 24.4 ((6.9 ± 0.22) %) in January.

The Regional Dependence
The ocean currents are complex in the North Pacific (Figure 1), thus seven subregions are defined (Appendix A) to investigate the detailed feature of the CW-ACEs [69].The mean radius, temperature anomalies, EKE, and amplitude of the CW-ACEs all have large standard deviations (Table 2), indicating less homogeneity.Thus, composite analysis cannot be used to compare the characteristics in the different subregions.Therefore, the boxplots of the CW-ACEs characteristics within the 25 % to 75 % are used to compare the differences in the seven subregions (Figure 9).

The Regional Dependence
The ocean currents are complex in the North Pacific (Figure 1), thus seven subregions are defined (Appendix A) to investigate the detailed feature of the CW-ACEs [69].The mean radius, temperature anomalies, EKE, and amplitude of the CW-ACEs all have large standard deviations (Table 2), indicating less homogeneity.Thus, composite analysis cannot be used to compare the characteristics in the different subregions.Therefore, the boxplots of the CW-ACEs characteristics within the 25 % to 75 % are used to compare the differences in the seven subregions (Figure 9).−0.18 ± 0.07 (−0.10, −0.37) The mean radius of the CW-ACEs varies in the different subregions (Figure 9a), and decreases with the increasing latitude.The CW-ACEs with the smallest radius appear in the OC region, while those with the largest radius are in the NEC region.This may be associated with the baroclinic Rossby deformation radius changing with latitude.Besides, the absolute value of the maximum temperature difference (|T |, Figure 9b) could be regarded as another indicator of the CW-ACEs activity.The largest |T | appears in the KE region, and the smallest |T | is in the NEC region.Correspondingly, the KE region is the strongest subregion of the CW-ACEs activity in the North Pacific.
EKE is considered as the third indicator to measure the CW-ACEs activity (Figure 9c).The EKE is highest in the KE region (especially for the CWEs), where it is several times higher than that in the AC region.The amplitude of the CW-ACEs, considered as the fourth indicator, is shown in Figure 9d.Similar to EKE, the largest and smallest amplitude also appear in the KE and AC region, respectively.To sum up, the KE region is the strongest subregion of the CW-ACEs activity in the North Pacific, in terms of the CW-ACE-induced temperature change, EKE, and amplitude.8.The abbreviations ("OC", "AC", "CC", "KE" "STCC", "NETP", and "NEC") are introduced in Appendix A.
EKE is considered as the third indicator to measure the CW-ACEs activity (Figure 9c).The EKE is highest in the KE region (especially for the CWEs), where it is several times higher than that in the AC region.The amplitude of the CW-ACEs, considered as the fourth indicator, is shown in Figure 9d.Similar to EKE, the largest and smallest amplitude also appear in the KE and AC region, respectively.To sum up, the KE region is the strongest subregion of the CW-ACEs activity in the North Pacific, in terms of the CW-ACE-induced temperature change, EKE, and amplitude.

Discussion
As presented in this study, CW-ACEs are widely distributed in the North Pacific, according to the definition.In order to deepen the understanding of this phenomenon, the cases of CW-ACEs from Figure 2a,d are chosen to discuss the generation mechanisms.As noted in Section 3, this eddy survived for 45 days, and appeared as a CWE from day 35 to day 43 (Figure 10).At the transition stage (i.e., day 35), the eddy's  T ,  p , and radius increase sharply.One possible theory is that the eddy loses its strong coherent structure and becomes very unstable when reaching the decaying stage [76].The eddy studied in this paper did appear that its radius suddenly increases when turning into CWE, and absorbs warm sea waters from the surrounding area, which in turn changes the water property within the eddy.8.The abbreviations ("OC", "AC", "CC", "KE" "STCC", "NETP", and "NEC") are introduced in Appendix A.

Discussion
As presented in this study, CW-ACEs are widely distributed in the North Pacific, according to the definition.In order to deepen the understanding of this phenomenon, the cases of CW-ACEs from Figure 2a,d are chosen to discuss the generation mechanisms.As noted in Section 3, this eddy survived for 45 days, and appeared as a CWE from day 35 to day 43 (Figure 10).At the transition stage (i.e., day 35), the eddy's T , γ p , and radius increase sharply.One possible theory is that the eddy loses its strong coherent structure and becomes very unstable when reaching the decaying stage [76].The eddy studied in this paper did appear that its radius suddenly increases when turning into CWE, and absorbs warm sea waters from the surrounding area, which in turn changes the water property within the eddy.
In order to support this assumption, Figure 11 gives the evolution of the SST and eddy boundary during before and after the transition time (day 35, i.e., 23 November 1998).It is clear that this eddy in Figure 11a is a "normal" eddy before the transition.The temperature in the surrounding area is generally lower to the west and north of the eddy and higher to the east.Compared with Figure 11a, the temperature of the sea water within the background area is rising over the period of observation, and this rising ambient temperature offers the chance to generate the CWE.On 24 November 1998 (Figure 11c), the eddy radius suddenly increases, while it absorbs warm sea waters from the surrounding area, and then turns to a CWE. Figure 11d shows the "first day" after it has turned into a CWE (25 November 1998).The eddy radius continues increasing and it moves to the north (the colder water area).This makes the mean temperature in the area surrounding the eddy decrease, and as a result, T increases.In short, at the decaying stage, the "normal" eddy radius suddenly increases and it absorbs warm sea waters from the surrounding area, and then turns into a CWE.Subsequently, the eddy moves to the colder water area and gets even stronger.In order to support this assumption, Figure 11 gives the evolution of the SST and eddy boundary during before and after the transition time (day 35, i.e., 23 November 1998).It is clear that this eddy in Figure 11a is a "normal" eddy before the transition.The temperature in the surrounding area is generally lower to the west and north of the eddy and higher to the east.Compared with Figure 11a, the temperature of the sea water within the background area is rising over the period of observation, and this rising ambient temperature offers the chance to generate the CWE.On 24 November 1998 (Figure 11c), the eddy radius suddenly increases, while it absorbs warm sea waters from the surrounding area, and then turns to a CWE. Figure 11d shows the "first day" after it has turned into a CWE (25 November 1998).The eddy radius continues increasing and it moves to the north (the colder water area).This makes the mean temperature in the area surrounding the eddy decrease, and as a result,  T increases.In short, at the decaying stage, the "normal" eddy radius suddenly increases and it absorbs warm sea waters from the surrounding area, and then turns into a CWE.Subsequently, the eddy moves to the colder water area and gets even stronger.Besides the CWE discussed above, the ACE depicted in Figure 2d is chosen to show its generation mechanism.Temporal evolution of T  , p  , and radius for the ACW is given in Figure 12.
Similar to the CWE depicted in Figure 10, the transition to ACE occurs in the final "decaying" stage.T  and p  have changed significantly during the transition process.However, the radius have not changed much then.This fact indicates that the generation mechanism of ACE is different from that of CWE.In other words, this transition to ACE is unlikely a result of its increased instability during its "decaying" stage.Besides the CWE discussed above, the ACE depicted in Figure 2d is chosen to show its generation mechanism.Temporal evolution of T , γ p , and radius for the ACW is given in Figure 12.Similar to the CWE depicted in Figure 10, the transition to ACE occurs in the final "decaying" stage.T and γ p have changed significantly during the transition process.However, the radius have not changed much then.This fact indicates that the generation mechanism of ACE is different from that of CWE.In other words, this transition to ACE is unlikely a result of its increased instability during its "decaying" stage.

T  and p
 have changed significantly during the transition process.However, the radius have not changed much then.This fact indicates that the generation mechanism of ACE is different from that of CWE.In other words, this transition to ACE is unlikely a result of its increased instability during its "decaying" stage.The sea conditions (temperature and current) of the area surrounding the ACE are given in Figure 13, for the period in which the transition into ACE occurs (11-14 June 2005), in an attempt to study its generation mechanism.Figure 13a,b are for the "normal" eddy, whereas Figure 13c,d are for the ACE. Figure 13a clearly shows that water to the north of the "normal" eddy is cold.There is another cyclonic eddy to the west, and a water mass of higher temperature exists between these two eddies.This warm water mass is dragged to the surroundings of the normal anticyclonic eddy by the strong current associated with sheared eddy.The cold water to the north is entrained into the interior of the eddy.As shown in Figure 13c, the water surrounding the anticyclonic eddy is warmed by the eddy strong shear processes, helping the transition from "normal" eddy to ACE.On the first day after the ACE's formation, the warm water mass is almost homogeneously distributed around the ACE, while the cold water from the north is entrained into the ACE interior, strengthening the ACE itself.In summary, it appears that the formation of this ACE is unlikely to be linked to instability.Rather, the interaction between adjacent eddies seem to have played an important role.It is this eddy-eddy interaction that transports warm water into the surrounding area and entrains the cold water into the ACE interior.
The generation mechanisms for CW-ACEs are discussed by studying one specific case for each kind of eddy.The two possible mechanisms found here are instability during the decay stage and eddy horizontal entrainment.The former mechanism is characterized with a sharp change in eddy radius, while in the latter case, little change in eddy radius is found.Other mechanisms, as discussed in the Introduction section, cannot be ruled out of course.The generation mechanism of CW-ACEs is still a very interesting scientific question.It appears insufficient to study this problem using only surface data, as the surface data can only reflect the surface conditions of oceanic eddies, while the information of the interior is absent.Further studies using 3-D observational data or numerical models are on the way.
the ACE, while the cold water from the north is entrained into the ACE interior, strengthening the ACE itself.In summary, it appears that the formation of this ACE is unlikely to be linked to instability.Rather, the interaction between adjacent eddies seem to have played an important role.It is this eddyeddy interaction that transports warm water into the surrounding area and entrains the cold water into the ACE interior.The generation mechanisms for CW-ACEs are discussed by studying one specific case for each kind of eddy.The two possible mechanisms found here are instability during the decay stage and eddy horizontal entrainment.The former mechanism is characterized with a sharp change in eddy radius, while in the latter case, little change in eddy radius is found.Other mechanisms, as discussed in the Introduction section, cannot be ruled out of course.The generation mechanism of CW-ACEs is still a very interesting scientific question.It appears insufficient to study this problem using only surface data, as the surface data can only reflect the surface conditions of oceanic eddies, while the information of the interior is absent.Further studies using 3-D observational data or numerical models are on the way.

Conclusions
The phenomenon of CW-ACEs is observed in the North Pacific.Following the definition of CW-ACEs presented in this study, this paper reveals the spatial distribution of the CW-ACEs in each 5° by 5° bin.The statistics indicate that the CWEs are mainly concentrated in the northwest and southeast of the North Pacific.Besides these two regions, the ACEs are also concentrated in the

Conclusions
The phenomenon of CW-ACEs is observed in the North Pacific.Following the definition of CW-ACEs presented in this study, this paper reveals the spatial distribution of the CW-ACEs in each 5 • by 5 • bin.The statistics indicate that the CWEs are mainly concentrated in the northwest and southeast of the North Pacific.Besides these two regions, the ACEs are also concentrated in the northeastern area.The highest γ b reaches 17% for the CWEs, appearing in the CC region, and 26% for the ACEs, appearing in the KE region.
The occurrence frequency (γ s ) of the ACEs is higher in the beginning and decaying stages, while lower in the mature stage.However, the occurrence frequency of the ACEs grow almost linearly in 0.1-0.9stages.It is also found that about (69.5 ± 0.8)% of the CWEs and (63.6 ± 1.1)% of the ACEs survive shorter than 10% of the whole eddy lifespan.Only 1.3 ± 0.6 out of ten thousand eddies for the CWEs and 2.9 ± 0.8 out of ten thousand eddies for the CWEs survive longer than 90% of the whole eddy lifespan.Therefore, most of the CW-ACEs are short-term phenomena and cannot survive for a long time like the "normal" eddies.
The annual mean number and the associated occurrence frequency (γ y ) decrease year by year for CW-ACEs, and the decreasing rate of the CWEs is almost twice as fast as that of the ACEs.The number of CW-ACEs (γ y ) reduces by about 162.1 ± 26.9 ((0.11 ± 0.03)%) each year for the CWEs, and about 80.3 ± 35.8 ((0.05 ± 0.03)%) for the ACEs.This decreasing trend is due to the SST gradient decreasing pattern.The CW-ACEs have a significant seasonal variation, being intense in summer and weak in winter.The maximum number (γ m ) of the CWEs appears in July, and the minimum number (γ m ) appears in February.The maximum number (γ m ) of the ACEs also appears in July, while the minimum number (γ m ) appears in January.
In view of the complex ocean currents and different dynamic environment, the North Pacific is divided into seven subregions, and the dynamic properties of the CW-ACEs (the mean radius, T , EKE, and amplitude) are discussed in each subregion.The KE region is the most active area of the CW-ACEs from the perspective of the CW-ACE-induced T , EKE, and amplitude.
In addition, generation mechanisms for CW-ACEs are discussed.The two possible mechanisms found here are (1) instability during the decay stage and (2) eddy-eddy interaction and horizontal entrainment.This study, though not exhaustive, paves the road to study CW-ACE generation mechanisms using three-dimensional data.
The statistical analysis of the CW-ACEs contributes to a comprehensive understanding of the mesoscale eddy activity, and it is beneficial to the selection of the buoy sites in the future.The analysis of the dynamic characteristics in the different subregions is conductive to comprehend the CW-ACEs activity in different environments.Although only sea surface remote sensing data is applied in this paper, it could be a foundation for further study about the three-dimensional structure and generation mechanisms of the CW-ACEs.
Remote Sens. 2018, 10, x FOR PEER REVIEW 3 of 23 the famous western boundary current--Kuroshio.The Kuroshio flows roughly northeast, past the south of Japan.At about 35°N, the Kuroshio Current turns to the east and re-enters the North Pacific, forming the Kuroshio Extension [67].To the north of the KE is the southwestward flowing Oyashio Current (OC).The OC splits into two branches at about 42°N; one turns east to join the Subarctic Current, and the other branch continues flowing southward and joins the KE.To the south of the KE is the Subtropical Countercurrent (STCC) flowing from west to east.At the eastern boundary, the eastward-flowing North Pacific Current also splits into two currents, the northward Alaska Current (AC) and the southward California Current (CC).The northernmost current in the North Pacific is the Alaska Stream, flowing southwestward along the Aleutian Islands.Besides, there are three energetic wind jets [68], also shown in Figure 1: the Panama Gulf (point A), Papagayo Gulf (point B), and Tehuantepec Gulf (point C).

Figure 2 .
Figure 2. Cases of the cyclonic warm-core eddies (a,b) and anticyclonic cold-core eddies (c,d) in the North Pacific.The shading represents the sea surface temperature (units: °C).The solid contour indicates the boundary of the cyclonic warm-core and anticyclonic cold-core eddies (CW-ACEs), which is also the inner boundary of the background region.The dashed contour refers to the outer boundary of the background region.The black asterisk is the center of the CW-ACE, and the vector denotes the sea surface geostrophic velocity anomaly (unit: m s −1 ).

Figure 2 .
Figure 2. Cases of the cyclonic warm-core eddies (a,b) and anticyclonic cold-core eddies (c,d) in the North Pacific.The shading represents the sea surface temperature (units: • C).The solid contour indicates the boundary of the cyclonic warm-core and anticyclonic cold-core eddies (CW-ACEs), which is also the inner boundary of the background region.The dashed contour refers to the outer boundary of the background region.The black asterisk is the center of the CW-ACE, and the vector denotes the sea surface geostrophic velocity anomaly (unit: m s −1 ). •

by 5 •
bin; γ b is defined as the ratio of N abn,bin to N tot,bin , where N abn,bin is the number of the CW-ACEs and N tot,bin is the total number of the same polarity eddies within the target bin, i.e., γ b = N abn,bin N tot,bin × 100%.Remote Sens. 2018, 10, x FOR PEER REVIEW 6 of 23 of the CW-ACEs and , tot bin Nis the total number of the same polarity eddies within the target bin, i.e.,

Figure 3 .
Figure 3. Distribution of the annual mean number (a) and the associated occurrence frequency (  b ) (b) of the CWEs in each 5° by 5° bin in the North Pacific.The red solid lines indicate the boundary of the active areas for the CWEs.

Figure 3 .
Figure 3. Distribution of the annual mean number (a) and the associated occurrence frequency (γ b ) (b) of the CWEs in each 5 • by 5 • bin in the North Pacific.The red solid lines indicate the boundary of the active areas for the CWEs.It is clearly observed that the CWEs are mostly concentrated in the northwest (north of 30 • N and west of 170 • W, mainly in the OC and the KE region) and southeast (north of 10 • N, south of 45 • N, and east of 135 • W, mainly in the CC region) of the North Pacific.Although the CWEs exist in other area as well, the number and γ b are significantly smaller.In the CC region, the maximum number (γ b ) reaches 194 (17%) per year, and the minimum number is only 5 (less than 1%) per year in the NEC region.Similarly, the annual mean number of the ACEs and the associated occurrence frequency are shown in Figure 4.The ACEs are mainly concentrated in the northwest (north of 25 • N, west of 180 • ,

Figure 4 .
Figure 4. Same as Figure 3, but for the ACEs.Distribution of the annual mean number (a) and the associated occurrence frequency (  b ) (b).The black solid rectangular region is the low-value band for the ACEs.

3. 3 .
General Characteristic of Cyclonic Warm-Core and Anticyclonic Cold-Core Eddy 3.3.1.Cyclonic Warm-Core and Anticyclonic Cold-Core Eddy Generation Time of newly generated CW-ACEs in the target normalized lifespan stage and tot N is the total

Figure 4 .
Figure 4. Same as Figure 3, but for the ACEs.Distribution of the annual mean number (a) and the associated occurrence frequency (γ b ) (b) of the ACEs in each 5 • by 5 • bin in the North Pacific.The black solid rectangular region is the low-value band for the ACEs.

(Figure 5 .
Figure 5. Distribution of the newly generated CW-ACEs ratio (  g ) in the different stages (a) and the annual mean occurrence frequency (  s ) of the CWEs (blue line) and ACEs (red line) within the normalized eddy lifespan (b).The error bar is defined as span is the time span during which the eddy is an CWE or ACE for each eddy

Figure 5 .
Figure 5. Distribution of the newly generated CW-ACEs ratio (γ g ) in the different stages (a) and the annual mean occurrence frequency (γ s ) of the CWEs (blue line) and ACEs (red line) within the normalized eddy lifespan (b).The error bar is defined as ±σ(t)/ √ N, where σ(t) is the standard deviation, and N is the number of years (N = 25).

Figure 6 .
Figure 6.Distribution of the annual mean  T of the CWEs (blue line) and ACEs (red line) within the eddy lifespan.Both  T and the associated standard error are in the logarithm (base 10) form.
(ACEs) in the target year, and, tot year Nis the total number of cyclonic (anticyclonic) eddies in the same year, i.e.,

Figure 6 .
Figure 6.Distribution of the annual mean γ T of the CWEs (blue line) and ACEs (red line) within the eddy lifespan.Both γ T and the associated standard error are in the logarithm (base 10) form.

Figure 7 .
Figure 7. Time series of the annual mean number (a) and the associated occurrence frequency (  y ) (b) of the CWEs (blue solid line) and ACEs (red solid line) SST gradient (c) in the North Pacific.The dashed line denotes the linear regression, the shadow indicates the confidence interval exceeding the 95 % significance.As a whole, the CW-ACE number and the associated occurrence frequency  y show a decreasing pattern year by year.The linear decreasing rate of the CWEs (ACEs) number is −162.1 ± 26.9 (−80.3 ± 35.8) (Figure 7a, dashed line) per year.It means that the CWEs decay about twice as fast as the ACEs.The linear slope of  y is −0.11 ± 0.03 for CWEs and −0.05 ± 0.03 for ACEs (Figure 7b, dashed line).Similarly, the decreasing rate of the CWEs occurrence frequency is about twice as fast

Figure 7 .
Figure 7. Time series of the annual mean number (a) and the associated occurrence frequency (γ y ) (b) of the CWEs (blue solid line) and ACEs (red solid line) SST gradient (c) in the North Pacific.The dashed line denotes the linear regression, the shadow indicates the confidence interval exceeding the 95 % significance.

23 Figure 8 .
Figure 8. Boxplots of the climatological monthly number (a,b) and the associated occurrence frequency (  m ) (c,d) of the CW-ACEs, in which (a) and (c) are for the CWEs, and (b) and (d) are for the ACEs.The line within the box means the median value and the box represents the interquartile range.The whiskers with 1.5 times of the interquartile range are lines extending from each end of the box to show the extent of the remaining data.If all data lie within the area of 1.5 times the interquartile range, the whiskers end with the minimum or maximum of the data.The symbol "+" indicates the abnormal value is out of the effective range.

Figure 8 .
Figure 8. Boxplots of the climatological monthly number (a,b) and the associated occurrence frequency (γ m ) (c,d) of the CW-ACEs, in which (a) and (c) are for the CWEs, and (b) and (d) are for the ACEs.The line within the box means the median value and the box represents the interquartile range.The whiskers with 1.5 times of the interquartile range are lines extending from each end of the box to show the extent of the remaining data.If all data lie within the area of 1.5 times the interquartile range, the whiskers end with the minimum or maximum of the data.The symbol "+" indicates the abnormal value is out of the effective range.

Figure 10 .
Figure 10.Times series of the  T (a),  p (b), and radius (c) for the CWE case.The dashed black lines are the threshold value of  T ,  p , and the mean radius.

Figure 10 . 23 Figure 11 .
Figure 10.Times series of the T (a, units: • C), γ p (b), and radius (c, units: km) for the CWE case.The dashed black lines are the threshold value of T', γ p , and the mean radius.Remote Sens. 2018, 10, x FOR PEER REVIEW 16 of 23

Figure 11 .
Figure 11.Evolution of the SST and the eddy boundary for the case of CWE on 17(a), 22 (b), 23 (c), and 25 (d) November 1998.The shading represents the SST (units: • C).The solid and dashed contours, black asterisk, and the vector are same meaning as in Figure 2.

Figure 12 .
Figure 12.Same as Figure 10, but for ACE.Figure 12. Same as Figure 10, but for ACE.(a-c) are the times series of T (units: • C), γ p and radius (units: km), respectively.

Figure 12 .
Figure 12.Same as Figure 10, but for ACE.Figure 12. Same as Figure 10, but for ACE.(a-c) are the times series of T (units: • C), γ p and radius (units: km), respectively.

Figure 13 .
Figure 13.Same as Figure 11, but for the ACE, which is from Figure 2d.

Figure 13 .
Figure 13.Same as Figure 11, but for the ACE, which is from Figure 2d.(a-d) are the evolutions of the SST and the eddy boundary on 11-14 June 2005, respectively.

Table 1 .
The number and the associated occurrence frequency (γ se ) of the CW-ACEs * in the four seasons (mean ± one standard deviation).
* Minimum and maximum values are in the brackets.

Table 2 .
Characteristics (mean radius,  T , amplitude and EKE) of the CW-ACEs * in the seven subregions of the North Pacific from 1993 to 2017 (mean ± one standard deviation).

Table 2 .
Characteristics (mean radius, T , amplitude and EKE) of the CW-ACEs * in the seven subregions of the North Pacific from 1993 to 2017 (mean ± one standard deviation).