Climate feedbacks in a general circulation model incorporating prognostic clouds

 This study performs a comprehensive feedback analysis on the Bureau of Meteorology Research Centre General Circulation Model, quantifying all important feedbacks operating under an increase in atmospheric CO₂. The individual feedbacks are analysed in detail, using an offline radiation perturbation method, looking at long- and shortwave components, latitudinal distributions, cloud impacts, non-linearities under 2xCO₂ and 4xCO₂ warmings and at interannual variability. The water vapour feedback is divided into terms due to moisture height and amount changes. The net cloud feedback is separated into terms due to cloud amount, height, water content, water phase, physical thickness and convective cloud fraction. Globally the most important feedbacks were found to be (from strongest positive to strongest negative) those due to water vapour, clouds, surface albedo, lapse rate and surface temperature. For the longwave (LW) response the most important term of the cloud ‘optical property’ feedbacks is due to the water content. In the shortwave (SW), both water content and water phase changes are important. Cloud amount and height terms are also important for both LW and SW. Feedbacks due to physical cloud thickness and convective cloud fraction are found to be relatively small. All cloud component feedbacks (other than height) produce conflicting LW/SW feedbacks in the model. Furthermore, the optical property and cloud fraction feedbacks are also of opposite sign. The result is that the net cloud feedback is the (relatively small) product of conflicting physical processes. Non-linearities in the feedbacks are found to be relatively small for all but the surface albedo response and some cloud component contributions. The cloud impact on non-cloud feedbacks is also discussed: greatest impact is on the surface albedo, but impact on water vapour feedback is also significant. The analysis method here proves to be a␣powerful tool for detailing the contributions from different model processes (and particularly those of the clouds) to the final climate model sensitivity. Received: 15 June 2000 / Accepted: 10 January 2001     

Sensitivity of simulated salinities in a three-dimensional ocean general circulation model to vertical mixing of destabilizing surface fluxes

 We present simulations performed with a three dimensional global ocean general circulation model which show that simulated salinities and amounts of convective mixing are very sensitive to vertical mixing of surface buoyancy fluxes. If, as usual, surface buoyancy fluxes are placed entirely in the topmost model level, our model produces excessive convective mixing in the Southern Ocean. This results in poor stimulated salinity in the Southern Ocean. In this simulation, we assume, as usual, that both surface buoyancy forcing and vertical mixing are homogeneous within each grid cell. If, on the other hand, destabilizing surface fluxes are instantaneously mixed into the subsurface ocean, the model produces much less convective mixing and much more realistic salinities. The vertical mixing of surface buoyancy fluxes performed in this simulation is equivalent to assuming that those fluxes affect only a small fraction of each grid cell, and cause vertical mixing only in that limited area. Our interpretation of these results is that the usual assumption that both surface buoyancy forcing and vertical mixing are uniform within each grid cell has a detrimental effect on model results; these results could be significantly improved by good parametrizations which treat the horizontal inhomogeneity of surface buoyancy forcing and of vertical mixing. Received: 25 February 1998 / Accepted: 9 September 1998     

Fluctuation dissipation in a general circulation model

This paper considers the climate response to step function changes in the solar constant in two versions of a general circulation model with simplified geography. The NCAR CCM0 model is applied to an all-land planet with no topography (Terra Blanda). In one version there is moisture in the air (as well as self-generated clouds) as evaporated from an ideal surface at a fixed 80% of saturation. In the other version there is no moisture in the atmosphere. We examine the decay of natural anomalies in the large-scale temperature field in each model and compare the time dependence of the ensemble average with the average temporal behavior of the response to step function changes in the solar constant. The fluctuation-dissipation theorem of statistical mechanics makes specific predictions about the relationship between the two curves. We conduct the experiments for both versions of the model since the sensitivity is quite different for each. The theorem is found to hold reasonably well in each case.     

A modelling study of aerosol processing by stratocumulus clouds and its impact on general circulation model parameterisations of cloud and aerosol

