Next Article in Journal
Application of OpenAir and AgDRIFT Models to Estimate Organophosphate Pesticide Spray Drift: A Case Study in Macon County, Alabama
Next Article in Special Issue
Multi-Sample Detection of Soil Nitrate Nitrogen Using a Digital Microfluidic Platform
Previous Article in Journal
The Efficacy of 5-Aminolevulinic Acid-Producing Luteovulum sphaeroides Strains on Saline Soil Fertility, Nutrient Uptakes, and Yield of Rice
Previous Article in Special Issue
Balancing Yield and Antioxidant Capacity in Basil Microgreens: An Exploration of Nutrient Solution Concentrations in a Floating System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Phosphorus Balance in Sandy Soil Subjected to 12 Years of Successive Applications of Animal Manure and Mineral Phosphate Fertilizer in Subtropical Climate

by
Carina Marchezan
1,†,
Paulo Ademar Avelar Ferreira
1,
Gustavo Boitt
2,
Natalia Moreira Palermo
1,
Ana Laura Thoma
1,
Ricardo Fagan Vidal
1,
Gustavo Scopel
1,
Cledimar Rogério Lourenzi
3,
Carlos Alberto Ceretta
1 and
Gustavo Brunetto
1,*
1
Department of Soil Science, Federal University of Santa Maria, Santa Maria 97105-900, RS, Brazil
2
Soil Science Building-Soils West, UWA School of Agriculture and Environment, University of Western Australia, Perth, WA 6009, Australia
3
Center of Agricultural Science, Rural Engineering Department, Federal University of Santa Catarina (UFSC), Florianópolis 88034-001, SC, Brazil
*
Author to whom correspondence should be addressed.
This work is part of the doctoral thesis Phosphorus balance in 12 years and annual crops strategies to access phosphorus reserves built over 15 years by organic and mineral sources, presented at the Federal University of Santa Maria.
Agriculture 2023, 13(9), 1762; https://doi.org/10.3390/agriculture13091762
Submission received: 20 July 2023 / Revised: 27 August 2023 / Accepted: 28 August 2023 / Published: 5 September 2023
(This article belongs to the Special Issue Advances in Nutrient Management in Soil-Plant System)

Abstract

:
Assessing the phosphorus (P) balance in agricultural soils is crucial for optimizing its use and reducing contamination risks. The objective of this study was to evaluate the impact of different animal wastes on the distribution of P in soil profiles subjected to 12 years of successive applications of animal wastes and mineral phosphate fertilizer, within a crop rotation system under no-till system. The study was conducted from 2004 to 2016 in the southern region of Brazil. The treatments were the applications of pig slurry (PS), cattle slurry (CS), pig deep-litter (PD), mineral fertilizer (MF), and a control treatment without application. The highest accumulation of P and its movement was observed in the 0–40 cm layer, in the soil submitted to applications of all P sources. The inputs of P via MF, CS, PS, and PD promoted the accumulation of 18, 42, 48, and 100 kg P ha−1 year−1. The P mass balance showed that between 77 and 98% of the P added by animal manure and MF was accounted for in grain exports (17–34%), soil storage (41–72%), and post-harvest residues (<1%), with the remaining 2–33%, unaccounted for, which was attributed, especially, to P transfer at the soil surface.

Graphical Abstract

1. Introduction

Phosphorus (P) is one of the most important nutrients for maintaining or increasing the productivity of plants in agroecosystems. The demand for this nutrient for food and bioenergy production has intensified over the years as the world’s population increases, raising concerns about its supply since P sources are limited and non-renewable [1,2]. Global swings in the price of phosphate fertilizers happen to accompany the depletion and scarcity of phosphate rock [3,4]. Therefore, it is critical to find viable alternatives that make it possible to meet food demand without relying solely on mineral phosphate fertilizers. An alternative could be the use of liquid animal waste, such as from pigs and cattle, or solid, such as pig deep-litter, as P sources in cropping systems [5,6].
The states in the southern region of Brazil (Rio Grande do Sul, Santa Catarina, and Paraná) are responsible for approximately 70% of Brazil’s swine production and export more than 84% of the meat from these animals, with 695,000 tons of meat (18.5% of total national production) produced in Rio Grande do Sul [7]. Moreover, Rio Grande do Sul contributed 13% of the national milk production (812 thousand liters in the first quarter of 2019) [8]. These activities, when conducted in confined system, which represents almost the totality for swine production, generate a large volume of waste rich in nutrients, including P. This large volume of waste is used in agricultural soils as an alternative of disposal, but also as a way to reuse nutrients. However, successive applications of manure as fertilizer, in amounts commonly greater than those exported by crops, cause accumulation of P in soils over time [9].
However, long-term studies quantifying the role of legacy P in sustaining agricultural systems for food production are still scarce, especially in soils with higher adsorption capacity, such as those in tropical and subtropical regions. A field trial was established in Rio Grande do Sul, Brazil, in 2004 to evaluate the agronomic responses of grain with different types of animal wastes, along with changes in soil properties and associated potential environmental impacts. The objective of this study was to quantify the accumulation, distribution, and indicators of potentially adverse environmental impacts following the application of different types of animal manure and mineral phosphate fertilizer in a long-term no-till soil. Due to the surface application and high adsorption capacity of the soil, most of the P added via animal manure and mineral fertilizer was retained in the topsoil layers.

2. Materials and Methods

2.1. Location of Experiment and Treatments

The study was performed in a long-term experiment, with 19 years of conduction, located in the experimental area of the Federal University of Santa Maria (UFSM), in Santa Maria, Rio Grande do Sul (RS), in the southern region of Brazil (29°43’12″S and 53°43’4″W). The climate of the region is humid subtropical (Cfa 2), according to Koppen’s classification, with annual averages of temperature, precipitation, and relative humidity of 19.3 °C, 1561 mm, and 82%, respectively. The soil is classified as Arenic Dystrophic Red Argissolo, according to Embrapa [10], corresponding to Typic Hapludalf [11], with 108 g kg−1 of clay, 183 g kg−1 of silt, and 709 g kg−1 of sand. The chemical properties of the soil before the installation of the experiment (2004) and at the time of the study (2016) are presented in Table 1.
The experiment was initiated in 2004, under a no-till farming system. The experimental design was randomized blocks, with four replicates. The plots had dimensions of 5 × 5 m The treatments were pig slurry (PS), cattle slurry (CS), pig deep-litter (PD), mineral fertilizer (urea + triple superphosphate—SFT + potassium chloride—KCl) (MF) and a control treatment, without the application of nutrients. The applications of PS, CS, PD, and mineral fertilizer were made on the surface, over the cultural residues of the previous crop and without incorporation. The treatments were applied only once a year and before summer crop sowing, in the period from 2004 to 2009. The treatments were not applied in winter. Starting in 2010, the treatments were applied twice each year, prior to winter and summer crop establishment. The crop succession used until August 2016 was black oats (Avena strigosa S.), corn (Zea mays), forage turnips (Raphanus sativus), black beans (Phaseolus vulgaris), and wheat (Triticum spp.). Further details about the experiment can be obtained in Lourenzi et al., Ferreira et al., Rodrigues et al., and Marhezan et al., [12,13,14].

2.2. Collecting Soil Samples

In June 2016, soil samples were collected in the layers 0–4, 4–10, 10–20, 20–30, 30–40, and 40–50 cm. For this, in each plot, two trenches with dimensions of 40 × 40 × 40 cmwere opened. In each trench, two soil samples per plot were collected, one in the longitudinal direction and the other transversal to the sowing line. To correct the soil attributes and concentrations in units of mass per area in hectares, in each sampled layer, the bulk density of the soil was determined by means of metal volumetric rings (78 cm3), as described by Embrapa [15].

2.3. P Fractionation in Soil

