Linking Soil Acidity to P Fractions and Exchangeable Base Cations under Increased N and P Fertilization of Mono and Mixed Plantations in Northeast China

: Atmospheric N deposition is increasing worldwide, especially in China, signiﬁcantly a ﬀ ecting soil health, i.e., increasing soil acidiﬁcation. The northern region of China is considered to be one of the N deposition points in Asia, ranging from 28.5 to 100.4 N ha − 1 yr − 1 . Phosphorus (P) is the limiting factor in the temperate ecosystem and an important factor that makes the ecosystem more susceptible to N-derived acidiﬁcation. However, it remained poorly understood how the soil acidiﬁcation process a ﬀ ects soil P availability and base cations in the temperate region to increased N deposition. To address this question, in May 2019, a factorial experiment was conducted under N and P additions with di ﬀ erent plantations in Maoershan Experimental Forest Farm, Northeast China, considering species and fertilization as variables. The e ﬀ ective acidity (EA) increased by N and NP fertilizations but was not signiﬁcantly a ﬀ ected by P fertilization. Similarly, the pH, base saturation percentage (BS%), calcium (Ca 2 + ), and magnesium (Mg 2 + ) were decreased under N addition, while the Al:Ca ratio increased, whereas NaHCO 3 inorganic phosphorus (Pi) and NaOH organic phosphorus (Po) signiﬁcantly decreased under N enrichments. However, NaOH Pi increased in N-enriched plots, while H 2 O Pi and NaHCO 3 Pi increased under the P addition. Thus, the results suggest that the availability of N triggers the P dynamics by increasing the P uptake by trees. The decrease in base cations, Ca 2 + , and Mg 2 + and increase in exchangeable Fe 3 + and Al 3 + ions are mainly responsible for soil acidiﬁcation and lead to the depletion of soil nutrients, which, ultimately, a ﬀ ects the vitality and health of forests, while the P addition showed a bu ﬀ ering e ﬀ ect but could not help to mitigate the soil acidity.


Introduction
The industrial revolution has caused an increase in atmospheric nitrogen (N) deposition around the globe [1]. According to the statistics, the deposition rate will be double by 2050 to that of 1995, i.e., 100 Tg N yr −1 -200 Tg N yr −1 [2]. Similarly, the mean wet N deposition rate in China has increased to 25% from 1990 to 2000 [3], and a dramatic increase is expected in the coming future [4]. Pan et al. [5] also documented that the northern part of China is considered one of the N deposition points in Asia, ranging from 28.5 to 100.4 N ha −1 yr −1 . According to the China Statsitcal Year Book, the fossil combustion cause a major contribution in NOx in the atmosphere and increased approximately six times in the last 38 years (1980-2018) [6]. The ecosystem's functions are greatly influenced by The soil is classified as Hap-Boric Luvisol in nature: high organic matter, up to 50-cm-deep, 1-10-cm loamy soil, and 10-20-cm sandy loam [13].
The study selected different plantations of Larix gmelinii (Rupr.) Rupr. (Dahurian larch; Larch), Juglans mandshurica Maxim. (Manchurian walnut; Walnut), and Fraxinus mandshurica Rupr. (Manchurian ash; Ash) in the monoculture and Ash× Lar and and Wal × Lar in the mixed culture (with 5 rows of larch with 3 rows of broadleaved species, alternately). The J. mandshurica and F. mandshurica are classified as angiosperms, having arbuscular mycorrhizal colonization, while the L. gmelinii was classified as a gymnosperm and having ectomycorrhizal colonization [50]. These plantations were planted with two-year-old seedlings at a spacing of 1.5 m × 2.0 m in 1987. Thus, a factorial experiment was conducted with five species and four treatments in a randomized fashion. The experiment was replicated three times, and each replication was considered as a big plot, i.e., 20 m × 20 m. Each big plot (replicate) was further subdivided into 3 subplots, and to each subplot, 20 sub-subpots were randomly allotted. Therefore, the experiment consisted of 180 experimental subunits. The application of nitrogen (N) as urea (CO(NH2)2) and phosphorus (P) as diammonium phosphate ((NH4)2 HPO4) fertilizers were carried out at 4 levels in pellet form per species. The subplots per species were designated as Control (C; no addition), nitrogen (N; 20-g N m −2 ), phosphorus (P; 10-g P m −2 ), and combination (NP; 20-g N m −2 and 10-g P m −2 ) [13,51,52]. The N was applied in three split doses (i.e., 30%, 40%, and 30%; May-Aug), preventing leaching, while P was applied at once [53]. Three sampling points in the sub-subplots were randomly selected to collect soil samples at 0-20-cm soil depth and made a composite sample. See Tables 1 and 2 for the preharvest data of stand characteristics and soil chemical properties.  (Manchurian ash; Ash) in the monoculture and Ash× Lar and and Wal × Lar in the mixed culture (with 5 rows of larch with 3 rows of broadleaved species, alternately). The J. mandshurica and F. mandshurica are classified as angiosperms, having arbuscular mycorrhizal colonization, while the L. gmelinii was classified as a gymnosperm and having ectomycorrhizal colonization [50]. These plantations were planted with two-year-old seedlings at a spacing of 1.5 m × 2.0 m in 1987. Thus, a factorial experiment was conducted with five species and four treatments in a randomized fashion. The experiment was replicated three times, and each replication was considered as a big plot, i.e., 20 m × 20 m. Each big plot (replicate) was further subdivided into 3 subplots, and to each subplot, 20 sub-subpots were randomly allotted. Therefore, the experiment consisted of 180 experimental subunits. The application of nitrogen (N) as urea (CO(NH 2 ) 2 ) and phosphorus (P) as diammonium phosphate ((NH 4 ) 2 HPO 4 ) fertilizers were carried out at 4 levels in pellet form per species. The subplots per species were designated as Control (C; no addition), nitrogen (N; 20-g N m −2 ), phosphorus (P; 10-g P m −2 ), and combination (NP; 20-g N m −2 and 10-g P m −2 ) [13,51,52]. The N was applied in three split doses (i.e., 30%, 40%, and 30%; May-Aug), preventing leaching, while P was applied at once [53]. Three sampling points in the sub-subplots were randomly selected to collect soil samples at 0-20-cm soil depth and made a composite sample. See Tables 1 and 2 for the preharvest data of stand characteristics and soil chemical properties.

