Impacts of long-term inorganic and organic fertilization on phosphorus adsorption and desorption characteristics in red paddies in southern China
DOI: 10.1371/journal.pone.0246428, Volume: 16, Issue: 1
Article Type: research-article, Article History

Soil phosphorus (P) adsorption and desorption occur in an important endogenous cycle linked with soil fertility problems and relevant to the environmental risk assessment of P. In our study, the effect of long-term inorganic and organic fertilization on P adsorption and desorption characteristics in relation to changes in soil properties was evaluated by selecting three long-term experimental sites in southern China. The selected treatments at each site were CK (unfertilized), NPK (synthetic nitrogen, phosphorus and potassium) and NPKM (synthetic NPK plus manure). The adsorption and desorption characteristics of P were evaluated using Langmuir and Freundlich isotherms. The results showed that long-term application of NPK plus manure significantly increased soil organic carbon (SOC), total P and available P at all three sites compared with the NPK and CK treatments. All three treatments fit these equations well. The maximum adsorption capacity (Qm) of P increased with NPKM treatment, and the binding energy of P (K) and the maximum buffering capacity (MBC) showed increasing trends. NPKM showed the highest Qm (2346.13 mg kg-1) at the Jinxian site, followed by Nanchang (221.16 mg kg-1) and Ningxiang (2219.36 mg kg-1). Compared to CK and NPK, the NPKM treatment showed a higher MBC as 66.64, 46.93 and 44.39 L kg-1 at all three sites. The maximum desorption capacity (Dm) of P in soil was highest with the NPKM treatment (157.58, 166.76, 143.13 mg kg-1), showing a better ability to release P in soil. The correlation matrix showed a significant positive correlation of SOC, total and available P with Qm, Dm and MBC. In conclusion, it is suggested that manure addition is crucial to improve P utilization in red paddy soils within the recommended range to avoid the risk of environmental pollution.

Ahmed, Jing, Kailou, Ali, Tianfu, Geng, Jin, Qaswar, Jiangxue, Mahmood, Akbar Maitlo, Haider Khan, Zhang, Chen, and Hui: Impacts of long-term inorganic and organic fertilization on phosphorus adsorption and desorption characteristics in red paddies in southern China


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 [1517]. 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 [1820]. 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.

Table 1
Locations, climate conditions and initial surface soil properties of the three long-term experimental sites.
Initiation year198419811986
Latitude (N)28.6228.7128.38
Longitude (E)115.99116.36112.72
Mean annual temperature (°C)17.818.516.2
Mean annual precipitation (mm)163215501570
Cropping systemRice-RiceRice-RiceRice-Rice
Soil classification in FAOEutric cambisolEutric cambisolEutric cambisol
Soil textureClay loamClay loamClay loam
Soil pH6.086.926.52
SOM (g kg-1)26.3529.3628.54
TN (g kg-1)1.371.522.05
AN (mg kg-1)82.36139.66142.36
TP (g kg-1)0.520.530.62
AP (mg kg-1)22.369.8113.38
TK (g kg-1)4.2613.2521.36
AK (mg kg-1)37.0583.3434.48
Fed (g kg-1)35.6655.2132.24
Ald (g kg-1)6.8110.545.46
*Abbreviations: MT: monsoon temperate, SM: subtropical monsoon, R-R: rice-rice, SOM: soil organic matter, TN: total nitrogen, AN: available nitrogen, TP: total phosphorus, AP: available phosphorus, TK: total potassium, AK: available potassium, Fed: dithionite-citrate bicarbonate-extractable Fe and Ald: dithionite-citrate bicarbonate-extractable Al.

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.

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

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.

