Chlorophyll Rings around Ocean Eddies in the North Pacific

Chlorophyll rings (CRs) are defined as elevated chlorophyll along eddy peripheries and have been observed in anticyclonic oceanic eddies occasionally. This study presents observations of CRs around both anticyclonic and cyclonic eddies from a large observational data set. An innovative algorithm is developed to identify CRs from satellite observations of sea level anomalies and near-surface chlorophyll concentration in the North Pacific Ocean between 2003 and 2010. The results show that only 1% of mesoscale eddies are associated with CRs, which implies the CRs are not ubiquitous. We propose two potential generation mechanisms for CRs: horizontal advection and wind-current interaction. The former dominates the formation of about two-thirds of the CRs. The CRs associated with both cyclones and anticyclones represents an important contribution to better understanding of mesoscale physical/biological coupled phenomena.


Results
The study is primarily based on chlorophyll concentration data and sea surface geostrophic currents velocities which are used to identify CRs. In addition, sea surface temperature (SST), and surface vector winds are also employed in the discussions of the potential mechanisms for CRs.
A total of 121,846 cyclonic and 119,534 anticyclonic realizations with lifetimes greater than 4 weeks are matched with concurrent chlorophyll concentration data between 2003 and 2010. Each of these realizations meets the criteria of more than 30% chlorophyll concentration coverage within the region inscribed by a circle whose radius is 2.5 times that of the eddy. Based on the criteria for CRs identification (see Methods for details), 1,286 cyclonic eddies and 1,506 anticyclonic eddies are associated with CRs. These account for about 1% of all the eddies detected. Anticyclones have a 17% higher likelihood of being surrounded by a ring when compared to cyclones.
To quantify the significance of CRs at a basin scale, area-integrated chlorophyll concentrations are calculated both over the area of all CRs and over the total area of the North Pacific basin over one climatological year. The averaged chlorophyll concentration over the unit area within CRs is 0.47 ± 0.07 mg m −3 , which is about twice of that over the whole North Pacific basin, 0.23 ± 0.01 mg m −3 . Considering that the total area covered by CRs   characterized by either high chlorophyll concentration (east and west boundary) or large eddy kinetic energy (Kuroshio extension). The lifetime of an individual eddy is separated into several stages, defined as proportion of total lifetime (ranging from 0 to 1) 2,21 : formation stage (0~0.1), intensification stage (0.1~0.3), mature stage (0.3~0.8) and aged stage (0.8~1). Specific life stages are counted for CRs within eddies, whose lifespans are longer than 10 weeks. This analysis reveals that CRs in eddies of both polarities are most common during the intensification and mature stage (Fig. 1b). In anticyclones, however, CRs occur more frequently in the late mature stage than that during eddy intensification. In contrast, CRs develop more frequently during eddy intensification than mature stage in cyclones.
Spatial distributions of the ensemble mean chlorophyll concentration of CRs around anticyclonic eddies and cyclonic eddies are shown in Fig. 1c,e. Normalized chlorophyll concentrations are presented in the figures, which

