Radiocarbon measurements of ecosystem respiration and soil pore-space CO2 in Utqiaġvik (Barrow), Alaska

Abstract. Radiocarbon measurements of ecosystem respiration and soil pore space
CO2 are useful for determining the sources of ecosystem respiration,
identifying environmental controls on soil carbon cycling rates, and
parameterizing and evaluating models of the carbon cycle. We measured flux
rates and radiocarbon content of
ecosystem respiration, as well as radiocarbon in soil profile CO2 in
Utqiaġvik (Barrow), Alaska, during the summers of 2012, 2013, and 2014.
We found that radiocarbon in ecosystem respiration
(Δ14CReco) ranged from +60.5 to −160 ‰ with a
median value of +23.3 ‰. Ecosystem respiration became more
depleted in radiocarbon from summer to autumn, indicating increased
decomposition of old soil organic carbon and/or decreased CO2
production from fast-cycling carbon pools. Across permafrost features,
ecosystem respiration from high-centered polygons was depleted in radiocarbon
relative to other polygon types. Radiocarbon content in soil pore-space
CO2 varied between −7.1 and −280 ‰, becoming more
negative with depth in individual soil profiles. These pore-space radiocarbon
values correspond to CO2 mean ages of 410 to 3350 years, based on a
steady-state, one-pool model. Together, these data indicate that
deep soil
respiration was derived primarily from old, slow-cycling carbon, but
that total CO2 fluxes depended largely on autotrophic respiration and
heterotrophic decomposition of fast-cycling carbon within the shallowest soil
layers. The relative contributions of these different CO2 sources
were highly variable across
microtopographic features and time in the sampling season. The highly
negative Δ14C values in soil pore-space CO2 and autumn
ecosystem respiration indicate that when it is not frozen, very old soil
carbon is vulnerable to decomposition. Radiocarbon data and associated
CO2 flux and temperature data are stored in the data repository of
the Next Generation Ecosystem Experiments (NGEE-Arctic) at
http://dx.doi.org/10.5440/1364062 and
https://doi.org/10.5440/1418853.



Introduction
The flux of CO 2 from ecosystems to the atmosphere is a critical component of the global carbon budget.This flux is highly heterogeneous in space and time, so extensive datasets are needed to evaluate the carbon balance within ecosystems.Many measurements have been made of soil surface CO 2 emissions using soil chambers and eddy covariance towers (e.g., Baldocchi, 2008;Bond-Lamberty and Thomson, 2010;Davidson et al., 2002;Norman et al., 1997;Xu and Baldocchi, 2004).Such measurements reveal spatial and temporal patterns in soil and ecosystem respiration that are important for scaling soil carbon emissions across the landscape or identifying drivers of respiration rates (Wainwright et al., 2015).Bulk fluxes of soil or ecosystem respiration, however, do not provide information on the cycling rates or relative contributions of different soil carbon pools to total CO 2 emissions.CO 2 emitted from the soil surface includes autotrophic respiration of rapidly cycling carbon as well as heterotrophic decomposition of carbon that cycles on a broad range of timescales.Shifts in these carbon pool distributions can have large long-term consequences for soil carbon stocks, but may be impossible to detect in bulk CO 2 flux rates (Hopkins et al., 2012;Schuur et al., 2009;Torn et al., 2009;Trumbore, 2000Trumbore, , 2009)).
Linking CO 2 emissions with carbon cycling rates requires a tracer of carbon dynamics that can differentiate between source carbon pools.Natural abundance radiocarbon provides such a tracer, as the radiocarbon content of CO 2 reflects the age (since photosynthesis) and decomposition rates of its component sources (Trumbore, 2000).Accordingly, radiocarbon in surface CO 2 emissions and soil pore gas can be used to quantify carbon cycling rates and assess their variability across space, time, carbon pools, and environmental factors such as soil moisture, temperature, and vegetation (Gaudinski et al., 2000;Trumbore, 2000).Within and across sites, variations in the radiocarbon abundance of respired CO 2 can indicate differences in substrate utilization by microbial decomposers (Borken et al., 2006;Chasar et al., 2000) and shifts between more fast-cycling and slow-cycling substrates (Hicks Pries et al., 2013;Hopkins et al., 2012;Schuur et al., 2009).Radiocarbon measurements of soil respiration may be particularly useful for models of the carbon cycle.Carbon pool structures and decomposition dynamics inferred from radiocarbon data may be used to parameterize pool-specific respiration rates; radiocarbon measurements in surface CO 2 emissions and soil pore gas can be used to benchmark model performance; and the radiocarbon signature of respiration can be used to constrain the terrestrial signal in top-down carbon cycle analyses (He et al., 2016;Randerson et al., 2002).
In the high latitudes in particular, increased decomposition rates from large, slow-cycling soil carbon stores have the potential to generate an important long-term climate change feedback (Schuur et al., 2015).There, an estimated 1300 Pg of soil carbon has been protected from decomposition by cold temperatures and often frozen or anoxic conditions (Hugelius et al., 2014).With climate change, these controls on decomposition rates are expected to change as soils warm, hydrological changes occur, and permafrost degradation intensifies (ACIA, 2004).In the lower depths of the seasonally thawed active layer, carbon that has historically cycled on millennial timescales may be highly decomposable under aerobic and thawed conditions (Elberling et al., 2013;Mueller et al., 2015;Strauss et al., 2015;Waldrop et al., 2010).Our current understanding of this old carbon's decomposability, however, is based primarily on CO 2 production rates from laboratory incubations, which may not accurately reflect in situ decomposition dynamics.To quantify in situ decomposition rates, field radiocarbon measurements can be used to differentiate between slowcycling and fast-cycling carbon and link decomposition dynamics to their environmental controls.
In spite of their use in tracing carbon cycle dynamics, only a limited set of radiocarbon data have been published for Arctic ecosystem respiration (e.g., Hardie et al., 2009;Hicks Pries et al., 2013;Lupascu et al., 2014a;Nowinski et al., 2010;Phillips et al., 2015;Schuur et al., 2009) or soil pore-space CO 2 (Czimczik and Welker, 2010;Lupascu et al., 2014a, b), from only a few sites and tundra types.To our knowledge, no data are currently available on radiocarbon in CO 2 from Arctic polygon tundra.Polygon tundra, which covers an estimated 250 000 km 2 throughout the Arctic (Donner et al., 2007), is characterized by microtopographic features ∼10-32 m in diameter known as ice wedge polygons, which form by the growth and subsequent degradation of subsurface ice formations (ice wedges) (Billings and Peterson, 1980;Brown et al., 1980).Polygons are separated by low-lying drainage channels called troughs, and can be classified according to the elevation of their centers into low-centered, flat-centered, and high-centered categories.Lowcentered polygons (LC) have raised, dry rims at their perimeter surrounding low, saturated or inundated centers; flat-centered polygons (FC) have flat, even surfaces with minimal rim delineation; and high-centered polygons (HC) have raised centers with dry surface soils that slope steeply into troughs (Fig. 1).Among LC, FC, and HC polygons, differences in active layer depth, oxygen availability, soil carbon distribution, and thermal conductivity (Liljedahl et al., 2016;Lipson et al., 2012;Ping et al., 1998) create strong differences in surface carbon fluxes (Vaughn et al., 2016;Wainwright et al., 2015).Given these steep spatial gradients in carbon cycling rates and their environmental controls, radiocarbon in respiration from polygon tundra ecosystems offers unique information on subsurface carbon dynamics across space, environmental factors, and the progression of seasonal thaw.
Here, we present measurements of radiocarbon in ecosystem respiration, soil pore-space CO 2 , and soil organic carbon from a polygon tundra ecosystem in Utqiaġvik, Alaska.We measured rates and radiocarbon contents of soil surface CO 2 emissions across three summer seasons, from a range of ice wedge polygon features.Using soil pore-space 14 CO 2 profiles, we quantified vertical distributions of carbon cycling rates within the soil profile.Together, these measurements reveal temporal patterns and spatial relationships in ecosystem respiration and its source carbon pools.

