Statistical upscaling of ecosystem CO <sub>2</sub> fluxes across the terrestrial tundra and boreal domain: Regional patterns and uncertainties

Anna‐Maria Virkkala, Juha Aalto, Brendan M. Rogers, Torbern Tagesson, Claire C. Treat, Susan M. Natali, Jennifer D. Watts, Stefano Potter, Aleksi Lehtonen, Marguerite Mauritz, Edward A. G. Schuur, John Kochendorfer, Donatella Zona, Walter C. Oechel, Hideki Kobayashi, Elyn Humphreys, Mathias Goeckede, Hiroyasu Iwata, Peter M. Lafleur, Eugénie Euskirchen, Stef Bokhorst, Maija E. Marushchak, Pertti J. Martikainen, Bo Elberling, Carolina Voigt, Christina Biasi, Oliver Sonnentag, Frans‐Jan W. Parmentier, Masahito Ueyama, Gerardo Celis, Vincent L. St. Louis, Craig A. Emmerton, Matthias Peichl, Jinshu Chi, Järvi Järveoja, Mats Nilsson, Steven F. Oberbauer, M. S. Torn, Sang Jong Park, Han Dolman, Ivan Mammarella, Namyi Chae, Rafael Poyatos, Efrèn López‐Blanco, Torben R. Christensen, Mi Hye Kwon, Torsten Sachs, David Holl, Miska Luoto

Abstract
The regional variability in tundra and boreal carbon dioxide (CO2) fluxes can be high, complicating efforts to quantify sink-source patterns across the entire region. Statistical models are increasingly used to predict (i.e., upscale) CO2 fluxes across large spatial domains, but the reliability of different modeling techniques, each with different specifications and assumptions, has not been assessed in detail. Here, we compile eddy covariance and chamber measurements of annual and growing season CO2 fluxes of gross primary productivity (GPP), ecosystem respiration (ER), and net ecosystem exchange (NEE) during 1990–2015 from 148 terrestrial high-latitude (i.e., tundra and boreal) sites to analyze the spatial patterns and drivers of CO2 fluxes and test the accuracy and uncertainty of different statistical models. CO2 fluxes were upscaled at relatively high spatial resolution (1 km2) across the high-latitude region using five commonly used statistical models and their ensemble, that is, the median of all five models, using climatic, vegetation, and soil predictors. We found the performance of machine learning and ensemble predictions to outperform traditional regression methods. We also found the predictive performance of NEE-focused models to be low, relative to models predicting GPP and ER. Our data compilation and ensemble predictions showed that CO2 sink strength was larger in the boreal biome (observed and predicted average annual NEE −46 and −29 g C m−2 yr−1, respectively) compared to tundra (average annual NEE +10 and −2 g C m−2 yr−1). This pattern was associated with large spatial variability, reflecting local heterogeneity in soil organic carbon stocks, climate, and vegetation productivity. The terrestrial ecosystem CO2 budget, estimated using the annual NEE ensemble prediction, suggests the high-latitude region was on average an annual CO2 sink during 1990–2015, although uncertainty remains high.
Cite:
Anna‐Maria Virkkala, Juha Aalto, Brendan M. Rogers, Torbern Tagesson, Claire C. Treat, Susan M. Natali, Jennifer D. Watts, Stefano Potter, Aleksi Lehtonen, Marguerite Mauritz, Edward A. G. Schuur, John Kochendorfer, Donatella Zona, Walter C. Oechel, Hideki Kobayashi, Elyn Humphreys, Mathias Goeckede, Hiroyasu Iwata, Peter M. Lafleur, et al.. 2021. Statistical upscaling of ecosystem CO 2 fluxes across the terrestrial tundra and boreal domain: Regional patterns and uncertainties. Global Change Biology, Volume 27, Issue 17, 27(17):4040–4059.
Copy Citation:
Cite Search