Soil Sampling Strategy
The soil samples were collected in October 2019 after removing the litter from the forest floor for each subplot, and randomly, three cores were selected, a soil sample taken from 0-20-cm soil depth mixed and put in Ziploc bags and transported to the laboratory. The soil samples were further subdivided into three parts. One subsample set was air-dried, grinded, and sieved through 2 mm to measure exchangeable cations and NH 4 + -N NO 3 − -N [24]. In contrast, the second set of subsamples was processed through a 0.5-mm sieve and used to measure P fractions and soil pH. The third set was simultaneously oven-dried and processed through a 0.25-mm sieve to estimate the soil organic carbon and total nitrogen [13].

Soil Chemical Analyses
The soil organic carbon and nitrogen were determined through the dry combustion method using a Vario EL III elemental analyzer, i.e., Elementar Analysensysteme GmbH, Hanau, Germany. The NH 4 + -N (ammonium) and NO 3 − -N (nitrate) were obtained using 2-M KCl, as determined by an auto analyzer III Bran, Luebbe GmbH, Norderstedt Germany [13]. One-step method i.e., 0.1-mol L −1 BaCl 2 (50:1, solution:soil) was used to examine the exchangeable cations H + , K + , Na + , Ca 2+ , Mg 2+ , Al 3+ , Fe 3+ , and Mn 2+ . The extract was filtered with an acetate filter. The H + ion was determined with a glass electrode, while the rest of the cations were measured by coupled plasma mass spectrometry (Agilent Technologies Co. Ltd., Santa Clara, CA, USA) [24]. To determine the soil pH, 1-g soil was mixed with distilled water at a ratio 1:2.5 and shaken on a mechanical shaker (30 min). The soil pH was measured using the calibrated pH meter (S220 Seven Compact pH Meter, Shanghai, China) [13]. Additionally, we calculated the effective acidity (EA), and it was the combination of Al and H. For the analysis of the P fraction, a method modified by Hedley et al. [54] was used, and soil samples (01g) were taken into a 50-mL centrifugal tube with an added 30 mL of distilled water, along with 05 drops of chloroform. The solute was shaken for 18 h mechanically and used for organic and inorganic H 2 O determinations. The resultant soils in tubes were again filled with 0.5-mol/L NaHCO 3 , with an added 05 drops of chloroform, and shaken for 18 h. The extracts were used for organic (Po) and inorganic (Pi) phoaphorus, NaHCO 3 determination. The tubes with remaining soil were again filled with 30 mL of distilled water, 0.1-mol/L NaOH, and an added 05 drops of chloroform. The solute was shaken for 18 h for the determination of Pi and Po. The same procedure was repeated for the determination of Pi and Po for 1-M HCl. While for the determination of residual P, the Tissen et al. [55] method was used, employing H 2 SO 4 and H 2 O 2 in the digestion chamber.

