Geo-Engineering to Confine Climate Change: Is it at all Feasible?

Paul Crutzen (2006) has suggested a research initiative to consider whether it would be feasible to artificially enhance the albedo of the planet Earth to counteract greenhouse warming. The enhancement of albedo would be achieved by intentionally injecting sulfur into the stratosphere. The rational for proposing the experiment is the observed cooling of the atmosphere following the recent major volcanic eruptions by El Chichon in 1984 and Mount Pinatubo in 1991 (Hansen et al., 1992). Although I am principally not against a research initiative to study such a potential experiment, I do have important reservations concerning its general feasibility. And its potential feasibility, I believe, must be the key motivation for embarking on such a study. Here I will bring up three major issues, which must be more thoroughly understood before any geo-engineering of climate could be considered, if at all. The three issues are (i) the lack of accuracy in climate prediction, (ii) the huge difference in timescale between the effect of greenhouse gases and the effect of aerosols and (iii) serious environmental problems which may be caused by high carbon dioxide concentration irrespective of the warming of the climate.

Origin of regional climate differences: role of boundary conditions and model formulation in two GCMs

Model differences in projections of extratropical regional climate change due to increasing greenhouse gases are investigated using two atmospheric general circulation models (AGCMs): ECHAM4 (Max Planck Institute, version 4) and CCM3 (National Center for Atmospheric Research Community Climate Model version 3). Sea-surface temperature (SST) fields calculated from observations and coupled versions of the two models are used to force each AGCM in experiments based on time-slice methodology. Results from the forced AGCMs are then compared to coupled model results from the Coupled Model Intercomparison Project 2 (CMIP2) database. The time-slice methodology is verified by showing that the response of each model to doubled CO2 and SST forcing from the CMIP2 experiments is consistent with the results of the coupled GCMs. The differences in the responses of the models are attributed to (1) the different tropical SST warmings in the coupled simulations and (2) the different atmospheric model responses to the same tropical SST warmings. Both are found to have important contributions to differences in implied Northern Hemisphere (NH) winter extratropical regional 500 mb height and tropical precipitation climate changes. Forced teleconnection patterns from tropical SST differences are primarily responsible for sensitivity differences in the extratropical North Pacific, but have relatively little impact on the North Atlantic. There are also significant differences in the extratropical response of the models to the same tropical SST anomalies due to differences in numerical and physical parameterizations. Differences due to parameterizations dominate in the North Atlantic. Differences in the control climates of the two coupled models from the current climate, in particular for the coupled model containing CCM3, are also demonstrated to be important in leading to differences in extratropical regional sensitivity.

Why is the global warming proceeding much slower than expected?

Upper air observations from radiosondes and microwave satellite instruments does not indicate any global warming during the last 19 years, contrary to surface measurements, where a warming trend is supposedly being found. This result is somewhat difficult to reconcile, since climate model experiments do indicate a reverse trend, namely, that upper tropospheric air should warm faster than the surface. To contribute toward an understanding of this difficulty, we have here undertaken some specific experiments to study the effect on climate due to the decrease in stratospheric ozone and the Mount Pinatubo eruption in 1991. The associated forcing was added to the forcing from greenhouse gases, sulfate aerosols (direct and indirect effect), and tropospheric ozone, which was investigated in a separate series of experiments. Furthermore, we have undertaken an ensemble study in order to explore the natural variability of an advanced climate model exposed to such a forcing over 19 years. The result shows that the reduction of stratospheric ozone cools not only the lower stratosphere but also the troposphere, in particular, the upper and middle part. In the upper troposphere the cooling from stratospheric ozone leads to a significant reduction of greenhouse warming. The modeled stratospheric aerosols from Mount Pinatubo generate a climate response (stratospheric warming and tropospheric cooling) in good agreement with microwave satellite measurements. Finally, analysis of a series of experiments with both stratospheric ozone and the Mount Pinatubo effect shows considerable variability in climate response, suggesting that an evolution having no warming in the period is as likely as another evolution showing modest warming. However, the observed trend of no warming in the midtroposphere and clear warming at the surface is not found in the model simulations.

Transient climate change simulations with a coupled atmosphere–ocean GCM including the tropospheric sulfur cycle

