high resolution surface current observations using high-frequency …efml.kaist.ac.kr › ppt ›...

1
ID: OS43A-02 High resolution surface current observations using high-frequency radar network on the U. S. West Coast Sung Yong Kim, Eric J. Terrill, and Bruce D. Cornuelle Scripps Institution of Oceanography, La Jolla, CA 92093-0213 USA [email protected] Abstract The newly-completed U. S. West Coast high-frequency radar (HFR) network provides an unprecedented capability to monitor and under- stand ocean dynamics and phenomenology through hourly surface current measurements at up to 1 km resolution. These novel obser- vations reveal coastal surface ocean variability having multiple tem- poral and spatial scales, poleward propagating alongshore currents, a clear seasonal transition forced by upwelling-favorable winds and their relaxation, and continuity across scales from sub-mesoscale to mesoscale. Observations Alongshore surface currents A time-alongcoast diagram of subinertial alongshore surface currents shows two types of poleward propagating signals with phase speeds of O(10) and O(100 – 300) km day 1 corresponding to typical shelf mode speeds [1]. The spring and fall transitions relating to the on- set and demise of seasonal upwelling-favorable winds are also visible, most prominently of central and northern California [2] (Figs. 1B and 1C). Variance of surface currents The spectra of surface vector currents off the USWC are characterized by variance peaks in the low frequency band [|σ |≤ 0.4 cycles per day (cpd)], in two bands centered at diurnal and semi-diurnal frequencies [tides (K 1 , M 2 , and S 2 ), diurnal wind (S 1 ) and its harmonics], and at in- ertial frequency (f c = 1.06 to 1.49 cpd for 32 N to 48 N) (Fig. 2A) [3]. The inertial variance of surface currents increases to the north, due in part to the response to increasingly persistent and strong winds and to energetic nonlinear interactions due to decreasing baroclinic Rossby deformation radius and increasing intermittency [4]. The vari- ance in the counter-clockwise inertial band implies elliptical motion of near-inertial surface currents rather than pure circular motions. This succinct view of surface current energy enables us to identify possi- ble driving forces for the surface circulation such as local and remote winds, tides, pressure gradients, and baroclinic motions due to coastal boundaries (e.g., capes, headlands, and bays) and changing bottom topography [5]. Scale continuity The continuity of the power spectra of ocean surface currents from sub-mesoscale to mesoscale is investigated by comparing their en- ergy in the wave-number (k =1/L) and frequency (σ =1/T ) domains. The wave-number spectra of HFR surface currents show a consistent and continuous variance distribution across three different resolutions (1, 6, and 20 km) (Fig. 3A). Resolved scales range from O(1000) km to O(1) km, and the spectra decay with k 2 at high wave-number in agreement with sub-mesoscale spectra [6]. Although their spectra can vary with location because of regional variations in driving forces and geostrophic contents, they have a robust k 2 decay. B SD LB SM VT SB SL RP MB SF PR PA SC TN CC CB WB NP LL -40 -30 -20 -10 0 10 20 30 40 C (cm s -1 ) 124° 120° 116° 32° 36° 40° 44° SD LB SM VT SB SL RP MB SF PR PA SC TN CC CB WB NP LL A Longitude (W) Latitude (N) Month (2008) a b c J J A S O SD LB SM VT SB SL RP MB SF PR PA SC TN CC CB WB NP LL Month (2007-2008) a b c J F M A M J J A S O N D J F M A M J J A S O N D J SD LB SM VT SB SL RP MB SF PR PA SC TN CC CB WB NP LL O D Month (2008) a b c S Yeardays (2008) 235 240 245 250 255 260 265 270 Figure 1: A. High-frequency radar (HFR) network on the USWC for the surface current observation (61 HFRs are marked as black dots as of January of 2009). The blue curve indicates the effective spatial cov- erage of surface current map over two years (2007 – 2008) and a red curve denotes the coastline axis. For reference, major coastal regions are denoted by abbreviated two letter names from south to north: San Diego (SD), Long Beach (LB), Santa Monica (SM), San BuenaVentura (VT), Santa Barbara Channel (SB), Port San Luis (SL), Ragged Point (RP), Monterey Bay (MB), San Francisco (SF), Point Reyes (PR), Point Arena (PA), Shelter Cove (SC), Trinidad (TN), Crescent City (CC), Cape Blanco (CB), Winchester Bay (WB), Newport (NP), and Loomis Lake (LL). B and C. Time-alongcoast diagram of subinertial along- shore surface currents (cm s 1 ) for two years (2007 – 2008) and about five months of 2008 (May 20 – October 10), respectively. Red and blue colors indicate the poleward (upcoast) and equatorward (down- coast) currents. The missing observations are presented with white space. Three black lines indicate the reference phase speeds of 10 (a), 100 (b), and 300 (c) km day 1 . D. Time-alongcoast diagram of hourly alongshore surface currents (cm s 1 ) on yeardays 232 – 271 of 2008 (August 20 – September 28). The first day of each month and the midnight of each day are labeled in B – D and D, respectively. Moreover, HFR observations resolve variability at scales smaller than L 100 km, where the noise in satellite altimeter-derived (ALT) cur- rents becomes dominant over oceanic signals [4]. For example, the wave-number spectra of HFR surface currents in a region (e.g., south- ern California) with minimum wind-driven components (Fig. 2B) are comparable to the spectra of along-track ALT geostrophic currents (Envisat and Jason-1) in the northeastern Pacific (30 N – 50 N, 114 W – 133 W). The mismatch at low wave-number (L> 500 km) is because of coverage area of each instrument (e.g., coastal regions and open ocean): HFR surface currents contain alongshore coastal signals with large wavelength (Figs. 1B – 1D). In the frequency domain (Fig. 3B), the energy at the major tidal con- stituents, inertial frequency, and diurnal harmonic frequencies, and their enhanced variance in nearby bands are clearly identified. Their energy decays with a factor between σ 1 and σ 2 at high frequency [7]. In the very-low-frequency band (σ< 2 × 10 2 cpd, i.e., T> 50 days), the coastal surface current spectra include mesoscale energy. Frequency (cpd) A -4 -3 -2 -1 0 1 2 3 4 SD LB SM VT SB SL RP MB SF PR PA SC TN CC CB WB NP LL -0.5 0 0.5 1 1.5 2 2.5 3 3.5 4 Frequency (cpd) α(σ) C K1 f c M2 0 0.5 1 1.5 2 2.5 3 3.5 4 -1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1 0 0.5 1 1.5 2 2.5 0 0.2 0.4 0.6 0.8 1 SD LB SM VT SB SL RP MB SF PR PA SC TN CC CB WB NP LL κ 2 B 3 (log 10 , cm 2 s -2 cpd -1 ) Figure 2: A. The rotary power spectra (log 10 scale, cm 2 s 2 cpd 1 ) of surface vector currents in the entire domain (Fig. 1A) as a function of frequency (cpd) and coastal regions on the USWC, averaged in each frequency bin and alongshore bin. A black curve indicates the inertial frequency. B. Wind skill (κ 2 ) – the fraction of variance explained by coastal surface winds at the NDBC buoys through linear regression. The wind skill between PA and CC is excluded due to lack of concur- rent observations of coastal winds and surface currents on the coast- line axis in Fig. 1A. C. A probability density function of polarization coefficients at each frequency bin. A range of the inertial frequency (f c ) and two tidal frequencies (K 1 and M 2 ) are indicated. See Fig. 1A for the abbreviated name of coastal regions. A 10 -3 10 -2 10 -1 10 0 10 1 10 0 10 1 10 2 10 3 10 4 10 5 Frequency (cpd) Power (cm 2 s -2 cpd -1 ) σ -2 σ -1 K1 M2 S2 SN LF B HFR (1 km) HFR (6 km) HFR (20 km) ALT (Envisat) ALT (Jason-1) SA6 SA1 SA3 S1 95% CI 10 -3 10 -2 10 -1 10 0 10 0 10 1 10 2 10 3 10 4 10 5 10 6 Wavenumber (km -1 ) Power (cm 2 s -2 km) k -2 k -5/3 k -1 1 km 10 km 50 km 500 km 200 km 95% CI HFR (1 km) HFR (6 km) HFR (20 km) ALT (Envisat) ALT (Jason-1) Figure 3: Power spectra of high-frequency radar-derived (HFR; 1, 6, and 20 km resolutions) surface currents and altimeter-derived geostrophic currents (ALT; along-track Envisat and Jason-1) for two years (2007 – 2008) in the wave-number domain (A) [Length scales (L) of 1, 5, 10, 50, 100, 200, 500, and 1000 km are marked] and fre- quency domain (B) [Six seasonal harmonics (SA 1 SA 6 ), spring-neap (SN , 14.765-day), lunar fortnightly (LF , 13.661-day), S 1 , K 1 , M 2 , and S 2 tidal frequencies are marked]. The auxiliary lines are denoted as k 1 , k 5/3 , and k 2 in the wave-number domain and σ 1 and σ 2 in the frequency domain. The 95% confidence interval (CI) of individual spectra is indicated. Concluding remarks High-resolution, large scale coastal surface current measurements en- abled by the USWC HFR network provide an observational basis to examine the coastal surface circulation at scales from sub-mesoscale through mesoscale, including poleward signals near the coast, tide- coherent (barotropic and baroclinic components), and wind-coherent surface currents. Coupled with other ongoing in-situ observational programs and satellite remote sensing missions, surface current ob- servations are crucial for monitoring continuous ocean variability by filling the gaps of in-situ instruments from offshore to nearshore and can provide timely input and fundamental scientific resource to the management of coastal waters. Acknowledgement Authors thank the State of California for the funding of the Coastal Ocean Currents Monitoring Program (COCMP) which established statewide current monitoring infrastructure, the Na- tional Oceanic and Atmospheric Administration (NOAA) for their con- tinued support of the Integrated Ocean Observing System (IOOS), and in particular, the leadership of Dr. Jack Harlan at NOAA in the establishment of a national HFR program, as well as the sponsorship of the National Science Foundation and Office of Naval Research who were early supporters of the development and use of HFR for ocean research. Data from this report originate from the following universi- ties and research organizations on the USWC: Scripps Institution of Oceanography (SIO) at University of California, San Diego, University of Southern California, Marine Science Institute at University of Cal- ifornia, Santa Barbara, California Polytechnic State University, Naval Postgraduate School, Romberg Tiburon Center at San Francisco State University, Humboldt State University, Bodega Marine Laboratory at University of California, Davis, and Oregon State University. The three USWC IOOS regional observing systems – SCCOOS, CeNCOOS, and NANOOS – are acknowledged for their continued support and advocacy of HFR current measurements for monitoring the oceans. The coastal wind data were provided by National Data Buoy Center (NDBC). The altimeter products were produced by SSALTO/DUACS and distributed by AVISO with support from CNES. Special thanks to Coastal Ocean Observing and Development Center (CORDC) staff at SIO. References [1]S. R. Ramp, L. K. Rosenfeld, T. D. Tisch, and M. R. Hicks. Moored observations of the current and temperature structure over the con- tinental shelf off central California 1. A basic description of the vari- ability. J. Geophys. Res., 102(C10):22877–22902, 1997. [2]P. T. Strub, J. S. Allen, A. Huyer, R. L. Smith, and R. C. Beard- sley. Seasonal cycles of currents, temperatures, winds, and sea level over the northeast Pacific continental shelf: 35 N to 48 N. J. Geophys. Res., 92(C2):1507–1526, 1987. [3]J. D. Paduan and L. K. Rosenfeld. Remotely sensed surface cur- rents in Monterey Bay from shore-based HF radar (Coastal Ocean Dynamics Application Radar). J. Geophys. Res., 101(C9):20669– 20686, 1996. [4] D. Stammer. Global characteristics of ocean variablity estimated from regional TOPEX/POSEIDON altimeter measurements. J. Phys. Oceanogr., 27(8):1743–1769, 1997. [5]L. K. Rosenfeld and R. C. Beardsley. Barotropic semidiurnal tidal currents off Northern California during the Coastal Ocean Dynam- ics Experiment (CODE). J. Geophys. Res., 92(C2):1721–1732, 1987. [6] J. C. McWilliams. Submesoscale, coherent vortices in the ocean. J. Geophys. Res., 23(2):162–182, 1985. [7] R. Ferrari and C. Wunsch. Ocean circulation kinetic energy: Reser- voirs, sources, and sinks. Annu. Rev. Fluid Mech., 41:253–282, 2009.

