Impacts of long-term inorganic and organic fertilization on phosphorus adsorption and desorption characteristics in red paddies in southern China

Abstract

Red soils, having a low or acidic soil pH, are widely distributed in the southern part of China. Because of the robust adsorption of soil P by aluminum-iron (Al-Fe) oxides, the bioavailability of P to crop plants is inadequate in red soils [5]. [...]approximately 80% of the applied P is fixed or adsorbed in the soil [7, 8]. [...]lessening P adsorption and enhancing the availability of P in red paddy soils is an imperative topic to be focused on in current environmental research. [...]the impacts of long-term organic and inorganic fertilizer amendment and changing soil properties on soil P adsorption and desorption in red paddy soils have not been well documented. [...]it is imperative to explore the effects of fertilization on the phytoavailability of P in these paddy soils.

About the Authors:

Waqas Ahmed

Roles Conceptualization, Data curation, Formal analysis, Methodology, Software, Validation, Visualization, Writing – original draft, Writing – review & editing

Affiliations National Engineering Laboratory for Improving Quality of Arable Land, Institute of Agricultural Resources and Regional Planning, Chinese Academy of Agricultural Sciences, Beijing, China, Guangdong Provincial Key Laboratory for Radionuclides Pollution Control and Resources, School of Environmental Science and Engineering, Guangzhou University, Guangzhou, China, School of Civil Engineering, Guangzhou University, Guangzhou, China

Huang Jing

Roles Project administration, Validation

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

Liu Kailou

Roles Project administration, Resources

Affiliation: National Engineering and Technology Research Center for Red Soil Improvement, Jiangxi Institute of Red Soil, Nanchang, China

Sehrish Ali

Roles Conceptualization, Software, Writing – original draft

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

Han Tianfu

Roles Data curation, Project administration

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

Sun Geng

Roles Investigation, Project administration

Affiliation: Hunan Institute of Soil and Fertilizer, Hunan Academy of Agricultural Sciences, Changsha, China

Chen Jin

Roles Project administration

Affiliation: Hunan Institute of Soil and Fertilizer, Hunan Academy of Agricultural Sciences, Changsha, China

Muhammad Qaswar

Roles Software, Visualization

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

Du Jiangxue

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

Sajid Mahmood

Roles Software, Writing – original draft

Affiliations Guangdong Provincial Key Laboratory for Radionuclides Pollution Control and Resources, School of Environmental Science and Engineering, Guangzhou University, Guangzhou, China, School of Civil Engineering, Guangzhou University, Guangzhou, China

Ali Akbar Maitlo

Roles Data curation, Formal analysis, Writing – original draft

Affiliations National Engineering Laboratory for Improving Quality of Arable Land, Institute of Agricultural Resources and Regional Planning, Chinese Academy of Agricultural Sciences, Beijing, China, Soil Fertility Research Institute, Agriculture Research Center, Tando Jam, Pakistan, Department of Agriculture, Government of Sindh, Karachi, Pakistan

Zulqarnain Haider Khan

Roles Software

Affiliations Agro-Environmental Protection Institute, Ministry of Agriculture of China, Tianjin, China, Chinese Academy of Agricultural Sciences, Beijing, China

Huimin Zhang

Roles Funding acquisition, Project administration, Supervision, Validation, Visualization, Writing – original draft

* E-mail: [email protected] (HZ); [email protected] (DYC)

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

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

Di-Yun Chen

Roles Conceptualization, Supervision, Validation, Writing – original draft

* E-mail: [email protected] (HZ); [email protected] (DYC)

Affiliations Guangdong Provincial Key Laboratory for Radionuclides Pollution Control and Resources, School of Environmental Science and Engineering, Guangzhou University, Guangzhou, China, School of Civil Engineering, Guangzhou University, Guangzhou, China

Introduction

Phosphorus (P) is an important macronutrient required for plant growth and therefore is a vital constituent of fertilizers applied to crops grown worldwide [1]. As the reserves of P are limited, attaining an adequate P supply is a large-scale environmental task in the 21st century [2]. P supply will limit increases in agricultural production worldwide in the coming decades. Excessive usage of fertilizers to overcome P deficiency has led to many problems, such as P accumulation in the soil system and surface water runoff [3, 4]. Red soils, having a low or acidic soil pH, are widely distributed in the southern part of China. Because of the robust adsorption of soil P by aluminum-iron (Al-Fe) oxides, the bioavailability of P to crop plants is inadequate in red soils [5]. In the past few decades, large amounts of organic and inorganic P have been applied to these soils to ensure maximum crop yields. This continuous long-term application of P fertilizer results in the gradual depletion of rock phosphate resources [6]. Moreover, approximately 80% of the applied P is fixed or adsorbed in the soil [7, 8]. The successive movement of P applied as fertilizer results in water eutrophication through runoff [9]. Hence, lessening P adsorption and enhancing the availability of P in red paddy soils is an imperative topic to be focused on in current environmental research.