The soil was air dried, ground, and passed through a 2 mm mesh sieve. It was then prepared and subjected to sequential P chemical fractionation with ammonium chloride (NH4 Cl 1 mol l−1), sodium bicarbonate (NaHCO3 0.5 mol l−1, pH 8.5), sodium hydroxide (NaOH 0.1 mol l−1), hydrochloric acid (HCl 1 mol l−1), and a second extraction with NaOH (0.1 mol l−1) [16,17,18]. Total soil P was determined in a soil subsample after digestion with concentrated sulfuric acid (H2SO4) and hydrogen peroxide (H2O2 at 30% v/v), according to Olsen and Sommers [19]. The inorganic P in the alkaline extracts of NaHCO3 and NaOH was assigned to the reactive P-molybdate analyzed by the method of Dick and Tabatabai [20]. In the alkaline extracts, total P was determined by digestion with ammonium persulfate + sulfuric acid in an autoclave. Unreacted P-molybdate was considered as organic P and was obtained by the difference between total P and inorganic P. The P in acid extracts was determined according to the method of Murphy and Riley [21].
Total organic P extracted was determined as the sum of the organic P extracted with NaHCO3 and the first and second NaOH extractions. Total extracted inorganic P was calculated as the sum of inorganic P in all extracted fractions. The difference between total soil P and total extracted P (inorganic + organic) was considered residual P. Soil P fractions were grouped into four P pools according to their relative lability: labile P (sum of NH4 Cl and NaHCO3 fractions (inorganic and organic)); moderately labile P (first NaOH (inorganic and organic)); stable P (HCl + second NaOH (inorganic and organic)); and residual P (non-extractable P) [22]. The amounts of total P, total inorganic P, and total organic P up to 50 cm were expressed in units of kg ha−1.
The potential risk of P transfer was evaluated by P concentration values (mg P kg−1) determined in the first extraction of sequential P fractionation (NH4 Cl-P). For this, we considered the most superficial layers (0–4 and 4–10 cm), following the recommendations of McDowell and Sharpley [23]. Easily extractable P in soil can be used as an indicator of environmental risk, because this P fraction has a high potential for transferring from the soil by erosion processes and surface runoff [24,25,26]. When the reactive functional groups of soil surfaces can no longer adsorb P, the equilibrium point is established and, as a consequence, large amounts of weakly bound P will be present in the soil solution (in this study represented by NH4 Cl-P), thus increasing potential P transfers [23].

2.4. P Added by Manure and Mineral Fertilizer

Liquid waste derived from pig and cattle farming was collected in open effluent tanks. The overlapping pig litter was collected in covered environments on the pig farms. The general physicochemical properties of the manure were analyzed on a representative sample before each application on all crops was performed. The applied dose of each organic waste was defined based on its N concentration and the N requirement needed by the crops, according to the regional recommendation of the Commission for Soil Chemistry and Fertility CQFS-RS/SC [27]. Therefore, the P inputs varied among the treatments. The total P inputs in the manure were calculated based on the P concentrations (kg m−3 for liquid manure and kg kg−1 for solid manure), present in the manure, multiplied by the amount of manure applied (m3 manure ha−1 for CS and PS and kg manure ha−1 for PD) and, therefore, were expressed in kg P ha−1. For MF, the P inputs were made according to the P contents present in the soil and the nutritional requirement of the crop. The P inputs (in kg P ha−1) in the manure and MF, by crop season/year are presented in Figure 1.

2.5. P in Residue and Exported by Crops

Plants were collected at flowering (winter and summer crops). For this, the shoot of the plants was cut close to the soil surface. Plants were dried in a forced air oven at ±65 °C to constant mass for dry matter determination. The tissue was ground, prepared, and subjected to nitric-perchloric digestion (ratio 4:1) [15]. P was determined in the extract according to Murphy and Riley [21]. This procedure was repeated for the determination of P present in grains and cereals. For the calculation of the P balance in the system, the Presidue, consists of the P present in the cultural residues after harvest. For this, we used the data of P contained in the dry matter of the oat crop established at the time of soil collection.

2.6. P Balance

The P balance in the soil–plant system (kg ha–1) was estimated considering the 50 cm layer, the amount of the nutrient added in fertilization (kg ha–1), the amount of P exported by grain crops (kg ha–1), and the change in the amount of P in the soil. To determine the P budget of the system, we used a similar approach to that described by Boitt et al. [28].
The P mass balance was calculated empirically as follows (units in kg P ha−1) (Equation (1)):
P a d d e d P s o i l P e x p o r t e d P r e s i d u e P u n a c c o u n t e d f o r = 0
where: Padded is the cumulative sum of P added by the treatments; ∆Psoil is the change in P storage in the soil profile (PinitialPfinal) calculated by summing the amounts of P determined in each soil layer (0 to 50 cm layer); Pexported is P exported by grains and cereals; Presidue is the P content of crop residues and Punaccounted for is the difference between Padded and the sum of Pcounted for pools, which was attributed to some accumulated uncertainties in P accounting in the mass balance. For the control treatment, Padded was zero and Punaccounted was considered insignificant, given the absence of P addition and minimal soil losses through erosion and runoff in this experiment [29]. Therefore, the soil from the control treatment was used to calculate Pinitial and was considered equal throughout the experiment area. This made it possible to determine the ΔPsoil of all treatments.
Total P inputs accounted for (Paccounted; Equation (2)) in the treatments fertilized with manure and mineral fertilizer are given by the difference in the sum of all P pools (given by P P o o l ; Equation (3)), minus the sum of the P pools of the control treatment (no P inputs) as follows:
P a c c o n t e d k g h a 1 = P p o o l s t r e a t . f e r t i l i z e d P p o o l s c o n t r o l
Being,
P p o o l s = P s o i l + P e x p o r t e d + P c r o p   r e s i d u e
Thus, the P accounted for in the system (for each treatment that received fertilization) was calculated according to Equation (4):
P a c c o n t e d % = P a c c o n t e d [ k g h a 1 ] P i n p u t [ k g h a 1 ] × 100
The cumulative agronomic efficiency of the P balance of the system was calculated according to Equation (5):
P   b a l a n c e   e f f i c i e n c y % = P e x p o r t e d P e x p o r t e d + P s o i l × 100
where ∆Psoil, in units of kg P ha−1, is the change in soil P storage in the 0 to 50 cm soil layer calculated using Equation (6):
P s o i l k g h a 1 = P f i n a l P i n i t i a l

2.7. Statistical Analysis

Differences in mean soil P fractions between treatments and depths, as well as crop yields, were determined using ANOVA and when significant were compared by the Scott–Knott 5% test. In addition, linear correlations (Pearson correlation coefficients) between the data were determined using R software [30].

3. Results

3.1. P in Soil

The additions of CS, PS, PD, and MF over 12 years increased the amounts of P in all fractions (Figure 2 and Supplementary Table S1). The highest values of labile P were observed in the superficial layers in soils submitted to animal manure applications (CS, PS, PD). However, among the manure sources, the highest values of labile P up to 40 cm layer were observed in the soil with PD applications (Figure 2a). The moderately labile P fractions increased up to 40 cm, especially in soils with animal manure applications, and the highest amounts were observed in soils with PD application (Figure 2b). The highest values of stable P were observed up to 30 cm, in soils submitted to residue and mineral fertilizer applications (Figure 2c). The highest values of residual P (not extractable by chemical P fractionation, Figure 2d) were observed up to 30 cm in the control and MF soils.
The highest values of inorganic P were observed in the soil fertilized with PD, followed by PS, CS, and MF, with increases of 73, 54, 50, and 40%, relative to the control soil (0–50 cm layer) (Figure 3a). The additions of PD, CS, PS, and MF increased surface inorganic P contents by 79, 69, 65, and 52%, respectively, relative to the control (0–10 cm layer), with observed amounts greater than 580 kg P ha−1 in the soil with PD addition. The contribution of inorganic P to total soil P varied from 26% in the control soil to 54, 41, 39, and 37% in the soils with PD, PS, CS, and MF applications, respectively (layer 0–50 cm, Figure 3a,c and Supplementary Table S2). On the other hand, the highest amounts of organic P were observed in soil with PD, followed by CS, PS, and MF, with increments of 36, 30, 29, and 6%, relative to the control (layer 0–50 cm Figure 3b). The additions of PD, PS, CS, and MF increased the amounts of organic P on the surface with values 50, 47, 43, and 42% higher than the control (layer 0–10 cm). The proportion of soil organic P varied from 37% in the soil with PD and MF additions, to 42% in the control soil and with PS and CS applications (layer 0–50 cm, Figure 3b,c and Supplementary Table S2).
The increase of total P and its distribution in the soil profile were proportional with the amounts of P added by the treatments. The highest amounts of total P were observed in the soil with PD applications, followed by PS, CS, and FM, with increases of 79, 41, and 37.16%, respectively, relative to the control soil (layer 0–50 cm, Figure 3c and Table 2). The highest values of total P in depth were seen in the soil with PD applications, with increment until the 40 cm layer. On average, applications of PD, PS, CS, and MF increased the amounts of total P by 91, 48, 42, and 18 kg P ha−1 year−1 in the soil up to the 50 cm layer, relative to the control.

