Next Article in Journal
A Robust and High-Precision Three-Step Positioning Method for an Airborne SAR Platform
Next Article in Special Issue
The Effect of Spatially Correlated Errors on Sea Surface Height Retrieval from SWOT Altimetry
Previous Article in Journal
Using Local Knowledge and Remote Sensing in the Identification of Informal Settlements in Riyadh City, Saudi Arabia
Previous Article in Special Issue
KaRIn Noise Reduction Using a Convolutional Neural Network for the SWOT Ocean Products
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Implementation of the Optical Flow to Estimate the Propagation of Eddies in the South Atlantic Ocean

by
Denis L. Volkov
1,2,* and
Shahriar Negahdaripour
3
1
Cooperative Institute for Marine and Atmospheric Studies, University of Miami, Miami, FL 33149, USA
2
NOAA Atlantic Oceanographic and Meteorological Laboratory, Miami, FL 33149, USA
3
College of Engineering, University of Miami, Coral Gables, FL 33146, USA
*
Author to whom correspondence should be addressed.
Remote Sens. 2023, 15(15), 3894; https://0-doi-org.brum.beds.ac.uk/10.3390/rs15153894
Submission received: 25 June 2023 / Revised: 27 July 2023 / Accepted: 5 August 2023 / Published: 7 August 2023
(This article belongs to the Special Issue Applications of Satellite Altimetry in Ocean Observation)

Abstract

:
The ocean is filled with mesoscale eddies that account for most of the oceanic kinetic energy. The importance of eddies in transporting properties and energy across the ocean basins has led to numerous efforts to track their motion. Here, we implement a computer vision technique—the optical flow—to map the pathways of mesoscale eddies in the South Atlantic Ocean. The optical flow is applied to the pairs of consecutive sea surface height maps produced from a nearly 30-year-long satellite altimetry record. In contrast to other methods to estimate the eddy propagation velocity, the optical flow can reveal the temporal evolution of eddy motion, which is particularly useful in the regions of strong currents. We present the time-dependent estimates of the speed and direction of eddy propagation in the Eulerian frame of reference. In an excellent agreement with earlier studies, the obtained pattern of eddy propagation reveals the interaction of eddies with the background flow and the bottom topography. We show that in the Antarctic Circumpolar Current, the variability of the eddy propagation velocity is correlated with the variability of the surface geostrophic velocity, demonstrating the robustness of the optical flow to detect the time-variable part of eddy motion.

1. Introduction

The ocean is always in turbulent, highly variable motion. Ocean circulation redistributes mass, heat, salt, nutrients, carbon, and other properties, and plays a crucial role in regulating the global and regional climate, weather, and ecosystems. The advent of satellite observations of sea surface temperature (SST) in the 1970s and sea surface height (SSH) in the 1980s has revolutionized our view on ocean variability. The nearly global coverage of satellite measurements has revealed the distribution and spectral characteristics of ocean variability at various time scales. It has been discovered that ocean variability is dominated by processes with spatial scales ranging from a few tens to a few hundreds of kilometers and a time scale on the order of 100 days [1,2,3]. This type of ocean variability is usually termed as transient mesoscale eddies, including individual coherent eddies (vortices characterized by a clear boundary), meanders, fronts, filaments, and planetary Rossby waves [4]. Transient mesoscale eddies can be found virtually everywhere in the ocean, but the largest and most energetic eddies emerge from the instabilities of the strongest currents, in particular the Antarctic Circumpolar Current (ACC) and western boundary currents, such as the Gulf Stream, Kuroshio, Agulhas Current, etc. Away from strong currents that can significantly influence their motion, eddies tend to propagate westward as Rossby waves. They carry properties horizontally and redistribute them vertically, contribute to the climatically meaningful meridional heat transport, drive regional sea level changes, influence shipping routes and fisheries, etc., e.g., refs. [5,6,7,8].
The importance of mesoscale eddies in ocean dynamics has stimulated numerous efforts to track the movement of individual eddies and to estimate their propagation velocities. The tracking of individual eddies consists of following their trajectories in a Lagrangian manner. Existing techniques include the Okubo–Weiss method [9,10,11,12], algorithms based on the geometry of velocity vectors and streamlines [13,14,15,16,17], a hybrid geometric Okubo–Weiss method [18], a wavelet-based approach [19,20], and an objective method based on geodesic transport theory [21]. Most procedures of eddy detection rely on using either surface geostrophic velocities derived from SSH or closed SSH contours. In altimetry applications, however, the use of geostrophic velocities and moreover the relative vorticity is not optimal, because it requires computation of the first and the second derivatives of SSH, respectively, which amplifies any noise present in satellite measurements and interpolation errors [15]. Different from individual eddy detection and tracking, the regional and global patterns of transient eddy propagation velocities in the Eulerian frame of reference have also been obtained by carrying out a space–time-lagged correlation analysis of SSH fields [22,23].
The objective of this study is to demonstrate the use of the optical flow (OF) technique to derive the time-dependent propagation velocities of transient mesoscale eddies in the Eulerian framework. The OF is a computer vision algorithm that estimates the apparent motion of brightness patterns in a sequence of images, e.g., [24]. Some examples of the earlier implementation of the OF in geosciences include the assessment of the quality of forecasts of meteorological fields [25,26,27], and the estimates of sea-ice motion [28,29] and surface currents from SST images [30,31]. The study reported in this paper is, to our knowledge, the first implementation of the OF method to satellite altimetry SSH data. Like the space–time-lagged correlation analysis of [22,23], the OF method cannot distinguish various forms of eddy variability, e.g., isolated cyclones and anticyclones, planetary waves, meanders, fronts, etc. On the other hand, while the space–time-lagged correlation analysis yields only the time–mean propagation velocities, the OF computes “instantaneous” velocities between the pairs of consecutive satellite images, so given the large number of images, it is well suited for estimating the temporal evolution of motion. In oceanographic applications, this might be particularly useful in regions where the propagation velocities of mesoscale eddies are strongly affected by the variable background ocean circulation.
In this work, we focus on the South Atlantic Ocean and the adjacent regions of the Southern Ocean (Figure 1). This region is characterized by complex dynamics and, therefore, it is particularly representative for testing the OF method. The same region was used by [22] to estimate eddy propagation velocities using a space–time-lagged correlation method, which we also use for validation. The mean dynamic topography contours displayed in Figure 1 indicate the streamlines of the surface flow (with the high SSH to the left of the flow direction in the Southern Hemisphere). The ocean circulation here is dominated by a system of western boundary currents (the Brazil and Malvinas Currents), the eddy-rich eastward-flowing ACC, and the Agulhas Retroflection, which produces prominent eddies that propagate into the South Atlantic. The ACC transports 130–140 Sv (1 Sv = 106 m3 s−1) of water through the Drake Passage [32] and supplies the Malvinas Current that carries 30–40 Sv northward along the Argentine coast [33,34,35]. The southward branch of the South Atlantic subtropical gyre along the Brazil coast is comprised of the Brazil Current with a mean strength of 14 Sv and a standard deviation of about 9 Sv at 34.5°S [36]. The Brazil and the Malvinas currents collide at approximately 39°S forming one of the most energetic regions of the World Ocean—the Brazil–Malvinas Confluence (Figure 1; ref. [37]). A very strong anticyclonic (anticlockwise) barotropic flow transporting approximately 100 Sv exists around a topographic anomaly known as the Zapiola Rise at about 45°S and 45°W [38,39]. This flow is also associated with the anticyclonic propagation of mesoscale eddies [22]. By applying the OF technique to SSH data in the South Atlantic Ocean, we provide a realistic description of eddy propagation velocity validated by several other methods and previous studies.