Table 3
The effect of long-term fertilization on soil pH, SOC, total P and available P concentrations and PAC.
SitesTreatmentspHSOC (mg kg-1)Total P (g kg-1)AP (mg kg-1)PAC (%)
NanchangCK5.71 ± 0.2 Ab16.97 ± 0.24 Ca0.36 ± 0.15 Cb7.30 ± 0.24 Ca2.03
NPK5.23 ± 0.12 Bb23.82 ± 0.13 Ba0.54 ± 0.05 Bb21.09 ± 1.35 Ba3.91
NPKM5.74 ± 0.15 Ab28.09 ± 1.33 Aa1.63 ± 0.07 Ac82.05 ± 2.14 Ab5.03
JinxianCK5.62 ± 0.15 ABb13.25 ± 0.52 Cb0.49 ± 0.04 Ca6.51 ± 0.55 Ca1.33
NPK5.50 ± 0.12 Bb18.87 ± 0.99 Bb0.93 ± 0.02 Ba30.14 ± 1.44 Bb3.24
NPKM5.94 ± 0.14 Ab21.76 ± 1.42 Ab1.63 ± 0.14 Ab64.78 ± 1.55 Ac3.97
NingxiangCK6.34 ± 0.11 Ba12.11 ± 0.74 Cc0.45 ± 0.05 Ca4.61 ± 0.39 Cb1.02
NPK6.51 ± 0.15 Aa14.35 ± 0.92 Bc0.96 ± 0.09 Ba14.82 ± 0.18 Bc1.54
NPKM6.71 ± 0.12 Aa21.64 ± 1.24 Ab2.21 ± 0.16 Aa90.71 ± 2.54 Aa4.10
*Abbreviations: PAC: phosphorus activation coefficient of soil.
Treatments: CK: unfertilized control; NPK: inorganic nitrogen, phosphorus and potassium; and NPKM: NPK plus manure.
Data (means ± SD, n = 3) followed by different uppercase letters denote significant differences (P≤0.05) between fertilization treatments at the same site (A, B, and C), and lowercase letters denote significant differences (P ≤ 0.05) between sites for the same fertilization treatment (a, b, and c).

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:

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:


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 [4446]. 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].

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

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).

Table 4
Isotherm parameters of the Langmuir and Freundlich models for soil P adsorption.
SitesTreatmentsLangmuir equationFreundlich equation
Qm (mg kg-1)KL (L mg-1)R2MBC (L kg-1)nKfR2
*Abbreviations: Qm: maximum P adsorption capacity, KL: P adsorption constant, R2: model fitting degree, MBC: maximum buffer capacity of soil P.

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.

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

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].

Table 5
P desorption isotherm parameters of the Langmuir and Freundlich models.
SitesTreatmentsLangmuir equationFreundlich equation
Dm (mg kg-1)KL (L mg-1)R2nKfR2
*Abbreviations: Dm: maximum P desorption capacity, KL: P desorption constant, R2: model fitting degree.

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.

Correlation analysis among soil properties, Qm, and Dm.
Fig 3
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.Correlation analysis among soil properties, Qm, and Dm.

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].

A proposed conceptual framework showing the linkage among soil pH, SOC, total P (TP), available P (AP), Qm, Dm and the MBC.
Fig 4
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.A proposed conceptual framework showing the linkage among soil pH, SOC, total P (TP), available P (AP), Qm, Dm and the MBC.
Curves of zeta potential vs. pH for red soils at three different experimental sites.
Fig 5
Curves of zeta potential vs. pH for red soils at three different experimental sites.


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.


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.


LiuS, MengJ, JiangL, YangX, LanY, ChengX, 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: 1222. doi: 10.1016/j.apsoil.2017.03.020

FinkJR, IndaA V., BavarescoJ, BarrónV, TorrentJ, BayerC. Adsorption and desorption of phosphorus in subtropical soils as affected by management system and mineralogy. Soil Tillage Res. 2016;155: 6268. doi: 10.1016/j.still.2015.07.017

van der SalmC, van MiddelkoopJC, EhlertPAI. Changes in soil phosphorus pools of grasslands following 17 yrs of balanced application of manure and fertilizer. Soil Use Manag. 2017;33: 212. doi: 10.1111/sum.12333

YangX, ChenX, YangX. Effect of organic matter on phosphorus adsorption and desorption in a black soil from Northeast China. Soil Tillage Res. 2019;187: 8591. doi: 10.1016/j.still.2018.11.016

MessigaAJ, ZiadiN, AngersDA, MorelC, ParentLE. Tillage practices of a clay loam soil affect soil aggregation and associated C and P concentrations. Geoderma. 2011;164: 225231. doi: 10.1016/j.geoderma.2011.06.014

GilbertN. 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: 716719. Available: doi: 10.1038/461716a

ZhuJ, LiM, WhelanM. Phosphorus activators contribute to legacy phosphorus availability in agricultural soils: A review. Sci Total Environ. 2018;612: 522537. doi: 10.1016/j.scitotenv.2017.08.095

WithersPJA, EdwardsAC, FoyRH. Phosphorus cycling in UK agriculture and implications for phosphorus loss from soil. Soil Use Manag. 2006;17: 139149. doi: 10.1111/j.1475-2743.2001.tb00020.x