3.2. P Mass Balance

Cumulative P inputs were 746, 938, 1431, and 538 kg P ha−1 in soils with CS, PS, PD, and MF applications, respectively (Table 2). On average, applications of CS, PS, PD, and MF corresponded to inputs of 62, 78, 119, and 45 kg P ha−1 year−1, respectively. The variation of P inputs among the crops (Figure 1) was correlated with the dry matter concentrations in the manure (Pearson’s correlation coefficient = 0.41, p value < 0.001) and also by the N content present in the manure, since the amounts of animal manure applied were established considering the N concentration.
Average grain and cereal yields (corn, beans, and wheat) over 12 years increased significantly in soils with PD, PS, CS, and FM applications, relative to the control (Figure 4a). On average, corn yield (10 crops) grown on soils with FM, PS, CS, and PD applications were 1.8, 2.2, 2.4, 2.0, and 2.5 times higher, respectively, relative to yield on the control soil. The yield of beans grown on soils with MF, PS, CS, and PD applications increased by 1.5, 2.7, 3.0, and 4.1, times respectively, relative to the control. The highest wheat yield was observed in soils grown with MF applications followed by PS, PD, and CS with increments of 2.6, 2.3, 2.3, and 2.2 times, respectively, relative to the control. More information on the grain yield responses of the crops in this study can be seen in Bacca et al., Marchezan et al., and Ferreira et al. [13]. The highest cumulative P exports (sum of corn, bean, and wheat exports) over 12 years were observed in crops grown on soils subjected to manure and mineral phosphate fertilizer applications (Figure 4b). Exports of P by crops grown on fertilized soils averaged 150, 143, 119, and 85 kg P ha−1, higher than that observed in the control soil (Figure 4b and Table 2).
The applications of MF, CS, PS, and PD promoted the accumulation of 18, 42, 48, and 100 kg P ha−1 year−1 (Table 2). The mass balance calculations showed that between 41 and 72% of the total P added by manure and MF was accumulated in the 50 cm layer. Furthermore, most of the P added in the cropping system can be explained by soil storage, exports by grains and crop residues (98, 89, 89, and 77%, for treatments CS, PS, PD, and MF, respectively). The P balance efficiency was highest in the soil with FM, at 40%, decreasing to 30, 29, and 19% in soils with CS, PS, and PD applications, respectively.

4. Discussion

4.1. Soil P Balance

Additions of CS, PS, PD, and MF in a no-till system for several years increased the amounts of P in the soil, with its vertical movement proportional to the amount of P applied by each source [31,32,33]. The difference in the increment in the amounts of total P compared to the control was observed in all treatments up to 40 cm depth, with a higher proportion in the upper soil layer (0–20 cm). This happened because, in the no-till farming system, the soil is not disturbed, and the fertilizer sources are applied on the surface. In addition, the rapid adsorption of the phosphate ion by the functional groups of reactive soil particles decreases the availability of P in solution, preventing it from moving to depth [34,35].
The fraction with which P was accumulated in the soil differed among treatments and in depth, with preferential accumulation of inorganic P in PD. In the soil layers with PS, CS, and MF applications, the accumulated amounts of organic and inorganic P were similar. The accumulation of P in the inorganic fraction in PD may have happened, especially, because of the higher concentration of P per kilogram of PD (Figure 1). This is because PD waste has animal feces and urine and food leftovers, which are deposited for several months in pig breeding stalls [36]. Therefore, PD has a high concentration of nutrients, among them, P. Moreover, the concentration of inorganic P in swine manure is higher in relation to the manure from other livestock farms, due to the absence of the enzyme phytase in the digestive tract, which prevents the degradation and utilization of compounds that contain P, excreting in the feces the vast majority of the P consumed via grains [37,38].
The higher amounts of organic P observed in the topsoil layers with CS, PS, PD, and MF applications can be explained by the higher biomass production of the soil cover plant species and grain crops [32]. Plants during their growth and development absorb P only in the orthophosphate (Pi) form and at the end of the crop cycle the P, after being incorporated into the tissues, returns to the soil in the form of organic P [32,39,40]. The replacement of P to the soil happens because of deposition and decomposition of the aerial part of the crops [35], as well as by root senescence [41].
The successive Inputs of P by CS, PS, and PD significantly increased the labile P fractions, especially the concentrations of weakly bound P (P-NH4 Cl) in the topsoil (0–4 and 4–10 cm), with the highest values observed in the soil with PD, because of the greater amount of P added (Figure 1 and Figure 5). The increase in P-NH4 Cl concentrations in the topsoil layers may potentiate possible transfers of soluble P or P adsorbed to functional groups of solid particles in the surface runoff solution [42,43]. Because of the high affinity of the phosphate ion to reactive soil particles and the small movement of P in the profile observed in this study, we believe that it is unlikely that P was transferred by leaching [23,44,45].
In a study conducted in southern Brazil, Ceretta et al. [44] observed losses ranging from 6.4 and 14.4 kg P ha−1 by surface runoff, in an experiment conducted on sandy soil (53% sand and 4% slope), subjected to applications of doses from 20 to 80 m3 ha−1 of liquid pig manure for 8 years. In the same experiment, Tiecher et al. [45] quantified during 6 crops (2.5 years) the forms of P transferred by leaching. The authors concluded that less than 0.2 kg P ha−1 year−1 applied at a dose of 80 m3 ha−1 of liquid pig manure was leached. When compared to P losses by surface runoff, leaching losses are economically and agronomically negligible. However, due to the numerous precipitations and high volume of rainfall in subtropical climate, which may increase the water flux in the profile, the leached P concentrations may exceed the limits recommended by the environmental legislation, especially in soils with high P contents. Thus, although the amounts of P transferred by leaching are small, these losses can have a significant environmental impact in groundwater contamination, especially when the water table is shallow [45,46,47,48].
Soil P accumulation in the inorganic fractions and P movement in the soil profile are similar to other studies that evaluated the effects of long-term P fertilization by phosphate fertilizers and different animal manure sources [28,49,50,51,52]. Studies conducted on soils and climatic conditions similar to those of the present study (<20% clay), with a history of swine manure application, also observed a movement of labile P greater than 30 cm [29,53,54]. While in more clayey soils (>40% clay), P movement observed in experiments with long application history was less than 20 cm [28,55]. The lower migration of P to deeper layers in clayey soils compared to sandy soils is directly related to the higher capacity of P adsorption in soils, especially in free functional groups on the surfaces of clays and Fe and Al oxides and hydroxides [25,56].
Most of the P added by CS, PS, PD, and MF was accounted for by mass balance calculations. However, it was not possible to calculate 0.8, 8.0, 12, and 10 kg P ha−1 year−1 in soils that received applications of CS, PS, PD, and MF, respectively (these values correspond respectively to 1, 13, 10, and 22% of the P added annually). The highest amount of unaccounted P in the PD soil (144 kg P ha−1) may be related to the highest concentration of weakly bound P among the treatments evaluated, which indicate potential P losses (Table 1 and Figure 5).
Accumulated P in the soil represented the main P compartment in this agricultural system, even in a sandy soil (10.8% clay), which has low P adsorption capacity [53,55]. In the soils subjected to the additions of CS, PS, and PD, this may have happened because of the constant additions of P many times higher than the demands of the crops. Added to this, the maintenance of cultural residues on the surface and the cultivation of cover crops over the years may have collaborated in the cycling and maintenance of P in this compartment. As discussed above, the transfer of P by leaching may be insignificant in this soil and, therefore, 2.5 to 23% of the P inputs, which could not be accounted for by the mass balance, may represent P lost by surface runoff.

