# Change in Soil and Forest Floor Carbon after Shelterwood Harvests in a New England Oak-Hardwood Forest, USA

Authors: Kayanna L. Warren; Mark S. Ashton
Journal: International Journal of Forestry Research
 (2014)

---

## Abstract

There has been effort worldwide to quantify how much carbon forests contain in order to designate appropriate offset credits to forest carbon climate mitigation. Carbon pools on or immediately below the soil surface are understood to be very active in response to environmental change but are not well understood. Our study focused on the effects of shelterwood regeneration harvests in New England on the carbon stored in litter, woody debris, and surface soil carbon. Results demonstrate significant difference in surface (0–10 cm) soil carbon between control (nonharvested) and harvested sites, with higher carbon percentage on control sites. Results showed a significant difference in coarse woody debris with higher amounts of carbon per area on harvested sites. No significant difference in litter mass was recorded between harvested and control sites. When coarse woody debris and litter are included with soil carbon, total carbon did not have a significant decline over 20 years following shelterwood treatment to the forest to secure regeneration, but there was considerable variability among sites. When taking all surface soil carbon measurements together, our results suggest that for accounting purposes the measurement of below-ground carbon after shelterwood harvests is not necessary for the southern New England region.


---

## Body

## 1. Introduction

The terrestrial carbon cycle takes up a considerable amount of atmospheric carbon, especially CO2, [1, 2], the increase of which has been cited as the main driver of global climate change [3, 4]. Maintenance of forest cover and reduction of deforestation rates have been proposed by numerous studies and policy processes as a means to reduce the magnitude of climate change [5–7]. It is therefore important that below-ground carbon pools are well understood.Persistent lack of clarity surrounds much of the efforts to tabulate, quantify, and account for the amount of carbon that is stored by forest ecosystems and thus makes its inclusion in policy protocols and its accounting expensive and challenging [8–11]. There is ongoing debate about how management affects overall carbon, with some studies suggesting that old growth forests store more carbon [12–16] while other studies find that continuous timber harvests reduce atmospheric carbon levels more [7, 17]. However, the effects of anthropogenic disturbance, such as harvest, on belowground carbon quantity—belowground pools—is still largely unknown and thinly researched [18–22]. Only a handful of studies on soil carbon pools represent all tropical forests combined, with most studies confined to temperate regions, and most of these are from tree plantations that examine the surface soils only [22]. For temperate forests primary obstacles to quantifying soil carbon include the slow rate of carbon accumulation [23] and the spatial variability of soil [7, 24–26].However, even though little research has focused on below-ground carbon, studies suggest forest soils, litter, and coarse woody debris may comprise half of terrestrial carbon storage and over two-thirds of forest carbon pools [1, 23, 27, 28]. Temperate forests in the eastern United States may be a significant carbon sink [29].Studies on the effects of timber harvest on soil carbon have yielded mixed results thus far [30–32]. For example, soil carbon that is initially lost may take a few years to hundreds of years to recover depending upon forest type, climate, and soil condition [30, 33–35]. Other studies report that soil carbon is not necessarily lost at all [27, 36, 37]. In a recent study Johnson et al. [38] found a consistent drop in soil carbon after a timber harvest for a northern hardwood forest in New England. Carbon accretion in soils of regenerating forests may be further limited by soil texture and other environmental conditions [39]. Information on second-growth hardwood forests is particularly important because conclusions from studies on carbon sink quantification conflict [22].Higher rates of carbon loss occur in the forest floor more than in the soil for hardwood forests [40, 41]. The reasons for changes in soil carbon relate strongly to litter inputs, decomposition, and respiration [42, 43]. Studies suggest that net soil carbon storage may be influenced by small changes in climate and soil condition [41, 44]. Over time and under the colder, wetter climate conditions the shift of live biomass into woody debris, litter, and soil humus may build up stored carbon in soils [7]. Studies conflict over whether forest floor and soil carbon increase or decrease after harvest, and it is unclear whether observed changes are more strongly correlated with the amount of debris left on-site or the change in efflux and respiration with disturbance of the soil [36, 45]. Other studies have shown weak ability to correlate litterfall and belowground carbon allocation in mature temperate hardwood forests [44] and weak conversion of litter carbon into long-term forest soil carbon storage [46]. Not only are the dynamics of decomposition from coarse woody debris and litter to soil carbon quantity poorly understood, but also much remains to be researched about the long-term dynamics of soil carbon after forest harvest.There are few studies which have examined the effects of silvicultural treatments (regeneration or thinning) on soil carbon stocks in the northeastern region of the United States [30, 47]. Some of the limitations of these studies are that they looked only at differences between thinnings [47], were modeled [38], or were focused on coniferous forests rather than the predominant hardwood type [48]. Oak-hardwood forests of the northeastern US have shown reduced soil carbon stocks, reduced carbon in litter, and increased woody debris with harvest [49], but this study was after clearcut harvests, rather than shelterwood. In another study coarse woody debris was reduced by 17% in the long term by a selective timber harvest [50]; but such harvests have shown dramatic increases in coarse woody debris in the short term [51]. Prescriptions that purposefully leave dead wood as snags and slash have been shown to negate any loss of coarse woody debris [52].Our study had two objectives. We wanted to first quantify the change in carbon for the uppermost layers of soil and surface horizons (coarse woody debris and litter) after shelterwood regeneration harvests in a second-growth oak-hardwood forest. Second, we wanted to document the change in soil carbon, coarse woody debris, and litter over time after shelterwood harvests. We sought to understand how soil carbon is affected by forest harvest and how surface carbon is partitioned between coarse woody debris and litter decomposition over a 20-year time horizon after harvest. Our study focuses on the most common type of forest regeneration harvest for oak forests in the northeastern United States: shelterwoods, a method of harvest that retains a limited canopy of masting nut trees (oak, hickory, and beech) to provide seed for dispersal and germination and shelter for newly established seedlings [53, 54]. We expected that surface carbon dynamics will show significant declines similar to findings by Covington [30] and to US Forest Service data from clearcuts. We also expected rate and amount of decline to be lower than in a clearcut because of the moderating shade from tree retention in shelterwood systems. Lastly we predicted that the surface carbon pools will decline with the most significant change in the surface litter.


