Description: Highlights: • Global mean sea level simulated in interannual CORE simulations. • Regional sea level patterns simulated in interannual CORE simulations. • Theoretical foundation for analysis of global mean sea level and regional patterns. Abstract: We provide an assessment of sea level simulated in a suite of global ocean-sea ice models using the interannual CORE atmospheric state to determine surface ocean boundary buoyancy and momentum fluxes. These CORE-II simulations are compared amongst themselves as well as to observation-based estimates. We focus on the final 15 years of the simulations (1993–2007), as this is a period where the CORE-II atmospheric state is well sampled, and it allows us to compare sea level related fields to both satellite and in situ analyses. The ensemble mean of the CORE-II simulations broadly agree with various global and regional observation-based analyses during this period, though with the global mean thermosteric sea level rise biased low relative to observation-based analyses. The simulations reveal a positive trend in dynamic sea level in the west Pacific and negative trend in the east, with this trend arising from wind shifts and regional changes in upper 700 m ocean heat content. The models also exhibit a thermosteric sea level rise in the subpolar North Atlantic associated with a transition around 1995/1996 of the North Atlantic Oscillation to its negative phase, and the advection of warm subtropical waters into the subpolar gyre. Sea level trends are predominantly associated with steric trends, with thermosteric effects generally far larger than halosteric effects, except in the Arctic and North Atlantic. There is a general anti-correlation between thermosteric and halosteric effects for much of the World Ocean, associated with density compensated changes.

Description: The upper ocean circulation of the Pacific and Indian Oceans is connected through both the Indonesian Throughflow north of Australia and the Tasman leakage around its south. The relative importance of these two pathways is examined using virtual Lagrangian particles in a high-resolution nested ocean model. The unprecedented combination of a long integration time within an eddy-permitting ocean model simulation allows the first assessment of the interannual variability of these pathways in a realistic setting. The mean Indonesian Throughflow, as diagnosed by the particles, is 14.3 Sv, considerably higher than the diagnosed average Tasman leakage of 4.2 Sv. The time series of Indonesian Throughflow agrees well with the Eulerian transport through the major Indonesian Passages, validating the Lagrangian approach using transport-tagged particles. While the Indonesian Throughflow is mainly associated with upper ocean pathways, the Tasman leakage is concentrated in the 400–900 m depth range at subtropical latitudes. Over the effective period considered (1968–1994), no apparent relationship is found between the Tasman leakage and Indonesian Throughflow. However, the Indonesian Throughflow transport correlates with ENSO. During strong La Niñas, more water of Southern Hemisphere origin flows through Makassar, Moluccas, Ombai, and Timor Straits, but less through Moluccas Strait. In general, each strait responds differently to ENSO, highlighting the complex nature of the ENSO-ITF interaction.

Description: The southwestern part of the subpolar North Atlantic east of the Grand Banks of Newfoundland and Flemish Cap is a crucial area for the Atlantic Meridional Overturning Circulation. Here the exchange between subpolar and subtropical gyre takes place, southward flowing cold and fresh water is replaced by northward flowing warm and salty water within the North Atlantic Current (NAC). As part of a long-term experiment, the circulation east of Flemish Cap has been studied by seven repeat hydrographic sections along 47 degrees N (2003-2011), a 2 year time series of current velocities at the continental slope (2009-2011), 19 years of sea surface height, and 47 years of output from an eddy resolving ocean circulation model. The structure of the flow field in the measurements and the model shows a deep reaching NAC with adjacent recirculation and two distinct cores of southward flow in the Deep Western Boundary Current (DWBC): one core above the continental slope with maximum velocities at mid-depth and the second farther east with bottom-intensified velocities. The western core of the DWBC is rather stable, while the offshore core shows high temporal variability that in the model is correlated with the NAC strength. About 30 Sv of deep water flow southward below a density of sigma=27.68 kg m(-3) in the DWBC. The NAC transports about 110 Sv northward, approximately 15 Sv originating from the DWBC, and 75 Sv recirculating locally east of the NAC, leaving 20 Sv to be supplied by the NAC from the south.

Description: Near the western boundary of the tropical North Atlantic, where the North Brazil Current (NBC) retroflects into the North Equatorial Countercurrent, large anticyclonic rings are shed. After separating from the retroflection region, the so-called NBC rings travel northwestward along the Brazilian coast, until they reach the island chain of the Lesser Antilles and disintegrate. These rings contribute substantially to the upper limb return flow of the Atlantic Meridional Overturning Circulation by carrying South Atlantic Water into the northern subtropical gyre. Their relevance for the northward transport of South Atlantic Water depends on the frequency of their generation as well as on their horizontal and vertical structure. The ring shedding and propagation and the complex interaction of the rings with the Lesser Antilles are investigated in the inline equation Family of Linked Atlantic Model Experiments (FLAME) model. The ring properties simulated in FLAME reach the upper limit of the observed rings in diameter and agree with recent observations on seasonal variability, which indicates a maximum shedding during the first half of the year. When the rings reach the shallow topography of the Lesser Antilles, they are trapped by the island triangle of St. Lucia, Barbados and Tobago and interact with the island chain. The model provides a resolution that is capable of resolving the complex topographic conditions at the islands and illuminates various possible fates for the water contained in the rings. It also reproduces laboratory experiments that indicate that both cyclones and anticyclones are formed after a ring passes through a topographic gap. Trajectories of artificial floats, which were inserted into the modeled velocity field, are used to investigate the pathways of the ring cores and their fate after they encounter the Lesser Antilles. The majority of the floats entered the Caribbean, while the northward Atlantic pathway was found to be of minor importance. No prominent pathway was found east of Barbados, where a ring could avoid the interaction with the islands and migrate toward the northern Lesser Antilles undisturbed.

