Terrestrial biosphere changes over the last 120 kyr and their impact on ocean δ 13C

A new global synthesis and biomization of long (>40 kyr) pollen-data records is presented, and used with simulations from the HadCM3 and FAMOUS climate models to analyse the dynamics of the global terrestrial biosphere and carbon storage over the last glacial–interglacial cycle. Global modelled (BIOME4) biome distributions over time generally agree well with those inferred from pollen data. The two climate models show good agreement in global net primary productivity (NPP). NPP is strongly influenced by atmospheric carbon dioxide (CO2) concentrations through CO2 fertilization. The combined effects of modelled changes in vegetation and (via a simple model) soil carbon result in a global terrestrial carbon storage at the Last Glacial Maximum that is 210–470 Pg C less than in pre-industrial time. Without the contribution from exposed glacial continental shelves the reduction would be larger, 330–960 Pg C. Other intervals of low terrestrial carbon storage include stadial intervals at 108 and 85 ka BP, and between 60 and 65 ka BP during Marine Isotope Stage 4. Terrestrial carbon storage, determined by the balance of global NPP and decomposition, influences the stable carbon isotope composition (δ13C) of seawater because terrestrial organic carbon is depleted in 13C. Using a simple carbon-isotope mass balance equation we find agreement in trends between modelled ocean δ13C based on modelled land carbon storage, and palaeo-archives of ocean δ13C, confirming that terrestrial carbon storage variations may be important drivers of ocean δ13C changes.

Recent lake-level and outflow variations at Lake Viljandi, Estonia: validation of a coupled lake-catchment modelling scheme for climate change studies

Changes in the depth of Lake Viljandi between 1940 and 1990 were simulated using a lake water and energy-balance model driven by standard monthly weather data. Catchment runoff was simulated using a one-dimensional hydrological model, with a two-layer soil, a single-layer snowpack, a simple representation of vegetation cover and similarly modest input requirements. Outflow was modelled as a function of lake level. The simulated record of lake level and outflow matched observations of lake-level variations (r = 0.78) and streamflow (r = 0.87) well. The ability of the model to capture both intra- and inter-annual variations in the behaviour of a specific lake, despite the relatively simple input requirements, makes it extremely suitable for investigations of the impacts of climate change on lake water balance.

Changes in fire regimes since the Last Glacial Maximum: an assessment based on a global synthesis and analysis of charcoal data

Fire activity has varied globally and continuously since the last glacial maximum (LGM) in response to long-term changes in global climate and shorter-term regional changes in climate, vegetation, and human land use. We have synthesized sedimentary charcoal records of biomass burning since the LGM and present global maps showing changes in fire activity for time slices during the past 21,000 years (as differences in charcoal accumulation values compared to pre-industrial). There is strong broad-scale coherence in fire activity after the LGM, but spatial heterogeneity in the signals increases thereafter. In North America, Europe and southern South America, charcoal records indicate less-than-present fire activity during the deglacial period, from 21,000 to ∼11,000 cal yr BP. In contrast, the tropical latitudes of South America and Africa show greater-than-present fire activity from ∼19,000 to ∼17,000 cal yr BP and most sites from Indochina and Australia show greater-than-present fire activity from 16,000 to ∼13,000 cal yr BP. Many sites indicate greater-than-present or near-present activity during the Holocene with the exception of eastern North America and eastern Asia from 8,000 to ∼3,000 cal yr BP, Indonesia and Australia from 11,000 to 4,000 cal yr BP, and southern South America from 6,000 to 3,000 cal yr BP where fire activity was less than present. Regional coherence in the patterns of change in fire activity was evident throughout the post-glacial period. These complex patterns can largely be explained in terms of large-scale climate controls modulated by local changes in vegetation and fuel load

Climate versus carbon dioxide controls on biomass burning: a model analysis of the glacial–interglacial contrast

