Uncertainties in air quality model predictions

As a result of several air quality model evaluation exercises involving a large number of source scenarios and types of models, it is becoming clear that the magnitudes of the uncertainties in model predictions are similar from one application to another. When considering continuous point sources and receptors at distances of about 0.1 km to 1 km downwind, the uncertainties in ground-level concentration predictions lead to typical mean biases of about ±20 to 40% and typical relative root-mean-square errors of about 60 to 80%. In fact, in two otherwise identical model applications at two independent sites, it is not unusual for the same model to overpredict by 50% at one site and underpredict by 50% at the second site. It is concluded that this fundamental level of model uncertainty is likely to exist due to data input errors and stochastic fluctuations, no matter how sophisticated a model becomes. The tracer studies that lead to these conclusions and have been considered in this study include: (1) tests of the Offshore and Coastal Dispersion (OCD) model at four coastal sites; (2) tests of the Hybrid Plume Dispersion Model (HPDM) at five power plants; (3) tests of a similarity model for near-surface point source releases at four sites; and (4) tests of 14 hazardous gas models at eight sites including six sets of experiments where dense gases were released.     

A 3-D model for the Antarctic ice sheet: a sensitivity study on the glacial-interglacial contrast

On the longer climatic time scales, changes in the elevation and extent of the Antarctic ice sheet have an important role in modulating global atmospheric and oceanographic processes, and contribute significantly to world-wide sea levels. In this paper, a 3-D time-dependent thermomechanical model for the entire ice sheet is presented, that is subsequently used to examine the effects of glacial-interglacial shifts in environmental boundary conditions on its geometry. The model takes into account a coupled ice shelf, grounding-line dynamics, basal sliding and isostatic bed adjustment and considers the fully coupled velocity and temperature fields. Ice flow is calculated on a fine mesh (40 km horizontal grid size and 10 layers in the vertical) for grounded and floating ice and a stress transition zone in between at the grounding line, where all stress components contribute in the effective stress in the flow law. There is free interaction between ice sheet and ice shelf, so that the entire geometry is internally generated. A simulation of the present ice sheet reveals that the model is able to yield realistic results. A series of sensitivity experiments are then performed, in which lower temperatures, reduced accumulation rates and lower global sea level stands are imposed, either singly or in combination. By comparing results of pairs of experiments, the effects of each of these environmental changes can be determined. In agreement with glacial-geological evidence, we found that the most pronounced changes show up in the West Antarctic ice sheet configuration. They appear to be essentially controlled by variations in eustatic sea level, whereas typical glacial-interglacial changes in temperature and ice deposition rates tend to balance one another. These findings support the hypothesis that the Antarctic ice sheet basically follows glacial episodes in the northern hemisphere by means of sea-level teleconnections. Grounding occurs more readily in the Weddell sea than in the Ross sea and long time scales appear to be involved: it may take up to 30–40000 years for these continental shelf areas to become completely grounded after an initial stepwise perturbation in boundary conditions. According to these reconstructions, a steady state Antarctic ice sheet may contribute some 16 m to global sea level lowering at maximum glaciation.     

Uncertainties in simulating regional climate of Southern Africa: sensitivity to physical parameterizations using WRF

This study aims at quantifying seasonal biases of regional climate model outputs during southern African summer, against a dense in situ measurement network (daily rain-gauge and surface air temperature records, and 12?h UTC radiosondes), and uncertainties associated with some physical parameterizations. Using the non-hydrostatic Advanced Research Weather Forecast (WRF) laterally forced by ERA40 reanalysis, twenty-seven experiments configured with three schemes of cumulus (CU), planetary boundary layer (PBL) and microphysics (MP), are performed at 35?km horizontal resolution during the core of a summer rainy season (December 1993 to February 1994 season) representative of the South African rainfall climatology. WRF simulates accurately seasonal large-scale rainfall patterns, as well as seasonal gradients of South African rainfall and 2-m temperature, and seasonal vertical profiles of the air temperature and humidity. However seasonal biases fluctuate strongly from an experiment to another, denoting considerable uncertainties generated by the physical package. Rainfall amounts are the most sensitive parameter to the tested schemes. Their geography, intensity, and intraseasonal characteristics are predominantly sensitive to CU schemes, and much less to PBL and MP schemes. Some CU-PBL combinations produce additive effects, which can dramatically either reduce or increase biases. Satisfactory configurations are found for South African climate, which would not have been possible without testing numerous physical parameterizations.     