Upload: others

Post on 05-Jul-2020

0 views

Category:

Documents


0 download

TRANSCRIPT

Page 1: High resolution surface current observations using high-frequency …efml.kaist.ac.kr › ppt › wcoija10_poster.pdf · 2019-01-04 · Scripps Institution of Oceanography, La Jolla,

ID: OS43A-02

High resolution surface current observationsusing high-frequency radar network on the U. S. West Coast

Sung Yong Kim,† Eric J. Terrill, and Bruce D. CornuelleScripps Institution of Oceanography, La Jolla, CA 92093-0213 USA

[email protected]

AbstractThe newly-completed U. S. West Coast high-frequency radar (HFR)network provides an unprecedented capability to monitor and under-stand ocean dynamics and phenomenology through hourly surfacecurrent measurements at up to 1 km resolution. These novel obser-vations reveal coastal surface ocean variability having multiple tem-poral and spatial scales, poleward propagating alongshore currents,a clear seasonal transition forced by upwelling-favorable winds andtheir relaxation, and continuity across scales from sub-mesoscale tomesoscale.

ObservationsAlongshore surface currentsA time-alongcoast diagram of subinertial alongshore surface currentsshows two types of poleward propagating signals with phase speedsof O(10) and O(100 – 300) km day−1 corresponding to typical shelfmode speeds [1]. The spring and fall transitions relating to the on-set and demise of seasonal upwelling-favorable winds are also visible,most prominently of central and northern California [2] (Figs. 1B and1C).

