Soil carbon (C), nitrogen (N) and phosphorus (P) stoichiometry drives phosphorus lability in paddy soil under long-term fertilization: A fractionation and path analysis study

Abstract

Soil C:N:P stoichiometry plays a vital role in nutrient cycling in ecosystems, but its importance to P transformation in paddy soil remains unclear. We investigated the effect of soil C:N:P stoichiometry on P mobility and uptake under long-term fertilization. Three treatments, CK (no fertilization), NPK (inorganic nitrogen, phosphorus and potassium fertilization) and NPKM (combined inorganic NPK fertilizer and manure application), were selected from two long-term experiments of paddy soil that were initiated in 1991 and 1982 in Chongqing and Suining, respectively. The results showed that in comparison the control treatment, under long-term fertilization, soil pH decreased. In comparison with the NPK and CK treatments, the NPKM treatment significantly increased soil nutrient contents, P uptake and phosphatase activities. In comparison to the CK treatment, the NPK and NPKM treatments significantly decreased soil C:N, C:P and N:P ratios. In comparison to NPK and CK treatments, the NPKM treatment decreased residual-P at both sites. Compared with CK treatment, the NPKM treatments increased labile-P and moderately labile-P by 987% and 144%, respectively, and NPK treatment increased these factors by 823% and 125%, respectively, at the Chongqing site. At the Suining site, with NPKM treatment, increases in labile-P and moderately labile-P were 706% and 73%, respectively, and with NPK treatment, the increases were 529% and 47%, respectively. In contrast, non-labile-P was significantly decreased with NPKM treatment in comparison to that with NPK and CK treatments. Moreover, increases in soil C:N and C:P ratios decreased the labile-P pools and increased non-labile-P pools. A path analysis indicated that soil C:N:P stoichiometry indirectly controlled P uptake by directly affecting P transformation from non-labile to labile-P pools. Moreover, the non-labile-P in soil with high SOM and P content directly affected P uptake, indicating that soil P transformation is mainly driven by soil C and P in paddy soil. In conclusion, understanding mechanism of P mobility influenced by soil C:N:P stoichiometry could be helpful to manage soil P fertility under long-term fertilization in paddy soils of these regions.

About the Authors:

Muhammad Qaswar

Roles Conceptualization, Data curation, Investigation, Methodology, Software, Writing – original draft, Writing – review & editing

Affiliation: National Engineering Laboratory for Improving Quality of Arable Land, Institute of Agricultural Resources and Regional Planning, Chinese Academy of Agricultural Sciences, Beijing, China

Waqas Ahmed

Roles Formal analysis, Writing – review & editing

Affiliation: National Engineering Laboratory for Improving Quality of Arable Land, Institute of Agricultural Resources and Regional Planning, Chinese Academy of Agricultural Sciences, Beijing, China

Huang Jing

Roles Data curation, Resources

Affiliations National Engineering Laboratory for Improving Quality of Arable Land, Institute of Agricultural Resources and Regional Planning, Chinese Academy of Agricultural Sciences, Beijing, China, National Observation Station of Qiyang Agri-ecology System, Institute of Agricultural Resources and Regional Planning, Chinese Academy of Agricultural Sciences, Qiyang, Hunan, China

Fan Hongzhu

Roles Data curation, Resources

Affiliation: Soil and Fertilizer Research Institute, Sichuan Academy of Agricultural Sciences, Chengdu, China

Shi Xiaojun

Roles Resources

Affiliation: College of Resources and Environment, Southwest University, Chongqing, China

Jiang Xianjun

Roles Resources

Affiliation: College of Resources and Environment, Southwest University, Chongqing, China

Liu Kailou

Roles Resources

Affiliations National Engineering Laboratory for Improving Quality of Arable Land, Institute of Agricultural Resources and Regional Planning, Chinese Academy of Agricultural Sciences, Beijing, China, Jiangxi Institute of Red Soil, National Engineering and Technology Research Center for Red Soil Improvement, Nanchang, China

Xu Yongmei

Roles Resources

Affiliation: Institute of Soil, Fertilizer and Agricultural Water Conservation, Xinjiang Academy of Agricultural Sciences, Urumqi, China

He Zhongqun

Roles Formal analysis, Resources

Affiliation: College of Horticulture, Sichuan Agricultural University, Chengdu, China

Waleed Asghar

Roles Formal analysis

Affiliation: School of Environment, Beijing Normal University, Beijing, China

Asad Shah

Roles Data curation, Formal analysis

Affiliation: National Engineering Laboratory for Improving Quality of Arable Land, Institute of Agricultural Resources and Regional Planning, Chinese Academy of Agricultural Sciences, Beijing, China

Huimin Zhang

Roles Funding acquisition, Project administration, Resources, Supervision

* E-mail: [email protected]

Affiliations National Engineering Laboratory for Improving Quality of Arable Land, Institute of Agricultural Resources and Regional Planning, Chinese Academy of Agricultural Sciences, Beijing, China, National Observation Station of Qiyang Agri-ecology System, Institute of Agricultural Resources and Regional Planning, Chinese Academy of Agricultural Sciences, Qiyang, Hunan, China

ORCID logo http://orcid.org/0000-0002-8509-731X

Introduction

Crop productivity is directly affected by P availability in agricultural soils [1]. However, in some agricultural soils, the concentration of total P is naturally high, which can limit plant growth due to the low solubility and rapid conversion of P compounds to unavailable or poorly available P after fertilizer application [2]. A synthetic chemical fertilizer application can ensure sufficient availability of P for plant uptake [3]. However, there has been growing confirmation that long-term excessive application of inorganic P fertilizers decreases P use efficiency that may result in excessive accumulation of P in soil [4].Specifically, in intensive cultivation areas, during surface runoff or erosion the P-rich soil particles may reach surface water and cause eutrophication [5,6].

The excessive use of inorganic fertilizers could be reduced by animal and green manure application in fields [7,8]. Organic amendments can increase the P recovery in soil by decreasing the formation of irreversible adsorption complexes in soil [9,10] or by increasing the P mobility in soil [11]. Additionally, organic amendments not only provide macronutrients and micronutrients, but also improve soil physical and chemical properties by increasing SOM [12]. However, in practice the rate of manure addition to fields is mostly calculated based on N input rates, which may cause excessive P addition with manure, especially under long-term fertilization [13]. Therefore, the combined application of inorganic fertilizers and manure could be beneficial to increasing the efficiency of applied P and thus improving crop productivity. The arable land in the Chongqing and the Suining (southern China) are mainly paddy soils, which are mostly cultivated under rice based cropping systems [14,15]. Paddy soils, especially in the southern regions of China, are inherently low in available P due to acidification. Therefore, to meet plant P requirements, excessive inorganic fertilization is a common practice in these regions, which causes environmental P losses and degraded soil fertility [16]. In these regions, previous studies mainly focused on long-term fertilization effect on crop yields and soil fertility [14,17], and recommended, integrated nutrient management through the combination of applying inorganic fertilizer and manure addition to improve soil P fertility [18]. Most existing studies have clearly shown that the combined application of inorganic and manure increases the availability of P and crop yield [19], but the mechanism that causes the redistribution of residual-P under long-term fertilization in paddy soils is not well studied [20]. There could be many factors regulating P in paddy soil under long-term fertilization. In a previous study, Ahmed et al. [21] found that the combined application of manure and inorganic fertilizers increased the labile-P fraction, and different soil properties (pH, SOM, total P and Olsen-P) showed different relationships with P fractions, however, the authors did not explain the mechanism of P pools transformation, mobility and uptake in paddy soils associated with soil nutrient stoichiometry. One of the most important factors for P mobility in paddy soil could be soil nutrient stoichiometry, coupling and interactions of the C, N and P cycles [22–24]. In recent years, great progress related to C:N:P stoichiometry in terrestrial ecosystems with focuses on microorganisms, plant leaves and litter, has been made [25,26]. Therefore, C:N:P stoichiometry could be used as a powerful tool to understand nutrient cycling and processes in soil [27]. Tian et al. [28] also reported that soil C:N:P stoichiometry could be a potential indicator for assessing soil nutrient status during soil development. However, the impact of soil C:N:P stoichiometry on P cycling, especially in paddy soil, is poorly understood.

