1 Introduction

The availability of nutrients is one of the leading environmental factors that limit the productivity of terrestrial ecosystems (Lambers and Oliveira 2019). In recent decades, nitrogen oxide (Nox) emissions into the atmosphere have changed nitrogen (N) availability in many forest ecosystems in Europe and North America (Crowley et al. 2012; Jonard et al. 2015; Talkner et al. 2015). An increased N pool in forests affects the mineral nutrition status of trees (Driscoll et al. 2003). Increased N availability stimulates tree growth, leading to phosphorus (P) deficiency (Akcelsson et al. 2008; Reich et al. 2006). Although total P resources in soils are relatively high, only a small proportion of this is available to plants (Hinsinger 2001). The content of N and P in leaves is perceived to be an indicator of the nutritional status of plants, and the N:P ratio can shed light on potential limitations—resulting from nutrient deficiencies—that affect the productivity of terrestrial ecosystems and plant growth (Reich and Oleksyn 2004). Numerous studies have shown that the content of N and P in leaves reflects the content of these elements in the soil (Aerts and Chapin 2000). One of the main stages in the circulation of nutrients between the plant and soil is leaf decomposition. In addition, the nutrient contents of the litter, soil properties, and climatic conditions are the primary factors influencing the accumulation rate of dead organic matter (Manzoni et al. 2008; Parton et al. 2007).

The accumulation of dead organic matter is crucial for soil formation in soilless areas (Woś and Pietrzykowski 2015). In mining areas, soil reconstruction and afforestation are potentially the best strategies for restoring post-mining areas to their previous state (Dutta and Agrawal 2003; Pietrzykowski and Krzaklewski 2007). From an ecological point of view, the goal of restoring post-mining areas is the development of a stable forest ecosystem (Bradshaw and Hüttl 2001). In order to alleviate nutrient deficiencies in reconstructed soils, tree species capable of binding atmospheric N are introduced, and thanks to symbiosis with bacteria, these species enrich the soil with N (Pietrzykowski et al. 2015; Vlachodimos et al. 2013; Zipper et al. 2011). In recent years, there has been increased interest in the selection of tree species for this purpose and in their adaptation to difficult habitat conditions in post-mining areas (Baumann et al. 2006; Chodak and Niklińska 2010; Pietrzykowski and Socha 2011).

The role of N-fixing species in alleviating P deficiencies has not been investigated in depth (Hinsinger et al. 2011). Plants play a central role in the transformation and circulation of P in the soil. Moreover, plants can take up P only in mineral forms and deposit it in organic forms, which become substrates for organic forms of P in the soil (Hinsinger et al. 2011). Although organic P may constitute up to 50% of the total P in the topsoil, it is not directly available to plants (Sulieman and Mühling 2021). In nature, however, several mechanisms allow plants to access organic forms of P. By increasing phosphatase activity and increasing P absorption by microorganisms, some tree species have improved their ability to absorb P from the soil (Richardson et al. 2009).

As global climate change can cause changes in the circulation of elements, it is vital to understand how soil and habitat conditions affect the mineral nutritional status of trees (Schaaf 2001). Such knowledge is essential at sites restored after severe disturbance since conditions are often harsh and nutrient supply is insufficient for healthy plant growth. Of particular concern are the low concentrations and availability of essential nutrients in the parent materials of technosols. N and P are considered the primary limiting elements for the development of vegetation planted on post-industrial barrens (Zipper et al. 2011; Manimel Wadu et al. 2017; Cross et al. 2019). The current knowledge on the impact of chemical properties and microbiological activity on forest tree N and P supply is mainly based on research material from intact areas. In contrast, reports from reclaimed sites are relatively scarce. Determining the factors influencing the supply of N and P to trees is crucial in monitoring stands mineral nutrition status. Therefore, this study aimed to determine the influence of microbiological and physicochemical properties of soils formed from different substrates on the mineral nutrition status of pioneer tree species (Scots pine and silver birch) and N-fixing tree species (black locust and black alder), expressed through the supply of nitrogen and phosphorus to their leaves. We hypothesized that the status of tree N and P supply in restored forest ecosystems depends on the chemical properties of parent materials and the microbial activity of developing soils.

