Western Boundary Current Extensions Research Articles

Low-frequency variability enhancement of the midlatitude climate in an eddy-resolving coupled ocean–atmosphere model—part II: ocean mechanisms

This paper investigates the spatial inhomogeneity of the time-averaged, quasigeostrophic, double-gyre circulation response to fixed, realistic, large-scale modes of wind-stress forcing. While the companion paper of this study focused on understanding the anatomy of low-frequency, midlatitude climate variability in an idealised, eddy-resolving coupled model, this paper looked at understanding the nature of the wind-induced ocean gyre response using an ocean-only configuration of the same model. Our analysis revealed two, time-averaged responses to an east–west dipole, wind-stress curl anomaly in the ocean basin. Firstly, wind-stress anomalies in the western ocean basin led to changes in relative strength of the inertial recirculation zones and jet-axis tilt. This is consistent with an advection-dominated, nonlinear adjustment of the ocean gyres to anomalous forcing. Secondly, wind-stress curl anomalies in the eastern ocean basin was found to induce a largely independent response involving meridional shifts of the western boundary current extension (WBCE). The effects of time-averaged advection in this region are weak and the discovery of westward-propagating Rossby waves along the WBCE revealed the response is more akin to a baroclinic Rossby wave adjustment.

Low-frequency variability enhancement of the midlatitude climate in an eddy-resolving coupled ocean–atmosphere model. Part I: anatomy

This study investigated the coupling of the wind-driven ocean gyres with the atmospheric westerly jet using an idealised, eddy-resolving, coupled model. An empirical orthogonal function analysis of the low-pass filtered data showed that the ocean gyre variability is dominated by meridional shifts of the western boundary current extension (WBCE) and changes in the strength of the subtropical inertial recirculation zone. On the other hand, the atmospheric potential vorticity (PV) variability is dominated by the growth of standing Rossby wave patterns, while its pressure variability is dominated by a zonally-asymmetric meridional shift of the atmospheric jet. Damping sea surface temperature (SST) variability in the atmosphere was shown to weaken its PV variability and reduce the zonal asymmetry of the jet-shift mode. Singular value decompositions revealed a positive feedback between meridional shifts of the WBCE and the growth of standing Rossby wave disturbances in the atmospheric jet. The atmosphere’s response is controlled by shifts in the meridional eddy heat flux over the SST front which triggers the growth of baroclinic instabilities. This instability growth eventually leads to a large-scale, barotropic pressure response over the eastern ocean basin, or an aforementioned meridional shift of the atmospheric jet. Reduction in the atmospheric resolution inhibits the ability of atmospheric eddies to resolve length scales associated with meridional shifts of the SST front and WBCE. The lack of resolution consequently weakens the influence of ocean gyre variability on the atmospheric jet and reduces the strength of the positive feedback.

Western boundary circulation and coastal sea-level variability in Northern Hemisphere oceans

Abstract. The northwest basins of the Atlantic and Pacific oceans are regions of intense western boundary currents (WBCs): the Gulf Stream and the Kuroshio. The variability of these poleward currents and their extensions in the open ocean is of major importance to the climate system. It is largely dominated by in-phase meridional shifts downstream of the points at which they separate from the coast. Tide gauges on the adjacent coastlines have measured the inshore sea level for many decades and provide a unique window on the past of the oceanic circulation. The relationship between coastal sea level and the variability of the western boundary currents has been previously studied in each basin separately, but comparison between the two basins is missing. Here we show for each basin that the inshore sea level upstream of the separation points is in sustained agreement with the meridional shifts of the western boundary current extension over the period studied, i.e. the past 7 (5) decades in the Atlantic (Pacific). Decomposition of the coastal sea level into principal components allows us to discriminate this variability in the upstream sea level from other sources of variability such as the influence of large meanders in the Pacific. Our result extends previous findings limited to the altimetry era and suggests that prediction of inshore sea-level changes could be improved by the inclusion of meridional shifts of the western boundary current extensions as predictors. Long-duration tide gauges, such as Key West, Fernandina Beach or Hosojima, could be used as proxies for the past meridional shifts of the western boundary current extensions.

Oceanic Circulation Drives the Deepest and Longest Marine Heatwaves in the East Australian Current System

AbstractAlthough the impacts of marine heatwaves (MHWs) can extend well below the ocean surface, little is known about how oceanic and atmospheric forcings control their vertical structure. Here, we relate the MHW drivers to their sub‐surface characteristics in different dynamical regimes including the East Australia Current. We detect MHWs in the depth‐dependent surface mixed layer over 30‐years and use a heat budget to identify the dominant mechanisms driving them. We show that MHWs in the Western Boundary Current (WBC) jet are predominantly driven by air‐sea heatflux whilst in the WBC extension, MHWs are advection‐driven. The deepest and longest MHWs are advection‐driven and are more prevalent in autumn and winter. Surface (latent) flux‐driven MHWs are shallower and occur predominantly in summer. Demonstrating how MHW characteristics are linked to their drivers facilitates their prediction through driver diagnosis, especially below the surface where observations are sparse and ecological impact is high.

Trends in the Agulhas Return Current

Recent observations report a poleward shift in several of the world’s western boundary current extensions in response to a warming climate. However, the nature and variability of the Agulhas Return Current (ARC) is still relatively unknown, leading to uncertainty in its response to gyre intensification. The objective of this work is to investigate meridional migration and eddy kinetic energy (EKE) trends in the ARC over the past 27 years using satellite-derived altimetry products. A difference is seen between the western and eastern regions of the ARC, the former exhibiting a predominantly equatorward migration and the latter moving predominantly poleward. The ARC is shown to be a region highly influenced by its bathymetry. Additionally, we show that the regional EKE peaks during austral summer months and that there has been a significant rise in summertime EKE levels between 1994–2011.

On the Upper-Ocean Vertical Eddy Heat Transport in the Kuroshio Extension. Part I: Variability and Dynamics

AbstractOceanic eddies play a crucial role in transporting heat from the subsurface to surface ocean. However, dynamics responsible for the vertical eddy heat transport QT have not been systematically understood, especially in the mixed layer of western boundary current extensions characterized by the coincidence of strong eddy activities and air–sea interactions. In this paper, the winter (December–March) QT in the Kuroshio Extension is simulated using a 1-km regional ocean model. An omega equation based on the geostrophic momentum approximation and generalized to include the viscous and diabatic effects is derived and used to decompose the contribution of QT from different dynamics. The simulated QT exhibits a pronounced positive peak around the center of the mixed layer (~60 m). The value of QT there exhibits multi-time-scale variations with irregularly occurring extreme events superimposed on a slowly varying seasonal cycle. The proposed omega equation shows good skills in reproducing QT, capturing its spatial and temporal variations. Geostrophic deformation and vertical mixing of momentum are found to be the two major processes generating QT in the mixed layer with the former and the latter accounting for its seasonal variation and extreme events, respectively. The mixed layer instability and the net effect of frontogenesis/frontolysis contribute comparably to the geostrophic deformation induced QT. The contribution of QT from vertical mixing of momentum can be understood on the basis of turbulent thermal wind balance.

Maintenance of mid-latitude oceanic fronts by mesoscale eddies.

Oceanic fronts associated with strong western boundary current extensions vent a vast amount of heat into the atmosphere, anchoring mid-latitude storm tracks and facilitating ocean carbon sequestration. However, it remains unclear how the surface heat reservoir is replenished by ocean processes to sustain the atmospheric heat uptake. Using high-resolution climate simulations, we find that the vertical heat transport by ocean mesoscale eddies acts as an important heat supplier to the surface ocean in frontal regions. This vertical eddy heat transport is not accounted for by the prevailing inviscid and adiabatic ocean dynamical theories such as baroclinic instability and frontogenesis but is tightly related to the atmospheric forcing. Strong surface cooling associated with intense winds in winter promotes turbulent mixing in the mixed layer, destructing the vertical shear of mesoscale eddies. The restoring of vertical shear induces an ageostrophic secondary circulation transporting heat from the subsurface to surface ocean.

Satellite-Driven Estimates of Water Mass Formation and Their Spatio-Temporal Evolution

