Eddy Statistics Research Articles

Reply to Comment by Chang on “Anticyclonic Suppression of the North Pacific Transient Eddy Activity in Midwinter”

AbstractChang (2024, https://doi.org/10.1029/2024gl110011) challenged the methodology proposed recently by Okajima et al. for evaluating cyclonic and anticyclonic contributions to Eulerian eddy statistics and atmospheric energetics based on the local flow curvature. He argued that using the local wind curvature to separate energetic contributions from cyclones and anticyclones is not physically meaningful. Here we argue that his claims are based on an unrealistic assumption of monopolar relative vorticity in an entire storm‐track domain and a meridionally uniform zonal background flow atypical to midlatitudes. We also demonstrate that the error in attributing eddy statistics to cyclones and anticyclones is significantly smaller than his estimation. Rather, we further demonstrate that the curvature‐based methodology effectively eliminates the shear influence to identify cyclonic and anticyclonic regions, which is dismissed in his argument. We conclude that the curvature‐based methodology is beneficial in evaluating distinct cyclonic and anticyclonic contributions to atmospheric energetics in realistic conditions.

Changes in the tropical upper-tropospheric zonal momentum balance due to global warming

Abstract. We study the zonal momentum budget of the deep upper tropics in the context of present and future climates. In the zonal mean, as is known, a robust balance exists between the acceleration by the horizontal eddy momentum flux convergence and the deceleration by the mean meridional momentum advection. During summer, climatological stationary Rossby waves over the Asian monsoon longitudes converge westerly momentum into the tropics and are the primary contributors to the eddy term. During winter, anomalous westerly winds over the tropical east Pacific allow extratropical waves to propagate into the deep tropics, where they tend to break and decelerate the flow. When integrated over all longitudes, we find that eddies from these two regions sum constructively in summer and destructively in winter, always yielding a net positive momentum forcing that balances the mean flow term. The state-of-the-art CMIP6 suite qualitatively captures these features in the historical run and shows that the momentum fluxes change due to global warming. In summer, stationary eddy circulations in the Asian monsoon zone weaken in the upper troposphere (UT) but strengthen in the lower stratosphere (LS). Greater upward mass flux from the troposphere forces a stronger divergence and a more intense circulation in the LS following a Sverdrup vorticity balance. We observe this strengthening of summertime tropical and subtropical stationary waves in the LS over all longitudes and verify it in an idealized aquaplanet general circulation model experiment. In winter, we find that eddy westerlies over the east Pacific longitudes decrease in strength due to the expected weakening of the subtropical stationary waves with warming. This causes a significant decrease in the propagation of extratropical waves into this region, along with a drop in eddy potential vorticity fluxes associated with these waves. Thus, apart from the mean meridional flux, which weakens due to the projected weakening of the Hadley cells, our analysis of warming simulations clearly suggests significant and robust changes in the eddy momentum fluxes in the deep tropics. Potential implications of these changes in the context of the zonal mean flow and regional circulations are discussed.

Dynamical reconstruction of the upper-ocean state in the central Arctic during the winter period of the MOSAiC expedition

Abstract. This paper presents a methodological tool for dynamic reconstruction of the state of the ocean, based, as an example, on observations from the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) experiment. The data used in this study were collected in the Amundsen Basin between October 2019 and January 2020. Analysing observational data to assess tracer field and upper-ocean dynamics is highly challenging when measurement platforms drift with the ice pack due to continuous drift speed and direction changes. We have equipped the new version of the coastal branch of the global Finite-volumE sea ice–Ocean Model (FESOM-C) with a nudging method. Model nudging was carried out assuming a quasi-steady state. Overall, the model can reproduce the lateral and vertical structure of the temperature, salinity, and density fields, which allows for projecting dynamically consistent features of these fields onto a regular grid. We identify two separate depth ranges of enhanced eddy kinetic energy located around two maxima in buoyancy frequency: the depth of the upper halocline and the depth of the warm (modified) Atlantic Water. Simulations reveal a notable decrease in surface layer salinity and density in the Amundsen Basin towards the north but no significant gradient from east to west. However, we find a mixed-layer deepening from east to west, with a 0.084 m km−1 gradient at 0.6 m km−1 standard deviation, compared to a weak deepening from south to north. The model resolves several stationary eddies in the warm Atlantic Water and provides insights into the associated dynamics. The model output can be used to further analyse the thermohaline structure and related dynamics associated with mesoscale and submesoscale processes in the central Arctic, such as estimates of heat fluxes or mass transport. The developed nudging method can be utilized to incorporate observational data from a diverse set of instruments and for further analysis of data from the MOSAiC expedition.

