NC State
BioResources
Aydın, M., and Rages, A. A. (2024). “Soil respiration and organic carbon changes along a chronosequence of Pinus nigra forest stands,” BioResources 19(3), 6095-6119.

Abstract

Understanding the trajectory of changes in soil respiration (Rs) and soil organic carbon (SOC) with stand ages of the black pine (Pinus nigra Arnold) forest is essential for forest management and carbon budget estimates. In this research, changes of Rs and SOC were studied with respect to stand age in a chronosequence of three age classes of P. nigra plantations consisting of young (0 to 10-year-olds), middle-aged (11- to 20-year-olds), and pre-mature (35- to 45-year-olds) forest stands. Rs rates, soil temperature, and soil moisture were measured using an automated dynamic survey chamber (Li-8100A) for a year, encompassing summer, fall, winter, and spring seasons. Mean Rs significantly increased from young- to middle-aged and then stabilized, with effluxes ranging from 2.46 to 2.94 µmol CO2 m-2 s-1. Forest litter significantly increased with stand age, but not the SOC in the mineral soil layers. The Rs showed a positive correlation with soil temperature (0.77) and air temperature (0.75) but not with soil moisture (-0.43). The present results highlight the importance of stand age in assessing carbon budget and provide essential information for forest managers and stakeholders in evaluating the potential of P. nigra forests as tools for carbon sequestration and mitigating global warming impacts.


Download PDF

Full Article

Soil Respiration and Organic Carbon Changes along a Chronosequence of Pinus nigra Forest Stands

Miraç Aydın * and Ashraf Anwar Rages

Understanding the trajectory of changes in soil respiration (Rs) and soil organic carbon (SOC) with stand ages of the black pine (Pinus nigra Arnold) forest is essential for forest management and carbon budget estimates. In this research, changes of Rs and SOC were studied with respect to stand age in a chronosequence of three age classes of P. nigra plantations consisting of young (0 to 10-year-olds), middle-aged (11- to 20-year-olds), and pre-mature (35- to 45-year-olds) forest stands. Rs rates, soil temperature, and soil moisture were measured using an automated dynamic survey chamber (Li-8100A) for a year, encompassing summer, fall, winter, and spring seasons. Mean Rs significantly increased from young- to middle-aged and then stabilized, with effluxes ranging from 2.46 to 2.94 µmol CO2 m-2 s-1. Forest litter significantly increased with stand age, but not the SOC in the mineral soil layers. The Rs showed a positive correlation with soil temperature (0.77) and air temperature (0.75) but not with soil moisture (-0.43). The present results highlight the importance of stand age in assessing carbon budget and provide essential information for forest managers and stakeholders in evaluating the potential of P. nigra forests as tools for carbon sequestration and mitigating global warming impacts.

DOI: 10.15376/biores.19.3.6095-6119

Keywords: Stand age; CO2 emissions; Soil temperature; Soil moisture; Black pine; Dynamic chamber

Contact information: Faculty of Forestry, Watershed Management Department, Kastamonu University, 37150, Kastamonu, Türkiye; *Corresponding author: maydin@kastamonu.edu.tr

INTRODUCTION

Under the Kyoto Protocol of the UNFCC, carbon sequestration in the afforestation and reforestation projects is one of the prescribed options to enhance the removal of atmospheric greenhouse gases (IPCC 2001). Article 12.5 of this Protocol requires that the certified emission reductions under the clean development mechanism (CDM) must be “environmentally additional” to what would have occurred without the project activities based on the tangible and measurable quantity of GHGs relative to the baseline (Baumert 2000). For this reason, it is essential to quantify the baseline GHG emissions accurately to prevent “over-crediting,” which could be done by estimating and validating the associated emissions and removals of GHGs (Pearson et al. 2013). Although there is no rigid rule for carbon sequestration accounting, the carbon budget method, which includes measuring the C input and output fluxes to and from the ecosystem, is considered the most accurate method of accounting for C sinks (Saint Andre et al. 2007).

Quantifying carbon budgets of a given forest ecosystem requires quantification of the different components of the C cycle, including the C inputs in the above- and below-ground biomass and soil organic carbon (SOC), as well as the fraction of gross primary production that is lost through plant respiration (Trenberth 2005; Zscheischler et al. 2014; Berg and Sheffield 2018; Xu et al. 2021). Plant respiration, the main pathway through which a large quantity of CO2 returns into the atmosphere, comprises soil respiration (roots and soil organisms) and respiration from the aboveground biomass (foliage, stems, and branches) (Wang et al. 2023; Pacaldo et al. 2024). Soil respiration (Rs) constitutes a large proportion of ecosystem respiration, contributing about 80% to 90% of the total plant respiration rates (Raich and Schlesinger 1992; Schimel et al. 2001). Some authors estimated that the CO2 emissions from the Rs are about 78 to 95 Pg of CO2 emissions back into the atmosphere annually (Bond-Lamberty and Thomson 2010; Hashimoto et al. 2015) and are considered the largest source of CO2 emissions (IPCC 2021; Nissan et al. 2023). Because of the large amount of CO2 emissions from Rs and SOC stocks (Raich and Schlesinger 1992; Schimel et al. 2001; Nissan et al. 2023) any small changes in these components of the C cycle could dramatically alter the C budget of a given ecosystem (Pacaldo et al. 2024).

The stand age of the forest plantation may exert a strong influence on Rs and SOC because as the forest develops and ages, some structural, morphological, and physiological changes occur, which are likely to affect the C cycling and other vital processes in the forest ecosystem (Yu et al. 2014). Many researchers have invested considerable efforts into understanding the age effects of forest ecosystems on Rs and SOC. However, these studies suggest different trajectories of stand age effects on Rs rates: (1) increased (Peichl et al. 2010; Peichl et al. 2014; Song et al. 2017; Yu et al. 2019); (2) decreased (Darenova et al. 2016); (3) no change (Chang et al. 2020); (4) decreased and then increased (Payeur et al. 2012); (5) increased and then stabilized (Uri et al. 2022); (6) and nonlinear response with high variability (Smith et al. 2010; Kukumägi et al. 2017). Similarly, the trajectory of SOC changes is less certain and may vary with stand age and other factors (Chen et al. 2010). Reported chronosequence studies conducted in forest ecosystems revealed inconsistent findings on stand age effects on the SOC: (1) increases with stand age (Li et al. 2013; Cheng et al. 2014; Zhao et al. 2014; Francis-Justine et al. 2015; Song et al. 2017; Smal et al. 2019; Lei et al. 2019; Chen et al. 2020; Zhu et al. 2020; De Marco et al. 2021); (2) decreases with stand age (Cao et al. 2012; Amir et al. 2018; Chen et al. 2010; ); (3) nonlinear response, i.e. a decrease and then increase with stand age (Covington 1981; Chen et al. 2013; Pacaldo et al. 2013), and (4) no significant pattern (Noh et al. 2010; Huang et al. 2021). The inconsistency of SOC changes with stand age may partly be explained by the complicated dynamics in the accumulation and decomposition of soil organic matter (Chen et al. 2013).

However, in the literature and to our knowledge, there has been no investigation conducted that compares and assesses the changes in Rs rates and SOC with stand ages of Pinus nigra stands. The deficiency of data precluded our understanding of the carbon sequestration potentials of different developmental stages of P. nigra stands, which is a widely distributed tree species in Turkiye with an estimated area of about 4.2 million hectares (Sakici et al. 2018; Pacaldo et al. 2024). Determining the trajectory of Rs and SOC changes with stand age is essential for accurate C budget estimates in forest ecosystems. Whether the Rs rates and SOC values significantly change as the forest ages is a critical question this study seeks to investigate. In this investigation, it was hypothesized that the Rs rates and SOC would increase with stand age. The objectives of this study were (1) to assess the Rs and SOC changes with stand age and (2) to determine the environmental factors affecting Rs. Quantifying the magnitude and trajectory of Rs and SOC changes with stand age provides valuable information to advance our understanding of the dynamics of the carbon budget in the P. nigra ecosystem, which is a critical factor in achieving accurate estimates of abated anthropogenic CO2 emissions by the reforestation and afforestation projects under the CDM of the Kyoto Protocol (Baumert 2000).

EXPERIMENTAL

Study Site

A field experiment was conducted in a black pine (Pinus nigra Arnold) forest at Kastamonu City, Türkiye, which is geographically located between 41°22’19.89 “N and 33°44’4.10″E, with a mid-latitude temperate climate under the Köppen classification system (Turkes 2020), and an annual temperature mean of 10 ℃, and annual precipitation mean of 538 mm (Turkish State Meteorological Service 2024). The forest plantation was established simultaneously, but restocking, replanting, and natural regeneration in canopy gaps resulted in the variability of age classes, in which the exact age was determined using an increment borer. The management history, soil properties, topography, and silvicultural treatments of the research site were more or less similar; hence, it is safe to assume homogeneity of the Rs and SOC before establishing the plantation forests. Based on the World Reference Base (WRB) classification system, the site’s soil is dominantly Lithic Leptosol (Özden et al. 2001; Aksoy et al. 2010), overlying calcareous and sedimentary rocks in the advanced stage of weathering. Soil properties of the study site are summarized in Table 1. The study site is a reforestation area containing a continuous track of homogenous and well-managed P. nigra plantations (35- to 45-year-olds) with considerable forest gaps in the periphery and inside the forest stands where the naturally regenerated forests of different age classes (5- to 10-year-olds and 11- to 15-year-olds) can be found. A timber inventory was carried out to determine tree stocking and sizes, which are summarized in Table 2.