Variance of surface currentsThe spectra of surface vector currents off the USWC are characterizedby variance peaks in the low frequency band [|σ| ≤ 0.4 cycles per day(cpd)], in two bands centered at diurnal and semi-diurnal frequencies[tides (K1, M2, and S2), diurnal wind (S1) and its harmonics], and at in-ertial frequency (fc = 1.06 to 1.49 cpd for 32◦N to 48◦N) (Fig. 2A) [3].The inertial variance of surface currents increases to the north, duein part to the response to increasingly persistent and strong windsand to energetic nonlinear interactions due to decreasing baroclinicRossby deformation radius and increasing intermittency [4]. The vari-ance in the counter-clockwise inertial band implies elliptical motion ofnear-inertial surface currents rather than pure circular motions. Thissuccinct view of surface current energy enables us to identify possi-ble driving forces for the surface circulation such as local and remotewinds, tides, pressure gradients, and baroclinic motions due to coastalboundaries (e.g., capes, headlands, and bays) and changing bottomtopography [5].

Scale continuityThe continuity of the power spectra of ocean surface currents fromsub-mesoscale to mesoscale is investigated by comparing their en-ergy in the wave-number (k = 1/L) and frequency (σ = 1/T ) domains.

The wave-number spectra of HFR surface currents show a consistentand continuous variance distribution across three different resolutions(1, 6, and 20 km) (Fig. 3A). Resolved scales range from O(1000)km to O(1) km, and the spectra decay with k−2 at high wave-number

