Description: Decadal variabilities in Indian Ocean subsurface ocean heat content (OHC; 50–300 m) since the 1950s are examined using ocean reanalyses. This study elaborates on how Pacific variability modulates the Indian Ocean on decadal time scales through both oceanic and atmospheric pathways. High correlations between OHC and thermocline depth variations across the entire Indian Ocean Basin suggest that OHC variability is primarily driven by thermocline fluctuations. The spatial pattern of the leading mode of decadal Indian Ocean OHC variability closely matches the regression pattern of OHC on the interdecadal Pacific oscillation (IPO), emphasizing the role of the Pacific Ocean in determining Indian Ocean OHC decadal variability. Further analyses identify different mechanisms by which the Pacific influences the eastern and western Indian Ocean. IPO-related anomalies from the Pacific propagate mainly through oceanic pathways in the Maritime Continent to impact the eastern Indian Ocean. By contrast, in the western Indian Ocean, the IPO induces wind-driven Ekman pumping in the central Indian Ocean via the atmospheric bridge, which in turn modifies conditions in the southwestern Indian Ocean via westward-propagating Rossby waves. To confirm this, a linear Rossby wave model is forced with wind stresses and eastern boundary conditions based on reanalyses. This linear model skillfully reproduces observed sea surface height anomalies and highlights both the oceanic connection in the eastern Indian Ocean and the role of wind-driven Ekman pumping in the west. These findings are also reproduced by OGCM hindcast experiments forced by interannual atmospheric boundary conditions applied only over the Pacific and Indian Oceans, respectively.

Description: The Eddy Kinetic Energy (EKE) associated with the Subtropical Countercurrent (STCC) in the western subtropical South Pacific is known to exhibit substantial seasonal and decadal variability. Using an eddy-permitting ocean general circulation model, which is able to reproduce the observed, salient features of the seasonal cycles of shear, stratification, baroclinic production and the associated EKE, we investigate the decadal changes of EKE. We show that the STCC region exhibits, uniquely among the subtropical gyres of the world’s oceans, significant, atmospherically forced, decadal EKE variability. The decadal variations are driven by changing vertical shear between the STCC in the upper 300 m and the South Equatorial Current below, predominantly caused by variations in STCC strength associated with a changing meridional density gradient. In the 1970s, an increased meridional density gradient results in EKE twice as large as in later decades in the model. Utilizing sensitivity experiments, decadal variations in the wind field are shown to be the essential driver. Local wind stress curl anomalies associated with the Interdecadal Pacific Oscillation (IPO) lead to up- and downwelling of the thermocline, inducing strengthening or weakening of the STCC and the associated EKE. Additionally, remote wind stress curl anomalies in the eastern subtropical South Pacific, which are not related to the IPO, generate density anomalies that propagate westward as Rossby waves and can account for up to 30–40 % of the density anomalies in the investigated region.

Description: Seasonal variability of the tropical Atlantic circulation is dominated by the annual cycle, but semi-annual variability is also pronounced, despite weak forcing at that period. Here we use multi-year, full-depth velocity measurements from the central equatorial Atlantic to analyze the vertical structure of annual and semi-annual variations of zonal velocity. A baroclinic modal decomposition finds that the annual cycle is dominated by the 4th mode and the semi-annual cycle by the 2nd mode. Similar local behavior is found in a high-resolution general circulation model. This simulation reveals that the annual and semi-annual cycles of the respective dominant baroclinic modes are associated with characteristic basin-wide structures. Using an idealized linear reduced-gravity model to simulate the dynamics of individual baroclinic modes, it is shown that the observed circulation variability can be explained by resonant equatorial basin modes. Corollary simulations of the reduced-gravity model with varying basin geometry (i.e. square basin versus realistic coastlines) or forcing (i.e. spatially uniform versus spatially variable wind) show a structural robustness of the simulated basin modes. A main focus of this study is the seasonal variability of the Equatorial Undercurrent (EUC) as identified in recent observational studies. Main characteristics of the observed EUC including seasonal variability of transport, core depth, and maximum core velocity can be explained by the linear superposition of the dominant equatorial basin modes as obtained from the reduced-gravity model.

