Modelling the impact of forest management and CO2-fertilisation on growth and demography in a Sitka spruce plantation.

Afforestation and reforestation to meet 'Net Zero' emissions targets are considered a necessary policy by many countries. Their potential benefits are usually assessed through forest carbon and growth models. The implementation of vegetation demography gives scope to represent forest management and other size-dependent processes within land surface models (LSMs). In this paper, we evaluate the impact of including management within an LSM that represents demography, using both in-situ and reanalysis climate drivers at a mature, upland Sitka spruce plantation in Northumberland, UK. We compare historical simulations with fixed and variable CO2 concentrations, and with and without tree thinning implemented. Simulations are evaluated against the observed vegetation structure and carbon fluxes. Including thinning and the impact of increasing CO2 concentration ('CO2 fertilisation') gave more realistic estimates of stand-structure and physical characteristics. Historical CO2 fertilisation had a noticeable effect on the Gross Primary Productivity seasonal-diurnal cycle and contributed to approximately 7% higher stand biomass by 2018. The net effect of both processes resulted in a decrease of tree density and biomass, but an increase in tree height and leaf area index.

The Montreal Protocol protects the terrestrial carbon sink.

The control of the production of ozone-depleting substances through the Montreal Protocol means that the stratospheric ozone layer is recovering1 and that consequent increases in harmful surface ultraviolet radiation are being avoided2,3. The Montreal Protocol has co-benefits for climate change mitigation, because ozone-depleting substances are potent greenhouse gases4-7. The avoided ultraviolet radiation and climate change also have co-benefits for plants and their capacity to store carbon through photosynthesis8, but this has not previously been investigated. Here, using a modelling framework that couples ozone depletion, climate change, damage to plants by ultraviolet radiation and the carbon cycle, we explore the benefits of avoided increases in ultraviolet radiation and changes in climate on the terrestrial biosphere and its capacity as a carbon sink. Considering a range of strengths for the effect of ultraviolet radiation on plant growth8-12, we estimate that there could have been 325-690 billion tonnes less carbon held in plants and soils by the end of this century (2080-2099) without the Montreal Protocol (as compared to climate projections with controls on ozone-depleting substances). This change could have resulted in an additional 115-235 parts per million of atmospheric carbon dioxide, which might have led to additional warming of global-mean surface temperature by 0.50-1.0 degrees. Our findings suggest that the Montreal Protocol may also be helping to mitigate climate change through avoided decreases in the land carbon sink.

Rainfall manipulation experiments as simulated by terrestrial biosphere models: Where do we stand?

Changes in rainfall amounts and patterns have been observed and are expected to continue in the near future with potentially significant ecological and societal consequences. Modelling vegetation responses to changes in rainfall is thus crucial to project water and carbon cycles in the future. In this study, we present the results of a new model-data intercomparison project, where we tested the ability of 10 terrestrial biosphere models to reproduce the observed sensitivity of ecosystem productivity to rainfall changes at 10 sites across the globe, in nine of which, rainfall exclusion and/or irrigation experiments had been performed. The key results are as follows: (a) Inter-model variation is generally large and model agreement varies with timescales. In severely water-limited sites, models only agree on the interannual variability of evapotranspiration and to a smaller extent on gross primary productivity. In more mesic sites, model agreement for both water and carbon fluxes is typically higher on fine (daily-monthly) timescales and reduces on longer (seasonal-annual) scales. (b) Models on average overestimate the relationship between ecosystem productivity and mean rainfall amounts across sites (in space) and have a low capacity in reproducing the temporal (interannual) sensitivity of vegetation productivity to annual rainfall at a given site, even though observation uncertainty is comparable to inter-model variability. (c) Most models reproduced the sign of the observed patterns in productivity changes in rainfall manipulation experiments but had a low capacity in reproducing the observed magnitude of productivity changes. Models better reproduced the observed productivity responses due to rainfall exclusion than addition. (d) All models attribute ecosystem productivity changes to the intensity of vegetation stress and peak leaf area, whereas the impact of the change in growing season length is negligible. The relative contribution of the peak leaf area and vegetation stress intensity was highly variable among models.

Stomatal optimization based on xylem hydraulics (SOX) improves land surface model simulation of vegetation responses to climate.

