Transitivity properties of surface temperature and ice cover in the CCM1

A frequently made assumption in simple models of long-term climatic behavior (e.g. ‘energy balance’ models) is that, owing to instability engendered by the ice-albedo feedback, the climatic system (surface temperature, in particular) can exhibit multiple steady states within the paleoclimatologically observed range of temperature. Here we show that, for a more comprehensive model than an energy balance model, the CCM1 general circulation model, such a bimodality does not exist for present values of the solar constant and atmospheric CO₂ even if one excludes the existence of a seasonal cycle of radiative forcing. Thus, no evidence is found to support the whole class of ice-age theories tuned to present CO₂ levels that depend on this bimodality. As a corollary, support is found for the idea that surficial temperature and snow-sea-ice cover are essentially ‘slaved’, fast-response, climatic variables that equilibrate uniquely with the prescribed external forcing and the slow-response climatic variables (e.g. the ice sheets and deep ocean state). It is also implied that, although care should be exercised in the choice of initial conditions to minimize computer time, the selection of initial conditions is not likely to affect the final outcome of general circulation model studies of climate for fixed (near-present) values of the solar constant and CO₂ forcing.     

Multiforced statistical assessments of greenhouse-gas-induced surface air temperature change 1890–1985

Based on univariate correlation and coherence analyses and considering the physical basis of the relationships, a simple multiforced (multiple) statistical concept is used which correlates observational climatic time series simultaneously with volcanic, solar, ENSO, and the anthropogenic greenhouse gases forcing. This is appropriate to remove some natural climate noise in the observed data and to evaluate the components (signals) possibly due to the anthropogenic greenhouse gas forcing (CO₂, or equivalent CO₂ implying additional gases) during industrial time. In this paper, we apply this technique to 100 global box data time series 1890–1985, of the surface air temperature, using observed data from Hansen and Lebedeff. The results are presented in terms of latitudinal-seasonal and regional trends, where the observed trend patterns are compared with the hypothetical signals (statistical assessments) possibly due to anthropogenic greenhouse forcing. These latter signals can be amplified to enable a comparison with corresponding results from general circulation model (GCM) CO₂ doubling experiments. These observed-statistical assessments lead to results which are, at least qualitatively and in respect to the zonal mean temperatures, very similar to some GCM experiments indicating the maximum CO₂ doubling signals (statistical assessment > 12 K) in the arctic winter. However, these signals are moderate in the tropics and in the Southern Hemisphere (global average 2.8–4.4 K). As far as the industrial signals are concerned (observed period) these signals are somewhat larger (maximum 7 K, global average 0.5–0.9 K) than the observed trends (maximum 5 K, global average 0.5 K). Phase shifts of cause and effect may amplify these signals but are very uncertain.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     

Modelling the equilibrium and transient responses of global temperature to past and future trace gas concentrations

Recent works with energy balance climate models and oceanic general circulation models have assessed the potential role of the world ocean for climatic changes on a decadal to secular time scale. This scientific challenge is illustrated by estimating the response of the global temperature to changes in trace gas concentration from the pre-industrial epoch to the middle of the next century. A simple energetic formulation is given to estimate the effect on global equilibrium temperature of a fixed instantaneous radiative forcing and of a time-dependent radiative forcing. An atmospheric energy balance model couple to a box-advection-diffusion ocean model is then used to estimate the past and future global climalic transient response to trace-gas concentration changes. The time-dependent radiative perturbation is estimated from a revised approximate radiative parameterization, and the recent reference set of trace gas scenarios proposed by Wuebbles et al. (1984) are adopted as standard scenarios. Similar computations for the past and future have recently been undertaken by Wigley (1985), but using a purely diffusive ocean and slightly different trace gas scenarios. The skill of the socalled standard experiment is finally assessed by examining the model sensitivity of different parameters such as the equilibrium surface air temperature change for a doubled CO₂ concentration [T _ae (2×CO₂)], the heat exchange with the deeper ocean and the trace gas scenarios. For T _ae (2×CO₂) between 1 K and 5 K, the following main results are obtained: (i) for a pre-industrial CO₂, concentration of 270 ppmv, the surface air warming between 1850 and 1980 ranges between 0.4 and 1.4 K (if a pre-industrial CO₂ concentration of 290 ppmv is chosen, the range is between 0.3 and 1 K); (ii) by comparison with the instantaneous equilibrium computations, the deeper ocean inertia induces a delay which amounts to between 6 years [for lower T _ae (2×CO₂)] and 23 years [for higher Tae(2×CO₂)] in 1980; (iii) for the standard future CO₂ and other trace gas scenarios of Wuebbles et al., the surface air warming between 1980 and 2050 is calculated to range between 0.9 and 3.4 K, with a delay amounting to between 7 years and 32 years in 2050 when compared to equilibrium computations.     

