Fingerprints of changes in the terrestrial carbon cycle in response to large reorganizations in ocean circulation

CO2 and carbon cycle changes in the land, ocean and atmosphere are investigated using the comprehensive carbon cycle-climate model NCAR CSM1.4-carbon. Ensemble simulations are forced with freshwater perturbations applied at the North Atlantic and Southern Ocean deep water formation sites under pre-industrial climate conditions. As a result, the Atlantic Meridional Overturning Circulation reduces in each experiment to varying degrees. The physical climate fields show changes qualitatively in agreement with results documented in the literature, but there is a clear distinction between northern and southern perturbations. Changes in the physical variables, in turn, affect the land and ocean biogeochemical cycles and cause a reduction, or an increase, in the atmospheric CO2 concentration by up to 20 ppmv, depending on the location of the perturbation. In the case of a North Atlantic perturbation, the land biosphere reacts with a strong reduction in carbon stocks in some tropical locations and in high northern latitudes. In contrast, land carbon stocks tend to increase in response to a southern perturbation. The ocean is generally a sink of carbon although large reorganizations occur throughout various basins. The response of the land biosphere is strongest in the tropical regions due to a shift of the Intertropical Convergence Zone. The carbon fingerprints of this shift, either to the south or to the north depending on where the freshwater is applied, can be found most clearly in South America. For this reason, a compilation of various paleoclimate proxy records of Younger Dryas precipitation changes are compared with our model results. The proxy records, in general, show good agreement with the model's response to a North Atlantic freshwater perturbation.

Diffusion-driven transport in clayrock formations

Clay mineral-rich sedimentary formations are currently under investigation to evaluate their potential use as host formations for installation of deep underground disposal facilities for radioactive waste (e.g. Boom Clay (BE), Opalinus Clay (CH), Callovo-Oxfordian argillite (FR)). The ultimate safety of the corresponding repository concepts depends largely on the capacity of the host formation to limit the flux towards the biosphere of radionuclides (RN) contained in the waste to acceptably low levels. Data for diffusion-driven transfer in these formations shows extreme differences in the measured or modelled behaviour for various radionuclides, e. g. between halogen RN (Cl-36, I-129) and actinides (U-238,U-235, Np-237, Th-232, etc.), which result from major differences between RN of the effects on transport of two phenomena: diffusion and sorption. This paper describes recent research aimed at improving understanding of these two phenomena, focusing on the results of studies carried out during the EC Funmig IP on clayrocks from the above three formations and from the Boda formation (HU). Project results regarding phenomena governing water, cation and anion distribution and mobility in the pore volumes influenced by the negatively-charged surfaces of clay minerals show a convergence of the modelling results for behaviour at the molecular scale and descriptions based on electrical double layer models. Transport models exist which couple ion distribution relative to the clay-solution interface and differentiated diffusive characteristics. These codes are able to reproduce the main trends in behaviour observed experimentally, e.g. D-e(anion) < D-e(HTO) < D-e(cation) and D-e(anion) variations as a function of ionic strength and material density. These trends are also well-explained by models of transport through ideal porous matrices made up of a charged surface material. Experimental validation of these models is good as regards monovalent alkaline cations, in progress for divalent electrostatically-interacting cations (e.g. Sr2+) and still relatively poor for 'strongly sorbing', high K-d cations. Funmig results have clarified understanding of how clayrock mineral composition, and the corresponding organisation of mineral grain assemblages and their associated porosity, can affect mobile solute (anions, HTO) diffusion at different scales (mm to geological formation). In particular, advances made in the capacity to map clayrock mineral grain-porosity organisation at high resolution provide additional elements for understanding diffusion anisotropy and for relating diffusion characteristics measured at different scales. On the other hand, the results of studies focusing on evaluating the potential effects of heterogeneity on mobile species diffusion at the formation scale tend to show that there is a minimal effect when compared to a homogeneous property model. Finally, the results of a natural tracer-based study carried out on the Opalinus Clay formation increase confidence in the use of diffusion parameters measured on laboratory scale samples for predicting diffusion over geological time-space scales. Much effort was placed on improving understanding of coupled sorption-diffusion phenomena for sorbing cations in clayrocks. Results regarding sorption equilibrium in dispersed and compacted materials for weakly to moderately sorbing cations (Sr2+, Cs+, Co2+) tend to show that the same sorption model probably holds in both systems. It was not possible to demonstrate this for highly sorbing elements such as Eu(III) because of the extremely long times needed to reach equilibrium conditions, but there does not seem to be any clear reason why such elements should not have similar behaviour. Diffusion experiments carried out with Sr2+, Cs+ and Eu(III) on all of the clayrocks gave mixed results and tend to show that coupled diffusion-sorption migration is much more complex than expected, leading generally to greater mobility than that predicted by coupling a batch-determined K-d and Ficks law based on the diffusion behaviour of HTO. If the K-d measured on equivalent dispersed systems holds as was shown to be the case for Sr, Cs (and probably Co) for Opalinus Clay, these results indicate that these cations have a D-e value higher than HTO (up to a factor of 10 for Cs+). Results are as yet very limited for very moderate to strongly sorbing species (e.g. Co(II), Eu(III), Cu(II)) because of their very slow transfer characteristics. (C) 2011 Elsevier Ltd. All rights reserved.

