Manuel Helbig

2023

DOI bib abs
Carbon uptake in Eurasian boreal forests dominates the high‐latitude net ecosystem carbon budget
Jennifer D. Watts, Mary Farina, John S. Kimball, Luke D. Schiferl, Zhihua Liu, Kyle A. Arndt, Donatella Zona, Ashley P. Ballantyne, E. S. Euskirchen, Frans‐Jan W. Parmentier, Manuel Helbig, Oliver Sonnentag, Torbern Tagesson, Janne Rinne, Hiroki Ikawa, Masahito Ueyama, Hideki Kobayashi, Torsten Sachs, Daniel F. Nadeau, John Kochendorfer, M. Jackowicz-Korczyński, Anna Virkkala, Mika Aurela, R. Commane, Brendan Byrne, Leah Birch, Matthew S. Johnson, Nima Madani, Brendan M. Rogers, Jinyang Du, Arthur Endsley, K. E. Savage, Benjamin Poulter, Zhen Zhang, L. M. Bruhwiler, Charles E. Miller, S. J. Goetz, Walter C. Oechel, Jennifer D. Watts, Mary Farina, John S. Kimball, Luke D. Schiferl, Zhihua Liu, Kyle A. Arndt, Donatella Zona, Ashley P. Ballantyne, E. S. Euskirchen, Frans‐Jan W. Parmentier, Manuel Helbig, Oliver Sonnentag, Torbern Tagesson, Janne Rinne, Hiroki Ikawa, Masahito Ueyama, Hideki Kobayashi, Torsten Sachs, Daniel F. Nadeau, John Kochendorfer, M. Jackowicz-Korczyński, Anna Virkkala, Mika Aurela, R. Commane, Brendan Byrne, Leah Birch, Matthew S. Johnson, Nima Madani, Brendan M. Rogers, Jinyang Du, Arthur Endsley, K. E. Savage, Benjamin Poulter, Zhen Zhang, L. M. Bruhwiler, Charles E. Miller, S. J. Goetz, Walter C. Oechel
Global Change Biology, Volume 29, Issue 7

Arctic-boreal landscapes are experiencing profound warming, along with changes in ecosystem moisture status and disturbance from fire. This region is of global importance in terms of carbon feedbacks to climate, yet the sign (sink or source) and magnitude of the Arctic-boreal carbon budget within recent years remains highly uncertain. Here, we provide new estimates of recent (2003-2015) vegetation gross primary productivity (GPP), ecosystem respiration (Reco ), net ecosystem CO2 exchange (NEE; Reco - GPP), and terrestrial methane (CH4 ) emissions for the Arctic-boreal zone using a satellite data-driven process-model for northern ecosystems (TCFM-Arctic), calibrated and evaluated using measurements from >60 tower eddy covariance (EC) sites. We used TCFM-Arctic to obtain daily 1-km2 flux estimates and annual carbon budgets for the pan-Arctic-boreal region. Across the domain, the model indicated an overall average NEE sink of -850 Tg CO2 -C year-1 . Eurasian boreal zones, especially those in Siberia, contributed to a majority of the net sink. In contrast, the tundra biome was relatively carbon neutral (ranging from small sink to source). Regional CH4 emissions from tundra and boreal wetlands (not accounting for aquatic CH4 ) were estimated at 35 Tg CH4 -C year-1 . Accounting for additional emissions from open water aquatic bodies and from fire, using available estimates from the literature, reduced the total regional NEE sink by 21% and shifted many far northern tundra landscapes, and some boreal forests, to a net carbon source. This assessment, based on in situ observations and models, improves our understanding of the high-latitude carbon status and also indicates a continued need for integrated site-to-regional assessments to monitor the vulnerability of these ecosystems to climate change.

DOI bib abs
Reduced Net CO2 Uptake During Dry Summers in a Boreal Shield Peatland
Renée M. McDonald, Paul Moore, Manuel Helbig, J. M. Waddington, Renée M. McDonald, Paul Moore, Manuel Helbig, J. M. Waddington, Renée M. McDonald, Paul Moore, Manuel Helbig, J. M. Waddington, Renée M. McDonald, Paul Moore, Manuel Helbig, J. M. Waddington
Journal of Geophysical Research: Biogeosciences, Volume 128, Issue 2

Abstract Peatlands are globally important long‐term sinks of atmospheric carbon dioxide (CO 2 ). However, there is concern that climate change‐mediated drying will reduce gross primary productivity (GPP) and increase ecosystem respiration (ER) making peatlands vulnerable to a weaker carbon sink function and potential net carbon loss. While large and deep peatlands are usually resilient to moderate summer drying, CO 2 exchange in shallow Boreal Shield peatlands is likely more sensitive to drying given the reduced groundwater connectivity and water storage potential. To better understand the carbon cycling responses of Boreal Shield peatlands to meteorological conditions, we examined ecohydrological controls on CO 2 fluxes using the eddy covariance technique at a shallow peatland during the summer season for 5 years, from 2016–2020. We found lower GPP in dry summer years. Mean summer water table depth (WTD) was found to be significantly correlated with summer total net ecosystem CO 2 exchange ( R 2 = 0.78; p value = 0.046) and GPP ( R 2 = 0.83; p value = 0.03), where wet summers with a WT close to the peat surface sequestered more than twice the amount of CO 2 than dry summers. Our findings suggest that shallow Boreal Shield peatland GPP may be sensitive to climate‐mediated drying as they may switch to a net CO 2 source in the summer season when WTDs exceed a critical ecohydrological threshold for a prolonged period of time.

DOI bib abs
Upscaling Wetland Methane Emissions From the FLUXNET‐CH4 Eddy Covariance Network (UpCH4 v1.0): Model Development, Network Assessment, and Budget Comparison
Gavin McNicol, Etienne Fluet‐Chouinard, Zutao Ouyang, Sara Knox, Zhen Zhang, Tuula Aalto, Sheel Bansal, Kuang‐Yu Chang, Min Chen, Kyle Delwiche, Sarah Féron, Mathias Goeckede, Jinxun Liu, Avni Malhotra, Joe R. Melton, W. J. Riley, Rodrigo Vargas, Kunxiaojia Yuan, Qing Ying, Qing Zhu, Pavel Alekseychik, Mika Aurela, David P. Billesbach, David I. Campbell, Jiquan Chen, Housen Chu, Ankur R. Desai, E. S. Euskirchen, Jordan P. Goodrich, Timothy J. Griffis, Manuel Helbig, Takashi Hirano, Hiroki Iwata, Gerald Jurasinski, John S. King, Franziska Koebsch, Randall K. Kolka, Ken W. Krauss, Annalea Lohila, Ivan Mammarella, Mats E Nilson, Asko Noormets, Walter C. Oechel, Matthias Peichl, Torsten Sachs, Ayaka Sakabe, Christopher Schulze, Oliver Sonnentag, Ryan C. Sullivan, Eeva‐Stiina Tuittila, Masahito Ueyama, Timo Vesala, Eric J. Ward, Christian Wille, Guan Xhuan Wong, Donatella Zona, L. Windham‐Myers, Benjamin Poulter, Robert B. Jackson
AGU Advances, Volume 4, Issue 5