Farm-Level Adaptation to Climatic Variability and Change: Crop Diversification in the Canadian Prairies

Among other foci, recent research on adaptation to climatic variability and change has sought to evaluate the merit of adaptation generally, as well as the suitability of particular adaptations. Additionally, there is a need to better understand the likely uptake of adaptations. For example, diversification is one adaptation that has been identified as a potential farm-level response to climatic variability and change, but its adoption by farmers for this reason is not well understood. This paper serves two purposes. The first is to document the adoption of crop diversification in Canadian prairie agriculture for the period 1994–2002, reflect upon its strengths and limitations for managing a variety of risks, including climatic ones, and gauge its likely adoption by producers in response to anticipated climate change. The second purpose is to draw on this case to refine our current understanding of climate change adaptation more generally. Based upon data from over 15 000 operations, it was determined that individual farms have become more specialized in their cropping patterns since 1994, and this trend is unlikely to change in the immediate future, notwithstanding anticipated climate change and the known risk-reducing benefits of crop diversification. More broadly, the analysis suggests that suitable and even possible climate change adaptations need to be more rigorously assessed in order to understand their wider strengths and limitations.     

Adapting stochastic weather generation algorithms for climate change studies

While large-scale climate models (GCMs) are in principle the most appropriate tools for predicting climate changes, at present little confidence can be placed in the details of their projections. Use of tools such as crop simulation models for investigation of potential impacts of climatic change requires daily data pertaining to small spatial scales, not the monthly-averaged and large-scale information typically available from the GCMs. A method is presented to adapt stochastic weather generation models, describing daily weather variations in the present-day climate at particular locations, to generate synthetic daily time series consistent with assumed future climates. These assumed climates are specified in terms of the commonly available monthly means and variances of temperature and precipitation, including time-dependent (so-called transient) climate changes. Unlike the usual practice of applying assumed changes in mean values to historically observed data, simulation of meteorological time series also exhibiting changes in variability is possible. Considerable freedom in climate change scenario construction is allowed. The results are suitable for investigating agricultural and other impacts of a variety of hypothetical climate changes specified in terms of monthly-averaged statistics.     

Monte Carlo climate change forecasts with a global coupled ocean-atmosphere model

Four time-dependent greenhouse warming experiments were performed with the same global coupled atmosphere-ocean model, but with each simulation using initial conditions from different snapshots of the control run climate. The radiative forcing — the increase in equivalent CO₂ concentrations from 1985–2035 specified in the Intergovernmental Panel on Climate Change (IPCC) scenario A — was identical in all four 50-year integrations. This approach to climate change experiments is called the Monte Carlo technique and is analogous to a similar experimental set-up used in the field of extended range weather forecasting. Despite the limitation of a very small sample size, this approach enables the estimation of both a mean response and the between-experiment variability, information which is not available from a single integration. The use of multiple realizations provides insights into the stability of the response, both spatially, seasonally and in terms of different climate variables. The results indicate that the time evolution of the global mean warming signal is strongly dependent on the initial state of the climate system. While the individual members of the ensemble show considerable variation in the pattern and amplitude of near-surface temperature change after 50 years, the ensemble mean climate change pattern closely resembles that obtained in a 100-year integration performed with the same model. In global mean terms, the climate change signals for near surface temperature, the hydrological cycle and sea level significantly exceed the variability among the members of the ensemble. Due to the high internal variability of the modelled climate system, the estimated detection time of the global mean temperature change signal is uncertain by at least one decade. While the ensemble mean surface temperature and sea level fields show regionally significant responses to greenhouse-gas forcing, it is not possible to identify a significant response in the precipitation and soil moisture fields, variables which are spatially noisy and characterized by large variability between the individual integrations.     