4.2. Environmental and Agronomic Implications

The applications, especially of liquid (CS and PS) and solid (PD) animal manure, significantly increased the P content of the soil compared to the control, causing concern about the continued use of these nutrient sources. Over 12 years, the amounts of animal manure were defined considering the efficiency index of each manure, the concentration of N, and the need of each crop, as established by the regional recommendation on the use of animal manure [27]. Probably, these variables are efficient in defining the amount of manure in the first years of applications. However, when the manure is applied repeatedly over several years, it can greatly increase the P contents in the soil, exceeding the extraction capacity of the plants, causing significant accumulation of this nutrient in the soil [28,32,49]. A study by Gatiboni et al. [25] established the environmental safety limits for P in soil in the southern region of Brazil. The method is based on the concentrations of bioavailable P in soil (Mehlich-1 extractor) in the 0–10 cm layer, soil P adsorption capacity, and clay content, according to the equation: P threshold (mg P kg−1) = 20 + clay content (%). When employing this prediction for the soil of the present study (clay content of 10.8%), the P threshold (30.8 mg P kg−1) was exceeded by 79, 135, 148, and 202% for MF, PS, CS, and PD applications, respectively. These values correspond to the significant increment in the concentrations of weakly bound P (NH4 Cl-P Figure 5) in the topsoil layer (0–4), representing potential risk of P losses by surface runoff [23,26]. In the absence of measures to mitigate excess soil P, ongoing applications could escalate environmental hazards, compromising surface water quality [24,56,57].
The grain yield data obtained for corn, beans, and wheat over the years showed the increase in productivity achieved through applications of CS, PS, and PD. The improvement in soil nutrient availability demonstrates the potential of using animal manure as a source of fertilizer [13,58]. However, ongoing monitoring of critical soil P levels is essential. This study’s P mass balance results indicate an imbalance in P additions among treatments, resulting in significant P accumulation in the soil compartment and low P efficiency for plant uptake [48,59,60]. In addition, the differences in the amounts of P added by the manure over the seasons may have been a consequence of the composition of the manure, especially the effect of dry matter concentration per kilogram of manure added. These factors contributed to the excessive addition of P in the PD treatment, thus decreasing the overall P mass balance efficiency of the culture system to unsatisfactory levels. These results raise concerns not only about the amount of animal waste that can be applied per area, but also about the chemical composition of the added manure in agricultural areas.
Continued additions of P and the accumulation of P legacy in soils after long periods of animal manure applications in the same areas can present a serious environmental risk, as was widely demonstrated in this long-term study. Similarly, applications of mineral phosphate fertilizers (MF) also implied a worrisome accumulation of P in the soil, especially in the topsoil. However, the amounts of plant-available P (Mehlich-1 extractor) is within the adequate range (30–50 mg P kg−1), requiring annual maintenance fertilization. We observed that the addition of MF increased the fractions of P not available to plants, with 51% of the total P in the stable and residual form in the surface layer (Figure 2d). This may have happened because of the high affinity of the phosphate ion that is rapidly adsorbed to the surface functional groups of reactive soil particles [34]. This effect was more pronounced in soil with MF applications, probably because of the lower pH value, caused by acidification, caused by the release of H+ coming from the solubilization of N fertilizers [61,62]. The decrease in pH causes the protonation of functional groups of reactive soil particles, thus increasing the binding energy of P with the soil solid phase, increasing the more stable P fractions [63]. Low efficiency and accumulation of P in the surface layers in soils with mineral phosphate fertilizer additions is very common in most highly weathered soils.
Considering the obtained results, research is needed to find strategies to mobilize recalcitrant P fractions present in the soil and decrease the dependence on constant phosphate fertilizations [40,64,65].

5. Conclusions

The additions of PD, PS, CS, and MF to the soil over 12 years caused the accumulation of 100, 48, 42, and 18 kg P ha−1 year−1, respectively.
All P sources caused significant P accumulation in the surface layers, especially in the inorganic P fraction. However, in the soil with PD, the P moved up to 40 cm and with applications of CS, PS, and MF, the movement was observed up to 30 cm.
Soils with animal manure applications, especially PD, have a worrisome potential for P transfer. Thus, threshold levels of P should be constantly monitored, and the fertilizer recommendation system should consider the history of nutrient source applications in the cropping system to establish the need and dose for a new application of a P-containing source.
The P mass balance revealed that 77 and 98% of P input by animal manure and MF was accounted for in grain exports (17–34%), soil storage (41–72%), and post-harvest residues (<1%), with the remaining 2–33% unaccounted for because of transfers by surface runoff. This reinforces the need for conservation farming practices to avoid runoff and soil losses.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/agriculture13091762/s1, Table S1. Distribution of P pools in the soil profile after 18 successive applications (2004–2016) of cattle slurry (CS), pig slurry (PS), pig deep-litter (PD), and mineral fertilizer (MF). Data presented are means (n = 4); Table S2. Quantities of inorganic, organic, and total P (kg P ha-1) accumulated in the soil profile after 18 successive applications (2004–2016) of cattle slurry (CS), pig slurry (PS), pig deep-litter (PD), and mineral fertilizer (MF). Data presented are means (n = 4); Table S3. P contents in the different forms of Hedley fractionation and sum of the P contents extracted from the soil profile after 18 successive applications (2004–2016) of cattle slurry (CS), pig slurry (PS), pig deep-litter (PD), and mineral fertilizer (MF). Data presented are means (n = 4).

Author Contributions

Conceptualization, C.M. and P.A.A.F.; methodology, R.F.V., G.S., N.M.P. and A.L.T.; investigation, C.M., G.B. (Gustavo Boitt) and C.R.L.; writing—review and editing, C.M., G.B. (Gustavo Boitt) and G.B. (Gustavo Brunetto); visualization, C.M.; supervision, G.B. (Gustavo Brunetto); project administration, C.A.C., P.A.A.F. and G.B. (Gustavo Brunetto). All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Foundation for Research Support of the State of Rio Grande do Sul (FAPERGS) (process number 19/2551-0001937-8), Agrisus foundation (PA 9491/18 and PA 2909/20), and the Brazilian National Council for Scientific and Technological Development (CNPq), process number 400982/2016-1.

Institutional Review Board Statement

