Vertical structure of mesoscale eddies in the eastern South Pacific

JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 116, C11025, doi:10.1029/2011JC007134, 2011
Vertical structure of mesoscale eddies in the eastern South Pacific
Ocean: A composite analysis from altimetry and Argo
profiling floats
Alexis Chaigneau,1,2,3 Marie Le Texier,4 Gérard Eldin,3 Carmen Grados,2
and Oscar Pizarro5
Received 16 March 2011; revised 7 September 2011; accepted 7 September 2011; published 17 November 2011.
[1] The mean vertical structure of mesoscale eddies in the Peru‐Chile Current System is
investigated by combining the historical records of Argo float profiles and satellite
altimetry data. A composite average of 420 (526) profiles acquired by Argo floats that
surfaced into cyclonic (anticyclonic) mesoscale eddies allowed constructing the mean
three‐dimensional eddy structure of the eastern South Pacific Ocean. Key differences in
their thermohaline vertical structure were revealed. The core of cyclonic eddies (CEs)
is centered at ∼150 m depth within the 25.2–26.0 kg m−3 potential density layer
corresponding to the thermocline. In contrast, the core of the anticyclonic eddies (AEs) is
located below the thermocline at ∼400 m depth impacting the 26.0–26.8 kg m−3 density
layer. This difference was attributed to the mechanisms involved in the eddy formation.
While intrathermocline CEs would be formed by instabilities of the surface equatorward
coastal currents, the subthermocline AEs are likely to be shed by the subsurface
poleward Peru‐Chile Undercurrent. In the eddy core, maximum temperature and salinity
anomalies are of ±1°C and ±0.1, with positive (negative) values for AEs (CEs). This
study also provides new insight into the potential impact of mesoscale eddies for the
cross‐shore transport of heat and salt in the eastern South Pacific. Considering only the
fraction of the water column associated with the fluid trapped within the eddies, each CE
and AE has a typical volume anomaly flux of ∼0.1 Sv and yields to a heat and salt
transport anomaly of ±1–3 × 1011 W and ±3–8 × 103 kg s−1, respectively.
Citation: Chaigneau, A., M. Le Texier, G. Eldin, C. Grados, and O. Pizarro (2011), Vertical structure of mesoscale eddies in the
eastern South Pacific Ocean: A composite analysis from altimetry and Argo profiling floats, J. Geophys. Res., 116, C11025,
doi:10.1029/2011JC007134.
1. Introduction
[2] The Peru‐Chile Current System (PCCS), also known as
the Humboldt Current System, is relatively complex, exhibiting several surface and subsurface currents. Its dynamics is
principally controlled by the atmospheric South Pacific
Anticyclone through Sverdrup dynamics. In the surface layers, the eastward flowing South Pacific Current feeds the
Chile‐Peru Current (CPC) between 30°S and 40°S, forming
the eastern branch of the South Pacific subtropical gyre [Strub
et al., 1998; Chaigneau and Pizarro, 2005a, 2005b]. Near the
South American coast, the persistent equatorward winds
1
Laboratoire d’Océanographie et de Climatologie: Expérimentation et
Analyse Numérique, UMR 7159, IRD, CNRS, UPMC, MNHN, Paris,
France.
2
Instituto del Mar de Perú, Callao, Peru.
3
Laboratoire d’Études en Géophysique et Océanographie Spatiale,
UMR 5566, IRD, CNRS, CNES, UPS, Toulouse, France.
4
Département de Formation en Hydraulique et Mécanique des Fluides,
ENSEEIHT, Toulouse, France.
5
Department of Geophysics, UDEC, Concepcion, Chile.
Copyright 2011 by the American Geophysical Union.
0148‐0227/11/2011JC007134
drive upwelling cells leading to the highest biological productivity of the world ocean in terms of fish [Chavez et al.,
2008]. Over the continental shelf, the upwelling of relatively cold deep water also gives rise to intense thermal fronts
which separate, over short distances, cold coastal water from
warmer and saltier subtropical water of the offshore ocean
[Wyrtki, 1967; Strub et al., 1998]. By geostrophic adjustment,
this relatively strong cross‐shore density gradient in the surface layers reinforces the equatorward Chile‐Coastal Current
and Peru‐Coastal Current (Figure 1a). The PCCS is also
characterized by two major subsurface poleward currents
(Figure 1a): the Peru‐Chile Countercurrent and the Peru‐
Chile Undercurrent (PCU) which have been both traced back
to the equatorial current system [Lukas, 1986; Tsuchiya,
1985; Strub et al., 1998; Montes et al., 2010]. The PCU
which transports relatively warm and salty equatorial subsurface water from the eastern tropical Pacific to at least 48°S
[Silva and Neshyba, 1979] along the continental slope and
shelf is a major source of the coastal upwelling off Peru and
northern Chile [Huyer et al., 1987; Montes et al., 2010;
Albert et al., 2010]. This subsurface current also supplies
low‐oxygenated water that helps maintaining the oxygen
C11025
1 of 16
C11025
CHAIGNEAU ET AL.: EDDY VERTICAL STRUCTURE IN THE ESP
Figure 1. Spatiotemporal distribution of the 4179 Argo profiles used in this study and typical water mass
properties observed in the PCCS. (a) Position of the 4179 Argo profiles and schematic circulation. Solid
blue lines indicate surface currents, and dashed red lines show subsurface currents. PCC, Peru Coastal
Current; CCC, Chile Coastal Current; PCCC, Peru‐Chile Countercurrent; and PCU, Peru‐Chile Undercurrent. (b) Meridional variation of the number of Argo profiles in 1° latitude bands. (c) Monthly variations of
the number of Argo profiles. (d) Mean temperature‐salinity diagrams observed in the two subregions
delimited by a black dashed line in Figure 1a. Grey curves represent potential density anomalies (s in
kg m−3). SPESTMW, South Pacific Eastern Subtropical Mode Water; ESPIW, Eastern South Pacific Intermediate Water; ESSW, Equatorial Subsurface Water; and AAIW, Antarctic Intermediate Water.
2 of 16
C11025
C11025
CHAIGNEAU ET AL.: EDDY VERTICAL STRUCTURE IN THE ESP
minimum zone of the eastern South Pacific [Wyrtki, 1963;
Fuenzalida et al., 2009; Stramma et al., 2010].
[3] The oceanic circulation along the South American
coast is also characterized by energetic mesoscale structures,
the oceanic cyclonic and anticyclonic eddies (CEs and AEs,
respectively). Such eddies have been mostly observed from
altimeter data [Chelton et al., 2007; Chaigneau et al., 2008,
2009], but they also have a signature in color satellite
images [Correa‐Ramirez et al., 2007] or in surface‐drifter
trajectories [Chaigneau and Pizarro, 2005c]. In the PCCS,
CEs and AEs, which have a typical diameter of 150–
300 km, are principally formed near the South American
coast where they locally impact the heat and salt budgets
through lateral turbulent fluxes [Chaigneau and Pizarro,
2005b; Colbo and Weller, 2007]. Then, CEs and AEs,
which acquire a water mass structure typical of their formation region, propagate seaward with translation velocities of
few cm s−1 owing to a combination of mean flow advection
and self‐propagation. In this new environment, eddies
appear as anomalous water masses with surface or subsurface temperature and salinity anomalies [Johnson and
McTaggart, 2010] which are progressively redistributed to
surrounding water during eddy decaying phase [e.g., Swart
et al., 2008]. The westward propagation of CEs can also
extend the area of high biological productivity offshore by
both eddy chlorophyll advection and eddy nutrient pumping
[Correa‐Ramirez et al., 2007].
[4] Although the main horizontal structure and kinematic
properties of the PCCS eddies were investigated during the
last decade, very little is known about their vertical structure
and their impact on the heat and salt transports. Chaigneau
and Pizarro [2005c] have briefly examined the vertical
structure of a particular CE sampled by the World Ocean
Circulation Experiment (WOCE) P19 hydrographic section
along 88°W, showing that the subsurface CE core have typical temperature and salinity anomalies of around −1°C and
−0.1, respectively. More recently, Johnson and McTaggart
[2010] used Argo float profile data to characterize AEs
of the PCCS, showing that their core is located in the
subthermocline and contains anomalous signature of the
Equatorial Subsurface Water originating from the PCU.
Although this analysis provides a first insight into the
vertical characteristics of AEs in the PCCS, some open
questions remain. For instance, do CEs exhibit the same
vertical structure as AEs? What is the respective role of such
eddies on the volume, heat and salt transports?
[5] The main objective of this work is thus to answer these
questions and extend the work of Johnson and McTaggart
[2010]. To achieve this goal, we propose a blended analysis of satellite altimetry data and Argo float profiles allowing
the three‐dimensional reconstruction of composite eddies.
The obtained average eddy fields are used to characterize the
mean vertical structure of both eddy types and to estimate
their contribution to the volume, heat and salt transports. The
paper is organized as follows. In section 2 we describe the
study region and the two data sets (altimetry and Argo data)
used in this work. Eddy identification algorithm from satellite data is also briefly presented in this section as well as the
methodology used to classify the Argo profiles. The methodology used for the construction of the three‐dimensional
composite eddies is presented in section 3. The mean thermohaline properties of the composite CE and AE are
C11025
described in section 4, as well as the eddy function in the
offshore transport of volume, heat and salt. Section 4 also
deals with the meridional variation of the eddy vertical
structure between the northern and southern areas of the
study region. The observed differences between cyclonic and
anticyclonic eddies allow a discussion in section 5 of the
most probable origin of both eddy types. Finally in section 6,
we summarize the results.
2. Data and Methods
[6] The study region extends from 10°S to 30°S and from
the South American coast to 100°W (Figure 1a). The
northern boundary is located south of the equatorial front
which separates the relatively fresh equatorial water from
the saltier subtropical water [Fiedler and Talley, 2006].
Similarly, the southern margin is located to the north of the
subtropical front, separating the subtropical water from the
relatively cold and fresh subantarctic water [Tomczak and
Godfrey, 1994; Stramma et al., 1995; Chaigneau and
Pizarro, 2005a].
2.1. Altimeter‐Derived Sea Level Anomaly Data
and Eddy Identification
[7] The presence and position of mesoscale eddies in the
study domain are determined by analyzing sea level
anomaly (SLA) maps from the multisatellite AVISO product
(http://www.aviso.oceanobs.com) between December 2003
and October 2009. This gridded multimission altimeter
product, produced by Ssalto/Duacs and distributed by CLS–
Space Oceanography Division (Toulouse, France), provides
the best available spatiotemporal resolution for observing
mesoscale features [Le Traon and Dibarboure, 1999; Pascual
et al., 2006]. SLA data, computed relatively to a 7 year mean
(1993–1999), are weekly mapped onto a 1/3° Mercator grid
and then bilinearly interpolated by AVISO onto a 0.25° ×
0.25° longitude/latitude grid (see Appendix A.2 of Chelton
et al. [2011] for a detailed description of the procedure
used by AVISO to produce SLA maps). This weekly SLA
product is further linearly interpolated in time to obtain SLA
maps on a 1 day basis.
[8] In each daily SLA map, mesoscale eddies were identified using the method recently developed by Chaigneau
et al. [2009]. The CE (AE) detection algorithm first
involves searching for eddy centers associated with local
SLA minima (maxima) in a moving window of 1° × 1°. Then,
for each possible cyclonic (anticyclonic) center, the algorithm
looks for closed SLA contours with an increment (decrement)
of 10−3 m. The outermost closed SLA contour, embedding
only the considered center, corresponds to the eddy edge.
Given the accuracy of satellite measurements and AVISO
product [Le Traon and Ogor, 1998; Ducet et al., 2000;
Chelton and Schlax, 2003], we only consider here eddies
having SLA amplitudes higher than 2 cm.
[9] A similar algorithm, based on closed SLA contours, has
also been used recently to describe mesoscale eddy characteristics at global scales [Chelton et al., 2011]. This eddy
identification method has been shown to be more accurate
than the commonly used Okubo Weiss‐based method [e.g.,
Isern Fontanet et al., 2003]. In particular, it allows to considerably reduce the number of “false” detected eddies
[Chaigneau et al., 2008]. Souza et al. [2011] have also
3 of 16
C11025
CHAIGNEAU ET AL.: EDDY VERTICAL STRUCTURE IN THE ESP
shown, through the comparison of three automatic detection
algorithms in the South Atlantic, that the geometric criterion
used in this study exhibits a better performance, mainly in
terms of the eddy detection number, their lifetimes and
propagation velocities.
2.2. Argo Data Set
[10] The vertical structure of mesoscale eddies was
investigated using the autonomous CTD profiling floats of
the Argo program (http://www.coriolis.eu.org) available in
the study region. This data set spans the same temporal
period as the altimetry product (December 2003 to October
2009). The original data set consists of 93 distinct floats
which acquired 6050 CTD profiles during the studied period.
Argo data are automatically preprocessed and quality controlled by the Argo data centers [Wong et al., 2003; Böhme
and Send, 2005; Owens and Wong, 2009]. Only pressure
(P), temperature (T), and salinity (S) data flagged as good
(Argo quality flag 1) were retained in the analyses. However,
some of these preprocessed profiles were still suspicious.
Thus, for our analysis we only considered profiles for which
(1) the shallowest data are located between the surface and
10 m depth and the deepest acquisition is below 1000 m;
(2) the depth difference between two consecutive data does
not exceed a given limit (Dzlim) which depends on the
considered depth (Dzlim = 25 m for the 0–100 m layer;
Dzlim = 50 m for the 100–300 m layer; and Dzlim = 100 m
for the 300–1000 m layer); and (3) the number of data levels
in the 0–1000 m layer is higher than 30.
[11] Once selected, every profile was visually checked,
and those which presented a suspicious T/S diagram were
systematically discarded. The final database is composed of
4179 profiles, corresponding to 69% of the initial data set
(Figure 1a). The PCCS was rather homogeneously sampled
by Argo floats except around 25°S where fewer profiles
were acquired (Figure 1b). The temporal distribution shows
that the number of CTD profiles linearly increased from
December 2003 (date of the first Argo float deployment in
the study region) to June 2008 (Figure 1c) and dramatically
decreased afterward.
[12] The mean T/S diagrams for the northern (10°S–20°S)
and southern (20°S–30°S) regions of the study domain
highlight the main water masses present in the PCCS
(Figure 1d). South Pacific Eastern Subtropical Mode Water
fills the major part of the surface ocean, with T > 20°C and
S > 35. In subsurface, the relatively fresh and cold Eastern
South Pacific Intermediate Water (T = 12–13°C, S < 34.8,
s ≈ 26 kg m−3) is located at the base of the thermocline and
surmounts the colder and saltier Equatorial Subsurface
Water (T ≈ 8–11°C, s ≈ 26.5 kg m−3). Finally, below the
Equatorial Subsurface Water is found the Antarctic Intermediate Water associated with a local salinity minimum
(T < 7°C, S < 34.6, s > 27 kg m−3). Note that both the
Eastern South Pacific Intermediate Water and the Antarctic
Intermediate Water become saltier and slightly warmer along
their northward/northwestward pathway in the PCCS,
whereas the Equatorial Subsurface Water of equatorial origin
gets fresher and colder toward the south (Figure 1d). Readers
interested in more details on the water mass distribution in
the eastern South Pacific should refer for instance to Wyrtki
[1967], Strub et al. [1998], Blanco et al. [2001], Wong and
Johnson [2003], and Schneider et al. [2003].
C11025
[13] The 4179 retained T/S profiles were then linearly
interpolated onto 101 regularly spaced vertical levels every
10 m from the surface to 1000 m, assuming that the shallowest values of T and S were representative of surface
values. Derived quantities were estimated from these vertically gridded values of T, S, such as potential temperature
() and density (s). At each vertical level, the dynamic
height (DH) relative to the 1000 m reference depth was also
computed. We considered that a 1000 m level of no motion
should not strongly impact our results since the average flow
field at 1000 m depth, estimated from Argo drift data has a
magnitude weaker than 1 cm s−1 in the PCCS [Davis, 2005;
Katsumata and Yoshinari, 2010].
[14] Finally, to better describe the mesoscale perturbations
captured by the floats, anomalies of every property (′, S′,
s′ , DH′) were computed by removing climatological profiles. Here, the climatological T/S profiles were obtained by
interpolating the CSIRO Atlas of Regional Seas (CARS)
[Dunn and Ridgway, 2002; Ridgway et al., 2002] to floats’
positions and times. Computed anomalies are expected to be
relatively weak for CTD profiles acquired outside eddies
and relatively strong inside eddies.
2.3. Classification of the Argo Profiles
[15] The 4179 CTD profiles were classified among three
distinct categories depending whether the floats surfaced
outside an eddy (OE) or inside a CE or AE identified from
daily SLA maps. Figure 2a shows the result of the classification for a given day in November 2005 with two floats
surfacing in CEs, one in an AE and two OEs. Over the entire
study period (December 2003 to October 2009), we obtained
a total of 420 profiles inside CEs, 526 inside AEs and
3233 OEs (Figures 2c–2e). Thus, AEs were slightly more
sampled than CEs and ∼23% of the profiles surfaced into
an eddy. This proportion is in agreement with the results of
Chaigneau et al. [2008], who estimated that mesoscale
structures cover ∼25% of the PCCS area.
[16] It is important to note that, apart from potential errors
in the automated eddy‐edge identification, three main factors
could cause errors in the classification of the Argo floats.
First, the accuracy of the CTD profiles locations depends on
the precision of the Argos positioning system which varies
from 150 m for a Class 3 Argos location to 1000 m for a
Class 1 Argos location. Second, the exact surfacing position
assigned to a CTD profile cannot be determined accurately
since an Argos positioning satellite does not necessarily fly
over a float immediately after it surfaces. A typical lag of
1 hour, during which Argo floats are advected by surface
currents, can be expected between the exact surfacing
location and the first reported position [Lebedev et al.,
2007]. However, considering maximum surface currents of
∼10 cm s−1 in the PCCS [Chaigneau and Pizarro, 2005b],
the error between reported CTD position and surfacing
position is relatively small (<500 m). Third, eddies were
identified from daily interpolated SLA maps and not at the
exact surfacing times reported by the Argo floats. On average, the 4179 surfacing times were recorded 1.3 h ahead the
associated SLA maps. Then, since long‐lived eddies propagate westward with typical velocities of 2–7 cm s−1 in the
study region [Chaigneau and Pizarro, 2005c; Chaigneau
et al., 2009], an additional small error of 100–300 m is
considered. We estimate that the total error introduced by
4 of 16
C11025
CHAIGNEAU ET AL.: EDDY VERTICAL STRUCTURE IN THE ESP
C11025
Figure 2. An illustrative example of the eddy detection algorithm and classification of the Argo profiles.
(a) Color shading corresponds to sea level anomalies (in cm) for the map of 15 November 2005, whereas the
edges of the automatically identified cyclonic and anticyclonic eddies are indicated by blue and red contours, respectively; Argo floats that surfaced on 15 November 2005 are shown as black dots. (b) Illustration
of the floats’ positions (M1 and M2) relative to the corresponding eddy centers (C1 and C2); this eddy‐
centered referential (Dx, Dy) is used to construct the composite eddies through objective interpolation
(see Figure 3 and text for details). (c) Position of the 420 Argo profiles that surfaced inside cyclonic
eddies between December 2003 and October 2009. (d) Same as Figure 2c but for the 526 profiles surfacing
inside anticyclonic eddies. (e) Same as Figure 2c but for the 3233 profiles surfacing outside eddies.
synchronization biases and positioning errors is of ∼1 km
on average, but can probably reach a distance as large as
∼5 km. This error can result in an erroneous classification of
a few Argo floats, in particular when they are close to the
identified eddy edges. However, only 0.1% (3.6%, respectively) of the CTD profile positions is closer than 1 km
(5 km) from a detected eddy edge. It is also worth to note that
the uncertainty of the eddy center and edge location determined from the SLA maps are probably much larger than the
synchronization biases and positioning errors of the Argo
floats.
3. Construction of Cyclonic and Anticyclonic
Composite Eddies
3.1. Objective Mapping of Eddy Properties
[17] CTD profiles surfacing into CEs or AEs can be
localized relative to their corresponding eddy centers identified as local maxima/minima in SLA (Figure 2b). Assuming first that all the CEs or AEs of the PCCS exhibit similar
3‐D structures, the spatial distribution of a given property (or
anomaly) inside CEs and AEs is thus investigated by a
composite analysis using a coordinate system (Dx, Dy) in
which each CTD profile is located respectively to the corresponding eddy center (Dx = Dy = 0). For instance, solid
dots in Figures 3a–3d show the temperature and salinity
anomaly distributions as a function of the distance (Dx, Dy)
to the eddy center. To illustrate the method, we chose a level
of 150 m for CEs (Figures 3a, 3c, and 3e) and 400 m for AEs
(Figures 3b, 3d, and 3f) which correspond to the vertical
location of the eddy cores (see section 4).
[18] The profiles inside the CEs and AEs are rather
homogeneously distributed around the composite eddy center. For example, at 150 m inside the CEs, the mean temperature and salinity anomalies obtained from the selected
420 Argo profiles are of −0.40°C and −0.05, respectively
(Table 1). At 400 m inside the AEs, the mean values obtained
are of 0.48°C and 0.04, respectively (Table 1). Although the
variability around these mean values is large, around 70%
(80%, respectively) of the 420 (526) ′ and S′ values are
negative (positive) for CEs (AEs). The mean ′ and S′ values
and the percentage of negative/positive values are strongly
different from the ones computed from the 3233 profiles
located OEs (Table 1). The nonparametric Mann‐Whitney
U‐test used to compare the different anomaly distributions
[Scherrer, 2007], confirmed that anomalies computed inside
eddies were significantly different (p < 0.05) from the ones
obtained OEs. Unlike the parametric t test, this nonparametric test makes no assumptions about the distribution of
the data (e.g., normality and equality of variance).
[19] At each depth level, anomalies that are more than
3 times the interquartile range from either the first or third
5 of 16
C11025
CHAIGNEAU ET AL.: EDDY VERTICAL STRUCTURE IN THE ESP
C11025
Figure 3. Objective interpolation, onto a regular grid (10 × 10 km2) of (a, b) potential temperature
anomalies and (c, d) salinity anomalies. Anomalies are shown at 150 m depth for cyclonic eddies
(Figures 3a and 3c) and 400 m depth for anticyclonic (Figures 3b and 3d) eddies. Solid dots in
Figures 3a–3d represent the anomalies estimated from Argo profiles in the eddy‐centered referential,
whereas color shadings correspond to the results of the objective interpolation. (e, f) Objectively interpolated dynamic height anomaly at 150 and 400 m depths relative to 1000 m (black contours) and horizontal
geostrophic speed (color shading, in cm s−1), respectively. White dots in Figures 3e and 3f show composite
eddy edges identified from the automatic algorithm, whereas black lines show eddy edges after fitting to an
ellipse (see text for details). Anomalies were computed relatively to the CSIRO Atlas of Regional Seas
(CARS) climatology interpolated to the positions and times of the Argo floats.
quartiles were considered as outliers and were discarded. The
remaining properties were then objectively mapped, assuming an isotropic Gaussian covariance decorrelation scale of
100 km. This scale allows filtering out the undesired small‐
scale variability and approximately corresponds to the mean
eddy radius previously observed in the study region
[Chaigneau et al., 2008, 2009]. Indeed from altimetry data,
the mean radius of the 971 eddies sampled by the profiling
floats is of 122 ± 30 km. Using the objective interpolation,
properties were mapped onto a regular 10 × 10 km2 grid.
This grid spacing was chosen since in the new referential,
80% of the minimum distances between two profiles are less
than 10 km (Figures 3a–3d). Color shading in Figures 3a–3d
shows the objectively mapped ′ and S′ at 150 m and 400 m
depth. At these levels, maximum ′ of 0.5°C–1°C are
observed near the composite center at Dx = Dy = 0. Although
negative (positive) anomalies persist over the interpolation
domain for CEs (AEs), the anomalies strongly weaken at
6 of 16
C11025
C11025
CHAIGNEAU ET AL.: EDDY VERTICAL STRUCTURE IN THE ESP
Table 1. Seawater Property Anomalies ′ and S′ Obtained From the 420 (526) Argo Floats That Surfaced Inside CEs (AEs)a
150 m Depth
N
′ (°C)
′ /max
′ (°C)
min
S′
′ /Smax
′
Smin
400 m Depth
Cyclonic Eddy
Outside Eddies
Anticyclonic Eddy
Outside Eddies
420
−0.40 ± 0.67
−2.96/1.88 (75%)
−0.05 ± 0.13
−0.61/0.35 (63%)
3233
0.10 ± 0.73
−2.75/2.85 (46%)
0.01 ± 0.14
−0.51/0.54 (48%)
526
0.48 ± 0.55
−0.69/2.58 (84%)
0.03 ± 0.05
−0.09/0.22 (76%)
3233
0.04 ± 0.33
−1.18/1.36 (53%)
0.00 ± 0.03
−0.12/0.13 (47%)
a
These values are compared to those obtained from the 3233 floats located outside eddies. A depth of 150 m (400 m) was chosen for CEs (AEs) since it
corresponds to the depth where density anomalies s′ are maximum. Numbers for ′ and S′ indicate the average ±1 standard deviation. Percentages
indicated correspond to the percentage of negative (positive) values observed at 150 m (400 m) depth in CEs (AEs) and OEs.
50–100 km from the eddy center. Thus at these levels, the
eddy core extends less than 100 km from the eddy center.
At 150 m depth in the core of the CE, the maximum
salinity anomalies are centered slightly more northward
than the associated temperature anomalies (Figure 3c). This
difference is induced by the objective mapping applied to
the much more variable salinity field (Figure 3c).
3.2. Horizontal Extent of the Composite Eddy Cores
[20] Figures 3e and 3f show the composite DH′ (black
contours) at 150 m and 400 m depth associated with CEs and
AEs and the corresponding geostrophic velocities computed
from the DH′ slopes (quivers and color shading). At each
vertical level, the center of the composite eddies correspond
to the local minimum/maximum of DH′ closest to the grid
center (Figures 3e and 3f). To determine the horizontal extent
of the eddy core, we averaged the geostrophic swirl velocities
along closed DH′ contours embedding the eddy center. The
closed DH′ contour associated with the strongest average
swirl velocity corresponds to the composite eddy‐core edge
(white dots in Figures 3e and 3f). To avoid spurious shapes
and vertical discontinuity, an ellipse is then least squares
fitted on this contour (black heavy ellipse in Figures 3e
and 3f). By definition, the composite eddy‐core edge coincides with a local maximum in geostrophic speed. This
maximum value can be highly variable from eddy to eddy
since it is proportional to eddy amplitudes for similar radii.
At 150 m depth (400 m, respectively), the maximum swirl
velocity value is of ∼7.8 cm s−1 (6.9 cm s−1) for the composite CE (AE). At these depths and considering only the
profiles located inside the ellipse, the mean temperature and
salinity anomalies of the eddy cores are on the order of
±0.8°C and ±0.05, respectively.
3.3. Vertical Extent of the Trapped Fluid
[21] One of the main objectives of this study is to estimate
the volume, heat and salt transports associated with the
PCCS eddies (see section 4.3). However, the water mass
anomalies in an eddy can only be maintained if the water
mass in the eddy is trapped for a considerable part, preventing surrounding water to enter the eddy when the eddy
is translated [Flierl, 1981; van Aken et al., 2003]. Thus, in
order to not overestimate the lateral transports we must only
consider the part of the water column which is effectively
trapped and transported by the eddies. On the basis of the
suggestions made by Flierl [1981], the amount of water
trapped inside a ring depends on the ratio of its drift velocity
to its tangential velocity. This ratio, which arises from the
comparison of nonlinear advection and acceleration, pro-
vides a measure of nonlinearity stating that the eddy
dynamics is nonlinear when this ratio exceeds 1 and maintains a coherent structure as it propagates [e.g., Flierl, 1981;
Chelton et al., 2007, 2011]. Thus, the composite eddies
were tracked with depth as long as their rotational speed
exceeded their translation speed. The rotational speeds of
our composite eddies were calculated at each depth by
averaging the geostrophic velocity along the eddy‐core
edge, and we defined their translation speed as being 4.3 cm
s−1. This value corresponds to the mean propagation speed
of long‐lived eddies in the study region [Chaigneau et al.,
2009]. Figure 4 shows the vertical profiles of the mean
swirl velocity at the eddy‐core edge for both the CE and AE.
The vertical extent of the trapped fluid in the composite CE
is of 240 m and of 530 m in the AE. Again, these values are
probably variable eddy to eddy and might be proportional to
the amplitude of the individual eddies. Hereafter, the term
“eddy” will refer to the eddy core extending from the surface to these “trapping depths.”
[22] Table 2 shows that on average along the vertical, the
composite AE is slightly smaller and of a more circular
shape than the composite CE. The mean equivalent radius
is of ∼60 km for both composite eddies corresponding to
typical areas of ∼10000 km2. These spatial scales correspond to the radii of the eddy cores which extend from the
eddy center to the closed DH′ contour associated with a
maximum rotational speed. The size of the eddy cores is
about half the total eddy size of 122 km determined previously from altimetry (see section 3.1). Table 2 also suggests
that both composites AE and CE are not perfectly circular
but slightly elongated along a southwestward/northeastward
direction. The mean ratios of the major to the minor axes of
the ellipses are of 1.3 for CE and 1.2 for AE, in agreement
with the values of 1.5–1.7 obtained at the surface from
altimetry [Chaigneau et al., 2008].
4. Eddy Vertical Structure and Associated
Transports
4.1. Mean Vertical Thermohaline Anomalies
of Composite Eddies
[23] In order to illustrate the differences between CE and
AE in terms of thermohaline structure, Figure 5 shows
vertical sections of the composite eddies along the zonal
direction at Dy = 0. CE (AE, respectively) shows a maximum
temperature anomaly of about −1°C (+1°C) centered at 150 m
(400 m). Although temperature anomalies greater than 0.5°C
extend between 50 m and 350 m for CE (Figure 5a) and
between 200 m and 600 m for AE (Figure 5b), we still observe
7 of 16
C11025
C11025
CHAIGNEAU ET AL.: EDDY VERTICAL STRUCTURE IN THE ESP
Figure 4. Vertical profiles of the swirl velocity averaged over the composite cyclonic (blue) and anticyclonic (red) eddy edges. The lower scale corresponds to the swirl velocity (in cm s−1), whereas the
upper scale corresponds to the nonlinearity parameter (ratio of the swirl velocity to the drift velocity)
using a mean drift velocity of 4.3 cm s−1. The indicated depths correspond to the vertical extent of the
trapped fluid within the composite eddies (240 m for CE and 530 m for AE).
a weak temperature anomaly of ±0.15°C at 1000 m depth.
Similarly, maximum salinity anomalies of ±0.1 are approximately centered at the same vertical levels than their temperature counterparts (Figures 5c and 5d). For the composite
AE, we note, however, that (1) maximum salinity anomalies
are slightly shifted upward compared to temperature
anomalies and (2) the eddy core is surmounted by relatively
fresh and cold water. This latter feature, also noted in model
simulations off Peru [Colas et al., 2011], is probably due to
the doming of the isotherms above the AE cores. The negative
temperature anomalies in the upper layers also suggest that
sea surface temperature alone would not be a good variable
for detection of AEs.
[24] Along the zonal section at Dy = 0, the combination of
′ and S′ leads to meridional geostrophic velocity anomalies
on the order of ±10 cm s−1, with maximum anomalies near
eddy‐core edges (Figures 5e and 5f). These geostrophic
velocities were computed considering the Coriolis parameter
at 20°S which corresponds to the mean latitude of the study
region. Again, there is a notable difference in terms of eddy
kinematics, since maximum velocity anomalies are observed
in the upper 0–200 m layer for CE and between 150 m and
400 m depth for AE (see also Figure 4). At the surface,
geostrophic velocities vary from a few cm s−1 near the
center to 7–10 cm s−1 at the eddy edges (Figures 4 and 5).
These later values are consistent with the eddy swirl speed
of 10–12 cm s−1 estimated from near‐surface drifter trajectories and altimetry measurements [Chaigneau and Pizarro,
2005c].
[25] Figures 6a–6c show the mean vertical temperature,
salinity and density anomalies obtained inside the ellipses
determining the composite eddy edges. These vertical
profiles are also compared with those obtained by averaging
all CTD profiles located either inside CEs, or inside AEs, or
OEs. In the eddy core of the composite eddies, we observed
mean maximum anomalies larger than ±0.7°C in temperature and ±0.06 in salinity. These /S anomalies lead to
maximum density anomalies on the order of +0.10 kg m−3
(−0.08 kg m−3, respectively) for CE (AE). Again, we note
the three main characteristics previously observed along the
zonal vertical section: (1) a core with maximum anomalies
centered at around 150 m for CE and 400 m depth for AE,
(2) a maximum salinity anomaly centered slightly above the
maximum temperature anomaly for AE, and (3) a layer of
relatively fresh and cold water overlaying the AE core. Thus,
Figures 5 and 6 show that on average, the thermohaline
structure of composite CEs and AEs and associated kinematics strongly differs. The core of the composite CE is centered at 150 m depth, a level corresponding to the seasonal
Table 2. Mean Geometrical Properties of the Composite Cyclonic
and Anticyclonic Eddiesa
Vertical extent (m)
Semimajor axis A (km)
Semiminor axis B (km)
Angle 8 (deg)
A/B
Radius (km)
Area (× 103 km2)
Cyclonic Eddy
Anticyclonic Eddy
240
72.3 ± 3.4
53.7 ± 2.3
14.1 ± 5.6
1.3 ± 0.0
62.3 ± 2.8
12.2 ± 1.1
530
63.4 ± 1.6
52.4 ± 1.2
44.7 ± 8.4
1.2 ± 0.0
57.6 ± 1.2
10.4 ± 0.4
a
Numbers indicate the vertical average ±1 standard deviation. A and B
represent the semimajor and semiminor axes of the fitted ellipses,
whereas 8 corresponds to the ellipse orientation relative to the eastward
direction.
8 of 16
C11025
CHAIGNEAU ET AL.: EDDY VERTICAL STRUCTURE IN THE ESP
C11025
Figure 5. Vertical section at Dy = 0 across the composite cyclonic and anticyclonic eddies. (a, b) Potential temperature anomaly (°C), (c, d) salinity anomaly (× 10−2), and (e, f) meridional (cross section) geostrophic speed (cm s−1) relative to 1000 m, indicating the clockwise (anticlockwise, respectively) rotation
for the composite cyclonic (anticyclonic) eddy. Eddy edges are denoted by black lines, whereas horizontal
dashed lines indicate the trapping depths.
thermocline/halocline/pycnocline characterized by a temperature of ∼14°C, a salinity of ∼34.7 and a density of
∼26 kg m−3 (Table 3). In contrast, the core of the composite
AE is centered below the seasonal thermocline at 400 m
depth, a level associated with Equatorial Subsurface Water
with a temperature of 9.5°C, a salinity of ∼34.6 and a
density of ∼26.7 kg m−3. The origin of this important difference between CE and AE morphology will be discussed
in section 5.
[26] Composite anomalies are consistent with those
obtained from all CTD profiles (thin lines in Figures 6a–6c)
but are stronger since they only correspond to the eddy core
located inside the ellipse. Outside the composite eddy core
and far from the eddy center, anomalies are weaker and
hence the average vertical profiles determined directly from
CTD measurements show weaker anomalies. Although
temperature and density anomalies are significantly different
from anomalies obtained outside eddies down to 1000 m
depth (Figures 6a and 6c; see also Figure 5), the fluid which
is more likely trapped and transported by the composite
eddies is limited to the upper 240 m for CE and 530 m for
AE. Table 4 shows the mean anomalies estimated from the
Argo profiles located within the trapped fluid. It confirms
that anomalies are relatively similar but of opposite signs
between CE and AE and are much larger than outside
eddies. Finally, as suggested by Figure 6 and Table 4, CTD
profiles located OEs show on average small positive biases
of 0.05°C in temperature and 0.5 × 10−2 in salinity compared to the CARS climatology. The positive salinity bias is
relatively strong in the upper 250 m and almost 0 in the
9 of 16
C11025
C11025
CHAIGNEAU ET AL.: EDDY VERTICAL STRUCTURE IN THE ESP
Figure 6. Mean vertical profiles of (a) temperature anomaly, (b) salinity anomaly, and (c) density
anomaly inside eddies. Thick lines correspond to the mean profiles obtained inside the composite CE (blue
line) and AE (red line). Thin lines correspond to the mean profiles obtained using all the Argo profiles
located inside the 420 CEs (blue line) and 526 AEs (red line). Black thick lines correspond to the mean
vertical profiles obtained from the 3233 Argo profilers located outside the eddies. Anomalies were computed relatively to the CARS climatology interpolated to the positions and times of the Argo floats.
deepest layers (black curve in Figure 6b). In contrast, the
positive temperature bias is observed from the surface to
1000 m depth (black curve in Figure 6a).
4.2. Meridional Variation of the Eddy Vertical
Structure
[27] In order to estimate the robustness of the previous
results and to study the potential meridional variation in the
eddy vertical structure, the same composite analysis was
repeated for the northern (10°S–20°S) and southern (20°S–
30°S) regions of the study domain (see Figure 1a). The
northern region includes 210–220 CTD profiles in both CEs
and AEs whereas the southern region contains a similar
number of profiles in CEs but 302 in AEs. In both regions,
the radius of composite eddies are similar and of ∼58 km
for AEs and ∼61 km for CEs.
[28] Figure 7 shows the eddy vertical structure along the
zonal direction at Dy = 0 for both regions. Vertical structures of eddies in the two regions are qualitatively similar to
the mean structure obtained for the whole PCCS. However,
Figure 7 shows that temperature and salinity anomalies
associated with CEs are larger in the northern region. In
contrast, the composite AE of the northern region shows
smaller ′ and S′ than in the southern region. The swirl
velocity of the composite eddies are maximum in the surface
for CEs (Figures 7i and 7k) and in subsurface for AEs
(Figures 7j and 7l), but eddies rotate more rapidly in the
northern region.
[29] Another clear difference between the northern and
southern regions is the vertical position of the eddy cores.
For both AEs and CEs the maximum ′ and S′ are shallower
in the northern region. For instance, maximum temperature
anomalies are found at 100 m depth in the composite CE of
the northern region and at 150 m depth in the southern
region (Figures 7a and 7c). Similarly, maximum ′ are found
at 350 m depth in the northern composite AE and at ∼500 m
depth in the southern AE (Figures 7b and 7d). Again, the
salinity anomaly cores in the AEs are found slightly above
the maximum ′ levels. Independently of the considered
subregion and as also noted in Figure 5, the relatively warm
and salty AE cores are overlaid by relatively fresh and cold
water. Finally, the vertical extent of the eddy core transporting trapped fluid can be estimated considering a mean
drift velocity of 5.6 cm s−1 between 10°S–20°S and of
3.2 cm s−1 between 20°S–30°S [Chaigneau et al., 2009].
For CEs, this vertical extent is of 200 m in the North and
280 m in the South. For AEs, the trapping depth also
increases southward, varying from 450 m in the northern
region to 570 m in the southern region.
4.3. Volume, Heat, and Salt Transports
[30] In sections 4.1 and 4.2, the combination of altimetry
data and Argo profiles has been shown to be ideally suited to
produce a good approximation of the three‐dimensional
structure of both cyclonic and anticyclonic eddies. We thus
now use this three‐dimensional structure to estimate the relative eddy contribution to fluxes of volume, heat, and salt in
the PCCS. In general, annual mean “eddy” transports can be
estimated by dividing the volume, heat or salt anomaly of one
eddy by 1 year and taking into account the number of eddies
per year [e.g., Gordon and Haxby, 1990; van Ballegooyen
et al., 1994; Doglioli et al., 2007]. The underlying assumption is that mesoscale eddies are sufficiently nonlinear so that
Table 3. Mean Seawater Properties (, S, s) and Their Anomalies
(′, S′, s′) Observed in the Core of the Composite Eddies at the
Depths Where Density Anomalies s′ Are Maximuma
Depth (m)
(°C)
S
s (kg m−3)
′ (°C)
S′
s′ (kg m−3)
Cyclonic Eddy
Anticyclonic Eddy
150
13.63 ± 0.21
34.68 ± 0.06
25.99 ± 0.05
−0.72 ± 0.14
0.05 ± 0.01
0.11 ± 0.03
400
9.53 ± 0.12
34.61 ± 0.02
26.73 ± 0.01
0.77 ± 0.10
0.06 ± 0.01
−0.08 ± 0.01
a
Here , S, and s are mean seawater properties and ′, S′, and s′ are
their anomalies. Density anomalies s′ are maximum at 150 m for CE
and at 400 m for AE. These values were obtained directly from the
objectively interpolated gridded fields. Numbers indicate the average
±1 standard deviation.
10 of 16
C11025
C11025
CHAIGNEAU ET AL.: EDDY VERTICAL STRUCTURE IN THE ESP
Table 4. Mean Temperature, Salinity, and Density Anomalies of the Eddy Cores From the Surface Down to the Trapping Depths of
240 m for CEs and 530 m for AEsa
Mean ′ (°C)
Mean S′ (× 10−2)
Mean s′ (× 10−2 kg m3)
Cyclonic Eddies
(0–240 m)
Outside Eddies
(0–240 m)
Anticyclonic Eddies
(0–530 m)
Outside Eddies
(0–530 m)
−0.47 ± 0.25
−3.7 ± 2.4
6.8 ± 3.9
0.07 ± 0.04
0.6 ± 0.5
−1.5 ± 0.6
0.41 ± 0.40
4.5 ± 3.8
−3.3 ± 4.9
0.06 ± 0.03
0.3 ± 0.4
−1.2 ± 0.5
a
Values were obtained considering only the Argo profiles located in the core of the composite eddies (93 inside CEs and 162 inside AEs) and can be
compared to the values obtained outside eddies in the same vertical ranges. Numbers indicate the vertical average ±1 standard deviation.
the anomalies are trapped inside their core, down to the
“trapping depth” defined using the Flierl [1981] criterion.
[31] On average over the PCCS, the volume of trapped
fluid transported by an eddy is estimated to be 3.4 × 1012 m3
for a CE and 5.1 × 1012 m3 for an AE (Table 5). Spread over
a 1 year period, the westward volume flux associated with
the trapped fluid of a single eddy is of 0.1–0.2 Sv (Table 5).
Table 5 also indicates that these volume transports are ∼40%
weaker in the northern region than in the south, owing to a
shallower trapping depth and thus a reduced volume. To
estimate how much warm and salty (cold and fresh,
respectively) water is being transported by AEs (CEs), we
calculated the available heat and salt content anomalies
(AHA and ASA) per meter on the vertical:
Z
AHA ¼
Cp ′dA;
ð1Þ
Z
ASA ¼ 0:001
S′dA;
ð2Þ
where r is the density (in kg m−3), Cp is the specific heat
capacity (4000 J kg−1 K−1), and ′, S′ are integrated over the
area (A) of the composite eddy delimited by the ellipse. The
Figure 7. Same as Figure 5 but for the northern and southern regions delimited by the black dashed line
in Figure 1a. The trapping depths related to the vertical extent of the trapped fluid were estimated using a
drift velocity of 5.6 and 3.2 cm s−1 in the northern and southern regions, respectively.
11 of 16
C11025
C11025
CHAIGNEAU ET AL.: EDDY VERTICAL STRUCTURE IN THE ESP
Figure 8. Mean available (a) heat anomaly and (b) salt anomaly, inside composite CE (blue lines) and
AE (red lines). Horizontal dashed lines indicate the trapping depths.
factor 0.001 converts salinity to salinity fraction (kg of salt
per kg of seawater). Figures 8a and 8b show the vertical
profiles of AHA and ASA within composite CE and AE. In
the core of the composite eddies, maximum AHA are larger
than ±0.3 × 1018 J m−1 whereas ASA are larger than ±0.8 ×
1010 kg m−1. We also note a secondary local maximum of
AHA and ASA at ∼300 m depth in CE. This local maximum, which is related to a larger eddy area in the deeper
layers (see Figure 5), is however probably not trapped and
thus not transported by the eddies.
[32] Integrated from the surface to the trapping depth
(240 m for CE and 530 m for AE), the total available heat
and salt content anomalies transported by a single CE (AE,
respectively) is of −5.5 × 1018 J (+8.7 × 1018 J) and −9.8 ×
1010 kg (+23.8 × 1010 kg) (Table 5). Spread over a 1 year
period, the anomalies of heat and salt transports associated
with the trapped fluid of a single CE (AE, respectively) are
of −1.7 × 1011 W (+2.8 × 1011 W) and −3.1 × 103 kg s−1
(+7.5 × 103 kg s−1), respectively (Table 5). The heat and
salt transports associated with CEs (AEs) are stronger
(weaker) in the northern region than in the southern region.
[33] Applying the eddy tracking algorithm developed by
Chaigneau et al. [2009] on the SLA data set, we determined
that around 310 eddies are formed each year near the coast
between 10°S and 30°S (not shown). Generally, these
eddies are quickly dissipated and thus only influence near‐
coastal regions. However, among these 310 eddies generated each year around 30 CEs and 25 AEs can be tracked for
more than 3 months and propagate at a mean distance of
700 km from the coast. These long‐lived eddies can thus
impact the offshore ocean of the PCCS. On the basis of the
average values shown in Table 5, we estimate that the
annual volume anomaly of trapped fluid transported by
these 55 long‐lived eddies is of ∼7.5 Sv. Similarly, the
annual heat (salt, respectively) transport anomaly associated
with 30 long‐lived CEs would be of around −50 × 1011 W
(−90 × 103 kg s−1) and of 70 × 1011 W (190.0 × 103 kg s−1)
for 25 AEs. Note that since anomalies associated with CEs
and AEs are of opposite signs (the negative anomalies
indicate heat and salt deficiencies with respect to the background water mass), the total contribution of the mesoscale
eddies on the heat and salt transport is probably weak.
However, as clearly observed in Figure 8, heat and salt flux
anomalies are transported in distinct depth layers depending
on the eddy rotation. On average, mesoscale eddies would
contribute to inject relatively warm and salty water in the
subthermocline water (below 200–300 m depth) and to cool
and freshen the upper ocean above 200 m depth.
[34] The volume, heat and salt transports estimated for the
PCCS are relatively weak compared to more energetic
regions of the Southern Hemisphere such as in the Antarctic
Circumpolar Current [Joyce et al., 1981; Peterson et al.,
1982; Morrow et al., 2004; Swart et al., 2008], the Agulhas
Current [Gordon and Haxby, 1990; van Ballegooyen et al.,
1994; van Aken et al., 2003], the Brazil‐Malvinas region
[Gordon, 1989; Lentini et al., 2002; de Souza et al., 2006], or
in the South Atlantic subtropical gyre [McCartney and
Woodgate‐Jones, 1991]. However, our volume flux estimates are in agreement with the fluxes obtained in the similar
California Current System where model simulations have
Table 5. Thermohaline Contents and Associated Transports
Integrated Over the Volume of the Composite Cyclonic and
Anticyclonic Eddiesa
PCCS
(10°S–30°S)
CE
Vertical extent (m)
Volume (× 1012 m3)
AHA (× 1018 J)
ASA (× 1010 kg)
Volume transport (Sv)
Heat transport (× 1011 W)
Salt transport (× 103 kg s−1)
AE
Northern
Region
(10°S–20°S)
CE
AE
Southern
Region
(20°S–30°S)
CE
0–230 0–540 0–200 0–450 0–280
3.1
5.5
2.6
4.9
3.4
−5.5
8.7
−5.9
6.5
−5.3
−9.8 23.8 −14.7 17.4 −7.5
0.10 0.18 0.08 0.15 0.11
−1.7
2.8
−1.9
2.1
−1.7
−3.1
7.5
−4.7
5.5
−2.4
AE
0–570
6.1
10.5
28.3
0.19
3.3
9.0
a
Estimates were computed in the entire PCCS (10°S–30°S) and in the
northern and southern regions delimited by black heavy lines in Figure 1d.
AHA, available heat anomaly; ASA, available salt anomaly.
12 of 16
C11025
CHAIGNEAU ET AL.: EDDY VERTICAL STRUCTURE IN THE ESP
C11025
Figure 9. Mean vertical profiles of (a, c) temperature anomaly and (b, d) salinity anomaly, inside CEs and
AEs as a function of their surface signature in dynamic height anomaly (relative to 1000 m). Anomalies are
projected either on depth (Figures 9a and 9b) or on s levels (Figures 9c and 9d). Grey lines in Figures 9a
and 9b correspond to s levels ranging from 24.8 to 27.2 kg m−3 with a contour interval of 0.2 kg m−3.
Black heavy lines correspond to the s levels delimiting the eddy cores where /S anomalies are maximum
(25.2–26.0 kg m−3 for CE and 26.0–26.8 kg m−3 for AE).
shown that the fluid trapped inside AEs extends from the
surface to 500 m depth and transport a volume flux of 0.25 Sv
[Cornuelle et al., 2000]. Note that our calculations only
represent part of the total eddy fluxes since we have not
included the possible large contribution of short‐lived eddies
to the volume, heat and salt budgets. Our results do indicate,
nonetheless, that mesoscale eddies can influence the zonal
heat and salt budgets of the PCCS.
5. Discussion
[35] In other eastern boundary current systems, and in
particular in the California Current System, it is now widely
accepted that nonlinear eddies are mainly generated by
instabilities of the near‐coastal currents [Batteen, 1997;
Huyer et al., 1998; Chereskin et al., 2000; Batteen et al.,
2003; Marchesiello et al., 2003; Capet et al., 2008; Kurian
et al., 2011]. In the California Current System, eddies generated by the surface California Current are predominantly
cyclonic with a surface core in the upper 150 m, whereas
those shed by the subsurface California Undercurrent are
mainly anticyclones with a subsurface core (400 m) [Simpson
and Lynn, 1990; Huyer et al., 1998; Garfield et al., 1999;
Chereskin et al., 2000; Cornuelle et al., 2000]. In the PCCS,
the instability of near‐coastal currents has been also considered as a key parameter for the generation of mesoscale
eddies that are mainly formed near the coast and propagate
westward [Leth and Middleton, 2004; Capet et al., 2008;
Johnson and McTaggart, 2010; Colas et al., 2011]. In fact,
two recent studies based on regional numerical simulations
have shown that modeled eddies in the California Current
System and the PCCS shared similar vertical characteristics
[Kurian et al., 2011; Colas et al., 2011]. In both these eastern
boundary current systems, modeled CEs were surface
intensified with their cores located in the upper 100–200 m
whereas AEs showed maximum subsurface /S anomalies
centered between 200 and 400 m depth [see Kurian et al.,
2011, Figures 16 and 17; Colas et al., 2011, Figures 14
and 15].
[36] In the present study, the composite analysis of Argo
CTD profiles corroborates the results obtained from regional
simulations, showing that in the PCCS the main core of the
CE is located above 200 m depth, whereas the core of the
AE is centered at 400 m depth (Figure 5). The vertical
position of the CE and AE cores is rather independent of the
eddy intensity, defined here as the surface DH′ value
(Figures 9a and 9b), and the maximum ′/S′ correspond to
density ranges of 25.2–26 kg m−3 for CEs and of 26–26.8 kg
m−3 for AEs (Figures 9c and 9d). The CE density class
(25.2–26 kg m−3) is associated with thermocline water
above the Eastern South Pacific Intermediate Water whereas
the AE density range (26.0–26.8 kg m−3) is related to subthermocline waters including both the Eastern South Pacific
Intermediate Water and the Equatorial Subsurface Water
13 of 16
C11025
CHAIGNEAU ET AL.: EDDY VERTICAL STRUCTURE IN THE ESP
(Figure 1d). Once generated near the coast, one can assume
that the fluid trapped in the westward propagating eddies
transports /S anomalies mainly along isopycnal levels. But
near the coast, the density range associated with CEs (25.2–
26 kg m−3) corresponds to the relatively cold and fresh
upwelled water suggesting that intrathermocline CEs may
arise from the meandering of the Chile and Peru surface
currents separating Cold Coastal Water from the warmer and
saltier offshore water. The destabilization of these currents
can trap Cold Coastal Water from their inshore side, forming
CEs which can inject during their offshore displacements
relatively cold and fresh coastal water inside the thermocline/halocline. In contrast, near the coast, the density range
associated with the core of AEs (26–26.8 kg m−3) corresponds to the subsurface PCU flowing poleward [e.g.,
Johnson and McTaggart, 2010]. Thus, these subthermocline
AEs are more likely to be shed by the PCU [Colas et al.,
2011; Johnson and McTaggart, 2010].
[37] The eddy generation by the surface and subsurface
currents are also supported by three other particular features
observed in this study. First, in AEs maximum salinity
anomalies are found above the core of maximum temperature
anomalies (Figures 5–9). But in the offshore ocean, the
Eastern South Pacific Intermediate Water, corresponding to
the 26.0–26.3 kg m−3 s layer, is fresher and warmer that the
underlying Equatorial Subsurface Water associated with the
26.3–26.8 kg m−3 s layer (Figure 1d). Thus, an AE transporting a rather homogeneous water mass from the PCU
toward the offshore ocean, would produce, far from the coast,
a stronger salinity (temperature, respectively) anomaly in
the upper (lower) layer of density 26.0–26.3 kg m−3 (26.3–
26.8 kg m−3). Second, in agreement with a southward
deepening of the PCU core [Silva and Neshyba, 1979; Colas
et al., 2011], the composite AEs also show a deeper core in
the southern region (Figure 7). Third, the temperature and
salinity anomalies were stronger in the core of CEs (AEs,
respectively) in the northern (southern) region. From the
CARS climatology (not shown), we determined that the
temperature and salinity differences between the Cold
Coastal Water (PCU water, respectively) and the offshore
intrathermocline (subthermocline) are also stronger in the
northern (southern) region. Thus, the observed vertical
structure of the composite eddies confirms that CEs are most
likely formed by instabilities of the surface currents whereas
AEs are shed by the PCU.
6. Summary and Future Works
[38] On the basis of altimetry data and available Argo
profiles, this study has characterized the mean vertical
structure of both cyclonic and anticyclonic eddies in the
Peru‐Chile Current System. Through a composite analysis of
the Argo profiles located inside eddies we have highlighted a
key difference between the thermohaline vertical structures
of CEs and AEs: On average, the core of CEs is located in the
25.2–26.0 kg m−3 s layer corresponding to the thermocline,
whereas the core of AEs is found in the subthermocline
layer characterized by s of 26.0–26.8 kg m−3. Within their
core, these eddies exhibit typical temperature and salinity
anomalies of ±1°C and ±0.1 respectively, associated with a
geostrophic velocity on the order of ±10 cm s−1. Intrathermocline CEs are likely to be formed by instabilities of the
C11025
near‐coastal surface currents and propagate offshoreward,
trapping in their cores recently upwelled cold coastal water.
In contrast and as also suggested by Johnson and McTaggart
[2010], subthermocline AEs seem to originate from the PCU
and to impact both the Eastern South Pacific Intermediate
Water and the Equatorial Subsurface water during their
propagation, with strongest salinity anomalies observed in
the layer of the relatively fresh Eastern South Pacific Intermediate Water and strongest temperature anomalies in the
underlying colder layer related to Equatorial Subsurface
Water. Spread over a 1 year period, the volume flux of near‐
coastal water trapped within each eddy is of 0.1–0.2 Sv. The
heat and salt transport anomalies associated with each of
these typical eddies are about ±2–3 × 1011 W and ±3–8 ×
103 kg s−1 depending on their rotation. Negative anomalies
are rather found in the upper thermocline level whereas
positive anomalies are located in subthermocline levels.
[39] The results may be useful for the validation of high‐
resolution regional models [e.g., Colas et al., 2011] and
could be of interest to the biogeochemical community to
investigate links between ecosystems and mesoscale eddies
in the highly productive Peru‐Chile Current System. In turn,
regional models could help to document the exact mechanisms involved in the formation of both the CEs and AEs in
the near‐coastal region. These simulations could also be used
to more precisely estimate the cross‐shore transports of heat
and salt associated with both the mesoscale eddies and the
large‐scale circulation.
[40] The results also raise a number of additional questions that require further investigation. For example, the core
of AEs is located in the Equatorial Subsurface Water layer
associated with the shallowest and most pronounced oxygen
minimum zone of the world ocean [Helly and Levin, 2004;
Fuenzalida et al., 2009; Stramma et al., 2010]. This hypoxic
layer influences biogeochemical cycling of elements and has
a strong impact not only on the local rich ecosystem but also
on the global climate, being a significant source of active
greenhouse gases (CO2 and N2O) [Paulmier and Ruiz‐Pino,
2009]. One can thus expect that AEs also impact the oxygen‐
minimum zone through the offshore propagation of dissolved oxygen anomalies from the near‐coastal region.
Similarly, CEs could also impact the biogeochemical tracer
budgets through the offshore advection of relatively nutrient‐
rich and oxygen‐poor cold coastal water into the thermocline of the offshore ocean. Thus, the impact of both CEs
and AEs on the nutrients and dissolved oxygen contents
should also be addressed in future research. These studies
should in particular focus on the decay of mesoscale eddies,
which is a key process for the redistribution of trapped
properties toward the surrounding water [e.g., Whitney and
Robert, 2002; Johnson et al., 2005; Swart et al., 2008;
van Sebille et al., 2010; Early et al., 2011].
[41] Acknowledgments. Float data used here were collected and
made freely available by Argo (http://www.argo.net/), a program of the
Global Ocean Observing System, and contributing national programs.
The altimeter products were produced by Ssalto/Duacs and distributed by
AVISO, with support from CNES. The preparation and development of this
work were made possible thanks to stays of A. Chaigneau and M. Le Texier
at IMARPE. Several Argo float deployments have been supported by
GMMC (Mercator‐Coriolis) through the Flotteurs du Pacifique Sud
(FLOPS) project. O.P. was supported by FONDECYT 1090791. Support
from ECOS‐CONICYT is also acknowledged. This work was partially
14 of 16
C11025
CHAIGNEAU ET AL.: EDDY VERTICAL STRUCTURE IN THE ESP
funded by the LMI “Dinamica del Sistema de la Corriente de Humboldt.”
The authors particularly thank F. Colas, D. Chelton, and E. van Sebille for
their constructive comments that considerably helped to improve a previous
version of the manuscript. We are also grateful to R. Samelson and
D. Chelton for interesting discussions on the definition of the “eddy trapping
depth.”
References
Albert, A., V. Echevin, M. Lévy, and O. Aumont (2010), Impact of nearshore wind stress curl on coastal circulation and primary productivity in
the Peru upwelling system, J. Geophys. Res., 115, C12033, doi:10.1029/
2010JC006569.
Batteen, M. (1997), Wind‐forced modeling studies of currents, meanders
and eddies in the California Current System, J. Geophys. Res., 102,
2199–2221.
Batteen, M. L., N. J. Cipriano, and J. T. Monroe (2003), A large‐scale seasonal modeling study of the California Current System, J. Oceanogr., 59,
545–562, doi:10.1023/B:JOCE.0000009585.24051.cc.
Blanco, J. L., A. C. Thomas, M.‐E. Carr, and P. T. Strub (2001), Seasonal
climatology of hydrographic conditions in the upwelling region off
northern Chile, J. Geophys. Res., 106, 11,451–11,467, doi:10.1029/
2000JC000540.
Böhme, L., and U. Send (2005), Objective analyses of hydrographic data
for referencing profiling float salinities in highly variable environments,
Deep Sea Res., Part II, 52, 651–664, doi:10.1016/j.dsr2.2004.12.014.
Capet, X., F. Colas, J. C. McWilliams, P. Penven, and P. Marchesiello
(2008), Eddies in eastern‐boundary subtropical upwelling systems, in
Eddy Resolving Ocean Modeling, Geophys. Monogr. Ser., vol. 177,
edited by M. Hecht and H. Hasumi, pp. 131–148, AGU, Washington, D. C.
Chaigneau, A., and O. Pizarro (2005a), Surface circulation and fronts of the
South Pacific Ocean, east of 120°W, Geophys. Res. Lett., 32, L08605,
doi:10.1029/2004GL022070.
Chaigneau, A., and O. Pizarro (2005b), Mean surface circulation and
mesoscale turbulent flow characteristics in the eastern South Pacific from
satellite tracked drifters, J. Geophys. Res., 110, C05014, doi:10.1029/
2004JC002628.
Chaigneau, A., and O. Pizarro (2005c), Eddy characteristics in the eastern
South Pacific, J. Geophys. Res., 110, C06005, doi:10.1029/
2004JC002815.
Chaigneau, A., A. Gizolme, and C. Grados (2008), Mesoscale eddies off
Peru in altimeter records: Identification algorithms and eddy spatio‐
temporal patterns, Prog. Oceanogr., 79, 106–119, doi:10.1016/j.pocean.
2008.10.013.
Chaigneau, A., G. Eldin, and B. Dewitte (2009), Eddy activity in the four
major upwelling systems from satellite altimetry (1992–2007), Prog.
Oceanogr., 83, 117–123, doi:10.1016/j.pocean.2009.07.012.
Chavez, F., A. Bertrand, R. Guevara‐Carrasco, P. Soler, and J. Csirke
(2008), The northern Humboldt Current System: Brief history, present
status and a view towards the future, Prog. Oceanogr., 79, 95–105,
doi:10.1016/j.pocean.2008.10.012.
Chelton, D. B., and M. G. Schlax (2003), The accuracies of smoothed sea
surface height fields constructed from tandem satellite altimeter datasets,
J. Atmos. Oceanic Technol., 20, 1276–1302, doi:10.1175/1520-0426(2003)
020<1276:TAOSSS>2.0.CO;2.
Chelton, D. B., M. G. Schlax, R. M. Samelson, and R. A. de Szoeke
(2007), Global observations of large oceanic eddies, Geophys. Res. Lett.,
34, L15606, doi:10.1029/2007GL030812.
Chelton, D. B., M. G. Schlax, and R. M. Samelson (2011), Global observations of nonlinear mesoscale eddies, Prog. Oceanogr., 91, 167–216,
doi:10.1016/j.pocean.2011.01.002.
Chereskin, T. K., M. Y. Morris, P. P. Niiler, P. M. Kosro, R. L. Smith,
S. R. Ramp, C. A. Collins, and D. L. Musgrave (2000), Spatial and
temporal characteristics of the mesoscale circulation of the California
current from eddy resolving moored and ship‐board measurements,
J. Geophys. Res., 105, 1245–1269, doi:10.1029/1999JC900252.
Colas, F., J. C. McWilliams, X. Capet, and J. Kurian (2011), Heat balance
and eddies in the Peru‐Chile current system, Clim. Dyn., doi:10.1007/
s00382-011-1170-6, in press.
Colbo, K., and R. Weller (2007), The variability and heat budget of the
upper ocean under the Chile‐Peru stratus, J. Mar. Res., 65, 607–637.
Cornuelle, B. D., T. K. Chereskin, P. P. Niiler, and M. Y. Morris
(2000), Observations and modeling of a California undercurrent eddy,
J. Geophys. Res., 105, 1227–1243.
Correa‐Ramirez, M. A., S. Hormazabal, and G. Yuras (2007), Mesoscale
eddies and high chlorophyll concentrations off central Chile (29°–
39°S), Geophys. Res. Lett., 34, L12604, doi:10.1029/2007GL029541.
Davis, R. E. (2005), Intermediate‐depth circulation of the Indian and South
Pacific Oceans measured by autonomous floats, J. Phys. Oceanogr., 35,
683–707, doi:10.1175/JPO2702.1.
C11025
de Souza, R. B., M. M. Mata, C. A. E. Garcia, M. Kampel, E. N. Oliveira,
and J. A. Lorenzzetti (2006), Multi‐sensor satellite and in situ measurements of a warm core ocean eddy south of the Brazil–Malvinas Confluence region, Remote Sens. Environ., 100, 52–66, doi:10.1016/j.rse.2005.
09.018.
Doglioli, A. M., B. Blanke, S. Speich, and G. Lapeyre (2007), Tracking
coherent structures in a regional ocean model with wavelet analysis:
Application to Cape Basin eddies, J. Geophys. Res., 112, C05043,
doi:10.1029/2006JC003952.
Ducet, N., P. Y. Le Traon, and G. Reverdin (2000), Global high‐resolution
mapping of ocean circulation from TOPEX/Poseidon and ERS‐1 and ‐2,
J. Geophys. Res., 105, 19,477–19,498, doi:10.1029/2000JC900063.
Dunn, J. R., and K. R. Ridgway (2002), Mapping ocean properties in
regions of complex topography, Deep Sea Res., Part I, 49, 591–604,
doi:10.1016/S0967-0637(01)00069-3.
Early, J. J., R. M. Samelson, and D. B. Chelton (2011), The evolution and
propagation of quasigeostrophic ocean eddies, J. Phys. Oceanogr., 41,
1535–1555, doi:10.1175/2011JPO4601.1.
Fiedler, P. F., and L. D. Talley (2006), Hydrography of the eastern tropical
Pacific: A review, Prog. Oceanogr., 69, 143–180, doi:10.1016/j.pocean.
2006.03.008.
Flierl, G. R. (1981), Particle motions in large‐amplitude wave fields,
Geophys. Astrophys. Fluid Dyn., 18, 39–74, doi:10.1080/
03091928108208773.
Fuenzalida, R., W. Schneider, J. Garcés‐Vargas, L. Bravo, and C. Lange
(2009), Vertical and horizontal extension of the oxygen minimum
zone in the eastern South Pacific Ocean, Deep Sea Res., Part II, 56,
992–1003, doi:10.1016/j.dsr2.2008.11.001.
Garfield, N., C. A. Collins, R. G. Paquette, and E. Carter (1999), Lagrangian
exploration of the California Undercurrent, 1992–95, J. Phys. Oceanogr.,
29, 560–583, doi:10.1175/1520-0485(1999)029<0560:LEOTCU>2.0.
CO;2.
Gordon, A. L. (1989), Brazil‐Malvinas Confluence: 1984, Deep Sea Res.,
Part A, 36, 359–384, doi:10.1016/0198-0149(89)90042-3.
Gordon, A. L., and W. F. Haxby (1990), Agulhas eddies invade the South
Atlantic: Evidence from Geosat altimeter and shipboard conductivity‐
temperature‐depth survey, J. Geophys. Res., 95, 3117–3125, doi:10.1029/
JC095iC03p03117.
Helly, J. J., and L. A. Levin (2004), Global distribution of naturally occurring marine hypoxia on continental margins, Deep Sea Res., Part I, 51,
1159–1168, doi:10.1016/j.dsr.2004.03.009.
Huyer, A., R. L. Smith, and T. Paluszkiewicz (1987), Coastal upwelling
off Peru during normal and El Niño times, J. Geophys. Res., 92,
14,297–14,307, doi:10.1029/JC092iC13p14297.
Huyer, A., J. A. Barth, P. M. Kosro, R. K. Shearman, and R. L. Smith
(1998), Upper‐ocean water mass characteristics of the California Current,
summer 1993, Deep Sea Res., Part II, 45, 1411–1442, doi:10.1016/
S0967-0645(98)80002-7.
Isern‐Fontanet, J., E. García‐Ladona, and J. Font (2003), Identification of
marine eddies from altimetry, J. Atmos. Oceanic Technol., 20, 772–778,
doi:10.1175/1520-0426(2003)20<772:IOMEFA>2.0.CO;2.
Johnson, G. C., and K. E. McTaggart (2010), Equatorial Pacific 13°C water
eddies in the eastern subtropical South Pacific Ocean, J. Phys. Oceanogr., 40, 226–236, doi:10.1175/2009JPO4287.1.
Johnson, W. K., L. A. Miller, N. E. Sutherland, and C. S. Wong (2005),
Iron transport by mesoscale Haida eddies in the Gulf of Alaska, Deep
Sea Res., Part II, 52, 933–953, doi:10.1016/j.dsr2.2004.08.017.
Joyce, T. M., S. L. Patterson, and R. C. Millard (1981), Anatomy of a cyclonic ring in the Drake Passage, Deep Sea Res., Part A, 28, 1265–1287,
doi:10.1016/0198-0149(81)90034-0.
Katsumata, K., and H. Yoshinari (2010), Uncertainties in global mapping
of Argo drift data at the parking level, J. Oceanogr., 66, 553–569,
doi:10.1007/s10872-010-0046-4.
Kurian, J., F. Colas, X. Capet, J. C. McWilliams, and D. B. Chelton (2011),
Eddy properties in the California Current System, J. Geophys. Res., 116,
C08027, doi:10.1029/2010JC006895.
Lebedev, K. V., H. Yoshinari, N. A. Maximenko, and P. W. Hacker (2007),
YoMaHa’07: Velocity data assessed from trajectories of Argo floats at
parking level and at the sea surface, Tech. Note 4(2), 16 pp., Int. Pac.
Res. Cent., Honolulu.
Lentini, C. A. D., D. B. Olson, and G. Podesta (2002), Statistics of Brazil
Current rings observed from AVHRR: 1993 to 1998, Geophys. Res. Lett.,
29(16), 1811, doi:10.1029/2002GL015221.
Leth, O., and J. Middleton (2004), A mechanism for enhanced upwelling
off central Chile: Eddy advection, J. Geophys. Res., 109, C12020,
doi:10.1029/2003JC002129.
Le Traon, P. Y., and G. Dibarboure (1999), Mesoscale mapping capabilities from multiple altimeter missions, J. Atmos. Oceanic Technol., 16,
15 of 16
C11025
CHAIGNEAU ET AL.: EDDY VERTICAL STRUCTURE IN THE ESP
1208–1223, doi:10.1175/1520-0426(1999)016<1208:MMCOMS>2.0.
CO;2.
Le Traon, P. Y., and F. Ogor (1998), ERS‐1/2 orbit improvement using T/P:
The 2 cm challenge, J. Geophys. Res., 103, 8045–8057, doi:10.1029/
97JC01917.
Lukas, R. (1986), The termination of the Equatorial Undercurrent in the
eastern Pacific, Prog. Oceanogr., 16, 63–90, doi:10.1016/0079-6611
(86)90007-8.
Marchesiello, P., J. C. McWilliams, and A. Shchepetkin (2003), Equilibrium
structure and dynamics of the California Current System, J. Phys.
Oceanogr., 33, 753–783, doi:10.1175/1520-0485(2003)33<753:ESADOT>
2.0.CO;2.
McCartney, M. S., and M. E. Woodgate‐Jones (1991), A deep‐reaching
anticyclonic eddy in the subtropical gyre of the eastern South Atlantic,
Deep Sea Res., 38, suppl. 1, S411–S443.
Montes, I., F. Colas, X. Capet, and W. Schneider (2010), On the pathways
of the equatorial subsurface currents in the eastern equatorial Pacific and
their contributions to the Peru‐Chile Undercurrent, J. Geophys. Res., 115,
C09003, doi:10.1029/2009JC005710.
Morrow, R. A., J.‐R. Donguy, A. Chaigneau, and S. Rintoul (2004), Cold‐
core anomalies at the subantarctic front, south of Tasmania, Deep Sea
Res., Part I, 51, 1417–1440.
Owens, W. B., and A. P. S. Wong (2009), An improved calibration method
for the drift of the conductivity sensor on autonomous CTD profiling floats
by ‐S climatology, Deep Sea Res., Part I, 56, 450–457, doi:10.1016/j.
dsr.2008.09.008.
Pascual, A., Y. Faugère, G. Larnicol, and P.‐Y. Le Traon (2006), Improved
description of the ocean mesoscale variability by combining four satellite
altimeters, Geophys. Res. Lett., 33, L02611, doi:10.1029/2005GL024633.
Paulmier, A., and D. Ruiz‐Pino (2009), Oxygen minimum zones (OMZs)
in the modern ocean, Prog. Oceanogr., 80, 113–128, doi:10.1016/j.
pocean.2008.08.001.
Peterson, R. G., W. D. Nowlin Jr., and T. Whitworth III (1982), Generation
and evolution of a cyclonic ring at Drake Passage in early 1979, J. Phys.
Oceanogr., 12, 712–719, doi:10.1175/1520-0485(1982)012<0712:
GAEOAC>2.0.CO;2.
Ridgway, K. R., J. R. Dunn, and J. L. Wilkin (2002), Ocean interpolation
by four‐dimensional weighted least squares: Application to the waters
around Australasia, J. Atmos. Oceanic Technol., 19, 1357–1375,
doi:10.1175/1520-0426(2002)019<1357:OIBFDW>2.0.CO;2.
Scherrer, B. (2007), Biostatistique, edited by G. Morin, Ed. de la Chenelière,
Montreal, Que., Canada.
Schneider, W., R. Fuenzalida, E. Rodríguez‐Rubio, J. Garcés‐Vargas, and
L. Bravo (2003), Characteristics and formation of eastern South Pacific
Intermediate Water, Geophys. Res. Lett., 30(11), 1581, doi:10.1029/
2003GL017086.
Silva, N., and S. Neshyba (1979), On the southernmost extension of the
Peru‐Chile Undercurrent, Deep Sea Res., Part A, 26, 1387–1393,
doi:10.1016/0198-0149(79)90006-2.
Simpson, J. J., and R. J. Lynn (1990), A mesoscale eddy dipole in the
offshore California Current, J. Geophys. Res., 95, 13,009–13,022,
doi:10.1029/JC095iC08p13009.
Souza, J. M. A., C. de Boyer Montégut, and P. Y. Le Traon (2011), Comparison between three implementations of automatic identification algorithms for the quantification and characterization of mesoscale eddies in
the South Atlantic Ocean, Ocean Sci. Discuss., 8, 483–531, doi:10.5194/
osd-8-483-2011.
C11025
Stramma, L., R. Peterson, and M. Tomczak (1995), The South Pacific
Current, J. Phys. Oceanogr., 25, 77–91, doi:10.1175/1520-0485(1995)
025<0077:TSPC>2.0.CO;2.
Stramma, L., G. C. Johnson, E. Firing, and S. Schmidtko (2010), Eastern
Pacific oxygen minimum zones: Supply paths and multidecadal changes,
J. Geophys. Res., 115, C09011, doi:10.1029/2009JC005976.
Strub, P. T., J. M. Mesias, V. Montecino, J. Rutlant, and S. Salinas (1998),
Coastal ocean circulation off western South America, in The Sea, vol. 11,
edited by A. R. Robinson and K. H. Brink, pp. 273–313, John Wiley,
Hoboken, N. J.
Swart, N. C., I. J. Ansorge, and J. R. E. Lutjeharms (2008), Detailed characterization of a cold Antarctic eddy, J. Geophys. Res., 113, C01009,
doi:10.1029/2007JC004190.
Tomczak, M., and J. Godfrey (1994), Regional Oceanography: An Introduction, 1st ed., 422 pp., Pergamon, New York.
Tsuchiya, M. (1985), The subthermocline phosphate distribution and circulation in the far eastern equatorial Pacific Ocean, Deep‐Sea Res. Part A,
32, 299–313, doi:10.1016/0198-0149(85)90081-0.
van Aken, H., A. K. van Veldhoven, C. Veth, W. P. M. de Ruijter, P. J. van
Leeuwen, S. S. Drijfhout, C. P. Whittle, and M. Rouault (2003), Observations of a young Agulhas ring, Astrid, during MARE in March 2000,
Deep Sea Res., Part II, 50, 167–195, doi:10.1016/S0967-0645(02)
00383-1.
van Ballegooyen, R. C., M. L. Grundlingh, and J. R. E. Lutjeharms
(1994), Eddy fluxes of heat and salt from the southwest Indian Ocean
into the southeast Atlantic Ocean: A case study, J. Geophys. Res., 99,
14,053–14,070.
van Sebille, E., P. J. van Leeuwen, A. Biastoch, and W. P. M. de Ruijter
(2010), On the fast decay of Agulhas rings, J. Geophys. Res., 115,
C03010, doi:10.1029/2009JC005585.
Whitney, F., and M. Robert (2002), Structure of Haida eddies and their
transport of nutrient from coastal margins into the NE Pacific Ocean,
J. Oceanogr., 58, 715–723, doi:10.1023/A:1022850508403.
Wong, A. P. S., and G. C. Johnson (2003), South Pacific Eastern Subtropical Mode Water, J. Phys. Oceanogr., 33, 1493–1509, doi:10.1175/15200485(2003)033<1493:SPESMW>2.0.CO;2.
Wong, A. P. S., G. C. Johnson, and W. B. Owens (2003), Delayed‐mode
calibration of autonomous CTD profiling float salinity data by ‐S climatology, J. Atmos. Oceanic Technol., 20, 308–318, doi:10.1175/15200426(2003)020<0308:DMCOAC>2.0.CO;2.
Wyrtki, K. (1963), The horizontal and vertical field of motion in the Peru
Current, Bull. Scripps Inst. Oceanogr., 8, 313–346.
Wyrtki, K. (1967), Circulation and water masses in the eastern equatorial
Pacific Ocean, Int. J. Oceanol. Limnol., 1(2), 117–147.
A. Chaigneau and G. Eldin, Laboratoire d’Études en Géophysique et
Océanographie Spatiale, UMR 5566, IRD, CNRS, CNES, UPS, 14 Av.
Edouard Belin, Toulouse F‐31400, France. ([email protected])
C. Grados, Instituto del Mar de Perú, Gamarra and Gral Valle Avenue,
Lima, Callao Callao 22, Peru.
M. Le Texier, Département de Formation en Hydraulique et Mécanique
des Fluides, ENSEEIHT, 2 Rue Charles Camichel, Toulouse F‐31071
CEDEX 07, France.
O. Pizarro, Department of Geophysics, UDEC, Avda. Esteban Iturra s/n ‐
Barrio Universitario, Casilla 160‐C, Concepcion, Chile.
16 of 16