Statistical Characteristics of Mesoscale Eddies in the North Pacific Derived from Satellite Altimetry

The sea level anomaly data derived from satellite altimetry are analyzed to investigate statistical characteristics of mesoscale eddies in the North Pacific. Eddies are detected by a free-threshold eddy identification algorithm. The results show that the distributions of size, amplitude, propagation speed, and eddy kinetic energy of eddy follow the Rayleigh distribution. The most active regions of eddies are the Kuroshio Extension region, the Subtropical Counter Current zone, and the Northeastern Tropical Pacific region. By contrast, eddies are seldom observed around the center of the eastern part of the North Pacific Subarctic Gyre. The propagation speed and kinetic energy of cyclonic and anticyclonic eddies are almost the same, but anticyclonic eddies possess greater lifespans, sizes, and amplitudes than those of cyclonic eddies. Most eddies in the North Pacific propagate westward except in the Oyashio region. Around the northeastern tropical Pacific and the California currents, cyclonic and anticyclonic eddies propagate westward with slightly equatorward (197° average azimuth relative to east) and poleward (165°) deflection, respectively. This implies that the background current may play an important role in formation of the eddy pathway patterns.


Introduction
Mesoscale eddies in the ocean profoundly influence the large-scale flux of momentum, heat, salt, and ocean-atmosphere systems [1,2].The kinetic energy of mesoscale eddies is one order greater than the mean energy of the global ocean motion [3,4].The mesoscale eddies of long lifespans enhance mixing in the upper ocean [5,6].Petersen et al. [7] applied a numerical model to investigate the three-dimensional structure of eddies in the global oceans and indicated that approximately one-third of eddies extended to at least 1000 m in depth.Meanwhile, 97% of eddies with a minimum lifespan of 4 weeks extend to the surface, so that it is feasible to detect eddies from the sea surface.For twenty years, satellite altimeter data have provided a powerful tool for detecting the mesoscale eddy fields in the global oceans.Analyzing sea surface height data from satellite altimetry, Chelton et al. [8] revealed many isolated eddy-like cyclonic and anticyclonic features and discovered that eddies propagated westward with the phase speed of long baroclinic Rossby waves.Furthermore, Chelton et al. [9] used merged altimeter data to improve the global mapping of nonlinear eddy observations and showed that the tracked eddies were generated nearly everywhere in the global oceans.Overall, the number of cyclonic eddies is greater than the number of anticyclonic eddies, but anticyclonic eddies possess greater lifespans and propagating distances.
Extensive observations by previous investigators [10,11] revealed current patterns in the North Pacific as schematically shown in Figure 1.One can see that the most remarkable feature is the western boundary current, the Kuroshio, which originates from east of Philippine Islands as a north branch of the North Equatorial Current.It enters the East China Sea through the passage of east Taiwan, and re-enters the North Pacific as the Kuroshio Extension after leaving from the Japanese coast.The Subtropical Countercurrent (STCC) is a shallow current that flows eastward between 22° and 25°N.The Alaskan Stream flows southwestward along the Aleutian Islands, whereas the Oyashio Current splits into two paths at approximately 42°N.One path veers offshore-ward and contributes to the Subarctic Current [11].The eastward-flowing North Pacific Current splits into two paths and contributes to the northward Alaskan Current and the southward California Current.In the regions with frequent eddy activities, previous investigators have analyzed the general characteristics of eddies.Liu et al. [12] constructed a dataset of the regional eddy dynamics in the STCC zone and area on the lee side of Hawaiian Islands.Yoshida et al. [13] investigated the eddy variability in the Hawaiian Lee Countercurrent (HLCC) and correlation to the time series of the Pacific Decadal Oscillation (PDO) index.Hwang et al. [14] used geostrophic current fields derived from satellite altimetry data to track mesoscale eddies in the STCC zone.Their results show that eddies propagate westward at about 8 km/day within 22°-24°N and reach the Kuroshio east of Taiwan.Itoh and Yasuda [15] analyzed characteristics of mesoscale eddies in the Kuroshio Extension/Oyashio region and revealed that eddies generally propagated at 4 km/day along the trenches south of 46°N and Subarctic Front (SAF).Giese et al. [16] indicated the existence of two bands of eddy activity off the coast of Central America in the Northeastern Tropical Pacific.Ladd et al. [17] studied eddies in the Northern Gulf of Alaska with satellite altimetry data and showed that eddies could persist for several years and moved southwestward along the Alaskan Peninsula to the Aleutian Islands.Ocean eddies can modify turbulences in the atmospheric boundary layer and locally affect near-surface wind, cloud properties, and rainfall [18].Thus, it is important to identify the distribution and characteristics of mesoscale eddies.In the past decades, many applications of satellite altimeter measurements have revolutionized studies of physical oceanography, such as estimating sea surface height and kinetic energy [19], documenting the drifting paths of mesoscale eddies [9].Previous investigators [12][13][14][15][16][17] have studied the regional characteristics of eddies in the respective active regions, using different eddy detection methods.However, little research has been conducted on the detailed properties of eddies in the North Pacific using the same detection method.This study aims to fill this gap.