2. Materials and Methods

2.1. Satellite Altimetry Data

Satellite altimetry has provided accurate, nearly global, sustained observations of SSH since the launch of the TOPEX/Poseidon mission in August 1992. Here, we use the daily SSH maps above the geoid (absolute dynamic topography) from January 1993 to June 2022 processed and distributed by the Copernicus Marine and Environment Monitoring Service (CMEMS; http://marine.copernicus.eu, accessed on 2 May 2023). The SSH above the geoid is the sum of the MDT-CNES-CLS22 mean dynamic topography [40] over the 1993–2012 period and the sea-level anomaly with respect to the mean sea surface. Prior to mapping, the along-track altimetry records are routinely corrected for instrumental noise, orbit determination error, atmospheric refraction, sea state bias, tides, as well as static and dynamic atmospheric-pressure effects [41].
The SSH maps are then produced on a global 1/4° × 1/4° (longitude × latitude) grid by the optimal interpolation of measurements from all the altimeter missions available at a given time. Due to the filtering properties of the optimal interpolation, the SSH maps have an effective spatial resolution ranging from about 400 km at 10°S to 100 km at 60°S and an effective temporal resolution varying between 14–45 days [42]. Although the number of simultaneously flying altimeter satellites has varied with time, it has been reported that the impact of satellite constellation on the effective spatial resolution is modest, with a globally averaged gain of resolution of ∼5% in maps constructed with three altimeters compared with two altimeters [42].
Because the smallest repeat period of altimeter satellites (excluding crossover points) is about 10 days and the effective temporal resolution of SSH maps is even longer, using every single daily SSH map for the OF computation is redundant. Therefore, to reduce the computational cost, we use every 10th SSH map (1076 maps in total). The SSH anomalies are computed with respect to the record-long mean SSH above the geoid (absolute dynamic topography). The global mean SSH time series is subtracted from each grid point to obtain the local dynamic SSH anomalies not related to the global mean sea-level rise. The seasonal cycle, which is mostly characterized by standing oscillations, is estimated by fitting both the annual and semi-annual harmonics in a least-squares sense and subtracted from the SSH fields. Furthermore, to focus on the mesoscale variability only, we high-pass filter the time series at each grid point by removing the signals with periods greater than 1.5 years.

2.2. Optical Flow

Here, we provide a brief description of the OF method as it is applied to calculate the perceived eddy propagation velocities between the pairs of consecutive SSH anomaly maps (images). Let η1(x, y, t) denote the SSH anomaly signal at a grid point or pixel (x, y) on the first map at time t. Suppose that during the time interval ∆t, the water parcel with the signal η1(x, y, t) has moved to a new location (x + ∆x, y + ∆y) on the second map and the SSH anomaly signal has been transformed to η2(x + ∆x, y + ∆y, t + ∆t) (Figure 2a). It is assumed that the SSH change over ∆t, i.e., between consecutive maps, is small, so that
η 2 ( x + Δ x , y + Δ y , t + Δ t ) η 1 ( x , y , t )
This is known as the brightness-constancy assumption in computer vision. Considering the effective resolution of SSH maps, it is reasonable to assume that the displacements ∆x and ∆y are relatively small over ∆t = 10 days. Thus, the left-hand side of (1) can be approximated by a Taylor series expansion up to the first-order terms:
η 2 ( x + Δ x , y + Δ y , t + Δ t ) η 1 ( x , y , t ) + η x Δ x + η y Δ y + η t Δ t
where (ηx, ηy, ηt) is a vector of the zonal, meridional, and temporal derivatives of η. Substituting (2) into (1) and dividing by ∆t yields the OF constraint equation [24]:
η x u + η y v + η t = 0
where u = ∆x/∆t and v = ∆y/∆t are the zonal and meridional propagation velocities (optical flow) of the SSH signal, respectively. While the vector derivative (ηx, ηy, ηt) can be estimated by spatio-temporal finite differencing at each pixel, the two unknowns (u and v) cannot be uniquely determined from (3) alone. Known as the aperture problem, this ambiguity translates to the unique computation of motion along the direction of the gradient (ηx, ηy), but not perpendicular to it. Additional constraints are required to arrive at a unique solution. Among many approaches to address this problem, the Lucas–Kanade method assumes that the OF within a small neighborhood W (e.g., an n × m rectangular patch) of each pixel (x, y) is constant [43]. This yields a linear system of equations for the n × m pixels within W with only two unknowns:
[ η x ( 1 ) η y ( 1 ) η x ( 2 ) η y ( 2 ) η x ( n m ) η y ( n m ) ] [ u v ] = [ η t ( 1 ) η t ( 2 ) η t ( n m ) ]
The least-squares optimization applied to the above over-constrained linear system of the form A n m × 2 u = B n m , where u(u, v) is the optical flow vector, yields the following solution:
u = ( A T A ) 1 A T B
Mathematically, we require a near-unity condition number of A T A for a robust solution. In our application, overcoming the aperture problem, this translates to having sufficient variation of η within W. Generally, SSH maps have sufficiently rich texture within a patch size of 5 × 5 pixels, which we assume in this work (Figure 2a). Any land and out-of-bound pixels falling within W are discarded, reducing the number of equations in (4). For the grid spacing of 1/4° and the patch size of 5 × 5 pixels, the spacing between the estimates of u(u, v) based on non-overlapping W is 1.25°. This is less than the effective spatial resolution of SSH maps everywhere within the domain.
Recall that the derivation of (3) assumes small displacements ∆x and ∆y, on the order of one pixel. For larger displacements, which may happen in dynamically active regions such as the ACC and western boundary currents, the Taylor series approximation (2) and the OF constraint Equation (3) are no longer valid. This is resolved by the recursive implementation of the Lucas–Kanade method applied to multi-resolution Gaussian pyramidal representation of SSH maps [44]. The OF is then computed recursively, starting from the lowest-resolution level of the pyramid, and propagating the resultant flow to and updating it at the next higher-resolution level until the original map (with the highest resolution) is reached (Figure 2b). The number of levels in the pyramid is set to achieve a motion of roughly one-pixel at the lowest resolution. In this work, we use three levels, where the original SSH map resolution is reduced by a factor of two relative to the previous level.
Unlike the space–time-lagged correlation analysis, the OF technique does not require spatial and temporal constraints to isolate the required scales of ocean variability. Instead, it estimates the motion of individual pixels that may belong to SSH patterns of any spatial scale. The scale separation can be performed at the stage of data preparation by applying spatial and temporal filtering, as was done in this work by removing the global mean SSH, the seasonal cycle, and the low-frequency signals with periods greater than 1.5 years. The residual variability is dominated by transient mesoscale eddies. As mentioned above, because u and v are computed between two consecutive images, applying the OF to a series of SSH maps yields the time series of the propagation velocities of SSH patterns.

3. Results

3.1. Time–Mean Eddy Propagation

The intensity of eddy variability in the South Atlantic is demonstrated by the standard deviation of SSH (Figure 1). The most energetic regions are associated with the Agulhas Retroflection, the Agulhas Return Current, and the Brazil–Malvinas Confluence. Elevated eddy energy is also observed along the ACC, especially near the Drake Passage and between 30–40°E. The time–mean speed and direction of the horizontal movement of SSH patterns associated with transient eddy variability estimated with the OF are displayed in Figure 3. Like the space–time correlation analysis, the OF provides a description of the transient eddy propagation velocity in the Eulerian frame of reference. The method reveals the relatively swift propagation of eddies along the Agulhas Current flowing along the southeastern coast of Africa with velocities of 5–7 km/day. Upon reaching the Agulhas Retroflection, the eddy motion becomes weak, and no eddies tend to follow the eastward Agulhas Return Current. Once new eddies are formed by the Agulhas Retroflection (around 20°E, 40°S), they start propagating westward-northwestward with a speed of 4–5 km/day. There is no significant northward eddy propagation along the relatively weak Benguela Current that flows along the southwestern coast of Africa.
Nevertheless, eddies are generated near the eastern boundary and propagate westward as Rossby waves in a vast region north of 35°S. This region represents the northern part of the South Atlantic subtropical gyre. The mesoscale variability towards the equator becomes weak (Figure 1), and westward-propagating signals have smaller amplitudes compared to signals at higher latitudes. The speeds of westward propagation increase from 2–3 km/day at 33–35°S to 5–7 km/day at 10–15°S (Figure 3). There is a remarkable transition between the coherent basin-wide westward propagation north of 33°S and the region south of this latitude and to the east of the Argentine Basin with a very weak westward flow. This meridional contrast is apparently the result of the mean-flow effect on eddy propagation. As illustrated by the contours of the mean dynamic topography (Figure 1), the latitude of 33°S divides the subtropical gyre in two parts: the northern part with the general northwestward direction of the surface flow and the southern part with the eastward direction of surface currents. This means that the mean flow reinforces the westward propagation north of 33°S and compensates for the westward propagation south of this latitude.
Near the western boundary, there is a well-defined band of southwestward velocities of up to 5 km/day between 25°S and 40°S associated with the Brazil Current. Further south, along the Argentine coast, the northward-flowing Malvinas Current is also a pathway of eddy movement, with eddy propagation velocities reaching 7 km/day. The bottom topography plays an important role in constraining the eddy motion, especially in the Southern Ocean where stratification is weak and the vertical extent of eddies can reach the bottom (e.g., [23]). The Argentine Basin represents another good example of eddy propagation steered by the bottom topography (Figure 4). The Southern Ocean eddies enter the basin from the Scotia Sea, mainly through channels in the North Scotia Ridge, and follow the continental shelf break. At about 40°S, eddies carried by the northward Malvinas Current collide with those carried by the southward Brazil Current. This process forms the highly energetic Brazil–Malvinas Confluence—the region with the high standard deviation of SSH (Figure 1) but very slow eddy propagation (Figure 3 and Figure 4).
The most remarkable feature of eddy dynamics in the Argentine Basin is the anticyclonic (anticlockwise) propagation of eddies around the Zapiola Rise (ZR). The ZR is a seamount that rises about 1500 m above the Argentine Abyssal Plain (Figure 4). A very strong anticyclonic recirculation gyre extending from surface to bottom and transporting about 100 Sv exists around this topographic anomaly [38]. The eddy propagation around the ZR mimics the recirculation gyre with velocities comparable to the mean flow. The eddy propagation velocities on the northern flank of the seamount are the strongest in the South Atlantic Ocean and often exceed 10 km/day. On the southern flank of the seamount, eddy propagation is slower with velocities of 1–3 km/day. It is interesting to note that the mean SSH slope is steeper to the south of the seamount crest than to the north (see contours in Figure 1), meaning that the eastward surface geostrophic velocity is greater south of the ZR than the westward surface geostrophic velocity north of the ZR. This suggests that eddy propagation velocities around the ZR are mostly affected by currents at depth rather than by the near-surface flow.
The ACC in the Southern Ocean has a sufficiently strong eastward component that reverses the inherent westward direction of eddy propagation and forces them to move in the direction of the mean flow. As illustrated in Figure 3, the zonal eddy propagation velocities within the ACC are much more variable in space than those within the South Atlantic subtropical gyre. Some regions of the ACC show very slow either eastward or westward eddy propagation, but the Drake Passage, the Scotia Sea, and the region east of 20°W are characterized by relatively swift eastward eddy propagation with speeds reaching 4–5 km/day. The latter segment of the ACC is zoomed in Figure 5, showing the contours of the time-mean SSH overlaid on the mean eastward eddy propagation velocity.
It appears that the intensification of the eastward eddy propagation happens in the locations where the ACC is also stronger, associated with steeper SSH gradients and, consequently, stronger surface geostrophic flows. There are two bands of steep SSH gradients or fronts within this segment of the ACC: (i) the Sub-Antarctic Front (SAF) centered approximately along −0.2 m mean SSH contour (red contour in Figure 5) and (ii) the Polar Front (PF) centered at about –0.6 m mean SSH contour (blue contour in Figure 5) [45,46].

3.2. Time-Variable Eddy Propagation

As noted earlier, the OF method estimates apparent motion between two consecutive images. By computing the OF between each pair of SSH maps (1075 pairs in total), we can investigate the temporal evolution of the eddy propagation velocity. Changes in the eddy propagation velocity can occur due to changes in the background flow and/or changes in the stratification. Chelton et al. [47] concluded that the effects of temporal variations of the stratification can be neglected for the purposes of estimating the phase speed of Rossby waves. They noted, however, that the lack of hydrographic time series at fixed locations make the quantitative estimates of the importance of temporal variability of the stratification difficult to obtain. In the following, we evaluate the temporal evolution of the eddy propagation velocity estimated with the OF method in the ACC region considered above (Figure 5). The region was selected because of the significant influence of the background flow that reverses the intrinsic westward eddy propagation. Therefore, by comparing the variability of the eddy propagation velocity with the variability of the SSH gradient or the geostrophic flow in the region, it is possible to analyze how robust the time-dependent part of the eddy propagation velocity is.
It appears that “instantaneous” velocities are rather noisy and some spatial and temporal filtering is required. The noise can be due to the measurement and optimal interpolation errors of SSH maps, as well as due to the local violation of the brightness-constancy assumption in highly energetic regions. Displayed in Figure 6a is the time–latitude diagram of the zonal component of the eddy propagation velocity averaged between 10°W and 10°E and low-pass filtered with a Lowess filter with a cutoff period of 3 months. Overlaid are the contours of the largest meridional SSH gradients. In the zonal segment considered, both the eastward eddy propagation velocity and the SSH gradient exhibit two distinct bands between about 44–46°S and 48–50°S associated with the SAF and the PF, respectively. Here, the eastward eddy propagation velocities are 2–4 km/day and the meridional SSH gradients often exceed 0.2 m per 100 km. Further downstream, these two fronts converge at 30–40°E (Figure 5).
The meridional position of the two fronts has remained rather stationary during the period of observations, but both the SAF and the PF have intensified over the observational period, as suggested by the SSH gradients. What is interesting is that in some periods the intensity of the fronts is correlated with the eastward eddy propagation velocity. The time series of the eddy propagation velocity and the meridional SSH gradient averaged zonally between 10°W–10°E and meridionally between 44–46°S for the SAF and 48–50°S for the PF are shown in Figure 6b and Figure 6c, respectively. The correlations between the eddy propagation velocity and the SSH gradient are 0.27 for the SAF and 0.46 for the PF, both significant at 95% confidence. Particularly good agreement between the eddy propagation velocity and the SAF intensity is observed in 2001–2009 and in 2015–2021 (Figure 6b). The agreement between the eddy propagation velocity and the PF intensity is more pronounced at interannual time scales, showing the general tendency for the intensification of both the SSH gradient and the eastward-propagation velocity (Figure 6c). Although there is no perfect match between the time series of the collocated eddy propagation velocity and the SSH gradient, the reasonable agreement between the two independent parameters suggests that the optical flow can detect the temporal evolution of the eddy propagation velocity, at least in the regions with strong background flow where the signal-to-noise ratio is sufficiently large.

4. Validation of the Optical Flow Estimates

The detection of westward-propagating eddies and Rossby waves in satellite altimetry measurements has been based on the analysis of time–longitude or Hovmöller diagrams of SSH anomalies. A time–longitude diagram is a two-dimensional plot that represents a time series of zonal sections. Westward-propagating signals (waves) are clearly seen as tilted crests and troughs of SSH anomalies moving linearly from right to left across the diagram. The speed of these signals can be easily calculated from the slope of the diagonal alignments of crests and troughs. Displayed in Figure 7 are the time–longitude diagrams of SSH anomalies along the zonal sections of 10°W–10°E at several latitudes that we use to validate the results obtained with the OF method. The diagrams reveal the westward propagation at 15°S, 25°S, and 33°S, i.e., within the South Atlantic subtropical gyre, and the eastward propagation at 45°S, where the Sub-Antarctic Front of the ACC is located.
We use the two-dimensional Radon transform (2D-RT; e.g., [48]) to objectively estimate the speed of propagating signals in the selected time–longitude plots (Figure 7). The 2D-RT has been widely used in satellite remote sensing for the global analysis of Rossby waves (e.g., [49,50,51,52]). The method consists of projecting the Hovmöller diagram onto a line at a range of angles from 0 to 179° (the time axis is at 90°). When the line is perpendicular to the alignment of crests and troughs of the Rossby waves in the diagram, the projected SSH anomalies have the maximum variance (energy). The angle of the perpendicular to the line with the maximum variance yields the speed of propagation. The 2D-RT is best applied to the regions where the speed is reasonably constant within the selected subdomain. This is why the longitudinal extent of the subdomain used in this study was limited to 20°. The zonal eddy propagation velocities between 10°W–10°E estimated with the 2D-RT are −5.3 km/day at 15°S, −3.6 km/day at 3.6°S, −2.8 km/day at 33°S, and 2.2 km/day at 45°S (positive eastward).
In addition, we estimate the eddy propagation velocities using the space–time-lagged cross-correlation analysis similar to [22,23]. In essence, this method is based on computing correlations between the SSH anomalies at each grid point with the SSH anomalies at all the neighboring locations and at various time lags. At each time lag, the location of the maximum correlation is determined, and a velocity is estimated from the time lag and the distance of the location from the origin. An average velocity vector weighted by the correlation coefficients is ultimately computed from the estimates at various time lags. To focus on the mesoscale eddies, the time lags are limited to less than 100 days. The zonal and meridional dimensions of the box, within which the correlations are computed, are set to latitude-dependent correlation scales from [53], ranging from 100 km at 60°S to about 270 km at 10°S.
Displayed in Figure 8 is a comparison of the eddy propagation velocities between 10°W–10°E estimated with the OF method (blue curve), the 2D-RT (red stars), and the space–time-lagged cross-correlation analysis (red curve). The 10°W–10°E segment clearly captures the westward propagation of eddies in the South Atlantic subtropical gyre and the eastward propagation of eddies within the ACC. There is an excellent agreement between the OF estimates and the other two methods, especially south of 20°S. The largest discrepancy is observed in the tropics, where the OF estimates have the strongest variability (standard deviation of the zonal velocity is about 3 km/day). The zonal phase speed of the first-mode nondispersive long baroclinic Rossby waves between 10°W–10°E in the South Atlantic Ocean, estimated from the global atlas of the first baroclinic Rossby radius of deformation [46], is also shown in Figure 8. There is a reasonable agreement between the observed propagation velocities and the ones predicted by the classical theory at mid-latitudes. At high latitudes, the eddy propagation velocity is governed by the ACC and nearly all eddies are advected eastward, which is not accounted for in the theoretical estimate. Equatorward of 25°S, the westward eddy propagation is much slower than Rossby waves, which is consistent with earlier studies [15,23,49].

5. Summary and Conclusions

The pathways of mesoscale eddies are an important informative characteristic of ocean dynamics. The objective of this paper is to introduce an application of the optical flow (OF) technique to satellite altimetry data to estimate the direction and speed of transient eddy propagation in the South Atlantic Ocean. The South Atlantic Ocean and the adjacent segment of the Southern Ocean feature several regions with specific eddy variability, such as the Brazil–Malvinas Confluence, the Zapiola Anticyclone, the Agulhas Retroflection, and the ACC and, therefore, represent an ideal polygon for testing the OF method. The method does not detect/track individual eddies or other SSH patterns. Instead, it provides a description of eddy motion in the Eulerian frame of reference, similar to the space–time-lagged correlation analysis of SSH fields implemented earlier in [22,23].
The obtained time–mean pattern of eddy propagation velocities (Figure 3) is in an excellent agreement with the results of [22], obtained for the same region (see Figure 1 in [22]). This pattern illustrates the latitude-dependent westward propagation of Rossby waves, as well as the interaction of eddies with ocean currents and with the bottom topography. For example, the westward propagation of eddies is faster in the northern part of the subtropical gyre, where it is reinforced by the northwestward background flow, and it is much slower in the southern part of the gyre, where the mean flow opposes and, thus, slows down the westward propagation of eddies. The eddy propagation in the Argentine Basin, including the Zapiola Anticyclone, provides a great illustration of the steering effect of bottom topography on eddy motion. Another example of the eddy–mean flow interaction is the generally eastward propagation of eddies in the ACC, where the strong background flow reverses the intrinsic westward propagation and forces eddies to move in the direction of the current. The eddy propagation velocities estimated with the OF method are validated by the 2D-RT of time–longitude diagrams at several latitudes and by carrying out the space–time-lagged correlation analysis of SSH fields (similar to [22]). The obtained estimates agree with each other very well, providing strong evidence for the utility of the OF algorithm in oceanographic applications to study the motion of eddies and other sea-surface patterns.
An advantage of the OF over the methods used in earlier studies, such as the Lagrangian tracking of eddies and the space–time-lagged correlation analysis, is that it allows the computation of the time-variable part of the eddy propagation velocity. While it has been reported that changes in the stratification probably have a negligible effect on the phase speed of westward propagating signals [47], the variability of strong currents, such as the ACC, does influence the eddy propagation velocity. Therefore, in the regions of strong background flow, the use of the OF to estimate the time-variable eddy motion might be particularly useful. Note that while changes in the surface geostrophic flow are effectively measured by satellite altimetry, changes in the eddy propagation velocity mainly reflect changes in the deep flow, which are not readily available from existing observations. We find that the “instantaneous” velocities calculated by the OF between two consecutive satellite images can be rather noisy, particularly in the regions of strong variability. Therefore, the removal of outliers and some low-pass spatial and temporal filtering of the resulting (u, v) fields might be required.
The validation of the time-variable part of the eddy propagation velocity is problematic because there are no relevant observations or estimates that can be used as the ground truth. In this work, to validate the result obtained with the OF, we use the meridional SSH gradient in the ACC region, which is proportional to the surface geostrophic flow. We show that although the match between the time series of the eddy propagation velocity and the SSH gradient cannot be perfect, a reasonable agreement between the time series indicates that the time-dependent part of the eddy propagation velocity is dominated by the physical signal rather than by noise due to data errors and/or the violation of the brightness-constancy assumption. Therefore, we conclude that the OF method can be used for the determination of time-variable eddy propagation patterns. The OF provides the potential to study changes in the mid-depth and deep circulation in the regions where strong background flow significantly affects the propagation of eddies. Further research is needed to characterize the sources of errors in the OF computation in more detail and to find an optimal method for noise reduction.

Author Contributions

Conceptualization, D.L.V. and S.N.; methodology, D.L.V. and S.N.; software, D.L.V. and S.N.; validation, D.L.V.; formal analysis, D.L.V.; investigation, D.L.V.; data curation, D.L.V.; writing—original draft preparation, D.L.V.; writing—review and editing, D.L.V. and S.N.; visualization, D.L.V.; funding acquisition, D.L.V. All authors have read and agreed to the published version of the manuscript.

Funding

D.L.V. was supported by the National Oceanic and Atmospheric Administration (NOAA) Climate Variability and Predictability program (grant number NA20OAR4310407) and by the NOAA Atlantic Oceanographic and Meteorological Laboratory. This research was carried out in part under the auspices of the Cooperative Institute for Marine and Atmospheric Studies, a cooperative institute of the University of Miami and NOAA (cooperative agreement number NA20OAR4320472).

Data Availability Statement

Publicly available datasets were used in this study. The daily maps of SSH were downloaded from the CMEMS portal (http://marine.copernicus.eu; accessed on 2 May 2023). The study also used the global atlas of the first baroclinic Rossby radius of deformation available at https://ceoas.oregonstate.edu/rossby_radius (accessed on 6 May 2023) [46].

Acknowledgments

The authors are grateful to Matthieu Le Henaff for reading the initial version of the manuscript and providing helpful suggestions. Comments and suggestions from three anonymous reviewers are greatly appreciated. Maps in Figure 1, Figure 3 and Figure 4 were plotted using the M_Map package for Matlab, available online at www.eoas.ubc.ca/~rich/map.html, accessed on 20 June 2023 [54].

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Stammer, D. Global characteristics of ocean variability estimated from regional TOPEX/POSEIDON altimeter measurements. J. Phys. Oceanogr. 1997, 27, 1743–1769. [Google Scholar] [CrossRef]
  2. Ducet, N.; Le Traon, P.Y.; Reverdin, G. Global high resolution mapping of ocean circulation from the combination of TOPEX/POSEIDON and ERS-1/2. J. Geophys. Res. 2000, 105, 19477–19498. [Google Scholar] [CrossRef]
  3. Le Traon, P.Y.; Morrow, R. Ocean currents and eddies. In Satellite Altimetry and Earth Sciences: A Handbook of Techniques and Applications; Fu, L.-L., Cazenave, A., Eds.; Academic: San Diego, CA, USA, 2001; pp. 171–215. [Google Scholar]
  4. Sun, W.; Yang, J.; Tan, W.; Liu, Y.; Zhao, B.; He, Y.; Dong, C. Eddy diffusivity and coherent mesoscale eddy analysis in the Southern Ocean. Acta Oceanol. Sin. 2021, 40, 1–16. [Google Scholar] [CrossRef]
  5. Dong, C.; McWilliams, J.; Liu, Y.; Chen, D. Global heat and salt transports by eddy movement. Nat. Commun. 2014, 5, 3294. [Google Scholar] [CrossRef] [Green Version]
  6. Zhang, Z.; Wang, W.; Qiu, B. Oceanic mass transport by mesoscale eddies. Science 2014, 345, 322–324. [Google Scholar] [CrossRef]
  7. Volkov, D.L.; Lee, T.; Fu, L.-L. Eddy-induced meridional heat transport in the ocean. Geophys. Res. Lett. 2008, 35, L20601. [Google Scholar] [CrossRef] [Green Version]
  8. Qiu, B.; Chen, S.; Wu, L.; Kida, S. Wind- versus eddy-forced regional sea level trends and variability in the North Pacific Ocean. J. Clim. 2015, 28, 1561–1577. [Google Scholar] [CrossRef] [Green Version]
  9. Okubo, A. Horizontal dispersion of floatable particles in the vicinity of velocity singularities such as convergences. Deep-Sea Res. 1970, 17, 445–454. [Google Scholar] [CrossRef]
  10. Weiss, J. The dynamics of enstrophy transfer in two-dimensional hydrodynamics. Phys. D 1991, 48, 273–294. [Google Scholar] [CrossRef]
  11. Isern-Fontanet, J.; Garcia-Ladona, E.; Font, J. Identification of marine eddies from altimetric maps. J. Atmos. Ocean. Technol. 2003, 20, 772–778. [Google Scholar] [CrossRef]
  12. Chelton, D.B.; Schlax, M.G.; Samelson, R.M.; de Szoeke, R.A. Global observations of large oceanic eddies. Geophys. Res. Lett. 2007, 34, L15606. [Google Scholar] [CrossRef]
  13. Sadarjoen, I.A.; Post, F.H. Detection, quantification and tracking of vortices using streamline geometry. Comput. Graph. 2000, 24, 333–341. [Google Scholar] [CrossRef]
  14. Chaigneau, S.; Gizolme, A.; Grados, C. Mesoscale eddies off Peru in altimeter records: Identification algorithms and eddy spatio-temporal patterns. Progr. Oceanogr. 2008, 79, 106–119. [Google Scholar] [CrossRef]
  15. Chelton, D.B.; Schlax, M.G.; Samelson, R.M. Global Observations of Nonlinear Mesoscale Eddies. Prog. Oceanogr. 2011, 91, 167–216. [Google Scholar] [CrossRef]
  16. Nencioli, F.; Dong, C.; Dickey, T.D.; Washburn, L.; 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. Ocean. Technol. 2010, 27, 564–579. [Google Scholar] [CrossRef]
  17. Liu, Y.; Dong, C.; Guan, Y.; Chen, D.; McWilliams, J.; Nencioli, F. Eddy analysis in the subtropical zonal band of the North Pacific Ocean. Deep Sea Res. 2012, 68, 54–67. [Google Scholar] [CrossRef]
  18. Halo, I.; Backeburg, B.; Penven, P.; Ansorge, I.; Reason, C.; Ullgren, J.E. Eddy properties in the Mozambique Channel: A comparison between observations and two numerical ocean circulation models. Deep-Sea Res. II 2014, 100, 38–53. [Google Scholar] [CrossRef]
  19. Doglioli, A.M.; Blanke, B.; Speich, S.; Lapeyre, G. Tracking coherent structures in a regional ocean model with wavelet analysis: Application to cape basin eddies. J. Geophys. Res. 2007, 112, C05043. [Google Scholar] [CrossRef] [Green Version]
  20. Rubio, A.; Blanke, B.; Speich, S.; Grima, N.; Roy, C. Mesoscale eddy activity in the southern Benguela upwelling system from satellite altimetry and model data. Prog. Oceanogr. 2009, 83, 288–295. [Google Scholar] [CrossRef]
  21. Beron-Vera, F.J.; Wang, Y.; Olascoaga, M.J.; Goni, G.J.; Haller, G. Objective detection of oceanic eddies and the Agulhas leakage. J. Phys. Oceanogr. 2013, 43, 1426–1438. [Google Scholar] [CrossRef]
  22. Fu, L.-L. Pathways of eddies in the South Atlantic Ocean revealed from satellite altimeter observations. Geophys. Res. Lett. 2006, 33, L14610. [Google Scholar] [CrossRef]
  23. Fu, L.-L. Pattern and velocity of propagation of the global ocean eddy variability. J. Geophys. Res. 2009, 114, C11017. [Google Scholar] [CrossRef] [Green Version]
  24. Szeliski, R. Computer Vision: Algorithms and Applications, 2nd ed.; Springer: Berlin/Heidelberg, Germany, 2022; p. 947. ISBN -10: 3030343715. [Google Scholar]
  25. Keil, C.; Craig, G.C. A Displacement-Based Error Measure Applied in a Regional Ensemble Forecasting System. Mon. Wea. Rev. 2007, 135, 3248–3259. [Google Scholar] [CrossRef] [Green Version]
  26. Marzban, C.; Sandgathe, S. Optical Flow for Verification. Weather. Forecast. 2010, 25, 1479–1494. [Google Scholar] [CrossRef] [Green Version]
  27. Bowler, N.E.H.; Pierce, C.E.; Seed, A. Development of a precipitation nowcasting algorithm based on optical flow techniques. J. Hydrol. 2004, 288, 74–91. [Google Scholar] [CrossRef]
  28. Petrou, Z.I.; Xian, Y.; Tian, Y. Towards breaking the spatial resolution barriers: An optical flow and super-resolution approach for sea ice motion estimation. ISPRS J. Photogramm. Remote Sens. 2018, 138, 164–175. [Google Scholar] [CrossRef]
  29. Li, H.; Ke, C.; Zhu, Q.; Shen, X.; Li, M. An improved optical flow method to estimate Arctic sea ice velocity (winter 2014−2016). Acta Oceanol. Sin. 2021, 40, 148–160. [Google Scholar] [CrossRef]
  30. Ghalenoei, E.; Hasanlou, M. Monitoring of sea surface currents by using sea surface temperature and satellite altimetry data in the Caspian Sea. Earth Obs. Geomat. Eng. 2017, 1, 36–46. [Google Scholar] [CrossRef]
  31. Sun, H.; Song, Q.; Shao, R.; Schlicke, T. Estimation of sea surface currents based on ocean color remote-sensing image analysis. Int. J. Remote Sens. 2016, 37, 5105–5121. [Google Scholar] [CrossRef]
  32. Cunningham, S.A.; Alderson, S.G.; King, B.A.; Brandon, M.A. Transport and variability of the Antarctic Circumpolar Current in Drake Passage. J. Geophys. Res. 2003, 108, 8084. [Google Scholar] [CrossRef] [Green Version]
  33. Combes, V.; Matano, R.P. A two-way nested simulation of the oceanic circulation in the Southwestern Atlantic. J. Geophys. Res. Ocean. 2014, 119, 731–756. [Google Scholar] [CrossRef]
  34. Frey, D.I.; Piola, A.R.; Krechik, V.A.; Fofanov, D.V.; Morozov, E.G.; Silvestrova, K.P.; Tarakanov, R.Y.; Gladyshev, S.V. Direct measurements of the Malvinas Current velocity structure. J. Geophys. Res. Ocean. 2021, 126, e2020JC016727. [Google Scholar] [CrossRef]
  35. Piola, A.R.; Franco, B.C.; Palma, E.D.; Saraceno, M. Multiple jets in the Malvinas Current. J. Geophys. Res. Ocean. 2013, 118, 2107–2117. [Google Scholar] [CrossRef] [Green Version]
  36. Chidichimo, M.P.; Piola, A.R.; Meinen, C.S.; Perez, R.C.; Campos, E.J.D.; Dong, S.; Lumpkin, R.; Garzoli, S. Brazil Current volume transport variability during 2009–2015 from a long-term moored array at 34.5°S. J. Geophys. Res. Ocean. 2021, 126, e2020JC017146. [Google Scholar] [CrossRef]
  37. Piola, A.R.; Figueroa, H.A.; Bianchi, A.A. Some aspects of the surface circulation south of 20°S revealed by First GARP Global Experiment Drifters. J. Geophys. Res. 1987, 92, 5101–5114. [Google Scholar] [CrossRef]
  38. Saunders, P.M.; King, B.A. Bottom currents derived from a shipborne ADCP on the WOCE Cruise A11 in the South Atlantic. J. Phys. Oceanogr. 1995, 25, 329–347. [Google Scholar] [CrossRef]
  39. Volkov, D.L.; Fu, L.-L. The role of vorticity fluxes in the dynamics of the Zapiola Anticyclone. J. Geophys. Res. 2008, 113, C11015. [Google Scholar] [CrossRef] [Green Version]
  40. Jousset, S.; Mulet, S.; Wilkin, J.; Greiner, E.; Dibarboure, G.; Picot, N. New global Mean Dynamic Topography CNES-CLS-22 combining drifters, hydrological profiles and high frequency radar data. Ocean. Surf. Topogr. Sci. Team 2022, 31. [Google Scholar] [CrossRef]
  41. Pujol, M.-I.; Faugeère, Y.; Taburet, G.; Dupuy, S.; Pelloquin, C.; Ablain, M.; Picot, N. DUACS DT2014: The new multi-mission altimeter data set reprocessed over 20 years. Ocean. Sci. 2016, 12, 1067–1090. [Google Scholar] [CrossRef] [Green Version]
  42. Ballarotta, M.; Ubelmann, C.; Pujol, M.-I.; Taburet, G.; Fournier, F.; Legeais, J.-F.; Faugère, Y.; Delepoulle, A.; Chelton, D.; Dibarboure, G.; et al. On the resolutions of ocean altimetry maps. Ocean. Sci. 2019, 15, 1091–1109. [Google Scholar] [CrossRef] [Green Version]
  43. Lucas, B.D.; Kanade, T. An iterative image registration technique with an application to stereo vision. Proc. Imaging Underst. Workshop 1981, 2, 121–130. [Google Scholar]
  44. Bouguet, J.-Y. Pyramidal Implementation of the Lucas Kanade Feature Tracker, Description of the Algorithm; Technical Report; Intel Corporation, Microprocessor Research Labs: 1999. Available online: http://robots.stanford.edu/cs223b04/algo_affine_tracking.pdf (accessed on 20 June 2023).
  45. Sokolov, S.; Rintoul, S.R. Circumpolar structure and distribution of the Antarctic circumpolar fronts: 1. Mean circumpolar paths. J. Geophys. Res. 2009, 114, C11018. [Google Scholar] [CrossRef]
  46. Volkov, D.L.; Zlotnicki, V. Antarctic Circumpolar Current fronts observed by GOCE and GRACE missions. Ocean. Dyn. 2012, 62, 893–905. [Google Scholar] [CrossRef]
  47. Chelton, D.B.; deSzoeke, R.A.; Schlax, M.G.; El Naggar, K.; Siwertz, N. Geographical Variability of the First Baroclinic Rossby Radius of Deformation. J. Phys. Oceanogr. 1998, 28, 433–460. [Google Scholar] [CrossRef]
  48. Challenor, P.G.; Cipollini, P.; Cromwell, D. Use of the 3D Radon Transform to Examine the Properties of Oceanic Rossby Waves. J. Atmos. Ocean. Technol. 2001, 18, 1558–1566. [Google Scholar] [CrossRef]
  49. Chelton, D.B.; Schlax, M.G. Global observations of oceanic Rossby waves. Science 1996, 272, 234–238. [Google Scholar] [CrossRef]
  50. Polito, P.S.; Cornillon, P. Long baroclinic Rossby waves detected by TOPEX/POSEIDON. J. Geophys. Res. 1997, 102, 3215–3235. [Google Scholar] [CrossRef]
  51. Cipollini, P.; Jones, M.S.; Quartly, G.D.; Challenor, P.G. Concurrent TOPEX/POSEIDON altimeter and Along-Track Scanning Radiometer observations of Rossby wave propagation near 34N in the northeast Atlantic. Geophys. Res. Lett. 1997, 24, 889–892. [Google Scholar] [CrossRef]
  52. Cipollini, P.; Cromwell, D.; Quartly, G.D. Observations of Rossby Wave propagation in the Northeast Atlantic with TOPEX/POSEIDON altimetry. Adv. Space Res. 1999, 22, 1553–1556. [Google Scholar] [CrossRef]
  53. Le Traon, P.Y.; Nadal, F.; Ducet, N. An Improved Mapping Method of Multisatellite Altimeter Data. J. Atmos. Ocean. Technol. 1998, 15, 522–534. [Google Scholar] [CrossRef]
  54. Pawlowicz, R. M_Map: A Mapping Package for MATLAB, Version 1.4m, [Computer Software]. 2020. Available online: www.eoas.ubc.ca/~rich/map.html (accessed on 6 May 2023).
Figure 1. The map of the study region with the color shading showing the standard deviation of the mesoscale sea surface height variability (in meters) and the contours representing the mean dynamic topography every 0.1 m. The general upper—ocean circulation is sketched by arrows. The Agulhas eddies are schematized by red circles. The white dashed lines show segments along which the time–longitude diagrams are analyzed in Section 4. The blank areas are places with data gaps due to seasonal sea-ice cover. Abbreviations: ACC—the Antarctic Circumpolar Current, BMC—Brazil–Malvinas Confluence, ZR—Zapiola Rise, ZA—Zapiola Anticyclone, AR—Agulhas Retroflection, ARC—Agulhas Return Current, DP—Drake Passage.
Figure 1. The map of the study region with the color shading showing the standard deviation of the mesoscale sea surface height variability (in meters) and the contours representing the mean dynamic topography every 0.1 m. The general upper—ocean circulation is sketched by arrows. The Agulhas eddies are schematized by red circles. The white dashed lines show segments along which the time–longitude diagrams are analyzed in Section 4. The blank areas are places with data gaps due to seasonal sea-ice cover. Abbreviations: ACC—the Antarctic Circumpolar Current, BMC—Brazil–Malvinas Confluence, ZR—Zapiola Rise, ZA—Zapiola Anticyclone, AR—Agulhas Retroflection, ARC—Agulhas Return Current, DP—Drake Passage.
Remotesensing 15 03894 g001
Figure 2. (a) An example of an anticyclonic eddy propagating westward over the 10-day time interval (∆t). The black and blue rectangles show 5 × 5 pixel patches, over which the flow is assumed to be constant. During ∆t, the patch shown by the black rectangle has moved in the direction shown by the black and blue arrows to the position outlined by the blue rectangle. The red rectangles outline individual pixels η1 and η2. (b) An illustration of the pyramid method with three layers. The optical flow (OF) is computed at each pyramid layer starting from the top layer at the lowest resolution down to the original image.
Figure 2. (a) An example of an anticyclonic eddy propagating westward over the 10-day time interval (∆t). The black and blue rectangles show 5 × 5 pixel patches, over which the flow is assumed to be constant. During ∆t, the patch shown by the black rectangle has moved in the direction shown by the black and blue arrows to the position outlined by the blue rectangle. The red rectangles outline individual pixels η1 and η2. (b) An illustration of the pyramid method with three layers. The optical flow (OF) is computed at each pyramid layer starting from the top layer at the lowest resolution down to the original image.
Remotesensing 15 03894 g002
Figure 3. The time–mean vectors of eddy propagation velocity obtained with the optical flow method. The color shading indicates the absolute speed in km per day.
Figure 3. The time–mean vectors of eddy propagation velocity obtained with the optical flow method. The color shading indicates the absolute speed in km per day.
Remotesensing 15 03894 g003
Figure 4. Bottom topography based on ETOPO2 Global Relief Model (color shading) and eddy propagation velocities (vectors) in the Argentine Basin. The black and blue contours denote the 200-m and 4000-m isobaths, respectively. Abbreviations: BMC—Brazil–Malvinas Confluence, AAP—Argentine Abyssal Plain.
Figure 4. Bottom topography based on ETOPO2 Global Relief Model (color shading) and eddy propagation velocities (vectors) in the Argentine Basin. The black and blue contours denote the 200-m and 4000-m isobaths, respectively. Abbreviations: BMC—Brazil–Malvinas Confluence, AAP—Argentine Abyssal Plain.
Remotesensing 15 03894 g004
Figure 5. The Southern Ocean segment of the study domain: (color shading) zonal (eastward) eddy propagation velocity (km/day) and the time–mean SSH (contours). The approximate locations of the SAF and the PF are shown by the red (–0.2 m) and blue (–0.6 m) contours, respectively.
Figure 5. The Southern Ocean segment of the study domain: (color shading) zonal (eastward) eddy propagation velocity (km/day) and the time–mean SSH (contours). The approximate locations of the SAF and the PF are shown by the red (–0.2 m) and blue (–0.6 m) contours, respectively.
Remotesensing 15 03894 g005
Figure 6. (a) The time–latitude diagram of (color shading) the zonal eddy propagation velocity (positive eastward) and (contours) the meridional SSH gradient averaged between 10°W–10°E. The black and blue contours show the meridional SSH gradients of 0.2 m per 100 km and 0.25 m per 100 km, respectively. (b,c) The time series of (blue) the zonal eddy propagation velocity and (red) the meridional SSH gradient averaged between 10°W–10°E and between (b) 44–46°S (Sub-Antarctic Front) and (c) 48–50°S (Polar Front).
Figure 6. (a) The time–latitude diagram of (color shading) the zonal eddy propagation velocity (positive eastward) and (contours) the meridional SSH gradient averaged between 10°W–10°E. The black and blue contours show the meridional SSH gradients of 0.2 m per 100 km and 0.25 m per 100 km, respectively. (b,c) The time series of (blue) the zonal eddy propagation velocity and (red) the meridional SSH gradient averaged between 10°W–10°E and between (b) 44–46°S (Sub-Antarctic Front) and (c) 48–50°S (Polar Front).
Remotesensing 15 03894 g006
Figure 7. Time–longitude (Hovmöller) diagrams of the high-pass filtered SSH for the zonal segments between 10°W–10°E at 15°S, 25°S, 33°S, and 45°S (the segments are shown by the white dashed lines in Figure 1). The color shading shows SSH anomalies in cm. The zonal eddy propagation velocities estimated using a Radon transform are shown on top of each panel.
Figure 7. Time–longitude (Hovmöller) diagrams of the high-pass filtered SSH for the zonal segments between 10°W–10°E at 15°S, 25°S, 33°S, and 45°S (the segments are shown by the white dashed lines in Figure 1). The color shading shows SSH anomalies in cm. The zonal eddy propagation velocities estimated using a Radon transform are shown on top of each panel.
Remotesensing 15 03894 g007
Figure 8. The time–mean zonal eddy propagation velocities estimated with the optical flow (OF) and averaged zonally between 10°W and 10°E (blue curve). The gray shading indicates the ±1 standard deviation. The red stars show the zonal propagation velocities of mesoscale eddies at 15°S, 25°S, 33°S, and 45°S estimated using the two-dimensional Radon transform of time–longitude (Hovmöller) diagrams at these latitudes for the zonal segments between 10°W–10°E. The corresponding time–longitude diagrams are shown in Figure 7 and the zonal segments are indicated by the white dashed lines in Figure 1. The red curve shows the zonal eddy propagation estimated with the space–time-lagged correlation analysis (STLC). The dotted black curve shows the zonally averaged (between 10°W and 10°E) propagation velocity of the first-mode long baroclinic Rossby waves (LBRW) computed using the Rossby radius of deformation from [47].
Figure 8. The time–mean zonal eddy propagation velocities estimated with the optical flow (OF) and averaged zonally between 10°W and 10°E (blue curve). The gray shading indicates the ±1 standard deviation. The red stars show the zonal propagation velocities of mesoscale eddies at 15°S, 25°S, 33°S, and 45°S estimated using the two-dimensional Radon transform of time–longitude (Hovmöller) diagrams at these latitudes for the zonal segments between 10°W–10°E. The corresponding time–longitude diagrams are shown in Figure 7 and the zonal segments are indicated by the white dashed lines in Figure 1. The red curve shows the zonal eddy propagation estimated with the space–time-lagged correlation analysis (STLC). The dotted black curve shows the zonally averaged (between 10°W and 10°E) propagation velocity of the first-mode long baroclinic Rossby waves (LBRW) computed using the Rossby radius of deformation from [47].
Remotesensing 15 03894 g008
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Volkov, D.L.; Negahdaripour, S. Implementation of the Optical Flow to Estimate the Propagation of Eddies in the South Atlantic Ocean. Remote Sens. 2023, 15, 3894. https://0-doi-org.brum.beds.ac.uk/10.3390/rs15153894

AMA Style

Volkov DL, Negahdaripour S. Implementation of the Optical Flow to Estimate the Propagation of Eddies in the South Atlantic Ocean. Remote Sensing. 2023; 15(15):3894. https://0-doi-org.brum.beds.ac.uk/10.3390/rs15153894

Chicago/Turabian Style

Volkov, Denis L., and Shahriar Negahdaripour. 2023. "Implementation of the Optical Flow to Estimate the Propagation of Eddies in the South Atlantic Ocean" Remote Sensing 15, no. 15: 3894. https://0-doi-org.brum.beds.ac.uk/10.3390/rs15153894

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop