Description: We report a study of a coastal frontal zone of the southeastern United States based on a field experiment and numerical modeling. The study was conducted in the spring of 1985 during weak to moderate wind stress and strong input of buoyancy from solar radiation and river discharge. The study confirms that the structure and slope of the frontal zone depends on a combination of wind stress and cross-shelf advection of buoyancy. A cross-shelf/depth two-dimensional (x, y), time-dependent numerical model illustrated the response of the frontal zone to the local wind stress regimes. A comparison of model results with field data showed that the model successfully predicted onsets of stratification and mixing. When alongshore wind stress was negative (southward), isopycnals in the frontal zone steepened due to a combination of horizontal advection and vertical convection. When stress was positive (northward), the offshore advection of low density water flattened the isopycnals and potential energy decreased, demonstrating that horizontal advection terms are important in the equation of conservation of buoyancy. The model predicts die offshore advection of lenses of less dense water during upwelling-favorable wind stress. These lenses are of the order of 20 km in cross-shelf scale and represent an efficient mechanism to export nearshore water. The lenses consist of a mixture of low-salinity coastal water and continental shelf water originating further offshore and advected onshore along the bottom. The mean flow inside the frontal zone opposed the mean alongshore wind stress. Part of the alongshore flow was in geostrophy with the cross-shore pressure gradient; the other part was due to an alongshore pressure gradient force (kinematic) of about 1 × 10−6 m s−2 (equivalent sea surface slope = 1 × 10−7), which was trapped along the coast with an offshore width scale of O(10 km). It is likely that the alongshore extent of this pressure gradient was governed by the scale at which freshwater is injected to the continental shelf, i.e., 20–30 km. The pressure gradient force immediately outside of the frontal zone was about −5 × 10−7 m s−2 in the direction of the mean alongshore wind stress. It is hypothesized that, as a result of wind setup and freshwater influx, the northward pressure gradient forced over outer shelf/slope by the Gulf Stream decreases in magnitude onshore, and can even change sign across a nearshore frontal zone of O(10 km). The implied flow field near the frontal zone is therefore highly three-dimensional with |∂v/∂y|≈|∂u/∂x|, where (u, v) are velocities in the cross-shore (x) and alongshore (y) directions, respectively.

Description: Rapid increases in upper 700‐m Indian Ocean heat content (IOHC) since the 2000s have focused attention on its role during the recent global surface warming hiatus. Here, we use ocean model simulations to assess distinct multidecadal IOHC variations since the 1960s and explore the relative contributions from wind stress and buoyancy forcing regionally and with depth. Multidecadal wind forcing counteracted IOHC increases due to buoyancy forcing from the 1960s to the 1990s. Wind and buoyancy forcing contribute positively since the mid‐2000s, accounting for the drastic IOHC change. Distinct timing and structure of upper ocean temperature changes in the eastern and western Indian Ocean are linked to the pathway how multidecadal wind forcing associated with the Interdecadal Pacific Oscillation is transmitted and affects IOHC through local and remote winds. Progressive shoaling of the equatorial thermocline—of importance for low‐frequency variations in Indian Ocean Dipole occurrence—appears to be dominated by multidecadal variations in wind forcing.

Description: We examine the diffusive behavior of the flow field in an eddy-resolving, primitive equation circulation model. Analysis of fluid particle trajectories illustrates the transport mechanisms, which are leading to uniform tracer and potential vorticity distributions in the interior of the subtropical thermocline. In contrast to the assumption of weak mixing in recent analytical theories, the numerical model indicates the alternative of tracer and potential vorticity homogenization on isopycnal surfaces taking place in a nonideal fluid with strong, along-isopycnal eddy mixing. The eastern, ventilated portion of the gyre appears to be sufficiently homogeneous to allow the concept of an eddy diffusivity to apply. A break in a random walk behavior of particle statistics occurs after about 100 days when along-flow dispersion sharply increases, indicative of mean shear effects. During the first months of particle spreading, eddy dispersal and mean advection are of similar magnitude. Eddy kinetic energy is of O(60–80 cm2 s−2) in the model thermocline, comparable to the pool of weak eddy intensity found in the eastern parts of the subtropical oceans. Eddy diffusivity in the model thermocline (Kxx = 8 × 107, Kyy = 3 × 107 cm2 s−1) seems to be higher by a factor of about 3 than oceanic values estimated for these area. Below the thermocline, model diffusivity decreases substantially and becomes much more anisotropic, with particle dispersal preferentially in the zonal direction. The strong nonisotropic behavior, prominent also in all other areas of water eddy intensity, appears as the major discrepancy when compared with the observed behavior of SOFAR floats and surface drifters in the ocean.