Description: Some studies of ocean climate model experiments suggest that regional changes in dynamic sea level could provide a valuable indicator of trends in the strength of the Atlantic meridional overturning circulation (MOC). This paper describes the use of a sequence of global ocean–ice model experiments to show that the diagnosed patterns of sea surface height (SSH) anomalies associated with changes in the MOC in the North Atlantic (NA) depend critically on the time scales of interest. Model hindcast simulations for 1958–2004 reproduce the observed pattern of SSH variability with extrema occurring along the Gulf Stream (GS) and in the subpolar gyre (SPG), but they also show that the pattern is primarily related to the wind-driven variability of MOC and gyre circulation on interannual time scales; it is reflected also in the leading EOF of SSH variability over the NA Ocean, as described in previous studies. The pattern, however, is not useful as a “fingerprint” of longer-term changes in the MOC: as shown with a companion experiment, a multidecadal, gradual decline in the MOC [of 5 Sv (1 Sv ≡ 106 m3 s−1) over 5 decades] induces a much broader, basin-scale SSH rise over the mid-to-high-latitude NA, with amplitudes of 20 cm. The detectability of such a trend is low along the GS since low-frequency SSH changes are effectively masked here by strong variability on shorter time scales. More favorable signal-to-noise ratios are found in the SPG and the eastern NA, where a MOC trend of 0.1 Sv yr−1 would leave a significant imprint in SSH already after about 20 years.

Description: The Indian Ocean has sustained robust surface warming in recent decades, but the role of multi-decadal variability remains unclear. Using ocean model hindcasts, characteristics of low-frequency Indian Ocean temperature variations are explored. Simulated upper-ocean temperature changes across the Indian Ocean in the hindcast are consistent with those recorded in observational products and ocean reanalyses. Indian Ocean temperatures exhibit strong warming trends since the 1950s limited to the surface and south of 30°S, while extensive subsurface cooling occurs over much of the tropical Indian Ocean. Previous work focused on diagnosing causes of these long-term trends in the Indian Ocean over the second half of the 20th Century. Instead, the temporal evolution of Indian Ocean subsurface heat content is shown here to reveal distinct multi-decadal variations associated with the Pacific Decadal Oscillation and the long-term trends are thus interpreted to result from aliasing of the low-frequency variability. Transmission of the multi-decadal signal occurs via an oceanic pathway through the Indonesian Throughflow and is manifest across the Indian Ocean centered along 12°S as westward propagating Rossby waves modulating thermocline and subsurface heat content variations. Resulting low-frequency changes in the eastern Indian Ocean thermocline depth are associated with decadal variations in the frequency of Indian Ocean Dipole (IOD) events, with positive IOD events unusually common in the 1960s and 1990s with a relatively shallow thermocline. In contrast, the deeper thermocline depth in the 1970s and 1980s is associated with frequent negative IOD and rare positive IOD events. Changes in Pacific wind forcing in recent decades and associated rapid increases in Indian Ocean subsurface heat content can thus affect the basin’s leading mode of variability, with implications for regional climate and vulnerable societies in surrounding countries.