The phosphorus that is applied to soil goes through many complex transformations (desorption-adsorption and dissolution-precipitation) controlled by soil properties such as pH, SOM and exchangeable iron-aluminum (Fe-Al) hydroxides [29,30]. The sequential procedure of P extraction developed by Hedley et al. [31], extracts the different fractions of P that is bound with complex organic and inorganic compounds using solutions of different strengths for the P extraction [32]. The Hedley method of P fractionation is one of the comprehensive methods that is used to investigate the dynamics of P lability and the P cycle in soil [33]. However, the sequential procedure of P fractionation alone does not provide insight into the transformation of P pools affected by different treatments and soil factors [34]. Path analysis could be used as a valuable tool to investigate the interrelationship between soil P pools in different soils and management practices [35]. Therefore, the aim of this study was to investigate the influence of soil nutrient stoichiometry on P pool transformation and P mobility in paddy soil after long-term fertilization using a structural equation modeling (SEM) pathway technique.

Materials and methods

Experimental design and site description

In this study, we selected two sites from long-term experiments located in Chongqing (106°24’33”E, 29°48’36”N) and Suining (105°03’26”E, 30°10’50”N) that were initiated in 1991 and 1981, respectively. These experimental sites belong to national public welfare research institutions. Research at these experiment sites was carried out with the permission of the Institute of Agricultural Resources and Agricultural Regional Planning, Chinese Academy of Agricultural Sciences, as cooperative research agreements were signed with these experimental sites. Climate at the Chongqing site is subtropical, and the climate at the Suining site is subtropical humid. The initial soil physical and chemical properties and climatic conditions of the experimental sites are given in Table 1. Soils in both sites are Purple-Udic Cambisols according to Chinese soil classification [36]. The cropping systems at both sites were rice-wheat double cropping systems. Hybrid indica rice was planted in the rice season at both sites, and we used data from only the rice season. Treatments at both sites were arranged in a randomized complete block design (RCBD) with three replicates. Each plot was separated by permanent cemented ridges. In this study, the three selected treatments were common at both sites. The treatments included no fertilization (CK), inorganic fertilization (NPK) and inorganic fertilizer plus manure (NPKM). The rates of fertilization are shown in Table 2. The chemical fertilizers applied at both sites were urea for N, calcium superphosphate for P and potassium sulfate for K. In the NPKM treatment, a half dose of the fertilizers was applied from an inorganic fertilizes, and the remaining half dose of nutrients was supplemented with pig fecal urine at the Chongqing site and fresh pig manure at the Suining site based on N content. The remaining P and K application rates were adjusted by inorganic fertilizers. In both experiments, phosphorus fertilizer and manure were applied as a basal fertilization before transplantation of the seedlings, and half of the N and K fertilizers were applied as a basal application. The remaining half were dressed in two splits at the tillering and panicle initiation stages of the crop. The nutrient content of the pig manure was 0.2% N, 0.17% P2O5, and 0.1% K2O, and the nutrient content of the pig fecal urine was 1.1% N, 0.9% P2O5, and 0.4% K2O experiment. The rice seedlings were transplanted in the middle of May, and the crop was harvested at the end of August at full maturity.

[Figure omitted. See PDF.]

Table 1. Experiment initiation date, climate and initial soil properties at the two long-term experimental sites.

https://doi.org/10.1371/journal.pone.0218195.t001

[Figure omitted. See PDF.]

Table 2. Annual fertilizer input rates (kg ha−1) for the crop at the long-term experimental sites of two typical croplands in China.

https://doi.org/10.1371/journal.pone.0218195.t002

Sampling and chemical analysis

After crop harvest, plant samples were air dried, and the yield was measured and the total P content in the plant samples was measured after oven drying at 60°C to calculate the P uptake. Soil samples were collected in the first week after the rice harvest in 2017. Composite soil samples were collected at a 0–20 cm soil depth from four randomly selected locations in each replicate block of treatments. Subsamples from each composite sample were stored at 4°C to prevent the moisture losses and to determine phosphatase activities within a week of sampling. Apart from the composite samples, the other samples were air dried and sieved < 2 mm to determine the soil chemical properties. Soil pH was determined in suspensions of water: soil (2.5: 1) [37]. Conventional routine methods [38] were used to analyze the soil organic carbon (SOC), total nitrogen (TN) and total phosphorus (TP). The wet oxidation method was used to determine SOC [39], and the Kjeldahl digestion-distillation method was used to determine soil total N. Soil total P was measured by the molybdenum blue colorimetric method using a HClO4-H2SO4 solution for digestion. Mineral N, also known as alkali hydrolysable N, was determined according to Lu [37], and available P was measured according to Olsen et al. [40]. Except for the refrigerated sample, the soil moisture contents of the samples were measured by oven drying at 105°C for 24 h before processing the soil samples for enzyme activities. The activities of acid phosphomonoesterase (AcP) and phosphodiesterase (PD) were determined according to Tabatabai [41]. The substrate solution used for AcP was p-nitrophenyl phosphate and the pH of the buffer solution was 6.5. The substrate solution used for PD was bis-p-nitrophenyl phosphate, and the pH of the buffer solution was 8. The activities of Acp and PD were expressed as mg p-nitrophenol kg−1 soil h−1.

Phosphorus fractionation

Organic and inorganic fractions of soil phosphorus were determined using the Hedley et al. [31] fractionation method modified by Tiessen and Moir [42]. Briefly, the soil samples were divided into three replicates each of 1 g soil and sequential extraction procedure was carried out as follows: NaHCO3-P was extracted by continuous shaking for 16 h after adding 0.5M NaHCO3 at pH 8.5 in 1 g soil. Then, shaking continued with 0.1M NaOH for 16 h to extract NaOH-P in solution. Next, shaking was continued with 1 M HCL for 16 consecutive hours to obtain diluted HCl-P and then heated with 10 ml of concentrated HCl at 80°C in a water bath for 10 minutes after adding 5 ml 12M HCl, a final volume of 50 ml was made with deionized water to obtain concentrated HCl-P (conc. HCl-P). Finally, to obtain residual-P, the residual soil samples were mineralized by adding 300 μl of concentrated H2SO4 in each soil sample of 30 mg followed by heating at 350°C for 3 h. The tubes were centrifuged at 25,000×g and 4°C for 10 minutes between two consecutive steps. The supernatant was passed through 0.45 μm cellulose nitrate filters, and extra soil particles were recovered by washing filters with extracted solution used in the following steps. Both organic and inorganic P (by taking the difference between total P after digestion with persulfate and inorganic P by colorimetry) levels were determined by measuring with 0.5 M NaHCO3, 0.1 NaOH and 12 M HCl extracts.