Climate controls fire regimes through its influence on the amount and types of fuel present and their dryness. CO2 concentration constrains primary production by limiting photosynthetic activity in plants. However, although fuel accumulation depends on biomass production, and hence on CO2 concentration, the quantitative relationship between atmospheric CO2 concentration and biomass burning is not well understood. Here a fire-enabled dynamic global vegetation model (the Land surface Processes and eXchanges model, LPX) is used to attribute glacial–interglacial changes in biomass burning to an increase in CO2, which would be expected to increase primary production and therefore fuel loads even in the absence of climate change, vs. climate change effects. Four general circulation models provided last glacial maximum (LGM) climate anomalies – that is, differences from the pre-industrial (PI) control climate – from the Palaeoclimate Modelling Intercomparison Project Phase~2, allowing the construction of four scenarios for LGM climate. Modelled carbon fluxes from biomass burning were corrected for the model's observed prediction biases in contemporary regional average values for biomes. With LGM climate and low CO2 (185 ppm) effects included, the modelled global flux at the LGM was in the range of 1.0–1.4 Pg C year-1, about a third less than that modelled for PI time. LGM climate with pre-industrial CO2 (280 ppm) yielded unrealistic results, with global biomass burning fluxes similar to or even greater than in the pre-industrial climate. It is inferred that a substantial part of the increase in biomass burning after the LGM must be attributed to the effect of increasing CO2 concentration on primary production and fuel load. Today, by analogy, both rising CO2 and global warming must be considered as risk factors for increasing biomass burning. Both effects need to be included in models to project future fire risks.

Orbital and millenial-scale environmental changes between 64 and 25 ka BP recorded in Black Sea sediments

High-resolution pollen and dinoflagellate cyst records from sediment core M72/5-25-GC1 were used to reconstruct vegetation dynamics in northern Anatolia and surface conditions of the Black Sea between 64 and 20 ka BP. During this period, the dominance of Artemisia in the pollen record indicates a steppe landscape and arid climate conditions. However, the concomitant presence of temperate arboreal pollen suggests the existence of glacial refugia in northern Anatolia. Long-term glacial vegetation dynamics reveal two major arid phases ~64–55 and 40–32 ka BP, and two major humid phases ~54–45 and 28–20 ka BP, correlating with higher and lower summer insolation, respectively. Dansgaard–Oeschger (D–O) cycles are clearly indicated by the 25-GC1 pollen record. Greenland interstadials are characterized by a marked increase in temperate tree pollen, indicating a spread of forests due to warm/wet conditions in northern Anatolia, whereas Greenland stadials reveal cold and arid conditions as indicated by spread of xerophytic biomes. There is evidence for a phase lag of ~500 to 1500 yr between initial warming and forest expansion, possibly due to successive changes in atmospheric circulation in the North Atlantic sector. The dominance of Pyxidinopsis psilata and Spiniferites cruciformis in the dinocyst record indicates brackish Black Sea conditions during the entire glacial period. The decrease of marine indicators (marine dinocysts, acritarchs) at ~54 ka BP and increase of freshwater algae (Pediastrum, Botryococcus) from 32 to 25 ka BP reveals freshening of the Black Sea surface water. This freshening is possibly related to humid phases in the region, to connection between Caspian Sea and Black Sea, to seasonal freshening by floating ice, and/or to closer position of river mouths due to low sea level. In the southern Black Sea, Greenland interstadials are clearly indicated by high dinocyst concentrations and calcium carbonate content, as a result of an increase in primary productivity. Heinrich events show a similar impact on the environment in the northern Anatolia/Black Sea region as Greenland stadials.

Palaeoenvironmental changes during the Valanginian: new insights from variations in phosphorus contents and bulk- and clay mineralogies in the western Tethys