Study site
Field sampling was conducted at the Barrow Environmental Observatory (BEO), ∼6 km east of Utqiaġvik (formerly Barrow), Alaska (71.3 ∘ N, 156.5 ∘ W), at the northern end of the Alaskan Arctic coastal plain.Utqiaġvik has a mean annual temperature of −9.8 ∘ C and mean annual precipitation of 146 mm, with long, dry winters and short, moist, cool summers (NOAA, 2018).The land surface has low topographic relief reaching a maximum elevation of 5 m (Brown et al., 1980;Hubbard et al., 2013).The seasonally thawed active layer ranges in depth from 20 to 60 cm, underlain by continuous ice-rich permafrost to depths greater than 400 m (Hinkel and Nelson, 2003).Formed from the late Pleistocene Gubic formation (Black, 1964), soils in the region are dominated by Typic Aquiturbels (53 %), Typic Histoturbels (22 %), and Typic Aquorthels (8.6 %) (Bockheim et al., 1999).Across the Utqiaġvik peninsula, 24 % of the ground surface is covered by LC polygons, 17 % by FC polygons, 11 % by HC polygons, and 5 % by their surrounding troughs (Lara et al., 2014).Within the specific study region, polygons cover the majority of the land surface, with the remainder comprised of a polygon-free drainage region.In the study area, FC polygons occupy the greatest percentage of the total ground surface, HC polygons occupy the smallest percentage, and LC polygons have intermediate percentage cover (Wainwright et al., 2015).Sample collection was distributed across centers, rims, and troughs of the three polygon types to capture variations in polygon microtopography and associated carbon cycle controls.
3 Field measurements and sample collection

