Description: The Joint Global Ocean Flux Study (JGOFS) has completed a decade of intensive process and time-series studies on the regional and temporal dynamics of biogeochemical processes in five diverse ocean basins. Its field program also included a global survey of dissolved inorganic carbon (DIC) in the ocean, including estimates of the exchange of carbon dioxide (CO2) between the ocean and the atmosphere, in cooperation with the World Ocean Circulation Experiment (WOCE). This report describes the principal achievements of JGOFS in ocean observations, technology development and modelling. The study has produced a comprehensive and high-quality database of measurements of ocean biogeochemical properties. Data on temporal and spatial changes in primary production and CO2 exchange, the dynamics of of marine food webs, and the availability of micronutrients have yielded new insights into what governs ocean productivity, carbon cycling and export into the deep ocean, the set of processes collectively known as the "biological pump." With large-scale, high-quality data sets for the partial pressure of CO2 in surface waters as well for other DIC parameters in the ocean and trace gases in the atmosphere, reliable estimates, maps and simulations of air-sea gas flux, anthropogenic carbon and inorganic carbon export are now available. JGOFS scientists have also obtained new insights into the export flux of particulate and dissolved organic carbon (POC and DOG), the variations that occur in the ratio of elements in organic matter, and the utilization and remineralization of organic matter as it falls through the ocean interior to the sediments. JGOFS scientists have amassed long-term data on temporal variability in the exchange of CO2 between the ocean and atmosphere, ecosystem dynamics, and carbon export in the oligotrophic subtropical gyres. They have documented strong links between these variables and large-scale climate patterns such as the El Nino-Southern Oscillation (ENSO) or the North Atlantic Oscillation (NAO). An increase in the abundance of organisms that fix free nitrogen (N-2) and a shift in nutrient limitation from nitrogen to phosphorus in the subtropical North Pacific provide evidence of the effects of a decade of strong El Ninos on ecosystem structure and nutrient dynamics. High-quality data sets, including ocean-color observations from satellites, have helped modellers make great strides in their ability to simulate the biogeochemical and physical constraints on the ocean carbon cycle and to extend their results from the local to the regional and global scales. Ocean carbon-cycle models, when coupled to atmospheric and terrestrial models, will make it possible in the future to predict ways in which land and ocean ecosystems might respond to changes in climate.

Description: The authors aim to acquaint the reader with the current state of ocean circulation models, their ability to model the present climate state and its variability, and their major shortcomings and uncertainties. They limit the discussion to three-dimensional models of the physical system. They begin by describing the basic structure of circulation models, and discussing various problems with their implementation. They give a brief overview of the types of observational data in oceanography, and the ways in which the data are used. Some results from models of the wind-driven circulation are discussed, with particular emphasis on the dynamics of mesoscale eddies. Considerable progress has been made in understanding short-term variability associated with ENSO, and the authors describe ocean-atmosphere interactions in the tropics as well as results from coupled ocean-atmosphere models for ENSO variability. Models of the thermohaline circulation are described and some emerging ideas regarding long-term changes are given.

Description: Observations from cruises in the Arabian Sea and data from satellites are interpreted using different realizations of a multi-level primitive equation model and an eddy-permitting reduced-gravity shallow water model of the Indian Ocean. The focus is on the interannual circulation variability of the Arabian Sea, and especially of the meridional location of the Great Whirl (GW). The results suggest that the variability in the western Arabian Sea is not only due to the interannual variability in the wind field, but that a substantial part is caused by the chaotic nature of the ocean dynamics. Decreasing the friction coefficient from 1000 to 500m2s-1 in a 19o numerical reduced-gravity model, the variance of the GW location increases dramatically, and the mean position moves southward by one degree. In the eddy-permitting experiments analyzed, both mechanisms appear to determine the GW location at the onset of the GW dynamics in late summer.

Description: The dependence of results from coarse-resolution models of the North Atlantic circulation on the numerical advection algorithm is studied. In particular, the sensitivity of parameters relevant for climate simulations as e.g., meridional transport of mass and heat and main thermocline thickness is investigated. Three algorithms were considered: (a) a central difference scheme with constant values for horizontal and vertical diffusion, (b) an upstream scheme with no explicit diffusion, and (c) a flux-corrected transport (FCT) scheme with constant and strictly isopycnal diffusion. The temporal evolution of the three models on time scales of centuries is markedly different, the upstream scheme resulting in much shorter adjustment time whereas the central difference scheme is slower and controlled by vertical diffusion rather than advection. In the steady state, the main thermocline structure is much less diffusive in the FCT calculation which also has much lower heat transport. Both horizontal circulation and overturning in the meridional-vertical plane are strongest in the upstream-model. The results are discussed in terms of the effective vertical (diapycnal) mixing in the different models. A significant increase in vertical resolution would be required to eliminate the high sensitivity due to the numerical algorithms, and allow physically motivated mixing formulations to become effective.

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.

Description: Assimilation experiments with data from the Bermuda Atlantic Time-series Study (BATS, 1989¯1993) were performed with a simple mixed-layer ecosystem model of dissolvedinorganic nitrogen (N), phytoplankton (P) and herbivorous zooplankton (H). Our aim is to optimize the biological model parameters, such that the misfits between model results andobservations are minimized. The utilized assimilation method is the variational adjoint technique, starting from a wide range of first-parameter guesses. A twin experiment displayedtwo kinds of solutions, when Gaussian noise was added to the model-generated data. The expected solution refers to the global minimum of the misfit model-data function, whereasthe other solution is biologically implausible and is associated with a local minimum. Experiments with real data showed either bottom-up or top-down controlled ecosystemdynamics, depending on the deep nutrient availability. To confine the solutions, an additional constraint on zooplankton biomass was added to the optimization procedure. Thisinclusion did not produce optimal model results that were consistent with observations. The modelled zooplankton biomass still exceeded the observations. From the model-datadiscrepancies systematic model errors could be determined, in particular when the chlorophyll concentration started to decline before primary production reached its maximum. Adirect comparision of measured 14C-production data with modelled phytoplankton production rates is inadequate at BATS, at least when a constant carbon to nitrogen C : N ratio isassumed for data assimilation.