An application of an efficient non-hydrostatic mesoscale model

This paper deals with a non-hydrostatic mesoscale model that achieves full vectorization on computers like the CYBER 205. The model formulation ensures the conservation of all fluxes and takes into account the terrain inhomogeneities by the aid of suitable transformations. The diagnostic equation for the pressure change is solved using a very efficient vectorized elliptic solver. By imposing appropriate boundary conditions no additional precautions at the boundaries are necessary to achieve meaningful results. As an application, the steady-state inviscid flow over a single mountain is simulated.     

Tropical convective responses to microphysical and radiative processes: a sensitivity study with a 2-D cloud resolving model

Summary Prognostic cloud schemes are increasingly used in weather and climate models in order to better treat cloud-radiation processes. Simplifications are often made in such schemes for computational efficiency, like the scheme being used in the National Centers for Environment Prediction models that excludes some microphysical processes and precipitation-radiation interaction. In this study, sensitivity tests with a 2-D cloud resolving model are carried out to examine effects of the excluded microphysical processes and precipitation-radiation interaction on tropical thermodynamics and cloud properties. The model is integrated for 10 days with the imposed vertical velocity derived from the Tropical Ocean Global Atmosphere Coupled Ocean–Atmosphere Response Experiment. The experiment excluding the depositional growth of snow from cloud ice shows anomalous growth of cloud ice and more than 20% increase of fractional cloud cover, indicating that the lack of the depositional snow growth causes unrealistically large mixing ratio of cloud ice. The experiment excluding the precipitation-radiation interaction displays a significant cooling and drying bias. The analysis of heat and moisture budgets shows that the simulation without the interaction produces more stable upper troposphere and more unstable mid and lower troposphere than does the simulation with the interaction. Thus, the suppressed growth of ice clouds in upper troposphere and stronger radiative cooling in mid and lower troposphere are responsible for the cooling bias, and less evaporation of rain associated with the large-scale subsidence induces the drying in mid and lower troposphere.     

A degradation approach to accelerate simulations to steady-state in a 3-D tracer transport model of the global ocean

 We have developed a new method to accelerate tracer simulations to steady-state in a 3-D global ocean model, run off-line. Using this technique, our simulations for natural ¹⁴C ran 17 times faster when compared to those made with the standard non-accelerated approach. For maximum acceleration we wish to initialize the model with tracer fields that are as close as possible to the final equilibrium solution. Our initial tracer fields were derived by judiciously constructing a much faster, lower-resolution (degraded), off-line model from advective and turbulent fields predicted from the parent on-line model, an ocean general circulation model (OGCM). No on-line version of the degraded model exists; it is based entirely on results from the parent OGCM. Degradation was made horizontally over sets of four adjacent grid-cell squares for each vertical layer of the parent model. However, final resolution did not suffer because as a second step, after allowing the degraded model to reach equilibrium, we used its tracer output to re-initialize the parent model (at the original resolution). After re-initialization, the parent model must then be integrated only to a few hundred years before reaching equilibrium. To validate our degradation-integration technique (DEGINT), we compared ¹⁴C results from runs with and without this approach. Differences are less than 10‰ throughout 98.5% of the ocean volume. Predicted natural ¹⁴C appears reasonable over most of the ocean. In the Atlantic, modeled Δ¹⁴C indicates that as observed, the North Atlantic Deep Water (NADW) fills the deep North Atlantic, and Antartic Intermediate Water (AAIW) infiltrates northward; conversely, simulated Antarctic Bottom Water (AABW) does not penetrate northward beyond the equator as it should. In the Pacific, in surface eastern equatorial waters, the model produces a north–south assymetry similar to that observed; other global ocean models do not, because their resolution is inadequate to resolve equatorial dynamics properly, particularly the intense equatorial undercurrent. The model’s oldest water in the deep Pacific (at −239‰) is close to that observed (−248‰), but is too deep. Surface waters in the Southern Ocean are too rich in natural ¹⁴C due to inadequacies in the OGCM’s thermohaline forcing. Received: 18 March 1997 / Accepted: 27 July 1997     

