A predictability study of simulated North Atlantic multidecadal variability

 The North Atlantic is one of the few places on the globe where the atmosphere is linked to the deep ocean through air–sea interaction. While the internal variability of the atmosphere by itself is only predictable over a period of one to two weeks, climate variations are potentially predictable for much longer periods of months or even years because of coupling with the ocean. This work presents details from the first study to quantify the predictability for simulated multidecadal climate variability over the North Atlantic. The model used for this purpose is the GFDL coupled ocean-atmosphere climate model used extensively for studies of global warming and natural climate variability. This model contains fluctuations of the North Atlantic and high-latitude oceanic circulation with variability concentrated in the 40–60 year range. Oceanic predictability is quantified through analysis of the time-dependent behavior of large-scale empirical orthogonal function (EOF) patterns for the meridional stream function, dynamic topography, 170 m temperature, surface temperature and surface salinity. The results indicate that predictability in the North Atlantic depends on three main physical mechanisms. The first involves the oceanic deep convection in the subpolar region which acts to integrate atmospheric fluctuations, thus providing for a red noise oceanic response as elaborated by Hasselmann. The second involves the large-scale dynamics of the thermohaline circulation, which can cause the oceanic variations to have an oscillatory character on the multidecadal time scale. The third involves nonlocal effects on the North Atlantic arising from periodic anomalous fresh water transport advecting southward from the polar regions in the East Greenland Current. When the multidecadal oscillatory variations of the thermohaline circulation are active, the first and second EOF patterns for the North Atlantic dynamic topography have predictability time scales on the order of 10–20 y, whereas EOF-1 of SST has predictability time scales of 5–7 y. When the thermohaline variability has weak multidecadal power, the Hasselmann mechanism is dominant and the predictability is reduced by at least a factor of two. When the third mechanism is in an extreme phase, the North Atlantic dynamic topography patterns realize a 10–20 year predictability time scale. Additional analysis of SST in the Greenland Sea, in a region associated with the southward propagating fresh water anomalies, indicates the potential for decadal scale predictability for this high latitude region as well. The model calculations also allow insight into regional variations of predictability, which might be useful information for the design of a monitoring system for the North Atlantic. Predictability appears to break down most rapidly in regions of active convection in the high-latitude regions of the North Atlantic. Received: 28 October 1996 / Accepted: 21 March 1997     

Contribution of glacier melt to sea-level rise since AD 1865: a regionally differentiated calculation

 The contribution of glacier melt, including the Greenland ice-sheet, to sea-level change since AD 1865 is estimated on the basis of modelled sensitivity of glacier mass balance to climate change and historical temperature data. Calculations are done in a regionally differentiated manner to overcome the inhomogeneity of the global distribution of glaciers. A distinction is made between changes in summer temperature and in temperature over the rest of the year. Our best estimate of the ice melt in the period 1865–1990 in terms of sea-level change equivalent is 5.7 cm (2.7 cm for glaciers and 3.0 cm for the Greenland ice-sheet). Additional calculations show that simpler methods, like using annual or even global mean temperature anomaly give estimates that differ by up to 55%. Consequently, a regionally differentiating approach is advised for making projections of glacier melt with GCM output. Received: 6 December 1996/Accepted: 30 May 1997     

Global dynamics of the Antarctic ice sheet

The total mass budget of the Antarctic ice sheet is studied with a simple axi-symmetrical model. The ice-sheet has a parabolic profile resting on a bed that slopes linearly downwards from the centre of the ice sheet into the ocean. The mean ice velocity at the grounding line is assumed to be proportional to the water depth. The accumulation rate is a linear function of the distance to the centre. Setting the total mass budget to zero yields a quadratic equation for the steady-state ice-sheet radius R. Analysis of the equilibrium states sheds light on the sensitivity of the ice-sheet radius to changes in sea level (S) and precipitation with respect to the present state (P_rel). For model parameters obtained by matching the analytical model to the present state of the Antarctic ice sheet, the sensitivity values are dR/dS = -2400 and dR/dP_rel = 4000 m/%. The model can also be used to study transient behaviour of the ice sheet. The characteristic relaxation time (e-folding time scale) is about 3500 years. Forcing the model with a sea-level and accumulation history over the past few hundred thousands of years yields Antarctic ice-volume curves that are similar to those obtained by comprehensive numerical modelling. The current imbalance predicted by the model corresponds to a sea-level rise of 0.25 mm yr^-1.     