Statistical analysis

The soil C:N, C:P and N:P ratios were calculated from SOC, total N and total P. We used one-way ANOVA followed by Tukey’s HSD test at the P = 0.05 level for each variable to evaluate the significant differences among the different treatment means [43]. The main effects of the fertilizer treatments, sites and their interactions were evaluated by two-way ANOVA. To determine the relationship between soil C:N:P stoichiometry and the P pools, linear regression analysis was performed using SPSS 16.0 software (SPSS, Chicago, IL, USA). Redundancy analysis (RDA) and variance partitioning analysis were performed using Canoco (version 5.0) software. Structure equation modeling (SEM) is a statistical procedure to determine the direct and indirect effects and estimates the partial contribution of explanatory variables [44,45]. We used SEM to explore the hypothesized linkage among C:N:P stoichiometry, P pools and P uptake using the AMOS package with SPSS 16.0 software (SPSS, Chicago, IL, USA).

Results

Effect of long-term fertilization on soil properties and soil C:N:P stoichiometry

The interaction of soil × fertilization significantly (P ≤ 0.05) affected soil properties and soil C:N:P stoichiometry except soil pH and mineral N. However, the main factors (site and fertilizer treatments) significantly (P ≤ 0.05) influenced pH and mineral N. Among both soils, soil pH ranged from 7.8 to 7.6 (Fig 1A). At the Suining, the pH in the control (CK) was significantly higher than the pH in the fertilizer treatments (NPK and NPKM). At both sites, in comparison to the CK and NPK treatments, the combined application of inorganic fertilizer plus manure (NPKM) significantly (P ≤ 0.05) increased soil organic matter (SOM), total N, total P, mineral N and Olsen-P. Soil total P at the NPKM treatment was more than double that in NPK and CK treatments at both sites. Compared with the CK treatment, in the NPKM treatment increased SOM, total N, total P, mineral N and Olsen-P by 74.4, 102, 343, 29 and 1307% at the Chongqing site and by 34.4, 104, 532, 57 and 1194% respectively at the Suining. Compared with the CK treatment, the NPKM treatment increased SOM, total N, total P, mineral N and Olsen-P by 34, 60.4, 42, 19, 965% respectively, at the Chongqing site and by 18, 82, 62, 21 and 597% respectively at the Suining site.

[Figure omitted. See PDF.]

Fig 1.

Effect of long-term fertilization on soil pH (A), organic matter (B), total nitrogen (C), total phosphorus (D), mineral nitrogen (E) and Olsen-P (F) in paddy soils. Error bars represent ± standard deviations; different lower-case letters over the bars indicate significant (P ≤ 0.05) differences according to Tukey's HSD test.

https://doi.org/10.1371/journal.pone.0218195.g001

Changes in soil total nutrient contents under long-term fertilization changed soil C:N:P stoichiometry in both soils. In comparison to no fertilizer and the inorganic fertilizer application, the long-term combined application of manure and inorganic fertilizer decreased the soil C:N, C:P and N:P ratios (Fig 2). However, the highest C:N and C:P ratio was in the CK treatment, and the highest N:P ratio was in NPK treatment at both sites.

[Figure omitted. See PDF.]

Fig 2.

Effect of long-term fertilization on soil C:N (A), C:P (B) and N:P (C) ratios in paddy soils. Error bars represent ± standard deviations; different lower-case letters over the bars indicate significant (P ≤ 0.05) differences according to Tukey's HSD test.

https://doi.org/10.1371/journal.pone.0218195.g002

Effect of long-term fertilization on phosphatase activities, phosphorus fractions and uptake

Interactive effect of fertilization × site significantly (P ≤ 0.05) influenced acid phosphomonoesterase (AcP) and phosphodiesterase (PD) activities (Fig 3). The activities of AcP were not significantly different between the CK and NPK treatments at both sites. Compared with the CK treatment, the NPKM and NPK treatments increased the AcP by 161 and 24%, respectively at the Chongqing site and 54 and 6.4%, respectively in Suining. The increases in PD with NPKM and NPK treatments compared the CK treatment were 294 and 121% at Chongqing site and by 75 and 69%, respectively at the Suining site. Averaged across the treatments, the activity of AcP at Chongqing was by 23% greater than that at Suining, and the activity of PD at the Chongqing by 32% lower than that at the Suining site.

[Figure omitted. See PDF.]

Fig 3.

Effect of long-term fertilization on acid phosphomonoesterase (A) and phosphodiesterase (B) activities in paddy soils. Error bars represent ± standard deviations; different lower-case letters over the bars indicate significant (P ≤ 0.05) differences according to Tukey's HSD test.

https://doi.org/10.1371/journal.pone.0218195.g003

The soil P fraction distribution and uptake were significantly (P ≤ 0.05) influenced by the site × fertilizer interaction. Regardless of the treatments, all organic and inorganic fractions at the Chongqing were greater than those at the Suining site except NaOH-Pi, HCl dil Pi and residual-Pi (Table 3). At the Chongqing site, compared with the CK treatment, the combined application of inorganic fertilizer and manure treatment increased NaHCO3-Pi, NaHCO3-Po, NaOH-Pi, NaOH-Po, HCl dil-Pi, HCl conc Pi and HCl conc-P fractions by 1481, 368, 16706, 244, 106, 31 and 81% respectively and the NPK increased these factors by 1322, 197, 11798, 225, 96, 29.5 and 86%, respectively. At the Suining site, compared with the CK treatment, the NPKM treatment increased the NaHCO3-Pi, NaHCO3-Po, NaOH-Pi, NaOH-Po, HCl dil-Pi, fractions by 888, 394, 5292, 77 and 37%, respectively and with the NPK treatment increased these fractions by 7.6, 1.4, 31, 0.3 and 0.28%. Residual-P was significantly decreased in the NPK and NPKM treatments compared with the CK treatment at both sites. In comparison with the NPK and CK treatments, long-term inorganic fertilizer plus manure application increased size of labile and moderately labile, while the non-labile pool was decreased in NPKM and NPK treatment compared with that in CK treatment (Fig 4).

[Figure omitted. See PDF.]

Fig 4.

Effect of long-term fertilization on labile-P (A), moderately labile-P and non-labile-P pools in paddy soils. Error bars represent ± standard deviations; different lower-case letters over the bars indicate significant (P ≤ 0.05) differences according to Tukey's HSD test.

https://doi.org/10.1371/journal.pone.0218195.g004

[Figure omitted. See PDF.]

Table 3. Concentrations of different fractions of phosphorus in the paddy soils affected by long-term fertilization at the two experimental sites.