Climate sensitivity due to increased CO2: experiments with a coupled atmosphere and ocean general circulation model

A version of the National Center for Atmospheric Research community climate model — a global, spectral (R15) general circulation model — is coupled to a coarse-grid (5° latitude-] longitude, four-layer) ocean general circulation model to study the response of the climate system to increases of atmospheric carbon dioxide (CO₂). Three simulations are run: one with an instantaneous doubling of atmospheric CO₂ (from 330 to 660 ppm), another with the CO₂ concentration starting at 330 ppm and increasing linearly at a rate of 1% per year, and a third with CO₂ held constant at 330 pm. Results at the end of 30 years of simulation indicate a globally averaged surface air temperature increase of 1.6° C for the instantaneous doubling case and 0.7°C for the transient forcing case. Inherent characteristics of the coarse-grid ocean model flow sea-surface temperatures (SSTs) in the tropics and higher-than-observed SSTs and reduced sea-ice extent at higher latitudes] produce lower sensitivity in this model after 30 years than in earlier simulations with the same atmosphere coupled to a 50-m, slab-ocean mixed layer. Within the limitations of the simulated meridional overturning, the thermohaline circulation weakens in the coupled model with doubled CO₂ as the high-latitude ocean-surface layer warms and freshens and westerly wind stress is decreased. In the transient forcing case with slowly increasing CO₂ (30% increase after 30 years), the zonal mean warming of the ocean is most evident in the surface layer near 30°–50° S. Geographical plots of surface air temperature change in the transient case show patterns of regional climate anomalies that differ from those in the instantaneous CO₂ doubling case, particularly in the North Atlantic and northern European regions. This suggests that differences in CO₂ forcing in the climate system are important in CO₂ response in regard to time-dependent climate anomaly regimes. This confirms earlier studies with simple climate models that instantaneous CO₂ doubling simulations may not be analogous in all respects to simulations with slowly increasing CO₂.A portion of this study is supported by the US Department of Energy as part of its Carbon Dioxide Research Program     

Regional climate model sensitivity to domain size

Regional climate models are increasingly used to add small-scale features that are not present in their lateral boundary conditions (LBC). It is well known that the limited area over which a model is integrated must be large enough to allow the full development of small-scale features. On the other hand, integrations on very large domains have shown important departures from the driving data, unless large scale nudging is applied. The issue of domain size is studied here by using the “perfect model” approach. This method consists first of generating a high-resolution climatic simulation, nicknamed big brother (BB), over a large domain of integration. The next step is to degrade this dataset with a low-pass filter emulating the usual coarse-resolution LBC. The filtered nesting data (FBB) are hence used to drive a set of four simulations (LBs for Little Brothers), with the same model, but on progressively smaller domain sizes. The LB statistics for a climate sample of four winter months are compared with BB over a common region. The time average (stationary) and transient-eddy standard deviation patterns of the LB atmospheric fields generally improve in terms of spatial correlation with the reference (BB) when domain gets smaller. The extraction of the small-scale features by using a spectral filter allows detecting important underestimations of the transient-eddy variability in the vicinity of the inflow boundary, which can penalize the use of small domains (less than 100 × 100 grid points). The permanent “spatial spin-up” corresponds to the characteristic distance that the large-scale flow needs to travel before developing small-scale features. The spin-up distance tends to grow in size at higher levels in the atmosphere.     