Paleoenvironmental and paleoclimatic changes during the Valanginian carbon isotopic excursion (CIE) have been investigated in the western Tethys. For this purpose, bulk-rock and clay mineralogies, as well as phosphorus (P) contents were evaluated in a selection of five sections located in the Vocontian Basin (Angles, SE France; Alvier, E Switzerland; Malleval, E France), and the Lombardian Basin (Capriolo, N Italy; Breggia, S Switzerland). Within the CIE interval, bulk-rock and clay mineralogies are inferred to reflect mostly climate change. The onset of the CIE (Busnardoites campylotoxus ammonite Zone) is characterized by higher detrital index (DI: sum of the detrital minerals divided by calcite contents) values and the presence of kaolinite in their clay-mineral assemblages. In the late Valanginian (from the Saynoceras verrucosum Zone up to the end of the Valanginian), the samples show relatively variable DI and lower values or the absence of kaolinite. The variation in the mineralogical composition is interpreted as reflecting a change from a climate characterized by optimal weathering conditions associated with an increase in terrigenous input on the southern European margin during the CIE towards an overall unstable climate associated with drier conditions in the late Valanginian. This is contrasted by a dissymmetry (proximal vs distal) along the studied transect, the northern Tethyan margin being more sensitive to changes in continental input compared to the distal environments. P accumulation rates (PAR) present similar features. In the Vocontian basin, P content variations are associated with changes in terrigenous influx, whereas in the Lombardian basin (i.e. Capriolo and Breggia), PAR values are less well correlated. This is mainly because the deeper part of the Tethys was less sensitive to changes in continental inputs. The onset of the CIE (top of the B. campylotoxus Zone) records a general increase in PAR suggesting an increase in marine nutrient levels. This is linked to higher continental weathering rates and the enhanced influx of nutrients into the ocean. In the period corresponding to the shift itself, P contents show a dissymmetry between the Vocontian and Lombardian basins (proximal vs distal). For the sections of Malleval, Alvier and Angles, a decrease in P concentrations associated to a decrease in detrital input is observed. In Capriolo and Breggia, PAR show maximum values during the plateau, indicating a more complex interaction between different P sources. The time interval including the top of S. verrucosum Zone up to the end of the Valanginian is characterized by variable PAR values, suggesting variable nutrient influxes. These changes are in agreement with an evolution towards seasonally contrasted conditions in the late Valanginian.

Two fast radiative transfer methods to improve the temporal sampling of clouds in numerical weather prediction and climate models

The high computational cost of calculating the radiative heating rates in numerical weather prediction (NWP) and climate models requires that calculations are made infrequently, leading to poor sampling of the fast-changing cloud field and a poor representation of the feedback that would occur. This paper presents two related schemes for improving the temporal sampling of the cloud field. Firstly, the ‘split time-stepping’ scheme takes advantage of the independent nature of the monochromatic calculations of the ‘correlated-k’ method to split the calculation into gaseous absorption terms that are highly dependent on changes in cloud (the optically thin terms) and those that are not (optically thick). The small number of optically thin terms can then be calculated more often to capture changes in the grey absorption and scattering associated with cloud droplets and ice crystals. Secondly, the ‘incremental time-stepping’ scheme uses a simple radiative transfer calculation using only one or two monochromatic calculations representing the optically thin part of the atmospheric spectrum. These are found to be sufficient to represent the heating rate increments caused by changes in the cloud field, which can then be added to the last full calculation of the radiation code. We test these schemes in an operational forecast model configuration and find a significant improvement is achieved, for a small computational cost, over the current scheme employed at the Met Office. The ‘incremental time-stepping’ scheme is recommended for operational use, along with a new scheme to correct the surface fluxes for the change in solar zenith angle between radiation calculations.

Vegetation and environmental changes in Northern Anatolia between 134 and 119ka recorded in Black Sea sediments

This multiproxy study on SE Black Sea sediments provides the first detailed reconstruction of vegetation and environmental history of Northern Anatolia between 134 and 119 ka. Here, the glacial–interglacial transition is characterized by several short-lived alternating cold and warm events preceding a meltwater pulse (~ 130.4–131.7 ka). The latter is reconstructed as a cold arid period correlated to Heinrich event 11. The initial warming is evidenced at ~ 130.4 ka by increased primary productivity in the Black Sea, disappearance of ice-rafted detritus, and spreading of oaks in Anatolia. A Younger Dryas-type event is not identifiable. The Eemian vegetation succession corresponds to the main climatic phases in Europe: i) the Quercus–Juniperus phase (128.7–126.4 ka) indicates a dry continental climate; ii) the Ostrya–Corylus–Quercus–Carpinus phase (126.4–122.9 ka) suggests warm summers, mild winters, and high year-round precipitation; iii) the Fagus–Carpinus phase (122.9–119.5 ka) indicates cooling and high precipitation; and iv) increasing Pinus at ~ 121 ka marks the onset of cooler/drier conditions. Generally, pollen reconstructions suggest altitudinal/latitudinal migrations of vegetation belts in Northern Anatolia during the Eemian caused by increased transport of moisture. The evidence for the wide distribution of Fagus around the Black Sea contrasts with the European records and is likely related to climatic and genetic factors

Quantifying the relative importance of land cover change from climate and land-use in the representative concentration pathways

