Single super-vortex as a proxy for ocean surface flow fields
- 1Department of Physics of Complex Systems, Eötvös Loránd University, Pázmány Péter s. 1/A, 1117 Budapest, Hungary
- 2Max Planck Institute for the Physics of Complex Systems, Nöthnitzer Str. 38, 01187 Dresden, Germany
- 3von Kármán Laboratory for Environmental Flows, Eötvös Loránd University, Pázmány Péter s. 1/A, 1117 Budapest, Hungary
- 4MTA-ELTE Theoretical Physics Research Group, Pázmány Péter s. 1/A, 1117 Budapest, Hungary
- 5Complexity Sciences Center, 9225 Collins Avenue Suite 1208, Surfside, FL 33154, USA
- 6Instituto de Altos Estudos da Paraíba, Rua Silvino Lopes 419-2502, 58039-190 João Pessoa, Brazil
Correspondence: Imre M. Jánosi (firstname.lastname@example.org)
Empirical flow field data evaluation in a well-studied ocean region along the US west coast revealed a surprisingly strong relationship between the surface integrals of kinetic energy and enstrophy (squared vorticity). This relationship defines a single isolated Gaussian super-vortex, whose fitted size parameter is related to the mean eddy size, and the square of the fitted height parameter is proportional to the sum of the square of all individual eddy amplitudes obtained by standard vortex census. This finding allows very effective coarse-grained eddy statistics with minimal computational efforts. As an illustrative example, the westward drift velocity of eddies is determined from a simple cross-correlation analysis of kinetic energy integrals.
Mesoscale eddies (MEs) are energetic, swirling, time-dependent circulatory flows on a characteristic scale of around 100 km (see Fig. 1), which are observed almost everywhere in satellite altimetry data of global sea surface height (Chelton et al., 2007, 2011). The total volume transport by drifting eddies is comparable in magnitude to that of the large-scale wind-driven and thermohaline circulations (Zhang et al., 2014); therefore, MEs play a crucial role in global material and heat transport and mixing of oceans. In spite of their importance, it is far from trivial to identify and characterize MEs from remote sensing data.
The vast majority of the ME studies is based on some automatic algorithm that identifies and tracks the eddies from gridded maps of sea level anomaly (SLA). Various Eulerian methods were developed and deployed in practice, such as detecting closed contours of SLA (Chelton et al., 2011; Mason et al., 2014; Li et al., 2016; Schütte et al., 2016; Pessini et al., 2018), evaluating the geometry of the velocity vectors (Nencioli et al., 2010), determining contours of the Okubo–Weiss parameter (Chelton et al., 2007; Kurian et al., 2011; Ubelmann and Fu, 2011; Schütte et al., 2016; Pessini et al., 2018) or using wavelet analysis to identify coherent eddy-like structures (Rubio et al., 2009; Pnyushkov et al., 2018). Critical comparisons show that none of the Eulerian methods are superior to the others (Souza et al., 2011; Escudier et al., 2016). The algorithms based on searching for finite-time Lagrangian coherent structures obey a better theoretical foundation (Haller, 2015; Beron-Vera et al., 2018; Haller et al., 2018); nevertheless, a recent test of 12 different approaches revealed that the various methods often produce very different predictions for coherent structures. In addition, false positives and negatives can be produced too (Hadjighasem et al., 2017). Apart from the difficulties of identifying MEs, Amores et al. (2018) pointed out that the spatial resolution of gridded fields is also a critical limiting factor. It is not surprising that small vortices are detected in large numbers at fine grid sizes. However, it is somewhat unexpected that many large eddies remain unidentified by close contour searching when the velocity field is represented at lower resolutions (Amores et al., 2018).
The original aim of our work was a detailed analysis of kinetic energy budget of the oceanic surface flow field along the US west coast. At the evaluation of integrated kinetic energy and enstrophy (squared vorticity), we found a non-trivial strong temporal correlation between these quantities. Since the dominating flow features are obviously mesoscale eddies (Fig. 1), it is rather straightforward to formulate an explanation related to the description of individual ocean vortices. One of the basic models is the Gaussian geostrophic vortex exhibiting the attractive features of finite total energy and total enstrophy over an infinite domain, and a simple closed relationship between them. We demonstrate here that a single Gaussian super-vortex properly describes the empirical energy/enstrophy ratio over an extended region; furthermore, the height and radius of such super-vortex are strongly related to the mean values over the same area obtained by classical vortex census.
As for the shape of ocean MEs, the common picture is that they are close to Gaussian humps or troughs (Hopfinger and van Heijst, 1993; Chelton et al., 2011). A detailed fitting procedure of about 5 million SLA profiles by Wang et al. (2015) revealed that around 50 % of MEs are indeed Gaussian, another ∼40 % are Gaussian over a sloping background or merger of two close Gaussian eddies, and the rest have a quadratic core resembling Rankine vortices.
An isolated Gaussian circular eddy in geostrophic equilibrium (where the hydrostatic pressure gradient force is balanced by the local Coriolis force) can be characterized by the following radial profiles of height η, tangential velocity v and vertical vorticity ξ (in cylindrical coordinates):
Here, η0 and R are the height and size parameters for the vortex, respectively, g is the gravitational acceleration, and f=2Ωsin (φ) is the local Coriolis parameter at latitude φ with s−1 for the Earth. The label “shielded” in the title of this section refers to the core of such a vortex being surrounded by a ring of opposite vorticity (Tóth and Jánosi, 2015); see Fig. 2c.
The simplest model of planetary-scale dynamics of the ocean is a single layer of homogeneous fluid, described by the two-dimensional (2-D) barotropic Navier–Stokes equations in a co-rotating frame of reference (Bracco et al., 2004). In the absence of dissipative processes, such a model conserves the total kinetic energy and total enstrophy . An appealing property of an isolated Gaussian vortex is that its total kinetic energy and enstrophy are finite over an infinite domain of integration:
Note that the total kinetic energy integral IKE depends only on the height parameter η0, reflecting self-similarity in the velocity field, and that the ratio of the two integrals is simply . This very relationship was utilized in a recent paper by Li et al. (2018), in a different context of studying viscous decay of individual MEs.
Simple visual inspection of a reconstructed geostrophic flow field (Fig. 1) reveals that MEs are indeed the dominating features. The area shown in Fig. 1 is an extremely well-studied region of the California Current System (CCS) both by observations and calibrated high-resolution numerical simulations (Kelly et al., 1998; Strub and James, 2000; Marchesiello et al., 2003; Castelao et al., 2006; Stegmann and Schwing, 2007; Capet et al., 2008a, b; Checkley and Barth, 2009; Matthews and Emery, 2009; Kurian et al., 2011; Molemaker et al., 2015; Yuan and Castelao, 2017). Openly available data compiled by Risien and Strub (2016) comprise a set of fields of sea level anomalies by combining gridded daily altimeter fields with coastal tide gauge data (Saraceno et al., 2008). The geographic area covers 32.0–48.5∘ N (latitude) and 135.0–111.25∘ W (longitude) with a spatial resolution of 0.25∘ × 0.25∘. Daily mean geostrophic velocity fields are produced for the period 1 January 1993–31 December 2014 (8035 d). The primary validation compares geostrophic velocities calculated from the SLA values and velocities measured at four mooring sites in the test region (Risien and Strub, 2016).
Figure 3 illustrates the total enstrophy (squared vorticity) and total kinetic energy (sum of squared velocity components) integrated over the offshore region (see the dashed frame in Fig. 4) for each day of the record. The correlation is strikingly strong, and it is not trivial. When the shore region is included, much larger differences appear, especially when the area of integration is restricted to a narrow band along the shoreline. Figure 5a clearly demonstrates that large correlation coefficients require large enough areas of integration; a value of 0.95 is reached around km2 (∼202 grid cells or 5∘ × 5∘). Nevertheless, the geometry of the area must not be a square. The red and black symbols in Fig. 5a belong to meridional stripes of width of 1 and 2∘ longitudes (smaller areas are stripes eastward from 125.0∘ W where the meridional length is restricted by the land). Their apparent scatter, however, is not random; the correlation coefficients in equal areas of integration (symbols lined up vertically in Fig. 5a) systematically increase with the distance from the shoreline.
By exploiting the strong correlations, the ratio of integrated kinetic energy and integrated enstrophy provides an effective size parameter of a hypothetical Gaussian super-vortex as . Results for the temporal mean values of this quantity are shown in Fig. 5b. Note that the obtained Reff≈50 km scale belongs to the 1σ width of a Gaussian profile given by Eq. (1). A visual contour of the super-vortex on a SLA map would have a radius closer to ∼2.5–3 Reff≈125–150 km (see Figs. 2a and 4).
As for the height parameter of the super-vortex, Eq. (4) is used for an estimate of shown in Fig. 5c. Since it is obtained from the total kinetic energy integrated over various areas A, an appropriate comparison requires a proper normalization. A practical choice somewhat correcting shape differences is the characteristic length scale . The error bars are much larger than the ones in Fig. 5b as a consequence of the marked annual oscillations shown in Fig. 3a. These oscillations are canceled when the ratio of strongly correlated kinetic energy and enstrophy is considered. Similarly to the correlation coefficients in Fig. 5a, the fitted height values of for the meridional stripes (red crosses and black squares) exhibit systematic changes with the distance from the shoreline, as discussed below.
The super-vortex fit makes only sense when the parameters have some relationship with the existing MEs. In order to make such a comparison, we implemented the eddy census procedure of Chelton et al. (2011) based on closed SLA contour searches. The methodology is described in Chelton et al. (2011) and Oliver et al. (2015); here, we emphasize three particular details. (i) The SLA fields in the data bank (Risien and Strub, 2016) exhibit marked annual oscillations; daily spatial mean values are changing between −8.6 and 10.1 cm. Since this range is comparable to the most common amplitude of the individual eddies (see below), we removed daily means before the eddy census. (ii) In order to avoid differences due to various definitions of the eddy amplitude, we adopted the following rule: when the algorithm identified the location (lat, long) of an eddy center, the amplitude value is imported directly from the (corrected) SLA field. (iii) We adopted the “equivalent radius” as scale parameter for an eddy (Chelton et al., 2011); that is, , where Atot is the total sum of grid cell areas identified inside closed SLA contours.
Figure 6 shows the results of eddy census. The histograms are very similar to previous statistics at the same spatial resolution (Stegmann and Schwing, 2007; Chelton et al., 2011; Kurian et al., 2011; Amores et al., 2018). Note that the eddy-scale histograms in Fig. 6a are sensitive to the level spacing parameter Δl of the closed contour search; fine-scale scans identify smaller eddies in a larger number. The oscillations at smaller eddy scales are due to the discretization error; the area of an eddy is composed of an integer number of grid cells. It is clear that the fitted super-vortex parameter Reff fluctuates around the mean values of eddy-scale histograms (black curve in Fig. 6a). We reiterate here that Reff is an 1σ radius of a Gaussian vortex, while S is closer to a “real” visual radius based on a closed contour estimate of zero height anomaly. As for the super-vortex height , Fig. 6b illustrates that it is much larger than the height of individual eddies, as expected, because it is related to the total kinetic energy over the test area (the offshore region, in the particular case). For this reason, we compare the square of eddy amplitudes in what follows.
The significant advantage of using the super-vortex picture emerges when the fits are performed over subregions of the test area. We have shown already results for meridional stripes of widths of 1 and 2∘ in Fig. 5b and c. Figure 7a illustrates local mean Reff values compared with local mean eddy scale 〈S〉 as a function of the mean distance from the shore. Both quantities exhibit very good agreement and a clear tendency of growth when eddies move away from the shore. The error bars reflect temporal fluctuations over the whole period of 8035 d which are much larger for the eddy census data, because their frequency fluctuates strongly from day to day in a given narrow meridional band.
Figure 7b is a comparison of the height parameters of the super-vortex and eddy census. It is reasonable to consider a relationship between and . The former measures the total kinetic energy (see Eq. 4), while the latter is proportional to the sum of kinetic energies of all individual eddies when we assume that all of them are Gaussian vortices. The sum of kinetic energies based on agrees pretty well with direct counting, when the kinetic energy is determined by adding up squared velocity components for each grid cell assigned to an eddy. Figure 7b illustrates that an empirical ratio of around 2 arises in each meridional stripe; that is, the long-term mean value of kinetic energy for individually identified eddies is ∼50 % of the total kinetic energy in the test region. Interestingly, Amores et al. (2018) reported on a partition ratio between 1 and 5 fluctuating strongly in time; however, they note that the total kinetic energy obtained for satellite altimetry accounts only for half of the real value. The tendency of initial growth up to ∼150 km (see Fig. 7b) might be related to the fact that eddies are generated mostly along the shore, and later they slowly decay during the drift in open water.
A well-known characteristic of eddy trajectories is the strong tendency for purely westward propagation (Cushman-Roisin et al., 1990; Chelton et al., 2007, 2011; Kurian et al., 2011). Chelton et al. (2007) found globally that only about 0.25 % of the eddies have mean drift directions that deviated by more than 10∘ from pure zonal; however, Kurian et al. (2011) and Stegmann and Schwing (2007) obtained stronger dispersion in the CCS study area. Together with the traditional eddy-tracking algorithm, we used our approach to evaluate the cross-correlations of total kinetic energy between neighboring meridional bands of width of a single grid cell (0.25∘):
where the time lag τ represents a temporal shift between the two time series by τ days, overbar denotes temporal mean, and σ is the standard deviation in the given band. Indeed, we find clear maxima at nonzero time lags (actual values are between 5 and 8 d) indicating that total kinetic energy and enstrophy are mostly advected in the offshore region; production or loss is almost negligible (considering geostrophic flow). The time lag and distance of neighboring bands permit an easy estimate of westward drift velocities; the results are shown in Fig. 7c. Drift velocity values in the literature are of the same order of magnitude (Stegmann and Schwing, 2007; Kurian et al., 2011; Chelton et al., 2007, 2011), similarly to our test. As for a direct validation, all individual eddy tracks are evaluated which had longer lifetime than 60 d (432 cyclonic and 422 anticyclonic MEs are identified). The cut at 60 d is somewhat arbitrary; however, we think that the detection error from both the limited spatial and temporal resolutions is larger for short living vortices (note that the typical westward traveling distance during 60 d is ∼155–200 km). Drift values estimated from vortex tracking belong to the centers of eddies, and as expected from a stable β drift, no spatial dependence in the zonal direction is revealed. Theoretical considerations suggest that anticyclonic eddies might drift faster than cyclonic ones (Cushman-Roisin et al., 1990); however, we could not detect statistically significant differences between the two subgroups of trajectories.
We proposed a simple description of geostrophic ocean surface flow fields by exploiting the following results. Firstly, a shielded Gaussian vortex has a finite total kinetic energy and finite total enstrophy; the ratio of them is proportional to the square of the radius of the vortex. Secondly, these two quantities determined from empirical velocity data are strongly correlated, and their ratio correlates with the mean eddy size obtained from traditional eddy census. Thirdly, the fitted amplitude parameter is strongly related to the sum of all squared eddy amplitudes. While this description cannot replace traditional eddy census algorithms, it is certainly able to extract coarse-grained eddy statistics in order to follow temporal and regional changes of eddy activity.
IMJ designed the research; IMJ and MV performed the research; GT and JACG contributed new numerical/analytical tools; IMJ and GT analyzed data; and IMJ, MV and JACG wrote the paper.
The authors declare no conflict of interest.
This work was supported by the Max-Planck Institute for the Physics of Complex Systems in the framework of an Advanced Study Group on “Forecasting with Lyapunov Vectors”. Jason A. C. Gallas was supported by CNPq, Brazil.
This research has been supported by the Hungarian National Research, Development and Innovation Office (grant nos. FK-125024 and K-125171).
The article processing charges for this open-access
publication were covered by the Max Planck Society.
This paper was edited by Eric J. M. Delhez and reviewed by two anonymous referees.
Amores, A., Jordà, G., Arsouze, T., and Le Sommer, J.: Up to what extent can we characterize ocean eddies using present-day gridded altimetric products?, J. Geophys. Res.-Oceans, 123, 7220–7236, https://doi.org/10.1029/2018JC014140, 2018. a, b, c, d
Beron-Vera, F. J., Hadjighasem, A., Xia, Q., Olascoaga, M. J., and Haller, G.: Coherent Lagrangian swirls among submesoscale motions, P. Natl. Acad. Sci. USA, 1–6, https://doi.org/10.1073/pnas.1701392115, 2018. a
Bracco, A., von Hardenberg, J., Provenzale, A., Weiss, J. B., and McWilliams, J. C.: Dispersion and mixing in quasigeostrophic turbulence, Phys. Rev. Lett., 92, 084501, https://doi.org/10.1103/PhysRevLett.92.084501, 2004. a
Capet, X., McWilliams, J. C., Molemaker, M. J., and Shchepetkin, A. F.: Mesoscale to submesoscale transition in the California Current System. Part I: Flow structure, eddy flux, and observational tests, J. Phys. Oceanogr., 38, 29–43, https://doi.org/10.1175/2007JPO3671.1, 2008a. a
Capet, X., McWilliams, J. C., Molemaker, M. J., and Shchepetkin, A. F.: Mesoscale to submesoscale transition in the California Current System. Part II: Frontal processes, J. Phys. Oceanogr., 38, 44–64, https://doi.org/10.1175/2007JPO3672.1, 2008b. a
Castelao, R. M., Mavor, T. P., Barth, J. A., and Breaker, L. C.: Sea surface temperature fronts in the California Current System from geostationary satellite observations, J. Geophys. Res.-Oceans, 111, C09026, https://doi.org/10.1029/2006JC003541, 2006. a
Chelton, D. B., Schlax, M. G., Samelson, R. M., and de Szoeke, R. A.: Global observations of large oceanic eddies, Geophys. Res. Lett., 34, L15606, https://doi.org/10.1029/2007GL030812, 2007. a, b, c, d, e
Chelton, D. B., Schlax, M. G., and Samelson, R. M.: Global observations of nonlinear mesoscale eddies, Progr. Oceanogr., 91, 167–216, https://doi.org/10.1016/j.pocean.2011.01.002, 2011. a, b, c, d, e, f, g, h, i, j, k
Cushman-Roisin, B., Tang, B., and Chassignet, E. P.: Westward motion of mesoscale eddies, J. Phys. Oceanogr., 20, 758–768, https://doi.org/10.1175/1520-0485(1990)020<0758:WMOME>2.0.CO;2, 1990. a, b
Escudier, R., Renault, L., Pascual, A., Brasseur, P., Chelton, D., and Beuvier, J.: Eddy properties in the Western Mediterranean Sea from satellite altimetry and a numerical simulation, J. Geophys. Res.-Oceans, 121, 3990–4006, https://doi.org/10.1002/2015JC011371, 2016. a
Hadjighasem, A., Farazmand, M., Blazevski, D., Froyland, G., and Haller, G.: A critical comparison of Lagrangian methods for coherent structure detection, Chaos, 27, 053104, https://doi.org/10.1063/1.4982720, 2017. a
Haller, G.: Lagrangian coherent structures, Annu. Rev. Fluid Mech., 47, 137–-162, https://doi.org/10.1146/annurev-fluid-010313-141322, 2015. a
Hopfinger, E. J. and van Heijst, G. J. F.: Vortices in rotating fluids, Annu. Rev. Fluid Mech., 25, 241–289, https://doi.org/10.1146/annurev.fl.25.010193.001325, 1993. a
Kelly, K. A., Beardsley, R. C., Brink, R. L. K. H., Paduan, J. D., and Chereskin, T. K.: Variability of the near-surface eddy kinetic energy in the California Current based on altimetric, drifter, and moored current data, J. Geophys. Res., 103, 13067–13083, https://doi.org/10.1029/97JC03760, 1998. a
Kurian, J., Colas, F., Capet, X., McWilliams, J. C., and Chelton, D. B.: Eddy properties in the California Current System, J. Geophys. Res., 116, C08027, https://doi.org/10.1029/2010JC006895, 2011. a, b, c, d, e, f
Marchesiello, P., McWilliams, J. C., and Shchepetkin, A.: Equilibrium structure and dynamics of the California Current System, J. Phys. Oceanogr., 33, 753–783, https://doi.org/10.1175/1520-0485(2003)33<753:ESADOT>2.0.CO;2, 2003. a
Mason, E., Pascual, A., and McWilliams, J. C.: A new sea surface height–based code for oceanic mesoscale eddy tracking, J. Atmos. Oceanic Technol., 31, 1181–1188, https://doi.org/10.1175/JTECH-D-14-00019.1, 2014. a
Molemaker, M. J., McWilliams, J. C., and Dewar, W. K.: Submesoscale instability and generation of mesoscale anticyclones near a separation of the California Undercurrent, J. Phys. Oceanogr., 45, 613–629, https://doi.org/10.1175/JPO-D-13-0225.1, 2015. a
Nencioli, F., Dong, C., Dickey, T., Washburn, L., and McWilliams, J. C.: A vector geometry–based eddy detection algorithm and its application to a high-resolution numerical model product and high-frequency radar surface velocities in the Southern California Bight, J. Atmos. Oceanic Technol., 27, 564–579, https://doi.org/10.1175/2009JTECHO725.1, 2010. a
Oliver, E. C. J., O'Kane, T. J., and Holbrook, N. J.: Projected changes to Tasman Sea eddies in a future climate, J. Geophys. Res.-Oceans, 120, 7150–7165, https://doi.org/10.1002/2015JC010993, 2015. a
Pessini, F., Olita, A., Cotroneo, Y., and Perilli, A.: Mesoscale eddies in the Algerian Basin: do they differ as a function of their formation site?, Ocean Science, 14, 669–688, https://doi.org/10.5194/os-14-669-2018, 2018. a, b
Pnyushkov, A., Polyakov, I. V., Padman, L., and Nguyen, A. T.: Structure and dynamics of mesoscale eddies over the Laptev Sea continental slope in the Arctic Ocean, Ocean Sci., 14, 1329–1347, https://doi.org/10.5194/os-14-1329-2018, 2018. a
Risien, C. M. and Strub, P. T.: Blended sea level anomaly fields with enhanced coastal coverage along the U.S. West Coast, Sci. Data, 3, 160013, https://doi.org/10.1038/sdata.2016.13, 2016. a, b, c, d, e, f
Rubio, A., Blanke, B., Speich, S., Grima, N., and Roy, C.: Mesoscale eddy activity in the southern Benguela upwelling system from satellite altimetry and model data, Progr. Oceanogr., 83, 288–295, https://doi.org/10.1016/j.pocean.2009.07.029, 2009. a
Saraceno, M., Strub, P. T., and Kosro, P. M.: Estimates of sea surface height and near-surface alongshore coastal currents from combinations of altimeters and tide gauges, J. Geophys. Res., 113, C11013, https://doi.org/10.1029/2008JC004756, 2008. a
Schütte, F., Brandt, P., and Karstensen, J.: Occurrence and characteristics of mesoscale eddies in the tropical northeastern Atlantic Ocean, Ocean Sci., 12, 663–685, https://doi.org/10.5194/os-12-663-2016, 2016. a, b
Souza, J. M. A. C., de Boyer Montégut, C., and Le Traon, P. Y.: Comparison between three implementations of automatic identification algorithms for the quantification and characterization of mesoscale eddies in the South Atlantic Ocean, Ocean Sci., 7, 317–334, https://doi.org/10.5194/os-7-317-2011, 2011. a
Strub, P. T. and James, C.: Altimeter-derived variability of surface velocities in the California Current System: 2. Seasonal circulation and eddy statistics, Deep-Sea Res. Pt. II, 47, 831–870, https://doi.org/10.1016/S0967-0645(99)00129-0, 2000. a
Wang, Z., Li, Q., Sun, L., Li, S., Yang, Y., and Liu, S.: The most typical shape of oceanic mesoscale eddies from global satellite sea level observations, Front. Earth Sci., 9, 202–208, https://doi.org/10.1007/s11707-014-0478-z, 2015. a
Yuan, Y. and Castelao, R. M.: Eddy-induced sea surface temperature gradients in Eastern Boundary Current Systems, J. Geophys. Res.-Oceans, 122, 4791–4801, https://doi.org/10.1002/2017JC012735, 2017. a