The relation of meridional pressure gradients to North Atlantic deep water volume transport in an ocean general circulation model

We use a coarse resolution ocean general circulation model to study the relation between meridional pressure and density gradients in the Southern Ocean and North Atlantic and the Atlantic meridional overturning circulation. In several experiments, we artificially modify the meridional density gradients by applying different magnitudes of the Gent–McWilliams isopycnal eddy diffusion coefficients in the Southern Ocean and in the North Atlantic and investigate the response of the simulated Atlantic meridional overturning to such changes. The simulations are carried out close to the limit of no diapycnal mixing, with a very small explicit vertical diffusivity and a tracer advection scheme with very low implicit diffusivities. Our results reveal that changes in eddy diffusivities in the North Atlantic affect the maximum of the Atlantic meridional overturning, but not the outflow of North Atlantic Deep Water into the Southern Ocean. In contrast, changes in eddy diffusivities in the Southern Ocean affect both the South Atlantic outflow of North Atlantic Deep Water and the maximum of the Atlantic meridional overturning. Results from these experiments are used to investigate the relation between meridional pressure gradients and the components of the Atlantic meridional overturning. Pressure gradients and overturning are found to be linearly related. We show that, in our simulations, zonally averaged deep pressure gradients are very weak between 20°S and about 30°N and that between 30°N and 60°N the zonally averaged pressure grows approximately linearly with latitude. This pressure difference balances a westward geostrophic flow at 30–40°N that feeds the southbound deep Atlantic western boundary current. We extend our analysis to a large variety of experiments in which surface freshwater forcing, vertical mixing and winds are modified. In all experiments, the pycnocline depth, assumed to be the relevant vertical scale for the northward volume transport in the Atlantic, is found to be approximately constant, at least within the coarse vertical resolution of the model. The model behaviour hence cannot directly be related to conceptual models in which changes in the pycnocline depth determine the strength of Atlantic meridional flow, and seems conceptually closer to Stommel’s box model. In all our simulations, the Atlantic overturning seems to be mainly driven by Southern Ocean westerlies. However, the actual strength of the Atlantic meridional overturning is not determined solely by the Southern Ocean wind stress but as well by the density/pressure gradients created between the deep water formation regions in the North Atlantic and the inflow/outflow region in the South Atlantic.     

Stability dependent parameterisations of vertical mixing in ocean general circulation models

Summary Parameterisations of mixing induced through shear instability, internal wave breaking, and double diffusion are investigated in simulations of ocean climate using a global ocean general circulation model (OGCM). Focus is placed on the sensitivity of the large scale circulation, water mass formation and transport of heat as measures of the model's ability to represent current climate. The model resolution is typical of OGCMs being coupled to atmospheric. GCMs in climate models and the parameterisations investigated are all computationally inexpensive enough to allow for integrations on long time scales. Under the assumption of constant vertical eddy coefficients (the control case), the model climatology displays acceptable values of North Atlantic Deep Water formation, Antarctic Circumpolar Current (ACC) transport, and Indonesian through-flow but an excessively deep and diffuse pycnocline structure with weak stratification in the deep ocean. It is found that various circulation and water mass properties are sensitive to the choice of parameterisation of vertical mixing and that determining a scheme which works satisfactorily over all regions (tropical, mid-latitude, and polar) of the domain is not straightforward. Parameterisations of internal wave breaking or upper ocean shear instability lead to some improvements in the model water mass formation. ACC and poleward heat transport when compared to the control case whereas parameterisations of double diffusive processes did not. Based on these and other results, various recommendations are made for mixing parameterisations in ocean climate models.With 8 Figures     

The influence of numerical advection schemes on the results of ocean general circulation models

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.This paper was presented at the International Conference on Modelling of Global Climate Change and Variability, held in Hamburg 11–15 September 1989 under the auspices of the Meteorological Institute of the University of Hamburg and the Max Planck Institute for Meteorology. Guest Editor for these papers is Dr. L. Dümenil     

On the control of subantarctic stratification by the ocean circulation

Climate Dynamics - A shallow mixed layer depth bias in Austral winter in the Subantarctic Zone is a common feature of Coupled Model Intercomparison Project (CMIP5) models, including the Community...     

A numerical world ocean general circulation model