Abstract Wetlands are responsible for 20%–31% of global methane (CH 4 ) emissions and account for a large source of uncertainty in the global CH 4 budget. Data‐driven upscaling of CH 4 fluxes from eddy covariance measurements can provide new and independent bottom‐up estimates of wetland CH 4 emissions. Here, we develop a six‐predictor random forest upscaling model (UpCH4), trained on 119 site‐years of eddy covariance CH 4 flux data from 43 freshwater wetland sites in the FLUXNET‐CH4 Community Product. Network patterns in site‐level annual means and mean seasonal cycles of CH 4 fluxes were reproduced accurately in tundra, boreal, and temperate regions (Nash‐Sutcliffe Efficiency ∼0.52–0.63 and 0.53). UpCH4 estimated annual global wetland CH 4 emissions of 146 ± 43 TgCH 4 y −1 for 2001–2018 which agrees closely with current bottom‐up land surface models (102–181 TgCH 4 y −1 ) and overlaps with top‐down atmospheric inversion models (155–200 TgCH 4 y −1 ). However, UpCH4 diverged from both types of models in the spatial pattern and seasonal dynamics of tropical wetland emissions. We conclude that upscaling of eddy covariance CH 4 fluxes has the potential to produce realistic extra‐tropical wetland CH 4 emissions estimates which will improve with more flux data. To reduce uncertainty in upscaled estimates, researchers could prioritize new wetland flux sites along humid‐to‐arid tropical climate gradients, from major rainforest basins (Congo, Amazon, and SE Asia), into monsoon (Bangladesh and India) and savannah regions (African Sahel) and be paired with improved knowledge of wetland extent seasonal dynamics in these regions. The monthly wetland methane products gridded at 0.25° from UpCH4 are available via ORNL DAAC ( https://doi.org/10.3334/ORNLDAAC/2253 ).

2022

DOI bib abs
Warming response of peatland CO2 sink is sensitive to seasonality in warming trends
Manuel Helbig, Tatjana Živković, Pavel Alekseychik, Mika Aurela, Tarek S. El‐Madany, E. S. Euskirchen, Lawrence B. Flanagan, Timothy J. Griffis, Paul J. Hanson, J. Hattakka, Carole Helfter, Takashi Hirano, Elyn Humphreys, Gerard Kiely, Randall K. Kolka, Tuomas Laurila, Paul Leahy, Annalea Lohila, Ivan Mammarella, Mats B. Nilsson, A. V. Panov, Frans‐Jan W. Parmentier, Matthias Peichl, Janne Rinne, D. Tyler Roman, Oliver Sonnentag, Eeva‐Stiina Tuittila, Masahito Ueyama, Timo Vesala, Patrik Vestin, Simon Weldon, Per Weslien, Sönke Zaehle
Nature Climate Change, Volume 12, Issue 8

Peatlands have acted as net CO2 sinks over millennia, exerting a global climate cooling effect. Rapid warming at northern latitudes, where peatlands are abundant, can disturb their CO2 sink function. Here we show that sensitivity of peatland net CO2 exchange to warming changes in sign and magnitude across seasons, resulting in complex net CO2 sink responses. We use multiannual net CO2 exchange observations from 20 northern peatlands to show that warmer early summers are linked to increased net CO2 uptake, while warmer late summers lead to decreased net CO2 uptake. Thus, net CO2 sinks of peatlands in regions experiencing early summer warming, such as central Siberia, are more likely to persist under warmer climate conditions than are those in other regions. Our results will be useful to improve the design of future warming experiments and to better interpret large-scale trends in peatland net CO2 uptake over the coming few decades.

DOI bib abs
The ABCflux database: Arctic–boreal CO2 flux observations and ancillary information aggregated to monthly time steps across terrestrial ecosystems
Anna‐Maria Virkkala, Susan M. Natali, Brendan M. Rogers, Jennifer D. Watts, K. E. Savage, Sara June Connon, Marguerite Mauritz, Edward A. G. Schuur, D. L. Peter, Christina Minions, Julia Nojeim, R. Commane, Craig A. Emmerton, Mathias Goeckede, Manuel Helbig, David Holl, Hiroki Iwata, Hideki Kobayashi, Pasi Kolari, Efrèn López‐Blanco, Maija E. Marushchak, Mikhail Mastepanov, Lutz Merbold, Frans‐Jan W. Parmentier, Matthias Peichl, Torsten Sachs, Oliver Sonnentag, Masahito Ueyama, Carolina Voigt, Mika Aurela, Julia Boike, Gerardo Celis, Namyi Chae, Torben R. Christensen, M. Syndonia Bret‐Harte, Sigrid Dengel, A. J. Dolman, Colin W. Edgar, Bo Elberling, E. S. Euskirchen, Achim Grelle, Juha Hatakka, Elyn Humphreys, Järvi Järveoja, Ayumi Kotani, Lars Kutzbach, Tuomas Laurila, Annalea Lohila, Ivan Mammarella, Yojiro Matsuura, Gesa Meyer, Mats B. Nilsson, Steven F. Oberbauer, Sang‐Jong Park, Roman Petrov, Anatoly Prokushkin, Christopher Schulze, Vincent L. St. Louis, Eeva‐Stiina Tuittila, Juha‐Pekka Tuovinen, W. L. Quinton, Andrej Varlagin, Donatella Zona, Viacheslav I. Zyryanov
Earth System Science Data, Volume 14, Issue 1

Abstract. Past efforts to synthesize and quantify the magnitude and change in carbon dioxide (CO2) fluxes in terrestrial ecosystems across the rapidly warming Arctic–boreal zone (ABZ) have provided valuable information but were limited in their geographical and temporal coverage. Furthermore, these efforts have been based on data aggregated over varying time periods, often with only minimal site ancillary data, thus limiting their potential to be used in large-scale carbon budget assessments. To bridge these gaps, we developed a standardized monthly database of Arctic–boreal CO2 fluxes (ABCflux) that aggregates in situ measurements of terrestrial net ecosystem CO2 exchange and its derived partitioned component fluxes: gross primary productivity and ecosystem respiration. The data span from 1989 to 2020 with over 70 supporting variables that describe key site conditions (e.g., vegetation and disturbance type), micrometeorological and environmental measurements (e.g., air and soil temperatures), and flux measurement techniques. Here, we describe these variables, the spatial and temporal distribution of observations, the main strengths and limitations of the database, and the potential research opportunities it enables. In total, ABCflux includes 244 sites and 6309 monthly observations; 136 sites and 2217 monthly observations represent tundra, and 108 sites and 4092 observations represent the boreal biome. The database includes fluxes estimated with chamber (19 % of the monthly observations), snow diffusion (3 %) and eddy covariance (78 %) techniques. The largest number of observations were collected during the climatological summer (June–August; 32 %), and fewer observations were available for autumn (September–October; 25 %), winter (December–February; 18 %), and spring (March–May; 25 %). ABCflux can be used in a wide array of empirical, remote sensing and modeling studies to improve understanding of the regional and temporal variability in CO2 fluxes and to better estimate the terrestrial ABZ CO2 budget. ABCflux is openly and freely available online (Virkkala et al., 2021b, https://doi.org/10.3334/ORNLDAAC/1934).

DOI bib abs
Causality guided machine learning model on wetland CH4 emissions across global wetlands
Kunxiaojia Yuan, Qing Zhu, Fa Li, W. J. Riley, Margaret Torn, Housen Chu, Gavin McNicol, Min Chen, Sara Knox, Kyle Delwiche, Huayi Wu, Dennis Baldocchi, Hongxu Ma, Ankur R. Desai, Jiquan Chen, Torsten Sachs, Masahito Ueyama, Oliver Sonnentag, Manuel Helbig, Eeva‐Stiina Tuittila, Gerald Jurasinski, Franziska Koebsch, David I. Campbell, Hans Peter Schmid, Annalea Lohila, Mathias Goeckede, Mats B. Nilsson, Thomas Friborg, Joachim Jansen, Donatella Zona, E. S. Euskirchen, Eric J. Ward, Gil Bohrer, Zhenong Jin, Licheng Liu, Hiroki Iwata, Jordan P. Goodrich, Robert B. Jackson
Agricultural and Forest Meteorology, Volume 324

Wetland CH4 emissions are among the most uncertain components of the global CH4 budget. The complex nature of wetland CH4 processes makes it challenging to identify causal relationships for improving our understanding and predictability of CH4 emissions. In this study, we used the flux measurements of CH4 from eddy covariance towers (30 sites from 4 wetlands types: bog, fen, marsh, and wet tundra) to construct a causality-constrained machine learning (ML) framework to explain the regulative factors and to capture CH4 emissions at sub-seasonal scale. We found that soil temperature is the dominant factor for CH4 emissions in all studied wetland types. Ecosystem respiration (CO2) and gross primary productivity exert controls at bog, fen, and marsh sites with lagged responses of days to weeks. Integrating these asynchronous environmental and biological causal relationships in predictive models significantly improved model performance. More importantly, modeled CH4 emissions differed by up to a factor of 4 under a +1°C warming scenario when causality constraints were considered. These results highlight the significant role of causality in modeling wetland CH4 emissions especially under future warming conditions, while traditional data-driven ML models may reproduce observations for the wrong reasons. Our proposed causality-guided model could benefit predictive modeling, large-scale upscaling, data gap-filling, and surrogate modeling of wetland CH4 emissions within earth system land models.

2021

DOI bib abs
The Boreal-Arctic Wetland and Lake Dataset (BAWLD)
David Olefeldt, Mikael Hovemyr, McKenzie A. Kuhn, David Bastviken, T. J. Bohn, John Connolly, P. M. Crill, E. S. Euskirchen, Sarah A. Finkelstein, Hélène Genet, Guido Grosse, Lorna I. Harris, Liam Heffernan, Manuel Helbig, Gustaf Hugelius, Ryan Hutchins, Sari Juutinen, Mark J. Lara, Avni Malhotra, Kristen Manies, A. David McGuire, Susan M. Natali, Jonathan A. O’Donnell, Frans‐Jan W. Parmentier, Aleksi Räsänen, Christina Schädel, Oliver Sonnentag, Maria Strack, Suzanne E. Tank, Claire C. Treat, R. K. Varner, Tarmo Virtanen, Rebecca K. Warren, Jennifer D. Watts, David Olefeldt, Mikael Hovemyr, McKenzie A. Kuhn, David Bastviken, T. J. Bohn, John Connolly, P. M. Crill, E. S. Euskirchen, Sarah A. Finkelstein, Hélène Genet, Guido Grosse, Lorna I. Harris, Liam Heffernan, Manuel Helbig, Gustaf Hugelius, Ryan Hutchins, Sari Juutinen, Mark J. Lara, Avni Malhotra, Kristen Manies, A. David McGuire, Susan M. Natali, Jonathan A. O’Donnell, Frans‐Jan W. Parmentier, Aleksi Räsänen, Christina Schädel, Oliver Sonnentag, Maria Strack, Suzanne E. Tank, Claire C. Treat, R. K. Varner, Tarmo Virtanen, Rebecca K. Warren, Jennifer D. Watts

Abstract. Methane emissions from boreal and arctic wetlands, lakes, and rivers are expected to increase in response to warming and associated permafrost thaw. However, the lack of appropriate land cover datasets for scaling field-measured methane emissions to circumpolar scales has contributed to a large uncertainty for our understanding of present-day and future methane emissions. Here we present the Boreal-Arctic Wetland and Lake Dataset (BAWLD), a land cover dataset based on an expert assessment, extrapolated using random forest modelling from available spatial datasets of climate, topography, soils, permafrost conditions, vegetation, wetlands, and surface water extents and dynamics. In BAWLD, we estimate the fractional coverage of five wetland, seven lake, and three river classes within 0.5 × 0.5° grid cells that cover the northern boreal and tundra biomes (17 % of the global land surface). Land cover classes were defined using criteria that ensured distinct methane emissions among classes, as indicated by a co-developed comprehensive dataset of methane flux observations. In BAWLD, wetlands occupied 3.2 × 106 km2 (14 % of domain) with a 95 % confidence interval between 2.8 and 3.8 × 106 km2. Bog, fen, and permafrost bog were the most abundant wetland classes, covering ~28 % each of the total wetland area, while the highest methane emitting marsh and tundra wetland classes occupied 5 and 12 %, respectively. Lakes, defined to include all lentic open-water ecosystems regardless of size, covered 1.4 × 106 km2 (6 % of domain). Low methane-emitting large lakes (> 10 km2) and glacial lakes jointly represented 78 % of the total lake area, while high-emitting peatland and yedoma lakes covered 18 and 4 %, respectively. Small (< 0.1 km2) glacial, peatland, and yedoma lakes combined covered 17 % of the total lake area, but contributed disproportionally to the overall spatial uncertainty of lake area with a 95 % confidence interval between 0.15 and 0.38 × 106 km2. Rivers and streams were estimated to cover 0.12 × 106 km2 (0.5 % of domain) of which 8 % was associated with high-methane emitting headwaters that drain organic-rich landscapes. Distinct combinations of spatially co-occurring wetland and lake classes were identified across the BAWLD domain, allowing for the mapping of “wetscapes” that will have characteristic methane emission magnitudes and sensitivities to climate change at regional scales. With BAWLD, we provide a dataset which avoids double-accounting of wetland, lake and river extents, and which includes confidence intervals for each land cover class. As such, BAWLD will be suitable for many hydrological and biogeochemical modelling and upscaling efforts for the northern Boreal and Arctic region, in particular those aimed at improving assessments of current and future methane emissions. Data is freely available at https://doi.org/10.18739/A2C824F9X (Olefeldt et al., 2021).

DOI bib abs
FLUXNET-CH4: a global, multi-ecosystem dataset and analysis of methane seasonality from freshwater wetlands
Kyle Delwiche, Sara Knox, Avni Malhotra, Etienne Fluet‐Chouinard, Gavin McNicol, Sarah Féron, Zutao Ouyang, Dario Papale, Carlo Trotta, Eleonora Canfora, You-Wei Cheah, Danielle Christianson, Ma. Carmelita Alberto, Pavel Alekseychik, Mika Aurela, Dennis Baldocchi, Sheel Bansal, David P. Billesbach, Gil Bohrer, Rosvel Bracho, Nina Buchmann, David I. Campbell, Gerardo Celis, Jiquan Chen, Weinan Chen, Housen Chu, Higo J. Dalmagro, Sigrid Dengel, Ankur R. Desai, Matteo Detto, A. J. Dolman, Elke Eichelmann, E. S. Euskirchen, D. Famulari, Kathrin Fuchs, Mathias Goeckede, Sébastien Gogo, Mangaliso J. Gondwe, Jordan P. Goodrich, Pia Gottschalk, Scott L. Graham, Martin Heimann, Manuel Helbig, Carole Helfter, Kyle S. Hemes, Takashi Hirano, David Y. Hollinger, Lukas Hörtnagl, Hiroki Iwata, Adrien Jacotot, Gerald Jurasinski, Minseok Kang, Kuno Kasak, John S. King, Janina Klatt, Franziska Koebsch, Ken W. Krauss, Derrick Y.F. Lai, Annalea Lohila, Ivan Mammarella, Luca Belelli Marchesini, Giovanni Manca, Jaclyn Hatala Matthes, Trofim C. Maximov, Lutz Merbold, Bhaskar Mitra, Timothy H. Morin, Eiko Nemitz, Mats B. Nilsson, Shuli Niu, Walter C. Oechel, Patricia Y. Oikawa, Keisuke Ono, Matthias Peichl, Olli Peltola, Michele L. Reba, Andrew D. Richardson, W. J. Riley, Benjamin R. K. Runkle, Youngryel Ryu, Torsten Sachs, Ayaka Sakabe, Camilo Rey‐Sánchez, Edward A. G. Schuur, Karina V. R. Schäfer, Oliver Sonnentag, Jed P. Sparks, Ellen Stuart‐Haëntjens, Cove Sturtevant, Ryan C. Sullivan, Daphne Szutu, Jonathan E. Thom, Margaret Torn, Eeva‐Stiina Tuittila, Jessica Turner, Masahito Ueyama, Alex Valach, Rodrigo Vargas, Andrej Varlagin, Alma Vázquez‐Lule, Joseph Verfaillie, Timo Vesala, George L. Vourlitis, Eric J. Ward, Christian Wille, Georg Wohlfahrt, Guan Xhuan Wong, Zhen Zhang, Donatella Zona, L. Windham‐Myers, Benjamin Poulter, Robert B. Jackson, Kyle Delwiche, Sara Knox, Avni Malhotra, Etienne Fluet‐Chouinard, Gavin McNicol, Sarah Féron, Zutao Ouyang, Dario Papale, Carlo Trotta, Eleonora Canfora, You-Wei Cheah, Danielle Christianson, Ma. Carmelita Alberto, Pavel Alekseychik, Mika Aurela, Dennis Baldocchi, Sheel Bansal, David P. Billesbach, Gil Bohrer, Rosvel Bracho, Nina Buchmann, David I. Campbell, Gerardo Celis, Jiquan Chen, Weinan Chen, Housen Chu, Higo J. Dalmagro, Sigrid Dengel, Ankur R. Desai, Matteo Detto, A. J. Dolman, Elke Eichelmann, E. S. Euskirchen, D. Famulari, Kathrin Fuchs, Mathias Goeckede, Sébastien Gogo, Mangaliso J. Gondwe, Jordan P. Goodrich, Pia Gottschalk, Scott L. Graham, Martin Heimann, Manuel Helbig, Carole Helfter, Kyle S. Hemes, Takashi Hirano, David Y. Hollinger, Lukas Hörtnagl, Hiroki Iwata, Adrien Jacotot, Gerald Jurasinski, Minseok Kang, Kuno Kasak, John S. King, Janina Klatt, Franziska Koebsch, Ken W. Krauss, Derrick Y.F. Lai, Annalea Lohila, Ivan Mammarella, Luca Belelli Marchesini, Giovanni Manca, Jaclyn Hatala Matthes, Trofim C. Maximov, Lutz Merbold, Bhaskar Mitra, Timothy H. Morin, Eiko Nemitz, Mats B. Nilsson, Shuli Niu, Walter C. Oechel, Patricia Y. Oikawa, Keisuke Ono, Matthias Peichl, Olli Peltola, Michele L. Reba, Andrew D. Richardson, W. J. Riley, Benjamin R. K. Runkle, Youngryel Ryu, Torsten Sachs, Ayaka Sakabe, Camilo Rey‐Sánchez, Edward A. G. Schuur, Karina V. R. Schäfer, Oliver Sonnentag, Jed P. Sparks, Ellen Stuart‐Haëntjens, Cove Sturtevant, Ryan C. Sullivan, Daphne Szutu, Jonathan E. Thom, Margaret Torn, Eeva‐Stiina Tuittila, Jessica Turner, Masahito Ueyama, Alex Valach, Rodrigo Vargas, Andrej Varlagin, Alma Vázquez‐Lule, Joseph Verfaillie, Timo Vesala, George L. Vourlitis, Eric J. Ward, Christian Wille, Georg Wohlfahrt, Guan Xhuan Wong, Zhen Zhang, Donatella Zona, L. Windham‐Myers, Benjamin Poulter, Robert B. Jackson
Earth System Science Data, Volume 13, Issue 7

Abstract. Methane (CH4) emissions from natural landscapes constitute roughly half of global CH4 contributions to the atmosphere, yet large uncertainties remain in the absolute magnitude and the seasonality of emission quantities and drivers. Eddy covariance (EC) measurements of CH4 flux are ideal for constraining ecosystem-scale CH4 emissions due to quasi-continuous and high-temporal-resolution CH4 flux measurements, coincident carbon dioxide, water, and energy flux measurements, lack of ecosystem disturbance, and increased availability of datasets over the last decade. Here, we (1) describe the newly published dataset, FLUXNET-CH4 Version 1.0, the first open-source global dataset of CH4 EC measurements (available at https://fluxnet.org/data/fluxnet-ch4-community-product/, last access: 7 April 2021). FLUXNET-CH4 includes half-hourly and daily gap-filled and non-gap-filled aggregated CH4 fluxes and meteorological data from 79 sites globally: 42 freshwater wetlands, 6 brackish and saline wetlands, 7 formerly drained ecosystems, 7 rice paddy sites, 2 lakes, and 15 uplands. Then, we (2) evaluate FLUXNET-CH4 representativeness for freshwater wetland coverage globally because the majority of sites in FLUXNET-CH4 Version 1.0 are freshwater wetlands which are a substantial source of total atmospheric CH4 emissions; and (3) we provide the first global estimates of the seasonal variability and seasonality predictors of freshwater wetland CH4 fluxes. Our representativeness analysis suggests that the freshwater wetland sites in the dataset cover global wetland bioclimatic attributes (encompassing energy, moisture, and vegetation-related parameters) in arctic, boreal, and temperate regions but only sparsely cover humid tropical regions. Seasonality metrics of wetland CH4 emissions vary considerably across latitudinal bands. In freshwater wetlands (except those between 20∘ S to 20∘ N) the spring onset of elevated CH4 emissions starts 3 d earlier, and the CH4 emission season lasts 4 d longer, for each degree Celsius increase in mean annual air temperature. On average, the spring onset of increasing CH4 emissions lags behind soil warming by 1 month, with very few sites experiencing increased CH4 emissions prior to the onset of soil warming. In contrast, roughly half of these sites experience the spring onset of rising CH4 emissions prior to the spring increase in gross primary productivity (GPP). The timing of peak summer CH4 emissions does not correlate with the timing for either peak summer temperature or peak GPP. Our results provide seasonality parameters for CH4 modeling and highlight seasonality metrics that cannot be predicted by temperature or GPP (i.e., seasonality of CH4 peak). FLUXNET-CH4 is a powerful new resource for diagnosing and understanding the role of terrestrial ecosystems and climate drivers in the global CH4 cycle, and future additions of sites in tropical ecosystems and site years of data collection will provide added value to this database. All seasonality parameters are available at https://doi.org/10.5281/zenodo.4672601 (Delwiche et al., 2021). Additionally, raw FLUXNET-CH4 data used to extract seasonality parameters can be downloaded from https://fluxnet.org/data/fluxnet-ch4-community-product/ (last access: 7 April 2021), and a complete list of the 79 individual site data DOIs is provided in Table 2 of this paper.

DOI bib abs
Substantial hysteresis in emergent temperature sensitivity of global wetland CH4 emissions
Kuang‐Yu Chang, W. J. Riley, Sara Knox, Robert B. Jackson, Gavin McNicol, Benjamin Poulter, Mika Aurela, Dennis Baldocchi, Sheel Bansal, Gil Bohrer, David I. Campbell, Alessandro Cescatti, Housen Chu, Kyle Delwiche, Ankur R. Desai, E. S. Euskirchen, Thomas Friborg, Mathias Goeckede, Manuel Helbig, Kyle S. Hemes, Takashi Hirano, Hiroki Iwata, Minseok Kang, Trevor F. Keenan, Ken W. Krauss, Annalea Lohila, Ivan Mammarella, Bhaskar Mitra, Akira Miyata, Mats B. Nilsson, Asko Noormets, Walter C. Oechel, Dario Papale, Matthias Peichl, Michele L. Reba, Janne Rinne, Benjamin R. K. Runkle, Youngryel Ryu, Torsten Sachs, Karina V. R. Schäfer, Hans Peter Schmid, Narasinha Shurpali, Oliver Sonnentag, Angela C. I. Tang, Margaret Torn, Carlo Trotta, Eeva‐Stiina Tuittila, Masahito Ueyama, Rodrigo Vargas, Timo Vesala, L. Windham‐Myers, Zhen Zhang, Donatella Zona, Kuang‐Yu Chang, W. J. Riley, Sara Knox, Robert B. Jackson, Gavin McNicol, Benjamin Poulter, Mika Aurela, Dennis Baldocchi, Sheel Bansal, Gil Bohrer, David I. Campbell, Alessandro Cescatti, Housen Chu, Kyle Delwiche, Ankur R. Desai, E. S. Euskirchen, Thomas Friborg, Mathias Goeckede, Manuel Helbig, Kyle S. Hemes, Takashi Hirano, Hiroki Iwata, Minseok Kang, Trevor F. Keenan, Ken W. Krauss, Annalea Lohila, Ivan Mammarella, Bhaskar Mitra, Akira Miyata, Mats B. Nilsson, Asko Noormets, Walter C. Oechel, Dario Papale, Matthias Peichl, Michele L. Reba, Janne Rinne, Benjamin R. K. Runkle, Youngryel Ryu, Torsten Sachs, Karina V. R. Schäfer, Hans Peter Schmid, Narasinha Shurpali, Oliver Sonnentag, Angela C. I. Tang, Margaret Torn, Carlo Trotta, Eeva‐Stiina Tuittila, Masahito Ueyama, Rodrigo Vargas, Timo Vesala, L. Windham‐Myers, Zhen Zhang, Donatella Zona
Nature Communications, Volume 12, Issue 1

Abstract Wetland methane (CH 4 ) emissions ( $${F}_{{{CH}}_{4}}$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:msub> <mml:mrow> <mml:mi>F</mml:mi> </mml:mrow> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi>C</mml:mi> <mml:mi>H</mml:mi> </mml:mrow> <mml:mrow> <mml:mn>4</mml:mn> </mml:mrow> </mml:msub> </mml:mrow> </mml:msub> </mml:math> ) are important in global carbon budgets and climate change assessments. Currently, $${F}_{{{CH}}_{4}}$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:msub> <mml:mrow> <mml:mi>F</mml:mi> </mml:mrow> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi>C</mml:mi> <mml:mi>H</mml:mi> </mml:mrow> <mml:mrow> <mml:mn>4</mml:mn> </mml:mrow> </mml:msub> </mml:mrow> </mml:msub> </mml:math> projections rely on prescribed static temperature sensitivity that varies among biogeochemical models. Meta-analyses have proposed a consistent $${F}_{{{CH}}_{4}}$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:msub> <mml:mrow> <mml:mi>F</mml:mi> </mml:mrow> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi>C</mml:mi> <mml:mi>H</mml:mi> </mml:mrow> <mml:mrow> <mml:mn>4</mml:mn> </mml:mrow> </mml:msub> </mml:mrow> </mml:msub> </mml:math> temperature dependence across spatial scales for use in models; however, site-level studies demonstrate that $${F}_{{{CH}}_{4}}$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:msub> <mml:mrow> <mml:mi>F</mml:mi> </mml:mrow> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi>C</mml:mi> <mml:mi>H</mml:mi> </mml:mrow> <mml:mrow> <mml:mn>4</mml:mn> </mml:mrow> </mml:msub> </mml:mrow> </mml:msub> </mml:math> are often controlled by factors beyond temperature. Here, we evaluate the relationship between $${F}_{{{CH}}_{4}}$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:msub> <mml:mrow> <mml:mi>F</mml:mi> </mml:mrow> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi>C</mml:mi> <mml:mi>H</mml:mi> </mml:mrow> <mml:mrow> <mml:mn>4</mml:mn> </mml:mrow> </mml:msub> </mml:mrow> </mml:msub> </mml:math> and temperature using observations from the FLUXNET-CH 4 database. Measurements collected across the globe show substantial seasonal hysteresis between $${F}_{{{CH}}_{4}}$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:msub> <mml:mrow> <mml:mi>F</mml:mi> </mml:mrow> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi>C</mml:mi> <mml:mi>H</mml:mi> </mml:mrow> <mml:mrow> <mml:mn>4</mml:mn> </mml:mrow> </mml:msub> </mml:mrow> </mml:msub> </mml:math> and temperature, suggesting larger $${F}_{{{CH}}_{4}}$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:msub> <mml:mrow> <mml:mi>F</mml:mi> </mml:mrow> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi>C</mml:mi> <mml:mi>H</mml:mi> </mml:mrow> <mml:mrow> <mml:mn>4</mml:mn> </mml:mrow> </mml:msub> </mml:mrow> </mml:msub> </mml:math> sensitivity to temperature later in the frost-free season (about 77% of site-years). Results derived from a machine-learning model and several regression models highlight the importance of representing the large spatial and temporal variability within site-years and ecosystem types. Mechanistic advancements in biogeochemical model parameterization and detailed measurements in factors modulating CH 4 production are thus needed to improve global CH 4 budget assessments.

DOI bib abs
The implications of permafrost thaw and land cover change on snow water equivalent accumulation, melt and runoff in discontinuous permafrost peatlands
Ryan F. Connon, L. Chasmer, Emily Haughton, Manuel Helbig, Chris Hopkinson, Oliver Sonnentag, W. L. Quinton, Ryan F. Connon, L. Chasmer, Emily Haughton, Manuel Helbig, Chris Hopkinson, Oliver Sonnentag, W. L. Quinton
Hydrological Processes, Volume 35, Issue 9

In the discontinuous permafrost zone of the Northwest Territories (NWT), Canada, snow covers the ground surface for half the year. Snowmelt constitutes a primary source of moisture supply for the short growing season and strongly influences stream hydrographs. Permafrost thaw has changed the landscape by increasing the proportional coverage of permafrost-free wetlands at the expense of permafrost-cored peat plateau forests. The biophysical characteristics of each feature affect snow water equivalent (SWE) accumulation and melt rates. In headwater streams in the southern Dehcho region of the NWT, snowmelt runoff has significantly increased over the past 50 years, despite no significant change in annual SWE. At the Fort Simpson A climate station, we found that SWE measurements made by Environment and Climate Change Canada using a Nipher precipitation gauge were more accurate than the Adjusted and Homogenized Canadian Climate Dataset which was derived from snow depth measurements. Here, we: (a) provide 13 years of snow survey data to demonstrate differences in end-of-season SWE between wetlands and plateau forests; (b) provide ablation stake and radiation measurements to document differences in snow melt patterns among wetlands, plateau forests, and upland forests; and (c) evaluate the potential impact of permafrost-thaw induced wetland expansion on SWE accumulation, melt, and runoff. We found that plateaus retain significantly (p < 0.01) more SWE than wetlands. However, the differences are too small (123 mm and 111 mm, respectively) to cause any substantial change in basin SWE. During the snowmelt period in 2015, wetlands were the first feature to become snow-free in mid-April, followed by plateau forests (7 days after wetlands) and upland forests (18 days after wetlands). A transition to a higher percentage cover of wetlands may lead to more rapid snowmelt and provide a more hydrologically-connected landscape, a plausible mechanism driving the observed increase in spring freshet runoff.

DOI bib abs
Gap-filling eddy covariance methane fluxes: Comparison of machine learning model predictions and uncertainties at FLUXNET-CH4 wetlands
Jeremy Irvin, Sharon Zhou, Gavin McNicol, Fred Lu, Vincent Liu, Etienne Fluet‐Chouinard, Zutao Ouyang, Sara Knox, Antje Lucas-Moffat, Carlo Trotta, Dario Papale, Domenico Vitale, Ivan Mammarella, Pavel Alekseychik, Mika Aurela, Anand Avati, Dennis Baldocchi, Sheel Bansal, Gil Bohrer, David I. Campbell, Jiquan Chen, Housen Chu, Higo J. Dalmagro, Kyle Delwiche, Ankur R. Desai, E. S. Euskirchen, Sarah Féron, Mathias Goeckede, Martin Heimann, Manuel Helbig, Carole Helfter, Kyle S. Hemes, Takashi Hirano, Hiroki Iwata, Gerald Jurasinski, Aram Kalhori, Andrew Kondrich, Derrick Y.F. Lai, Annalea Lohila, Avni Malhotra, Lutz Merbold, Bhaskar Mitra, Andrew Y. Ng, Mats B. Nilsson, Asko Noormets, Matthias Peichl, Camilo Rey‐Sánchez, Andrew D. Richardson, Benjamin R. K. Runkle, Karina VR Schäfer, Oliver Sonnentag, Ellen Stuart‐Haëntjens, Cove Sturtevant, Masahito Ueyama, Alex Valach, Rodrigo Vargas, George L. Vourlitis, Eric J. Ward, Guan Xhuan Wong, Donatella Zona, Ma. Carmelita Alberto, David P. Billesbach, Gerardo Celis, A. J. Dolman, Thomas Friborg, Kathrin Fuchs, Sébastien Gogo, Mangaliso J. Gondwe, Jordan P. Goodrich, Pia Gottschalk, Lukas Hörtnagl, Adrien Jacotot, Franziska Koebsch, Kuno Kasak, Regine Maier, Timothy H. Morin, Eiko Nemitz, Walter C. Oechel, Patricia Y. Oikawa, Keisuke Ono, Torsten Sachs, Ayaka Sakabe, Edward A. G. Schuur, Robert Shortt, Ryan C. Sullivan, Daphne Szutu, Eeva‐Stiina Tuittila, Andrej Varlagin, Joeseph G Verfaillie, Christian Wille, L. Windham‐Myers, Benjamin Poulter, Robert B. Jackson, Jeremy Irvin, Sharon Zhou, Gavin McNicol, Fred Lu, Vincent Liu, Etienne Fluet‐Chouinard, Zutao Ouyang, Sara Knox, Antje Lucas-Moffat, Carlo Trotta, Dario Papale, Domenico Vitale, Ivan Mammarella, Pavel Alekseychik, Mika Aurela, Anand Avati, Dennis Baldocchi, Sheel Bansal, Gil Bohrer, David I. Campbell, Jiquan Chen, Housen Chu, Higo J. Dalmagro, Kyle Delwiche, Ankur R. Desai, E. S. Euskirchen, Sarah Féron, Mathias Goeckede, Martin Heimann, Manuel Helbig, Carole Helfter, Kyle S. Hemes, Takashi Hirano, Hiroki Iwata, Gerald Jurasinski, Aram Kalhori, Andrew Kondrich, Derrick Y.F. Lai, Annalea Lohila, Avni Malhotra, Lutz Merbold, Bhaskar Mitra, Andrew Y. Ng, Mats B. Nilsson, Asko Noormets, Matthias Peichl, Camilo Rey‐Sánchez, Andrew D. Richardson, Benjamin R. K. Runkle, Karina VR Schäfer, Oliver Sonnentag, Ellen Stuart‐Haëntjens, Cove Sturtevant, Masahito Ueyama, Alex Valach, Rodrigo Vargas, George L. Vourlitis, Eric J. Ward, Guan Xhuan Wong, Donatella Zona, Ma. Carmelita Alberto, David P. Billesbach, Gerardo Celis, A. J. Dolman, Thomas Friborg, Kathrin Fuchs, Sébastien Gogo, Mangaliso J. Gondwe, Jordan P. Goodrich, Pia Gottschalk, Lukas Hörtnagl, Adrien Jacotot, Franziska Koebsch, Kuno Kasak, Regine Maier, Timothy H. Morin, Eiko Nemitz, Walter C. Oechel, Patricia Y. Oikawa, Keisuke Ono, Torsten Sachs, Ayaka Sakabe, Edward A. G. Schuur, Robert Shortt, Ryan C. Sullivan, Daphne Szutu, Eeva‐Stiina Tuittila, Andrej Varlagin, Joeseph G Verfaillie, Christian Wille, L. Windham‐Myers, Benjamin Poulter, Robert B. Jackson
Agricultural and Forest Meteorology, Volume 308-309

• We evaluate methane flux gap-filling methods across 17 boreal-to-tropical wetlands • New methods for generating realistic artificial gaps and uncertainties are proposed • Decision tree algorithms perform slightly better than neural networks on average • Soil temperature and generic seasonality are the most important predictors • Open-source code is released for gap-filling steps and uncertainty evaluation Time series of wetland methane fluxes measured by eddy covariance require gap-filling to estimate daily, seasonal, and annual emissions. Gap-filling methane fluxes is challenging because of high variability and complex responses to multiple drivers. To date, there is no widely established gap-filling standard for wetland methane fluxes, with regards both to the best model algorithms and predictors. This study synthesizes results of different gap-filling methods systematically applied at 17 wetland sites spanning boreal to tropical regions and including all major wetland classes and two rice paddies. Procedures are proposed for: 1) creating realistic artificial gap scenarios, 2) training and evaluating gap-filling models without overstating performance, and 3) predicting half-hourly methane fluxes and annual emissions with realistic uncertainty estimates. Performance is compared between a conventional method (marginal distribution sampling) and four machine learning algorithms. The conventional method achieved similar median performance as the machine learning models but was worse than the best machine learning models and relatively insensitive to predictor choices. Of the machine learning models, decision tree algorithms performed the best in cross-validation experiments, even with a baseline predictor set, and artificial neural networks showed comparable performance when using all predictors. Soil temperature was frequently the most important predictor whilst water table depth was important at sites with substantial water table fluctuations, highlighting the value of data on wetland soil conditions. Raw gap-filling uncertainties from the machine learning models were underestimated and we propose a method to calibrate uncertainties to observations. The python code for model development, evaluation, and uncertainty estimation is publicly available. This study outlines a modular and robust machine learning workflow and makes recommendations for, and evaluates an improved baseline of, methane gap-filling models that can be implemented in multi-site syntheses or standardized products from regional and global flux networks (e.g., FLUXNET).