Adsorption and desorption mechanisms are vital aspects of P behavior with respect to the solid phases in soils. Adsorption processes can limit P phyto-availability, where desorption can allow P to migrate or vanish from soils [9, 10]. P adsorption includes a continuous mechanism of adsorption and precipitation. Adsorption is a reversible but rapid process, whereas the process of precipitation is irreversible and slow, and the difference between these processes is unclear [11, 12]. P desorption, the reverse process of adsorption, has an influential role on soil P availability [13, 14]. It has been stated that the soil P adsorption and desorption process is primarily affected by the soil pH, SOC and Al-Fe oxides [15–17]. Al and Fe oxides are the primary sources of fluctuating charge in red soils. Many studies have reported that SOC results in competitive adsorption by occupying some P adsorption sites [18–20]. In contrast, [21] revealed that SOC can boost the adsorption capability of P in some sandy textured soils.

Previous studies have reported that the adsorption and desorption reactions of soil P directly or indirectly affect P availability and mobility in upland soils [22, 23]. The adsorption and desorption processes of soil P can be analyzed by using several isotherms, i.e., the Langmuir and Freundlich models [24, 25]. The key factors are the adsorption capacity of soil (Qm), the adsorption constant (K), the maximum buffer capacity of soil P (MBC), and many others [26]. Continuous fertilizer application can alter soil properties including soil pH, SOC, total P, available P and soil biological indexes [27], which might influence the ability of soil to adsorb P.

To the best of our knowledge, most studies have focused on the adsorption and desorption behavior of P in upland soils with a short experimental duration. However, the impacts of long-term organic and inorganic fertilizer amendment and changing soil properties on soil P adsorption and desorption in red paddy soils have not been well documented. Therefore, it is imperative to explore the effects of fertilization on the phytoavailability of P in these paddy soils. The objectives of this study were to perform batch experiments to examine the adsorption and desorption behavior of P in red paddy soil under long-term organic and inorganic fertilization and the mechanism by which changes in soil properties affect P availability in typical red paddy soils. We hypothesized that the combined application of NPK plus manure would enhance the maximum adsorption capacity of P, P storage capacity, maximum buffering capacity and P binding energy. The study will contribute to an improved understanding of the processes P cycling driven by changes of sorption and adsorption resulted from fertilization practices with intensively managed rice-rice cropping system.

Materials and methods

Study sites and sample collection

The three sites under long-term fertilization selected for this study are located in Nanchang (NC) which belongs to Soil, Fertilizer & Resources and Environment Institute, Jiangxi Academy of Agricultural Sciences, Jinxian (JX) and Ningxiang (NX) that belongs to Red Soil Institute of Jiangxi province, China. Research at these experimental sites was carried out with the consent of the Institute of Agricultural Resources and Agricultural regionalization, Chinese Academy of Agricultural Sciences, Beijing, China, under the cooperative research agreement. The geographic locations, climatic details and physicochemical properties of these three sites are given in Table 1.

[Figure omitted. See PDF.]

Table 1. Locations, climate conditions and initial surface soil properties of the three long-term experimental sites.

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

The agricultural soil of these regions is classified as red soil based on the soil classification system of China [28]. Three treatments were applied in each long-term experiment involving different inorganic and organic fertilization approaches under a rice-based cropping system: (1) CK (unfertilized soil); (2) NPK (synthetic nitrogen, phosphorus and potassium); and (3) NPKM (synthetic NPK plus organic manure). The soil had a heavy clayey texture, containing 38% clay. All the treatments were set in a randomized complete block design (RCBD). Triplicate treatments were separated from each other by cemented ridges to prevent water percolation. Inorganic amendments were applied in the form of urea as nitrogen (N), calcium phosphate as phosphorus (P), and potassium chloride as potassium (K) at each site. The respective annual fertilizer inputs are presented in Table 2.

[Figure omitted. See PDF.]

Table 2. Fertilizer input rates (kg ha−1) at the three long-term experimental sites.

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

For the NPKM treatment, 50% inorganic NPK was applied, and the remaining 50% was added with pig manure based on the N content, while the P and K rates were modified with inorganic fertilizers. Organic manure was applied at rates of 16700 kg ha-1, 15000 kg ha-1, and 4050 kg ha-1 at Nanchang, Jinxian and Ningxiang, respectively. The manure and P fertilizers were applied as a basal dose before rice seedling transplantation, for N and K, 50% was applied as a basal dose, 25% was top-dressed at the tillering stage, and another 25% was added at the panicle stage.

Soil samples were collected from the topsoil (0–20 cm) from all three sites immediately after the early rice harvest using a steel auger. The samples were air-dried, crumbled gently, and passed through a 2 mm sieve before analysis. The representative samples were subjected to analysis for organic C [29], total P [30] and available P [31]. The soil pH was assessed in a 1:2.5 soil and water suspension [32]. The soil properties are given in Table 3.

[Figure omitted. See PDF.]

Table 3. The effect of long-term fertilization on soil pH, SOC, total P and available P concentrations and PAC.

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

