The albedo of snow for partially cloudy skies

The albedo of snow for different cloudiness conditions is an important parameter in the Earth's radiation budget analysis and in the study of snowpack's thermal conditions. In this study an efficient approximate method is derived to calculate the incident spectral solar flux and snow-cover albedo in terms of different atmospheric, cloud, and snow parameters. The global flux under partially cloudy skies is expressed in terms of the clear sky flux and a coefficient which models the effect of scattering and absorption by cloud patches and multiple reflections between the cloud base and snowcover. The direct and the diffuse components of the clear sky flux are obtained using the spectral flux outside the atmosphere and the spectral transmission coefficients for absorption and scattering by molecules and aerosols.The spectral snow reflectance model considers both specular surface reflection and volumetric multiple scattering. The surface reflection is calculated by using a crystal-shape-dependent bidirectional reflectance distribution function; the volumetric multiple scattering is calculated by using a crystal-size-dependent approximate solution in the radiative transfer equation. The input parameters to the model are atmospheric precipitable water, ozone content, turbidity, cloud optical thickness, the size and shape of ice crystals of snow and surface pressure. The model yields spectral and integrated solar flux and snow reflectance as a function of solar elevation and fractional cloudcover.The model is illustrated using representative parameters for the Antarctic coastal regions. The albedo for a clear sky depends inversely on the solar elevation. At high elevations the albedo depends primarily upon the grain size; at low elevation the albedo depends on grain size and shape. The gradient of the albedo-elevation curve increases as the grains become larger and faceted. The albedo for a densely overcast sky is a few percent higher than the clear-sky albedo at high elevations. A simple relationship between grain size and the overcast albedo is obtained. For a set of grain size and shape, the albedo as a function of solar elevation and fractional cloud cover is tabulated.     

Global-scale climate impact functions: the relationship between climate forcing and impact

Although there is a strong policy interest in the impacts of climate change corresponding to different degrees of climate change, there is so far little consistent empirical evidence of the relationship between climate forcing and impact. This is because the vast majority of impact assessments use emissions-based scenarios with associated socio-economic assumptions, and it is not feasible to infer impacts at other temperature changes by interpolation. This paper presents an assessment of the global-scale impacts of climate change in 2050 corresponding to defined increases in global mean temperature, using spatially-explicit impacts models representing impacts in the water resources, river flooding, coastal, agriculture, ecosystem and built environment sectors. Pattern-scaling is used to construct climate scenarios associated with specific changes in global mean surface temperature, and a relationship between temperature and sea level used to construct sea level rise scenarios. Climate scenarios are constructed from 21 climate models to give an indication of the uncertainty between forcing and response. The analysis shows that there is considerable uncertainty in the impacts associated with a given increase in global mean temperature, due largely to uncertainty in the projected regional change in precipitation. This has important policy implications. There is evidence for some sectors of a non-linear relationship between global mean temperature change and impact, due to the changing relative importance of temperature and precipitation change. In the socio-economic sectors considered here, the relationships are reasonably consistent between socio-economic scenarios if impacts are expressed in proportional terms, but there can be large differences in absolute terms. There are a number of caveats with the approach, including the use of pattern-scaling to construct scenarios, the use of one impacts model per sector, and the sensitivity of the shape of the relationships between forcing and response to the definition of the impact indicator.     

On tropospheric adjustment to forcing and climate feedbacks