DOI bib abs
Soil respiration strongly offsets carbon uptake in Alaska and Northwest Canada
Jennifer D. Watts, Susan M. Natali, Christina Minions, D. A. Risk, Kyle A. Arndt, Donatella Zona, E. S. Euskirchen, Adrian V. Rocha, Oliver Sonnentag, Manuel Helbig, Aram Kalhori, Walter C. Oechel, Hiroki Ikawa, Masahito Ueyama, Rikie Suzuki, Hideki Kobayashi, Gerardo Celis, Edward A. G. Schuur, Elyn Humphreys, Yongwon Kim, Bang‐Yong Lee, S. J. Goetz, Nima Madani, Luke D. Schiferl, R. Commane, John S. Kimball, Zhihua Liu, Margaret Torn, Stefano Potter, Jonathan Wang, M. Torre Jorgenson, Jingfeng Xiao, Xing Li, Colin W. Edgar, Jennifer D. Watts, Susan M. Natali, Christina Minions, D. A. Risk, Kyle A. Arndt, Donatella Zona, E. S. Euskirchen, Adrian V. Rocha, Oliver Sonnentag, Manuel Helbig, Aram Kalhori, Walter C. Oechel, Hiroki Ikawa, Masahito Ueyama, Rikie Suzuki, Hideki Kobayashi, Gerardo Celis, Edward A. G. Schuur, Elyn Humphreys, Yongwon Kim, Bang‐Yong Lee, S. J. Goetz, Nima Madani, Luke D. Schiferl, R. Commane, John S. Kimball, Zhihua Liu, Margaret Torn, Stefano Potter, Jonathan Wang, M. Torre Jorgenson, Jingfeng Xiao, Xing Li, Colin W. Edgar
Environmental Research Letters, Volume 16, Issue 8