https://doi.org/10.1371/journal.pone.0218195.t003

In comparison with the NPK and CK treatments, the long-term inorganic plus manure application treatment significantly (P ≤ 0.05) increased P uptake (Fig 5) by increasing P lability. Compared with the CK treatment, the NPKM and NPK treatments increased P uptake by 345 and 287%, respectively, at the Chongqing site and by 246 and 204%, respectively at the Suining site. Regardless of the treatments, the total P uptake at the Chongqing site was 46% greater than that at the Suining site.

[Figure omitted. See PDF.]

Fig 5. Effect of long-term fertilization on phosphorus uptake in rice crop in paddy soils.

Error bars represent ± standard deviations; different lower-case letters over the bars indicate significant (P ≤ 0.05) differences according to Tukey's HSD test.

https://doi.org/10.1371/journal.pone.0218195.g005

Relationship between soil C:N:P stoichiometry and phosphorus mobility

The results from the linear regression showed that the soil C:N, C:P and N:P ratios were negatively correlated with labile-P and moderately labile-P pools (Fig 6). The C:N ratio had a highly significant (P ≤ 0.05) positive correlation (R2 = 0.73) with non-labile-P. The soil C:P and N:P ratios had a weak correlation with soil non-labile-P pool. Redundancy analysis (RDA) showed the effect of soil C:N:P stoichiometry on P fractions (Fig 7). Non-labile-P fractions especially residual-P, showed a close relationship with the N:P and C:P ratios, while, HCl conc-Pi and HCl conc-Po showed a strong positive relation with the C:N ratio. RDA-1 explained 55% of the total variation, and RDA-2 explained 8.6% of the total variation.

[Figure omitted. See PDF.]

Fig 6. Linear regression indicating the relationship between soil C:N:P stoichiometry and phosphorus pools affected by long-term fertilization in paddy soil.

https://doi.org/10.1371/journal.pone.0218195.g006

[Figure omitted. See PDF.]

Fig 7. Redundancy analysis of the correlations between the soil C:N:P ratio, phosphorus fractions and phosphatase activities.

Red lines indicate the effect of C:N, C:P and N:P ratios on phosphorus fractions and phosphatase enzyme activities. Blue dashed lines indicate the explanatory variables. Abbreviations: AcP: acid phosphomonoesterase activity; DP: phosphodiesterase activity.

https://doi.org/10.1371/journal.pone.0218195.g007

The variance portioning analysis (VPA) showed that soil properties, soil C:N:P stoichiometry and phosphatase totally explained 97.5% of the variation in the soil P fraction distribution (Fig 8). Soil properties included soil pH, SOM, total N, total P, and soil C:N:P stoichiometry (C:N, C:P and N:P ratios), and phosphatase activities included AcP and PD. Out of total variance the soil properties, soil C:N:P stoichiometry and phosphatase explained 30.3, 4.9 and 4.6% of the total variation, respectively, indicating the relative dominance of soil properties as a major factor.

[Figure omitted. See PDF.]

Fig 8.

The variance partitioning analysis showing the proportional contribution (%) of soil properties (S), soil C:N:P stoichiometry (C), phosphatase (F) and their interactions on phosphorus fractions.

https://doi.org/10.1371/journal.pone.0218195.g008

We also investigated the direct, indirect and total effects of soil C:N:P stoichiometry on P pools and P uptake by a SEM pathway analysis (Fig 9A and 9B). The results showed that at the Chongqing site, SEM explained 99, 95, 90 and 99% of the variance in labile-P, moderately labile-P, non-labile-P and P uptake, respectively. N:P stoichiometry explained that the C:N and C:P ratios directly affected the non-labile P pool and indirectly affected the moderately labile and labile-P pools (Fig 9A). However, the C:N, C:P and N:P ratios indirectly affected P uptake by directly influencing P pools. Non-labile and labile-P both directly influenced P uptake, while moderately labile-P indirectly affected P uptake by directly affecting labile-P. At the Suining site, SEM explained 99, 99, 89 and 99% of the variance in labile-P, moderately labile-P, non-labile-P and P uptake, respectively. C:N, C:P and N:P ratios directly affected non-labile and moderately labile-P and indirectly affected labile-P (Fig 9B). However, P uptake was directly affected only by labile-P. Non-labile had no significant effect on moderately labile-P, but moderately labile-P had a direct effect on labile-P and an indirect effect on P uptake.

[Figure omitted. See PDF.]

Fig 9.

Structural equation modeling results for the direct and indirect effects of soil C:N:P stoichiometry on P lability and uptake at the Chongqing site (A) and the Suining (B) site. The numbers adjacent to the arrows are standardized path coefficients, which are analogous to relative regression weights, and indicative of the effect size of the relationship. The continuous and dashed arrows indicate direct and indirect relationships, respectively. Gray arrows indicate a non-significant effect. The proportion of explained variance (R2) appears above every response variable in the model. The goodness-of-fit statistics for each model are shown in the lower right corner. Chi/df = 0.86. RMSEA: Root mean square error of approximation.

https://doi.org/10.1371/journal.pone.0218195.g009

Discussion

Effect of long-term fertilization on soil properties, soil C:N:P stoichiometry and enzyme activities

Long-term inorganic and organic fertilizer application significantly influence soil properties [46]. In this study, we found a significant (P ≤ 0.05) effect of long-term chemical fertilizer and manure application on soil chemical properties (Fig 1) that affected soil C:N:P stoichiometry and enzyme activities. In this study, soil pH was significantly decreased by long-term chemical fertilization. Our results are consistent with those of Chen et al. [47] and Cui et al. [48], who investigated the long-term effects of organic and inorganic fertilization on soil acidification and found that the soils receiving different rates of chemical fertilizer inputs showed accumulation of exchangeable acidity and a reduction in soil pH. Soil pH may also influenced by the removal of base cations from soil with long-term cultivation [49]. Some studies have suggested that long-term manure application may increase the soil pH due to alkalinity of manure [48], a result is consistent with our results. In comparison with no fertilization and inorganic fertilization, long-term fertilization, especially combined application of NPK plus manure, highly increased SOM (Fig 1B). These results were also supported by the results of previous studies under long-term fertilization [49]. As a long-term combined application of manure and inorganic fertilizer enhances the productivity, an increase in carbon input in soil is also expected because it is a key driver of increased SOM [50,51]. The results were also consistent with findings of Ladha et al. [52], who found that long-term fertilization increases the soil organic matter, especially in paddy soil, due to a slower decomposition of SOM under flooded conditions [53]. However, due to different soil properties, the SOM content in both soils was significantly different, which could have been due to different soil properties and different N input rates (Table 2). SOM and nutrient availability are directly affected by fertilizer management [54]. In our study, the long-term combined application of inorganic fertilizer and manure also significantly increased the total and available N and P in at both experimental sites (Fig 1). The increases in nutrient availability with a continuous manure addition were observed in previous studies [46,55,56], and this increase could be due to the contribution of nutrient reserves in SOM-rich soils [47].

