The role of dust in climate changes today, at the last glacial maximum and in the future.

Natural mineral aerosol (dust) is an active component of the climate system and plays multiple roles in mediating physical and biogeochemical exchanges between the atmosphere, land surface and ocean. Changes in the amount of dust in the atmosphere are caused both by changes in climate (precipitation, wind strength, regional moisture balance) and changes in the extent of dust sources caused by either anthropogenic or climatically induced changes in vegetation cover. Models of the global dust cycle take into account the physical controls on dust deflation from prescribed source areas (based largely on soil wetness and vegetation cover thresholds), dust transport within the atmospheric column, and dust deposition through sedimentation and scavenging by precipitation. These models successfully reproduce the first-order spatial and temporal patterns in atmospheric dust loading under modern conditions. Atmospheric dust loading was as much as an order-of-magnitude larger than today during the last glacial maximum (LGM). While the observed increase in emissions from northern Africa can be explained solely in terms of climate changes (colder, drier and windier glacial climates), increased emissions from other regions appear to have been largely a response to climatically induced changes in vegetation cover and hence in the extent of dust source areas. Model experiments suggest that the increased dust loading in tropical regions had an effect on radiative forcing comparable to that of low glacial CO2 levels. Changes in land-use are already increasing the dust loading of the atmosphere. However, simulations show that anthropogenically forced climate changes substantially reduce the extent and productivity of natural dust sources. Positive feedbacks initiated by a reduction of dust emissions from natural source areas on both radiative forcing and atmospheric CO2 could substantially mitigate the impacts of land-use changes, and need to be considered in climate change assessments.

Intercomparison of simulated global vegetation distribution in response to 6 kyr BP orbital forcing

The response of ten atmospheric general circulation models to orbital forcing at 6 kyr BP has been investigated using the BIOME model, which predicts equilibrium vegetation distribution, as a diagnostic. Several common features emerge: (a) reduced tropical rain forest as a consequence of increased aridity in the equatorial zone, (b) expansion of moisture-demanding vegetation in the Old World subtropics as a consequence of the expansion of the Afro–Asian monsoon, (c) an increase in warm grass/shrub in the Northern Hemisphere continental interiors in response to warming and enhanced aridity, and (d) a northward shift in the tundra–forest boundary in response to a warmer growing season at high northern latitudes. These broadscale features are consistent from model to model, but there are differences in their expression at a regional scale. Vegetation changes associated with monsoon enhancement and high-latitude summer warming are consistent with palaeoenvironmental observations, but the simulated shifts in vegetation belts are too small in both cases. Vegetation changes due to warmer and more arid conditions in the midcontinents of the Northern Hemisphere are consistent with palaeoenvironmental data from North America, but data from Eurasia suggests conditions were wetter at 6 kyr BP than today. The models show quantitatively similar vegetation changes in the intertropical zone, and in the northern and southern extratropics. The small differences among models in the magnitude of the global vegetation response are not related to differences in global or zonal climate averages, but reflect differences in simulated regional features. Regional-scale analyses will therefore be necessary to identify the underlying causes of such differences among models.

Climate and biome simulations for the past 21,000 years

This paper reports on a set of paleoclimate simulations for 21, 16, 14, 11 and 6 ka (thousands of years ago) carried out with the Community Climate Model, Version 1 (CCM1) of the National Center for Atmospheric Research (NCAR). This climate model uses four interactive components that were not available in our previous simulations with the NCAR CCM0 (COHMAP, 1988Science, 241, 1043–1052; Wright et al., 1993Global Climate Since the Last Glocial Maximum, University of Minnesota Press, MN): soil moisture, snow hydrology, sea-ice, and mixed-layer ocean temperature. The new simulations also use new estimates of ice sheet height and size from ( Peltier 1994, Science, 265, 195–201), and synchronize the astronomically dated orbital forcing with the ice sheet and atmospheric CO2 levels corrected from radiocarbon years to calendar years. The CCM1 simulations agree with the previous simulations in their most general characteristics. The 21 ka climate is cold and dry, in response to the presence of the ice sheets and lowered CO2 levels. The period 14–6 ka has strengthened northern summer monsoons and warm mid-latitude continental interiors in response to orbital changes. Regional differences between the CCM1 and CCM0 simulations can be traced to the effects of either the new interactive model components or the new boundary conditions. CCM1 simulates climate processes more realistically, but has additional degrees of freedom that can allow the model to ‘drift’ toward less realistic solutions in some instances. The CCM1 simulations are expressed in terms of equilibrium vegetation using BIOME 1, and indicate large shifts in biomes. Northern tundra and forest biomes are displaced southward at glacial maximum and subtropical deserts contract in the mid-Holocene when monsoons strengthen. These vegetation changes could, if simulated interactively, introduce additional climate feedbacks. The total area of vegetated land remains nearly constant through time because the exposure of continental shelves with lowered sea level largely compensates for the land covered by the expanded ice sheets.

