Elsevier

Forest Ecology and Management

Volume 319, 1 May 2014, Pages 150-160
Forest Ecology and Management

Spatial heterogeneity of soil CO2 efflux after harvest and prescribed fire in a California mixed conifer forest

https://doi.org/10.1016/j.foreco.2014.02.012Get rights and content

Highlights

  • We assessed the effects of forest management on spatial variability of soil CO2 efflux.

  • We compared a control, a harvested and a burned stand in a mixed conifer forest.

  • Spatial variability was increased by disturbance, especially by the clear-cut.

  • Harvest reduced the ability to quantify soil CO2 flux and detect treatment effects.

Abstract

Spatial variability is a key factor when quantifying soil CO2 efflux and punctual measurements need to be extended to larger stand, ecosystem, or regional scales. Spatial variation also affects comparisons among ecosystems, as when quantifying effects of disturbances on ecosystem carbon dynamics. However, spatial variability of soil CO2 efflux is still unknown and difficult to predict. We quantified the effects of silvicultural practices (prescribed fire and harvesting) on spatial variability of soil CO2 efflux in a mixed conifer forest from the central Sierra Nevada in California, USA. Soil CO2 efflux was measured using a portable chamber system, on 20–29 locations in four treatment sites: an untreated control, a prescribed fire site (burned in 2002 and 2009) and two clear cut sites harvested in 2010. In one of the harvested sites the soils were mechanically ripped to reduce soil compaction, a common practice done on industrial timber forest lands in the Sierra Nevada. Results showed that disturbance increased spatial variability of soil CO2 efflux. Coefficient of variability increased from an annual average of 32% at the control site to 37% at the burned site, and 49–51% at the harvested sites (without and with soil ripping, respectively), mirroring post-disturbance increases in spatial variability of soil temperature and soil water content. Because of the post-harvest increase in spatial variability, the ability to detect differences became lower, and the number of samples needed to obtain a value representative of the full population mean (within a 10% range) increased by 100%, from 60 to 120 samples. To reduce uncertainty in our soil CO2 efflux treatment estimates, more than 10–15 randomly selected locations per study site were necessary. Spatial variability of soil CO2 efflux at our sites was not affected by distance between measurement locations, was correlated to fine root and litter biomass at the control site, negatively correlated to soil bulk density at the fire site, and un-correlated to soil temperature and water content at all sites. The increase of spatial variability in soil CO2 efflux after disturbance and the requirement for a sufficient number of measurement locations should be considered when quantifying carbon dynamics of disturbed ecosystems, or assessing effects of different forest management practices.

Introduction

Soil-surface CO2 efflux, commonly referred to as soil respiration, is one of the main carbon fluxes in the global carbon cycle (IPCC, 1996), and the largest carbon source to the atmosphere (60–90% of total ecosystem respiration; Liang et al., 2004) in forest ecosystems. Forest ecosystems act generally as sink of carbon, however disturbances can switch forests from carbon sinks to sources (Amiro et al., 2010, Dore et al., 2012, Thornton et al., 2002).

Disturbances can vary in type and intensity. They can range from high intensity events, such as land-use changes, stand replacing fires or clear-cut harvest, to low intensity events, such as thinning and low intensity fires. Thus forest management practices, that include removal of biomass by mechanical methods or by prescribed fire, represent a disturbance to forest ecosystems. The strength and persistence of disturbance effects are greatly determined by their impacts on soil CO2 efflux, because even if vegetation recovers promptly after the disturbance, the carbon lost in decomposition of old and new material can exceed the newly restored carbon sink (Restaino and Peterson, 2013).

High spatial and temporal variability in soil CO2 efflux has been found in numerous studies (Han et al., 2007, Ngao et al., 2012, Tedeschi et al., 2006, Vincent et al., 2006) and control of soil CO2 efflux has been attributed mainly to soil temperature and soil water content (Hanson et al., 1993, Raich and Schlesinger, 1992). While these variables can explain most of the temporal variability of soil CO2 efflux, they are often unable to explain its spatial variability (Tedeschi et al., 2006, Xu and Qi, 2001, Yim et al., 2003). Temporal variation of soil CO2 efflux is relatively easy to quantify, especially with the development of techniques based on continuous measurements made using chambers (Edwards and Riggs, 2003), below canopy eddy covariance (Law et al., 1999) and CO2 soil profiles (Jassal et al., 2005). However, spatial variability is still understudied and mostly uncertain (Ngao et al., 2012). Spatial variability of soil CO2 efflux has been linked to several biotic and abiotic factors, such as species composition, leaf area, fine root biomass, litter depth, soil bulk density, and soil carbon content (Ngao et al., 2012). Nonetheless, the contribution of these factors is highly variable across ecosystems, sites, and seasonally within the same sites. This is because soil CO2 efflux is the result of both heterotrophic and autotrophic processes, and these processes are controlled independently by different factors (Tang and Baldocchi, 2005).

