Global biomass burning: a synthesis and review of Holocene paleofire records and their controls.

We synthesize existing sedimentary charcoal records to reconstruct Holocene fire history at regional, continental and global scales. The reconstructions are compared with the two potential controls of burning at these broad scales – changes in climate and human activities – to assess their relative importance on trends in biomass burning. Here we consider several hypotheses that have been advanced to explain the Holocene record of fire, including climate, human activities and synergies between the two. Our results suggest that 1) episodes of high fire activity were relatively common in the early Holocene and were consistent with climate changes despite low global temperatures and low levels of biomass burning globally; 2) there is little evidence from the paleofire record to support the Early Anthropocene Hypothesis of human modification of the global carbon cycle; 3) there was a nearly-global increase in fire activity from 3 to 2 ka that is difficult to explain with either climate or humans, but the widespread and synchronous nature of the increase suggests at least a partial climate forcing; and 4) burning during the past century generally decreased but was spatially variable; it declined sharply in many areas, but there were also large increases (e.g., Australia and parts of Europe). Our analysis does not exclude an important role for human activities on global biomass burning during the Holocene, but instead provides evidence for a pervasive influence of climate across multiple spatial and temporal scales.

Climatic control of the biomass-burning decline in the Americas after AD 1500

The significance and cause of the decline in biomass burning across the Americas after AD 1500 is a topic of considerable debate. We synthesized charcoal records (a proxy for biomass burning) from the Americas and from the remainder of the globe over the past 2000 years, and compared these with paleoclimatic records and population reconstructions. A distinct post-AD 1500 decrease in biomass burning is evident, not only in the Americas, but also globally, and both are similar in duration and timing to ‘Little Ice Age’ climate change. There is temporal and spatial variability in the expression of the biomass-burning decline across the Americas but, at a regional–continental scale, ‘Little Ice Age’ climate change was likely more important than indigenous population collapse in driving this decline.

How do plants respond to nutrient shortage by biomass allocation?

Plants constantly sense the changes in their environment; when mineral elements are scarce, they often allocate a greater proportion of their biomass to the root system. This acclimatory response is a consequence of metabolic changes in the shoot and an adjustment of carbohydrate transport to the root. It has long been known that deficiencies of essential macronutrients (nitrogen, phosphorus, potassium and magnesium) result in an accumulation of carbohydrates in leaves and roots, and modify the shoot-to-root biomass ratio. Here, we present an update on the effects of mineral deficiencies on the expression of genes involved in primary metabolism in the shoot, the evidence for increased carbohydrate concentrations and altered biomass allocation between shoot and root, and the consequences of these changes on the growth and morphology of the plant root system.

Building the niche through time: using 13,000 years of data to predict the effects of climate change on three tree species in Europe

Aim  Species distribution models (SDMs) based on current species ranges underestimate the potential distribution when projected in time and/or space. A multi-temporal model calibration approach has been suggested as an alternative, and we evaluate this using 13,000 years of data. Location  Europe. Methods  We used fossil-based records of presence for Picea abies, Abies alba and Fagus sylvatica and six climatic variables for the period 13,000 to 1000 yr bp. To measure the contribution of each 1000-year time step to the total niche of each species (the niche measured by pooling all the data), we employed a principal components analysis (PCA) calibrated with data over the entire range of possible climates. Then we projected both the total niche and the partial niches from single time frames into the PCA space, and tested if the partial niches were more similar to the total niche than random. Using an ensemble forecasting approach, we calibrated SDMs for each time frame and for the pooled database. We projected each model to current climate and evaluated the results against current pollen data. We also projected all models into the future. Results  Niche similarity between the partial and the total-SDMs was almost always statistically significant and increased through time. SDMs calibrated from single time frames gave different results when projected to current climate, providing evidence of a change in the species realized niches through time. Moreover, they predicted limited climate suitability when compared with the total-SDMs. The same results were obtained when projected to future climates. Main conclusions  The realized climatic niche of species differed for current and future climates when SDMs were calibrated considering different past climates. Building the niche as an ensemble through time represents a way forward to a better understanding of a species' range and its ecology in a changing climate.

Vertical distribution and radiative effects of mineral dust and biomass burning aerosol over West Africa during DABEX

This paper presents measurements of the vertical distribution of aerosol extinction coefficient over West Africa during the Dust and Biomass-burning Aerosol Experiment (DABEX)/African Monsoon Multidisciplinary Analysis dry season Special Observing Period Zero (AMMA-SOP0). In situ aircraft measurements from the UK FAAM aircraft have been compared with two ground-based lidars (POLIS and ARM MPL) and an airborne lidar on an ultralight aircraft. In general, mineral dust was observed at low altitudes (up to 2 km), and a mixture of biomass burning aerosol and dust was observed at altitudes of 2–5 km. The study exposes difficulties associated with spatial and temporal variability when intercomparing aircraft and ground measurements. Averaging over many profiles provided a better means of assessing consistent errors and biases associated with in situ sampling instruments and retrievals of lidar ratios. Shortwave radiative transfer calculations and a 3-year simulation with the HadGEM2-A climate model show that the radiative effect of biomass burning aerosol was somewhat sensitive to the vertical distribution of aerosol. In particular, when the observed low-level dust layer was included in the model, the absorption of solar radiation by the biomass burning aerosols increased by 10%. We conclude that this absorption enhancement was caused by the dust reflecting solar radiation up into the biomass burning aerosol layer. This result illustrates that the radiative forcing of anthropogenic absorbing aerosol can be sensitive to the presence of natural aerosol species.

The influence of vegetation, fire spread and fire behaviour on biomass burning and trace gas emissions: results from a process-based model

A process-based fire regime model (SPITFIRE) has been developed, coupled with ecosystem dynamics in the LPJ Dynamic Global Vegetation Model, and used to explore fire regimes and the current impact of fire on the terrestrial carbon cycle and associated emissions of trace atmospheric constituents. The model estimates an average release of 2.24 Pg C yr−1 as CO2 from biomass burning during the 1980s and 1990s. Comparison with observed active fire counts shows that the model reproduces where fire occurs and can mimic broad geographic patterns in the peak fire season, although the predicted peak is 1–2 months late in some regions. Modelled fire season length is generally overestimated by about one month, but shows a realistic pattern of differences among biomes. Comparisons with remotely sensed burnt-area products indicate that the model reproduces broad geographic patterns of annual fractional burnt area over most regions, including the boreal forest, although interannual variability in the boreal zone is underestimated.

Effective LAI and CHP of a single tree from small-footprint full-waveform LiDAR

This letter has tested the canopy height profile (CHP) methodology as a way of effective leaf area index (LAIe) and vertical vegetation profile retrieval at a single-tree level. Waveform and discrete airborne LiDAR data from six swaths, as well as from the combined data of six swaths, were used to extract the LAIe of a single live Callitris glaucophylla tree. LAIe was extracted from raw waveform as an intermediate step in the CHP methodology, with two different vegetation-ground reflectance ratios. Discrete point LAIe estimates were derived from the gap probability using the following: 1) single ground returns and 2) all ground returns. LiDAR LAIe retrievals were subsequently compared to hemispherical photography estimates, yielding mean values within ±7% of the latter, depending on the method used. The CHP of a single dead Callitris glaucophylla tree, representing the distribution of vegetation material, was verified with a field profile manually reconstructed from convergent photographs taken with a fixed-focal-length camera. A binwise comparison of the two profiles showed very high correlation between the data reaching R2 of 0.86 for the CHP from combined swaths. Using a study-area-adjusted reflectance ratio improved the correlation between the profiles, but only marginally in comparison to using an arbitrary ratio of 0.5 for the laser wavelength of 1550 nm.

How body mass and lifestyle affect juvenile biomass production in placental mammals