Experimental Design

A field experiment in a complete randomized block design (10 m × 10 m) was established across regeneration and plantation forest stands. Three age classes of P. nigra plantations and one control were selected in this study. Age classes were based on the age of stands during the 2023 growing season. These included young forest stands (0- to 10-year-olds), middle-aged stands (11- to 20-year-olds), pre-mature stands (35- to 45-year-olds), and control (treeless undisturbed sites located along forest borders). Each of the three age classes, including the control, was replicated four times, represented by 16 sampling plots (i.e., three age classes + control × 4 replications) distributed across four blocks.

In this experiment, the following parameters were evaluated: Rs, SOC, soil temperature, air temperature, soil moisture, and stand age. The key parameters include the Rs and SOC, which are the main focus of the investigation. These parameters constitute the main components of output and input in the C cycle and C budget estimates; hence, any small changes in these components can dramatically alter the C balance of the ecosystem. Soil temperature, air temperature, and soil moisture affect Rs rates and soil organic matter decomposition.

Soil Respiration Measurement

An automated soil CO2 efflux measurement system (LI- 8100A) equipped with a 10-cm survey chamber (LI-8100-103), soil temperature probe (6000-09TC Omega), and EC-5 soil moisture sensor (Li-COR, NE, USA) was used to measure the Rs, soil temperature, air temperature, and soil moisture. The analyzer unit (LI-8100) houses the infrared gas analyzer (IRGA) and stores the data. The survey chamber is equipped with a pressure vent on its top, alleviating the errors due to differences in pressures between inside the chamber and the ambient environment (Liang et al. 2004; Xu et al. 2006). The base of the chamber is fitted with a rubber seal, which prevents air leakage in and outside the chamber’s headspace during measurements. During measurements, the chamber moves automatically into the soil collar by the control of the analyzer unit, which pumps air into the chamber. Measurement of Rs rates was done by mounting the survey chamber on the top of the soil collar with a total duration of 240 seconds for each measurement, consisting of 30 seconds of equilibration/deadband (i.e., length of time when chamber closes completely and mixes with air before measurement begins), 150 seconds observation length, and 60 seconds of purge time (Pacaldo et al. 2024). Simultaneous measurements were made of soil temperature and soil volumetric moisture content with the Rs at 5 cm soil depth.

In each plot, one polyvinyl chloride (PVC) soil collar (5 cm diameter and 7 cm height) was inserted into 4 cm soil depth, leaving a 3 cm soil collar. The soil collars were inserted a few days before the first measurement to minimize the effects of soil disturbance and artifacts. The Rs was measured for one year, encompassing the four seasons (summer, fall, winter, and spring). Bi-monthly Rs measurements were conducted, except during the winter months (January to March), wherein a monthly measurement was done. In winter, the snowpacks were removed carefully from the soil collars before mounting the survey chamber on the top of the soil collars. The Rs between 10:00 and 16:00 were measured based on previous studies using continuous and unattended Rs measurements in the field (Pacaldo et al. 2014). The annual cumulative Rs was evaluated using the mean values and scaled to metric tons of CO2 per hectare per year (Mt CO2 ha-1yr-1) using a conversion factor 12.59 (Pacaldo et al. 2024). Calculations were done also for the sensitivity of Rs to soil temperature using the Lloyd and Taylor (1994) model based on an exponential relationship,

Rs = Rref expE0T       (1)

Q10 = exp10E0           (2)

where Rref (µmol m-2 s-1) is the basal respiration at the reference temperature, E0 (°C) is the parameter of temperature sensitivity, T (°C) is the soil temperature, and Q10 refers to temperature sensitivity, representing the response of Rs to a 10 °C increase in temperature (Han et al. 2023).

Soil and Forest Litter Sampling

Composite soil samples were collected at two 20 cm increment depths: 0 to 20 cm as topsoil and 21 to 40 cm as subsoil. A metal cylindrical bulk density corer (5 cm dia. x 5 cm ht) was used to collect soil samples in the field. Before collecting soil samples, organic matter (undecomposed and partially decomposed) was removed to the depth of the mineral soil surface. Then, the bulk density soil corer was hammered into the middle section of each soil layer depth.

Composite forest litter and other organic debris on the soil surface (O horizon) were collected in each plot using a forest litter sampler (20 cm dia. × 20 cm ht). The samples were transported to the laboratory and weighed per plot to determine the total green (fresh). The dry matter weight was determined randomly by oven-drying the samples at 105 ℃ to constant weight.

Soil Laboratory Analysis

We measured the soil pH and electrical conductivity (EC) using a Hanna HI 9812-5 pH/EC/TDS/Temperature meter in a 1:2.5 soil-water mixture, and EC was then converted to EC of saturated paste extract using the equation ECe = 4.34 × EC1:2.5 – 0.17 (Sonmez et al. 2008) for sandy soils in Türkiye. We determined the soil organic carbon content using the Loss-On-Ignition method by burning the soil samples at 450 °C for 12 hours. The values of SOM were converted to SOC using a conversion factor of 1.72 (Post et al. 1982; Nelson and Sommers 1983). We used an established procedure for the hydrometer method (Bouyoucos 1962; Gangwar and Baskar 2019) to determine the different proportions of soil particles (clay, silt, and sand). We used a textural triangle to determine the textural classification of the soil.

The physical and chemical analyses revealed that the site’s soils, both upper and subsurface layers, are moderately dense, with values ranging from 1.47 to 1.75 g cm-3 and 1.42 to 1.89 g cm-3, respectively. The soil reactions (pH) are nearly neutral, ranging from 7.53 to 7.65 in the upper 10 cm depth and 7.60 to 7.75 in the subsurface layer. The soil EC is slightly saline, ranging from 3.2 to 4.72 dSm-1 in the upper 10 cm depth and 2.85 to 3.63 dSm-1 in the subsurface layer. Generally, the site has a sandy loam soil texture (Table 1).

Table 1. Soil pH, EC, BD, and Soil Texture at the Upper and Lower Layers of Different Forest Types (Mean ± MSE; n = 4)

Aboveground Biomass Inventory

A tree inventory was carried out to determine the tree density and estimate the size of aboveground biomass. In each stand age, four 20 m x 20 m sampling plots were delineated for a total sampling area of 800 m2 or eight percent (8%) sampling intensity per hectare. In each sampling plot, the diameter breast height (DBH) and height of all trees within the plot’s borders were measured. A steel increment borer was used to determine the tree’s exact age. The biomass was estimated based on the formula CAG =0.010dbh2h, developed by Sakici et al. (2018) for P. nigra in the study site. The results of the inventory are summarized below (Table 2).

Table 2. Characteristics of the Aboveground Standing Biomass of P. nigra Forest in the study site (Mean ± MSE)

*Calculated using the equation: Biomass = 0.020 DBH2 x Height (Sakici et al. 2018)

Statistical Analyses

The Rs rates of different forest types and controls were analyzed using the general linear (PROC GLM) model, in which treatment was considered as the main effect, block as the random effect, and time as a second qualitative factor to test if the Rs vary among different measurement times. The significant differences in Rs and SOC among treatments were tested using one-way analysis of variance (ANOVA) with p < 0.05 considered a significant value. Multiple-wise comparisons were done with Tukey’s test to separate significant differences among treatment means. Pearson correlation analysis determined the relationship between soil respiration, soil temperature, soil moisture, and air temperature. The relationships among soil respiration, soil temperature, air temperature, and soil moisture were analyzed using regression analyses based on collected data throughout the study. Values of r-square, Mallows’ Cp Statistics, Akaike Information Criterion (AIC), and mean standard error (MSE) were used to select the number of independent variables in the multiple regression model (Pacaldo et al. 2024). All statistical analyses were performed using a SAS Statistical Package (SAS 9.1 SAS Institute).

RESULTS

Soil Respiration Rates

Soil respiration rates significantly differed across all sources of variation. The Rs showed highly significant differences among treatment (age classes) (p = 0.009) and time (p < 0.0001). However, the combined effects of time and treatment on Rs rates did not show significant differences (p = 0.089) (Table 3).

 

Table 3. Results of ANOVA Test for Rs Rates Among Treatments, Time, and Interaction Effects Between Treatment and Time (n = 240; p = 0.05)

Mean separation using Tukey’s test revealed that the middle-aged stand (11- to 20-year-olds) exhibited the highest Rs rates, which is significantly different from young stand (0- to 10-year-olds) and the control, but not with the Rs of the 35-45-year-olds (pre-mature) plantation. The Rs rates, ranging from 2.36 to 2.94 µmol m-2 s-1, agreed well with some reported values in literature (e.g., Wiseman and Seiler 2004; Payeur-Poirier et al. 2012; Luan et al. 2012; Pang et al. 2013; Wei et al. 2022; Amarille et al. 2023; Pacaldo and Aydın 2023; Tong et al. 2023). The cumulative annual Rs rates ranged from 29.72 to 37.06 Mt CO2 ha-1 yr-1 (Table 4).