## 2. Methods

### 2.1. Site Description of the Study Area

This study was conducted at the 3173-ha Yale-Myers Forest near Eastford, CT, in the northeastern United States (41°57′N, 72°07′W). The Yale-Myers Forest is a research and demonstration forest managed by the School of Forestry and Environmental Studies and owned by Yale University. The forest overstory is predominantly oak-maple-pine. The climate is characterized as cool temperate and humid. Mean annual summer temperature is 20°C, and mean annual winter temperature is –4°C. Precipitation is evenly distributed throughout the year with a mean of 110 cm [55, 56]. This is a moist temperate forest underlain by a metamorphic bedrock of schist/gneiss and overlain with soils that originated from till and fluvial sediments of the last Wisconsin glacial period 20,000 b.p. [55]. The soils are mostly coarse-loamy, mesic Typic Dystrudepts [57]. The topography is characterized by undulating hills with broad ridges and narrow valleys. The elevation ranges from a low point of 170 m above sea level to a high point of 300 m above sea level [56].Starting in the early 18th century, the forest, as with many areas of New England, had a long period of agricultural use with land used for various kinds of improved and unimproved pasture, tilled for crops, or relict forest patches that were used as woodlots. Since the 1820s this land has largely been reforested naturally [55].


### 2.2. Treatment Description