The role of vegetation feedbacks on Greenland glaciation

The role of vegetation feedbacks for the process of ice-sheet evolution could potentially be important in realistically modeling the past and future evolution of the Greenland ice-sheet. We use a fully coupled atmosphere–ocean model to assess the response of the climate when the Greenland ice-sheet is replaced with a number of fixed vegetation types (bare soil, broadleaf and needleleaf trees, C3 and C4 grasses and shrubs) in conjunction with loaded and unloaded bedrock orography. These sensitivity experiments show that albedo changes dominate the climate response during the summer months while temperature changes during winter are attributed to altitude change and changes in atmospheric circulation over Greenland. Snow-free summers occur for all fixed vegetation types, except for high altitude eastern regions for bare soil. We perform further simulations with dynamic vegetation resulting in dominant shrub coverage over central and southern Greenland with grasses supported in the north. Ice-sheet modeling shows significant regrowth of the Greenland ice-sheet can occur for a bare soil surface type, dependent on ice-sheet model parameters, while Greenland remains almost ice-free for needleleaf tree coverage. Furthermore, a realistically vegetated Greenland can only support a small amount of ice-sheet regrowth implying multi-stability of the Greenland ice-sheet under a preindustrial climate. This study highlights the importance of considering vegetation climate ice-sheet interactions, and uncertainty in ice-sheet model parameters.     

GCM simulations of the Last Glacial Maximum surface climate of Greenland and Antarctica

 The LMDz variable grid GCM was used to simulate the Last Glacial Maximum (LGM, 21 ky Bp.) climate of Greenland and Antarctica at a spatial resolution of about 100 km.The high spatial resolution allows to investigate the spatial variability of surface climate change signals, and thus to address the question whether the sparse ice core data can be viewed as representative for the regional scale climate change. This study addresses primarily surface climate parameters because these can be checked against the, limited, ice core record. The changes are generally stronger for Greenland than for Antarctica, as the imposed changes of the forcing boundary conditions (e.g., sea surface temperatures) are more important in the vicinity of Greenland. Over Greenland, and to a limited extent also in Antarctica, the climate shows stronger changes in winter than in summer. The model suggests that the linear relationship between the surface temperature and inversion strength is modified during the LGM. The temperature dependency of the moisture holding capacity of the atmosphere alone cannot explain the strong reduction in snowfall over central Greenland; atmospheric circulation changes also play a crucial role. Changes in the high frequency variability of snowfall, atmospheric pressure and temperature are investigated and possible consequences for the interpretation of ice core records are discussed. Using an objective cyclone tracking scheme, the importance of changes of the atmospheric dynamics off the coasts of the ice sheets, especially for the high frequency variability of surface climate parameters, is illustrated. The importance of the choice of the LGM ice sheet topography is illustrated for Greenland, where two different topographies have been used, yielding results that differ quite strongly in certain nontrivial respects. This means that the paleo-topography is a significant source of uncertainty for the modelled paleoclimate. The sensitivity of the Greenland LGM climate to the prescribed sea surface conditions is examined by using two different LGM North Atlantic data sets. Received: 23 October 1997 / Accepted: 17 March 1998     

Snowpack model estimates of the mass balance of the Greenland ice sheet and its changes over the twentyfirst century