In comparison with no fertilization, long-term addition of chemical fertilization and manure in paddy soils changed the soil C:N:P stoichiometry (Fig 2). Surplus P accumulation in paddy soil under long-term combined inorganic fertilizer and manure application decreased the soil C:N, C:P and N:P ratio in the NPKM treatment in comparison with the inorganic fertilization and no fertilization treatments. This result was also due to the low supply of N and P in the CK and NPK treatments. Previous studies have also suggested that land-use changes such as cultivation practices and long-term fertilization may produce soil nutrient stoichiometric shifts, especially in paddy soils [57]. In our study, in both soils, the C:N and C:P ratios in the topsoil (0–20 cm) were approximately similar to the average ratios for China (12.3 and 52.7, respectively) [28] and to global average ratios [27]. The N:P ratio in our results was lower than those found by Wojciech et al. [58] in a forestland and by Wang et al. [57] in subtropical wetlands in China. In another study, an N:P ratio similar to our results was reported by Zhang et al. [59] during the afforestation on loess Plateau of China. Regardless of the treatments, the overall lower N:P may be due to the higher solubility of N than of P and the consequently higher losses of N than of P due to continuous flooding [57]. In addition, a lower N:P ratio under long-term combined application of inorganic fertilizer and manure supports the hypothesis of a more productive ecosystem with high P accumulation and a lower N:P ratio [60].

Enzyme activities play a vital role in the biochemical reactions of C, N and P cycling [61–63] and have different responses to different type of fertilization, soil and environment [64]. The changes in soil nutrient availability and stoichiometry under different types of long-term fertilization changed phosphatase activities (Fig 3). In both soils, the highest activities AcP and PD were found in combined inorganic and manure application treatment (Fig 3), which might be due to long-term manure addition and high C inputs compared with those in the NPK and CK treatments. Nayak et al. [65] also suggested that phosphatase activities increased by increasing soil organic matter. In another study Ohm et al. [66] also found positive relationship between soil carbon and phosphatase activities. Therefore, studying phosphatase activities could be helpful to understand the lability of P under long-term fertilization regulated by soil C:N:P stoichiometry.

Phosphorus fractions, uptake and its relationship with soil C:N:P stoichiometry

Long-term fertilization significantly influenced P fractions distribution (Table 3) and P uptake (Fig 5). The changes in P fractions could be associated with soil properties. Many studies have explained the effect of soil pH, SOM, and P input rates on P fractions [67,68]. In our results, in contrast to previous studies, in comparison to no fertilization, long-term combined application of inorganic fertilizers and manure greatly increased the inorganic and organic fractions of P, except for the concentrated HCl organic and inorganic fractions of P (HCl conc-Pi and HCl conc-Po) and residual-P (Table 3). Yan [69] found that soil inorganic P fractions significantly increased following long-term inorganic fertilization but no corresponding differences were found in organic P fractions. In comparison to no fertilization, both the NPK and NPKM treatments decreased residual-P by increasing the labile-P. Because continuous manure application increased the lability of P by increasing both organic and inorganic fractions in labile pool which decrease the residual-P [20]. The P characteristics of organic material may also influence soil P content. For example, Li et al. [70] reported that in comparison to chemical fertilization, pig manure contained higher PKCl content, therefore, the major available proportion of P in pig manure could be H2O-P or NaHCO3-P, which is readily available P for P uptake [70] Similarly, in our results, in comparison with inorganic and no fertilization, long-term combined application of inorganic fertilizer plus manure significantly increased P uptake in rice crops more by increasing the labile-P pool due to long-term manure addition (Fig 4). Our results are also consistent with results from previous studies [71]. Dobermann et al. [72] also found that P application with or without manure increased the soluble inorganic P fractions with little or no effect on residual-P. Different studies have reported that, soil containing high organic C content showed less P sorption and retention capacities for P [73], which indicates the high P availability for uptake.

Soil nutrient stoichiometry is more complex than in plants in the response to nutrient availability, because nutrient mobility in soil is influenced by both fertilization and plant uptake. In our results, soil C:N:P stoichiometry influenced P pools in paddy soils (Figs 6 and 7). Many studies have reported poor P fertility in flooded soil conditions [69]. However, there could be many soil and environmental factors that control the P mobility in soil [74,75], but most important factor controlling P mobility in paddy soil could be soil C:N:P stoichiometry. Because under long-term fertilization, the increase in carbon inputs by long-term combined inorganic fertilizer and manure addition shift the nutrient cycling in cropland, the decoupling of soil C:N:P in previous studies is well documented [76], however, the effect of soil C:N:P stoichiometry on P lability specially in paddy soil has not been studied before. The negative relationship of C:N and C:P ratios with labile-P and the positive relationship of C:N and C:P ratios with non-labile-P (Fig 6) indicate that long-term manure addition increased carbon input [51] and increased the non-labile fractions of P in paddy soils. The non-labile-P fractions contained organic P with different characteristics [69,74,77], which could be increased by long-term manure addition. Lehmann et al. [78] also reported that organic amendments mainly increase the organic fractions of P. In another study, González Jiménez et al. [77] also found that organic soils contained relatively higher concentrations of labile and non-labile organic proportion of P. Moreover, some studies have explained that manure addition in paddy soil decrease the decomposition rate of organic matter, especially under flooded conditions [53], which may affect the P release mechanism [20]. Therefore, in the present study high C:N and C:P ratios showed negative relationship with labile-P.

The SEM showed that soil C:N:P stoichiometry indirectly affected the P uptake by controlling P transformation and mobility in paddy soils (Fig 9). In a previous study, González Jiménez et al. [77] explained the pathways of P pools. However, to our knowledge no study before explained the soil C:N:P stoichiometric effect on P mobility in paddy soils by using pathway analysis. The C:N and C:P ratios directly affected the non-labile-P, possibly due to high carbon input under long-term manure addition. Soil organic matter and soil total P could be the main factors affecting P mobility in paddy soils [77]. Labile-P had a direct effect on P uptake (Fig 9A and 9B), because a labile-P pool is the sum of the all readily available P fractions, which have a direct effect on P uptake [79,80].The indirect effect of the soil N:P ratio on P uptake could be due to a positive direct effect of the N:P ratio on moderately labile-P. Previous studies have explained the synergistic effect of N on P uptake [23,57,81]. N application also affects P uptake by increasing biomass production and soil carbon inputs [82]. Interestingly, at the Chongqing site, non-labile-P had a direct effect on P uptake, and at the Suining site, non-labile-P had an indirect effect on P uptake (Fig 9B), because the soil P at the Chongqing site was relatively higher than that at the Suining site, which directly affected P uptake. These results are consistent with the results of previous studies, for example, in a path analysis study by González Jiménez et al. [77], soil residual-P directly affected the labile-P and therefore, increased the P uptake. In this study we found that the soil C:N:P stoichiometry controlled P lability in paddy soil. Therefore, our results emphasized the significance of soil C:N:P stoichiometric effect on P mobility and uptake in paddy soil. SEM explained the unique relationship between soil C:N:P stoichiometry and P mobility in paddy soils affected by long-term fertilization, and it could be helpful to understand the mechanism of P transformation in paddy soil affected by soil C, N and P cycle interaction. It also could be helpful to manage soil P fertility by understanding the mechanism of P mobility under long-term fertilization in paddy soils of the Chongqing and the Suining.