A model of the aqueous phase processing of an aerosol population undergoing multiple cycling through a stratocumulus (Sc) cloud layer is presented. Results indicate that a significant modification of the aerosol properties is achieved following the first cycle through cloud. In a polluted atmosphere, further modification in subsequent cycles is seen to be hydrogen peroxide limited unless there is a flux of ammonia entering the system through cloud base (CB). The modification of the aerosol population is seen to have little effect on the microphysics (specifically the cloud droplet concentration and effective radius) of the processing cloud. However, it enables processed aerosols to subsequently act as efficient cloud condensation nuclei (CCN) in less vigorous clouds (as a result of reducing the critical supersaturation required to activate them). The effects of variations in the internal mixture of soluble components of aerosols on the microphysics of clouds forming on them are also investigated using the cloud model. A (K2) parameterisation of the effects of variations in internally mixed nitrate loadings on the cloud droplet number concentration is presented. The effects of applying this K2 correction to the droplet number (derived from a parameterisation based on sulphate) for the presence of nitrate in aerosol have been investigated using the HadAM3 version of the Hadley Centre General Circulation Model (GCM). The effect on global annual mean simulations of the indirect forcing and effective radius is small, but more pronounced regionally. Suggestions (based on model results and observations) for parameterising the size distribution and in-cloud growth of aerosols for use in GCMs are presented.     

A look at the ocean in the EC-Earth climate model

EC-Earth is a newly developed global climate system model. Its core components are the Integrated Forecast System (IFS) of the European Centre for Medium Range Weather Forecasts (ECMWF) as the atmosphere component and the Nucleus for European Modelling of the Ocean (NEMO) developed by Institute Pierre Simon Laplace (IPSL) as the ocean component. Both components are used with a horizontal resolution of roughly one degree. In this paper we describe the performance of NEMO in the coupled system by comparing model output with ocean observations. We concentrate on the surface ocean and mass transports. It appears that in general the model has a cold and fresh bias, but a much too warm Southern Ocean. While sea ice concentration and extent have realistic values, the ice tends to be too thick along the Siberian coast. Transports through important straits have realistic values, but generally are at the lower end of the range of observational estimates. Exceptions are very narrow straits (Gibraltar, Bering) which are too wide due to the limited resolution. Consequently the modelled transports through them are too high. The strength of the Atlantic meridional overturning circulation is also at the lower end of observational estimates. The interannual variability of key variables and correlations between them are realistic in size and pattern. This is especially true for the variability of surface temperature in the tropical Pacific (El Ni?o). Overall the ocean component of EC-Earth performs well and helps making EC-Earth a reliable climate model.     

Multi-year predictability in a coupled general circulation model

Multi-year to decadal variability in a 100-year integration of a BMRC coupled atmosphere-ocean general circulation model (CGCM) is examined. The fractional contribution made by the decadal component generally increases with depth and latitude away from surface waters in the equatorial Indo-Pacific Ocean. The relative importance of decadal variability is enhanced in off-equatorial “wings” in the subtropical eastern Pacific. The model and observations exhibit “ENSO-like” decadal patterns. Analytic results are derived, which show that the patterns can, in theory, occur in the absence of any predictability beyond ENSO time-scales. In practice, however, modification to this stochastic view is needed to account for robust differences between ENSO-like decadal patterns and their interannual counterparts. An analysis of variability in the CGCM, a wind-forced shallow water model, and a simple mixed layer model together with existing and new theoretical results are used to improve upon this stochastic paradigm and to provide a new theory for the origin of decadal ENSO-like patterns like the Interdecadal Pacific Oscillation and Pacific Decadal Oscillation. In this theory, ENSO-driven wind-stress variability forces internal equatorially-trapped Kelvin waves that propagate towards the eastern boundary. Kelvin waves can excite reflected internal westward propagating equatorially-trapped Rossby waves (RWs) and coastally-trapped waves (CTWs). CTWs have no impact on the off-equatorial sub-surface ocean outside the coastal wave guide, whereas the RWs do. If the frequency of the incident wave is too high, then only CTWs are excited. At lower frequencies, both CTWs and RWs can be excited. The lower the frequency, the greater the fraction of energy transmitted to RWs. This lowers the characteristic frequency (reddens the spectrum) of variability off the equator relative to its equatorial counterpart. At low frequencies, dissipation acts as an additional low pass filter that becomes more effective, as latitude increases. At the same time, ENSO-driven off-equatorial surface heating anomalies drive mixed layer temperature responses in both hemispheres. Both the eastern boundary interactions and the accumulation of surface heat fluxes by the surface mixed layer act to low pass filter the ENSO-forcing. The resulting off-equatorial variability is therefore more coherent with low pass filtered (decadal) ENSO indices [e.g. NINO3 sea-surface temperature (SST)] than with unfiltered ENSO indices. Consequently large correlations between variability and NINO3 extend further poleward on decadal time-scales than they do on interannual time-scales. This explains why decadal ENSO-like patterns have a broader meridional structure than their interannual counterparts. This difference in appearance can occur even if ENSO indices do not have any predictability beyond interannual time-scales. The wings around 15–20°S, and sub-surface variability at many other locations are predictable on interannual and multi-year time-scales. This includes westward propagating internal RWs within about 25° of the equator. The slowest of these take up to 4 years to reach the western boundary. This sub-surface predictability has significant oceanographic interest. However, it is linked to only low levels of SST variability. Consequently, extrapolation of delayed action oscillator theory to decadal time-scales might not be justified.     

