Next Article in Journal
Sight versus Sound: Do Visual Assessments of Dead Standing Trees Reflect Acoustic Nondestructive Evaluations of Wood Quality?
Previous Article in Journal
Distribution Pattern and Factors Influencing Spontaneous Plant Diversity in Different Wetland Habitats
Previous Article in Special Issue
The Unabated Atmospheric Carbon Losses in a Drowning Wetland Forest of North Carolina: A Point of No Return?
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Impact of Thinning and Clear Cut on the Ecosystem Carbon Storage of Scots Pine Stands under Maritime Influence in Flanders, Belgium

1
PLECO (Plants and Ecosystems), Department of Biology, University of Antwerp, 2610 Antwerp, Belgium
2
Computational Applied Vegetation Ecology, Department of Environment, Faculty of Bioscience Engineering, Ghent University, 9000 Gent, Belgium
3
Isotope Bioscience Laboratory—ISOFYS, Department of Green Chemistry and Technology, Faculty of Bioscience Engineering, Ghent University, 9000 Ghent, Belgium
4
Forest and Nature Lab, Department Environment, Faculty of Bioscience Engineering, Ghent University, 9090 Melle, Belgium
*
Author to whom correspondence should be addressed.
Forests 2022, 13(10), 1679; https://doi.org/10.3390/f13101679
Submission received: 31 August 2022 / Revised: 26 September 2022 / Accepted: 4 October 2022 / Published: 12 October 2022

Abstract

:
A shift in management to improve the ecological function of mature plantations of exotic species can have important effects on the ecosystem climate mitigation potential. This study investigated the effect of two common forest management strategies for Scots pine (Pinus sylvestris L.) stands on the C storage after 15 years of management. Two pairs of forest stands on poor sandy soil and under the maritime influence in Brasschaat, Belgium, were observed as case studies. The observed forest management strategies were (i) thinning and group planting of oak saplings (Quercus robur L.) and (ii) clear cut, followed by replanting of young oak. For each stand, all forest C pools (aboveground biomass, belowground biomass, litter, and mineral soil) were determined. Results showed, surprisingly, no significant difference in the whole ecosystem C stock for both forest management strategies after 15 years of management. However, after the clear cut and the new plantation, the C in the top 30 cm layer of the mineral soil increased, while it decreased on the forest floor. For thinning with group planting, the C stocks reduced within the 10–30 cm soil layer without impact on the total soil C. Therefore, the shift in management did result in a different allocation of the belowground C, particularly after a clear cut. The results are not only relevant for the study region but also for managed Scots pine forests in neighboring regions of the Atlantic zone of Western Europe.

1. Introduction

Forests store about 50% of the terrestrial carbon (C), making them a very important component within the global C cycle [1,2]. Annually, about 40–60% of the C assimilated by forests through photosynthesis is allocated to biomass [3,4,5]. Part of this biomass dies, falls as litter, and is then partly broken down by decomposers, which return C to the atmosphere (heterotrophic respiration). The remaining part of the litter is incorporated in a recalcitrant litter C pool and soil organic matter. A fraction of the C in the soil organic matter (in particular the microbial products) also enters the mineral soil. The recalcitrant litter, the soil organic matter and the mineral soil store the majority of the belowground C of forest ecosystems [3,5,6,7].
The amount of C stored in a forest depends on the environmental (e.g., climate) and site conditions (e.g., soil hydrology). This explains the high C stocks in tropical forests (climate promoting forest productivity) and boreal forests (climate and hydrology delaying decomposition) [8]. However, temperate forests also store an important proportion of the C globally present in forest ecosystems, i.e., ca. 10–14% [9,10]. In temperate forests, the C storage crucially depends on factors related to the forest management, such as the tree species planted, the amount of wood and other materials that is removed from the site or the age of the forest (rotation time) [1,5,11].
For a long time, forest management in the Atlantic and Central European temperate zone focused mainly on the economic function of commercial wood production [12]. In Flanders, Belgium, monocultures of exotic evergreen species (e.g., Pinus sylvestris L., Pinus nigra subsp. laricio (Poiret) Maire, Picea abies (L.) H. Karst.) often dominated the forest cover because of their high production and adaptability to suboptimal site conditions [5,12,13,14,15]. The classical management of these types of stands was based on a short rotation time (ca. 40–50 years) and a final clear cut [12]. In the past three decades however, other forest functions, such as biodiversity conservation, recreation, air quality control, and climate mitigation, received more attention, resulting in a change of forest management policy [12]. The conversion of evergreen monocultures to mixed or deciduous stands of native species (e.g., pedunculate oak) became one of the applied interventions aiming to improve the forest services. For Flanders, the forest inventory of the Agency for Nature and Forest [16] indicated indeed an increase in the share of hardwood (+3%) and mixed stands (+6%) between 1997–1999 and 2009–2019 at the expense of coniferous stands (−10%). However, we still do not know the impact of interventions improving the ecological function of forests on other ecosystem services such as climate mitigation. Studies focusing on the indirect effects of forest management shifts on the ecosystem C cycles are therefore timely and needed.
This study examines the effect of two typical management strategies of Scots pine (Pinus sylvestris L.) stands on the ecosystem C stocks. The first management strategy combined classical thinning regimes (every 5–8 years) with planting of saplings of pedunculate oak (Quercus robur L.) in understory gaps (group planting). A gradual transition was therefore established by harvesting pine for wood production and replacing it with a native deciduous species. The light and space created by the thinning allow the understory to develop, giving the planted saplings, but also other trees that grow spontaneously, more opportunities to grow and thereby becoming seed trees for the newly developed broadleaved stand in the long term [17,18,19,20]. According to previous studies [21,22], an increasing intensity of thinning results, as expected, in an increasing loss of C in the aboveground biomass, but it has generally no significant impact on the C storage in the mineral soil while it can have a negative impact on the forest floor C. However, in the literature, the reported effect of thinning on the belowground C storage is not always consistent, making further research on this topic necessary [11]. The second management strategy was clear cut followed by installing a new saplings plantation of oak. Clear cut is an intensive forest management technique to harvest wood, that can also be used to realize a fast forest conversion [11,12,23]. However, due to the abrupt changes in the soil (with large C losses) in site micro-climate and in the landscape following the clear cut, the application of this management method is questioned [11,23]. Studies showed that clear cut can result in a decrease in C storage mainly in the forest floor, but also in the mineral soil [11,24,25]. Scots pine is a widely cultivated species in Europe. Sites with similar management and environmental conditions (e.g., soil, climate) are commonly found not only in the studied region (Flanders), but also in other vast areas in Europe with maritime climate and coastal influence.
We put forward two hypotheses. The first is that the aboveground part of both reference pine stands store more C than the stands that underwent the studied management regimes, since the C storage increases with time, and the new management was associated with larger wood removal than the reference management [5,14]. The second hypothesis is that the belowground part of the reference pine stands would also store more C than the other stands. In fact, with the new management, the soil and the forest litter floor are expected to lose more C, due to the thinning (decreasing the input of litter) and the clear cut (stimulating the degradation of organic materials due to disturbance and increased irradiation) [11,20,25,26].

2. Materials and Methods

2.1. Site Description

The study area is located in Brasschaat, Antwerp, in the Campine region of Belgium. The climate is temperate maritime, with an average annual temperature of 11 °C and an average annual precipitation of 938 mm (data for period 1991–2020). The coldest and warmest months during this period have a mean temperature of, respectively, 4 °C and 19 °C [27]. The study area is part of the Atlantic biogeographic region of Europe, originally dominated by broadleaved deciduous forests, but today extensively covered also by introduced evergreen coniferous species [28]. In the study region, the soils are strongly leached, nutrient-poor podzol types, originated under heather ecosystems [29].
Two study cases were used to investigate the impact of the two management regimes described above (i.e., thinning and group planting, clear cut, and new plantation) on the ecosystem C stocks when applied to Scots pine stands on sandy soils for ca. 15 years. In the first, the ecosystem C stocks of an extensively managed Scots pine stand (reference stand 1) were compared to the C stocks of a Scots pine stand with regular thinning and group planting of oak (thinned stand). In the second, the C stocks of a Scots pine stand under regular thinning (reference stand 2) were compared to the C stocks of a young oak stand planted after clear cut of a pine stand very similar to the former one (clear cut stand).
First study case—thinning with group planting. The two stands investigated in this regard are situated in De Inslag (51°18′ N, 4°31′ O, Figure 1a) and were planted in 1929 with Scots pine. After a similar (extensive) management and a last thinning in 1999 (removing about 30% of the trees [29]), one stand was further extensively managed with minimal intervention (De Inslag Ref., 1.76 ha, [30]). On the other hand, the other stand (1.80 ha) underwent a more intensive management with three thinnings in 2005, 2012 and 2015 (De Inslag Thin.). Following the yield tables, it was estimated that each thinning removed about 8% of the standing trees [31]. The thinning targeted the less vigorous trees following the “future trees” strategy [12]. Furthermore, in the latter stand, in 2012, pedunculate oak saplings were planted in understory gaps (group planting) generally in groups of 19–27 trees [19]. Rhododendron (Rhododendron ponticum L.) was also removed to decrease competition with the young, planted oaks [18,32]. In both stands the understory mainly consisted of young trees and shrubs. The soil in De Inslag is classified as a Endogleyic Hypobrunic Albic Hypoluvic Arenosol (Hyperdystric) according to the World Reference Base (WRB) for Soil Resources or medium wet (Zdg) to wet sandy podzol (Zeg) following the local Belgian classification [30,33]. De Inslag Ref. is intensively investigated within ICOS Belgium [30,34].
Second study case—clear cut with replanting. The two stands investigated here are a mature Scots pine stand and a young pedunculate oak stand planted after a clear cut, in De Mik (51°18′ N 4°33′ E, Figure 1b). In the first stand (Mik Ref.), pine was planted in 1943 on rows of artificially elevated ground (ca. 0.5 m) on a parcel of 3 ha. For the study only 1.5 ha was considered. Six thinnings followed until 2015. The second stand (Mik Oak; 1.52 ha) was planted in 1922 with pine and was clear cut in 2006 at an age of about 80 years (similar to the age of Mik Ref. at the time of measurements). Subsequently, young oak trees were planted with a density of 2500 trees ha−1. In 2017, pre-commercial thinning halved the tree density. In the Mik Ref. the understory consisted of grasses (mainly purple moor-grass, Molina caerulea L. Moench), while in the oak stand little understory was present. The soil is an Abrubtic, Gleyic Arenosol (Dystric) following the WRB, also classified as a medium wet to wet sandy podzols with a clay-sand substrate (w-Zdg, w-Zeg to w-Zfg) following the Belgian classification [33].
An overview of the different characteristics of all four studied stands is presented in Table 1. Note that De Inslag and De Mik are less than 3 km apart from each other at bird’s eye view.