Conclusion

We concluded that in comparison with no fertilization and chemical fertilization, long-term combined application of manure and chemical fertilizers significantly increased soil total and available nutrient content. In comparison to the CK and NPK treatments, the NPKM treatment significantly increased phosphatase activities in both soils. The soil C:N and C:P ratios was highest in the CK treatment, and no significant difference was observed between the NPK and NPKM treatments for the C:N ratio at both sites. While, the lowest soil C:P and N:P ratios were observed in the NPKM treatment at both sites. The NPKM treatment significantly increased P mobility and uptake, and reduced residual-P at both sites via increasing SOM. SEM pathways provided new insights into P lability mediated by the soil C:N, C:P and N:P ratios. The C:N, C:P and N:P ratios had a direct effect on moderately labile-P and non-labile-P, thus indirectly affecting mobile-P and P uptake. The soil organic matter and soil P content at the Chongqing site were relatively higher than those at the Suining site. Therefore, non-labile-P at the Chongqing site directly affected P uptake, and at the Suining site indirectly affected P uptake, indicating that soil C and P were two most important driving factors of P lability in paddy soil. Our findings suggest that soil nutrient stoichiometry is very important for managing the soil P fertility in paddy soil under long-term fertilization in these regions.

Acknowledgments

We acknowledge the staff at the two long-term experimental sites (Chongqing and Suining) for technical support to collect soil samples. We also appreciate Mr. Normatov Yodgor Mamanazarovich (Department of Soil Science, Gulistan State University, Uzbekistan) for his help during the sample analysis. This research was supported by the National Key Research and Development Program of China (2016YFD0300901, 2016YFD0300902 and 2017YFD0800101) and the National Natural Science Foundation of China (Nos. 41671301, 41371293 and 41561070).

Citation: Qaswar M, Ahmed W, Jing H, Hongzhu F, Xiaojun S, Xianjun J, et al. (2019) Soil carbon (C), nitrogen (N) and phosphorus (P) stoichiometry drives phosphorus lability in paddy soil under long-term fertilization: A fractionation and path analysis study. PLoS ONE 14(6): e0218195. https://doi.org/10.1371/journal.pone.0218195

References

1. Williams A, Börjesson G, Hedlund K. Soil Biology & Biochemistry The effects of 55 years of different inorganic fertiliser regimes on soil properties and microbial community composition. Soil Biol Biochem. Elsevier Ltd; 2013;67: 41–46.

2. Malik MA, Marschner P, Khan KS. Addition of organic and inorganic P sources to soil e Effects on P pools and microorganisms. Soil Biol Biochem. Elsevier Ltd; 2012;49: 106–113.

3. Ayaga G, Todd A, Brookes PC. Enhanced biological cycling of phosphorus increases its availability to crops in low-input sub-Saharan farming systems. Soil Biol Biochem. 2006;38: 81–90.

4. Liu Y, Shi G, Mao L, Cheng G, Jiang S, Ma X, et al. Direct and indirect influences of 8 yr of nitrogen and phosphorus fertilization on Glomeromycota in an alpine meadow ecosystem. New Phytol. 2012; 523–535. pmid:22292929

5. Higgs B, Johnston AE, Salter JL, Dawson CJ. Some Aspects of Achieving Sustainable Phosphorus Use in Agriculture. J Environ Qual. 2000; 80–87.

6. Huang J, Xu C, Ridoutt BG, Wang X, Ren P. Nitrogen and phosphorus losses and eutrophication potential associated with fertilizer application to cropland in China. J Clean Prod. Elsevier Ltd; 2017;159: 171–179.

7. Reddy DD, Rao SA, Singh M. Changes in P fractions and sorption in an Alfisol following crop residues application. J Plant Nutr Soil Sci. 2005;168: 241–247.

8. Zhao Y, Wang P, Li J, Chen Y, Ying X, Liu S. The effects of two organic manures on soil properties and crop yields on a temperate calcareous soil under a wheat–maize cropping system. Eur J Agron. 2009;31: 36–42.

9. Delgado A, Madrid A, Kassem S, Andreu L. Phosphorus fertilizer recovery from calcareous soils amended with humic and fulvic acids. Plant Soil. 2002;245: 277–286.

10. Siddique MT, Robinson JS. Phosphorus Sorption and Availability in Soils Amended with Animal Manures and Sewage Sludge. J Environ Qual. 2003;32: 1114–1121. pmid:12809313

11. Mohanty S, Kumar N, Rajan AR. Availability and uptake of phosphorus from organic manures in groundnut (Arachis hypogea L.)–corn (Zea mays L.) sequence using radio tracer technique. Geoderma. 2006;133: 225–230.

12. Agbede TM. Soil & Tillage Research Tillage and fertilizer effects on some soil properties, leaf nutrient concentrations, growth and sweet potato yield on an Alfisol in southwestern Nigeria. Soil Tillage Res. Elsevier B.V.; 2010;110: 25–32.

13. Sharpley AN, Herron S, Daniel T. Overcoming the challenges of phosphorus-based management in poultry farming. J soil water Conserv. 2007;62: 375–389.

14. Zhang H, Xu M, Shi X, Li Z, Huang Q, Wang X. Rice yield, potassium uptake and apparent balance under long-term fertilization in rice-based cropping systems in southern China. Nutr Cycl Agroecosystems. 2010;88: 341–349.

15. Fan HZ, Chen QR, Qin YS, Chen K, Tu SH, Xu MG, et al. Soil carbon sequestration under long-term rice-based cropping systems of purple soil in Southwest China. J Integr Agric. Chinese Academy of Agricultural Sciences; 2015;14: 2417–2425.

16. Xiao Y, De-Jian W, Gang Z, Lu-Ji B, Xiao-Lan P. Soil phosphorous accumulation in long-term P fertilization paddy field and its environmental effects. Chinese J Eco-Agriculture. 2013;4: 393–400 (In Chinese).

17. An N, Wei W, Qiao L, Zhang F, Christie P, Jiang R, et al. Agronomic and environmental causes of yield and nitrogen use efficiency gaps in Chinese rice farming systems. Eur J Agron. Elsevier; 2018;93: 40–49.

18. Xu ZX. Effect of long-term located fertilization on the yields of rice and wheat and soil nutrients. Acta Agric Zhejiangensis. 2009;5: 485–489 (in Chinese).

19. Xin X, Qin S, Zhang J, Zhu A, Yang W, Zhang X. Yield, phosphorus use efficiency and balance response to substituting long-term chemical fertilizer use with organic manure in a wheat-maize system. F Crop Res. Elsevier; 2017;208: 27–33.

20. Wu Q, Zhang S, Zhu P, Huang S, Wang B, Zhao L, et al. Characterizing differences in the phosphorus activation coefficient of three typical cropland soils and the influencing factors under long- term fertilization. PLoS One. 2017;12: 1–15. pmid:28467425

21. Ahmed W, Jing H, Kaillou L, Qaswar M, Khan MN, Jin C, et al. Changes in phosphorus fractions associated with soil chemical properties under long-term organic and inorganic fertilization in paddy soils of southern China. Bui E, editor. PLoS One. 2019;14: e0216881. pmid:31075143

