Full effects of land use change in the representative concentration pathways

Future land use change (LUC) is an important component of the IPCC representative concentration pathways (RCPs), but in these scenarios' radiative forcing targets the climate impact of LUC only includes greenhouse gases. However, climate effects due to physical changes of the land surface can be as large. Here we show the critical importance of including non-carbon impacts of LUC when considering the RCPs. Using an ensemble of climate model simulations with and without LUC, we show that the net climate effect is very different from the carbon-only effect. Despite opposite signs of LUC, all the RCPs assessed here have a small net warming from LUC because of varying biogeophysical effects, and in RCP4.5 the warming is outside of the expected variability. The afforestation in RCP4.5 decreases surface albedo, making the net global temperature anomaly over land around five times larger than RCPs 2.6 and 8.5, for around twice the amount of LUC. Consequent changes to circulation in RCP4.5 in turn reduce Arctic sea ice cover. The small net positive temperature effect from LUC could make RCP4.5's universal carbon tax, which incentivizes retaining and growing forest, counter productive with respect to climate. However, there are spatial differences in the balance of impacts, and potential climate gains would need to be assessed against other environmental aims.

El Nino-Southern Oscillation

This 20th anniversary edition of PAGES news explores the elusive El Niño-Southern Oscillation (ENSO) from a paleo-perspective. The initiative for this special section came out of a 2011 PMIP workshop. In a mini section, a newly introduced format in the PAGES newsletter, four articles on data assimilation address methods of combining observations and model simulations. Finally, this newsletter also features several reports on the 2013 2nd Young Scientists Meeting held in Goa, India.

The climate of Europe 6000 years ago

Palaeoclimates across Europe for 6000 y BP were estimated from pollen data using the modern pollen analogue technique constrained with lake-level data. The constraint consists of restricting the set of modern pollen samples considered as analogues of the fossil samples to those locations where the implied change in annual precipitation minus evapotranspiration (P–E) is consistent with the regional change in moisture balance as indicated by lakes. An artificial neural network was used for the spatial interpolation of lake-level changes to the pollen sites, and for mapping palaeoclimate anomalies. The climate variables reconstructed were mean temperature of the coldest month (T c ), growing degree days above 5  °C (GDD), moisture availability expressed as the ratio of actual to equilibrium evapotranspiration (α), and P–E. The constraint improved the spatial coherency of the reconstructed palaeoclimate anomalies, especially for P–E. The reconstructions indicate clear spatial and seasonal patterns of Holocene climate change, which can provide a quantitative benchmark for the evaluation of palaeoclimate model simulations. Winter temperatures (T c ) were 1–3 K greater than present in the far N and NE of Europe, but 2–4 K less than present in the Mediterranean region. Summer warmth (GDD) was greater than present in NW Europe (by 400–800 K day at the highest elevations) and in the Alps, but >400 K day less than present at lower elevations in S Europe. P–E was 50–250 mm less than present in NW Europe and the Alps, but α was 10–15% greater than present in S Europe and P–E was 50–200 mm greater than present in S and E Europe.

Vegetation and soils feedbacks on the response of the African monsoon response to orbital forcing in early to middle Holocene

Fossil pollen, ancient lake sediments and archaeological evidence from Africa indicate that the Sahel and Sahara regions were considerably wetter than today during the early to middle Holocene period, about 12,000 to 5,000 years ago1–4. Vegetation associated with the modern Sahara/Sahel boundary was about 5° farther north, and there were more and larger lakes between 15 and 30° N. Simulations with climate models have shown that these wetter conditions were probably caused by changes in Earth's orbital parameters that increased the amplitude of the seasonal cycle of solar radiation in the Northern Hemisphere, enhanced the land-ocean temperature contrast, and thereby strengthened the African summer monsoon5–7. However, these simulations underestimated the consequent monsoon enhancement as inferred from palaeorecords4. Here we use a climate model to show that changes in vegetation and soil may have increased the climate response to orbital forcing. We find that replacing today's orbital forcing with that of the mid-Holocene increases summer precipitation by 12% between 15 and 22° N. Replacing desert with grassland, and desert soil with more loamy soil, further enhances the summer precipitation (by 6 and 10% respectively), giving a total precipitation increase of 28%. When the simulated climate changes are applied to a biome model, vegetation becomes established north of the current Sahara/Sahel boundary, thereby shrinking the area of the Sahara by 11% owing to orbital forcing alone, and by 20% owing to the combined influence of orbital forcing and the prescribed vegetation and soil changes. The inclusion of the vegetation and soil feedbacks thus brings the model simulations and palaeovegetation observations into closer agreement.