In Search of Warmer Climates? The Impact of Climate Change on Flows of British Tourists

This paper investigates the impact of climate change on the chosendestinations of Britishtourists. Destinations are characterised in terms of attractors includingclimate variables, traveland accommodation costs. These variables are used to explain the currentobserved pattern ofoverseas travel in terms of a model based upon the idea of utilitymaximisation. The approachpermits the trade-offs between climate and holiday expenditure to be analysedand effectivelyidentifies the optimal climate for generating tourism. The findings are usedto predict the impactof various climate change scenarios on popular tourist destinations.     

General circulation model experiments with surface albedo changes

In 1975, Charney proposed a biogeophysical feedback mechanism to partly explain the droughts that occur in desert border areas. He showed that a perturbation of albedo (due to a natural or anthropogenic decrease of vegetation) can be unstable and lead to a variation of precipitation in the region where albedo is changed.Several numerical experiments have been achieved with general circulation models to study the sensitivity of climate to surface albedo. We compare the GLAS and LMD model results for the Sahel. For all models, rainfall decreases when albedo increases and net radiative heating of soil decreases. We show the variations of circulation simulated by the LMD model that we obtain when albedo is increased. These changes are compared to the weakening of Easterly Jet at 200 mb observed during dry years.     

Can market mechanisms ameliorate the effects of long-term climate change?

Land prices, insurance rates, future markets, mortgage terms, and other market mechanisms may be expected to guide short-term economic response to climatic change as its effects become apparent to investors. On the other hand, the pervasive influence of discount rates on investment decisions makes it unlikely that the market will give satisfactory guidance to investments that must be undertaken long before the appearance of the climatic effect they are intended to mitigate. For this reason, only government is likely to undertake such long-term investments as large civil works intended to modify hydrology (irrigation, seawall dikes), and research and development in agriculture or other technologies which need to be adapted to new climatic conditions, while the effects of climate change are still distant and uncertain.Decisions regarding very large macroprojects, as well as decisions which determine the siting of installations that have long-term consequences for the environment (e.g., dump sites for disposal of long-term hazardous wastes) should carefully consider the effects of climate change regardless of current market signals. Strategic planners for government and industry should take steps now to identify those decisions for which planning is now beginning, and which need to take into account the effects of long-term climate change.     

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     

Sensitivity of glaciation to initial snow cover,CO2, snow albedo,and oceanic roughness in the NCAR CCM

We present results from numerical experiments made with a GCM, the NCAR CCM1, that were designed to estimate the annual balance between snow-fall accumulation and ablation for geographically important land regions for a variety of conditions. We also attempt to assess the reliability of these results by investigating model sensitivity to changes in prescribed physical parameters. Experiments were run with an initial imposition of 1 m of (midwinter) snowcover over all northern hemisphere land points. Over Alaska, western Canada, Siberia, and the Tibetan Plateau the model tended to retain this snow cover through the summer and in some cases increase its depth as well. We define these regions as glaciation sensitive and note some correspondence between them and source regions for the Pleistocene ice sheets. An experiment with greatly reduced CO₂ (100 ppm) showed a tendency towards spontaneous glaciation, i.e., the model remained snow-covered throughout the summer over the same geographic regions noted above. With 200 ppm CO₂ (roughly equal to values at the last glacial maximum), snow cover over these regions did not quite survive the summer on a consistent basis. Combining 200 ppm CO₂ and 1 m of initial northern hemisphere snow cover yielded glaciation-sensitive conditions, agreeing remarkably well with locations undergoing glaciation during the Pleistocene. To assess the reliability of these results, we have determined minimal model uncertainty by varying two of the empirical coefficients in the model within physically plausible ranges. In one case surface roughness of all ocean gridpoints was reduced by an order of magnitude, leading to local 10% reductions in precipitation (snowfall), a change hard to distinguish from inherent model variability. In the other case, the fraction of a land grid square assumed to be occupied by snow cover for albedo purposes was varied from one-half to unity. Large changes occurred in the degree of summer melting, and in some cases the sign of the net balance changed as fractional snow cover was changed. We conclude that the model may be able to reveal regions sensitive to glaciation, but that it cannot yield a reliable quantitative computation of the magnitude of the net snow accumulation that can be implicitly or explicitly integrated through time.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. Dilmenil     