2 Materials and Methods

2.1 Study Sites and Field Study

The research was carried out in stands growing on reclaimed post-mining soils (RMS) composed of various types of parent materials representing sands (sand quarry, Szczakowa; open-pit lignite mine, Bełchatów), clays (open-pit lignite mine, Turów; open-pit sulfur mine, Piaseczno), and ashes (combustion waste disposal site, Lubień). The age of forest stands growing on the ashes was 18–20 years, those growing on the sands 30–35 years, and those growing on clays 38–44 years. Detailed characteristics of the examined objects are presented in Table 1.

Table 1 Characteristics of the tested objects

2.2 Soil Sampling

To determine the factors influencing the state of N and P supply to trees in the designated parent material variants, for stands of age classes I and II (class I, stands up to 20 years old; class II, 21–40 years), permanent research plots with an area of 100 m2 have been created, within which there are single-species stands of black locust, black alder, common birch, and Scots pine. The experiment was fully randomized. At the turn of August and September, leaf samples were collected from the sunlit part of the crown. The samples were placed in PVC bags, transported to the laboratory, and dried. Five soil samples distributed in the envelope system were taken on the surface of each of the 0–5 cm and 5–20 cm layers, from which one composite soil sample was prepared. After collection, the samples were transported to a laboratory and separated into samples for physicochemical and biochemical analyses. The samples for physicochemical and chemical analysis were dried to constant weight and passed through a 2-mm mesh before analysis. The samples for biochemical analysis were stored in a refrigerator at 4 °C.

2.3 Determination of Soil Physicochemical Properties and Biochemical Analyses

Soil texture was measured spectrally using the Fritsch GmbH Laser Particle Sizer ANALYSETTE 22. The pH was measured potentiometrically in H2O and 1 M KCl. The contents of N and carbon (C) were measured using the dry combustion method with a TruMac® CNS analyzer. Exchangeable K+ and Mg2+ were extracted with 1 M CH3COONH4 and determined using an iCAP 6000 Series ICP OES spectrometer. The total Mg, K, and P were measured using the same instrument after extraction with a mixture of HNO3 and HClO4 at a ratio of 3:1.

Acid and alkaline phosphatase (phosphomonoesterase) activities (AcdPhoSP and AlkoPhoSP) and inorganic pyrophosphatase (PyroPho) were measured as described by Margesin (1996). Briefly, the soil samples (1 g d.w.) were mixed with disodium p-nitrophenyl phosphate solution (115 mM) and a buffer solution (pH = 6.5 for AcdPhoSP; pH = 11.0 for AlkPhoSP) and incubated at 37 °C for 1 h. The p-nitrophenol released by phosphatase activity was extracted, reacted with NaOH, and determined photometrically at 400 nm. Activities of AcdPhoSP and AlkPhoSP were expressed per dry soil mass in units of μg p-NP·g−1·h−1. Total phosphatase activity (PhoSuma) was expressed as the sum of AcdPhoSP and AlkPhoSP. For PyroPho measurement, the soil samples (1 g d.w.) were mixed with buffered pyrophosphate solution and incubated for 1 h at 37 °C. Orthophosphate released by inorganic pyrophosphatase activity was extracted with sulfuric acid, reacted with ammonium molybdate, and determined photometrically at 700 nm. PyroPho activity was expressed per dry soil mass in units of μg P-PO43-·g−1·h−1.

2.4 Statistical Procedures