This paper describes a numerical model of the world ocean based on the fully primitive equations. A “Standard” ocean state is introduced into the equations of the model and the perturbed thermodynamic variables are used in the modle’s calculations. Both a free upper surface and a bottom topography are included in the model and a sigma coordinate is used to normalize the model’s vertical component. The model has four unevenly-spaced layers and 4 × 5 horizontal resolution based on C-grid system. The finite-difference scheme of the model is designed to conserve the gross available energy in order to avoid fictitious energy generation or decay.The model has been tested in response to the annual mean surface wind stress, sea level air pressure and sea level air temperature as a preliminary step to its further improvement and its coupling with a global atmospheric general circulation model. Some of results, including currents, temperature and sea surface elevation simulated by the model are presented.     

The axial angular momentum balance of a global ocean general circulation model

In this study we examine the axial angular momentum balance of a non-eddy-resolving global ocean general circulation model, from the perspective of the geographical and seasonal variability of angular momentum and from the perspective of the torques acting on the ocean through its surfaces. Our purpose is to provide an estimate of the magnitude of the seasonal storage of angular momentum in the ocean and hence the oceanic excitation of variability in length of day, and to elucidate the role of the ocean in transferring angular momentum between the atmosphere and the Earth's crust. We provide an assessment of the reliability of the model results by examining the sensitivity of the angular momentum and torque distributions to several model parameters.Although the Southern Ocean region containing the Antarctic Circumpolar Current (ACC) makes the largest contribution to both the annual mean oceanic angular momentum and its seasonal variability, inclusion of the rest of the world ocean reduces both of these quantities to about two-thirds of the value of the Southern Ocean alone. The annual, global mean angular momentum is found to be insensitive to most model choices except for the isopycnal diffusivity. The seasonal variability, on the other hand, is insensitive to the isopycnal diffusivity, but sensitive to the smoothness of the representation of topography and moderately sensitive to horizontal and vertical friction parameterizations. The torque balance at all latitudes, including within the Antarctic circumpolar belt, is between wind stress and bottom pressure torques. Horizontal friction torques are small but non-negligible. Bottom friction and storage of angular momentum are negligible in angular momentum budgets on seasonal time scales. Two commonly used wind stress climatologies, one based on historical marine meteorological observations and the other based on operational weather analyses, differ in the sign of the globally integrated wind stress torque.     

Internal secular variability in an ocean general circulation model

We describe results of an experiment in which the Hamburg Large-Scale Geostrophic Ocean General Circulation Model was driven by a spatially correlated white-noise freshwater flux superimposed on the climatological fluxes. In addition to the red-noise character of the oceanic response, the model exhibits pronounced variability in a frequency band around 320 years. The centers of action of this oscillation are the Southern Ocean and the Atlantic.This paper was presented at the International Conference on Modelling of Global Climate Change and Variability, held in Hamburg 11–15 September 1989 under the auspices of the Meteorological Institute of the University of Hamburg and the Max Planck Institute for Meteorology. Guest Editor for these papers is Dr. L. Dümenil.     

The seasonal cycle in coupled ocean-atmosphere general circulation models

We examine the seasonal cycle of near-surface air temperature simulated by 17 coupled ocean-atmosphere general circulation models participating in the Coupled Model Intercomparison Project (CMIP). Nine of the models use ad hoc “flux adjustment” at the ocean surface to bring model simulations close to observations of the present-day climate. We group flux-adjusted and non-flux-adjusted models separately and examine the behavior of each class. When averaged over all of the flux-adjusted model simulations, near-surface air temperature falls within 2?K of observed values over the oceans. The corresponding average over non-flux-adjusted models shows errors up to ～6?K in extensive ocean areas. Flux adjustments are not directly applied over land, and near-surface land temperature errors are substantial in the average over flux-adjusted models, which systematically underestimates (by ～5?K) temperature in areas of elevated terrain. The corresponding average over non-flux-adjusted models forms a similar error pattern (with somewhat increased amplitude) over land. We use the temperature difference between July and January to measure seasonal cycle amplitude. Zonal means of this quantity from the individual flux-adjusted models form a fairly tight cluster (all within ～30% of the mean) centered on the observed values. The non-flux-adjusted models perform nearly as well at most latitudes. In Southern Ocean mid-latitudes, however, the non-flux-adjusted models overestimate the magnitude of January-minus-July temperature differences by ～5?K due to an overestimate of summer (January) near-surface temperature. This error is common to five of the eight non-flux-adjusted models. Also, over Northern Hemisphere mid-latitude land areas, zonal mean differences between July and January temperatures simulated by the non-flux-adjusted models show a greater spread (positive and negative) about observed values than results from the flux-adjusted models. Elsewhere, differences between the two classes of models are less obvious. At no latitude is the zonal mean difference between averages over the two classes of models greater than the standard deviation over models. The ability of coupled GCMs to simulate a reasonable seasonal cycle is a necessary condition for confidence in their prediction of long-term climatic changes (such as global warming), but it is not a sufficient condition unless the seasonal cycle and long-term changes involve similar climatic processes. To test this possible connection, we compare seasonal cycle amplitude with equilibrium warming under doubled atmospheric carbon dioxide for the models in our data base. A small but positive correlation exists between these two quantities. This result is predicted by a simple conceptual model of the climate system, and it is consistent with other modeling experience, which indicates that the seasonal cycle depends only weakly on climate sensitivity.     

