Soil carbon and nutrient stocks under Scots pine plantations in comparison to European beech forests: a paired-plot study across forests with different management history and precipitation regimes

Background: Organic carbon stored in forest soils (SOC) represents an important element of the global C cycle. It is thought that the C storage capacity of the stable pool can be enhanced by increasing forest productivity, but empirical evidence in support of this assumption from forests differing in tree species and productivity, while stocking on similar substrate, is scarce. Methods: We determined the stocks of SOC and macro-nutrients (nitrogen, phosphorus, calcium, potassium and magnesium) in nine paired European beech/Scots pine stands on similar Pleistocene sandy substrates across a precipitation gradient (560–820 mm∙yr ) in northern Germany and explored the influence of tree species, forest history, climate, and soil pH on SOC and nutrient pools. Results: While the organic layer stored on average about 80% more C under pine than beech, the pools of SOC and total N in the total profile (organic layer plus mineral soil measured to 60 cm and extrapolated to 100 cm) were greater under pine by about 40% and 20%, respectively. This contrasts with a higher annual production of foliar litter and a much higher fine root biomass in beech stands, indicating that soil C sequestration is unrelated to the production of leaf litter and fine roots in these stands on Pleistocene sandy soils. The pools of available P and basic cations tended to be higher under beech. Neither precipitation nor temperature influenced the SOC pool, whereas tree species was a key driver. An extended data set (which included additional pine stands established more recently on former agricultural soil) revealed that, besides tree species identity, forest continuity is an important factor determining the SOC and nutrient pools of these stands. Conclusion: We conclude that tree species identity can exert a considerable influence on the stocks of SOC and macronutrients, which may be unrelated to productivity but closely linked to species-specific forest management histories, thus masking weaker climate and soil chemistry effects on pool sizes.