Table 4. Mean Annual Rs and Cumulative Annual Rs Rates of All Age Classes and Control

*Values with the same letters are not significantly different at a 95% confidence level, based on Tukey’s Test

Seasonal Soil Respiration Rates

In the winter season, the Rs rates showed no significant differences in all stand ages, including the control. The middle-aged stand in springtime demonstrated significantly higher Rs than other treatments. From summer to fall periods, the control indicated the highest Rs. However, it did not significantly differ from the middle-aged and pre-mature stands (Fig. 1). Proportional seasonal contributions to the total cumulative annual soil CO2 emissions showed that winter contributed only about 8.84 %. In comparison, summer contributed about twice as much (42.60%) as spring (21.66%) and fall (26.90%) (Fig. 2), which is consistent with some previous reports (e.g., Groffman et al. 2001; Liptzin et al. 2009; Pacaldo et al. 2012, 2024).

Fig. 1. Mean seasonal Rs rates in all age classes and control (n = 16; Mean ± MSE)

Fig. 2. Proportion of soil respiration (Rs) in all stand ages and control (Mt CO2 ha-1 per season) and seasonal contribution of Rs to the cumulative annual CO2 emission rates (%)

Relationship between Soil Respiration, Temperature, and Soil Moisture

Pearson’s correlation analysis showed a strong significant relationship between Rs and soil temperature (0.77) and air temperature (0.75), indicating that it tends to increase with increasing soil and air temperatures. In contrast, the Rs showed a weak negative significant relationship with soil moisture (-0.43), suggesting that the Rs tend to decrease with increasing soil moisture contents. Notably, the soil temperature and air temperature showed a robust positive correlation (0.91), suggesting that the air temperature provides a good approximation of the soil temperature (Pacaldo et al. 2024). The soil moisture demonstrated a stronger negative relationship with the air temperature (-0.58) than with soil temperature (-0.45), indicating that the soil moisture tends to dry faster with an increasing air temperature than with soil temperature (Table 5). The Rs significantly increased in dryer than saturated soil conditions (Figs. 3 and 4). The sensitivity analysis showed higher sensitivity of Rs to temperature below 20 °C, with sensitivity Q10 values ranging from 2.8 to 7.3. The sensitivity of Rs with soil temperature more or less stabilizes with increasing temperature above 30 °C, with Q10 values ranging from 1.28 to 1.39.

Table 5. Estimated Pearson’s Correlation Coefficients of the Relationship between Soil Respiration, Soil Temperature, Air Temperature, and Soil Moisture (n = 240)

Abbreviations: Rs, soil respiration; SoilT, soil temperature; AirT, air temperature; SM, soil moisture

The highest Rs rates were observed in the summer months (i.e., June to August), with soil and air temperatures ranging from 17.47 to 22.81 ℃ and soil moisture from 18.4% to 18.5%. The trend follows a gradual decrease in Rs rates at the start of wet months in September to the end of winter in February. Following the melting of snow, gradual increase in soil and air temperatures, and increased soil moisture in spring (i.e., March through May), the Rs rates gradually increased until they peaked in June. The Rs showed significant positive relationships with soil temperature and air temperature (0.56) but indicated a negative relation with soil moisture (Figs. 3 and 4). These findings suggest that soil respiration increases with an increase in soil and air temperatures, while it decreases as soil moisture increases.

Fig. 3. Soil respiration rates as a function of soil and air temperatures. The correlations show a significant strong positive relationship of Rs rates with soil (0.77) and air temperature (0.75) (n = 240).

Fig. 4. Soil respiration rates as a function of soil moisture. The correlation shows a weak negative relationship of Rs with soil moisture (n = 240).

The calculated values were R2, Mallow’s Cp statistics, Akaike Information Criterion (AIC), and mean standard error (MSE) to determine the variables useful for predicting Rs rates. Variables with the highest R2 but low values of Mallow’s Cp, AIC, and MSE are considered as best candidate variables for the multiple regression model (Pacaldo et al. 2024). Table 6 summarizes the analysis results, which showed that the soil temperature and soil moisture provide a good approximation of Rs, as indicated by a high R2 value but low Cp, AIC, and MSE. Including air temperature in the model did not significantly improve the precision of the model in predicting Rs. Estimated parameter estimates of the regression line of the multiple regression model are summarized in Table 7, which shows that the intercept was not significantly different from zero (p = 0.0898). The soil temperature is shown to be the only variable with a significant probability value (p <0.0001), suggesting that the soil temperature is a significant parameter in predicting Rs rates.

Table 6. Estimated r2 values, Mallows’ Cp Statistics, Akaike Information Criterion (AIC), and Mean Standard Error (MSE)

Abbreviations: SoilT, Soil temperature: AirT, Air temperature; SM, Soil moisture.

Table 7. Estimated Line Intercepts and Constants Values for Independent Variables in the Multiple Regression Model

Changes in Soil Organic Carbon

The ANOVA revealed no significant differences in SOC contents among age classes and interaction effects (Table 8). By contrast, the SOC significantly differed between soil depths (p = 0.0002), with the upper 15-cm soil depth containing significantly lower SOC contents in the upper 10 cm (18.60 to 25.68 Mt C ha-1) than the subsurface layer (31.31 to 45.30 Mt CO2 ha-1) of the mineral soil. The soil organic matter contents on the soil surface (O-horizon) significantly differed in all stand ages, with increasing accumulation of SOM as the forest ages (11.14 to 78.57 Mt SOM ha-1) (Table 9). As expected, the control contained the smallest volume of litter because it receives SOM inputs from grasses only.

Table 8. ANOVA Showing Differences in SOC among Stand Age, Soil Depths (0 to 10 cm and 11 to 20 cm), and the Interaction Effects between Age and Soil Depth (n = 4; p = 0.05)

Table 9. Soil Organic Carbon of the Study Sites’ Different Soil Layers and Forest Litter Deposits

Subscripts in the same letters denote non-significant relationships between forest types of the same soil layer at a 95% confidence level based on Tukey’s test (Mean ± MSE; n = 4).

DISCUSSION

Stand Age Influence on Soil Respiration Rates

The chronosequence approach provides a means to investigate processes that may take decades to develop by utilizing sites of different ages as treatment or basis for describing patterns attributable to individual sites as they age (Yanai et al. 2000). In this study, the Rs increases from young (0- to 10-year-olds) to middle-aged (11- to 20-year-olds) stands and then stabilizes onwards, as indicated by the non-significant differences of Rs rates between middle-aged and pre-mature (35- to 45- year-olds) stands (Table 4). At a young age, Rs increases because of high tree density per square meter, rapid growth, and high fine root production (Litton et al. 2003; Montagnoli et al. 2012; Pregitzer et al. 2000), until it reaches the stability period, which usually occurs at the middle-aged when the canopy fully occupies the available space (Helmisaari and Hallbacken 1999; Vanninen and Makela 1999; Makkonen and Helmisaari 2001; Børja et al. 2008; Tang et al. 2009; Claus and George 2011; Konôpka et al. 2011). When the forest matures, the root and shoot growth ratio stabilizes and reaches an equilibrium in which there is a balance between ecosystem production (inputs) and ecosystem respiration (outputs) (Van Noordwijk and De Willigen 1987; Vogt et al. 1987). The present finding agrees well with some previous studies, reporting Rs increases at a young age (Makkonen and Helmisaari 2001; Sulzman et al. 2005; Sayer et al. 2007; Prévost-Bouré et al. 2010; Zhuang et al. 2023), peaks at middle ages, and stabilizes at the time of canopy closure (Law et al. 2003; Bond- Lamberty et al. 2004; Wiseman and Seiler 2004; Saurette et al. 2006; Tang et al. 2009; Arevalo et al. 2010; Chang et al. 2020).

However, some studies found lower or higher Rs values due to stand density, canopy gaps, root biomass, forest litter production, soil organic matter inputs, and other environmental factors (Irvine and Law 2002; Mayer et al. 2017; Chin et al. 2023). Varik et al. (2015) reported an increasing trend of Rs from young to middle-aged and then stabilized to pre-maturity age due to equilibrium between SOC input and heterotrophic respiration at the maturity period. Similarly, Peichl et al. (2010) observed an increased Rs with stand age in the Pinus strobus forest in Canada due to a steady increase of SOM from aboveground biomass and roots. Other authors also reported a similar pattern in the loblolly pine forest in Virginia, U.S.A (Wiseman et al. 2004), boreal jack pine forest in Canada (Smith et al. 2010), and Scots pine in Estonia (Uri et al. 2022). Payeur-Poirier et al. (2012) observed a decreasing trend of Rs after harvest and then increased with further stand development. In the Scots pine forest in Canada, Uri et al. (2022) reported that Rs showed similar patterns following seasonal fluctuations in soil temperature irrespective of stand age. Other authors reported no change in Rs with stand age in the hybrid poplar plantations in Canada (Chang et al. 2020), Norway Spruce forest in Estonia (Kukumägi et al. 2017), and White pine (Pinus strobus) in Canada (Peichl et al. 2014).