2.2. Data Collection

The data were collected December–March 2020–2021. To determine the ecosystem C stocks at each studied stand, the following main components of the ecosystem were considered: the phytomass, the litter, and the mineral soil [30]. Phytomass was further divided into aboveground biomass of the overstory, understory, and belowground biomass. Trees were considered as part of the understory when their stem diameter (measured at 5 cm aboveground) was smaller than 5 cm. Litter was divided into forest floor and dead, woody debris. The mineral soil was sampled in three layers (0–10 cm, 10–30 cm, and 30–60 cm) [30].
In every stand, four center points were randomly chosen to mark the plots. Using these center points, three concentric circles with different radii were marked. In the first circle with a radius of 15 m (R15 plot), we determined the C pools (i) in the aboveground biomass of the overstory and (ii) the coarse roots with a diameter greater than 5 mm. The second circle had a radius of 5.6 m (R5.6 plot). Coarse, woody debris with a diameter greater than 10 cm was determined in these plots. In the third circle with a radius of 2.8 m (R2.8 plot) we determined the C storage in the understory and in the fine woody debris. To do so, in each wind direction and in the center of the plot, the vegetation was harvested in subplots (n = 5) of 0.3 × 0.3 m or 1 m2 depending on the homogeneity and abundance of the understory present. For each site, in one of the four R5.6 plot, 16 points were randomly selected to collect the samples for the C pools in (i) forest floor, (ii) the fine and medium roots (with a diameter less than 5 mm) and (iii) the mineral soil. An overview of the plots in which the components were observed is presented in Figure 2.

2.3. Measurements of Phytomass

2.3.1. Aboveground Biomass of the Overstory

The aboveground biomass of the overstory was determined using allometric relations (see Section 2.7.1) using data of height (H), diameter at breast height (DBH, at 1.30 m height) and species of each tree in the 15 m radius plots [30,36]. To determine the tree height (m), Field-map was used (IFER, Jilove u Prahy, Czech Republic). This was equipped with a rangefinder and an inclinometer to measure the height of the trees using trigonometry. The DBH (cm) was calculated by measuring the circumference manually with a measuring tape.

2.3.2. Understory

The understory consisted of grasses, mosses, ferns, shrubs and young trees with a stem less than 5 cm in diameter, Ø < 5 cm, measured 5 cm above ground. The biomass of young trees with a stem diameter between 1 cm and 5 cm was determined using allometric relations derived from the literature (Appendix A, Table A1). The height, stem diameter measured 5 cm above ground were determined as described above for the overstory. The biomass of the Rhododendron was determined using a relationship between biomass and volume determined in this study (Appendix B).
The biomass of the grasses, mosses, ferns, shrubs, and young trees with a diameter of less than 1 cm, was determined using destructive measurements [30,36]. Living and dead standing biomass was collected and further separated into (i) grasses, mosses and ferns and (ii) shrubs and young trees (Ø < 1 cm). The dry biomass was determined after drying at 70 °C for at least 48 h. For the data analysis, the average of the five subplots was used for each plot (see A2.8m in Figure 2).

2.3.3. Belowground Biomass

The belowground biomass consisted of the roots of the overstory and of the understory. These were separated in coarse, medium, and fine roots. The coarse roots were roots with a diameter larger than 5 mm and were estimated for the overstory and the young trees (1 < Ø < 5 cm) in each R15 and R2.8 plot, respectively. For the overstory, allometric relations (see Section 2.7.1) were used [30,36]. Due to a lack of allometric relations for young trees, for each species, their coarse root biomass was derived from data of aboveground biomass and the coarse root/shoot ratio of the mature trees (the ratios used are presented in Appendix A, Table A1). Bias was introduced by using these ratios (that change with age), but the relative importance of the C stock in the coarse roots of young trees was so small (see Section 3.1.3 in Results) that this approximation had little impact [37].
The fine and medium roots consisted of roots with diameter less than 5 mm and were measured in 16 points of one R5.6 plot in each site [30,38]. After the removal of the forest floor (see Section 2.4.1), the samples (n = 16) were taken with a root auger (Ø 7.85 cm, height 14.9 cm) for two layers at depth 0–10 cm and 10–30 cm (> 90% of fine roots reside above this depth [29]). At four of the 16 points, the 30–60 cm layer was also sampled. In the laboratory, the roots were manually picked out of the soil samples [39]. While sorting, dead and live roots were not separated. Thus, fine and medium root biomass also included their litter component. After washing, the roots were sorted into two classes: (a) medium roots, with a diameter between 2 mm and 5 mm and (b) fine roots, with a diameter less than 2 mm. The roots were dried at 70 °C for at least 48 h, after which they were weighed.

2.4. Litter

2.4.1. Forest Floor

The forest floor was considered as the thick layer (ca. 5 cm) consisting of fallen leaves, needles, pinecones, acorns, and partly decomposed organic material. This excludes the dark decomposed humus, which was considered the A-horizon of the mineral soil. Woody debris was collected separately (see Section 2.4.2).
For each stand, at each of the 16 points in the R5.6 plot, a square of 0.25 × 0.25 m of the forest floor was collected. The material was then dried at 40 °C for 7 days and weighed [30,38].

2.4.2. Woody Debris

The woody debris was separated in coarse (Ø > 10 cm) and fine (Ø < 10 cm) woody debris [36]. Fine woody debris (or fine wood, e.g., small branches) was collected for each R2.8 plot, in the same five subplots (of 0.09–1 m2) used to collect the understory. It was dried at 70 °C and weighed [36]. The biomass of the coarse woody debris (coarse wood, e.g., large branches, logs) was determined in the four R5.6 plots. The circumference and length of each wood debris within the plot were measured. To determine the biomass of each piece of coarse wood, the stadium of decomposition was accounted for [30,36]: 1—solid, not rotten; 2—solid in the center, rotting starting on the outside ; 3—rotten to the center. For each stadium, 3 to 5 pieces of wood were collected and their volume was determined. After drying at 70 °C, these coarse wood pieces were weighed and the density of the three stadia of decomposition was calculated. By using the length, circumference, and estimated density of the decomposition level, the dry biomass of the coarse wood in each plot was calculated.
In De Mik also (i) broken dead standing trees, (ii) stumps and (iii) dead wood left behind after the pre-commercial thinning of oaks, were also present and thus also considered as coarse woody debris. The biomass was determined as above, from measurements of volume and wood density. For standing dead trees and stumps, the biomass of the (dead) coarse root system (determined with allometric relationships assuming a decomposition level of 1) was also added to the coarse wood.

2.5. Mineral Soil

Soil samples of layers 0–10 cm and 10–30 cm were taken at the same 16 points of one R5.6 plot measured for the forest floor and fine root biomass. After removing the forest floor, the samples to measure the soil C were taken with a soil auger. Furthermore, for 4 sampling points (one in each quadrant of R5.6 plot) samples from the layer 30–60 cm were also collected.
For each layer and measuring point, two soil samples were taken: one to determine the bulk density and the other to determine the soil C content. To determine the bulk density, a sample was taken using a bulk density auger with Kopecky rings. These samples were dried for 48 h at 105 °C and weighed. The density of the soil was calculated using the dry mass value and the volume of the auger (0.00010 m3) [38]. By multiplying by the thickness of the layer, the mass of soil per surface unit could be calculated. Samples for soil C concentration were collected close to the sampling point of the bulk density with a smaller auger (Ø 2.45 cm). The soil samples for C analysis were dried at 40 °C for about 7 days [38].

2.6. C Content and C Stocks

For the phytomass and the woody debris a C content of 0.5 g C per g dry biomass was assumed [30]. The C content of the samples of the forest floor and mineral soil was analyzed using dry combustion. The dried samples were first sieved to pass a 2 mm mesh to remove roots and stones. Next, soil samples were ground using a ball mill (Fisher Bioblock Scientific Retsch), whereas dried forest floor samples were ground using a centrifugal mill (Model ZM 200, Retsch GmbH, Haan, Germany) with a sieve of 0.25 mm mesh size. The ground samples were analyzed for C content by dry combustion, based on the Dumas-method using an elemental analyzer (Model FLASH 2000, Interscience, Louvain-la-Neuve, Belgium).
Using the bulk density of the soil and the soil C content (g C g dry mass−1), the C stored in each layer of the mineral soil was calculated (tC ha−1). The C stocks (tC ha−1) in phytomass, forest floor, and woody debris were calculated using the dry mass (t ha−1) and the C content.

2.7. Data Analysis

2.7.1. Allometric Relationships

For Scots pine, allometric equations were available for the study area [40]. For the other species, we used equations determined for other Atlantic or Central European regions and for trees of comparable age as the investigated trees. In case equations for a representative region or age was not available, the relationships of a closely related species were used (these approximations were only applied for young trees of less common species). For all deciduous species, the use of allometric relationships introduced some bias, as some (dead) leaves of the current year were also considered within the forest floor pool. However, litter traps were not available for more accurate estimation and the bias is assumed minor. All details information about the relationships used are reported in Appendix A (Table A1) and Appendix B.

2.7.2. Statistics

