Archived Soil Incubations Manuscript Figures

Notes

• This workbook is intended to load and prepare the key data for analysis for the archive incubation project.
• This is an updated version of RMD file “./src/arc_inc_data-wrangling_2020-04-20.Rmd” with the explanatory text removed
• all code chunk options are set to “echo = FALSE”; see raw .Rmd file for data wrangling code.

Methods and materials

Radiocarbon corrections

Although differences were small, radiocarbon data for archived samples were first corrected for decay since the year of collection before assessing treatment effects.

Decay correction formula: \[1000 \cdot \left( (FM \cdot e^{\frac{-year_{sampled} + 1950}{8267}}) - 1 \right)\]

where FM is the fraction modern, and 8267 is the inverse of the product of the natural log of two and the true half life of ¹⁴C (5730 y).

Results

Respiration rates

note that the text below was left because statistics were calculated on the fly

Experiment 1 (air-dry + storage treatment)

Among the air-dry + storage samples, respiration rates were more than twice as high in grassland soils than in forest soils, reaching a maximum of 7.1 mg CO₂ g soil C^-1 d^-1 after 0.7 d, followed by a sharp decline. Mean respiration rates in forest sites peaked at 1.5 mg CO₂ g soil C^-1 d^-1 after 6.9 d, followed by a much more gradual decline than in grassland sites. Control samples responded more weakly and more gradually to rewetting, although as in the treatment samples respiration was greater in grassland soils than in forest soils. Peak respiration rates for control incubations were 1.9 and 0.6 mg CO₂ g soil C^-1 d^-1 after 4.8 d for grassland and forest soils, respectively.

Experiment 2 (air-dry treatment, 2019 samples)

Peak respiration rates were not significantly different (p > 0.05) between forest and grassland soils in Experiment 2, peaking at 3.0 and 3.3 mg CO₂ g soil C^-1 d^-1 after 95 h for grassland and forest soils, respectively (Fig. 1).

Fig. 1. Respiration rates for Experiment 1 (air-dry + storage treatment, 2011) and Experiment 2 (air-dry only treatment, 2019)

Caption: Top panel shows data from samples collected in 2011 for Experiment 1 (air-dry + storage treatment), bottom panel shows data from samples collected in 2019 for Experiment 2 (air-dry only treatment). Vertical gray line at day 4 demarcates the end of the pre-incubation period and the start of the equilibrium respiration period. Points show measurements and lines show trends in mean respiration rate. Shaded ribbons represent one standard error. The final measurement points for a few samples which took >18 days to reach CO₂ targets are excluded for display reasons; respiration rates for those samples remained flat. Note that headspace CO₂ concentrations for Experiment 1 control samples were only measured once during the pre-incubation period (day 4) in contrast to daily measurements for all other samples. Consequently the respiration rate for those samples is the cumulative average rate over the first 4 d.

Experiment 3 (storage duration)

Grassland soils from the storage duration treatment group responded rapidly to rewetting and reached target CO₂ levels after just 72 h of incubation. Only a single observation was made for the grassland treament samples due to the rapid respiration rates, which peaked at 5.6 mg CO₂ g soil C^-1 d^-1. The mean peak respiration rate for forest treatment samples was lower and lagged in comparison to the grassland soils, reaching 2.2 mg CO₂ g soil C^-1 d^-1 after 97 h.

Control-3 respiration rates peaked after the pre-incubation period (115 h), at 2.3 mg CO₂ g soil C^-1 d^-1. Forest control samples were pre-incubated under various conditions, but respiration rates were only measured during the pre-incubation period for two samples from the Sierra Nevada mountains (USA). In general, respiration rates for forest control samples in Experiment 3 were much lower than in the treatment incubations, peaking at 0.6 mg CO₂ g soil C^-1 d^-1 after 120 h.

Radiocarbon data

no non-missing arguments to max; returning -Infno non-missing arguments to max; returning -Inf

Pre-incubation versus equilibrium respiration ¹⁴C-CO₂

Fig. 2. \(\Delta\)¹⁴C-CO₂ of the rewetting pulse and the equilibrium respiration period

Caption: Points are means of laboratory duplicates and error bars are the min and max (except for Experiment 1 control samples, which were not replicated). Note that rewetting pulse \(\Delta\)¹⁴C was not measured for control-1 samples; additionally samples from three of the forest plots of the air-dry + storage samples from Experiment 1 failed to accumulate enough CO₂ during the pre-incubation period to measure \(\Delta\)¹⁴C. The outlier point with the substantially depleted pre-incubation \(\Delta\)¹⁴C is from Experiment 2 (control).

Treatment effect on \(\Delta\)¹⁴C-CO₂ for all samples (Experiments 1, 2, and 3)

no non-missing arguments to min; returning Infno non-missing arguments to min; returning Infno non-missing arguments to min; returning Infno non-missing arguments to min; returning Inf

