Search

Accurate representation of the turbulent exchange of carbon, water, and heat between the land surface and the atmosphere is critical for modelling global energy, water, and carbon cycles in both future climate projections and weather forecasts. Evaluation of models' ability to do this is performed in a wide range of simulation environments, often without explicit consideration of the degree of observational constraint or uncertainty and typically without quantification of benchmark performance expectations. We describe a Model Intercomparison Project (MIP) that attempts to resolve these shortcomings, comparing the surface turbulent heat flux predictions of around 20 different land models provided with in situ meteorological forcing evaluated with measured surface fluxes using quality-controlled data from 170 eddy-covariance-based flux tower sites. Predictions from seven out-of-sample empirical models are used to quantify the information available to land models in their forcing data and so the potential for land model performance improvement. Sites with unusual behaviour, complicated processes, poor data quality, or uncommon flux magnitude are more difficult to predict for both mechanistic and empirical models, providing a means of fairer assessment of land model performance. When examining observational uncertainty, model performance does not appear to improve in low-turbulence periods or with energy-balance-corrected flux tower data, and indeed some results raise questions about whether the energy balance correction process itself is appropriate. In all cases the results are broadly consistent, with simple out-of-sample empirical models, including linear regression, comfortably outperforming mechanistic land models. In all but two cases, latent heat flux and net ecosystem exchange of CO2 are better predicted by land models than sensible heat flux, despite it seeming to have fewer physical controlling processes. Land models that are implemented in Earth system models also appear to perform notably better than stand-alone ecosystem (including demographic) models, at least in terms of the fluxes examined here. The approach we outline enables isolation of the locations and conditions under which model developers can know that a land model can improve, allowing information pathways and discrete parameterisations in models to be identified and targeted for future model development.

Land biosphere processes are of central importance to the climate system. Specifically, ecosystems interact with the atmosphere through a variety of feedback loops that modulate energy, water, and CO2 fluxes between the land surface and the atmosphere across a wide range of temporal and spatial scales. Human land use and land cover modification add a further level of complexity to land–atmosphere interactions. Dynamic global vegetation models (DGVMs) attempt to capture land ecosystem processes and are increasingly incorporated into Earth system models (ESMs), which makes it possible to study the coupled dynamics of the land biosphere and the climate. In this work we describe a number of modifications to the LPJ-GUESS DGVM, aimed at enabling direct integration into an ESM. These include energy balance closure, the introduction of a sub-daily time step, a new radiative transfer scheme, and improved soil physics. The implemented modifications allow the model (LPJ-GUESS/LSM) to simulate the diurnal exchange of energy, water, and CO2 between the land ecosystem and the atmosphere and thus provide surface boundary conditions to an atmospheric model over land. A site-based evaluation against FLUXNET2015 data shows reasonable agreement between observed and modelled sensible and latent heat fluxes. Differences in predicted ecosystem function between standard LPJ-GUESS and LPJ-GUESS/LSM vary across land cover types. We find that the emerging ecosystem composition and carbon fluxes are sensitive to both the choice of stomatal conductance model and the response of plant water uptake to soil moisture. The new implementation described in this work lays the foundation for using the well-established LPJ-GUESS DGVM as an alternative land surface model (LSM) in coupled land–biosphere–climate studies, where an accurate representation of ecosystem processes is essential.

Understanding uncertainties and sensitivities of projected ecosystem dynamics under environmental change is of immense value for research and climate change policy. Here, we analyze sensitivities (change in model outputs per unit change in inputs) and uncertainties (changes in model outputs scaled to uncertainty in inputs) of vegetation dynamics under climate change, projected by a state-of-the-art dynamic vegetation model (LPJ-GUESS v4.0) across European forests (the species Picea abies, Fagus sylvatica and Pinus sylvestris), considering uncertainties of both model parameters and environmental drivers. We find that projected forest carbon fluxes are most sensitive to photosynthesis-, water-, and mortality-related parameters, while predictive uncertainties are dominantly induced by environmental drivers and parameters related to water and mortality. The importance of environmental drivers for predictive uncertainty increases with increasing temperature. Moreover, most of the interactions of model inputs (environmental drivers and parameters) are between environmental drivers themselves or between parameters and environmental drivers. In conclusion, our study highlights the importance of environmental drivers not only as contributors to predictive uncertainty in their own right but also as modifiers of sensitivities and thus uncertainties in other ecosystem processes. Reducing uncertainty in mortality-related processes and accounting for environmental influence on processes should therefore be a focus in further model development.

The Earth system model EC-Earth3 for contributions to CMIP6 is documented here, with its flexible coupling framework, major model configurations, a methodology for ensuring the simulations are comparable across different high-performance computing (HPC) systems, and with the physical performance of base configurations over the historical period. The variety of possible configurations and sub-models reflects the broad interests in the EC-Earth community. EC-Earth3 key performance metrics demonstrate physical behavior and biases well within the frame known from recent CMIP models. With improved physical and dynamic features, new Earth system model (ESM) components, community tools, and largely improved physical performance compared to the CMIP5 version, EC-Earth3 represents a clear step forward for the only European community ESM. We demonstrate here that EC-Earth3 is suited for a range of tasks in CMIP6 and beyond.