The time-dependent climate response to changing concentrations of greenhouse gases and sulfate aerosols is studied using a coupled general circulation model of the atmosphere and the ocean (ECHAM4/OPYC3). The concentrations of the well-mixed greenhouse gases like CO2, CH4, N2O, and CFCs are prescribed for the past (1860–1990) and projected into the future according to International Panel on Climate Change (IPCC) scenario IS92a. In addition, the space–time distribution of tropospheric ozone is prescribed, and the tropospheric sulfur cycle is calculated within the coupled model using sulfur emissions of the past and projected into the future (IS92a). The radiative impact of the aerosols is considered via both the direct and the indirect (i.e., through cloud albedo) effect. It is shown that the simulated trend in sulfate deposition since the end of the last century is broadly consistent with ice core measurements, and the calculated radiative forcings from preindustrial to present time are within the uncertainty range estimated by IPCC. Three climate perturbation experiments are performed, applying different forcing mechanisms, and the results are compared with those obtained from a 300-yr unforced control experiment. As in previous experiments, the climate response is similar, but weaker, if aerosol effects are included in addition to greenhouse gases. One notable difference to previous experiments is that the strength of the Indian summer monsoon is not fundamentally affected by the inclusion of aerosol effects. Although the monsoon is damped compared to a greenhouse gas only experiment, it is still more vigorous than in the control experiment. This different behavior, compared to previous studies, is the result of the different land–sea distribution of aerosol forcing. Somewhat unexpected, the intensity of the global hydrological cycle becomes weaker in a warmer climate if both direct and indirect aerosol effects are included in addition to the greenhouse gases. This can be related to anomalous net radiative cooling of the earth’s surface through aerosols, which is balanced by reduced turbulent transfer of both sensible and latent heat from the surface to the atmosphere.

Impact of global warming on the interannual and interdecadal climate modes in a coupled GCM

In this study, we investigated the impact of global warming on the variabilities of large-scale interannual and interdecadal climate modes and teleconnection patterns with two long-term integrations of the coupled general circulation model of ECHAM4/OPYC3 at the Max-Planck-Institute for Meteorology, Hamburg. One is the control (CTRL) run with fixed present-day concentrations of greenhouse gases. The other experiment is a simulation of transient greenhouse warming, named GHG run. In the GHG run the averaged geopotential height at 500 hPa is increased significantly, and a negative phase of the Pacific/North American (PNA) teleconnection-like distribution pattern is intensified. The standard deviation over the tropics (high latitudes) is enhanced (reduced) on the interdecadal time scales and reduced (enhanced) on the interannual time scales in the GHG run. Except for an interdecadal mode related to the Southern Oscillation (SO) in the GHG run, the spatial variation patterns are similar for different (interannual + interdecadal, interannual, and interdecadal) time scales in the GHG and CTRL runs. Spatial distributions of the teleconnection patterns on the interannual and interdecadal time scales in the GHG run are also similar to those in the CTRL run. But some teleconnection patterns show linear trends and changes of variances and frequencies in the GHG run. Apart from the positive linear trend of the SO, the interdecadal modulation to the El Niño/SO cycle is enhanced during the GHG 2040 ∼ 2099. This is the result of an enhancement of the Walker circulation during that period. La Niña events intensify and El Niño events relatively weaken during the GHG 2070 ∼ 2090. It is interesting to note that with increasing greenhouse gas concentrations the relation between the SO and the PNA pattern is reversed significantly from a negative to a positive correlation on the interdecadal time scales and weakened on the interannual time scales. This suggests that the increase of the greenhouse gas concentrations will trigger the nonstationary correlation between the SO and the PNA pattern both on the interdecadal and interannual time scales.

The effect of nonlinearity in CO2 heating rates on the attribution of stratospheric ozone and temperature changes