Gérard-MarchantP, HivelyWD, SteenhuisTS. Distributed hydrological modelling of total dissolved phosphorus transport in an agricultural landscape, part I: Distributed runoff generation. Hydrol Earth Syst Sci. 2006;10: 245261. doi: 10.5194/hess-10-245-2006


LairGJ, ZehetnerF, KhanZH, GerzabekMH. Phosphorus sorption-desorption in alluvial soils of a young weathering sequence at the Danube River. Geoderma. 2009;149: 3944. doi: 10.1016/j.geoderma.2008.11.011


AriasM, Da Silva-CarballalJ, García-RíoL, MejutoJ, NúñezA. Retention of phosphorus by iron and aluminum-oxides-coated quartz particles. J Colloid Interface Sci. 2006;295: 6570. doi: 10.1016/j.jcis.2005.08.001


GouX, CaiY, WangC, LiB, ZhangY, TangX, et al Effects of different long-term cropping systems on phosphorus adsorption and desorption characteristics in red soils. J Soils Sediments. 2020;20: 13711382. doi: 10.1007/s11368-019-02493-2


HeidariS, ReyhanitabarA, OustanS. Kinetics of phosphorus desorption from calcareous soils using DGT technique. Geoderma. 2017;305: 275280. doi: 10.1016/j.geoderma.2017.06.012


McDonaldGJ, NortonSA, FernandezIJ, HoppeKM, DennisJ, AmirbahmanA. Chemical controls on dissolved phosphorus mobilization in a calcareous agricultural stream during base flow. Sci Total Environ. 2019;660: 876885. doi: 10.1016/j.scitotenv.2019.01.059


QuesadaCA, LloydJ, AndersonLO, FyllasNM, SchwarzM, CzimczikCI. Soils of Amazonia with particular reference to the RAINFOR sites. Biogeosciences. 2011;8: 14151440. doi: 10.5194/bg-8-1415-2011


FangH, CuiZ, HeG, HuangL, ChenM. Phosphorus adsorption onto clay minerals and iron oxide with consideration of heterogeneous particle morphology. Sci Total Environ. 2017;605606: 357–367. doi: 10.1016/j.scitotenv.2017.05.133


BornøML, Müller-StöverDS, LiuF. Contrasting effects of biochar on phosphorus dynamics and bioavailability in different soil types. Sci Total Environ. 2018;627: 963974. doi: 10.1016/j.scitotenv.2018.01.283


LiuY, ZhuZQ, HeXS, YangC, DuYQ, HuangYD, et al Mechanisms of rice straw biochar effects on phosphorus sorption characteristics of acid upland red soils. Chemosphere. 2018;207: 267277. doi: 10.1016/j.chemosphere.2018.05.086


ZhangY, HuangS, GuoD, ZhangS, SongX, YueK, et al Phosphorus adsorption and desorption characteristics of different textural fluvo-aquic soils under long-term fertilization. J Soils Sediments. 2019;19: 13061318. doi: 10.1007/s11368-018-2122-0


AhmedW, JingH, KaillouL, QaswarM, KhanMN, JinC, 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 doi: 10.1371/journal.pone.0216881


DebickaM, KocowiczA, WeberJ, JamrozE. Organic matter effects on phosphorus sorption in sandy soils. Arch Agron Soil Sci. 2016;62: 840855. doi: 10.1080/03650340.2015.1083981


BeraR, SealA, BhattacharyyaP, MukhopadhyayK, GiriR. Phosphate sorption desorption characteristics of some ferruginous soils of tropical region in Eastern India. Environ Geol. 2006;51: 399407. doi: 10.1007/s00254-006-0335-9


AhmedW, KailouL, QaswarM, JingH, QinghaiH, YongmeiX, 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 doi: 10.3390/agronomy9120784


BhattacharyyaP, NayakAK, ShahidM, TripathiR, MohantyS, KumarA, 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: 387395. doi: 10.1016/j.cj.2015.03.009


GuedesRS, MeloLCA, VergützL, Rodríguez-VilaA, CoveloEF, FernandesAR. Adsorption and desorption kinetics and phosphorus hysteresis in highly weathered soil by stirred flow chamber experiments. Soil Tillage Res. 2016;162: 4654. doi: 10.1016/j.still.2016.04.018


LiRY, QiuYX, LiuCY, WangYL, HuangLD. Adsorption-desorption behaviors of phosphorus under different silicon concentrations in paddy soils. Chin J Soil Sci. 2013;44: 11341139.