All stands for our study were on glacial till soils of ablation origin and comprising the Charlton-Chatfield series [57]. This soil series is a dominant one for the oak-hardwood of southern New England making our study results generalizable to a much larger region than the forest. All study stands were originally cleared for marginal pasture but were never plowed. After abandonment in the 1850s these sites transformed to eastern white pine (Pinus strobus L.) and were then cutover in the early 1900s for the pine boxwood industry. After cutting the pine, understory hardwood regeneration comprising oaks, hickories, birches, and maples was released. Today this second-growth hardwood forest varies in age between 90 and 110 years of age. This age class and forest history are very representative of the southern New England area.Using available GIS maps of soil type and harvests within Yale-Myers, a map of shelterwood harvests and reserves of similar land use history (as described above) overlain with the Charlton-Chatfield soil type was created to identify sites that had suitable treated and control plots adjacent to each other (Figure1). Controls were therefore selected as an uncut reference of similar cover, soil, land use history, and topographic position adjacent to each shelterwood harvest treatment. We did paired comparisons of treated versus control plots to determine if there is a statistical difference (at the 0.05 level or greater) between shelterwood harvest and no harvest. Harvested sites were selected to represent a range of ages from 0 to 20 years since harvest and so can be regarded as surrogates to investigate effects of age.Figure 1
Locations of each of the 15 shelterwood treatment and adjacent control sites in the Yale-Myers Forest.At the Yale-Myers Forest, shelterwood harvests reduce the canopy in a single cutting from a fully stocked basal area of about 37 m2/ha to about 11 m2/ha, largely comprising about forty 50–60 cm dbh well-spaced trees per ha. The aim is to maintain species diversity among regenerating trees and seed sources by insuring the heavy seeded masting trees (oaks, hickories, and beech) adequately disperse seeds and that the partial shade of the canopy will moderate the harsh open conditions during germination and seedling establishment. Seven to eight years later half the standing trees (basal area = 5 m2/ha) are removed in a second and final cut leaving the rest as irregular structures to grow within the new regenerating stand. Harvest operations are generally undertaken during winter months, when the ground is frozen and covered with snow. Skid trails and landings are designated ahead of operations, often utilizing previous trails and historic roads.


### 2.3. Sampling Design

We selected fifteen study sites, comprising pairs of treated and control plots. Based on a prior study on soils within this forest, the range was expected to lie between 1.0 and 7.5 kg m2 of carbon in the top 30 cm [58], with a world average around 11.7 kg m2 of carbon for all soils to 100 cm depth [59]. This study estimated that northeastern temperate soils, less carbon-rich than boreal soils in Canada, would have a spread toward the lower end of that range. The sample size for each treatment block was selected with the University of Iowa’s online sample size calculator [60] with settings for two-way ANOVA. Using this estimation, a sample of five and three plots was determined for the treated and control sites, respectively. Individual plots on each site were selected using a random number generator and a GPS coordinate grid.Plots were centered at the GPS coordinates. Plot corners were defined as three meters north, east, south, and west of the plot center. Two 2.54 cm wide soil cores were taken at each point (center, north, south, east, and west) and pooled into one plot sample (Figure2). Because the upper soil horizon has more carbon than other increments of the mineral soil at deeper levels [58] and because most of the observable change 20 years after harvest was expected to be in the upper 10 cm [58, 61], we cored to a depth of 10 cm at each core point from the mineral soil surface and below the litter and organic soil horizon. This meant that we carefully scraped of the organic horizon to carefully identify the mineral soil surface layer. The forest floor organic horizon was sampled using a 15 cm × 15 cm square template that allowed removal of the complete incorporated and unincorporated litter from the center of the plots. This included all detritus and dead wood less than 2.5 cm in diameter. This material was bagged and taken to the lab for processing. In addition all coarse woody debris greater than 2.5 cm in diameter was measured for length and diameter but only for material that was within the 4.24 m × 4.24 m plot. Woody debris volume was calculated based on the volume of a cylinder. Lastly, tree species was recorded and an estimation of basal area was taken using a variable radius plot protocol and a 2 BAF metric angle gauge to characterize stand structure and composition and to insure there were no differences across sites and treatments [62].Figure 2
Diagram of plot design. Two soil cores were taken at the center and at each corner.


### 2.4. Soil Measurements and Analysis