Motivated by findings that major components of so-called cloud ??feedbacks?? are best understood as rapid responses to CO₂ forcing (Gregory and Webb in J Clim 21:58?C71, 2008), the top of atmosphere (TOA) radiative effects from forcing, and the subsequent responses to global surface temperature changes from all ??atmospheric feedbacks?? (water vapour, lapse rate, surface albedo, ??surface temperature?? and cloud) are examined in detail in a General Circulation Model. Two approaches are used: applying regressions to experiments as they approach equilibrium, and equilibrium experiments forced separately by CO₂ and patterned sea surface temperature perturbations alone. Results are analysed using the partial radiative perturbation (??PRP??) technique. In common with Gregory and Webb (J Clim 21:58?C71, 2008) a strong positive addition to ??forcing?? is found in the short wave (SW) from clouds. There is little evidence, however, of significant global scale rapid responses from long wave (LW) cloud, nor from surface albedo, SW water vapour or ??surface temperature??. These responses may be well understood to first order as classical ??feedbacks????i.e. as a function of global mean temperature alone and linearly related to it. Linear regression provides some evidence of a small rapid negative response in the LW from water vapour, related largely to decreased relative humidity (RH), but the response here, too, is dwarfed by subsequent response to warming. The large rapid SW cloud response is related to cloud fraction changes??and not optical properties??resulting from small cloud decreases ranging from the tropical mid troposphere to the mid latitude lower troposphere, in turn associated with decreased lower tropospheric RH. These regions correspond with levels of enhanced heating rates and increased temperatures from the CO₂ increase. The pattern of SW cloud fraction response to SST changes differs quite markedly to this, with large positive radiation responses originating in the upper troposphere, positive contributions in the lowest levels and patterns of positive/negative contributions in mid latitude low levels. Overall SW cloud feedback was diagnosed as negative, due to the substantial negative SW feedback in cloud optical properties more than offsetting these. This study therefore suggests the rapid response to CO₂ forcing is (apart from a possible small negative response from LW water vapour) essentially confined to cloud fraction changes affecting SW radiation, and further that significant feedbacks with temperature occur in all cloud components (including this one), and indeed in all other classically understood ??feedbacks??.     

On summing the components of radiative forcing of climate change

 Radiative forcing is a useful concept in determining the potential influence of a particular mechanism of climate change. However, due to the increased number of forcing agents identified over the past decade, the total radiative forcing is difficult to assess. By assigning a range of probability distribution functions to the individual radiative forcings and using a Monte-Carlo approach, we estimate the total radiative forcing since pre-industrial times including all quantitative radiative forcing estimates to date. The resulting total radiative forcing has a 75–97% probability of being positive (or similarly a 3–25% probability of being negative), with mean radiative forcing ranging from +0.68 to +1.34 W m⁻², and median radiative forcing ranging from +0.94 to +1.39 W m⁻². Received: 14 March 2001 / Accepted: 1 June 2001     

Aerosol and ozone changes as forcing for climate evolution between 1850 and 2100

Global aerosol and ozone distributions and their associated radiative forcings were simulated between 1850 and 2100 following a recent historical emission dataset and under the representative concentration pathways (RCP) for the future. These simulations were used in an Earth System Model to account for the changes in both radiatively and chemically active compounds, when simulating the climate evolution. The past negative stratospheric ozone trends result in a negative climate forcing culminating at ?0.15 W m^?2 in the 1990s. In the meantime, the tropospheric ozone burden increase generates a positive climate forcing peaking at 0.41 W m^?2. The future evolution of ozone strongly depends on the RCP scenario considered. In RCP4.5 and RCP6.0, the evolution of both stratospheric and tropospheric ozone generate relatively weak radiative forcing changes until 2060–2070 followed by a relative 30 % decrease in radiative forcing by 2100. In contrast, RCP8.5 and RCP2.6 model projections exhibit strongly different ozone radiative forcing trajectories. In the RCP2.6 scenario, both effects (stratospheric ozone, a negative forcing, and tropospheric ozone, a positive forcing) decline towards 1950s values while they both get stronger in the RCP8.5 scenario. Over the twentieth century, the evolution of the total aerosol burden is characterized by a strong increase after World War II until the middle of the 1980s followed by a stabilization during the last decade due to the strong decrease in sulfates in OECD countries since the 1970s. The cooling effects reach their maximal values in 1980, with ?0.34 and ?0.28 W m^?2 respectively for direct and indirect total radiative forcings. According to the RCP scenarios, the aerosol content, after peaking around 2010, is projected to decline strongly and monotonically during the twenty-first century for the RCP8.5, 4.5 and 2.6 scenarios. While for RCP6.0 the decline occurs later, after peaking around 2050. As a consequence the relative importance of the total cooling effect of aerosols becomes weaker throughout the twenty-first century compared with the positive forcing of greenhouse gases. Nevertheless, both surface ozone and aerosol content show very different regional features depending on the future scenario considered. Hence, in 2050, surface ozone changes vary between ?12 and +12 ppbv over Asia depending on the RCP projection, whereas the regional direct aerosol radiative forcing can locally exceed ?3 W m^?2.     