in agreement with sub-mesoscale spectra [6]. Although their spectracan vary with location because of regional variations in driving forcesand geostrophic contents, they have a robust k−2 decay.

B

SDLBSMVTSBSLRPMBSFPRPA

SCTNCC

CBWBNP

LL

−40 −30 −20 −10 0 10 20 30 40

C

(cm s-1)

124° 120° 116°

32°

36°

40°

44°

SD

LBSM

VTSBSLRP

MB

SFPR

PA

SCTNCC

CBWBNP

LL

A

Longitude (W)

Latitude (N)

Month (2008)

a

b

c

J J A S O

SDLBSMVTSB

SLRPMBSFPRPA

SCTNCC

CBWBNP

LL

Month (2007−2008)

a

b

c

J F M A M J J A S O N D J F M A M J J A S O N D J

SDLBSMVTSBSLRPMBSFPR

PA

SCTNCC

CBWBNP

LL

O

D

Month (2008)

a

b

c

S

Yeardays (2008)

235 240 245 250 255 260 265 270

Figure 1: A. High-frequency radar (HFR) network on the USWC for thesurface current observation (61 HFRs are marked as black dots as ofJanuary of 2009). The blue curve indicates the effective spatial cov-erage of surface current map over two years (2007 – 2008) and a redcurve denotes the coastline axis. For reference, major coastal regionsare denoted by abbreviated two letter names from south to north: SanDiego (SD), Long Beach (LB), Santa Monica (SM), San BuenaVentura(VT), Santa Barbara Channel (SB), Port San Luis (SL), Ragged Point(RP), Monterey Bay (MB), San Francisco (SF), Point Reyes (PR), PointArena (PA), Shelter Cove (SC), Trinidad (TN), Crescent City (CC),Cape Blanco (CB), Winchester Bay (WB), Newport (NP), and LoomisLake (LL). B and C. Time-alongcoast diagram of subinertial along-shore surface currents (cm s−1) for two years (2007 – 2008) and aboutfive months of 2008 (May 20 – October 10), respectively. Red andblue colors indicate the poleward (upcoast) and equatorward (down-coast) currents. The missing observations are presented with whitespace. Three black lines indicate the reference phase speeds of 10(a), 100 (b), and 300 (c) km day−1. D. Time-alongcoast diagram ofhourly alongshore surface currents (cm s−1) on yeardays 232 – 271 of2008 (August 20 – September 28). The first day of each month andthe midnight of each day are labeled in B – D and D, respectively.