Trends in Atmospheric Heat Transport Since 1980

Abstract We investigate the linear trends in meridional atmospheric heat transport (AHT) since 1980 in atmospheric reanalysis datasets, coupled climate models, and atmosphere-only climate models forced with historical sea surface temperatures. Trends in AHT are decomposed into contributions from three components of circulation: (i) transient eddies, (ii) stationary eddies, and (iii) the mean meridional circulation. All reanalyses and models agree on the pattern of AHT trends in the Southern Ocean, providing confidence in the trends in this region. There are robust increases in transient-eddy AHT magnitude in the Southern Ocean in the reanalyses, which are well replicated by the atmosphere-only models, while coupled models show smaller magnitude trends. This suggests that the pattern of sea surface temperature trends contributes to the transient-eddy AHT trends in this region. In the tropics, we find large differences between mean-meridional circulation AHT trends in models and the reanalyses, which we connect to discrepancies in tropical precipitation trends. In the Northern Hemisphere, we find less evidence of large-scale trends and more uncertainty, but note several regions with mismatches between models and the reanalyses that have dynamical explanations. Throughout this work we find strong compensation between the different components of AHT, most notably in the Southern Ocean where transient-eddy AHT trends are well compensated by trends in the mean-meridional circulation AHT, resulting in relatively small total AHT trends. This highlights the importance of considering AHT changes holistically, rather than each AHT component individually.

Application of a Simple Diffusivity Formulation to Examine Regime Transition and Jet–Eddy Energy Partitioning in Quasi-Geostrophic Turbulence

Abstract This study uses a simple diffusivity formulation to examine flow regime transition and jet–eddy energy partitioning in two-layer quasigeostrophic turbulence. Guided by simulations, the formulation is empirically constructed so that the diffusivity is bounded by a f-plane asymptote (Df) in the limit of vanishing β (termed drag-controlled) while reduced to a drag-independent scaling (Dβ) of Lapeyre and Held toward large β (termed β-controlled). Good agreement is found for diffusivities diagnosed from simulations with both quadratic and linear drag and in 2D turbulence. From the formulation, a regime diagram is readily constructed, with Df/Dβ = 1 separating the drag-controlled and β-controlled regimes. The diagram also sets the parameter range where an eddy velocity scaling is applicable. The quantitative representations of eddy variables then enable a reasonably skillful theory for zonal jet speed to be developed from energy balance. It is shown that, using Df/Dβ ≥ 10, a state where eddy statistics are approximately drag insensitive could be identified and interpreted using wave-damping competitions in slowing an inverse cascade. However, contrary to an existing hypothesis, the energy dissipation in such a state is not dominated by zonal jets. A modest revision for a way to maintain balance while keeping eddies drag insensitive is proposed. In the regime diagram, a subspace of zonostrophic condition, defined as jet dissipation surpassing eddy, is further quantified. It is demonstrated that a rough scaling could help interpret how the relative importance of jet and eddy dissipation varies across the parameter space.

Anticyclonic Suppression of the North Pacific Transient Eddy Activity in Midwinter

AbstractDynamical understandings of midlatitude transient eddy activity, especially its midwinter minimum over the North Pacific, are still limited, partly because conventional Eulerian eddy statistics are incapable of separating cyclonic and anticyclonic contributions. Here we evaluate the two contributions separately based on local curvature of instantaneous flow fields to compare their seasonality between the North Pacific and North Atlantic storm‐tracks. The anticyclonic contribution is found crucial for the midwinter minimum of the North Pacific transient eddy activity. Eddy energetics reveals that the net efficiency of the anticyclonic contribution in replenishing total transient eddy energy over the North Pacific exhibits a pronounced midwinter minimum leading to net energy loss, while that of its cyclonic counterpart does not, in harmony with a precipitation peak around midwinter. This study suggests that more attention should be paid to anticyclones in studying midlatitude storm‐track dynamics.

A Detailed Limited-Area Atmospheric Energy Cycle for Climate and Weather Studies

Lorenz’ seminal work on global atmospheric energetics improved our understanding of the general circulation. With the advent of Regional Climate Models (RCMs), it is important to have a limited-area energetic budget available that is applicable for both weather and climate, analogous to Lorenz’ global atmospheric energetics. A regional-scale energetic budget is obtained in this study by applying Reynolds decomposition rules to quadratic forms of the kinetic energy K and the available enthalpy A, to obtain time mean and time deviation contributions. According to the employed definition, the time mean energy contributions are decomposed in a component associated with the time-averaged atmospheric state and a component due to the time-averaged statistics of transient eddies; these contributions are suitable for the study of the climate over a region of interest. Energy fluctuations (the deviations of instantaneous energies from their climate value) that are appropriate for weather studies are split into quadratic and linear contributions. The sum of all the contributions returns exactly to the total primitive kinetic energy and available enthalpy equations.

Distinct Roles of Cyclones and Anticyclones in Setting the Midwinter Minimum of the North Pacific Eddy Activity: A Lagrangian Perspective

Abstract The North Pacific storm-track activity is suppressed substantially under the excessively strong westerlies to form a distinct minimum in midwinter, which seems inconsistent with linear baroclinic instability theory. This “midwinter minimum” of the storm-track activity has been intensively investigated for decades as a test case for storm-track dynamics. However, the mechanisms controlling it are yet to be fully unveiled and are still under debate. Here we investigate the detailed seasonal evolution of the climatological density of surface migratory anticyclones over the North Pacific, in comparison with its counterpart for cyclones, based on a Lagrangian tracking algorithm. We demonstrate that the frequency of surface cyclones over the North Pacific maximizes in midwinter, whereas that of anticyclones exhibits a distinct midwinter minimum under the upstream influence, especially from the Japan Sea region. In midwinter, it is only on such a rare occasion that prominent weakening of the East Asian winter monsoon allows a migratory surface anticyclone to form over the Japan Sea, despite the unfavorable climatological-mean conditions due to persistent monsoonal cold-air outbreaks and the excessively strong upper-tropospheric westerlies. The midwinter minimum of the North Pacific anticyclone density suggests that anticyclones are likely the key to understanding the midwinter minimum of the North Pacific storm-track activity as measured by Eulerian eddy statistics.

Spatio‐Temporal Coarse‐Graining Decomposition of the Global Ocean Geostrophic Kinetic Energy

AbstractWe expand on a recent determination of the first global energy spectrum of the ocean's surface geostrophic circulation (Storer et al., 2022, https://doi.org/10.1038/s41467-022-33031-3) using a coarse‐graining (CG) method. We compare spectra from CG to those from spherical harmonics by treating land in a manner consistent with the boundary conditions. While the two methods yield qualitatively consistent domain‐averaged results, spherical harmonics spectra are too noisy at gyre‐scales (>1,000 km). More importantly, spherical harmonics are inherently global and cannot provide local information connecting scales with currents geographically. CG shows that the extra‐tropics mesoscales (100–500 km) have a root‐mean‐square (rms) velocity of ∼15 cm/s, which increases to ∼30–40 cm/s locally in the Gulf Stream and Kuroshio and to ∼16–28 cm/s in the ACC. There is notable hemispheric asymmetry in mesoscale energy‐per‐area, which is higher in the north due to continental boundaries. We estimate that ≈25%–50% of total geostrophic energy is at scales smaller than 100 km, and is un(der)‐resolved by pre‐SWOT satellite products. Spectra of the time‐mean circulation show that most of its energy (up to 70%) resides in stationary eddies with characteristic scales smaller than (<500 km). This highlights the preponderance of “standing” small‐scale structures in the global ocean due to the temporally coherent forcing by boundaries. By coarse‐graining in space and time, we compute the first spatio‐temporal global spectrum of geostrophic circulation from AVISO and NEMO. These spectra show that every length‐scale evolves over a wide range of time‐scales with a consistent peak at ≈200 km and ≈2–3 weeks.

Antarctica Slope Front bifurcation eddy: A stationary feature influencing CO2 dynamics in the northern Antarctic Peninsula

The Southern Ocean is a key region for analyzing environmental drivers that regulate sea-air CO2 exchanges. These CO2 fluxes are influenced by several mesoscale structures, such as meanders, eddies and other mechanisms responsible for energy dissipation. Aiming to better understand sea-air CO2 dynamics in the northern Antarctica Peninsula, we investigated an anticyclonic stationary eddy located south of Clarence Island, in the eastern basin of Bransfield Strait – named the Antarctica Slope Front bifurcation (ASFb) eddy. Physical, chemical and biological data were sampled, and remote sensing measurements taken, in the region during late summer conditions in February 2020. The eddy’s core consisted of cold (0.31 °C), salty (34.38) and carbon-rich (2247 μmol kg−1) waters with dissolved oxygen depletion (337 μmol kg−1). The core retains a mixture of local surface waters with waters derived from Circumpolar Deep Water (i.e., Warm Deep Water from the Weddell Sea and modified Circumpolar Deep Water from the Bransfield Strait) and Dense Shelf Water. The ASFb eddy acts as a CO2 outgassing structure that reaches a CO2 emission to the atmosphere of ∼1.5 mmol m−2 d–1 in the eddy’s core, mostly due to enhanced dissolved inorganic carbon (DIC). The results suggest that surface variation in DIC in the eddy’s core is modulated by (i) the entrainment of CO2-rich intermediate waters at ∼500 m, (ii) low primary productivity, associated with small phytoplankton cells such as cryptophytes and green flagellates, and (iii) respiration processes caused by heterotrophic organisms (i.e., zooplankton community). By providing a comprehensive view of these physical and biogeochemical properties of this stationary eddy, our results are key to adding new insights to a better understanding of the behavior of mesoscale features influencing sea-air CO2 exchanges in polar environments.

Seasonal Superrotation in Earth’s Troposphere

Abstract Although Earth’s troposphere does not superrotate in the annual mean, for most of the year—from October to May—the winds of the tropical upper troposphere are westerly. We investigate this seasonal superrotation using reanalysis data and a single-layer model for the winds of the tropical upper troposphere. We characterize the temporal and spatial structures of the tropospheric superrotation, and quantify the relationships between the superrotation and the leading modes of tropical interannual variability. We also find that the strength of the superrotation has remained roughly constant over the past few decades, despite the winds of the tropical upper troposphere decelerating (becoming more easterly) in other months. We analyze the monthly zonal-mean zonal momentum budget and use numerical simulations with an axisymmetric, single-layer model of the tropical upper troposphere to study the underlying dynamics of the seasonal superrotation. Momentum flux convergence by stationary eddies accelerates the superrotation, while cross-equatorial easterly momentum transport associated with the Hadley circulation decelerates the superrotation. The seasonal modulations of these two competing factors shape the superrotation. The single-layer model is able to qualitatively reproduce the seasonal progression of the winds in the tropical upper troposphere, and highlights the northward displacement of the intertropical convergence zone in the annual mean as a key factor responsible for the annual cycle of the tropical winds.

Responses of Global Atmospheric Energy Transport to Idealized Groundwater Conditions in a General Circulation Model

Abstract The representation of groundwater dynamics in land surface models and their roles in global precipitation variations has received attention in recent years. Studies have revealed the overall higher soil moisture but rather diverse precipitation changes after incorporating the groundwater component in climate models. However, groundwater effects on large-scale atmospheric energy transport, the fundamental atmospheric variable regulating Earth’s climate, have not been explored thoroughly. In this study, a pair of idealized experiments corresponding to contrast globally fixed water table depths by AMIP-type simulations in the Community Earth System Model was conducted. In the wet (shallow water table) experiments, an increased meridional surface temperature gradient makes the mean meridional energy transports and Hadley circulation stronger than dry (deep water table) experiments over the tropics. Such energy transport changes are primarily attributed to the dynamic contribution (intensified Hadley circulation). The wet experiments make the simulated world be like an aquaplanet simulation with less land–sea temperature contrast and the enhancement (reduction) of mean meridional circulation (stationary eddies) energy transports. Furthermore, the South Asian monsoon circulation in the wet experiment shows a southward shift in the premonsoon season (April–June) and slight weakening in the mature phase (July and August). This study explores the impacts of the soil conditions caused by various water table depths on global energy transport and has further implications for climate model developments and experiment designs.