For the statistical analyses, the software program IBM SPSS Statistics (Version 28.0, IBM Corp., Armonk, NY, USA) was used, while R (Version 4.1.1, R Foundation for Statistical Computing, Vienna, Austria) was used for data visualization [41,42]. The C stocks of all components characterizing phytomass, litter, and mineral soil were compared between the two stands at De Inslag and the two stands at De Mik separately. For C stocks in the fine roots and the mineral soil, the soil layers were considered both separately, as well as aggregated.
When data presented a normal distribution (Shapiro–Wilk test) and homoscedasticity (Levene’s test), an independent, two-tailed t-test was performed. In case of non-normal distribution, the non-parametric Mann–Whitney U-test was applied. When the homoscedasticity was not met, the Welch test was applied. The difference was considered significant when the p < 0.05.
This analysis could not be applied to the total C stock of the woody debris, the mineral soil (0–60 cm), total C stock in fine and medium roots (for each layer and in total depth) and of the total ecosystem C stock, due to different sample size or plot size of the components making up the aggregated values. Testing was therefore based on comparing the means and their standard deviation (SD) [43]. For example, the following equations were used to determine the standard deviation of the mineral soil C stock (SDmineral soil) and of the total ecosystem C stock (SDtotal ecosystem).
S D C   m i n e r a l   s o i l = S D C   0 10   c m 2 + S D C   10 30   c m 2 + S D C   30 60   c m 2
S D C   e c o s y s t e m = S D C   p h y t o m a s s 2 + S D C   f o r e s t   f l o o r 2 + S D C   w o o d y   d e b r i s 2 + S D C   m i n e r a l   s o i l 2
With the mean values of an aggregated C stock (e.g., mineral soil, total roots, total ecosystem) of two stands being m1 and m2 (with m1 < m2), and SD1 and SD2 their respective standard deviations, the values were considered to be significantly different at p < 0.05 when m1 + 1.96 × SD1 < m2 − 1.96 × SD2 [43].

3. Results

3.1. Phytomass

3.1.1. Overstory

In De Inslag, the overstory stored at the reference stand 93.3 tC ha−1 and at the thinned stand 75.3 tC ha−1, which were not significantly different (p = 0.15; Table 2). In both stands, Scots pine contributed to about 85% of the overstory biomass. In the reference stand, the most important species next to Scots pine, was birch (Betula pendula Roth; accounting for 10% of the overstorey biomass), whereas in the thinned stand, it was pedunculate oak (7%) (Figure 3).
In De Mik, the C stock in the aboveground tree biomass was 69.1 and 64.2 tC ha−1 in the reference and oak stand, respectively (not significantly different; Table 2). At Mik Ref., Scots pine made up 97% of the biomass. At Mik Oak, pedunculate oak made up 77% of the biomass.

3.1.2. Understory

In De Inslag, thinning and new planting did not result in a significant difference in C stock in the young trees and shrubs (1 < Ø < 5 cm) with 4.1 and 4.6 tC ha−1 in the reference and thinned stand, respectively. In the understory, the largest amount of C was stored in this compartment. Young trees and shrubs with Ø < 1 cm (< 0.1 tC ha−1) and grasses, mosses, and ferns (0.2 tC ha−1), accounted for a negligible part of the ecosystem C stock.
However, the grasses, mosses, and ferns represented the largest share of the understory biomass in De Mik, with values at Mik Ref. (1.38 tC ha−1) significantly larger than at Mik Oak (0.35 tC ha−1) (p < 0.001; Table 2). Young trees and shrubs with Ø < 1 cm and 1 < Ø < 5 cm accounted for a negligible portion of the ecosystem C stock at both stands (< 0.1 tC ha−1).

3.1.3. Belowground Biomass

In De Inslag, C stock in the coarse roots of the overstory trees (ca. 9–11 tC ha−1) and the young trees (ca. 1 tC ha−1) were both not significantly different between stands (Table 2). In De Mik, the C stock in the coarse roots of the overstory trees was 7.8 and 10.9 tC ha−1 in the reference and oak stand, respectively, with a significant difference (p < 0.001; Table 2). Understory trees were largely absent at De Mik, and thus did not contribute to the belowground biomass.
An overview of the C stored in the fine and medium roots for the different soil layers in both stands of De Inslag and De Mik is presented in Table 3. Overall, for the aggregated fine and medium roots in De Inslag, the total sampled layer (0–60 cm) stored 4.7 ± 2.6 and 7.4 ± 3.3 tC ha−1 (mean ± SD) in the thinned and reference stand, respectively, without significant difference between stands. The fine and medium roots accounted for about 50% of the total root C stock. In De Inslag, the C storage in the fine roots was only significantly different for the layer 30–60 cm, with larger C stock in the reference than in the thinned stand (p = 0.03). For the medium roots, the C stock of the separated layers was not significantly different between stands. Though, in the layer 0–10 cm, there was a trend toward a larger C stock in De Inslag Ref. (p = 0.06; Table 3). In the aggregated fine and medium roots, the C stock showed only a significant difference in the layer 30–60 cm, as in the fine roots, with larger values in the reference stand (3.0 ± 0.9 tC ha−1) than in the thinned stand (0.7 ± 0.2 tC ha−1) (p = 0.040; Table 2).
The total C in the aggregated fine and medium roots in De Mik was 4.6 ± 2.0 and 6.9 ± 2.8 tC ha−1 (mean ± SD) in the reference and oak stand, respectively, without significant difference between stands. However, in De Mik, the C stock in the fine roots in the layer 0–10 cm was lower in the reference than in the oak stand (0.8 ± 0.1 and 1.3 ± 0.1 tC ha−1, respectively) (p < 0.001; Table 3) and a similar trend (p = 0.05) was observed in the same layer for the aggregated fine and medium root C stock (Table 2). There were no significant differences in the root C stocks in the layers of 10–30 cm (1–3 tC ha−1 ) and 30–60 cm (1–2 tC ha−1) (Table 2).

3.2. Litter

3.2.1. Forest Floor

In De Inslag, the amount of C stored in the forest floor in the extensively managed stand and in the thinned stand was about 42–45 tC ha−1 (the difference was not significant; Table 2). In De Mik, the amount of C stored in the forest floor of the pine stand was twice as large as the amount of C stored in the forest floor of the oak stand (p < 0.001; Table 2).

3.2.2. Woody Debris

Regarding the C stored in the coarse woody debris, there was no significant difference between both stands in De Inslag, neither between both stands in De Mik (Table 2). In De Inslag and De Mik, the C stocks of fine woody debris also did not show a significant difference between stands, but Mik Oak presented marginally significant (p = 0.08) larger fine woody debris than Mik Ref (Table 2). Overall, in De Inslag, a total of 5–6 tC ha−1 was stored in the dead wood, whereas in De Mik, this amount was about 4–7 tC ha−1 (without significant differences between stands).

3.3. Mineral Soil

In De Inslag, there was no significant difference in C stocks in the total sampled soil (0–60 cm) between the extensively managed stand (107 ± 22 tC ha−1) and the thinned stand (105 ± 31 tC ha−1). In addition, in De Mik, the total C stored in the soil did not differ significantly between the pine stand (185 ± 31 tC ha−1) and the oak stand (205 ± 56 tC ha−1). The results are given with their corresponding SD.
In De Inslag, there was, however, a difference in the layer 10–30 cm, with larger C stock in the extensively managed stand than in the thinned stand (p < 0.001; Table 2). This difference was related to C content (C%) and bulk density, which tended to be larger in the reference stand (Table 4). The layers 0–10 cm (40-55 tC ha−1) and 30–60 cm (22 tC ha−1) did not present a significant difference between stands. Note that the aggregated layer of 0–30 cm showed no significant difference in C stocks between the extensively managed stand (84.6 ± 5.9 tC ha−1) and the thinned stand (82.5 ± 8.4 tC ha−1) (p = 0.84).
In De Mik, the layer 0–10 cm presented larger C stock in the oak stand than in the pine stand (p = 0.036; Table 2). The C stock in the layer 10–30 cm also tended to be larger in the oak stand than the pine stand (Table 2). In fact, when considering the aggregated layer 0–30 cm, the C stored in the pine stand (118.1 ± 7.5 tC ha−1) was significantly lower than in the oak stand (142.6 ± 6.5 tC ha−1) with p = 0.020. In both layers, 0–10 cm and 10–30 cm, the difference in C stocks between stands was associated with a marginal difference (p = 0.05–0.06) in C%, with Mik Oak presenting larger values than Mik Ref. (Table 4). In the layer 30–60 cm the C stock did not differ significantly between the two stands (60–65 tC ha−1), though the bulk density was marginally larger at Mik Ref. (Table 4).

3.4. Total Ecosystem

The total C stock of the study stands is reported in Figure 4. For both De Inslag and De Mik, the total C stock in the ecosystem was not significantly different between stands with different management (Figure 4).
In De Inslag, the aboveground C plus woody debris C represented about 30–38% of the C stored in the ecosystem. On the other hand, soil and roots represented 43–44% of the ecosystem C stock. In the thinned stand, the lower proportion of C in the soil layer 10–30 cm than the extensively managed stand (13% vs. 15%, respectively) was partly compensated by a larger amount of C (even if not significant) in the 0–10 cm layer, so that no difference in C stock in the top 30 cm of soil was recorded. Differences between stands in the relative proportion of C in forest floor, coarse roots and the soil layer at 30–60 cm were also not relevant (Table 2).
At De Mik, only 20–25% of the C was aboveground, whereas the majority (60–70%) was belowground. The relative proportion of C stocks in the aboveground components, in coarse roots and in the deepest soil layer did not vary between stands (Table 2). However, there were differences between stands in the proportion of C in the forest floor and the soil layers until 30 cm (Figure 4). In fact, the pine stand had largest proportion of stored C in the forest floor (15% vs. 7%, for Mik Ref. and Mik Oak, respectively), whereas the oak stand had the largest proportion of C in the top 30 cm of soil (47% vs. 39%, for Mik Oak and Mik Ref., respectively), consisting of mineral soil and fine and medium roots.

4. Discussion

4.1. Aboveground Biomass