The data presented in this study is available upon request to the corresponding author. The data is not publicly available as it is part of a larger dataset for which the results have not yet been published, but will be made available soon.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Scholz, R.W.; Ulrich, A.E.; Eilittä, M.; Roy, A. Sustainable Use of Phosphorus: A Finite Resource. Sci. Total Environ. 2013, 461–462, 799–803. [Google Scholar] [CrossRef] [PubMed]
  2. Scholz, R.W.; Wellmer, F.W. Approaching a Dynamic View on the Availability of Mineral Resources: What We May Learn from the Case of Phosphorus? Glob. Environ. Change 2013, 23, 11–27. [Google Scholar] [CrossRef]
  3. Cordell, D.; Drangert, J.-O.; White, S. The Story of Phosphorus: Global Food Security and Food for Thought. Glob. Environ. Change 2009, 19, 292–305. [Google Scholar] [CrossRef]
  4. Geissler, B.; Mew, M.C.; Steiner, G. Phosphate Supply Security for Importing Countries: Developments and the Current Situation. Sci. Total Environ. 2019, 677, 511–523. [Google Scholar] [CrossRef] [PubMed]
  5. Cordell, D.; Rosemarin, A.; Schröder, J.J.; Smit, A.L. Towards Global Phosphorus Security: A Systems Framework for Phosphorus Recovery and Reuse Options. Chemosphere 2011, 84, 747–758. [Google Scholar] [CrossRef]
  6. Schneider, K.D.; Thiessen Martens, J.R.; Zvomuya, F.; Reid, D.K.; Fraser, T.D.; Lynch, D.H.; O’Halloran, I.P.; Wilson, H.F. Options for Improved Phosphorus Cycling and Use in Agriculture at the Field and Regional Scales. J. Environ. Qual. 2019, 48, 1247–1264. [Google Scholar] [CrossRef] [PubMed]
  7. ABPA Relatório Anual; Associação Brasileira de Proteína Animal: São Paulo, Brazil, 2020; Volume 160.
  8. IBGE Pesquisa Trimestral Do Leite—Primeiros Resultados Dados Da Quantidade de Leite Cru Processado Pela Indústria; IBGE: Rio de Janeiro, Brazil, 2019.
  9. Haygarth, P.M.; Jarvie, H.P.; Powers, S.M.; Sharpley, A.N.; Elser, J.J.; Shen, J.; Peterson, H.M.; Chan, N.-I.; Howden, N.J.K.; Burt, T.; et al. Sustainable Phosphorus Management and the Need for a Long-Term Perspective: The Legacy Hypothesis. Environ. Sci. Technol. 2014, 48, 8417–8419. [Google Scholar] [CrossRef]
  10. Embrapa Sistema Brasileiro de Classificação de Solos; 5a.; Embrapa: Rio de Janeiro, Brazil, 2018; ISBN 978-85-7035-800-4.
  11. Soil Survey Staff Keys to Soil Taxonomy. Soil. Conserv. Serv. 2014, 12, 410. [CrossRef]
  12. Lourenzi, C.R.; Ceretta, C.A.; Ciancio, N.H.R.; Tiecher, T.L.; da Silva, L.O.S.; De Conti, L.; Girotto, E.; Ferreira, P.A.A.; Vidal, R.F.; Scopel, G.; et al. Forms of Nitrogen and Phosphorus Transfer by Runoff in Soil under No-Tillage with Successive Organic Waste and Mineral Fertilizers Applications. Agric. Water Manag. 2021, 248, 106779. [Google Scholar] [CrossRef]
  13. Ferreira, P.A.A.; Ceretta, C.A.; Lourenzi, C.R.; De Conti, L.; Marchezan, C.; Girotto, E.; Tiecher, T.L.; Palermo, N.M.; Parent, L.É.; Brunetto, G. Long-Term Effects of Animal Manures on Nutrient Recovery and Soil Quality in Acid Typic Hapludalf under No-Till Conditions. Agronomy 2022, 12, 243. [Google Scholar] [CrossRef]
  14. Marchezan, C.; Ferreira, P.A.A.; Silva, L.S.; Bacca, A.; Krug, A.V.; Nicoloso, F.T.; Tarouco, C.P.; Tiecher, T.L.; Brunetto, G.; Ceretta, C.A. Nitrogen Availability and Physiological Response of Corn After 12 Years with Organic and Mineral Fertilization. J. Soil. Sci. Plant Nutr. 2020, 20, 979–989. [Google Scholar] [CrossRef]
  15. Manual de Métodos de Análise de Solo; Embrapa: Brasília, Brazil, 1997; Volume 2, ISBN 8585864036.
  16. Hedley, M.J.; Stewart, J.W.B.; Chauhan, B.S. Changes in Inorganic and Organic Soil Phosphorus Fractions Induced by Cultivation Practices and by Laboratory Incubations1. Soil Sci. Soc. Am. J. 1982, 46, 970. [Google Scholar] [CrossRef]
  17. Condron, L.M.; Davis, M.R.; Newman, R.H.; Cornforth, I.S. Influence of Conifers on the Forms of Phosphorus in Selected New Zealand Grassland Soils. Biol. Fertil. Soils 1996, 21, 37–42. [Google Scholar] [CrossRef]
  18. Condron, L.M.; Newman, S. Revisiting the Fundamentals of Phosphorus Fractionation of Sediments and Soils. J. Soils Sediments 2011, 11, 830–840. [Google Scholar] [CrossRef]
  19. Olsen, S.R.; Sommers, L.E. Phosphorus. Methods of Soil Analysis, Part 2. In Chemical and Microbiological Properties; Page, A.L., Miller, R.H., Keeney, D.R., Eds.; Soil Science Society of America: Madison, WI, USA, 1982; pp. 403–430. [Google Scholar]
  20. Dick, W.A.; Tabatabai, M.A. Determination of Orthophosphate in Aqueous Solutions Containing Labile Organic and Inorganic Phosphorus Compounds 1. J. Environ. Qual. 1977, 6, 82. [Google Scholar] [CrossRef]
  21. Murphy, J.; Riley, J.P. A Modified Single Solution Method for the Determination of Phosphate in Natural Waters. Anal. Chim. Acta 1962, 27, 31–36. [Google Scholar] [CrossRef]
  22. Cross, A.F.; Schlesinger, W.H. A Literature Review and Evaluation of the. Hedley Fractionation: Applications to the Biogeochemical Cycle of Soil Phosphorus in Natural Ecosystems. Geoderma 1995, 64, 197–214. [Google Scholar] [CrossRef]
  23. McDowell, R.W.; Sharpley, A.N. Approximating Phosphorus Release from Soils to Surface Runoff and Subsurface Drainage. J. Environ. Qual. 2001, 30, 508. [Google Scholar] [CrossRef] [PubMed]
  24. Sharpley, A.N.; Mcdowell, R.W.; Kleinman, P.J.A. Phosphorus Loss from Land to Water: Integrating Agricultural and Environmental Management. Proc. Plant Soil 2001, 237, 287–307. [Google Scholar] [CrossRef]
  25. Gatiboni, L.C.; Smyth, T.J.; Schmitt, D.E.; Cassol, P.C.; de Oliveira, C.M.B. de Soil Phosphorus Thresholds in Evaluating Risk of Environmental Transfer to Surface Waters in Santa Catarina, Brazil. Rev. Bras. Cienc. Solo 2015, 39, 1225–1234. [Google Scholar] [CrossRef]
  26. Gatiboni, L.C.; Nicoloso, R.d.S.; Mumbach, G.L.; de Souza Junior, A.A.; Dall’orsoletta, D.J.; Schmitt, D.E.; Smyth, T.J. Establishing Environmental Soil Phosphorus Thresholds to Decrease the Risk of Losses to Water in Soils from Rio Grande Do Sul, Brazil. Rev. Bras. Cienc. Solo 2020, 44, 1–14. [Google Scholar] [CrossRef]
  27. Manual of Liming and Fertilization for the States of Rio Grande do Sul and Santa Catarina; Brazilian Soil Science Society: Viçosa, Brazil, 2016; 376p, ISBN 978-85-66301-80-9.
  28. Boitt, G.; Schmitt, D.E.; Gatiboni, L.C.; Wakelin, S.A.; Black, A.; Sacomori, W.; Cassol, P.C.; Condron, L.M. Fate of Phosphorus Applied to Soil in Pig Slurry under Cropping in Southern Brazil. Geoderma 2018, 321, 164–172. [Google Scholar] [CrossRef]
  29. Lourenzi, C.R.; Ceretta, C.A.; Cerini, J.B.; Ferreira, P.A.A.; Lorensini, F.; Girotto, E.; Tiecher, T.L.; Schapanski, D.E.; Brunetto, G. Available Content, Surface Runoff and Leaching of Phosphorus Forms in a Typic Hapludalf Treated with Organic and Mineral Nutrient Sources. Rev. Bras. Cienc. Solo 2014, 38, 544–556. [Google Scholar] [CrossRef]
  30. RStudio Team. 2015. RStudio: Integrated Development for R. RStudio, PBC, Boston, MA. Available online: https://www.rstudio.com/ (accessed on 27 August 2023).
  31. Filho, J.d.S.O.; Ferrari, A.C.; Pereira, M.G.; Pinto, L.A.d.S.R.; Pinto, L.A.d.S.R.; Matos, T.S. Phosphorus Accumulation in Soil after Successive Applications of Swine Manure: A Long-Term Study in Brazil. Environ. Earth Sci. 2020, 79, 62. [Google Scholar] [CrossRef]
  32. Oyeleke Azeez, M.; Therese Christensen, J.; Ravnskov, S.; Johann Heckrath, G.; Labouriau, R.; Christensen, B.T.; Holton Rubaek, G. Phosphorus in an Arable Coarse Sandy Soil Profile after 74 Years with Different Lime and P Fertilizer Applications. Geoderma 2020, 376, 114555. [Google Scholar] [CrossRef]
  33. Parfitt, R.L. Phosphate Reactions with Natural Allophane, Ferrihydrite and Goethite. J. Soil. Sci. 1989, 40, 359–369. [Google Scholar] [CrossRef]
  34. Barrow, N.J. Landmark Papers A Mechanistic Model for Describing the Sorption and Desorption of Phosphate by Soil. Eur. J. Soil. Sci. 2015, 66, 9–18. [Google Scholar] [CrossRef]
  35. Oliveira, P.; Soares, W. Produção e Manejo de Dejetos de Suínos. Produção Anim. Visão Bras. 2001, d, 72–90. [Google Scholar]
  36. Baxter, C.A.; Joern, B.C.; Ragland, D.; Sands, J.S.; Adeola, O. Phytase, High-Available-Phosphorus Corn, and Storage Effects on Phosphorus Levels in Pig Excreta. J. Environ. Qual. 2003, 32, 1481. [Google Scholar] [CrossRef] [PubMed]
  37. Wienhold, B.J.; Miller, P.S. Phosphorus Fractionation in Manure from Swine Fed Traditional and Low-Phytate Corn Diets. J. Environ. Qual. 2004, 33, 389–393. [Google Scholar] [CrossRef]
  38. Pagliari, P.H.; Laboski, C.A.M. Investigation of the Inorganic and Organic Phosphorus Forms in Animal Manure. J. Environ. Qual. 2012, 41, 901–910. [Google Scholar] [CrossRef]
  39. Hallama, M.; Pekrun, C.; Lambers, H.; Kandeler, E. Hidden Miners—The Roles of Cover Crops and Soil Microorganisms in Phosphorus Cycling through Agroecosystems. Plant Soil. 2019, 434, 7–45. [Google Scholar] [CrossRef]
  40. de Oliveira, R.A.; José Comin, J.; Tiecher, T.; Piccin, R.; Michelon Somavilla, L.; Loss, A.; Rogério Lourenzi, C.; Brunetto, G. Release of Phosphorus Forms from Cover Crop Residues in Agroecological No-Till Onion Production. Artic. Rev. Bras. Cienc. Solo 2017, 41, 160272. [Google Scholar] [CrossRef]
  41. Pausch, J.; Kuzyakov, Y. Carbon input by roots into the soil: Quantification of rhizodeposition from root to ecosystem scale. Glob. Chang. Biol. 2018, 24, 1–12. [Google Scholar] [CrossRef]
  42. Pierzynski, G.M.; McDowell, R.W.; Sims, J.T. Chemistry, Cycling, and Potatial Movemet of Inorganic Phosphorus in Soil. Phosphorus Agric. Environ. 2005, 46, 53–86. [Google Scholar]
  43. Mcdowell, R.; Mcdowell, R.W.; Sharpley, A.N. Soil Phosphorus Fractions in Solution: Influence of Fertiliser and Manure, Filtration and Method of Determination. Chemosphere 2014, 45, 737–748. [Google Scholar] [CrossRef]
  44. Ceretta, C.A.; Girotto, E.; Lourenzi, C.R.; Trentin, G.; Vieira, R.C.B.; Brunetto, G. Nutrient Transfer by Runoff under No Tillage in a Soil Treated with Successive Applications of Pig Slurry. Agric. Ecosyst. Environ. 2010, 139, 689–699. [Google Scholar] [CrossRef]
  45. Tiecher, T.L.; Lourenzi, C.R.; Girotto, E.; Tiecher, T.; De Conti, L.; Marques, A.C.R.; Silva, L.O.S.; Marchezan, C.; Brunetto, G.; Ceretta, C.A. Phosphorus Forms Leached in a Sandy Typic Hapludalf Soil under No-Tillage with Successive Pig Slurry Applications. Agric. Water Manag. 2020, 242, 106406. [Google Scholar] [CrossRef]
  46. Girotto, E.; Ceretta, C.A.; Lourenzi, C.R.; Lorensini, F.; Tiecher, T.L.; Vieira, R.C.B.; Trentin, G.; Basso, C.J.; Miotto, A.; Brunetto, G. Nutrient Transfers by Leaching in a No-Tillage System through Soil Treated with Repeated Pig Slurry Applications. Nutr. Cycl. Agroecosyst. 2013, 95, 115–131. [Google Scholar] [CrossRef]
  47. Toor, G.S.; Sims, J.T. Phosphorus Leaching in Soils Amended with Animal Manures Generated from Modified Diets. J. Environ. Qual. 2016, 45, 1385–1391. [Google Scholar] [CrossRef] [PubMed]
  48. Khan, A.; Lu, G.; Ayaz, M.; Zhang, H.; Wang, R.; Lv, F.; Yang, X.; Sun, B.; Zhang, S. Phosphorus Efficiency, Soil Phosphorus Dynamics and Critical Phosphorus Level under Long-Term Fertilization for Single and Double Cropping Systems. Agric. Ecosyst. Environ. 2018, 256, 1–11. [Google Scholar] [CrossRef]
  49. Guardini, R.; Comin, J.J.; Schmitt, D.E.; Tiecher, T.; Bender, M.A.; dos Santos, D.R.; Mezzari, C.P.; Oliveira, B.S.; Gatiboni, L.C.; Brunetto, G. Accumulation of Phosphorus Fractions in Typic Hapludalf Soil after Long-Term Application of Pig Slurry and Deep Pig Litter in a No-Tillage System. Nutr. Cycl. Agroecosyst. 2012, 93, 215–225. [Google Scholar] [CrossRef]
  50. Vanden Nest, T.; Vandecasteele, B.; Ruysschaert, G.; Cougnon, M.; Merckx, R.; Reheul, D. Effect of Organic and Mineral Fertilizers on Soil P and C Levels, Crop Yield and P Leaching in a Long Term Trial on a Silt Loam Soil. Agric. Ecosyst. Environ. 2014, 197, 309–317. [Google Scholar] [CrossRef]
  51. Requejo, M.I.; Eichler-Löbermann, B. Organic and Inorganic Phosphorus Forms in Soil as Affected by Long-Term Application of Organic Amendments. Nutr. Cycl. Agroecosyst. 2014, 100, 245–255. [Google Scholar] [CrossRef]
  52. Tian, J.; Boitt, G.; Black, A.; Wakelin, S.; Condron, L.M.; Chen, L. Accumulation and Distribution of Phosphorus in the Soil Profile under Fertilized Grazed Pasture. Agric. Ecosyst. Environ. 2017, 239, 228–235. [Google Scholar] [CrossRef]
  53. De Conti, L.; Ceretta, C.A.; Ferreira, P.A.A.; Lorensini, F.; Lourenzi, C.R.; Vidal, R.F.; Tassinari, A.; Brunetto, G. Effects of Pig Slurry Application and Crops on Phosphorus Content in Soil and the Chemical Species in Solution. Rev. Bras. Ciência Solo 2015, 39, 774–787. [Google Scholar] [CrossRef]
  54. Cassol, P.C.; da Costa, A.C.; Ciprandi, O.; Pandolfo, C.M.; Ernani, P.R. Disponibilidade de Macronutrientes e Rendimento de Milho Em Latossolo Fertilizado Com Dejeto Suíno. Rev. Bras. Cienc. Solo 2012, 36, 1911–1923. [Google Scholar] [CrossRef]
  55. Rheinheimer, D.S.; Anghinoni, I.; Conte, E.; Kaminski, J.; Gatiboni, L.C. Dessorção de Fósforo Avaliada Por Extrações Sucessivas Em Amostras de Solo Provenientes Dos Sistemas Plantio Direto e Convencional. Ciência Rural. 2003, 33, 1053–1059. [Google Scholar] [CrossRef]
  56. Chen, M.; Graedel, T.E. A Half-Century of Global Phosphorus Flows, Stocks, Production, Consumption, Recycling, and Environmental Impacts. Glob. Environ. Change 2016, 36, 139–152. [Google Scholar] [CrossRef]
  57. Penn, C.; Bowen, J.M. Design and Construction of Phosphorus Removal Structures for Improving Water Quality. In Design and Construction of Phosphorus Removal Structures for Improving Water Quality; Springer: Berlin/Heidelberg, Germany, 2017; pp. 1–228. [Google Scholar] [CrossRef]
  58. Ferreira, P.A.A.; Coronas, M.V.; Dantas, M.K.L.; Somavilla, A.; Brunetto, G.; Ceretta, C.A.; Giacomini, S.J.; Gubiani, P.I.; Boitt, G.; Soares, C.R.F.S.; et al. Repeated Manure Application for Eleven Years Stimulates Enzymatic Activities and Improves Soil Attributes in a Typic Hapludalf. Agronomy 2021, 11, 2467. [Google Scholar] [CrossRef]
  59. Sharpley, A.; Jarvie, H.P.; Buda, A.; May, L.; Spears, B.; Kleinman, P. Phosphorus Legacy: Overcoming the Effects of Past Management Practices to Mitigate Future Water Quality Impairment. J. Environ. Qual. 2013, 42, 1308–1326. [Google Scholar] [CrossRef] [PubMed]
  60. Wironen, M.B.; Bennett, E.M.; Erickson, J.D. Phosphorus Flows and Legacy Accumulation in an Animal-Dominated Agricultural Region from 1925 to 2012. Glob. Environ. Change 2018, 50, 88–99. [Google Scholar] [CrossRef]
  61. Peryea, F.J.; Burrows, R.L. Soil Acidification Caused by Four Commercial Nitrogen Fertilizer Solutions and Subsequent Soil PH Rebound. Commun. Soil. Sci. Plant Anal. 1999, 30, 525–533. [Google Scholar] [CrossRef]
  62. Schroder, J.L.; Zhang, H.; Girma, K.; Raun, W.R.; Penn, C.J.; Payton, M.E. Soil Acidification from Long-Term Use of Nitrogen Fertilizers on Winter Wheat. Soil. Sci. Soc. Am. J. 2011, 75, 957–964. [Google Scholar] [CrossRef]
  63. Barrow, N.J. The Effects of PH on Phosphate Uptake from the Soil. Plant Soil. 2017, 410, 401–410. [Google Scholar] [CrossRef]
  64. Haygarth, P.M.; Bardgett, R.D.; Condron, L.M. Nitrogen and Phosphorus Cycles and Their Management; Gregory, P.J., Nortcliff, S., Eds.; Soil Conditions and Plant Growth; Blackwell Publishing Ltd.: Oxford, UK, 2013. [Google Scholar]
  65. Sattari, S.Z.; Bouwman, A.F.; Giller, K.E.; van Ittersum, M.K. Residual Soil Phosphorus as the Missing Piece in the Global Phosphorus Crisis Puzzle. Proc. Natl. Acad. Sci. USA 2012, 109, 6348–6353. [Google Scholar] [CrossRef]