Improved agricultural management plays a vital role in protecting soils from degradation in eastern Africa. Changing practices such as reducing tillage, fertilizer use, or cover crops are expected to enhance soil organic carbon (SOC) storage, with climate change mitigation co-benefits, while increasing crop production. However, the quantification of cropland management effects on agricultural ecosystems remains inadequate in this region. Here, we explored seven management practices and their potential effects on soil carbon (C) pools, nitrogen (N) losses, and crop yields under different climate scenarios, using the dynamic vegetation model LPJ-GUESS. The model performance is evaluated against observations from two long-term maize field trials in western Kenya and reported estimates from published sources. LPJ-GUESS generally produces soil C stocks and maize productivity comparable with measurements and mostly captures the SOC decline under some management practices that is observed in the field experiments. We found that for large parts of Kenya and Ethiopia, an integrated conservation agriculture practice (no-tillage, residue and manure application, and cover crops) increases SOC levels in the long term (+11 % on average), accompanied by increased crop yields (+22 %) in comparison to the conventional management. Planting nitrogen-fixing cover crops in our simulations is also identified as a promising individual practice in eastern Africa to increase soil C storage (+4 %) and crop production (+18 %), with low environmental cost of N losses (+24 %). These management impacts are also sustained in simulations of three future climate pathways. This study highlights the possibilities of conservation agriculture when targeting long-term environmental sustainability and food security in crop ecosystems, particularly for those with poor soil conditions in tropical climates.

Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere in a changing climate is critical to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe and synthesize datasets and methodology to quantify the five major components of the global carbon budget and their uncertainties. Fossil CO2 emissions (EFOS) are based on energy statistics and cement production data, while emissions from land-use change (ELUC), mainly deforestation, are based on land use and land-use change data and bookkeeping models. Atmospheric CO2 concentration is measured directly, and its growth rate (GATM) is computed from the annual changes in concentration. The ocean CO2 sink (SOCEAN) is estimated with global ocean biogeochemistry models and observation-based data products. The terrestrial CO2 sink (SLAND) is estimated with dynamic global vegetation models. The resulting carbon budget imbalance (BIM), the difference between the estimated total emissions and the estimated changes in the atmosphere, ocean, and terrestrial biosphere, is a measure of imperfect data and understanding of the contemporary carbon cycle. All uncertainties are reported as ±1σ. For the first time, an approach is shown to reconcile the difference in our ELUC estimate with the one from national greenhouse gas inventories, supporting the assessment of collective countries' climate progress. For the year 2020, EFOS declined by 5.4 % relative to 2019, with fossil emissions at 9.5 ± 0.5 GtC yr−1 (9.3 ± 0.5 GtC yr−1 when the cement carbonation sink is included), and ELUC was 0.9 ± 0.7 GtC yr−1, for a total anthropogenic CO2 emission of 10.2 ± 0.8 GtC yr−1 (37.4 ± 2.9 GtCO2). Also, for 2020, GATM was 5.0 ± 0.2 GtC yr−1 (2.4 ± 0.1 ppm yr−1), SOCEAN was 3.0 ± 0.4 GtC yr−1, and SLAND was 2.9 ± 1 GtC yr−1, with a BIM of −0.8 GtC yr−1. The global atmospheric CO2 concentration averaged over 2020 reached 412.45 ± 0.1 ppm. Preliminary data for 2021 suggest a rebound in EFOS relative to 2020 of +4.8 % (4.2 % to 5.4 %) globally. Overall, the mean and trend in the components of the global carbon budget are consistently estimated over the period 1959–2020, but discrepancies of up to 1 GtC yr−1 persist for the representation of annual to semi-decadal variability in CO2 fluxes. Comparison of estimates from multiple approaches and observations shows (1) a persistent large uncertainty in the estimate of land-use changes emissions, (2) a low agreement between the different methods on the magnitude of the land CO2 flux in the northern extra-tropics, and (3) a discrepancy between the different methods on the strength of the ocean sink over the last decade. This living data update documents changes in the methods and datasets used in this new global carbon budget and the progress in understanding of the global carbon cycle compared with previous publications of this dataset (Friedlingstein et al., 2020, 2019; Le Quéré et al., 2018b, a, 2016, 2015b, a, 2014, 2013). The data presented in this work are available at https://doi.org/10.18160/gcp-2021 (Friedlingstein et al., 2021).