Statistical Analysis
ANOVA (Analysis of variance) of factorial design for species, fertilization, and their interactions were carried out to examine whether the N and P fertilizers affected the phosphorus fractions, exchangeable cations, or otherwise. The Shapiro-Wilk test was used to check the normality, while the post-hoc Tukey honest significant difference (HSD) test was applied to compare treatments. Heatmap for Pearson's correlation coefficients was performed to determine positive or negative correlations within soil chemical properties and soil P forms using the "ggplot2" and "reshape2" packages in "R" v. 3.6.1 [57]. A redundancy analysis (RDA) was used to identify the soil properties (i.e., soil acidity) that predicted the variations in P fractions and base cations. The RDA was performed with the "rda" function by using the "vegan" package in R v. 3.6.1 [57]. All statistical analyses were performed in the program R v. 3.6.1 [57]. Sigma Plot v. 12.5 (Systat software Inc., San Jose, CA, USA) was used for graphical representations. Mean data were displayed as mean ± SE. If not noted otherwise, a significance level was reported at 0.05.

Soil Chemical Properties
Different elements showed variations concerning species and fertilizer application ( Figure 2, Table 3, and Supplementary Table S1). However, of the selected soil chemical properties, i.e., the H + ion, exchangeable ions Al 3+ and Fe 3+ , EA (the combination of Al and H ions), Al:Ca, soil organic carbon (SOC), C:N ratios, and NO 3 − -N significantly increased (p < 0.05) by N and NP fertilizations but remained stable under P application relative to the control in all plantations (Figure 2b-f and Tables 3  and 4), while the Mg 2+ and Ca 2+ ions were significantly decreased (p < 0.05) by N addition but remained stable under P and NP fertilization in all plantations (Tables 2 and 3). The base saturation percentage (BS%) was also affected (p < 0.05) in N and NP addition plots when compared to that of the control, while remaining stable under the P addition across the plantations (Figure 2a). Similarly, the soil pH significantly decreased by the soil N enrichment but was not significantly affected by either P or NP fertilizations across all plantations (Tables 3 and 4), whereas the total P was not significantly affected by any of the fertilization's treatments in all plantations (Supplementary Table S1 and Figure S1f).

Soil P Fractions
All phosphorus fractions were significantly affected by species and fertilization ( Table 3). The H 2 O Pi significantly increased (p < 0.05) by P fertilization but remained stable under N and NP fertilizations across all plantations (Figure 3a). However, NaHCO 3 Pi, the labile P, significantly decreased (p < 0.05) under N enrichment while increasing under P fertilization but remained stable under NP application across all plantations (Figure 3b). Similarly, the NaOH Pi increased significantly (p < 0.05) by N enrichment while remaining stable under P and NP fertilizations in all plantations (Figure 3c), whereas NaOH-Po significantly decreased (p < 0.05) by N and NP additions but remain unaffected under P applications across all plantations (Figure 3d). However, the H 2 O Po, NaHCO 3 Po, HCl Pi, HCl Po, and residual P were not significantly affected by any of the fertilization treatments in all plantations (Supplementary Figure S1a-e and Table S1).

Soil P Fractions
All phosphorus fractions were significantly affected by species and fertilization ( Table 3). The H2O Pi significantly increased (p < 0.05) by P fertilization but remained stable under N and NP fertilizations across all plantations (Figure 3a). However, NaHCO3 Pi, the labile P, significantly decreased (p < 0.05) under N enrichment while increasing under P fertilization but remained stable under NP application across all plantations (Figure 3b). Similarly, the NaOH Pi increased significantly (p < 0.05) by N enrichment while remaining stable under P and NP fertilizations in all plantations (Figure 3c), whereas NaOH-Po significantly decreased (p < 0.05) by N and NP additions but remain unaffected under P applications across all plantations (Figure 3d). However, the H2O Po, NaHCO3 Po, HCl Pi, HCl Po, and residual P were not significantly affected by any of the fertilization treatments in all plantations (Supplementary Figure S1a-e and Table S1).