Thinning can significantly decrease the C stock in overstory biomass, depending on its intensity and the recovery time [21,22,44,45]. However, the results do not show an impact of thinning on the overstory C stock at De Inslag. This lack of significant effect could be attributed to a less intensive thinning, than in other studies [21,22,44]. Even though the exact data of the trees or wood volume removed was not available, following the standard yield tables it was estimated that only about 18% of the standing trees was thinned in the 15 years before measurements [29]. Another factor partly compensating for the loss of biomass due to the tree harvest could have been the improved growth of the trees left after thinning. This was shown by the difference in the DBH of the pines in De Inslag with an average of 37.3 and 31.1 cm in the thinned and reference stand, respectively. This was observed despite a considerable amount of younger and smaller pine trees in the thinned stand, due to natural regeneration, which reduced the mean DBH and increased the data uncertainty (Figure 3). The oaks planted in 2012 at De Inslag Thin. did not contribute substantially to the overall C storage (ca. 7% of the overstory biomass), but their relevance can be assumed to increase in the next decades as the conversion to an oak forest continues.
Additionally, in De Mik, there was no significant difference in the aboveground C stock between both stands. Generally, younger trees store less C than older ones, and pioneer coniferous evergreen trees grow faster than late successional broadleaved deciduous trees [5,46,47]. However, the standard applied guidelines for the planting and the first thinning assured that, after 15 years, the young oak stand already stored as much C as the ca. 80-year-old pine stand.
Both in De Inslag and De Mik, the understory accounted for a minor contribution to the C storage in the ecosystem. In De Inslag, young trees and shrubs represented the most important C stock in the understory. In De Mik, grasses, mosses, and ferns represented the largest understory C stock in the open Mik Ref. stand, whereas the closed Mik Oak stand had practically no understory. For De Inslag Ref., C storage in grasses, mosses, and ferns in the peak of the 2021 growing season was 0.18 tC ha−1 according to ICOS Belgium, which compares well to our results (0.21 tC ha−1). This might indicate that the standing dead understory biomass measured here outside of the growing season was an acceptable proxy of the peak understory biomass. On the other hand, it is possible that the spatial heterogeneity of the understory confounds the results.

4.2. Roots

Both in De Inslag and De Mik, there was no effect of forest management on the C stored in the fine and medium roots in the total sampled depth (0–60 cm). However, in De Inslag, in the layer 30–60 cm, the C stock of the fine and medium roots decreased alongside management intensity. It is likely that the amount of (deeper) pine roots decreased with thinning, and this could not have been compensated by the young oaks. On the other hand, in De Mik, in the layer 0–10 cm, there was a significantly higher amount of roots in the oak stand, which could be attributed to the high tree density. Regarding coarse roots, the significantly lower C storage in the pine stand in De Mik could be the result of a lower coarse root/shoot ratio in the pine stand compared to the oak stand (0.11 and 0.17, respectively). In fact, the root/shoot ratio can decrease with age and with lower stand density ([37] for root/shoot ratio considering total roots).

4.3. Forest Floor

In De Inslag, the effect of different forest management was not noticeable on the forest floor C stock. This could be explained by the similar composition of the litter in both stands, each dominated (for ca. 80–85%) by Scots pine.
In Mik Oak, the litter, and as a result, the forest floor composition, changed with the sudden conversion after the clear cut. The litter of a broadleaved deciduous stand has better quality than the one of an evergreen coniferous stand. As a result, the litter in the former is broken down faster by soil biota than in the latter [6,11]. Pine litter is recalcitrant, partly because of the low pH and low calcium content [46]. Secondly, the soil of pine stands typically contains fewer earthworms, which are important for the breaking down of litter [46]. Overall, due to the change in litter and soil ecology following the clear cut and the new planting, the C stock of the forest floor in Mik Oak was lower than in Mik Ref. In addition, the clear cut had likely promoted additional breakdown of forest floor organic matter through the increase in insolation, and thus soil temperature, that it created [11].

4.4. Mineral Soil

The amount of C stored in the soil can strongly depend on the intensity and kind of forest management [1,5,48,49]. For example, intensive disturbance of the soil can foster the breakdown of organic matter and, thereby, the loss of soil C [11,48]. However, the effect of thinning on the soil C stocks can be ambiguous, and literature shows that thinning does not necessarily influence the organic C stored in the mineral soil [11,21,50,51]. This corresponds to our observations at De Inslag. In the soil layer of 10–30 cm, however, the C stock in the reference stand was larger than in the thinned stand. A possible reason for this could be a smaller belowground input of C from the trees in the thinned stand, which is corroborated by the smaller amount of roots in the soil layer of 30–60 cm. On the other hand, we note that the C stock in the 0–10 cm layer was, although non-significant, larger in the thinned stand, balancing out the difference in soil C stock when the top 30 cm of soil was considered. Therefore, the different management could induce a modification in the C allocation pattern among pools rather than a net C loss, which warrants further investigation over longer time scales.
In De Mik, the C stock in the topsoil layer (0–30 cm) was larger in the oak than in the pine stand. This is likely the result of a combination of a better quality of the litter (that is more easily broken down) and of better-developed soil ecology, allowing the transfer of the organic matter from the forest floor to the mineral soil [6,46]. Another factor contributing to the larger topsoil C stock in Mik Oak could have been the increased fine rootstock in the layer 0–10 cm. Due to their turnover, more fine roots could result in larger C input into the soil [52]. The relocation of C from litter to the mineral soil did not reach the deepest measured layer of 30–60 cm, which stored the same amount of C in both stands. This could be explained by the young age of the oak stand. Literature reports that modifications in C content in mineral soil layers can occur over long-time scales and that there are still many unknown aspects of these dynamics [11,53].
The total amount of soil C stock was remarkably higher at De Mik when compared to De Inslag. This difference could be attributed to the soil C content (C%), which was larger at De Mik, particularly for the deeper layers of the soil (10–60 cm). Due to the larger clay fraction, it is likely that more water is retained and less soil organic matter is decomposed at De Mik than at De Inslag [54].

4.5. Woody Debris

It was expected that thinning would result in less woody debris due to the lower density of trees, as well as the lower tree competition [12]. This was not reflected in the results for De Inslag. This might be due partly to the fact that, in the reference stand, the forest canopy was very open, so competition between the trees was avoided even in the stand that was not thinned. Additionally, woody debris (crown wood) was left on site after thinning.
In addition, in De Mik, the management strategy did not have a major influence on the C stored in the woody debris. The slightly larger C stock in the fine woody debris in the oak stand could be attributed to the higher amount of branches and the higher tree density at the oak stand, both resulting in higher branch mortality [55].

4.6. Total Ecosystem

For De Inslag Ref. (the only stand with previous measurements), our data match previous C stock assessments. For the mineral soil, the C stock was reported to be 115 tC ha−1 in layer 0–100 cm in 1999 [29] and 68 tC ha−1 in layer 0–30 cm in the period 2002–2010 [30]. This is similar to the 85 tC ha−1 and 107 tC ha−1 in the 0–30 cm and 0–60 cm layers, respectively, measured in De Inslag Ref. in the current study. The C stock in the overstory biomass was similar in 1999 [29] (90 tC ha−1), 2002–2010 [30] (76–90 tC ha−1), and in this study (93 tC ha−1). Lastly, the litter/forest floor stored 26, 30–34, and 42 tC ha−1 in 1999, 2002–2010, and 2021, respectively. Overall, these observations confirm that our measurement protocol was of satisfactorily quality.
In De Inslag, as well as in De Mik, there was no noticeable effect of the different forest management on the C stock of the total ecosystem. This contrasts with the expectation that the forest C stock would decrease after the new management strategies, i.e., thinning or clear cut. According to literature, clear cut can have a large influence on the C stored in the mineral soil and, mainly, the forest floor, with recovering time depending on multiple factors (such as soil type, climate, tree species, etc.) and up to five decades [11,21,25]. In this case study, the C stock has recovered in, at maximum, 15 years.
The management seems to have modified the C allocation pattern in the forest floor and soil rather than the absolute values of the total ecosystem C stock. This was evident at De Mik, where the clear cut and new planting resulted in less C on the forest floor (−7% considering the total ecosystem C stock) and more C in the top 30 cm of soil (+9%). Although less obvious, this effect was also present in De Inslag, where management seemed to influence the C distribution among the topsoil layers.
Based on the similarity in soil and maritime climate influence, the results and insights obtained here are also of relevance for the widespread Scots pine forests in other areas of the Atlantic region (e.g., The Netherlands, the western part of Germany, and Norway) [56,57], but also in the coastal areas of more continental European regions (e.g., southern Baltic coast) [58].

4.7. Methodological Limitations

We need to stress that the belowground sampling was performed only at one plot of about 100 m2 at each stand. Therefore, the belowground sampling was not suited to capture the intra-stand belowground spatial variability. Previous research indicated, however, that the variability in C stocks of soil and litter is relatively low at the De Inslag (confidence interval about 6–13% of the mean), but it is larger for fine roots (confidence interval about 30–38% of the mean) [30]. Fine roots data should therefore be treated with caution, as well as the soil data of the 30–60 cm layer, which were analyzed via four replicas only. For the roots, there was no separation made between live and dead roots. Different live vs. dead roots ratios across stands were possible if the management strategy affected the root turnover rate, but this is unknown.
Finally, note that, in our analysis, the four plots per stand from which the aboveground data were derived were considered independent. A more refined approach would have been achieved by selecting multiple stands for each management strategy. However, such stand replicas were not present in the study area.

5. Conclusions

The effect of typical management strategies on Scots pine stands on sandy soil was clarified by investigating gradual thinning with concurrent planting of pedunculate oak, on the one hand, and clear cut followed by a new planting of pedunculate oak, on the other hand. Contrary to the hypotheses, both management methods did not show a significant effect on the total C stock of the ecosystem or above- and belowground C stocks taken separately, but the distribution of C within the ecosystem C pools did change. In particular, after the clear cut and new plantation, important changes in the C stock of the mineral soil (increase) and the forest floor (decrease) could have occurred. The effect of thinning and concurrent planting had less impact on the ecosystem C storage, although modifications in the distribution of the C in the top 30 cm soil layer were observed. Because of the significant ecological and economic value of Scots pine, more information on the impact of forest management on the ecosystem C stocks of maritime Scots pine forests on poor sandy soils was needed. Therefore, the data presented in this study are of importance.

Author Contributions

Conceptualization, Methodology, Funding acquisition, M.C.; Formal analysis, F.V.D., H.M., and T.H.; Investigation, F.V.D., H.M., L.L., T.D.M., and M.P.-E.; Visualization, F.V.D., Writing—original draft, F.V.D., H.M., and M.C., Writing—review and editing, F.V.D., H.M., T.H., L.L., T.D.M., M.P.-E., M.R., B.G., I.A.J., K.V., and M.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Agentschap voor Natuur en Bos (ANB, AGENCY FOR NATURE AND FOREST) and Fonds Wetenschappelijk Onderzoek (FWO, research Foundation-Flanders).

Data Availability Statement

Data can be retrieved from the first and last authors, F.V.D. and M.C.