Zeta potential measurement

Zeta potential is the potential at the shear plane of the electric double layer on colloidal particles. The value and sign (positive or negative) of the zeta potential depend on the surface charges of the soil particles. The isoelectric point (IEP) indicates the pH value at which the zeta potential is 0 mV [33]. The zeta potentials of soil colloids were measured by electrophoresis using a Zetasizer Nano ZS particle size analyzer (ZEN3600, Malvern, Worcestershire, UK). Briefly, 0.05 g of soil sample (0.05 mm) was weighed into a 250-mL Erlenmeyer flask in duplicate followed by the addition of 200 mL of KCl solution (0.01M). Each suspension was dispersed under ultrasound and then dispensed into four 50-mL plastic bottles (~ 30 mL each). The pH was adjusted to 3, 4, 5, 6, or 7 with HCl or KOH. The zeta potential was measured after 2 h of equilibration [34].

P adsorption and desorption experiments

Adsorption isotherms were assessed using the batch equilibrium method [35]. Briefly, 2.5 g of soil was added to a 50 ml centrifuge tube accompanied by the addition of 25 ml of P working solution, i.e., 0.01 M CaCl2. A known concentration of P was maintained: 0 mg L-1, 10 mg L-1, 20 mg L-1, 40 mg L-1, 60 mg L-1, 100 mg L-1 or 150 mg L-1. Two to three drops of phenol were further added to inhibit any microbial growth. The mixture was shaken overnight at room temperature (180 r.p.m at 25°C) and then centrifuged (4000g) for approximately 10 minutes. Afterwards, 5 mL of the supernatant was filtered and transferred into a 25 mL tube, and molybdenum blue spectroscopy was used to analyze the P concentration. The soil samples used for adsorption were then washed two times with saturated NaCl solution, centrifuged and subjected to the desorption experiment. After supernatant removal, 25 mL of CaCl2 and 3 drops of phenol were added to each sample, followed by overnight shaking. After centrifugation for 10 minutes, 5 mL of the supernatant was taken for the measurement of P concentration by the same molybdenum blue method. This solution P concentration was defined as the desorbed P.

P adsorption and desorption models

The Langmuir and Freundlich models are widely used for the quantitative description of the P adsorption characteristics of soil particle surfaces [36]. The adsorption isotherms are expressed as:(1)and(2)where C is the concentration of P in solution at the equilibrium stage (mg L-1); Q is the adsorption capacity (mg kg-1); Qm is the maximum P adsorption (mg kg-1); K is the equation constant, defined as the factor indicating the intensity of P adsorption, with a maximum quantity consistent with higher soil P adsorption [37]; and n is the constant of heterogeneity related to the intensity of adsorption.

The maximum buffer capacity of P in soil (MBC, mg kg-1) was calculated from the series of indexes derived from the Langmuir and Freundlich equations. MBC is an index that merges the P adsorption capacity (Qm) with its intensity (K). The resilience of soil to variations in solution P concentration is explained by MBC [38]. This parameter is expressed as:(3)

Statistical analysis

All statistical analyses were performed using SPSS 20.0 (IBM SPSS, Somers, NY, USA). Data fitting and mapping were conducted in Origin Pro 19.0 (Origin Lab Corp, MA, USA). Data were subjected to one-way ANOVA for comparison of treatments and two-way ANOVA for comparison of sites, followed by the least significant difference (LSD) test at the P ≤ 0.05 level to calculate significant differences among the mean values. Pearson’s correlation coefficients (r) were computed using R 3.6.1 (R Foundation for Statistical Computing, Vienna, Austria). Correlation coefficients were evaluated to determine the relationships between different soil properties, Qm, and Dm.

Results and discussion

SOC and total and available P contents in soil

Long-term continuous fertilization significantly (P ≤ 0.05) affected the soil chemical properties at all three sites (Table 3). Soil pH showed a nonsignificant trend among the treatments, and it ranged from 5.23 in Nanchang to 6.91 in Ningxiang. Compared with the CK treatments, NPK and NPKM significantly increased the SOC, total P and available P concentrations in the soil. NPKM showed maximum SOC contents of 28.09 mg kg-1, 21.76 mg kg-1 and 21.64 mg kg-1 at Nanchang, Jinxian and Ningxiang, respectively. These results could be due to the manure addition, which supplies enough nutrients to the soil, primarily by directly increasing carbon inputs into the soil, which assists in carbon sequestration. Similar results were reported in previous long-term studies [39, 40]. Long-term NPK plus manure addition significantly increased the soil total P and available P concentrations in the soil at all three sites (Table 3). The maximum values were observed at Ningxiang, with values of 2.21 g kg-1 and 90.71 mg kg-1, respectively. [41] also reported that long-term fertilization increased soil nutrient concentrations. Similar studies reported that increases in SOC, total P and available P concentrations could be associated with long-term soil mineral and organic inputs that build SOM in the topsoil [42, 43].