Actinic fluxes: The role of clouds in photodissociation

The role of clouds in photodissociation is examined by both modeling and observations. It is emphasized that the photodissociation rate is proportional to the actinic flux rather than to the irradiance. (The actinic flux concerns the energy that is incident on a molecule, irrespective of the direction of incidence. The irradiance concerns the energy that is incident on a plane.) A 3-layer model is used to calculate the actinic flux above and below a cloud, relative to the incident flux, in terms of cloud albedo, zenith angle and the albedo of the underlying and overlying atmosphere. Cloud albedo is mainly determined by cloud optical thickness. An expression for the in-cloud actinic flux is given as a function of in-cloud optical thickness. The 3-layer model seems to be an useful model for estimation of photodissociation rates in dispersion models. Further, a multi-layer delta-Eddington model is used to calculate irradiances, actinic fluxes and photodissociation rates of nitrogen dioxide J(NO₂) as a function of height in inhomogeneous atmospheres. For the considered wavelength interval [290–420 nm], Rayleigh scattering, ozone absorption and Mie scattering and absorption by cloud drops and aerosols should be taken into account. It is stressed that both models are one-dimensional and as such are unable to deal with partial cloudiness. It is shown that if no clouds are present, the actinic flux depends primarily on the solar zenith angle. The actinic flux usually increases with height. For cloudy atmospheres, another important parameter with respect to the actinic flux is added: cloud optical thickness, which determines cloud albedo. It can be shown that in-cloud characteristics and cloud height are less important in describing the effect of a cloud on the actinic flux (outside the cloud). The in-cloud values of the actinic flux can exceed the values outside the cloud. Finally, using the photostationary state relationship, good agreement is found between model results and aircraft measurements.     

Simulation of thermohaline circulation with a twenty-layer oceanic general circulation model

Summary This paper presents the basic configuration and preliminary performance of a twenty-layer oceanic general circulation model which represents a portion of the recent progress in developing coupled ocean-atmosphere general circulation models made by the authors. The model uses latitude/depthdependent thermohaline-stratification subduction, -coordinate, three-dimensional implicit diffusion, complete convective adjustment, separating and coupling of external and internal modes and Asselin temporal filter, and thermodynamic sea-ice calculation. With seasonally varying climatological forcing at the surface and enhanced surface salinities in the region adjacent Antarctica, the model has been integrated for one thousand years to reach a quasiequilibrium state. Preliminary verification shows that the model is capable of simulating successfully not only many aspects of the upper ocean circulation but also an acceptable thermohaline circulation. The modelled overturning rate of the North Atlantic Deep Water (NADW) is greater than 15Sv. The simulated overturning rate of the Antarctic Bottom Water (AABW) is about 20Sv. The southward outflow of NADW can be identified from not only the meridional overturning streamfunction but also the current fields at four deeper levels from 1455m to 2475m. The AABW northward outflow exists at some bottom levels below 2600m, and mainly flows towards the Pacific basin.Major problems in the present simulation include the underestimate of the NADW outflow, the failure to simulate the Antarctic Intermediate Water (AAIW), the too fresh bottom water and the too diffuse thermocline of the model. A sensitivity experiment has revealed that the model diffusion process has an important impact on the simulation of both the thermocline and the NADW outflow.With 16 Figures     