Abstract Soil respiration (i.e. from soils and roots) provides one of the largest global fluxes of carbon dioxide (CO 2 ) to the atmosphere and is likely to increase with warming, yet the magnitude of soil respiration from rapidly thawing Arctic-boreal regions is not well understood. To address this knowledge gap, we first compiled a new CO 2 flux database for permafrost-affected tundra and boreal ecosystems in Alaska and Northwest Canada. We then used the CO 2 database, multi-sensor satellite imagery, and random forest models to assess the regional magnitude of soil respiration. The flux database includes a new Soil Respiration Station network of chamber-based fluxes, and fluxes from eddy covariance towers. Our site-level data, spanning September 2016 to August 2017, revealed that the largest soil respiration emissions occurred during the summer (June–August) and that summer fluxes were higher in boreal sites (1.87 ± 0.67 g CO 2 –C m −2 d −1 ) relative to tundra (0.94 ± 0.4 g CO 2 –C m −2 d −1 ). We also observed considerable emissions (boreal: 0.24 ± 0.2 g CO 2 –C m −2 d −1 ; tundra: 0.18 ± 0.16 g CO 2 –C m −2 d −1 ) from soils during the winter (November–March) despite frozen surface conditions. Our model estimates indicated an annual region-wide loss from soil respiration of 591 ± 120 Tg CO 2 –C during the 2016–2017 period. Summer months contributed to 58% of the regional soil respiration, winter months contributed to 15%, and the shoulder months contributed to 27%. In total, soil respiration offset 54% of annual gross primary productivity (GPP) across the study domain. We also found that in tundra environments, transitional tundra/boreal ecotones, and in landscapes recently affected by fire, soil respiration often exceeded GPP, resulting in a net annual source of CO 2 to the atmosphere. As this region continues to warm, soil respiration may increasingly offset GPP, further amplifying global climate change.