YanX, WangD, ZhangH, ZhangG, WeiZ. Organic amendments affect phosphorus sorption characteristics in a paddy soil. Agric Ecosyst Environ. 2013;175: 4753. doi: 10.1016/j.agee.2013.05.009


ZhangWL, XuAG, ZhangRL, JiHJ, WuSX. Review of soil classification and revision of China soil classification system. Sci Agric Sin. 2014;47: 32143230. doi: 10.3864/j.issn.0578-1752.2014.16.009


WalkleyA. An examination of methods for determining organic carbon and nitrogen in soils. (with one text-figure). J Agric Sci. 1935;25: 598609. doi: 10.1017/S0021859600019687


MurphyJ, RileyJP. A modified single solution method for the determination of phosphate in natural waters. Anal Chim Acta. 1962;27: 3136. doi: 10.1016/S0003-2670(00)88444-5


OlsenSR, ColeC V, WatandbeF, DeanL. Estimation of Available Phosphorus in Soil by Extraction with sodium Bicarbonate. Journal of Chemical Information and Modeling. United States Department Of Agriculture; Washington; 1954 doi: 10.1017/CBO9781107415324.004


LuRK. Soil agricultural chemical analysis method. China Agric Sci Technol Press Beijing 2000; 1315.


XuRK, QafokuNP, Van RanstE, LiJY, JiangJ. Adsorption Properties of Subtropical and Tropical Variable Charge Soils: Implications from Climate Change and Biochar Amendment. Advances in Agronomy. 2016 doi: 10.1016/bs.agron.2015.09.001


JiangJ, XuR, WangY, ZhaoA. The mechanism of chromate sorption by three variable charge soils. Chemosphere. 2008 doi: 10.1016/j.chemosphere.2007.12.012


SuiY, ThompsonML. Phosphorus Sorption, Desorption, and Buffering Capacity in a Biosolids-Amended Mollisol. Soil Sci Soc Am J. 2000;64: 164169. doi: 10.2136/sssaj2000.641164x


LuH, YangL, ShabbirS, WuY. The adsorption process during inorganic phosphorus removal by cultured periphyton. Environ Sci Pollut Res. 2014;21: 87828791. doi: 10.1007/s11356-014-2813-z


BACHEBW, WILLIAMSEG. a Phosphate Sorption Index for Soils. J Soil Sci. 1971;22: 289301. doi: 10.1111/j.1365-2389.1971.tb01617.x


HolfordICR. Soil phosphorus: Its measurement, and its uptake by plants. Australian Journal of Soil Research. 1997 pp. 227239. doi: 10.1071/S96047


QaswarM, JingH, AhmedW, DongchuL, ShujunL, AliS, 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 doi: 10.3390/agronomy9120780


CaiA, ZhangW, XuM, WangB, WenS, ShahSAA. Soil fertility and crop yield after manure addition to acidic soils in South China. Nutr Cycl Agroecosystems. 2018;111: 6172. doi: 10.1007/s10705-018-9918-6


LanZM, LinXJ, WangF, ZhangH, ChenCR. Phosphorus availability and rice grain yield in a paddy soil in response to long-term fertilization. Biol Fertil Soils. 2012;48: 579588. doi: 10.1007/s00374-011-0650-5


FageriaNK, Dos SantosAB, MoraesMF. Influence of urea and ammonium sulfate on soil acidity indices in lowland rice production. Commun Soil Sci Plant Anal. 2010;41: 15651575. doi: 10.1080/00103624.2010.485237


WuQ, ZhangS, ZhuP, HuangS, WangB, ZhaoLP, 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 doi: 10.1371/journal.pone.0176437


WangL, LiangT. Effects of exogenous rare earth elements on phosphorus adsorption and desorption in different types of soils. Chemosphere. 2014;103: 148155. doi: 10.1016/j.chemosphere.2013.11.050


SiddiqueMT, RobinsonJS. Differences in Phosphorus Retention and Release in Soils Amended with Animal Manures and Sewage Sludge. Soil Sci Soc Am J. 2004;68: 14211428. doi: 10.2136/sssaj2004.1421


AbdalaDB, GhoshAK, da SilvaIR, de NovaisRF, Alvarez VenegasVH. Phosphorus saturation of a tropical soil and related P leaching caused by poultry litter addition. Agric Ecosyst Environ. 2012;162: 1523. doi: 10.1016/j.agee.2012.08.004


BhadhaJH, DaroubSH, LangTA. Effect of kinetic control, soil: Solution ratio, electrolyte cation, and others, on equilibrium phosphorus concentration. Geoderma. 2012;173174: 209–214. doi: 10.1016/j.geoderma.2011.12.027