The stability of the North Atlantic thermohaline circulation in a coupled ocean-atmosphere general circulation model

 The stability of the Atlantic thermohaline circulation against meltwater input is investigated in a coupled ocean-atmosphere general circulation model. The meltwater input to the Labrador Sea is increased linearly for 250 years to a maximum input of 0.625 Sv and then reduced again to 0 (both instantaneously and linearly decreasing over 250 years). The resulting freshening forces a shutdown of the formation of North Atlantic deepwater and a subsequent reversal of the thermohaline circulation of the Atlantic, filling the deep Atlantic with Antarctic bottom water. The change in the overturning pattern causes a drastic reduction of the Atlantic northward heat transport, resulting in a strong cooling with maximum amplitude over the northern North Atlantic and a southward shift of the sea-ice margin in the Atlantic. Due to the increased meridional temperature gradient, the intertropical convergence zone over the Atlantic is displaced southward and the westerlies in the Northern Hemisphere gain strength. We identify four main feedbacks affecting the stability of the thermohaline circulation: the change in the overturning circulation of the Atlantic leads to longer residence times of the surface water in high-northern latitudes, which allows them to accumulate more precipitation and runoff from the continents. As a consequence the stratification in the North Atlantic becomes more stable. This effect is further amplified by an enhanced northward atmospheric water vapour transport, which increases the freshwater input into the North Atlantic. The reduced northward oceanic heat transport leads to colder sea-surface temperatures and an intensification of the atmospheric cyclonic circulation over the Norwegian Sea. The associated Ekman transports cause increased upwelling and increased freshwater export with the East Greenland Current. Both the cooling and the wind-driven circulation changes largely compensate for the effects of the first two feedbacks. The wind-stress feedback destabilizes modes without deep water formation in the North Atlantic, but has been neglected in almost all studies so far. After the meltwater input stops, the North Atlantic deepwater formation resumed in all experiments and the meridional overturning returned within 200 years to a conveyor belt pattern. This happened although the formation of North Atlantic deep water was suppressed in one experiment for more than 300 years and the Atlantic overturning had settled into a circulation pattern with Antarctic bottom water as the only source of deep water. It is a clear indication that cooling and wind-stress feedback are more effective, at least in our model, than advection feedback and increased atmospheric water vapour transport. We conclude that the conveyor belt-type thermohaline circulation seems to be much more stable than hitherto assumed from experiments with simpler models. Received 31 January 1996/Accepted 22 August 1996     

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.     

A diagnostic study of the southern hemisphere summer circulation of the CCC general circulation model

Abstract

The medium‐scale wave regime, consisting largely of zonal wavenumbers 5–7, frequently dominates the summer Southern Hemisphere tropospheric circulation. We perform a diagnostic study of this circulation as simulated by the Canadian Climate Centre (CCC) general circulation model (GCM). The analysis of Hövmöller diagrams, space‐time and zonal wavenumber spectra shows that the CCC GCM is able to simulate the observed medium‐scale wave regime.

The zonally averaged meridional eddy heat and momentum transports and the associated baroclinic and barotropic energy conversions are also examined. The distributions of the transports on the vertical plane agree well with the observations. After comparison with the observed December‐January‐February 1979 distributions, some quantitative differences remain: the heat transport is too weak aloft and too large near the surface, whereas the momentum transport tends to be too weak. The baroclinic and barotropic conversions show a maximum in the medium‐scale waves. The time evolution of the Richardson number of the mean flow suggests that the medium‐scale wave is due to a baroclinic instability.     

Comparison of model predicted cloud parameters and surface radiative fluxes with observations on the 100 km scale