DOI bib abs
Representativeness of Eddy-Covariance flux footprints for areas surrounding AmeriFlux sites
Housen Chu, Xiangzhong Luo, Zutao Ouyang, Stephen Chan, Sigrid Dengel, Sébastien Biraud, Margaret Torn, Stefan Metzger, Jitendra Kumar, M. Altaf Arain, T. J. Arkebauer, Dennis Baldocchi, Carl J. Bernacchi, D. P. Billesbach, T. Andrew Black, Peter D. Blanken, Gil Bohrer, Rosvel Bracho, S. M. Brown, N. A. Brunsell, Jiquan Chen, Xingyuan Chen, Kenneth L. Clark, Ankur R. Desai, Tomer Duman, David Durden, Silvano Fares, Inke Forbrich, John A. Gamon, Christopher M. Gough, Timothy J. Griffis, Manuel Helbig, David Y. Hollinger, Elyn Humphreys, Hiroki Ikawa, Hiroki Iwata, Yang Ju, John F. Knowles, Sara Knox, Hideki Kobayashi, Thomas E. Kolb, B. E. Law, Xuhui Lee, M. E. Litvak, Heping Liu, J. William Munger, Asko Noormets, Kim Novick, Steven F. Oberbauer, Walter C. Oechel, Patricia Y. Oikawa, S. A. Papuga, Elise Pendall, Prajaya Prajapati, John H. Prueger, W. L. Quinton, Andrew D. Richardson, Eric S. Russell, Russell L. Scott, Gregory Starr, R. M. Staebler, Paul C. Stoy, Ellen Stuart‐Haëntjens, Oliver Sonnentag, Ryan C. Sullivan, Andy Suyker, Masahito Ueyama, Rodrigo Vargas, Jeffrey D. Wood, Donatella Zona, Housen Chu, Xiangzhong Luo, Zutao Ouyang, Stephen Chan, Sigrid Dengel, Sébastien Biraud, Margaret Torn, Stefan Metzger, Jitendra Kumar, M. Altaf Arain, T. J. Arkebauer, Dennis Baldocchi, Carl J. Bernacchi, D. P. Billesbach, T. Andrew Black, Peter D. Blanken, Gil Bohrer, Rosvel Bracho, S. M. Brown, N. A. Brunsell, Jiquan Chen, Xingyuan Chen, Kenneth L. Clark, Ankur R. Desai, Tomer Duman, David Durden, Silvano Fares, Inke Forbrich, John A. Gamon, Christopher M. Gough, Timothy J. Griffis, Manuel Helbig, David Y. Hollinger, Elyn Humphreys, Hiroki Ikawa, Hiroki Iwata, Yang Ju, John F. Knowles, Sara Knox, Hideki Kobayashi, Thomas E. Kolb, B. E. Law, Xuhui Lee, M. E. Litvak, Heping Liu, J. William Munger, Asko Noormets, Kim Novick, Steven F. Oberbauer, Walter C. Oechel, Patricia Y. Oikawa, S. A. Papuga, Elise Pendall, Prajaya Prajapati, John H. Prueger, W. L. Quinton, Andrew D. Richardson, Eric S. Russell, Russell L. Scott, Gregory Starr, R. M. Staebler, Paul C. Stoy, Ellen Stuart‐Haëntjens, Oliver Sonnentag, Ryan C. Sullivan, Andy Suyker, Masahito Ueyama, Rodrigo Vargas, Jeffrey D. Wood, Donatella Zona
Agricultural and Forest Meteorology, Volume 301-302