22. Marklein AR, Houlton BZ. Nitrogen inputs accelerate phosphorus cycling rates across a wide variety of terrestrial ecosystems. New Phytol. 2012;193: 696–704. pmid:22122515

23. Yuan ZY, Chen HYH. Decoupling of nitrogen and phosphorus in terrestrial plants associated with global changes. Nat Clim Chang. 2015;5: 465–469.

24. Wang YP, Law RM, Pak B. A global model of carbon, nitrogen and phosphorus cycles for the terrestrial biosphere. Biogeosciences. 2010;7: 2261–2282.

25. Anzoni STM, Rofymow JOHNAT, Ackson ROBJ. Stoichiometric controls on carbon, nitrogen, and phosphorus dynamics in decomposing litter. Ecol Monogr. 2010;80: 89–106.

26. Liu Z, Fu B, Zheng X, Liu G. Soil Biology & Biochemistry Plant biomass, soil water content and soil N: P ratio regulating soil microbial functional diversity in a temperate steppe: A regional scale study. Soil Biol Biochem. Elsevier Ltd; 2010;42: 445–450.

27. Cleveland CC, Liptzin D. S Y N T H E S IS AND E M E R G IN G ID EA S C: N: P stoichiometry in soil: is there a “Redfield ratio” for the microbial biomass? Biogeochemistry. 2007;85: 235–252.

28. Tian H, Chen G, Zhang C, Melillo JM, Hall CAS. Pattern and variation of C: N: P ratios in China ‘ s soils: a synthesis of observational data. Biogeochemistry. 2010;98: 139–151.

29. Gue´rin JE, Parent L-E, Si BC. Spatial and seasonal variability of phosphorus risk indexes in cultivated organic Spatial and seasonal variability of phosphorus risk indexes in cultivated organic soils. Can J Soil Sci. 2011;91: 291?302.

30. Janardhanan L, Daroub SH. Phosphorus Sorption in Organic Soils in South Florida. Soil Sci Soc Am J. 2010;74.

31. Hedley MJ, Stewart JWB, Chauhan BS. Changes in Inorganic and Organic Soil Phosphorus Fractions Induced by Cultivation Practices and by Laboratory Incubations. Soil Sci Soc Am J. 1982;46: 970.

32. Cross AF, Schlesinger WH. A literature review and evaluation of the Hedley fractionation: Applications to the biogeochemical cycle of soil phosphorus in natural ecosystems. Geoderma. 1995;64: 197–214.

33. Condron LM, Newman S. Revisiting the fundamentals of phosphorus fractionation of sediments and soils. J Soils Sediments. 2011;11: 830–840.

34. Gama-Rodrigues AC, Sales MVS, Silva PSD, Comerford NB, Cropper WP, Gama-Rodrigues EF. An exploratory analysis of phosphorus transformations in tropical soils using structural equation modeling. Biogeochemistry. 2014;118: 453–469.

35. Tiecher T, Gomes MV, Ambrosini VG, Amorim MB, Bayer C. Assessing linkage between soil phosphorus forms in contrasting tillage systems by path analysis. Soil Tillage Res. Elsevier; 2018;175: 276–280.

36. Baxter S. World Reference Base for Soil Resources. World Soil Resources Report 103. Rome: Food and Agriculture Organization of the United Nations (2006), pp. 132, US$22.00 (paperback). ISBN 92-5-10511-4. Exp Agric. 2007;43: 264. https://doi.org/10.1017/S0014479706394902

37. Lu RK. Analytical Methods of Soil Agricultural Chemistry. Beijing: China Agricultural Science and Technology Press; 2000.

38. Sparks DL, Page AL, Helmke PA, Loeppert RH, Soltanpour PN, Tabatabai MA, et al. Methods of Soil Analysis. Part 3-Chemical Methods. Madison, USA: Soil Science Society of America; 1996.

39. Walkley A, Black IA. An examination of the Degtjareff method for determining soil organic matter and a proposed modification of the chromic acid titration method. Soil Sci. 1934;37: 29–34.

40. Olsen SR, Sommers LE, Page AL. Methods of soil analysis. Part 2. Chem Microbiol Prop Phosphorus ASA Monogr. 1982; 403–430.

41. Tabatabai MA. No Title. In: S. H. Mickelson JMB, editor. Methods of Soil Analysis Part 2: Microbiological and Biochemical Properties. Madison: Soil Science Society of America; 1994. pp. 775–833.

42. Tiessen H, Moir JO. Characterization of Available P by Sequential Extraction. In: Carter M., Gregorich EG, editors. Soil Sampling and Methods of Analysis. 2nd ed. Taylor & Francis Group, LLC.; 1993. pp. 5–229.

43. Steel RGD, Torrie JH, Dickey DA. Principles and procedures of statistics: A biological approach. McGraw-Hill; 1997.

44. Luo Z, Feng W, Luo Y, Baldock J, Wang E. Soil organic carbon dynamics jointly controlled by climate, carbon inputs, soil properties and soil carbon fractions. Glob Chang Biol. 2017;23: 4430–4439. pmid:28544252

45. Han T, Cai A, Liu K, Huang J, Wang B, Li D, et al. The links between potassium availability and soil exchangeable calcium, magnesium, and aluminum are mediated by lime in acidic soil. J Soils Sediments. Journal of Soils and Sediments; 2018;

46. Mi W, Wu L, Brookes PC, Liu Y, Zhang X, Yang X. Changes in soil organic carbon fractions under integrated management systems in a low-productivity paddy soil given different organic amendments and chemical fertilizers. Soil Tillage Res. Elsevier B.V.; 2016;163: 64–70.

47. Chen D, Yuan L, Liu Y, Ji J, Hou H. Long-term application of manures plus chemical fertilizers sustained high rice yield and improved soil chemical and bacterial properties. Eur J Agron. Elsevier; 2017;90: 34–42.

48. Cui X, Zhang Y, Gao J, Peng F, Gao P. Long-term combined application of manure and chemical fertilizer sustained higher nutrient status and rhizospheric bacterial diversity in reddish paddy soil of Central. Sci Rep. Springer US; 2018; 1–11.

49. Ye G, Lin Y, Liu D, Chen Z, Luo J, Bolan N, et al. Long-term application of manure over plant residues mitigates acidi fi cation, builds soil organic carbon and shifts prokaryotic diversity in acidic Ultisols. Appl Soil Ecol. Elsevier; 2019;133: 24–33.

50. Tian J, Lu S, Fan M, Li X, Kuzyakov Y. Integrated management systems and N fertilization: effect on soil organic matter in rice-rapeseed rotation. Plant Soil. 2013;372: 53–63.

51. Wei W, Yan Y, Cao J, Christie P, Zhang F, Fan M. Effects of combined application of organic amendments and fertilizers on crop yield and soil organic matter: An integrated analysis of long-term experiments. Agric Ecosyst Environ. Elsevier B.V.; 2016;225: 86–92.

52. Ladha JK, Reddy CK, Padre AT, van Kessel C. Role of Nitrogen Fertilization in Sustaining Organic Matter in Cultivated Soils. J Environ Qual. 2011;40: 1756. pmid:22031558

53. Kögel-Knabner I, Amelung W, Cao Z, Fiedler S, Frenzel P, Jahn R, et al. Biogeochemistry of paddy soils. Geoderma. Elsevier B.V.; 2010;157: 1–14.

54. Domagała-S´wia˛tkiewicz , Ga˛stoł M. Soil chemical properties under organic and conventional crop management systems in south Poland. Biol Agric Hortic. 2013;

55. Steiner C, Teixeira WG, Lehmann J, Nehls T, de Macêdo JLV, Blum WEH, et al. Long term effects of manure, charcoal and mineral fertilization on crop production and fertility on a highly weathered Central Amazonian upland soil. Plant Soil. 2007;291: 275–290.

56. Sadeghpour A, Ketterings QM, Godwin GS, Czymmek KJ. Shifting from N-based to P-based manure management maintains soil test phosphorus dynamics in a long-term corn and alfalfa rotation. Agron Sustain Dev. Agronomy for Sustainable Development; 2017;37: 1–11.

57. Wang W, Wang C, Sardans J, Min Q, Zeng C, Tong C, et al. Agricultural land use decouples soil nutrient cycles in a subtropical riparian wetland in China. Catena. Elsevier B.V.; 2015;133: 171–178.

58. Wojciech P, Jaros L, Martin L. Catena A comparison of C: N: P stoichiometry in soil and deadwood at an advanced decomposition stage. Catena. 2019;179: 1–5.

59. Zhang W, Liu W, Xu M, Deng J, Han X, Yang G. Geoderma Response of forest growth to C: N: P stoichiometry in plants and soils during Robinia pseudoacacia a ff orestation on the Loess Plateau, China. Geoderma. Elsevier; 2019;337: 280–289.

60. Sterner RW, Elser JJ. Ecological stoichiometry: the biology of elements from molecules to the biosphere. Princeton University Press; 2002.

61. Acosta-Martinez V, Cruz L, Sotomayor-Rmirz D, Perez-Alegria L. Enzyme activities as affected by soil properties and land use in a tropical watershed. Appl soil Ecol. 2007;35: 35–45.

62. Böhme L, Langer U, Böhme F. Microbial biomass, enzyme activities and microbial community structure in two European long-term field experiments. Agric Ecosyst Environ. 2005;109: 141–152.

63. Fließbach A, Oberholzer H-R, Gunst L, Mäder P. Soil organic matter and biological soil quality indicators after 21 years of organic and conventional farming. Agric Ecosyst Environ. 2007;118: 273–284.

64. Margalef O, Sardans J, Fernández-Martínez M, Molowny-Horas R, Janssens IA, Ciais P, et al. Global patterns of phosphatase activity in natural soils. Sci Rep. 2017;7: 1337. pmid:28465504

65. Nayak DR, Babu YJ, Adhya TK. Long-term application of compost influences microbial biomass and enzyme activities in a tropical Aeric Endoaquept planted to rice under flooded condition. Soil Biol Biochem. 2007;39: 1897–1906.

66. Ohm M, Paulsen HM, Moos JH, Eichler-löbermann B. Long-term negative phosphorus budgets in organic crop rotations deplete plant-available phosphorus from soil. Agron Sustain Dev. Agronomy for Sustainable Development; 2017;37: 17.

67. Hinsinger P. Bioavailability of soil inorganic P in the rhizosphere as affected by root-induced chemical changes: a review. Plant Soil. 2001;237: 173–195.

68. Cao N, Chen X, Cui Z. Change in soil available phosphorus in relation to the phosphorus budget in China. Nutr Cycl Agroecosystems. 2012;94: 161–170.

69. Yan X, Wei Z, Hong Q, Lu Z, Wu J. Phosphorus fractions and sorption characteristics in a subtropical paddy soil as influenced by fertilizer sources. Geoderma. Elsevier B.V.; 2017;295: 80–85.

70. Li G, Li H, Leffelaar PA, Shen J, Zhang F. Characterization of Phosphorus in Animal Manures Collected from Three (Dairy, Swine, and Broiler) Farms in China. Hui D, editor. PLoS One. 2014;9: e102698. pmid:25051245

71. Negassa W, Leinweber P. How does the Hedley sequential phosphorus fractionation reflect impacts of land use and management on soil phosphorus: A review. J Plant Nutr Soil Sci. 2009; 305–325.

72. Dobermann A, George T, Thevs N. Phosphorus Fertilizer Effects on Soil Phosphorus Pools in Acid Upland Soils. Soil Sci Soc Am J. 2002;66: 652–660.

73. Guppy CN, Menzies NW, Moody PW, Blamey FPC. Competitive sorption reactions between phosphorus and organic matter in soil: a review. Soil Res. 2005;43: 189.

74. Adhami E, Owliaie H., Molavi R, Rezaei Rashti M, Esfandbod M. Effects of soil properties on phosphorus fractions in subtropical soils of Iran. J soil Sci plant Nutr. 2013;13: 11–21.

75. Simonsson M, Östlund A, Renfjäll L, Sigtryggsson C, Börjesson G, Kätterer T. Pools and solubility of soil phosphorus as affected by liming in long-term agricultural field experiments. Geoderma. Elsevier; 2018;315: 208–219.

76. Goll DS, Brovkin V, Parida BR, Reick CH, Kattge J, Reich PB, et al. Nutrient limitation reduces land carbon uptake in simulations with a model of combined carbon, nitrogen and phosphorus cycling. Biogeosciences. 2012;9: 3547–3569.

77. González Jiménez JL, Healy MG, Daly K. Effects of fertiliser on phosphorus pools in soils with contrasting organic matter content: A fractionation and path analysis study. Geoderma. Elsevier; 2019;338: 128–135.

78. Lehmann J, Lan Z, Hyland C, Sato S, Solomon D, Ketterings QM. Long-Term Dynamics of Phosphorus Forms and Retention in Manure-Amended Soils. Environ Sci Technol. 2005;39: 6672–6680. pmid:16190226

79. Zhang ZS, Cao CG, Cai ML, Li CF. Crop yield, P uptake and soil organic phosphorus fractions in response to short-term tillage and fertilization under a rape- rice rotation in central China. J Soil Sci Plant Nutr. 2013;13: 871–882.

80. Chen CR, Hou EQ, Condron LM, Bacon G, Esfandbod M, Olley J, et al. Soil phosphorus fractionation and nutrient dynamics along the Cooloola coastal dune chronosequence, southern Queensland, Australia. Geoderma. Elsevier B.V.; 2015;257–258: 4–13.

81. De Groot CC, Marcelis LFM, van den Boogaard R, Kaiser WM, Lambers H. Interaction of nitrogen and phosphorus nutrition in determining growth. Plant Soil. 2003;248: 257–268.

82. Joshi R V, Patel BJ, Patel KM. Effect of Nitrogen Levels and Time of Application on Growth, Yield, Quality, Nitrogen, Phosphorus Content and Uptake for Seed Production of Oat (Avena Sativa L.). Forage Res. 2015;41: 104–108.