In the field, soil core depth was measured to calculate total soil volume of the collected cores. In the lab, the wet and dry weights, the percent moisture and percent carbon were determined. While there are difficulties with assessing total carbon losses by using percent carbon, other studies have demonstrated that it can be used as a proxy [63]. Soil was sieved in the lab to 2 mm using the protocol of Bradford et al. [64]. The dried soil samples were sent to the Ecology Lab at the University of Georgia, where carbon content of each sample was determined using a micro-Dumas flash combustion, CHN Analyzer, essentially combusting soil components and measuring carbon via gas (CO2) release [64]. Carbon per unit forest area was calculated with the following equation:
(1)
C
(
Mg
ha
)
=
[
soil
bulk
density
(
g
/
m
3
)
×
soil
depth
(
m
)
×
C
]
×
100
×
10000
m
2
1
ha
×
1
Mg
1000000
g
.See [65]. After initial tests of bulk density in our plots, our measurements corroborated those of Kulmatiski et al. [58]. Using both our initial tests and a penetrometer we found no evidence of differences in bulk densities between our treatment and control for each our sites or across all our sites more generally. We therefore used 0.81, from an average surface soil bulk density for the Yale-Myers Forest [58]. We used the Kulmatiski et al. [58] measure of soil bulk density because it was based on a much more substantial soil survey of the forest soils at Yale-Myers than our study.Litter mass was measured wet and then dried in the laboratory at 80°C for 48 hours [66]. Carbon content was found by converting a ratio of carbon to dried material for litter and, through additional steps, for coarse woody debris. Biomass was assumed to comprise 49.8% of the wood volume for coarse woody debris and 50% of the mass for dry litter for oak-hickory stands in the northeastern US [67]. Specific wood density was estimated at 0.636 [67] in this study. There are lower estimates from US Forest Inventory and Analysis (FIA) data but we selected a higher value given the larger proportion of slower growing denser woods within our forest [68].Based on field observation, decay classes were assigned an average decay class of two for harvested plots younger than six years, a decay class of three for harvested plots six years or greater, and a decay class of three for control sites. Our approach was a simplified protocol using the Brown et al. [65] decay classification and decay coefficients from Waddell [69]. The equation to convert volume of coarse woody debris per plot to carbon mass per area was
(2)
C
(
g
)
area
(
m
2
)
=
[
volume
,
(
m
3
)
*
specific
gravity
1
*
biomass
ratio
*
decay
class
(
m
3
)
]
.


### 2.5. Data Analysis

The fifteen sites were treated as paired replicates (control and shelterwood harvest) using a two-way ANOVA using R, version 2.12.2 [70]. The model comprised treatment (harvest, control), sites over time (1–20 years), and an interaction term (treatment × site). Variables analyzed included (i) percent soil carbon in the top 10 cm; (ii) carbon in the litter and forest floor; (iii) estimated carbon in coarse woody debris; and (iv) all pools combined. All soil carbon and litter data was nonnormally distributed and was log-transformed while coarse woody data was cube-root-transformed prior to analysis. We used Tukey’s studentizedt-test to compare levels of significance among treatments for all variables we tested. In addition, regression analysis was used to explore change in carbon over time to determine if a relationship exists in surface and soil carbon and time since harvest treatment.


## 2.1. Site Description of the Study Area

This study was conducted at the 3173-ha Yale-Myers Forest near Eastford, CT, in the northeastern United States (41°57′N, 72°07′W). The Yale-Myers Forest is a research and demonstration forest managed by the School of Forestry and Environmental Studies and owned by Yale University. The forest overstory is predominantly oak-maple-pine. The climate is characterized as cool temperate and humid. Mean annual summer temperature is 20°C, and mean annual winter temperature is –4°C. Precipitation is evenly distributed throughout the year with a mean of 110 cm [55, 56]. This is a moist temperate forest underlain by a metamorphic bedrock of schist/gneiss and overlain with soils that originated from till and fluvial sediments of the last Wisconsin glacial period 20,000 b.p. [55]. The soils are mostly coarse-loamy, mesic Typic Dystrudepts [57]. The topography is characterized by undulating hills with broad ridges and narrow valleys. The elevation ranges from a low point of 170 m above sea level to a high point of 300 m above sea level [56].Starting in the early 18th century, the forest, as with many areas of New England, had a long period of agricultural use with land used for various kinds of improved and unimproved pasture, tilled for crops, or relict forest patches that were used as woodlots. Since the 1820s this land has largely been reforested naturally [55].


## 2.2. Treatment Description