Acknowledgments

We thank the Agency for Nature and Forest (ANB) of the Flemish Community for support to Freke Van Damme and Hana Mertens. Hereby, we want to give special thanks to two ANB officers: Mario De Block and Dirk Leyssens. The authors also acknowledge the financial support of the Research Foundation-Flanders to the ICOS infrastructure.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Allometric Relations

Table A1. Allometric relations used to calculate the aboveground biomass and biomass of the coarse roots of the overstory trees and understory trees.
Table A1. Allometric relations used to calculate the aboveground biomass and biomass of the coarse roots of the overstory trees and understory trees.
Species Aboveground Biomass (AGB) (kg/Tree) Biomass Coarse Roots (BGB) (kg/Tree) Remarks
Overstory(D > 5 cm)
Acer pseudoplantanus L.AGBleave s= 10(−1.1619)+1.2237 × LOG(DBH)
(Q. robur) [59]
AGBtotal = e(−2.7018)+2.5751 × ln(DBH) + AGBleaves
[60]
BGB = 0.040113 × DBH2.227842
(Q. robur) [61]
a
Amelanchier arborea FernaldAGBLeaves = 10(−1.1619)+1.2237 × LOG(DBH)
(Q. robur) [59]
AGBTotal = 0.0927 × DBH2.5097+ AGBLeaves
(Q. robur) [62]
BGB = 0.040113 × DBH2.227842
(Q. robur) [61]
a
Betula pendulaAGBLeaves, Inslag = 0.80 kg/tree
AGBLeaves, Mik = 1.06 kg/tree
[63]
AGBBranches = 0.0742 × DBH2,24
AGBStem = 0.193 × DBH2,25
AGBTotal = 0.2511 × DBH2.29 + AGBleaves
[64]
BGB = 0.04582 × DBH2.23951
[61]
AGBTotal: De Inslag
AGBBranches and AGBStem: Mik

a
Castanea sativa Mill.AGBLeaves = 10(−1.1619)+1.2237 × LOG(DBH)
(Q. robur) [59]
AGBTotal = 0.0927 × DBH2.5097+ AGBLeaves
(Q. robur) [62]
BGB = 0.040113 × DBH2.227842
(Q. robur) [61]
a
Corylus avellanaAGBTotal = 0.068 × DBH2.745
[65]
BGB = 0.040113 × DBH2.227842
(Q. robur) [61]
a
Fagus sylvaticaAGBLeaves = 0.0167 × DBH2.951 × H1.101
AGBBranches = 0.0114 × DBH3.682 × H1.031
AGBStem = 0.0109 × DBH1.951 × H1.262
AGBTotal = 0.0306 × DBH2.347 × H0.590
[66]
BGB = e(−4.1302)+2.6099 × ln(DBH)
[67]
AGBTotal: De Inslag
AGBLeaves AGBBranches and AGBStem: Mik
Frangula alnus Mill.AGBTotal = 0.121 × DBH2.480
[65]
BGB = 0.040113 × DBH2.227842
(Q. robur) [61]
a
Ilex aquifoliumAGBLeaves = 10(−1.1619)+1.2237 × LOG(DBH)
(Q. robur) [59]
AGBTotal = 0.0927 × DBH2.5097+ AGBLeaves
(Q. robur) [62]
BGB = 0.040113 × DBH2.227842
(Q. robur) [61]
a
Pinus sylvestris L.AGBbranches = 0.0022 × DBH2.9123
AGBLeaves = 0.0045 × DBH2.2372
AGBStem = 0.1227 × DBH2.3272
AGBtotal = AGBbranches + AGBLeaves + AGBStem
[40]
BGB = 0.3399 × DBH1.4728
[40]
Prunus serotina Ehrh.AGBLeaves = 10(−1.1619)+1.2237 × LOG(DBH)
(Q. robur) [59]
AGBTotal = 0.0927 × DBH2.5097+ AGBLeaves
(Q. robur) [62]
BGB = 0.040113 × DBH2.227842
(Q. robur) [61]
a
Pseudotsuga menziesii (Mirb.) FrancoAGB = e(1.620)+2.410 × ln(DBH)
[66]
BGB = 0.040113 × DBH2.227842
(Q. robur) [61]
a
Quercus roburAGBLeaves = 10(−1.1619)+1.2237 × LOG(DBH)
[59]
AGBBranches = 0.0149 × DBH2.5994
AGBstem = 0.0722 × DBH2.5135
AGBTotal = 0.0927 × DBH2.5097+ AGBLeaves
[62]
BGB = 0.040113 × DBH2.227842
[61]
AGBTota: De Inslag
AGBBranches, AGBstem: De Mik
a
Quercus rubraAGBLeaves = 10(−1.1619)+1.2237 × LOG(DBH)
(Q. robur) [59]
AGBTotal = 0.0927 × DBH2.5097+ AGBLeaves
(Q. robur) [62]
BGB = 0.040113 × DBH2.227842
(Q. robur) [61]
a
Salix caprea L.AGBLeaves = 10(−1.1619)+1.2237 × LOG(DBH)
(Q. robur) [59]
AGBTotal = 0.0927 × DBH2.5097+ AGBLeaves
(Q. robur) [62]
BGB = 0.040113 × DBH2.227842
(Q. robur) [61]
a
Sorbus aucupariaAGBTotal = 0.164 × DBH2.047
(De Inslag) [65]
AGBLeaves = 10(−1.1619)+1.2237 × LOG(DBH)
(Mik, Q. robur) [59]
AGBBranches = 0.0149 × DBH2.5994
AGBstem = 0.0722 × DBH2.5135
(Mik, Q. robur) [62]
BGB = 0.040113 × DBH2.227842
(Q. robur) [61]
a
Taxus baccata L.AGBbranches = 0.0022 × DBH2.9123
AGBLeaves = 0.0045 × DBH2.2372
AGBStem = 0.1227 × DBH2.3272
AGBtotal = AGBbranches + AGBLeaves + AGBStem
(P. sylvestris) [40]
BGB = 0.3399 × DBH1.4728
(P. sylvestris) [40]
Young trees (D > 5 cm)
Betula pendula AGBLeaves = e(−1.6938)+1.7455 × ln(D)
(Q. robur) [68]
AGBTotal = (0.259 × D2.132)/1000 + AGBLeaves
[69]
BGB = 0.16 × AGB
Frangula alnusAGBLeaves = e(−1.6938)+1.7455 × ln(D)
[68]
AGBTotal = (0.027 × (D)2.769)/1000 + AGBLeaves
[69]
BGB = 0.21 × AGB
Quercus roburAGBLeaves = e(−1.6938)+1.7455 × ln(D)
[68]
AGBTotal = (0.027 × (D)2.769)/1000 + AGBLeaves
[69]
BGB = 0.23 × AGB
Pinus sylvestris L.AGBLeaves = e(−5.478)+2,494 × ln(DBH)
[70]
AGBTotal = (0.015 × D2.881)/1000 + AGBLeaves
[69]
BGB = 0.17 × AGBNeedles: Eq. from literature uses DBH; here, D was used, resulting in a small error
Sorbus aucupariaAGBLeaves = e(−1.6938)+1.7455 × ln(D)
(Q. robur) [68]
AGBTotal = (0.143 × D2.064)/1000+AGBLeaves
[69]
BGB = 0.23 × AGB
RhododendronBiomass per volume-unit (Appendix B):
AGBTotal = 68,131.70 g/m3
Root/shoot = 0.68
[71]
BGB: root/shoot ratio of Phillyrea latifolia L
Dry aboveground biomass (AGB) and coarse roots (BGB) in kg, DBH in cm, stem diameter measured 5 cm above ground (D) in mm. a: Coarse roots: allometric relations of roots with diameter > 2 mm. The diameter class of 2–5 mm is counted twice (as coarse roots and as medium roots); compared to the total amount of C in the roots, this error is small.

Appendix B. Calculations of the Allometric Relationships for Rhododendron

The amount of dry weight (DW) of stems (or branches) plus leaves of Rhododendron was calculated from measurements of volume and density of the stem. Six full stems were harvested from a Rhododendron located outside the measured plots. The length and diameter at 5 cm above ground were determined for each stem so that the stem volume could be calculated. The fresh weight was also determined. Next, sub-samples of each stem were taken and weighed, after which they were dried and weighed again. By extrapolation, the DW of the total stem was calculated. Using the DW (g) and the volume (m3) of the stem, the stem density (g m−3) was calculated. The average density of Rhododendron was 68.132 kg m−3. Using this density, DW was determined from measurements of diameter and length for each of the several stems making up the Rhododendron shrubs in the measurements plots.