In mammals, the mass-specific rate of biomass production during gestation and lactation, here called maternal productivity, has been shown to vary with body size and lifestyle. Metabolic theory predicts that post-weaning growth of offspring, here termed juvenile productivity, should be higher than maternal productivity, and juveniles of smaller species should be more productive than those of larger species. Furthermore because juveniles generally have similar lifestyles to their mothers, across species juvenile and maternal productivities should be correlated. We evaluated these predictions with data from 270 species of placental mammals in 14 taxonomic/lifestyle groups. All three predictions were supported. Lagomorphs, perissodactyls and artiodactyls were very productive both as juveniles and as mothers as expected from the abundance and reliability of their foods. Primates and bats were unproductive as juveniles and as mothers, as expected as an indirect consequence of their low predation risk and consequent low mortality. Our results point the way to a mechanistic explanation for the suite of correlated life-history traits that has been called the slow–fast continuum.

Classification of seed storage behaviour of 67 Amazonian tree species

Information was collated on the seed storage behaviour of 67 tree species native to the Amazon rainforest of Brazil; 38 appeared to show orthodox, 23 recalcitrant and six intermediate seed storage behaviour. A double-criteria key based on thousand-seed weight and seed moisture content at shedding to estimate likely seed storage behaviour, developed previously, showed good agreement with the above classifications. The key can aid seed storage behaviour identification considerably.

Simulation of tree-ring widths with a model for primary production, carbon allocation, and growth

We present a simple, generic model of annual tree growth, called "T". This model accepts input from a first-principles light-use efficiency model (the "P" model). The P model provides values for gross primary production (GPP) per unit of absorbed photosynthetically active radiation (PAR). Absorbed PAR is estimated from the current leaf area. GPP is allocated to foliage, transport tissue, and fine-root production and respiration in such a way as to satisfy well-understood dimensional and functional relationships. Our approach thereby integrates two modelling approaches separately developed in the global carbon-cycle and forest-science literature. The T model can represent both ontogenetic effects (the impact of ageing) and the effects of environmental variations and trends (climate and CO2) on growth. Driven by local climate records, the model was applied to simulate ring widths during the period 1958–2006 for multiple trees of Pinus koraiensis from the Changbai Mountains in northeastern China. Each tree was initialised at its actual diameter at the time when local climate records started. The model produces realistic simulations of the interannual variability in ring width for different age cohorts (young, mature, and old). Both the simulations and observations show a significant positive response of tree-ring width to growing-season total photosynthetically active radiation (PAR0) and the ratio of actual to potential evapotranspiration (α), and a significant negative response to mean annual temperature (MAT). The slopes of the simulated and observed relationships with PAR0 and α are similar; the negative response to MAT is underestimated by the model. Comparison of simulations with fixed and changing atmospheric CO2 concentration shows that CO2 fertilisation over the past 50 years is too small to be distinguished in the ring-width data, given ontogenetic trends and interannual variability in climate.

Validation of Canopy Height Profile methodology for small-footprint full-waveform airborne LiDAR data in a discontinuous canopy environment