Surprisingly, the control showed significantly higher Rs than the 0 to 10 and 35 to 45 age classes in the summer (Fig. 1), suggesting that treeless areas could contribute a higher CO2 into the atmosphere during warmer periods. The high CO2 emission rates in treeless areas covered by grasses (control) could probably be explained not only by the direct exposure of the soil surface to solar radiation and high temperature but also by the influence of litter quality (Han et al. 2015; Petraglia et al. 2019). Grasses produce a highly decomposable organic matter, hence a rich source of labile carbon, which drives microbial activities and Rs rates (Post and Kwon 2000; Saurette et al. 2006; Teklay and Chang 2008; Petraglia et al. 2019; Chang et al. 2020; Pacaldo et al. 2024). In contrast, the Rs rates under pine forests did not dramatically increase during summer, probably due to the high acid-nonhydrolyzable residue (AUR) or lignin contents of needles, bark, and other residues of pines that inhibit decomposition rates and microbial activities (Prescott 2010; Hasbullah and Marschner 2015).

Environmental Factors Affecting Soil Respiration

In this study, the Rs showed a strong positive correlation with soil temperature (0.77) and air temperatures (0.75), indicating the tendency of Rs to increase with temperatures, which agrees well with some previous reports with the strength of the relationship ranging from moderate to a strong relationship (Dinca et al. 2018; Cui et al. 2020; Pacaldo et al. 2023, 2024; Amarille et al. 2023; Pacaldo et al. 2024). The Rs rates increase as the temperature rises because autotrophic and heterotrophic activities respond positively to high temperatures, particularly at times of abundant soil moisture contents, which could drive decomposition rates (Salah and Scholes 2011; Petraglia et al. 2019). This pattern conforms with our calculations on the sensitivity of Rs to temperature, in which we found that Rs is highly sensitive to increasing temperature with Q10 values ranging from 2.8 to 7.3 at temperatures ranging from 5 °C to 19 °C, consistent with some previous reports (Han et al. 2023).

In contrast, the Rs was negatively correlated with soil moisture, indicating an inversely proportional relationship pattern (Amarille et al. 2023; Pacaldo et al. 2023; Pacaldo et al. 2024). However, findings on the relationship between Rs and soil moisture are not consistent, with some authors reporting a positive relationship (Raich and Schlesinger 1992; Wood et al. 2013; Fei et al. 2015), negative relationship (Adachi et al. 2005; Yanni et al. 2020), and no significant relationship (Borken et al. 2006; Bréchet et al. 2009). Sealing of soil pores, which occurs when soil pores are thoroughly saturated with water, and freezing effects of cold temperature during winter results in suppression of Rs due to the reduced diffusion of CO2 and decreased microbial and root activities (Du et al. 2013; Chang et al. 2014; Pacaldo et al. 2024). In this study, it was also observed that the Rs significantly decreased when the soil moisture dropped to 5.74%, despite the warm temperature (>15 °C), suggesting that, during dry periods, the soil moisture mainly regulates Rs, not the temperature, particularly at times when the soil water becomes the limiting factor of microbial and root activities.

In this study, a combination of soil temperature and moisture factors could strongly predict Rs rates, as indicated by high r-square but low values of Mallow’s Cp, AIC, and MSE, suggesting that these factors influence each other and Rs rates (Dinca et al. 2018; Pacaldo et al. 2024), probably due to their direct influence on root and microbial activities (Subke and Bahn 2010; Chang et al. 2014), ecosystem productivity, and hydrological processes (Kanmani et al. 2023).

Stand Age Effects on Soil Organic Carbon

Contrary to the hypothesis, the statistical analysis failed to detect significant differences in SOC with stand age. The lack of statistically detectable significant changes in the upper 10 cm of the mineral soil suggests that changes in SOC in the pine forest ecosystem occur at a prolonged process, and these changes could not be statistically detected not only due to high spatial variations but also because of slow transformation rate of forest litter into soil organic carbon. Recalcitrant organic matter, such as roots, usually decays very slowly, which can be detected only after a few years (Prescott 2010; Hasbullah and Marschner 2015; Aydin et al. 2018). In contrast, the significant differences in forest litter on the soil surface across different age classes demonstrate that the decomposition rate of organic matter from P. nigra occurs at a very slow pace resulting in its steady accumulation on the forest floor as the forest ages (Table 6). Lignin (AUR), abundant in the forest litter of pines, provides an effective shield against a rapid decomposition process and inhibits microbial activities (Prescott 2010).

Reported chronosequence studies in other forest ecosystem types suggested different trajectories of SOC changes with stand age. A study in the northern hardwood, New Hampshire, predicted a 50 % loss of organic matter in the first 20 years before the disturbed site slowly recovered in the next 50 years and then stabilized onwards. The author attributed the decreased SOM to the rapid loss of organic matter in young stands due to increased decomposition and reduced litter inputs (Covington 1981). Chen et al. (2013) reported a similar pattern of SOC with stand ages in the Chinese fir (Cunninghamia lanceolate), due to management regimes, climate, and edaphic conditions interacting with the SOC. In Pakistan, a chronosequence of Chir Pine (Pinus roxburghii) revealed a decreasing trend of SOC with stand age due to disturbances of forest management operations (Amir et al. 2018, 2019). Chen et al. (2010) reported a similar pattern for the Mongolian pine (Pinus sylvestris var. mongolica Litv.) in China, wherein the SOC decreased from 12 to 40 years due to disturbance caused by wind erosion.

In contrast, Smal et al. (2019) reported an increasing SOC with stand age because litter production exceeded decomposition, gradually increasing SOC in the organic layer. De Marco et al. (2021) observed a similar pattern in Stone pine (Pinus pinea) in Italy, wherein the SOC increases with stand age with about 80 % accumulation in the organic horizon due to high litter production. Other authors also reported a similar pattern of increased SOC with stand age in Chinese pine (Pinus tabulaeformis) (Zhao et al. 2014), lacebark pine (Li et al. 2013), Mongolian pine (Pinus sylvestris var. mongolica) (Lei et al. 2019), Pinus massoniana (Song et al. 2017), and boreal larch forest (Larix gmelinii) (Zhu et al. 2020). Although Yanai et al. (2000) found no pattern in the change in organic matter with time in the Northern Hardwood, New Hampshire, they observed the highest organic matter accumulation in oldest stands and the least in young and middle-aged stands.

Surprisingly, the subsurface layer (11 to 20 cm depth) contained significantly higher SOC stocks compared to the upper 10 cm depth, which is somewhat inconsistent with the authors’ expectations and findings in some previous reports (Smal et al. 2008; Bayramin et al. 2009; Pacaldo 2012; Pacaldo et al. 2013). The higher SOC contents in the subsurface may be associated with the cultural treatment received by the soil before the reforestation project. The plowing and disking of the soil, part of the site preparation, resulted in mixing the soil organic matter in the plow layer (Ap) into the subsurface layer, which remained in place and protected against the attack of decomposers. In contrast, the decomposition of the SOM in the upper 10 cm is expected to occur much faster than the subsurface layer because it is directly exposed to ambient conditions and is the optimum depth of rapid microbial and root activities (Wang et al. 2019; Hao et al. 2021).

CONCLUSIONS

The present findings failed to support the hypothesis. The Rs rates in the P. nigra forest ecosystem significantly increased from the young (0 to 10-year-olds) to middle-aged (11- to 20-year-olds) stands and then stabilized onwards. In contrast, the SOC in the mineral soil showed no statistically significant differences at two 10-cm soil depth increments, with higher SOC contents in the subsurface layer than the upper 10-cm depth. Forest litter on the soil surface significantly increased with stand age and accumulated in large quantities as the forest ages. Although these results may not be conclusive to other types of forest ecosystems, these findings provide additional information regarding the trajectory of Rs and SOC changes in the forest ecosystem with stand age. To our knowledge, this is the first chronosequence study assessing the Rs and SOC changes for the P. nigra forest, a vital piece of information in the assessment of the carbon budget in this type of forest ecosystem. The present findings could also be used in alleviating the uncertainty regarding the critical question of whether to “freeze” the reference baseline over the project’s lifetime in the accounting of carbon sequestration of registered CDM reforestation projects, particularly for P. nigra forest. The results also imply that periodic recalculation of the baseline, e.g., every five years, is needed to increase the reliability and accuracy of the “baseline” estimates, as highlighted by significant increases of Rs rates from young to middle-aged stands and continued accumulation of forest litter as the forest ages. The SOC contents in the mineral soil layers are C neutral, indicating that this C pool could be frozen over the forest’s lifetime, at least in the P. nigra plantations.