References

  1. Fahey, T.J.; Woodbury, P.B.; Battles, J.J.; Goodale, C.L.; Hamburg, S.P.; Ollinger, S.V.; Woodall, C.W. Forest carbon storage: Ecology, management, and policy. Front. Ecol. Environ. 2010, 8, 245–252. [Google Scholar] [CrossRef] [Green Version]
  2. Hui, D.; Deng, Q.; Tian, H.; Luo, Y. Climate change and carbon sequestration in forest ecosystems. In Handbook of Climate Change Mitigation and Adaptation, 2nd ed.; Chen, W.-Y., Suzuki, T., Lackner, M., Eds.; Springer International Publishing: Cham, Switserland, 2017; Volume 555, p. 594. [Google Scholar] [CrossRef]
  3. Falkowski, P.; Scholes, R.J.; Boyle, E.E.A.; Canadell, J.; Canfield, D.; Elser, J.; Gruber, N.; Hibbard, K.; Högberg, P.; Linder, S.; et al. The global carbon cycle: A test of our knowledge of earth as a system. Science 2000, 290, 291–296. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Campioli, M.; Vicca, S.; Luyssaert, S.; Bilcke, J.; Ceschia, E.; Iii, F.S.C.; Ciais, P.; Fernández-Martínez, M.; Malhi, Y.; Obersteiner, M.; et al. Biomass production efficiency controlled by management in temperate and boreal ecosystems. Nat. Geosci. 2015, 8, 843–846. [Google Scholar] [CrossRef]
  5. Anderson-Teixeira, K.J.; Herrmann, V.; Morgan, R.B.; Bond-Lamberty, B.; Cook-Patton, S.C.; Ferson, A.E.; Muller-Landau, H.C.; Wang, M.M. Carbon cycling in mature and regrowth forests globally. Environ. Res. Lett. 2021, 16, 053009. [Google Scholar] [CrossRef]
  6. De Deyn, G.B.; Cornelissen, J.H.; Bardgett, R.D. Plant functional traits and soil carbon sequestration in contrasting biomes. Ecol. Lett. 2008, 11, 516–531. [Google Scholar] [CrossRef]
  7. Cotrufo, M.F.; Wallenstein, M.D.; Boot, C.M.; Denef, K.; Paul, E. The Microbial Efficiency-Matrix Stabilization (MEMS) framework integrates plant litter decomposition with soil organic matter stabilization: Do labile plant inputs form stable soil organic matter? Glob. Chang. Biol. 2013, 19, 988–995. [Google Scholar] [CrossRef] [Green Version]
  8. Ameray, A.; Bergeron, Y.; Valeria, O.; Montoro Girona, M.; Cavard, X. Forest Carbon Management: A Review of Silvicultural Practices and Management Strategies Across Boreal, Temperate and Tropical Forests. Curr. For. Rep. 2021, 7, 245–266. [Google Scholar] [CrossRef]
  9. Bonan, G.B. Forests and climate change: Forcings, feedbacks, and the climate benefits of forests. Science 2008, 320, 1444–1449. [Google Scholar] [CrossRef] [Green Version]
  10. Pan, Y.; Birdsey, R.A.; Fang, J.; Houghton, R.; Kauppi, P.E.; Kurz, W.A.; Phillips, O.L.; Shvidenko, A.; Lewis, S.L.; Canadell, J.G.; et al. A large and persistent carbon sink in the world’s forests. Science 2011, 333, 988–993. [Google Scholar] [CrossRef]
  11. Mayer, M.; Prescott, C.E.; Abaker, W.E.; Augusto, L.; Cécillon, L.; Ferreira, G.W.; James, J.; Jandl, R.; Katzensteiner, K.; Laclau, J.-P.; et al. Tamm Review: Influence of forest management activities on soil organic carbon stocks: A knowledge synthesis. For. Ecol. Manag. 2020, 466, 118127. [Google Scholar] [CrossRef]
  12. den Ouden, J.; Muys, B.; Mohren, G.M.J.; Verheyen, K. Bosecologie en Bosbeheer; Acco: Leuven, Belgium, 2010. [Google Scholar]
  13. Van der Meijden, R. Heukels Flora van Nederland 23e Druk; Wolters-Noordhoff: Groningen, The Netherlands, 2005. [Google Scholar]
  14. Alonso, I.; Weston, K.; Gregg, R.; Morecroft, M. Carbon storage by habitat: Review of the evidence of the impacts of management decisions and condition of carbon stores and sources. In Natural England Research Reports; Number NERR043; Natural England: Pitlochry, UK, 2012. [Google Scholar]
  15. Adriaens, T.; Cartuyvels, E.; Denys, L.; Devisscher, S.; Oldoni, D.; Packet, J.; Provoost, S.; Scheers, K.; Soors, J.; Vandevoorde, B.; et al. Invasieve Exoten in Vlaanderen: Toestand en beleidsaanbevelingen. Uitgebreid achtergrondrapport bij het Natuurrapport 2020. In Rapporten van Het Instituut voor Natuur- en Bosonderzoek 2020; Instituut voor Natuur- en Bosonderzoek: Belgium, Brussel, 2020; Volume 41. [Google Scholar] [CrossRef]
  16. Agency for Nature and Forest (ANB). Meetvraag 3: Bestandsopbouw. Vlaamse Bosinventaris (1997–1999, 2009–2019). Available online: https://www.natuurenbos.be/vlaamse-bosinventaris/oppervlakteperbestandstype.html (accessed on 21 December 2021).
  17. Verheyen, K.; De Schrijver, A.; Wuyts, K.; Gielis, M.; Van Gossum, P.; Geudens, G.; Van Herzele, A.; De Boever, L.; Vanhellemont, M. Van dennenplantages naar een beloofd land?! Theoretische en praktische aspecten van bosomvorming. Silva Belg. 2007, 114, 20–26. [Google Scholar]
  18. Agency for Nature and Forest (ANB). Bosomvorming [Brochure]; Agency for Nature and Forest: Brussels, Belgium, 2008. [Google Scholar]
  19. Saha, S.; Kuehne, C.; Kohnle, U.; Brang, P.; Ehring, A.; Geisel, J.; Leder, B.; Muth, M.; Peterson, R.; Peter, J.; et al. Growth and quality of young oaks (Quercus robur and Quercus petraea) grown in cluster plantings in central Europe: A weighted meta-analysis. For. Ecol. Manag. 2012, 283, 106–118. [Google Scholar] [CrossRef]
  20. del Río, M.; Bravo-Oviedo, A.; Pretzsch, H.; Löf, M.; Ruiz-Peinado, R. A review of thinning effects on Scots pine stands: From growth and yield to new challenges under global change. For. Syst. 2017, 26, eR03S. [Google Scholar] [CrossRef]
  21. Zhou, D.; Zhao, S.Q.; Liu, S.; Oeding, J. A meta-analysis on the impacts of partial cutting on forest structure and carbon storage. Biogeosciences 2013, 10, 3691–3703. [Google Scholar] [CrossRef] [Green Version]
  22. Ruiz-Peinado, R.; Bravo-Oviedo, A.; Montero, G.; Del Río, M. Carbon stocks in a Scots pine afforestation under different thinning intensities management. Mitig. Adapt. Strateg. Glob. Chang. 2016, 21, 1059–1072. [Google Scholar] [CrossRef]
  23. Seedre, M.; Felton, A.; Lindbladh, M. What is the impact of continuous cover forestry compared to clearcut forestry on stand-level biodiversity in boreal and temperate forests? A systematic review protocol. Environ. Evid. 2018, 7, 28. [Google Scholar] [CrossRef] [Green Version]
  24. Davis, S.C.; Hessl, A.E.; Scott, C.J.; Adams, M.B.; Thomas, R.B. Forest carbon sequestration changes in response to timber harvest. For. Ecol. Manag. 2009, 258, 2101–2109. [Google Scholar] [CrossRef]
  25. Nave, L.E.; Vance, E.D.; Swanston, C.W.; Curtis, P.S. Harvest impacts on soil carbon storage in temperate forests. For. Ecol. Manag. 2010, 259, 857–866. [Google Scholar] [CrossRef]
  26. Hughes, J.W.; Fahey, T.J. Litterfall dynamics and ecosystem recovery during forest development. For. Ecol. Manag. 1994, 63, 181–198. [Google Scholar] [CrossRef]
  27. Royal Meteorological Institute of Belgium (KMI). Klimaatstatistieken van de Belgische Gemeenten: Brasschaat (NIS 11008); Royal Meteorological Institute of Belgium (KMI): Brussels, Belgium, 2020. [Google Scholar]
  28. Polunin, O.; Walters, M. Guide to the Vegetation of Britain and Europe; Oxford University Press: Oxford, UK, 1985. [Google Scholar]
  29. Janssens, I.A.; Sampson, D.A.; Cermak, J.; Meiresonne, L.; Riguzzi, F.; Overloop, S.; Ceulemans, R. Above-and belowground phytomass and carbon storage in a Belgian Scots pine stand. Ann. For. Sci. 1999, 56, 81–90. [Google Scholar] [CrossRef]
  30. Gielen, B.; De Vos, B.; Campioli, M.; Neirynck, J.; Papale, D.; Verstraeten, A.; Ceulemans, R.; Janssens, I.A. Biometric and eddy covariance-based assessment of decadal carbon sequestration of a temperate Scots pine forest. Agric. For. Meteorol. 2013, 174, 135–143. [Google Scholar] [CrossRef]
  31. Jansen, J.J.; Sevenster, J.A.; Faber, J. Opbrengsttabellen voor Belangrijke Boomsoorten in Nederland; IBN-Rapport 221, Hinkeloord rapport No. 17; Institute for Forestry and Nature Research: Wageningen, The Netherlands, 1996; pp. 1–240. [Google Scholar]
  32. Carrara, A.; Kowalski, A.S.; Neirynck, J.; Janssens, I.A.; Yuste, J.C.; Ceulemans, R. Net ecosystem CO2 exchange of mixed forest in Belgium over 5 years. Agric. For. Meteorol. 2003, 119, 209–227. [Google Scholar] [CrossRef]
  33. Databank Ondergrond Vlaanderen—Vlaamse Overheid, Departement Omgeving, Vlaams Planbureau voor Omgeving (VPO). Digitale bodemkaart van het Vlaams Gewest: Bodemtypes. Available online: https://www.dov.vlaanderen.be (accessed on 5 October 2021).
  34. Integrated Carbon Observation System (ICOS) Research Infrastructure. ICOS Handbook: Knowledge through Observations, 2nd ed.; ICOS ERIC: Helsinki, Finland, 2020. [Google Scholar]
  35. Geopunt. Orthofotomozaïek, Middenschalig, Zomeropnamen. 2018, Vlaanderen. Available online: www.geopunt.be (accessed on 21 October 2021).
  36. Gielen, B.; Acosta, M.; Altimir, N.; Buchmann, N.; Cescatti, A.; Ceschia, E.; Fleck, S.; Hörtnagl, L.; Klumpp, K.; Kolari, P.; et al. Ancillary vegetation measurements at ICOS ecosystem stations. Int. Agrophys. 2018, 32, 645–664. [Google Scholar] [CrossRef]
  37. Mokany, K.; Raison, R.J.; Prokushkin, A.S. Critical analysis of root: Shoot ratios in terrestrial biomes. Glob. Chang. Biol. 2006, 12, 84–96. [Google Scholar] [CrossRef]
  38. Sleutel, S.; D’Hose, T.; Lettens, S.; Ruysschaert, G.; De Vos, B. Monitoring Koolstofgehalte in Vlaamse Bodems in Openbaar Domein en Particuliere Tuinen (Opdracht VPO-OMG_VPO_2018_15-F02)–Eindrapport; Vlaams Planbureau voor Omgeving: Belgium, Brussel, 2020. [Google Scholar]
  39. Berhongaray, G.; King, J.S.; Janssens, I.A.; Ceulemans, R. An optimized fine root sampling methodology balancing accuracy and time investment. Plant Soil 2013, 366, 351–361. [Google Scholar] [CrossRef]
  40. Xiao, C.W.; Yuste, J.C.; Janssens, I.A.; Roskams, P.; Nachtergale, L.; Carrara, A.; Sanchez, B.Y.; Ceulemans, R. Above- and belowground biomass and net primary production in a 73-year-old Scots pine forest. Tree Physiol. 2003, 23, 505–516. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. IBM Corp. IBM SPSS Statistics for Windows; Version 28.0; IBM Corp: Armonk, NY, USA, 2021. [Google Scholar]
  42. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2020. [Google Scholar]
  43. Campioli, M.; Gielen, B.; Göckede, M.; Papale, D.; Bouriaud, O.; Granier, A. Temporal variability of the NPP-GPP ratio at seasonal and interannual time scales in a temperate beech forest. Biogeosciences 2011, 8, 2481–2492. [Google Scholar] [CrossRef] [Green Version]
  44. Moreno-Fernandez, D.; Diaz-Pines, E.; Barbeito, I.; Sanchez-Gonzalez, M.; Montes, F.; Rubio, A.; Canellas, I. Temporal carbon dynamics over the rotation period of two alternative management systems in Mediterranean mountain Scots pine forests. For. Ecol. Manag. 2015, 348, 186–195. [Google Scholar] [CrossRef]
  45. Collalti, A.; Trotta, C.; Keenan, T.F.; Ibrom, A.; Bond-Lamberty, B.; Grote, R.; Vicca, S.; Reyer, C.P.; Migliavacca, M.; Veroustraete, F.; et al. Thinning can reduce losses in carbon use efficiency and carbon stocks in managed forests under warmer climate. J. Adv. Model. Earth Syst. 2018, 10, 2427–2452. [Google Scholar] [CrossRef]
  46. Augusto, L.; De Schrijver, A.; Vesterdal, L.; Smolander, A.; Prescott, C.; Ranger, J. Influences of evergreen gymnosperm and deciduous angiosperm tree species on the functioning of temperate and boreal forests. Biol. Rev. 2015, 90, 444–466. [Google Scholar] [CrossRef]
  47. Pretzsch, H.; del Río, M.; Ammer, C.; Avdagic, A.; Barbeito, I.; Bielak, K.; Brazaitis, G.; Coll, L.; Dirnberger, G.; Drössler, L.; et al. Growth and yield of mixed versus pure stands of Scots pine (Pinus sylvestris L.) and European beech (Fagus sylvatica L.) analysed along a productivity gradient through Europe. Eur. J. For. Res. 2015, 134, 927–947. [Google Scholar] [CrossRef] [Green Version]
  48. Jandl, R.; Lindner, M.; Vesterdal, L.; Bauwens, B.; Baritz, R.; Hagedorn, F.; Johnson, D.W.; Minkkinen, K.; Byrne, K.A. How strongly can forest management influence soil carbon sequestration? Geoderma 2007, 137, 253–268. [Google Scholar] [CrossRef]
  49. Noormets, A.; Epron, D.; Domec, J.C.; McNulty, S.G.; Fox, T.; Sun, G.; King, J.S. Effects of forest management on productivity and carbon sequestration: A review and hypothesis. For. Ecol. Manag. 2015, 355, 124–140. [Google Scholar] [CrossRef] [Green Version]
  50. Ruiz-Peinado, R.; Bravo-Oviedo, A.; López-Senespleda, E.; Montero, G.; Río, M. Do thinnings influence biomass and soil carbon stocks in Mediterranean maritime pinewoods? Eur. J. For. Res. 2013, 132, 253–262. [Google Scholar] [CrossRef]
  51. Zhang, X.; Guan, D.; Li, W.; Sun, D.; Jin, C.; Yuan, F.; Wang, A.; Wu, J. The effects of forest thinning on soil carbon stocks and dynamics: A meta-analysis. For. Ecol. Manag. 2018, 429, 36–43. [Google Scholar] [CrossRef]
  52. Brunner, I.; Godbold, D.L. Tree roots in a changing world. J. For. Res. 2007, 12, 78–82. [Google Scholar] [CrossRef] [Green Version]
  53. Buchholz, T.; Friedland, A.J.; Hornig, C.E.; Keeton, W.S.; Zanchi, G.; Nunery, J. Mineral soil carbon fluxes in forests and implications for carbon balance assessments. Glob. Chang. Biol. Bioenergy 2014, 6, 305–311. [Google Scholar] [CrossRef]
  54. Gruba, P.; Socha, J. Exploring the effects of dominant forest tree species, soil texture, altitude, and pHH2O on soil carbon stocks using generalized additive models. For. Ecol. Manag. 2019, 447, 105–114. [Google Scholar] [CrossRef]
  55. Van Damme, F. Invloed van de Omschakeling van Groenblijvend Naaldblad naar Bladverliezende Breedbladerige Bossen op de Koolstofopslag. Master’s Thesis, Antwerp University, Antwerp, Belgium, 2021. [Google Scholar]
  56. Wardenaar, E.C.; Sevink, J. A comparative study of soil formation in primary stands of Scots pine (planted) and poplar (natural) on calcareous dune sands in the Netherlands. Plant Soil 1992, 140, 109–120. [Google Scholar] [CrossRef]
  57. Øyen, B.H.; Blom, H.H.; Gjerde, I.; Myking, T.; Sætersdal, M.; Thunes, K.H. Ecology, history and silviculture of Scots pine (Pinus sylvestris L.) in western Norway–a literature review. Forestry 2006, 79, 319–329. [Google Scholar] [CrossRef]
  58. Janecka, K.; Harvey, J.E.; Trouillier, M.; Kaczka, R.J.; Metslaid, S.; Metslaid, M.; Buras, A.; Wilmking, M. Higher winter-spring temperature and winter-spring/summer moisture availability increase scots pine growth on coastal dune microsites around the South Baltic Sea. Front. For. Glob. Chang. 2020, 3, 578912. [Google Scholar] [CrossRef]
  59. Schroeder, J.; Klinner, S.; Koerner, M. A New Set of Biomass Functions for Quercus petraea in Western Pomerania. Balt. For. 2017, 23, 449–462. [Google Scholar]
  60. Zianis, D.; Muukkonen, P.; Mäkipää, R.; Mencuccini, M. Biomass and stem volume equations for tree species in Europe. Silva Fenn. Monogr. 2005, 4, 63. [Google Scholar] [CrossRef]
  61. Röhling, S.; Demant, B.; Dunger, K.; Neubauer, M.; Oehmichen, K.; Riedel, T.; Stümer, W. Equations for estimating belowground biomass of Silver Birch, Oak and Scots Pine in Germany. iForest 2019, 12, 166. [Google Scholar] [CrossRef]
  62. Suchomel, C.; Pyttel, P.; Becker, G.; Bauhus, J. Biomass equations for sessile oak (Quercus petraea (Matt.) Liebl.) and hornbeam (Carpinus betulus L.) in aged coppiced forests in southwest Germany. Biomass Bioenergy 2012, 46, 722–730. [Google Scholar] [CrossRef]
  63. Uri, V.; Varik, M.; Aosaar, J.; Kanal, A.; Kukumägi, M.; Lõhmus, K. Biomass production and carbon sequestration in a fertile silver birch (Betula pendula Roth) forest chronosequence. For. Ecol. Manag. 2012, 267, 117–126. [Google Scholar] [CrossRef]
  64. Hughes, M.K. Tree biocontent, net production and litter fall in a deciduous woodland. Oikos 1971, 22, 62–73. [Google Scholar] [CrossRef]
  65. Škėma, M.; Mikšys, V.; Aleinikovas, M.; Šilinskas, B.; Varnagirytė-Kabašinskienė. Biomass Structure and morphometric parameters for non-destructive biomass estimation of common forest underbrush species in Lithuania. Pol. J. Environ. Stud. 2018, 27, 325–333. [Google Scholar] [CrossRef]
  66. Bartelink, H.H. Allometric relationship for biomass and leaf area of beech (Fagus sylvatica L). Ann. For. Sci. 1997, 54, 39–50. [Google Scholar] [CrossRef] [Green Version]
  67. le Goff, N.; Ottorini, J.M. Root biomass and biomass increment in a beech (Fagus sylvatica L.) stand in North-East France. Ann. For. Sci. 2001, 58, 1–13. [Google Scholar] [CrossRef]
  68. Blujdea, V.N.B.; Pilli, R.; Dutca, I.; Ciuvat, L.; Abrudan, I.V. Allometric biomass equations for young broadleaved trees in plantations in Romania. For. Ecol. Manag. 2012, 264, 172–184. [Google Scholar] [CrossRef]
  69. Annighöfer, P.; Ameztegui, A.; Ammer, C.; Balandier, P.; Bartsch, N.; Bolte, A.; Coll, L.; Collet, C.; Ewald, J.; Frischbier, N.; et al. Species-specific and generic biomass equations for seedlings and saplings of European tree species. Eur. J. For. Res. 2016, 135, 313–329. [Google Scholar] [CrossRef] [Green Version]
  70. Oleksyn, J.; Reich, P.B.; Chalupka, W.; Tjoelker, M.G. Differential above-and below-ground biomass accumulation of European Pinus sylvestris populations in a 12-year-old provenance experiment. Scand. J. For. Res. 1999, 14, 7–17. [Google Scholar] [CrossRef]
  71. Marziliano, P.A.; Lafortezza, R.; Medicamento, U.; Lorusso, L.; Giannico, V.; Colangelo, G.; Sanesi, G. Estimating belowground biomass and root/shoot ratio of Phillyrea latifolia L. in the Mediterranean forest landscapes. Ann. For. Sci. 2015, 72, 585–593. [Google Scholar] [CrossRef]