Surface CO 2 emissions
In August and October 2012, July and September 2013, and September 2014, we collected soil surface CO 2 emissions for radiocarbon analysis using a static chamber method modified from Hahn et al. (2006).Samples were collected from a total of 19 locations within 11 polygons (4 LC, 4 FC, 3 HC).Opaque chambers (25 cm diameter) were seated on circular PVC chamber bases extending to depths of ∼10 cm below the soil surface, with aboveground chamber height varying between ∼15 and 20 cm.We installed all bases at least 2 days prior to sampling to limit the influence of disturbance on gas flux rates and radiocarbon values.Chambers blocked light transmission and were tall enough to enclose surface vegetation, so CO 2 emissions were equivalent to ecosystem respiration.For sample collection, the chamber was placed in a 3 cm deep channel on the top rim of each base, which was filled with water to create an airtight seal.We then circulated chamber gas through soda lime for 20 min at a flow rate of 1 L min −1 to remove ambient CO 2 .CO 2 was allowed to accumulate in the chamber over 2 to 48 h, depending on the rate of CO 2 accumulation, which we monitored periodically by injecting a 30 mL sample of chamber gas into a LI-820 CO 2 gas analyzer (LI-COR) at a flow rate of ∼1 L min −1 .For all samples, the final chamber CO 2 concentration was more than twice its initial concentration, important for accuracy in chamber radiocarbon measurements (Egan et al., 2014).Based on this concentration measurement, a volume of chamber gas sufficient for radiocarbon analysis was collected in one or more 500-1000 mL evacuated stainless steel canisters connected by capillary tubing to a chamber sampling port.High-concentration samples were collected with a syringe and needle through a septum in the sampling port and immediately injected into evacuated glass vials sealed with 14 mm thick chlorobutyl septa (Bellco Glass, Inc.).To correct chamber gas samples for atmospheric contamination, we collected local air samples in 3000 mL stainless steel canisters on 12 August 2012, 13 July 2013, and 2 and 7 September 2014.
In July and September 2013 and September 2014, we measured rates of ecosystem respiration from radiocarbon sampling locations.CO 2 fluxes were measured within 2 days of radiocarbon sample collection, using opaque static chambers seated on bases described above and vented according to Xu et al. (2006) to minimize pressure changes due to the Venturi effect.For each measurement, we measured CO 2 concentrations within the chamber over a period of 4-8 min using a Los Gatos Research, Inc. Portable Greenhouse Gas Analyzer.We calculated the CO 2 flux rate (equivalent to ecosystem respiration) as the slope of the linear portion of its concentration vs. time curve, converted to units of µmol m −2 s −2 according to chamber volume and temperature.Endpoints of this linear region were determined manually for each curve, and curves lacking a clear linear range were not included in the dataset.Corresponding to each CO 2 flux measurement, we measured thaw depth with a tile probe.

Soil pore-space CO 2
In August 2012 and July 2013, we collected soil pore gas from a total of six soil profiles within five polygons (one LC, two FC, two HC).Samples were collected with a method similar to Czimczik and Welker (2010).Briefly, 1/4'' diameter stainless steel probes were inserted into the soil at 45 ∘ angles to vertical depths of 10, 20, and 30 cm, or to 2 cm above the frost table if thaw depth was less than 30 cm.Probes were capped with gastight septa and allowed to remain in place throughout the sampling season.Before collecting each sample, we purged 10 mL of gas from the probe and measured the CO 2 concentration with a LI-820 CO 2 gas analyzer as described above.Radiocarbon samples were collected by connecting evacuated 500-1000 mL stainless steel canisters to probes via flow-restricting tubing (Upchurch scientific, 0.01'' ID ×10 cm length), which allowed canisters to fill slowly over 4 h with minimal disturbance to the soil CO 2 concentration gradient (Gaudinski et al., 2000).In-line Drierite water traps were used during sample collection to prevent moisture accumulation in canisters.As with surface respiration samples, high-concentration samples were collected from probes with a syringe and needle and immediately injected into evacuated glass vials.Due to water-saturated soils or clogged soil probes, we were unable to obtain samples from all profiles and depths.For this reason, the final sample set represents only a subset of depths and sampling locations and is not a random sample of the entire landscape.

Soil organic matter
On 8 September 2014, we collected 1'' diameter soil cores from the centers of the nine polygons used for September 2014 soil surface CO 2 sampling (3 LC, 3 FC, 3 HC).Soil cores were collected from within 0.5 m of each static chamber.Cores were collected using a manual soil recovery probe (AMS Samplers Inc.) to the base of the thawed soil layer (16-44 cm depth) and stored intact for 8 days at 5 ∘ C before further processing.Following this 8-day period, cores from FC and HC polygons were divided into 2 equal depth increments, which ranged from 8 to 15.5 cm in length.In the case of LC polygons, where standing water was present at the time of sampling, the shallow and deep increments were defined as surface water and sediment, respectively.For a separate experiment, core increments were incubated at 3 ∘ C for 387 days, during which CO 2 production rates were monitored with a LI-820 CO 2 analyzer (LI-COR).Following this 387-day period, we collected evolved CO 2 for radiocarbon analysis using stainless steel air sampling canisters or glass serum vials.Incubated soil, water, and sediment samples were freeze-dried and ground for carbon content and radiocarbon analyses.