An analysis of the attribution of past and future changes in stratospheric ozone and temperature to anthropogenic forcings is presented. The analysis is an extension of the study of Shepherd and Jonsson (2008) who analyzed chemistry-climate simulations from the Canadian Middle Atmosphere Model (CMAM) and attributed both past and future changes to changes in the external forcings, i.e. the abundances of ozone-depleting substances (ODS) and well-mixed greenhouse gases. The current study is based on a new CMAM dataset and includes two important changes. First, we account for the nonlinear radiative response to changes in CO2. It is shown that over centennial time scales the radiative response in the upper stratosphere to CO2 changes is significantly nonlinear and that failure to account for this effect leads to a significant error in the attribution. To our knowledge this nonlinearity has not been considered before in attribution analysis, including multiple linear regression studies. For the regression analysis presented here the nonlinearity was taken into account by using CO2 heating rate, rather than CO2 abundance, as the explanatory variable. This approach yields considerable corrections to the results of the previous study and can be recommended to other researchers. Second, an error in the way the CO2 forcing changes are implemented in the CMAM was corrected, which significantly affects the results for the recent past. As the radiation scheme, based on Fomichev et al. (1998), is used in several other models we provide some description of the problem and how it was fixed.

Dynamics, stratospheric ozone, and climate change

Dynamics affects the distribution and abundance of stratospheric ozone directly through transport of ozone itself and indirectly through its effect on ozone chemistry via temperature and transport of other chemical species. Dynamical processes must be considered in order to understand past ozone changes, especially in the northern hemisphere where there appears to be significant low-frequency variability which can look “trend-like” on decadal time scales. A major challenge is to quantify the predictable, or deterministic, component of past ozone changes. Over the coming century, changes in climate will affect the expected recovery of ozone. For policy reasons it is important to be able to distinguish and separately attribute the effects of ozone-depleting substances and greenhouse gases on both ozone and climate. While the radiative-chemical effects can be relatively easily identified, this is not so evident for dynamics — yet dynamical changes (e.g., changes in the Brewer-Dobson circulation) could have a first-order effect on ozone over particular regions. Understanding the predictability and robustness of such dynamical changes represents another major challenge. Chemistry-climate models have recently emerged as useful tools for addressing these questions, as they provide a self-consistent representation of dynamical aspects of climate and their coupling to ozone chemistry. We can expect such models to play an increasingly central role in the study of ozone and climate in the future, analogous to the central role of global climate models in the study of tropospheric climate change.

Multimodel projections of stratospheric ozone in the 21st century

Simulations from eleven coupled chemistry-climate models (CCMs) employing nearly identical forcings have been used to project the evolution of stratospheric ozone throughout the 21st century. The model-to-model agreement in projected temperature trends is good, and all CCMs predict continued, global mean cooling of the stratosphere over the next 5 decades, increasing from around 0.25 K/decade at 50 hPa to around 1 K/ decade at 1 hPa under the Intergovernmental Panel on Climate Change (IPCC) Special Report on Emissions Scenarios (SRES) A1B scenario. In general, the simulated ozone evolution is mainly determined by decreases in halogen concentrations and continued cooling of the global stratosphere due to increases in greenhouse gases (GHGs). Column ozone is projected to increase as stratospheric halogen concentrations return to 1980s levels. Because of ozone increases in the middle and upper stratosphere due to GHGinduced cooling, total ozone averaged over midlatitudes, outside the polar regions, and globally, is projected to increase to 1980 values between 2035 and 2050 and before lower stratospheric halogen amounts decrease to 1980 values. In the polar regions the CCMs simulate small temperature trends in the first and second half of the 21st century in midwinter. Differences in stratospheric inorganic chlorine (Cly) among the CCMs are key to diagnosing the intermodel differences in simulated ozone recovery, in particular in the Antarctic. It is found that there are substantial quantitative differences in the simulated Cly, with the October mean Antarctic Cly peak value varying from less than 2 ppb to over 3.5 ppb in the CCMs, and the date at which the Cly returns to 1980 values varying from before 2030 to after 2050. There is a similar variation in the timing of recovery of Antarctic springtime column ozone back to 1980 values. As most models underestimate peak Cly near 2000, ozone recovery in the Antarctic could occur even later, between 2060 and 2070. In the Arctic the column ozone increase in spring does not follow halogen decreases as closely as in the Antarctic, reaching 1980 values before Arctic halogen amounts decrease to 1980 values and before the Antarctic. None of the CCMs predict future large decreases in the Arctic column ozone. By 2100, total column ozone is projected to be substantially above 1980 values in all regions except in the tropics.

Detection of anthropogenic climate change using a fingerprint method