Climate model sensitivity to sea ice albedo parameterization

Summary Three one-year experimental simulations with the National Center for Atmospheric Research Community Climate Model (NCAR CCM) were performed with three sea ice albedo parameterizations and compared with control run results to examine their impact on polar surface temperature, planetary albedo and clouds. The first integration utilized sea ice albedos of the Arctic Basin for the spring and summer of 1977 derived from defence Meteorological Satellite Imagery (DMSP). The second simulation employed prescribed lead and melt pond fractions and an albedo weighting scheme. The third simulation involved the coupling of an interactive sea ice/snow albedo parameterization made a function of surface state.Results show that prescribed, and assumed true satellite sea ice albedos produced higher planetary albedos than those calculated with the standard CCM sea ice albedo scheme in the control run. As a result, lower temperatures (up to 0.5 K) and increased cloudiness are generated for the Arctic region. The standard CCM sea ice albedo scheme is used as an adjustment to maintain normal temperatures for the polar oceans. The radiative impact of leads and melt ponds warmed sea ice regions only for short time periods. The third scheme generated markedly lower planetary albedos (reductions of 0.07 to 0.17) and higher surface temperatures (up to 2.0 K) than control values.The CCM simulates a gradual decrease in spring and summer Arctic cloud cover whereas observations show a sharp spring increase. Examination of the CCM code, particularly the cloud parameterization, is required to address this problem.With 12 Figures     

Millennial timescale carbon cycle and climate change in an efficient Earth system model

A new Earth system model, GENIE-1, is presented which comprises a 3-D frictional geostrophic ocean, phosphate-restoring marine biogeochemistry, dynamic and thermodynamic sea-ice, land surface physics and carbon cycling, and a seasonal 2-D energy-moisture balance atmosphere. Three sets of model climate parameters are used to explore the robustness of the results and for traceability to earlier work. The model versions have climate sensitivity of 2.8–3.3°C and predict atmospheric CO₂ close to present observations. Six idealized total fossil fuel CO₂ emissions scenarios are used to explore a range of 1,100–15,000 GtC total emissions and the effect of rate of emissions. Atmospheric CO₂ approaches equilibrium in year 3000 at 420–5,660 ppmv, giving 1.5–12.5°C global warming. The ocean is a robust carbon sink of up to 6.5 GtC year⁻¹. Under ‘business as usual’, the land becomes a carbon source around year 2100 which peaks at up to 2.5 GtC year⁻¹. Soil carbon is lost globally, boreal vegetation generally increases, whilst under extreme forcing, dieback of some tropical and sub-tropical vegetation occurs. Average ocean surface pH drops by up to 1.15 units. A Greenland ice sheet melt threshold of 2.6°C local warming is only briefly exceeded if total emissions are limited to 1,100 GtC, whilst 15,000 GtC emissions cause complete Greenland melt by year 3000, contributing 7 m to sea level rise. Total sea-level rise, including thermal expansion, is 0.4–10 m in year 3000 and ongoing. The Atlantic meridional overturning circulation shuts down in two out of three model versions, but only under extreme emissions including exotic fossil fuel resources.     

A 3-D sea breeze model of the PBL including pollutant dispersion

A 3-D sea breeze model has been developed to study the sea breeze and the pollutant dispersion over a coastal area. Topography is also included. Similarity theory is adopted while a turbulent kinetic energy equation is used to study turbulence and to parameterize eddy diffusion. The integration is repeated in diurnal cycles and the results — meteorological and pollutant — are compared with data from the Athens area.     

CO2 climate sensitivity and snow-sea-ice albedo parameterization in an atmospheric GCM coupled to a mixed-layer ocean model