DENITRIFICATION IN SOILS: FROM GENES TO ENVIRONMENTAL OUTCOMES

Denitrification is an important process of global nitrogen cycle as it removes reactive nitrogen from the biosphere, and acts as the primary source of nitrous oxide (N2O). This thesis seeks to gain better understanding of the biogeochemistry of denitrification by investigating the process from four different aspects: genetic basis, enzymatic kinetics, environmental interactions, and environmental consequences. Laboratory and field experiments were combined with modeling efforts to unravel the complexity of denitrification process under microbiological and environmental controls. Dynamics of denitrification products observed in laboratory experiments revealed an important role of constitutive denitrification enzymes, whose presence were further confirmed with quantitative analysis of functional genes encoding nitrite reductase and nitrous oxide reductase. A metabolic model of denitrification developed with explicit denitrification enzyme kinetics and representation of constitutive enzymes successfully reproduced the dynamics of N2O and N2 accumulation observed in the incubation experiments, revealing important regulatory effect of denitrification enzyme kinetics on the accumulation of denitrification products. Field studies demonstrated complex interaction of belowground N2O production, consumption and transport, resulting in two pulse pattern in the surface flux. Coupled soil gas diffusion/denitrification model showed great potential in simulating the dynamics of N2O below ground, with explicit representation of the activity of constitutive denitrification enzymes. A complete survey of environmental variables showed distinct regulation regimes on the denitrification activity from constitutive enzymes and new synthesized enzymes. Uncertainties in N2O estimation with current biogeochemical models may be reduced as accurate simulation of the dynamics of N2O in soil and surface fluxes is possible with a coupled diffusion/denitrification model that includes explicit representation of denitrification enzyme kinetics. In conclusion, denitrification is a complex ecological function regulated at cellular level. To assess the environmental consequences of denitrification and develop useful tools to mitigate N2O emissions require a comprehensive understanding of the regulatory network of denitrification with respect to microbial physiology and environmental interactions.

Net CO 2 surface emissions at Bern, Switzerland inferred from ambient observations of CO 2 , δ(O 2 /N 2 ), and 222 Rn using a customized radon tracer inversion

The 222Radon tracer method is a powerful tool to estimate local and regional surface emissions of, e.g., greenhouse gases. In this paper we demonstrate that in practice, the method as it is commonly used, produces inaccurate results in case of nonhomogeneously spread emission sources, and we propose a different approach to account for this. We have applied the new methodology to ambient observations of CO2 and 222Radon to estimate CO2 surface emissions for the city of Bern, Switzerland. Furthermore, by utilizing combined measurements of CO2 and δ(O2/N2) we obtain valuable information about the spatial and temporal variability of the main emission sources. Mean net CO2 emissions based on 2 years of observations are estimated at (11.2 ± 2.9) kt km−2 a−1. Oxidative ratios indicate a significant influence from the regional biosphere in summer/spring and fossil fuel combustion processes in winter/autumn. Our data indicate that the emissions from fossil fuels are, to a large degree, related to the combustion of natural gas which is used for heating purposes.

Current Regulating and Supporting Services: Nutrient Cycles