A fingerprint method for detecting anthropogenic climate change is applied to new simulations with a coupled ocean-atmosphere general circulation model (CGCM) forced by increasing concentrations of greenhouse gases and aerosols covering the years 1880 to 2050. In addition to the anthropogenic climate change signal, the space-time structure of the natural climate variability for near-surface temperatures is estimated from instrumental data over the last 134 years and two 1000 year simulations with CGCMs. The estimates are compared with paleoclimate data over 570 years. The space-time information on both the signal and the noise is used to maximize the signal-to-noise ratio of a detection variable obtained by applying an optimal filter (fingerprint) to the observed data. The inclusion of aerosols slows the predicted future warming. The probability that the observed increase in near-surface temperatures in recent decades is of natural origin is estimated to be less than 5%. However, this number is dependent on the estimated natural variability level, which is still subject to some uncertainty.

Multimodel assessment of the factors driving stratospheric ozone evolution over the 21st century

The evolution of stratospheric ozone from 1960 to 2100 is examined in simulations from 14 chemistry‐climate models, driven by prescribed levels of halogens and greenhouse gases. There is general agreement among the models that total column ozone reached a minimum around year 2000 at all latitudes, projected to be followed by an increase over the first half of the 21st century. In the second half of the 21st century, ozone is projected to continue increasing, level off, or even decrease depending on the latitude. Separation into partial columns above and below 20 hPa reveals that these latitudinal differences are almost completely caused by differences in the model projections of ozone in the lower stratosphere. At all latitudes, upper stratospheric ozone increases throughout the 21st century and is projected to return to 1960 levels well before the end of the century, although there is a spread among models in the dates that ozone returns to specific historical values. We find decreasing halogens and declining upper atmospheric temperatures, driven by increasing greenhouse gases, contribute almost equally to increases in upper stratospheric ozone. In the tropical lower stratosphere, an increase in upwelling causes a steady decrease in ozone through the 21st century, and total column ozone does not return to 1960 levels in most of the models. In contrast, lower stratospheric and total column ozone in middle and high latitudes increases during the 21st century, returning to 1960 levels well before the end of the century in most models.

Issues in stratosphere-troposphere coupling

An overview is given of current issues concerning the coupling between the stratosphere and troposphere. The tropopause region, more generally the upper troposphere/lower stratosphere, is the region of direct contact where exchange of material takes place. Dynamical coupling through angular momentum transfer by waves occurs nonlocally, and provides a generally negative torque on the stratosphere which drives an equator-to-pole circulation (i.e., towards the Earth’s axis of rotation). This wave-driven circulation is the principal mechanism for intraseasonal and interannual variability in the extratropical stratosphere. Although such variability is generally dynamical in origin, there are important chemical and radiative feedbacks. The location of the tropopause has implications for radiative forcing of climate, through its effect on the distribution of relatively short-lived greenhouse gases (ozone and water vapour). Some outstanding puzzles in our current understanding are identified. Attention is focused on possible climate sensitivities, and how these may be tested and constrained. Results from the Canadian Middle Atmosphere Model (CMAM), a fully interactive radiative-chemical-dynamical general circulation model, are used to illustrate some of the points.

Determination of the gas diffusion coefficient of a peat grassland soil

Peatland habitats are important carbon stocks that also have the potential to be significant sources of greenhouse gases, particularly when subject to changes such as artificial drainage and application of fertilizer. Models aiming to estimate greenhouse gas release from peatlands require an accurate estimate of the diffusion coefficient of gas transport through soil (Ds). The availability of specific measurements for peatland soils is currently limited. This study measured Ds for a peat soil with an overlying clay horizon and compared values with those from widely available models. The Ds value of a sandy loam reference soil was measured for comparison. Using the Currie (1960) method, Ds was measured between an air-filled porosity (ϵ) range of 0 and 0.5 cm3 cm−3. Values of Ds for the peat cores ranged between 3.2 × 10−4 and 4.4 × 10−3 m2 hour−1, for loamy clay cores between 0 and 4.7 × 10−3 m2 hour−1 and for the sandy reference soil they were between 5.4 × 10−4 and 3.4 × 10−3 m2 hour−1. The agreement of measured and modelled values of relative diffusivity (Ds/D0, with D0 the diffusion coefficient through free air) varied with soil type; however, the Campbell (1985) model provided the best replication of measured values for all soils. This research therefore suggests that the use of the Campbell model in the absence of accurately measured Ds and porosity values for a study soil would be appropriate. Future research into methods to reduce shrinkage of peat during measurement and therefore allow measurement of Ds for a greater range of ϵ would be beneficial.