The snow-sea-ice albedo parameterization in an atmospheric general circulation model (GCM), coupled to a simple mixed-layer ocean and run with an annual cycle of solar forcing, is altered from a version of the same model described by Washington and Meehl (1984). The model with the revised formulation is run to equilibrium for 1 × CO₂ and 2 × CO₂ experiments. The 1 ×CO₂ (control) simulation produces a global mean climate about 1° warmer than the original version, and sea-ice extent is reduced. The model with the altered parameterization displays heightened sensitivity in the global means, but the geographical patterns of climate change due to increased carbon dioxide (CO₂) are qualitatively similar. The magnitude of the climate change is affected, not only in areas directly influenced by snow and ice changes but also in other regions of the globe, including the tropics where sea-surface temperature, evaporation, and precipitation over the oceans are greater. With the less-sensitive formulation, the global mean surface air temperature increase is 3.5 °C, and the increase of global mean precipitation is 7.12%. The revised formulation produces a globally averaged surface air temperature increase of 4.04 °C and a precipitation increase of 7.25%, as well as greater warming of the upper tropical troposphere. Sensitivity of surface hydrology is qualitatively similar between the two cases with the larger-magnitude changes in the revised snow and ice-albedo scheme experiment. Variability of surface air temperature in the model is comparable to observations in most areas except at high latitudes during winter. In those regions, temporal variation of the sea-ice margin and fluctuations of snow cover dependent on the snow-ice-albedo formulation contribute to larger-than-observed temperature variability. This study highlights an uncertainty associated with results from current climate GCMs that use highly parameterized snow-sea-ice albedo schemes with simple mixed-layer ocean models.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

Uncertainties linked to land-surface processes in climate change simulations

 The impact of climate change on the hydrology of continental surfaces is critical for human activities but the response of the surface to this perturbation may also affect the sensitivity of the climate. This complex feedback is simulated in general circulation models (GCMs) used for climate change predictions by their land-surface schemes. The present study attempts to quantify the uncertainty associated with these schemes and what impact it has on our confidence in the simulated climate anomalies. Four GCMs, each coupled to two different land-surface schemes, are used to explore the spectrum of uncertainties. It is shown that, in this sample, surface processes have a significant contribution to our ability to predict surface temperature changes and perturbations of the hydrological cycle in an environment with doubled greenhouse gas concentration. The results reveal that the uncertainty introduced by land-surface processes in the simulated climate is different from its impact on the sensitivity of GCMs to climate change, indeed an alteration of the surface parametrization with little impact on model climate can affect sensitivity significantly. This result leads us to believe that the validation of land-surface schemes should not be limited to the current climate but should also cover their sensitivity to variations in climatic forcing. Received: 24 June 1999 / Accepted: 20 April 2000     

The 3-D model characteristics of a Cb cloud which moves along a valley

Summary The three-dimensional cloud-resolving mesoscale model is used to simulate an individual Cb cloud in condition of real orography. We have conducted our numerical experiments over an area known as hail bearing clouds source. Once formed such clouds regenerate and propagate along the valley if the shallow layer of strong wind shear exists. The orographic effects on model Cb cloud are recognized through comparison of simulated cloud characteristics with those calculated for the flat terrain. Sensitivity experiments with respect to the wind shear layer depth are also conducted.Our results demonstrate the model capability to simulate well some observational Cb cloud characteristics. It is shown that the river valley is of essential importance for Cb cloud development. The most prominent features of the model cloud in this case are as follows: fast propagation along the river valley; considerable depth of the cold air nose with pronounced pulsation mechanism and intense cell regeneration at the leading edge of cold air outflow. Model and observed radar reflectivities in the vertical cross-section are in agreement. Some characteristics of surface cumulative rain precipitation are also well reproduced by the model. In contrast with the real orography case, the model cloud is more intense and it propagates freely in lateral direction for the flat terrain. The cell regeneration associated with forced updraft above the cold air nose is not pronounced in contrast with earlier considerations. Reflectivity pattern near the ground, having mainly bat-like wings, encircles much larger area with altitude compared to the case of real orography. Finally, the model cloud characteristics depend strongly on shear layer.Received June 2002; revised August 22, 2002; accepted October 1, 2002 Published online: April 10, 2003     