Figure 1. Total P applied (expressed in kg P ha−1) to field plots receiving cattle slurry (CS), pig slurry (PS), pig deep-litter (PD), and mineral fertilizer (MF) from 2004 to 2012 under no-tillage system.
Figure 1. Total P applied (expressed in kg P ha−1) to field plots receiving cattle slurry (CS), pig slurry (PS), pig deep-litter (PD), and mineral fertilizer (MF) from 2004 to 2012 under no-tillage system.
Agriculture 13 01762 g001
Figure 2. Average P quantities (kg P ha−1) determined in P fractions for soils sampled from different depths after 12 years under applications cattle slurry (CS), pig slurry (PS), pig deep-litter (PD), and mineral fertilizer (MF): labile P (PiAM + Pibic + Pobic) (a), moderately labile P (PiOH_I + PoOH_I) (b), stable P (PiHCl + PiOH_II + PoOH_II) (c), and residual P (Pres) (d). * p < 0.05; ** p < 0.01 and ns: not significant (Scott–Knott test at 5%). Standard deviation is presented on the surface of the bars. Data presented are means (n = 4).
Figure 2. Average P quantities (kg P ha−1) determined in P fractions for soils sampled from different depths after 12 years under applications cattle slurry (CS), pig slurry (PS), pig deep-litter (PD), and mineral fertilizer (MF): labile P (PiAM + Pibic + Pobic) (a), moderately labile P (PiOH_I + PoOH_I) (b), stable P (PiHCl + PiOH_II + PoOH_II) (c), and residual P (Pres) (d). * p < 0.05; ** p < 0.01 and ns: not significant (Scott–Knott test at 5%). Standard deviation is presented on the surface of the bars. Data presented are means (n = 4).
Agriculture 13 01762 g002
Figure 3. Average quantities (kg P ha−1) of total extracted inorganic P (Pi), total extracted organic P (Po), and total P determined for soils sampled from different depths after 12 years under applications cattle slurry (CS), pig slurry (PS), pig deep-litter (PD), and mineral fertilizer (MF) (total extracted inorganic P = sum of PiAM, Pibic, PiOH_I, PiHCl, PiOH_II (a), total extracted organic P = sum of Pobic + PoOH_I + PoOH_II) (b) and Total P (c). ** p < 0.01 and ns: not significant (Scott–Knott test at 5%). Standard deviation is presented on the surface of the bars. Data presented are means (n = 4).
Figure 3. Average quantities (kg P ha−1) of total extracted inorganic P (Pi), total extracted organic P (Po), and total P determined for soils sampled from different depths after 12 years under applications cattle slurry (CS), pig slurry (PS), pig deep-litter (PD), and mineral fertilizer (MF) (total extracted inorganic P = sum of PiAM, Pibic, PiOH_I, PiHCl, PiOH_II (a), total extracted organic P = sum of Pobic + PoOH_I + PoOH_II) (b) and Total P (c). ** p < 0.01 and ns: not significant (Scott–Knott test at 5%). Standard deviation is presented on the surface of the bars. Data presented are means (n = 4).
Agriculture 13 01762 g003
Figure 4. Response of the production of grains of corn, beans, and wheat (a) and dry mass maize, beans, wheat, oats, and turnips (b), depending on the successive applications cattle slurry (CS), pig slurry (PS), pig deep-litter (PD), and mineral fertilizer (MF). The results represent an average of 12 years (2004–20016) of crop rotation management with 18 treatment applications. Corn was cultivated in 10 seasons, while beans and wheat were grown in the 2012/2013 and 2014 cropping seasons, respectively. For cover crops, oats were grown for 8 cropping seasons and O. Radish in three cropping seasons (2005, 2007, and 2009). Different letters between treatments at the same depth indicate significant difference (Scott–Knott test at 5%). Standard deviation is presented on the surface of the bars.
Figure 4. Response of the production of grains of corn, beans, and wheat (a) and dry mass maize, beans, wheat, oats, and turnips (b), depending on the successive applications cattle slurry (CS), pig slurry (PS), pig deep-litter (PD), and mineral fertilizer (MF). The results represent an average of 12 years (2004–20016) of crop rotation management with 18 treatment applications. Corn was cultivated in 10 seasons, while beans and wheat were grown in the 2012/2013 and 2014 cropping seasons, respectively. For cover crops, oats were grown for 8 cropping seasons and O. Radish in three cropping seasons (2005, 2007, and 2009). Different letters between treatments at the same depth indicate significant difference (Scott–Knott test at 5%). Standard deviation is presented on the surface of the bars.
Agriculture 13 01762 g004
Figure 5. Representation of weakly bound P extracted with 1 mol L−1 NH4Cl after long-term animal slurry and mineral fertilizer additions in the topmost soil layers (0–10 cm). Substantial increases in this labile fraction pose potential environmental risks, mainly following soil erosion by runoff. Vertical bars represent standard deviation of averages (McDowell et al. [46]). Different letters between treatments at the same depth indicate significant difference (Scott–Knott test at 5%).
Figure 5. Representation of weakly bound P extracted with 1 mol L−1 NH4Cl after long-term animal slurry and mineral fertilizer additions in the topmost soil layers (0–10 cm). Substantial increases in this labile fraction pose potential environmental risks, mainly following soil erosion by runoff. Vertical bars represent standard deviation of averages (McDowell et al. [46]). Different letters between treatments at the same depth indicate significant difference (Scott–Knott test at 5%).
Agriculture 13 01762 g005
Table 1. Chemical properties of soil collected at 0–10 cm before experiment implant (2004) and after 18 applications of organic and mineral N sources (August 2016). The values are the averages of the four plots.
Table 1. Chemical properties of soil collected at 0–10 cm before experiment implant (2004) and after 18 applications of organic and mineral N sources (August 2016). The values are the averages of the four plots.
Chemical Properties Treatments (2016)
Initial Soil (2004) *ControlCSPSPDFM
pHH2O4.95.25.34.95.74.5
SOM (g kg−1) (1)1914.021.916.920.915.9
Available P (mg kg−1) (2)20.39.776.572.993.255.1
Available K (mg kg−1) (2)6029.0117.071.8116.875.0
Exchangeable Ca (cmolc dm−3) (3)0.82.23.32.14.41.4
Exchangeable Mg (cmolc dm−3) (3)0.31.82.61.72.81.0
Exchangeable Al (cmolc dm−3) (3)0.030.50.20.60.11.3
H+Al (cmolc kg−1)3.73.73.73.31.94.1
CECeffective (cmolc kg−1)1.34.66.44.57.63.9
CECpH 7.0 (cmolc kg−1)57.89.97.29.56.7
Base Saturation (%)25.453.763.556.380.038.8
Aluminum Saturation (%)211.02.912.91.333.2
* Data adapted from Lourenzi et al. [12]. (1) Determined using modified Walkley–Black method + Van Bemmelen factor; (2) Extracted by Mehlich-1; (3) Extracted by KCl 1 mol∙L−1; Control = no fertilization; CS = cattle slurry; PS = pig slurry; PD = pig deep-litter; FM = mineral fertilizer.
Table 2. P mass balance in the no-tillage cropping system receiving long-term P inputs in cattle slurry (CS), pig slurry (PS), pig deep-litter (PD), and mineral fertilizer (MF). The P pools presented are the cumulative amounts (in units of kg P ha−1) after 18 successive applications of the treatments from 2004 to 2016. Data presented are means (n = 4).
Table 2. P mass balance in the no-tillage cropping system receiving long-term P inputs in cattle slurry (CS), pig slurry (PS), pig deep-litter (PD), and mineral fertilizer (MF). The P pools presented are the cumulative amounts (in units of kg P ha−1) after 18 successive applications of the treatments from 2004 to 2016. Data presented are means (n = 4).
P Budget in the SystemTreatments
ControlCSPSPDFM
kg P ha−1
InputsP added via treatments a0747.5938.21431.1538.71
PoolsSoil P pools b1383.9 e *1892.9 c1971.1 b2417.1 a1605.5 d
Crop P residue c4.8 d10.3 a10.3 a9.5 b8.7 c
OutputsGrain P export d98.88 d218.08 b242.28 a243.87 a184.22 c
Total P accounted e1487.6 e2121.3 c2223.7 b2670.5 a1798.4 d
P unaccounted f010.1 b98.4 a144.4 a124.2 a
Accounted P in the system (%) g098.6 a89.5 b89.9 b76.9 c
Added P accumulated in soil (%) h068.1 b62.6 c72.2 a41.1 d
P balance efficiency (%) i030.0 b29.2 b19.1 c45.5 a
* Letters identify significant differences in the means between treatments for a given P pool (Scott–Knott test at 5%). a Total P added via pig slurry after 12 years (2004–2016). b Total soil P in the layer 0–50 cm depth. c Amounts of P stored in crop residues. d Total amount of P exported in grain (corn and beans [in 2016]) in 12 years as a function of P concentration in grain (% P) and total quantities of produced grain (kg ha−1). e Sum of P accounted for in pools (soil and residues) and P outputs (P exported in grain). f Punaccounted (kg P ha−1) = Pinput − (ΔPsoil + Poutput + Presidue). g Accounted P (%) = [1 − (Punaccounted:Pinput)] × 100. h Accumulated P in soil (%) = (Δpsoil:Pinput) × 100. i P balance efficiency (%) = [Poutput:(Poutput + ΔPsoil)] × 100.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Marchezan, C.; Ferreira, P.A.A.; Boitt, G.; Palermo, N.M.; Thoma, A.L.; Vidal, R.F.; Scopel, G.; Lourenzi, C.R.; Ceretta, C.A.; Brunetto, G. Phosphorus Balance in Sandy Soil Subjected to 12 Years of Successive Applications of Animal Manure and Mineral Phosphate Fertilizer in Subtropical Climate. Agriculture 2023, 13, 1762. https://doi.org/10.3390/agriculture13091762