The performance of a snow cover model in capturing the ablation on the Greenland ice sheet is evaluated. This model allows an explicit calculation of the formation of melt water, of the fraction of melt water which re-freezes, and of runoff in the ablation region. The input climate variables to the snowpack model come from two climate models. While the higher resolution general circulation model (ECHAM 4), is closest to observations in its estimate of accumulation, it fails to give accurate results in its predictions of runoff, primarily in the southern half of the ice sheet. The two-dimensional low-resolution climate model (MIT 2D LO) produces estimates of runoff from the Greenland ice sheet within the range of uncertainty of the Inter governmental Panel on Climate Change (IPCC1) 1995 estimates. Both models reproduce some of the characteristics of the extent of the wet snow zone observed with satellite remote sensing; the MIT model is closer to observations in terms of areal extent and intensity of the melting in the southern half of the ice-sheet in July and August while the ECHAM model reproduces melting in the northern half of the ice sheet well. Changes in runoff from Greenland and Antarctica are often cited as one of the major concerns linked to anthropogenic changes in climate. Because it is based on physical principles and relies on the surface energy balance as input, the snow cover model can respond to the current climatic forcing as well as to future changes in climate on the century time scale without the limitations inherent in empirical parametrizations. For a reference climate scenario similar to the IPCC's IS92a, the model projects that the Greenland ice sheet does not contribute significantly to changes in the level of the ocean over the twenty-first century. Increases in accumulation over the central portion of the ice sheet offset most of the increase in melting and runoff, which takes place along the margins of the ice sheet. The range of uncertainty in the predictions of sea-level rise is estimated by repeating the calculation with the MIT model for seven climate change scenarios. The range is –0.5 to 1.7 cm.     

Thresholds for irreversible decline of the Greenland ice sheet

The Greenland ice sheet will decline in volume in a warmer climate. If a sufficiently warm climate is maintained for a few thousand years, the ice sheet will be completely melted. This raises the question of whether the decline would be reversible: would the ice sheet regrow if the climate cooled down? To address this question, we conduct a number of experiments using a climate model and a high-resolution ice-sheet model. The experiments are initialised with ice sheet states obtained from various points during its decline as simulated in a high-CO₂ scenario, and they are then forced with a climate simulated for pre-industrial greenhouse gas concentrations, to determine the possible trajectories of subsequent ice sheet evolution. These trajectories are not the reverse of the trajectory during decline. They converge on three different steady states. The original ice-sheet volume can be regained only if the volume has not fallen below a threshold of irreversibility, which lies between 80 and 90% of the original value. Depending on the degree of warming and the sensitivity of the climate and the ice-sheet, this point of no return could be reached within a few hundred years, sooner than CO₂ and global climate could revert to a pre-industrial state, and in that case global sea level rise of at least 1.3 m would be irreversible. An even larger irreversible change to sea level rise of 5 m may occur if ice sheet volume drops below half of its current size. The set of steady states depends on the CO₂ concentration. Since we expect the results to be quantitatively affected by resolution and other aspects of model formulation, we would encourage similar investigations with other models.     