Spectral surface albedo over Morocco and its impact on radiative forcing of Saharan dust

In May–June 2006, airborne and ground-based solar (0.3–2.2 μm) and thermal infrared (4–42 μm) radiation measurements have been performed in Morocco within the Saharan Mineral Dust Experiment (SAMUM). Upwelling and downwelling solar irradiances have been measured using the Spectral Modular Airborne Radiation Measurement System (SMART)-Albedometer. With these data, the areal spectral surface albedo for typical surface types in southeastern Morocco was derived from airborne measurements for the first time. The results are compared to the surface albedo retrieved from collocated satellite measurements, and partly considerable deviations are observed. Using measured surface and atmospheric properties, the spectral and broad-band dust radiative forcing at top-of-atmosphere (TOA) and at the surface has been estimated. The impact of the surface albedo on the solar radiative forcing of Saharan dust is quantified. In the SAMUM case of 19 May 2006, TOA solar radiative forcing varies by 12 W m⁻² per 0.1 surface-albedo change. For the thermal infrared component, values of up to +22 W m⁻² were derived. The net (solar plus thermal infrared) TOA radiative forcing varies between −19 and +24 W m⁻² for a broad-band solar surface albedo of 0.0 and 0.32, respectively. Over the bright surface of southeastern Morocco, the Saharan dust always has a net warming effect.     

Climate sensitivity and climate change under strong forcing

A version of the National Centre for Atmospheric Research (NCAR) coupled climate model is integrated under current climate conditions and in a series of experiments with climate forcings ranging from modest to very strong. The purpose of the experiments is to investigate the nature and behaviour of the climate feedback/sensitivity of the model, its evolution with time and climate state, the robustness of model parameterizations as forcing levels increase, and the possibility of a “runaway” warming under strong forcing. The model is integrated for 50 years, or to failure, after increasing the solar constant by 2.5, 10, 15, 25, 35, and 45% of its control value. The model successfully completes 50 years of integration for the 2.5, 10, 15, and 25% solar constant increases but fails for increases of 35% and 45%. The effective global climate sensitivity evolves with time and analysis indicates that a new equilibrium will be obtained for the 2.5, 10, and 15% cases but that runaway warming is underway for the 25% increase in solar constant. Feedback processes are analysed both locally and globally in terms of longwave and shortwave, clear-sky/surface, and cloud forcing components. Feedbacks in the system must be negative overall and of sufficient strength to balance the positive forcing if the system is to attain a new equilibrium. Longwave negative feedback processes strengthen in a reasonably linear fashion as temperature increases but shortwave feedback processes do not. In particular, solar cloud feedback becomes less negative and, for the 25% forcing case, eventually becomes positive, resulting in temperatures that “run away”. The conditions under which a runaway climate warming might occur have previously been investigated using simpler models. For sufficiently strong forcing, the greenhouse effect of increasing water vapour in a warmer atmosphere is expected to overwhelm the negative feedback of the longwave cooling to space as temperature increases. This is not, however, the reason for the climate instability experienced in the GCM. Instead, the model experiences a “cloud feedback” warming whereby the decrease in cloudiness that occurs when temperature increases beyond a critical value results in an increased absorption of solar radiation by the system, leading to the runaway warming.     