Biological nitrogen fixation (BNF) from grain legumes is of significant importance in global agricultural ecosystems. Crops with BNF capability are expected to support the need to increase food production while reducing nitrogen (N) fertilizer input for agricultural sustainability, but quantification of N fixing rates and BNF crop yields remains inadequate on a global scale. Here we incorporate two legume crops (soybean and faba bean) with BNF into a dynamic vegetation model LPJ-GUESS (Lund–Potsdam–Jena General Ecosystem Simulator). The performance of this new implementation is evaluated against observations from a range of water and N management trials. LPJ-GUESS generally captures the observed response to these management practices for legume biomass production, soil N uptake, and N fixation, despite some deviations from observations in some cases. Globally, simulated BNF is dominated by soil moisture and temperature, as well as N fertilizer addition. Annual inputs through BNF are modeled to be 11.6±2.2 Tg N for soybean and 5.6±1.0 Tg N for all pulses, with a total fixation of 17.2±2.9 Tg N yr−1 for all grain legumes during the period 1981–2016 on a global scale. Our estimates show good agreement with some previous statistical estimates but are relatively high compared to some estimates for pulses. This study highlights the importance of accounting for legume N fixation process when modeling C–N interactions in agricultural ecosystems, particularly when it comes to accounting for the combined effects of climate and land-use change on the global terrestrial N cycle.

In 2018 and 2019, central Europe was affected by two consecutive extreme dry and hot summers (DH18 and DH19). The DH18 event had severe impacts on ecosystems and likely affected vegetation activity in the subsequent year, for example through depletion of carbon reserves or damage from drought. Such legacies from drought and heat stress can further increase vegetation susceptibility to additional hazards. Temporally compound extremes such as DH18 and DH19 can, therefore, result in an amplification of impacts due to preconditioning effects of past disturbance legacies. Here, we evaluate how these two consecutive extreme summers impacted ecosystems in central Europe and how the vegetation responses to the first compound event (DH18) modulated the impacts of the second (DH19). To quantify changes in vegetation vulnerability to each compound event, we first train a set of statistical models for the period 2001–2017, which are then used to predict the impacts of DH18 and DH19 on enhanced vegetation index (EVI) anomalies from MODIS. These estimates correspond to expected EVI anomalies in DH18 and DH19 based on past sensitivity to climate. Large departures from the predicted values can indicate changes in vulnerability to dry and hot conditions and be used to identify modulating effects by vegetation activity and composition or other environmental factors on observed impacts. We find two regions in which the impacts of the two compound dry and hot (DH) events were significantly stronger than those expected based on previous climate–vegetation relationships. One region, largely dominated by grasslands and crops, showed much stronger impacts than expected in both DH events due to an amplification of their sensitivity to heat and drought, possibly linked to changing background CO2 and temperature conditions. A second region, dominated by forests and grasslands, showed browning from DH18 to DH19, even though dry and hot conditions were partly alleviated in 2019. This browning trajectory was mainly explained by the preconditioning role of DH18 on the impacts of DH19 due to interannual legacy effects and possibly by increased susceptibility to biotic disturbances, which are also promoted by warm conditions. Dry and hot summers are expected to become more frequent in the coming decades, posing a major threat to the stability of European forests. We show that state-of-the-art process-based models could not represent the decline in response to DH19 because they missed the interannual legacy effects from DH18 impacts. These gaps may result in an overestimation of the resilience and stability of temperate ecosystems in future model projections.

Quantifying the net carbon flux from land use and land cover changes (fLULCC) is critical for understanding the global carbon cycle and, hence, to support climate change mitigation. However, large-scale fLULCC is not directly measurable and has to be inferred from models instead, such as semi-empirical bookkeeping models and process-based dynamic global vegetation models (DGVMs). By definition, fLULCC estimates are not directly comparable between these two different model types. As an important example, DGVM-based fLULCC in the annual global carbon budgets is estimated under transient environmental forcing and includes the so-called loss of additional sink capacity (LASC). The LASC results from the impact of environmental changes on land carbon storage potential of managed land compared to potential vegetation and accumulates over time, which is not captured in bookkeeping models. The fLULCC from transient DGVM simulations, thus, strongly depends on the timing of land use and land cover changes mainly because LASC accumulation is cut off at the end of the simulated period. To estimate the LASC, the fLULCC from pre-industrial DGVM simulations, which is independent of changing environmental conditions, can be used. Additionally, DGVMs using constant present-day environmental forcing enable an approximation of bookkeeping estimates. Here, we analyse these three DGVM-derived fLULCC estimations (under transient, pre-industrial, and present-day forcing) for 12 models within 18 regions and quantify their differences as well as climate- and CO2-induced components and compare them to bookkeeping estimates. Averaged across the models, we find a global fLULCC (under transient conditions) of 2.0±0.6 PgC yr−1 for 2009–2018, of which ∼40 % are attributable to the LASC (0.8±0.3 PgC yr−1). From 1850 onward, the fLULCC accumulated to 189±56 PgC with 40±15 PgC from the LASC. Around 1960, the accumulating nature of the LASC causes global transient fLULCC estimates to exceed estimates under present-day conditions, despite generally increased carbon stocks in the latter. Regional hotspots of high cumulative and annual LASC values are found in the USA, China, Brazil, equatorial Africa, and Southeast Asia, mainly due to deforestation for cropland. Distinct negative LASC estimates in Europe (early reforestation) and from 2000 onward in the Ukraine (recultivation of post-Soviet abandoned agricultural land), indicate that fLULCC estimates in these regions are lower in transient DGVM compared to bookkeeping approaches. Our study unravels the strong dependence of fLULCC estimates on the time a certain land use and land cover change event happened to occur and on the chosen time period for the forcing of environmental conditions in the underlying simulations. We argue for an approach that provides an accounting of the fLULCC that is more robust against these choices, for example by estimating a mean DGVM ensemble fLULCC and LASC for a defined reference period and homogeneous environmental changes (CO2 only).