Correlation between Soil P Fractions and Chemical Properties
Among plantations, the RDA showed the interrelationship of soil P fractions and base cations with soil acidity and the responses to additional N and P fertilizations (Figure 4). The axis 1 and axis 2 of the RDA accounted for 72.6% and 5.6% of the total variations, respectively. Based on our aim, the

Correlation between Soil P Fractions and Chemical Properties
Among plantations, the RDA showed the interrelationship of soil P fractions and base cations with soil acidity and the responses to additional N and P fertilizations (Figure 4) Figure 4). Besides, the SOC and CN ratio were strongly positively related to the Al:Ca ratio. Based on the RDA, up to 77.7% of the total variation in soil P fractions and other chemical properties were explained by N and P availability. In contrast, Al 3+ , NO 3 − -N, and EA showed weaker responses.
Forests 2020, 11, x FOR PEER REVIEW 10 of 19 factors responsible for soil acidification, like Fe 3+ , H + , NH4 + -N, Mn, TN, and pH, strongly correlated to NaOH Pi, HCl Pi, NaHCO3 Pi, Resi P, total phosphorious (TP), NaOH Po, H2O Po, NaHCO3 Po, and H2O Pi while weakly correlating to HCl Po, Mg, C, BS%, and Na ( Figure 4). Besides, the SOC and CN ratio were strongly positively related to the Al:Ca ratio. Based on the RDA, up to 77.7% of the total variation in soil P fractions and other chemical properties were explained by N and P availability. In contrast, Al 3+ , NO3 − -N, and EA showed weaker responses. Looking at the trait interrelations between the mono-and mixed ( Figure 5 and Supplementary  Figures S2 and S3) cultures, the total P was found to be significantly positively related to all P fractions in the mixed culture. Still, in the monoculture, HCl Po was not significantly related to total P. All the P fractions were significantly related to each other. Still, HCl was not significantly related to H2O Po, Looking at the trait interrelations between the mono-and mixed ( Figure 5 and Supplementary  Figures S2 and S3) cultures, the total P was found to be significantly positively related to all P fractions in the mixed culture. Still, in the monoculture, HCl Po was not significantly related to total P. All the P fractions were significantly related to each other. Still, HCl was not significantly related to H 2 O Po, NaHCO 3 Pi, NaOH Pi, and NaOH Po in the monoculture, while, in the mixed culture, NaHCO 3 Pi was not significantly related to H 2 O Po or NaOH Pi to NaHCO 3 Pi or HCl Po to H 2 O Po, H 2 O Pi, NaOH Po, and NaOH Pi, respectively. However, K was significantly related to all P fractions in the monoculture but not in the mixed culture, except for HCl Po.
was not significantly related to H2O Po or NaOH Pi to NaHCO3 Pi or HCl Po to H2O Po, H2O Pi, NaOH Po, and NaOH Pi, respectively. However, K was significantly related to all P fractions in the monoculture but not in the mixed culture, except for HCl Po.
Similarly, Mg was significantly positively related to NaHCO3 Pi and HCl Po in the monoculture. In contrast, in the mixed, it was found to be significantly related to NaOH Po. Interestingly, Na was significantly positively related to H2O Pi, H2O Po, Total P, NaHCO3 Po, and HCl Po in the monoculture but not significantly related to any P fractions in the mixed culture. Ca 2+ was found to be significantly positively related to H2O Pi, H2O Po, NaHCO3 Pi, NaOH Po, and total P in the monoculture, but in the mixed culture, it was found to be significantly related to H2O Pi, H2O Po, NaHCO3 Pi-Po, and total P. Looking at the trait interrelations between the base cations, K + was found to be significantly positively related to Mg 2+ and Na + only in the monoculture and, in the mixed culture, significantly related to Mg 2+ and Ca 2+ . The BS% was significantly related to NaOH Pi, HCl Pi, and Ca 2+ in the monoculture and, in the mixed culture, significantly related to NaOH Pi and Ca 2+ . Interestingly, it was found that both the cultures showed a more or less similar pattern of correlations to their trait's interrelations. See Supplementary Figures S2 and S3 for details. Heat map for Pearson's correlation coefficients of soil acidification, P fractions, and other soil chemical properties in the mono-and mixed cultures at three replications. The P fractions and other chemical properties were identified in each sample by colours deduced from red to blue (−1.0 to +1.0) and the cross indicates not significant differences between parameters. Abbreviations: phosphorus (Pi and Po) fractions (extracted by H2O, NaHCO3, NaOH, and HCl); residual P (Resi P); hydrogen ions (H + ); potassium ions (K + ); magnesium ions (Mg 2+ ); manganese ions (Mn 2+ ); sodium ions (Na + ); calcium ions (Ca 2+ ); aluminum ions (Al 3+ ); iron ions (Fe 3+ ); effective cation exchange capacity (CEC); base saturation percentage (BS%); effective acidity (EA); aluminum calcium ratio (Al:Ca); ammonium (NH4 + -N); nitrate (NO3 + -N); total nitrogen (TN); soil organic carbon (SOC); carbon nitrogen ratio (C:N ratio); soil pH (pH); and total phosphorus (TP). Similarly, Mg was significantly positively related to NaHCO 3 Pi and HCl Po in the monoculture. In contrast, in the mixed, it was found to be significantly related to NaOH Po. Interestingly, Na was significantly positively related to H 2 O Pi, H 2 O Po, Total P, NaHCO 3 Po, and HCl Po in the monoculture but not significantly related to any P fractions in the mixed culture. Ca 2+ was found to be significantly positively related to H 2 O Pi, H 2 O Po, NaHCO 3 Pi, NaOH Po, and total P in the monoculture, but in the mixed culture, it was found to be significantly related to H 2 O Pi, H 2 O Po, NaHCO 3 Pi-Po, and total P. Looking at the trait interrelations between the base cations, K + was found to be significantly positively related to Mg 2+ and Na + only in the monoculture and, in the mixed culture, significantly related to Mg 2+ and Ca 2+ . The BS% was significantly related to NaOH Pi, HCl Pi, and Ca 2+ in the monoculture and, in the mixed culture, significantly related to NaOH Pi and Ca 2+ . Interestingly, it was found that both the cultures showed a more or less similar pattern of correlations to their trait's interrelations. See Supplementary Figures S2 and S3 for details.