Partial Control of the Gulf of Mexico Dynamics by the Current Feedback to the Atmosphere

Abstract The surface oceanic current feedback (CFB) to the atmosphere has been shown to correct long-lasting biases in the representation of ocean dynamics by providing an unambiguous energy sink mechanism. However, its effects on the Gulf of Mexico (GoM) oceanic circulation are not known. Here, twin ocean–atmosphere eddy-rich coupled simulations, with and without CFB, are performed for the period 1993–2016 over the GoM to assess to which extent CFB modulates the GoM dynamics. CFB, through the eddy killing mechanism and the associated transfer of momentum from mesoscale currents to the atmosphere, damps the mesoscale activity by roughly 20% and alters eddy statistics. We furthermore show that the Loop Current (LC) extensions can be classified into three categories: a retracted LC, a canonical LC, and an elongated LC. CFB, by damping the mesoscale activity, enhance the occurrence of the elongated category (by about 7%). Finally, by increasing the LC extension, CFB plays a key role in determining LC eddy separations and statistics. Taking into account CFB improves the representation of the GoM dynamics, and it should be taken into account in ocean models.

Global coarse-grained mesoscale eddy statistics based on integrated kinetic energy and enstrophy correlations

Abstract. Recently, Jánosi et al. (2019) introduced the concept of a “vortex proxy” based on an observation of strong correlations between integrated kinetic energy and integrated enstrophy over a large enough surface area. When mesoscale vortices are assumed to exhibit a Gaussian shape, the two spatial integrals have particularly simple functional forms, and a ratio of them defines an effective radius of a “proxy vortex”. In the original work, the idea was tested over a restricted area in the Californian Current System. Here we extend the analysis to global scale by means of 25 years of AVISO altimetry data covering the (ice-free) global ocean. The results are compared with a global vortex database containing over 64 million mesoscale eddies. We demonstrate that the proxy vortex representation of surface flow fields also works globally and provides a quick and reliable way to obtain coarse-grained vortex statistics. Estimated mean eddy sizes (effective radii) are extracted in very good agreement with the data from the vortex census. Recorded eddy amplitudes are directly used to infer the kinetic energy transported by the mesoscale vortices. The ratio of total and eddy kinetic energies is somewhat higher than found in previous studies. The characteristic westward drift velocities are evaluated by a time-lagged cross-correlation analysis of the kinetic energy fields. While zonal mean drift speeds are in good agreement with vortex trajectory evaluation in the latitude bands 30–5∘ S and 5–30∘ N, discrepancies are exhibited mostly at higher latitudes on both hemispheres. A plausible reason for somewhat different drift velocities obtained by eddy tracking and cross-correlation analysis is the fact that the drift of mesoscale eddies is only one component of the surface flow fields. Rossby wave activities, coherent currents, and other propagating features on the ocean surface apparently contribute to the zonal transport of kinetic energy.

Near Invariance of Poleward Atmospheric Heat Transport in Response to Midlatitude Orography