Although this study increased our understanding of the stand age effects of SOC and Rs rates in the P. nigra forests, there is no absolute certainty whether the observed trajectory is due to changes of age or other factors because of the high spatial variation of soil organic matter stocks, roots, and microclimate conditions. A chronosequence approach is subject to several potential sources of error due to site variability and other environmental factors affecting Rs and SOC that could be unrelated to the time since the establishment of the forest or the occurrence of disturbance (Yanai et al. 2000). Furthermore, it is also beyond the scope of this study to assess the microclimate effects on Rs and SOC due to differences of canopy sizes and tree height or trees within the borders of different stand age classes. Moreover, climate change impacts and other environmental stresses would also create uncertainties about this observed trajectory of Rs and SOC changes in the P. nigra forest, particularly in the event of forest fires. These subjects are interesting topics for future investigations, particularly in light of changing climate patterns.

ACKNOWLEDGMENTS

The authors are deeply grateful to Prof. Dr. Renato S. Pacaldo of the Mindanao State University-Marawi City and a visiting professor at the Kastamonu for his guidance in preparing the manuscript and conducting this study. Grateful appreciation also goes to Randell Keith Amarille, a graduate student of Kastamonu University, for his invaluable assistance in data collection and analysis, laboratory analysis, and manuscript preparation.

FUNDING SUPPORT

This study was supported by the Turkey Scientific and Technological Research Council (Türkiye Bilimsel ve Teknolojik Araştırma Kurumu) (TÜBİTAK), the Science and Fellowship Grant No. 121C066, under the CoCirculation2 with funding from the European Union’s Horizon 2020 research and innovation program under the Marie Sktodowska-Curie grant agreement No. 801509.

DECLARATION OF COMPETING INTEREST

As part of our ethical obligation as a researcher, we declare that no competing interest is involved in this study, and the research outcome, either financial, commercial, or non-financial conflicts.

REFERENCES CITED

Adachi, M., Bekku, Y. S., Konuma, A., Kadir, W. R., Okuda, T., and Koizumi, H. (2005). “Required sample size for estimating soil respiration rates in large areas of two tropical forests and two types of plantations in Malaysia,” Forest Ecology and Management 210(1-3), 455-459. DOI: 10.1016/j.foreco.2005.02.011

Aksoy, E., Panagos, P., Montanarella, L., and Jones, A. (2010). Integration of the Soil Database of Turkey into European Soil Database 1: 1.000. 000 (European Commission JRC Research Report: EUR 24295EN), European Commission Joint Research Centre, Ispra, Italy.

Amarille, R., Pacaldo, R., and Aydın, M. (2023). “Effects of forest litter reduction on soil respiration rates across a chronosequence of black pine forest,” in: 3rd International Congress on Engineering and Life Science, Trabzon, Turkiye, pp. 114. DOI: 10.61326/icelis.2023.40

Amir, M., Liu, X., Ahmad, A., Saeed, S., Mannan, A., and Muneer, M. A. (2018). “Patterns of biomass and carbon allocation across chronosequence of chir pine (Pinus roxburghii) Forest in Pakistan: Inventory‐based estimate,” Advances in Meteorology 2018(1), article 3095891. DOI: 10.1155/2018/3095891

Arevalo, C. B. M., Bhatti, J. S., Chang, S. X., Jassal, R. S., and Sidders, D. (2010). “Soil respiration in four different land use systems in north-central Alberta, Canada,” Journal of Geophysical Research 115(G1), article ID G01003. DOI: 10.1029/2009JG001006

Aydın, M., Pacaldo, R., and Volk, T. (2018). “Soil respiration in shrub willow (Salix x dasyclados) biomass crop increased in the third year after removal,” International Journal of Global Warming 15, 54-66. DOI: 10.1504/IJGW.2018.091953

Baumert, K. A. (2000). “The clean development mechanism: Understanding additionality,” World Resources Institute, (www.Academia.edu), Accessed 15 March 2024.

Bayramin, I., Basaran, M., Erpul, G., Dolarslan, M., and Canga, M. R. (2009). “Comparison of soil organic carbon content, hydraulic conductivity, and particle size fractions between a grassland and a nearby black pine plantation of 40 years in two surface depths,” Environmental Geology 56, 1563-1575. DOI: 10.1007/s00254-008-1254-8

Berg, A., and Sheffield, J. (2018). “Climate change and drought: The soil moisture perspective,” Current Climate Change Reports 4, 180-191. DOI: 10.1007/s40641-018-0095-0

Bond-Lamberty B., Wang C. K., and Gower S. T. (2004). “Contribution of root respiration to soil surface CO2 flux in a boreal black spruce chronosequence,” Tree Physiology 24(12), 1387-1395. DOI: 10.1093/treephys/24.12.1387

Bond-Lamberty, B., and Thomson, A. (2010). “Temperature-associated increases in the global soil respiration record,” Nature 464, 579-582. DOI: 10.1038/nature08930

Børja, I., de Wit, H. A., Steffenrem, A., and Majdi, H. (2008). “Stand age and fine root biomass, distribution, and morphology in a Norway spruce chronosequence in southeast Norway,” Tree Physiology 28(5), 773-784. DOI: 10.1093/treephys/28.5.773

Borken, W., Savage, K., Davidson, E. A., and Trumbore, S. E. (2006). “Effects of experimental drought on soil respiration and radiocarbon efflux from temperate forest soil,” Global Change Biology 12(2), 177-193. DOI: 10.1111/j.1365-2486.2005.001058.x

Bouyoucos, G. J. (1962). Hydrometer method improved for making particle size analyses of soils. 1,” Agronomy Journal 54(5), 464. DOI: 10.2134/agronj1962.000219.

Bréchet, L., Ponton, S., Roy, J., Freycon, V., Couteaux, M. M., Bonal, D., Epron, D. (2009). “Do tree species characteristics influence soil respiration in tropical forests? A test based on 16 tree species planted in monospecific plots,” Plant and Soil 319, 235-246. DOI: 10.1007/s11104-008-9866-z

Cao, J., Wang, X., Tian, Y., Wen, Z., and Zha, T. (2012). “Pattern of carbon allocation across three different stages of stand development of a Chinese pine (Pinus tabulaeformis) forest,” Ecological Research 27, 883-892. DOI: 10.1007/s11284-012-0965-1

Chang, C. T., Sabate, S., Sperlich, D., Poblador, S., Sabater, F., and Gracia, C. (2014). “Does soil moisture overrule temperature dependence on soil respiration in Mediterranean riparian forests?,” Biogeosciences 11(21), 6173-6185. DOI: 10.5194/bg-11-6173-2014

Chang, S. X., Shi, Z., and Thomas, B. R. (2020). “Soil respiration and net ecosystem productivity in a chronosequence of hybrid poplar plantations,” Canadian Journal of Soil Science 100(4), 488-502. DOI: 10.1139/cjss-2020-0006

Chen, A., Wang, Z., Lin, Y., Wang, X., Li, Y., Zhang, Y., Zhang Y., Tao Z., Gao Q., and Tang, G. (2020). “Temporal variation of soil organic carbon pools along a chronosequence of reforested land in Southwest China,” Catena 194, article 104650. DOI: 10.1016/j.catena.2020.104650

Chen, F. S., Zeng, D. H., Fahey, T. J., and Liao, P. F. (2010). “Organic carbon in soil physical fractions under different-aged plantations of Mongolian pine in semi-arid region of Northeast China,” Applied Soil Ecology 44(1), 42-48. DOI: 10.1016/j.apsoil.2009.09.003

Chen, G. S., Yang, Z. J., Gao, R., Xie, J. S., Guo, J. F., Huang, Z. Q., and Yang, Y. S. (2013). “Carbon storage in a chronosequence of Chinese fir plantations in southern China,” Forest Ecology and Management 300, 68-76. DOI: 10.1016/j.foreco.2012.07.046

Cheng, X., Han, H., Kang, F., Song, Y., and Liu, K. (2014). “Variation in biomass and carbon storage by stand age in pine (Pinus tabulaeformis) planted ecosystem in Mt. Taiyue, Shanxi, China,” Journal of Plant Interactions 9(1), 521-528. DOI: 10.1080/17429145.2013.862360

Chin, M. Y., Lau, S. Y. L., Midot, F., Jee, M. S., Lo, M. L., Sangok, F. E., and Melling, L. (2023). “Root exclusion methods for partitioning of soil respiration: Review and methodological considerations,” Pedosphere 33(5), 683-699. DOI: 10.1016/j.pedsph.2023.01.015

Claus, A., and George, E. (2011). “Effect of stand age on fine-root biomass and biomass distribution in three European forest chronosequences,” Canadian Journal of Forest Research 35, 1617-1625. DOI: 10.1139/x05-079

Covington, W. W. (1981). “Changes in forest floor organic matter and nutrient content following clear cutting in northern hardwoods,” Ecology 62(1), 41-48.

Cui, Y. B., Feng, J. G., Liao, L. G., Yu, R., Zhang, X., Liu, Y. H., Yang, L. Y., Zhao, J. F., and Tan, Z. H. (2020). “Controls of temporal variations on soil respiration in a tropical lowland rainforest in Hainan Island, China,” Tropical Conservation Science 13, 1-14. DOI: 10.1177/1940082920914902