LaiDYF, LamKC. Phosphorus sorption by sediments in a subtropical constructed wetland receiving stormwater runoff. Ecol Eng. 2009;35: 735743. doi: 10.1016/j.ecoleng.2008.11.009


GérardF. Clay minerals, iron/aluminum oxides, and their contribution to phosphate sorption in soils—A myth revisited. Geoderma. 2016;262: 213226. doi: 10.1016/j.geoderma.2015.08.036


WangY, ZhangY, HeY. Effect of soil matrix components on phosphate sorption index in red soil. Acta Ped0L0Gica Sin. 2012;49: 552559.


MoshiAO, WildA, GreenlandDJ. Effect of organic matter on the charge and phosphate adsorption characteristics of kikuyu red clay from kenya. Geoderma. 1974;11: 275285. doi: 10.1016/0016-7061(74)90054-8


LinJ, ZhangZ, ZhanY. Effect of humic acid preloading on phosphate adsorption onto zirconium-modified zeolite. Environ Sci Pollut Res. 2017;24: 1219512211. doi: 10.1007/s11356-017-8873-0


YanJ, JiangT, YaoY, LuS, WangQ, WeiS. Preliminary investigation of phosphorus adsorption onto two types of iron oxide-organic matter complexes. J Environ Sci (China). 2016;42: 152162. doi: 10.1016/j.jes.2015.08.008


YokoyamaD, MoriT, WagaiR, HiradateS, KitayamaK. 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: 218221. doi: 10.1080/00380768.2017.1421018


ZhouQ, ZhuY. Potential pollution and recommended critical levels of phosphorus in paddy soils of the southern Lake Tai area, China. Geoderma. 2003;115: 4554. doi: 10.1016/S0016-7061(03)00074-0


ZhongX, ZhaoX, BaoH, LiH, LiG, LinQ. The evaluation of phosphorus leaching risk of 23 Chinese soils Ⅰ. Leaching criterion. Acta Ecol Sin. 2004;24: 22752280.


WILLIAMSEG, SAUNDERSWMH. Distribution of Phosphorus in Profiles and Particle‐Size Fractions of Some Scottish Soils. J Soil Sci. 1956;7: 90109. doi: 10.1111/j.1365-2389.1956.tb00866.x


WangJB, ChenZH, ChenLJ, ZhuAN, WuZJ. Surface soil phosphorus and phosphatase activities affected by tillage and crop residue input amounts. Plant, Soil Environ. 2011;57: 251257. doi: 10.17221/437/2010-pse


WangS, WangE, others. Desorption characteristics of phosphorus from different used sandy soil in Western Liao River Basin. Res Environ Sci. 2011;24: 756762.


TellenVA, YerimaBPK. Effects of land use change on soil physicochemical properties in selected areas in the North West region of Cameroon. Environ Syst Res. 2018 doi: 10.1186/s40068-018-0106-0
This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. of long-term inorganic and organic fertilization on phosphorus adsorption and desorption characteristics in red paddies in southern China&author=Waqas Ahmed,Huang Jing,Liu Kailou,Sehrish Ali,Han Tianfu,Sun Geng,Chen Jin,Muhammad Qaswar,Du Jiangxue,Sajid Mahmood,Ali Akbar Maitlo,Zulqarnain Haider Khan,Huimin Zhang,Di-Yun Chen,Dafeng Hui,&keyword=&subject=Research Article,Physical Sciences,Chemistry,Physical Chemistry,Sorption,Adsorption,Physical Sciences,Chemistry,Physical Chemistry,Sorption,Desorption,Biology and Life Sciences,Agriculture,Agrochemicals,Fertilizers,Biology and Life Sciences,Agriculture,Agricultural Soil Science,Earth Sciences,Soil Science,Agricultural Soil Science,Earth Sciences,Geography,Cartography,Isotherms,Physical sciences,Chemistry,Environmental chemistry,Soil chemistry,Soil pH,Ecology and environmental sciences,Environmental chemistry,Soil chemistry,Soil pH,Earth sciences,Soil science,Soil chemistry,Soil pH,Physical Sciences,Chemistry,Environmental Chemistry,Soil Chemistry,Ecology and Environmental Sciences,Environmental Chemistry,Soil Chemistry,Earth Sciences,Soil Science,Soil Chemistry,Biology and Life Sciences,Agriculture,Crop Science,Crops,