Sea-ice anomalies observed in the Greenland and Labrador seas during 1901–1984 and their relation to an interdecadal Arctic climate cycle

Two independent ice data sets from the Greenland and Labrador Seas have been analyzed for the purpose of characterizing interannual and decadal time scale sea-ice extent anomalies during this century. Sea-ice concentration data for the 1953–1984 period revealed the presence of a large positive anomaly in the Greenland Sea during the 1960s which coincided with the great salinity anomaly, an upper-ocean low-salinity water mass that was observed to travel cyclonically around the northern North Atlantic during 1968–1982. This ice anomaly as well as several smaller ones propagated into the Labrador Sea and then across to the Labrador and east Newfoundland coast, over a period of 3 to 5 years. A complex empirical orthogonal function analysis of the same data also confirmed this propagation phenomenon. An inverse relation between sea-ice and salinity anomalies in the Greenland-Labrador Sea region was also generally found. An analysis of spring and summer ice-limit data obtained from Danish Meteorological Institute charts for the period 1901–1956 indicated the presence of heavy ice conditions (i.e., positive ice anomalies) in the Greenland Sea during 1902–1920 and in the late 1940s, and generally negative ice anomalies during the 1920s and 1930s. Only limited evidence of the propagation of Greenland Sea ice anomalies into the Labrador Sea was observed, however, probably because the data were from the ice-melt seasons. On the other hand, several large ice anomalies in the Greenland Sea occurred 2–3 years after large runoffs (in the early 1930s and the late 1940s) from northern Canada into the western Arctic Ocean. Similarly, a large runoff into the Arctic during 1964–1966 preceded the large Greenland Sea ice anomaly of the 1960s. These facts, together with recent evidence of climatic jumps in the Northern Hemisphere tropospheric circulation, suggest the existence of an interdecadal self-sustained climate cycle in the Arctic. In the Greenland Sea, this cycle is characterized by a state of large sea-ice extent overlying an upper layer of cool, relatively fresh water that does not convectively overturn, which alternates every 10–15 years with a state of small sea-ice extent and relatively warm saline surface water that frequently overturns.Dedicated to Robert W. Stewart on the occasion of his retirement     

Impacts on Global Ozone and Climate from Use and Emission of 2,2-Dichloro-1,1,1-Trifluoroethane (HCFC-123)

Analyses of emissions, and consequent chlorine loading, show that projected use of 2,2-dichloro-1,1,1-trifluoroethane (HCFC-123) will result in a virtually indiscernible impact on stratospheric ozone. Parametric scenarios uphold this conclusion, even for extreme levels of emissions far exceeding those of current technologies and practices. Additional scenarios reaffirm the conclusion for continued use – beyond the scheduled phaseout date – as a refrigerant in closed systems. By contrast, use of this compound offers unique opportunities to reduce global warming. Moreover, time-dependent analyses show that the minimal contribution to stratospheric chlorine from HCFC-123 emissions will not peak until more than a decade after the residual peaks of chlorine and bromine, from prior chlorofluorocarbon and halon releases, subside. While no single index exists to compare the relative demerits of ozone depletion and climate change, three conclusions are clear. First, reversal of the buildup of bromine and chlorine (i.e., healing of the ozone layer) is underway and progressing on target, while sufficient practical remedies for global climate change are far more difficult. Second, the analyses show that phaseout of all chlorinated, and conceptually – but much less probably – all brominated, compounds of anthropogenic origin targets some compounds that provide environmental benefits. Most chlorinated and brominated compounds do warrant phaseout; the exceptions are those with very short atmospheric lifetimes, and consequent low ozone depletion potential (ODP), that also offer offsetting environmental benefits. And third, since new global environmental concerns may, and probably will, be identified in the future, a more scientific approach is needed to determine environmental acceptability or rejection.     

A study of the climatic regimes of the Pleistocene using a stochastic resonance model