Summary Cloud parameters and surface radiative fluxes predicted by regional atmospheric models are directly compared with observations for a 10-day period in late summer 1995 characterized by predominantly large-scale synoptic conditions. Observations of total cloud cover and vertical cloud structure are inferred from measurements with a ground-based network of Lidar ceilometers and IR-radiometers and from satellite observations on a 100 kilometer scale. Ground-based observations show that at altitudes below 3 km, implying liquid water clouds, there is a considerable portion of optically non-opaque clouds. Vertical distributions of cloud temperatures simultaneously inferred from the ground-based infrared radiometer network and from satellite can only be reconciled if the occurrence of optically thin cloud structures at mid- and high tropospheric levels is assumed to be frequent. Results of three regional atmospheric models, i.e. the GKSS-REMO, SMHI-HIRLAM, and KNMI-RACMO, are quantitatively compared with the observations. The main finding is that all models predict too much cloud amount at low altitude below 900 hPa, which is then compensated by an underestimation of cloud amount around 800 hPa. This is likely to be related with the finding that all models tend to underestimate the planetary boundary layer height. All models overpredict the high-level cloud amount albeit it is difficult to quantify to what extent due to the frequent presence of optically thin clouds. Whereas reasonably alike in cloud parameters, the models differ considerably in radiative fluxes. One model links a well matching incoming solar radiation to a radiatively transparent atmosphere over a too cool surface, another model underpredicts incoming solar radiation at the surface due to a too strong cloud feedback to radiation, the last model represents all surface radiative fluxes quite well on average, but underestimates the sensitivity of atmospheric transmissivity to cloud amount. Received August 31, 2000 Revised March 15, 2001     

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.     

On the vertical profile of surface radiative fluxes in southwest Germany

Summary Vertical profile of surface radiative fluxes in an area of heterogeneous terrain in south-west Germany is presented. Main data sets utilized for the study were recorded during the REgio KLIma Projekt (REKLIP). Supporting observational data were provided by the German weather service and German geophysical consultant service. Elevation of the study sites ranged from 212 m a.s.l. to 1489 m a.s.l. From May to September, monthly mean albedo was generally low at the study sites, ranging from 19% to 24%. For the other months, monthly mean albedo lie between 22% and 25% at the lowland site but extended between 27% and 71% at the highly elevated mountain site. Following the altitudinal increase in surface albedo, net radiative flux and radiation efficiency declined with elevation at an annual mean of 1.15 Wm⁻²/100 m and 0.008/100 m respectively. Absorbed shortwave radiation and effective terrestrial radiation showed mean decline of 1.54 Wm⁻²/100 m and 0.34 Wm⁻²/100 m, respectively, with the mean sky-to-earth radiation deficit amounting to about 52 Wm⁻² for the lowland site and 73 Wm⁻² for the highest elevated site. Some empirical models which express shortwave and longwave radiative fluxes in terms of meteorological variables have been validated for the lowland and mountain sites. Monthly mean daily total estimates of solar radiation obtained from ?ngst?m-Prescott relation were quite consistent with observed values. Parameterisation of downward atmospheric radiation under all sky condition was achieved by extending Brutsaert clear sky atmospheric model. Relationship between outgoing longwave radiation and screen temperature at the study sites was best described by an exponential function unlike the linear relationship proposed by Monteith and Unsworth. Net radiative flux for the lowland and mountain sites has been expressed in terms of absorbed shortwave radiation, cloud amount and screen temperature. Received March 5, 2001 Revised October 29, 2001     

Sensitivity of Atlantic meridional overturning circulation to the dynamical framework in an ocean general circulation model

The horizontal coordinate systems commonly used in most global ocean models are the spherical latitude–longitude grid and displaced poles, such as a tripolar grid. The effect of the horizontal coordinate system on Atlantic meridional overturning circulation (AMOC) is evaluated by using an OGCM (ocean general circulation model). Two experiments are conducted with the model—one using a latitude–longitude grid (referred to as Lat_1) and the other using a tripolar grid (referred to as Tri). The results show that Tri simulates a stronger North Atlantic deep water (NADW) than Lat_1, as more saline water masses enter the Greenland–Iceland–Norwegian (GIN) seas in Tri. The stronger NADW can be attributed to two factors. One is the removal of the zonal filter in Tri, which leads to an increasing of the zonal gradient of temperature and salinity, thus strengthening the north geostrophic flow. In turn, it decreases the positive subsurface temperature and salinity biases in the subtropical regions. The other may be associated with topography at the North Pole, because realistic topography is applied in the tripolar grid while the latitude–longitude grid employs an artificial island around the North Pole. In order to evaluate the effect of the filter on AMOC, three enhanced filter experiments are carried out. Compared to Lat_1, an enhanced filter can also augment NADW formation, since more saline water is suppressed in the GIN seas, but accumulated in the Labrador Sea, especially in experiment Lat_2_S, which is the experiment with an enhanced filter on salinity.     