Climate change is projected to cause substantial alterations in vegetation distribution, but these have been given little attention in comparison to land-use in the Representative Concentration Pathway (RCP) scenarios. Here we assess the climate-induced land cover changes (CILCC) in the RCPs, and compare them to land-use land cover change (LULCC). To do this, we use an ensemble of simulations with and without LULCC in earth system model HadGEM2-ES for RCP2.6, RCP4.5 and RCP8.5. We find that climate change causes an expansion poleward of vegetation that affects more land area than LULCC in all of the RCPs considered here. The terrestrial carbon changes from CILCC are also larger than for LULCC. When considering only forest, the LULCC is larger, but the CILCC is highly variable with the overall radiative forcing of the scenario. The CILCC forest increase compensates 90% of the global anthropogenic deforestation by 2100 in RCP8.5, but just 3% in RCP2.6. Overall, bigger land cover changes tend to originate from LULCC in the shorter term or lower radiative forcing scenarios, and from CILCC in the longer term and higher radiative forcing scenarios. The extent to which CILCC could compensate for LULCC raises difficult questions regarding global forest and biodiversity offsetting, especially at different timescales. This research shows the importance of considering the relative size of CILCC to LULCC, especially with regard to the ecological effects of the different RCPs.

Climate responses to anthropogenic emissions of short-lived climate pollutants

Policies to control air quality focus on mitigating emissions of aerosols and their precursors, and other short-lived climate pollutants (SLCPs). On a local scale, these policies will have beneficial impacts on health and crop yields, by reducing particulate matter (PM) and surface ozone concentrations; however, the climate impacts of reducing emissions of SLCPs are less straightforward to predict. In this paper we consider a set of idealised, extreme mitigation strategies, in which the total anthropogenic emissions of individual SLCP emissions species are removed. This provides an upper bound on the potential climate impacts of such air quality strategies. We focus on evaluating the climate responses to changes in anthropogenic emissions of aerosol precursor species: black carbon (BC), organic carbon (OC) and sulphur dioxide (SO2). We perform climate integrations with four fully coupled atmosphere-ocean global climate models (AOGCMs), and examine the effects on global and regional climate of removing the total land-based anthropogenic emissions of each of the three aerosol precursor species. We find that the SO2 emissions reductions lead to the strongest response, with all three models showing an increase in surface temperature focussed in the northern hemisphere high latitudes, and a corresponding increase in global mean precipitation and run-off. Changes in precipitation and run-off patterns are driven mostly by a northward shift in the ITCZ, consistent with the hemispherically asymmetric warming pattern driven by the emissions changes. The BC and OC emissions reductions give a much weaker forcing signal, and there is some disagreement between models in the sign of the climate responses to these perturbations. These differences between models are due largely to natural variability in sea-ice extent, circulation patterns and cloud changes. This large natural variability component to the signal when the ocean circulation and sea-ice are free-running means that the BC and OC mitigation measures do not necessarily lead to a discernible climate response.

Impacts of uplift of northern Tibetan Plateau and formation of Asian inland deserts on regional climate and environment

Based on the geological evidence that the northern Tibetan Plateau (NTP) had an uplift of a finite magnitude since the Miocene and the major Asian inland deserts formed in the early Pliocene, a regional climate model (RegCM4.1) with a horizontal resolution of 50 km was used to explore the effects of the NTP uplift and the related aridification of inland Asia on regional climate. We designed three numerical experiments including the control experiment representing the present-day condition, the high-mountain experiment representing the early Pliocene condition with uplifted NTP but absence of the Asian inland deserts, and the low-mountain experiment representing the mid-Miocene condition with reduced topography in the NTP (by as much as 2400 m) and also absence of the deserts. Our simulation results indicated that the NTP uplift caused significant reductions in annual precipitation in a broad region of inland Asia north of the Tibetan Plateau (TP) mainly due to the enhanced rain shadow effect of the mountains and changes in the regional circulations. However, four mountainous regions located in the uplift showed significant increases in precipitation, stretching from the Pamir Plateau in the west to the Qilian Mountains in the east. These mountainous areas also experienced different changes in the rainfall seasonality with the greatest increases occurring during the respective rainy seasons, predominantly resulted from the enhanced orographically forced upwind ascents. The appearance of the major deserts in the inland Asia further reduced precipitation in the region and led to increased dust emission and deposition fluxes, while the spatial patterns of dust deposition were also changed, not only in the regions of uplift-impacted topography, but also in the downwind regions. One major contribution from this study is the comparison of the simulation results with 11 existing geological records representing the moisture conditions from Miocene to Pliocene. The comparisons revealed good matches between the simulation results and the published geological records. Therefore, we conclude that the NTP uplift and the related formation of the major deserts played a controlling role in the evolution of regional climatic conditions in a broad region in inland Asia since the Miocene.