Satellite Data
In this study, 20 years of sea level anomaly (SLA) are used to investigate ocean eddies in the North Pacific from October 1992 to December 2012.The SLA data is a merged product of satellite altimetry and in situ measurements provided by the Archiving Validation and Interpretation of Satellite Data in Oceanography (AVISO), the Centre National d'Études Spatiales (CNES) of France.Absolute dynamic topography (ADT) also is used to reveal background currents in the North Pacific.ADT is defined as the sea surface height above the geoid and is obtained by summing sea level anomaly (SLA) and mean dynamic topography (MDT), which is derived from multivariate objective analysis.This analysis is based on a 4.5-year Gravity Recovery and Climate Experiment (GRACE), and 15-year altimetry and in situ data (Argo floats and Conductivity-Temperature-Depth measurements).The SLA data, possessing a spatial resolution of 1°/4 by 1°/4 in a 7-day interval, are merged by multiple satellite altimeters onboard ERS-1/2, Envisat, TOPEX/Poseidon, and Jason-1/2 satellites.Ray [20] indicated that because tidal corrections are typically inadequate in shelf areas, the altimetry data in shallow-water applications is contaminated by full sea-level aliasing problems, even though tidal and inverted barometer height corrections have been incorporated.Hence, this study excluded the altimetry data in shelf areas where the water depth was less than 200 m.