All stands for our study were on glacial till soils of ablation origin and comprising the Charlton-Chatfield series [57]. This soil series is a dominant one for the oak-hardwood of southern New England making our study results generalizable to a much larger region than the forest. All study stands were originally cleared for marginal pasture but were never plowed. After abandonment in the 1850s these sites transformed to eastern white pine (Pinus strobus L.) and were then cutover in the early 1900s for the pine boxwood industry. After cutting the pine, understory hardwood regeneration comprising oaks, hickories, birches, and maples was released. Today this second-growth hardwood forest varies in age between 90 and 110 years of age. This age class and forest history are very representative of the southern New England area.Using available GIS maps of soil type and harvests within Yale-Myers, a map of shelterwood harvests and reserves of similar land use history (as described above) overlain with the Charlton-Chatfield soil type was created to identify sites that had suitable treated and control plots adjacent to each other (Figure1). Controls were therefore selected as an uncut reference of similar cover, soil, land use history, and topographic position adjacent to each shelterwood harvest treatment. We did paired comparisons of treated versus control plots to determine if there is a statistical difference (at the 0.05 level or greater) between shelterwood harvest and no harvest. Harvested sites were selected to represent a range of ages from 0 to 20 years since harvest and so can be regarded as surrogates to investigate effects of age.Figure 1
Locations of each of the 15 shelterwood treatment and adjacent control sites in the Yale-Myers Forest.At the Yale-Myers Forest, shelterwood harvests reduce the canopy in a single cutting from a fully stocked basal area of about 37 m2/ha to about 11 m2/ha, largely comprising about forty 50–60 cm dbh well-spaced trees per ha. The aim is to maintain species diversity among regenerating trees and seed sources by insuring the heavy seeded masting trees (oaks, hickories, and beech) adequately disperse seeds and that the partial shade of the canopy will moderate the harsh open conditions during germination and seedling establishment. Seven to eight years later half the standing trees (basal area = 5 m2/ha) are removed in a second and final cut leaving the rest as irregular structures to grow within the new regenerating stand. Harvest operations are generally undertaken during winter months, when the ground is frozen and covered with snow. Skid trails and landings are designated ahead of operations, often utilizing previous trails and historic roads.


## 2.3. Sampling Design

We selected fifteen study sites, comprising pairs of treated and control plots. Based on a prior study on soils within this forest, the range was expected to lie between 1.0 and 7.5 kg m2 of carbon in the top 30 cm [58], with a world average around 11.7 kg m2 of carbon for all soils to 100 cm depth [59]. This study estimated that northeastern temperate soils, less carbon-rich than boreal soils in Canada, would have a spread toward the lower end of that range. The sample size for each treatment block was selected with the University of Iowa’s online sample size calculator [60] with settings for two-way ANOVA. Using this estimation, a sample of five and three plots was determined for the treated and control sites, respectively. Individual plots on each site were selected using a random number generator and a GPS coordinate grid.Plots were centered at the GPS coordinates. Plot corners were defined as three meters north, east, south, and west of the plot center. Two 2.54 cm wide soil cores were taken at each point (center, north, south, east, and west) and pooled into one plot sample (Figure2). Because the upper soil horizon has more carbon than other increments of the mineral soil at deeper levels [58] and because most of the observable change 20 years after harvest was expected to be in the upper 10 cm [58, 61], we cored to a depth of 10 cm at each core point from the mineral soil surface and below the litter and organic soil horizon. This meant that we carefully scraped of the organic horizon to carefully identify the mineral soil surface layer. The forest floor organic horizon was sampled using a 15 cm × 15 cm square template that allowed removal of the complete incorporated and unincorporated litter from the center of the plots. This included all detritus and dead wood less than 2.5 cm in diameter. This material was bagged and taken to the lab for processing. In addition all coarse woody debris greater than 2.5 cm in diameter was measured for length and diameter but only for material that was within the 4.24 m × 4.24 m plot. Woody debris volume was calculated based on the volume of a cylinder. Lastly, tree species was recorded and an estimation of basal area was taken using a variable radius plot protocol and a 2 BAF metric angle gauge to characterize stand structure and composition and to insure there were no differences across sites and treatments [62].Figure 2
Diagram of plot design. Two soil cores were taken at the center and at each corner.


## 2.4. Soil Measurements and Analysis