Description: The Makassar Strait, the main passageway of the Indonesian Throughflow (ITF), is an important component of Indo-Pacific climate through its inter-basin redistribution of heat and freshwater. Observational studies suggest that wind-driven freshwater advection from the marginal seas into the Makassar Strait modulates the strait's surface transport. However, direct observations are too short (〈15 years) to resolve variability on decadal timescales. Here we use a series of global ocean simulations to assess the advected freshwater contributions to ITF transport across a range of timescales. The simulated seasonal and interannual freshwater dynamics are consistent with previous studies. On decadal timescales, we find that wind-driven advection of South China Sea (SCS) waters into the Makassar Strait modulates upper-ocean ITF transport. Atmospheric circulation changes associated with Pacific decadal variability appear to drive this mechanism via Pacific lower-latitude western boundary current interactions that affect the SCS circulation. Key Points: - A global ocean model is used to show how freshwater impacts the decadal variability of transport through the main Indonesian Throughflow pathway - Wind-driven advection of South China Sea freshwater induces an upstream pressure gradient that reduces transport - Freshwater input is modulated by atmospheric circulation changes associated with Pacific decadal variability

Description: Deep convection and associated deep water formation are key processes for climate variability, since they impact the oceanic uptake of heat and trace gases and alter the structure and strength of the global overturning circulation. For long, deep convection in the subpolar North Atlantic was thought to be confined to the central Labrador Sea in the western subpolar gyre (SPG). However, there is increasing observational evidence that deep convection also has occurred in the eastern SPG south of Cape Farewell and in the Irminger Sea, in particular, in 2015–2018. Here we assess this recent event in the context of the temporal evolution of spatial deep convection patterns in the SPG since the mid-twentieth century, using realistic eddy-rich ocean model simulations. These reveal a large interannual variability with changing contributions of the eastern SPG to the total deep convection volume. Notably, in the late 1980s to early 1990s, the period with highest deep convection intensity in the Labrador Sea related to a persistent positive phase of the North Atlantic Oscillation, the relative contribution of the eastern SPG was small. In contrast, in 2015–2018, deep convection occurred with an unprecedented large relative contribution of the eastern SPG. This is partly linked to a smaller north-westward extent of deep convection in the Labrador Sea compared to previous periods of intensified deep convection, and may be a first fingerprint of freshening trends in the Labrador Sea potentially associated with enhanced Greenland melting and the oceanic advection of the 2012–2016 eastern North Atlantic fresh anomaly.

Description: Marine heatwaves along the coast ofWestern Australia, referred to as Ningaloo Niño, have had dramatic impacts on the ecosystem in the recent decade. A number of local and remote forcing mechanisms have been put forward, however little is known about the depth structure of such temperature extremes. Utilizing an eddy-active global Ocean General Circulation Model, Ningaloo Niño and the corresponding cold Ningaloo Niña events are investigated between 1958-2016, with focus on their depth structure. The relative roles of buoyancy and wind forcing are inferred from sensitivity experiments. Composites reveal a strong symmetry between cold and warm events in their vertical structure and associated large-scale spatial patterns. Temperature anomalies are largest at the surface, where buoyancy forcing is dominant and extend down to 300m depth (or deeper), with wind forcing being the main driver. Large-scale subsurface anomalies arise from a vertical modulation of the thermocline, extending from the western Pacific into the tropical eastern Indian Ocean. The strongest Ningaloo Niños in 2000 and 2011 are unprecedented compound events, where long-lasting high temperatures are accompanied by extreme freshening, which emerges in association with La Niñas, more common and persistent during the negative phase of the Interdecadal Pacific Oscillation. It is shown that Ningaloo Niños during La Nina phases have a distinctively deeper reach and are associated with a strengthening of the Leeuwin Current, while events during El Niño are limited to the surface layer temperatures, likely driven by local atmosphere-ocean feedbacks, without a clear imprint on salinity and velocity.

Description: Enhanced Southern Ocean ventilation in recent decades has been suggested to be a relevant modulator of the observed changes in ocean heat and carbon uptake. This study focuses on the Southern Ocean midlatitude ventilation changes from the 1960s to the 2010s. A global 1/4° configuration of the NEMO–Louvain-la-Neuve sea ice model, version 2 (LIM2), including the inert tracer CFC-12 (a proxy of ocean ventilation) is forced with the CORE, phase II (CORE-II), and JRA-55 driving ocean (JRA55-do) atmospheric reanalyses. Sensitivity experiments, where the variability of wind stress and/or the buoyancy forcing is suppressed on interannual time scales, are used to unravel the mechanisms driving ventilation changes. Ventilation changes are estimated by comparing CFC-12 interior inventories among the different experiments. All simulations suggest a multidecadal fluctuation of Southern Ocean ventilation, with a decrease until the 1980s–90s and a subsequent increase. This evolution is related to the atmospheric forcing and is caused by the (often counteracting) effects of wind stress and buoyancy forcing. Until the 1980s, increased buoyancy gains caused the ventilation decrease, whereas the subsequent ventilation increase was driven by strengthened wind stress causing deeper mixed layers and a stronger meridional overturning circulation. Wind stress emerges as the main driver of ventilation changes, even though buoyancy forcing modulates its trend and decadal variability. The three Southern Ocean basins take up CFC-12 in distinct density intervals but overall respond similarly to the atmospheric forcing. This study suggests that Southern Ocean ventilation is expected to increase as long as the effect of increasing Southern Hemisphere wind stress overwhelms that of increased stratification.