• Large-scale eddy-covariance flux datasets need to be used with footprint-awareness • Using a fixed-extent target area across sites can bias model-data integration • Most sites do not represent the dominant land-cover type at a larger spatial extent • A representativeness index provides general guidance for site selection and data use Large datasets of greenhouse gas and energy surface-atmosphere fluxes measured with the eddy-covariance technique (e.g., FLUXNET2015, AmeriFlux BASE) are widely used to benchmark models and remote-sensing products. This study addresses one of the major challenges facing model-data integration: To what spatial extent do flux measurements taken at individual eddy-covariance sites reflect model- or satellite-based grid cells? We evaluate flux footprints—the temporally dynamic source areas that contribute to measured fluxes—and the representativeness of these footprints for target areas (e.g., within 250–3000 m radii around flux towers) that are often used in flux-data synthesis and modeling studies. We examine the land-cover composition and vegetation characteristics, represented here by the Enhanced Vegetation Index (EVI), in the flux footprints and target areas across 214 AmeriFlux sites, and evaluate potential biases as a consequence of the footprint-to-target-area mismatch. Monthly 80% footprint climatologies vary across sites and through time ranging four orders of magnitude from 10 3 to 10 7 m 2 due to the measurement heights, underlying vegetation- and ground-surface characteristics, wind directions, and turbulent state of the atmosphere. Few eddy-covariance sites are located in a truly homogeneous landscape. Thus, the common model-data integration approaches that use a fixed-extent target area across sites introduce biases on the order of 4%–20% for EVI and 6%–20% for the dominant land cover percentage. These biases are site-specific functions of measurement heights, target area extents, and land-surface characteristics. We advocate that flux datasets need to be used with footprint awareness, especially in research and applications that benchmark against models and data products with explicit spatial information. We propose a simple representativeness index based on our evaluations that can be used as a guide to identify site-periods suitable for specific applications and to provide general guidance for data use.

DOI bib abs
The Boreal–Arctic Wetland and Lake Dataset (BAWLD)
David Olefeldt, Mikael Hovemyr, McKenzie A. Kuhn, David Bastviken, T. J. Bohn, John Connolly, P. M. Crill, E. S. Euskirchen, Sarah A. Finkelstein, Hélène Genet, Guido Grosse, Lorna I. Harris, Liam Heffernan, Manuel Helbig, Gustaf Hugelius, Ryan Hutchins, Sari Juutinen, Mark J. Lara, Avni Malhotra, Kristen Manies, A. David McGuire, Susan M. Natali, Jonathan A. O’Donnell, Frans‐Jan W. Parmentier, Aleksi Räsänen, Christina Schädel, Oliver Sonnentag, Maria Strack, Suzanne E. Tank, Claire C. Treat, R. K. Varner, Tarmo Virtanen, Rebecca K. Warren, Jennifer D. Watts
Earth System Science Data, Volume 13, Issue 11

Abstract. Methane emissions from boreal and arctic wetlands, lakes, and rivers are expected to increase in response to warming and associated permafrost thaw. However, the lack of appropriate land cover datasets for scaling field-measured methane emissions to circumpolar scales has contributed to a large uncertainty for our understanding of present-day and future methane emissions. Here we present the Boreal–Arctic Wetland and Lake Dataset (BAWLD), a land cover dataset based on an expert assessment, extrapolated using random forest modelling from available spatial datasets of climate, topography, soils, permafrost conditions, vegetation, wetlands, and surface water extents and dynamics. In BAWLD, we estimate the fractional coverage of five wetland, seven lake, and three river classes within 0.5 × 0.5∘ grid cells that cover the northern boreal and tundra biomes (17 % of the global land surface). Land cover classes were defined using criteria that ensured distinct methane emissions among classes, as indicated by a co-developed comprehensive dataset of methane flux observations. In BAWLD, wetlands occupied 3.2 × 106 km2 (14 % of domain) with a 95 % confidence interval between 2.8 and 3.8 × 106 km2. Bog, fen, and permafrost bog were the most abundant wetland classes, covering ∼ 28 % each of the total wetland area, while the highest-methane-emitting marsh and tundra wetland classes occupied 5 % and 12 %, respectively. Lakes, defined to include all lentic open-water ecosystems regardless of size, covered 1.4 × 106 km2 (6 % of domain). Low-methane-emitting large lakes (>10 km2) and glacial lakes jointly represented 78 % of the total lake area, while high-emitting peatland and yedoma lakes covered 18 % and 4 %, respectively. Small (<0.1 km2) glacial, peatland, and yedoma lakes combined covered 17 % of the total lake area but contributed disproportionally to the overall spatial uncertainty in lake area with a 95 % confidence interval between 0.15 and 0.38 × 106 km2. Rivers and streams were estimated to cover 0.12 × 106 km2 (0.5 % of domain), of which 8 % was associated with high-methane-emitting headwaters that drain organic-rich landscapes. Distinct combinations of spatially co-occurring wetland and lake classes were identified across the BAWLD domain, allowing for the mapping of “wetscapes” that have characteristic methane emission magnitudes and sensitivities to climate change at regional scales. With BAWLD, we provide a dataset which avoids double-accounting of wetland, lake, and river extents and which includes confidence intervals for each land cover class. As such, BAWLD will be suitable for many hydrological and biogeochemical modelling and upscaling efforts for the northern boreal and arctic region, in particular those aimed at improving assessments of current and future methane emissions. Data are freely available at https://doi.org/10.18739/A2C824F9X (Olefeldt et al., 2021).

