Biosolids Soil Application: Agronomic and Environmental ImplicationsView this Special Issue
Phosphorus Sorption and Redistribution on Soil Solid Phase in a Brazilian Haplorthox Amended with Biosolids
Land application of biosolids (SS) can cause a buildup of phosphorus (P) in the top soil. The changes in the soil P characteristics may be assessed by the sorption isotherm and the sequential fractionation techniques. Samples of Haplorthox were collected from a field experiment where maize was cultivated for two years, after two applications of SS originated from two cities of São Paulo State, Brazil. SS applications added a total of 125, 250, 500, 1000 and 2000 kg ha−1 of P in the area. To perform the sorption isotherms and obtain P maximum sorption capacity () and the binding energy, soil samples were submitted to increasing P concentration solutions until equilibrium was reached. Sequential fractionation was done by a sequential extraction with CaCl2, NaHCO3, NaOH, HCl, and HNO3 + HClO4 (residual). Addition of biosolids from both cities to the soil decreased and the binding energy obtained by the Langmuir equation. SS additions changed the P fractions distribution in the soil by increasing the labile fractions (P-CaCl2 and P-NaHCO3) and the moderately labile fraction (P-NaOH) by 11.2% and 20.3%, respectively, in detriment of the most resistant P fraction.
A sustainable option to dispose biosolids (SS) originated from sewage treatment plants is the use as soil conditioner or as a source of nutrients to plants. Nevertheless, high SS rates might add to soil excessive amounts of phosphorus (P). One of the consequences is the increase on the risk of environment contamination by P transport to superficial water reservoirs, since P can be weakly sorbed in the superficial soil layer and carried over by rainfall runoff. Furthermore, topsoil P concentrations have been related to P movement below the plough layer and tile flow . Recently, Gebrim et al.  observed that P leaching through the profile can be significant in highly weathered soils, where residual P is accumulated by successive crops, particularly when poultry litter is applied as fertilizer.
Among the useful parameters to monitor the increase in soil P availability is the P sorption capacity and P distribution among the main soil fractions. The application of organic residues to soils might implicate in change in soil P sorption capacity. According to Reddy et al. , application of organic wastes increases soluble P and decreases the P sorption capacity of the soil, and these parameters were directly related to the loading rates of animal wastes. The positive effects of manure on the solubility of P are due to the addition of soluble inorganic P and mineralization of organic P during decomposition of waste, which probably saturated the adsorption sites. Or else, some organic acids released during residue decomposition might compete with P for the same sorption sites in the soil solid phase, decreasing the element immobilization or forming stable complexes with Fe and Al and consequently blocking P retention by them . Nevertheless, it is supposed that frequent SS addition can increase the soil organic matter content which has a positive correlation with the capacity of soil in retaining P .
The sorption isotherms can be useful to detect such alterations. This technique consists of mixing a P-enriched solution with a soil sample, at constant temperature and ionic strength, and graphically plotting the P quantity sorbed by soil versus the P remained in solution. The curve yielded, which although asymptotic, can be described by equations such as the Freundlich or Langmuir. The Langmuir equation can better provide estimates of P sorption maximum and the sorption strength of the soil, parameters inherent in determining the rate of P movement through the soil . The Langmuir equation in its simplest form assumes only monolayer sorption and hence tends to greatly underestimate the P sorption maximum but yields an estimate of P sorption strength that represents an average of the total population of sorption sites . Holford et al.  proposed that to measure a reasonable P sorption maximum, a “two-surface” Langmuir equation should be used. This equation assumes that the total number of sorption sites can be separated into two populations, one with an affinity much greater than the other. According to these authors, the application of this equation to studies of P sorption in soil has revealed much in terms of the relative effect of soil constituents, including P fertilizer and manure application rates and soil fertility.
After SS addition to the soil, the P species present in the material will react in different ways, depending on the soil chemical, physical, or biological characteristics, until the steady state equilibrium. This dynamics can be evaluated through the element redistribution in the main soil fractions. Therefore, the sequential fractionation is an important tool once it is based on the extractor affinity with a specific-soil P species. The method of Heddley et al. , for instance, is able to identify four main P fractions: (a) labile P (P extractable by anionic resin plus P extractable by NaHCO3); (b) moderately labile P ( extractable by NaOH); (c) P bound to Ca (extractable by HCl) and (d) residual P (determined after nitric-perchloric acid digestion).
In Brazil, there is little information about the SS effects on soil P dynamics. Thus, any new available information quantifying and/or identifying some of the soil P dynamic steps will be relevant to improve the efficacy and safety of actual recommendations on SS application in agricultural lands.
The objectives of this work were to evaluate, by means of sorption isotherms, the possible changes in soil P sorption capacity and also to assess P redistribution among the main P fractions in an oxisol that received two successive applications of sewage sludge.
2. Material and Methods
Soil samples were collected from a long-term experiment under field conditions, on an Haplorthox area (dystrophic Red Yellow Oxisol) in the experiment station of Embrapa, at Jaguariúna (22°41′S; 47°06′WGr.; altitude 570 m), State of São Paulo, Brazil, after two SS additions. Some of the soil chemical and physical characteristics before SS additions are given in Table 1. Organic matter (OM) was determined by a colorimetric procedure, after digestion of samples with sodium dichromate solution in sulfuric acid and comparing the readings with a standard curve prepared with soil samples previously analyzed by Walkley-Black’s method . Soil pH (1 : 2.5 solid/solution ratio) was determined in a CaCl2 0.01 mol L−1 solution. In the same suspension, the SMP buffer pH was measured, and with this value, the potential acidity was evaluated . Phosphorus (P), calcium (Ca), magnesium (Mg), and potassium (K) were extracted from the soil by an ion-exchange resin procedure, after a 16-hour shaking period. P was determined by colorimetric method, K by flame photometry, and Ca and Mg by atomic absorption spectrophotometry . The values of sum of bases (S), cation exchange capacity (CEC), and percent base saturation (V) were calculated. Oxalate-extractable Fe and Al were obtained by extraction with ammonium oxalate-oxalic acid  to determine P associated with amorphous oxides of Fe and Al. Physical analysis was performed according to Camargo et al. . The hydrometer method was used for particle-size analysis after dispersion with a mixture of 0.015 mol L−1 (NaPO3)6·Na2O and 0.1 mol L−1 NaOH by overnight shaking. Bulk density was determined on 76 mm × 76 mm diameter undisturbed soil cores .
The main objective was to evaluate the agroenvironmental effects caused by successive SS applications before maize (Zea mays L., CATI AL 30 cultivar) cropping. Hence, in order to evaluate the environmental impacts, such as heavy-metal phytoavailability and accumulation in soil and nitrate leaching, biosolids were also applied in higher amounts than the recommended for optimal crop growth. Biosolids from the cities of Barueri and Franca, both from the State of São Paulo, were used in the experiment. Some SS characteristics are presented in Table 2. Nitrogen was determined in the original SS-samples using the Kjeldahl method (vapor distillation), and the organic-C was determined by titrimetry after dichromate digestion in digester block, according to Raij et al. . Phosphorus, sulfur, calcium, magnesium, and potassium were analyzed according to the recommendations of US-EPA (SW-846) method no. 3051  K by flame photometry, and other elements by ICP-AES. Humidity was determined by mass loss at 60°C and SS sample pH was determined in water extract (1 : 5).
The experiment consisted of twelve treatments in a randomized complete block design with three replications, in a total of 36 plots, each plot measuring 20 m length versus 10 m width. The study treatments were: two control (no N, P and K applied); two mineral fertilization recommended for the crop (46 kg N ha−1 for the first crop and 90 kg N ha−1 for the second crop; 37 kg P ha−1 and 56 kg K ha−1 on both crops) ; four treatments with the SS from the city of Franca, four treatments with the SS from the city of Barueri. The sewage sludge treatments were applied based on the N concentration that provides the same amount of N as in the mineral fertilization recommended for corn; two, four, and eight times the N recommended dosage for corn crop. Calculations of sludge rates were performed as a function of the N available for plants and were calculated by the formula where Kjeldahl-N: total nitrogen (mg kg−1), NH3-N: nitrogen as ammonia (mg kg−1), -N: nitrogen as nitrate (mg kg−1), -N: nitrogen as nitrite (mg kg−1), and MF: mineralization fraction (30%), for aerobically digested sludge.
Supplementary K was applied for treatments with sewage sludge, when necessary. The wet sludge was incorporated on April 1999 and on December 1999. It was toss-distributed in the total area of the experimental plots and incorporated to a depth of 20 cm with a rotary harrow, 3-4 days before sowing.
The P quantity in each SS rate was calculated based on the SS total P content added after two SS applications. Treatments applied were: FC = control without SS (for the plots with SS from Franca); BC = control without SS (for the plots with SS from Barueri); FP1, FP2, FP4, FP8 = 125, 250, 500, and 1000 kg ha−1 of P from Franca-SS applications; BP2, BP4, BP8, and BP16 = 250, 500, 1000, and 2000 kg ha−1 of P from Barueri-SS applications; FNPK and BNPK = mineral fertilized plots.
Soil samples were collected at 0–20 cm depth, after the second maize harvesting, air dried and passed through 2 mm sieve. For the P sorption essay, 3 g of soil samples of each SS- treated plot with 0, 125 (FP1), 250 (BP2), 500 (FP4), 1000 (FP8), and 2000 (BP16) kg ha−1 of P were equilibrated in a 50 mL centrifuge tubes with 30 mL of 0.01 mol L−1 CaCl2 solution containing 0, 5, 10, 15, 30, 60, 100 and 250 mg L−1 of P. Two drops of toluene were added to preserve solution from microorganism contamination. During six days, solutions were shaken for 30 min thrice a day. Finally, the tubes were centrifuged at 577 ×g for 15 min to separate supernatant from the solid material, and P was determined in the suspension . The difference between the P in the equilibrium solution and P added is the quantity of P adsorbed by the soil. The adsorption isotherms, which represent the linear fit of the relationship between adsorbed P and the equilibrium solution P, were built with the results obtained with all P concentrations. Sorption strengths were calculated with the simple (monolayer) Langmuir equation, which gives an average sorption strength of the P sorbed , where is the P sorbed (mg kg−1), is the sorption strength (L mg−1 P), is the maximum quantity of P sorbed for monolayer P sorption (mg kg−1), and is the concentration of P in solution (mg L−1).
The total P sorption maximum (: mg kg−1) was calculated as the sum of P sorption maxima ( and ) derived from the two-surface Langmuir equation . For the two-surface Langmuir equation, where is the P sorbed (mg kg−1), and are the sorption strengths of the first and the second sorption layers (L mg−1 P), and are the maximum quantities of P sorbed for two sorption layers (mg kg−1), and is the concentration of P in solution (mg L−1).
The sequential fractionation was based on the method of Heddley et al. , modified by Sui , and also with the following modifications: 0.01 mol L−1 CaCl2 solution instead of H2O, centrifugation at 1073 ×g at 5°C, and filtering through Whatman 41 paper (20–25 μm) instead of 0.22 μm. Samples were selected only from Barueri SS treatments that added 0, 250 (BP2), 1000 (BP8), and 2000 (BP16) kg P ha−1. An aliquot of 0.5 g of each sample was submitted to nitric-perchloric digestion to determine the total P (Pt) concentration by Inductively Coupled Plasma Atomic Emission Spectroscopy (ICP-AES). Another equal aliquot sample was transferred to 50 mL centrifuge tubes and submitted to P extraction using a sequence of P extracting solutions: 30 mL of 0.01 mol L−1 CaCl2 solution; 30 mL of 0.5 mol L−1 NaHCO3 (pH 8.2); 30 mL 0.1 mol L−1 NaOH; 30 mL of 1.0 mol L−1 HCl. Each soil-solution mixture was shaken for 16 hours, centrifuged at 3,000 g for 15 minutes, and filtered through Whatman 41 filter paper. The inorganic phosphorus (Pi) was determined by colorimetric method described by Murphy and Riley , and the organic phosphorus (Po) was determined by the difference between Pt and Pi in the same extract. After the last extracting solution, the soil residue remained in each centrifuge tube was submitted to nitric-perchloric digestion, and P content was determined by ICP-AES.
This procedure of sequential fractionation allowed identifying four types of P bindings to the solid phase: (a) labile or biologically available P, extracted by CaCl2 solution and by NaHCO3 solution; (b) chemisorbed or moderately labile Pi and Po (bound to Fe and Al oxides) extracted by NaOH solution; (c) P (bound to Ca) extracted by HCl solution; (d) residual P.
Langmuir monolayer equations were evaluated on the basis of statistical significance of the regression coefficients.
3. Results and Discussion
Sewage sludge addition from both cities (Franca and Barueri) decreased soil P sorption capacity in all rates under investigation (Figures 1(a) and 1(b)). Such behavior was expected since SS usually contains water soluble P , and the SS tested here presented more than 0.3% total P  and a C : P ratio less than 100  which also favors the decrease in P sorption by releasing more P into soil solution as a result of organic P mineralization. Berton and Pratt  also used the isotherms to evaluate the effects of increasing rates of four different organic materials (including SS) on the P adsorption capacity of two soils and concluded that SS decreased P adsorption more effectively compared to the other materials.
The linear form of Langmuir regression equations presented a significant for all SS rates applied indicating that the sorption phenomena was adequately described by this isotherm (Table 3). P sorption maxima ranged from 1130 to 1008 mg kg−1 and the binding energy also varied from 0.148 to 0.072 L mg−1 for control and treatment 2000 kg ha−1, respectively (Table 3). Both and the binding energy values decreased as a result of soil SS incorporation, but this effect was linear only for the residue derived from Barueri. Although soil usually presents a positive correlation with soil organic matter content , the addition of organic residues to soil, such as manures, causes an inverse effect by decreasing P sorption with increase in P availability to plants . Mechanisms for the reduction of sorption maxima and the affinity constants by organic amendments include the competition of organic acids produced during mineralization for the same sites of P fixation and the complexing of exchangeable Al by organic acids [4, 25, 26]. Furthermore, both sewage sludges contain inorganic ligands such as which can complex exchangeable Al and metals such as Ca and Mg which may form soluble complexes with P in the soil solution, thus preventing P from being fixed by Al or Fe .
Only the treatment with the highest P input (2000 kg ha−1) reached the threshold value of 0.07 L mg−1 proposed by McDowell et al. , whereby soils less than this value for binding strength are subject to enhance loss of P in subsurface flow. Since the recommended SS rate is eight times less than the highest rate tested (125 kg ha−1 and 250 kg ha−1 for Franca and Barueri cities, resp.), it indicates that the amount of SS recommended as source for corn crop may be considered safe for potential P losses in subsurface flow, but it should be monitored after successive SS applications.
The amount of P (mg kg−1) in equilibrium with 0.2 mg P L−1 (P0.2) calculated from each isotherm equation (Table 3) ranged from 32,49 to 14,31 mg kg−1 for control and treatment 2000 kg ha−1, respectively. The amount of P sorbed at P0.2 in the equilibrium solution was decreased by increasing rates of SS from both cities and markedly reduced by the highest rates of these amendments. Sewage sludge from Franca required the double amount of added P to attain P0.2 in equilibrium solution compared to the sludge from Barueri (500 to 250 kg ha−1), suggesting that the latter was more efficient in maintaining P in soil solution.
The SS from Barueri city applied to soil changed the P distribution within the studied P fractions (Table 4). Phosphorus extracted by CaCl2 was markedly increased only at P rate of 2000 kg ha−1. The P rates applied as SS had little effect over PoNaHCO3. In relative values, this P fraction decreased from 10.8 to 4.9% as long as P rates increased from 250 kg ha−1 to 2000 kg ha−1 (Table 5). This apparent little variation of Po-NaHCO3 fraction might probably be attributed to an intense organic P mineralization which could have contributed to increase Pi-NaHCO3 fraction from 9.1% (control plot) to 18.1% (highest P rate-2000 kg ha−1). Similar results were found by Sui et al.  who suggested that Po-NaHCO3 fraction might be responsible for an extra increase in Pi-NaHCO3 fraction. Nevertheless, the Po from the Po-NaHCO3 fraction might also be adsorbed to Al and Fe oxides, favoring the increasing Po-NaOH fraction from 1.5 % to 11.3 % (Table 5).
Compared with any other P fractions, the contribution of P added as SS to the NaOH fraction was most evident and illustrates the P affinity to Fe and Al oxides, a relationship that predominates in oxisols. The increases of chemisorbed P observed on all SS treatments were likely the result of soluble inorganic P added with the organic residue and P mineralized from added organic P . Since the NaOH extractant removes a P fraction which is associated with amorphous and crystalline Al and Fe phosphates [29, 30], the addition of SS appears to have a significant impact on soil P sorption, as showed by the sorption isotherms (Figure 1(b)).
Addition of SS increased P-HCl fraction, but it was not proportional to the rates of SS applied. Novais and Smyth  suggested that the low variation of PHCl fraction with the increasing P rates applied to soil evidenced that P bound to Ca did not have any important role in P dynamics of more weathered acid soils and not treated with natural phosphates as P source. In these soils, the poor calcium phosphate content was probably the main factor determining the HCl-extractable P quantity. In a similar way, Rheinheimer  observed in two oxisols that P extracted with 1.0 mol L−1 HCl (after resin, 0.5 mol L−1 NaHCO3 and 0.1 mol L−1 NaOH extractions) did not change with the successive extractions.
The effect of SS on soil P species distribution, increasing the element availability, can be observed by the variation occurred in three P fractions: (a) the labile fraction represented by PCaCl2 plus PNaHCO3 increased from 13.2% (control) to 24.4% of total P; (b) the chemisorbed P or moderately labile fraction suggested by Bowman and Cole , represented by P forms extracted by 0.1 mol L−1 NaOH increased from 9.6 % (control) to 29.9 % of total P; consequently, (c) the most resistant residual P fraction decreased from 73.4 % (control) to 41.4 % (Figure 2). The redistribution of residual P to more labile forms was also observed by Barcellos  on an eutroferric red latosol amended with successive applications of SS.
(1) The biosolids application to soil decreased soil P sorption which resulted in an increase of P in equilibrium solution.
(2) The positive correlation of P extracted by NaOH with total P added as SS suggested that the decrease in soil P sorption can be mainly attributed to the P chemisorbed to amorphous and crystalline Al and Fe.
(3) The addition of biosolids changed the P fractions distribution in the soil, favoring P availability to plants and its geoavailability by increasing the labile fractions (P-CaCl2 and P-NaHCO3) and the moderately labile fraction (P-NaOH) by 11.2 % and 20.3 %, respectively, in detriment of the most resistant P fraction.
R. McDowell, A. Sharpley, and P. Withers, “Indicator to predict the movement of phosphorus from soil to subsurface flow,” Environmental Science and Technology, vol. 36, no. 7, pp. 1505–1509, 2002.View at: Publisher Site | Google Scholar
F. O. Gebrim, R. F. Novais, I. R. Silva et al., “Mobility of inorganic and organic phosphorus forms under different levels of phosphate and poultry litter fertilization in soils,” Revista Brasileira de Ciência do Solo, vol. 34, no. 4, pp. 1195–1205, 2010.View at: Google Scholar
K. R. Reddy, M. R. Overcash, R. Khaleel, and P. W. Westerman, “Phosphorus adsorption-desorption characteristics of two soils utilized for disposal of animal wastes,” Journal of Environmental Quality, vol. 9, no. 1, pp. 86–92, 1980.View at: Google Scholar
S. Nagarajah, A. M. Posner, and J. P. Quirk, “Competitive adsorption of phosphate with polygalacturonate and other organic anions on kaolinite and oxide surfaces,” Nature, vol. 228, no. 5266, pp. 83–85, 1970.View at: Publisher Site | Google Scholar
S. K. Sanyal and S. K. De Datta, “Chemistry of phosphorus transformations in soil,” Advances in Soil Sciences, vol. 16, pp. 1–120, 1991.View at: Google Scholar
I. C. R. Holford, C. Hird, and R. Lawrie, “Effects of animal effluents on the phosphorus sorption characteristics of soils,” The Australian Journal of Soil Research, vol. 35, no. 2, pp. 365–373, 1997.View at: Google Scholar
I. C. R. Holford, R. W. M. Wedderburn, and G. E. G. Mattingly, “A Langmnuir two-surface equation as a model of phosphate adsorption by soils,” Journal of Soil Science, vol. 25, pp. 242–254, 1974.View at: Google Scholar
M. J. Heddley, J. W. B. Stewart, and B. S. Chauhuan, “Changes in organic and organic soil phosphorus fractions induced by cultivation practices and laboratory incubations,” Soil Science Society of America Journal, vol. 46, pp. 970–976, 1982.View at: Google Scholar
A. Walkley and I. A. Black, “An examination of the Degtjareff method for determining organic carbon in soils: effect of variations in digestion conditions and of inorganic soil constituents,” Soil Science, vol. 63, pp. 251–263, 1934.View at: Google Scholar
H. E. Shoemaker, E. O. McLean, and P. F. Pratt, “Buffer methods for determining lime requirement of soils with appreciable amounts of extractable aluminum,” Soil Science Society of America Proceedings, vol. 25, pp. 274–277, 1961.View at: Google Scholar
B. van Raij, J. A. Quaggio, and N. M. Silva, “Extraction of phosphorus, potassium, calcium and magnesium from soils by an ion-exchange resin procedure,” Communications in Soil Science and Plant Analysis, vol. 14, no. 5, pp. 547–566, 1986.View at: Google Scholar
L. J. Evans and W. G. Wilson, “Extractable Fe, Al, Si and C in B-horizons of podzolic and brunisolic soils from Ontario,” The Canadian Journal of Soil Science, vol. 65, no. 3, pp. 489–496, 1985.View at: Google Scholar
O. A. Camargo, A. C. Moniz, T. A. Jorge, and J. M. A. S. Valadares, Métodos de Análise Química, Mineralógica e Física de Solos do Instituto Agronômico de Campinas, Instituto Agronômico, Campinas, Brazil, 1986.
R. E. Uhland, “Physical properties of soils as modified by crops and management,” Soil Science Society of America Proceedings, vol. 14, pp. 361–366, 1950.View at: Google Scholar
B. van Raij, J. A. Quaggio, H. Cantarella, and C. A. Abreu, “Os métodos de análise química do sistema I.A.C. de análise de solo no contexto nacional,” in Análise Química para Avaliação da Fertilidade de Solos Tropicais, B. van Raij, J. C. Andrade, H. Cantarella, and J. A. Quaggio, Eds., pp. 5–39, Instituto Agronômico, Campinas, Brazil, 2001.View at: Google Scholar
US-EPA, Test Methods for Evaluating Solid Waste: Physical/Chemical Methods, U. S. Environmental Protection Agency/Office of Solid Waste and Emergency Response/U. S. Government Printing Office, Washington, DC, USA, 3rd edition, 1983, (SW-846, Method 3051).
B. van Raij, H. Cantarella, J. A. Quaggio, and A. M. C. Furlani, Recomendações de Adubação e Calagem para o Estado de São Paulo, Instituto Agronômico/ Fundação IAC, Campinas, Brazil, 2nd edition, 1997.
J. Murphy and J. P. Riley, “A modified single solution method for the determination of phosphate in natural waters,” Analytica Chimica Acta, vol. 27, pp. 125–151, 1962.View at: Google Scholar
I. C. R. Holford and G. E. G. Mattingley, “The high- and low-energy phosphate adsorbing surfaces in calcareous soils,” Journal of Soil Science, vol. 26, pp. 407–417, 1975.View at: Google Scholar
Y. Sui, M. L. Thompson, and C. Shang, “Fractionation of phosphorus in a Mollisol amended with biosolids,” Soil Science Society of America Journal, vol. 63, no. 5, pp. 1174–1180, 1999.View at: Google Scholar
R. S. Berton and P. F. Pratt, “Evaluation of phosphorus requirement by the sorption isotherm technique in soils amended with organic materials,” Revista Brasileira de Ciência do Solo, vol. 21, no. 2, pp. 199–206, 1997.View at: Google Scholar
B. B. Singh and J. P. Jones, “Phosphorus sorption and desorption characteristics of soil as affected by organic residues,” Soil Science Society of America Proceedings, vol. 40, pp. 389–395, 1976.View at: Google Scholar
F. Iyamuremye and R. P. Dick, “Organic amendments and phosphorus sorption by soils,” Advances in Agronomy, vol. 56, pp. 139–185, 1996.View at: Publisher Site | Google Scholar
C. A. M. Laboski and J. A. Lamb, “Impact of manure application on soil phosphorus sorption characteristics and subsequent water quality implications,” Soil Science, vol. 169, no. 6, pp. 440–448, 2004.View at: Google Scholar
N. V. Hue, “Correcting soil acidity of a highly weathered ultisol with chicken manure and sewage sludge,” Communications in Soil Science and Plant Analysis, vol. 23, no. 3-4, pp. 241–264, 1992.View at: Google Scholar
A. H. Bumaya and D. V. Naylor, “Phosphorus sorption and extractability in Andic soil incubated with plant residues of variable P content,” Plant and Soil, vol. 112, no. 1, pp. 77–81, 1988.View at: Publisher Site | Google Scholar
F. Iyamuremye, R. P. Dick, and J. Baham, “Organic amendments and phosphorus dynamics: I. Phosphorus chemistry and sorption,” Soil Science, vol. 161, no. 7, pp. 426–435, 1996.View at: Publisher Site | Google Scholar
F. Iyamuremye, R. P. Dick, and J. Baham, “Organic amendments and phosphorus dynamics: II. Distribution of soil phosphorus fractions,” Soil Science, vol. 161, no. 7, pp. 436–443, 1996.View at: Publisher Site | Google Scholar
J. R. McLaughlin, J. C. Ryden, and J. K. Syers, “Development and evaluation of a kinetic model to describe phosphate sorption by hydrous ferric oxide gel,” Geoderma, vol. 18, no. 4, pp. 295–307, 1977.View at: Google Scholar
R. A. Bowman and C. V. Cole, “An exploratory method for fractionation of organic phosphorus from grassland soil,” Soil Science, vol. 125, pp. 95–101, 1978.View at: Google Scholar
R. F. Novais and T. J. Smyth, Fósforo em Solo e Planta em Condições Tropicais, Universidade Federal de Viçosa, Viçosa, Brazil, 1999.
D. S. Rheinheimer, I. Anghinoni, and J. Kaminski, “Depleção do fósforo inorgânico de diferentes frações provocada pela extração sucessiva com resina em diferentes solos e manejos,” Revista Brasileira de Ciência do Solo, vol. 24, pp. 345–354, 2000.View at: Google Scholar
I. F. Barcellos, Partição de fósforo (P) em sedimento de solo tratado com lodo de esgoto e cultivado com milho, M.S. thesis, IAC, Campinas, Brazil, 2006.