Moreover, HFR observations resolve variability at scales smaller thanL ≈ 100 km, where the noise in satellite altimeter-derived (ALT) cur-

rents becomes dominant over oceanic signals [4]. For example, thewave-number spectra of HFR surface currents in a region (e.g., south-ern California) with minimum wind-driven components (Fig. 2B) arecomparable to the spectra of along-track ALT geostrophic currents(Envisat and Jason-1) in the northeastern Pacific (30◦N – 50◦N, 114◦W– 133◦W). The mismatch at low wave-number (L > 500 km) is becauseof coverage area of each instrument (e.g., coastal regions and openocean): HFR surface currents contain alongshore coastal signals withlarge wavelength (Figs. 1B – 1D).

In the frequency domain (Fig. 3B), the energy at the major tidal con-stituents, inertial frequency, and diurnal harmonic frequencies, andtheir enhanced variance in nearby bands are clearly identified. Theirenergy decays with a factor between σ−1 and σ−2 at high frequency[7]. In the very-low-frequency band (σ < 2 × 10−2 cpd, i.e., T > 50days), the coastal surface current spectra include mesoscale energy.

Frequency (cpd)

A

−4 −3 −2 −1 0 1 2 3 4

SDLB

SMVTSB

SLRP

MBSFPR

PA

SC

TNCC

CBWB

NP

LL

−0.5 0 0.5 1 1.5 2 2.5 3 3.5 4

Frequency (cpd)

α(σ

)

C

K1 fcM2

0 0.5 1 1.5 2 2.5 3 3.5 4−1

−0.8

−0.6

−0.4

−0.2

0

0.2

0.4

0.6

0.8

1

0 0.5 1 1.5 2 2.5

0 0.2 0.4 0.6 0.8 1

SDLB

SMVTSB

SLRP

MBSFPR

PA

SC

TNCC

CBWB

NP

LL

κ2

B

3

(log10 , cm2 s-2 cpd-1)

Figure 2: A. The rotary power spectra (log10 scale, cm2 s−2 cpd−1) ofsurface vector currents in the entire domain (Fig. 1A) as a function offrequency (cpd) and coastal regions on the USWC, averaged in eachfrequency bin and alongshore bin. A black curve indicates the inertialfrequency. B. Wind skill (κ2) – the fraction of variance explained bycoastal surface winds at the NDBC buoys through linear regression.The wind skill between PA and CC is excluded due to lack of concur-rent observations of coastal winds and surface currents on the coast-line axis in Fig. 1A. C. A probability density function of polarizationcoefficients at each frequency bin. A range of the inertial frequency(fc) and two tidal frequencies (K1 and M2) are indicated. See Fig. 1Afor the abbreviated name of coastal regions.

A

10−3

10−2

10−1

100

101

100

101

102

103

104

105

Frequency (cpd)

Power (cm

2 s−2 cpd−1)

σ−2

σ−1

K1 M2S2

SN LF

B

HFR (1 km)

HFR (6 km)

HFR (20 km)

ALT (Envisat)

ALT (Jason−1)

SA6SA1 SA3 S1

95% CI

10−3

10−2

10−1

100

100

101

102

103

104

105

106

Wavenumber (km−1)

Power (cm

2 s−2 km)

k−2

k−5/3

k−1

1 km10 km50 km500 km 200 km

95% CI

HFR (1 km)

HFR (6 km)

HFR (20 km)

ALT (Envisat)

ALT (Jason−1)