Abstract Total poleward atmospheric heat transport (AHT) is similar in both magnitude and latitudinal structure between the Northern and Southern Hemispheres. These similarities occur despite more major mountain ranges in the Northern Hemisphere, which help create substantial stationary eddy AHT that is largely absent in the Southern Hemisphere. However, this hemispheric difference in stationary eddy AHT is compensated by hemispheric differences in other dynamic components of AHT so that total AHT is similar between hemispheres. To understand how AHT compensation occurs, we add midlatitude mountain ranges in two different general circulation models that are otherwise configured as aquaplanets. Even when midlatitude mountains are introduced, total AHT is nearly invariant. We explore the near invariance of total AHT in response to orography through dynamic, energetic, and diffusive perspectives. Dynamically, orographically induced changes to stationary eddy AHT are compensated by changes in both transient eddy and mean meridional circulation AHT. This creates an AHT system with three interconnected components that resist large changes to total AHT. Energetically, the total AHT can only change if the top-of-the-atmosphere net radiation changes at the equator-to-pole scale. Midlatitude orography does not create large-enough changes in the equator-to-pole temperature gradient to alter outgoing longwave radiation enough to substantially change total AHT. In the zonal mean, changes to absorbed shortwave radiation also often compensate for changes in outgoing longwave radiation. Diffusively, the atmosphere smooths anomalies in temperature and humidity created by the addition of midlatitude orography, such that total AHT is relatively invariant. Significance Statement The purpose of this study is to better understand how orography influences heat transport in the atmosphere. Enhancing our understanding of how atmospheric heat transport works is important, as heat transport helps moderate Earth’s surface temperatures and influences precipitation patterns. We find that the total amount of atmospheric heat transport does not change in the presence of mountains in the midlatitudes. Different pieces of the heat transport change, but they change in compensatory ways, such that the total heat transport remains roughly constant.

Revisiting the physical mechanisms of East Asian summer monsoon precipitation changes during the mid-Holocene: a data–model comparison

The mid-Holocene (MH; 6 ka) is one of the benchmark periods for the Paleoclimate Modeling Intercomparison Project (PMIP) and provides a unique opportunity to study monsoon dynamics and orbital forcing (i.e., mostly precession) that differ significantly from the present day. We conducted a data–model comparison along with a dynamic analysis to investigate monsoonal (i.e., East Asian summer monsoon; EASM) precipitation changes over East Asia during the MH. We used the three phases of the PMIP simulations for the MH, and quantitatively compare model results with pollen-based climate records. The data–model comparison shows an overall increase in the summer monsoon precipitation, except a local decrease during the MH. Decomposition of the moisture budget into thermodynamic and dynamic components allows us to assess their relative role in controlling EASM precipitation during the MH, and to investigate the precipitation changes obtained from pollen records in terms of physical processes. We show that the dynamic effect, rather than the thermodynamic effect, is the dominant control in increased EASM precipitation during the MH in both the proxy records and models. The dynamic increase in precipitation results mainly from the enhancement of horizontal monsoonal moisture transport that is caused by intensified stationary eddy horizontal circulation over East Asia. In addition, a cloud-related cooling effect reduced the thermodynamic contribution to the increase in EASM precipitation during the MH.

Interdecadal Variation in Available Potential Energy of Stationary Eddies in the Midlatitude Northern Hemisphere in Response to the North Pacific Gyre Oscillation

AbstractThis study investigates the variation in available potential energy of stationary eddies in the mid‐troposphere over the midlatitude Northern Hemisphere (30°–60°N) from 1950 to 2021, based on the ERA5 reanalysis and its preliminary back extension. The variation is attributed mainly to changes in the wavenumber‐1 eddy. It is also noted that potential energy conversion from the zonal mean flow in early winter plays an important role in the 20‐year variation. In addition, statistical analysis and numerical simulation of an ICTP model demonstrate that the potential energy conversion is controlled by the North Pacific Gyre Oscillation, which induces anomalous meridional temperature advection in the mid‐troposphere over Northeast Asia.

Modulation of Coupled Modes of Tibetan Plateau Heating and Indian Summer Monsoon on Summer Rainfall over Central Asia

Abstract It has been suggested that summer rainfall over central Asia (CA) is significantly correlated with the summer thermal distribution of the Tibetan Plateau (TP) and the Indian summer monsoon (ISM). However, relatively few studies have investigated their synergistic effects of different distribution. This study documents the significant correlations between precipitation in CA and the diabatic heating of TP and the ISM in summer based on the results of statistical analysis and numerical simulation. Summer precipitation in CA is dominated by two water vapor transport branches from the south that are related to the two primary modes of anomalous diabatic heating distribution contributed by the TP and ISM precipitation, that is, the “+−” dipole mode in the southeastern TP and the Indian subcontinent (IS), and the “+−+” tripole mode in the southeastern TP, the IS, and southern India. Both modes exhibit obvious midlatitude Silk Road pattern (SRP) wave trains with cyclone anomalies over CA, but with different transient and stationary eddies over South Asia. The different locations of anomalous anticyclones over India govern two water vapor transport branches to CA, which are from the Arabian Sea and the Bay of Bengal. The water vapor flux climbs while being transported northward and can be transported to CA with the cooperation of cyclonic circulation. The convergent water vapor and ascending motion caused by cyclonic anomalies favor the precipitation in CA. Further analysis corroborates the negative south Indian Ocean dipole in February could affect the tripole mode distribution of TP heating and ISM via the atmospheric circulation, water vapor transport, and anomalous Hadley cell circulation. The results indicate a reliable prediction reference for summer precipitation in CA.