The proportion of available P to total P is defined as the P activation coefficient, i.e., PAC, which characterizes the degree of difficulty with which transformation occurs among total P and available P. PAC is a crucial soil fertility indicator, and a higher PAC indicates a maximum P that could be available for plant growth [4, 43]. Similar to SOC, total P and available P, PAC increased significantly with NPKM treatment compared with NPK and CK (Table 3). The PAC values ranged from 1.02 to 5.03 among the treatments at all three sites.

Characteristics of P adsorption

P adsorption isotherms.

The sorption isotherms presented in Fig 1 show that the soil samples adsorbed extra P as the exogenous P concentration increased; however, there was a gradual decline in the increasing rate of P adsorption capacity (Fig 1). The variations in the P adsorption capacity between different treatments became increasingly pronounced with increasing added P concentration. Among the treatments, the P adsorption capacity at an added P concentration of 150 mg L-1 was maximum for NPKM, with values of 1786.54 mg kg-1, 1684.39 mg kg-1, and 1663.05 mg kg-1 at Jinxian, Ningxiang and Nanchang, respectively. Overall, the P adsorption capacity was higher with NPKM and lower in the CK treatment across all three sites. These results were in line with previous studies [44–46]. This observation was possibly due to the adsorption capacity of the solid phase and the residence time (contact time between the soil and solution) [47]. The P adsorption mechanism comprises two different phases, i.e., chemical and physical adsorption. The process of chemical adsorption dominated at comparatively low exogenous P intensities and was completed rapidly. Ligand and ion exchange are the probable dominant mechanisms that contribute to higher adsorption rates [48]. This phenomenon is commonly termed the relatively quick adsorption phase. At higher P concentrations, the process of chemical adsorption slows rapidly because the existing adsorption sites become saturated so quickly, and the P present in the liquid phase is physicochemically adsorbed to the surface of the soil at a relatively slow rate. This process is termed the slower adsorption phase [44, 48].

[Figure omitted. See PDF.]

Fig 1. The adsorption isotherms of P in red paddy soils under three different long-term experimental sites.

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

Al and Fe oxides are vital for P adsorption, ultimately affecting its availability in soil [49]. SOC alters soil P adsorption by blocking the adsorption sites on Fe/Al oxides or by the formation of Fe/Al-SOC-P multicomplexes with relatively variable structures [50, 51]. Significant differences were found in the P adsorbed amount in soils under different treatments that were subjected to solution with the same exogenous content. The NPKM treatment had the maximum P retention capacity, followed by NPK and CK. This result was probably due to the higher SOC content in NPKM than in the other treatments. These results were also supported by the study of [44].

Adsorption equations.

Various models have been developed and used for the quantitative description of P adsorption isotherms. The Langmuir and Freundlich models are widely used for P adsorption [10]. The data presented in Table 4 show that the adsorption equations fit the P adsorption isotherms well for different treatments at different sites. The correlation coefficients (R2) showed that all the correlations were quite significant. Therefore, either of these equations could be selected to illustrate the characteristics of P adsorption in soils with different treatments. These results were consistent with previous studies [10, 52, 53]. Stated that the P adsorption capacity and its availability are commonly described by the MBC calculated from the Langmuir isotherm, the maximum adsorption capacity (Qm), and the binding energy constant (K).

[Figure omitted. See PDF.]

Table 4. Isotherm parameters of the Langmuir and Freundlich models for soil P adsorption.

https://doi.org/10.1371/journal.pone.0246428.t004

Parameters of P adsorption.

The maximum P adsorption capacity (Qm) represents the quantity of available adsorption sites per unit soil weight and is commonly used for the determination of soil P adsorption capacity [27, 54]. The Qm values for different treatments across the sites ranged from 549.15 mg kg-1 to 2219.36 mg kg-1 (Table 4). Long-term inorganic fertilization and manure addition clearly influenced Qm. Qm increased with NPKM treatment compared to the other treatments, and the relationships were statistically significant (P ≤ 0.05). This difference could be due to increased SOM concentrations because of manure application in the NPKM treatments across the sites, demonstrating that the capacity of soil P storage was enhanced with increased SOC content. Similar results were reported by [55]. In contrast [27], revealed that the SOC content does not directly influence Qm. These conflicting findings might have been determined by other soil factors, such as soil pH, soil type, Fe and Al forms, etc., that could have an effect on soil P adsorption capacity. Further studies are recommended to improve the understanding of such influencing factors.

In addition, it has previously been reported that increased SOC content enhances the P adsorption capacity of soils at pH values greater than 6 but decreases the P adsorption capacity at pH values lower than 6 [56]. The P binding energy (K) is important in defining the affinity of soil for P. A maximum K value indicates much stronger adsorption of P [44]. An uneven trend in K values was observed among the treatments (Table 4). This variation could be explained by the SOC content available to react with phosphate anions, as mentioned earlier. MBC is a unified index that merges K and Qm [44]. A high MBC suggests that there will be more P adsorption. The NPKM treatment showed a higher MBC than did the NPK and CK treatments. Similar results were reported in previous studies [57, 58]. This finding led to the assumption that soil P adsorption was mainly controlled by the SOM content in the soil.