Towards a parametrization of river discharges into ocean general circulation models: a closure through energy conservation

Diagnostic methods are defined in order to compare two numerical simulations of ocean dynamics in a region of freshwater influence. The first one is a river plume simulation based on a high resolution numerical configuration of the POM coastal ocean model in which mixing parametrizations have been previously defined. The second one is a simulation based on the NEMO Global Ocean Model used for climate simulations in its half-a-degree configuration in which a river inflow is represented as precipitation on two coastal grid cells. Both simulations are forced with the same freshwater inflows and wind stresses. The divergence of volumetric fluxes above and below the halocline are compared. Results show that when an upwelling wind blows, the two models display similar behavior although the impact of lack of precision can be observed in the NEMO configuration. When a downwelling wind blows, the NEMO Global Ocean configuration can not reproduce the coastally trapped baroclinic dynamics because its grid resolution is too coarse. To find a parametrization to help represent these dynamics in ocean general circulation models, a method based on energy conservation is investigated. This method shows that it is possible to link the energy fluxes provided by river inflows to the divergence of energy fluxes integrated over the grid cells of ocean general circulation models. A parametrization of the dynamics created by freshwater inflows is deduced from this method. This enabled creation of a box model that proved to have the same behavior as the fluxes previously computed from the high resolution configuration.     

Effect of seasonal forcing on global circulation in a world ocean general circulation model

 Effects of the seasonal variation in thermohaline and wind forcing on the abyssal circulation are investigated by using an ocean general circulation model. To isolate effects of the seasonality in the thermohaline forcing from those in the wind forcing, we carry out three experiments with (1) annual-mean wind forcing and perpetual-winter thermohaline forcing, (2) annual-mean wind forcing and seasonal thermohaline forcing, and (3) seasonal wind forcing and seasonal thermohaline forcing. The deep water under the seasonal thermohaline forcing becomes warmer than under the perpetual-winter thermohaline forcing. Although the perpetual-winter thermohaline forcing is widely used and believed to reproduce the deep water better than the annual-mean forcing, the difference between the results of the perpetual-winter and the seasonal thermohaline forcing is significant. The seasonal variation of the Ekman convergence and divergence produces meridional overturning cells extending to the bottom because the period of seasonal cycle is shorter than the adjustment timescale by baroclinic Rossby waves. The heat transport owing to those Ekman flows and temperature anomalies makes the upper water (0–200 m) colder at low to mid-latitudes (40S–40N) and warmer at high latitudes. Also the deep water becomes warmer owing to the warming of the northern North Atlantic, the main source region of North Atlantic Deep Water. The model is also synchronously (i.e., without acceleration) integrated with seasonal forcing for 5400 y. A past study suggested that under seasonal forcing, a sufficient equilibrium state can be achieved after only decades of synchronous integration following more than 10 000 y of accelerated integration. Here, the result so obtained is compared with that of the 5400-y synchronous integration. The difference in the global average temperature is as small as 0.12 °C, and most of the difference is confined to the Southern Ocean. Received: 1 May 1998 / Accepted: 5 January 1999     

A global ocean biogeochemistry general circulation model and its simulations