Greenhouse-gas versus aerosol forcing and African climate response

There are many indicators that human activity may change climate conditions all around the globe through emissions of greenhouse gases. In addition, aerosol particles are emitted from various natural and anthropogenic sources. One important source of aerosols arises from biomass burning, particularly in low latitudes where shifting cultivation and land degradation lead to enhanced aerosol burden. In this study the counteracting effects of greenhouse gases and aerosols on African climate are compared using climate model experiments with fully interactive aerosols from different sources. The consideration of aerosol emissions induces a remarkable decrease in short-wave solar irradiation near the surface, especially in winter and autumn in tropical West Africa and the Congo Basin where biomass burning is mainly prevailing. This directly leads to a modification of the surface energy budget with reduced sensible heat fluxes. As a consequence, temperature decreases, compensating the strong warming signal due to enhanced trace gas concentrations. While precipitation in tropical Africa is less sensitive to the greenhouse warming, it tends to decrease, if the effect of aerosols from biomass burning is taken into account. This is partly due to the local impact of enhanced aerosol burden and partly to modifications of the large-scale monsoon circulation in the lower troposphere, usually lagging behind the season with maximum aerosol emissions. In the model equilibrium experiments, the greenhouse gas impact on temperature stands out from internal variability at various time scales from daily to decadaland the same holds for precipitation under the additional aerosol forcing. Greenhouse gases and aerosols exhibit an opposite effect on daily temperature extremes, resulting in an compensation of the individual responses under the combined forcing. In terms of precipitation, daily extreme events tend to be reduced under aerosol forcing, particularly over the tropical Atlantic and the Congo basin. These results suggest that the simulation of the multiple aerosol effects from anthropogenic sources represents an important factor in tropical climate change, hence, requiring more attention in climate modelling attempts.     

Quantifying global climate feedbacks, responses and forcing under abrupt and gradual CO2 forcing

Experiments with abrupt CO₂ forcing allow the diagnosis of the response of global mean temperature and precipitation in terms of fast temperature independent adjustments and slow, linear temperature-dependent feedbacks. Here we compare responses, feedbacks and forcings in experiments performed as part of version 5 of the coupled model inter-comparison project (CMIP5). The experiments facilitate, for the first time, a comparison of fully coupled atmosphere-ocean general circulation models (GCM’s) under both linearly increasing and abrupt radiative forcing. In the case of a 1 % per year compounded increase in CO₂ concentration, we find that the non-linear evolution of surface air temperature in time, when combined with the linear evolution of the radiative balance at the top of the atmosphere, results in a feedback parameter and effective climate sensitivity having an offset compared to values computed from abrupt 4× CO₂ forcing experiments. The linear evolution of the radiative balance at the top of the atmosphere also contributes to an offset between the global mean precipitation response predicted in the 1 % experiment using linear theory and that diagnosed from the experiments themselves, and a potential error between the adjusted radiative forcing and that produced using a standard linear formula. The non-linear evolution of temperature and precipitation responses are also evident in the RCP8.5 scenario and have implications for understanding, quantifying and emulating the global response of the CMIP5 climate GCMs.     

Origins of differences in climate sensitivity, forcing and feedback in climate models