Fig. 3. Overall treatment effect on \(\Delta\)¹⁴C-CO₂

Caption: Points show data from Experiments 1 and 3 (air-dry + storage treatment) and Experiment 2 (air-dry only treatment). Points are the mean of laboratory replicates (for replicated samples); error bars are 2x standard error. Solid line is 1:1. For context, the dashed and dotted lines show differences of ±20‰ and ±$40‰, equivalent to the decline in \(\Delta\)¹⁴C in atmospheric CO₂ over 4 and 8 y respectively, during the period of 2000 to 2020 (Graven et al. 2017).

Storage duration effect on ¹⁴C-CO₂ (Experiment 3)

Fig. 4. Change in ¹⁴C-CO₂ in relation to storage duration

Caption: Data are from both Experiment 1 (in black) and Experiment 3 (all other points), averaged by site and ecosystem type. Points are the mean, error bars are 2x standard error. For context, the dashed and dotted lines are the same is in Fig. 4 and show a difference of 20‰ and 40‰, equivalent to the decline in \(\Delta\)¹⁴C in atmospheric CO₂ over 4 and 8 y respectively, during the period of 2000 to 2020 (Graven et al. 2017). Position of points jittered to avoid overplotting; storage duration has been rounded down to the nearest whole year.

Fig 5. Time series of control and treatment \(\Delta\)¹⁴C-CO₂ (Experiments 1 and 2) Caption: Filled circles show \(\Delta\)¹⁴C-CO₂ observed for both control-1 and control-2 samples (2011 and 2019 points, respectively). Open symbols show \(\Delta\)¹⁴C-CO₂ observed for treament samples: open squares = air-dry + storage treatment, Experiment 1; open circles = air-dry only treatment, Experiment 2. Points are means and error bars show the pooled standard deviation. The black line shows \(\Delta\)¹⁴C of the atmosphere.

Supplemental Fig 3. Time series of control and treatment \(\delta\)¹³C-CO₂ (Experiments 1 and 2)

Caption: Filled circles show \(\delta\)¹³C-CO₂ observed for control samples, while open symbols show \(\delta\)¹³C-CO₂ observed for treament samples (open squares = air-dry + storage treatment, Experiment 1, 2011; open circles = air-dry only treatment, Experiment 2, 2019). Points are means and error bars show standard deviations.

Effect of cumulative respired carbon on ¹⁴C-CO₂

Supplementary Fig. 4. Change in ¹⁴C-CO₂ in relation to cumulative soil carbon respired

Caption: Note that pre-incubation \(\Delta\)¹⁴C was not measured for the control-1 samples in 2011. Limits exclude outlier point (HEW22 control-2, pre-incubation) for improved legibility. Points are means, error bars show min and max of duplicate samples.

Discussion

Fig 6a. Trajectories of \(\Delta\)¹⁴C over time in soil C and respired CO₂ for a 2-pool model fit to the Hainich-Dün forest sites in relation to atmospheric \(\Delta\)¹⁴C

Caption: Modeled curves derived from a two-pool parallel model parameterized with data from Schrumpf et al. (2015): (k_fast = 1/4) and a more slowly cycling pool (k_slow = 1/115). Gold points show observed \(\Delta\)¹⁴C-CO₂ from this study. Solid symbols show control samples, open symbols show treatment samples. Atmospheric \(\Delta\)¹⁴C data up to the year 2015 are from Graven et al. (2017), while data points beyond 2015 use the extrapolation method from Sierra (2018). All atmospheric radiocarbon data is for the northern hemisphere (zone 2).

Fig 6b. Potential shifts in \(\Delta\)¹⁴C of respired CO₂ in response to treatment

Caption: Zooming in on the study period, we can see that the treatment sample \(\Delta\)¹⁴C-CO₂ shifts in the direction of the slowly cycling soil C pool (blue line), indicating increased contribution from this pool to respiration following air-drying and rewetting. Note that the increased contribution from the slow cycling pool leads to depletion in \(\Delta\)¹⁴C-CO₂ relative to the controls in 2011, but enrichment in 2019 due to the crossing of the slow and fast cycling pool curves in 2016. Shaded ribbon around the respiration curve shows the model parameterized with the 75th quartile estimates for k_fast (1/7) and k_slow (1/200).

Change in soil moisture content with moisture adjustment

Supplemental Fig 5. Change in \(\Delta\)¹⁴C-CO₂ (control - treatment) relative to field moisture

Caption: Data are from Experiment 1 (“arc”) and Experiment 2 (“rewet”). All samples were moisture-adjusted prior to incubation, but control samples were adjusted from field moisture, “whcFresh” (percent of WHC), whereas treatment samples were moisture adjusted after air-drying, i.e. at approximately 0% of WHC.