| Literature DB >> 30305744 |
Sebastiaan Luyssaert1,2, Guillaume Marie3, Aude Valade4,5, Yi-Ying Chen6,7, Sylvestre Njakou Djomo8, James Ryder6,9, Juliane Otto6,10, Kim Naudts6,11, Anne Sofie Lansø6, Josefine Ghattas4, Matthew J McGrath6.
Abstract
The Paris Agreement promotes forest management as a pathway towards halting climate warming through the reduction of carbon dioxide (CO2) emissions1. However, the climate benefits from carbon sequestration through forest management may be reinforced, counteracted or even offset by concurrent management-induced changes in surface albedo, land-surface roughness, emissions of biogenic volatile organic compounds, transpiration and sensible heat flux2-4. Consequently, forest management could offset CO2 emissions without halting global temperature rise. It therefore remains to be confirmed whether commonly proposed sustainable European forest-management portfolios would comply with the Paris Agreement-that is, whether they can reduce the growth rate of atmospheric CO2, reduce the radiative imbalance at the top of the atmosphere, and neither increase the near-surface air temperature nor decrease precipitation by the end of the twenty-first century. Here we show that the portfolio made up of management systems that locally maximize the carbon sink through carbon sequestration, wood use and product and energy substitution reduces the growth rate of atmospheric CO2, but does not meet any of the other criteria. The portfolios that maximize the carbon sink or forest albedo pass only one-different in each case-criterion. Managing the European forests with the objective of reducing near-surface air temperature, on the other hand, will also reduce the atmospheric CO2 growth rate, thus meeting two of the four criteria. Trade-off are thus unavoidable when using European forests to meet climate objectives. Furthermore, our results demonstrate that if present-day forest cover is sustained, the additional climate benefits achieved through forest management would be modest and local, rather than global. On the basis of these findings, we argue that Europe should not rely on forest management to mitigate climate change. The modest climate effects from changes in forest management imply, however, that if adaptation to future climate were to require large-scale changes in species composition and silvicultural systems over Europe5,6, the forests could be adapted to climate change with neither positive nor negative climate effects.Entities:
Mesh:
Substances:
Year: 2018 PMID: 30305744 PMCID: PMC6277009 DOI: 10.1038/s41586-018-0577-1
Source DB: PubMed Journal: Nature ISSN: 0028-0836 Impact factor: 49.962
Figure 1Surface areas (x 10,000 km2) under forest management by the year 2100 for portfolios that target maximising the carbon sink, continue present-day management, and reduce the near-surface air temperature. Forest management distinguishes between tree species composition and silvicultural systems. The inset presents the mean values for all of Europe. Regional difference are shown for three geographical regions, each shown in a different shade of grey.
Changes in surface area (km2) by 2100 for six different forest management portfolios over Europe. Note that the total surface area of forests was held constant at 2,000,000 km2 between 2010 and 2100, for reasons described in “Simulation experiment”.
| Change in surface area (km2) | Business as usual (BAU) | Maximise carbon sink | Maximise albedo | Minimise carbon sink | Minimise albedo | Reduce near-surface temperature |
|---|---|---|---|---|---|---|
| Deciduous to conifers | 0 | 475,000 | 30,000 | 6,000 | 516,000 | 41,000 |
| Conifers to deciduous | 0 | 266,000 | 590,000 | 236,000 | 26,000 | 534,000 |
| Net increase conifers | 0 | 209,000 | -560,000 | -230,000 | 490,000 | -493,000 |
| Net increase thin and fell | 0 | -280,000 | -330,000 | -390,000 | -230,000 | -680,000 |
| Net increase coppice | 0 | -20,000 | 130,000 | -130,000 | -210,000 | 600,000 |
| Net increase unmanaged | 0 | 300,000 | 200,000 | 520,000 | 440,000 | 80,000 |
Biogeochemical and biophysical effects in 2100 for four different forest management portfolios over Europe. The business as usual simulation which served as a control, was repeated three times with slightly different initial atmospheric conditions (see Methods “Equilibrium climate for the management portfolios”). The variability between these three repetitions was considered the minimal model noise of the climate model. The reported noise was taken to be the definition of one standard deviation. TOA denotes the radiative imbalance at the top of the atmosphere.
| Variable name (units) | Business as usual (BAU) | Maximise carbon sink | Maximise albedo | Reduce near-surface temperature |
|---|---|---|---|---|
| Global average TOA (W m-2) | 4.31 ± 0.01 | 4.31 | 4.33 | 4.32 |
| Δ2100-2010 CO2 sink & avoided emissions (Pg C) | 4.7 | 12.8 | 5.0 | 8.1 |
| Δ2100-2010 net cumulated atmospheric CO2 (Pg C) | 2.7 | 7.0 | 2.8 | 4.5 |
| Atmospheric CO2 (ppm) | 934.6 | 932.6 | 934.6 | 933.8 |
| Near surface temperature (K) | 283.84 ± <0.001 | 283.84 | 283.83 | 283.81 |
| Annual precipitation (mm) | 734.7 ± 0.1 | 732.6 | 730.0 | 730.9 |
| Summer precipitation (mm) | 166.1 ± 0.1 | 165.2 | 163.7 | 165.0 |
| Wood harvest (Tg C y-1) | 203.2 | 179.5 | 144.5 | 151.6 |
| Surface albedo (-) | 0.113 ± <0.0001 | 0.113 | 0.128 | 0.126 |
| Evapotranspiration (mm) | 555.5 ± 0.1 | 552.8 | 546.4 | 549.2 |
| Latent heat (W m-2) | 44.35 ± <0.01 | 44.13 | 43.60 | 43.82 |
| Sensible heat (W m-2) | 26.67 ± <0.01 | 26.82 | 27.28 | 27.00 |
| Total cloud cover (%) | 46.8 ± <0. 1 | 46.7 | 46.7 | 46.6 |
Figure 2Changes in, and main drivers of, near-surface air temperature (ΔT; K) in February and March by the turn of the 21st-century for a forest management portfolio that reduces the near-surface air temperature. (a) Spatially explicit changes in near-surface air temperature (K) in February and March. (b) Drivers of the changes in springtime near-surface air temperature for 0.5 degree latitudinal bands. In subplot (a) temperature changes less than 1.96 times the standard deviations are shown in white. Where, the standard deviation represents the minimal noise of LMDzORCAN (see Methods “Equilibrium climate for the management portfolios”). The change in near-surface air temperature (T) due to changes in atmospheric emissivity (ε) is written as ΔT. By analogy ΔT is the change in air temperature due to change in the ground heat flux, ΔT due to changes in turbulent fluxes, ΔT due to changes in shortwave incoming radiation which in this simulation experiment is a proxy for cloud cover, ΔT due to changes in surface albedo, and ΔT due to changes in atmospheric circulation.
Extended Data Figure 1Drivers of the changes in mean bimonthly near-surface air temperature (ΔT; K) for 0.5 degree latitudinal bands. The change in near-surface air temperature (T) due to changes in atmospheric emissivity ε is written as ΔT|ε. By analogy ΔT|G is the change in air temperature due to change in the ground heat flux, ΔT|LE+H due to changes in turbulent fluxes, ΔT|R due to changes in shortwave incoming radiation (which in this simulation experiment is a proxy for cloud cover), ΔT|α due to changes in surface albedo, and ΔT|circ due to changes in atmospheric circulation. Although all the components contribute to the near-surface air temperature, changes in emissivity always result in a cooling and changes in shortwave incoming radiation always result in warming. Consequently, emissivity and incoming shortwave radiation cannot explain the seasonal variation in the changes in near-surface air temperature. The other components are in some months positively correlated with near-surface air temperature whereas they are negatively correlated for other months, excluding them from being the main driver of changes in near-surface air temperature. Suggesting the net effect is the outcome of the interplay between the different components.
Extended Data Figure 2Relationship between changes in springtime near-surface air temperature (K) and changes in fractional cover of deciduous forest (km2) for 0.5 degree latitudinal bands over Europe. Locations where the tree species were maintained between 2010 and 2100 (i.e. Δ deciduous area on the X-axis equals zero) could experience similar air temperature changes as neighbouring locations where one species was converted into another, especially in Scandinavia, suggesting advection of heat and moisture. Nevertheless, at lower latitudes the spatial scale of this advection was limited to a few pixels (e.g., Fig. 2A) corresponding to a range of 50 to 200 km. Furthermore, the temperature effect quickly saturated with the fractional cover change and showed a strong dependence on geographical location91. Whether the apparent geographical dependency was the outcome of climatic differences and/or differences between northern and southern European deciduous species could not be established by the experimental setup used in this study.
Extended Data Figure 3Setup of the simulation experiments following the description in “Simulation experiment”. Simulations with ORCHIDEE-CAN are shown in black, and simulations with LMDzORCAN are shown in red. Blue boxes show intermediate calculations making use of simulation results (see Methods “Spatially optimised management portfolios” and “Equilibrium climate for the management portfolios”). The labels of the different simulations shown in this figure are the same labels as those used to run and archive the simulations. Note that in the main text the results of BBESTT2M were presented as “reduced near-surface air temperature”, BESTALBEDO as “maximise surface albedo”, BWORSTALBEDO as “minimise surface albedo, BBESTLCA as “maximise C-sink”, BWORSTLCA as “minimise C-sink”, and BWAC as “business as usual”. BWAC, BWAC-P1 and BWAC-P2 were used to calculate the minimal model noise as explained in Methods “Simulation experiment”.
Biogeochemical and biophysical effects in 2100 for two additional –compared to Table 1– forest management portfolios over Europe.
| Variable name (units) | Minimise carbon sink | Minimise albedo |
|---|---|---|
| TOA (W m-2) | 4.32 | 4.32 |
| Δ2100-2010 CO2 sink & avoided emissions (Pg C) | 0.7 | 10.5 |
| Δ2100-2010 net cumulated atmospheric CO2 (Pg C) | 0.5 | 5.7 |
| Atmospheric CO2 (ppm) | 935.7 | 933.2 |
| Near surface temperature (K) | 283.85 | 283.86 |
| Annual precipitation (mm) | 733.1 | 734.2 |
| Summer precipitation (mm) | 164.0 | 165.4 |
| Wood harvest (Tg C y-1) | 122.9 | 176.2 |
| Surface albedo (-) | 0.119 | 0.107 |
| Evapotranspiration (mm) | 550.0 | 553.9 |
| Latent heat (W m-2) | 43.90 | 44.23 |
| Sensible heat (W m-2) | 27.12 | 26.81 |
| Total cloud cover (%) | 46.8 | 46.8 |
Key characteristics of the individual model runs in the simulation experiment. The model runs are presented in the same order as Extended Data Fig. 3. For each model run the time period is described by its start year, end year and the length of the simulation in years (Years) together with the simulation used to initialize key characteristics of the biosphere (Initial state). The atmospheric CO2, CH4, N2O, CFC11, and CFC12 concentrations at the end of the simulation were reported and their values were extracted from Refs. 86 and 95. For the portfolio model runs atmospheric CO2 concentrations were adjusted for the simulated carbon sink after discounting for ocean uptake as outlined in “Atmospheric composition” and “Life cycle analysis”. Climate and other forcing agents including sea surface temperature, sea ice extent, and atmospheric aerosol concentrations were retrieved from the RCP8.5 simulation with the IPSL-CM5 model59 as part of the AR5 model inter-comparison. In this study, forest management consisted of two activities: species changes (Species) and a silvicultural treatment (Silviculture). For historical model runs a forest management reconstruction was used47, and a single year indicates the reconstruction for that specific year was used. For future simulations, species distribution and/or silvicultural management was either fixed to the 2010 distribution or was changed towards deciduous or conifers for the species and/or conservation, high-stand, or coppice for the silvicultural system (see Methods “Forest cover and forest silvicultural reconstruction”, “Future species”, and “Future silviculture”). The labels of the different simulations shown in this table are the same labels as those used to run and archive the simulations. Note that in the main text the results of BBESTT2M were presented as “reduced near-surface air temperature”, BESTALBEDO as “maximise surface albedo”, BWORSTALBEDO as “minimise surface albedo, BBESTLCA as “maximise C-sink”, BWORSTLCA as “minimise C-sink”, and BWAC as “business as usual”. BWAC, BWAC-P1 and BWAC-P2 were used to calculate the minimal model noise as explained in “Simulation experiment”.
| Simulation label | Period | Years | Initial state | Climate | CO2 (ppm) | CH4 (ppb) | N2O (ppb) | CFC11 (ppt) | CFC12 (ppt) | Other | Species | Silviculture |
|---|---|---|---|---|---|---|---|---|---|---|---|---|
| SPIN1 | 1600/1600 | 260 | n.a. | 1901/1920 | 277.9 | n.a. | n.a. | n.a. | n.a. | n.a. | 1600 | 1600 |
| SPIN2 | 1600/1600 | 40 | SPIN1 | 1901/1920 | 277.9 | n.a. | n.a. | n.a. | n.a. | n.a. | 1600 | 1600 |
| TRANS1 | 1601/1900 | 300 | SPIN2 | 1901/1930 | 295.8 | n.a. | n.a. | n.a. | n.a. | n.a. | Recon. | Recon. |
| TRANS2 | 1901/2010 | 110 | TRANS1 | 1901/2010 | 395.8 | n.a. | n.a. | n.a. | n.a. | n.a. | Recon. | Recon. |
| CWAC | 2011/2100 | 90 | TRANS2 | RCP8.5 | 935.8 | n.a. | n.a. | n.a. | n.a. | n.a. | 2010 | 2010 |
| CWA1 | 2011/2100 | 90 | TRANS2 | RCP8.5 | 935.8 | n.a. | n.a. | n.a. | n.a. | n.a. | 2010 | Conser. |
| CWA2 | 2011/2100 | 90 | TRANS2 | RCP8.5 | 935.8 | n.a. | n.a. | n.a. | n.a. | n.a. | 2010 | Thin&F. |
| CWC2 | 2011/2100 | 90 | TRANS2 | RCP8.5 | 935.8 | n.a. | n.a. | n.a. | n.a. | n.a. | Coni. | Thin&F. |
| CWD2 | 2011/2100 | 90 | TRANS2 | RCP8.5 | 935.8 | n.a. | n.a. | n.a. | n.a. | n.a. | Deci. | Thin&F. |
| CWD3 | 2011/2100 | 90 | TRANS2 | RCP8.5 | 935.8 | n.a. | n.a. | n.a. | n.a. | n.a. | Deci. | Coppice |
| BWAC | 2101/2101 | 10 | CWAC | n.a. | 934.6 | 3751 | 435 | 26 | 167 | RCP8.5 | 2010 | 2010 |
| BWA1 | 2101/2101 | 10 | CWA1 | n.a. | 934.6 | 3751 | 435 | 26 | 167 | RCP8.5 | 2010 | Conser. |
| BWA2 | 2101/2101 | 10 | CWA2 | n.a. | 934.6 | 3751 | 435 | 26 | 167 | RCP8.5 | 2010 | Thin&F. |
| BWC2 | 2101/2101 | 10 | CWC2 | n.a. | 934.6 | 3751 | 435 | 26 | 167 | RCP8.5 | Coni. | Thin&F. |
| BWD2 | 2101/2101 | 10 | CWD2 | n.a. | 934.6 | 3751 | 435 | 26 | 167 | RCP8.5 | Deci. | Thin&F. |
| BWD3 | 2101/2101 | 10 | CWD3 | n.a. | 934.6 | 3751 | 435 | 26 | 167 | RCP8.5 | Deci. | Coppice |
| CBESTT2M | 2011/2100 | 90 | Optimised | RCP8.5 | 935.8 | n.a. | n.a. | n.a. | n.a. | n.a. | Optim. | Optim. |
| CBESTLCA | 2011/2100 | 90 | Optimised | RCP8.5 | 935.8 | n.a. | n.a. | n.a. | n.a. | n.a. | Optim. | Optim. |
| CWORSTLCA | 2011/2100 | 90 | Optimised | RCP8.5 | 935.8 | n.a. | n.a. | n.a. | n.a. | n.a. | Optim. | Optim. |
| CBESTALBEDO | 2011/2100 | 90 | Optimised | RCP8.5 | 935.8 | n.a. | n.a. | n.a. | n.a. | n.a. | Optim. | Optim. |
| CWORSTALBEDO | 2011/2100 | 90 | Optimised | RCP8.5 | 935.8 | n.a. | n.a. | n.a. | n.a. | n.a. | Optim. | Optim. |
| BWAC | 2101/2101 | 20 | CWAC | n.a. | 934.6 | 3751 | 435 | 26 | 167 | RCP8.5 | 2010 | 2010 |
| BWAC-P1 | 2101/2101 | 20 | CWAC | n.a. | 934.6 | 3751 | 435 | 26 | 167 | RCP8.5 | 2010 | 2010 |
| BWAC-P2 | 2101/2101 | 20 | CWAC | n.a. | 934.6 | 3751 | 435 | 26 | 167 | RCP8.5 | 2010 | 2010 |
| BBESTT2M | 2101/2101 | 20 | CBESTT2M | n.a. | 933.8 | 3751 | 435 | 26 | 167 | RCP8.5 | 2010 | 2010 |
| BBESTLCA | 2101/2101 | 20 | CBESTLCA | n.a. | 932.6 | 3751 | 435 | 26 | 167 | RCP8.5 | 2010 | 2010 |
| BWORSTLCA | 2101/2101 | 20 | CWORSTLCA | n.a. | 935.7 | 3751 | 435 | 26 | 167 | RCP8.5 | 2010 | 2010 |
| BBESTALBEDO | 2101/2101 | 20 | CBESTALBEDO | n.a. | 934.6 | 3751 | 435 | 26 | 167 | RCP8.5 | 2010 | 2010 |
| BWORSTALBEDO | 2101/2101 | 20 | CWORSTALBEDO | n.a. | 933.2 | 3751 | 435 | 26 | 167 | RCP8.5 | 2010 | 2010 |
Emission coefficients, conversion factors, and assumptions used in the European wide life cycle analysis.
| Component | Unit | Value | Source |
|---|---|---|---|
| Carbon in biomass | g g-1 | 0.5 | Assumed |
| Transport distance roundwood | Km | 100 | Assumed |
| Transport distance fuelwood | Km | 50 | Assumed |
| Transport by truck | kg CO2 tkm-1 | 1.12 | Ref. 100 |
| Mechanized harvest | kg CO2 ha-1 | 233 | Ref. 101 |
| Mechanized planting | kg CO2 ha-1 | 93 | Ref. 101 |
| Mechanized thinning | kg CO2 ha-1 | 69 | Ref. 101 |
| Product substitution | kg CO2 kg-1 CO2 | 1.1 | Ref. 101 |
| Energy density of biomass | GJ t-1 | 19.3 | Ref. 102 |
| Conversion efficiency | % | 34 | Ref. 103 |
| Energy from biomass-based electricity | GJ, t-1 | 6.6 | Energy density of biomass multiplied with the conversion efficiency |
| Emissions from biomass-based electricity | kg CO2 kg-1 CO2 | 1.05 | Assuming that drying consumes 0.05 kg CO2 kg-1 CO2 and burning or gasifying woody biomass produces 1 kg CO2 kg-1 CO2 |
Country based CO2 emission factors (g CO2 eq kWh-1) for the current non-renewable electricity mix of energy production based on ref. 100.
| Country | Emission factor |
|---|---|
| Albania, Belarus, Kosovo, Macedonia, Moldova, | 1020 |
| Montenegro & Ukraine | |
| Andorra, France & Monaco | 810 |
| Austria & Liechtenstein | 777 |
| Belgium | 687 |
| Bosnia & Herzegovina | 1017 |
| Bulgaria | 981 |
| Croatia | 812 |
| Cyprus, Iceland & Malta | 868 |
| Czech Republic | 1010 |
| Denmark | 904 |
| Estonia | 1014 |
| Finland | 853 |
| Germany | 927 |
| Greece | 894 |
| Hungary | 780 |
| Ireland | 766 |
| Italy | 744 |
| Latvia | 615 |
| Lithuania | 591 |
| Luxembourg | 614 |
| Netherlands | 748 |
| Norway | 641 |
| Poland | 1000 |
| Portugal | 840 |
| Romania | 907 |
| Serbia | 1012 |
| Slovakia | 842 |
| Slovenia | 982 |
| Spain | 797 |
| Sweden | 857 |
| Switzerland | 628 |
| United Kingdom | 854 |