Detection Algorithm
Several automated algorithms have been developed to identify eddies by utilizing satellite data or ocean model data.In this study, a free-threshold eddy identification algorithm with the SLA data is employed.This algorithm is based on the vector geometry method [21] and Okubo-Weiss method [22].We assume that the ocean is in a steady state; therefore, it is examined by considering only the balance of the pressure gradient and Coriolis force.The two-dimensional flow fields can be derived from SLA data through the geostrophic relations as: and where u' and v' are zonal and meridional velocity components, g is gravitational acceleration, f is the Coriolis parameter, R T is the mean earth radius (6371 km), ϕ and λ are latitude and longitude, respectively, and η is the SLA.The Okubo-Weiss parameter (W) is defined as [23,24]. ( where s n and s s are the normal and shear strain, respectively, and ζ is the relative vorticity defined as and Previous investigators used the W-based method [25] to define an eddy as a connected region, in which the W parameter falls below a specified negative threshold value.For example, Pasquero et al. [26] divided the flow field into a vorticity-dominated region (W < −W T ), a strain-dominated region (W > W T ), and a background field (|W| ≤ W T ), where W T = 0.2σ w , in which σ w is the standard deviation of W in the study region.Isern-Fontanet et al. [5] and Morrow et al. [27] demonstrated that a threshold value is specified according to regional characteristics.Specifying the threshold value is crucial but also subjective.Nevertheless, the Okubo-Weiss parameter is still significant regarding its physical meaning, which implies a struggle between the deformation field (|s n | + |s s |) and vorticity field (|ζ|).To objectively detect eddies based on the Okubo-Weiss method, we consider that the eddy has an isolated rotary structure exhibiting variations in the sea surface height.Six constraints are applied to the SLA to detect an eddy: (1) a vorticity-dominated region at the eddy center (W < 0) must exist; (2) the zonal gradient of SLA must reverse in sign across the eddy center; (3) the meridional gradient of SLA must reverse in sign across the eddy center and the sense of reversed sign has to be the same as for the zonal component; (4) the SLA magnitude has a local extreme value (minimum or maximum); (5) closed contours of SLA around the eddy center must exist; (6) the eddy amplitude must be larger than the sea level accuracy value which is provided by AVISO.In this study, the amplitude is defined as the absolute value of the SLA difference between the eddy center and the SLA along the eddy edge.Continuous contouring schemes of SLA every 1 cm (Figure 2b) are used to find the eddy edge, and the outermost closed contour line within which there is the same sign of relative vorticity around the eddy center is defined as the eddy edge.Figure 2 shows an example of the results of the eddy detection algorithm using SLA data.One can see that ocean eddies with various shapes and signs of vorticity can easily be detected.Furthermore, based on eddy centers and edges, eddy properties can be estimated.For example, eddy radius (R) is defined as the radius of an equivalent circular with the same area (A) as the closed contour For the characterization of eddies, the eddy kinetic energy (EKE) per unit mass is defined as

Detectability of the Detection Algorithm
Because eddy is detected by SLA data, several eddies might not extend to the surface and cannot be detected by the eddy detection method referred above.Petersen et al. [7] indicated that 97% of eddies with a minimum lifespan of 4 weeks extended to the surface, thus most mesoscale eddies could be detected from the sea surface.Detectability of the eddy detection algorithm is strongly affected by the spatial and temporal resolutions of satellite data.The spatial and temporal resolutions of SLA are 25 km and 7 days, respectively, thus an eddy diameter smaller than 50 km or a lifespan shorter than one week could not be detected in this study.The sea level accuracy used in this study is provided by AVISO.It represents errors induced by the constellation sampling capability and consistency with the spatial/temporal scales.More detailed descriptions can be found in Ducet et al. [28].The mean value and standard deviation in the North Pacific are 0.49 cm and 0.66 cm, respectively.Moreover, the shape of eddy would be affected by strong background flow or other eddies.Following previous studies [8,10,12,23], this study also uses SLA to avoid the effect of mean flow.Several eddies might be absorbed into the strong background flow, such as the Kuroshio Extension, and then be shed from it after some days [29]; this is considered as two different eddies in this study because the relevant processes cannot be fully resolved solely from SLA maps [15].

Tracking Algorithm
In addition, we develop an automated eddy tracking scheme based on the method by Chelton et al. [9].After the eddy center is detected, an eddy can be tracked at successive points in time; each eddy is tracked from a given time step t to the next step t + 1 by the closest eddy center with the same sign of relative vorticity (cyclonic/anticyclonic).The range of the searching window affects the accuracy of the eddy tracking.To avoid the risk of error tracking, such as jumping from one track to another or splitting a continuous track into multiple tracks, the search window at the latter time t should be restricted depending on the spatial and temporal resolutions of satellite data.In this study, the search radius around the eddy center is restricted within 150 km referred to the study of Chelton et al. [9].If the closest eddy is found at the latter time t + 1, it is associated with the track of the reference eddy at the time t.The entire track of each eddy is found at each subsequent time step using the above process.When no center is detected within the search window, the eddy is considered dissipated and the track is terminated.The start (end) position of an eddy track is defined as eddy generation (termination).To reduce the risk of spurious eddies arisen by noise in the SLA field, eddies with a lifespan shorter than 4 weeks are not considered in this study.

Statistical Characteristics
This study examines eddies with a lifespan longer than 4 weeks in the North Pacific.The total number of eddies is 19,836, including 10,398 cyclonic eddies and 9438 anticyclonic eddies from October 1992 to December 2012.Figure 3a shows the histogram of eddy radius distribution.One can see that it is close to the Gaussian distribution, but the peak value skews a little toward the smaller size.This result agrees with that in the STCC by Liu et al. [12].The mean radius value is 86 km for anticyclonic eddies and 85 km for cyclonic eddies.Figure 3b shows the distribution of lifespans.The mean value of lifespan is about 18 weeks, and some eddies are able to survive for a couple of years.Figure 3c shows the distribution of travel distances, which looks like the distribution of lifespans, because eddies with longer lifespans certainly propagate for longer distances.On the other hand, some eddies can have high EKE (>1000 cm 2 /s 2 ), though the mean value of EKE is about 176 cm 2 /s 2 .The kinetic energy of mesoscale eddies in the North Pacific is one order greater than the mean kinetic ocean energy [3,4].Figure 3   The eddy characteristics in the North Pacific are displayed in Figure 4. Based on the maps of EKE (Figure 4d) and the amplitude (Figure 4e), the figures show that there are several active regions of eddy, including the Oyashio, the Kuroshio Extension, the STCC, the lee side of Hawaiian Islands, the northeastern tropical Pacific, Californian currents, and Alaskan currents.Using SLA data, Liu et al. [12] showed that the regions of Kuroshio Extension and the STCC were more active with high EKE in the western North Pacific.On the other side, eddies are hardly observed around the center of the eastern part of the North Pacific Subarctic Gyre (discussed in Section 4.2).The radius, travel distance, and speed mainly vary with latitude.At lower latitudes, eddies possess greater radii, speeds, and travel distances.Due to the 25-km spatial resolution of SLA, we assume that travel distances shorter than 75 km are stationary eddies.Figure 4c reveals that most of the eddies in the North Pacific are moving eddies; meanwhile a sporadic distribution of stationary eddies appears at higher latitudes.Chelton et al. [9] indicated that eddy speeds were mainly affected by the phase speed of long baroclinic Rossby waves decreasing with latitude.The speed was smaller than 0.8 km/day at 40°N, thus stationary eddies potentially appear at higher latitudes.Furthermore, seven more active regions in the North Pacific are drawn out and discussed below (Figure 5).The Kuroshio Extension/Oyashio region encompasses the Kuroshio Extension, the Oyashio, and the Subarctic Current.In this region, eddies can persist from several months to longer than two years.The existence of eddies has a strong influence on the local climate, hydrography, and fishery [30,31].Qiu et al. [32] indicated that eddies in this region could also affect the North Pacific subtropical mode waters.The behavior of these mesoscale eddies are related to the bottom topography and variability of currents, such as meanders of the Kuroshio Extension, and southward intrusions and the eastward extension of Oyashio [33].The number of eddies in this region is around 200 per 1° latitude × 1° longitude in the past 20 years.It is more concentrated than in other places in the North Pacific.The average radius of eddies in the Kuroshio Extension and in the Oyashio is 80 km and 70 km, respectively.These results are close to those in the study of Itoh and Yasuda [15].
The lee side of Hawaiian Islands and the STCC region compose the second zonal band with strong eddy activity in the western Pacific [34].The higher density of number is similar to the Kuroshio Extension region.Qiu [35] suggested that eddies within this region could impact the subtropical gyre.The average size and propagation speed are about 90 km and 8 km/day, which is in agreement with the study of Hwang et al. [14].In this region, eddies can also persist for longer than one year.
The equatorial current zone is an area with inactivity of eddy and a lack of investigation.A total of 691 cyclonic eddies and 412 anticyclonic eddies have been detected in the last two decades.The number of cyclonic eddy is larger than that of anticyclonic eddy.The phenomenon is different from other regions iwhich the numbers of cyclonic eddies and anticyclonic eddies are similar.Velocity shear between the westward North Equatorial Current and the eastward Equatorial Countercurrent generating positive vorticity could be the reason; hence the number of cyclonic eddies is a little larger than that of anticyclonic eddies.However, further studies are needed in this area.Bauer et al. [36] investigated the eddy-mean flow decomposition in the tropical Pacific and indicated that the eddy energy was about 380 cm 2 /s 2 and the mean radius was 100 km.The lifespan of eddies is a few months, and it is the shortest in the North Pacific.
In the northern Gulf of Alaska region, mesoscale eddies are an important mechanism for cross-shelf exchange [37].As described in the study of Crawford et al. [38], eddies in this region have radii of around 80 km and can live up to 3 years.In this study we also observe that the average radius of eddies is 74 km and the lifespan could be longer than 3 years.The EKE is about 50 cm 2 /s 2 , which is smaller than in other regions.The generation mechanisms of these eddies could be reversals of downwelling winds or the interaction of the Alaska currents with the topography [39].
In the California Current region, the average eddy radius and EKE are about 85 km and 80 cm 2 /s 2 , respectively.The propagation speed is about 4 km/day.This region is the third active region of eddies in the North Pacific with an eddy number of 170/degree 2 .The distance that eddies can propagate ranges from an average of 300 km to more than 1600 km.It is believed that instabilities of the California Current and wind forcing create those energetic mesoscale eddies [40].Kurian et al. [41] revealed that most eddies had radii within a range of 25-100 km and propagation speeds from 0 to 10 km/day.
The Northeastern Tropical Pacific region is defined as the ocean from 5° to 15°N and from 120°W to the coast of Central America.Palacios and Bograd [42] used a 12-year time series of satellite altimetry data to characterize eddies.They showed that the average eddy lifespan was around 84-143 days and the radius could range from 50 to over 200 km, but cyclonic eddies dissipated more quickly [43].

Uncertainties and Errors
As described in Section 3.2, the limitations of the eddy detection method used in this study are mainly influenced by the spatial/temporal resolution of SLA.In this study, eddies can only be detected if their diameters and lifespans are greater than 50 km and 7 day, respectively.Moreover, eddies with lifespans longer than 4 weeks are only considered in this study to avoid the risk of spurious eddies that arise by noise in the SLA field.In this section, uncertainties and errors of statistical analysis are discussed.Table 1 shows the mean, standard deviation, and maximum (minimum) value of eddy characteristics in the seven active regions (Figure 5).As can be seen, lifespan, EKE, distance, and amplitude of eddies have large standard deviations.A large standard deviation implies that the values of these characters are more uncertain.Although the mean lifespan at higher latitudes is around 20 weeks, sometimes eddies can live for more than 3 years.As shown in Figure 3, the distributions of eddy kinematic properties in the North Pacific are positively skewed which is indicated by the fact that the tail on the high value side is longer than on the low value side.Because the mean values of lifespan, EKE, distance, and amplitude are not highly significant, they are not indicated for comparing the differences of eddy characters in the seven regions.Therefore, boxplots of eddy kinematics with the 25th and 75th percentile are used to compare the differences (Figure 6).The radius, lifespan, travel distance, and speed are mainly affected by latitude.Although anticyclonic eddies possess greater radii, amplitudes, and lifespans than cyclonic eddies, EKE of cyclonic eddies is not smaller than that of anticyclonic eddies.Compared to other regions, the lifespan of eddies is the greatest, but EKE is the weakest in the Alaska Current region.Moreover, the travel distance of eddies in this region can be greater than 600 km.Ladd et al. [44] a Minimum and maximum values are in brackets.

Shape Error of Eddy
The shape of eddy might be affected by a strong background flow or other eddies, hence the eddy shape would not be as an ideal circle.A test of shape error of eddies detected by the algorithm is discussed in this section.As shown in Figure 7a, closed contours of the eddy edge would have gap and bulge areas within and outside the fitted ideal circle, respectively.Kurian et al. [41] defined the shape error of eddy as the ratio between the area of deviations to the area of the ideal circle.A histogram of the shape error reveals that the eddy detected by the algorithm in the North Pacific is 24% of the mean value with 11% of one standard deviation (Figure 7b).The shape error of less than 50% is account for by 98% of the total eddies in the North Pacific.This implies that the eddy detection algorithm used in this study has successfully reduced the effect of mean flow on eddy shape.

Eddy Generation and Termination
The formation and the termination of eddies are defined by the first and the last record of time series of each eddy lifespan, respectively.Figure 8 shows the spatial distribution of eddy formation and termination.The results indicate that eddies can generate and terminate in most regions, except for the center of the eastern part of the North Pacific Subarctic Gyre.The west coast of North America, the coast of Central America, the Kuroshio Extension, the Oyashio, and the lee side of Hawaiian Islands are the major generation locations.In the Alaska Current region, eddies are often formed near the coast of Yakutat and Sitka [17].Kurian et al. [41] observed eddies in the California Current system and indicated that many of eddies formed near the coast with high eddy kinetic energy (0.018 m 2 /s 2 ).In the northeastern tropical Pacific, the ocean is strongly modulated by three powerful wind jets over the Tehuantepec, Papagayo, and Panama Gulfs [45].Therefore, there are three bands of strong mesoscale eddy variability off the coast of Central America from the Tehuantepec, Papagayo, and Panama Gulfs [42].Itoh and Yasuda [15] revealed that many eddies were generated along the Kuroshio Extension.Anticyclonic and cyclonic eddies shed off from the northern and southern Kuroshio Extension, respectively.On the lee side of Hawaiian Islands and in the STCC region, eddies are generated by different mechanisms.Qiu [35] suggested that eddies in the STCC region could be generated by the baroclinic instability of the vertical current shear.On the lee side of Hawaiian Islands, eddies are generated mainly due to wind curls and wakes of the Hawaiian Islands [46].Because of the nearly westward propagation of the eddies, the Kuroshio and the Oyashio currents play a role in terminating eddies [9].Contrarily, eddies are seldom observed in vast areas around the center of the eastern part of the North Pacific Subarctic Gyre (Figure 8a), because the eddy sizes there are too small to be resolved [8].From the kinetic energy distribution of currents in the study area derived from altimetry ADT (Figure 9), a relatively small kinetic energy is found in the North Pacific subarctic Gyre.Based on the fact that the formation of eddies requires energy supplies provided by the unstable oceanic kinetics, this may be the reason that eddies rarely form in this region.

Eddy Path
Chelton et al. [9] revealed that a striking feature of eddy trajectories in global oceans is obviously traveling westward with small deflections.In this study the results of eddy tracking similarly show that eddies propagate westward with various meridional deflections.As shown in Figure 10, we apply the linear regression method to eddy trajectories and define the angle of azimuth in relation to the east direction.The map of eddy azimuth shows that most eddies propagate nearly westward.Nevertheless, eddies propagate southwestward near the coast of eastern Japan and northeastward around the Subarctic Front.Similar results were also found in the study of Itoh and Yasuda [15].In the Kuroshio Extension and the STCC regions, the trajectories of eddies exhibit a strong westward propagating tendency (180°).Eddies are generated by the frontal instability associated with the STCC and can be carried by the westward propagating Rossby waves [47].Eddies in the southern Japan region propagate eastward, which means that they may be affected by the strong eastward Kuroshio Extension.Although this region is around the large meander in the Kuroshio south of central Japan, we could detect mesoscale eddies but not the Kuroshio meanders.Near the eastern coast of Taiwan and Luzon, eddies propagating with about 100° azimuth might be affected by the northward Kuroshio.Moreover, Eddies purely propagate along the coast around the Alaska Stream.Eddies are strongly influenced by the Alaskan Stream system in this region.The Gulf of Alaska turns southwestward following the self-break and forms the Alaskan Stream, which may carry eddies [37].Xiu et al. [48] demonstrated that eddies detached and propagated away from the coast along the Alaskan Stream.On the other hands, the propagation of anticyclonic and cyclonic eddies differs more around the California currents.As noted in previous studies [9,41], eddies show a tendency for westward propagation with cyclonic and anticyclonic eddies exhibiting a poleward and equatorward deflection.This suggests that nonlinear eddy propagation may be influenced by the β-effect and self-advection [49].
Furthermore, circular histograms of the mean azimuth in seven active regions are used to display the difference (Figure 11).Eddies propagate southwestward with 210° and northeastward with 30° nearly parallel to the coastal line of eastern Japan around the Oyashio currents.In the North Equatorial Current region, regardless of whether the eddies are cyclonic or anticyclonic, the northwestward propagation near an 170° azimuth may be affected by the North Equatorial Current.In the Alaska Current region, eddies are strongly influenced by the Alaskan Stream that propagates westward at 208° azimuth for anticyclonic eddies and 203° for cyclonic eddies.The differences of deflection between anticyclonic and cyclonic eddies may be affected by the combined influence of the β effect and the advection effect of the Alaskan Stream.Therefore, the equatorward deflections of anticyclonic eddies are more noticeable than the deflections exhibited by cyclonic eddies.In the regions of the California Current and the Northeastern Tropical Pacific, eddies propagate westward with cyclonic eddies exhibiting a slight deflection toward to the pole (165°); and anticyclonic eddies tending toward to the equator (197°).However, this kind of distribution of the eddy pathways in the other regions is not remarkable.The eddy trajectories may be affected mainly by eddy-eddy and eddy-mean flow interactions.Based on the aforementioned results, the effects of the background current may not be crucial for eddy trajectories around the California currents compared with other regions in the North Pacific.

Conclusions
The weekly SLA data derived from satellite altimetry are used to investigate the mesoscale eddies in the North Pacific.To extract eddies from the ocean surface, we use the free-threshold eddy identification algorithm.This study considers eddies with lifespans longer than 4 weeks and investigates their characters.According to results of analyses of the SLA between November 1992 and December 2012, histograms of radius, propagation speed, amplitude, and eddy kinetic energy are approximated by the Rayleigh distribution.From maps of EKE and amplitude, eddies are more active in the regions of Kuroshio Extension, STCC, the lee side of the Hawaiian Islands and the northeastern tropical Pacific.Although eddies can generate and terminate in most regions, eddies are rarely observed around the center of the eastern part of the North Pacific Subarctic Gyre.The radius, travel distance, and speed of eddies mainly decrease with increase in latitude.On the whole, anticyclonic eddies possess greater radii, amplitudes, and lifespans than cyclonic eddies, but the EKE of cyclonic eddies is not smaller than that of anticyclonic eddies.Lifespans of eddies which can live up to 3 years in the Alaska Current region are the greatest in the North Pacific.The largest amplitude and EKE are revealed in the Kuroshio Extension and Northeastern Tropical Pacific regions, hence these regions possess a bigger potential to carry water mass and influence the large-scale flux of momentum, heat, and salt.Regarding of whether the eddies are found in the regions of California Current or the Northeastern Tropical Pacific, most eddies propagate westward with slight equatorward and poleward deflections for cyclonic and anti-cyclonic eddies, respectively.This result suggests that eddy propagation may be influenced by the β-effect.However, this divergence of the eddy pathways is not remarkable in other regions.Around the Oyashio currents, most eddies propagate nearly parallel to the coastal line of eastern Japan.This study gives a statistical analysis and premilinary results of eddies in the North Pacific.Nevertheless, more detailed studies of the dynamic processes and eddy generation mechanisms constitute a challenge for future research.

Figure 1 .
Figure 1.Schematic current patterns in the North Pacific.Geostrophic currents are calculated from altimeter absolute dynamic topography data and averaged over a period from 1992 to 2012.Background color is the speed of geostrophic currents.

Figure 2 .
Figure 2. Sea level anomaly (SLA) and eddies on 14 July 2010.Red and blue dots are the centers of anticyclonic and cyclonic eddies, respectively.The eddy edges are enclosed by red (blue) lines.(a) Continuous contour lines are displayed every 4 cm; (b) Zoom in of the region of black box in (a).Continuous contour lines are displayed every 1 cm.
shows the histogram of eddy characteristics in the North Pacific, and we use the probability density function of the Rayleigh distribution to fit them.The results reveal that histograms of radius, propagation speed, amplitude, and eddy kinetic energy are approximated by the Rayleigh distribution.

Figure 3 .
Figure 3. Histogram of eddy radius (a); lifespan (b); travel distance (c); eddy kinetic energy (EKE) (d); amplitude (e); and speed (f).Red and blue bars stand for anticyclonic and cyclonic eddies, respectively.The dotted line is the fitting line of a Rayleigh distribution.
also showed that water had been transported for hundreds of kilometers by eddies in the Gulf of Alaska region.The largest amplitude is revealed in the Kuroshio Extension and Northeastern Tropical Pacific regions.The smallest lifespan and amplitude appear in the North Equatorial Current region.Although the greatest EKE in the western North Pacific appears in the Northeastern Tropical Pacific region, the number of eddy in the Northeastern Tropical Pacific region (about 100/degree 2 ) is smaller than that in the California Current region (about 250/degree 2 ) in the last two decades.Eddies in the Kuroshio Extension and the Northeastern Tropical Pacific regions possess a bigger potential to carry water mass and influence the large-scale flux of momentum, heat, and salt, because they possess greater a EKE in those regions.

Figure 5 .
Figure 5. Map of active regions in the North Pacific.

Figure 6 .
Figure 6.Boxplot of eddy kinematic characteristics in the active regions in Figure 5.The green line is the median value and the box represents the interquartile range.The whiskers are lines extending from each end of the box to show the extent of the rest of the data.The whiskers' length is 1.5 times that of the interquartile range.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.(a) Radius; (b) Lifespan; (c) Distance; (d) EKE; (e) Amplitude; (f) Speed.Red and blue represent the anticyclonic and cyclonic eddies, respectively.

Figure 7 .
Figure 7. Sketch (a) and histogram (b) of shape error.The closed contour of the eddy edge and the ideal circle are represented by the red and black lines, respectively.Area deviations are the sum of the shaded regions, and the shape area is the area deviation expressed as the percentage of the area of the fitted circle.

Figure 8 .
Figure 8. Statistics for eddies with lifespan more than or equal to 4 weeks showing the numbers of (a) eddy generations and (b) eddy terminations for each 1° × 1° region over a 20-year period.

Figure 9 .
Figure 9. Spatial distribution of kinetic energy in the North Pacific calculated from altimeter absolute dynamic topography (ADT) data and averaged over the 20-year period.

Figure 10 .
Figure 10.The maps of the mean azimuth relative to the eastern anticyclonic (a) and cyclonic eddies (b).

Figure 11 .
Figure 11.The propagation characteristics of cyclonic (blue) and anticyclonic eddies (red) with lifespans longer than or equal to 4 weeks.Circular histogram of the mean azimuth relative to the east in seven active regions in Figure 5.

Table 1 .
Statistics of eddy kinematic properties in seven regions of the North Pacific a from 1992 to 2012 (Mean ± One Standard Deviation).