Effects of Additional N on the Soil Chemical Properties
The study revealed that soil pH and BS% significantly decreased in N addition plots and increased the EA (Figure 2a,e). Our results are confirmed by the previous study of Shi et al. [58], which reported that the pH and BS% significantly decreased by the N addition in two contrasting warm temperate forests. Some further studies suggested that N is limited in temperate ecosystems, but due to an increase in N's atmospheric deposition, it can lead to N saturation [59,60]. The leaching down of nitrate due to its high mobility and weak bind to soil particles (i.e., negative charge on clay particle) results in a decrease in base cations and soil pH [61]. The increase in soil acidity might be the release of H + at the time of nitrification processes [62]. The pattern of the cation pool also indicates that all plantations might be approaching towards the Al-buffering stage reported in previous studies [60,63]. However, we did not measure the precipitation rate and chemical composition, but Al buffering in the study area might be because of two factors, i.e., precipitation and N deposition. Precipitation plays a vital role in soil biochemistry, as reported in a previous study that the mean precipitation rate would rise in the 21st century [64]. Secondly, in nature, a N wet deposition never occurs alone; it will occur with precipitation [25]. We also observed a significant increase in the Al:Ca ratio, while a significant decrease in BS% (Figure 2b,d). Our findings confirmed the previous findings of Mao et al. [24] and Malek [12], which suggests that these plantations are probable in Al stress; in this regard, further long-term study is needed.
We also observed that NO 3 − -N significantly increased while NH 4 + -N remained stable in N additional plots (Table 4 and Supplementary Table S2). Yan et al. [65] reported that, in Eastern China's terrestrial ecosystem, the N's input and its forms create a balance between H + ion consumption and production and plays a significant role in soil N mineralization with increased N application. However, the nonsignificant effect of NH 4 + -N and increase of NO 3 − -N also indicate that the decrease in ammonification and increase in nitrification suggest N saturation [66]. Similarly, the long-term N addition increases the N mineralization rate in the soil due to microorganism adaptability to increased soil acidification and the N content [67]. We found an increase in soil organic carbon (SOC) in N additional plots, which confirmed the findings of Tafazoli et al. [66] and Nave et al. [15], Since the global C cycle is affected by additional atmospheric N and essential to monitor the forest soil C:N ratio, it plays a significant role in forest health, as the plant uptake of available N is made through the C and microbial activities in soil [68]. The decrease in base cations like Ca 2+ and Mg 2+ in N additional plots (Supplementary Table S2) is attributed to the increase in Fe 3+ and Al 3+ ions, mainly responsible for soil acidification, and leads to the depletion of soil nutrients, which ultimately affect the vitality and health of forests [12,69,70]. Lucas et al. [9] reported that the emission of NOx and SO 2 cause soil acidification, resulting in weathering and leaching down of base cations. The selective weathering of certain base cations like Ca > Na > Mg > K would cause a reduction in base cation budgets and an imbalance of metal ions in soils [71].