Characteristics of P desorption

P desorption isotherms.

Soil desorption is thought to be an inverse process of sorption. It is considered more important than adsorption due to the significance of immobilized P in soil, which can be reused, and due to environment-related issues that could be caused by P released from soil [4, 44]. The concentration of P desorbed from soil was lower than the amount of P adsorbed in CK, NPK and NPKM (Fig 2). This result indicates that the P adsorbed to the soil could easily be desorbed and be released back into the soil system but that the adsorption was not fully reversible. With increasing exogenous P concentration, the desorbed P amount increased gradually in each treatment, and the trend was more obvious in the case of CK. The rate of P desorption was lower initially when the exogenous P content in solution was lower (< 40 mg L-1) but was approximately 20% higher at added P concentrations > 60 mg L-1. This result shows that there were abundant adsorption sites on soil colloids when the exogenous P content was lower and there was a higher chemical adsorption binding capacity for soil P, which resulted in a higher adsorption degree and lower desorption. The adsorption sites on soil colloids gradually became covered and saturated as the concentration of exogenous P increased, which resulted in a decreased adsorption binding capacity of the soil. This process is generally described as a physical adsorption level at which physically bound P can be easily desorbed. Consequently, a high P application rate could increase the desorbed P amount.

[Figure omitted. See PDF.]

Fig 2. Desorption isotherms of P in red paddy soils at three different long-term experimental sites.

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

P desorption parameters.

The desorption equations fit the isotherms of the CK, NPK and NPKM treatments well, showing significant (P ≤ 0.05) correlation coefficients (Table 5). Either the Langmuir or Freundlich equation could be used to describe the desorption characteristics, as was observed in a previous study by Wang and Wang [59]. The amount of P desorbed when the adsorbent was saturated with exogenous P, termed the maximum desorption capacity (Dm), reflected the maximum amount of P provided to the soil. The Dm values ranged from 78.62 mg kg-1 to 157.58 mg kg-1 across the three sites under CK, NPK and NPKM. The Dm values were significantly higher in the NPKM treatment than in the NPK and CK treatments. This difference could be attributed to the higher SOM content in NPKM due to manure addition, showing competition for adsorption sites between SOM and P. This result implies that more SOC entered the soil solution and enhanced the desorbed P amount. However, from an environmental pollution risk assessment perspective, the addition of manure may increase the risk of P loss, primarily by surface water eutrophication. Similar results were previously reported [4].

[Figure omitted. See PDF.]

Table 5. P desorption isotherm parameters of the Langmuir and Freundlich models.

https://doi.org/10.1371/journal.pone.0246428.t005

The results and findings in this study suggest practical implications in soil fertility management, primarily in the use of manure addition to control soil SOM content. These findings will assist in regulating P adsorption and desorption in soil to regulate P nutrition in crop plants. By increasing the SOC content in soil by fertilizer application, the P fixation capacity can be increased. However, when there is a surplus SOC content, sufficient fertilizer must be added to sustain the maximum P supply under long-term experiments.

Relationships between different soil properties and parameters of P adsorption/desorption

The correlation matrix showed correlations between soil pH, SOC, total P, available P, Qm, MBC and Dm (Fig 3). SOC content, total P and available P showed a strong positive correlation with Qm, MBC and Dm.

[Figure omitted. See PDF.]

Fig 3. Correlation analysis among soil properties, Qm, and Dm.

The upper diagonal represents the correlation coefficient (the overall value is shown in black; the treatments are differentiated by color); the lower diagonal plots show the distribution of data along the axis; the box plots show the overall difference between the treatments.

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

Organic carbon in soil regulates nutrient supply, and it also interacts with and boosts the supply of other essential nutrients in soil. P adsorption and desorption processes in soil are greatly influenced by SOM content. A similar study was reported by [4]. Fig 4 shows a conceptual framework showing the linkage among soil pH, SOC, total P (TP), available P (AP), Qm, Dm and MBC. This linkage is based on the relationships between different variables. The variations in the changes in Qm and Dm were mainly due to SOC, TP and AP. SOC played a dominant role among these factors. Subsequently, adsorption and desorption were affected by the SOC content in the soil. The absolute value of zeta potential increased with increasing soil pH and ranked in the following order for the three soils: NX< NC <JX (Fig 5). This result indicates that the P adsorption capacity of NX soil was the lowest, while that of JX soil was the highest among the three sites. The high P adsorption capacity of JX soil was confirmed by its high IEP, which indicates that the soil surface was positively charged and advantageous to P adsorption. This is probably due to the naturally acidic parent material which is characteristic of red soils. The weathering of parent material could result in abundant Fe-Al oxides, leading to substantial leaching of soluble soil minerals and basic cations [60].

[Figure omitted. See PDF.]