Differences in the mean values of the basic soil properties, and results of biochemical analyses between different types of parent materials, were tested using one-way ANOVA followed by an RIR-Tukey multiple comparison procedure (at p < 0.05). Additionally, two-way ANOVA was used to test for differences in the N, P, K, and Mg contents, and N-to-P ratios, in the leaves of the studied tree species growing on different substrates. Before statistical analyses, test assumptions were checked using the Shapiro-Wilk test (for normal distribution) and the Brown-Forsyth test (for homogeneity of variances). Statistical analysis included the implementation of generalized additive models (GAM). The relationships between N and P content in leaves and the predictor variables Ntot and Mg2+ and the activity of AcdPhoSP were modeled using GAM models, a non-parametric extension of the generalized linear model (GLM). An identity join function was specified for normally distributed response data. Predictor variables were entered into the models one at a time. The selection of the best model (i.e., the model containing independent variables explaining the largest amount of observed variance) was made based on the adjusted model determination/efficiency coefficient (R2adj). Smoothing spline function plots of GAM models were used to illustrate the influence of individual predictor variables on dependent variables. The remaining analyzed variables (pH KCl, sand, silt, clay, Ctot, Ptot, Ktot, Mgtot, AlkPhoSP, PhoSuma, and PyroPho) were not significant or showed collinearity. Statistical analysis was performed using Statistica 13.3 software (StatSoft Inc.) (Table 2).

Table 2 Approximate significance of seven predictor variables on the periodic annual volume increment described using the GAM (generalized additive model)

3 Results

3.1 Basic Soil Properties

There were significant differences in texture between the tested substrates in the 0–5 cm soil layer. The tested substrates were characterized by their reactivity as very strongly acid (SS), moderately acid (CS), or slightly alkaline (AS). The sand content in the SS, AS, and CS substrates was 81%, 70%, and 29%, respectively. In addition, there were differences in the content of Ctot and Ntot between the tested substrates. The highest content of Ctot and Ntot was found in the CS substrate and the lowest in the SS substrate.

The AS and CS substrates contained more Ptot and Mgtot. The Mgtot content in the SS variant was ten times lower than in the AS and CS variants. In addition, higher concentrations of K+ cmol·kg−1 and Mg2+ cmol·kg−1 were found in AS and CS soils than in SS soils. AcidPhoSP and PyroPho activities were higher in the CS substrate than in the AS and SS substrates. The activity of AlkPhoSP was lower in the SS substrate than in the AS and CS substrates.

There were also significant differences in texture between the substrates in the 5–20 cm layer. The sand content in the SS substrate was 80%, in the AS substrate 85%, and in the CS substrate only 21%. There were differences in the content of Ctot and Ntot between the tested substrates. The highest content of Ctot and Ntot was found in the CS substrate and the lowest in the SS substrate. The CS substrate contained more Ptot and Mgtot. The Mgtot content in the SS variant was ten times lower than in the CS variants. In addition, CS soils contained higher concentrations of K+ cmol∙kg−1 and Mg2+ cmol∙kg−1 than SS soils. AcidPhoSP, PyroPho, and AlkPhoSP activities were higher in the CS substrate than in the AS and SS substrates (Table 3).

Table 3 One-way ANOVA for the basic soil properties (mean ± standard errors) in two soil layers (0–5 cm and 5–20 cm)

3.2 Mineral Nutritional Status and Quantification in Leaves

In substrate variants AS and CS, the N content in alder and black locust leaves was higher than in birch and pine leaves. In the SS substrate variant, the highest content of N was found in black locust leaves and the lowest in pine leaves. The type of substrate did not affect the N content in the leaves of birch, pine, and alder. A higher N content was found in the black locust leaves in the SS substrate variant than in the AS and CS substrate variants. Furthermore, no differences in P content were found in the birch, alder, pine, and black locust leaves in the AS and CS variants. In the SS substrate variant, pine leaves contained less P than alder and black locust leaves. The substrate type also affected the K and Mg content in the leaves of alder, birch, and black locust. K in birch and black locust leaves was higher in variants CS and SS than in variant AS. A higher content of Mg was found in pine, black locust, and birch leaves in the AS and SS substrate variants compared with the CS substrate. The lowest N:P ratio was found in pine leaves and the highest in black locust leaves (Table 4).