4 We derive water mass transformation and formation rates using satellite-derived datasets of 5 salinity, temperature and fluxes of heat and freshwater over the North Atlantic, North Pacific 6 and Southern Ocean. The formation rates are expressed in three coordinate systems: (1) density, 7 (2) temperature-salinity and (3) latitude-longitude. In the North Atlantic and North Pacific, peak 8 formation occurs south of the western boundary current extensions during the winter months of 9 the study period. In the Southern Ocean, wintertime peak formation occurs just north of the sub- 10 Antarctic Front. The satellite-derived water mass properties and formation areas agree well with 11 previous estimates from literature. The location of peak Mode Water formation varies slightly with 12 time in all coordinate systems. We assess seasonal and inter-annual variability in all three basins 13 from 2012 to 2014. We assess the impact of satellite uncertainties on final estimates of formation 14 rates and areas with Monte-Carlo simulations. The simulations provide insights on the associated 15 uncertainty of formation estimates. They also provide information on the geographic spread of the 16 water mass formation area subject to the satellite errors. We find that the total uncertainty is 17 dominated by the uncertainty in the sea surface salinity dataset. This stresses the need for frequent 18 and increasingly accurate sea surface salinity data for reliable estimates of water mass formation 19 rates and areas. Our study highlights the feasibility of providing satellite-based estimates of water 20 mass formation rates and areas. The good spatio-temporal coverage of satellite data further adds 21 to the utility of the approach.

Inhomogeneous Turbulent Dispersion across the Nearshore Induced by Surfzone Eddies

AbstractIn various oceanic regions, drifter-derived diffusivities reach a temporal maximum and subsequently decrease. Often, these are regions of inhomogeneous eddies, however, the effect of inhomogeneous turbulence on dispersion is poorly understood. The nearshore region (spanning from the surfzone to the inner shelf) also has strong cross-shore inhomogeneous turbulence. Nearshore Lagrangian statistics are estimated from drifter trajectories simulated with a wave-resolving two-dimensional Boussinesq model with random, normally incident, and directionally spread waves. The simulation is idealized and does not include other (wind, tidal, Coriolis) processes. The eddy field cross-shore inhomogeneity affects both the mean cross-shore drift and cross- and alongshore diffusivities. Short-time diffusivities are locally ballistic, and the mean drift is toward the eddy velocity variance maximum. The diffusivities reach a maximum and subsequently decrease, that is, are subdiffusive. The diffusivity maximum and time to maximum are parameterized taking into account the eddy field inhomogeneity. At long times, the cross- and alongshore diffusivities scale as t−1/2 and t−1/4, respectively, which is related to the offshore decay of the eddy intensity. A diffusion equation, with a space-dependent Fickian diffusivity related to the eddy velocity variance, reproduced the short-, intermediate-, and long-time behavior of the mean drift and cross-shore diffusivity. The small Middleton parameter, indicating fixed float dispersion, suggests the Eulerian time scale can parameterize the Lagrangian time scale in this region. Although this idealized simulation had no mean currents, and thus no shear dispersion or mixing suppression, inhomogeneous turbulence effects may be relevant in other regions such as the Antarctic Circumpolar Current (ACC) and western boundary current extensions.

Subtropical Mode Water and Permanent Pycnocline Properties in the World Ocean

AbstractA global reference state of the subtropical mode water and permanent pycnocline properties for the 2000–2015 period is presented. The climatology is obtained from a pattern recognition algorithm applied to stratification profiles from the Argo global array. The stratification features are identified as permanent upper ocean pycnostad and pycnocline even when the seasonal pycnocline is developed. The climatology shows that both Northern Hemisphere subtropical gyres have a qualitatively very similar stratification structure. The permanent pycnoclines in the North Atlantic and North Pacific show two deep centers colocated with thick subtropical and subpolar mode waters. These centers coincide with modes in the density and stratification space. These deep pycnocline centers are separated by a region with a shallower and thinner permanent pycnocline that is located downstream of Western Boundary Current Extensions and upstream of Eastern Subtropical Fronts. This feature creates a remarkable double‐bowl pattern at the basin scale. In the subtropical gyres of the South Atlantic and South Pacific Oceans, the mode water and permanent pycnocline structures are characterized by two modes in the density and stratification space that, unlike in the Northern Hemisphere, do not necessarily correspond to deep and thick permanent pycnocline regions. In the subtropical gyre of the South Indian Ocean, a single mode is found to correspond to a single center in the western part of the gyre. This study also shows that away from these deep centers where the pycnocline depth almost follows isopycnals, the permanent pycnocline experiences significant thermohaline gradients that are not density compensated.

Enhanced warming of the subtropical mode water in the North Pacific and North Atlantic

Warming of surface ocean waters is well known, but how the subsurface waters are changing is less clear. This study shows that subtropical mode water in the North Atlantic and North Pacific is warming at twice the rate of the surface waters. Over the past six decades, the subtropical surface ocean has warmed at rates close to those of global mean surface ocean temperature1 except in western boundary current regions where the surface warming is locally enhanced by a factor of two2. Changes in the subsurface ocean, however, remain unclear because of lack of data. Compiling historical temperature measurements—some available for the first time—here we show that the subtropical mode water has warmed over the past six decades in both the North Pacific and North Atlantic. The rate of the warming is twice as large in the mode waters than at the surface. Subtropical mode waters are important water masses of vertically uniform temperature that are a few hundred metres thick and distributed widely in the main thermocline of the subtropical oceans3. The enhanced warming of subtropical mode waters can be traced back to the surface warming in the formation regions along the western boundary current extensions. Furthermore, we detect increased temperature stratification and decreased dissolved oxygen in the subtropical mode waters. The latter change has clear implications for predicting biogeochemical responses to climate warming.

Decadal Variability of Upper-Ocean Heat Content Associated with Meridional Shifts of Western Boundary Current Extensions in the North Pacific

Generation and propagation processes of upper-ocean heat content (OHC) in the North Pacific are investigated using oceanic subsurface observations and an eddy-resolving ocean general circulation model hindcast simulation. OHC anomalies are decomposed into physically distinct dynamical components (OHC ρ) due to temperature anomalies that are associated with density anomalies and spiciness components (OHC χ) due to temperature anomalies that are density compensating with salinity. Analysis of the observational and model data consistently shows that both dynamical and spiciness components contribute to interannual–decadal OHC variability, with the former (latter) component dominating in the subtropical (subpolar) North Pacific. OHC ρ variability represents heaving of thermocline, propagates westward, and intensifies along the Kuroshio Extension, consistent with jet-trapped Rossby waves, while OHC χ variability propagates eastward along the subarctic frontal zone, suggesting advection by mean eastward currents. OHC χ variability tightly corresponds in space to horizontal mean spiciness gradients. Meanwhile, area-averaged OHC χ anomalies in the western subarctic frontal zone closely correspond in time to meridional shifts of the subarctic frontal zone. Regression coefficient of the OHC χ time series on the frontal displacement anomalies quantitatively agree with the area-averaged mean spiciness gradient in the region, and suggest that OHC χ is generated via frontal variability in the subarctic frontal zone.

A global probabilistic study of the ocean heat content low‐frequency variability: Atmospheric forcing versus oceanic chaos

AbstractA global 1/4° ocean/sea ice 50‐member ensemble simulation is used to disentangle the low‐frequency imprints of the atmospherically forced oceanic variability and of the chaotic intrinsic oceanic variability (IOV) on the large‐scale (10° × 10°) ocean heat content (OHC) between 1980 and 2010. The IOV explains most of the interannual‐to‐decadal large‐scale OHC variance over substantial fractions of the global ocean area that increase with depth: 9%, 22%, and 31% in the 0–700 m, 700–2000 m and 2000 m bottom layers, respectively. Such areas concern principally eddy‐active regions, mostly found in the Southern Ocean and in western boundary current extensions, and also concern the subtropical gyres at intermediate and deep levels. The oceanic chaos may also induce random multidecadal fluctuations so that large‐scale regional OHC trends computed on the 1980–2010 period cannot be unambiguously attributed to the atmospheric forcing in several oceanic basins at various depths. These results are likely to raise detection and attribution issues from real observations.

Eddy stirring and horizontal diffusivity from Argo float observations: Geographic and depth variability

AbstractStirring along isopycnals is a significant factor in determining the distribution of tracers within the ocean. Salinity anomalies on density surfaces from Argo float profiles are used to investigate horizontal stirring and estimate eddy mixing lengths. Eddy mixing length and velocity fluctuations from the ECCO2 global state estimate are used to estimate horizontal diffusivity at a 300 km scale in the upper 2000 m with near‐global coverage. Diffusivity varies by over two orders of magnitude with latitude, longitude, and depth. In all basins, diffusivity is elevated in zonal bands corresponding to strong current regions, including western boundary current extension regions, the Antarctic Circumpolar Current, and equatorial current systems. The estimated mixing lengths and diffusivities provide an observationally based data set that can be used to test and constrain predictions and parameterizations of eddy stirring.

Ocean currents generate large footprints in marine palaeoclimate proxies.

Fossils of marine microorganisms such as planktic foraminifera are among the cornerstones of palaeoclimatological studies. It is often assumed that the proxies derived from their shells represent ocean conditions above the location where they were deposited. Planktic foraminifera, however, are carried by ocean currents and, depending on the life traits of the species, potentially incorporate distant ocean conditions. Here we use high-resolution ocean models to assess the footprint of planktic foraminifera and validate our method with proxy analyses from two locations. Results show that foraminifera, and thus recorded palaeoclimatic conditions, may originate from areas up to several thousands of kilometres away, reflecting an ocean state significantly different from the core site. In the eastern equatorial regions and the western boundary current extensions, the offset may reach 1.5 °C for species living for a month and 3.0 °C for longer-living species. Oceanic transport hence appears to be a crucial aspect in the interpretation of proxy signals.

Focusing and vertical mode scattering of the first mode internal tide by mesoscale eddy interaction

Numerical experiments are performed using the MITgcm to investigate the interaction of a mode‐one internal tide with barotropic and baroclinic mode‐one mesoscale eddies. Results show that after a mode‐one internal tide passes through a barotropic eddy, spatial hot and cold spots of energy flux are produced in beam‐like patterns. The magnitude of the energy flux in the hot spots can exceed twice the incident flux while in the cold spots can reach nearly zero. Passing a mode‐one internal tide through a mode‐one baroclinic eddy results in the scattering of energy from the incident mode‐one to modes two and higher. The higher mode waves are produced in beam‐like patterns. For the parameter regime explored here, we find conversion efficiencies that reach 13% for eddies of diameter 120 km. The Rossby numbers for our experiments are order one, corresponding to energetic mesoscale eddies that are typically found in western boundary current extensions and in the southern ocean. These eddies have length scales comparable to those of low‐mode internal tides, and we expect that interaction between the two will be easily formed in locations where these phenomena coexist.

Ocean bottom pressure seasonal cycles and decadal trends from GRACE Release‐05: Ocean circulation implications

Ocean mass variations are important for diagnosing sea level budgets, the hydrological cycle, the global energy budget, and ocean circulation variability. Here seasonal cycles and decadal trends of ocean mass from January 2003 to December 2012, both global and regional, are analyzed using GRACE Release‐05 data. The trend of global flux of mass into the ocean approaches 2 cm decade−1 in equivalent sea level rise. Regional trends are of similar magnitude, with the North Pacific, South Atlantic, and South Indian oceans generally gaining mass and other regions losing mass. These trends suggest a spin‐down of the North Pacific western boundary current extension and the Antarctic Circumpolar Current in the South Atlantic and South Indian oceans. The global average seasonal cycle of ocean mass is about 1 cm in amplitude, with a maximum in early October and volume fluxes in and out of the ocean reaching 0.5 Sv (1 Sv = 1 × 106 m3 s−1) when integrated over the area analyzed here. Regional patterns of seasonal ocean mass change have typical amplitudes of 1–4 cm, and include maxima in the subtropics and minima in the subpolar regions in hemispheric winters. The subtropical mass gains and subpolar mass losses in the winter spin‐up both subtropical and subpolar gyres, hence the western boundary current extensions. Seasonal variations in these currents are order 10 Sv, but since the associated depth‐averaged current variations are only order 0.1 cm s−1, they would be difficult to detect using in situ oceanographic instruments.

Eddy Shape, Orientation, Propagation, and Mean Flow Feedback in Western Boundary Current Jets

AbstractThis manuscript revisits a study of eddy–mean flow interactions in an idealized model of a western boundary current extension jet using properties of the horizontal velocity correlation tensor to diagnose characteristics of average eddy shape, orientation, propagation, and mean flow feedback. These eddy characteristics are then used to provide a new description of the eddy–mean flow interactions observed in terms of different ingredients of the eddy motion. The diagnostics show patterns in average eddy shape, orientation, and propagation that are consistent with the signatures of jet instability in the upstream region and wave radiation in the downstream region. Together they give a feedback onto the mean flow that gives the downstream character of the jet and drives the jet's recirculation gyres. A breakdown of the eddy forcing into contributions from individual terms confirms the expected role of cross-jet gradients in meridional eddy tilt in stabilizing the jet to its barotropic instability; however, it also reveals important roles played by the along-jet evolution of eddy zonal–meridional elongation. It is the mean flow forcing derived from these patterns that acts to strengthen and extend the jet downstream and forces the time-mean recirculation gyres. This understanding of the dependence of mean flow forcing on eddy structural properties suggests that failure to adequately resolve eddy elongation could underlie the weakened jet strength, extent, and changed recirculation structure seen in this idealized model for reduced spatial resolutions. Further, it may suggest new ideas for the parameterization of this forcing.

Response of the Inertial Recirculation to Intensified Stratification in a Two-Layer Quasigeostrophic Ocean Circulation Model

Abstract Previous observation and model studies show that the upper-ocean stratification is enhanced under global warming (Capotondi et al.; Cravatte et al.; Deser et al., etc.). The response of the recirculation, which is associated with the western boundary current (WBC) jet extension and significantly increases its transport, to the intensified stratification, is studied in a two-layer quasigeostrophic ocean circulation model. It is found that the barotropic transport of the circulation first increases with stratification but then decreases as a result of saturation of the surface-layer circulation intensity when the stratification exceeds a threshold. PV budget analysis indicates that the saturation is caused by the increased intergyre transport of relative potential vorticity resulting from the intensified variability of the jet location. Both the barotropic instability and bifurcation mechanisms contribute to the intensified variability of the jet location. Because of barotropic instability, eddies are generated in the confluence region of the WBCs and advected eastward, causing the variability of the jet location. With increased stratification, the surface-layer circulation is strengthened and the barotropic instability is intensified. As a result, the surface flow becomes more variable with excessive eddies and intense variability of the jet. With the increasing stratification, three regimes, each marked by its own variation of the jet location, emerge owing to the successive system bifurcations. In the last two regimes, variability of the jet location is further enhanced by frequent switches among the different dynamic states on multidecadal time scales.

Sea Level Expression of Intrinsic and Forced Ocean Variabilities at Interannual Time Scales

Abstract This paper evaluates in a realistic context the local contributions of direct atmospheric forcing and intrinsic oceanic processes on interannual sea level anomalies (SLAs). A ¼° global ocean–sea ice general circulation model, driven over 47 yr by the full range of atmospheric time scales, is quantitatively assessed against altimetry and shown to reproduce most observed features of the interannual SLA variability from 1993 to 2004. Comparing this simulation with a second driven only by the climatological annual cycle reveals that the intrinsic part of the total interannual SLA variance exceeds 40% over half of the open-ocean area and exceeds 80% over one-fifth of it. This intrinsic contribution is particularly strong in eddy-active regions (more than 70%–80% in the Southern Ocean and western boundary current extensions) as predicted by idealized studies, as well as within the 20°–35° latitude bands. The atmosphere directly forces most of the interannual SLA variance at low latitudes and in most midlatitude eastern basins, in particular north of about 40°N in the Pacific. The interannual SLA variance is almost entirely due to intrinsic processes south of the Antarctic Circumpolar Current in the Indian Ocean sector, while half of this variance is forced by the atmosphere north of it. The same simulations were performed and analyzed at 2° resolution as well: switching to this laminar regime yields a comparable forced variability (large-scale distribution and magnitude) but almost suppresses the intrinsic variability. This likely explains why laminar ocean models largely underestimate the interannual SLA variance.