Articles | Volume 19, issue 24
Biogeosciences, 19, 5953–5972, 2022
Biogeosciences, 19, 5953–5972, 2022
Research article
 | Highlight paper
22 Dec 2022
Research article  | Highlight paper | 22 Dec 2022

Using atmospheric observations to quantify annual biogenic carbon dioxide fluxes on the Alaska North Slope

Using atmospheric observations to quantify annual biogenic carbon dioxide fluxes on the Alaska North Slope
Luke D. Schiferl1,2, Jennifer D. Watts3, Erik J. L. Larson4, Kyle A. Arndt3,5,6, Sébastien C. Biraud7, Eugénie S. Euskirchen8, Jordan P. Goodrich5,9, John M. Henderson10, Aram Kalhori5,11, Kathryn McKain12,13, Marikate E. Mountain10, J. William Munger2, Walter C. Oechel5,14, Colm Sweeney12, Yonghong Yi15,16, Donatella Zona5,17, and Róisín Commane1,18 Luke D. Schiferl et al.
  • 1Lamont-Doherty Earth Observatory, Columbia University, Palisades, New York, USA
  • 2Harvard John A. Paulson School of Engineering and Applied Sciences, Cambridge, Massachusetts, USA
  • 3Woodwell Climate Research Center, Falmouth, Massachusetts, USA
  • 4Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, Massachusetts, USA
  • 5Department of Biology, San Diego State University, San Diego, California, USA
  • 6Earth Systems Research Center, Institute for the Study of Earth, Oceans, and Space, University of New Hampshire, Durham, New Hampshire, USA
  • 7Lawrence Berkeley National Laboratory, Berkeley, California, USA
  • 8Institute of Arctic Biology, University of Alaska Fairbanks, Fairbanks, Alaska, USA
  • 9Ministry for the Environment, Wellington, New Zealand
  • 10Atmospheric and Environmental Research, Inc., Lexington, Massachusetts, USA
  • 11GFZ German Research Centre for Geosciences, Potsdam, Germany
  • 12Global Monitoring Laboratory, Earth System Research Laboratories, NOAA, Boulder, Colorado, USA
  • 13Cooperative Institute for Research in Environmental Sciences, University of Colorado, Boulder, Colorado, USA
  • 14Department of Geography, University of Exeter, Exeter, UK
  • 15Joint Institute for Regional Earth System Science and Engineering, University of California, Los Angeles, California, USA
  • 16College of Surveying and Geo-Informatics, Tongji University, Shanghai, China
  • 17Department of Animal and Plant Sciences, University of Sheffield, Western Bank, Sheffield, UK
  • 18Department of Earth and Environmental Sciences, Columbia University, New York, New York, USA