Improved understanding of the mechanisms and quantification of the spatial heterogeneity of soil CO2 efflux is essential to scale up from point measurements to the stand and all the way to the global scale, or to verify and/or integrate flux measurements obtained using different techniques (Dore et al., 2003). For example, when using eddy covariance, efflux of CO2 from soil (though scaled up to the same spatial scale of eddy covariance) can replace low quality night ecosystem respiration measurements (Wohlfahrt et al., 2005), explain heterogeneity in the footprint of measured ecosystem fluxes (Ngao et al., 2012), and help partition the respiratory flux between aboveground and soil fluxes (Baldocchi, 2003). Still, determination of spatial variability is difficult and methods are limited (Tang and Baldocchi, 2005). High spatial variability requires a high number of samples to be taken to obtain meaningful results, especially when studies are aimed to detect differences among ecosystems or treatments. Common use of a relatively small number of measurements could explain in part the lack of agreement among studies quantifying effects of disturbances on soil CO2 efflux (Kobziar and Stephens, 2006).

In addition to the high number of replicates needed, results are confounded by practical difficulties in the measurement techniques, such as the disturbance of the soil–air interface, effect of leakage, effect of pressure on diffusivity of CO2 from soil, and the time and labor needed to take measurements in several locations. Also, because most of the factors controlling soil CO2 efflux lie underground and are difficult to quantify, it is not possible to visually estimate spatial variability and easily evaluate the number and position of soil CO2 efflux measurement locations.

In our study we analyzed the effect of forest management practices on spatial variability of soil CO2 efflux. We compared four different treatment types: un-manipulated control, prescribed fire, clear cut harvest, and clear cut harvest followed by mechanical soil ripping (sub-soiling). Clear cut harvests intensely disturbed the ecosystem by totally removing aboveground biomass and affecting the soil during mechanical operations, prescribed fire was a less intense disturbance. Our aim was to: (1) quantify and characterize spatial variability of soil CO2 efflux and its temporal changes in a mixed conifer forest in the central Sierra Nevada of California; (2) determine if spatial variability of soil CO2 efflux was affected by the most commonly used fuel treatments and commercial harvesting methods used in these mixed conifer forests; (3) and to analyze the effects of samples size on soil CO2 efflux estimates. To explore this we compared uncertainty of soil CO2 efflux estimates using a different number of random subsamples compared to the full measured dataset. Second, we developed a protocol aimed to select a smallest number of soil CO2 efflux measurement locations among the locations measured initially. Even if our protocol implies an a priori assessment of the spatial variability of soil CO2 efflux, it could be useful for long term monitoring of fewer selected locations to complement/validate continuous soil CO2 efflux or eddy covariance. Alternatively, it could be used when selecting the location for a permanent, continuous soil CO2 efflux system, because these systems usually have a low number of chambers and thus a high risk of measuring locations that do not represent a larger area.

We had the opportunity to characterize the spatial and temporal variability of soil CO2 efflux in an area subject to a range of disturbance severity but with generally the same climate, vegetation, management history, and soil type. In addition, the treatments created conditions ranging from the high autotrophic contribution of an undisturbed, dense forest, to the single heterotrophic contribution of a clear cut harvested stand.

Section snippets

Methods

This study was conducted at Blodgett Forest (38°54′N, 120°39′W), a University of California Research Station in the central Sierra Nevada near Georgetown, California, that is actively managed as a commercial timberland. Total annual precipitation in the mixed-conifer forest located between 1100 and 1410 m a.s.l. averages about 1600 mm, falling between September and May, and almost absent in the summer (Stephens and Collins, 2004). The average minimum daily temperature in January is 0.6 °C and the

Results

Soil CO2 efflux varied in space, and this spatial variability differed among treatments (p < 0.001) and in time. Coefficient of variation (CV) was lowest at the UND site, was slightly higher at the FIRE site and highest at the HARV sites, mirroring the intensity of disturbance (Fig. 3). CV in soil CO2 efflux was highest at the beginning of both spring and fall.

We can be confident that the higher CV at the HARV site was not caused by a difference in the site characteristics, such as a different

Discussion