Use of an upwelling-diffusion energy balance climate model to simulate and diagnose A/OGCM results

 We demonstrate that a hemispherically averaged upwelling-diffusion energy-balance climate model (UD/EBM) can emulate the surface air temperature change and sea-level rise due to thermal expansion, predicted by the HadCM2 coupled atmosphere-ocean general circulation model, for various scenarios of anthropogenic radiative forcing over 1860–2100. A climate sensitivity of 2.6 °C is assumed, and a representation of the effect of sea-ice retreat on surface air temperature is required. In an extended experiment, with CO₂ concentration held constant at twice the control run value, the HadCM2 effective climate sensitivity is found to increase from about 2.0 °C at the beginning of the integration to 3.85 °C after 900 years. The sea-level rise by this time is almost 1.0 m and the rate of rise fairly steady, implying that the final equilibrium value (the `commitment') is large. The base UD/EBM can fit the 900-year simulation of surface temperature change and thermal expansion provided that the time-dependent climate sensitivity is specified, but the vertical profile of warming in the ocean is not well reproduced. The main discrepancy is the relatively large mid-depth warming in the HadCM2 ocean, that can be emulated by (1) diagnosing depth-dependent diffusivities that increase through time; (2) diagnosing depth-dependent diffusivities for a pure-diffusion (zero upwelling) model; or (3) diagnosing higher depth-dependent diffusivities that are applied to temperature perturbations only. The latter two models can be run to equilibrium, and with a climate sensitivity of 3.85 °C, they give sea-level rise commitments of 1.7 m and 1.3 m, respectively. Received: 27 April 1999 / Accepted: 13 September 2000     

Greenland��s contribution to global sea-level rise by the end of the 21st century

The Greenland ice sheet holds enough water to raise the global sea level with ??7 m. Over the last few decades, observations manifest a substantial increase of the mass loss of this ice sheet. Both enhanced melting and increase of the dynamical discharge, associated with calving at the outlet-glacier fronts, are contributing to the mass imbalance. Using a dynamical and thermodynamical ice-sheet model, and taking into account speed up of outlet glaciers, we estimate Greenland??s contribution to the 21st-century global sea-level rise and the uncertainty of this estimate. Boundary fields of temperature and precipitation extracted from coupled climate-model projections used for the IPCC Fourth Assessment Report, are applied to the ice-sheet model. We implement a simple parameterization for increased flow of outlet glaciers, which decreases the bias of the modeled present-day surface height. It also allows for taking into account the observed recent increase in dynamical discharge, and it can be used for future projections associated with outlet-glacier speed up. Greenland contributes 0?C17?cm to global sea-level rise by the end of the 21st century. This range includes the uncertainties in climate-model projections, the uncertainty associated with scenarios of greenhouse-gas emissions, as well as the uncertainties in future outlet-glacier discharge. In addition, the range takes into account the uncertainty of the ice-sheet model and its boundary fields.     

A 20-year coupled ocean-sea ice-atmosphere variability mode in the North Atlantic in an AOGCM

In order to understand potential predictability of the ocean and climate at the decadal time scales, it is crucial to improve our understanding of internal variability at this time scale. Here, we describe a 20-year mode of variability found in the North Atlantic in a 1,000-year pre-industrial simulation of the IPSL-CM5A-LR climate model. This mode involves the propagation of near-surface temperature and salinity anomalies along the southern branch of the subpolar gyre, leading to anomalous sea-ice melting in the Nordic Seas, which then forces sea-level pressure anomalies through anomalous surface atmospheric temperatures. The wind stress associated to this atmospheric structure influences the strength of the East Greenland Current across the Denmark Strait, which, in turn, induces near-surface temperature and salinity anomalies of opposite sign at the entrance of the Labrador Sea. This starts the second half of the cycle after approximatively 10 years. The time scale of the cycle is thus essentially set by advection of tracers along the southern branch of the subpolar gyre, and by the time needed for anomalous East Greenland Current to accumulate heat and freshwater anomalies at the entrance of the Labrador Sea. The Atlantic meridional overturning circulation (AMOC) does not play a dominant role in the mode that is confined in the subpolar North Atlantic, but it also has a 20-year preferred timescale. This is due to the influence of the propagating salinity anomalies on the oceanic deep convection. The existence of this preferred timescale has important implications in terms of potential predictability of the North Atlantic climate in the model, although its realism remains questionable and is discussed.     

Climate-model induced differences in the 21st century global and regional glacier contributions to sea-level rise

The large uncertainty in future global glacier volume projections partly results from a substantial range in future climate conditions projected by global climate models. This study addresses the effect of global and regional differences in climate input data on the projected twenty-first century glacier contribution to sea-level rise. Glacier volume changes are calculated with a surface mass balance model combined with volume-area scaling, applied to 89 glaciers in different climatic regions. The mass balance model is based on a simplified energy balance approach, with separated contributions by net solar radiation and the combined other fluxes. Future mass balance is calculated from anomalies in air temperature, precipitation and atmospheric transmissivity, taken from eight global climate models forced with the A1B emission scenario. Regional and global sea-level contributions are obtained by scaling the volume changes at the modelled glaciers to all glaciers larger than 0.1 km² outside the Greenland and Antarctic ice sheets. This results in a global value of 0.102 ± 0.028 m (multi-model mean and standard deviation) relative sea-level equivalent for the period 2012–2099, corresponding to 18 ± 5 % of the estimated total volume of glaciers. Glaciers in the Antarctic, Alaska, Central Asia and Greenland together account for 65 ± 4 % of the total multi-model mean projected sea-level rise. The projected sea-level contribution is 35 ± 17 % larger when only anomalies in air temperature are taken into account, demonstrating an important compensating effect by increased precipitation and possibly reduced atmospheric transmissivity. The variability in projected precipitation and atmospheric transmissivity changes is especially large in the Arctic regions, making the sea-level contribution for these regions particularly sensitive to the climate model used. Including additional uncertainties in the modelling procedure and the input data, the total uncertainty estimate for the future projections becomes ±0.063 m.     

On the Response of the Greenland Ice Sheet to Greenhouse Climate Change

Numerical computations are performed with the three-dimensional polythermal ice-sheet model SICOPOLIS in order to investigate the possible impact of a greenhouse-gas-induced climate change on the Greenland ice sheet. The assumed increase of the mean annual air temperature above the ice covers a range from T = 1°C to 12°C, and several parameterizations for the snowfall and the surface melting are considered. The simulated shrinking of the ice sheet is a smooth function of the temperature rise, indications for the existence of critical thresholds of the climate input are not found. Within 1000 model years, the ice-volume decrease is limited to 10% of the present volume for T 3°C, whereas the most extreme scenario, T = 12°C, leads to an almost entire disintegration, which corresponds to a sea-level equivalent of 7 m. The different snowfall and melting parameterizations yield an uncertainty range of up to 20% of the present ice volume after 1000 model years.     

A continuous simulation of global ice volume over the past 1 million years with 3-D ice-sheet models

Sea-level records show large glacial-interglacial changes over the past million years, which on these time scales are related to changes of ice volume on land. During the Pleistocene, sea-level changes induced by ice volume are largely caused by the waxing and waning of the large ice sheets in the Northern Hemisphere. However, the individual contributions of ice in the Northern and Southern Hemisphere are poorly constrained. In this study, for the first time a fully coupled system of four 3-D ice-sheet models is used, simulating glaciations on Eurasia, North America, Greenland and Antarctica. The ice-sheet models use a combination of the shallow ice and shelf approximations to determine sheet, shelf and sliding velocities. The framework consists of an inverse forward modelling approach to derive a self-consistent record of temperature and ice volume from deep-sea benthic δ¹⁸O data over the past 1 million years, a proxy for ice volume and temperature. It is shown that for both eustatic sea level and sea water δ¹⁸O changes, the Eurasian and North American ice sheets are responsible for the largest part of the variability. The combined contribution of the Antarctic and Greenland ice sheets is about 10 % for sea level and about 20 % for sea water δ¹⁸O during glacial maxima. However, changes in interglacials are mainly caused by melt of the Greenland and Antarctic ice sheets, with an average time lag of 4 kyr between melt and temperature. Furthermore, we have tested the separate response to changes in temperature and sea level for each ice sheet, indicating that ice volume can be significantly influenced by changes in eustatic sea level alone. Hence, showing the importance of a simultaneous simulation of all four ice sheets. This paper describes the first complete simulation of global ice-volume variations over the late Pleistocene with the possibility to model changes above and below present-day ice volume, constrained by observations of benthic δ¹⁸O proxy data.     

Predictability of SST forced climate signals in two atmospheric general circulation models

 The predictability of atmospheric responses to global sea surface temperature (SST) anomalies is evaluated using ensemble simulations of two general circulation models (GCMs): the GENESIS version 1.5 (GEN) and the ECMWF cycle 36 (ECM). The integrations incorporate observed SST variations but start from different initial land and atmospheric states. Five GEN 1980–1992 and six ECM 1980–1988 realizations are compared with observations to distinguish predictable SST forced climate signals from internal variability. To facilitate the study, correlation analysis and significance evaluation techniques are developed on the basis of time series permutations. It is found that the annual mean global area with realistic signals is variable dependent and ranges from 3 to 20% in GEN and 6 to 28% in ECM. More than 95% of these signal areas occur between 35 °S–35 °N. Due to the existence of model biases, robust responses, which are independent of initial condition, are identified over broader areas. Both GCMs demonstrate that the sensitivity to initial conditions decreases and the predictability of SST forced responses increases, in order, from 850 hPa zonal wind, outgoing longwave radiation, 200 hPa zonal wind, sea-level pressure to 500 hPa height. The predictable signals are concentrated in the tropical and subtropical Pacific Ocean and are identified with typical El Ni?o/ Southern Oscillation phenomena that occur in response to SST and diabatic heating anomalies over the equatorial central Pacific. ECM is less sensitive to initial conditions and better predicts SST forced climate changes. This results from (1) a more realistic basic climatology, especially of the upper-level wind circulation, that produces more realistic interactions between the mean flow, stationary waves and tropical forcing; (2) a more vigorous hydrologic cycle that amplifies the tropical forcing signals, which can exceed internal variability and be more efficiently transported from the forcing region. Differences between the models and observations are identified. For GEN during El Ni?o, the convection does not carry energy to a sufficiently high altitude, while the spread of the tropospheric warming along the equator is slower and the anomaly magnitude smaller than observed. This impacts model ability to simulate realistic responses over Eurasia and the Indian Ocean. Similar biases exist in the ECM responses. In addition, the relationships between upper and lower tropospheric wind responses to SST forcing are not well reproduced by either model. The identification of these model biases leads to the conclusion that improvements in convective heat and momentum transport parametrizations and basic climate simulations could substantially increase predictive skill. Received: 25 April 1996 / Accepted: 9 December 1996     

Effects of a melted greenland ice sheet on climate, vegetation, and the cryosphere

This paper investigates the possible implications for the earth-system of a melting of the Greenland ice-sheet. Such a melting is a possible result of increased high latitude temperatures due to increasing anthropogenic greenhouse gas emissions. Using an atmosphere-ocean general circulation model (AOGCM), we investigate the effects of the removal of the ice sheet on atmospheric temperatures, circulation, and precipitation. We find that locally over Greenland, there is a warming associated directly with the altitude change in winter, and the altitude and albedo change in summer. Outside of Greenland, the largest signal is a cooling over the Barents sea in winter. We attribute this cooling to a decrease in poleward heat transport in the region due to changes to the time mean circulation and eddies, and interaction with sea-ice. The simulated climate is used to force a vegetation model and an ice-sheet model. We find that the Greenland climate in the absence of an ice sheet supports the growth of trees in southern Greenland, and grass in central Greenland. We find that the ice sheet is likely to regrow following a melting of the Greenland ice sheet, the subsequent rebound of its bedrock, and a return to present day atmospheric CO₂ concentrations. This regrowth is due to the high altitude bedrock in eastern Greenland which allows the growth of glaciers which develop into an ice sheet.     

On the role of flux adjustments in an idealized coupled climate model

 The effect of employing flux adjustments on the climatic response of an idealized coupled model to an imposed radiative forcing is investigated with two coupled models, one of which employs flux adjustments. A linear reduction (to the planetary longwave flux) of 4 W/m² is applied over a 70 y period and held constant thereafter. Similar model responses are found (during the initial 70 y period) for global-scale diagnostics of hemispheric air temperature due to the nearly linear surface-air temperature response to the radiative forcing. Significant regional scale differences do exist, however. As the perturbation away from the present climate grows, basin-scale diagnostics (such as meridional overturning rates) begin to diverge between flux adjusted and non-flux adjusted models. Once the imposed radiative forcing is held constant, differences in global mean air temperature of up to 0.5 °C are found, with large regional-scale differences in air temperature and overturning rates within the North Atlantic and Southern Ocean. Two additional experiments with the flux adjusted model (beginning from points further along the control integration) suggest that the elimination of much of the coupling shock before the radiative forcing is applied leads to results slightly closer to the non-flux adjusted case, although large differences still persist. In particular a dipole structure indicating an enhanced warming within the Pacific sector of the Southern Ocean, and cooling within the Atlantic sector is not reproduced by the flux adjusted models. This disparity is intimately linked to the Southern Ocean overturning cell along with the flux adjustments employed as well as the drift arising from coupling shock. If a similar form of sensitivity exists in more realistic coupled models, our results suggest: (1) perturbation experiments should not be undertaken until after the coupled model control experiment is carried out for several hundred years (thereby minimizing the coupling shock); (2) care should be exercised in the interpretation of regional-scale results (over the ocean) in coupled models which employ flux adjustments; (3) care should also be taken in interpreting even global-scale diagnostics in flux adjusted models for large perturbations about the present climate. Received: 15 November 1996 / Accepted: 4 June 1997     

Global and regional variability in a coupled AOGCM

 The variability of near surface temperature on global and regional spatial scales and interannual time scales from a 1000 year control integration of the Hadley Centre coupled model (HADCM2-CTL) are compared with the observational record of surface temperature. The model succeeds in reproducing the observed patterns of natural variability, with high variability over the northern continents and low variability over much of the tropics. The model global mean variability has similar strength to observed global mean variability on time scales less than 20 years. The warming seen in the historical record is outside the range of natural variability as simulated in HADCM2-CTL. The model has El-Ni?o/Southern Oscillation (ENSO)-like behaviour with a central Pacific, peak to peak, strength of approximately 3 K. Changes in near surface temperature in the central Pacific are strongly correlated with changes in near surface temperature over most of the tropics, large regions of the extra-tropics and changes in tropical ocean upper 250 m heat content. Tropospheric temperature changes and tropical surface pressure changes are also strongly correlated with changes in the central Pacific surface temperature. Oceanic regions show significant departures from an AR1 or first order Markov behaviour in the Northwest Atlantic, Northwest Pacific and Arctic oceans. The Northwest Atlantic region has large amounts of variability over periods greater than 50 years. This variability is associated with a jump in the strength of North Atlantic meridional stream function. The spectra of the Western European and Continental US land regions are not significantly different from an AR1 process. The flow through the Drake Passage has an interannual standard deviation of approximately 2.5 Sv with significant departures from an AR1 process at time scales greater than 40 years. Winter northern hemispheric 500 hPa geopotential height shows some evidence of multiple regimes but no year to year persistence of these regimes. Received: 31 January 1996/Accepted: 22 July 1996     

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.     

Intercomparison of interannual variability and potential predictability: an AMIP diagnostic subproject

 The inter-annual variability and potential predictability of 850 hPa temperature (T ₈₅₀), 500 hPa geopotential (φ₅₀₀) and 300 hPa stream function (ψ₃₀₀) simulated by the models participating in the Atmospheric Model Intercomparison Project (AMIP) are examined. The total inter-annual variability is partitioned into a potentially predictable component which arises from the forcing implied by the prescribed SST and sea-ice evolution, or from sources internal to the simulated climate, and an unpredictable low frequency component induced by “weather noise”. There is wide variation in the ability to simulate observed inter-annual variability, both total and weather-noise induced. A majority of models under simulate seasonal mean φ₅₀₀ variability in DJF and JJA and over simulate ψ₃₀₀ variability in JJA. All but one model simulates less T ₈₅₀ total inter-annual variability than in the analysed data. There is little apparent connection between gross model characteristics and the corresponding ability to simulate observed variability, with the possible exceptions of resolution. Received: 7 July 1996 / Accepted: 8 January 1998     

Linear response functions to project contributions to future sea level

We propose linear response functions to separately estimate the sea-level contributions of thermal expansion and solid ice discharge from Greenland and Antarctica. The response function formalism introduces a time-dependence which allows for future rates of sea-level rise to be influenced by past climate variations. We find that this time-dependence is of the same functional type, R(t) ～ t ^α, for each of the three subsystems considered here. The validity of the approach is assessed by comparing the sea-level estimates obtained via the response functions to projections from comprehensive models. The pure vertical diffusion case in one dimension, corresponding to α =  ?0.5, is a valid approximation for thermal expansion within the ocean up to the middle of the twenty first century for all Representative Concentration Pathways. The approximation is significantly improved for α =  ? 0.7. For the solid ice discharge from Greenland we find an optimal value of α =  ?0.7. Different from earlier studies we conclude that solid ice discharge from Greenland due to dynamic thinning is bounded by 0.42 m sea-level equivalent. Ice discharge induced by surface warming on Antarctica is best captured by a positive value of α = 0.1 which reflects the fact that ice loss increases with the cumulative amount of heat available for softening the ice in our model.