Correspondence: Luke D. Schiferl (


The continued warming of the Arctic could release vast stores of carbon into the atmosphere from high-latitude ecosystems, especially from thawing permafrost. Increasing uptake of carbon dioxide (CO2) by vegetation during longer growing seasons may partially offset such release of carbon. However, evidence of significant net annual release of carbon from site-level observations and model simulations across tundra ecosystems has been inconclusive. To address this knowledge gap, we combined top-down observations of atmospheric CO2 concentration enhancements from aircraft and a tall tower, which integrate ecosystem exchange over large regions, with bottom-up observed CO2 fluxes from tundra environments and found that the Alaska North Slope is not a consistent net source nor net sink of CO2 to the atmosphere (ranging from −6 to +6Tg C yr−1 for 2012–2017). Our analysis suggests that significant biogenic CO2 fluxes from unfrozen terrestrial soils, and likely inland waters, during the early cold season (September–December) are major factors in determining the net annual carbon balance of the North Slope, implying strong sensitivity to the rapidly warming freeze-up period. At the regional level, we find no evidence of the previously reported large late-cold-season (January–April) CO2 emissions to the atmosphere during the study period. Despite the importance of the cold-season CO2 emissions to the annual total, the interannual variability in the net CO2 flux is driven by the variability in growing season fluxes. During the growing season, the regional net CO2 flux is also highly sensitive to the distribution of tundra vegetation types throughout the North Slope. This study shows that quantification and characterization of year-round CO2 fluxes from the heterogeneous terrestrial and aquatic ecosystems in the Arctic using both site-level and atmospheric observations are important to accurately project the Earth system response to future warming.

1 Introduction

The Arctic surface air temperature is warming at twice the rate of the global average (Box et al., 2019; Meredith et al., 2019). Continued thawing of Arctic permafrost has the potential to release vast stores of carbon into the atmosphere, thereby further accelerating warming (Schuur et al., 2015; Hugelius et al., 2014). In the biosphere, the net CO2 flux is the balance between the uptake of CO2 by vegetation through photosynthesis (negative net CO2 flux indicates removal from the atmosphere) and the release of CO2 into the atmosphere by plant and microbial respiration (positive net CO2 flux indicates a source to the atmosphere). Arctic growing seasons are short (∼3 months), and the long cold season dominates the seasonal cycle. The transition between the growing and cold seasons is marked by the soil zero-curtain period, when belowground temperatures of the active layer above frozen permafrost remain near freezing; the active layer is insulated by snow and ice at the surface and warmed by the latent heat release of freezing water (Outcalt et al., 1990). During the zero-curtain period, soil respiration can remain active in deeper soils for weeks to months after the end of the growing season (Zona et al., 2016; Romanovsky and Osterkamp, 2000). As the climate warms, the active layer above permafrost deepens, thawed soils become wetter, a larger volume of soil remains unfrozen for a longer period of time, and the duration of the zero-curtain period plays an increasingly important role in determining the net carbon exchange in Arctic ecosystems (Kim et al., 2012; Arndt et al., 2019). Recent work has shown a significant cold-season source of CO2 from Arctic ecosystems, including more than a 70 % increase in October–December CO2 concentration enhancements in the past 40 years, consistent with an increase in cold-season respiration, which is not well represented in Earth system models (Commane et al., 2017; Natali and Watts et al., 2019). Neglecting these processes could lead to a large underestimation of CO2 emissions, biassing current and future climate projections.

Figure 1Alaska North Slope study region, eddy flux site locations, area sampled by aircraft and the tower, and example results from the eddy flux site measurement–model comparison. (a) North Slope region (red box) within Alaska and northwestern Canada; tundra areas are shown in purple, and boreal forest areas are shown in green (Luus et al., 2017). (b) Location of eddy flux measurement sites on the Alaska North Slope used in this analysis. (c) The 10 d WRF-STILT (Weather Research and Forecasting–Stochastic Time-Inverted Lagrangian Transport) footprints used to sample CO2 flux models, summed for all aircraft and tall-tower CO2 observations used in this analysis; colors represent values greater than 0 and are saturated at 60 ppm(µmolm-2s-1)-1, and the maximum value near Utqiaġvik, Alaska, is 324 ppm(µmolm-2s-1)-1. (d) Time series of the observed (black dots) and simulated (colored lines) site-level daily mean net CO2 flux for 2014 at the Ivotuk (IVO; left) and Climate Monitoring and Diagnostics Laboratory (CMDL; right) eddy flux measurement sites, where site-level Tundra Vegetation Photosynthesis and Respiration Model (TVPRM) net CO2 flux simulations are driven by North American Regional Reanalysis (NARR) meteorology and the solar-induced chlorophyll fluorescence (SIF) product from the contiguous SIF (CSIF) dataset. Positive net CO2 flux values indicate CO2 fluxes into the atmosphere throughout this study. A comparison for all eight eddy flux sites is provided in Fig. S1 in the Supplement.

Tundra ecosystems, characterized by frozen soils covered in low shrubs, sedges, grasses, and mosses, make up approximately 50 % of the Arctic landscape (Raynolds et al., 2019). Due to the lack of trees, the magnitude of net CO2 uptake in tundra during the growing season is relatively small and may be offset by emissions from respiration that can continue well into the cold season (Watts et al., 2021). In the past, year-round CO2 flux measurements from tundra ecosystems were rare due to difficulties in maintaining instrumentation under remote and extreme cold conditions (Euskirchen et al., 2017; Kittler et al., 2017; Goodrich et al., 2016). Long-term year-round CO2 concentration measurements have been made in the Arctic at a small number of tall towers, which have been situated to sample clean marine air off the ocean (Jeong et al., 2018; Worthy et al., 2009). While aircraft provide greater spatial coverage over land than these towers, they tend to operate for short durations, and their temporal coverage is limited by weather and visibility during the cold season (Chang et al., 2014; Commane et al., 2017; Miller et al., 2016). However, the recent increase in the availability of observations of gas fluxes and concentrations within a particular tundra region, the Alaska North Slope (Fig. 1a), is making it possible to better conduct year-round multi-scale assessments of tundra ecosystems, with the aim of improving our understanding of CO2 sink/source activity and carbon budgets in these environments.

Currently, observations and models do not agree on the sign of the annual net CO2 flux across the Alaska North Slope region. Site-level measurements and atmospheric observations suggest that this region is a net CO2 source (Commane et al., 2017; Oechel et al., 2014; Euskirchen et al., 2017). However, a comparison of process-based models of the North Slope found large variability in the sign and magnitude of the net CO2 flux with an approximately neutral regional annual net CO2 flux multi-model mean of -3.5±67Tg C yr−1 (Fisher et al., 2014). In a more recent study, Tao et al. (2021) found an annual net CO2 flux range of −9 to 12 Tg C yr−1 for the years 2010–2016, with only 2014 being an annual net CO2 source. Extrapolating from site-level CO2 flux measurements to regional budgets is difficult due to the extreme heterogeneity of tundra ecosystems in the North Slope region as well as a lack of spatial and seasonal representativeness by existing flux monitoring sites (Pallandt et al., 2022).

In this study, we compare bottom-up flux estimates with top-down atmospheric observations from aircraft and a tall tower using an integrated modeling approach to quantify the CO2 budget sign and magnitude of the Alaska North Slope. Our framework first applies a bottom-up approach to understand Arctic tundra ecosystem CO2 fluxes, constrained by site-level observations, using an empirical model ensemble of CO2 fluxes derived from eddy flux measurements representing varied tundra ecosystems within the region. We then apply top-down information gained from regional CO2 concentration enhancement observations measured by a tall tower and aircraft, which sample the atmosphere–biosphere exchange throughout the Alaska North Slope, to evaluate the range of potential CO2 fluxes identified by the bottom-up model ensemble for 2012–2017. This evaluation also identifies the ecosystem parameterizations, vegetation distributions, and environmental drivers that best characterize the observed spatial and temporal distribution of biogenic CO2 in the atmosphere across the region. By developing regional CO2 budgets constrained by both atmospheric observations and ecosystem environmental responses, we can better project how Arctic tundra ecosystems will respond to climate change on annual and decadal timescales.

2 Materials and methods

2.1 Observed CO2 concentrations and fluxes on the Alaska North Slope

2.1.1 Atmospheric CO2 concentration observations

We use a suite of CO2 concentration observations from various sources on the North Slope for our analysis. The United States (US) National Oceanic and Atmospheric Administration (NOAA) Barrow Atmospheric Baseline Observatory (BRW) tall tower near Utqiaġvik, Alaska, has made continuous in situ CO2 concentration measurements since 1973 (Sweeney et al., 2016). The US Department of Energy (DOE) Atmospheric Radiation Measurement Climate Research Facility Airborne Carbon Measurements V (ARM-ACME V) airborne campaign measured CO2 concentrations sub-weekly from June to September 2015 over the North Slope (Biraud et al., 2016; Tadić et al., 2021). The US National Aeronautics and Space Administration (NASA) Arctic-Boreal Vulnerability Experiment (ABoVE) Arctic Carbon Atmospheric Profiles (Arctic-CAP) airborne campaign flew throughout Alaska and northwestern Canada approximately every month from May to November 2017 (Sweeney and McKain, 2019; Sweeney et al., 2022). Carbon dioxide concentration observations from the NASA Carbon in Arctic Reservoirs Vulnerability Experiment (CARVE) flights for 2012–2014 are incorporated into the Commane et al. (2017) optimized CO2 fluxes used in our analysis below. The NOAA/US Coast Guard collaborative Alaska Coast Guard (ACG) flights have also made aircraft CO2 concentration measurements in the region, but these coastal flights observe only limited spatial coverage of the North Slope, and we do not use them here.

For the NOAA BRW tower, we use hourly CO2 concentration observations with wind direction from the land (135–202.5 clockwise with respect to north) and ocean sectors (0–45), avoiding Utqiaġvik anthropogenic activity, with wind speeds > 2.5 m s−1 (Fig. S2) (Commane et al., 2017; Sweeney et al., 2016). We only use land-sector observations from the cold season (defined here as September–April) because seasonal wind patterns do not favor transport from those directions during the growing season (defined here as May–August). For the ARM-ACME V and ABoVE Arctic-CAP aircraft campaign observations, we group averaged sampling points into 50 m vertical bins after removing data influenced by combustion sources such as anthropogenic activity and biomass burning events. These combustion sources of CO2 are expected to be small (<1Tg C yr−1 on the North Slope; see Table S1 in the Supplement) during our study period. They are not accounted for in biogenic CO2 flux models, however, and must be removed from our analysis when observed. We remove time periods with an elevated carbon monoxide (CO) concentration (above 150 ppb), as in Chang et al. (2014) and Commane et al. (2017), which indicates local combustion sources. Time periods with highly variable CO concentrations (ΔCO>40ppb) indicate complex mixing of more remote combustion sources and are also removed (Chang et al., 2014). The remaining grouped sampling points correspond to the available Weather Research and Forecasting–Stochastic Time-Inverted Lagrangian Transport (WRF-STILT) modeling system simulations (Henderson et al., 2015; see below): ARM-ACME V points are calculated every 50 m vertically below 1 km, every 100 m vertically above 1 km, and every 10 km horizontally from 1 s observations, and ABoVE Arctic-CAP points are matched every 20 s from averaged 10 s observations. To ensure these points observe the Alaska North Slope, we only use points with at least 70 % of the total 10 d WRF-STILT-simulated surface influence occurring in our regional domain.

2.1.2 Eddy covariance CO2 flux tower observations

We also use up to 5 years (2013–2017) of year-round observations of net CO2 flux from eight eddy covariance tower sites (for 32 total site-years) representing an array of tundra ecosystems throughout the Alaska North Slope (Figs. 1b, S1; Table S2). These half-hourly eddy flux measurements of net CO2 flux are not gap filled to avoid introducing additional uncertainties. Three of the sites are located near Imnavait Creek along a wetness gradient from valley to hilltop: wet sedge tundra (Imnavait Creek sedge – ICS), moist acidic tussock tundra (Imnavait Creek tussock – ICT), and dry heath tundra (Imnavait Creek heath – ICH) (Euskirchen et al., 2017, 2012). The other sites include tussock tundra at Ivotuk (IVO), wet polygonized tundra at Atqasuk (ATQ), and three sites near Utqiaġvik: wetland tundra (Biocomplexity Experiment, South – BES), wet polygonized tundra (Barrow Environmental Observatory – BEO), and moist tundra (Climate Monitoring and Diagnostics Laboratory – CMDL) (Zona et al., 2016; Arndt et al., 2020).

2.2 Observed atmospheric CO2 concentration enhancement calculation

We calculate the observed top-down atmospheric CO2 concentration enhancement (ΔCO2) for the North Slope region for every land-sector hour at the NOAA BRW tower and for every 50 m of vertical distance transited during the airborne campaigns (ARM-ACME V and ABoVE Arctic-CAP). The observed ΔCO2 (in units of ppm) generated by the North Slope ecosystem is calculated relative to the background concentration without influence from this region such that

(1) observed  Δ CO 2 = observed  [ CO 2 ] - background  [ CO 2 ]

following previous work (Sweeney et al., 2016; Commane et al., 2017; Jeong et al., 2018).

The background CO2 concentrations at the NOAA BRW tower are determined by smoothing the 10 d mean of the observed ocean-sector concentrations using spline fitting to produce a daily CO2 background concentration. We calculate the uncertainty of these background concentrations by both (1) varying the starting hour of the 10 d mean calculation prior to spline fitting and (2) randomly sub-selecting 50 % of the ocean-sector concentrations 1000 times. The interval that contains 95 % of these 240 000 fits represents our daily background uncertainty. Figure S2 shows the ocean-sector concentrations, the resulting background concentration, and the uncertainty described here.

To determine the background CO2 concentrations for the ARM-ACME V and ABoVE Arctic-CAP aircraft campaigns, we isolate aircraft observations without surface influence from the North Slope using the WRF-STILT footprints, as done for larger regions in Chang et al. (2014) and Commane et al. (2017). These observed CO2 concentrations represent the state of the air before it interacts with the surface in the study region. The regional backgrounds vary by the direction from which the air enters the domain. For example, the backgrounds from the south and from over land generally experience CO2 drawdown prior to those from over the Arctic Ocean. The time- and direction-dependent backgrounds that we use are shown in Fig. S3. We apply the uncertainty from the NOAA BRW tower background to the aircraft backgrounds as a reasonable representation of the variability associated with available background CO2 concentration data.

2.3 Simulated atmospheric CO2 concentration enhancement calculation

To understand how landscape interactions with the atmosphere (through CO2 flux) influenced the observed CO2 concentrations across space and time, we calculate the corresponding simulated ΔCO2 (in units of ppm) by transporting bottom-up biogenic CO2 fluxes to each observation site such that

(2) simulated  Δ CO 2 = simulated  CO 2  flux  ×  simulated footprint .

In this calculation, we multiply the hourly simulated CO2 flux (in units of µmolCO2m-2s-1) by the footprint (in units of ppm(µmolCO2m-2s-1)-1) for that hour starting at the observation point, backward in time for each hour up to 10 d, where the footprint quantifies the influence of the land surface on the concentration observed at a measurement point. The simulated ΔCO2 is then the sum of these hours.

We use expected CO2 fluxes based on a variety of bottom-up model approaches which represent North Slope ecosystems. Year-round bottom-up estimates of net CO2 fluxes (defined by the models as net ecosystem exchange, NEE) are obtained from the Tundra Vegetation Photosynthesis and Respiration Model (TVPRM) ensemble as well as from existing model output from Luus et al. (2017) and Commane et al. (2017). Independent bottom-up estimates of belowground CO2 emissions (equal to the NEE) for the cold season (net CO2 uptake = 0) were obtained from Natali and Watts et al. (2019) and Watts et al. (2021). The TVPRM model ensemble development process is described in Sect. 2.4, and the other CO2 flux models, including their native spatial and temporal resolutions, are listed in Table S3.

The footprints are generated by the WRF-STILT atmospheric transport modeling system (Henderson et al., 2015). In this system, WRF meteorological fields are first generated for the study region and time period (v3.5.1 for ARM-ACME V and NOAA BRW tower footprints used here, and v3.9.1 for ABoVE Arctic-CAP footprints). STILT then uses the WRF meteorology to estimate the contribution of surface fluxes to the atmospheric concentration at a specified time and place, called a receptor, by calculating the amount of time that air (represented by a distribution of particles) spends in the lower half of the boundary layer at a given location. The WRF-STILT model configurations from Henderson et al. (2015) have been extensively employed in numerous previous papers to study greenhouse gas fluxes using observations from aircraft and towers in Alaska, including on the North Slope (e.g., Chang et al., 2014; Miller et al., 2016; Zona et al., 2016; Commane et al., 2017; Karion et al., 2015; Hartery et al., 2018). An evaluation by Henderson et al. (2015) for WRF v.3.4.1 and v3.5.1 showed that their polar WRF configuration performs well against surface observations of air temperature and wind speed in Alaska and that WRF-STILT can capture the shape and approximate depth of greenhouse gases in the column. Zona et al. (2016) note that WRF planetary boundary layer ventilation rates may be biased in the fall (and winter) when heat fluxes are low, but this error is difficult to assess quantitatively. For this study, we use receptors set to correspond to the tower and aircraft CO2 concentration observations. The footprints (and their corresponding measurements) for these receptors sample air from throughout the North Slope but are concentrated more heavily toward the area around the NOAA BRW tower (Fig. 1c).

For calculating simulated ΔCO2 from the TVPRM ensemble, we grid the distribution of WRF-STILT particles and their corresponding surface influence to the spatial resolution of the meteorological reanalysis products driving the model. The CO2 flux models used for comparison to the TVPRM ensemble are similarly treated using 0.5 gridded 10 d WRF-STILT footprints, which are available on a circumpolar grid poleward of 30 N. The simulated CO2 fluxes from Luus et al. (2017), Natali and Watts et al. (2019), and Watts et al. (2021) are regridded to a 0.5 spatial resolution. For the models by Natali and Watts et al. (2019) and Watts et al. (2021), which only estimate monthly CO2 fluxes, we apply a constant flux for that month. As the ends of our defined cold season (September–April) include transitional periods when some biogenic plant activity does occur (hence belowground CO2 emissions  NEE), we add in estimates of photosynthesis and plant respiration fluxes from the TVPRM ensemble for April and September for the Natali and Watts et al. (2019) and Watts et al. (2021) bottom-up scenarios.

2.4 Empirically simulated biogenic CO2 fluxes from tundra ecosystems

We develop the TVPRM as an ensemble of ecosystem-resolved models that represent a more extensive range of potential tundra ecosystem functional relationships, environmental drivers, and scaling assumptions than available from other CO2 flux models. For this study, TVPRM generates a set of spatially and temporally varying CO2 flux maps for a 6-year period (2012–2017) at a 30 km× 30 km spatial and a 1 h temporal resolution for the Alaska North Slope.

TVPRM is driven by parameterized functional relationships for soil respiration (Rsoil), plant respiration (Rplant), and photosynthesis (gross primary productivity, GPP), which are described by the following respective equations:


The simulated hourly CO2 fluxes (in units of µmolCO2m-2s-1) are determined as responses to light and heat: Rsoil is a function of near-surface soil temperature (Ts; C); Rplant is a function of air temperature (Ta; C); and GPP is a function of a temperature scalar (Tscale) and photosynthetically active radiation (PAR; µmolphotonm-2s-1), with solar-induced chlorophyll fluorescence (SIF; mWm-2nm-1sr-1) used to define the seasonal cycle of photosynthetic capacity. Ts depths are determined by reanalysis product and listed in Table S4. Tscale ranges from zero to one based on the position of Ta on the continuum between minimum temperature (Tmin= 0 C), maximum temperature (Tmax= 40 C), and optimal temperature (Topt= 15 C). NEE is then calculated as

(7) NEE = R soil + R plant - GPP .

Positive NEE values indicate a net source of CO2 into the atmosphere, and negative NEE values indicate net movement of CO2 into the biosphere. We use NEE to be synonymous with net CO2 flux. Using SIF, which correlates to photosynthetic activity (Porcar-Castell et al., 2014; Yang et al., 2015), in the modeling framework provides an advantage over indices such as the enhanced vegetation index (EVI) due to the limited canopy and evergreen nature of tundra ecosystems (Luus et al., 2017).

The parameter values (αs, βs, αa, βa, λ, and PAR0) for the site-level relationships used by TVPRM are determined first using the observed net CO2 fluxes from the eddy flux sites (see Sect. S1 in the Supplement). We determine the site-level parameters separately for each combination of reanalysis product (North American Regional Reanalysis – NARR, Mesinger et al., 2006, and the fifth-generation European Centre for Medium-Range Weather Forecasts atmospheric reanalysis – ERA5, Hersbach et al., 2020), which provide Ta, Ts, and PAR, and the SIF product (Global Ozone Monitoring Experiment-2 – GOME-2, Joiner et al., 2016; Global OCO-2, where OCO-2 is the Orbiting Carbon Observatory-2, SIF – GOSIF, Li and Xiao, 2019; and contiguous SIF – CSIF, Zhang et al., 2018) that will later be used to generate the regional TVPRM ensemble (Tables S4 and S5; see Sects. S2 and S3). Additional αs and βs parameters are determined using Ts from the remote-sensing-driven permafrost model (RS-PM; Yi et al., 2019, 2018) to test its implementation in TVPRM. RS-PM uses tailored input for Alaska permafrost zones, such as downscaled snow depth and aircraft-observed soil dielectric constants, and was developed and tested using Ts and active-layer thickness measurements from the North Slope. RS-PM also produces Ts at a higher vertical resolution in the near-surface region than the reanalysis products in order to capture subsurface heterogeneity in unfrozen soil, which is important to represent the zero-curtain throughout the freezing and thawing periods in Alaska.

Using the median parameter value sets for each site, we simulate the TVPRM net CO2 flux for our study period at every site location to perform a cross-site evaluation (Fig. S1). These simulated net CO2 fluxes perform well against the net CO2 flux observations at their corresponding sites (Figs. 1d, S4; see Sect. S4). This process also identifies two distinct ecosystem groups – “inland”, predominately graminoid and shrub tundra (ICS, ICT, ICH, and IVO), and “coastal”, predominately wetland tundra (ATQ, BES, BEO, and CMDL) – based on the similar simulated CO2 flux responses to the meteorology- and SIF-determined functional relationships within each group demonstrated by the cross-site evaluation (Fig. S1).

The net CO2 flux for each meteorological grid box in our study domain is then calculated using the site-level functional relationships for both tundra groups. These fluxes are weighted by the spatial distribution of inland and coastal tundra from three different vegetation maps (Circumpolar Arctic Vegetation Map – CAVM, Walker et al., 2005; Raster Circumpolar Arctic Vegetation Map – RasterCAVM, Raynolds et al., 2019; and ABoVE land cover – ABoVE LC, Wang et al., 2020; the reader is also referred to Fig. S5, Table S6, and Sect. S5 in the Supplement) to produce the regionally scaled TVPRM net CO2 flux. By varying the choice of representative inland and coastal tundra sites, meteorological reanalysis product, vegetation map, and SIF product, we generate 288 different simulations (members) of net CO2 flux (referred to here as the unconstrained TVPRM ensemble) for each grid box across the region for each of the 6 study years. Monthly and annual regional net CO2 flux values are calculated as the area-weighted sum of all grid boxes simulated in our domain. Notable changes since the previous iteration of this empirical CO2 flux model (Commane et al., 2017; Luus et al., 2017) include the expansion of the model to include multiple ensemble members to account for variability and uncertainty in model formulation, the use of additional site-years of CO2 flux observations (with increased data coverage over the cold season), more inclusive data filtering methods, and much higher temporal (1, 4, and 8 d rather than monthly) and spatial (0.01 and 0.05 rather than 0.5) resolution SIF datasets. We compare TVPRM to the previous model version by Luus et al. (2017) and its CARVE-informed optimization by Commane et al. (2017) in Sect. 3.3.

2.5 Evaluation framework

We use the atmospheric CO2 concentration observations to evaluate the many tundra ecosystem parameterizations, vegetation distributions, and environmental drivers that represent the net CO2 flux on the North Slope over various spatial and temporal scales. For this assessment, we compare the observed ΔCO2 values, which are the observed CO2 concentration changes driven by regional CO2 fluxes, with the simulated ΔCO2 values determined by combining the regional biogenic CO2 flux models with the atmospheric transport model.

To compare the regional observed ΔCO2 and simulated ΔCO2, we calculated the coefficient of determination (R2) as the square of the Pearson correlation coefficient for all points. The slope is determined by ordinary least squares using the median of each 10 % bin of ordered observed and corresponding simulated net CO2 flux. The normalized mean bias (NMB) of all points is defined as (simulated-observed)observed. The root-mean-square error (RMSE) of all points is defined as (simulated-observed)2.

Figure 2Aircraft and tower CO2 concentration measurements constrain year-round simulated CO2 fluxes on the Alaska North Slope. (a–c) Comparison of observed and simulated ΔCO2 during the ARM-ACME V flight campaign (a), during the ABoVE Arctic-CAP flight campaign (b), and at the NOAA BRW tower (c) for air over the Alaska North Slope. Horizontal lines indicate the range of uncertainty in the NOAA BRW tower ocean-sector background calculation. Vertical boxes colored by month of the year represent 50 % and whiskers represent 95 % of the ΔCO2 values from all members of the unconstrained TVPRM ensemble (see Sect. 2.4) from all binned points. Black points show values from the constrained TVPRM member with additional zero-curtain (ZC) emissions and inland water (IW) fluxes (see Sect. 3.4). For panels (a) and (b), observed values are vertically binned medians, and vertical lines contain the middle 95 % of ΔCO2 values from all binned points for the constrained TVPRM member + ZC and IW. (d) Combined comparison of observed and simulated ΔCO2 for all aircraft and tower points using the constrained TVPRM member + ZC and IW. The linear best fit (red line), the slope determined by ordinary least squares, and the coefficient of determination (R2) of all points (n=455) are shown. The 1:1 is line shown in dark gray.


These comparisons enable us to constrain the regional net CO2 flux on the Alaska North Slope. First, we identify the year-round empirically driven net CO2 fluxes from the TVPRM ensemble (TVPRM Unconstrained) which are most consistent with the CO2 concentration observations from the two aircraft campaigns and at the tower (TVPRM Constrained) (Sect. 3.1 and 3.2). Then, noting the large range of potential cold-season CO2 fluxes, we compare our constrained TVPRM member with CO2 fluxes from previous studies (Sect. 3.3). Finally, we suggest and quantify sources of the missing CO2 flux observed during the early cold season (defined here as September–December) and incorporate those fluxes into our net CO2 budget (TVPRM Constrained + Additional Zero-Curtain Emissions and Inland Water Fluxes) (Sect. 3.4). This analysis provides a unique regional net CO2 flux quantification for the North Slope that is verified using atmospheric observations and can also be explained from an ecological and physical perspective.

3 Results

3.1 Evaluation of the unconstrained empirical net CO2 flux model ensemble

3.1.1 Using aircraft-observed CO2 enhancements

The observed ΔCO2 during the ARM-ACME V (June–September 2015) and ABoVE Arctic-CAP (May–November 2017) airborne campaigns show a strong seasonal uptake pattern throughout the growing season (Fig. 2a, b). The frequent flights during ARM-ACME V (multiple flights per week) observe the transition from early to peak growing season uptake (observed ΔCO2=11 ppm) and on into cold-season respiration, which results in net CO2 source conditions in September (+5 ppm). While less frequent, the ABoVE Arctic-CAP flights begin at the end of the cold season, extend later into following cold season, and cover a larger area of the North Slope. Peak growing season uptake observed by the ABoVE Arctic-CAP flights (14 ppm) is slightly stronger than for during ARM-ACME V, and the ABoVE Arctic-CAP flights observe a strong CO2 source throughout the North Slope (+10 ppm) by November. The difference in observed ΔCO2 during peak growing season uptake between 2015 and 2017 is likely similar to the uncertainty in the respective values and could be due to differences in areas of the North Slope sampled between years.

Figure 3Cold-season CO2 emissions for the inland tundra site parameterizations and comparison to tower observations. (a) Time series of simulated daily mean Alaska North Slope net CO2 flux for the median of all unconstrained TVPRM ensemble members using each of the three inland tundra site parameterizations: ICS (red), ICT (orange), and IVO (green). Yearly colored lines shown for September–April beginning in September 2012 and ending in April 2017. The same is shown for all eight eddy flux sites in Fig. S9. (b, c) Comparison of observed and simulated ΔCO2 at the NOAA BRW tower for air over the North Slope using the median of all unconstrained TVPRM ensemble members for the inland tundra site parameterizations at ICS (b) and IVO (c). All points are colored by day of year. The slope determined by ordinary least squares and the coefficient of determination (R2) of all points (n=191) are shown. The 1:1 line is also shown in dark gray.


The magnitude and timing of the observed net CO2 uptake throughout the growing season is generally well represented by the empirical net CO2 flux model ensemble (TVPRM Unconstrained; Figs. 2a, b, and S6). The median coefficients of determination (R2) and ordinary least squares slopes between the observed and simulated ΔCO2 for this time are 0.54 and 0.41 for ARM-ACME V and 0.82 and 0.72 for ABoVE Arctic-CAP, respectively. Only for the July observations during the ABoVE Arctic-CAP campaign do many members of the CO2 flux trend toward an underestimate of net CO2 uptake, with all points showing a much larger range of simulated values compared to ARM-ACME V. The net CO2 release tends to be overestimated by the TVPRM ensemble during the ABoVE Arctic-CAP seasonal transitions in May and September, but the observed Rsoil is consistently underestimated during November.

Given the large range of unconstrained representations of the regional CO2 flux, the accuracy in simulating the aircraft-observed ΔCO2 varies between TVPRM ensemble members. For example, members using the RasterCAVM vegetation map, which places less coastal tundra area cover in the south (Fig. S5), produce a smaller mean July net CO2 uptake flux (by  1 µmolm-2s-1, Fig. S7a) throughout the southern North Slope than members using other vegetation maps (CAVM and ABoVE LC), and this placement consistently underestimates the net ΔCO2 uptake during the growing season by 5–10 ppm compared with the aircraft observations (Fig. S8). Also, members driven by SIF products that integrate additional remote sensing and/or meteorological data (GOSIF and CSIF) better reflect the timing and magnitude of the peak season carbon uptake in tundra ecosystems than members produced by interpolated SIF retrievals (GOME-2 SIF product), which underestimate the observed CO2 uptake during July (Fig. S8).

Using these comparisons, we identify less-representative ensemble members that generally underestimate the observed ΔCO2 uptake during the growing season (RasterCAVM vegetation map and GOME-2 SIF product members). Removing these members from the TVPRM ensemble improves the collective performance of the remaining members during the growing season (Fig. S6), brings the median slope of agreement closer to 1 for both campaigns (improves from 0.53 to 0.64 and from 0.71 to 0.94 for ARM-ACME V and ABoVE Arctic-CAP, respectively), and reduces the median NMB (−0.34 to −0.03) and median RMSE (3.12 to 2.73) for ABoVE Arctic-CAP.

3.1.2 Using tower-observed CO2 enhancements

As seen with the September–November aircraft data, the observed ΔCO2 at the NOAA BRW tower (Fig. 2c) indicate that the CO2 source to the atmosphere increases substantially from September to peak in October and November (+12 ppm) before decreasing to near zero throughout the late cold season (January–April).

Most of the TVPRM ensemble members substantially underestimate the observed ΔCO2 in the early cold season (September–December) as the soils freeze, and some simulations produce too much CO2 in the late cold season when the soils are frozen (Fig. 2c). The cold-season CO2 flux differs greatest in magnitude and spatial extent between the ensemble members parameterized for the ICS and ICT inland tundra sites (Figs. 3a, S9, S10), with a net CO2 flux difference of  0.2 µmolm-2s-1 throughout the region (Fig. S7b).

While the magnitude of CO2 flux from ICS members better matches the observed ΔCO2 in the early cold season than that from other sites (Figs. 3b, c, and S11), the response to Ts at ICS shows only a modest decrease in CO2 flux between the early and late cold season (32 % decrease between October and March; Fig. 3a), resulting in an overestimate of the regional ΔCO2 in the late cold season. The CO2 flux response to Ts for ICT members is similar to that for ICS but lower in magnitude, and the simulated ΔCO2 from members of neither site performs well against the observations in both the early and late cold season. Therefore, ICS and ICT inland tundra responses to Ts are not representative of the regional ΔCO2 observed at the NOAA BRW tower throughout the entire cold season, and we remove those members from our TVPRM ensemble.

The observed net CO2 fluxes at the IVO inland tundra and CMDL coastal tundra sites both show prolonged zero-curtain emissions (Fig. S1) and respond strongly to Ts in the early cold season (Fig. S9). The stronger response of CO2 fluxes to Ts from the early to late cold season at IVO (70 % decrease by January; Fig. 3a) compared with at the Imnavait Creek sites produces TVPRM members that better represent the large regional decrease in ΔCO2 observed on the North Slope (Fig. 3c). While all coastal tundra sites respond similarly to Ts during the cold season, we determine that the CO2 flux magnitude at CMDL is most consistent with the regional observations (Fig. S11). Ts values from ERA5 remain warmer throughout the late cold season compared with those from NARR, which causes simulations using ERA5 Ts to overestimate CO2 release during that time (Fig. S11). Unlike during the growing season, cold-season CO2 fluxes are not sensitive to the vegetation distribution nor the SIF products.

Figure 4Tall-tower atmospheric observations of the Alaska North Slope support early-cold-season emissions not driven by soil temperature (Ts) and present no evidence of elevated late-cold-season emissions. (a, b) Comparison of hourly cold-season (September–April) observed and simulated ΔCO2 at the NOAA BRW tower for the constrained TVPRM member, where soil respiration (Rsoil) is determined only by Ts (a) and for the constrained TVPRM member + additional zero-curtain (ZC) emissions and inland water (IW) fluxes (b). Horizontal segments indicate the range of uncertainty in the NOAA BRW tower ocean-sector background calculation. The slope determined by ordinary least squares and the coefficient of determination (R2) of all points (n=191) are shown. The 1:1 line is also shown in dark gray. (c) Monthly cold-season total Alaska North Slope net CO2 fluxes for various CO2 flux models. TVPRM-based simulations and Natali and Watts et al. (2019) show values for 2012–2017, Luus et al. (2017) show 2012–2014, and Watts et al. (2021) show September 2016–April 2017. Ribbons represent the range of all years, where applicable. The area of the North Slope domain used to calculate regional totals is 3.537×105km2.


Finally, we identify the TVPRM member that best matches the observed ΔCO2: parameterized by IVO inland tundra and CMDL coastal tundra site responses, distributed by the ABoVE LC vegetation map, and driven by NARR reanalysis and the CSIF SIF product (referred to here as TVPRM Constrained; Figs. S6, S12). This constrained simulation estimates a mean regional CO2 flux of 0.05 µmolm-2s-1 for the late cold season in 2012–2015 and reproduces the observed ΔCO2 during this time well (Fig. 4a). The late-cold-season NMB and RMSE against the observations at the NOAA BRW tower are reduced from 4.91 to 2.04 and from 1.94 to 1.30, respectively, for the constrained simulation compared with the median of the entire TVPRM ensemble (Fig. S12). However, the early-cold-season CO2 emissions, with a mean regional CO2 flux of 0.25 µmolm-2s-1 for September–December (Fig. S13a), are still underestimated, with the simulated ΔCO2 lower than the observed ΔCO2 by ∼5ppm (Fig. 4a).

3.2 Alternative Ts products and Rsoil parameterizations

To test the impact of reanalysis Ts on the early-cold-season CO2 fluxes, we implement Ts values that are more specifically developed to represent Alaska tundra permafrost soils during freeze–thaw processes than the reanalysis products driving our constrained TVPRM member. A single layer of Ts at 8 cm depth from RS-PM (Fig. S14a) captures the magnitude and temporal behavior of the observed early-cold-season CO2 fluxes slightly better than the constrained member (Figs. 4a, S12), which uses NARR reanalysis Ts and does not incorporate permafrost-model-derived Ts. The RS-PM Ts extends CO2 emission fluxes further into the cold season by up to a month, which is consistent with a better representation of the zero-curtain period; however, emissions remain higher throughout the late cold season than our atmospheric-observation-constrained CO2 fluxes (Fig. S15). We also test the implementation of a multilayer fit driven by soil column temperature from RS-PM, but neither of these instances of remote-sensing-informed Ts substantially improve the agreement of the ΔCO2 at the NOAA BRW tower during the early cold season. Attempts to use alternative Rsoil formulations based on Ts, including Q10 relationships, also fail to reproduce the observed elevated CO2 fluxes during the cold season.

3.3 Evaluation of other CO2 flux models during the cold season

More early-cold-season (September–December) CO2 flux into the atmosphere is observed at the NOAA BRW tower than is emitted by our constrained empirical simulation member, and these observations also indicate low late-cold-season (January–April) CO2 emissions. We compare our constrained CO2 fluxes to several other representations of gridded CO2 flux on the North Slope (Table S3) and find that difficulty in simulating the magnitude and timing of the observed ΔCO2 throughout the cold season is not unique to the constrained fluxes from our study.

The net CO2 fluxes from Luus et al. (2017) are similar to the constrained TVPRM member during the growing season (Fig. S16) but release more than 3 times as much CO2 into the atmosphere throughout the late cold season (Fig. 4c). This large late-cold-season CO2 flux leads to a large overestimate compared with the observed ΔCO2 (Fig. S14b). The optimization employed by Commane et al. (2017) increases the September–October CO2 flux to a range that matches our observations at the NOAA BRW tower. However, Commane et al. (2017) did not optimize the cold-season fluxes from November to March, instead reverting to Luus et al. (2017) fluxes during this time, and thus produced late-cold-season fluxes that are too large. Overall, Commane et al. (2017) projected a regional total cold-season CO2 source of 37–40 Tg C for 2012–2014, which is more than twice as high as our constrained TVPRM member CO2 flux (15–18 Tg C) for those years.

Carbon dioxide fluxes from work by Natali and Watts et al. (2019), a cold-season model developed for the global high-latitude permafrost region, are similar to our constrained TVPRM member in September, but the fluxes remain high throughout the cold season (Fig. 4c), similarly to Luus et al. (2017), for a range of total cold-season CO2 flux of 40–43 Tg C for 2012–2017. This sustained CO2 release also leads to an overestimation in the ΔCO2 in the late cold season for this region (Fig. S14c). Tao et al. (2021) also show that the cold-season CO2 fluxes of Natali and Watts et al. (2019) are high compared with their model. More recent work by Watts et al. (2021), using observations from new Soil Respiration Station monitoring sites in Alaska, produces cold-season CO2 fluxes more similar to our constrained CO2 fluxes, with an underestimate in the simulated ΔCO2 during the early cold season (Fig. S14d), for a total cold-season CO2 flux of 18 Tg C for September 2016 to April 2017.

3.4 Sources of missing CO2 fluxes

None of the flux products discussed above, including our TVPRM ensemble, account for any potential CO2 fluxes during the zero-curtain period that are not driven by Ts or are from areas on the terrestrial–aquatic interface. To account for these processes, we first add an additional CO2 flux with zero-curtain timing to our constrained CO2 flux (TVPRM) member from both inland and coastal tundra areas that consists of 0.25 µmolm-2s-1 for October with a reduction to 0 by the end of December. This peak additional CO2 flux is within the daily variability in the observed CO2 flux at the IVO and CMDL eddy flux sites during the zero-curtain period (Fig. S9), and the reduction into December is consistent with these observations. The additional zero-curtain flux improves the ability of the model to reproduce the observed ΔCO2 at the NOAA BRW tower (slope = 0.46, R2= 0.41). We also apply the coastal tundra site ecosystem parameterization used in our constrained TVPRM member to all areas of inland water on the North Slope, which account for 4 % of the domain according to the ABoVE LC map (Fig. S5) and were previously set to zero CO2 flux. Representing these aquatic areas with biogenic CO2 fluxes consistent with coastal tundra ecosystems is one simple way to bridge the terrestrial–aquatic gap in tundra ecosystem models, where portions of aquatic systems on the land–water gradient (i.e., the edges) may be more likely to respond to the environment as coastal tundra than with the zero-flux assumed by water area. The ice phenology for areas of inland water producing CO2 flux is then considered to be similar to that of the freeze–thaw timing in coastal tundra soils. Adding these coastal tundra fluxes to inland water areas also improves the performance of our model (slope = 0.32 and R2= 0.29 against NOAA BRW tower observations). The magnitude of additional zero-curtain flux suggested here and the portion of inland water represented with coastal tundra site parameterizations produce the best statistical comparison for a range of choices tested (Fig. S17).

Together, adding these zero-curtain (ZC) and inland water (IW) CO2 fluxes to our constrained simulation (referred to as TVPRM Constrained + ZC and IW) increases the mean regional CO2 flux in the early cold season by 70 % (0.18 µmolm-2s-1; Fig. S13b) and results in a large improvement to our comparison of ΔCO2 at the NOAA BRW tower (slope =0.54, R2=0.40; Figs. 4b, S12) and across the region using airborne data, especially during the November ABoVE Arctic-CAP flights (Figs. 2, S6). The year-round comparison using all available aircraft and tower observations shows that these net CO2 fluxes are now representative of the region (slope =0.90, R2=0.80; Fig. 2d). As a result, the North Slope regional total cold-season CO2 flux increases by 6 Tg C (∼38 %) to 20–24 Tg C for 2012–2017 compared with the constrained empirical CO2 flux model member.

3.5 Alaska North Slope annual net CO2 flux

The median Alaska North Slope annual net CO2 flux from the TVPRM ensemble (5 Tg C yr−1) for 2012–2017 is consistent with the previous multi-model comparison (Fisher et al., 2014), but we find a much smaller range of regional CO2 flux values (26 to 29 Tg C yr−1 for 95 % of TVPRM members; Fig. S18). The largest contribution to this ensemble range comes from the difference in parameterizations determined for the ICS and ICT inland tundra sites, with TVPRM members using ICS trending toward a net CO2 source, while ICT trends toward net CO2 uptake. The distribution of inland and coastal tundra throughout the region represented by the vegetation maps also has a noticeable impact on the sign of the net CO2 flux, with members using the RasterCAVM more likely to release net CO2 into the atmosphere than members using the other maps. There is also little interannual variability in the unconstrained TVPRM ensemble, with only 2014 moving toward a net CO2 source, consistent with Tao et al. (2021) for these years.

Figure 5Annual and seasonal Alaska North Slope net CO2 flux constrained by aircraft and tower observations. (a) Annual, (b) late-cold-season (January–April), (c) growing-season (May–August), and (d) early-cold-season (September–December) total Alaska North Slope net CO2 fluxes for various CO2 flux models for 2012–2017, as in Fig. 4. Purple squares indicate the middle 95 % of all TVPRM ensemble members.


Our best quantification of the annual net CO2 flux for the North Slope informed by atmospheric observations, TVPRM Constrained + ZC and IW, indicates that the region is a small net sink for 2013 (5 Tg C yr−1) and 2015 (6 Tg C yr−1) and a small net source for 2012 (+6 Tg C yr−1), 2014 (+6 Tg C yr−1), 2016 (+2 Tg C yr−1), and 2017 (+2 Tg C yr−1) (Fig. 5a). We estimate a 10 % uncertainty in the net annual CO2 flux based on the slope from our final comparison with the year-round observations (Fig. 2d). The year-round net CO2 fluxes from Luus et al. (2017) (driven with NARR meteorology, monthly GOME-2 SIF, and the CAVM vegetation map) indicate the North Slope to be a strong annual net CO2 source for 2012–2014 (+9 to +15 Tg C yr−1; Fig. S18) and are inconsistent with our results. Our results are more consistent with Tao et al. (2021), but we find a smaller range in the magnitude of net CO2 flux over the same years and more years trending toward a net CO2 source.

We find that the regional net growing season CO2 uptake and the cold-season emissions on the North Slope are comparable in magnitude, so the net balance could depend on small perturbations in either flux. However, the regional cold-season CO2 emissions for these years were relatively similar from year to year: 18–21 Tg C for the early cold season (Fig. 5d), diminishing to only 2–3 Tg C for the late cold season (Fig. 5b). Therefore, the interannual variability in the regional carbon balance is largely driven by fluctuating net growing season CO2 fluxes during these years: greater net growing season uptake in 2013 and 2015 than in 2012, 2014, 2016, and 2017 (Fig. 5c).

4 Discussion

4.1 Tundra ecosystem growing season net CO2 fluxes

The good performance of the TVPRM ensemble against the atmospheric observations during the growing season indicates that the tundra ecosystems of the Alaska North Slope respond to light and heat, as quantified by PAR, Ts, and Ta, and that the net CO2 flux is largely controlled by the simple Rsoil, Rplant, and GPP relationships in the empirical model over this time.

The growing season of each year determines the sign of the regional annual net CO2 flux during our study period, with 2013 and 2015 being strong net sinks and 2014 being the strongest net source. The relative magnitude of each component of the net CO2 flux during the growing season (i.e., Rsoil, Rplant, and GPP) varies from year to year (Table S7) and helps explain the interannual variability in the net source or sink status of the North Slope. The growing season in 2015 was very warm, dry, and sunny in Alaska and resulted in extreme biomass burning activity outside of the North Slope (Table S1). High regional mean Ta and PAR (Table S8) and low accumulated precipitation (Table S9) in NARR confirm this was the case for North Slope as well, with high Ta and PAR contributing to a very high GPP. The growing season SIF signal from the CSIF product, which determines the seasonal cycle and relative magnitude of photosynthetic activity, is also large in 2015 (Table S8), further enhancing GPP. This year and others with a larger GPP component of NEE correspond to growing seasons with stronger SIF signals, which is an indicator of increased productivity and consistent with previous studies (e.g., Magney et al., 2019; Sun et al., 2017). While fairly high Ta and Ts in 2015 also result in high Rsoil and Rplant, respectively, this elevated respiration is not enough to offset the very high GPP and results in a large net CO2 sink. In contrast, the summer of 2014 was cool, wet, and cloudy, and the North Slope experienced a very low Ta, PAR, and SIF signal, producing very low GPP. Lower-than-normal Ta also results in very low Rplant, but (as for 2015) this is not enough to offset the extremely low uptake, resulting in a large net CO2 source for 2014. In 2013, the other growing season with a strong net CO2 sink, moderately high GPP combines with moderately low Rplant and very low Rsoil. Extremely low Ts causes this very low Rsoil, which, relative to moderate Ta and PAR, is likely a result of an above-average lingering snowpack into May (Table S9). This lingering snowpack is perhaps surprising given that the mean snowpack for the proceeding cold season was not particularly deep. The important impact that snow cover and the timing of snowmelt has on Ts and carbon response in tundra ecosystems has recently been emphasized (e.g., Kim et al., 2021) and is also supported by our work, which shows that the prevalence of snow in the spring may determine the sign of the regional net CO2 for an entire year.

The regional net CO2 flux is highly sensitive, however, to the distribution of tundra vegetation types (upland vs. coastal) throughout the North Slope during the growing season. Coastal tundra takes up more CO2 for a given unit PAR compared with inland tundra, based on the relationships between observed site-level net CO2 flux and PAR in this study (TVPRM parameters; Fig. S1), which could be evidence of an adaptation to lower light levels. This difference is consistent with Luus et al. (2017), who calculated greater uptake at “wetland” sites like Atqasuk and Barrow than at “graminoid tundra” sites like Ivotuk and Imnavait when all driver inputs are constant, and with Mbufong et al. (2014), who also found that peak growing season net uptake for constant light is greater at Barrow than at Ivotuk. The stronger CO2 uptake response of coastal tundra to light is important to consider due to the fact that the vegetation distributions assessed here, with more coastal tundra to the south (CAVM, Walker et al., 2005; ABoVE LC, Wang et al., 2020), better agree with the atmospheric observations. When considering the ability of coastal tundra to take up CO2 when moved toward the south, Patankar et al. (2013) saw that tundra plants exposed to additional intense light did not respond with additional uptake. Therefore, while the ecosystem response of the southern North Slope is more consistent with coastal ecosystems, it seems possible that these areas are misclassified in either our simplified two-tundra-type scheme or in the vegetation maps themselves. The large variability in net CO2 flux calculated using the different maps supports the importance of accurate ecosystem type locations in upscaling eddy flux measurements and highlights the need for improved vegetation mapping and classification schemes in the Arctic ecology research community.

4.2 Regional-scale cold-season CO2 emissions

Observations across scales, at the in situ eddy flux towers, the NOAA BRW tower, and from aircraft, consistently show signs of large early-cold-season CO2 emissions from ecosystems on the Alaska North Slope. However, there is no evidence of widespread elevated emissions in this region during the late cold season, contrary to other studies (Commane et al., 2017; Natali and Watts et al., 2019). The TVPRM ensemble parameterizations using terrestrial eddy flux sites and the fluxes from other terrestrial CO2 models cannot reproduce both the observed magnitude and across-season timing of these cold-season CO2 emissions.

The largest differences in the net CO2 flux between TVPRM ensemble members result from the contrasting site conditions driving the ICS and ICT Rsoil parameterizations during the cold season. When taken separately by cold-season segment, ICS members perform quite well against observations at the NOAA BRW tower for the early cold season and ICT members perform well for the late cold season. The contrasting performance between site parameterizations is due to the topographic and hydrologic conditions, which are quite heterogeneous over a short distance and influence the plant communities and carbon storage, at each site. The ecosystems sampled by the ICS tower are seasonally inundated and retain a deep layer of organic soil that can be respired in greater amounts longer into the early cold season, whereas the well-drained hillslope at ICT does not allow for the accumulation of organic matter in the same way (Euskirchen et al., 2017; Larson et al., 2021). While varying topography and soil inundation throughout the North Slope means that each of these site relationships is likely to be representative of many other locations in the region with similar conditions, the early- to late-cold-season reduction in CO2 fluxes at these sites is not consistent with the observed regional atmospheric trend, and we remove the members parameterized by them from the ensemble. Individual eddy flux site parameterizations may reproduce regional CO2 fluxes for a given season, but it is important to consider their response to drivers across multiple seasons when scaling from the site-level to regional domains.

The observed cold-season CO2 flux pattern on the North Slope may be unique to tundra ecosystems of this region. For example, the CO2 fluxes from Natali and Watts et al. (2019) and Watts et al. (2021) both incorporate measurements from the North Slope. However, Natali and Watts et al. (2019) used boosted regression trees trained on belowground respiration measurements from across the pan-Arctic tundra and boreal zones, which may not be representative of our study region. The fluxes from Watts et al. (2021) are based on respiration measurements from throughout only Alaska and northwest Canada and conform better to local conditions. The evaluation of these CO2 fluxes against atmospheric CO2 measurements also produces results that are more consistent with our TVPRM ensemble determined by North Slope eddy flux tower measurements.

We find that the atmospheric observations are best matched by biogenic CO2 fluxes that include an additional CO2 source from tundra ecosystems during the zero-curtain period that are independent of Ts variability and year-round net CO2 fluxes from areas of inland water. The additional zero-curtain flux represents large-scale emission events that are not directly timed to microbial activity and root respiration controlled by Ts but that could be related to the delayed physical release of previously produced CO2 from soil through the snowpack as the soil layers remain unfrozen (Bowling and Massman, 2011). The Alaska North Slope also has many waterbodies distributed throughout the coastal tundra region, and the extent to which carbon cycles between small, shallow ponds and their surrounding terrestrial components is unclear (Magnússon et al., 2020). The biogenic CO2 fluxes in these areas are likely driven by ecosystem-scale CO2 fluxes from both coastal tundra and small ponds (Holgerson and Raymond, 2016; Tan et al., 2017), and their impact on the regional net CO2 flux, via both emissions and uptake, may be significant (Elder et al., 2018; Beckebanze et al., 2022). Only by adding fluxes that match observed zero-curtain CO2 emission pulses and by approximating net CO2 fluxes in aquatic areas can we reproduce the observed ΔCO2 magnitude in both the early and late cold season. The resulting seasonal change between the early and late cold season is consistent with the extended duration of the observed regional-scale zero curtain. The simplistic approximations suggested here are not inconsistent with the existing uncertainties in tundra CO2 flux modeling and demonstrate the importance of considering these additional CO2 fluxes and their mechanisms for future study.

4.3 Future state of net CO2 flux on the Alaska North Slope

As the Arctic warms rapidly, the competition between the growing and cold-season Arctic CO2 fluxes will determine the net biogenic CO2 flux into the atmosphere. Warming Ta warms soils, thaws permafrost, increases the active-layer thickness, and has extended the duration of the zero curtain from weeks to over 100 d (Romanovsky and Osterkamp, 2000; Schuur et al., 2015; Zona et al., 2016), all of which increase cold-season CO2 emissions. The warming may also increase net growing season uptake, but the severe light limitation at high northern latitudes limits the extent of the growing season, especially on the North Slope (Zhang et al., 2020). The future of CO2 fluxes from inland waters and wetlands in the Arctic is uncertain, but some studies suggest CO2 emissions from lakes may increase (Bayer et al., 2019). The culmination of these effects will likely push the North Slope into being a consistent net source in the future. However, observations at the NOAA BRW tower during our study period do not show elevated late-cold-season CO2 emissions; therefore, the North Slope was not a consistent net source through 2017. Accordingly, care must be taken to accurately represent CO2 fluxes from Arctic ecosystems during both the early and late cold season when calculating the annual net CO2 budget. TVPRM could be used with projections of meteorology and SIF to calculate the future net CO2 balance for this region, but we caution against overuse of the model using current parameters, as the flux–driver relationships in the rapidly warming Arctic ecosystems are changing so quickly that we would not assume accuracy into the future. While we can constrain the annual net CO2 budget with existing data, the Arctic is rapidly changing and needs constant monitoring. The following recommendations would provide more detailed spatial and seasonal constraints and up-to-date information on the processes driving CO2 fluxes across the region.

4.3.1 Future observation efforts

Our results motivate the need for a more extensive network of CO2 eddy flux towers operating year-round, alongside sensors for soil moisture and Ts profiles throughout the active layer to better understand the mechanisms driving year-round and especially early-cold-season CO2 fluxes. Noting that automated or semiautomated monitoring systems for aquatic environments currently do not exist for the North Slope or other high-latitude regions, this sensor network should be distributed throughout poorly sampled ecosystem types, particularly along wetness gradients that span mixed terrestrial–aquatic environments. The results of this study also support the need for additional continuous CO2 concentration measurements at tall towers across the North Slope (including away from the coast) to increase coverage of observed ΔCO2 during all seasons and to better constrain the regional background. Airborne measurements of both CO2 concentrations and CO2 fluxes remain valuable to sample areas less accessible via ground-based measurements, but a large-scale flight campaign in the region has not occurred since 2017. Any additional flights should be targeted as early before, and as late after, the growing season as possible. Satellites that rely on reflected sunlight to detect CO2 have increasingly been used to constrain CO2 budgets in the northern latitudes (e.g., Byrne et al., 2022), but data are very limited in the cold season, especially in far northern regions like the North Slope.

4.3.2 Future modeling efforts

The large initial range of potential regional net CO2 flux values that we found for the Alaska North Slope indicates a large sensitivity to the choices and assumptions made when scaling eddy flux observations from the site to the regional scale. The most important of these choices are the representation of the upland tundra, particularly for the response of Rsoil to Ts during the cold season, and the distribution of vegetation types throughout the domain. Future tundra CO2 modeling efforts should focus on using site-level data that are the most consistent with regional-scale fluxes, rather than incorporating data from all available sites. Consistency and accuracy in classification schemes used in vegetation maps must also be addressed. As we have shown with the atmospheric observations, not all model scenarios have equal likelihood to be true, and the mean of the model ensemble is not necessarily the most likely or most consistent with the atmosphere. Using these atmospheric observations is uncertain, however, due to potential errors in the transport modeling, which are difficult to quantify. Atmospheric modeling of remote areas such as the Alaska North Slope requires further evaluation and improvement. Moreover, increasing the model temporal resolution should be considered, as the importance of the zero-curtain and snow cover to the net CO2 flux of tundra ecosystems is recognized, both of which vary on the order of days and weeks, rather than months.

5 Conclusions

Observed atmospheric concentrations from aircraft and towers are a powerful tool that provide a regional constraint on the many combinations of possible CO2 flux parameterizations and distributions of tundra ecosystems on the North Slope of Alaska. We find that the annual regional net CO2 flux on the North Slope in not a consistent net source nor sink but instead varies between 6 and +6 Tg C yr−1 for 2012–2017. We can also identify ecosystem relationships and driver combinations that best represent both local CO2 flux patterns and regional atmospheric CO2 enhancements. The simulated regional net CO2 flux is highly sensitive to the assumptions made while scaling up eddy flux observations, especially the ecosystem response to Ts of tundra during the cold season and the spatial distribution of tundra types across the North Slope. Additionally, scaling methods that average observations from multiple eddy covariance flux sites should consider which sites are most representative of the regional impact of the biosphere on the atmosphere using integrative top-down observations.

This work shows that year-round measurements of atmospheric CO2 concentrations and fluxes across heterogeneous terrestrial and aquatic ecosystems are needed to represent the drivers of CO2 fluxes from Arctic regions. Arctic ecosystems have the potential to accelerate warming if vast stores of carbon are released or to buffer warming if increasing carbon uptake from vegetation occurs. All components of Arctic tundra ecosystems must be fully incorporated into Earth system models to improve projections of future climate warming and associated carbon cycle feedbacks.

Data availability

Data that support the findings of this study are listed below:


The supplement related to this article is available online at:

Author contributions

LDS and RC designed the study. KAA, ESE, JPG, AK, WCO, and DZ provided eddy covariance flux tower data. SCB, KM, and CS provided aircraft concentration data. JMH and MEM provided WRF-STILT particle files and footprints. YY provided RS-PM Ts data. JDW provided Watts et al. (2021) cold-season belowground CO2 fluxes. LDS developed and evaluated TVPRM net CO2 fluxes against observations. RC, EJLL, JWM, and JDW assisted the analysis. LDS wrote the paper. All co-authors contributed to the preparation of the manuscript.

Competing interests

The contact author has declared that none of the authors has any competing interests.


Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.


We would like to acknowledge that the Alaskan North Slope is home to multiple Alaska Native nations, including the Nunamiut, Gwich'in, Koyukuk, and Iñupiaq peoples. We support and honor the place-based knowledge of Indigenous Peoples and recognize their ancestral and contemporary stewardship of their homelands that we research. Luke D. Schiferl and Róisín Commane are supported by research funding from the Department of Earth and Environmental Sciences at Columbia University and the NASA ABoVE grant no. NNX17AC61A. Luke D. Schiferl is additionally supported by the National Science Foundation (NSF) Office of Polar Programs grant no. 1848620. Erik J. L. Larson and J. William Munger are supported by NASA ABoVE grant no. NNX17AE75G. Jennifer D. Watts is supported by NASA ABoVE grant no. 80NSSC19M0209 and NASA grant no. NNH17ZDA001N-NIP. Part of the research was carried out at the Jet Propulsion Laboratory, California Institute of Technology, under a contract with NASA (contract no. 80NM0018D0004). Imnavait Creek flux towers are funded under grants from the NSF Office of Polar Programs (grant nos. 1503912 and 0632264). Resources supporting John M. Henderson and WRF-STILT modeling were provided by NASA grant nos. NNX17AE75G and NNX17AC61A as well as by the NASA High-End Computing (HEC) program through the NASA Advanced Supercomputing (NAS) Division at Ames Research Center. We thank the R Project community for analysis and plotting tools, especially the ggplot2, ggpattern, magick, anytime, lubridate, raster, and cowplot packages. NCEP Reanalysis data were provided by the NOAA/OAR/ESRL PSL, Boulder, Colorado, USA. Some of the data products used in this paper were acquired for CARVE, a NASA Earth Ventures Suborbital (EV-S1) investigation.

Financial support

This research has been supported by the National Aeronautics and Space Administration (grant nos. NNX17AC61A, NNX17AE75G, 80NSSC19M0209, NNH17ZDA001N-NIP, and 80NM0018D0004) and the National Science Foundation, Office of Polar Programs (grant nos. 1848620, 1503912, and 0632264).

Review statement

This paper was edited by Paul Stoy and reviewed by two anonymous referees.


Arndt, K. A., Oechel, W. C., Goodrich, J. P., Bailey, B. A., Kalhori, A., Hashemi, J., Sweeney, C., and Zona, D.: Sensitivity of Methane Emissions to Later Soil Freezing in Arctic Tundra Ecosystems, J. Geophys. Res.-Biogeo., 124, 2595–2609,, 2019. 

Arndt, K. A., Lipson, D. A., Hashemi, J., Oechel, W. C., and Zona, D.: Snow melt stimulates ecosystem respiration in Arctic ecosystems, Glob. Change Biol., 26, 5042–5051,, 2020. 

Bayer, T. K., Gustafsson, E., Brakebusch, M., and Beer, C.: Future Carbon Emission From Boreal and Permafrost Lakes Are Sensitive to Catchment Organic Carbon Loads, J. Geophys. Res.-Biogeo., 124, 1827–1848,, 2019. 

Beckebanze, L., Rehder, Z., Holl, D., Wille, C., Mirbach, C., and Kutzbach, L.: Ignoring carbon emissions from thermokarst ponds results in overestimation of tundra net carbon uptake, Biogeosciences, 19, 1225–1244,, 2022. 

Biraud, S., Mei, F., Flynn, C., Hubbe, J., Long, C., Matthews, A., Pekour, M., Sedlacek, A., Springston, S., Tomlinson, J., and Chand, D.: Campaign datasets for ARM Airborne Carbon Measurements (ARM-ACME-V), [data set],, last access: 15 October 2019, 2016. 

Bowling, D. R. and Massman, W. J.: Persistent wind-induced enhancement of diffusive CO2 transport in a mountain forest snowpack, J. Geophys. Res.-Biogeo., 116, G04006,, 2011. 

Box, J. E., Colgan, W. T., Christensen, T. R., Schmidt, N. M., Lund, M., Parmentier, F.-J. W., Brown, R., Bhatt, U. S., Euskirchen, E. S., Romanovsky, V. E., Walsh, J. E., Overland, J. E., Wang, M., Corell, R. W., Meier, W. N., Wouters, B., Mernild, S., Mård, J., Pawlak, J., and Olsen, M. S.: Key indicators of Arctic climate change: 1971–2017, Environ. Res. Lett., 14, 045010,, 2019. 

Byrne, B., Liu, J., Yi, Y., Chatterjee, A., Basu, S., Cheng, R., Doughty, R., Chevallier, F., Bowman, K. W., Parazoo, N. C., Crisp, D., Li, X., Xiao, J., Sitch, S., Guenet, B., Deng, F., Johnson, M. S., Philip, S., McGuire, P. C., and Miller, C. E.: Multi-year observations reveal a larger than expected autumn respiration signal across northeast Eurasia, Biogeosciences, 19, 4779–4799,, 2022. 

CARVE Science Team: CARVE: L4 Gridded Particle Trajectories for WRF-STILT model, 2012–2016, ORNL DAAC [data set], Oak Ridge, Tennessee, USA,, 2017. 

CAVM Team: Circumpolar Arctic Vegetation Map (1:7,500,000 scale), Conservation of Arctic Flora and Fauna (CAFF) Map No. 1, U.S. Fish and Wildlife Service [data set], Anchorage, Alaska, ISBN: 0-9767525-0-6, ISBN-13: 978-0-9767525-0-9, 2003 (data available at:, last access: 5 August 2018), 2003. 

Chang, R. Y.-W., Miller, C. E., Dinardo, S. J., Karion, A., Sweeney, C., Daube, B. C., Henderson, J. M., Mountain, M. E., Eluszkiewicz, J., Miller, J. B., Bruhwiler, L. M. P., and Wofsy, S. C.: Methane emissions from Alaska in 2012 from CARVE airborne observations, P. Natl. Acad. Sci. USA, 111, 16694–16699,, 2014. 

Commane, R., Benmergui, J., Lindaas, J. O. W., Miller, S., Luus, K. A., Chang, R. Y.-W., Daube, B. C., Euskirchen, S., Henderson, J., Karion, A., Miller, J. B., Parazoo, N. C., Randerson, J. T., Sweeney, C., Tans, P., Thoning, K., Veraverbeke, S., Miller, C. E., and Wofsy, S. C.: CARVE: Net Ecosystem CO2 Exchange and Regional Carbon Budgets for Alaska, 2012–2014, ORNL DAAC [data set], Oak Ridge, Tennessee, USA,, 2017a. 

Commane, R., Lindaas, J., Benmergui, J., Luus, K. A., Chang, R. Y.-W., Daube, B. C., Euskirchen, E. S., Henderson, J. M., Karion, A., Miller, J. B., Miller, S. M., Parazoo, N. C., Randerson, J. T., Sweeney, C., Tans, P., Thoning, K., Veraverbeke, S., Miller, C. E., and Wofsy, S. C.: Carbon dioxide sources from Alaska driven by increasing early winter respiration from Arctic tundra, P. Natl. Acad. Sci. USA, 114, 5361–5366,, 2017b. 

Elder, C. D., Xu, X., Walker, J., Schnell, J. L., Hinkel, K. M., Townsend-Small, A., Arp, C. D., Pohlman, J. W., Gaglioti, B. V., and Czimczik, C. I.: Greenhouse gas emissions from diverse Arctic Alaskan lakes are dominated by young carbon, Nat. Clim. Change, 8, 166–171,, 2018. 

Euskirchen, E. S. and Edgar, C.: Arctic Observatory Network (AON) Data [data set],, last access: 15 August 2019 

Euskirchen, E. S., Bret-Harte, M. S., Scott, G. J., Edgar, C., and Shaver, G. R.: Seasonal patterns of carbon dioxide and water fluxes in three representative tundra ecosystems in northern Alaska, Ecosphere, 3, art4,, 2012. 

Euskirchen, E. S., Bret-Harte, M. S., Shaver, G. R., Edgar, C. W., and Romanovsky, V. E.: Long-Term Release of Carbon Dioxide from Arctic Tundra Ecosystems in Alaska, Ecosystems, 20, 960–974,, 2017. 

Fisher, J. B., Sikka, M., Oechel, W. C., Huntzinger, D. N., Melton, J. R., Koven, C. D., Ahlström, A., Arain, M. A., Baker, I., Chen, J. M., Ciais, P., Davidson, C., Dietze, M., El-Masri, B., Hayes, D., Huntingford, C., Jain, A. K., Levy, P. E., Lomas, M. R., Poulter, B., Price, D., Sahoo, A. K., Schaefer, K., Tian, H., Tomelleri, E., Verbeeck, H., Viovy, N., Wania, R., Zeng, N., and Miller, C. E.: Carbon cycle uncertainty in the Alaskan Arctic, Biogeosciences, 11, 4271–4288,, 2014. 

Goodrich, J. P., Oechel, W. C., Gioli, B., Moreaux, V., Murphy, P. C., Burba, G., and Zona, D.: Impact of different eddy covariance sensors, site set-up, and maintenance on the annual balance of CO2 and CH4 in the harsh Arctic environment, Agr. Forest Meteorol., 228–229, 239–251,, 2016. 

Hartery, S., Commane, R., Lindaas, J., Sweeney, C., Henderson, J., Mountain, M., Steiner, N., McDonald, K., Dinardo, S. J., Miller, C. E., Wofsy, S. C., and Chang, R. Y.-W.: Estimating regional-scale methane flux and budgets using CARVE aircraft measurements over Alaska, Atmos. Chem. Phys., 18, 185–202,, 2018. 

Henderson, J. M., Eluszkiewicz, J., Mountain, M. E., Nehrkorn, T., Chang, R. Y.-W., Karion, A., Miller, J. B., Sweeney, C., Steiner, N., Wofsy, S. C., and Miller, C. E.: Atmospheric transport simulations in support of the Carbon in Arctic Reservoirs Vulnerability Experiment (CARVE), Atmos. Chem. Phys., 15, 4093–4116,, 2015. 

Henderson, J., Miller, J. B., Nehrkorn, T., Chang, R. Y.-W., Sweeney, C., Steiner, N., Wofsy, S. C., and Miller, C. E.: CARVE: L4 Gridded Footprints from WRF-STILT model, 2012–2016, ORNL DAAC [data set], Oak Ridge, Tennessee, USA,, 2017. 

Henderson, J., Mountain, M., Dayalu, A., McKain, K., Hu, L., and Nehrkorn, T: ABoVE: Level-4 WRF-STILT Footprint Files for Circumpolar Receptors, 2016–2019, ORNL DAAC [data set], Oak Ridge, Tennessee, USA,, 2021a. 

Henderson, J., Mountain, M., Dayalu, A., McKain, K., Hu, L., and Nehrkorn, T: ABoVE: Level-4 WRF-STILT Particle Trajectories for Circumpolar Receptors, 2016–2019, ORNL DAAC [data set], Oak Ridge, Tennessee, USA,, 2021b. 

Hersbach, H., Bell, B., Berrisford, P., Hirahara, S., Horányi, A., Muñoz‐Sabater, J., Nicolas, J., Peubey, C., Radu, R., Schepers, D., Simmons, A., Soci, C., Abdalla, S., Abellan, X., Balsamo, G., Bechtold, P., Biavati, G., Bidlot, J., Bonavita, M., De Chiara, G., Dahlgren, P., Dee, D., Diamantakis, M., Dragani, R., Flemming, J., Forbes, R., Fuentes, M., Geer, A., Haimberger, L., Healy, S., Hogan, R.J., Hólm, E., Janisková, M., Keeley, S., Laloyaux, P., Lopez, P., Lupu, C., Radnoti, G., de Rosnay, P., Rozum, I., Vamborg, F., Villaume, S., and Thépaut, J.-N.: Complete ERA5 from 1979: Fifth generation of ECMWF atmospheric reanalyses of the global climate. Copernicus Climate Change Service (C3S) Data Store (CDS) [data set], (last access: 6 February 2019), 2017. 

Hersbach, H., Bell, B., Berrisford, P., Hirahara, S., Horányi, A., Muñoz-Sabater, J., Nicolas, J., Peubey, C., Radu, R., Schepers, D., Simmons, A., Soci, C., Abdalla, S., Abellan, X., Balsamo, G., Bechtold, P., Biavati, G., Bidlot, J., Bonavita, M., Chiara, G. D., Dahlgren, P., Dee, D., Diamantakis, M., Dragani, R., Flemming, J., Forbes, R., Fuentes, M., Geer, A., Haimberger, L., Healy, S., Hogan, R. J., Hólm, E., Janisková, M., Keeley, S., Laloyaux, P., Lopez, P., Lupu, C., Radnoti, G., de Rosnay, P., Rozum, I., Vamborg, F., Villaume, S., and Thépaut, J.-N.: The ERA5 global reanalysis, Q. J. Roy. Meteor. Soc., 146, 1999–2049,, 2020. 

Holgerson, M. A. and Raymond, P. A.: Large contribution to inland water CO2 and CH4 emissions from very small ponds, Nat. Geosci., 9, 222–226,, 2016. 

Hugelius, G., Strauss, J., Zubrzycki, S., Harden, J. W., Schuur, E. A. G., Ping, C.-L., Schirrmeister, L., Grosse, G., Michaelson, G. J., Koven, C. D., O'Donnell, J. A., Elberling, B., Mishra, U., Camill, P., Yu, Z., Palmtag, J., and Kuhry, P.: Estimated stocks of circumpolar permafrost carbon with quantified uncertainty ranges and identified data gaps, Biogeosciences, 11, 6573–6593,, 2014. 

Jeong, S.-J., Bloom, A. A., Schimel, D., Sweeney, C., Parazoo, N. C., Medvigy, D., Schaepman-Strub, G., Zheng, C., Schwalm, C. R., Huntzinger, D. N., Michalak, A. M., and Miller, C. E.: Accelerating rates of Arctic carbon cycling revealed by long-term atmospheric CO2 measurements, Sci. Adv., 4, eaao1167,, 2018. 

Joiner, J., Yoshida, Y., Guanter, L., and Middleton, E. M.: New methods for the retrieval of chlorophyll red fluorescence from hyperspectral satellite instruments: simulations and application to GOME-2 and SCIAMACHY, Atmos. Meas. Tech., 9, 3939–3967,, 2016. 

Joiner, J., Yoshida, Y., Koehler, P., Frankenberg, C., and Parazoo, N. C.: L2 Daily Solar-Induced Fluorescence (SIF) from MetOp-A GOME-2, 2007–2018, ORNL DAAC [data set], Oak Ridge, Tennessee, USA,, 2022. 

Karion, A., Sweeney, C., Kort, E. A., Shepson, P. B., Brewer, A., Cambaliza, M., Conley, S. A., Davis, K., Deng, A., Hardesty, M., Herndon, S. C., Lauvaux, T., Lavoie, T., Lyon, D., Newberger, T., Pétron, G., Rella, C., Smith, M., Wolter, S., Yacovitch, T. I., and Tans, P.: Aircraft-Based Estimate of Total Methane Emissions from the Barnett Shale Region, Environ. Sci. Technol., 49, 8124–8131,, 2015. 

Kim, J., Kim, Y., Zona, D., Oechel, W., Park, S.-J., Lee, B.-Y., Yi, Y., Erb, A., and Schaaf, C. L.: Carbon response of tundra ecosystems to advancing greenup and snowmelt in Alaska, Nat. Commun., 12, 6879,, 2021. 

Kim, Y., Kimball, J. S., Zhang, K., and McDonald, K. C.: Satellite detection of increasing Northern Hemisphere non-frozen seasons from 1979 to 2008: Implications for regional vegetation growth, Remote Sens. Environ., 121, 472–487,, 2012. 

Kittler, F., Eugster, W., Foken, T., Heimann, M., Kolle, O., and Göckede, M.: High-quality eddy-covariance CO2 budgets under cold climate conditions, J. Geophys. Res.-Biogeo., 122, 2064–2084,, 2017. 

Larson, E. J. L., Schiferl, L. D., Commane, R., Munger, J. W., Trugman, A. T., Ise, T., Euskirchen, E. S., Wofsy, S., and Moorcroft, P. M.: The changing carbon balance of tundra ecosystems: results from a vertically-resolved peatland biosphere model, Environ. Res. Lett., 17, 014019,, 2021. 

Li, X. and Xiao, J.: GOSIF – Global, OCO-2 based SIF product [data set],, last access: 24 April 2019. 

Luus, K. A. and Lin, J. C.: CARVE Modeled Gross Ecosystem CO2 Exchange and Respiration, Alaska, 2012–2014, ORNL DAAC [data set], Oak Ridge, Tennessee, USA,, 2017. 

Luus, K. A., Commane, R., Parazoo, N. C., Benmergui, J., Euskirchen, E. S., Frankenberg, C., Joiner, J., Lindaas, J., Miller, C. E., Oechel, W. C., Zona, D., Wofsy, S., and Lin, J. C.: Tundra photosynthesis captured by satellite-observed solar-induced chlorophyll fluorescence, Geophys. Res. Lett., 44, 2016GL070842,, 2017. 

Magney, T. S., Bowling, D. R., Logan, B. A., Grossmann, K., Stutz, J., Blanken, P. D., Burns, S. P., Cheng, R., Garcia, M. A., K?hler, P., Lopez, S., Parazoo, N. C., Raczka, B., Schimel, D., and Frankenberg, C.: Mechanistic evidence for tracking the seasonality of photosynthesis with solar-induced fluorescence, P. Natl. Acad. Sci. USA, 116, 11640–11645,, 2019. 

Magnússon, R. Í., Limpens, J., van Huissteden, J., Kleijn, D., Maximov, T. C., Rotbarth, R., Sass-Klaassen, U., and Heijmans, M. M. P. D.: Rapid Vegetation Succession and Coupled Permafrost Dynamics in Arctic Thaw Ponds in the Siberian Lowland Tundra, J. Geophys. Res.-Biogeo., 125, 2019JG005618,, 2020. 

Mbufong, H. N., Lund, M., Aurela, M., Christensen, T. R., Eugster, W., Friborg, T., Hansen, B. U., Humphreys, E. R., Jackowicz-Korczynski, M., Kutzbach, L., Lafleur, P. M., Oechel, W. C., Parmentier, F. J. W., Rasse, D. P., Rocha, A. V., Sachs, T., van der Molen, M. K., and Tamstorf, M. P.: Assessing the spatial variability in peak season CO2 exchange characteristics across the Arctic tundra using a light response curve parameterization, Biogeosciences, 11, 4897–4912,, 2014. 

Meredith, M., Sommerkorn, M., Cassotta, S., Derksen, C., Ekaykin, A., Hollowed, A., Kofinas, G., Mackintosh, A., Melbourne-Thomas, J., Muelbert, M. M. C., Ottersen, G., Pritchard, H., and Schuur, E. A. G.: Polar Regions, in: IPCC Special Report on the Ocean and Cryosphere in a Changing Climate, edited by: Pörtner, H.-O., Roberts, D. C., Masson-Delmotte, V., Zhai, P., Tignor, M., Poloczanska, E., Mintenbeck, K., Alegría, A., Nicolai, M., Okem, A., Petzold, J., Rama, B., and Weyer, N. M., Cambridge University Press, Cambridge, UK and New York, NY, USA,, 2019. 

Mesinger, F., DiMego, G., Kalnay, E., Mitchell, K., Shafran, P. C., Ebisuzaki, W., Jović, D., Woollen, J., Rogers, E., Berbery, E. H., Ek, M. B., Fan, Y., Grumbine, R., Higgins, W., Li, H., Lin, Y., Manikin, G., Parrish, D., and Shi, W.: North American Regional Reanalysis, B. Am. Meteorol. Soc., 87, 343–360,, 2006. 

Miller, S. M., Miller, C. E., Commane, R., Chang, R. Y.-W., Dinardo, S. J., Henderson, J. M., Karion, A., Lindaas, J., Melton, J. R., Miller, J. B., Sweeney, C., Wofsy, S. C., and Michalak, A. M.: A multiyear estimate of methane fluxes in Alaska from CARVE atmospheric observations, Global Biogeochem. Cy., 30, 1441–1453,, 2016. 

Natali, S. M., Watts, J. D., Rogers, B. M., Potter, S., Ludwig, S. M., Selbmann, A.-K., Sullivan, P. F., Abbott, B. W., Arndt, K. A., Birch, L., Björkman, M. P., Bloom, A. A., Celis, G., Christensen, T. R., Christiansen, C. T., Commane, R., Cooper, E. J., Crill, P., Czimczik, C., Davydov, S., Du, J., Egan, J. E., Elberling, B., Euskirchen, E. S., Friborg, T., Genet, H., Göckede, M., Goodrich, J. P., Grogan, P., Helbig, M., Jafarov, E. E., Jastrow, J. D., Kalhori, A. A. M., Kim, Y., Kimball, J. S., Kutzbach, L., Lara, M. J., Larsen, K. S., Lee, B.-Y., Liu, Z., Loranty, M. M., Lund, M., Lupascu, M., Madani, N., Malhotra, A., Matamala, R., McFarland, J., McGuire, A. D., Michelsen, A., Minions, C., Oechel, W. C., Olefeldt, D., Parmentier, F.-J. W., Pirk, N., Poulter, B., Quinton, W., Rezanezhad, F., Risk, D., Sachs, T., Schaefer, K., Schmidt, N. M., Schuur, E. A. G., Semenchuk, P. R., Shaver, G., Sonnentag, O., Starr, G., Treat, C. C., Waldrop, M. P., Wang, Y., Welker, J., Wille, C., Xu, X., Zhang, Z., Zhuang, Q., and Zona, D.: Large loss of CO2 in winter observed across the northern permafrost region, Nat. Clim. Change, 9, 852–857,, 2019. 

NOAA PSL: NCEP North American Regional Reanalysis (NARR) data provided by the NOAA PSL [data set], Boulder, Colorado, USA,, last access: 26 April 2018. 

Oechel, W. C., Laskowski, C. A., Burba, G., Gioli, B., and Kalhori, A. A. M.: Annual patterns and budget of CO2 flux in an Arctic tussock tundra ecosystem, J. Geophys. Res.-Biogeo., 119, 323–339,, 2014. 

Outcalt, S. I., Nelson, F. E., and Hinkel, K. M.: The zero-curtain effect: Heat and mass transfer across an isothermal region in freezing soil, Water Resour. Res., 26, 1509–1516,, 1990. 

Pallandt, M. M. T. A., Kumar, J., Mauritz, M., Schuur, E. A. G., Virkkala, A.-M., Celis, G., Hoffman, F. M., and Göckede, M.: Representativeness assessment of the pan-Arctic eddy covariance site network and optimized future enhancements, Biogeosciences, 19, 559–583,, 2022. 

Patankar, R., Mortazavi, B., Oberbauer, S. F., and Starr, G.: Diurnal patterns of gas-exchange and metabolic pools in tundra plants during three phases of the arctic growing season, Ecol. Evol., 3, 375–388,, 2013. 

Porcar-Castell, A., Tyystjärvi, E., Atherton, J., van der Tol, C., Flexas, J., Pfündel, E. E., Moreno, J., Frankenberg, C., and Berry, J. A.: Linking chlorophyll a fluorescence to photosynthesis for remote sensing applications: mechanisms and challenges, J. Exp. Bot., 65, 4065–4095,, 2014. 

Raynolds, M. and Walker, D.: Raster Circumpolar Arctic Vegetation Map, V1, Mendeley Data [data set],, 2019. 

Raynolds, M. K., Walker, D. A., Balser, A., Bay, C., Campbell, M., Cherosov, M. M., Daniëls, F. J. A., Eidesen, P. B., Ermokhina, K. A., Frost, G. V., Jedrzejek, B., Jorgenson, M. T., Kennedy, B. E., Kholod, S. S., Lavrinenko, I. A., Lavrinenko, O. V., Magnússon, B., Matveyeva, N. V., Metúsalemsson, S., Nilsen, L., Olthof, I., Pospelov, I. N., Pospelova, E. B., Pouliot, D., Razzhivin, V., Schaepman-Strub, G., Šibík, J., Telyatnikov, M. Yu., and Troeva, E.: A raster version of the Circumpolar Arctic Vegetation Map (CAVM), Remote Sens. Environ., 232, 111297,, 2019. 

Romanovsky, V. E. and Osterkamp, T. E.: Effects of unfrozen water on heat and mass transport processes in the active layer and permafrost, Permafrost Periglac., 11, 219–239,<219::AID-PPP352>3.0.CO;2-7, 2000. 

Schiferl, L. and Commane, R.: ABoVE: TVPRM Simulated Net Ecosystem Exchange, Alaskan North Slope, 2008–2017, ORNL DAAC [data set], Oak Ridge, Tennessee, USA,, 2022. 

Schuur, E. A. G., McGuire, A. D., Schädel, C., Grosse, G., Harden, J. W., Hayes, D. J., Hugelius, G., Koven, C. D., Kuhry, P., Lawrence, D. M., Natali, S. M., Olefeldt, D., Romanovsky, V. E., Schaefer, K., Turetsky, M. R., Treat, C. C., and Vonk, J. E.: Climate change and the permafrost carbon feedback, Nature, 520, 171–179,, 2015. 

Sun, Y., Frankenberg, C., Wood, J. D., Schimel, D. S., Jung, M., Guanter, L., Drewry, D. T., Verma, M., Porcar-Castell, A., Griffis, T. J., Gu, L., Magney, T. S., Köhler, P., Evans, B., and Yuen, K.: OCO-2 advances photosynthesis observation from space via solar-induced chlorophyll fluorescence, Science, 358, eaam5747,, 2017. 

Sweeney, C. and McKain, K.: ABoVE: Atmospheric Profiles of CO, CO2 and CH4 Concentrations from Arctic-CAP, 2017, ORNL DAAC [data set], Oak Ridge, Tennessee, USA,, 2019. 

Sweeney, C., Dlugokencky, E., Miller, C. E., Wofsy, S., Karion, A., Dinardo, S., Chang, R. Y.-W., Miller, J. B., Bruhwiler, L., Crotwell, A. M., Newberger, T., McKain, K., Stone, R. S., Wolter, S. E., Lang, P. E., and Tans, P.: No significant increase in long-term CH4 emissions on North Slope of Alaska despite significant increase in air temperature, Geophys. Res. Lett., 43, 6604–6611,, 2016. 

Sweeney, C., Chatterjee, A., Wolter, S., McKain, K., Bogue, R., Conley, S., Newberger, T., Hu, L., Ott, L., Poulter, B., Schiferl, L., Weir, B., Zhang, Z., and Miller, C. E.: Using atmospheric trace gas vertical profiles to evaluate model fluxes: a case study of Arctic-CAP observations and GEOS simulations for the ABoVE domain, Atmos. Chem. Phys., 22, 6347–6364,, 2022. 

Tadić, J. M., Miller, S., Yadav, V., and Biraud, S. C.: Greenhouse gas fluxes from Alaska's North Slope inferred from the Airborne Carbon Measurements Campaign (ACME-V), Atmos. Environ., 248, 118239,, 2021. 

Tan, Z., Zhuang, Q., Shurpali, N. J., Marushchak, M. E., Biasi, C., Eugster, W., and Anthony, K. W.: Modeling CO2 emissions from Arctic lakes: Model development and site-level study, J. Adv. Model. Earth Sy., 9, 2190–2213,, 2017. 

Tao, J., Zhu, Q., Riley, W. J., and Neumann, R. B.: Warm-season net CO2 uptake outweighs cold-season emissions over Alaskan North Slope tundra under current and RCP8.5 climate, Environ. Res. Lett., 16, 055012,, 2021. 

Thoning, K. W., Kitzis, D. R., and Crotwell, A.: Atmospheric Carbon Dioxide Dry Air Mole Fractions from quasi-continuous measurements at Barrow, Alaska, 1973–2017, Version 2018-10, National Oceanic and Atmospheric Administration (NOAA), Earth System Research Laboratory (ESRL) [data set], Boulder, Colorado, USA,, last access: 16 November 2018. 

Walker, D. A., Raynolds, M. K., Daniëls, F. J. A., Einarsson, E., Elvebakk, A., Gould, W. A., Katenin, A. E., Kholod, S. S., Markon, C. J., Melnikov, E. S., Moskalenko, N. G., Talbot, S. S., Yurtsev, B. A., and Team, T. other members of the C.: The Circumpolar Arctic vegetation map, J. Veg. Sci., 16, 267–282,, 2005. 

Wang, J. A., Sulla-Menashe, D., Woodcock, C. E., Sonnentag, O., Keeling, R. F., and Friedl, M. A.: ABoVE: Landsat-derived Annual Dominant Land Cover Across ABoVE Core Domain, 1984–2014, ORNL DAAC [data set], Oak Ridge, Tennessee, USA,, 2019. 

Wang, J. A., Sulla-Menashe, D., Woodcock, C. E., Sonnentag, O., Keeling, R. F., and Friedl, M. A.: Extensive land cover change across Arctic–Boreal Northwestern North America from disturbance and climate forcing, Glob. Change Biol., 26, 807–822,, 2020. 

Watts, J. D., Natali, S., Potter, S., and Rogers, B. M.: Gridded Winter Soil CO2 Flux Estimates for pan-Arctic and Boreal Regions, 2003–2100, ORNL DAAC [data set], Oak Ridge, Tennessee, USA,, 2019. 

Watts, J. D., Natali, S. M., Minions, C., Risk, D., Arndt, K., Zona, D., Euskirchen, E. S., Rocha, A. V., Sonnentag, O., Helbig, M., Kalhori, A., Oechel, W., Ikawa, H., Ueyama, M., Suzuki, R., Kobayashi, H., Celis, G., Schuur, E. A. G., Humphreys, E., Kim, Y., Lee, B.-Y., Goetz, S., Madani, N., Schiferl, L. D., Commane, R., Kimball, J. S., Liu, Z., Torn, M. S., Potter, S., Wang, J. A., Jorgenson, M. T., Xiao, J., Li, X., and Edgar, C.: Soil respiration strongly offsets carbon uptake in Alaska and Northwest Canada, Environ. Res. Lett., 16, 084051,, 2021. 

Watts, J. D., Natali, S., and Minions, C.: Soil Respiration Maps for the ABoVE Domain, 2016–2017, ORNL DAAC [data set], Oak Ridge, Tennessee, USA,, 2022. 

Worthy, D. E. J., Chan, E., Ishizawa, M., Chan, D., Poss, C., Dlugokencky, E. J., Maksyutov, S., and Levin, I.: Decreasing anthropogenic methane emissions in Europe and Siberia inferred from continuous carbon dioxide and methane observations at Alert, Canada, J. Geophys. Res., 114, D10301,, 2009.  

Yang, X., Tang, J., Mustard, J. F., Lee, J.-E., Rossini, M., Joiner, J., Munger, J. W., Kornfeld, A., and Richardson, A. D.: Solar-induced chlorophyll fluorescence that correlates with canopy photosynthesis on diurnal and seasonal scales in a temperate deciduous forest, Geophys. Res. Lett., 42, 2977–2987,, 2015. 

Yi, Y., Kimball, J. S., Chen, R. H., Moghaddam, M., Reichle, R. H., Mishra, U., Zona, D., and Oechel, W. C.: Characterizing permafrost active layer dynamics and sensitivity to landscape spatial heterogeneity in Alaska, The Cryosphere, 12, 145–161,, 2018. 

Yi, Y., Kimball, J. S., Chen, R. H., Moghaddam, M., and Miller, C. E.: Sensitivity of active-layer freezing process to snow cover in Arctic Alaska, The Cryosphere, 13, 197–218,, 2019. 

Zhang, Y.: CSIF, figshare [data set],, 2018. 

Zhang, Y., Joiner, J., Alemohammad, S. H., Zhou, S., and Gentine, P.: A global spatially contiguous solar-induced fluorescence (CSIF) dataset using neural networks, Biogeosciences, 15, 5779–5800,, 2018. 

Zhang, Y., Commane, R., Zhou, S., Williams, A. P., and Gentine, P.: Light limitation regulates the response of autumn terrestrial carbon uptake to warming, Nat. Clim. Change, 10, 739–743,, 2020. 

Zona, D.: Greenhouse gas flux measurements at the zero curtain, North Slope, Alaska, 2012–2019, Arctic Data Center [data set],, 2019. 

Zona, D., Gioli, B., Commane, R., Lindaas, J., Wofsy, S. C., Miller, C. E., Dinardo, S. J., Dengel, S., Sweeney, C., Karion, A., Chang, R. Y.-W., Henderson, J. M., Murphy, P. C., Goodrich, J. P., Moreaux, V., Liljedahl, A., Watts, J. D., Kimball, J. S., Lipson, D. A., and Oechel, W. C.: Cold season emissions dominate the Arctic tundra methane budget, P. Natl. Acad. Sci. USA, 113, 40–45,, 2016. 

Schiferl and colleagues combine top-down and bottom-up estimates of carbon dioxide emissions from the North Slope of Alaska to find that CO2 efflux from terrestrial and aquatic ecosystems during the early cold season (September – December) are critical for explaining its carbon balance. Fluxes during the late cold season (January through April) were muted in comparison. Despite the importance of cold-season efflux, growing season CO2 uptake and release processes dominated its interannual variability. This study helps clarify how the carbon cycle of complex tundra ecosystems across large Arctic regions respond to ongoing climate changes.
Short summary
As the Arctic rapidly warms, vast stores of thawing permafrost could release carbon dioxide (CO2) into the atmosphere. We combined observations of atmospheric CO2 concentrations from aircraft and a tower with observed CO2 fluxes from tundra ecosystems and found that the Alaskan North Slope in not a consistent source nor sink of CO2. Our study shows the importance of using both site-level and atmospheric measurements to constrain regional net CO2 fluxes and improve biogenic processes in models.
Final-revised paper