Climate predictability experiments with a general circulation model

The atmospheric response to the evolution of the global sea surface temperatures from 1979 to 1992 is studied using the Max-Planck-Institut 19 level atmospheric general circulation model, ECHAM3 at T 42 resolution. Five separate 14-year integrations are performed and results are presented for each individual realization and for the ensemble-averaged response. The results are compared to a 30-year control integration using a climate monthly mean state of the sea surface temperatures and to analysis data. It is found that the ECHAM3 model, by and large, does reproduce the observed response pattern to El Nino and La Niña. During the El Nino events, the subtropical jet streams in both hemispheres are intensified and displaced equatorward, and there is a tendency towards weak upper easterlies over the equator. The Southern Oscillation is a very stable feature of the integrations and is accurately reproduced in all experiments. The inter-annual variability at middle- and high-latitudes, on the other hand, is strongly dominated by chaotic dynamics, and the tropical SST forcing only modulates the atmospheric circulation. The potential predictability of the model is investigated for six different regions. Signal to noise ratio is large in most parts of the tropical belt, of medium strength in the western hemisphere and generally small over the European area. The ENSO signal is most pronounced during the boreal spring. A particularly strong signal in the precipitation field in the extratropics during spring can be found over the southern United States. Western Canada is normally warmer during the warm ENSO phase, while northern Europe is warmer than normal during the ENSO cold phase. The reason is advection of warm air due to a more intense Pacific low than normal during the warm ENSO phase and a more intense Icelandic low than normal during the cold ENSO phase, respectively.     

Parameterization of the shortwave radiative properties of water clouds for use in GCMs

Summary The influence of the micro- and macrophysics of water clouds on the scattering and radiative properties of clouds is investigated using versatile cloud drop size distributions (DSDs) and Mie theory for single scattering and the delta-Eddington approximation for multiple scattering. A new parameterization scheme for the shortwave radiative properties of water clouds is presented. As for single- scattering properties, a new parameterization for cloud optical thickness () is proposed. This is based upon the seperation of the dependence of on the total number of DSDs, the cloud thickness, and the liquid water content, combined with equivalent radius. The cloud bulk radiative properties, i.e., the flux reflectance, transmittance, and absorptance, are uniquely fitted by a dimensionless parameter (X) or the optical thickness. The parameterization is compared with other schemes. Finally, the features and potential use of the scheme are discussed.With 5 Figures     

A comparison of solar radiative flux above clouds from MODIS with BALTIMOS regional climate model simulations

A comparison study for the solar radiative flux above clouds is presented between the regional climate model system BALTEX integrated model system (BALTIMOS) and Moderate Resolution Imaging Spectroradiometer (MODIS) satellite observations. For MODIS, an algorithm has been developed to retrieve reflected shortwave fluxes over clouds. The study area is the Baltic Sea catchment area during an 11-month period from February to December 2002. The intercomparison focuses on the variations of the daily and seasonal cycle and the spatial distributions. We found good agreement between the observed and the simulated data with a bias of the temporal mean of 13.6 W/m² and a bias of the spatial mean of 35.5 W/m². For summer months, BALTIMOS overestimates the solar flux with up to 90 W/m² (20%). This might be explained by the insufficient representation of cirrus clouds in the regional climate model.     

Isopycnal mixing and convective adjustment in an ocean general circulation model

Abstract

Convective adjustment is examined in an ocean general circulation model which uses an isopycnal mixing parametrization. It is found that the use of an explicit convective adjustment scheme is not needed in a variety of equilibria and climate change scenario simulations. A numerical mechanism is proposed to explain this as well as the localized appearance of ‘negative’ diffusion.