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.     

A comparison of the climatology of a tropospherestratosphere-mesosphere model with observations

Mean fields from a perpetual January simulation of a GCM extending from the surface to 0.01 hPa (near 80 km) are compared to observations. The zonal mean temperature and wind fields correspond quite well with reality; the low stratosphere, especially in the polar night, is too cold, but warmer than in the original version of the model, with an upper boundary at 25 hPa. Mean fields at standard levels show that the major features of the troposphere are represented by the model, but rather over emphasised; the stratospheric winter polar vortex is too strong, too cold, and too barotropic; it resembles an `undisturbed' January rather than the climatology. Differences in the stationary eddy activity between the extended and orginal versions of the model are noted, and used to explain some differences between the two simulations.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 für Meteorology. Guest Editor for these papers is Dr. L. Dümenil     

Review of the development and applications of the Adem thermodynamic climate model

Recent advances in the development and applications of the author's Hemispheric Thermodynamic Climate Model are presented. The model has been adapted to simulate the climates from 18 kyr BP to the present, and to study the effect of the ice sheets, the insolation anomalies and the atmospheric CO₂ content on such climates. The surface ocean temperature anomaly is also simulated in the model, and comparison with values of CLIMAP (1981) for 18 kyr BP shows some agreement. A long series of numerical experiments have lead to the improvement in prediction of the monthly surface temperature anomalies. Verification of 93 predictions over the contiguous United States of America shows a useful skill in the predictions. The model is being adapted for forecasting in the Mexican Republic. Experiments to improve the skill in prediction of surface ocean temperature anomalies in the Northern Hemisphere have been carried out, and using a fine resolution grid, the model has been used to simulate the annual cycle of the normal sea surface temperatures in the Gulf of Mexico, that agrees well with observations.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     

A model of late Pleistocene ice sheet growth with realistic geography and simplified cryodynamics and geodynamics

A global two-dimensional one-level seasonal energy-balance model is asynchronously coupled to vertically integrated ice-flow models (which depend both on latitude and longitude) to study the response of the atmosphere-ocean-cryosphere-lithosphere system to solar forcing for the last ice age cycle of the late Pleistocene. The model simulates the position of the North American and European ice sheet complexes at the last glacial maximum satisfactorily. Both the geographic distributions of the ice volumes delivered by the model and their masses are a reasonable approximation to those inferred on the basis of relative sea level data (Tushingham and Peltier 1990). The sensitivity of the coupled model over the last glacial-interglacial cycle to solar forcing is nevertheless low, which suggests that further physical mechanisms will have to be added to the model (such as explicit basal sliding and ice shelves which would respond to sea-level variations and therefore permit marine incursions), if it is to adequately simulate the terminations that control the 10⁵ year ice age cycle. One should also incorporate long-term variations of the greenhouse gases (Manabe et al. 1985b).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     

Validation of an annual cycle simulation with a T42-L20 GCM

An annual cycle of an atmospheric general circulation model (AGCM) is presented. The winter and summer zonal averages of the atmospheric fields are compared with an observed climatology. The main features of the observed seasonal means are well reproduced by the model. One of the main discrepancies is that the simulated atmosphere is too cold, particularly in its upper part. Some other discrepancies might be explained by the interannual variability. The AGCM surface fluxes are directly compared to climatological estimates. On the other hand, the calculation of meridional heat transport by the ocean, inferred from the simulated energy budget, can be compared to transport induced from climatologies. The main result of this double comparison is that AGCM fluxes generally are within the range of climatological estimates. The main deficiency of the model is poor partitioning between solar and non-solar heat fluxes in the tropical belt. The meridional heat transport also reveals a significant energy-loss by the Northern Hemisphere ocean north of 45° N. The possible implications of model surface flux deficiencies on coupling with an oceanic model are discussed.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     

Refinement of dynamically downscaled precipitation and temperature scenarios

A method for adjusting dynamically downscaled precipitation and temperature scenarios representing specific sites is presented. The method reproduces mean monthly values and standard deviations based on daily observations. The trend obtained in the regional climate model both for temperature and precipitation is maintained, and the frequency of modelled and observed rainy days shows better agreement. Thus, the method is appropriate for tailoring dynamically downscaled temperature and precipitation values for climate change impact studies. One precipitation and temperature scenario dynamically downscaled with HIRHAM from the Atmospheric-Ocean General Circulation Model at the Max-Planck Institute in Hamburg, ECHAM4/OPYC4 GSDIO with emission scenario IS92a, is chosen to illustrate the adjustment method.     

Arctic radiation deficit and climate variability

One of the generally accepted climatic effects of stratospheric aerosol injection is the reduction of the global radiation in high latitudes by an order of 5% during El Chichon type eruptions. To test the effect of a high-latitude radiation deficit on global climate, a GCM experiment was performed with the ECMWF T21 atmosphere general circulation model (AGCM). The results provide physically-consistent evidence that this radiation deficit is a possible external forcing factor for severe climatic anomalies not only in the area directly affected by the reduced radiation, but also in the tropics. The most important factor is the creation of enhanced snow cover in regions of Asia that are distant from the location of the introduced radiation anomaly. The simulated results show certain features that are well known from observations in weak monsoon years, i.e. the weakened easterly jet in the upper troposphere over northern India, prolonged winter monsoon conditions, and prevailing anticyclonic vorticity anomalies over the entire Indian summer monsoon region. Over the western Pacific at the end of boreal winter (May), increased convective activity leads to a negative Walker circulation anomaly with westerly wind anomalies near the surface and easterly anomalies in the upper troposphere. This is known as one of the most important anomalies at the beginning of an El Niño/Southern Oscillation (ENSO) event.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     

A comparative study of the observations of high clouds and simulations by an atmospheric general circulation model

The importance of clouds in the upper troposphere (cirrus) for the sensitivity of the Earth's climate e.g., requires that these clouds be modeled accurately in general circulation model (GCM) studies of the atmosphere. Bearing in mind the lack of unambiguous quantitative information on the geographical distribution and properties of high clouds, the simulated distribution of upper tropospheric clouds in a spectral GCM is compared with several satellite-derived data-sets that pertain to high clouds only, for both winter and summer seasons. In the model, clouds are assumed to occupy an entire gridbox whenever the relative humidity exceeds 99%: otherwise the grid box is assumed to be free of cloud. Despite the simplicity of the cloud prediction scheme, the geographical distribution of the maxima in the model's upper tropospheric cloud cover coincides approximately with the regions of the observed maxima in the high cloud amount and their frequency of occurrence (e.g., intertropical convergence zone and the monsoon areas). These areas exhibit a minimum in the outgoing longwave radiation (OLR; Nimbus-7) and are also coincident with regions of heavy precipitation. The model, with its relatively simple cloud formation scheme, appears to capture the principal large-scale features of the tropical convective processes that are evident in the satellite and precipitation datasets, wherein the intense, upward motion is accompanied by condensation and the spreading of thick upper tropospheric layers of high relative humidity and cloudiness in the vicinity of the tropical rainbelt regions.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     

A sensitivity experiment for the removal of Arctic sea ice with the French spectral general circulation model

Sea ice has a major influence on climate in high latitudes. In this paper we analyzed the impact of removal of Arctic sea-ice cover on the climate simulated by a T42 20-level version of the French spectral model Emeraude. The control experiment was the second winter of an annual cycle simulation of the present climate. In the perturbed simulation the Arctic sea-ice cover was replaced by open ocean maintained at the freezing temperature of sea water. The zonal mean patterns of the model response were found to be in good agreement with earlier simulations of Fletcher et al. and Warshaw and Rapp. The atmospheric warming, caused by the increase of upward fluxes of sensible and latent heat and of longwave radiation from the ice-free ocean surface, is largely limited to the high latitudes poleward of 70° N and the lower half of the troposphere and leads to a surface pressure decrease and a precipitation increase over this area. We also analyze the geographical distribution of the response and the mechanisms that can explain the simulated cooling over Eurasia in relation to the energy budget at the surface. Finally, we discuss the reduction of cloud cover over the ice-free Arctic, which was an unexpected result of our simulation, and conclude that further studies are necessary to resolve the question of cumulus convection and cloud process parameterization in high latitudes.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     

A first-order global model of late Cenozoic climatic change II. Further analysis based on a simplification of CO2 dynamics

A model developed recently for the long-term variations of global ice mass, carbon dioxide, and mean ocean temperature through the late Cenozoic is simplified by hypothesizing a new equation for the CO₂ variations containing one less adjustable parameter, but retaining the essential physical content of the previous equation (including nonlinearity and the potential for instability). By assuming plausible time constants for the glacial ice mass and global mean ocean temperature, and setting the values of six adjustable parameters (rate constants), a solution for the last 5 My is obtained displaying many of the features observed over this period, including the transition to the near-100 ky major ice-age oscillations of the late Pleistocene. In obtaining this solution it is also assumed that variations in tectonic forcing lead to a reduction of the equilibrium CO₂ concentration (perhaps due to increased weathering of rapidly uplifted mountain ranges over this period). As a consequence of this CO₂ reduction, the model dynamical system can bifurcate to a free oscillatory ice-age regime that is under the pacemaker influence of earthorbital (Milankovitch) forcing. Expanded discussions are given of the surface temperature variations accompanying the evolution of ice, CO₂, and ocean temperature, and of the bifurcation properties of the model from both mathematical and physical viewpoints.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. Dumenil     

Three basic problems of paleoclimatic modeling: a personal perspective and review

The development of a theory of the evolution of the climate of the earth over millions of years can be subdivided into three fundamental, nested, problems: Firstly, to establish by equilibrium climate models (e.g., general circulation models) the diagnostic relations, valid at any time, between the fast-response climate variables (i.e., the weather statistics) and both the prescribed external radiative forcing and the prescribed distribution of the slow-response variables (e.g., the ice sheets and shelves, the deep ocean state, and the atmospheric CO₂ concentration). Secondly, to construct, by an essentially inductive process, a model of the time-dependent evolution of the slow-response climatic variables over time scales longer than the damping times of these variables but shorter than the time scale of ultra-slow tectonic and astronomical changes in the boundary conditions (e.g., altered geography and elevation of the continents, slow outgassing and weathering and and solar radiative output). Thirdly, to determine the nature of these ultra-slow processes and their effects on the evolution of the equilibrium state of the climatic system about which the previously mentioned time-dependent variations occur. In this review we discuss the basis for this resolution, and give a broad overview of the contributions that have been made thus far in each area, emphasizing the work of the Yale climate group.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     

Simulation of the response of the North Pacific Ocean to the anomalous atmospheric circulation associated with El Niño

During El Niño events when positive sea surface temperature (SST) anomalies form in the equatorial Pacific, SST anomalies also tend to develop in the North Pacific. This study attempts to model and explain the large-scale features of the observed SST anomaly field in the North Pacific during the fall and winter of the El Niño year. The experiment design consists of a mixed layer ocean model of the North Pacific which is forced by atmospheric surface fields from two sets of Community Climate Model (CCM) integrations: the El Niño set with prescribed positive SST anomalies in the tropical Pacific; and the control set which is obtained from an extended CCM integration with prescribed climatological SSTs. The response of the midlatitude ocean to atmospheric surface fields associated with El Niño is obtained by compositing each set of model integrations (El Niño and Control) and then taking the difference between the composites. The ocean model is able to reproduce the general features of the observed midlatitude SST anomaly pattern: warm water in the northeast Pacific and an elliptically shaped cold pool in the central Pacific. In these regions, a large fraction of the temperature anomalies are significant at the 95% level as indicated by a two tailed t-test. The ocean temperature anomalies simulated by the model are primarily caused by changes in the sensible and latent heat flux and to a lesser extent the longwave radiation flux. Entrainment of cold water from below the mixed layer also influences ocean temperatures. However, the entrainment anomaly pattern has a complex spatial structure which does not always coincide with the simulated mixed layer temperature anomalies.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     

Design and performance of a new cloud microphysics scheme developed for the ECHAM general circulation model

A new cloud microphysics scheme including a prognostic treatment of cloud ice (PCI) is developed to yield a more physically based representation of the components of the atmospheric moisture budget in the general circulation model ECHAM. The new approach considers cloud water and cloud ice as separate prognostic variables. The precipitation formation scheme for warm clouds distinguishes between maritime and continental clouds by considering the cloud droplet number concentration, in addition to the liquid water content. Based on several observational data sets, the cloud droplet number concentration is derived from the sulfate aerosol mass concentration as given from the sulfur cycle simulated by ECHAM. Results obtained with the new scheme are compared to satellite observations and in situ measurements of cloud physical and radiative properties. In general, the standard model ECHAM4 and also PCI capture the overall features, and the simulated results usually lie within the range of observed uncertainty. As compared to ECHAM4, only slight improvements are achieved with the new scheme. For example, the overestimated liquid water path and total cloud cover over convectively active regions are reduced in PCI. On the other hand, some shortcomings of the standard model such as underestimated shortwave cloud forcing over the extratropical oceans of the respective summer hemisphere are more pronounced in PCI.This paper was presented at the Third International Conference on Modelling of Global Climate Change and Variability, held in Hamburg 4–9 Sept. 1995 under the auspices of the Max Planck Institute for Meteorology, Hamburg. Editor for these papers is L. Dümenil.     

Specifying the atmospheric pollutant distribution forecast using the mathematical modeling and measurement data

Considered is a possibility of specifying the forecast of volcanic aerosol distribution in the atmosphere using the data of ground-based actinometric measurements. To simulate the aerosol distribution, the FlexPart software (Norwegian Institute for Air Research) was used. Using the FlexGraph software (Planeta Research Center for Space Hydrometeorology), the visualization and interactive analysis were carried out of aerosol concentration fields at different time moments. The prognostic concentration of volcanic ash from Grímsvötn volcano (May 2011) is bound to the results of measurements of aerosol optical depth in Hamburg at one of the points of AERONET global network. The comparison of the corrected aerosol concentration with the measurements at other AERONET points demonstrated the competence of the proposed approach. To specify the aerosol cloud distribution forecast on the territory of Russia, the Roshydromet actinometric network data can be used.     

A technique for verification of weather forecast and climate simulation with fuzzy sets

In the light of the idea of fuzzy neartude, a fuzzy procedure of verifying weather forecast and climate simula-tion is demonstrated. Several verification indices used in weather forecast, such as hit ratio, Heidke score, Brier score, correlation coefficient and information entropy score, are unified under fuzzy neartude and a new verifica-tion index is proposed on the basis of Hamming distance neartude. Further, D2-statistics used in climate sim-ulation verification is derived from a weighted Euclidean distance neartude. The “January climate” generated by general circulation model (GCM) is also numerically verified.     

Global warming and the growth of ice sheets

Recent research has suggested that warmer conditions, that may result from increased levels of CO₂ in the atmosphere, may induce the growth of the Northern Hemisphere ice sheets (Miller and de Vernal 1992) through the impact of warmer temperature on the water carrying capacity of the atmosphere and thus on precipitation. In this study we examine this possibility by using a coupled energy balance climate-thermo-dynamic sea ice model. Results indicate that if summer ice albedo is high enough, and there is some mechanism for initially maintaining ice through the summer season, then it may be possible to have ice sheet growth under the conditions of CO₂ induced warming.This paper was presented at the Second International Conference on Modelling of Global Climate Variability, held in Hamburg 7–11 September 1992 under the auspices of the Max Planck Institute for Meteorology. Guest Editor for these papers is L. Dümenil     

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     

Response of the Météo-France climate model to changes in CO2 and sea surface temperature

The climate response to an increase in carbon dioxide and sea surface temperatures is examined using the Météo-France climate model. This model has a high vertical resolution in the stratosphere and predicts the evolution of the ozone mixing ratio. This quantity is fully interactive with radiation and photochemical production and loss rates are accounted for. Results from a 5-year control run indicate a reasonable agreement with observed climatologies. A 5-year simulation is performed with a doubled CO₂ concentration using, as lower boundary conditions, mean surface temperatures anomalies and sea ice limits predicted for the years 56–65 of a 100-year transient simulation performed at Hamburg with a global coupled atmosphere-ocean model. The perturbed simulation produces a global mean surface air warming of 1.4 K and an increase in global mean precipitation rate of 4%. Outside the high latitudes in the Northern Hemisphere, the model simulates a strong cooling in the stratosphere reaching 10 K near the stratopause. Temperature increases are noticed in the lower polar stratosphere of the Northern Hemisphere caused by an intensification in the frequency of sudden warmings in the perturbed simulation. The low and mid-latitude stratospheric cooling leads to an ozone column enhancement of about 5%. Other features present in similar studies are exhibited in the troposphere such as the stronger surface warming over polar regions of the Northern Hemisphere, the summer time soil moisture drying in mid-latitudes and the increase in high convective cloudiness in tropical regions.This paper was presented at the Second International Conference on Modelling of Global Climate Variability, held in Hamburg 7–11 September 1992 under the auspices of the Max Planck Institute for Meteorology. Guest Editor for these papers is L. Dümenil Correspondence to: JF Mahfouf     

Water Resources of the South Asian Region in a Warmer Atmosphere

The global mean surface temperature may rise by about 0.3oC per decade during the next Few decades as a result of anthropogenic greenhouse gas emissions in the earth’s atmosphere. The data generated in the greenhouse warming simulations (Business-as-Usual scenario of IPCC) with the climate models developed at Max Planck Institute for Meteorology, Hamburg have been used to assess future plausible hydrological scenario for the South Asian region. The model results indicate enhanced surface warming (2.7oC for summer and 3.6oC for winter) over the land regions of South Asia during the next hundred years. While there is no significant change in the precipitation over most of the land regions during winter, substantial increase in precipitation is likely to occur during summer. As a result, an increase in soil moisture is likely over central India, Bangladesh and South China during summer but a statistically sig-nificant decline in soil moisture is expected over central China in winter. A moderate decrease in surface runoff may occur over large areas of central China during winter while the flood prone areas of NE-India. Bangladesh and South China are likely to have an increase in surface runoff during summer by the end of next century.     

Validation of GCM control simulations using indices of daily airflow types over the British Isles

In this study, the control simulations of two general circulation model (GCM) experiments are assessed in terms of their ability to reproduce realistic real world weather. The models examined are the UK Meteorological Office high-resolution atmospheric model (UKHI) and a coupled ocean/atmosphere model of the Max Planck Institut für Meteorologic, Hamburg (MPI). An objective classification of daily airflow patterns over the British Isles is used as a basis for comparing the frequencies of model-generated weather types with the frequencies derived from 110 years of observed mean-sea-level pressure (MSLP) fields. The weather-type frequencies generated by the GCMs, and their relationships with simulated monthly mean temperatures and total precipitation over the UK, are compared, season by season, with similar results derived using the observational data. An index of gale frequencies over the British Isles, derived from a similar objective analysis of daily MSLP fields, is used to evaluate the ability of the GCMs to simulate the observed frequency of storm events. One advantage of using 110 years of observational data is that the observed decadal-scale variability of climate can be introduced into this type of validation exercise. Both the GCMs assessed here are too cyclonic in winter. The seasonality of both anticyclonic and cyclonic types is much too strong in MPI and summer precipitation in this model is greatly underestimated. MPI simulates the annual cycle of temperature well, while UKHI successfully reproduces the annual cycle of precipitation. The analysis also indicates that the summer temperature variability of the two models is not driven by circulation changes.This paper was presented at the Second International Conference on Modelling of Global Climate Variability, held in Hamburg 7–11 September 1992 under the auspices of the Max Planck Institute for Meteorology. Guest Editor for these papers is L. Dümenil