Land surface models (LSMs) typically use empirical functions to represent vegetation responses to soil drought. These functions largely neglect recent advances in plant ecophysiology that link xylem hydraulic functioning with stomatal responses to climate. We developed an analytical stomatal optimization model based on xylem hydraulics (SOX) to predict plant responses to drought. Coupling SOX to the Joint UK Land Environment Simulator (JULES) LSM, we conducted a global evaluation of SOX against leaf- and ecosystem-level observations. SOX simulates leaf stomatal conductance responses to climate for woody plants more accurately and parsimoniously than the existing JULES stomatal conductance model. An ecosystem-level evaluation at 70 eddy flux sites shows that SOX decreases the sensitivity of gross primary productivity (GPP) to soil moisture, which improves the model agreement with observations and increases the predicted annual GPP by 30% in relation to JULES. SOX decreases JULES root-mean-square error in GPP by up to 45% in evergreen tropical forests, and can simulate realistic patterns of canopy water potential and soil water dynamics at the studied sites. SOX provides a parsimonious way to incorporate recent advances in plant hydraulics and optimality theory into LSMs, and an alternative to empirical stress factors.

Surface-Atmosphere Coupling Scale, the Fate of Water, and Ecophysiological Function in a Brazilian Forest.

Tropical South America plays a central role in global climate. Bowen ratio teleconnects to circulation and precipitation processes far afield, and the global CO2 growth rate is strongly influenced by carbon cycle processes in South America. However, quantification of basin-wide seasonality of flux partitioning between latent and sensible heat, the response to anomalies around climatic norms, and understanding of the processes and mechanisms that control the carbon cycle remains elusive. Here, we investigate simulated surface-atmosphere interaction at a single site in Brazil, using models with different representations of precipitation and cloud processes, as well as differences in scale of coupling between the surface and atmosphere. We find that the model with parameterized clouds/precipitation has a tendency toward unrealistic perpetual light precipitation, while models with explicit treatment of clouds produce more intense and less frequent rain. Models that couple the surface to the atmosphere on the scale of kilometers, as opposed to tens or hundreds of kilometers, produce even more realistic distributions of rainfall. Rainfall intensity has direct consequences for the "fate of water," or the pathway that a hydrometeor follows once it interacts with the surface. We find that the model with explicit treatment of cloud processes, coupled to the surface at small scales, is the most realistic when compared to observations. These results have implications for simulations of global climate, as the use of models with explicit (as opposed to parameterized) cloud representations becomes more widespread.

Decadal biomass increment in early secondary succession woody ecosystems is increased by CO2 enrichment.

Increasing atmospheric CO2 stimulates photosynthesis which can increase net primary production (NPP), but at longer timescales may not necessarily increase plant biomass. Here we analyse the four decade-long CO2-enrichment experiments in woody ecosystems that measured total NPP and biomass. CO2 enrichment increased biomass increment by 1.05 ± 0.26 kg C m-2 over a full decade, a 29.1 ± 11.7% stimulation of biomass gain in these early-secondary-succession temperate ecosystems. This response is predictable by combining the CO2 response of NPP (0.16 ± 0.03 kg C m-2 y-1) and the CO2-independent, linear slope between biomass increment and cumulative NPP (0.55 ± 0.17). An ensemble of terrestrial ecosystem models fail to predict both terms correctly. Allocation to wood was a driver of across-site, and across-model, response variability and together with CO2-independence of biomass retention highlights the value of understanding drivers of wood allocation under ambient conditions to correctly interpret and predict CO2 responses.

Land-use emissions play a critical role in land-based mitigation for Paris climate targets.

Scenarios that limit global warming to below 2 °C by 2100 assume significant land-use change to support large-scale carbon dioxide (CO2) removal from the atmosphere by afforestation/reforestation, avoided deforestation, and Biomass Energy with Carbon Capture and Storage (BECCS). The more ambitious mitigation scenarios require even greater land area for mitigation and/or earlier adoption of CO2 removal strategies. Here we show that additional land-use change to meet a 1.5 °C climate change target could result in net losses of carbon from the land. The effectiveness of BECCS strongly depends on several assumptions related to the choice of biomass, the fate of initial above ground biomass, and the fossil-fuel emissions offset in the energy system. Depending on these factors, carbon removed from the atmosphere through BECCS could easily be offset by losses due to land-use change. If BECCS involves replacing high-carbon content ecosystems with crops, then forest-based mitigation could be more efficient for atmospheric CO2 removal than BECCS.

Large sensitivity in land carbon storage due to geographical and temporal variation in the thermal response of photosynthetic capacity.