Background
The sequestration of atmospheric carbon dioxide in soils as stable organic matter is discussed as a means of mitigating anthropogenic climate warming (Spalding et al. 2012). Forest soils store one third (c. 720 Pg) of the global organic carbon (C) (Schlesinger 1991, Price et al. 2012) and thus are an important element of the global C cycle (Lal 2005). Soil organic C (SOC) storage depends not only on climate (Lorenz & Lal 2010) and mineral properties, notably the potential of soils to stabilize organic C on reactive mineral surfaces (Vogel et al. 2015), but it may also be influenced by tree species identity and forest management (Binkley & Giardina 1998;Jandl et al. 2007;Grüneberg et al. 2019). Tree species influence the C input to the soil by differences in the amounts of aboveground litter (leaves, fruits, coarse woody debris) and root litter (Roy et al. 2001;Finér et al. 2011;Schlesinger & Bernhardt 2013) as well as in root exudation rates and carbohydrate transfer to the mycorrhiza (Godbold et al. 2006;Pausch & Kuzyakov 2018;Akatsuki & Makita 2020). Litter chemistry, notably the content of lignin, nitrogen (N) and calcium (Ca) (and possibly exudate chemistry) are additional important factors that affect C turnover in the soil (Berg & McClaugherty 2003;Reich et al. 2005).
Tree species differences in litter properties should be particularly large between angiosperm and gymnosperm trees, as the latter produce needles that are generally poorer in N and Ca than leaves of broad-leaved trees (Berg & McClaugherty 2003;Hobbie et al. 2006), and conifers are often more shallow-rooted than many angiosperm trees (Polomski & Kuhn 1998), with possible consequences for root-borne C input. Species differences in the amount and seasonality of litter production may increase the dissimilarity between these groups. While it is generally agreed that tree species often have a profound effect on the C storage in the organic layer on the forest floor (e.g. Binkley & Valentine 1991;Vesterdal & Raulund-Rasmussen 1998;Fischer et al. 2002), the species effect on mineral soil SOC pools is less clear (e.g. Heinsdorf 2002;Ladegaard-Pedersen et al. 2005;Jandl et al. 2007;Schulp et al. 2008;Grüneberg et al. 2019). Thus, which forest type in general sequesters more SOC cannot be simply be stated, because C sequestration also depends on soil and climate. Hence, with respect to current efforts of quantifying C stock changes in forests and their development under altered climatic conditions, a better understanding of how the forest type and the forest composition influence C sequestration is needed. In the debate about how temperate production forests of the future should look like, this understanding is particularly needed for comparing the impact of conifers, which are widely planted in monocultures and valued in forestry for their timber productivity and economic value, and native hardwood trees, which may have advantages in terms of biodiversity conservation, lower soil acidification potential and lower susceptibility to insect attack (Ammer et al. 2008;Gessler et al. 2006).
In the past, Scots pine (Pinus sylvestris L.) plantations have widely replaced the natural broadleaf forests in the lowlands of northern Central Europe, i.e. in the north of the Netherlands, Germany and Poland. At least in Germany and the Netherlands, most of the pine stands have been planted during the last 250 years on soils that once were stocked with European beech (Fagus sylvatica L.) and oak forests (Quercus petraea Matt. Liebl and Q. robur L.) (Kremser 1990;Milnik 2007). Pine-dominated stands nowadays cover 52.1% of the current forest area in the northern German lowlands (deviating from the region defined as North German Plain Fig. 1), while the natural occurrence of Scots pine is only estimated at 2.1% of the area (Förster et al., 2021). In absolute terms, a rough calculation of the current pine forest area gives for Poland c. 6.5 Mio ha (1990), for Germany about 2.4 Mio ha (2012), and for the Netherlands c. 0.2 Mio ha (Kuusela 1994;Leuschner et al. 2013;BMEL 2020). These numbers demonstrate the enormous extension of pine plantations in this region. However, while Scots pine is valued as a fast-growing, non-demanding timber species, it has several disadvantages, notably a high susceptibility to fire and pathogen attack, the build-up of a relatively dry and infertile organic layer, and the tendency to acidify the soil (Bublinec 1974;Reich et al. 2005;Prietzel et al. 2006). In terms of C sequestration in the soil, it is not fully clear, how such pine plantations differ from the naturally prevailing beech forests in this area.
Several regional studies that compared Scots pine and the natural broad-leaved forests (beech) on similar soil found in most cases larger organic layer SOC stocks under pine (Fischer et al. 2002;Leuschner et al. 2013), while the information for the mineral soil is contradictory (e.g. Heinsdorf 2002;Jandl et al. 2007;Schulp et al. 2008;Leuschner et al. 2013). The nation-wide Forest Soil Inventory (NFSI) of Germany, which analyzed data from 1900 soil pits, gave roughly three times lager C stores in the organic layer for pine stands (mostly Scots pine) than for beech stands, while the mineral soil SOC pool was significantly smaller than under beech (Grüneberg et al. 2019). Regional studies from northern Germany also found larger C pools in the organic layer and lower pools in the mineral soil under pine compared to beech (e.g. Wördehoff et al. 2017aWördehoff et al. , 2017b. However, the tree species effect is in these large-scale inventories overlaid by variation in geological substrate, topography, and soil moisture regimes. Moreover, stand history, which can strongly influence SOC pools (Compton and Boone 2000), was not considered in these inventories. Many pine plantations in the lowlands of northern Germany have been established on soils that were heathlands or arable fields for long periods, which may imprint on soil C storage (von Oheimb et al. 2008;Laganière et al. 2010;Leuschner et al. 2013). The comparison of paired pine and beech stands on similar substrate and on soils with long forest continuity can minimize these confounding factors. This approach has been adopted on a local scale (e.g. in the Veluwe/the Netherlands and Brandenburg/ Germany, Heinsdorf 2002;Schulp et al. 2008;Leuschner et al. 2013), but not yet on a supra-regional level.
Here, we present the results of a country-wide study on the soil C and nutrient stocks (organic layer plus mineral soil measured to 60 cm and extrapolated to 100 cm) in paired beech and pine stands along a climate gradient (temperate-oceanic to temperate sub-continental) in the Pleistocene lowlands of northern Germany, spanning from the North Sea coast to the Polish border. As many pine stands stock on soils that were arable land or heathland in the more distant past, we also address the importance of stand history on SOC stocks.
The goal of the study was to examine the influence of tree species (beech vs. pine), stand history (ancient vs. recent forest), and climate (560-820 mm•yr − 1 precipitation) on (i) the SOC pool and SOC distribution in the profile, and (ii) the soil pools of the macro-nutrients N, P, Ca, Mg and K. To reduce the influence of soil chemical factors, the stands were selected on fairly uniform geological substrate, i.e. base-poor unconsolidated glacial deposits which dominate in the region. The influence of forest continuity ('ancient forest sites': > 230 years vs. 'recent forest sites': c. 50-130 years) was addressed by comparing our results (mostly ancient forest sites) to data from 7 recent pine forest sites provided by the study of Leuschner et al. (2013) in the study region.

Study sites
The study was conducted in the lowlands of northern Germany in nine pairs of European beech and Scots pine stands on glacial unconsolidated sandy to loamy sediments from the last (Weichselian) or penultimate Ice Age (Saalian). The nine study sites were located between Bremen in the west with a cool-temperate oceanic climate to the river Oder in the east with a cooltemperate sub-continental climate (Fig. 1). Mean annual temperature varied between 9.0 and 9.6°C, mean annual precipitation decreased from 712 to 823 mm•yr − 1 in the west to 564-657 mm•yr − 1 in the east ( Table 1). The air temperature (2 m) and precipitation data for the nine sites were derived from the 1 km × 1 km grid data of temperature and precipitation interpolated by Deutscher Wetterdient (DWD) from its dense climate station net (DWD Climate Data Center, v1.0). The geological substrates comprise sandy to loamy moraine deposits, glacio-fluviatile sandy deposits and sandy substrates covered by aeolian sandy loess. All substrates led to the formation of acidic, nutrient-poor soils without Fig. 1 Location of the study sites. Location of the nine study sites in Germany (inlet) and location of the 18 stands in the North German Plain (grey-shaded, Be = near-natural beech forests, Pi = pine forests). The white triangles mark the location of the additionally analyzed seven recent pine stands in western Brandenburg adopted from Leuschner et al. (2013) groundwater influence. The soil types were spododystric Cambisols and Podzols.
Information given in Glaser and Hauke (2004) and consultation of local sources (Meyer et al. 2006) indicated that most of the forest stands (all beech forests and 6 out of 9 pine stands) grew on sites with long forest continuity, i.e. uninterrupted forest cover during at least the last 200 years. It is likely that these 15 sites were also covered by forest during the Middle Ages and probably have never been cleared for agricultural use. The remaining three pine stands were either heathland or farmland 200 years ago and thus have a relatively short continuity. The paired beech and pine stands were selected in short distance to each other on similar geological substrate. In most pairs, the distance between the beech and pine stands was less than 1 km. Only in case of site 1 (Weichel), the distance was greater (20 km).

Study design and soil sampling
In each stand, ten circular plots were established at random positions. The plots had a size of 0.1 ha (35.68 m in diameter), and the minimum distance between the plot centers was 66 m. Thus, at least 30 m separated the circle edges from each other. To exclude edge effects, the distance to forest paths and forest edges was at least 20 m.
One soil profile each per plot was dug at a random location to 60 cm depth and soil samples were taken at four depths (0-10, 10-20, 20-40 and 40-60 cm) and in the organic layer (L, Of and Oh horizon according to Arbeitskreis Standortkartierung 2016). Three samples per depth were collected at the three walls of the pit and combined to one composite sample for each depth and profile. We chose a sampling depth to 60 cm according to fine root studies on sandy soils by Meier et al. (2018) in the North German Lowlands and investigations by Anders et al. (2004) in a forest area in the north-east German lowlands, which show that beech and pine have the largest amount of root mass in the uppermost 60 cm of the mineral soil. Thus, a very large proportion of the rooted soil volume was covered by our soil chemical sampling scheme. Soil was sampled with a metal cylinder of 100 cm 3 volume. The organic layer was sampled with a root corer (5.5 cm in diameter) at three locations close to the soil profile; the samples were combined to a single sample as well. Corresponding to the number of circular sampling plots, the number of analytical replicates was 10 at the stand level (one sample per depth per plot) and thus 900 in total (2 species × 9 stands × 5 depths × 10 stand-level replicate plots). The soil samples were stored in plastic bags and kept at 4°C in the laboratory. Each study site is presented by a pair of beech (Be) and pine stands (Pi). Given are stand age, forest continuity, mean annual precipitation (MAP), mean annual temperature (MAT), stand basal area, stem density, and stand area. Forest stands were classified as "ancient" for > 200 years of forest continuity, and "recent" for < 150 years of continuity

Soil analysis
Prior to analysis, the mineral soil was sieved through a 2-mm sieve and roots were picked out by hand from the mineral soil and organic layer. The pH was measured in a suspension of 10 g fresh sieved soil (or 2.5 g of organic layer material) in 25 mL deionized water (pH(H 2 O)) or in 1 M KCl (pH (KCl)). Total C and total N were measured in ground soil that was dried at 70°C for 48 h. C and N concentrations were determined by gas chromatography (vario EL III, elementar, Hanau, Germany). While all C in the organic layer was assumed to be organic C, inorganic C was determined in the mineral soil as the difference between gas chromatographic C analyses before and after combustion at 640°C. To obtain organic C content, inorganic C content (determined as the C content remaining after combustion) was subsequently subtracted from total C content (determined as the C content before combustion). The total content of K, Mg and Ca in the organic layer was determined by ICP-OES analysis (Perkin Elmer Optima 5300 DV). In the mineral soil, only the salt-exchangeable base cations were determined. 2.5 g of fresh soil was percolated with 0.2 M BaCl 2 solution and the concentration of K, Ca, Mg, Al, Fe and Mn determined in the solution with ICP-OES. Total phosphorus (P) concentration was measured in the organic layer samples with ICP-OES, while in the mineral soil, resin-extractable P was determined (resinbag method). 1 g of fresh soil was suspended in 30 mL of water and P was extracted with the anion exchanger resin Dowex 1 × 8-50. P was then re-exchanged from the resin with NaCl and NaOH solutions, and the P concentration determined in a photometer at 712 nm against water using the colorimetric molybdate-ascorbic acid method. To obtain element pools for the mineral soil depth levels, the concentration data were converted to volumetric data using the bulk soil density of the depth levels, that were determined separately in all profiles and depths. The organic layer data were multiplied with organic layer mass per area to obtain element stocks per ground area.

Statistical analysis
The data analysis was conducted with R, version 3.6.3 (R Core Team 2020) in the framework of the tidyverse package (Wickham et al. 2019). Stand-level averages were calculated from the 10 samples per stand. The means (and standard errors) per tree species used in the graphs are derived from the mean values calculated at the stand level for each of the nine pine and beech stands. In order to obtain element pools for a standard 1-m soil profile, we fitted a Michaelis-Menten curve to the element pools of the 0-10, 10-20, 20-40, and 40-60 cm soil depth levels and extrapolated the pools down to a depth of 1 m based on the resulting concentration-depth curve fitted for each soil plot. We applied ANOVA (F-test) over stand-level averages of all study sites and all soil depth levels to test for significant differences between stand types (beech vs. pine). Here, we also accounted for study site (stand) and soil depth level as additional model variables and considered different trends across soil depth levels (interaction forest type x soil depth). Significantly different means were identified post-hoc with Tukey's HSD test. As recent pine stands, established on sites with former agricultural or heathland use, are more widespread in the study region than ancient pine stands with longer forest continuity, we extended our stand sample by including data from 7 recent pine forests on similar soil in western Brandenburg (study of Leuschner et al. 2013), and tested for the combined effect of stand type and forest continuity (ancient beech stands vs. recent pine stands or ancient pine stands) in a similar testing procedure (recent beech forests were not available). Significantly different means were identified post-hoc with Tukey's HSD test. Further, the influence of important climatic and edaphic variables (temperature, precipitation, pH of the subsoil) on the soil C, N, P and base cation pools in the 0-100 cm profile (C and N: organic layer, mineral soil, and total profile; base cations: only mineral soil considered in the analysis) was tested in linear models, that of management factors (stand age, stand type and history) in multiple-regression analyses. In case of the climatic variables, we considered 30-year averages of annual and summer temperature means and precipitation sums. Before we built the maximal models, we tested the explanatory variables for autocorrelation using the R package cormorant v. 0.0.0.9007 (Link 2020), which revealed a close correlation of summer temperature and precipitation to annual temperature and precipitation, respectively ( Fig. A5 in the Appendix). Subsequently, we applied Spearman's rank correlation analysis (R package Hmisc v.4.4-1, Harrell Jr, 2020) to investigate direct univariate relationships between mineral and organic layer C, N, P and base cation pools with stand age, temperature and precipitation variables, and pH of the subsoil (Table 2). Accordingly, annual temperatures explained the trends in SOC and nutrient pools across the study sites better than summer temperatures, while summer precipitation had a similar or slightly larger influence than annual precipitation. Consequently, we excluded summer temperature and annual precipitation from the following multiple-regression analyses, also in view of the high autocorrelation of the variables related to annual and summer climatic conditions (Fig. A5 in the Appendix). This variable selection matches the ecology of the tree species, as growth is most closely related to summer precipitation, while growing season length may depend more on spring and autumn temperature, which is better reflected in mean annual temperature. Thus, starting from maximal models that contained stand age, pH (subsoil), mean summer precipitation (MSP) and mean annual temperature in addition to stand type and forest continuity as possible explanatory variables of the SOC and nutrient pools, the variables in each model were backward-selected by iterative reduction of the least significant variable until only significant variables remained.

Results
Influence of tree species on the soil C and nutrient pools The organic layers tended to be thicker under pine stands by about 6 mm and stored about 80% more C than the organic layer under beech stands (52.0 vs. 28.6 Mg C•ha − 1 ; difference significant; Fig. 2a). In particular, the slightly decomposed Of layer was markedly thicker under pine than beech, while the fresh litter layer (L) was thinner (Fig. A1 in the Appendix). The SOC pool in the mineral soil to 100 cm was about 27% higher in the pine stands (89.1 Mg C•ha − 1 ) than in the beech stands (70.4 Mg C•ha − 1 ). The total SOC pool (organic layer and mineral soil) was on average about 40% greater under pine than beech (141.2 vs. 99.1 Mg C•ha − 1 ; difference significant; Fig. 2a). Yet, variation between the nine investigated sites was considerable (Fig. 3) and in two of the 9 pairs, the SOC pool was larger in the beech than the pine stand (sites 4 and 9). The soil N content was closely related to the SOC content but the difference between beech and pine stands was smaller. The pine stands stored about 50% more N in the organic layer than the beech stands (2007 vs. 1330 kg N•ha − 1 ; Fig. 2b), but the difference was not significant. While the mineral topsoil under pine was slightly richer in N than the equivalent layers under beech, the subsoil tended to have slightly smaller N pools, resulting in similar mineral soil N pools (3897 vs. 3638 kg N•ha − 1 ). The total soil profile stored about 20% more N under pine than under beech (5903 vs. 4968 kg N•ha − 1 ; difference significant; Fig. 2b). The C:N ratio was by about 3 g g − 1 higher under pine in the entire soil profile (Fig. 4c).
With 87 kg P•ha − 1 , the pools of total phosphorus in the organic layer were about 25% larger under pine than under beech (70 kg P•ha − 1 , Fig. 5a; difference not significant). In the mineral soil to 100 cm, the pool of plantavailable P (resin-P) was slightly (but not significantly) higher under beech than pine (62 vs. 58 kg P resin •ha − 1 ; Fig. 5b). Due to the higher C org content, the C:P ratio in the mineral topsoil (> 30 cm depth) was somewhat higher under pine than beech (Fig. A2b in the Appendix), while the N:P ratio was similar among the forest types (Fig. 4d).
In contrast to N, we found somewhat larger pools of basic cations (Ca, K, Mg) in the organic layer (total pool) and in the mineral soil (BaCl 2 -exchangeable pools) under beech than pine (17.6 vs. 13.9 kmol c •ha − 1 to 100 cm, Fig. 6a and b), but the differences were not significant. Cation exchange capacity (CEC) and base saturation were generally low in the subsoil (< 20 μmol c •g − 1 and < 10%), with a somewhat higher CEC at 50 cm and a slightly higher base saturation in the topsoil under beech (Figs. A3a and A3b in the Appendix). Acidity was very similar in the mineral soil of beech and pine stands (Fig. A4 in the Appendix). Only in the subsoil (50 cm depth), pH (KCl) was slightly (not significantly) lower under beech, while it was much lower (by 0.5 units) under Analysis of the dependence of C org and N pools of the organic layers (OL), and C org , N, P and base cation pools (BC) of the mineral soil (MS) (0-100 cm), on stand age (age), mean annual temperature (MAT), mean summer temperature (MST), mean annual precipitation (MAP), mean summer precipitation (MSP), and subsoil pH (KCl) (40-60 cm soil depth) in the beech and pine stands. Given are the ρ values of the relationships. Significant values are printed in bold. Climatological means refer to the last 30 years Fig. 3 Mean stocks of soil organic carbon in the mineral soil and organic layer. SOC stocks in the mineral soil (0-100 cm, data for 60-100 cm were computed by extrapolation) and the organic layers of the each nine beech (black bars) and pine stands (white bars) at the sites 1 to 9 (given are the means of 10 profiles per stand). Significantly different pool sizes between the stand types are marked with an asterisk (paired t-test) Fig. 2 Soil pool of organic carbon and total nitrogen. Pool of organic C (a) and total N (b) in soil profiles to 100 cm depth (data for 60-100 cm were computed by extrapolation) under beech (black bars) and pine stands (white bars). Total profile refers to mineral soil plus organic layer. Significantly different pool sizes between the stand types are marked with an asterisk (ANOVA). Given are means and standard error of each nine stands, in which each 10 profiles were investigated pine in the organic layer (3.1 vs. 3.6; difference significant).

Differences in SOC and soil nutrient pools between ancient and recent forest stands
The analysis of an extended forest stand sample completed with data from recent pine forests (Leuschner et al. 2013) showed that ancient pine stands had a total SOC stock (organic layer and mineral soil to 100 cm) of 148.9 Mg•ha − 1 , which is about 27% higher than in recent pine stands (116.8 Mg•ha − 1 ) and 50% higher than in ancient beech stands (99.1 Mg•ha − 1 ) (Fig. 7a). While the difference in the mineral soil was not significant, ancient pine stands stored almost twice as much C in the organic layer than recent pine stands (60.5 vs. 33.9 Mg C•ha − 1 , Fig. 7c). Total N pools were also higher in the Fig. 4 Soil concentrations of organic carbon, total nitrogen, C:N ratio, and N:P resin ratio. Mean concentration (and standard error) of (a) organic C, (b) total N, (c) C:N ratio and (d) N:P resin ratio in mineral soil profiles under beech (solid line) and pine stands (dotted line). Different capital letters mark significant differences between stand types for the whole profile (ANOVA). Given are means and standard error of each nine stands, in which each 10 profiles were investigated. Depth-wise differences were not significant Significantly different pool sizes between the stand types are marked with an asterisk (ANOVA). Given are means and standard error of each nine stands, in which each 10 profiles were investigated organic layer of ancient than recent pine stands, while the mineral soil N pool tended to be larger in recent pine stands ( Fig. 7e and f). Consequently, the profile total N pools differed not significantly among the stand types, but were in tendency larger in recent and ancient pine than beech stands. The resin-P pool in the mineral soil was significantly larger in recent than in ancient pine stands and higher than in ancient beech stands (Fig. 8a). For the mineral soil base cation pool (BaCl 2 -exchangeable fraction), a tendency for lower stocks under ancient pine stands was found, but the difference to recent pine and beech stands was not significant (Fig. 8b).
Influence of climatic and edaphic factors, and stand structure combined with stand history on SOC and soil nutrient pools The univariate relationships of SOC with the covariates stand age, subsoil pH and temperature and precipitation differed between beech and pine stands ( Table 2). The organic layer C pool in beech stands tended to increase with mean annual temperature (MAT) and to decrease with stand age, but both relationships were not significant (Table 2). In contrast, the organic layer C pool in pine forests significantly increased with a pH decrease, while the influence of climate and stand age was very weak ( Table 2). The mineral soil SOC pool was not related to any of the tested climatic, edaphic or stand structural (species, forest history, age) variables, neither in beech nor in pine stands (Tables 2 and 3). Modeling pool size variation with a SOC model that accounts for both the environmental and stand structural (species) factors indicated in case of the total profile a significant effect only for stand type (beech vs. pine) and stand history, in case of the organic layer only for stand age, stand type and stand history, while climate and soil pH were insignificant factors Fig. 7 Soil pools of C and N of different forest history stages. Mean pools of organic C in the soil under ancient beech (BeAnc), ancient pine (PiAnc) and recent pine stands (PiRec) for the total profile to 100 cm (including organic layer, data for 60-100 cm were computed by extrapolation) (a), the mineral soil (0-100 cm) (b), and the organic layer (c), and mean pools of total N in the three forest categories for the total profile to 100 cm (including organic layer) (d), the mineral soil (0-100 cm) (e), and the organic layer (f). Given is the population mean (black dot) and the standard error, and the individual means of the stands (open circles). Significantly different population means between the stand types are marked with different small letters. The stands comprise those analyzed in this study and seven additional recent pine stands from western Brandenburg according to Leuschner et al. (2013) ( Table 3), despite a soil acidity increase with stand age (Fig. A5 in the Appendix). The N t pool in the organic layer significantly increased with a pH decrease in the pine stands; in the beech stands, only an insignificant trend toward a N pool increase with MAT was visible (Table 2). In the mineral soil, the N t pool increased with mean annual precipitation in the beech, but not in the pine stands. However, in a similar manner as was observed for the C pools, these direct effects of climatic and edaphic factors disappeared, when the influence of stand type and history was included in linear models. Stand type and history remained as the only significant explanatory variable for the organic layer N pool (Table 3). None of the tested variables had a significant influence on the mineral soil C and N pools in the linear models. Our results thus suggest that due the close linkage of N accumulation to C sequestration in the soil, tree species and stand history effects are overriding effects of temperature, precipitation and also atmospheric N deposition patterns in our study region.
The resin-P pool of the mineral soil of both beech and pine stands decreased towards drier sites (significant only in pine), and significantly decreased with increasing soil acidity in pine (Table 2). In the linear models with all factors included, these effects where of minor Fig. 8 Soil pools of P and base cations of different forest history stages. a Mean pools of available P (resin P) and (b) BaCl 2 -exchangeable base cations in the mineral soil (0-100 cm, data for 60-100 cm were computed by extrapolation) under ancient beech (BeAnc), ancient pine (PiAnc) and recent pine stands (PiRec). Given is the population mean (black dot) and the standard error, and the individual means of the stands (open circles). Significantly different population means between the stand types are marked with different small letters. The stands comprise those analyzed in this study and seven additional recent pine stands from western Brandenburg according to Leuschner et al. (2013) Variables influencing the pools of C org and N t (total soil profile, mineral soil 0-100 cm (MS) and organic layer (OL)) and resin-P (mineral soil 0-100 cm) and exchangeable base cations (BC, mineral soil 0-100 cm) in the beech and pine stands according to linear models. The initial models included mean annual temperature (MAT, last 30 years), mean summer precipitation (MSP, last 30 years), stand age, soil pH at 40-60 cm and stand type in combination with forest history (ancient beech, ancient pine and recent pine stands); insignificant variables (−) were subsequently removed until only significant explanatory variables (×) remained importance, and only summer precipitation and stand type and history remained as the principal explanatory variables for the resin-P pool (Table 3). The exchangeable pool of basic cations in the mineral soil in the beech stands showed a strong and significant relationship to mean annual temperature (MAT) and an insignificant influence of mean summer precipitation (MSP), with an increasing pool of basic cations towards warmer and moister sites (Table 2). However, MAT remained as the only significant explanatory variable, when exploring the combined effect of all explanatory variables in our model analysis (Table 3).

Discussion
Tree species effects on the SOC pool Our paired-plot study revealed under pine an on average about 40% greater SOC pool than under beech (organic layer and mineral soil combined). Only two of the nine investigated plot pairs showed a higher C content under beech. Our results are generally in line with the conclusions of Boča et al. (2014) from a meta-analysis, who found conifers to have on average 26% and 19% more SOC in the forest floor and mineral soil, respectively, compared to beech-dominated stands. It has been discussed, whether there is a difference in the vertical allocation of SOC stocks between conifers and broadleaf trees rather than a difference in total SOC stocks. Various authors suggested that less C in the mineral soil offsets higher C stocks in the forest floor and vice versa (Vesterdal et al. 2008;Prescott and Vesterdal 2013). This is not the case in our study, where SOC stocks were higher under pine in both the organic layer and mineral soil. Other studies directly comparing pine and beech forests on similar soil in the Pleistocene lowlands of northern Central Europe came to partly contradicting results: Fischer et al. (2002) reported 30% higher total SOC pools under pine compared to beech in eastern Brandenburg (north-eastern Germany), Heinsdorf (2002) and Schulp et al. (2008) found similar pools under pine and beech in Brandenburg and the Veluwe (eastern Netherlands), and Leuschner et al. (2013) measured 25% higher pools under beech than pine in western Brandenburg. It appears that the effect of tree species on SOC pools varies from one region to another. This is underpinned by the great variation in SOC pools among the nine stand pairs of our study (beech: 60.8-160.9 Mg C•ha − 1 , pine 89.2-191.7 Mg C•ha − 1 , 0-100 cm). Nevertheless, stand characteristics, i.e. tree species and stand history, were the most important predictors of soil C, N and P pools, while climatic and edaphic factors had only a weak or no influence in the global models. Climate was an additional explanatory factor only for the mineral soil P-resin pool, and it was the only factor in case of the base cation pool.
Lower SOC pools under beech than pine match the observation of Prietzel (2004) that planting beech in pine stands reduced C storage in the organic layer over 16 to 83 years and resulted in a 5% smaller total SOC pool. Beech and pine forests differ in many aspects, which are relevant for SOC accumulation. First, mature beech forests produce substantially more aboveground litter than pine plantations. Litter bucket measurements of leaf litter production in the stands of this study gave for the years 2015 and 2016 by 25 to 48% higher annual litter amounts in the beech stands than the pine stands (Förster et al. 2021). Moreover, fine root biomass is much higher in beech than pine stands, indicating a higher C transfer from the beech root system to the SOC pool (A. Förster, pers. comm.). Anders et al. (2004) reported a total root biomass of 17.9 Mg•ha − 1 in beech stands compared to 10.0 Mg•ha − 1 in pine stands in north-eastern Germany, i.e. a difference by 80%. Thus, the annual organic matter input with aboveground and belowground litter is considerably higher in beech than pine stands. In addition, the beech fine root system typically penetrates to greater depths than the often shallow fine root system of pine (Polomski & Kuhn 1998), which should lead to higher C inputs particularly in the subsoil under beech.
There is a striking mismatch between the larger SOC pool under pine and the relatively low annual C input with above-and belowground litter. This sheds some light on the relative importance of C input on the one side, and organic matter processing and C stabilization in the soil on the other side, for C sequestration in the soil (Angst et al. 2019). While the SOC accumulation of agricultural soils generally increases with the biomass input (Jenkinson 1988), this seems not to be the case in these forests on acidic, clay-poor soil. It has been found that the net accumulation of new tree-derived C can be greater under a low-productive stand on loamy soil than under a high-productive stand on fertile sandy soil (Hagedorn et al. 2004). According to biomass inventories in our stands, the pine forests stored only half of the aboveground biomass C that was recorded in the beech forests (Förster et al. 2021). Summing up the C pools in aboveground biomass and in the soil (without roots) gives on average an about 15% higher ecosystem C storage in the beech than the pine forests (246 Mg C•ha − 1 vs. 213 Mg C•ha − 1 ). This result may help to quantify the consequences for C sequestration of future changes in Scots pine plantation area in the lowlands of northern Central Europe. According to current forestry planning schemes, the area of pine plantations will shrink in northern Germany in future times to give way to more broad-leaved forests (predominantly beech) and Douglas fir (Pseudotsuga menziesii (Mirb.) Franco) plantations (Anders et al. 2004), which would increase ecosystem C storage. In the German state of Lower Saxony, for example, the forest area planted with pine is expected to decrease from 21% in 2015 to 15% in 2055 (NW-FVA & NL 2019).

Litter quality effects on the SOC pool
Litter quality is probably much more decisive for explaining the elevated SOC pools under pine than the amount of C transferred from biomass. The chemistry of root tissues may play a key role (Angst et al. 2018(Angst et al. , 2019). The C:N ratio of fresh foliage litter is typically by 15 to 40 g•g − 1 higher in pine than beech with N concentrations of 8-9 mg•g − 1 in beech and 5-7 mg•g − 1 in pine litter (Hertel 1995;Lorenz et al. 2004;Hobbie et al. 2006). Similar to needle litter, pine root litter has a low N content and high C:N ratio, and the lignin content is also high (Silver & Miya 2001). Correspondingly, the C: N ratio of soil organic matter was in our study by~3 g•g − 1 higher under pine than beech throughout the profile. As the C:N ratio is a good predictor of decomposition rate across different litter types (Enriquez et al. 1993;Pérez-Harguindeguy et al. 2000), one would assume that litter decomposition should proceed faster under beech than pine, which would explain the smaller SOC stocks under beech. Moreover, pine needle litter has a much lower Ca content than beech litter, which may hamper forest floor C dynamics (Hobbie et al. 2006). However, litter bag experiments with pine and beech foliage in neighboring stands showed that pine needles decomposed significantly faster in the first and also the second year in stands on sandy soil than beech litter (Hertel 1995;Lorenz et al. 2004;Hobbie et al. 2006). The higher decomposition rate of pine litter is unexpected not only from the higher C:N ratio, but also from the often more pronounced drying of the organic layer in pine forests as compared to beech forests (Leuschner 2001;Anders et al. 2004). One explanation could be the lignin and tannin content of fresh beech litter, which is considerably higher in beech than in pine leaf litter (Lorenz et al. 2004;Hobbie et al. 2006). Whatsoever the reason, the turnover of the labile and intermediate SOC pools is thus unrelated to the flux of C into the stable soil C pool in these forests (Poeplau et al. 2018).
It is astonishing that soil N concentration, soil C:P and N:P ratios, and base saturation at the cation exchangers all were similar in soil profiles under beech and pine in our study, even though the foliage and leaf litter of beech contains significantly more N and other nutrients, which are returned to the soil. The 'soil-ameliorating' effect of the more nutrient-rich beech litter (Anders et al. 2004;Prietzel 2004,) seems to be limited to the organic layer, affecting the mineral soil only marginally.
Our results suggest that the P-resin content may primarily depend on climatic factors, notably growing season precipitation. According to the meta-analysis of Hou et al. (2018), available P in forest soils generally decreases with increasing MAT and MAP, matching our results of a significant negative relationship between MSP and P-resin pools in the mineral soil of the pine stands. Possible drivers of this relationship are a higher P leaching when precipitation is high, and decreased P uptake in drier soil. Our data do not show a clear temperature effect on P-resin pools.
As beech leaves are typically richer in lignin than pine needles (Berg & McClaugherty 2003;Hobbie et al. 2006), it is probably not lignin that fosters the transfer of pine litter components to the stable SOC pool. It is also unlikely that differences in the stabilization of SOM (e.g. the abundance of reactive clay mineral surfaces) have caused the higher SOC pools under pine, given that the geological substrates and associated soil mineral contents were similar among the paired beech and pine sites. This is supported by the fact that the profiles under pine and beech had very similar cation exchange capacities. Other, probably local, factors must be responsible for the higher soil C accumulation under pine than beech. One factor could be stand history.
The combined influence of forest history, climate, and soil Due to century-or millennia-long human impact on the forests, land use history is often a main determinant of the actual SOC pools in temperate forests (Verheyen et al. 1999;Compton & Boone 2000;Leuschner et al. 2013). Many pine plantations in northern Germany have been established on former heathland or arable land during the last 250 years, causing re-accumulation of C especially in the organic layer, but also in the mineral soil (Heinsdorf 2002;Vesterdal et al. 2002). The linear models developed to test for the combined effects of climate, soil pH, stand age, and tree species in conjunction with stand continuity (ancient beech, ancient pine, recent pine) revealed no consistent effects of soil (pH) and climate (temperature and precipitation) on SOC stores. Rather, the models showed a dominant influence of tree species together with stand continuity, which masked any significant relationship between SOC pool size and environmental factors that was visible in the beech and the pine data sets in correlation analyses. Our data set, completed by seven pine stands with less than 150 years of forest continuity ('recent pine stands'), thus shows a large increase in organic layer C org (by on average 30 Mg C•ha − 1 ) from recent to ancient pine stands (> 230 years of continuity) and a smaller (non-significant) increase in the mineral soil (Figs. 7a-c). The C accumulation corresponds to a parallel accumulation of N mainly in the organic layers (Fig. 7).
The extended data set (16 pine and 9 beech stands) further indicates significantly higher plant-available P stocks (P resin ) in the mineral soil of the recent pine stands (Fig. 8a), which may point at agricultural use of the sites in the distant past. Thus, legacy effects seem to play an important role for P availability in the studied forests, besides the dependence of available P on precipitation. The significant positive relation between the pool of exchangeable base cations and MAT in the beech stands and the positive effect of MAT on the base cation pool in the model suggest that a warming climate might affect base cation availability in these forests, perhaps through altered weathering rates (Akselsson et al. 2016;Houle et al. 2020). Despite an assumed role for temperature and precipitation for the availability of P and base cations in our stands, the temperature and precipitation gradients covered in our study (8.4°C-9.1°C and 560-820 mm•yr − 1 ) may be too short to significantly influence the soil stocks of C and N in our sample in relation to the influence of the other factors. On the national level, climate was also found to be only of secondary importance for soil C stocks in the soil survey of German forests (Grüneberg et al. 2019). Thus, attempts to explain SOC stocks by edaphic and climatic factors without considering forest history may in our study region be not comprehensive enough.

Conclusions
Our study suggests that the tree species effect on soil C and N accumulation can hardly be deduced from forest productivity and litter decomposition rate in the studied forest ecosystems. Regarding the climate warming mitigation potential of pine and beech forests in the study region, it is an important finding that the SOC stocks are on average larger under pine than beech. However, this must be put in relation to the C stock in the aboveground biomass, which is considerably lower in pine than beech. Although covering large parts of Northern Germany, it must be kept in mind that this study is a case study and the results should be interpreted under the constraints that local variation in geological substrate, forest management and forest history may have significantly affected current SOC stocks. In agreement with other studies, it is evident from this study that tree species effects on soil C relations often are closely linked to forest histories and management regimes, which makes it difficult to disentangle biological effects related to productivity, litter chemistry, and root system properties from management effects. This is even more so in large-scale forest soil inventories, where a strong interaction between tree species and substrate type likely is obscuring tree species effects on soil C storage.