The study forest regulates nutrient cycles as a supporting ecosystem service mainly via retention in the biosphere and the soil organic layer. How tight the nutrient cycles are depends on environmental conditions. In this chapter, we focus on the roles of (1) deposition from the atmosphere, (2) soil moisture regime, and (3) conversion to pasture in the nutrient cycle. Between 1998 and 2010, there were a seasonal deposition of salpetric acid, an episodic deposition of Ca and Mg from Sahara dusts, and a continuous increase in reactive N inputs related to Amazonian forest fires, the El Niño Southern Oscillation cycle, and the economic development, respectively. Simultaneously, soils became increasingly drier enhancing nutrient release by mineralization. An increasing number of rain storms could considerably increase the export of N and base metals (K, Ca, Mg) via fast surface-near lateral transport in soil. Land-use change from forest to pasture introduces alkaline ashes and grass-derived organic matter. The resulting increases in soil pH and nutrient and substrate supply increase nutrient cycling rates because of enhanced microbial activity.

Transient simulations of the carbon and nitrogen dynamics in northern peatlands: from the Last Glacial Maximum to the 21st century

The development of northern high-latitude peatlands played an important role in the carbon (C) balance of the land biosphere since the Last Glacial Maximum (LGM). At present, carbon storage in northern peatlands is substantial and estimated to be 500 ± 100 Pg C (1 Pg C = 1015 g C). Here, we develop and apply a peatland module embedded in a dynamic global vegetation and land surface process model (LPX-Bern 1.0). The peatland module features a dynamic nitrogen cycle, a dynamic C transfer between peatland acrotelm (upper oxic layer) and catotelm (deep anoxic layer), hydrology- and temperature-dependent respiration rates, and peatland specific plant functional types. Nitrogen limitation down-regulates average modern net primary productivity over peatlands by about half. Decadal acrotelm-to-catotelm C fluxes vary between −20 and +50 g C m−2 yr−1 over the Holocene. Key model parameters are calibrated with reconstructed peat accumulation rates from peat-core data. The model reproduces the major features of the peat core data and of the observation-based modern circumpolar soil carbon distribution. Results from a set of simulations for possible evolutions of northern peat development and areal extent show that soil C stocks in modern peatlands increased by 365–550 Pg C since the LGM, of which 175–272 Pg C accumulated between 11 and 5 kyr BP. Furthermore, our simulations suggest a persistent C sequestration rate of 35–50 Pg C per 1000 yr in present-day peatlands under current climate conditions, and that this C sink could either sustain or turn towards a source by 2100 AD depending on climate trajectories as projected for different representative greenhouse gas concentration pathways.

Central European hardwood trees in a high-CO 2 future: synthesis of an 8-year forest canopy CO 2 enrichment project

Rapidly increasing atmospheric CO2 is not only changing the climate system but may also affect the biosphere directly through stimulation of plant growth and ecosystem carbon and nutrient cycling. Although forest ecosystems play a critical role in the global carbon cycle, experimental information on forest responses to rising CO2 is scarce, due to the sheer size of trees. Here, we present a synthesis of the only study world-wide where a diverse set of mature broadleaved trees growing in a natural forest has been exposed to future atmospheric CO2 levels (c. 550ppm) by free-air CO2 enrichment (FACE). We show that litter production, leaf traits and radial growth across the studied hardwood species remained unaffected by elevated CO2 over 8years. CO2 enrichment reduced tree water consumption resulting in detectable soil moisture savings. Soil air CO2 and dissolved inorganic carbon both increased suggesting enhanced below-ground activity. Carbon release to the rhizosphere and/or higher soil moisture primed nitrification and nitrate leaching under elevated CO2; however, the export of dissolved organic carbon remained unaltered.Synthesis. Our findings provide no evidence for carbon-limitation in five central European hardwood trees at current ambient CO2 concentrations. The results of this long-term study challenge the idea of a universal CO2 fertilization effect on forests, as commonly assumed in climate-carbon cycle models.

Climate extremes and the carbon cycle

The terrestrial biosphere is a key component of the global carbon cycle and its carbon balance is strongly influenced by climate. Continuing environmental changes are thought to increase global terrestrial carbon uptake. But evidence is mounting that climate extremes such as droughts or storms can lead to a decrease in regional ecosystem carbon stocks and therefore have the potential to negate an expected increase in terrestrial carbon uptake. Here we explore the mechanisms and impacts of climate extremes on the terrestrial carbon cycle, and propose a pathway to improve our understanding of present and future impacts of climate extremes on the terrestrial carbon budget.