The interpretation and use of biases in decadal climate predictions

Decadal climate predictions exhibit large biases, which are often subtracted and forgotten. However, understanding the causes of bias is essential to guide efforts to improve prediction systems, and may offer additional benefits. Here the origins of biases in decadal predictions are investigated, including whether analysis of these biases might provide useful information. The focus is especially on the lead-time-dependent bias tendency. A “toy” model of a prediction system is initially developed and used to show that there are several distinct contributions to bias tendency. Contributions from sampling of internal variability and a start-time-dependent forcing bias can be estimated and removed to obtain a much improved estimate of the true bias tendency, which can provide information about errors in the underlying model and/or errors in the specification of forcings. It is argued that the true bias tendency, not the total bias tendency, should be used to adjust decadal forecasts. The methods developed are applied to decadal hindcasts of global mean temperature made using the Hadley Centre Coupled Model, version 3 (HadCM3), climate model, and it is found that this model exhibits a small positive bias tendency in the ensemble mean. When considering different model versions, it is shown that the true bias tendency is very highly correlated with both the transient climate response (TCR) and non–greenhouse gas forcing trends, and can therefore be used to obtain observationally constrained estimates of these relevant physical quantities.

A simple framework for assessing the trade-off between the climate impact of aviation carbon dioxide emissions and contrails for a single flight

Persistent contrails are an important climate impact of aviation which could potentially be reduced by re-routing aircraft to avoid contrailing; however this generally increases both the flight length and its corresponding CO emissions. Here, we provide a simple framework to assess the trade-off between the climate impact of CO emissions and contrails for a single flight, in terms of the absolute global warming potential and absolute global temperature potential metrics for time horizons of 20, 50 and 100 years. We use the framework to illustrate the maximum extra distance (with no altitude changes) that can be added to a flight and still reduce its overall climate impact. Small aircraft can fly up to four times further to avoid contrailing than large aircraft. The results have a strong dependence on the applied metric and time horizon. Applying a conservative estimate of the uncertainty in the contrail radiative forcing and climate efficacy leads to a factor of 20 difference in the maximum extra distance that could be flown to avoid a contrail. The impact of re-routing on other climatically-important aviation emissions could also be considered in this framework.

Mean climate and variability of the atmosphere and ocean on an aquaplanet

Numerical experiments are described that pertain to the climate of a coupled atmosphere–ocean–ice system in the absence of land, driven by modern-day orbital and CO2 forcing. Millennial time-scale simulations yield a mean state in which ice caps reach down to 55° of latitude and both the atmosphere and ocean comprise eastward- and westward-flowing zonal jets, whose structure is set by their respective baroclinic instabilities. Despite the zonality of the ocean, it is remarkably efficient at transporting heat meridionally through the agency of Ekman transport and eddy-driven subduction. Indeed the partition of heat transport between the atmosphere and ocean is much the same as the present climate, with the ocean dominating in the Tropics and the atmosphere in the mid–high latitudes. Variability of the system is dominated by the coupling of annular modes in the atmosphere and ocean. Stochastic variability inherent to the atmospheric jets drives variability in the ocean. Zonal flows in the ocean exhibit decadal variability, which, remarkably, feeds back to the atmosphere, coloring the spectrum of annular variability. A simple stochastic model can capture the essence of the process. Finally, it is briefly reviewed how the aquaplanet can provide information about the processes that set the partition of heat transport and the climate of Earth.

Sensitivity of simulated regional Arctic climate to the choice of coupled model domain