Darenova, E., Fabiánek, T., and Pavelka, M. (2016). “Efflux of CO2 from soil in Norway spruce stands of different ages: A case study,” European Journal of Environmental Sciences 6(2). DOI: 10.14712/23361964.2016.14

De Marco, A., Berg, B., Zarrelli, A., and De Santo, A. V. (2021). “Shifts in soil chemical and microbial properties across forest chronosequence on recent volcanic deposits,” Applied Soil Ecology 161, article 103880. DOI: 10.1016/j.apsoil.2021.103880

Dinca, L., Badea, O., Guiman, G., Braga, C., Crisan, V., Greavu, V., Murariu, G., and Georgescu, L. (2018). “Monitoring of soil moisture in Long-Term Ecological Research (LTER) sites of Romanian Carpathians,” Annals of Forest Research 61(2), 171-188. DOI: 10.1016/j.soilbio.2011.08.012

Du, E., Zhou, Z., Li, P., Jiang, L., Hu, X., and Fang, F. (2013). “Winter soil respiration during the soil-freezing process in a boreal forest in Northeast China,” Journal of Plant Ecology 6(5), 349-357. DOI: 10.1093/jpe/rtt012

Fei, P., Manhou, X., Quangang, Y., Xuhui, Z., Tao, W., and Xian, X. (2015). “Different responses of soil respiration and its components to experimental warming with contrasting soil water content,” Arctic, Antarctic, and Alpine Research 47(2), 359-368. DOI: 10.1657/AAAR0014-018

Francis Justine, M., Yang, W., Wu, F., Tan, B., Naeem Khan, M., and Zhao, Y. (2015). “Biomass stock and carbon sequestration in a chronosequence of Pinus massoniana plantations in the upper reaches of the Yangtze River,” Forests 6(10), 3665-3682. DOI: 10.3390/f6103665

Gangwar, D. P., and Baskar, M. (2019). Texture Determination of Soil by Hydrometer Method for Forensic Purpose, Central Forensic Science Laboratory, Chandigarh, India.

Groffman, P. M., Driscoll, C. T., Fahey, T. J., Hardy, J. P., Fitzhugh, R. D., and Tierney, G. L. (2001). “Colder soils in a warmer world: a snow manipulation study in a northern hardwood forest ecosystem,” Biogeochemistry 56, 135-150. DOI: 10.1023/A:1013039830323

Han, T., Huang, W., Liu, J., Zhou, G., and Xiao, Y. (2015). “Different soil respiration responses to litter manipulation in three subtropical successional forests,” Scientific Reports 5(1), article 18166. DOI: 10.1038/srep18166

Han, Y., Wang, G., Zhou, S., Li, W., and Xiong, L. (2023). “Day–night discrepancy in soil respiration varies with seasons in a temperate forest,” Functional Ecology 37(7), 2002-2013. DOI: 10.1111/1365-2435.14358

Hao, J., Chai, Y. N., Lopes, L. D., Ordonez, R., Wright, E., Archontoulis, S., and Schachtman, D. (2021). “The effects of soil depth on the structure of microbial communities in agricultural soils in Iowa (United States),” Applied and Environmental Microbiology 87(4), article ID e02673-20. DOI: 10.1128/AEM.02673-20

Hasbullah, H., and Marschner, P. (2015). “Residue properties influence the impact of salinity on soil respiration,” Biology and Fertility of Soils 51, 99-111. DOI: 10.1007/s00374-014-0955-2

Hashimoto, S., Carvalhais, N., Ito, A., Migliavacca, M., Nishina, K., and Reichstein, M. (2015). “Global spatiotemporal distribution of soil respiration modeled using a global database,” Biogeosciences 12, 4331-4364. DOI: 10.5194/bg-12-4121-2015

Helmisaari, H. S., and Hallbäcken, L. (1999). “Fine-root biomass and necromass in limed and fertilized Norway spruce (Picea abies (L.) Karst.) stands,” Forest Ecology and Management 119(1-3), 99-110. DOI: 10.1016/S0378-1127(98)00514-3

Huang, Z., Cui, Z., Liu, Y., and Wu, G. L. (2021). “Carbon accumulation by Pinus sylvestris forest plantations after different periods of afforestation in a semiarid sandy ecosystem,” Land Degradation & Development 32(6), 2094-2104. DOI: 10.1002/ldr.3858

IPCC (2001). Climate Change 2001: The Scientific Basis, Cambridge University Press, Cambridge, UK.

IPCC (2021). Climate Change 2021: The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change, Cambridge University Press, Cambridge, UK.

IPCC. (2006) Forest lands. Intergovernmental Panel on Climate Change Guidelines for National Greenhouse Gas Inventories; Institute for Global Environmental Strategies (IGES): Hayama, Japan; Volume 4, p. 83

Irvine, J., and Law, B. (2002). “Contrasting soil respiration in young and old-growth ponderosa pine forests,” Global Change Biology 8, 1183-1194. DOI: 10.1046/j.1365-2486.2002.00544.x

Kanmani, K., Vasanthi, P., Pari, P., and Shafeer Ahamed, N. S. (2023). “Estimation of soil moisture for different crops using SAR polarimetric data,” Civ. Eng. J. 9(6), 1402-11. DOI: 10.28991/CEJ-2023-09-06-08

Konôpka, B., Pajtík, J., Moravčík, M., and Lukac, M. (2011). “Biomass partitioning and growth efficiency in four naturally regenerated forest tree species,” Basic and Applied Ecology 2(3), 234-243. DOI: 10.1016/j.baae.2010.02.004

Kukumägi, M., Ostonen, I., Uri, V., Helmisaari, H. S., Kanal, A., Kull, O., and Lõhmus, K. (2017). “Variation of soil respiration and its components in hemiboreal Norway spruce stands of different ages,” Plant and Soil 414, 265-280. DOI: 10.1007/s11104-016-3133-5

Law, B. E., Sun, O. J., Campbell, J., Van Tuyl, S., and Thornton, P. E. (2003). “Changes in carbon storage and fluxes in a chronosequence of ponderosa pine,” Global Change Biology 9, 510-524. DOI: 10.1046/j.1365-2486.2003.00624.x

Lei, Z., Yu, D., Zhou, F., Zhang, Y., Yu, D., Zhou, Y., and Han, Y. (2019). “Changes in soil organic carbon and its influencing factors in the growth of Pinus sylvestris var. mongolica plantation in Horqin Sandy Land, Northeast China,” Scientific Reports 9(1), article 16453. DOI: 10.1038/s41598-019-52945-5

Li, C., Zha, T., Liu, J., and Jia, X. (2013). “Carbon and nitrogen distribution across a chronosequence of secondary lacebark pine in China,” The Forestry Chronicle 89(2), 192-198. DOI: 10.5558/tfc2013-037

Liang N, Nakadai T, Hirano T, Qu L, Koike T, Fujinuma Y, and Inoue G. (2004). “In situ comparison of four approaches to estimating soil CO2 efflux in a northern larch (Larix kaempferi Sarg.) forest,” Agricultural and Forest Meteorology 123, 97-117. DOI: 10.1016/j.agrformet.2003.10.002

Liptzin, D., Williams, M. W., Helmig, D., Seok, B., Filippa, G., Chowanski, K., and Hueber, J. (2009). “Process-level controls on CO2 fluxes from a seasonally snow-covered subalpine meadow soil, Niwot Ridge, Colorado,” Biogeochemistry 95, 151-166. DOI: 10.1007/s10533-009-9303-2

Litton, C. M., Ryan, M. G., Tinker, D. B. and Knight, D. H. (2003). “Belowground and aboveground biomass in young postfire lodge pole pine forests of contrasting tree density,” Canadian Journal of Forest Research 33, 351-353. DOI: 10.1139/x02-181

Lloyd, J., and Taylor, J. A. (1994). “On the temperature dependence of soil respiration,” Functional Ecology 8(3), 315-323. DOI: 10.2307/2389824

Luan, J., Liu, S., Zhuf, X., Wang, J., Liu, K. (2012). “Roles of biotic and abiotic variables in determining the spatial variation of soil respiration in secondary oak and planted pine forests,” Soil Biology and Biochemistry 44, 143-150. DOI: 10.1016/j.soilbio.2011.08.012

Makkonen, K., and Helmisaari, H. S. (2001). “Fine root biomass and production in Scots pine stands in relation to stand age,” Tree Physiology 21(2-3), 193-198. DOI: 10.1093/treephys/21.2-3.193

Mayer, M., Matthews, B., Rosinger, C., Sandén, H., Godbold, D. L., and Katzensteiner, K. (2017). “Tree regeneration retards decomposition in temperate mountain soil after forest gap disturbance,” Soil Biology and Biochemistry 115, 490–498. DOI: 10.1016/j.soilbio.2017.09.010

Montagnoli, A., Terzaghi, M., di Iorio, A., Scippa, G., and Chiatante, D. (2012). “Fine-root seasonal pattern, production, and turnover rate of European beech (Fagus sylvatica L.) stands in Italy Prealps: Possible implications of coppice conversion to high forest,” Plant Biosystems 146, 1012-1022. DOI: 10.1080/11263504.2012.741626