AMA Style

Marchezan C, Ferreira PAA, Boitt G, Palermo NM, Thoma AL, Vidal RF, Scopel G, Lourenzi CR, Ceretta CA, Brunetto G. Phosphorus Balance in Sandy Soil Subjected to 12 Years of Successive Applications of Animal Manure and Mineral Phosphate Fertilizer in Subtropical Climate. Agriculture. 2023; 13(9):1762. https://doi.org/10.3390/agriculture13091762

Chicago/Turabian Style

Marchezan, Carina, Paulo Ademar Avelar Ferreira, Gustavo Boitt, Natalia Moreira Palermo, Ana Laura Thoma, Ricardo Fagan Vidal, Gustavo Scopel, Cledimar Rogério Lourenzi, Carlos Alberto Ceretta, and Gustavo Brunetto. 2023. "Phosphorus Balance in Sandy Soil Subjected to 12 Years of Successive Applications of Animal Manure and Mineral Phosphate Fertilizer in Subtropical Climate" Agriculture 13, no. 9: 1762. https://doi.org/10.3390/agriculture13091762

APA Style

Marchezan, C., Ferreira, P. A. A., Boitt, G., Palermo, N. M., Thoma, A. L., Vidal, R. F., Scopel, G., Lourenzi, C. R., Ceretta, C. A., & Brunetto, G. (2023). Phosphorus Balance in Sandy Soil Subjected to 12 Years of Successive Applications of Animal Manure and Mineral Phosphate Fertilizer in Subtropical Climate. Agriculture, 13(9), 1762. https://doi.org/10.3390/agriculture13091762

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop