Divergence and dispersion of global eddy propagation from satellite altimetry

It is well understood that isolated eddies are presumed to propagate westward intrinsically at the speed of the annual baroclinic Rossby wave. This classic description, however, is known to be frequently violated in both propagation speed and its direction in the real ocean. Here, we present a systematic analysis on the divergence of eddy propagation direction (i.e., global pattern of departure from due west) and dispersion of eddy propagation speed (i.e., zonal pattern of departure from Rossby wave phase speed). Our main findings include the following: 1) A global climatological phase map (the first of its kind to our knowledge) indicating localized direction of most likely eddy propagation has been derived from twenty-eight years (1993-2020) of satellite altimetry, leading to a leaf-like full-angle pattern in its overall divergence. 2) A meridional deflection map of eddy motion is created with prominent equatorward/poleward deflecting zones identified, revealing that it is more geographically correlated rather than polarity determined as previously thought (i.e., poleward for cyclonic eddies and equatorward for anticyclonic ones). 3) The eddy-Rossby wave relationship has a duality nature (waves riding by eddies) in five subtropical bands centered around 27°N and 26°S in the two hemispheres, outside which their relationship has a dispersive nature with dominant waves (eddies) propagating faster in the tropical (extratropical) oceans. Current, wind and topographic effects are major external forcings responsible for the observed divergence and dispersion of eddy propagations. These results are expected to make a significant contribution to eddy trajectory prediction using physically based and/or data-driven models.

Overlapping boundary layers in coastal oceans

Boundary layer turbulence in coastal regions differs from that in deep ocean because of bottom interactions. In this paper, we focus on the merging of surface and bottom boundary layers in a finite-depth coastal ocean by numerically solving the wave-averaged equations using a large eddy simulation method. The ocean fluid is driven by combined effects of wind stress, surface wave, and a steady current in the presence of stable vertical stratification. The resulting flow consists of two overlapping boundary layers, i.e. surface and bottom boundary layers, separated by an interior stratification. The overlapping boundary layers evolve through three phases, i.e. a rapid deepening, an oscillatory equilibrium and a prompt merger, separated by two transitions. Before the merger, internal waves are observed in the stratified layer, and they are excited mainly by Langmuir turbulence in the surface boundary layer. These waves induce a clear modulation on the bottom-generated turbulence, facilitating the interaction between the surface and bottom boundary layers. After the merger, the Langmuir circulations originally confined to the surface layer are found to grow in size and extend down to the sea bottom (even though the surface waves do not feel the bottom), reminiscent of the well-organized Langmuir supercells. These full-depth Langmuir circulations promote the vertical mixing and enhance the bottom shear, leading to a significant enhancement of turbulence levels in the vertical column.

Unified Wind Wave Growth and Spectrum Functions for All Water Depths: Field Observations and Model Results

Wind wave development is governed by the fetch- or duration-limited growth principle that is expressed as a pair of similarity functions relating the dimensionless elevation variance (wave energy) and spectral peak frequency to fetch or duration. Combining the pair of similarity funtions the fetch or duration variable can be removed to form a dimensionless function of elevation variance and spectral peak frequency, which is interepreated as the wave enegry evolution with wave age. The relationship is initially developed for quasi-neural stability and quasi-steady wind forcing conditions. Further analyses show that the same fetch, duration, and wave age similarity functions are applicable to unsteady wind forcing conditions, including rapidly accelerating and decelerating mountain gap wind episodes and tropical cyclone (TC) wind fields. Here it is shown that with the dimensionless frequency converted to dimensionless wavenumber using the surface wave dispersion relationship, the same similarity function is applicable in all water depths. Field data collected in shallow to deep waters and mild to TC wind conditions, and synthetic data generated by spectrum model computations are assembled to illustrate the applicability. For the simulation work, the finite-depth wind wave spectrum model and its shoaling function are formulated for variable spectral slopes. Given wind speed, wave age, and water depth, the measrued and spectrum-computed significant wave heights and the associated growth parameters are in good agreement in forcing conditions from mild to TC winds and in all depths from deep ocean to shallow lake.

The 1994 positive Indian Ocean Dipole event as investigated by the transfer routes of oceanic wave energy

The present study investigates the interannual variability of the tropical Indian Ocean (IO) based on the transfer routes of wave energy in a set of 61-year hindcast experiments using a linear ocean model. To understand the basic feature of the IO Dipole mode, this paper focuses on the 1994 pure positive event. Two sets of westward transfer episodes in the energy flux associated with Rossby waves (RWs) are identified along the equator during 1994. One set represents the same phase speed as the linear theory of equatorial RWs, while the other set is slightly slower than the theoretical phase speed. The first set originates from the reflection of equatorial Kelvin waves at the eastern boundary of the IO. On the other hand, the second set is found to be associated with off-equatorial RWs generated by southeasterly winds in the southeastern IO, which may account for the appearance of the slower group velocity. A combined empirical orthogonal function (EOF) analysis of energy-flux streamfunction and potential reveals the intense westward signals of energy flux are attributed to off-equatorial RWs associated with predominant wind input in the southeastern IO corresponding to the positive IO Dipole event.

Indian Ocean Subtropical Underwater and the Interannual Variability in its Annual Subduction Rate Associated with the Southern Annular Mode

In this study, the Indian Ocean subtropical underwater (IOSTUW) was investigated as a subsurface salinity maximum using Argo floats (2000–2020) for the first time. It has mean salinity, potential temperature and potential density values of 35.54 ± 0.29 psu, 17.91 ± 1.66 °C, and 25.56 ± 0.35 kg m−3, respectively, and mainly extends between 10°S and 30°S along the isopycnal surface in the subtropical south Indian Ocean. The annual subduction rate of the IOSTUW during the period of 2004-2019 was investigated based on a gridded Argo dataset. The results revealed a mean value of 4.39 Sv (1 Sv=106 m3s−1) with an interannual variability that is closely related to the Southern Annular Mode (SAM). The variation in the annual subduction rate of the IOSTUW is dominated by the lateral induction term, which largely depends on the winter mixed layer depth (MLD) in the sea surface salinity (SSS) maximum region. The anomalies of winter MLD is primarily determined by SAM-related air-sea heat flux and zonal wind anomalies through modulation of the buoyancy. As a result, the annual subduction rate of the IOSTUW generally increased when the SAM index showed negative anomalies and decreased when the SAM index showed positive anomalies. Exceptional cases occurred when the wind anomaly within the SSS maximum region was weak or was dominated by its meridional component.

Export of ice sheet meltwater from Upernavik Fjord, West Greenland

Meltwater from Greenland is an important freshwater source for the North Atlantic Ocean, released into the ocean at the head of fjords in the form of runoff, submarine melt and icebergs. The meltwater release gives rise to complex in-fjord transformations that result in its dilution through mixing with other water masses. The transformed waters, which contain the meltwater, are exported from the fjords as a new water mass “Glacially Modified Water” (GMW). Here we use summer hydrographic data collected from 2013 to 2019 in Upernavik, a major glacial fjord in northwest Greenland, to describe the water masses that flow into the fjord from the shelf and the exported GMWs. Using an Optimum Multi-Parameter technique across multiple years we then show that GMW is composed of 57.8 ±8.1% Atlantic Water, 41.0 ±8.3% Polar Water, 1.0 ±0.1% subglacial discharge and 0.2 ±0.2% submarine meltwater. We show that the GMW fractional composition cannot be described by buoyant plume theory alone since it includes lateral mixing within the upper layers of the fjord not accounted for by buoyant plume dynamics. Consistent with its composition, we find that changes in GMW properties reflect changes in the AW and PW source waters. Using the obtained dilution ratios, this study suggests that the exchange across the fjord mouth during summer is on the order of 50 mSv (compared to a freshwater input of 0.5 mSv). This study provides a first order parameterization for the exchange at the mouth of glacial fjords for large-scale ocean models.

Reduction in meridional heat export contributes to recent Indian Ocean warming

Since 2000, the Indian Ocean has warmed more rapidly than the Atlantic or Pacific. Air-sea fluxes alone cannot explain the rapid Indian Ocean warming, which has so far been linked to an increase in temperature transport into the basin through the Indonesian Throughflow (ITF). Here, we investigate the role that the heat transport out of the basin at 36°S plays in the warming. Adding the heat transport out of the basin to the ITF temperature transport into the basin, we calculate the decadal mean Indian Ocean heat budget over the 2010s. We find that heat convergence increased within the Indian Ocean over 2000-2019. The heat convergence over the 2010s is the same order as the warming rate, and thus the net air-sea fluxes are near zero. This is a significant change from previous analyses using trans-basin hydrographic sections from 1987, 2002, and 2009, which all found divergences of heat. A two year time series shows that seasonal aliasing is not responsible for the decadal change. The anomalous ocean heat convergence over the 2010s compared to previous estimates is due to changes in ocean currents at both the southern boundary (33%) and the ITF (67%). We hypothesize that the changes at the southern boundary are linked to an observed broadening of the Agulhas Current, implying that temperature and velocity data at the western boundary are crucial to constrain heat budget changes.

On the Origins of the Oceanic Ultraviolet Catastrophe

We provide a first-principles analysis of the energy fluxes in the oceanic internal wavefield. The resulting formula is remarkably similar to the renowned phenomenological formula for the turbulent dissipation rate in the ocean which is known as the Finescale Parameterization. The prediction is based on the wave turbulence theory of internal gravity waves and on a new methodology devised for the computation of the associated energy fluxes. In the standard spectral representation of the wave energy density, in the two-dimensional vertical wavenumber – frequency (m – w) domain, the energy fluxes associated with the steady state are found to be directed downscale in both coordinates, closely matching the Finescale-Parameterization formula in functional form and in magnitude. These energy transfers are composed of a ‘local’ and a ‘scale-separated’ contributions; while the former is quantified numerically, the latter is dominated by the Induced Diffusion process and is amenable to analytical treatment. Contrary to previous results indicating an inverse energy cascade from high frequency to low, at odds with observations, our analysis of all non-zero coefficients of the diffusion tensor predicts a direct energy cascade. Moreover, by the same analysis fundamental spectra that had been deemed ‘no-flux’ solutions are reinstated to the status of ‘constant-downscale-flux’ solutions. This is consequential for an understanding of energy fluxes, sources and sinks that fits in the observational paradigm of the Finescale Parameterization, solving at once two long-standing paradoxes that had earned the name of ‘Oceanic Ultraviolet Catastrophe’.

Spatial and subannual variability of the Antarctic Slope Current in an eddying ocean-sea ice model

The Antarctic Slope Current (ASC) circumnavigates the Antarctic continent following the continental slope and separating the waters on the continental shelf from the deeper offshore Southern Ocean. Water mass exchanges across the continental slope are critical for the global climate as they impact the global overturning circulation and the mass balance of the Antarctic ice sheet via basal melting. Despite the ASC’s global importance, little is known about its spatial and subannual variability, as direct measurements of the velocity field are sparse. Here, we describe the ASC in a global eddying ocean-sea ice model and reveal its large-scale spatial variability by characterising the continental slope using three regimes: the surface-intensified ASC, the bottom-intensified ASC and the reversed ASC. Each ASC regime corresponds to a distinct classification of the density field as previously introduced in the literature, suggesting that the velocity and density fields are governed by the same leading-order dynamics around the Antarctic continental slope. Only the surface-intensified ASC regime has a strong seasonality. However, large temporal variability at a range of other timescales occurs across all regimes, including frequent reversals of the current. We anticipate our description of the ASC’s spatial and subannual variability to be helpful to guide future studies of the ASC aiming to advance our understanding of the region’s response to a changing climate.

Winter Convective Mixing in the northern Arabian Sea during Contrasting Monsoons

Along-track Argo observations in the northern Arabian Sea during 2017 – 19 showed by far the most contrasting winter convective mixing; 2017 – 18 was characterized by less intense convective mixing resulting in a mixed layer depth of 110 m, while 2018 – 19 experienced strong and prolonged convective mixing with the mixed layer deepening to 150 m. The response of the mixed layer to contrasting atmospheric forcing and the associated formation of Arabian Sea High Salinity Water (ASHSW) in the northeastern Arabian Sea are studied using a combination of Argo float observations, gridded observations, a data assimilative general circulation model and a series of 1-D model simulations. The 1-D model experiments show that the response of winter mixed layer to atmospheric forcing is not only influenced by winter surface buoyancy loss, but also by a preconditioned response to freshwater fluxes and associated buoyancy gain by the ocean during the summer that is preceding the following winter. A shallower and short-lived winter mixed layer occurred during 2017 – 18 following the exceptionally high precipitation over evaporation during the summer monsoon in 2017. The precipitation induced salinity stratification (a salinity anomaly of -0.7 psu) during summer inhibited convective mixing in the following winter resulting in a shallow winter mixed layer (103 m). Combined with weak buoyancy loss due to weaker surface heat loss in the northeastern Arabian Sea, this caused an early termination of the convective mixing (February 26, 2018). In contrast, the winter convective mixing during 2018 – 19 was deeper (143 m) and long-lived. The 2018 summer, by comparison, was characterized by normal or below normal precipitation which generated a weakly stratified ocean pre-conditioned to winter mixing. This combined with colder and drier air from the land mass to the north with low specific humidity lead to strong buoyancy loss, and resulted in prolonged winter convective mixing through March 25, 2019.