Potential role of vegetation feedback in the climate sensitivity of high-latitude regions: a case study at 6000 years B.P.

Previous climate model simulations have shown that the configuration of the Earth's orbit during the early to mid-Holocene (approximately 10–5 kyr) can account for the generally warmer-than-present conditions experienced by the high latitudes of the northern hemisphere. New simulations for 6 kyr with two atmospheric/mixed-layer ocean models (Community Climate Model, version 1, CCMl, and Global ENvironmental and Ecological Simulation of Interactive Systems, version 2, GENESIS 2) are presented here and compared with results from two previous simulations with GENESIS 1 that were obtained with and without the albedo feedback due to climate-induced poleward expansion of the boreal forest. The climate model results are summarized in the form of potential vegetation maps obtained with the global BIOME model, which facilitates visual comparisons both among models and with pollen and plant macrofossil data recording shifts of the forest-tundra boundary. A preliminary synthesis shows that the forest limit was shifted 100–200 km north in most sectors. Both CCMl and GENESIS 2 produced a shift of this magnitude. GENESIS 1 however produced too small a shift, except when the boreal forest albedo feedback was included. The feedback in this case was estimated to have amplified forest expansion by approximately 50%. The forest limit changes also show meridional patterns (greatest expansion in central Siberia and little or none in Alaska and Labrador) which have yet to be reproduced by models. Further progress in understanding of the processes involved in the response of climate and vegetation to orbital forcing will require both the deployment of coupled atmosphere-biosphere-ocean models and the development of more comprehensive observational data sets

Middle and Late Pleistocene humid periods recorded in palaeolake deposits of the Nafud desert, Saudi Arabia

Present climate in the Nafud desert of northern Saudi Arabia is hyper-arid and moisture brought by north-westerly winds scarcely reaches the region. The existence of abundant palaeolake sediments provides evidence for a considerably wetter climate in the past. However, the existing chronological framework of these deposits is solely based on radiocarbon dating of questionable reliability, due to potential post-depositional contamination with younger 14C. By using luminescence dating, we show that the lake deposits were not formed between 40 and 20 ka as suggested previously, but approximately ca 410 ka, 320 ka, 200 ka, 125 ka, and 100 ka ago. All of these humid phases are in good agreement with those recorded in lake sediments and speleothems from southern Arabia. Surprisingly, no Holocene lake deposits were identified. Geological characteristics of the deposits and diatom analysis suggest that a single, perennial lake covered the entire south-western Nafud ca 320 ka ago. In contrast, lakes of the 200 ka, 125 ka, and 100 ka humid intervals were smaller and restricted to interdune depressions of a pre-existing dune relief. The concurrent occurrence of humid phases in the Nafud, southern Arabia and the eastern Mediterranean suggests that moisture in northern Arabia originated either from the Mediterranean due to more frequent frontal depression systems or from stronger Indian monsoon circulation, respectively. However, based on previously published climate model simulations and palaecolimate evidence from central Arabia and the Negev desert, we argue that humid climate conditions in the Nafud were probably caused by a stronger African monsoon and a distinct change in zonal atmospheric circulation.

Probabilistic downscaling of remote sensing data with applications for multi-scale biogeochemical flux modeling

Upscaling ecological information to larger scales in space and downscaling remote sensing observations or model simulations to finer scales remain grand challenges in Earth system science. Downscaling often involves inferring subgrid information from coarse-scale data, and such ill-posed problems are classically addressed using regularization. Here, we apply two-dimensional Tikhonov Regularization (2DTR) to simulate subgrid surface patterns for ecological applications. Specifically, we test the ability of 2DTR to simulate the spatial statistics of high-resolution (4 m) remote sensing observations of the normalized difference vegetation index (NDVI) in a tundra landscape. We find that the 2DTR approach as applied here can capture the major mode of spatial variability of the high-resolution information, but not multiple modes of spatial variability, and that the Lagrange multiplier (γ) used to impose the condition of smoothness across space is related to the range of the experimental semivariogram. We used observed and 2DTR-simulated maps of NDVI to estimate landscape-level leaf area index (LAI) and gross primary productivity (GPP). NDVI maps simulated using a γ value that approximates the range of observed NDVI result in a landscape-level GPP estimate that differs by ca 2% from those created using observed NDVI. Following findings that GPP per unit LAI is lower near vegetation patch edges, we simulated vegetation patch edges using multiple approaches and found that simulated GPP declined by up to 12% as a result. 2DTR can generate random landscapes rapidly and can be applied to disaggregate ecological information and compare of spatial observations against simulated landscapes.

Efficient incorporation of channel cross-section geometry uncertainty into regional and global scale flood inundation models

This paper investigates the challenge of representing structural differences in river channel cross-section geometry for regional to global scale river hydraulic models and the effect this can have on simulations of wave dynamics. Classically, channel geometry is defined using data, yet at larger scales the necessary information and model structures do not exist to take this approach. We therefore propose a fundamentally different approach where the structural uncertainty in channel geometry is represented using a simple parameterization, which could then be estimated through calibration or data assimilation. This paper first outlines the development of a computationally efficient numerical scheme to represent generalised channel shapes using a single parameter, which is then validated using a simple straight channel test case and shown to predict wetted perimeter to within 2% for the channels tested. An application to the River Severn, UK is also presented, along with an analysis of model sensitivity to channel shape, depth and friction. The channel shape parameter was shown to improve model simulations of river level, particularly for more physically plausible channel roughness and depth parameter ranges. Calibrating channel Manning’s coefficient in a rectangular channel provided similar water level simulation accuracy in terms of Nash-Sutcliffe efficiency to a model where friction and shape or depth were calibrated. However, the calibrated Manning coefficient in the rectangular channel model was ~2/3 greater than the likely physically realistic value for this reach and this erroneously slowed wave propagation times through the reach by several hours. Therefore, for large scale models applied in data sparse areas, calibrating channel depth and/or shape may be preferable to assuming a rectangular geometry and calibrating friction alone.

Modelling primary producer interaction and composition: an example of a UK lowland river

Nutrient enrichment and drought conditions are major threats to lowland rivers causing ecosystem degradation and composition changes in plant communities. The controls on primary producer composition in chalk rivers are investigated using a new model and existing data from the River Frome (UK) to explore abiotic and biotic interactions. The growth and interaction of four primary producer functional groups (suspended algae, macrophytes, epiphytes, sediment biofilm) were successfully linked with flow, nutrients (N, P), light and water temperature such that the modelled biomass dynamics of the four groups matched that of the observed. Simulated growth of suspended algae was limited mainly by the residence time of the river rather than in-stream phosphorus concentrations. The simulated growth of the fixed vegetation (macrophytes, epiphytes, sediment biofilm) was overwhelmingly controlled by incoming solar radiation and light attenuation in the water column. Nutrients and grazing have little control when compared to the other physical controls in the simulations. A number of environmental threshold values were identified in the model simulations for the different producer types. The simulation results highlighted the importance of the pelagic–benthic interactions within the River Frome and indicated that process interaction defined the behaviour of the primary producers, rather than a single, dominant driver. The model simulations pose interesting questions to be considered in the next iteration of field- and laboratory based studies.

Interhemispheric dynamics of the African rainbelt during the late Quaternary

The spatial pattern of precipitation variability in tropical and subtropical Africa over the late Quaternary has long been debated. Prevailing hypotheses variously infer (1) insolation-controlled asymmetry of wet phases between hemispheres, (2) symmetric contraction and expansion of the tropical rainbelt, and (3) independent control on moisture available in Southern Africa via sea surface temperatures in the Indian Ocean. In this study we use climate-model simulations covering the last glacial cycle (120 kyr) with HadCM3 and the multi-model ensembles from PMIP3 (the Palaeoclimate Model Intercomparison Project) to investigate the long-term behaviour of the African rainbelt, and test these simulations against existing empirical palaeohydrological records. Through regional model-data comparisons we find evidence for the validity of several hypotheses, with various proposed processes occurring concurrently but with different regional emphasis (e.g. asymmetric shifts at the seasonal extremes and symmetric expansions/ contractions towards West equatorial regions). Crucially, variations in rainfall are associated with multiple forcing mechanisms that vary in their dominance both spatially and temporally over the glacial cycle; an important consideration when interpreting and extrapolating from often relatively short palaeoenvironmental records.

Extratropical cyclones in a warmer, moister climate: a recent Atlantic analogue

Current climate model projections do not exhibit a large change in the intensity of extratropical cyclones. However, there are concerns that current models represent moist processes poorly, and this provides motivation for investigating observational evidence for how cyclones behave in warmer climates. In the North Atlantic in particular, recent decades provide a clear contrast between warm and cold climates due to Atlantic Multidecadal Variability. In this paper we investigate these periods as analogues which may provide a guide to future cyclone behavior. While temperature and moisture rise in recent warm periods as in the projections, differences in energetics and temperature gradients imply that these periods are only partial analogues. The main result from current reanalyses is that while increased cyclone-associated precipitation is seen in the recent warm periods, there is no robust evidence of an increase in cyclone intensity by other measures, such as maximum wind speed or vorticity. A set of low- and high-resolution model simulations are also studied, suggesting that changes in cyclone intensity may be different in higher-resolution reanalyses.

Combining satellite observations and reanalysis energy transports to estimate global net surface energy fluxes 1985-2012

Two methods are developed to estimate net surface energy fluxes based upon satellite-based reconstructions of radiative fluxes at the top of atmosphere and the atmospheric energy tendencies and transports from the ERA-Interim reanalysis. Method 1 applies the mass adjusted energy divergence from ERA-Interim while method 2 estimates energy divergence based upon the net energy difference at the top of atmosphere and the surface from ERA-Interim. To optimise the surface flux and its variability over ocean, the divergences over land are constrained to match the monthly area mean surface net energy flux variability derived from a simple relationship between the surface net energy flux and the surface temperature change. The energy divergences over the oceans are then adjusted to remove an unphysical residual global mean atmospheric energy divergence. The estimated net surface energy fluxes are compared with other data sets from reanalysis and atmospheric model simulations. The spatial correlation coefficients of multi-annual means between the estimations made here and other data sets are all around 0.9. There are good agreements in area mean anomaly variability over the global ocean, but discrepancies in the trend over the eastern Pacific are apparent.

Evaluating the climate and air quality impacts of short-lived pollutants

This paper presents a summary of the work done within the European Union's Seventh Framework Programme project ECLIPSE (Evaluating the Climate and Air Quality Impacts of Short-Lived Pollutants). ECLIPSE had a unique systematic concept for designing a realistic and effective mitigation scenario for short-lived climate pollutants (SLCPs; methane, aerosols and ozone, and their precursor species) and quantifying its climate and air quality impacts, and this paper presents the results in the context of this overarching strategy. The first step in ECLIPSE was to create a new emission inventory based on current legislation (CLE) for the recent past and until 2050. Substantial progress compared to previous work was made by including previously unaccounted types of sources such as flaring of gas associated with oil production, and wick lamps. These emission data were used for present-day reference simulations with four advanced Earth system models (ESMs) and six chemistry transport models (CTMs). The model simulations were compared with a variety of ground-based and satellite observational data sets from Asia, Europe and the Arctic. It was found that the models still underestimate the measured seasonality of aerosols in the Arctic but to a lesser extent than in previous studies. Problems likely related to the emissions were identified for northern Russia and India, in particular. To estimate the climate impacts of SLCPs, ECLIPSE followed two paths of research: the first path calculated radiative forcing (RF) values for a large matrix of SLCP species emissions, for different seasons and regions independently. Based on these RF calculations, the Global Temperature change Potential metric for a time horizon of 20 years (GTP20) was calculated for each SLCP emission type. This climate metric was then used in an integrated assessment model to identify all emission mitigation measures with a beneficial air quality and short-term (20-year) climate impact. These measures together defined a SLCP mitigation (MIT) scenario. Compared to CLE, the MIT scenario would reduce global methane (CH4) and black carbon (BC) emissions by about 50 and 80 %, respectively. For CH4, measures on shale gas production, waste management and coal mines were most important. For non-CH4 SLCPs, elimination of high-emitting vehicles and wick lamps, as well as reducing emissions from gas flaring, coal and biomass stoves, agricultural waste, solvents and diesel engines were most important. These measures lead to large reductions in calculated surface concentrations of ozone and particulate matter. We estimate that in the EU, the loss of statistical life expectancy due to air pollution was 7.5 months in 2010, which will be reduced to 5.2 months by 2030 in the CLE scenario. The MIT scenario would reduce this value by another 0.9 to 4.3 months. Substantially larger reductions due to the mitigation are found for China (1.8 months) and India (11–12 months). The climate metrics cannot fully quantify the climate response. Therefore, a second research path was taken. Transient climate ensemble simulations with the four ESMs were run for the CLE and MIT scenarios, to determine the climate impacts of the mitigation. In these simulations, the CLE scenario resulted in a surface temperature increase of 0.70 ± 0.14 K between the years 2006 and 2050. For the decade 2041–2050, the warming was reduced by 0.22 ± 0.07 K in the MIT scenario, and this result was in almost exact agreement with the response calculated based on the emission metrics (reduced warming of 0.22 ± 0.09 K). The metrics calculations suggest that non-CH4 SLCPs contribute ~ 22 % to this response and CH4 78 %. This could not be fully confirmed by the transient simulations, which attributed about 90 % of the temperature response to CH4 reductions. Attribution of the observed temperature response to non-CH4 SLCP emission reductions and BC specifically is hampered in the transient simulations by small forcing and co-emitted species of the emission basket chosen. Nevertheless, an important conclusion is that our mitigation basket as a whole would lead to clear benefits for both air quality and climate. The climate response from BC reductions in our study is smaller than reported previously, possibly because our study is one of the first to use fully coupled climate models, where unforced variability and sea ice responses cause relatively strong temperature fluctuations that may counteract (and, thus, mask) the impacts of small emission reductions. The temperature responses to the mitigation were generally stronger over the continents than over the oceans, and with a warming reduction of 0.44 K (0.39–0.49) K the largest over the Arctic. Our calculations suggest particularly beneficial climate responses in southern Europe, where surface warming was reduced by about 0.3 K and precipitation rates were increased by about 15 (6–21) mm yr−1 (more than 4 % of total precipitation) from spring to autumn. Thus, the mitigation could help to alleviate expected future drought and water shortages in the Mediterranean area. We also report other important results of the ECLIPSE project.