A global energy balance model employing the stochastic resonance mechanism, previously used to explain the climatic variability of the late Pleistocene, has now been extended to account for the climatic variations over the full Pleistocene. The possibility that extremely long-term changes (of the order of millions of years) in the boundary conditions of the climate system have altered the response of the Pleistocene climate to the external orbital forcing has been investigated. It is shown that, by slowly changing the only free parameter of the model, the system can undergo a pitchfork bifurcation. The bifurcation point separates a linear regime (identified with the early Pleistocene climate) from a strongly nonlinear regime (the late Pleistocene) where the stochastic resonance mechanism produces rapid and symmetric transitions between the two stable steady states of the system. The main differences in the dynamic features of the two regimes are the change in amplitude of the oscillations, the relative importance of the stochastic forcing, the change in shape of the probability distribution, and the corresponding change in the power centered around the 100000 year cycle: in qualitative agreement with the observed geological record. With the introduction of the external orbital forcing, now spectrally complete and included without requiring any additional hypothesis, the model reproduces the previous results, namely the good correlation with the isotopic record, the appearance of the dominant spectral peaks, as well as the redness of the power spectrum. In particular, it is shown that the orbital forcing in eccentricity acts as a pacemaker of the major glacial cycles of the late Pleistocene through the mechanism of stochastic resonance. A stochastic sensitivity analysis is then applied to validate the significance of the results and to investigate the predictability of the climate system over the time-scales of the orbital cycles.     

Reproducibility of seasonal ensemble integrations with ECMWF GCM and its association with ENSO

Summary This study evaluates seasonal climate potential predictability with a dataset of nine-member ensemble seasonal integrations produced by the ECMWF GCM for the ERA-15 period (1979–1993). The methodology used here is the measure of the ensemble reproducibility for a particular season defined by Yang et al (1998). High reproducibility reflects the dominant role of the lower boundary forcing in seasonal climate anomaly, indicating good potential predictability. Spatial patterns of the reproducibility for selected variables are documented, which exhibit obvious regionality and seasonality. Such variables are always highly reproducible over most of the tropical regions. Over the northern extratropics, primary reproducible information, taking the 500hPa geopotential height for example, is found over the PNA region during winter while over most of Asia during summer. Winter has the largest reproducible area; autumn has the lowest, while summer and spring are in between. Association of the reproducibility with ENSO events was examined regionally. Internal variances due to ensemble spread were broken down for individual years, and the reproducibility was computed for four categories: El Niño, La Niña, ENSO, and non-ENSO years. The reproducibility during winter especially over the PNA region is insensitive to ENSO events, with exceptions over the tropical western Pacific, central Siberia and Western Europe. Contrarily, ENSO events have significant impacts on the reproducibility over the southwest USA and most of Asian monsoon region during summer. These results suggest that ENSO events may not be helpful to the seasonal climate predictability over the PNA region during winter, but they may increase predictable information over many regions of the northern extratropical continents during summer.     

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     

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     

On Trajectory Curvature as a Selection Criterion for Valid Lagrangian Stochastic Dispersion Models

Recently Wilson and Flesch (Boundary-Layer Meteorology, 84, 411-426, 1997) suggested that the average increment d z to the orientation = arctan(w/u) of the Lagrangian velocity-fluctuation vector can be used to distinguish the better Lagrangian stochastic models within the well-mixed class. Here it is demonstrated that the specification of d z constitutes neither a sufficient or universally applicable criterion to distinguish the better Lagrangian stochastic models within the well-mixed class. The hypothesis made by Wilson and Flesch that Lagrangian stochastic models with /P_E irrotational are zero-spin models, having d z=0, is proven     

Mechanisms of shrubland expansion: land use,climate or CO2?

Encroachment of trees and shrubs into grasslands and the thicketization of savannas has occurred worldwide over the past century. These changes in vegetation structure are potentially relevant to climatic change as they may be indicative of historical shifts in climate and as they may influence biophysical aspects of land surface-atmosphere interactions and alter carbon and nitrogen cycles. Traditional explanations offered to account for the historic displacement of grasses by woody plants in many arid and semi-arid ecosystems have centered around changes in climatic, livestock grazing and fire regimes. More recently, it has been suggested that the increase in atmospheric CO₂ since the industrial revolution has been the driving force. In this paper we evaluate the CO₂ enrichment hypotheses and argue that historic, positive correlations between woody plant expansion and atmospheric CO₂ are not cause and effect.Please direct all correspondence to the senior author.