In the field, soil core depth was measured to calculate total soil volume of the collected cores. In the lab, the wet and dry weights, the percent moisture and percent carbon were determined. While there are difficulties with assessing total carbon losses by using percent carbon, other studies have demonstrated that it can be used as a proxy [63]. Soil was sieved in the lab to 2 mm using the protocol of Bradford et al. [64]. The dried soil samples were sent to the Ecology Lab at the University of Georgia, where carbon content of each sample was determined using a micro-Dumas flash combustion, CHN Analyzer, essentially combusting soil components and measuring carbon via gas (CO2) release [64]. Carbon per unit forest area was calculated with the following equation:
(1)
C
(
Mg
ha
)
=
[
soil
bulk
density
(
g
/
m
3
)
×
soil
depth
(
m
)
×
C
]
×
100
×
10000
m
2
1
ha
×
1
Mg
1000000
g
.See [65]. After initial tests of bulk density in our plots, our measurements corroborated those of Kulmatiski et al. [58]. Using both our initial tests and a penetrometer we found no evidence of differences in bulk densities between our treatment and control for each our sites or across all our sites more generally. We therefore used 0.81, from an average surface soil bulk density for the Yale-Myers Forest [58]. We used the Kulmatiski et al. [58] measure of soil bulk density because it was based on a much more substantial soil survey of the forest soils at Yale-Myers than our study.Litter mass was measured wet and then dried in the laboratory at 80°C for 48 hours [66]. Carbon content was found by converting a ratio of carbon to dried material for litter and, through additional steps, for coarse woody debris. Biomass was assumed to comprise 49.8% of the wood volume for coarse woody debris and 50% of the mass for dry litter for oak-hickory stands in the northeastern US [67]. Specific wood density was estimated at 0.636 [67] in this study. There are lower estimates from US Forest Inventory and Analysis (FIA) data but we selected a higher value given the larger proportion of slower growing denser woods within our forest [68].Based on field observation, decay classes were assigned an average decay class of two for harvested plots younger than six years, a decay class of three for harvested plots six years or greater, and a decay class of three for control sites. Our approach was a simplified protocol using the Brown et al. [65] decay classification and decay coefficients from Waddell [69]. The equation to convert volume of coarse woody debris per plot to carbon mass per area was
(2)
C
(
g
)
area
(
m
2
)
=
[
volume
,
(
m
3
)
*
specific
gravity
1
*
biomass
ratio
*
decay
class
(
m
3
)
]
.


## 2.5. Data Analysis

The fifteen sites were treated as paired replicates (control and shelterwood harvest) using a two-way ANOVA using R, version 2.12.2 [70]. The model comprised treatment (harvest, control), sites over time (1–20 years), and an interaction term (treatment × site). Variables analyzed included (i) percent soil carbon in the top 10 cm; (ii) carbon in the litter and forest floor; (iii) estimated carbon in coarse woody debris; and (iv) all pools combined. All soil carbon and litter data was nonnormally distributed and was log-transformed while coarse woody data was cube-root-transformed prior to analysis. We used Tukey’s studentizedt-test to compare levels of significance among treatments for all variables we tested. In addition, regression analysis was used to explore change in carbon over time to determine if a relationship exists in surface and soil carbon and time since harvest treatment.


## 3. Results

### 3.1. Comparisons in Carbon Pools between Shelterwood Harvest and Control Sites

Harvested sites were shown to have significantly lower amounts of surface soil carbon (carbon in the top 30 cm of the soil) (6.01 kg m2) than control sites (8.36 kg m2) that were left intact and there was significant difference among sites and in interaction between treatments and sites (Figure 3, Table 1).Table 1
ANOVA of soil, litter, and coarse woody debris (CWD) carbon mass per area by treatment.
df
Surface soil
Litter
CWD
Total
F value
P
F value
P
F value
P
F value
P
Treatment
1
28.273
0.000***
5.5742
0.020*
25.314
0.000***
0.024
0.877
Site
14
5.188
0.000***
1.5889
0.099
3.868
0.000***
2.977
0.001***
Treatment × site
14
4.649
0.000***
1.2263
0.272
2.041
0.024*
3.887
0.000***
Significance codes are as follows: ***0.001, **0.01, and *0.05. Treatment refers to shelterwood harvest versus control. Site refers to replicate sites examined that span a 20-year period since harvest.Mass per area of each carbon pool by treatment. Tukey pairedt-tests between the variables indicate level of significance at 0.05 level whereby a
>
b. Different letters indicate significant difference.
(a)
(b)
(c)Shelterwood harvest also had an impact on coarse woody debris. Harvest treatments had a significantly greater amount of carbon in coarse woody debris (3.16 kg m2) as compared to unharvested controls (0.62 kg m2) (Figure 3, Table 1), but again significant differences were shown among sites and in interactions between site and treatment. Total carbon in litter per unit area was found to have a significant difference between harvested (0.59 kg m2) and unharvested controls (0.47 kg m2) (Figure 3, Table 1). But no difference was shown among sites and interactions between site and treatment.When surface soil carbon was combined with carbon in litter and coarse woody debris on a mass per unit area, no significant difference between shelterwood harvested treatments and unharvested controls could be demonstrated (Figure3, Table 1). However, there was a significant difference among sites and in an interaction between treatment and site, indicating that sites (as surrogates for time since shelterwood cut) not only had different starting carbon stocks but also responded differently to harvests over time.