Effects of Additional P on the Soil Chemical Properties
In our study, we observed that the H 2 O Pi and NaHCO 3 Pi significantly increased with the additional P (Figure 3a,b), while the rest of the P fractions remained stable in all plantations. The stability in other fractions might be refer to adding P, which increased the N uptake by plants and subsequently reduced the excess N from the ecosystem. This finding is an agreement with previous research, wherein P significantly increased the litter N and microbial biomass and reduced N mineralization [72], which reduced the nitrate leach-down effect and the emission of nitrous oxide [73]. In the acidification mitigation, nitrate uptakes by trees also provide hydroxyl ion in the soil [62]. An increase in H 2 O Pi and NaHCO 3 Pi also suggests that P biotic assimilation might be augmented, i.e., an increase in soil hydroxyl ions [23,62]. Our current study was also confirmed by the previous study of Zarif et al. (unpublished data) that the P increase in leaf P contents, and an increase in the microbial biomass carbon was also in-line with previous studies conducted in the temperate and tropical ecosystems [13,74]. Mao et al. [24] also reported that the amount of N was relatively higher in the winter than summer, which might depict that the continuous N deposition from the atmosphere increased the soil acidity, which indicates that a buffering effect by the additional P keeps the soil pH stable [13,24].
The exchangeable soil cations like K + , Ca 2+ , and Mg 2+ were not significantly affected by the P application (Table 4 and Supplementary Table S2). Our results align with the previous study of Mori et al. [62], who found that the response of exchangeable cations and base saturation percentage to the increased P deposition disappeared after a long term of P addition, contrary to Zhu et al. [75], who reported that the two-year P fertilization significantly increased the total soil exchangeable base cations by 60%. The Al concentration was also not significantly affected by the P fertilization ( Figure 1c). P's nonsignificant effect on Al might be exchangeable with Al's concentration on the soil pH [76]. The soil also maintains constant Al 3+ ions, resulting in a balance in soil acidity immobilization that might occur due to phosphate availability [24].
Last, we hypothesized that the combination of N and P might increase the process of acidification, i.e., a decrease in soil pH, increase in effective acidity (EA), and a reduction in BS. P added to the plots is the main evidence of the buffering effect. The EA in NP additional plots showed a significant increase (Figure 2a), which suggests that the N application might significantly increase the nitrification process. This also describes Al and Fe's solubilization, which decreased P binding considerably [13,77,78]. Similarly, Guo et al. [79] reported that the soil acidification in a major cropland of China was influenced by the N cycle during 1980-2000. Another research spread over seven plantations of Australia found that the combined applications of N and P fertilizers significantly increased the soil acidification in topsoil at five experimental sites [80]. However, more comprehensive monitoring would be necessary to understand P's role in buffering N-deprived soil acidification in temperate ecosystems.