Figure 3: Power spectra of high-frequency radar-derived (HFR; 1,6, and 20 km resolutions) surface currents and altimeter-derivedgeostrophic currents (ALT; along-track Envisat and Jason-1) for twoyears (2007 – 2008) in the wave-number domain (A) [Length scales(L) of 1, 5, 10, 50, 100, 200, 500, and 1000 km are marked] and fre-quency domain (B) [Six seasonal harmonics (SA1 – SA6), spring-neap(SN , 14.765-day), lunar fortnightly (LF , 13.661-day), S1, K1, M2, andS2 tidal frequencies are marked]. The auxiliary lines are denoted ask−1, k−5/3, and k−2 in the wave-number domain and σ−1 and σ−2 inthe frequency domain. The 95% confidence interval (CI) of individualspectra is indicated.

Concluding remarksHigh-resolution, large scale coastal surface current measurements en-abled by the USWC HFR network provide an observational basis toexamine the coastal surface circulation at scales from sub-mesoscalethrough mesoscale, including poleward signals near the coast, tide-coherent (barotropic and baroclinic components), and wind-coherentsurface currents. Coupled with other ongoing in-situ observationalprograms and satellite remote sensing missions, surface current ob-servations are crucial for monitoring continuous ocean variability byfilling the gaps of in-situ instruments from offshore to nearshore and

can provide timely input and fundamental scientific resource to themanagement of coastal waters.

Acknowledgement Authors thank the State of California for thefunding of the Coastal Ocean Currents Monitoring Program (COCMP)which established statewide current monitoring infrastructure, the Na-tional Oceanic and Atmospheric Administration (NOAA) for their con-tinued support of the Integrated Ocean Observing System (IOOS),and in particular, the leadership of Dr. Jack Harlan at NOAA in theestablishment of a national HFR program, as well as the sponsorshipof the National Science Foundation and Office of Naval Research whowere early supporters of the development and use of HFR for oceanresearch. Data from this report originate from the following universi-ties and research organizations on the USWC: Scripps Institution ofOceanography (SIO) at University of California, San Diego, Universityof Southern California, Marine Science Institute at University of Cal-ifornia, Santa Barbara, California Polytechnic State University, NavalPostgraduate School, Romberg Tiburon Center at San Francisco StateUniversity, Humboldt State University, Bodega Marine Laboratory atUniversity of California, Davis, and Oregon State University. The threeUSWC IOOS regional observing systems – SCCOOS, CeNCOOS,and NANOOS – are acknowledged for their continued support andadvocacy of HFR current measurements for monitoring the oceans.The coastal wind data were provided by National Data Buoy Center(NDBC). The altimeter products were produced by SSALTO/DUACSand distributed by AVISO with support from CNES. Special thanks toCoastal Ocean Observing and Development Center (CORDC) staff atSIO.

References[1] S. R. Ramp, L. K. Rosenfeld, T. D. Tisch, and M. R. Hicks. Moored

observations of the current and temperature structure over the con-tinental shelf off central California 1. A basic description of the vari-ability. J. Geophys. Res., 102(C10):22877–22902, 1997.

[2] P. T. Strub, J. S. Allen, A. Huyer, R. L. Smith, and R. C. Beard-sley. Seasonal cycles of currents, temperatures, winds, and sealevel over the northeast Pacific continental shelf: 35◦N to 48◦N. J.Geophys. Res., 92(C2):1507–1526, 1987.

[3] J. D. Paduan and L. K. Rosenfeld. Remotely sensed surface cur-rents in Monterey Bay from shore-based HF radar (Coastal OceanDynamics Application Radar). J. Geophys. Res., 101(C9):20669–20686, 1996.

[4] D. Stammer. Global characteristics of ocean variablity estimatedfrom regional TOPEX/POSEIDON altimeter measurements. J.Phys. Oceanogr., 27(8):1743–1769, 1997.

[5] L. K. Rosenfeld and R. C. Beardsley. Barotropic semidiurnal tidalcurrents off Northern California during the Coastal Ocean Dynam-ics Experiment (CODE). J. Geophys. Res., 92(C2):1721–1732,1987.

[6] J. C. McWilliams. Submesoscale, coherent vortices in the ocean.J. Geophys. Res., 23(2):162–182, 1985.

[7] R. Ferrari and C. Wunsch. Ocean circulation kinetic energy: Reser-voirs, sources, and sinks. Annu. Rev. Fluid Mech., 41:253–282,2009.