Land-use models and integrated assessment models provide scenarios of land-use and land-cover (LULC) changes following pathways or storylines related to different socioeconomic and environmental developments. The large diversity of available scenario projections leads to a recognizable variability in impacts on land ecosystems and the levels of services provided. We evaluated 16 projections of future LULC until 2040 that reflected different assumptions regarding socioeconomic demands and modeling protocols. By using these LULC projections in a state-of-the-art dynamic global vegetation model, we simulated their effect on selected ecosystem service indicators related to ecosystem productivity and carbon sequestration potential, agricultural production and the water cycle. We found that although a common trend for agricultural expansion exists across the scenarios, where and how particular LULC changes are realized differs widely across models and scenarios. They are linked to model-specific considerations of some demands over others and their respective translation into LULC changes and also reflect the simplified or missing representation of processes related to land dynamics or other influencing factors (e.g., trade, climate change). As a result, some scenarios show questionable and possibly unrealistic features in their LULC allocations, including highly regionalized LULC changes with rates of conversion that are contrary to or exceed rates observed in the past. Across the diverging LULC projections, we identified positive global trends of net primary productivity (+10.2 % ± 1.4 %), vegetation carbon (+9.2 % ± 4.1 %), crop production (+31.2 % ± 12.2 %) and water runoff (+9.3 % ± 1.7 %), and a negative trend of soil and litter carbon stocks (−0.5 % ± 0.4 %). The variability in ecosystem service indicators across scenarios was especially high for vegetation carbon stocks and crop production. Regionally, variability was highest in tropical forest regions, especially at current forest boundaries, because of intense and strongly diverging LULC change projections in combination with high vegetation productivity dampening or amplifying the effects of climatic change. Our results emphasize that information on future changes in ecosystem functioning and the related ecosystem service indicators should be seen in light of the variability originating from diverging projections of LULC. This is necessary to allow for adequate policy support towards sustainable transformations.

Changes in forest cover have a strong effect on climate through the alteration of surface biogeophysical and biogeochemical properties that affect energy, water and carbon exchange with the atmosphere. To quantify biogeophysical and biogeochemical effects of deforestation in a consistent setup, nine Earth system models (ESMs) carried out an idealized experiment in the framework of the Coupled Model Intercomparison Project, phase 6 (CMIP6). Starting from their pre-industrial state, models linearly replace 20×106 km2 of forest area in densely forested regions with grasslands over a period of 50 years followed by a stabilization period of 30 years. Most of the deforested area is in the tropics, with a secondary peak in the boreal region. The effect on global annual near-surface temperature ranges from no significant change to a cooling by 0.55 ∘C, with a multi-model mean of -0.22±0.21 ∘C. Five models simulate a temperature increase over deforested land in the tropics and a cooling over deforested boreal land. In these models, the latitude at which the temperature response changes sign ranges from 11 to 43∘ N, with a multi-model mean of 23∘ N. A multi-ensemble analysis reveals that the detection of near-surface temperature changes even under such a strong deforestation scenario may take decades and thus longer than current policy horizons. The observed changes emerge first in the centre of deforestation in tropical regions and propagate edges, indicating the influence of non-local effects. The biogeochemical effect of deforestation are land carbon losses of 259±80 PgC that emerge already within the first decade. Based on the transient climate response to cumulative emissions (TCRE) this would yield a warming by 0.46 ± 0.22 ∘C, suggesting a net warming effect of deforestation. Lastly, this study introduces the “forest sensitivity” (as a measure of climate or carbon change per fraction or area of deforestation), which has the potential to provide lookup tables for deforestation–climate emulators in the absence of strong non-local climate feedbacks. While there is general agreement across models in their response to deforestation in terms of change in global temperatures and land carbon pools, the underlying changes in energy and carbon fluxes diverge substantially across models and geographical regions. Future analyses of the global deforestation experiments could further explore the effect on changes in seasonality of the climate response as well as large-scale circulation changes to advance our understanding and quantification of deforestation effects in the ESM frameworks.