Plant temperature responses vary geographically, reflecting thermally contrasting habitats and long-term species adaptations to their climate of origin. Plants also can acclimate to fast temporal changes in temperature regime to mitigate stress. Although plant photosynthetic responses are known to acclimate to temperature, many global models used to predict future vegetation and climate-carbon interactions do not include this process. We quantify the global and regional impacts of biogeographical variability and thermal acclimation of temperature response of photosynthetic capacity on the terrestrial carbon (C) cycle between 1860 and 2100 within a coupled climate-carbon cycle model, that emulates 22 global climate models. Results indicate that inclusion of biogeographical variation in photosynthetic temperature response is most important for present-day and future C uptake, with increasing importance of thermal acclimation under future warming. Accounting for both effects narrows the range of predictions of the simulated global land C storage in 2100 across climate projections (29% and 43% globally and in the tropics, respectively). Contrary to earlier studies, our results suggest that thermal acclimation of photosynthetic capacity makes tropical and temperate C less vulnerable to warming, but reduces the warming-induced C uptake in the boreal region under elevated CO2 .

Land use change and El Niño-Southern Oscillation drive decadal carbon balance shifts in Southeast Asia.

An integrated understanding of the biogeochemical consequences of climate extremes and land use changes is needed to constrain land-surface feedbacks to atmospheric CO2 from associated climate change. Past assessments of the global carbon balance have shown particularly high uncertainty in Southeast Asia. Here, we use a combination of model ensembles to show that intensified land use change made Southeast Asia a strong source of CO2 from the 1980s to 1990s, whereas the region was close to carbon neutral in the 2000s due to an enhanced CO2 fertilization effect and absence of moderate-to-strong El Niño events. Our findings suggest that despite ongoing deforestation, CO2 emissions were substantially decreased during the 2000s, largely owing to milder climate that restores photosynthetic capacity and suppresses peat and deforestation fire emissions. The occurrence of strong El Niño events after 2009 suggests that the region has returned to conditions of increased vulnerability of carbon stocks.

Implications of improved representations of plant respiration in a changing climate.

Land-atmosphere exchanges influence atmospheric CO2. Emphasis has been on describing photosynthetic CO2 uptake, but less on respiration losses. New global datasets describe upper canopy dark respiration (R d) and temperature dependencies. This allows characterisation of baseline R d, instantaneous temperature responses and longer-term thermal acclimation effects. Here we show the global implications of these parameterisations with a global gridded land model. This model aggregates R d to whole-plant respiration R p, driven with meteorological forcings spanning uncertainty across climate change models. For pre-industrial estimates, new baseline R d increases R p and especially in the tropics. Compared to new baseline, revised instantaneous response decreases R p for mid-latitudes, while acclimation lowers this for the tropics with increases elsewhere. Under global warming, new R d estimates amplify modelled respiration increases, although partially lowered by acclimation. Future measurements will refine how R d aggregates to whole-plant respiration. Our analysis suggests R p could be around 30% higher than existing estimates.

Challenging terrestrial biosphere models with data from the long-term multifactor Prairie Heating and CO2 Enrichment experiment.

Multifactor experiments are often advocated as important for advancing terrestrial biosphere models (TBMs), yet to date, such models have only been tested against single-factor experiments. We applied 10 TBMs to the multifactor Prairie Heating and CO2 Enrichment (PHACE) experiment in Wyoming, USA. Our goals were to investigate how multifactor experiments can be used to constrain models and to identify a road map for model improvement. We found models performed poorly in ambient conditions; there was a wide spread in simulated above-ground net primary productivity (range: 31-390 g C m-2  yr-1 ). Comparison with data highlighted model failures particularly with respect to carbon allocation, phenology, and the impact of water stress on phenology. Performance against the observations from single-factors treatments was also relatively poor. In addition, similar responses were predicted for different reasons across models: there were large differences among models in sensitivity to water stress and, among the N cycle models, N availability during the experiment. Models were also unable to capture observed treatment effects on phenology: they overestimated the effect of warming on leaf onset and did not allow CO2 -induced water savings to extend the growing season length. Observed interactive (CO2  × warming) treatment effects were subtle and contingent on water stress, phenology, and species composition. As the models did not correctly represent these processes under ambient and single-factor conditions, little extra information was gained by comparing model predictions against interactive responses. We outline a series of key areas in which this and future experiments could be used to improve model predictions of grassland responses to global change.

Gross primary production responses to warming, elevated CO2 , and irrigation: quantifying the drivers of ecosystem physiology in a semiarid grassland.