Discussion
Two potential mechanisms are proposed for the CR generation: horizontal advection and wind-current interaction.
(1) Horizontal advection of chlorophyll in the ambient water of eddies is reported to be one of the dominant mechanisms by which eddies influence chlorophyll distribution 2,7 . The generation of CRs by horizontal advection can be seen in the example shown in Fig. 2a~d where an anticyclonic eddy in the eastern boundary of the North Pacific is observed to generate a CR. The rotating flow of the eddy results in a chlorophyll "tongue" south of the eddy on Aug 16, 2006. The chlorophyll tongue moves to the north with the development of eddy. Eventually higher chlorophyll concentration appears as a ring structure. The chlorophyll concentration within the region 0.5 Ls~1.5 Ls of the anticyclonic eddy is enhanced by 0.12~0.21 mg m −3 during the formation, relative to the climatological value. The SST shows additional evidence of horizontal advection (Fig. 2e~h). A cold tongue propagates from the shore to the offshore while the anticyclonic eddy rotates. Consequently, a ring structure with cold water is developed along the periphery of the eddy. The anticyclonic eddy survives for 44 weeks and the analysis of the evolution of the CRs suggests that the ring occurs at a normalized lifetime of 0.2~0.35, defined here as the eddy intensification stage.
Analogously, a CR is formed around a cyclonic eddy at the boundary of the Extension of Kuroshio in Fig. 3a~d. The local chlorophyll concentrations within this CR increase by about 0.01 mg m −3 . In SST images Comparing the averaged chlorophyll concentration along the eddies' peripheries (between 0.5 Ls and 1.5 Ls) and outside the eddies (greater than 1.5 Ls), we distinguish whether CRs are generated by the horizontal transport of chlorophyll. If the averaged chlorophyll concentration outside the eddies is higher than that within a CR along the periphery, the CR is judged to be associated with the contribution of horizontal advection. 2123 CRs (about 66%) are found to be formed by this mechanism.
(2) Wind-current interaction via eddy-induced Ekman pumping could generate vertical transport, possibly resulting in CRs observed around eddies. In regions of intense wind stress and relatively strong horizontal current gradients (e.g., at the edges of mesoscale features), Ekman-induced vertical velocities are generated not only by the curl of the wind stress but also by the interaction of the surface wind stress with the geostrophic relative vorticity gradient 19,[22][23][24][25] . With the Ekman transport modified by the surface geostrophic vorticity ζ g , the total eddy-induced Ekman pumping W tot is where ζ g is the vertical component of the geostrophic relative vorticity, ρ 0 =1020 kg m −3 is the (assumed constant) density of sea surface water, f is the Coriolis parameter, τ x and τ y are zonal and meridional components of the surface wind stress  τ , respectively. Small magnitude of the meridional derivative of the Coriolis parameter contribution to the total Ekman pumping in (1) is ignored. Surface stress curl-induced Ekman pumping, W c (the first term on the right side of (1)), is determined by the wind tress curl because that ζ + f g in the expression of W c is always positive for mesoscale eddies: when the wind stress curl is positive (negative), it generates the upwelling (downwelling). Vorticity gradient-induced Ekman pumping, ζ W (the second term on the right side of (1)), generates a dipole of Ekman upwelling and downwelling within the interior of ocean eddies. Figure 4a shows a CR in the anticyclonic eddy on Aug 31, 2005. The interaction of intense wind stress with sea surface currents results in upward vertical transport (1.5 m day −1 ) in the east part of the eddy (Fig. 4b). The portion resulting from the curl of the surface stress, in Fig. 4c, generates upwelling within the anticyclone, while the contribution resulting the interaction of the surface stress with surface current vorticity gradient induces a dipole of upwelling and downwelling in the eddy, in Fig. 4d.
Similarly, eddy-induced Ekman pumping could also generate along the periphery of cyclonic eddy (Fig. 5). Upward vertical transport (0.6 m day −1 ) resulting from the total Ekman pumping appeares in the southeastern side of eddy edge. Additionally, the eddy is located south of a high-chlorophyll region, hence horizontal advection can also play a role in the CR formation observed. This is the first endeavor to investigate CRs systematically. It shows that the CRs can occur in both cyclonic and anticyclonic eddies. An innovative algorithm to identify CRs from a large set of eddies is developed. CRs are more likely to be observed along the periphery of anticyclones when compared to cyclones. Two potential CR generation mechanisms are proposed: horizontal advection which transports the chlorophyll from adjacent areas; while vertical transports caused by eddy-induced Ekman pumping can bring nutrient and phytoplankton into the surface layer. Satellite observations and a reanalysis numerical product are applied to explore mechanisms that form observed CRs. We found that horizontal advection explains the formation of about two-thirds of CRs identified in the North Pacific.
Two potential mechanisms are proposed in this study, however, there are also various mechanisms of vertical transport of nutrients, chlorophyll and phytoplankton below the surface, which might be at play at any given time, for examples: local upwelling at submesoscale induced by symmetric instability could bring waters rich in nutrients into the euphotic layer 15,26 ; internal waves could also uplift isopycnals high in nutrients within eddies 27 ; ageostrophic secondary circulation could contribute to the vertical transport upward and downward as well 25,[27][28][29][30] . Currently, there is still a challenge as the spatial-temporal resolution of observations is not available to resolve the submesoscale processes. Because of the interdependencies among various processes and scales, it is difficult to ascribe the vertical advection nutrients to individual mechanisms 30 , which requires further study in future.

