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: The repercussions of surface ocean currents for the near-surface wind and the air-sea momentum flux are investigated in two versions of a global climate model with eddying ocean. The focus is on the effect of mesoscale ocean current features at scales of less than 150 km, by considering high-pass filtered, monthly-mean model output fields. We find a clear signature of a mesoscale oceanic imprint in the wind fields over the energetic areas of the oceans, particularly along the extensions of the western boundary currents and the Antarctic Circumpolar Current. These areas are characterized by a positive correlation between mesoscale perturbations in the curl of the surface currents and the wind curl. The coupling coefficients are spatially non-uniform and show a pronounced seasonal cycle. The positive feedback of mesoscale current features on the near-surface wind acts in opposition to their damping effect on the wind stress. A tentative incorporation of this feedback in the surface stress formulation of an eddy-permitting global ocean-only model leads to a gain in the kinetic energy of up to 10 %, suggesting a fundamental shortcoming of present ocean model configurations.

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: We revisit the challenges and prospects for ocean circulation models following Griffies et al. (2010). Over the past decade, ocean circulation models evolved through improved understanding, numerics, spatial discretization, grid configurations, parameterizations, data assimilation, environmental monitoring, and process-level observations and modeling. Important large scale applications over the last decade are simulations of the Southern Ocean, the Meridional Overturning Circulation and its variability, and regional sea level change. Submesoscale variability is now routinely resolved in process models and permitted in a few global models, and submesoscale effects are parameterized in most global models. The scales where nonhydrostatic effects become important are beginning to be resolved in regional and process models. Coupling to sea ice, ice shelves, and high-resolution atmospheric models has stimulated new ideas and driven improvements in numerics. Observations have provided insight into turbulence and mixing around the globe and its consequences are assessed through perturbed physics models. Relatedly, parameterizations of the mixing and overturning processes in boundary layers and the ocean interior have improved. New diagnostics being used for evaluating models alongside present and novel observations are briefly referenced. The overall goal is summarizing new developments in ocean modeling, including: how new and existing observations can be used, what modeling challenges remain, and how simulations can be used to support observations.

Description: Satellite observations and output from a high-resolution ocean model are used to investigate how the Loop Current in the Gulf of Mexico affects the Gulf Stream transport through the Florida Straits. We find that the expansion (contraction) of the Loop Current leads to lower (higher) transports through the Straits of Florida. The associated surface velocity anomalies are coherent from the southwestern tip of Florida to Cape Hatteras. A simple continuity-based argument can be used to explain the link between the Loop Current and the downstream Gulf Stream transport: as the Loop Current lengthens (shortens) its path in the Gulf of Mexico, the flow out of the Gulf decreases (increases). Anomalies in the surface velocity field are first seen to the southwest of Florida and within 4 weeks propagate through the Florida Straits up to Cape Hatteras and into the Gulf Stream Extension. In both the observations and the model this propagation can be seen as pulses in the surface velocities. We estimate that the Loop Current variability can be linked to a variability of several Sverdrups (1Sv = 10(6) m(3) s(-1)) through the Florida Straits. The exact timing of the Loop Current variability is largely unpredictable beyond a few weeks and its variability is therefore likely a major contributor to the chaotic/intrinsic variability of the Gulf Stream. However, the time lag between the Loop Current and the flow downstream of the Gulf of Mexico means that if a lengthening/shortening of the Loop Current is observed this introduces some predictability in the downstream flow for a few weeks.

Description: Oceanic eddies are an important component in preconditioning the central Labrador Sea (LS) for deep convection and in restratifying the convected water. This study investigates the different sources and impacts of Eddy Kinetic Energy (EKE) and its temporal variability in the LS with the help of a 52-year long hindcast simulation of a 1/20° ocean model. Irminger Rings (IR) are generated in the West Greenland Current (WGC) between 60 and 62°N, mainly affect preconditioning and limit the northward extent of the convection area. The IR exhibit a seasonal cycle and decadal variations linked to the WGC strength, varying with the circulation of the subpolar gyre. The mean and temporal variations of IR generation can be attributed to changes in deep ocean baroclinic and upper ocean barotropic instabilities at comparable magnitudes. The main source of EKE and restratification in the central LS are Convective Eddies (CE). They are generated by baroclinic instabilities near the bottom of the mixed layer during and after convection. The CE have a mid-depth core and reflect the hydrographic properties of the convected water mass with a distinct minimum in potential vorticity. Their seasonal to decadal variability is tightly connected to the local atmospheric forcing and the associated air-sea heat fluxes. A third class of eddies in the LS are the Boundary Current Eddies shed from the Labrador Current (LC). Since they are mostly confined to the vicinity of the LC, these eddies appear to exert only minor influence on preconditioning and restratification.

Description: The Agulhas Current, the western boundary current of the South Indian Ocean, has been shown to play an important role in the connectivity between the Indian and Atlantic oceans. The greater Agulhas Current system is highly dominated by mesoscale dynamics. To investigate their influence on the regional and global circulations, a family of high-resolution ocean general circulation model configurations based on the NEMO code has been developed. Horizontal resolution refinement is achieved by embedding “nests” covering the South Atlantic and the western Indian oceans at 1/10∘ (INALT10) and 1/20∘ (INALT20) within global hosts with coarser resolutions. Nests and hosts are connected through two-way interaction, allowing the nests not only to receive boundary conditions from their respective host but also to feed back the impact of regional dynamics onto the global ocean. A double-nested configuration at 1/60∘ resolution (INALT60) has been developed to gain insights into submesoscale processes within the Agulhas Current system. Large-scale measures such as the Drake Passage transport and the strength of the Atlantic meridional overturning circulation are rather robust among the different configurations, indicating the important role of the hosts in providing a consistent embedment of the regionally refined grids into the global circulation. The dynamics of the Agulhas Current system strongly depend on the representation of mesoscale processes. Both the southward-flowing Agulhas Current and the northward-flowing Agulhas Undercurrent increase in strength with increasing resolution towards more realistic values, which suggests the importance of improving mesoscale dynamics as well as bathymetric slopes along this narrow western boundary current regime. The exploration of numerical choices such as lateral boundary conditions and details of the implementation of surface wind stress forcing demonstrates the range of solutions within any given configuration.