Sample purification and radiocarbon analysis
CO 2 from gas samples was cryogenically purified under vacuum, divided for 14 C and 13 C analysis, and sealed in 9 mm quartz tubes.For radiocarbon analysis, we sent samples to Lawrence Livermore National Laboratory's Center for Accelerator Mass Spectrometry (CAMS) or the Carbon, Water, and Soils Research Lab at the USDA-FS Northern Research Station, where CO 2 was reduced to graphite on iron powder under H 2 . 14C abundance was then measured at CAMS using an HVEC FN Tandem Van de Graaff accelerator mass spectrometer or at UC Irvine's Keck Carbon Cycle AMS facility. 13C∕ 12 C in CO 2 splits was analyzed on the UC Davis Stable Isotope Laboratory GVI Optima Stable Isotope Ratio Mass Spectrometer.
Two quality control measures were employed to test for contamination by non-sample carbon during gas sample storage and preparation.First, we evaluated the air sampling canisters with repeated leak testing and an ethanol-derived certified CO 2 reference gas standard with empirically determined 14 C (Airgas, USA) to confirm that no contamination was introduced during sample storage.Second, the vacuum line used for CO 2 purification was tested with USGS coal and oxalic acid II (NIST SRM 4990C) standards for contamination or leakage.Standards were well characterized with respect to 14 C. Subsamples of bulk soil from the 2014 incubation were sent to the Carbon, Water, and Soils Research Lab, where they were combusted to CO 2 and analyzed for radiocarbon as above.Carbon concentrations and 13 C abundances were measured on a Costech analytical Technologies ECS 4010 elemental analyzer coupled to a Thermo Fischer Scientific Delta V plus isotope ratio mass spectrometer.Given the acidic pH range of soils at this site (Zona et al., 2011), we consider bulk soil carbon to be equivalent to soil organic carbon.
Following the conventions of Stuiver and Polach (1977), radiocarbon results are presented as fraction modern relative to the NBS Oxalic Acid I (OX1) standard (F 14 C), and deviations in parts per thousand (‰) from the absolute (decay-corrected) OX1 standard (Δ 14 C).All results have been corrected for mass-dependent isotopic fractionation using 13 C measurements.
4 Calculations and data quality control 4.1 Correction of surface CO 2 emissions for atmospheric contamination We corrected surface-chamber radiocarbon measurements for atmospheric contamination using the method described in Schuur and Trumbore (2006).Briefly, we determined fractional contributions of background atmosphere and ecosystem respiration to total chamber gas using 13 C values in a twopool mixing model: where f Reco and f atm are the fractional contributions of ecosystem respiration and background atmosphere, 13 C S and 13 C atm are the measured 13 C abundances in the sample and background atmosphere in units of at.%, and 13 C Reco is the 13 C abundance in ecosystem respiration, approximated separately for each polygon type as the mean 13 C of chamber CO 2 samples with [CO 2 ] > 4000 ppm.Mean δ 13 C Reco values calculated this way were −24.6, −26.5, and −26.2‰ for LC, FC, and HC polygons.To minimize error due to large proportions of atmospheric CO 2 , we omitted samples with f Reco <0.5.For each sample, we calculated Δ 14 C of ecosystem respiration (Δ 14 C Reco ) according to Eq. ( 3): where Δ 14 C S and Δ 14 C atm are the measured Δ 14 C values of the sample and background atmosphere, and f Reco and f atm were calculated from Eqs. ( 1) and (2).Errors due to analytical precision and variations in source isotopic signatures were propagated through this correction according to the formulation in Phillips and Gregg (2001).
Soil surface CO 2 flux measurements were assessed for quality using two criteria, the standard error of the slope (SE) and the per cent relative standard error (PRSE), defined as 100×SEslope/Estimateslope. Flux measurements with SE > 0.05 and PRSE >5 were omitted from the dataset.This set of dual criteria avoided biasing the dataset toward low fluxes (if SE alone were used) or high fluxes (if PRSE alone or R 2 were used).

Modeled mean age of soil pore-space CO 2
With soil pore-space radiocarbon data, we omitted measurements with CO 2 yields below the expected yield for atmospheric measurements due to possible leakage and atmospheric contamination during sampling, with the exception of one sample from 31 cm depth with highly negative Δ 14 C, indicating a low proportion of atmospheric CO 2 .At this field site, 13 C abundances in pore-space CO 2 vary greatly due to isotopic fractionation from methane production and consumption (Vaughn et al., 2016).For this reason, we could not correct subsurface samples for atmospheric CO 2 .Reported radiocarbon values thus represent the total CO 2 present in the soil pore space, sourced from heterotrophic respiration, root respiration, and downward atmospheric diffusion.
With soil-profile CO 2 samples, we used radiocarbon measurements to model the mean age of carbon in respired CO 2 using the time-dependent steadystate turnover time model described in Torn et al. (2009), modified to account for carbon residence time in plant tissues: where of the given carbon pool, equal to F ′ of the CO 2 sample, F′atm=F′ of CO 2 in the local atmosphere, I= input rate of carbon from the atmosphere to the given carbon pool (g C y −1 ), C= stock of carbon in the given carbon pool (g), τ= turnover time of the given carbon pool, equivalent to the mean age of carbon in its decomposition flux (y) at steady state, λ= radioactive decay rate of 14 C (1∕8267 y), and T R = mean residence time of carbon in plants before entering soil organic matter (y).
At steady state, C t =C t −1=I×τ, so Eq. ( 4) reduces to Following Eq. ( 5), the Δ 14 C value of CO 2 at time t thus depends on the turnover time of carbon in the decomposing carbon pool, the mean residence time of carbon in plant material, and the Δ 14 C of atmospheric CO 2 in the current and previous years, which has changed continuously since the release of radiocarbon into the atmosphere from nuclear weapons testing between 1950 and the mid 1960s (Trumbore, 2000).This model assumes that CO 2 is derived from a homogeneous pool of decomposing carbon, such that the turnover time of this pool is equal to the mean age of its decomposition flux at steady state (Sierra et al., 2017).
The mean residence time of carbon in plants (T R ) reflects a mixture of materials with varying transfer rates.In the Arctic vegetation dominant at our site, some photosynthates enter the soil within 1 day of fixation (Loya et al., 2002); leaf tissues with transit times of <1 year represent the majority of annual production (Shaver and Kummerow, 1991); roots and shoots are estimated to live for ranges of 3-7 and 1-8 years, respectively (Chapin III et al., 1980;Shaver and Billings, 1975), and rhizomes and stem bases, the longest-lived of belowground tissues, have estimated lifetimes of 2.7-15.6 years (Dennis, 1977).Overall, T R represents the flux-weighted average of these various pools.Given the large percentage of carbon flux dedicated to annual leaf production and rapid transport to soils, we estimate that across plant organs, plant species, and seasons, the mean value of T R lies between 0 and 5 years.
Annual atmospheric Δ 14 C values were compiled from our data and three other sources: the IntCal13 dataset (Reimer et al., 2013), measurements from Fruholmen, Norway, between 1962and 1991(Nydal and Lövseth, 1996), measurements from Utqiaġvik between 1999 and 2007 (Graven et al., 2012), and our measurements from the BEO in 2012, 2013, and 2014.From roughly the same latitude as Utqiaġvik, Fruholmen Δ 14 C measurements provide a close approximation for missing Utqiaġvik data (Meijer et al., 2008).Because ecosystem CO 2 uptake occurs primarily during the growing season, we averaged June-August Δ 14 C values to produce an annualized summer dataset.Data gaps from 1992-1998 and 2008-2011 were filled using an exponential interpolation constrained by the available data from the 10 years surrounding each gap.
Using our measured Δ 14 C values and annually resolved atmospheric Δ 14 C data, we iteratively solved for the mean age of each CO 2 sample.We performed this calculation twice, using T R values of 0 and 5 to bracket the likely T R range.Samples containing a large percentage of recently fixed carbon yielded two possible solutions for each T R value (Trumbore, 2000).In such cases, we chose the appropriate solution either by comparing the two values to other (unique) mean age values within the same profile, or by comparing the measured carbon stock with the carbon stocks calculated from the CO 2 production rate and the candidate solutions (Torn et al., 2009).