What does global mean temperature tell us about local climate?

The subject of climate feedbacks focuses attention on global mean surface air temperature (GMST) as the key metric of climate change. But what does knowledge of past and future GMST tell us about the climate of specific regions? In the context of the ongoing UNFCCC process, this is an important question for policy-makers as well as for scientists. The answer depends on many factors, including the mechanisms causing changes, the timescale of the changes, and the variables and regions of interest. This paper provides a review and analysis of the relationship between changes in GMST and changes in local climate, first in observational records and then in a range of climate model simulations, which are used to interpret the observations. The focus is on decadal timescales, which are of particular interest in relation to recent and near-future anthropogenic climate change. It is shown that GMST primarily provides information about forced responses, but that understanding and quantifying internal variability is essential to projecting climate and climate impacts on regional-to-local scales. The relationship between local forced responses and GMST is often linear but may be nonlinear, and can be greatly complicated by competition between different forcing factors. Climate projections are limited not only by uncertainties in the signal of climate change but also by uncertainties in the characteristics of real-world internal variability. Finally, it is shown that the relationship between GMST and local climate provides a simple approach to climate change detection, and a useful guide to attribution studies.

Atmosphere-ocean coupled processes in the Madden-Julian oscillation

The Madden-Julian oscillation (MJO) is a convectively coupled 30-70 day (intraseasonal) tropical atmospheric mode that drives variations in global weather, but which is poorly simulated in most atmospheric general circulation models. Over the past two decades, field campaigns and modeling experiments have suggested that tropical atmosphere-ocean interactions may sustain or amplify the pattern of enhanced and suppressed atmospheric convection that defines the MJO, and encourage its eastward propagation through the Indian and Pacific Oceans. New observations collected during the past decade have advanced our understand of the ocean response to atmospheric MJO forcing and the resulting intraseasonal sea surface temperature (SST) fluctuations. Numerous modeling studies have revealed a considerable impact of the mean state on MJO ocean-atmosphere coupled processes, as well as the importance of resolving the diurnal cycle of atmosphere--upper-ocean interactions. New diagnostic methods provide insight to atmospheric variability and physical processes associated with the MJO, but offer limited insight on the role of ocean feedbacks. Consequently, uncertainty remains concerning the role of the ocean in MJO theory. Our understanding of how atmosphere-ocean coupled processes affect the MJO can be improved by collecting observations in poorly sampled regions of MJO activity, assessing oceanic and atmospheric drivers of surface fluxes, improving the representation of upper-ocean mixing in coupled-model simulations, designing model experiments that minimize mean-state differences, and developing diagnostic tools to evaluate the nature and role of coupled ocean-atmosphere processes over the MJO cycle.