Antarctic Climate Change and the Environment

The Antarctic climate system varies on timescales from orbital, through millennial to sub-annual, and is closely coupled to other parts of the global climate system. We review these variations from the perspective of the geological and glaciological records and the recent historical period from which we have instrumental data (similar to the last 50 years). We consider their consequences for the biosphere, and show how the latest numerical models project changes into the future, taking into account human actions in the form of the release of greenhouse gases and chlorofluorocarbons into the atmosphere. In doing so, we provide an essential Southern Hemisphere companion to the Arctic Climate Impact Assessment.

Water–montmorillonite systems: Neutron scattering and tracer throughdiffusion studies

Designs for deep geological respositories of nuclear waste include bentonite as a hydraulic and chemisorption buffer material to protect the biosphere from leakage of radionuclides. Bentonite is chosen because it is a cheap, naturally occurring material with the required properties. It consists essentially of montmorillonite, a swelling clay mineral. Upon contact with groundwater such clays can seal the repository by incorporating water in the interlayers of their crystalline structure. The intercalated water exhibits significantly different properties to bulk water in the surrounding interparticle pores, such as lower diffusion coefficients (González Sánchez et. al. 2008). This doctoral thesis presents water distribution and diffusion behavior on various time and space scales in montmorillonite. Experimental results are presented for Na- and Cs-montmorillonite samples with a range of bulk dry densities (0.8 to 1.7 g/cm3). The experimental methods employed were neutron scattering (backscattering, diffraction, time-of-flight), adsorption measurements (water, nitrogen) and tracer-through diffusion. For the tracer experiments the samples were fully saturated via the liquid phase under volume-constrained conditions. In contrast, for the neutron scattering experiments, the samples were hydrated via the vapor phase and subsequently compacted, leaving a significant fraction of interparticle pores unfilled with water. Owing to these differences in saturation, the water contents of the samples for neutron scattering were characterized by gravimetry whereas those for the tracer experiments were obtained from the bulk dry density. The amount of surface water in interlayer pores could be successfully discriminated from the amount of bulk-like water in interparticle pores in Na- and Csmontmorillonite using neutron spectroscopy. For the first time in the literature, the distribution of water between these two pore environments was deciphered as a function of gravimetric water content. The amount was compared to a geometrical estimation of the amount of interlayer and interparticle water determined by neutron diffraction and adsorption measurements. The relative abundances of the 1 to 4 molecular water layers in the interlayer were determined from the area ratios of the (001)-diffraction peaks. Depending on the characterization method, different fractions of surface water and interlayer water were obtained. Only surface and interlayer water exists in amontmorillonite with water contents up to 0.18 g/g according to spectroscopic measurements and up to 0.32 g/g according to geometrical estimations, respectively. At higher water contents, bulk-like and interparticle water also exists. The amounts increase monotonically, but not linearly, from zero to 0.33 g/g for bulk-like water and to 0.43 g/g for interparticle water. It was found that water most likely redistributes between the surface and interlayer sites during the spectroscopic measurements and therefore the reported fraction is relevant only below about -10 ºC (Anderson, 1967). The redistribution effect can explain the discrepancy in fractions between the methods. In a novel approach the fractions of water in different pore environments were treated as a fixed parameter to derive local diffusion coefficients for water from quasielastic neutron scattering data, in particular for samples with high water contents. Local diffusion coefficients were obtained for the 1 to 4 molecular water layers in the interlayer of 0.5·10–9, 0.9·10–9, 1.5·10–9 and 1.4·10–9 m²/s, respectively, taking account of the different water fractions (molecular water layer, bulk-like water). The diffusive transport of 22Na and HTO through Na-montmorillonite was measured on the laboratory experimental scale (i.e. cm, days) by tracer through-diffusion experiments. We confirmed that diffusion of HTO is independent of the ionic strength of the external solution in contact with the clay sample but dependent on the bulk dry density. In contrast, the diffusion of 22Na was found to depend on both the ionic strength of the pore solution and on the bulk dry density. The ratio of the pore and surface diffusion could be experimentally determined for 22Na from the dependence of the diffusion coefficient on the ionic strength. Activation energies were derived from the temperaturedependent diffusion coefficients via the Arrhenius relation. In samples with high bulk dry density the activation energies are slightly higher than those of bulk water whereas in low density samples they are lower. The activation energies as a function of ionic strengths of the pore solutions are similar for 22Na and HTO. The facts that (i) the slope of the logarithmic effective diffusion coefficients as a function of the logarithmic ionic strength is less than unity for low bulk dry densities and (ii) two water populations can be observed for high gravimetric water contents (low bulk dry densities) support the interlayer and interparticle porosity model proposed by Glaus et al. (2007), Bourg et al. (2006, 2007) and Gimmi and Kosakowski (2011).