Correction of soil organic carbon for respiration losses
To calculate Δ 14 C values of soil organic carbon, we corrected bulk soil Δ 14 C measurements for loss of CO 2 respired during the 387-day incubation.We used the post-incubation masses of non-respired soil carbon (C soil , calculated from measured soil carbon contents and sample masses) and carbon respired as CO 2 (CCO2, calculated from CO 2 production rates) to calculate pre-incubation Δ 14 C values of bulk soil organic carbon (Δ 14 C SOC ) according to where Δ 14 C soil and Δ 14 C CO2 are the radiocarbon abundances measured from bulk soil and respired CO 2 following the incubation.

Radiocarbon in ecosystem respiration
Radiocarbon contents of soil surface CO 2 emissions, equivalent to ecosystem respiration (Δ 14 C Reco ), ranged from +60.5 to −160 ‰, with a median value of +23.3 ‰ (Table S1 in the Supplement).The positive Δ 14 C Reco values measured in 28 of the 37 samples indicate high proportions of carbon fixed since 1950, likely sourced from both autotrophic respiration and decomposition of rapidly cycling soil carbon in shallow soil layers.In contrast, the negative values measured in the other nine samples show that the sources of CO 2 at those times and places were dominated by carbon that cycles on centennial to millennial timescales.Δ 14 C Reco followed a left-skewed distribution, with notably negative Δ 14 C Reco values measured from LC1-center in October 2012 (−159.5 ‰) and HC1-center in September 2013 (−115.5 ‰).The variations in these Δ 14 C Reco values reflect differences in the relative CO 2 production and transport rates from both autotrophic and heterotrophic source pools, whose decomposition rates and radiocarbon contents vary spatially and temporally (Nowinski et al., 2010).Together, seasonal variations in soil thaw depths and spatial variations in soil organic carbon 14 C depth profiles are known to affect the radiocarbon signature of surface CO 2 emissions.
From data compiled over the 3 measurement years, the lowest Δ 14 C Reco values were measured in September, when the depth of soil thaw was near its annual maximum (Fig. 2).This observation is consistent with other data from high-latitude sites (Hicks Pries et al., 2013;Trumbore, 2000) and reflects seasonal changes in thaw depth, soil temperature, and vegetation activity.Surface soils at this site begin to thaw in June, typically reaching their maximum temperatures in July (Hinkel et al., 2001;Torn, 2015).As shallow soils warm and plant activity increases in this early summer period, ecosystem respiration includes high proportions of 14 C-enriched CO 2 from autotrophic respiration and heterotrophic decomposition of shallow, rapidcycling soil carbon (Fig. 2).Autotrophic respiration peaks in July or August, and decreases substantially into the fall after plants senesce (Hicks Pries et al., 2013), shifting the balance of respiration toward soil carbon decomposition.As thaw depth increases into October (Table S1) (Zona et al., 2016), decomposition of deep, 14 C-depleted soil carbon can become an increasing fraction of total soil respiration.The effect of these changes is a seasonal shift in respiration from primarily shallow, fast-cycling source carbon pools to deeper, more 14 C-depleted soil organic matter.
Whereas September Δ 14 C Reco values were notably low, this was not the case with October measurements (Fig. 2).Although it is possible that this observation reflects the relative carbon mineralization rates of the different soil depths, the relatively enriched Δ 14 C Reco values measured in October may instead be due to physical limitations to gas transport within the soil profile or bias in the distribution of samples across polygon types and features.Surface soils were frozen at the time of October sample collection, possibly trapping deeper, 14 C-depleted CO 2 below the surface soil layers.This effect has been proposed as a mechanism underlying pulses of CO 2 release during the autumn freeze-up and spring thaw periods (Bubier et al., 2002;Oechel et al., 1997;Raz-Yaseef et al., 2017).In the absence of such a pulse, surface CO 2 fluxes may be derived largely from shallow soil decomposition, yielding relatively enriched Δ 14 C Reco values late in the autumn season.Additionally, because of snow and frozen surface soils, we were unable to obtain CO 2 samples from a number of locations.Notably, measurements from the centers of HC and FC polygons are absent from the October dataset.If the seasonal decrease in Δ 14 C Reco were greatest in polygon centers, the lack of 14 C-depleted October measurements may thus be due to this sampling bias.Time-series measurements from individual profiles or polygon types indicate that Δ 14 C Reco varied not only with time in the sampling season, but also with polygon type (Fig. 2, Table S1).HC polygons displayed particularly strong seasonal trends in Δ 14 C Reco ; in each HC polygon from which multiple measurements were made, Δ 14 C Reco decreased throughout the summer months to a minimum value in September or October.The magnitude of this seasonal decrease varied greatly among individual HC polygons and microtopographic positions within the polygons, reaching a minimum Δ 14 C Reco value of −115 ‰ from the center of polygon HC1, but only +13 ‰ from the trough of HC3 (Table S1).In contrast with HC polygons, Δ 14 C Reco from FC and LC polygons remained closer to atmospheric values throughout the sampling season (Fig. 2).One exception was one highly negative Δ 14 C Reco value measured in October from the center of polygon LC1.In this instance, with October thaw depths close to their annual maximum (Zona et al., 2016), thaw may have penetrated into the transition layer at the top of the permafrost, exposing old, previously frozen carbon to decomposition.Alternatively, this isolated measurement of 14 C-depleted CO 2 may reflect heterogeneity within the active layer due to cryoturbation, the vertical transport of soil material due to freeze-thaw processes (Kaiser et al., 2007;Ping et al., 1998).
In September 2014, we measured Δ 14 C Reco , Δ 14 C SOC , and ecosystem respiration rates from the centers of three polygons of each type (Fig. 3, Table S2).The influence of microtopography on old carbon emissions was particularly apparent in this measurement set, which offers an even distribution of measurements across polygon types.At this time, decomposition of old carbon consistently dominated the respiration flux from HC polygon centers (Δ 14 C Reco =−51.6±32.4‰), whereas fast-cycling pools dominated respiration from FC and LC polygon centers (13.0±10.9 and 5.9±5.6 ‰, respectively; Fig. 3c).Ecosystem respiration rates were comparable across the three polygon types (Fig. 3b), which indicates that absolute rates of old carbon decomposition were greater from HC polygons than from LC or FC polygons.The spatial pattern in old carbon decomposition was not tightly linked to Δ 14 C SOC , which was 14 C-depleted relative to the respiration flux in both shallow and deep samples (Fig. 3c, d, Tables S1, S2).Differences in 14 C between bulk soil carbon and the respired fraction are commonly observed across ecosystems and soils, as soil carbon is comprised of components with a range of cycling rates (Dioumaeva et al., 2002;Hicks Pries et al., 2013;Schuur and Trumbore, 2006;Trumbore, 2000).Additionally, spatial and temporal variations in autotrophic respiration obscure relationships between soil-respired CO 2 and ecosystem respiration.
We saw no correlation in September 2014 between Δ 14 C Reco and thaw depth, which was greatest at this time in FC polygons (Fig. 3a).In contrast, repeated measurements from individual soil locations showed a general decrease in Δ 14 C Reco as thaw depth increased and exposed deeper soil layers to unfrozen conditions (Table S1).These findings suggest that the relationship between the depth of thaw and old carbon mineralization depends on the spatial and temporal scales of observation.At the scale of an individual soil profile, seasonal variations in Δ 14 C Reco correspond with changes in thaw depth.Δ 14 C of soil organic matter tends to increase with depth in the soil, particularly in permafrost environments where the brevity of summer thaw strongly limits carbon turnover at depth (Hicks Pries et al., 2013;McFarlane et al., 2013;Torn et al., 2002;Trumbore, 2000;Fig. 3d).Accordingly, increasingly 14 C-depleted soil organic matter becomes available to decomposers as soil thaw deepens over the summer season.At the site scale, however, thaw depth may not be a useful predictor of spatial variations in Δ 14 C Reco .Other differences -including variations in soil temperature profiles, organic layer thickness, vegetation productivity and rooting depth, the presence of cryoturbation, and oxygen availability to decomposers (Newman et al., 2015;Ping et al., 2015;Sloan et al., 2014;Vaughn et al., 2016) -may lead to microtopographic differences in the shape of Δ 14 C SOC depth profiles and the relative contributions of different carbon source pools to the decomposition flux.
In July and September 2013 and September 2014, ecosystem respiration ranged between 0.32 and 2.5 µmol CO 2 m −2 s −1 (Fig. 4).The relationship between the rate and radiocarbon content of ecosystem respiration was highly variable, but with two notable patterns across all data.First, as ecosystem respiration increased, the variance of Δ 14 C Reco tended to decrease.Second, there were no negative Δ 14 C Reco measurements associated with fluxes above 1 µmol CO 2 m −2 D −1 (Fig. 4).These patterns reflect the sensitivity of ecosystem respiration to the most dynamic carbon sources, autotrophic respiration and shallow soil carbon mineralization.Autotrophic respiration can contribute as much as 70 % of ecosystem respiration at the peak of the growing season (Hicks Pries et al., 2013), declining markedly as plants senesce.Similarly, heterotrophic decomposition rates are both highest and most temporally variable in shallow soil layers (Hicks Pries et al., 2013), where soil temperatures exhibit a large seasonal range.Supporting our interpretation, previous studies documenting seasonal variations in CO 2 exchange at this site (Oechel et al., 1995;Vaughn et al., 2016;Zona et al., 2014) indicate that ecosystem respiration, gross primary production, and soil surface temperatures peak in July or August, then decline into the late summer and fall.During the early summer period when respiration fluxes are high, high rates of autotrophic respiration and shallow soil respiration dominate the overall respiration flux.As a result, old, slow-cycling carbon respired from deep soils comprises a large percentage of ecosystem respiration only when respiration rates are low from autotrophic and shallow (fast-cycling) soil sources.