Fig 4. A proposed conceptual framework showing the linkage among soil pH, SOC, total P (TP), available P (AP), Qm, Dm and the MBC.

The linkage is based on the relationships between different variables. The color of each line refers to the degree of relationship between the variables. The variations in the changes in Qm and Dm were mainly due to SOC, TP and AP. SOC played a dominant role among these factors. Consequently, adsorption and desorption were affected by the SOC content in the soil.

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

[Figure omitted. See PDF.]

Fig 5. Curves of zeta potential vs. pH for red soils at three different experimental sites.

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

Conclusions

The relationships between the soil properties of red paddy soils from southern China under long-term fertilization experiments and their adsorption-desorption characteristics were assessed by different batch experiments. The Langmuir and Freundlich equations fit the isotherms well for every treatment at each site. The combined application of NPK plus manure enhanced the maximum adsorption capacity of P, P storage capacity, maximum buffering capacity and P binding energy. Manure addition significantly altered the MBC and Dm values and increased the P desorption process, suggesting an effective supply of P to crops from the soil. In short, it is suggested that manure addition is important to improve P utilization in red paddies within the recommended range to avoid the environmental risk of pollution.

Acknowledgments

We are thankful to all managers of long-term experiments at Nanchang, Jinxian and Ningxiang, China for their valuable work associated with the long-term Monitoring Network of Soil Fertility and Fertilizer Effects in China.

Citation: Ahmed W, Jing H, Kailou L, Ali S, Tianfu H, Geng S, et al. (2021) Impacts of long-term inorganic and organic fertilization on phosphorus adsorption and desorption characteristics in red paddies in southern China. PLoS ONE 16(1): e0246428. https://doi.org/10.1371/journal.pone.0246428

References

1. Liu S, Meng J, Jiang L, Yang X, Lan Y, Cheng X, et al. Rice husk biochar impacts soil phosphorous availability, phosphatase activities and bacterial community characteristics in three different soil types. Appl Soil Ecol. 2017;116: 12–22.

2. Fink JR, Inda A V., Bavaresco J, Barrón V, Torrent J, Bayer C. Adsorption and desorption of phosphorus in subtropical soils as affected by management system and mineralogy. Soil Tillage Res. 2016;155: 62–68.

3. van der Salm C, van Middelkoop JC, Ehlert PAI. Changes in soil phosphorus pools of grasslands following 17 yrs of balanced application of manure and fertilizer. Soil Use Manag. 2017;33: 2–12.

4. Yang X, Chen X, Yang X. Effect of organic matter on phosphorus adsorption and desorption in a black soil from Northeast China. Soil Tillage Res. 2019;187: 85–91.

5. Messiga AJ, Ziadi N, Angers DA, Morel C, Parent LE. Tillage practices of a clay loam soil affect soil aggregation and associated C and P concentrations. Geoderma. 2011;164: 225–231.

6. Gilbert N. The disappearing nutrient: phosphate-based fertilizers have helped spur agricultural gains in the past century, but the world may soon run out of them. Natasha Gilbert. Nature. 2009;461: 716–719. Available: http://go.galegroup.com/ps/i.do?p=AONE&sw=w&issn=00280836&v=2.1&it=r&id=GALE%7CA210224903&sid=googleScholar&linkaccess=fulltext pmid:19812648

7. Zhu J, Li M, Whelan M. Phosphorus activators contribute to legacy phosphorus availability in agricultural soils: A review. Sci Total Environ. 2018;612: 522–537. pmid:28865270

8. Withers PJA, Edwards AC, Foy RH. Phosphorus cycling in UK agriculture and implications for phosphorus loss from soil. Soil Use Manag. 2006;17: 139–149.

9. Gérard-Marchant P, Hively WD, Steenhuis TS. Distributed hydrological modelling of total dissolved phosphorus transport in an agricultural landscape, part I: Distributed runoff generation. Hydrol Earth Syst Sci. 2006;10: 245–261.

10. Lair GJ, Zehetner F, Khan ZH, Gerzabek MH. Phosphorus sorption-desorption in alluvial soils of a young weathering sequence at the Danube River. Geoderma. 2009;149: 39–44.

11. Arias M, Da Silva-Carballal J, García-Río L, Mejuto J, Núñez A. Retention of phosphorus by iron and aluminum-oxides-coated quartz particles. J Colloid Interface Sci. 2006;295: 65–70. pmid:16125184

12. Gou X, Cai Y, Wang C, Li B, Zhang Y, Tang X, et al. Effects of different long-term cropping systems on phosphorus adsorption and desorption characteristics in red soils. J Soils Sediments. 2020;20: 1371–1382.

13. Heidari S, Reyhanitabar A, Oustan S. Kinetics of phosphorus desorption from calcareous soils using DGT technique. Geoderma. 2017;305: 275–280.

14. McDonald GJ, Norton SA, Fernandez IJ, Hoppe KM, Dennis J, Amirbahman A. Chemical controls on dissolved phosphorus mobilization in a calcareous agricultural stream during base flow. Sci Total Environ. 2019;660: 876–885. pmid:30743973