Time-dependent response of the overturning circulation and pycnocline depth to Southern Ocean surface wind stress changes

Abstract Changes in the Southern Ocean (SO) surface wind stress influence both the meridional overturning circulation (MOC) and stratification not only in the SO but in the global oceans, which can take multiple millennia to fully equilibrate. We use a hierarchy of models to investigate the time-dependent response of the MOC and low-latitude pycnocline depth (which quantifies the stratification) to SO wind stress changes: a two-layer analytical theory, a multi-column model (PyMOC), and an idealized general circulation model (GCM). We find that in both the GCM and PyMOC, the MOC has a multi-decadal adjustment timescale while the pycnocline depth has a multi-centennial timescale. The two-layer theory instead predicts the MOC and pycnocline depth to adjust on the same, multi-decadal timescale. We argue that this discrepancy arises because the pycnocline depth depends on the bulk stratification, while the MOC amplitude is sensitive mostly to isopycnals within the overturning cell. We can reconcile the discrepancy by interpreting the “pycnocline depth” in the theory as the depth of a specific isopycnal near the maximum of the MOC. We also find that SO stationary eddies respond very quickly to a sudden wind stress change, compensating for most of the change in the Ekman-driven MOC. This effect is missing in the theory, where the eddy-induced MOC only follows the adjustment of the pycnocline depth. Our results emphasize the importance of depth-dependence in the oceans’ transient response to changes in surface boundary conditions, and the distinct role played by stationary eddies in the SO.

Bathymetric Control of Subpolar Gyres and the Overturning Circulation in the Southern Ocean

AbstractThe subpolar gyres of the Southern Ocean form an important dynamical link between the Antarctic Circumpolar Current (ACC) and the coastline of Antarctica. Despite their key involvement in the production and export of bottom water and the poleward transport of oceanic heat, these gyres are rarely acknowledged in conceptual models of the Southern Ocean circulation, which tend to focus on the zonally averaged overturning across the ACC. To isolate the effect of these gyres on the regional circulation, we carried out a set of numerical simulations with idealized representations of the Weddell Sea sector in the Southern Ocean. A key result is that the zonally oriented submarine ridge along the northern periphery of the subpolar gyre plays a fundamental role in setting the stratification and circulation across the entire region. In addition to sharpening and strengthening the horizontal circulation of the gyre, the zonal ridge establishes a strong meridional density front that separates the weakly stratified subpolar gyre from the more stratified circumpolar flow. Critically, the formation of this front shifts the latitudinal outcrop position of certain deep isopycnals such that they experience different buoyancy forcing at the surface. Additionally, the zonal ridge modifies the mechanisms by which heat is transported poleward by the ocean, favoring heat transport by transient eddies while suppressing that by stationary eddies. This study highlights the need to characterize how bathymetry at the subpolar gyre–ACC boundary may constrain the transient response of the regional circulation to changes in surface forcing.Significance StatementThis study explores the impact of seafloor bathymetry on the dynamics of subpolar gyres in the Southern Ocean. The subpolar gyres are major circulation features that connect the Antarctic Circumpolar Current (ACC) and the coastline of Antarctica. This work provides deeper insight for how the submarine ridges that exist along the northern periphery of these gyres shape the vertical distribution of tracers and overturning circulation in these regions. These findings highlight an underappreciated yet fundamentally important topographical constraint on the three-dimensional cycling of heat and carbon in the Southern Ocean—processes that have far-reaching implications for the global climate. Future work should explore how the presence of these ridges affect the time-evolving response of the Southern Ocean to changes in surface conditions.