A future of increasing atmospheric carbon dioxide concentrations, changing climate, growing human populations, and shifting socioeconomic conditions means that the global agricultural system will need to adapt in order to feed the world. These changes will affect not only agricultural land but terrestrial ecosystems in general. Here, we use the coupled land use and vegetation model LandSyMM (Land System Modular Model) to quantify future land use change (LUC) and resulting impacts on ecosystem service indicators relating to carbon sequestration, runoff, biodiversity, and nitrogen pollution. We additionally hold certain variables, such as climate or land use, constant to assess the relative contribution of different drivers to the projected impacts. Some ecosystem services depend critically on land use and management: for example, carbon storage, the gain in which is more than 2.5 times higher in a low-LUC scenario (Shared Socioeconomic Pathway 4 and Representative Concentration Pathway 6.0; SSP4-60) than a high-LUC one with the same carbon dioxide and climate trajectory (SSP3-60). Other trends are mostly dominated by the direct effects of climate change and carbon dioxide increase. For example, in those two scenarios, extreme high monthly runoff increases across 54 % and 53 % of land, respectively, with a mean increase of 23 % in both. Scenarios in which climate change mitigation is more difficult (SSPs 3 and 5) have the strongest impacts on ecosystem service indicators, such as a loss of 13 %–19 % of land in biodiversity hotspots and a 28 % increase in nitrogen pollution. Evaluating a suite of ecosystem service indicators across scenarios enables the identification of tradeoffs and co-benefits associated with different climate change mitigation and adaptation strategies and socioeconomic developments.

FLUXNET comprises globally distributed eddy-covariance-based estimates of carbon fluxes between the biosphere and the atmosphere. Since eddy covariance flux towers have a relatively small footprint and are distributed unevenly across the world, upscaling the observations is necessary to obtain global-scale estimates of biosphere–atmosphere exchange. Based on cross-consistency checks with atmospheric inversions, sun-induced fluorescence (SIF) and dynamic global vegetation models (DGVMs), here we provide a systematic assessment of the latest upscaling efforts for gross primary production (GPP) and net ecosystem exchange (NEE) of the FLUXCOM initiative, where different machine learning methods, forcing data sets and sets of predictor variables were employed. Spatial patterns of mean GPP are consistent across FLUXCOM and DGVM ensembles (R2>0.94 at 1∘ spatial resolution) while the majority of DGVMs show, for 70 % of the land surface, values outside the FLUXCOM range. Global mean GPP magnitudes for 2008–2010 from FLUXCOM members vary within 106 and 130 PgC yr−1 with the largest uncertainty in the tropics. Seasonal variations in independent SIF estimates agree better with FLUXCOM GPP (mean global pixel-wise R2∼0.75) than with GPP from DGVMs (mean global pixel-wise R2∼0.6). Seasonal variations in FLUXCOM NEE show good consistency with atmospheric inversion-based net land carbon fluxes, particularly for temperate and boreal regions (R2>0.92). Interannual variability of global NEE in FLUXCOM is underestimated compared to inversions and DGVMs. The FLUXCOM version which also uses meteorological inputs shows a strong co-variation in interannual patterns with inversions (R2=0.87 for 2001–2010). Mean regional NEE from FLUXCOM shows larger uptake than inversion and DGVM-based estimates, particularly in the tropics with discrepancies of up to several hundred grammes of carbon per square metre per year. These discrepancies can only partly be reconciled by carbon loss pathways that are implicit in inversions but not captured by the flux tower measurements such as carbon emissions from fires and water bodies. We hypothesize that a combination of systematic biases in the underlying eddy covariance data, in particular in tall tropical forests, and a lack of site history effects on NEE in FLUXCOM are likely responsible for the too strong tropical carbon sink estimated by FLUXCOM. Furthermore, as FLUXCOM does not account for CO2 fertilization effects, carbon flux trends are not realistic. Overall, current FLUXCOM estimates of mean annual and seasonal cycles of GPP as well as seasonal NEE variations provide useful constraints of global carbon cycling, while interannual variability patterns from FLUXCOM are valuable but require cautious interpretation. Exploring the diversity of Earth observation data and of machine learning concepts along with improved quality and quantity of flux tower measurements will facilitate further improvements of the FLUXCOM approach overall.

Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere – the “global carbon budget” – 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 methodology to quantify the five major components of the global carbon budget and their uncertainties. Fossil CO2 emissions (EFF) are based on energy statistics and cement production data, while emissions from land use change (ELUC), mainly deforestation, are based on land use and land use change data and bookkeeping models. Atmospheric CO2 concentration is measured directly and its growth rate (GATM) is computed from the annual changes in concentration. The ocean CO2 sink (SOCEAN) and terrestrial CO2 sink (SLAND) are estimated with global process models constrained by observations. The resulting carbon budget imbalance (BIM), the difference between the estimated total emissions and the estimated changes in the atmosphere, ocean, and terrestrial biosphere, is a measure of imperfect data and understanding of the contemporary carbon cycle. All uncertainties are reported as ±1σ. For the last decade available (2009–2018), EFF was 9.5±0.5 GtC yr−1, ELUC 1.5±0.7 GtC yr−1, GATM 4.9±0.02 GtC yr−1 (2.3±0.01 ppm yr−1), SOCEAN 2.5±0.6 GtC yr−1, and SLAND 3.2±0.6 GtC yr−1, with a budget imbalance BIM of 0.4 GtC yr−1 indicating overestimated emissions and/or underestimated sinks. For the year 2018 alone, the growth in EFF was about 2.1 % and fossil emissions increased to 10.0±0.5 GtC yr−1, reaching 10 GtC yr−1 for the first time in history, ELUC was 1.5±0.7 GtC yr−1, for total anthropogenic CO2 emissions of 11.5±0.9 GtC yr−1 (42.5±3.3 GtCO2). Also for 2018, GATM was 5.1±0.2 GtC yr−1 (2.4±0.1 ppm yr−1), SOCEAN was 2.6±0.6 GtC yr−1, and SLAND was 3.5±0.7 GtC yr−1, with a BIM of 0.3 GtC. The global atmospheric CO2 concentration reached 407.38±0.1 ppm averaged over 2018. For 2019, preliminary data for the first 6–10 months indicate a reduced growth in EFF of +0.6 % (range of −0.2 % to 1.5 %) based on national emissions projections for China, the USA, the EU, and India and projections of gross domestic product corrected for recent changes in the carbon intensity of the economy for the rest of the world. Overall, the mean and trend in the five components of the global carbon budget are consistently estimated over the period 1959–2018, but discrepancies of up to 1 GtC yr−1 persist for the representation of semi-decadal variability in CO2 fluxes. A detailed comparison among individual estimates and the introduction of a broad range of observations shows (1) no consensus in the mean and trend in land use change emissions over the last decade, (2) a persistent low agreement between the different methods on the magnitude of the land CO2 flux in the northern extra-tropics, and (3) an apparent underestimation of the CO2 variability by ocean models outside the tropics. This living data update documents changes in the methods and data sets used in this new global carbon budget and the progress in understanding of the global carbon cycle compared with previous publications of this data set (Le Quéré et al., 2018a, b, 2016, 2015a, b, 2014, 2013). The data generated by this work are available at https://doi.org/10.18160/gcp-2019 (Friedlingstein et al., 2019).