Nelson, D. A., and Sommers, L. (1983). “Total carbon, organic carbon, and organic matter,” in: Methods of Soil Analysis: Part 2 Chemical and Microbiological Properties 9, 539-579. DOI: 10.2134/agronmonogr9.2.2ed.c29

Nissan, A., Alcolombri, U., Peleg, N., Galili, N., Jimenez-Martinez, J., Molnar, P., and Holzner, M. (2023). “Global warming accelerates soil heterotrophic respiration,” Nature Communications 14(1), article 3452. DOI: 10.1038/s41467-023-38981-w

Noh, N. J., Son, Y., Lee, S. K., Seo, K. W., Heo, S. J., Yi, M. J., … and Lee, K. H. (2010). “Carbon and nitrogen storage in an age-sequence of Pinus densiflora stands in Korea,” Science China Life Sciences 53, 822-830. DOI: 10.1007/s11427-010-4018-0

Özden, D. M., Keskin, S., Dinç, U., Kapur, S., Akça, E., Şenol, S., and Dinç, O. (2001). “Soil geographical database of Turkey at a scale of 1: 1.000. 000,” World Soil Survey Archive and Catalogue (WOSSAC), (http://www.wossac.com/downloads/19811_turkeysoilmap.pdf), Accessed 06 Feb 2024.

Pacaldo, R. S. (2012). Carbon Balances in Shrub Willow Biomass Crops Along A 19-year Chronosequence as Affected by Continuous Production and Crop Removal (Tear-Out) Treatments, SUNY College of Environmental Science and Forestry, Syracuse, NY, USA.

Pacaldo, R. S., Aydin, M., and Amarille, R. K. (2024). “Soil respiration and controls in warmer winter: A snow manipulation study in postfire and undisturbed black pine forests,” Ecology and Evolution 14(3), article e11075. DOI: 10.1002/ece3.11075

Pacaldo, R. S., Volk, T. A., and Briggs, R. D. (2013). “No significant differences in soil organic carbon contents along a chronosequence of shrub willow biomass crop fields,” Biomass and Bioenergy 58, 136-142. DOI: 10.1016/j.biombioe.2013.10.018

Pacaldo, R. S., Volk, T. A., Briggs, R. D., Abrahamson, L. P., Bevilacqua, E., and Fabio, E. S. (2014). “Soil CO2 effluxes, temporal and spatial variations, and root respiration in shrub willow biomass crop fields along a 19‐year chronosequence as affected by regrowth and removal treatments,” GCB Bioenergy 6(5), 488-498. DOI: 10.1111/gcbb.12108

Pacaldo, R., and Aydın, M. (2023). “Soil respiration in a natural forest and a plantation during a dry period in the Philippines,” Journal of Forestry Research 34, 1975-1983. DOI: 10.1007/s11676-023-01636-z

Pang, X., Bao, W., Zhu, B., and Cheng, W. (2013). “Responses of soil respiration and its temperature sensitivity to thinning in a pine plantation,” Agricultural and Forest Meteorology 171–172, 57-64. DOI: 10.1016/J.AGRFORMET.2012.12.001

Payeur-Poirier, J. L., Coursolle, C., Margolis, H. A., and Giasson, M. A. (2012). “CO2 fluxes of a boreal black spruce chronosequence in eastern North America,” Agricultural and Forest Meteorology 153, 94-105. DOI: 10.1016/j.agrformet.2011.07.009

Pearson, T., Walker, S., and Brown, S. (2013). “Sourcebook for land use, land-use change, and forestry projects,” Winrock International, (https://winrock.org/wp-content/uploads/2016/03/Winrock-BioCarbon_Fund_Sourcebook-compressed.pdf), Accessed 03 April 2024.

Peichl, M., Arain, A. M., Moore, T. R., Brodeur, J. J., Khomik, M., Ullah, S., … and Pejam, M. R. (2014). “Carbon and greenhouse gas balances in an age sequence of temperate pine plantations,” Biogeosciences 11(19), 5399-5410. DOI: 10.5194/bg-11-5399-2014

Peichl, M., Arain, M. A., and Brodeur, J. J. (2010). “Age effects on carbon fluxes in temperate pine forests,” Agricultural and Forest Meteorology 150(7-8), 1090-1101. DOI: 10.1016/j.agrformet.2010.04.008

Petraglia, A., Cacciatori, C., Chelli, S., Fenu, G., Calderisi, G., Gargano, D., Abeli, T., Orsenigo, S., and Carbognani, M. (2019). “Litter decomposition: Effects of temperature driven by soil moisture and vegetation type,” Plant and Soil 435, 187-200. DOI: 10.1007/s11104-018-3889-x

Post, W. M., and Kwon, K. C. (2000). “Soil carbon sequestration and land‐use change: Processes and potential,” Global Change Biology 6(3), 317-327. DOI: 10.1046/j.1365-2486.2000.00308.x

Post, W. M., Emanuel, W. R., Zinke, P. J., and Stangenberger, A. G. (1982). “Soil carbon pools and world life zones,” Nature 298, 156-159. DOI: 10.1038/298156a0

Pregitzer, K., Zak, D., Maziasz, J., Deforest, J., Curtis, P., and Lussenhop, J. (2000). “Interactive effects of atmospheric CO2 and soil‐N availability on fine roots of Populus tremuloides,” Ecological Applications 10, 18-33. DOI: 10.1890/1051-0761(2000)010[0018:IEOACA]2.0.CO;2

Prescott, C. E. (2010). “Litter decomposition: What controls it, and how can we alter it to sequester more carbon in forest soils?,” Biogeochemistry 101, 133-149. DOI: 10.1007/s10533-010-9439-0

Prévost-Bouré, N. C., Soudani, K., Damesin, C., Berveiller, D., Lata, J. C., and Dufrêne, E. (2010). “An increase in aboveground fresh litter quantity over-stimulates soil respiration in a temperate deciduous forest,” Applied Soil Ecology 46(1), 26-34. DOI: 10.1016/j.apsoil.2010.06.004

Raich, J. W., and Schlesinger, W. H. (1992). “The global carbon dioxide flux in soil respiration and its relationship to vegetation and climate,” Tellus B 44(2), 81-99. DOI: 10.1034/j.1600-0889.1992.t01-1-00001.x

Saint André, L., Roupsard, O., Marsden, C., Thongo M’Bou, A., D’Annunzio, R., De Grandcourt, A., Jourdan C., Derrien D., Picard N., Zeller B., …. and Laclau, J. P. (2007). “Literature review on current methodologies to assess C balance in CDM afforestation/reforestation projects and a few relevant alternatives for assessing water and nutrient balance, as a complement to carbon sequestration assessments,” Project no 037132. CARBOAFRICA. Quantification, Understanding and Prediction of Carbon Cycle, and other GHG Gases, in Sub-Saharan Africa. Work page N° 6.

Sakici, O. E., Seki, M., and Saglam, F. (2018). “Above-ground biomass and carbon stock equations for Crimean pine stands in Kastamonu region of Turkey,” Fresenius Environmental Bulletin 27(10), 7079-7089.

Salah, Y. M. S., and Scholes, M. C. (2011). “Effect of temperature and litter quality on the decomposition rate of Pinus patula needle litter,” Procedia Environmental Sciences 6, 180-193. DOI: 10.1016/j.proenv.2011.05.019

Saurette, D. D., Chang, S. X., and Thomas, B. R. (2006). “Some characteristics of soil respiration in hybrid poplar plantations in northern Alberta,” Canadian Journal of Soil Science 86, 257-268. DOI: 10.4141/S05-083

Sayer, E. J., Powers, J. S., and Tanner, E. V. J. (2007). “Increased litterfall in tropical forests boosts the transfer of soil CO2 to the atmosphere,” PLoS One 2(12), article e1299. DOI: 10.1371/journal.pone.0001299

Schimel, D., House, J., Hibbard, K., Bousquet, P., Ciais, P., Peylin, P., Braswell, B., Apps, M. J., Baker, D., Bondeau, A., et al. (2001). “Recent patterns and mechanisms of C exchange by terrestrial ecosystems,” Nature 414, 169-172. DOI: 10.1038/35102500

Smal, H., Ligęza, S., Pranagal, J., Urban, D., and Pietruczyk-Popławska, D. (2019). “Changes in the stocks of soil organic carbon, total nitrogen and phosphorus following afforestation of post-arable soils: A chronosequence study,” Forest Ecology and Management 451, and 117536. DOI: 10.1016/j.foreco.2019.117536

Smith, D. R., Kaduk, J. D., Balzter, H., Wooster, M. J., Mottram, G. N., Hartley, G., … and Stocks, B. J. (2010). “Soil surface CO2 flux increases with successional time in a fire scar chronosequence of Canadian boreal jack pine forest,” Biogeosciences 7(5), 1375-1381. DOI: 10.5194/bg-7-1375-2010

Song, X., Kimberley, M. O., Zhou, G., and Wang, H. (2017). “Soil carbon dynamics in successional and plantation forests in subtropical China,” Journal of Soils and Sediments 17, 2250-2256. DOI: 10.1007/s11368-016-1421-6

Sonmez, S., Buyuktas, D., Okturen Asri, F., and Citak, S. (2008). “Assessment of different soil-to-water ratios (1:1, 1:2.5, 1:5) in soil salinity studies,” Geoderma 144, 361-369. DOI: 10.1016/j.geoderma.2007.12.005

Subke, J.-A., and Bahn, M. (2010). “On the ‘temperature sensitivity’ of soil respiration: Can we use the immeasurable to predict the unknown?” Soil Biology and Biochemistry 42, 1653-1656. DOI: 10.1016/j.soilbio.2010.05.026

Sulzman, E. W., Brant, J. B., Bowden, R. D., and Lajtha, K. (2005). “Contribution of aboveground litter, belowground litter and rhizosphere respiration to total soil CO2 efflux in an old growth coniferous forest,” Biogeochemistry 73, 231-256. DOI: 10.1007/s10533-004-7314-6

Tang, J., Bolstad, P., and Martin, J. (2009). “Soil carbon fluxes and stocks in a Great Lakes forest chronosequence,” Global Change Biology 15(1), 145-155. DOI: 10.1111/j.1365-2486.2008.01741.x

Teklay, T., and Chang, S. X. (2008). “Temporal changes in soil carbon and nitrogen storage in a hybrid poplar chronosequence in northern Alberta,” Geoderma 144, 613-619. DOI: 10.1016/j.geoderma.2008.01.023

Tong, X., Xiao, J., Liu, P., Zhang, J., Zhang, J., Yu, P., Meng, P., and Li, J. (2023). “Carbon exchange of forest plantations: global patterns and biophysical drivers,” Agricultural and Forest Meteorology 336, article ID 109379. DOI: 10.1016/j.agrformet.2023.109379

Trenberth, K. E. (2005). “The impact of climate change and variability on heavy precipitation, floods, and droughts,” Encyclopedia of Hydrological Sciences 17, 1-11.

Turkes, M. (2020). “Climate and drought in Turkey,” Water Resources of Turkey, 85-125. DOI: 10.1007/978-3-030-11729-0_4.

Turkish State Meteorological Service. (2024). “State of the Climate of Turkey in 2023”. www.mgm.gov.tr/eng/Yearly-climate/State_of_the_Climate_in_Turkey_in_2023.pdf

Uri, V., Kukumägi, M., Aosaar, J., Varik, M., Becker, H., Aun, K., … and Padari, A. (2022). “The dynamics of the carbon storage and fluxes in Scots pine (Pinus sylvestris) chronosequence,” Science of the Total Environment 817, article 152973. DOI: 10.1016/j.scitotenv.2022.152973

Van Noordwijk, M., and De Willigen, P. (1987). “Agricultural concepts of roots: From morphogenetic to functional equilibrium between root and shoot growth,” Netherlands Journal of Agricultural Science 35(4), 487-496. DOI: 10.18174/njas.v35i4.16707

Vanninen, P., and Mäkelä, A. (1999). “Fine root biomass of Scots pine stands differing in age and soil fertility in southern Finland,” Tree Physiology 19, 823-830. DOI: 10.1093/treephys/19.12.823

Varik, M., Kukumägi, M., Aosaar, J., Becker, H., Ostonen, I., Lõhmus, K., and Uri, V. (2015). “Carbon budgets in fertile silver birch (Betula pendula Roth) chronosequence stands,” Ecological Engineering 77, 284-296. DOI: 10.1016/j.ecoleng.2015.01.041

Vogt, K. A., Vogt, D. J., Moore, E. E., Fatuga, B. A., Redlin, M. R., and Edmonds, R. L. (1987). “Conifer and angiosperm fine-root biomass in relation to stand age and site productivity in Douglas-fir forests,” The Journal of Ecology 75(3), 857-870. DOI: 10.2307/2260210

Wang, J., Liu, H., Hu, M., Du, Y., Liu, Y., Lu, L., and Han S. (2023). “Effects of decreased precipitation and thinning on soil respiration in a temperate forest: A one-year field experiment in Central China,” Catena 229, article ID 107239. DOI: 10.1016/j.catena.2023.107239

Wang, S., Sun, L., Ling, N., Zhu, C., Chi, F., Li, W., Hao, X., Zhang, W., Bian, J., Chen, L., et al. (2019). “Exploring soil factors determining composition and structure of the bacterial communities in saline-alkali soil of Songnen plain,” Frontiers in Microbiology 10, article 2902. DOI: 10.3389/fmicb.2019.02902

Wei, Z., Lin, C., Xu, C., Xiong, D., Liu, X., Chen, S., Lin, T., Yang, Z., and Yang, Y. (2022). “Soil respiration in planted and naturally regenerated Castanopsis carlesii forests during three years post-establishment,” Forests 13(6), article 931. DOI: 10.3390/f13060931

Wiseman, P. E., and Seiler, J. R. (2004). “Soil CO2 efflux across four age classes of plantation loblolly pine (Pinus taeda L.) on the Virginia Piedmont,” Forest Ecology and Management 192(2-3), 297-311. DOI: 10.1016/j.foreco.2004.01.017

Wiseman, P. E., and Seiler, J. R. (2004). “Soil CO2 efflux across four age classes of plantation loblolly pine (Pinus taeda L.) on the Virginia Piedmont,” Forest Ecology and Management 192(2-3), 297-311. DOI: 10.1016/j.foreco.2004.01.017

Wood, T. E., Detto, M., and Silver, W. L. (2013). “Sensitivity of soil respiration to variability in soil moisture and temperature in a humid tropical forest,” PLoS ONE 8(12), article e80965. DOI: 10.1371/journal.pone.0080965

Xu, L., Furtaw, M. D., Madsen, R. A., Garcia, R. L., Anderson, D. J., and McDermitt, D. K. (2006). “On maintaining pressure equilibrium between a soil CO2 flux chamber and the ambient air,” Journal of Geophysical Research: Atmospheres 111(D8). DOI: 10.1029/2005JD006435

Xu, M., Zhang, T., Zhang, Y., Chen, N., Zhu, J., He, Y., Zhao, T., and Yu, G. (2021). “Drought limits alpine meadow productivity in northern Tibet,” Agricultural and Forest Meteorology 303, article ID 108371. DOI: 10.1016/j.agrformet.2021.108371

Yanai, R. D., Arthur, M. A., Siccama, T. G., and Federer, C. A. (2000). “Challenges of measuring forest floor organic matter dynamics: Repeated measures from a chronosequence,” Forest Ecology and Management 138(1-3), 273-283. DOI: 10.1016/S0378-1127(00)00402-3

Yanni, S. F., Helgason, B. L., Janzen, H. H., Ellert, B. H., and Gregorich, E. G. (2020). “Warming effects on carbon dynamics and microbial communities in soils of diverse texture,” Soil Biology and Biochemistry 140, article ID 107631. DOI: 10.1016/j.soilbio.2019.107631

Yu, K., Yao, X., Deng, Y., Lai, Z., Lin, L., and Liu, J. (2019). “Effects of stand age on soil respiration in Pinus massoniana plantations in the hilly red soil region of Southern China,” Catena 178, 313-321. DOI: 10.1016/j.catena.2019.03.038

Yu, S., Wang, D., Dai, W., and Li, P. (2014). “Soil carbon budget in different-aged Chinese fir plantations in south China,” Journal of Forestry Research 25(3), 621-626. DOI: 10.1007/s11676-014-0500-z

Zhao, J., Kang, F., Wang, L., Yu, X., Zhao, W., Song, X., Zhang Y., Chen F., Sun Y., He T., and Han, H. (2014). “Patterns of biomass and carbon distribution across a chronosequence of Chinese pine (Pinus tabulaeformis) forests,” PLoS One 9(4), article e94966. DOI: 10.1371/journal.pone.0094966

Zhu, J., Wang, C., Zhou, Z., Zhou, G., Hu, X., Jiang, L., Li Y., Liu G., Ji C., .,…. and Fang, J. (2020). “Increasing soil carbon stocks in eight permanent forest plots in China,” Biogeosciences 17(3), 715-726. DOI: 10.5194/bg-17-715-2020

Zhuang, W., Liu, M., Wu, Y., Ma, J., Zhang, Y., Su, L., Liu, Y., Zhao, C., and Fu, S. (2023). “Litter inputs exert greater influence over soil respiration and its temperature sensitivity than roots in a coniferous forest in a north-south transition zone,” Science of The Total Environment 886, article ID 164009. DOI: 10.1016/j.scitotenv.2023.164009

Zscheischler, J., Mahecha, M. D., Von Buttlar, J., Harmeling, S., Jung, M., Rammig, A., Randerson, J. T., Schölkopf, B., Seneviratne, S., Tomelleri, E., et al. (2014). “A few extreme events dominate global interannual variability in gross primary production,” Environmental Research Letters 9(3), article ID 035001. DOI: 10.1088/1748-9326/9/3/035001

Article submitted: April 19, 2024; Peer review completed: June 16, 2024; Revised version received: July 7, 2024; Accepted: July 8, 2024; Published: July 18, 2024.

DOI: 10.15376/biores.19.3.6095-6119