Radiocarbon in soil pore-space CO 2
At a subset of locations and sampling dates, we measured the radiocarbon content of CO 2 in soil pore gas (Δ 14 C CO2p ).Δ 14 CCO2p varied widely among profiles and sampling dates from −7.1 to −280 ‰ (Table 1, Fig. 5).These consistently negative values indicate that pore-space CO 2 was derived primarily from older soil organic matter, with minimal contributions from plant-respired carbon or fast-cycling soil organic carbon.In contrast, ecosystem respiration was generally enriched in radiocarbon relative to soil pore-space CO 2 , even at only 10 cm depth (Fig. 5).This observation suggests that although old, slow-cycling carbon dominates soil pore-space CO 2 , the shallowest CO 2 sources -autotrophic respiration and/or heterotrophic decomposition of fast-cycling (annual-decadal) organic carbon -contribute large proportions of the total respiration flux, even late in the season when plants have largely senesced.Detecting and characterizing the decomposition of older, deeper soil organic carbon requires direct measurements of soil pore-space CO 2 .Such measurements provide a qualitative indicator of old carbon decomposition; as with surface CO 2 effluxes, proportional or absolute contributions from distinct carbon pools cannot be calculated without well-resolved vertical distributions of 14 C and 13 C source pools, as well as gas transport rates within the profile.
In most soil profiles, Δ 14 C CO2p became increasingly negative with depth in the soil.This depth trend is similar to that seen in Δ 14 C SOC profiles measured from this site (Fig. 3d, Table S2), indicating that the cycling rates of both bulk soil organic matter and the decomposing carbon fraction tend to decrease with depth in the soil.Because frozen or near-freezing temperatures slow decomposition from the deep active layer throughout the majority of summer, this pattern is what would be expected in the absence of vertical mixing or rapid contributions of fast-cycling carbon at depth.In contrast with this general trend, soil pore-space CO 2 from three soil profiles (HC3-trough, HC3-center, and FC2-center) became enriched in radiocarbon near the permafrost table (Fig. 5, Table 1), a pattern that has been previously observed in the Arctic (Lupascu et al., 2014a).In these three profiles, concentrations of CO 2 in deep pore-space samples were 6 to 25 times higher than background atmosphere, so we infer that the higher Δ 14 C CO2p values at depth were not caused by downward transport of atmospheric CO 2 .Instead, this observation indicates significant CO 2 production from fast-cycling carbon near the permafrost table.As in other Arctic sites, this depth trend is likely due either to cryoturbation (Bockheim and Tarnocai, 1998) or DOC leaching (Lupascu et al., 2014a), both of which can transport recently fixed, relatively decomposable carbon to the deep section of the active layer.
With Δ 14 C CO2p values, we used a single pool turnover time model to calculate the mean age of pore-space CO 2 (Table 1).This model assumes that CO 2 is respired from a homogeneous pool of decomposing carbon (Trumbore, 2000) such that the turnover time equals the mean age of the respired carbon (Sierra et al., 2017).Across all measurement dates and depths, we found the mean age of respired CO 2 ranged from 410 years in one 10 cm sample to 3350 years at only 20 cm depth, indicating decomposition of old and/or slowcycling carbon, even in relatively shallow soils.HC polygon centers produced particularly old CO 2 , with ages greater than 3000 years from HC3-center in July 2013 and HC1-center in August 2012 and September 2013.This observation aligns with the spatial patterns we observed in surface CO 2 emissions, which were particularly 14 C-depleted in HC polygons (Fig. 2, Table S1).Together, these findings suggest that low Δ 14 C Reco from HC polygon centers was due to mineralization of very old carbon at depth.Because of unknown diffusion rates within the soil profile, it is challenging to quantitatively use pore-space CO 2 measurements to link soil carbon cycling rates to soil surface Δ 14 C Reco .Low diffusion rates can lead to the accumulation of high concentrations of 14 C-depleted, slow-cycling CO 2 deep in the soil profile (Lee et al., 2010), and vertical mixing can transport CO 2 away from the site of production prior to its collection as soil pore gas.
Our findings of consistently negative Δ 14 C values from soil pore-space CO 2 are similar to those from the other Arctic sites where it has been measured (Czimczik and Welker, 2010;Lupascu et al., 2014b, a).In contrast, studies in temperate or tropical sites find Δ 14 C CO2p values that are more similar to atmospheric values (Borken et al., 2006;Gaudinski et al., 2000;Trumbore, 2000), where soil respiration is dominated by root respiration and decomposition of rapidly cycling soil organic matter.This difference may be due in part to differences among biomes in rooting distributions.Because of strong vertical gradients in soil temperature and nutrient and water availability, rooting distributions in Arctic tundra are extremely shallow (Iversen et al., 2015), limiting most root respiration to the shallow organic layer.Additionally, such negative Δ 14 C CO2p values suggest that old, slowcycling organic matter in Arctic soils is readily decomposable under thawed conditions.Cold soil temperatures and frozen conditions limit microbial activity throughout much of the year, allowing otherwise-decomposable soil organic matter to accumulate.During the short season in which Arctic soils are thawed, these large pools of 14 C-depleted soil carbon may be largely unprotected against microbial mineralization.Compared with tropical or temperate soils, existing Arctic soil carbon stores may thus be particularly vulnerable to warming.
Due to sampling limitations, 14 C data from soil pore-space CO 2 are available from only a subset of polygon types, positions, sampling dates, and depths.No samples were collected from polygon rims, and only few samples were obtained where soils were saturated (Table 1).For this reason, our dataset does not capture the full range of microtopographic variations in deep soil decomposition rates and controls.Soil temperature profiles, soil pore-space oxygen availability, soil organic matter concentration and chemistry, and vegetation composition, rooting distribution, and productivity influence microbial activity and vary among profiles and with time in the thawed season (Lipson et al., 2012;Olivas et al., 2011).In particular, these controls on decomposition dynamics likely differ between saturated and unsaturated soils, due both to the direct effects of saturation and to covarying effects of microtopography.For this reason, radiocarbon values presented in this manuscript cannot be scaled across the full range of environmental variation present at the site.To better characterize the relationships between these variables and soil carbon decomposition rates, further measurements are needed of 14 C in soil pore-space DIC and CO 2 , across spatial, seasonal, and hydrological gradients.
6 Data availability Data are stored in the Next-Generation Ecosystem Experiments (NGEE-Arctic) data repository (Vaughn et al., 2018;Vaughn and Torn, 2018) and may be accessed at http://dx.doi.org/10.5440/1364062and https://doi.org/10.5440/1418853.The code used to generate plots and correct chamber samples for atmospheric contamination can be found at https://github.com/lydiajsvaughn/Radiocarbon_field_Barrow.The code used for turnover time modeling can be found at https://github.com/lydiajsvaughn/Radiocarbon_inc_2012.

