Description: A set of experiments utilizing different implementations of the global ORCA-LIM model with horizontal resolutions of 2°, 0.5° and 0.25° is used to investigate tropical and extra-tropical influences on equatorial Pacific SST variability at interannual to decadal time scales. The model experiments use a bulk forcing methodology building on the global forcing data set for 1958 to 2000 developed by Large and Yeager (2004) that is based on a blend of atmospheric reanalysis data and satellite products. Whereas representation of the mean structure and transports of the (sub-) tropical Pacific current fields is much improved with the enhanced horizontal resolution, there is only little difference in the simulation of the interannual variability in the equatorial regime between the 0.5° and 0.25° model versions, with both solutions capturing the observed SST variability in the Niño3-region. The question of remotely forced oceanic contributions to the equatorial variability, in particular, the role of low-frequency changes in the transports of the Subtropical Cells (STCs), is addressed by a sequence of perturbation experiments using different combinations of fluxes. The solutions show the near-surface temperature variability to be governed by wind-driven changes in the Equatorial Undercurrent. The relative contributions of equatorial and off-equatorial atmospheric forcing differ between interannual and longer, (multi-) decadal timescales: for the latter there is a significant impact of changes in the equatorward transport of subtropical thermocline water associated with the lower branches of the STCs, related to variations in the off-equatorial trade winds. A conspicuous feature of the STC variability is that the equatorward transports in the interior and along the western boundary partially compensate each other at both decadal and interannual time scales, with the strongest transport extrema occurring during El Niño episodes. The behaviour is rationalized in terms of a wobbling in the poleward extents of the tropical gyres, which is manifested also in a meridional shifting of the bifurcation latitudes of the North and South Equatorial Current systems.

Description: The relation between the Agulhas Current retroflection location and the magnitude of Agulhas leakage, the transport of water from the Indian to the Atlantic Ocean, is investigated in a high-resolution numerical ocean model. Sudden eastward retreats of the Agulhas Current retroflection loop are linearly related to the shedding of Agulhas rings, where larger retreats generate larger rings. Using numerical Lagrangian floats a 37 year time series of the magnitude of Agulhas leakage in the model is constructed. The time series exhibits large amounts of variability, both on weekly and annual time scales. A linear relation is found between the magnitude of Agulhas leakage and the location of the Agulhas Current retroflection, both binned to three month averages. In the relation, a more westward location of the Agulhas Current retroflection corresponds to an increased transport from the Indian Ocean to the Atlantic Ocean. When this relation is used in a linear regression and applied to almost 20 years of altimetry data, it yields a best estimate of the mean magnitude of Agulhas leakage of 13.2 Sv. The early retroflection of 2000, when Agulhas leakage was probably halved, can be identified using the regression.

Description: The skill of numerical Lagrangian drifter trajectories in three numerical models is assessed by comparing these numerically obtained paths to the trajectories of drifting buoys in the real ocean. The skill assessment is performed using the two-sample Kolmogorov-Smirnov statistical test. To demonstrate the assessment procedure, it is applied to three different models of the Agulhas region. The test can either be performed using crossing positions of one-dimensional sections in order to test model performance in specific locations, or using the total two-dimensional data set of trajectories. The test yields four quantities: a binary decision of model skill, a confidence level which can be used as a measure of goodness-of-fit of the model, a test statistic which can be used to determine the sensitivity of the confidence level, and cumulative distribution functions that aid in the qualitative analysis. The ordering of models by their confidence levels is the same as the ordering based on the qualitative analysis, which suggests that the method is suited for model validation. Only one of the three models, a 1/10 degree two-way nested regional ocean model, might have skill in the Agulhas region. The other two models, a 1/2 degree global model and a 1/8 degree assimilative model, might have skill only on some sections in the region.

Description: Coordinated Ocean-ice Reference Experiments (COREs) are presented as a tool to explore the behaviour of global ocean-ice models under forcing from a common atmospheric dataset. We highlight issues arising when designing coupled global ocean and sea ice experiments, such as difficulties formulating a consistent forcing methodology and experimental protocol. Particular focus is given to the hydrological forcing, the details of which are key to realizing simulations with stable meridional overturning circulations. The atmospheric forcing from [Large, W., Yeager, S., 2004. Diurnal to decadal global forcing for ocean and sea-ice models: the data sets and flux climatologies. NCAR Technical Note: NCAR/TN-460+STR. CGD Division of the National Center for Atmospheric Research] was developed for coupled-ocean and sea ice models. We found it to be suitable for our purposes, even though its evaluation originally focussed more on the ocean than on the sea-ice. Simulations with this atmospheric forcing are presented from seven global ocean-ice models using the CORE-I design (repeating annual cycle of atmospheric forcing for 500 years). These simulations test the hypothesis that global ocean-ice models run under the same atmospheric state produce qualitatively similar simulations. The validity of this hypothesis is shown to depend on the chosen diagnostic. The CORE simulations provide feedback to the fidelity of the atmospheric forcing and model configuration, with identification of biases promoting avenues for forcing dataset and/or model development.

Description: The greater Agulhas Current system has several components with high mesoscale turbulence. The phytoplankton distribution in the southwest Indian Ocean reflects this activity. We have used a regional eddy-permitting, coupled physical–biological model to study the physical–biological interactions and to address the main processes responsible for phytoplankton distribution in three different biogeochemical provinces: the southwest Subtropical Indian Gyre (SWSIG), the subtropical convergence zone (SCZ) and the subantarctic waters (SAW) south of South Africa. The biological model with four compartments (Nitrate–Phytoplankton–Zooplankton–Detritus) adequately reproduces the observed field of chlorophyll a. The phase of the strong modelled seasonality in the SWSIG is opposite to that of the SCZ that forms the southern boundary of the subtropical gyre. Phytoplankton concentrations are governed by the source-minus-sink terms, which are one order of magnitude greater than the dynamical diffusion and advection terms. North of 35°S, in the SWSIG, phytoplankton growth is limited by nutrients supply throughout the year. However, deeper stratification, enhanced cross-frontal transport and higher detritus remineralization explain the simulated higher concentrations of phytoplankton found in winter in the SWSIG. The region between 35° and 40°S constitutes a transition zone between the SCZ and the oligotrophic subtropical province. Horizontal advection is the main process bringing nutrients for phytoplankton growth. The front at 34°S represents a dynamical barrier to an extension further to the north of this advection of nutrients. Within the SCZ, primary production is high during spring and summer. This high productivity depletes the nutrient standing stock built up during winter time. In winter, nutrients supply in the convergence zone is indeed large, but the deep mixing removes phytoplankton from the euphotic zone and inhibits photosynthesis, yielding lower surface chlorophyll a concentrations. Waters south of the Subantarctic Front have a summer biomass close to that of frontal waters and higher than for subtropical waters. However, these simulated concentrations are slightly higher than the observed ones suggesting that limitation by iron and/or silica may play a role