We diagnose climate feedback parameters and CO₂ forcing including rapid adjustment in twelve atmosphere/mixed-layer-ocean (“slab”) climate models from the CMIP3/CFMIP-1 project (the AR4 ensemble) and fifteen parameter-perturbed versions of the HadSM3 slab model (the PPE). In both ensembles, differences in climate feedbacks can account for approximately twice as much of the range in climate sensitivity as differences in CO₂ forcing. In the AR4 ensemble, cloud effects can explain the full range of climate sensitivities, and cloud feedback components contribute four times as much as cloud components of CO₂ forcing to the range. Non-cloud feedbacks are required to fully account for the high sensitivities of some models however. The largest contribution to the high sensitivity of HadGEM1 is from a high latitude clear-sky shortwave feedback, and clear-sky longwave feedbacks contribute substantially to the highest sensitivity members of the PPE. Differences in low latitude ocean regions (30°N/S) contribute more to the range than those in mid-latitude oceans (30–55°N/S), low/mid latitude land (55°N/S) or high latitude ocean/land (55–90°N/S), but contributions from these other regions are required to account fully for the higher model sensitivities, for example from land areas in IPSL CM4. Net cloud feedback components over the low latitude oceans sorted into percentile ranges of lower tropospheric stability (LTS) show largest differences among models in stable regions, mainly due to their shortwave components, most of which are positive in spite of increasing LTS. Differences in the mid-stability range are smaller, but cover a larger area, contributing a comparable amount to the range in climate sensitivity. These are strongly anti-correlated with changes in subsidence. Cloud components of CO₂ forcing also show the largest differences in stable regions, and are strongly anticorrelated with changes in estimated inversion strength (EIS). This is qualitatively consistent with what would be expected from observed relationships between EIS and low-level cloud fraction. We identify a number of cases where individual models show unusually strong forcings and feedbacks compared to other members of the ensemble. We encourage modelling groups to investigate unusual model behaviours further with sensitivity experiments. Most of the models fail to correctly reproduce the observed relationships between stability and cloud radiative effect in the subtropics, indicating that there remains considerable room for model improvements in the future.     

The influence of inter-annually varying albedo on regional climate and drought

Albedo plays an important role in land–atmosphere interactions and local climate. This study presents the impact on simulating regional climate, and the evolution of a drought, when using the default climatological albedo as is usually done in regional climate modelling, or using the actual observed albedo which is rarely done. Here, time-varying satellite derived albedo data is used to update the lower boundary condition of the Weather Research and Forecasting regional climate model in order to investigate the influence of observed albedo on regional climate simulations and also potential changes to land–atmosphere feedback over south-east Australia. During the study period from 2000 to 2008, observations show that albedo increased with an increasingly negative precipitation anomaly, though it lagged precipitation by several months. Compared to in-situ observations, using satellite observed albedo instead of the default climatological albedo provided an improvement in the simulated seasonal mean air temperature. In terms of precipitation, both simulations reproduced the drought that occurred from 2002 through 2006. Using the observed albedo produced a drier simulation overall. During the onset of the 2002 drought, albedo changes enhanced the precipitation reduction by 20 % on average, over locations where it was active. The area experiencing drought increased 6.3 % due to the albedo changes. Two mechanisms for albedo changes to impact land–atmosphere drought feedback are investigated. One accounts for the increased albedo, leading to reduced turbulent heat flux and an associated decrease of moist static energy density in the planetary boundary layer; the other considers that enhanced local radiative heating, due to the drought, favours a deeper planetary boundary layer, subsequently decreasing the moist static energy density through entrainment of the free atmosphere. Analysis shows that drought related large-scale changes in the regional climate favour a strengthening of the second mechanism. That is, the second mechanism is stronger in a drought year compared to a normal year and this difference is larger than for the first mechanism. When both mechanisms are active, the second mechanism tends to dominate across the model domain, particularly during the 2002 drought period. The introduction of observed inter-annual variations in albedo produces an enhancement of the first mechanism and a weakening of the second mechanism during the onset of the drought.     

Solar irradiance forcing of centennial climate variability during the Holocene