Abstract Land use and climate changes both affect terrestrial ecosystems. Here, we used three combinations of Shared Socioeconomic Pathways and Representative Concentration Pathways (SSP1xRCP26, SSP3xRCP60, and SSP5xRCP85) as input to three dynamic global vegetation models to assess the impacts and associated uncertainty on several ecosystem functions: terrestrial carbon storage and fluxes, evapotranspiration, surface albedo, and runoff. We also performed sensitivity simulations in which we kept either land use or climate (including atmospheric CO2) constant from year 2015 on to calculate the isolated land use versus climate effects. By the 2080–2099 period, carbon storage increases by up to 87 ± 47 Gt (SSP1xRCP26) compared to present day, with large spatial variance across scenarios and models. Most of the carbon uptake is attributed to drivers beyond future land use and climate change, particularly the lagged effects of historic environmental changes. Future climate change typically increases carbon stocks in vegetation but not soils, while future land use change causes carbon losses, even for net agricultural abandonment (SSP1xRCP26). Evapotranspiration changes are highly variable across scenarios, and models do not agree on the magnitude or even sign of change of the individual effects. A calculated decrease in January and July surface albedo (up to −0.021 ± 0.007 and −0.004 ± 0.004 for SSP5xRCP85) and increase in runoff (+67 ± 6 mm/year) is largely driven by climate change. Overall, our results show that future land use and climate change will both have substantial impacts on ecosystem functioning. However, future changes can often not be fully explained by these two drivers and legacy effects have to be considered.

Abstract Ecosystems continuously adapt to interacting environmental drivers that change over time. Consequently, the carbon balance of terrestrial ecosystem may presently still be affected by past anthropogenic disturbances (e.g., deforestation) and other environmental changes (e.g., climate change). However, even though such so‐called “legacy effects” are implicitly included in many carbon cycle modeling studies, they are typically not explicitly quantified and therefore scientists might not be aware of their long‐term importance. Here, we use the ecosystem model LPJ‐GUESS to quantify legacy effects for the 21st century and the respective contributions of the following environmental drivers: climate change, CO2 fertilization, land use change, wood harvest, nitrogen deposition, and nitrogen fertilization. According to our simulations, the combined legacy effects of historical (1850–2015) environmental changes result in a land carbon uptake of +126 Gt C over the future (2015–2099) period. This by far exceeds the impacts of future environmental changes (range −53 Gt C to +16 Gt C for three scenarios) and is comparable in magnitude to historical carbon losses (−154 Gt C). Legacy effects can mainly be attributed to ecosystems still adapting to historical increases in atmospheric CO2 (+65 Gt C) and nitrogen deposition (+33 Gt C), but long‐term vegetation regrowth following agricultural abandonment (+8 Gt C) and wood harvest (+19 Gt C) also play a role. The response of the biosphere to historical environmental changes dominates future terrestrial carbon cycling at least until midcentury. Legacy effects persist many decades after environmental changes occurred and need to be considered when interpreting changes and estimating terrestrial carbon uptake potentials.

