Constraining Uncertainty in Projected Gross Primary Production With Machine Learning

The terrestrial biosphere is currently slowing down global warming by absorbing about 30% of human emissions of carbon dioxide (CO2). The largest flux of the terrestrial carbon uptake is gross primary production (GPP) defined as the production of carbohydrates by photosynthesis. Elevated atmospheric CO2 concentration is expected to increase GPP (“CO2 fertilization effect”). However, Earth system models (ESMs) exhibit a large range in simulated GPP projections. In this study, we combine an existing emergent constraint on CO2 fertilization with a machine learning approach to constrain the spatial variations of multimodel GPP projections. In a first step, we use observed changes in the CO2 seasonal cycle at Cape Kumukahi to constrain the global mean GPP at the end of the 21st century (2091–2100) in Representative Concentration Pathway 8.5 simulations with ESMs participating in the Coupled Model Intercomparison Project Phase 5 (CMIP5) to 171 ± 12 Gt C yr, compared to the unconstrainedmodel range of 156–247 Gt C yr. In a second step, we use a machine learning model to constrain gridded future absolute GPP and gridded fractional GPP change in two independent approaches. For this, observational data are fed into the machine learning algorithm that has been trained on CMIP5 data to learn relationships between present‐day physically relevant diagnostics and the target variable. In a leave‐one‐model‐out cross‐validation approach, the machine learning model shows superior performance to the CMIP5 ensemble mean. Our approach predicts an increased GPP change in northern high latitudes compared to regions closer to the equator. Plain Language Summary About a quarter of human emissions of carbon dioxide (CO2) is absorbed by vegetation and soil on the Earth's surface and hence does not contribute to global warming caused by CO2 in the atmosphere. Thus, in order to better define the remaining carbon budgets left to meet particular warming targets like the 1.5°C of the Paris Agreement, it is important to accurately quantify the carbon uptake by plants in the future. Currently, this is modeled by Earth system models yet with great uncertainties. In this work, we present an alternative machine learning approach to reduce spatial uncertainties of vegetation carbon uptake in future climate projections using observations of today's conditions.

[1]  H. Damon Matthews,et al.  The proportionality of global warming to cumulative carbon emissions , 2009, Nature.

[2]  M. Salvucci,et al.  Rubisco activase constrains the photosynthetic potential of leaves at high temperature and CO2. , 2000, Proceedings of the National Academy of Sciences of the United States of America.

[3]  J. Gregory,et al.  Quantifying Carbon Cycle Feedbacks , 2009 .

[4]  B. Dawson,et al.  UNITED NATIONS FRAMEWORK CONVENTION ON CLIMATE CHANGE (UNFCCC) , 2008 .

[5]  Alexander J. Winkler,et al.  Earth system models underestimate carbon fixation by plants in the high latitudes , 2019, Nature Communications.

[6]  R. Betts,et al.  High sensitivity of future global warming to land carbon cycle processes , 2012 .

[7]  G. Danabasoglu,et al.  The Community Climate System Model Version 4 , 2011 .

[8]  P. Cox,et al.  Projected land photosynthesis constrained by changes in the seasonal cycle of atmospheric CO2 , 2016, Nature.

[9]  Atul K. Jain,et al.  Contrasting effects of CO2 fertilization, land-use change and warming on seasonal amplitude of Northern Hemisphere CO2 exchange , 2019, Atmospheric Chemistry and Physics.

[10]  C. Jones,et al.  Development and evaluation of an Earth-System model - HadGEM2 , 2011 .

[11]  Ramond Eeli Aa Ndrenelaprise,et al.  Forecasting Skill Limits of Nested, Limited-Area Models: A Perfect-Model Approach , 2002 .

[12]  Alexander J. Winkler,et al.  Investigating the applicability of emergent constraints , 2019, Earth System Dynamics.

[13]  D. Maraun,et al.  Improving Antarctic Total Ozone Projections by a Process-Oriented Multiple Diagnostic Ensemble Regression , 2013 .

[14]  Atul K. Jain,et al.  Role of CO2, climate and land use in regulating the seasonal amplitude increase of carbon fluxes in terrestrial ecosystems: a multimodel analysis , 2016, Biogeosciences.