Conclusions
We measured the radiocarbon contents of ecosystem respiration and soil pore-space CO 2 between 2012 and 2014 in Utqiaġvik, Alaska.As a naturally occurring tracer of soil carbon dynamics, radiocarbon in CO 2 and soil carbon reflects the rates at which carbon cycles through plant and soil pools (Trumbore, 2000).Radiocarbon provides a powerful tool to test and parameterize bottom-up models of the carbon cycle (He et al., 2016), constrain source signatures in inverse modeling studies (Randerson et al., 2002), and evaluate how environmental variables, seasonal changes, and disturbance influence decomposition rates from fast-and slow-cycling carbon pools.In Arctic sites, where cycling rates can be very slow, radiocarbon can be used to detect changes in decomposition from millennial-cycling pools (Hicks Pries et al., 2013;Schuur et al., 2009), which strongly influence longterm carbon dynamics but may be impossible to detect with other methods.Only limited data are available of radiocarbon in ecosystem respiration or soil pore-space CO 2 , particularly from Arctic sites where samples can be difficult to obtain due to low rates of CO 2 production.To characterize seasonal and spatial variations in the age of carbon pools contributing to total CO 2 efflux, we measured radiocarbon in ecosystem respiration and soil pore-space CO 2 across three summer seasons in Arctic polygon tundra.Δ 14 C in ecosystem respiration varied between +60.5 and −160 ‰ across the 3 years, with a strong seasonal trend (Fig. 2).July and August Δ 14 C Reco measurements were generally close to the Δ 14 CO 2 values of the local atmosphere, which declined from 22.2 ‰ in 2012 to 17.7 ‰ in 2014.In September, in contrast, over half the Δ 14 CO 2 measurements differed greatly from the local atmosphere (Fig. 2), contributing 14 C-depleted CO 2 to the atmospheric pool.These seasonal variations in Δ 14 C Reco and respiration rates contribute to the strong seasonal cycles of atmospheric Δ 14 CO 2 observed in the high latitudes (Graven et al., 2012).By quantifying these variationsrespiration flux rates and radiocarbon values over space and time -this dataset offers useful information for atmospheric budgeting and inversions.
Measurements of radiocarbon in late-season ecosystem respiration indicate that carbon that cycles on millennial timescales contributes substantially to soil respiration.When thaw depth approached its maximum in September and October, highly depleted 14 C in respiration indicated that carbon older than 1000 years was a major source of heterotrophic respiration.Decomposition of old, slow-cycling soil carbon, however, may be missed when measuring surface Δ 14 C Reco alone.In the soil pore space, Δ 14 C CO2p declined steeply with depth, particularly in high-centered polygons (Fig. 4); at only 20 cm below the soil surface, the mean age of carbon in the decomposition flux was as old as 3000 years (Table 1).Together, these observations suggest that during the short summer thaw season, ancient carbon stores in the deep active layer become available to decomposition, stabilized otherwise by cold, often frozen, and often anaerobic conditions (Mikan et al., 2002;Schmidt et al., 2011;Schuur et al., 2015).
As climate change alters these environmental controls and soils warm and thaw, a key question is how decomposition rates will change.Where permafrost degradation occurs, either through gradual deepening of the active layer or rapid thaw events, old soil carbon that has remained frozen for years to millennia is exposed to thawed conditions.A particularly important -but unknown -factor is the decomposition rate of this carbon released from thawing permafrost (Hicks Pries et al., 2013;Koven et al., 2015;Kuhry et al., 2013).Our measurements cannot differentiate between such newly thawed soil organic matter and carbon that has historically experienced an annual thaw.Our data do show, however, that when it is unfrozen, even very old Arctic soil organic carbon will readily decompose.As thaw depth progressed to the permafrost boundary, we consistently documented CO 2 production via decomposition of centennial carbon pools.As permafrost thaw progresses and Arctic soils warm, further measurements of soil-respired 14 CO 2 across spatial and temporal gradients will provide critical information on soil carbon vulnerability.