Table 4 Two-way ANOVA for the N, P, K, and Mg contents and N-to-P ratios (mean values ± standard errors) in leaves studied tree species growing on different substrates (N, contents of nitrogen; P, contents of phosphorus; K, contents of potassium; Mg, contents of magnesium)

Analysis of individual soil properties showed that N content was strongly influenced by the contents of Ntot and Mg2+ and by the activity of AcdPhoSP in the soil (Fig. 1). Increasing Ntot in the topsoil layer (0–5 cm) caused a substantial increase in the N content of the leaves. However, further growth was much slower above 2.5 mg·g−1 Ntot (Fig. 1a). A relatively steep decrease in the N content of leaves was observed with the increase of Mg2+ from 2 to 110 cmol·kg−1. Furthemore, the same change in Mg2+ caused a decrease of the N supply to trees by 20 mg·g−1 (Fig. 1b). The increase in the activity of AcidPhoSP in the range 55 to 500 ug NP·g-1·dm·h-1 also contributed to the improvement of the N supply in trees (Fig. 1c). However, a further increase in the activity of AcidPhoSP did not significantly increase the content of N in the leaves. The analysis also revealed the influence of soil properties on the P content in leaves. A positive effect of the AcidPhoSP activity on P content in the leaves was apparent up to 500 ug NP·g-1·dm·h-1; further increase in AcidPhoSP activity did not affect P content in leaves (Fig. 1d).

Fig. 1
figure 1

a Varied nitrogen (N) content in leaves depending on nitrogen (N) content in the soil, b varied nitrogen (N) content in leaves depending on exchangeable magnesium (Mg2+) concentration in the soil, c Varied nitrogen (N) content in leaves depending on acid phosphatase (AcidPhoSP) activity in the soil, d varied phosphorus (P) content in leaves depending on acid phosphatase (AcidPhoSP) activity in the soil (shown using generalized additive models). An identity join function was specified for normally distributed response data. R2adj models = 80%

4 Discussion

As expected, the leaves of alder and black locust contained more N than those of birch and pine, on each tested substrate. Compared to the available literature, the concentration of N in the leaves of the examined tree species was similar or high in relation to the values obtained on natural and reclaimed post-mining areas. In the forests of Europe, the concentration of N in alder leaves varies from 20 to 40 mg·g−1 (Kuznetsova et al. 2010; Lorenc-Płucińska et al. 2013; Uri et al. 2002). The N content in the leaves of black alder found in the combustion waste dump of the Bełchatów power plant was 24 mg·g−1, and the N content in the leaves of alder found on extremely sterile sands was 29 mg·g−1. The N content of the birch leaves was high compared to data in previous studies (Hytonen et al. 1995; Saarsalmi 1995). Studies in Estonia show that, with high N soil concentration, the N concentration in birch leaves can be similar to or even higher than in alder leaves. In the forests of the temperate climate zone, the N content in pine needles ranges from 11.7 to 29.1 mg·g−1 (Fober 1993). Similar results were obtained for Scots pine stands on extremely poor sands reclaimed for forestry, where the N concentration was 12.1 mg·g−1 (Pietrzykowski et al. 2013). N content in black locust leaves under natural conditions is about 32 mg·g−1, compared with 30 mg·g−1 in the Bełchatów furnace waste landfill (Woś et al. 2020). However, the 25-year-old black locust stands found in unpolluted areas of eastern Bulgaria were characterized by lower N leaf content, amounting to 21 mg·g−1 (Tzvetkova and Petkova 2015). Furthermore, a similar N content in black locust leaves was obtained in a reclaimed opencast lignite mine in Appalachia (Brinks et al. 2011).