[15]  J Elith,et al.  A working guide to boosted regression trees. , 2008, The Journal of animal ecology.

[16]  J. Thepaut,et al.  The ERA‐Interim reanalysis: configuration and performance of the data assimilation system , 2011 .

[17]  C. Tebaldi,et al.  Long-term Climate Change: Projections, Commitments and Irreversibility , 2013 .

[18]  I. C. Prentice,et al.  Climatic Control of the High-Latitude Vegetation Greening Trend and Pinatubo Effect , 2002, Science.

[19]  R. Knutti,et al.  Skill and independence weighting for multi-model assessments , 2016 .

[20]  Ivar A. Seierstad,et al.  The Norwegian Earth System Model, NorESM1-M – Part 2: Climate response and scenario projections , 2012 .

[21]  Radford M. Neal Pattern Recognition and Machine Learning , 2007, Technometrics.

[22]  K. Denman,et al.  Carbon emission limits required to satisfy future representative concentration pathways of greenhouse gases , 2011 .

[23]  Gaël Varoquaux,et al.  Scikit-learn: Machine Learning in Python , 2011, J. Mach. Learn. Res..

[24]  Karl E. Taylor,et al.  An overview of CMIP5 and the experiment design , 2012 .

[25]  Stephen Sitch,et al.  A roadmap for improving the representation of photosynthesis in Earth system models. , 2017, The New phytologist.

[26]  A M Michalak,et al.  Uncertainty in the response of terrestrial carbon sink to environmental drivers undermines carbon-climate feedback predictions , 2017, Scientific Reports.

[27]  A. Arneth,et al.  Global patterns of land-atmosphere fluxes of carbon dioxide, latent heat, and sensible heat derived from eddy covariance, satellite, and meteorological observations , 2011 .

[28]  Axel Lauer,et al.  Earth System Model Evaluation Tool (ESMValTool) v2.0 – diagnostics for emergent constraints and future projections from Earth system models in CMIP , 2020, Geoscientific Model Development.

[29]  P. Cox,et al.  Evaluating the Land and Ocean Components of the Global Carbon Cycle in the CMIP5 Earth System Models , 2013 .

[30]  E. A. Kort,et al.  Enhanced Seasonal Exchange of CO2 by Northern Ecosystems Since 1960 , 2013, Science.

[31]  Nuno Carvalhais,et al.  Enhanced seasonal CO2 exchange caused by amplified plant productivity in northern ecosystems , 2016, Science.

[32]  Philippe Ciais,et al.  On the causes of trends in the seasonal amplitude of atmospheric CO2 , 2018, Global change biology.

[33]  Klaus Zimmermann,et al.  Earth System Model Evaluation Tool (ESMValTool) v2.0 – technical overview , 2020, Geoscientific Model Development.

[34]  P. Ciais,et al.  Partitioning global land evapotranspiration using CMIP5 models constrained by observations , 2018, Nature Climate Change.

[35]  Yoav Freund,et al.  Experiments with a New Boosting Algorithm , 1996, ICML.

[36]  P. Cox,et al.  Emergent constraints on climate‐carbon cycle feedbacks in the CMIP5 Earth system models , 2014 .

[37]  Robert B. Jackson,et al.  Global patterns of terrestrial nitrogen and phosphorus limitation , 2020, Nature Geoscience.

[38]  P. Cox,et al.  Sensitivity of tropical carbon to climate change constrained by carbon dioxide variability , 2013, Nature.

[39]  G. De’ath Boosted trees for ecological modeling and prediction. , 2007, Ecology.

[40]  Corinne Le Quéré,et al.  Carbon and Other Biogeochemical Cycles , 2014 .

[41]  K.,et al.  Carbon–Concentration and Carbon–Climate Feedbacks in CMIP5 Earth System Models , 2012 .

[42]  William Woodgate,et al.  Higher than expected CO2 fertilization inferred from leaf to global observations , 2020, Global change biology.

[43]  P. Jones,et al.  Updated high‐resolution grids of monthly climatic observations – the CRU TS3.10 Dataset , 2014 .

[44]  Ruth Lorenz,et al.  A climate model projection weighting scheme accounting for performance and interdependence , 2017 .