To support the assessments of the Intergovernmental Science-Policy Platform on Biodiversity and Ecosystem Services (IPBES), the IPBES Expert Group on Scenarios and Models is carrying out an intercomparison of biodiversity and ecosystem services models using harmonized scenarios (BES-SIM). The goals of BES-SIM are (1) to project the global impacts of land-use and climate change on biodiversity and ecosystem services (i.e., nature's contributions to people) over the coming decades, compared to the 20th century, using a set of common metrics at multiple scales, and (2) to identify model uncertainties and research gaps through the comparisons of projected biodiversity and ecosystem services across models. BES-SIM uses three scenarios combining specific Shared Socio-economic Pathways (SSPs) and Representative Concentration Pathways (RCPs) – SSP1xRCP2.6, SSP3xRCP6.0, SSP5xRCP8.6 – to explore a wide range of land-use change and climate change futures. This paper describes the rationale for scenario selection, the process of harmonizing input data for land use, based on the second phase of the Land Use Harmonization Project (LUH2), and climate, the biodiversity and ecosystem services models used, the core simulations carried out, the harmonization of the model output metrics, and the treatment of uncertainty. The results of this collaborative modeling project will support the ongoing global assessment of IPBES, strengthen ties between IPBES and the Intergovernmental Panel on Climate Change (IPCC) scenarios and modeling processes, advise the Convention on Biological Diversity (CBD) on its development of a post-2020 strategic plans and conservation goals, and inform the development of a new generation of nature-centred scenarios.

Land management for carbon storage is discussed as being indispensable for climate change mitigation because of its large potential to remove carbon dioxide from the atmosphere, and to avoid further emissions from deforestation. However, the acceptance and feasibility of land-based mitigation projects depends on potential side effects on other important ecosystem functions and their services. Here, we use projections of future land use and land cover for different land-based mitigation options from two land-use models (IMAGE and MAgPIE) and evaluate their effects with a global dynamic vegetation model (LPJ-GUESS). In the land-use models, carbon removal was achieved either via growth of bioenergy crops combined with carbon capture and storage, via avoided deforestation and afforestation, or via a combination of both. We compare these scenarios to a reference scenario without land-based mitigation and analyse the LPJ-GUESS simulations with the aim of assessing synergies and trade-offs across a range of ecosystem service indicators: carbon storage, surface albedo, evapotranspiration, water runoff, crop production, nitrogen loss, and emissions of biogenic volatile organic compounds. In our mitigation simulations cumulative carbon storage by year 2099 ranged between 55 and 89 GtC. Other ecosystem service indicators were influenced heterogeneously both positively and negatively, with large variability across regions and land-use scenarios. Avoided deforestation and afforestation led to an increase in evapotranspiration and enhanced emissions of biogenic volatile organic compounds, and to a decrease in albedo, runoff, and nitrogen loss. Crop production could also decrease in the afforestation scenarios as a result of reduced crop area, especially for MAgPIE land-use patterns, if assumed increases in crop yields cannot be realized. Bioenergy-based climate change mitigation was projected to affect less area globally than in the forest expansion scenarios, and resulted in less pronounced changes in most ecosystem service indicators than forest-based mitigation, but included a possible decrease in nitrogen loss, crop production, and biogenic volatile organic compounds emissions.

Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere – the “global carbon budget” – 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 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 and consistency within and among components, alongside methodology and data limitations. CO2 emissions from fossil fuels and industry (EFF) are based on energy statistics and cement production data, respectively, 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 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. We compare the mean land and ocean fluxes and their variability to estimates from three atmospheric inverse methods for three broad latitude bands. 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 (2006–2015), EFF was 9.3 ± 0.5 GtC yr−1, ELUC 1.0 ± 0.5 GtC yr−1, GATM 4.5 ± 0.1 GtC yr−1, SOCEAN 2.6 ± 0.5 GtC yr−1, and SLAND 3.1 ± 0.9 GtC yr−1. For year 2015 alone, the growth in EFF was approximately zero and emissions remained at 9.9 ± 0.5 GtC yr−1, showing a slowdown in growth of these emissions compared to the average growth of 1.8 % yr−1 that took place during 2006–2015. Also, for 2015, ELUC was 1.3 ± 0.5 GtC yr−1, GATM was 6.3 ± 0.2 GtC yr−1, SOCEAN was 3.0 ± 0.5 GtC yr−1, and SLAND was 1.9 ± 0.9 GtC yr−1. GATM was higher in 2015 compared to the past decade (2006–2015), reflecting a smaller SLAND for that year. The global atmospheric CO2 concentration reached 399.4 ± 0.1 ppm averaged over 2015. For 2016, preliminary data indicate the continuation of low growth in EFF with +0.2 % (range of −1.0 to +1.8 %) based on national emissions projections for China and USA, and projections of gross domestic product corrected for recent changes in the carbon intensity of the economy for the rest of the world. In spite of the low growth of EFF in 2016, the growth rate in atmospheric CO2 concentration is expected to be relatively high because of the persistence of the smaller residual terrestrial sink (SLAND) in response to El Niño conditions of 2015–2016. From this projection of EFF and assumed constant ELUC for 2016, cumulative emissions of CO2 will reach 565 ± 55 GtC (2075 ± 205 GtCO2) for 1870–2016, about 75 % from EFF and 25 % from ELUC. This living data update documents changes in the methods and data sets used in this new carbon budget compared with previous publications of this data set (Le Quéré et al., 2015b, a, 2014, 2013). All observations presented here can be downloaded from the Carbon Dioxide Information Analysis Center (doi:10.3334/CDIAC/GCP_2016).

Southeast Asia is a region known for active land-use changes (LUC) over the past 60 years; yet, how trends in net CO2 uptake and release resulting from LUC activities (net LUC flux) have changed through past decades remains uncertain. The level of uncertainty in net LUC flux from process-based models is so high that it cannot be concluded that newer estimates are necessarily more reliable than older ones. Here, we examined net LUC flux estimates of Southeast Asia for the 1980s−2010s from older and newer sets of Dynamic Global Vegetation Model simulations (TRENDY v2 and v7, respectively), and forcing data used for running those simulations, along with two book-keeping estimates (H&N and BLUE). These estimates yielded two contrasting historical LUC transitions, such that TRENDY v2 and H&N showed a transition from increased emissions from the 1980s to 1990s to declining emissions in the 2000s, while TRENDY v7 and BLUE showed the opposite transition. We found that these contrasting transitions originated in the update of LUC forcing data, which reduced the loss of forest area during the 1990s. Further evaluation of remote sensing studies, atmospheric inversions, and the history of forestry and environmental policies in Southeast Asia supported the occurrence of peak emissions in the 1990s and declining thereafter. However, whether LUC emissions continue to decline in Southeast Asia remains uncertain as key processes in recent years, such as conversion of peat forest to oil-palm plantation, are yet to be represented in the forcing data, suggesting a need for further revision.

The present study examined whether gender, ethnicity, instructional medium and school category differences manifest in science performance and motivation to learn science among secondary school students in Sri Lanka. The mean of five successive term test scores was used as the measure of science performance. Level of motivation in terms of six dimensions was measured by using the Science Motivation Questionnaire. A sample of 1316 grade 11 students representing Sinhala, Tamil, and Muslim ethnic groups from the three categories of public schools, which provide instructions in the vernacular languages of Sinhala or Tamil, participated in the study. Girls showed significantly higher performance in science compared to boys, and there was a significant gender difference in the levels of motivational dimensions in favor of girls. Although Tamil medium students possessed a higher level of motivation to learn science, Sinhala medium students outperformed their Tamil medium counterparts in science performance. Significant differences in science performance between Sinhala, Tamil, and Muslim students were also observed. However, motivation towards learning science between Tamils and Muslims was not significantly different. Highly significant differences in both motivation to learn science and performance in science were found between three categories of schools. The present study provides information to education officials who have to achieve equity across gender, ethnicity, medium of instruction, and school category, teachers who deliver the subject and school principals who design academic support programs.

"Swing away from science" is a dilemma that researchers are seeking solutions over 50 years. The dilemma is becoming more complicated due to the developmental decline in students' motivation toward learning science around the world, particularly at the secondary level. While the parental involvement becomes more vital in the secondary school years, it also dramatically declines in this period. The study investigated the level of parental involvement in terms of parents' socioeconomic factors and students' demographic factors. The association between parental involvement measured in terms of four dimensions and Sri Lankan secondary public school students' motivation toward learning science measured in terms of six dimensions was also studied. A random sample of 689 students and their parents participated in this quantitative survey. While the parental involvement showed significant differences in relation to the parents' income level and the school category that the child attends, the same showed no significant difference in relation to parents' ethnicity and child's gender. Parental involvement showed a strong association with students' intrinsically and extrinsically motivated science learning and self-efficacy. Home-based parental involvements had a stronger impact on students' motivation compared to the school-based involvements. Measures that can be taken by the stakeholders to strengthen the parental involvement activities, which enhance students' motivation to learn science and mitigate swing away from science were discussed.

This paper aims to address the impact of jet impingement on nano titanium-based ENi-P composite coatings deposited through electroless deposition on AH36 steel through Taguchi design of experiments. L9 Taguchi DOE was utilised to synthesise nine coatings by varying the four process parameter combinations suggested by the technique. Jet impingement combined with AFS 50/70 quartz sand erosion experimental run was conducted on the acidic, alkaline and nine ENi-P-TiO2 composite coatings administered through a waterjet erosion tester. The influence of C14-SB type Zwitterionic surfactant and nano TiO2 were discussed along with temperature and bath pH to achieve better resistance against the jet impingement on the surface of the coatings. Analysis of variance, the S/N ratio, and interaction plots were used to determine the optimal bath input parameters to deliver the least weight loss from the erosive environment and identify its use in the marine environment. The proposed optimal coating was validated using a confirmation test, and the weight loss was reduced by more than 35% compared to the initial condition parameter (A2B2C2D2). In addition to mechanical properties, the quality of the deposits was analysed through morphological and chemical composition probed through scanning electron microscope, energy dispersive X-Ray analysis and powder X-Ray diffraction methods. As a result, the Composite Modified Resilience (CMR) value of ENi-P-TiO2 optimal coating was improved by over 41.17%, and its corresponding weight loss was reduced to 35% compared to the initial condition. This experiment paves the way for further study on electroless deposits and their characteristics on AH36 steel generally used for shipbuilding in commercial ships, naval, and marine associated defence applications. [ABSTRACT FROM AUTHOR]