Seeding agent dispersion within convective cloud as simulated by a 3-D numerical model

Summary The three-dimensional mesoscale cloud-resolving model ARPS (Advanced Regional Prediction System) was used to investigate the dispersal of an inert seeding agent within a cumulonimbus (Cb) cloud developing from two different initial states. In this paper, we stress the influence vortices in the cloud have on seeding agent dispersion. If a strong directional ambient wind shear is present in the lowest layer, a vortex pair formed at the flanks of the simulated cloud. Following the velocity field, a considerable amount of the injected seeding agent would be thrown out to the rear of the cloud, where both updrafts associated with vortices and downdrafts occurred. After a short time the agent was present only in the cloud periphery. If the Cb cloud developed under conditions where directional ambient wind shear did not exist, seeding agent dispersion would be quite different. In this case, almost all the seeding agent was transported into the main updraft, while the residence time of the agent within the cloud was longer due to the weaker cloud dynamics. Therefore, we must pay attention to whether or not the cloud contains vortices when we make the decision where to seed. This is necessary in order to minimize the loss of seeding material.     

Quantifying the sensitivity of simulated climate change to model configuration

This study used “factor separation” to quantify the sensitivity of simulated present and future surface temperatures and precipitation to alternative regional climate model physics components. The method enables a quantitative isolation of the effects of using each physical component as well as the combined effect of two or more components. Simulation results are presented from eight versions of the Mesoscale Modeling System Version 5 (MM5), one-way nested within one version of the Goddard Institute for Space Studies Atmosphere-Ocean Global Climate Model (GISS AOGCM). The MM5 simulations were made at 108 km grid spacing over the continental United States for five summers in the 1990s and 2050s. Results show that the choice of cumulus convection parameterization is the most important “factor” in the simulation of contemporary surface summer temperatures and precipitation over both the western and eastern USA. The choice of boundary layer scheme and radiation package also increases the range of model simulation results. Moreover, the alternative configurations give quite different results for surface temperature and precipitation in the 2050s. For example, simulated 2050s surface temperatures by the scheme with the coolest 1990s surface temperatures are comparable to 1990s temperatures produced by other schemes. The study analyzes the spatial distribution of 1990s to 2050s projected changes in the surface temperature for the eight MM5 versions. The predicted surface temperature change at a given grid point, averaged over all eight model configurations, is generally about twice the standard deviation of the eight predicted changes, indicating relative consensus among the different model projections. Factor separation analysis indicates that the choice of cumulus parameterization is the most important modeling factor amongst the three tested contributing to the computed 1990s to 2050s surface temperature change, although enhanced warming over many areas is also attributable to synergistic effects of changing all three model components. Simulated ensemble mean precipitation changes, however, are very small and generally smaller than the inter-model standard deviations. The MM5 versions therefore offer little consensus regarding 1990s to 2050s changes in precipitation rates.     

Climate sensitivity to changes in solar insolation in a simple coupled climate model

A simple coupled ocean, atmosphere and sea-ice model is presented. The idealised model consists of a zonally averaged land and ocean strip of constant angular width extending from pole to pole. The meridional energy transport in the ocean is modelled by contributions from the large scale thermohaline overturning cells and from horizontal diffusive fluxes. The atmospheric meridional energy transports are parametrised as diffusive fluxes in addition to advective transports in the Hadley domain. This parametrisation resolves the equatorward moisture transport as well as the poleward transport of potential energy in the upper branch of the Hadley circulation. The model reproduces the annual averaged meridional energy transports in the climate system with a small number of free model parameters. The basic feedbacks between the three climatic components are studied by investigating the model's sensitivity towards reductions in the solar insolation. It is found that the meridional energy transport in the ocean does not amplify the ice albedo feedback. This has important implications for modelling the climate sensitivity in atmosphere-only models, as these would exaggerate the sensitivity to changes in the solar insolation if their parametrisations of the meridional energy transport are constrained by surface temperatures. The role of the dependence of the atmospheric transports on the meridional temperature gradient is shown to have a significant influence on the sensitivity on the coupled model, and the inclusion of seasonal cycles greatly increase the models sensitivity. The Hadley circulation does significantly alter the strength of the ice-albedo feedback in the coupled model. The idealised configuration of the model makes it a useful tool for studying the feedbacks in the ocean-atmosphere-sea ice system in the context of the "Snowball Earth" hypothesis.     