2020

DOI bib abs
Increasing contribution of peatlands to boreal evapotranspiration in a warming climate
Manuel Helbig, J. M. Waddington, Pavel Alekseychik, B. D. Amiro, Mika Aurela, Alan Barr, T. Andrew Black, Peter D. Blanken, Sean K. Carey, Jiquan Chen, Jinshu Chi, Ankur R. Desai, Allison L. Dunn, E. S. Euskirchen, Lawrence B. Flanagan, Inke Forbrich, Thomas Friborg, Achim Grelle, Silvie Harder, Michal Heliasz, Elyn Humphreys, Hiroki Ikawa, Pierre‐Erik Isabelle, Hiroki Iwata, Rachhpal S. Jassal, Mika Korkiakoski, J. Kurbatova, Lars Kutzbach, Anders Lindroth, Mikaell Ottosson Löfvenius, Annalea Lohila, Ivan Mammarella, Philip Marsh, Trofim C. Maximov, Joe R. Melton, Paul Moore, Daniel F. Nadeau, Erin M. Nicholls, Mats B. Nilsson, Takeshi Ohta, Matthias Peichl, Richard M. Petrone, Roman Petrov, Anatoly Prokushkin, W. L. Quinton, David E. Reed, Nigel T. Roulet, Benjamin R. K. Runkle, Oliver Sonnentag, Ian B. Strachan, Pierre Taillardat, Eeva‐Stiina Tuittila, Juha‐Pekka Tuovinen, Jessica Turner, Masahito Ueyama, Andrej Varlagin, Martin Wilmking, Steven C. Wofsy, Vyacheslav Zyrianov
Nature Climate Change, Volume 10, Issue 6

The response of evapotranspiration (ET) to warming is of critical importance to the water and carbon cycle of the boreal biome, a mosaic of land cover types dominated by forests and peatlands. The effect of warming-induced vapour pressure deficit (VPD) increases on boreal ET remains poorly understood because peatlands are not specifically represented as plant functional types in Earth system models. Here we show that peatland ET increases more than forest ET with increasing VPD using observations from 95 eddy covariance tower sites. At high VPD of more than 2 kPa, peatland ET exceeds forest ET by up to 30%. Future (2091–2100) mid-growing season peatland ET is estimated to exceed forest ET by over 20% in about one-third of the boreal biome for RCP4.5 and about two-thirds for RCP8.5. Peatland-specific ET responses to VPD should therefore be included in Earth system models to avoid biases in water and carbon cycle projections.

DOI bib abs
The biophysical climate mitigation potential of boreal peatlands during the growing season
Manuel Helbig, J. M. Waddington, Pavel Alekseychik, B. D. Amiro, Mika Aurela, Alan Barr, T. Andrew Black, Sean K. Carey, Jiquan Chen, Jinshu Chi, Ankur R. Desai, Allison L. Dunn, E. S. Euskirchen, Lawrence B. Flanagan, Thomas Friborg, Michelle Garneau, Achim Grelle, Silvie Harder, Michal Heliasz, Elyn Humphreys, Hiroki Ikawa, Pierre‐Erik Isabelle, Hiroki Iwata, Rachhpal S. Jassal, Mika Korkiakoski, J. Kurbatova, Lars Kutzbach, Е. Д. Лапшина, Anders Lindroth, Mikaell Ottosson Löfvenius, Annalea Lohila, Ivan Mammarella, Philip Marsh, Paul Moore, Trofim C. Maximov, Daniel F. Nadeau, Erin M. Nicholls, Mats B. Nilsson, Takeshi Ohta, Matthias Peichl, Richard M. Petrone, Anatoly Prokushkin, W. L. Quinton, Nigel T. Roulet, Benjamin R. K. Runkle, Oliver Sonnentag, Ian B. Strachan, Pierre Taillardat, Eeva‐Stiina Tuittila, Juha‐Pekka Tuovinen, Jessica Turner, Masahito Ueyama, Andrej Varlagin, Timo Vesala, Martin Wilmking, Vyacheslav Zyrianov, Christopher Schulze
Environmental Research Letters, Volume 15, Issue 10

Peatlands and forests cover large areas of the boreal biome and are critical for global climate regulation. They also regulate regional climate through heat and water vapour exchange with the atmosphere. Understanding how land-atmosphere interactions in peatlands differ from forests may therefore be crucial for modelling boreal climate system dynamics and for assessing climate benefits of peatland conservation and restoration. To assess the biophysical impacts of peatlands and forests on peak growing season air temperature and humidity, we analysed surface energy fluxes and albedo from 35 peatlands and 37 evergreen needleleaf forests - the dominant boreal forest type - and simulated air temperature and vapour pressure deficit (VPD) over hypothetical homogeneous peatland and forest landscapes. We ran an evapotranspiration model using land surface parameters derived from energy flux observations and coupled an analytical solution for the surface energy balance to an atmospheric boundary layer (ABL) model. We found that peatlands, compared to forests, are characterized by higher growing season albedo, lower aerodynamic conductance, and higher surface conductance for an equivalent VPD. This combination of peatland surface properties results in a ∼20% decrease in afternoon ABL height, a cooling (from 1.7 to 2.5 °C) in afternoon air temperatures, and a decrease in afternoon VPD (from 0.4 to 0.7 kPa) for peatland landscapes compared to forest landscapes. These biophysical climate impacts of peatlands are most pronounced at lower latitudes (∼45°N) and decrease toward the northern limit of the boreal biome (∼70°N). Thus, boreal peatlands have the potential to mitigate the effect of regional climate warming during the growing season. The biophysical climate mitigation potential of peatlands needs to be accounted for when projecting the future climate of the boreal biome, when assessing the climate benefits of conserving pristine boreal peatlands, and when restoring peatlands that have experienced peatland drainage and mining. © 2020 The Author(s). Published by IOP Publishing Ltd. (Less)

2019

DOI bib abs
Increased high‐latitude photosynthetic carbon gain offset by respiration carbon loss during an anomalous warm winter to spring transition
Zhihua Liu, John S. Kimball, N. Parazoo, Ashley P. Ballantyne, Wen J. Wang, Nima Madani, Caleb G. Pan, Jennifer D. Watts, Rolf H. Reichle, Oliver Sonnentag, Philip Marsh, Miriam Hurkuck, Manuel Helbig, W. L. Quinton, Donatella Zona, Masahito Ueyama, Hideki Kobayashi, E. S. Euskirchen
Global Change Biology, Volume 26, Issue 2

Arctic and boreal ecosystems play an important role in the global carbon (C) budget, and whether they act as a future net C sink or source depends on climate and environmental change. Here, we used complementary in situ measurements, model simulations, and satellite observations to investigate the net carbon dioxide (CO2 ) seasonal cycle and its climatic and environmental controls across Alaska and northwestern Canada during the anomalously warm winter to spring conditions of 2015 and 2016 (relative to 2010-2014). In the warm spring, we found that photosynthesis was enhanced more than respiration, leading to greater CO2 uptake. However, photosynthetic enhancement from spring warming was partially offset by greater ecosystem respiration during the preceding anomalously warm winter, resulting in nearly neutral effects on the annual net CO2 balance. Eddy covariance CO2 flux measurements showed that air temperature has a primary influence on net CO2 exchange in winter and spring, while soil moisture has a primary control on net CO2 exchange in the fall. The net CO2 exchange was generally more moisture limited in the boreal region than in the Arctic tundra. Our analysis indicates complex seasonal interactions of underlying C cycle processes in response to changing climate and hydrology that may not manifest in changes in net annual CO2 exchange. Therefore, a better understanding of the seasonal response of C cycle processes may provide important insights for predicting future carbon-climate feedbacks and their consequences on atmospheric CO2 dynamics in the northern high latitudes.

DOI bib abs
Large loss of CO2 in winter observed across the northern permafrost region
Susan M. Natali, Jennifer D. Watts, Brendan M. Rogers, Stefano Potter, S. Ludwig, A. K. Selbmann, Patrick F. Sullivan, Benjamin W. Abbott, Kyle A. Arndt, Leah Birch, Mats P. Björkman, A. Anthony Bloom, Gerardo Celis, Torben R. Christensen, Casper T. Christiansen, R. Commane, Elisabeth J. Cooper, P. M. Crill, C. I. Czimczik, S. P. Davydov, Jinyang Du, J. E. Egan, Bo Elberling, E. S. Euskirchen, Thomas Friborg, Hélène Genet, Mathias Göckede, Jordan P. Goodrich, Paul Grogan, Manuel Helbig, Elchin Jafarov, Julie Jastrow, Aram Kalhori, Yongwon Kim, John S. Kimball, Lars Kutzbach, Mark J. Lara, Klaus Steenberg Larsen, Bang-Yong Lee, Zhihua Liu, M. M. Loranty, Magnus Lund, Massimo Lupascu, Nima Madani, Avni Malhotra, Roser Matamala, Jack W. McFarland, A. David McGuire, Anders Michelsen, Christina Minions, Walter C. Oechel, David Olefeldt, Frans‐Jan W. Parmentier, Norbert Pirk, Benjamin Poulter, W. L. Quinton, Fereidoun Rezanezhad, David Risk, Torsten Sachs, Kevin Schaefer, Niels Martin Schmidt, Edward A. G. Schuur, Philipp Semenchuk, Gaius R. Shaver, Oliver Sonnentag, Gregory Starr, Claire C. Treat, Mark P. Waldrop, Yihui Wang, J. M. Welker, Christian Wille, Xiaofeng Xu, Zhen Zhang, Qianlai Zhuang, Donatella Zona
Nature Climate Change, Volume 9, Issue 11