Determining whether the terrestrial biosphere will be a source or sink of carbon (C) under a future climate of elevated CO2 (eCO2 ) and warming requires accurate quantification of gross primary production (GPP), the largest flux of C in the global C cycle. We evaluated 6 years (2007-2012) of flux-derived GPP data from the Prairie Heating and CO2 Enrichment (PHACE) experiment, situated in a grassland in Wyoming, USA. The GPP data were used to calibrate a light response model whose basic formulation has been successfully used in a variety of ecosystems. The model was extended by modeling maximum photosynthetic rate (Amax ) and light-use efficiency (Q) as functions of soil water, air temperature, vapor pressure deficit, vegetation greenness, and nitrogen at current and antecedent (past) timescales. The model fits the observed GPP well (R2  = 0.79), which was confirmed by other model performance checks that compared different variants of the model (e.g. with and without antecedent effects). Stimulation of cumulative 6-year GPP by warming (29%, P = 0.02) and eCO2 (26%, P = 0.07) was primarily driven by enhanced C uptake during spring (129%, P = 0.001) and fall (124%, P = 0.001), respectively, which was consistent across years. Antecedent air temperature (Tairant ) and vapor pressure deficit (VPDant ) effects on Amax (over the past 3-4 days and 1-3 days, respectively) were the most significant predictors of temporal variability in GPP among most treatments. The importance of VPDant suggests that atmospheric drought is important for predicting GPP under current and future climate; we highlight the need for experimental studies to identify the mechanisms underlying such antecedent effects. Finally, posterior estimates of cumulative GPP under control and eCO2 treatments were tested as a benchmark against 12 terrestrial biosphere models (TBMs). The narrow uncertainties of these data-driven GPP estimates suggest that they could be useful semi-independent data streams for validating TBMs.

Patterns of Use, Acute Subjective Experiences, and Motivations for Using Synthetic Cathinones ("Bath Salts") in Recreational Users.

Given the variety and potential toxicity of synthetic cathinones, clinicians and educators would benefit from information about patterns of and motivations for use, frequency of psychosocial consequences, and experience of acute subjective effects. We administered a comprehensive, web-based survey to 104 recreational users of synthetic cathinones. Sixty percent of respondents consumed synthetic cathinones once or more per month, usually snorting or swallowing these drugs, typically at home, usually with others, customarily during the evening and nighttime hours, and often in combination with another drug such as alcohol or marijuana. Acute subjective effects attributed to synthetic cathinones were similar to those of other psychostimulants, including increased energy, rapid heartbeat, racing thoughts, difficulty sleeping, euphoria, decreased appetite, open-mindedness, and increased sex drive. Reported reasons for using synthetic cathinones included its stimulating effects, curiosity, substitution for another drug, and being at a party/music event. Respondents had experienced an average of six negative consequences of using synthetic cathinones during the previous year (e.g., tolerance, neglecting responsibilities, personality change). In combination with previously published investigations, these findings increase our understanding of the reported rationales and outcomes of recreational use of synthetic cathinones.

Confronting model predictions of carbon fluxes with measurements of Amazon forests subjected to experimental drought.

Considerable uncertainty surrounds the fate of Amazon rainforests in response to climate change. Here, carbon (C) flux predictions of five terrestrial biosphere models (Community Land Model version 3.5 (CLM3.5), Ecosystem Demography model version 2.1 (ED2), Integrated BIosphere Simulator version 2.6.4 (IBIS), Joint UK Land Environment Simulator version 2.1 (JULES) and Simple Biosphere model version 3 (SiB3)) and a hydrodynamic terrestrial ecosystem model (the Soil-Plant-Atmosphere (SPA) model) were evaluated against measurements from two large-scale Amazon drought experiments. Model predictions agreed with the observed C fluxes in the control plots of both experiments, but poorly replicated the responses to the drought treatments. Most notably, with the exception of ED2, the models predicted negligible reductions in aboveground biomass in response to the drought treatments, which was in contrast to an observed c. 20% reduction at both sites. For ED2, the timing of the decline in aboveground biomass was accurate, but the magnitude was too high for one site and too low for the other. Three key findings indicate critical areas for future research and model development. First, the models predicted declines in autotrophic respiration under prolonged drought in contrast to measured increases at one of the sites. Secondly, models lacking a phenological response to drought introduced bias in the sensitivity of canopy productivity and respiration to drought. Thirdly, the phenomenological water-stress functions used by the terrestrial biosphere models to represent the effects of soil moisture on stomatal conductance yielded unrealistic diurnal and seasonal responses to drought.