Centennial climate variability during the Holocene has been simulated in two 10,000 year experiments using the intermediate-complexity ECBilt model. ECBilt contains a dynamic atmosphere, a global 3-D ocean model and a thermodynamic sea-ice model. One experiment uses orbital forcing and solar irradiance forcing, which is based on the Stuiver et al. residual ¹⁴C record spliced into the Lean et al. reconstruction. The other experiment uses orbital forcing alone. A glacier model is coupled off-line to the climate model. A time scale analysis shows that the response in atmospheric parameters to the irradiance forcing can be characterised as the direct response of a system with a large thermal inertia. This is evident in parameters like surface air temperature, monsoon precipitation and glacier length, which show a stronger response for longer time scales. The oceanic response, on the other hand, is strongly modified by internal feedback processes. The solar irradiance forcing excites a (damped) mode of the thermohaline circulation (THC) in the North Atlantic Ocean, similar to the loop-oscillator modes associated with random-noise freshwater forcing. This results in a significant peak (at time scales 200–250 year) in the THC spectrum which is absent in the reference run. The THC response diminishes the sea surface temperature response at high latitudes, while it gives rise to a signal in the sea surface salinity. A comparison of the model results with observations shows a number of encouraging similarities.     

Late quaternary volcanic modulation of Milankovitch climate forcing

Summary A preliminary global volcanicity index time series was developed from the total available record of radiocarbon dated eruptions. This was combined with a previously developed simple model of Pleistocene glaciation based on the integral of the seasonal contrast of hemispheric irradiance. The combined integral gave a modeled glacial volume time-series which matched both the general variation of glacial volume in the last forty thousand years and the short-term events as well. Particularly well-modeled were the Alleröd-Younger Dryas-Two Creeks complex of Late Glacial events and the Cochrane stillstand.From this time-series was derived an ice albedo series which was combined with the calculated irradiance to drive a hemispheric heat budget model to obtain time series of monthly hemispheric surface temperatures. The climatic episodes depicted by these time series match the known field data quite well. Indeed the modeling of the Little Ice Age and other features suggest that most of the century-to-millennia events are driven by variations in volcanicity.

---

Zusammenfassung Von der Gesamtheit der erhältlichen Aufzeichnungen radiokarbon-datierter Ausbrüche wurde eine vorläufige globale Vulkanitätsindex-Zeitreihe erstellt. Diese wurde mit einem früher entwickelten, einfachen Modell der Pleistozänvereisung, das auf dem Integral des jahreszeitlichen Kontrasts der hemisphärischen Einstrahlung basiert, kombiniert. Das so kombinierte Integral ergab ein Modell der Zeitreihe des Eisvolumens, das sowohl den allgemeinen Schwankungen (des Eisvolumens) in den letzten 40 000 Jahren als auch kurzfristigen Ereignissen Rechnung trug. Besonders gut wiedergegeben waren der Allerod-Younger-Dryas-Two-Creeks-Komplex der jüngeren Eiszeit und der Cochrane-Stillstand.Von diesen Zeitreihen wurde eine Eisalbedo-Zeitreihe abgeleitet, die in Kombination mit der errechneten Einstrahlung die Aufstellung eines hemisphärischen Wärmebilanzmodells erlaubte und somit Zeitreihen des Monatsmittels hemisphärischer Oberflächentemperatur. Die klimatischen Episoden, die in diesen Zeitserien dargestellt sind, passen recht gut zu den bereits bekannten Felddaten. Tatsächlich läßt die Modellierung der Kleinen Eiszeit und anderer Erscheinungen vermuten, daß die meisten Ereignisse im Zeitraum von Jahrhunderten bis Jahrtausenden durch Schwankungen der Vulkanität gesteuert werden.

---

---

With 8 Figures     

Intercomparison and validation of snow albedo parameterization schemes in climate models

Snow albedo is known to be crucial for heat exchange at high latitudes and high altitudes, and is also an important parameter in General Circulation Models (GCMs) because of its strong positive feedback properties. In this study, seven GCM snow albedo schemes and a multiple linear regression model were intercompared and validated against 59 years of in situ data from Svalbard, the French Alps and six stations in the former Soviet Union. For each site, the significant meteorological parameters for modeling the snow albedo were identified by constructing the 95% confidence intervals. The significant parameters were found to be: temperature, snow depth, positive degree day and a dummy of snow depth, and the multiple linear regression model was constructed to include these. Overall, the intercomparison showed that the modeled snow albedo varied more than the observed albedo for all models, and that the albedo was often underestimated. In addition, for several of the models, the snow albedo decreased at a faster rate or by a greater magnitude during the winter snow metamorphosis than the observed albedo. Both the temperature dependent schemes and the prognostic schemes showed shortcomings.     