Recent warming in the Arctic, which has been amplified during the winter1-3, greatly enhances microbial decomposition of soil organic matter and subsequent release of carbon dioxide (CO2)4. However, the amount of CO2 released in winter is highly uncertain and has not been well represented by ecosystem models or by empirically-based estimates5,6. Here we synthesize regional in situ observations of CO2 flux from arctic and boreal soils to assess current and future winter carbon losses from the northern permafrost domain. We estimate a contemporary loss of 1662 Tg C yr-1 from the permafrost region during the winter season (October through April). This loss is greater than the average growing season carbon uptake for this region estimated from process models (-1032 Tg C yr-1). Extending model predictions to warmer conditions in 2100 indicates that winter CO2 emissions will increase 17% under a moderate mitigation scenario-Representative Concentration Pathway (RCP) 4.5-and 41% under business-as-usual emissions scenario-RCP 8.5. Our results provide a new baseline for winter CO2 emissions from northern terrestrial regions and indicate that enhanced soil CO2 loss due to winter warming may offset growing season carbon uptake under future climatic conditions.

DOI bib abs
Monthly gridded data product of northern wetland methane emissions based on upscaling eddy covariance observations
Olli Peltola, Timo Vesala, Yao Gao, Olle Räty, Pavel Alekseychik, Mika Aurela, Bogdan H. Chojnicki, Ankur R. Desai, A. J. Dolman, E. S. Euskirchen, Thomas Friborg, Mathias Göckede, Manuel Helbig, Elyn Humphreys, Robert B. Jackson, Georg Jocher, Fortunat Joos, Janina Klatt, Sara Knox, Natalia Kowalska, Lars Kutzbach, Sebastian Lienert, Annalea Lohila, Ivan Mammarella, Daniel F. Nadeau, Mats B. Nilsson, Walter C. Oechel, Matthias Peichl, Thomas G. Pypker, W. L. Quinton, Janne Rinne, Torsten Sachs, Mateusz Samson, Hans Peter Schmid, Oliver Sonnentag, Christian Wille, Donatella Zona, Tuula Aalto
Earth System Science Data, Volume 11, Issue 3

Abstract. Natural wetlands constitute the largest and most uncertain source of methane (CH4) to the atmosphere and a large fraction of them are found in the northern latitudes. These emissions are typically estimated using process (“bottom-up”) or inversion (“top-down”) models. However, estimates from these two types of models are not independent of each other since the top-down estimates usually rely on the a priori estimation of these emissions obtained with process models. Hence, independent spatially explicit validation data are needed. Here we utilize a random forest (RF) machine-learning technique to upscale CH4 eddy covariance flux measurements from 25 sites to estimate CH4 wetland emissions from the northern latitudes (north of 45∘ N). Eddy covariance data from 2005 to 2016 are used for model development. The model is then used to predict emissions during 2013 and 2014. The predictive performance of the RF model is evaluated using a leave-one-site-out cross-validation scheme. The performance (Nash–Sutcliffe model efficiency =0.47) is comparable to previous studies upscaling net ecosystem exchange of carbon dioxide and studies comparing process model output against site-level CH4 emission data. The global distribution of wetlands is one major source of uncertainty for upscaling CH4. Thus, three wetland distribution maps are utilized in the upscaling. Depending on the wetland distribution map, the annual emissions for the northern wetlands yield 32 (22.3–41.2, 95 % confidence interval calculated from a RF model ensemble), 31 (21.4–39.9) or 38 (25.9–49.5) Tg(CH4) yr−1. To further evaluate the uncertainties of the upscaled CH4 flux data products we also compared them against output from two process models (LPX-Bern and WetCHARTs), and methodological issues related to CH4 flux upscaling are discussed. The monthly upscaled CH4 flux data products are available at https://doi.org/10.5281/zenodo.2560163 (Peltola et al., 2019).

DOI bib abs
Modelling the effects of permafrost loss on discharge from a wetland‐dominated, discontinuous permafrost basin
Lindsay E. Stone, Xing Fang, Kristine M. Haynes, Manuel Helbig, John W. Pomeroy, Oliver Sonnentag, W. L. Quinton
Hydrological Processes, Volume 33, Issue 20

Permafrost degradation in the peat‐rich southern fringe of the discontinuous permafrost zone is catalysing substantial changes to land cover with expansion of permafrost‐free wetlands (bogs and fens) and shrinkage of forest‐dominated permafrost peat plateaux. Predicting discharge from headwater basins in this region depends upon understanding and numerically representing the interactions between storage and discharge within and between the major land cover types and how these interactions are changing. To better understand the implications of advanced permafrost thaw‐induced land cover change on wetland discharge, with all landscape features capable of contributing to drainage networks, the hydrological behaviour of a channel fen sub‐basin in the headwaters of Scotty Creek, Northwest Territories, Canada, dominated by peat plateau–bog complexes, was modelled using the Cold Regions Hydrological Modelling platform for the period of 2009 to 2015. The model construction was based on field water balance observations, and performance was deemed adequate when evaluated against measured water balance components. A sensitivity analysis was conducted to assess the impact of progressive permafrost loss on discharge from the sub‐basin, in which all units of the sub‐basin have the potential to contribute to the drainage network, by incrementally reducing the ratio of wetland to plateau in the modelled sub‐basin. Simulated reductions in permafrost extent decreased total annual discharge from the channel fen by 2.5% for every 10% decrease in permafrost area due to increased surface storage capacity, reduced run‐off efficiency, and increased landscape evapotranspiration. Runoff ratios for the fen hydrological response unit dropped from 0.54 to 0.48 after the simulated 50% permafrost area loss with a substantial reduction of 0.47 to 0.31 during the snowmelt season. The reduction in peat plateau area resulted in decreased seasonal variability in discharge due to changes in the flow path routing, with amplified low flows associated with small increases in subsurface discharge, and decreased peak discharge with large reductions in surface run‐off.

2018

DOI bib abs
Minor contribution of overstorey transpiration to landscape evapotranspiration in boreal permafrost peatlands
Rebecca K. Warren, Christoforos Pappas, Manuel Helbig, L. Chasmer, Aaron Berg, Jennifer L. Baltzer, W. L. Quinton, Oliver Sonnentag
Ecohydrology, Volume 11, Issue 5

Evapotranspiration (ET) is a key component of the water cycle, whereby accurate partitioning of ET into evaporation and transpiration provides important information about the intrinsically coupled carbon, water, and energy fluxes. Currently, global estimates of partitioned evaporative and transpiration fluxes remain highly uncertain, especially for high‐latitude ecosystems where measurements are scarce. Forested peat plateaus underlain by permafrost and surrounded by permafrost‐free wetlands characterize approximately 60% (7.0 × 107 km2) of Canadian peatlands. In this study, 22 Picea mariana (black spruce) individuals, the most common tree species of the North American boreal forest, were instrumented with sap flow sensors within the footprint of an eddy covariance tower measuring ET from a forest–wetland mosaic landscape. Sap flux density (JS), together with remote sensing data and in situ measurements of canopy structure, was used to upscale tree‐level JS to overstorey transpiration (TBS). Black spruce trees growing in nutrient‐poor permafrost peat soils were found to have lower mean JS than those growing in mineral soils. Overall, TBS contributed less than 1% to landscape ET. Climate‐change‐induced forest loss and the expansion of wetlands may further minimize the contributions of TBS to ET and increase the contribution of standing water.

2017

DOI bib abs
Warmer spring conditions increase annual methane emissions from a boreal peat landscape with sporadic permafrost
Manuel Helbig, W. L. Quinton, Oliver Sonnentag
Environmental Research Letters, Volume 12, Issue 11

About a fifth of the global wetland methane emissions originate from boreal peatlands, which represent an important land cover type in boreal landscapes in the sporadic permafrost zone. There, rising air temperatures could lead to warmer spring and longer growing seasons, changing landscape methane emissions. To quantify the effect of warmer spring conditions on methane emissions of a boreal peat landscape in the sporadic permafrost zone of northwestern Canada, we analyzed four years (2013 – 2016) of methane fluxes measured with the eddy covariance technique and long-term (1951-2016) meteorological observations from a nearby climate station. In May, after snowmelt was complete, mean air temperatures were more than 2 °C warmer in 2013, 2015, and 2016 than in 2014. Mean growing season (May-August) air temperatures, in contrast, differed by less than 1 °C over the four years. Warmer May air temperatures caused earlier wetland soil warming, with temperatures rising from ~0 °C to g12 °C 25 to 40 days earlier and leading to ~6 °C warmer mean soil temperatures between May and June. However, from July to August, soil temperatures were similar among years. Mean May to August and annual methane emissions (6.4 g CH4 m-2 and 9.4 g CH4 m-2, respectively) of years with warmer spring (i.e., May) temperatures exceeded emissions during the cooler year by 20-30 % (4.5 g CH4 m-2 and 7.2 g CH4 m-2, respectively). Among years with warmer springs, growing season methane emissions varied little (0.5 g CH4 m-2). The observed interannual differences are most likely caused by a strong soil temperature control on methane fluxes and large soil temperature differences during the spring. Thus, in a warming climate, methane emissions from waterlogged boreal peat landscapes at the southern limit of permafrost are likely to increase in response to more frequent occurrences of warm springs.