Historical and idealized climate model experiments: an intercomparison of Earth system models of intermediate complexity

Both historical and idealized climate model experiments are performed with a variety of Earth system models of intermediate complexity (EMICs) as part of a community contribution to the Intergovernmental Panel on Climate Change Fifth Assessment Report. Historical simulations start at 850 CE and continue through to 2005. The standard simulations include changes in forcing from solar luminosity, Earth's orbital configuration, CO2, additional greenhouse gases, land use, and sulphate and volcanic aerosols. In spite of very different modelled pre-industrial global surface air temperatures, overall 20th century trends in surface air temperature and carbon uptake are reasonably well simulated when compared to observed trends. Land carbon fluxes show much more variation between models than ocean carbon fluxes, and recent land fluxes appear to be slightly underestimated. It is possible that recent modelled climate trends or climate–carbon feedbacks are overestimated resulting in too much land carbon loss or that carbon uptake due to CO2 and/or nitrogen fertilization is underestimated. Several one thousand year long, idealized, 2 × and 4 × CO2 experiments are used to quantify standard model characteristics, including transient and equilibrium climate sensitivities, and climate–carbon feedbacks. The values from EMICs generally fall within the range given by general circulation models. Seven additional historical simulations, each including a single specified forcing, are used to assess the contributions of different climate forcings to the overall climate and carbon cycle response. The response of surface air temperature is the linear sum of the individual forcings, while the carbon cycle response shows a non-linear interaction between land-use change and CO2 forcings for some models. Finally, the preindustrial portions of the last millennium simulations are used to assess historical model carbon-climate feedbacks. Given the specified forcing, there is a tendency for the EMICs to underestimate the drop in surface air temperature and CO2 between the Medieval Climate Anomaly and the Little Ice Age estimated from palaeoclimate reconstructions. This in turn could be a result of unforced variability within the climate system, uncertainty in the reconstructions of temperature and CO2, errors in the reconstructions of forcing used to drive the models, or the incomplete representation of certain processes within the models. Given the forcing datasets used in this study, the models calculate significant land-use emissions over the pre-industrial period. This implies that land-use emissions might need to be taken into account, when making estimates of climate–carbon feedbacks from palaeoclimate reconstructions.

Identification of micro-scale anthropogenic CO2, heat and moisture sources – processing eddy covariance fluxes for a dense urban environment

Anthropogenic emissions of heat and exhaust gases play an important role in the atmospheric boundary layer, altering air quality, greenhouse gas concentrations and the transport of heat and moisture at various scales. This is particularly evident in urban areas where emission sources are integrated in the highly heterogeneous urban canopy layer and directly linked to human activities which exhibit significant temporal variability. It is common practice to use eddy covariance observations to estimate turbulent surface fluxes of latent heat, sensible heat and carbon dioxide, which can be attributed to a local scale source area. This study provides a method to assess the influence of micro-scale anthropogenic emissions on heat, moisture and carbon dioxide exchange in a highly urbanized environment for two sites in central London, UK. A new algorithm for the Identification of Micro-scale Anthropogenic Sources (IMAS) is presented, with two aims. Firstly, IMAS filters out the influence of micro-scale emissions and allows for the analysis of the turbulent fluxes representative of the local scale source area. Secondly, it is used to give a first order estimate of anthropogenic heat flux and carbon dioxide flux representative of the building scale. The algorithm is evaluated using directional and temporal analysis. The algorithm is then used at a second site which was not incorporated in its development. The spatial and temporal local scale patterns, as well as micro-scale fluxes, appear physically reasonable and can be incorporated in the analysis of long-term eddy covariance measurements at the sites in central London. In addition to the new IMAS-technique, further steps in quality control and quality assurance used for the flux processing are presented. The methods and results have implications for urban flux measurements in dense urbanised settings with significant sources of heat and greenhouse gases.