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