Coupling an ocean wave model to an atmospheric general circulation model

A coupled model, consisting of an ocean wave model and an atmospheric general circulation model (AGCM), is integrated under permanent July conditions. The wave model is forced by the AGCM wind stress, whereas the wind waves modify the AGCM surface fluxes of momentum, sensible and latent heat. We investigate the following aspects of the coupled model: how realistic are the wave fields, how strong is the coupling, and how sensitive is the atmospheric circulation to the spatially and temporally varying wave field. The wave climatology of the coupled model compares favorably with observational data. The interaction between the two models is largest (although weak) in the storm track in the Southern Hemisphere. Young windsea, which is associated with enhanced surface fluxes is generated mostly in the equatorward frontal area of an individual cyclone. However, the enhancement of the surface fluxes is too small to significantly modify the climatological mean atmospheric circulation.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     

Long-term climate changes due to increased CO2 concentration in the coupled atmosphere-ocean general circulation model ECHAM3/LSG

 The long-term adjustment processes of atmosphere and ocean in response to gradually increased atmospheric CO₂ concentration have been analysed in two 850-year integrations with a coupled atmosphere-ocean general circulation model (AOGCM). In these experiments the CO₂ concentration has been increased to double and four times the initial concentration, respectively, and is kept fixed thereafter. Three characteristic time scales have been identified: a very fast response associated with processes dominated by the atmospheric adjustment, an intermediate time scale of a few decades connected with processes in the upper ocean, and adjustment processes with time scales of centuries and longer due to the inertia of the deep ocean. The latter in particular is responsible for a still ongoing adjustment of the atmosphere-ocean system at the end of the integrations after 850 years. After 60 years, at the time of CO₂ doubling, the global mean near-surface air temperature rises by 1.4 K. In spite of the constant CO₂ concentration during the following centuries the warming continues to 2.6 K after 850 years. The behaviour of the quadrupling run is similar: global mean near-surface air temperature increases by 3.8 K at the time of CO₂ quadrupling and by 4.8 K at the end of the simulation. The thermohaline circulation undergoes remarkable changes. Temporarily, the North Atlantic overturning circulation weakens by up to 30% in the CO₂ doubling experiment and up to 50% in the CO₂ quadrupling experiment. After reaching the minimum the North Atlantic overturning slowly recovers in both experiments. Received: 23 August 1999 / Accepted: 27 April 2000     

Storm track sensitivity to sea surface temperature resolution in a regional atmosphere model

A high resolution regional atmosphere model is used to investigate the sensitivity of the North Atlantic storm track to the spatial and temporal resolution of the sea surface temperature (SST) data used as a lower boundary condition. The model is run over an unusually large domain covering all of the North Atlantic and Europe, and is shown to produce a very good simulation of the observed storm track structure. The model is forced at the lateral boundaries with 15–20 years of data from the ERA-40 reanalysis, and at the lower boundary by SST data of differing resolution. The impacts of increasing spatial and temporal resolution are assessed separately, and in both cases increasing the resolution leads to subtle, but significant changes in the storm track. In some, but not all cases these changes act to reduce the small storm track biases seen in the model when it is forced with low-resolution SSTs. In addition there are several clear mesoscale responses to increased spatial SST resolution, with surface heat fluxes and convective precipitation increasing by 10–20% along the Gulf Stream SST gradient.