15. Quesada CA, Lloyd J, Anderson LO, Fyllas NM, Schwarz M, Czimczik CI. Soils of Amazonia with particular reference to the RAINFOR sites. Biogeosciences. 2011;8: 1415–1440.

16. Fang H, Cui Z, He G, Huang L, Chen M. Phosphorus adsorption onto clay minerals and iron oxide with consideration of heterogeneous particle morphology. Sci Total Environ. 2017;605–606: 357–367. pmid:28668747

17. Bornø ML, Müller-Stöver DS, Liu F. Contrasting effects of biochar on phosphorus dynamics and bioavailability in different soil types. Sci Total Environ. 2018;627: 963–974. pmid:29426221

18. Liu Y, Zhu ZQ, He XS, Yang C, Du YQ, Huang YD, et al. Mechanisms of rice straw biochar effects on phosphorus sorption characteristics of acid upland red soils. Chemosphere. 2018;207: 267–277. pmid:29803158

19. Zhang Y, Huang S, Guo D, Zhang S, Song X, Yue K, et al. Phosphorus adsorption and desorption characteristics of different textural fluvo-aquic soils under long-term fertilization. J Soils Sediments. 2019;19: 1306–1318.

20. 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. PLoS One. 2019;14: e0216881. pmid:31075143

21. Debicka M, Kocowicz A, Weber J, Jamroz E. Organic matter effects on phosphorus sorption in sandy soils. Arch Agron Soil Sci. 2016;62: 840–855.

22. Bera R, Seal A, Bhattacharyya P, Mukhopadhyay K, Giri R. Phosphate sorption desorption characteristics of some ferruginous soils of tropical region in Eastern India. Environ Geol. 2006;51: 399–407.

23. Ahmed W, Kailou L, Qaswar M, Jing H, Qinghai H, Yongmei X, et al. Long-term mineral fertilization improved the grain yield and phosphorus use efficiency by changing soil P fractions in ferralic Cambisol. Agronomy. 2019;9: 784.

24. Bhattacharyya P, Nayak AK, Shahid M, Tripathi R, Mohanty S, Kumar A, et al. Effects of 42-year long-term fertilizer management on soil phosphorus availability, fractionation, adsorption–desorption isotherm and plant uptake in flooded tropical rice. Crop J. 2015;3: 387–395.

25. Guedes RS, Melo LCA, Vergütz L, Rodríguez-Vila A, Covelo EF, Fernandes AR. Adsorption and desorption kinetics and phosphorus hysteresis in highly weathered soil by stirred flow chamber experiments. Soil Tillage Res. 2016;162: 46–54.

26. Li RY, Qiu YX, Liu CY, Wang YL, Huang LD. Adsorption-desorption behaviors of phosphorus under different silicon concentrations in paddy soils. Chin J Soil Sci. 2013;44: 1134–1139.

27. Yan X, Wang D, Zhang H, Zhang G, Wei Z. Organic amendments affect phosphorus sorption characteristics in a paddy soil. Agric Ecosyst Environ. 2013;175: 47–53.

28. Zhang WL, Xu AG, Zhang RL, Ji HJ, Wu SX. Review of soil classification and revision of China soil classification system. Sci Agric Sin. 2014;47: 3214–3230.

29. Walkley A. An examination of methods for determining organic carbon and nitrogen in soils. (with one text-figure). J Agric Sci. 1935;25: 598–609.

30. Murphy J, Riley JP. A modified single solution method for the determination of phosphate in natural waters. Anal Chim Acta. 1962;27: 31–36.

31. Olsen SR, Cole C V, Watandbe F, Dean L. Estimation of Available Phosphorus in Soil by Extraction with sodium Bicarbonate. Journal of Chemical Information and Modeling. United States Department Of Agriculture; Washington; 1954.

32. Lu RK. Soil agricultural chemical analysis method. China Agric Sci Technol Press Beijing. 2000; 1–315.

33. Xu RK, Qafoku NP, Van Ranst E, Li JY, Jiang J. Adsorption Properties of Subtropical and Tropical Variable Charge Soils: Implications from Climate Change and Biochar Amendment. Advances in Agronomy. 2016.

34. Jiang J, Xu R, Wang Y, Zhao A. The mechanism of chromate sorption by three variable charge soils. Chemosphere. 2008. pmid:18291439

35. Sui Y, Thompson ML. Phosphorus Sorption, Desorption, and Buffering Capacity in a Biosolids-Amended Mollisol. Soil Sci Soc Am J. 2000;64: 164–169.

36. Lu H, Yang L, Shabbir S, Wu Y. The adsorption process during inorganic phosphorus removal by cultured periphyton. Environ Sci Pollut Res. 2014;21: 8782–8791. pmid:24728572

37. BACHE BW, WILLIAMS EG. a Phosphate Sorption Index for Soils. J Soil Sci. 1971;22: 289–301.