Description: Model drift in the Labrador Sea in eddy permitting model simulations is examined using a series of configurations based on the NEMO numerical framework. There are two phases of the drift that we can identify, beginning with an initial rapid 3-year period, associated with the adjustment of the model from its initial conditions followed by an extended model drift/adjustment that continued for at least another decade. The drift controlled the model salinity in the Labrador Sea, over-riding the variability. Thus, during this initial period, similar behavior was observed between the inter-annually forced experiments as with perpetual year forcing. The results also did not depend on whether the configuration was global, or regional North Atlantic Ocean. The inclusion of an explicit sea-ice component did not seem to have a significant impact on the interior drift. Clear cut evidence for the drift having an advective nature was shown, based on two separate currents/flow regimes. We find, as expected, the representation of freshwater in the sub-polar gyre’s boundary currents important. But this study also points out another, equally important process and pathway: the input of high salinity mode water from the subtropical North Atlantic. The advective regime is dependent on the details of the model, such as the representation of the freshwater transport in the model’s East Greenland Current being very sensitive to the strength of the local sea surface salinity restoring (and the underlying field that the model is being restored to).

Description: Representing correctly the distribution of biogeochemical tracers in the interior ocean, such as oxygen or phosphate, is hampered by large biases in the representation of circulation in the coarse resolution models. Here we assess the oxygen and phosphate budget in two configurations of a coupled circulation biogeochemical model (NEMO - NPZD), focusing on the Atlantic Ocean. These two configurations have been integrated using realistic atmospheric forcings for the period 1948-2007. While a coarse (0.5°) configuration displays the common bias of too low oxygen associated with too high phosphate concentration, particularly at intermediate depth in the eastern side of the basin, the values are closer to the observations in an eddying (0.1°) configuration. The improvement in the representation of oxygen and phosphate is traced to a stronger transport by a more realistic representation of the equatorial and off-equatorial undercurrents. The biogeochemical fluxes are less sensitive to the current strength as the phytoplankton growth is mainly limited by the available light in the two configurations. This study emphasizes the need of high resolution models to tackle coupled biogeochemical problematics, such as the extension of oxygen minimum zones or variability in the eastern boundary upwelling system productivity.

Description: A sequence of global ocean circulation models, with horizontal mesh sizes of 0.5°, 0.25° and 0.1°, are used to estimate the long-term dispersion by ocean currents and mesoscale eddies of a slowly decaying tracer (half-life of 30 years, comparable to that of 137Cs) from the local waters off the Fukushima Dai-ichi Nuclear Power Plants. The tracer was continuously injected into the coastal waters over some weeks; its subsequent spreading and dilution in the Pacific Ocean was then simulated for 10 years. The simulations do not include any data assimilation, and thus, do not account for the actual state of the local ocean currents during the release of highly contaminated water from the damaged plants in March–April 2011. An ensemble differing in initial current distributions illustrates their importance for the tracer patterns evolving during the first months, but suggests a minor relevance for the large-scale tracer distributions after 2–3 years. By then the tracer cloud has penetrated to depths of more than 400 m, spanning the western and central North Pacific between 25°N and 55°N, leading to a rapid dilution of concentrations. The rate of dilution declines in the following years, while the main tracer patch propagates eastward across the Pacific Ocean, reaching the coastal waters of North America after about 5–6 years. Tentatively assuming a value of 10 PBq for the net 137Cs input during the first weeks after the Fukushima incident, the simulation suggests a rapid dilution of peak radioactivity values to about 10 Bq m−3 during the first two years, followed by a gradual decline to 1–2 Bq m−3 over the next 4–7 years. The total peak radioactivity levels would then still be about twice the pre-Fukushima values.

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: Previous studies found a negative trend in oxygen concentrations in tropical regions during the last decades. Employing a biogeochemical circulation model, we highlight the importance of wind driven ocean transport associated with the Subtropical-Tropical Cells (STCs) in setting the oxygen levels in the tropical ocean. The observed and simulated slowdown of the STCs by 30 percent from the 1960s to the 1990s caused a decrease in oxygen transport to the tropics. Transport of phosphate was similarly reduced, decreasing export production and respiration. The effects of physical transport and biological consumption partly compensate, damping oxygen interannual and decadal variability. Our results suggest that the observed residual oxygen trend in the tropical Pacific is mainly driven by changes in oxygen transport. Accordingly, the observed recent strengthening of the STCs leads us to expect a pause in the oxygen decrease or even an increase of tropical Pacific oxygen values in the near future.