Figure 1. Aerial photos of the study sites. (a): De Inslag with the reference pine stand (dashed line) and the thinned pine stand (solid line); (b): De Mik with the reference pine stand (dashed line) and oak stand planted after the clear cut (solid line). Screenshot from [35].
Figure 1. Aerial photos of the study sites. (a): De Inslag with the reference pine stand (dashed line) and the thinned pine stand (solid line); (b): De Mik with the reference pine stand (dashed line) and oak stand planted after the clear cut (solid line). Screenshot from [35].
Forests 13 01679 g001
Figure 2. Overview of the different measuring plots and the components measured. Ø indicates the diameter.
Figure 2. Overview of the different measuring plots and the components measured. Ø indicates the diameter.
Forests 13 01679 g002
Figure 3. C stocks (tC ha−1) in the aboveground biomass of the overstory of the investigated stands with differentiation between the contribution of Scots pine (dark green) and of the other species (light green).
Figure 3. C stocks (tC ha−1) in the aboveground biomass of the overstory of the investigated stands with differentiation between the contribution of Scots pine (dark green) and of the other species (light green).
Forests 13 01679 g003
Figure 4. Total ecosystem C storage (tC ha−1; mean and SD) for the investigated stands, with a separation between components. Note that (i) the aboveground C consists of the aboveground biomass and the woody debris, and (ii) for each soil layer, the C stock of the mineral soil and the fine and medium roots were aggregated.
Figure 4. Total ecosystem C storage (tC ha−1; mean and SD) for the investigated stands, with a separation between components. Note that (i) the aboveground C consists of the aboveground biomass and the woody debris, and (ii) for each soil layer, the C stock of the mineral soil and the fine and medium roots were aggregated.
Forests 13 01679 g004
Table 1. Stand characteristics in winter 2020–2021.
Table 1. Stand characteristics in winter 2020–2021.
De Inslag Ref.De Inslag Thin.Mik Ref.Mik Oak
Total basal area (m2 ha1)31.026.523.223.9
Dominant speciesP. sylvestrisP. sylvestrisP. sylvestrisQ. robur
Basal area dominant species (m2 ha1)26.820.222.518.9
Stand age (years)91917714
Density (trees ha1)3751842551225
Tree height (m, mean ± SE)17.9 ± 0.720.5 ± 0.718.1 ± 0.512.4 ± 0.2
Diameter breast height (cm, mean ± SE)31.1 ± 1.437.3 ± 1.733.2 ± 1.014.0 ± 0.3
Table 2. Overview of the mean values of C stocks (tC ha−1) with standard error (SE; tC ha−1) at the study stands. ‘%’ indicates the percentage relative to the total C storage in the ecosystem for each component. The p-value shows the statistical significance between pairs of stands (De Inslag Ref. vs. De Inslag Thin. and Mik Ref. vs. Mik Thin.); the star (*) indicates when p < 0.05, the plus (+) indicates when 0.05 < p < 0.10. NA indicates not applicable.
Table 2. Overview of the mean values of C stocks (tC ha−1) with standard error (SE; tC ha−1) at the study stands. ‘%’ indicates the percentage relative to the total C storage in the ecosystem for each component. The p-value shows the statistical significance between pairs of stands (De Inslag Ref. vs. De Inslag Thin. and Mik Ref. vs. Mik Thin.); the star (*) indicates when p < 0.05, the plus (+) indicates when 0.05 < p < 0.10. NA indicates not applicable.
De Inslag Ref.De Inslag Thin. Mik Ref.Mik Oak
Mean ± SE (tC ha1)%Mean ± SE (tC ha1)%pMean ± SE (tC ha1)%Mean ± SE (tC ha1)%p
Aboveground biomass
 Overstory93.3 ± 3.4 3475.3 ± 10.2300.1569.1 ± 1.12264.2 ± 3.6200.27
 Understory
  Grass, moss, and fern0.21 ± 0.05 <10.21 ± 0.1111.001.38 ± 0.12<10.35 ± 0.06<1<0.001 *
  Young trees and shrubs a4.1 ± 2.024.6 ± 0.720.810000NA