The climate over the Arctic has undergone changes in recent decades. In order to evaluate the coupled response of the Arctic system to external and internal forcing, our study focuses on the estimation of regional climate variability and its dependence on large-scale atmospheric and regional ocean circulations. A global ocean–sea ice model with regionally high horizontal resolution is coupled to an atmospheric regional model and global terrestrial hydrology model. This way of coupling divides the global ocean model setup into two different domains: one coupled, where the ocean and the atmosphere are interacting, and one uncoupled, where the ocean model is driven by prescribed atmospheric forcing and runs in a so-called stand-alone mode. Therefore, selecting a specific area for the regional atmosphere implies that the ocean–atmosphere system can develop ‘freely’ in that area, whereas for the rest of the global ocean, the circulation is driven by prescribed atmospheric forcing without any feedbacks. Five different coupled setups are chosen for ensemble simulations. The choice of the coupled domains was done to estimate the influences of the Subtropical Atlantic, Eurasian and North Pacific regions on northern North Atlantic and Arctic climate. Our simulations show that the regional coupled ocean–atmosphere model is sensitive to the choice of the modelled area. The different model configurations reproduce differently both the mean climate and its variability. Only two out of five model setups were able to reproduce the Arctic climate as observed under recent climate conditions (ERA-40 Reanalysis). Evidence is found that the main source of uncertainty for Arctic climate variability and its predictability is the North Pacific. The prescription of North Pacific conditions in the regional model leads to significant correlation with observations, even if the whole North Atlantic is within the coupled model domain. However, the inclusion of the North Pacific area into the coupled system drastically changes the Arctic climate variability to a point where the Arctic Oscillation becomes an ‘internal mode’ of variability and correlations of year-to-year variability with observational data vanish. In line with previous studies, our simulations provide evidence that Arctic sea ice export is mainly due to ‘internal variability’ within the Arctic region. We conclude that the choice of model domains should be based on physical knowledge of the atmospheric and oceanic processes and not on ‘geographic’ reasons. This is particularly the case for areas like the Arctic, which has very complex feedbacks between components of the regional climate system.

Patterns of maximum body size evolution in Cenozoic land mammals: eco-evolutionary processes and abiotic forcing

There is accumulating evidence that macroevolutionary patterns of mammal evolution during the Cenozoic follow similar trajectories on different continents. This would suggest that such patterns are strongly determined by global abiotic factors, such as climate, or by basic eco-evolutionary processes such as filling of niches by specialization. The similarity of pattern would be expected to extend to the history of individual clades. Here, we investigate the temporal distribution of maximum size observed within individual orders globally and on separate continents. While the maximum size of individual orders of large land mammals show differences and comprise several families, the times at which orders reach their maximum size over time show strong congruence, peaking in the Middle Eocene, the Oligocene and the Plio-Pleistocene. The Eocene peak occurs when global temperature and land mammal diversity are high and is best explained as a result of niche expansion rather than abiotic forcing. Since the Eocene, there is a significant correlation between maximum size frequency and global temperature proxy. The Oligocene peak is not statistically significant and may in part be due to sampling issues. The peak in the Plio-Pleistocene occurs when global temperature and land mammal diversity are low, it is statistically the most robust one and it is best explained by global cooling. We conclude that the macroevolutionary patterns observed are a result of the interplay between eco-evolutionary processes and abiotic forcing

Feedbacks and mechanisms affecting the global sensitivity of glaciers to climate change

Mass loss by glaciers has been an important contributor to sea level rise in the past, and is projected to contribute a substantial fraction of total sea level rise during the 21st century. Here, we use a model of the world's glaciers to quantify equilibrium sensitivities of global glacier mass to climate change, and to investigate the role of changes in glacier hypsometry for long-term mass changes. We find that 21st century glacier-mass loss is largely governed by the glacier's response to 20th century climate change. This limits the influence of 21st century climate change on glacier-mass loss, and explains why there are relatively small differences in glacier-mass loss under greatly different scenarios of climate change. The projected future changes in both temperature and precipitation experienced by glaciers are amplified relative to the global average. The projected increase in precipitation partly compensates for the mass loss caused by warming, but this compensation is negligible at higher temperature anomalies since an increasing fraction of precipitation at the glacier sites is liquid. Loss of low-lying glacier area, and more importantly, eventual complete disappearance of glaciers, strongly limit the projected sea level contribution from glaciers in coming centuries. The adjustment of glacier hypsometry to changes in the forcing strongly reduces the rates of global glacier-mass loss caused by changes in global mean temperature compared to rates of mass loss when hypsometric changes are neglected. This result is a second reason for the relatively weak dependence of glacier-mass loss on future climate scenario, and helps explain why glacier-mass loss in the first half of the 20th century was of the same order of magnitude as in the second half of the 20th century, even though the rate of warming was considerably smaller.

Broad range of 2050 warming from an observationally constrained large climate model ensemble

Incomplete understanding of three aspects of the climate system—equilibrium climate sensitivity, rate of ocean heat uptake and historical aerosol forcing—and the physical processes underlying them lead to uncertainties in our assessment of the global-mean temperature evolution in the twenty-first century1,2. Explorations of these uncertainties have so far relied on scaling approaches3,4, large ensembles of simplified climate models1,2, or small ensembles of complex coupled atmosphere–ocean general circulation models5,6 which under-represent uncertainties in key climate system properties derived from independent sources7–9. Here we present results from a multi-thousand-member perturbed-physics ensemble of transient coupled atmosphere–ocean general circulation model simulations. We find that model versions that reproduce observed surface temperature changes over the past 50 years show global-mean temperature increases of 1.4–3 K by 2050, relative to 1961–1990, under a mid-range forcing scenario. This range of warming is broadly consistent with the expert assessment provided by the Intergovernmental Panel on Climate Change Fourth Assessment Report10, but extends towards larger warming than observed in ensemblesof-opportunity5 typically used for climate impact assessments. From our simulations, we conclude that warming by the middle of the twenty-first century that is stronger than earlier estimates is consistent with recent observed temperature changes and a mid-range ‘no mitigation’ scenario for greenhouse-gas emissions.

Meteorological forcing data for urban outdoor thermal comfort models from a coupled convective boundary layer and surface energy balance scheme

Site-specific meteorological forcing appropriate for applications such as urban outdoor thermal comfort simulations can be obtained using a newly coupled scheme that combines a simple slab convective boundary layer (CBL) model and urban land surface model (ULSM) (here two ULSMs are considered). The former simulates daytime CBL height, air temperature and humidity, and the latter estimates urban surface energy and water balance fluxes accounting for changes in land surface cover. The coupled models are tested at a suburban site and two rural sites, one irrigated and one unirrigated grass, in Sacramento, U.S.A. All the variables modelled compare well to measurements (e.g. coefficient of determination = 0.97 and root mean square error = 1.5 °C for air temperature). The current version is applicable to daytime conditions and needs initial state conditions for the CBL model in the appropriate range to obtain the required performance. The coupled model allows routine observations from distant sites (e.g. rural, airport) to be used to predict air temperature and relative humidity in an urban area of interest. This simple model, which can be rapidly applied, could provide urban data for applications such as air quality forecasting and building energy modelling, in addition to outdoor thermal comfort.

Global analysis of radiative forcing from fire-induced shortwave albedo change

Land surface albedo, a key parameter to derive Earth's surface energy balance, is used in the parameterization of numerical weather prediction, climate monitoring and climate change impact assessments. Changes in albedo due to fire have not been fully investigated on a continental and global scale. The main goal of this study, therefore, is to quantify the changes in instantaneous shortwave albedo produced by biomass burning activities and their associated radiative forcing. The study relies on the MODerate-resolution Imaging Spectroradiometer (MODIS) MCD64A1 burned-area product to create an annual composite of areas affected by fire and the MCD43C2 bidirectional reflectance distribution function (BRDF) albedo snow-free product to compute a bihemispherical reflectance time series. The approximate day of burning is used to calculate the instantaneous change in shortwave albedo. Using the corresponding National Centers for Environmental Prediction (NCEP) monthly mean downward solar radiation flux at the surface, the global radiative forcing associated with fire was computed. The analysis reveals a mean decrease in shortwave albedo of −0.014 (1σ = 0.017), causing a mean positive radiative forcing of 3.99 Wm−2 (1σ = 4.89) over the 2002–20012 time period in areas affected by fire. The greatest drop in mean shortwave albedo change occurs in 2002, which corresponds to the highest total area burned (378 Mha) observed in the same year and produces the highest mean radiative forcing (4.5 Wm−2). Africa is the main contributor in terms of burned area, but forests globally give the highest radiative forcing per unit area and thus give detectable changes in shortwave albedo. The global mean radiative forcing for the whole period studied (~0.0275 Wm−2) shows that the contribution of fires to the Earth system is not insignificant.

Climate change projections and selected impacts for Vietnam

In September 2013, the 5th Assessment Report (5AR) of the International Panel on Climate Change (IPCC) has been released. Taking the 5AR cli-mate change scenarios into account, the World Bank published an earli-er report on climate change and its impacts on selected hot spot re-gions, including Southeast Asia. Currently, dynamical and statistical-dynamical downscaling efforts are underway to obtain higher resolution and more robust regional climate change projections for tropical South-east Asia, including Vietnam. Such initiatives are formalized under the World Meteorological Organization (WMO) Coordinated Regional Dynamic Downscaling Experiment (CORDEX) East Asia and Southeast Asia and also take place in climate change impact projects such as the joint Vietnam-ese-German project “Environmental and Water Protection Technologies of Coastal Zones in Vietnam (EWATEC-COAST)”. In this contribution, the lat-est assessments for changes in temperature, precipitation, sea level, and tropical cyclones (TCs) under the 5AR Representative Concentration Pathway (RCP) scenarios 4.5 and 8.5 are reviewed. Special emphasis is put on changes in extreme events like heat waves and/or heavy precipita-tion. A regional focus is Vietnam south of 16°N. A continued increase in mean near surface temperature is projected, reaching up to 5°C at the end of this century in northern Vietnam un-der the high greenhouse-gas forcing scenario RCP8.5. Overall, project-ed changes in annual precipitation are small, but there is a tendency of more rainfall in the boreal winter dry season. Unprecedented heat waves and an increase in extreme precipitation events are projected by both global and regional climate models. Globally, TCs are projected to decrease in number, but an increase in intensity of peak winds and rain-fall in the inner core region is estimated. Though an assessment of changes in land-falling frequency in Vietnam is uncertain due to difficul-ties in assessing changes in TC tracks, some work indicates a reduction in the number of land-falling TCs in Vietnam. Sea level may rise by 75-100 cm until the end of the century with the Vietnamese coastline experienc-ing 10-15% higher rise than on global average. Given the large rice and aquaculture production in the Mekong and Red River Deltas, that are both prone to TC-related storm surges and flooding, this poses a challenge to foodsecurity and protection of coastal population and assets.

Climate and carbon cycle changes from 1850 to 2100 in MPI-ESM simulations for the Coupled Model Intercomparison Project phase 5

The new Max-Planck-Institute Earth System Model (MPI-ESM) is used in the Coupled Model Intercomparison Project phase 5 (CMIP5) in a series of climate change experiments for either idealized CO2-only forcing or forcings based on observations and the Representative Concentration Pathway (RCP) scenarios. The paper gives an overview of the model configurations, experiments related forcings, and initialization procedures and presents results for the simulated changes in climate and carbon cycle. It is found that the climate feedback depends on the global warming and possibly the forcing history. The global warming from climatological 1850 conditions to 2080–2100 ranges from 1.5°C under the RCP2.6 scenario to 4.4°C under the RCP8.5 scenario. Over this range, the patterns of temperature and precipitation change are nearly independent of the global warming. The model shows a tendency to reduce the ocean heat uptake efficiency toward a warmer climate, and hence acceleration in warming in the later years. The precipitation sensitivity can be as high as 2.5% K−1 if the CO2 concentration is constant, or as small as 1.6% K−1, if the CO2 concentration is increasing. The oceanic uptake of anthropogenic carbon increases over time in all scenarios, being smallest in the experiment forced by RCP2.6 and largest in that for RCP8.5. The land also serves as a net carbon sink in all scenarios, predominantly in boreal regions. The strong tropical carbon sources found in the RCP2.6 and RCP8.5 experiments are almost absent in the RCP4.5 experiment, which can be explained by reforestation in the RCP4.5 scenario.

Deconstructing the climate change response of the Northern Hemisphere wintertime storm tracks

There are large uncertainties in the circulation response of the atmosphere to climate change. One manifestation of this is the substantial spread in projections for the extratropical storm tracks made by different state-of-the-art climate models. In this study we perform a series of sensitivity experiments, with the atmosphere component of a single climate model, in order to identify the causes of the differences between storm track responses in different models. In particular, the Northern Hemisphere wintertime storm tracks in the CMIP3 multi-model ensemble are considered. A number of potential physical drivers of storm track change are identified and their influence on the storm tracks is assessed. The experimental design aims to perturb the different physical drivers independently, by magnitudes representative of the range of values present in the CMIP3 model runs, and this is achieved via perturbations to the sea surface temperature and the sea-ice concentration forcing fields. We ask the question: can the spread of projections for the extratropical storm tracks present in the CMIP3 models be accounted for in a simple way by any of the identified drivers? The results suggest that, whilst the changes in the upper-tropospheric equator-to-pole temperature difference have an influence on the storm track response to climate change, the large spread of projections for the extratropical storm track present in the northern North Atlantic in particular is more strongly associated with changes in the lower-tropospheric equator-to-pole temperature difference.