### 3.2. Comparisons in Carbon Pools over Time for Shelterwood Harvest and Control Sites

Downward trends in carbon content for sites over time since harvest were seen for all three pools but with a weakR
2 value when analyzed with a linear regression model (Figure 4). However, P values were significant for the harvested trend lines for both the coarse woody debris and total pooled carbon, indicating a downward trend in both of these pools over time and supporting the contention that time since harvest is the interacting effect in the results for the ANOVA. Coarse woody debris was the most significant with a decline from 4.5 kg m2 to 0.40 kg m2 over the twenty-year period suggesting that it had reached control baseline levels of unharvested sites within this time frame. Soil carbon and litter suggested a weak downward trend over time relative to the control but both were not significant.Regressions depicting mass per area change of each carbon pool over time for control and shelterwood harvested sites: (a) soil carbon change, (b) litter carbon change, (c) coarse woody debris carbon change, and (d) total carbon change for three combined pools.
(a)
(b)
(c)
(d)


## 3.1. Comparisons in Carbon Pools between Shelterwood Harvest and Control Sites

Harvested sites were shown to have significantly lower amounts of surface soil carbon (carbon in the top 30 cm of the soil) (6.01 kg m2) than control sites (8.36 kg m2) that were left intact and there was significant difference among sites and in interaction between treatments and sites (Figure 3, Table 1).Table 1
ANOVA of soil, litter, and coarse woody debris (CWD) carbon mass per area by treatment.
df
Surface soil
Litter
CWD
Total
F value
P
F value
P
F value
P
F value
P
Treatment
1
28.273
0.000***
5.5742
0.020*
25.314
0.000***
0.024
0.877
Site
14
5.188
0.000***
1.5889
0.099
3.868
0.000***
2.977
0.001***
Treatment × site
14
4.649
0.000***
1.2263
0.272
2.041
0.024*
3.887
0.000***
Significance codes are as follows: ***0.001, **0.01, and *0.05. Treatment refers to shelterwood harvest versus control. Site refers to replicate sites examined that span a 20-year period since harvest.Mass per area of each carbon pool by treatment. Tukey pairedt-tests between the variables indicate level of significance at 0.05 level whereby a
>
b. Different letters indicate significant difference.
(a)
(b)
(c)Shelterwood harvest also had an impact on coarse woody debris. Harvest treatments had a significantly greater amount of carbon in coarse woody debris (3.16 kg m2) as compared to unharvested controls (0.62 kg m2) (Figure 3, Table 1), but again significant differences were shown among sites and in interactions between site and treatment. Total carbon in litter per unit area was found to have a significant difference between harvested (0.59 kg m2) and unharvested controls (0.47 kg m2) (Figure 3, Table 1). But no difference was shown among sites and interactions between site and treatment.When surface soil carbon was combined with carbon in litter and coarse woody debris on a mass per unit area, no significant difference between shelterwood harvested treatments and unharvested controls could be demonstrated (Figure3, Table 1). However, there was a significant difference among sites and in an interaction between treatment and site, indicating that sites (as surrogates for time since shelterwood cut) not only had different starting carbon stocks but also responded differently to harvests over time.


## 3.2. Comparisons in Carbon Pools over Time for Shelterwood Harvest and Control Sites