As the deficient component, the content of P was lower than the adopted limit numbers for the supply of this element for the tested tree species. The optimal P content for pine, black locust, alder, and birch is 1.8, 2.0, 3.6, and 3.7 mg·g−1, respectively. The P content in the leaves of alder found in the furnace waste landfill was 1.1 mg·g−1 and in the oil shale extraction area in northeastern Estonia, 1.9 mg·g−1 (Kuznetsova et al. 2010). P content in birch leaves was lower than on post-agricultural lands in Estonia (3890 mg·g−1) (Uri et al. 2007). Furthermore, the P content in birch leaves in the oil shale mining area was 3.2 mg·g−1 (Kuznetsova et al. 2010). The P content in black locust leaves was similar to that determined in the USA waste landfill (Sopper 1992). Moreover, the concentration of P in the leaves of black locust growing in areas affected by industrial emissions ranged from 1680 to 2530 mg·g−1 (Tzvetkova and Petkova 2015). A much higher P content (3.8 mg·g−1) in black locust leaves was found in the areas of a reclaimed coal mine in the USA (Marschner 2012). The P content of pine needles found in natural forests of Europe ranges from 1.52 mg·g−1 for the central population to 1.56 mg·g−1 for the northern population (Oleksyn et al. 2002). A similarly low P content has been found in pine needles on Poland’s reclaimed land, ranging from 1.03 mg·g−1 on coal-bearing shales to 1.12 mg·g−1 on neogene sands with loam and clay, both carbonated and carbonated sulfurized (Pietrzykowski et al. 2013).

The content of individual elements in the assimilation apparatus and the mutual proportions between them are important in assessing the mineral nutrition status of trees (Baule and Fricker 1971; Pietrzykowski et al. 2013). A N:P ratio above 14 was determined to limit plant growth due to P deficiency, while a value below 14 limited growth due to N deficiency. Atmospheric N-fixing trees (alder and black locust) contained a much higher N:P ratio in leaves than birch and pine. Low concentrations of P in leaves, and high N:P ratios, indicate an insufficient supply of P for all species except pine and birch in the SS substrate. A decrease in the P concentration in leaves may cause health problems for trees: a N:P ratio above 14.8 for deciduous trees and above 7.3 for conifers can cause an increase in defoliation (Verosoglou et al. 2014). Many studies show that an unbalanced increase in the concentration of C and N in relation to P significantly impacts ecosystem functioning (Sardans et al. 2012).

The content of N in the soil affected the content of N in the leaves of the examined tree species across the entire range of analyzed concentrations. Studies conducted on a global scale have also shown that soil N content is a good predictor of phosphatase activity (Margalef et al. 2017). Therefore, acid phosphatase concentration can be a good predictor of N and P concentration in leaves. An increase in AcdPhoSP activity in the soil to about 500 ug NP·g-1·dm·h-1 favors an increase in the concentration of N and P in leaves (see Fig. 1c, d). The stimulating effect of N on the AcdPhoSP content is particularly visible in soils poor in this element (Olander and Vitousek 2000). The relationship between AcdPhoSP activity and N and P concentrations in leaves may also be due to intensive accumulation of N in soils forested with N-fixing trees increasing acid phosphatase activity, which results from the greater availability of N necessary for the synthesis of this enzyme (Olander and Vitousek 2000).

5 Conclusions

Our research showed differences in the content of nitrogen in the leaves of black locust, black alder, common birch, and Scots pine occurring on the same substrate. As expected, the alder and black locust leaves contained more nitrogen than the leaves of birch and pine on each tested substrate. However, there were no differences in the phosphorus content in the leaves of these tree species. The results confirm that nitrogen content was strongly influenced by the content of nitrogen total and exchangeable magnesium and the activity of acid phosphatase in the soil. Moreover, acid phosphatase activity can be a good predictor of nitrogen and phosphorus concentration in leaves. An increase in the activity of acid phosphatase in the soil to about 500 ug NP·g-1 dm·h-1 favors an increase of nitrogen and phosphorus concentrations in the leaves.