A Canopy Height Profile (CHP) procedure presented in Harding et al. (2001) for large footprint LiDAR data was tested in a closed canopy environment as a way of extracting vertical foliage profiles from LiDAR raw-waveform. In this study, an adaptation of this method to small-footprint data has been shown, tested and validated in an Australian sparse canopy forest at plot- and site-level. Further, the methodology itself has been enhanced by implementing a dataset-adjusted reflectance ratio calculation according to Armston et al. (2013) in the processing chain, and tested against a fixed ratio of 0.5 estimated for the laser wavelength of 1550nm. As a by-product of the methodology, effective leaf area index (LAIe) estimates were derived and compared to hemispherical photography-derived values. To assess the influence of LiDAR aggregation area size on the estimates in a sparse canopy environment, LiDAR CHPs and LAIes were generated by aggregating waveforms to plot- and site-level footprints (plot/site-aggregated) as well as in 5m grids (grid-processed). LiDAR profiles were then compared to leaf biomass field profiles generated based on field tree measurements. The correlation between field and LiDAR profiles was very high, with a mean R2 of 0.75 at plot-level and 0.86 at site-level for 55 plots and the corresponding 11 sites. Gridding had almost no impact on the correlation between LiDAR and field profiles (only marginally improvement), nor did the dataset-adjusted reflectance ratio. However, gridding and the dataset-adjusted reflectance ratio were found to improve the correlation between raw-waveform LiDAR and hemispherical photography LAIe estimates, yielding the highest correlations of 0.61 at plot-level and of 0.83 at site-level. This proved the validity of the approach and superiority of dataset-adjusted reflectance ratio of Armston et al. (2013) over a fixed ratio of 0.5 for LAIe estimation, as well as showed the adequacy of small-footprint LiDAR data for LAIe estimation in discontinuous canopy forests.

A model analysis of climate and CO2 controls on tree growth in a semi-arid woodland

We used a light-use efficiency model of photosynthesis coupled with a dynamic carbon allocation and tree-growth model to simulate annual growth of the gymnosperm Callitris columellaris in the semi-arid Great Western Woodlands, Western Australia, over the past 100 years. Parameter values were derived from independent observations except for sapwood specific respiration rate, fine-root turnover time, fine-root specific respiration rate and the ratio of fine-root mass to foliage area, which were estimated by Bayesian optimization. The model reproduced the general pattern of interannual variability in radial growth (tree-ring width), including the response to the shift in precipitation regimes that occurred in the 1960s. Simulated and observed responses to climate were consistent. Both showed a significant positive response of tree-ring width to total photosynthetically active radiation received and to the ratio of modeled actual to equilibrium evapotranspiration, and a significant negative response to vapour pressure deficit. However, the simulations showed an enhancement of radial growth in response to increasing atmospheric CO2 concentration (ppm) ([CO2]) during recent decades that is not present in the observations. The discrepancy disappeared when the model was recalibrated on successive 30-year windows. Then the ratio of fine-root mass to foliage area increases by 14% (from 0.127 to 0.144 kg C m-2) as [CO2] increased while the other three estimated parameters remained constant. The absence of a signal of increasing [CO2] has been noted in many tree-ring records, despite the enhancement of photosynthetic rates and water-use efficiency resulting from increasing [CO2]. Our simulations suggest that this behaviour could be explained as a consequence of a shift towards below-ground carbon allocation.

Metabolic theory predicts whole-ecosystem properties

Understanding the effects of individual organisms on material cycles and energy fluxes within ecosystems is central to predicting the impacts of human-caused changes on climate, land use, and biodiversity. Here we present a theory that integrates metabolic (organism-based bottom-up) and systems (ecosystem-based top-down) approaches to characterize how the metabolism of individuals affects the flows and stores of materials and energy in ecosystems. The theory predicts how the average residence time of carbon molecules, total system throughflow (TST), and amount of recycling vary with the body size and temperature of the organisms and with trophic organization. We evaluate the theory by comparing theoretical predictions with outputs of numerical models designed to simulate diverse ecosystem types and with empirical data for real ecosystems. Although residence times within different ecosystems vary by orders of magnitude—from weeks in warm pelagic oceans with minute phytoplankton producers to centuries in cold forests with large tree producers—as predicted, all ecosystems fall along a single line: residence time increases linearly with slope = 1.0 with the ratio of whole-ecosystem biomass to primary productivity (B/P). TST was affected predominantly by primary productivity and recycling by the transfer of energy from microbial decomposers to animal consumers. The theory provides a robust basis for estimating the flux and storage of energy, carbon, and other materials in terrestrial, marine, and freshwater ecosystems and for quantifying the roles of different kinds of organisms and environments at scales from local ecosystems to the biosphere.

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.