Downward trends in carbon content for sites over time since harvest were seen for all three pools but with a weakR
2 value when analyzed with a linear regression model (Figure 4). However, P values were significant for the harvested trend lines for both the coarse woody debris and total pooled carbon, indicating a downward trend in both of these pools over time and supporting the contention that time since harvest is the interacting effect in the results for the ANOVA. Coarse woody debris was the most significant with a decline from 4.5 kg m2 to 0.40 kg m2 over the twenty-year period suggesting that it had reached control baseline levels of unharvested sites within this time frame. Soil carbon and litter suggested a weak downward trend over time relative to the control but both were not significant.Regressions depicting mass per area change of each carbon pool over time for control and shelterwood harvested sites: (a) soil carbon change, (b) litter carbon change, (c) coarse woody debris carbon change, and (d) total carbon change for three combined pools.
(a)
(b)
(c)
(d)


## 4. Discussion

The results indicate that soil carbon pools significantly declined after shelterwood harvest but that total soil carbon mass per area in the top 10 cm showed no significant difference as compared to uncut controls. Carbon mass per area of coarse woody debris was found to increase significantly after harvest, as expected. Furthermore, the difference between more recently cut sites and older sites appeared to narrow, as expected with the decomposition of woody material. It could be concluded that over the 20-year period of shelterwood harvests differences in amount of coarse woody debris could be related to changes in site treatment (e.g., prescribed burning, whole-tree chipping, etc.). We can attest that no changes occurred. All site treatments were the same with all stems and tops less than 20 cm dbh being left purposefully scattered across the harvest sites. However, this trend may have been more accurately depicted if woody debris volume had been individually classified to a decay class rather than pooled to particular years since harvest. Litter carbon was not shown to be significantly different between shelterwood harvest and control plots, although it was shown to be generally higher than US Forest Service estimates for postharvest oak-hickory forests in the northeast [49]. Coarse woody debris, however, was shown to be less for those forest service estimates [49], which in turn were less than those found by Covington [30]. The difference between our results and those of Covington [30] and Smith et al. [49] could be attributed to method of harvest (shelterwood versus clearcut) or to differences in measurement protocols and sampling design.Our study also showed higher on average soil carbon content as compared to a study at the same forest by Kulmatiski et al. [71] whose objective was to quantify total soil carbon at depth. The Kulmatiski et al. [71] study was forest wide, irrespective of soil type. On average, it also sampled soil three times deeper than our study. Our soil cores went to a 10 cm depth because our focus was the effects of harvesting on carbon in surface horizons. They also found topographic position accounted for 18% of soil carbon variation, with wetter lower lying areas containing more carbon than drier upland soils. Furthermore, carbon content, bulk density, and depth are correlated [72] and further research needs to be done to measure each of these variables to verify and refine more accurate calculations for soil carbon content unique to each site.The difference in total carbon stocks between control and harvested sites, however, was not found to be significant when pooling together all sites of varying time since harvest. The potential reason for the lack of significance is the dynamic change among the different carbon pools, some increasing and some decreasing over time. While coarse woody debris represented a smaller carbon pool than mineral soil carbon, the large increase counterbalanced the soil carbon decrease. Litter was the smallest carbon pool and also experienced no significant change. However, when examining trends over the 20-year time period total carbon among all pools (soil, litter, and debris) showed a significant decline. This was primarily due to coarse woody debris which showed the only significant downward trend over time. The trend line for mineral soil carbon and litter suggested a downward trend over time relative to the control but was not significant. Obtaining data from older harvests may have strengthened this relationship but no sites exist.Overall, these results indicate that these three pools are a significant portion of total carbon stocks but that the changes that may occur due to harvest represent a relatively small portion of total site carbon stocks. Coarse woody debris showed the strongest significant change with shelterwood harvest, but it was a smaller pool than soil carbon. Surface soil carbon significantly declined in the top 10 cm, but further studies that extend over a longer time period will be important to do.


## 5. Conclusions

Taking all sites together, total pools of carbon (litter, woody debris, and surface soil carbon) showed no significant decline with harvest primarily because declines in surface soil carbon and litter were more than made up for by increases in coarse woody debris. However, when comparing sites with time since shelterwood harvest, there are significant declines in total carbon. Our results show this is strongly driven by decline in coarse woody debris. There are also significant declines in surface soil carbon with harvest, but there is considerable variability between sites, and this trend is weak with time since harvest. These results suggest that for the purposes of forest carbon accounting auditors do not need to monitor changes in below-ground carbon stocks with shelterwood harvests for this region. We believe there is not enough of a significant decline or change in carbon to merit the time and cost in measurement.



---

*Source: 527236-2014-05-06.xml*