38. Holford ICR. Soil phosphorus: Its measurement, and its uptake by plants. Australian Journal of Soil Research. 1997. pp. 227–239.

39. Qaswar M, Jing H, Ahmed W, Dongchu L, Shujun L, Ali S, et al. Long-term green manure rotations improve soil biochemical properties, yield sustainability and nutrient balances in acidic paddy soil under a rice-based cropping system. Agronomy. 2019;9: 780.

40. Cai A, Zhang W, Xu M, Wang B, Wen S, Shah SAA. Soil fertility and crop yield after manure addition to acidic soils in South China. Nutr Cycl Agroecosystems. 2018;111: 61–72.

41. Lan ZM, Lin XJ, Wang F, Zhang H, Chen CR. Phosphorus availability and rice grain yield in a paddy soil in response to long-term fertilization. Biol Fertil Soils. 2012;48: 579–588.

42. Fageria NK, Dos Santos AB, Moraes MF. Influence of urea and ammonium sulfate on soil acidity indices in lowland rice production. Commun Soil Sci Plant Anal. 2010;41: 1565–1575.

43. Wu Q, Zhang S, Zhu P, Huang S, Wang B, Zhao LP, et al. Characterizing differences in the phosphorus activation coefficient of three typical cropland soils and the influencing factors under longterm fertilization. PLoS One. 2017;12. pmid:28467425

44. Wang L, Liang T. Effects of exogenous rare earth elements on phosphorus adsorption and desorption in different types of soils. Chemosphere. 2014;103: 148–155. pmid:24342358

45. Siddique MT, Robinson JS. Differences in Phosphorus Retention and Release in Soils Amended with Animal Manures and Sewage Sludge. Soil Sci Soc Am J. 2004;68: 1421–1428.

46. Abdala DB, Ghosh AK, da Silva IR, de Novais RF, Alvarez Venegas VH. Phosphorus saturation of a tropical soil and related P leaching caused by poultry litter addition. Agric Ecosyst Environ. 2012;162: 15–23.

47. Bhadha JH, Daroub SH, Lang TA. Effect of kinetic control, soil: Solution ratio, electrolyte cation, and others, on equilibrium phosphorus concentration. Geoderma. 2012;173–174: 209–214.

48. Lai DYF, Lam KC. Phosphorus sorption by sediments in a subtropical constructed wetland receiving stormwater runoff. Ecol Eng. 2009;35: 735–743.

49. Gérard F. Clay minerals, iron/aluminum oxides, and their contribution to phosphate sorption in soils—A myth revisited. Geoderma. 2016;262: 213–226.

50. Wang Y, Zhang Y, He Y. Effect of soil matrix components on phosphate sorption index in red soil. Acta Ped0L0Gica Sin. 2012;49: 552–559.

51. Moshi AO, Wild A, Greenland DJ. Effect of organic matter on the charge and phosphate adsorption characteristics of kikuyu red clay from kenya. Geoderma. 1974;11: 275–285.

52. Lin J, Zhang Z, Zhan Y. Effect of humic acid preloading on phosphate adsorption onto zirconium-modified zeolite. Environ Sci Pollut Res. 2017;24: 12195–12211. pmid:28353102

53. Yan J, Jiang T, Yao Y, Lu S, Wang Q, Wei S. Preliminary investigation of phosphorus adsorption onto two types of iron oxide-organic matter complexes. J Environ Sci (China). 2016;42: 152–162. pmid:27090706

54. Yokoyama D, Mori T, Wagai R, Hiradate S, Kitayama K. Characteristics of phosphorus fractions in the soils derived from sedimentary and serpentinite rocks in lowland tropical rain forests, Borneo. Soil Sci Plant Nutr. 2018;64: 218–221.

55. Zhou Q, Zhu Y. Potential pollution and recommended critical levels of phosphorus in paddy soils of the southern Lake Tai area, China. Geoderma. 2003;115: 45–54.

56. Zhong X, Zhao X, Bao H, Li H, Li G, Lin Q. The evaluation of phosphorus leaching risk of 23 Chinese soils Ⅰ. Leaching criterion. Acta Ecol Sin. 2004;24: 2275–2280.

57. WILLIAMS EG, SAUNDERS WMH. Distribution of Phosphorus in Profiles and Particle‐Size Fractions of Some Scottish Soils. J Soil Sci. 1956;7: 90–109.

58. Wang JB, Chen ZH, Chen LJ, Zhu AN, Wu ZJ. Surface soil phosphorus and phosphatase activities affected by tillage and crop residue input amounts. Plant, Soil Environ. 2011;57: 251–257.

59. Wang S, Wang E, others. Desorption characteristics of phosphorus from different used sandy soil in Western Liao River Basin. Res Environ Sci. 2011;24: 756–762.

60. Tellen VA, Yerima BPK. Effects of land use change on soil physicochemical properties in selected areas in the North West region of Cameroon. Environ Syst Res. 2018.