An ocean biogeochemistry model was developed and incorporated into a global ocean general circulation model (LICOM) to form an ocean biogeochemistry general circulation model (OBGCM). The model was used to study the natural carbon cycle and the uptake and storage of anthropogenic CO2 in the ocean. A global export production of 12.5 Pg C yr-1 was obtained. The model estimated that in the pre-industrial era the global equatorial region within 15o of the equator released 0.97 Pg C yr-1 to the atmosphere, which was balanced by the gain of CO2 in other regions. The post-industrial air-sea CO2 flux indicated the oceanic uptake of CO2 emitted by human activities. An increase of 20－50 mol kg-1 for surface dissolved inorganic carbon (DIC) concentrations in the 1990s relative to pre-industrial times was obtained in the simulation, which was consistent with data-based estimates. The model generated a total anthropogenic carbon inventory of 105 Pg C as of 1994, which was within the range of estimates by other researchers. Various transports of both natural and anthropogenic DIC as well as labile dissolved organic carbon (LDOC) were estimated from the simulation. It was realized that the Southern Ocean and the high-latitude region of the North Pacific are important export regions where accumulative air-sea CO2 fluxes are larger than the DIC inventory, whereas the subtropical regions are acceptance regions. The interhemispheric transport of total natural carbon (DIC+LDOC) was found to be northward (0.11 Pg C yr-1), which was just balanced by the gain of carbon from the atmosphere in the Southern Hemisphere.     

Error evolution in the dynamics of an ocean general circulation model

The problem of error propagation is considered for spatially uncorrelated errors of the barotropic stream function in an oceanic general circulation model (OGCM). Such errors typically occur when altimetric data from satellites are assimilated into ocean models. It is shown that the error decays at first due to the dissipation of the smallest scales in the error field. The error then grows exponentially before it saturates at the value corresponding to the difference between independent realizations. A simple analytic formula for the error behavior is derived; it matches the numerical results documented for the present primitive-equation ocean model, and other models in the literature.     

The seasonal cycle of mixed layer temperatures in a global ocean general circulation model

We examine mixed layer temperatures in a global ocean general circulation model subjected to seasonally varying climatological forcing. Harmonic analysis of monthly mixed layer temperatures and climatological sea surface temperatures (SSTs) shows that, on the average, the annual harmonic accounts for 90% of the total seasonal variance in both fields, while the semiannual harmonic accounts for about 8%. The semiannual signal is mostly confined to equatorial and high-latitude regions. The model mixed layer temperatures underestimate the mean amplitude of the annual harmonic in middle latitudes (65°||10°) by about 26%, while lagging climatological SSTs by 22 days, on average. In several parameter sensitivity experiments, these differences could be reduced to as little as 12% and 12.5 days, respectively, though most of this gain occurred when the mixed layer was unrealistically shallow (mean depth less than 65 m). At least part of the differences in amplitudes and phases of the annual harmonic is linked to the uncoupled formulation of the surface heat flux, which is computed using specified and seasonally varying climatological air temperatures. In ice-free areas, seasonal amplitudes and phases of air temperatures are almost identical to those of climatological SSTs. Thus, differences between model mixed layer temperatures and climatological SSTs give rise to Newtonian relaxation to SSTs, which then leads to amplitude damping and time lags in mixed layer temperatures relative to the SSTs.     

The impact of a downslope water-transport parametrization in a global ocean general circulation model

The impact of a downslope water-transport parametrization on the circulation and water mass characteristics of a global depth-level ocean general circulation model is investigated. The spreading of dense water from the formation regions into the deep ocean is known to be poorly represented in depth-level models with no bottom boundary layer resolved or attached. The new scheme is simple and intends to parametrize the effects of various oceanographic processes (rather than the processes themselves) that help dense water to descend topographic slopes by which the formation regions are separated from the world ocean. The new scheme significantly improves the large scale properties of the North Atlantic Deep Water. Changes in the North Atlantic circulation, however, are rather small. In the Southern Ocean, the exchange between the dense water formation regions on the continental shelves and the deep ocean is strengthened at the expense of deep water mass formation by open ocean convection. In all three ocean basins, the density of the deep and bottom water is higher with the new parametrization, which brings the simulations closer to observations in the Atlantic and Indian Oceans. In the Pacific Ocean, however, where the density has already been well reproduced without the downslope transport, it becomes slightly too high. The results are in agreement with those from other model studies.