The global carbon budget 1959–2011

Accurate assessments of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere is important to better understand the global carbon cycle, support the climate policy process, and project future climate change. Present-day analysis requires the combination of a range of data, algorithms, statistics and model estimates and their interpretation by a broad scientific community. Here we describe datasets and a methodology developed by the global carbon cycle science community to quantify all major components of the global carbon budget, including their uncertainties. We discuss changes compared to previous estimates, consistency within and among components, and methodology and data limitations. CO2 emissions from fossil fuel combustion and cement production (EFF) are based on energy statistics, while emissions from Land-Use Change (ELUC), including deforestation, are based on combined evidence from land cover change data, fire activity in regions undergoing deforestation, and models. The global atmospheric CO2 concentration is measured directly and its rate of growth (GATM) is computed from the concentration. The mean ocean CO2 sink (SOCEAN) is based on observations from the 1990s, while the annual anomalies and trends are estimated with ocean models. Finally, the global residual terrestrial CO2 sink (SLAND) is estimated by the difference of the other terms. For the last decade available (2002–2011), EFF was 8.3 ± 0.4 PgC yr−1, ELUC 1.0 ± 0.5 PgC yr−1, GATM 4.3 ± 0.1 PgC yr−1, SOCEAN 2.5 ± 0.5 PgC yr−1, and SLAND 2.6 ± 0.8 PgC yr−1. For year 2011 alone, EFF was 9.5 ± 0.5 PgC yr−1, 3.0 percent above 2010, reflecting a continued trend in these emissions; ELUC was 0.9 ± 0.5 PgC yr−1, approximately constant throughout the decade; GATM was 3.6 ± 0.2 PgC yr−1, SOCEAN was 2.7 ± 0.5 PgC yr−1, and SLAND was 4.1 ± 0.9 PgC yr−1. GATM was low in 2011 compared to the 2002–2011 average because of a high uptake by the land probably in response to natural climate variability associated to La Niña conditions in the Pacific Ocean. The global atmospheric CO2 concentration reached 391.31 ± 0.13 ppm at the end of year 2011. We estimate that EFF will have increased by 2.6% (1.9–3.5%) in 2012 based on projections of gross world product and recent changes in the carbon intensity of the economy. All uncertainties are reported as ±1 sigma (68% confidence assuming Gaussian error distributions that the real value lies within the given interval), reflecting the current capacity to characterise the annual estimates of each component of the global carbon budget. This paper is intended to provide a baseline to keep track of annual carbon budgets in the future.

The global carbon budget 1959–2011

Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere is important to better understand the global carbon cycle, support the climate policy process, and project future climate change. Present-day analysis requires the combination of a range of data, algorithms, statistics and model estimates and their interpretation by a broad scientific community. Here we describe datasets and a methodology developed by the global carbon cycle science community to quantify all major components of the global carbon budget, including their uncertainties. We discuss changes compared to previous estimates, consistency within and among components, and methodology and data limitations. Based on energy statistics, we estimate that the global emissions of CO2 from fossil fuel combustion and cement production were 9.5 ± 0.5 PgC yr−1 in 2011, 3.0 percent above 2010 levels. We project these emissions will increase by 2.6% (1.9–3.5%) in 2012 based on projections of Gross World Product and recent changes in the carbon intensity of the economy. Global net CO2 emissions from Land-Use Change, including deforestation, are more difficult to update annually because of data availability, but combined evidence from land cover change data, fire activity in regions undergoing deforestation and models suggests those net emissions were 0.9 ± 0.5 PgC yr−1 in 2011. The global atmospheric CO2 concentration is measured directly and reached 391.38 ± 0.13 ppm at the end of year 2011, increasing 1.70 ± 0.09 ppm yr−1 or 3.6 ± 0.2 PgC yr−1 in 2011. Estimates from four ocean models suggest that the ocean CO2 sink was 2.6 ± 0.5 PgC yr−1 in 2011, implying a global residual terrestrial CO2 sink of 4.1 ± 0.9 PgC yr−1. All uncertainties are reported as ±1 sigma (68% confidence assuming Gaussian error distributions that the real value lies within the given interval), reflecting the current capacity to characterise the annual estimates of each component of the global carbon budget. This paper is intended to provide a baseline to keep track of annual carbon budgets in the future.