Effects of Fertilizations on P Fractions
The total P is the sum of organic and inorganic P, which relates to the net input (litter) and net output (nutrients) of the plant uptake. In our results, we observed that the total P was not significantly affected by additional N and P fertilizations (Supplementary Figure S1f), which reflects that there might be balance in P's input and output by below-and aboveground [13]. However, we found a significant decrease in the NaHCO 3 Pi concentration under N application and a significant increase under the P addition, while no considerable effect was observed under NP application in all plantations (Figure 3b). Our findings are not in-line with the previous findings of Weand et al. [81], who stated that the abundance of available P and biotic P varied from species to species and remained stable with fertilizer application. The N addition could stimulate the potential net N mineralization rate and possible net nitrification processes [81], resulting in an increase in the P biotic demand. Similarly, the study conducted in a larch plantation revealed that N and P concentrations might be increased in the roots, stem, branches, and leaves [82]. However, we did not measure any root traits, but the decrease in NaHCO 3 Pi concentration under the N additional plot depicted that the plants consumed for growth and similar findings were also reported by Yang et al. [13] and that the root biomass of the Larch plantation was significantly increased by 20% after the long-term addition (nine years) of N fertilization. In contrast, the increase in H 2 O Pi and NaHCO 3 Pi (Figure 3a,b) showed that the P biotic assimilation might be increased under P addition by providing more hydroxyl ions in the soil, which justifies the previous results of Zarif et al. (unpublished data) and Mao et al. [24], who reported an increase in P contents in leaves. The decrease in NaHCO 3 Pi in the N additional plots could be due to the microbial activity. Khan et al. [83] reported a linear relationship between NaHCO 3 Pi (Olsen-P) and the microbial biomass under different management activities. However, our results contradict the previous findings of Sardans et al. [84], who suggested that the P uptake decreased in the evergreen Mediterranean forest due to the decrease of soil available P. Furthermore, we observed that the overall P inorganic concentration in the monoculture showed significant variations (i.e., Larch > Ash > Walnut) to that of the mixed culture, which suggests a decrease in the microbial and enzymatic activities.
The amount of labile Pi indicates that the mycorrhizae play an essential role in the nutrient uptake, especially under a stressed condition. In accordance to the findings of Cross and Schlesinger [85] and Yang et al. [13], who reported that the intermediate P fraction, i.e., NaOH, was not readily available to plants, made chalets with Fe and Al oxides and may represent P desorption sites. However, this study observed that the increase in NaOH Pi under N deposition attributed the significant increase in the nitrification process, i.e., to soil acidification by the microbial mediation process of ammonia oxidation H + ion. Furthermore, the soil pH decreasing by N addition plots while remaining stable in other plots suggests that the increase in Al and Fe solubilization significantly increased the soil P binding [77,78]. However, we observed that the N deposition might reduce the available P in soil by changing the inorganic P solubilization, P uptake, and microbial and enzymatic activities. Our study further confirmed the findings of Liu et al. [86], and Yang et al. [13], who reported that the long-term N application decreased the NaOH Po concentration of Pinus radiata and Larix gmelinii in subtropical and temperate forests. It was suggested that the NaOH Po pool can be a significant sink for available P in soil and may maintain the levels of plant available P through P mineralization (Yang et al. [13] and Beck and Sanchez [87]).

Conclusions
Our results indicate a significant increase in effective acidity under N and NP fertilizations but were not significantly affected by P fertilization in all plantations. Similarly, under N addition, the significant decrease in pH, base cations, Ca 2+ , and Mg 2+ was attributed to the increase in Fe 3+ and Al 3+ ions, mainly responsible for soil acidification and leading to the depletion of soil nutrients, which, ultimately, affects the vitality and health of forests. Similarly, NaHCO 3 Pi and NaOH Pi pointed out that N deposition may stimulate plant P uptake. In contrast, the increase in P deposition plots of H 2 O Pi and NaHCO 3 Pi depicts that P biotic assimilation might be increased. The study further concluded that, although the P deposition showed a buffering effect, it could not help to mitigate the soil acidity. The P concentration was ranked as Larch >Ash > Walnut, respectively, in the monoculture, while, in the mixed, it more or less remained the same. This study gave detailed insight and understanding about the potential factors affecting P fractions and soil acidification in a temperate ecosystem. However, we think it will be necessary to conduct long-term research on soil aggregates in relation to fertilization, atmospheric N deposition, and precipitation rates of different forest stands to assess the generality of our results.

Supplementary Materials:
The following are available online at http://www.mdpi.com/1999-4907/11/12/1274/s1, Figure S1. The organic P fractions (extracted by H 2 O, NaHCO 3 , and HCl) Residual and P Total P. Figure S2. Heat map for Pearson's correlation coefficients of soil acidification, P fractions, and other soil chemical properties in mono cultures at three replications. Figure S3. Heat map for Pearson's correlation coefficients of soil acidification, P fractions, and other soil chemical properties in mixed cultures at three replications. Table S1. ANOVA results of species, fertilizations and their interactions in mono and mixed cultures. Table S2. Soil chemical properties of mono culture and mixed culture.