Description: Experiments with a suite of North Atlantic general circulation models are used to examine the sources of eddy kinetic energy (EKE) in the Labrador Sea. A high-resolution model version (112°) quantitatively reproduces the observed signature. A particular feature of the EKE in the Labrador Sea is its pronounced seasonal cycle, with a maximum intensity in early winter, as already found in earlier studies based on altimeter data. In contrast to a previously advanced hypothesis, the seasonally varying eddy field is not related to a forcing by high-frequency wind variations but can be explained by a seasonally modulated instability of the West Greenland Current (WGC). The main source of EKE in the Labrador Sea is an energy transfer due to Reynolds interaction work (barotropic instability) in a confined region near Cape Desolation where the WGC adjusts to a change in the topographic slope: Geostrophic contours tend to converge upstream of Cape Desolation, such that the topographically guided WGC narrows as well and becomes barotropically unstable. The eddies spawned from the WGC instability area, dominating the EKE in the interior Labrador Sea, are predominantly anticyclonic with warm and saline cores in the upper kilometer of the water column, while the few cyclones originating as well from the instability area show a more depth-independent structure. Companion experiments with a ⅓° model exhibit the strength of the WGC, influenced by either changes in the wind stress or heat flux forcing, as a leading factor determining seasonal to interannual changes of EKE in the Labrador Sea

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: A high-resolution model of the wind-driven and thermohaline circulation in the North and equatorial Atlantic Ocean is used to study the structure and variability of the boundary current system at 26°N, including the Florida Current, the Antilles Current, and the Deep Western Boundary Current (DWBC). The model was developed by Bryan and Holland as a Community Modeling Effort of the World Ocean Circulation Experiment. Subsequent experiments have been performed at IfM Kiel, with different friction coefficients, and different climatologies of monthly mean wind stress: Hellerman–Rosenstein (HR) and Isemer–Hasse (IH). The southward volume transports in the upper 1000 m of the interior Atlantic, at 26°N, are 25.0 Sv (Sv ≡ 106m3s−1) for HR, and 34.9 Sv for IH forcing, in good agreement with the transport from the integrated Sverdrup balance at this latitude (23.9 Sv for HR, 35.6 Sv for IH). The return flow of this wind-driven transport, plus the southward transport of the DWBC (6–8 Sv), is partitioned between the Florida Current and Antilles Current. With HR forcing, the transport through the Straits of Florida is 23.2 Sv; this increases to 29.1 Sv when the wind stresses of IH are used. The annual variation of the simulated Florida Current is very similar to previous, coarse-resolution models when using the same wind-stress climatology (HR); the annual range (3.4 Sv) obtained with HR forcing is strongly enhanced (6.3 Sv) with IH forcing. The meridional heat transport at 26°N, zonally integrated across the basin, is in phase with the Florida Current; its annual range increases from 0.44 PW (HR) to 0.80 PW (IH). The annual signal east of the Bahamas is masked by strong transport fluctuations on a time scale of O(100 days), caused by an instability of the Antilles Current. By averaging over several model years, an annual cycle is extracted, which is in phase with the wind stress curl over the western part of the basin.

Description: To avoid an explicit simulation of the overflows across the Greenland-Scotland ridge, many models of the large-scale ocean circulation seek to include the net effect of the inflowing dense water masses by restoring temperature and salinity near the ridge to observed conditions. In this paper the authors examine the effect of different datasets for the northern restoring condition in two versions, eddy resolving and non-eddy resolving, of the model of the North and equatorial Atlantic that has been developed in recent years as a Community Modeling Effort for WOCE. It is shown that the use of smoothed climatological fields of temperature and salinity south of the Denmark Strait leads to strong deficiencies in the simulation of the deep flow field in the basin. A switch to actual hydrographic data from the Denmark Strait ignites a rapid dynamic response throughout the North Atlantic, affecting the transport and vertical structure of the deep western boundary current and, by virtue of the JEBAR efffect, the transport of the horizontal gyres. Meridional overturning and northward heat transport too weak in the cases with climatological boundary conditions, increase to more realistic levels in the subtropical North Atlantic. The initial response to switches in the high-latitude thermohaline forcing is mediated by fast waves along the westurn boundary, leading to changes in the deep western boundary current in low latitudes after about two years in the non-eddy-resolving cast. The initial timescale depends on the horizontal grid spacing of the model; in the high-resolution case, the first signal reaches the equator in a few months. The adjustment to a new, dynamic quasi equilibrium involves Kelvin waves along the equator and Rossby wave in the interior and is attained in less than two decades throughout the North Atlantic. It is suggested that these fast dynamic adjustment processes could play an important role in possible fluctuations of the thermohaline circulation, or transitions between different equilibrium states of the coupled ocean–atmosphere system, and may have determined the timescale of the observed climatic transitions before and during the last deglaciation.

Description: Global mean and eddy fields from a four-year experiment with a 1/6° × 1/5° horizontal resolution implementation of the CME North Atlantic model are presented. The time-averaged wind-driven and thermohaline circulation in the model is compared to the results of a 1/3° × 2/5° model run in very similar configuration. In general, the higher resolution results are found to confirm that the resolution of previous CME experiments is sufficient to describe many features of the large-scale circulation and water mass distribution quite well. While the increased resolution does not lead to large changes in the mean flow patterns, the variability in the model is enhanced significantly. On the other hand, however, not all aspects of the circulation have improved with resolution. The Azores Current Frontal Zone with its variability in the eastern basin is still represented very poorly. Particular attention is also directed toward the unrealistic stationary anticyclones north of Cape Hatteras and in the Gulf of Mexico.