Global carbon budget 2013

Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere is important to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe data sets and a methodology to quantify all major components of the global carbon budget, including their uncertainties, based on the combination of a range of data, algorithms, statistics and model estimates and their interpretation by a broad scientific community. We discuss changes compared to previous estimates, consistency within and among components, alongside methodology and data limitations. CO2 emissions from fossil-fuel combustion and cement production (EFF) are based on energy statistics, while emissions from land-use change (ELUC), mainly deforestation, are based on combined evidence from land-cover change data, fire activity associated with deforestation, and models. The global atmospheric CO2 concentration is measured directly and its rate of growth (GATM) is computed from the annual changes in concentration. The mean ocean CO2 sink (SOCEAN) is based on observations from the 1990s, while the annual anomalies and trends are estimated with ocean models. The variability in SOCEAN is evaluated for the first time in this budget with data products based on surveys of ocean CO2 measurements. The global residual terrestrial CO2 sink (SLAND) is estimated by the difference of the other terms of the global carbon budget and compared to results of independent dynamic global vegetation models forced by observed climate, CO2 and land cover change (some including nitrogen–carbon interactions). All uncertainties are reported as ± 1 σ, reflecting the current capacity to characterise the annual estimates of each component of the global carbon budget. For the last decade available (2003–2012), EFF was 8.6 ± 0.4 GtC yr − 1, ELUC 0.9 ± 0.5 GtC yr − 1, GATM 4.3 ± 0.1 GtC yr − 1, S OCEAN 2.5 ± 0.5 GtC yr − 1, and S LAND 2.8 ± 0.8 GtC yr − 1. For year 2012 alone, EFF grew to 9.7 ± 0.5 GtC yr − 1, 2.2 % above 2011, reflecting a continued growing trend in these emissions, GATM was 5.1 ± 0.2 GtC yr − 1, SOCEANwas 2.9 ± 0.5 GtC yr −1, and assuming an ELU Cof 1.0 ± 0.5 GtC yr − 1 (based on the 2001–2010 average), SLAND was 2.7 ± 0.9 GtC yr − 1. GATM was high in 2012 compared to the 2003–2012 average, almost entirely reflecting the high EFF. The global atmospheric CO2 con- centration reached 392.52 ± 0.10 ppm averaged over 2012. We estimate that EFF will increase by 2.1 % (1.1–3.1 %) to 9.9 ± 0.5 GtC in 2013, 61 % above emissions in 1990, based on projections of world gross domestic product and recent changes in the carbon intensity of the economy. With this projection, cumulative emissions of CO2 will reach about 535 ± 55 GtC for 1870–2013, about 70 % from EFF (390 ± 20 GtC) and 30 % from ELUC (145 ± 50 GtC). This paper also documents any changes in the methods and data sets used in this new carbon budget from previous budgets (Le Quéré et al., 2013). All observations presented here can be downloaded from the Carbon Dioxide Information Analysis Center.

Mapping tree density at a global scale

The global extent and distribution of forest trees is central to our understanding of the terrestrial biosphere. We provide the first spatially continuous map of forest tree density at a global scale. This map reveals that the global number of trees is approximately 3.04 trillion, an order of magnitude higher than the previous estimate. Of these trees, approximately 1.39 trillion exist in tropical and subtropical forests, with 0.74 trillion in boreal regions and 0.61 trillion in temperate regions. Biome-level trends in tree density demonstrate the importance of climate and topography in controlling local tree densities at finer scales, as well as the overwhelming effect of humans across most of the world. Based on our projected tree densities, we estimate that over 15 billion trees are cut down each year, and the global number of trees has fallen by approximately 46% since the start of human civilization.