1. Introduction
The main source of water for the Gulf of Mexico is the Yucatan Current, which enters through the Yucatan Channel from the Caribbean Sea and forms, within the Gulf, an anticyclonic retroflection or “loop” current known as the Loop Current (LC). The Yucatan Current is a baroclinic jet with the bulk of the transport above 800 m (Sheinbaum et al. 2002) flowing along the eastern coast of the Yucatan Peninsula at speeds of 1 to 2 m s−1 (Badan et al. 2005) before it enters the Gulf to form the LC. Unlike the Yucatan Current, which is confined to a nearly fixed geographic position, the LC can intrude far into the Gulf, reaching on average to 26.2°N and a maximum of 28°N in the north–south direction and 85.8°W and a maximum of 92°W in the east–west direction based on altimeter-derived LC metrics (Leben 2005). During each intrusion cycle, the LC sheds a large warm core eddy that ranges in diameter from 300 to 400 km; after separating from the LC the eddy moves westward at speeds of 1–8 km day−1 (Elliot 1982; Vukovich and Crissman 1986; Vukovich 2007). These eddies are commonly referred to as LC eddies (LCEs). Separation, therefore, refers to the detachment and ultimate separation of an LCE into the western Gulf away from the LC. The LC retracts by an amount equal to the diameter of the eddy that has been shed, but in no instance has it been observed to retreat southward of 24°N (Leben 2005). These large intrusions and rapid retreats are a unique time-dependent characteristic of LC dynamics.
The time between eddy separation events is commonly called the eddy separation period and ranges from 0.5 to 18.5 months (Vukovich 1995, 2007; Sturges and Leben 2000; Leben 2005). Eddy separation is irregular, but not chaotic (Lugo-Fernández 2007). While the cause of eddy separation remains largely unexplained, we do know that the baroclinic Yucatan Current inflow (Hurlburt and Thompson 1980) and the associated vorticity and transport fluctuations affect eddy shedding (Candela et al. 2002; Oey et al. 2003; Lugo-Fernández and Badan 2007). Pichevin and Nof (1997) proposed that eddy separation results from a northward current entering an ocean and turning eastward, because of vorticity conservation, in order to balance the momentum of the system. Another proposed mechanism for eddy separation is the shearing of the LC by cyclonic eddies that have been generated by the interaction of the current with topography (Cochrane 1972; Zavala-Hidalgo et al. 2003; Schmitz 2005; Chérubin et al. 2006). A comprehensive theory capable of predicting the highly variable eddy separation periods exhibited by individual LC intrusion events has not been found.
Leben (2005) reported a correlation between the retreat latitude of LC after eddy separation and the eddy separation period (see his Fig. 5). This author found that for retreat latitudes ≥26°N there is a strong linear relationship, but for retreat latitudes <26°N there is more scatter. In spite of the scatter, the correlation between the retreat latitude and eddy separation period was −0.83 for 16 events. Leben (2005) hypothesized that this linear relationship reflects the sensitivity of LC intrusion to the initial condition set at the time of eddy separation.
The main objectives of this study are to provide a semitheoretical justification for the observed relationship between retreat latitude and separation period and to verify the functional form of the relationship. A simple statistical model based on the observed relationship could provide a forecast tool for predicting the period of LC intrusion and eddy separation from LC retreat. The ability to accurately predict the time scale of LC intrusion and eddy separation has large practical implications for offshore oil and gas activities in the Gulf of Mexico, scientific investigations, navigation, and hurricane forecasting (Lewis et al. 1991; Shay and Uhlhorn 2006).
2. Derivation
At this point, two aspects of Eq. (1) need to be discussed. First, we assumed that H and W are constants. The justification for this assumption comes from observations and theoretical considerations: 1) Sheinbaum et al. (2002) show that H ∼ 800 m and is set by the channel depth off Miami; 2) observations in the LC confirm the two-layer approximation with an interface close to 800 m (Inoue et al. 2008; Welsh et al. 2009); and 3) analysis of hydrographic data in the Gulf also indicates that H ≈ 800 m (Sturges 2005). The constancy of W, about twice the width of the Yucatan Channel, is also set by the geometrical constraint on the LC (Leben 2005; Vukovich 2007) and velocity cross sections through the LC that show a width of ∼200 km (Nowlin and Hubertz 1972). The second aspect is whether Eq. (1) is a lower bound (Lugo-Fernández and Badan 2007). An examination of Eq. (11) in Lugo-Fernández and Badan (2007) shows that it is the amplitude of a sinusoidal variation, and as such is not a lower bound. Second, a comparison of the observed separation periods and the predictions of Eq. (11) [see Fig. 5 in Lugo-Fernández and Badan (2007)] reveals that this equation tends to overestimate T in general, which suggests again that it is not a lower bound.
3. Data and methods
4. Results
Figure 1 shows θs versus T for the 25 separation events observed from 1993 to 2009, which updates Fig. 5 from Leben (2005). First, observe that the new data display a clear linear trend also and include events around θs ∼ 25°N that were absent in Leben (2005). The correlation between θs versus T is −0.84, comparable to that in Leben (2005). Thus, the assumptions and the equality sign in Eq. (5) are supported by the data.
Plots of T versus θs and θs versus T were constructed with the data in Table 1 (Fig. 3). Note that values T and θs (upper panel of Fig. 3) are anticorrelated, as predicted by Eq. (5). The slopes obtained from the regression equations are mT = −153° day−1 and mθ = −218° day−1 (see Table 2). The geometric mean slope estimate is −183° day−1. The slope from Eq. (6), −168° day−1, is in very good agreement with the estimated geometric mean slope estimate, −183° day−1. The y intercept estimated from the regressions equals 27.3° versus the 28° assumed in deriving Eq. (5). Note, also, that the maximum observed intrusion latitude is 28.4° (Table 1), which is a little bit larger than our initial assumed value. The error estimates for the slope are ±25° day−1and ±2° for the y intercept. The agreement between the altimeter-derived results and the semiempirical relationship regression parameters is very good, with discrepancies of 9% and 2% for the slope and y intercept, respectively. Based on the error estimates, the regression parameters are statistically significant.
5. Discussion and conclusions
First, we point out that Eq. (1) represents an integration over an LC intrusion cycle and that each cycle is independent of the previous cycle; thus, we can apply Eq. (1) to each individual cycle. The independence of the LC cycle from previous cycles reflects the short-term memory of the LC system, which was discussed in more detail in Lugo-Fernández (2007). As a result, each intrusion cycle is a unique event; nevertheless, the fundamental balance of vorticity and mass conservation represented by Eq. (1) still holds. Second, the deceptively simple Eq. (5) warrants a physical interpretation. To start, κ represents the constraint imposed by the Gulf’s geometry on the northward intrusion of the LC. The slope represents the effects of the conditions at the Yucatan Channel on the LC motion and vorticity conservation. The slope also accounts for the dependence of the LC on the initial conditions as hypothesized by Leben (2005). Third, the estimated slope and y intercept from the semiempirical theory are in good agreement with the linear regression values based on the altimeter-derived LC metrics, with respectively 9% and 2% discrepancies between the estimated and regression derived values. The agreement between theory and observations suggests that this linear relationship is a fundamental physical property of the LC and arises from a straightforward conservation of mass and vorticity and the geometrical constraint imposed by the Gulf’s deepwater basin. Finally, the correlation between T and θs (upper panel of Fig. 3) provides the recommended forecasting tool for LC eddy separation in the Gulf of Mexico.
Physically, the mechanism embodied in this work is as follows. During an intrusion cycle, the Loop Current penetrates northward into the Gulf, displacing water and gaining anticyclonic vorticity, until halted by the northern Gulf continental slope. Simultaneously, the displaced water forces flow back into the Caribbean Sea through the Yucatan Channel and under the Loop Current. This outflow advects anticyclonic vorticity out of the Gulf. Thus, the time scale is determined by the balance between accumulation of vorticity in the Gulf and outflow of vorticity from the Gulf. Since the final state is imposed by the geometry of the Gulf, the intrusion cycle depends primarily on the initial state of the Loop or its southern retreat position.
Acknowledgments
A. Lugo-Fernández appreciates the support of the U.S. Department of the Interior, Bureau of Ocean Energy Management, Regulation and Enforcement (BOEMRE), Gulf of Mexico Region during the preparation of this manuscript. R. Leben acknowledges support from BOEMRE Contract M08PC20043 to Science Applications International Corp. and from NASA Ocean Surface Topography Mission Science Team Grant NNX08AR60G. The opinions expressed by the authors are their own and do not necessarily reflect the opinion or policy of the U.S. Government. Thanks to two anonymous reviewers whose comments and suggestions helped to improve this work.
REFERENCES
Badan, A., J. Candela, J. Sheinbaum, and J. Ochoa, 2005: Upper-layer circulation in the approaches to Yucatan Channel. Circulation in the Gulf of Mexico: Observations and Models, Geophys. Monogr., Vol. 161, Amer. Geophys. Union, 57–69.
Candela, J., J. Sheinbaum, J. Ochoa, A. Badan, and R. Leben, 2002: The potential vorticity flux through the Yucatan Channel and the Loop Current in the Gulf of Mexico. Geophys. Res. Lett., 29 , 2059–2062.
Candela, J., S. Tanahara, M. Crepon, B. Barnier, and J. Sheinbaum, 2003: Yucatan Channel flow: Observations versus CLIPPER ATL6 and MERCATOR PAM models. J. Geophys. Res., 108 , 3385. doi:10.1029/2003JC001961.
Chérubin, L. M., Y. Morel, and E. P. Chassignet, 2006: Loop Current ring shedding: The formation of cyclones and the effect of topography. J. Phys. Oceanogr., 36 , 569–591.
Cochrane, J. D., 1972: Separation of an anticyclone and subsequent developments in the Loop Current (1969). Contributions on the Physical Oceanography of the Gulf of Mexico, Vol. II, L. R. A. Capurro and J. L. Reid, Eds., Gulf Publishing, 91–106.
Elliot, B. A., 1982: Anticyclonic rings in the Gulf of Mexico. J. Phys. Oceanogr., 12 , 1292–1309.
Hurlburt, H. E., and J. D. Thompson, 1980: A numerical study of Loop Current intrusions and eddy shedding. J. Phys. Oceanogr., 10 , 1611–1651.
Inoue, M., S. E. Welsh, L. J. Rouse Jr., and E. Weeks, 2008: Deepwater currents in the Eastern Gulf of Mexico: Observations at 25.5°N and 87°W. OCS Study MMS 2008-001, 95 pp. [Available online at http://www.gomr.boemre.gov/PI/PDFImages/ESPIS/4/4305.pdf].
Laws, E., 1997: Mathematical Methods for Oceanographers: An Introduction. John Wiley & Sons, 343 pp.
Leben, R. R., 2005: Altimeter-derived Loop Current metrics. Circulation in the Gulf of Mexico: Observations and Models, Geophys. Monogr., Vol. 161, Amer. Geophys. Union, 181–201.
Lewis, D. B., J. B. Adams, F. E. Shanks, and D. Szabo, 1991: The Loop Current experience: Ewing Bank 871. J. Pet. Technol., 43 , 1038–1044.
Lugo-Fernández, A., 2007: Is the Loop Current a chaotic oscillator? J. Phys. Oceanogr., 37 , 1455–1469.
Lugo-Fernández, A., and A. Badan, 2007: On the vorticity cycle of the Loop Current. J. Mar. Res., 65 , 471–489.
Nowlin, W. D., and J. M. Hubertz, 1972: Contrasting summer circulation patterns for the eastern Gulf. Contributions on the Physical Oceanography of the Gulf of Mexico, Vol. II, L. R. A. Capurro and J. L. Reid, Eds., Gulf Publishing, 119–137.
Oey, L-Y., H-C. Lee, and W. J. Schmitz, 2003: Effects of winds and Caribbean eddies on the Loop Current eddy shedding: A numerical model study. J. Geophys. Res., 108 , 3324. doi:10.1029/2002JC001698.
Pedlosky, J., 1998: Ocean Circulation Theory. 2nd ed. Springer-Verlag, 453 pp.
Pichevin, T., and D. Nof, 1997: The momentum imbalance paradox. Tellus, 49A , 298–319.
Schmitz, W. J., 2005: Cyclones and westward propagation in the shedding of anticyclonic rings from the Loop Current. Circulation in the Gulf of Mexico: Observations and Models, Geophys. Monogr., Vol. 161, Amer. Geophys. Union, 241–261.
Shay, L. K., and E. W. Uhlhorn, 2006: Loop Current interactions during Hurricanes Isidore and Lilli. Proc. 27th Conf. on Hurricanes and Tropical Meteorology, Monterey, CA, Amer. Meteor. Soc., 3C.3.
Sheinbaum, J., J. Candela, A. Badan, and J. Ochoa, 2002: Flow structure and transport in the Yucatan Channel. Geophys. Res. Lett., 29 , 1040. doi:10.1029/2001GL013990.
Sturges, W., 2005: Deep-water exchange between the Atlantic, Caribbean, and Gulf of Mexico. Circulation in the Gulf of Mexico: Observations and Models, Geophys. Monogr., Vol. 161, Amer. Geophys. Union, 263–278.
Sturges, W., and R. Leben, 2000: Frequency of ring separations from the Loop Current in the Gulf of Mexico: A revised estimate. J. Phys. Oceanogr., 30 , 1814–1818.
Vukovich, F. M., 1995: An updated evaluation of the Loop Current eddy-shedding frequency. J. Geophys. Res., 100 , (C5). 8655–8659.
Vukovich, F. M., 2007: Climatology of ocean features in the Gulf of Mexico using satellite remote sensing data. J. Phys. Oceanogr., 37 , 689–707.
Vukovich, F. M., and B. W. Crissman, 1986: Aspects of warm rings in the Gulf of Mexico. J. Geophys. Res., 91 , 2645–2660.
Welsh, S. E., M. Inoue, L. J. Rouse Jr., and E. Weeks, 2009: Observation of the deepwater manifestation of the Loop Current and Loop Current rings in the eastern Gulf of Mexico. OCS Study MMS 2009-050, 110 pp.
Zavala-Hidalgo, J., S. L. Morey, and J. J. O’Brien, 2003: Cyclonic eddies northeast of the Campeche Bank from altimetry data. J. Phys. Oceanogr., 33 , 623–629.
Correlation between separation period T and retreat latitude θs for the Loop Current based on altimetry observations from September 1993 to March 2009. The correlation between θs and T is −0.84 and the linear fit has an R2 of 0.70.
Citation: Journal of Physical Oceanography 40, 12; 10.1175/2010JPO4354.1
(top) Lower-layer transport in the Yucatan Channel and (bottom) its integrated transport anomaly from 10 Sep 1999 through 31 May 2001. Data from Candela et al. (2003).
Citation: Journal of Physical Oceanography 40, 12; 10.1175/2010JPO4354.1
Regressions of (top) eddy separation period T vs retreat latitude θs and (bottom) θs vs T using the CCAR altimetry data in Table 1.
Citation: Journal of Physical Oceanography 40, 12; 10.1175/2010JPO4354.1
Updated information on eddy separation period and retreat latitude from the Colorado Center for Astrodynamics Research (CCAR) altimetry data between September 1993 to March 2009.
Regression results for T vs θs and θs vs T for the CCAR data. The asterisk (*) indicates the geometric mean.