Description: Many models of the large-scale thermohaline circulation in the ocean exhibit strong zonally integrated upwelling in the midlatitude North Atlantic that significantly decreases the amount of deep water that is carried from the formation regions in the subpolar North Atlantic toward low latitudes and across the equator. In an analysis of results from the Community Modeling Effort using a suite of models with different horizontal resolution, wind and thermohaline forcing, and mixing parameters, it is shown that the upwelling is always concentrated in the western boundary layer between roughly 30° and 40°N. The vertical transport across 1000 m appears to be controlled by local dynamics and strongly depends on the horizontal resolution and mixing parameters of the model. It is suggested that in models with a realistic deep-water formation rate in the subpolar North Atlantic, the excessive upwelling can be considered as the prime reason for the typically too low meridional overturning rates and northward heat transports in the subtropical North Atlantic. A new isopycnal advection and mixing parameterization of tracer transports by mesoscale eddies yield substantial improvements in these integral measures of the circulation.

Description: The monthly mean wind stress climatology of Hellerman and Rosenstein (HR) is compared with the climatology of Isemer and Hasse (IH), which represents a version of the Bunker atlas (BU) for the North Atlantic based on revised parameterizations. The drag coefficients adopted by IH are 21% smaller than the values of BU and HR, and the calculation of wind speed from marine estimates of Beaufort force (Bft) is based on a revised Beaufort equivalent scale similar to the scientific scale recommended by WMO. The latter choice significantly increases wind speed below Bft 8, and effectively counteracts the reduction of the drag coefficients. Comparing the IH stresses with HR reveals substantially enhanced magnitudes in the trade wind region throughout the year. At 15°N the mean easterly stress increases from about 0.9 (HR) to about 1.2 dyn cm−1 (IH). Annual mean differences are smaller in the region of the westerlies. In winter, the effect due to the reduced drag coefficient dominates and leads to smaller stress values in IH; during summer season the revision of the Beaufort equivalents is more effective and leads to increased stresses. Implications of the different wind stress climatologies for forcing the large-scale ocean circulation are discussed by means of the Sverdrup transport streamfunction (ψs): Throughout the subtropical gyre a significant intensification of ψs takes place with IH. At 27°N, differences of more than 10 Sv (1 Sv ≡ 106 m3 s−1) are found near the western boundary. Differences in the seasonality of ψs are more pronounced in near-equatorial regions where IH increase the amplitude of the annual cycle by about 50%. An eddy-resolving model of the North Atlantic circulation is used to examine the effect of the different wind stresses on the seasonal cycle of the Florida Current. The transport predicted by the numerical model is in much better agreement with observations when the circulation is forced by IH than by HR, regarding both the annual mean (29.1 Sv vs 23.2 Sv) and the seasonal range (6.3 Sv vs 3.4 Sv).

Description: This paper describes, and establishes the dynamical mechanisms responsible for, the large-scale, time-mean, midlatitude circulation in a high-resolution model of the North Atlantic basin. The model solution is compared with recently proposed transport schemes and interpretations of the dynamical balances operating in the sub-tropical gyre. In particular, the question of the degree to which Sverdrup balance holds for the subtropical gyre is addressed. At 25°N, thermohaline-driven bottom flows cause strong local departures from the Sverdrup solution for the vertically integrated meridional mass transport, but these nearly integrate to zero across the interior of the basin. In the northwestern region of the subtropical gyre, in the vicinity of the Gulf Stream, higher-order dynamics become important, and linear vorticity dynamics is unable to explain the model's vertically integrated transport. In the subpolar gyre, the model transport bears little resemblance to the Sverdrup prediction, and higher-order dynamics are important across the entire longitudinal extent of the basin. The sensitivity of the model transport amplitudes, patterns, and dynamical balances are estimated by examining the solutions under a range of parameter choices and for four different wind stress forcing specifications. Taking into account a deficit of 7–10 Sv (Sv ≡ 106 m3 s−1) in the contribution of the model thermohaline circulation to the meridional transports at 25°N, the wind stress climatology of Isemer and Hasse appears to yield too strong of a circulation, while that derived from the NCAR Community Climate Model yields too weak of a circulation. The Hellerman and Rosenstein and ECMWF climatologies result in wind-driven transports close to observational estimates at 25°N. The range between cases for the annual mean southward transport in the interior above 1000 m is 14 Sv, which is 40%–70% of the mean transport itself. There is little sensitivity to the model closure parameters at this latitude. At 55°N, in the subpolar gyre, there is little sensitivity of the model solution to the choice of either closure parameters or wind climatology, despite large differences in the Sverdrup transports implied by the different wind stress datasets. Large year to year variability of the meridional transport east of the Bahamas makes it difficult to provide robust estimates of the sensitivity of the Antilles and deep western boundary current systems to forcing and parameter changes.