Methods
Observational Data. The data used in the present study include several observational data: satellite remote sensing ocean color data (chlorophyll concentration), sea surface geostrophic currents velocities, sea surface temperature and wind vectors.
Eddies identification is based on weekly-mean sea surface geostrophic currents with a spatial resolution of 1/4°×1/4° which are obtained from AVISO (www.aviso.oceanobs.com). An automated algorithm based on the geometry of the velocity vectors of the flow field 31 is used to identify and track mesoscale eddies. Our research is based on eddies with lifetimes of four weeks and longer. Chlorophyll concentration data used in this paper are downloaded from http://wiki.icess.ucsb.edu/measures/ index.php/GSM. It is a merged product of SeaWiFS, Meris, and MODIS-Aqua with a resolution of 9 km every 8 days, generated using the Garver-Siegel-Maritorena (GSM) model 32,33 .
SST data are used to illustrate the developments of eddies and to present the process of advection. The SST data used are provided by Group for High Resolution Sea Surface Temperature (GHRSST) Level 4 sea temperature analysis production based on a global 0.01 degree grid with the temporal resolution of 1 day. The Multiscale Ultrahigh Resolution (MUR) L4 analysis is based on SST observations from AMSRE, MODIS, WindSat, AVHRR and in situ observations from NASA.
Wind data used to calculate the nonlinear Ekman pumping combined with geostrophic currents. Wind data are the Cross-Calibrated Multi-Platform (CCMP) gridded surface vector winds which are produced by satellite, moored buoy, and model wind data. This daily production covers the global regions with the spatial resolution of 0.25°.
Chlorophyll rings identification. Chlorophyll observations are linearly interpolated onto a normalized grid spanning +/−2.5Ls allowing missing data to be filled (Ls denotes the eddy radius). The normalized, eddy-centric chlorophyll concentrations are then transformed from a Cartesian grid (Fig. 6a) into polar coordinates (Fig. 6b). The average chlorophyll concentrations in different radial bins are calculated (Fig. 6c) and if the local extrema of averaged chlorophyll concentration are located between 0.5 Ls and 1.5 Ls, the feature is identified as a possible CR. In order to avoid that the chlorophyll concentration is extremely high locally, a similar process is applied to chlorophyll concentration along the 16 directions in the polar coordinates which are evenly distributed from 0° to 360°, checking whether the peak of chlorophyll concentration is between 0.5 Ls and 1.5 Ls in more than 9 directions. To all the possible rings, a threshold is applied to the range of chlorophyll between the maximum and minimum of 3 times of the minimum value. Figure 6d-f demonstrate the application of the method in a cyclonic eddy on Apr 27, 2005. Eddies identification. An automated algorithm based on the geometry of the velocity vectors of the flow field 31 is used to identify and track mesoscale eddies based on the geostrophic current velocity which is calculated from the dataset of Merged Sea Level Anomalies (MSLA) acquired from AVISO. Eddy centers are determined by the spatial characters of the velocity fields, eddy sizes are computed from closed contours of the stream function field, and eddy tracks are retrieved by comparing the distribution of eddy centers at successive time steps.
Each eddy is matched with chlorophyll concentration data if it has a minimum lifetime of 4 weeks and the proportion of missing chlorophyll data do not exceed 70% for total number of pixels within an area defined as 2.5 time the radius scale of the eddy.