In this study we quantified the effects of silvicultural practices (prescribed fire and harvesting) on spatial variability of soil CO2 efflux in a mixed conifer forest in the central Sierra Nevada. Spatial variability increased after disturbance, particularly after harvest, despite the harvested site being less complex and heterogeneous than the undisturbed stand. The harvested areas had minimal vegetation cover, consisting of one year old conifer seedling planted on a 2.5 m grid, and thus

Acknowledgements

This project was funded by The US Joint Fire Sciences Program. We thank Rob York, Jen York, and Blodgett staff for collecting data and maintaining instruments during this project.

References (46)

  • J.C. Restaino et al.

    Wildfire and fuel treatment effects on forest carbon dynamics in the western United States

    For. Ecol. Manage.

    (2013)
  • S.L. Stephens et al.

    Experimental fuel treatment impacts on forest structure, potential fire behavior, and predicted tree mortality in a mixed conifer forest

    Forest Ecol. Manag.

    (2005)
  • B.W. Sullivan et al.

    Thinning reduces soil carbon dioxide but not methane flux from southwestern USA ponderosa pine forests

    Forest Ecol. Manag.

    (2008)
  • P.E. Thornton et al.

    Modeling and measuring the effects of disturbance history and climate on carbon and water budgets in evergreen needleleaf forests

    Agric. For. Meteorol.

    (2002)
  • G. Vincent et al.

    Spatial and seasonal variations in soil respiration in a temperate deciduous forest with fluctuating water table

    Soil Biol. Biochem.

    (2006)
  • G. Wohlfahrt et al.

    Quantifying nighttime ecosystem respiration of a meadow using eddy covariance, chambers and modelling

    Agric. For. Meteorol.

    (2005)
  • M.H. Yim et al.

    Spatial variability of soil respiration in a larch plantation: estimation of the number of sampling points required

    Forest Ecol. Manag.

    (2003)
  • B.D. Amiro et al.

    Ecosystem carbon dioxide fluxes after disturbance in forests of North America

    J. Geophys. Res.

    (2010)
  • D.D. Baldocchi

    Assessing the eddy covariance technique for evaluating carbon dioxide exchange rates of ecosystems: past, present and future

    Global Change Biol.

    (2003)
  • Y. Chen

    New approaches for calculating Moran’s index of spatial autocorrelation

    PLoS ONE

    (2013)
  • A. Concilio et al.

    Soil respiration response to prescribed burning and thinning in mixed-conifer and hardwood forests

    Can. J. For. Res.

    (2005)
  • S. Dore et al.

    Cross validation of open-top chamber and eddy covariance measurements of ecosystem CO2 exchange in a Florida scrub-oak ecosystem

    Global Change Biol.

    (2003)
  • S. Dore et al.

    Recovery of ponderosa pine ecosystem carbon and water fluxes from thinning and stand replacing fire

    Global Change Biol.

    (2012)
  • Cited by (31)

    • Spatial heterogeneity of soil carbon exchanges and their drivers in a boreal forest

      2022, Science of the Total Environment
      Citation Excerpt :

      The reason is not always well understood, demonstrating that there is still a lack of understanding of how the C balance of boreal forests is regulated. Soil respiration has two main components: 1) microbial decomposition of SOM and plant litter, and 2) respiration by living plant roots (Basiliko et al., 2012; Dore et al., 2014). Microbial growth and activity are controlled by multiple abiotic factors, such as soil moisture, pH, temperature, SOM quality, and nutrient availability (Bengtson et al., 2005; Fenn et al., 2010; Heinemeyer et al., 2007; Yu et al., 2015).

    • Harvesting intensity and tree species affect soil respiration in uneven-aged Dinaric forest stands

      2021, Forest Ecology and Management
      Citation Excerpt :

      Seasonal variations in Rs, observed in almost all ecosystems, have often been associated with changes in temperature, moisture, photosynthetic production, root growth and their combination (Yan et al., 2011; Acosta et al., 2018; Makita et al., 2018; Zhang et al., 2018). On the other hand, the spatial variability in Rs results from a large variability in soil physical properties, soil chemistry, fine-root biomass, fungi and bacteria, nutrient availability and others (disturbance and weathering) (Hanson et al., 1993; Fang et al., 1998; La Scala et al., 2000; Xu and Qi, 2001; Merbold et al., 2011; Allaire et al., 2012; Dore et al., 2014; De Carlo et al., 2019; D'Andrea et al., 2020). The differences in Rs among stands result also from the prevailing tree species, age and management practices (Peng et al., 2008; Akburak and Makineci, 2013; Wang et al., 2013).

    View all citing articles on Scopus
    View full text