Figures

Results

Respiration rates

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 second enclosure 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 CO2 targets are excluded for display reasons; respiration rates for those samples remained flat. Note that headspace CO2 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.

Radiocarbon data

Pre-incubation versus equilibrium respiration 14C-CO2

Fig. 2. \(\Delta\)14C-CO2 of the rewetting pulse compared to that of the second enclosure period

Caption: Points are means of laboratory duplicates and error bars are the min and max. Note that rewetting pulse \(\Delta\)14C 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 CO2 during the pre-incubation period to measure \(\Delta\)14C. The outlier point with the substantially depleted pre-incubation \(\Delta\)14C is from Experiment 2 (control).

Treatment effect on 14C-CO2 for all samples (Experiments 1, 2, and 3)

Fig. 3. Overall effect of air-dry/rewet treatment with and without storage on 14C-CO2

Caption: Points are from all three experiments: the treatment for Experiments 1 and 3 was air-dry/rewet + storage, while Experiment 2 samples underwent air-drying and rewetting without storage. Points are the mean of laboratory replicates (for replicated samples); error bars on points 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/\)14C in atmospheric CO2 over 4 and 8 y respectively, during the period of 2000 to 2020 (Graven et al. 2017).

Fig 4. Time series of control and treatment \(\Delta\)14C-CO2 (Experiments 1 and 2)

Caption: Filled circles show \(\Delta\)14C-CO2 observed for both control-1 and control-2 samples (2011 and 2019 points, respectively). Open symbols show \(\Delta\)14C-CO2 observed for treament samples: open squares = air-dry + storage treatment, Experiment 1; open circles = air-dry only treatment, Experiment 2. Arrows show the direction of change in \(\Delta\)14C-CO2 relative to the controls. Points are means and error bars show 2x standard error. The gray line shows \(\Delta\)14C of the atmosphere.

Storage duration effect on 14C-CO2 (Experiment 3)

Fig. 5. Change in 14C-CO2 in relation to storage duration

Caption: Points show the mean differences of paired treatment and control samples, averaged by site and ecosystem type. All data are from Experiment 3, except for the Central Germany points, where data from both Experiment 1 and Experiment 3 were included. 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/\)14C in atmospheric CO2 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.

Discussion

Fig 6a. Modeled trajectories of \(\Delta\)14C over time for a hypothetical 2-pool soil carbon model in relationship to observed atmospheric \(\Delta\)14C

Caption: Modeled curves derived from a hypothetical two-pool parallel model system in which inputs are partitioned between a fast cycling soil carbon pool (k = 1/6) and a more slowly cycling pool (k = 1/100) without any tranfers between the pools. Carbon stocks and pool sizes are based on density fraction data for the Hainich Forest (Schrumpf et al. 2013); decomposition rates are chosen to be realistic but are arbitrary. Atmospheric \(\Delta\)14C 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\)14C of respired CO2 in response to treatment

Caption: Blue and magenta lines show modeled trajectories of \(\Delta\)14C in fast and slow cycling soil carbon pools (same model as Fig 1a), while \(\Delta\)14C of respired CO2 is shown in yellow. Filled circles show hypothetical observations of \(\Delta\)14C-CO2 from control incubations, while the open symbols represent two possible scenarios in which an air-drying and rewetting treatment leads to shifts in \(\Delta\)14C-CO2. In the first scenario (open squares), an increased contribution from the slow pool shifts \(\Delta\)14C-CO2 toward the slow pool curve (blue line), while in the second scenario (open circles), an increased contribution from the fast pool shifts \(\Delta\)14C-CO2 toward the fast pool curve (magenta line). Due to the crossing of the slow and fast pool curves in 2009, an increased contribution of the slow pool to respiration following treatment leads to relative depletion of \(\Delta\)14C-CO2 in 1991, but relative enrichment of \(\Delta\)14C-CO2 in 2019, while the opposite is observed for the increased fast pool contribution scenario.