Relative roles of climate sensitivity and forcing in defining the ocean circulation response to climate change

The response of the ocean’s meridional overturning circulation (MOC) to increased greenhouse gas forcing is examined using a coupled model of intermediate complexity, including a dynamic 3-D ocean subcomponent. Parameters are the increase in CO₂ forcing (with stabilization after a specified time interval) and the model’s climate sensitivity. In this model, the cessation of deep sinking in the north “Atlantic” (hereinafter, a “collapse”), as indicated by changes in the MOC, behaves like a simple bifurcation. The final surface air temperature (SAT) change, which is closely predicted by the product of the radiative forcing and the climate sensitivity, determines whether a collapse occurs. The initial transient response in SAT is largely a function of the forcing increase, with higher sensitivity runs exhibiting delayed behavior; accordingly, high CO₂-low sensitivity scenarios can be assessed as a recovering or collapsing circulation shortly after stabilization, whereas low CO₂-high sensitivity scenarios require several hundred additional years to make such a determination. We also systemically examine how the rate of forcing, for a given CO₂ stabilization, affects the ocean response. In contrast with previous studies based on results using simpler ocean models, we find that except for a narrow range of marginally stable to marginally unstable scenarios, the forcing rate has little impact on whether the run collapses or recovers. In this narrow range, however, forcing increases on a time scale of slow ocean advective processes results in weaker declines in overturning strength and can permit a run to recover that would otherwise collapse.     

Satellite radiances and clear air turbulence probabilities

Radiance gradients recorded by Nimbus III and IV, over the United States and the Atlantic Ocean, respectively, are well-related to probability of clear air turbulence. Making use of satellite data only, it is possible to avoid CAT almost entirely. In general, radiosonde-measured temperature gradients are not as well correlated with CAT as satellite radiance gradients; but if the satellite track is parallel to the isotherms, radiance gradients along the satellite track are negligible; and temperature gradients at constant pressure are better indicators of CAT unless the satellite instrument scans at right angles to the track.     

Sub-Saharan Northern African climate at the end of the twenty-first century: forcing factors and climate change processes

A regional climate model, the Weather Research and Forecasting (WRF) Model, is forced with increased atmospheric CO₂ and anomalous SSTs and lateral boundary conditions derived from nine coupled atmosphere–ocean general circulation models to produce an ensemble set of nine future climate simulations for northern Africa at the end of the twenty-first century. A well validated control simulation, agreement among ensemble members, and a physical understanding of the future climate change enhance confidence in the predictions. The regional model ensembles produce consistent precipitation projections over much of northern tropical Africa. A moisture budget analysis is used to identify the circulation changes that support future precipitation anomalies. The projected midsummer drought over the Guinean Coast region is related partly to weakened monsoon flow. Since the rainfall maximum demonstrates a southward bias in the control simulation in July–August, this may be indicative of future summer drying over the Sahel. Wetter conditions in late summer over the Sahel are associated with enhanced moisture transport by the West African westerly jet, a strengthening of the jet itself, and moisture transport from the Mediterranean. Severe drought in East Africa during August and September is accompanied by a weakened Indian monsoon and Somali jet. Simulations with projected and idealized SST forcing suggest that overall SST warming in part supports this regional model ensemble agreement, although changes in SST gradients are important over West Africa in spring and fall. Simulations which isolate the role of individual climate forcings suggest that the spatial distribution of the rainfall predictions is controlled by the anomalous SST and lateral boundary conditions, while CO₂ forcing within the regional model domain plays an important secondary role and generally produces wetter conditions.