[45]  King-Sun Fu,et al.  Pattern Recognition and Machine Learning , 2012 .

[46]  Nicholas C. Parazoo,et al.  Large and projected strengthening moisture limitation on end-of-season photosynthesis , 2020, Proceedings of the National Academy of Sciences.

[47]  Victor O. Leshyk,et al.  Integrating the evidence for a terrestrial carbon sink caused by increasing atmospheric CO2. , 2020, The New phytologist.

[48]  Veronika Eyring,et al.  Constraining Future Summer Austral Jet Stream Positions in the CMIP5 Ensemble by Process-Oriented Multiple Diagnostic Regression* , 2016 .

[49]  C. Tucker,et al.  Increased plant growth in the northern high latitudes from 1981 to 1991 , 1997, Nature.

[50]  N. Zeng,et al.  Continued increase in atmospheric CO 2 seasonal amplitude in the 21st century projected by the CMIP5 Earth system models , 2014 .

[51]  J. Friedman Greedy function approximation: A gradient boosting machine. , 2001 .

[52]  R. Schnur,et al.  Climate-carbon cycle feedback analysis: Results from the C , 2006 .

[53]  Nuno Carvalhais,et al.  Earth System Model Evaluation Tool (ESMValTool) v2.0 – an extended set of large-scale diagnostics for quasi-operational and comprehensive evaluation of Earth system models in CMIP , 2020, Geoscientific Model Development.

[54]  J. Friedman Stochastic gradient boosting , 2002 .

[55]  A. Lauer,et al.  Constraining Uncertainties in CMIP5 Projections of September Arctic Sea Ice Extent with Observations , 2020, Journal of Climate.

[56]  Atul K. Jain,et al.  Global Carbon Budget 2019 , 2019, Earth System Science Data.

[57]  Carlos Guestrin,et al.  "Why Should I Trust You?": Explaining the Predictions of Any Classifier , 2016, ArXiv.

[58]  B. Stevens,et al.  Climate and carbon cycle changes from 1850 to 2100 in MPI‐ESM simulations for the Coupled Model Intercomparison Project phase 5 , 2013 .

[59]  N. Nakicenovic,et al.  RCP 8.5—A scenario of comparatively high greenhouse gas emissions , 2011 .

[60]  Alessandro Anav,et al.  Global Data Sets of Vegetation Leaf Area Index (LAI)3g and Fraction of Photosynthetically Active Radiation (FPAR)3g Derived from Global Inventory Modeling and Mapping Studies (GIMMS) Normalized Difference Vegetation Index (NDVI3g) for the Period 1981 to 2011 , 2013, Remote. Sens..

[61]  Tianqi Chen,et al.  XGBoost: A Scalable Tree Boosting System , 2016, KDD.

[62]  C. D. Keeling,et al.  Increased activity of northern vegetation inferred from atmospheric CO2 measurements , 1996, Nature.

[63]  Peter J. Rayner,et al.  Uncertainty in temperature projections reduced using carbon cycle and climate observations , 2013 .

[64]  C. D. Keeling,et al.  Atmospheric CO2 and 13CO2 Exchange with the Terrestrial Biosphere and Oceans from 1978 to 2000: Observations and Carbon Cycle Implications , 2005 .

[65]  Bettina K. Gier,et al.  Taking climate model evaluation to the next level , 2019, Nature Climate Change.

[66]  Samuel B. Williams,et al.  ASSOCIATION FOR COMPUTING MACHINERY , 2000 .

[67]  M. Wahlen,et al.  Interannual extremes in the rate of rise of atmospheric carbon dioxide since 1980 , 1995, Nature.

[68]  Ronald,et al.  GFDL’s ESM2 Global Coupled Climate–Carbon Earth System Models. Part I: Physical Formulation and Baseline Simulation Characteristics , 2012 .

[69]  Mark A. Friedl,et al.  Direct human influence on atmospheric CO2 seasonality from increased cropland productivity , 2014, Nature.

[70]  Anja Rammig,et al.  Amazon forest response to CO2 fertilization dependent on plant phosphorus acquisition , 2019, Nature Geoscience.

[71]  R. Tibshirani Regression Shrinkage and Selection via the Lasso , 1996 .