Radiative forcing and climate change

Aviation causes climate change as a result of its emissions of CO2, oxides of nitrogen, aerosols, and water vapor. One simple method of quantifying the climate impact of past emissions is radiative forcing. The radiative forcing due to changes in CO2 is best characterized, but there are formidable difficulties in estimating the non-CO2 forcings – this is particularly the case for possible aviation-induced changes in cloudiness (AIC). The most recent comprehensive assessment gave a best estimate of the 2005 total radiative forcing due to aviation of about 55–78 mW m−2 depending on whether AIC was included or not, with an uncertainty of at least a factor of 2. The aviation CO2 radiative forcing represents about 1.6% of the total CO2 forcing from all human activities. It is estimated that, including the non-CO2 effects, aviation contributes between 1.3 and 14% of the total radiative forcing due to all human activities. Alternative methods for comparing the future impact of present-day aviation emissions are presented – the perception of the relative importance of the non-CO2 emissions, relative to CO2, depends considerably on the chosen method and the parameters chosen within those methods.

Attribution of climate extreme events

There is a tremendous desire to attribute causes to weather and climate events that is often challenging from a physical standpoint. Headlines attributing an event solely to either human-induced climate change or natural variability can be misleading when both are invariably in play. The conventional attribution framework struggles with dynamically driven extremes because of the small signal-to-noise ratios and often uncertain nature of the forced changes. Here, we suggest that a different framing is desirable, which asks why such extremes unfold the way they do. Specifically, we suggest that it is more useful to regard the extreme circulation regime or weather event as being largely unaffected by climate change, and question whether known changes in the climate system's thermodynamic state affected the impact of the particular event. Some examples briefly illustrated include 'snowmaggedon' in February 2010, superstorm Sandy in October 2012 and supertyphoon Haiyan in November 2013, and, in more detail, the Boulder floods of September 2013, all of which were influenced by high sea surface temperatures that had a discernible human component.

A mechanism of internal decadal atlantic ocean variability in a high-resolution coupled climate model

The North Atlantic Ocean subpolar gyre (NA SPG) is an important region for initialising decadal climate forecasts. Climate model simulations and palaeo climate reconstructions have indicated that this region could also exhibit large, internally generated variability on decadal timescales. Understanding these modes of variability, their consistency across models, and the conditions in which they exist, is clearly important for improving the skill of decadal predictions — particularly when these predictions are made with the same underlying climate models. Here we describe and analyse a mode of internal variability in the NA SPG in a state-of-the-art, high resolution, coupled climate model. This mode has a period of 17 years and explains 15–30% of the annual variance in related ocean indices. It arises due to the advection of heat content anomalies around the NA SPG. Anomalous circulation drives the variability in the southern half of the NA SPG, whilst mean circulation and anomalous temperatures are important in the northern half. A negative feedback between Labrador Sea temperatures/densities and those in the North Atlantic Current is identified, which allows for the phase reversal. The atmosphere is found to act as a positive feedback on to this mode via the North Atlantic Oscillation which itself exhibits a spectral peak at 17 years. Decadal ocean density changes associated with this mode are driven by variations in temperature, rather than salinity — a point which models often disagree on and which we suggest may affect the veracity of the underlying assumptions of anomaly-assimilating decadal prediction methodologies.

Nonlinearity of ocean heat uptake during warming and cooling in the FAMOUS climate model

Atmospheric CO2 concentration is expected to continue rising in the coming decades, but natural or artificial processes may eventually reduce it. We show that, in the FAMOUS atmosphere-ocean general circulation model, the reduction of ocean heat content as radiative forcing decreases is greater than would be expected from a linear model simulation of the response to the applied forcings. We relate this effect to the behavior of the Atlantic meridional overturning circulation (AMOC): the ocean cools more efficiently with a strong AMOC. The AMOC weakens as CO2 rises, then strengthens as CO2 declines, but temporarily overshoots its original strength. This nonlinearity comes mainly from the accumulated advection of salt into the North Atlantic, which gives the system a longer memory. This implies that changes observed in response to different CO2 scenarios or from different initial states, such as from past changes, may not be a reliable basis for making projections.