Belowground biomass
 Coarse roots of overstory10.8 ± 0.348.8 ± 1.040.107.8 ± 0.2210.9 ± 0.43<0.001
 Coarse roots of young trees and shrubs 0.88 ± 0.4<11.1 ± 0.2<10.620000NA
 Fine and medium roots
  0–10 cm2.3 ± 0.411.8 ± 0.310.371.5 ± 0.3<11.8 ± 0.210.05 +
  10–30 cm2.1 ± 0.612.2 ± 0.610.792.1 ± 0.413.2 ± 0.710.24
  30–60 cm3.0 ± 0.910.71 ± 0.23<10.040 *1.0 ± 0.1<11.9 ± 0.410.25
Forest floor41.9 ± 2.71544.7 ± 3.7180.5446.2 ± 2.11521.7 ± 1.07<0.001 *
Woody debris
 Fine woody debris2.5 ± 1.212.6 ± 1.210.570.80 ± 0.16<11.36 ± 0.30<10.08 +
 Coarse woody debris2.1 ± 1.713.0 ± 2.810.893.2 ± 1.415.9 ± 1.720.27
Mineral soil
 0–10 cm44.9 ± 4.81752.7 ± 7.3210.4750.8 ± 3.51663.8 ± 4.8200.036 *
 10–30 cm39.7 ± 2.11529.7 ± 2.212<0.001 *68.1 ± 5.82180.3 ± 3.6250.08 +
 30–60 cm22.9 ± 3.1822.1 ± 3.890.8865.8 ± 7.72160.4 ± 25.6190.85
a refers to individuals with a diameter between 1 and 5 cm at 5 cm aboveground (young trees and shrubs with a diameter < 1 cm were not reported as with a negligible stock of < 0.1 tC ha−1).
Table 3. Mean values of C stocks (tC ha−1) with standard error (SE; tC ha−1) of the fine (with less than 2 mm in diameter) and medium roots (with diameter between 2 and 5 mm) in different soil layers at the study stands. The star (*) indicates a statistically significant difference between stands at p < 0.05, the plus (+) indicates 0.05 < p < 0.10.
Table 3. Mean values of C stocks (tC ha−1) with standard error (SE; tC ha−1) of the fine (with less than 2 mm in diameter) and medium roots (with diameter between 2 and 5 mm) in different soil layers at the study stands. The star (*) indicates a statistically significant difference between stands at p < 0.05, the plus (+) indicates 0.05 < p < 0.10.
De Inslag Ref.De Inslag Thin. Mik Ref.Mik Oak
Mean ± SE
(tC ha1)
Mean ± SE
(tC ha1)
pMean ± SE
(tC ha1)
Mean ± SE
(tC ha1)
p
0–10 cm
Ø < 2 mm1.12 ± 0.131.08 ± 0.100.820.78 ± 0.071.26 ± 0.09<0.001 *
2 < Ø < 5 mm1.18 ± 0.300.67 ± 0.310.06 +0.69 ± 0.300.52 ± 0.120.57
10–30 cm
Ø < 2 mm0.71 ± 0.080.66 ± 0.110.841.27 ± 0.201.38 ± 0.220.61
2 < Ø < 5 mm1.40 ± 0.561.54 ± 0.510.830.85 ± 0.341.86 ± 0.530.18
30–60 cm
Ø < 2 mm1.91 ± 0.750.53 ± 0.150.029 *0.70 ± 0.271.28 ± 0.320.15
2 < Ø < 5 mm1.12 ± 0.470.18 ± 0.140.110.31 ± 0.270.60 ± 0.170.25
Table 4. Characteristics of the mineral soil for each soil layer in each investigated stand. The star (*) indicates the statistically significant difference between stands at p < 0.05, the plus (+) indicates 0.05 < p < 0.10.
Table 4. Characteristics of the mineral soil for each soil layer in each investigated stand. The star (*) indicates the statistically significant difference between stands at p < 0.05, the plus (+) indicates 0.05 < p < 0.10.
De Inslag Ref.De Inslag Thin. Mik Ref.Mik Oak
Mean ± SEMean ± SEpMean ± SEMean ± SE p
0–10 cm
bulk density (g cm−3)1.18 ± 0.031.23 ± 0.020.191.20 ± 0.041.23 ± 0.040.37
C content (%)3.80 ± 0.384.36 ± 0.650.764.30 ± 0.285.26 ± 0.390.05 +
10–30 cm
bulk density (g cm−3)1.28 ± 0.011.22 ± 0.030.07 +1.23 ± 0.051.21 ± 0.030.60
C content (%)1.55 ± 0.091.22 ± 0.090.004 *2.82 ± 0.233.39 ± 0.190.06 +
30–60 cm
bulk density (g cm−3)1.29 ± 0.021.17 ± 0.050.781.45 ± 0.041.33 ± 0.040.07 +
C content (%)0.59 ± 0.080.64 ± 0.130.09 +1.53 ± 0.221.51 ± 0.620.97
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Van Damme, F.; Mertens, H.; Heinecke, T.; Lefevre, L.; De Meulder, T.; Portillo-Estrada, M.; Roland, M.; Gielen, B.; Janssens, I.A.; Verheyen, K.; et al. The Impact of Thinning and Clear Cut on the Ecosystem Carbon Storage of Scots Pine Stands under Maritime Influence in Flanders, Belgium. Forests 2022, 13, 1679. https://doi.org/10.3390/f13101679

AMA Style

Van Damme F, Mertens H, Heinecke T, Lefevre L, De Meulder T, Portillo-Estrada M, Roland M, Gielen B, Janssens IA, Verheyen K, et al. The Impact of Thinning and Clear Cut on the Ecosystem Carbon Storage of Scots Pine Stands under Maritime Influence in Flanders, Belgium. Forests. 2022; 13(10):1679. https://doi.org/10.3390/f13101679

Chicago/Turabian Style

Van Damme, Freke, Hana Mertens, Thilo Heinecke, Lodewijk Lefevre, Tim De Meulder, Miguel Portillo-Estrada, Marilyn Roland, Bert Gielen, Ivan A. Janssens, Kris Verheyen, and et al. 2022. "The Impact of Thinning and Clear Cut on the Ecosystem Carbon Storage of Scots Pine Stands under Maritime Influence in Flanders, Belgium" Forests 13, no. 10: 1679. https://doi.org/10.3390/f13101679

APA Style

Van Damme, F., Mertens, H., Heinecke, T., Lefevre, L., De Meulder, T., Portillo-Estrada, M., Roland, M., Gielen, B., Janssens, I. A., Verheyen, K., & Campioli, M. (2022). The Impact of Thinning and Clear Cut on the Ecosystem Carbon Storage of Scots Pine Stands under Maritime Influence in Flanders, Belgium. Forests, 13(10), 1679. https://doi.org/10.3390/f13101679

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop