Phosphate Fertilizer in Soybean-Wheat Cropping System Under No-Till Management

Moreira, Adônis; Moraes, Larissa A. C.; Petineli, Rafael

doi:10.3389/fsufs.2020.564586

ORIGINAL RESEARCH article

Front. Sustain. Food Syst., 15 October 2020

Sec. Crop Biology and Sustainability

Volume 4 - 2020 | https://doi.org/10.3389/fsufs.2020.564586

This article is part of the Research Topic Intercropping Systems in Sustainable Agriculture View all 17 articles

Phosphate Fertilizer in Soybean-Wheat Cropping System Under No-Till Management

$\nAdnis Moreira$ Adônis Moreira¹^*

Larissa A. C. Moraes¹

Rafael Petineli²

¹Department of Plant Nutrition, Brazilian Agricultural Research Corporation (Embrapa Soja), Londrina, Brazil
²Departmant of Crop Science, Londrina State University, Londrina, Brazil

Improvements in soybean (Glycine max [L.] Merrill) and wheat (Triticum aestivum L.) yields in cropping systems under no-till management (NTM) in the subtropics have been obtained through advances in phosphorus (P) utilization, cultivar selection, and planting and harvesting strategies. This fact, along with the P application in band application in consolidated planting, has resulted in adequate P utilization efficiency and small depleted soil P levels. The aim of this study was to evaluate the efficiency of P₂O₅ rates and application methods in soybean-wheat cropping under NTM during two growing seasons (2014–2016) in a Typical Eutrorthox. Four P₂O₅ rates (0, 30, 60, and 120 kg ha⁻¹) of triple superphosphate (45% P₂O₅) were applied using two application methods (broadcast and band application at sowing) before the soybean and wheat crops were cultivated. Both P application methods resulted in a significant yield response to P₂O₅ rates for the soybean and wheat crops. Under NTM, broadcast P₂O₅ application was more effective than band application in soybean, and maximum grain yields were obtained with 129.3 and 88.1 kg P₂O₅ ha⁻¹, respectively. Maximum wheat grain yields were obtained with 91.8 and 99.7 kg P₂O₅ ha⁻¹ for broadcast and band application, respectively. Except for total P in leaves, nutrient uptake and yield components were not affected by the P application methods and rates. The results suggested that in soils with adequate available P levels of >15 mg kg⁻¹ (Mehlich 1 extractant), the grain yield is equivalent or superior with broadcast P application as compared with band application for both soybean and wheat crops when cultivated in a Typical Eutrorthox under subtropical conditions. However, broadcast application is considered more effective for large areas of grain cultivation, where it is necessary to sow crops uniformly within the shortest possible time.

Highlights

- Phosphorus (P) application increased soybean and wheat grain yield independently of application method.

- Broadcast P₂O₅ application in a consolidated no-till management was efficient in making P available to plants.

- The P₂O₅ application methods (broadcast and band application) did not influence the nutritional status of soybean and wheat plants.

- Our study provides a comprehensive survey of the response of application type (broadcast and band application) and P₂O₅ rates in soybean-wheat cropping system in no-till management and provides novel insights to improve the sustainability and quality of P (P₂O₅) fertilization under subtropical conditions.

Introduction

Most soils in Brazil are deficient in phosphorus (P). Associated crop plant P deficiencies are more evident in grain crops, and result in reduced plant size, low pod insertion height, and small ears (Sfredo, 2008). In no-till management (NTM), crop rotation and/or cropping system in tropical and subtropical climates became viable only with the development of technologies to increase the use efficiency of fertilizers and correctives in different soil types in these regions. In these places, the use of large amounts of fertilizer, combined with the absence of abundant reserves of phosphate rocks of sedimentary origin and the need to maximize the time between one cultivation and the next, justifies studies to investigate optimization of phosphate fertilizer use efficiency in agriculture.

Crop rotation is feasible for several crops as it provides continuous soil protection, a good amount of soil organic matter (SOM), and nutrient accumulation which are of fundamental importance for maintaining high yields over the years. No-till management involves the implantation of crops in unturned soil and their protection by cover with straw residues. In NTM, the soil is not de-structured; however, this approach is restricted to band application (Fidelis et al., 2003).

The management of phosphate fertilizer influences the availability of nutrients in the soil and in the plants (Hansel, 2019). The application method can change the speed and capacity of the fertilizer to react in the soil, which in turn can change the solubilization and P availability for the plant, influencing the use efficiency (Lacerda et al., 2015). Fertilizer application must be planned from a long-term perspective, since the cost of fertilization and the responses in productivity are subject to many uncertainties and may vary from year to year (Fixen and Halvorson, 1991; Resende et al., 2006). Among the systems adopted, there are controversies about the ways of supplying the soil with adequate P levels since soluble phosphates, while promptly making P available, have shown good results in different application forms (Borges and Mallarino, 2003).

The use of soybean (Glycine max [L.] Merril) and wheat (Triticum aestivum L.) in a single cropping system is common in southern Brazil. Both crops benefit from this management system: the higher soil fertility requirement for wheat cultivation results in greater soil nutrient uptake in the straw, mainly P, and contributes to a reduction of fertilization needed for soybean crops, while wheat and other cereals benefit from the provision of significant amounts of N, P, and K left by soybean crop residues (Stainer et al., 2000; Rossi et al., 2013).

There are controversies about the efficiency of P₂O₅ application in relation to the method used (broadcast or band application) and the appropriate rate. Therefore, the objective of this study was to determine if the method and rate of P₂O₅ application has an influence on soil fertility, grain yield (GY), yield components, and nutritional status of soybean and wheat crops grown under subtropical conditions in a Kaolinitic Typical Eutrorthox.

Materials and Methods

Site and Experimental Design

The experiment was carried out in the field under rainfed conditions in an NTM (operational for 10 years). The study was conducted across two growing seasons (2014–2015 and 2015–2016) in the same field rotation located in Londrina, Paraná State, in southern Brazil (23° 23′ 30″ LS and 51° 11′ 05″ LW) (Figure 1). The soil, classified as a loamy (710 g kg⁻¹ clay) Kaolinitic Typical Eutrorthox (USDA, soil taxonomy classification), had the following chemical properties (determined at a depth of 0–0.1 m) prior to soybean cultivation: pH (CaCl₂) = 4.9, soil organic matter (SOM) = 31.3 g kg⁻¹, P (Mehlich 1 extractant) = 3.7 mg kg⁻¹, P (Resin extractant) = 9.1 mg kg⁻¹, potassium (K⁺) (Mehlich 1) = 0.7 cmol_c kg⁻¹, calcium (Ca²⁺) (KCl, 1.0 mol L⁻¹) = 5.1 cmol_c kg⁻¹, magnesium (Mg²⁺) (KCl, 1.0 mol L⁻¹) = 1.9 cmol_c kg⁻¹, sulfur ( ${S-SO}_{4}^{2 -}$ ) = 79.1 mg kg⁻¹, aluminum (Al³⁺) = 0.1 cmol_c ‘kg⁻¹, potential acidity (H+Al) = 3.8 cmol_c kg⁻¹, cation exchange capacity (CEC) = 11.5 cmol_c kg⁻¹, base saturation (V) = 67%, boron (B) available = 0.6 mg kg⁻¹, copper (Cu) available (DTPA-TEA) = 18.3 mg kg⁻¹, iron (Fe) available (DTPA-TEA) = 101.6 mg kg⁻¹, manganese (Mn) available (DTPA-TEA) = 173.8 mg kg⁻¹, and zinc (Zn) available (DTPA-TEA) = 7.9 mg kg⁻¹. The region has a humid subtropical climate (Köppen classification: Cfa) and sequential water balance and temperature (Figure 2) in the soybean and wheat crop seasons (2014–2016) according to Thornthwaite and Mather (1955).

FIGURE 1

Figure 1. Overview of the experimental site in Londrina County, Paraná State, Brazil.

FIGURE 2

Figure 2. Sequential water balance in the soybean and wheat crop seasons (2014–2016) according to Thornthwaite and Mather (1955). PET, potential evapotranspiration withdraw, H₂O removal system, and air temperature. Londrina, Paraná State, Brazil.

The experiment was set up as a randomized block design in a 2 × 4 factorial scheme [two application methods (broadcast, total area of plot, and band application together with the seed) and four P₂O₅ rates (0, 30, 60, and 120 kg ha⁻¹, equivalent to 0, 13.1, 26.2, and 52.4 kg ha⁻¹ of P)] in a split-plot arrangement (growing season), with four replicates. For the first cultivation, the soybean cv. ‘M5947 IPRO' crop was grown in plots of 4.0 × 8.0 m with a row spacing of 0.5 m. After each soybean crop, wheat cv. ‘Pardela' was cultivated in succession with a spacing of 0.175 m to quantify the residual effects of phosphate application to the soybean crop. Phytosanitary treatments and weed control of the soybean and wheat crops were carried out according to the recommendations described in TPS (2013) and ITTT (2011).

The micronutrients (B, Cu, Fe, Mn, and Zn) were applied in the form of salts mixed with 2.0 Mg ha⁻¹ of gypsum (CaSO₄·2H₂O). In the second cultivation, fertilization of the soybean crop was based on soil chemical analysis sampled after the first wheat crop and carried out according to Moreira et al. (2019b). In both soybean cultivations, seeds were inoculated with Bradyrhizobium elkanii SEMIA 587 and SEMIA 5019 (4.0 × 10⁹ viable cells g⁻¹) and treated with a solution containing 20 g ha⁻¹ of molybdenum (Mo), 2.0 g ha⁻¹ of cobalt (Co), and 10 g ha⁻¹ of nickel (Ni). For the second soybean crop, in addition to the P₂O₅ treatments, the experimental area received 183 kg ha⁻¹ of potassium chloride [KCl (60% K₂O)] (Moreira et al., 2019b). In both cultivations, the wheat crop received 180 kg ha⁻¹ of urea (45% N) in sowing.

Evaluation of the Experiments

At the R2 reproductive stage, leaves 3 and 4 were collected at random from the apex of 30 plants to determine the nutritional status and chlorophyll level of the plants with hand-held chlorophyll meters (SPAD-502, Minolta, Japan) of the plants, afterwards, the data were transformed in chlorophyll contents (mg cm⁻²) using the equation y = 16.033 + (7.5774 × SPAD) using the equation y = 16.033 + (7.5774 × SPAD) (Moreira et al., 2020). The leaves were collected and dried under forced ventilation oven at 65 ± 5 °C until constant weight. The total N was extracted by sulfuric digestion and determined by the micro-Kjeldahl method (Nelson and Sommers, 1972). The total P, K, Ca, Mg, S, B, Cu, Fe, Mn, and Zn concentrations were extracted by digestion in a nitric-perchloric solution. The total P concentration was determined by spectrophotometry with molybdenum blue and the S concentration by turbidimetry. The total K, Ca, Mg, Cu, Fe, Mn, and Zn concentrations were analyzed by atomic absorption spectrophotometry, according to the methods described by Malavolta et al. (1997). Soil chemical analyses (pH, C, P, K, Ca, Mg, S, B, Cu, Fe, Mn, and Zn) were carried out according to the methodologies described in Teixeira et al. (2017).

At the R6 growth stage (Fehr et al., 1971), 30 soybean plants were collected from the four central rows in each plot to determine the number of pods (NP), number of grains per pod (NGP), and plant height (PH). Soybean grain yield (GY) and 100 grain mass (R8 stage) were determined from an area of 3.0 × 7.0 m, leaving a 0.5-m border on all sides of each plot.

For the wheat crop, at the 10.1 stage, leaves were collected from 30 plants at random within each treatment to determine the nutritional status (N, P, K, Ca, Mg, S, Cu, Fe, Mn, and Zn—as per Malavolta et al., 1997), and chlorophyll level (SPAD-502, Minolta, Japan), and transformed with the equation described by Uddling et al. (2007). The GY, PH, hectoliter volume (HV), number of spikelet (NS), number of spikelets per clump (NSC), and 1000 grain mass were quantified. At the end of the crop cycle (harvest), soils were sampled from each treatment at a depth of 0–0.1 m to quantify available P with Mehlich 1 and resin extractants as described by Teixeira et al. (2017).

Statistical Analysis

Normality of the distribution of the traits was tested via Shapiro–Wilk's normality test (Shapiro and Wilk, 1965), testing the null hypothesis that the sample belongs to a population with normal distribution with the statistic X (0 < X < 1); if X is equal to 1, the data perfectly fit normal distribution, whereas small values of X are evidence of deviations from normality. The value of the statistic and the associated p-value were obtained; if this p-value was less than the nominal value of significance p ≤ 0.05, the null hypothesis of normality was rejected. Normality of residuals was examined using the UNIVARIATE procedure (p ≤ 0.05). Squared and absolute values of residuals were examined with Levene's Test to confirm homogeneity of variances (p ≤ 0.05). After the normality test, the treatment variance data (soybean and wheat) were subjected to analysis of variance (ANOVA) and F-test using the standard least squares procedure of JMP by SAS. Regression and Pearson product-moment correlations were derived using the REG procedure of SAS to investigate the relation between P₂O₅ rates with grain yield (GY), yield components, nutrient concentrations, and chlorophyll level of soybean and wheat, and the regression equations for GY were chosen according to the best fit of the points within the curve. Where there was interaction of treatments × growing seasons (p ≤ 0.05), the data were separated for each growing season.

Results and Discussion

There was an interaction of P₂O₅ rate and application method on GY for both the soybean (Figure 3) and wheat (Figure 4) crops. In the two growing seasons, the yield of soybeans responded to the P₂O₅ application with a quadratic effect in the first and a linear effect in the second cultivation cycle. In addition, the results from the GY regression analysis under the effect of P₂O₅ rates indicated that, when grown in a Typic Eutrorthox, the broadcast application for the soybean crop showed an average value for the two harvests that was 5.4% higher for the maximum estimated rate (MER) of 99.8 kg ha⁻¹ of P₂O₅ for the band application at sowing method and 129.3 kg ha⁻¹ for the broadcast method with a MER of 115.0 kg ha⁻¹ for the mean of the two application methods (Figure 2). These values are beyond the official recommendation for yield 3.0 to 4.0 Mg ha⁻¹ (Moreira et al., 2019b). The wheat crop GY also increased with increasing P₂O₅ rates, but there was no effect of the application method (Figure 4). However, there was less yield in the second crop, probably due to less rain at the beginning of flowering (Figure 2). The mean GY of the two crops varied from 2,975.5 kg ha⁻¹ for the MER of 91.8 kg ha⁻¹ of P₂O₅ applied using the band application method to 2,977.7 kg ha⁻¹ for the MER of 99.7 kg ha⁻¹ of P₂O₅ applied using the broadcast approach.

FIGURE 3

Figure 3. Relationship of the effect of P₂O₅ rates applied through broadcast or band application in soybean crop, and soybean grain yield in two growing season. MER, maximum estimated rate. *Significant (F-test, p ≤ 0.05). Londrina, Paraná State, Brazil.

FIGURE 4

Figure 4. Relationship of the residual effect of P₂O₅ rates applied through broadcast or band application in soybean crop and wheat grain yield in two growing seasons. MER, maximum estimated rate. *Significant (F-test, p ≤ 0.05). Londrina, Paraná State, Brazil.

Similar results reflecting the positive effects of P₂O₅ rates on GY of soybeans and wheat have been reported by Moreira et al. (2014) and Blue et al. (1990). Regarding the P₂O₅ rates and the two application methods, Heckman and Kamprath (1992) and Borges and Mallarino (2000) obtained positive responses in soybean GY with an increase in P rates application in the soil with both application methods. The band application or broadcast method showed no difference, even with available P levels in the soil within ranges considered low and very low (Moreira et al., 2019b). A low soil P level was a characteristic of the present study, where the average P content available in the soil before planting was 3.7 mg kg⁻¹ (determined using the Mehlich 1 extractant). In a cultivation study carried out under “Cerrado” conditions, in a Typical Oxisol with a clayey texture and low P concentration, Broch and Chueiri (2005) obtained similar soybean yields when applying P fertilizer using the broadcast or the band application method. Regarding yield components, only the chlorophyll levels and number of grains per plant (NGP) in soybean crop and the number of spikelets (NS), and chlorophyll levels in wheat were significantly influenced by the P₂O₅ rates (Table 1). Fageria et al. (2011, 2013) reported that in phosphate fertilization studies, the NGP and NS have a greater influence in increasing GY as compared with other yield components, while P in the plants acts directly as an energy donor in the process of photosynthesis (Malavolta et al., 1997).

TABLE 1

Table 1. Yield components and chlorophyll level as influenced by P₂O₅ rates, under broadcast and band application, on the average of two growing seasons of a soybean–wheat cropping system, Londrina, Paraná State, Brazil.

Studies by Anghinoni (1992) and Sousa et al. (2002) have shown that rates higher than 100 kg ha⁻¹ of P₂O₅ can be applied using the broadcast method without influencing the GY. Resende et al. (2006) and Borges and Mallarino (2003) investigated P sources, P rates, and application methods in corn (Zea mays L.) and soybean crop and observed that fertilization by broadcast and/or band application with soluble P sources yielded similar productivity. Prado et al. (2001) point out that in soils with a low P concentration, the adsorption of this nutrient is maximized when the phosphate fertilizer is applied as a broadcast with incorporation. This application methods lead to a greater contact of the P with the solid phase of the soil and, at the same time, the P contact with the plant root system is reduced. According to Lana et al. (2003), where there is non-incorporation of fertilizer in the soil, P fertilization can reduce losses due to fixation, being provisionally adsorbed on SOM, which becomes available to plants after mineralization (Anderson, 1980) and provides greater root development in the topsoil (Barber, 1995). Another factor to be highlighted is that in NTM, characterized by greater SOM accumulation in the upper layers of the soil, the restrictive effect described by Novais and Smyth (1999) was not observed. This restrictive effect is that more weathered tropical soils with high fixation power represent a drain preferential of P, and the broadcast application of P₂O₅ proves inadequate to meet plant demands. Although the results of the present study indicate better plant responses with broadcast application of P₂O₅, it is worth mentioning that due to the slope of the land and excessive rain (runoff effect), the band application at sowing method can reduce fertilizer losses by laminar leaching (Borges and Mallarino, 2003).

In this study, the available P levels in the soil under soybean and wheat crops showed linear increases in Mehlich 1 (P-M1) and ion exchange resin (P-RE) in relation to P₂O₅ rates (Table 2). In Brazil, these two methods are considered official; however, the interpretations must be in accordance with the corresponding calibrations made in each region of the country. The P-M1 extractant available P-values were lower than the P-RE values extracted until 5.7 times as much available P as that of P-M1, a result similar to that observed by Moreira and Malavolta (2001). In addition, both had a high correlation coefficient (soybean: y (P-RE) = 6.628 + 3.202 (P-M1), r = 0.90, and wheat: y (P-RE) = 1.127 + 4.264 (P-M1), r = 0.86, p ≤ 0.05). These results confirm the findings of Moreira et al. (1997) and Moreira and Malavolta (2001) for these two extractants, indicating that both are efficient ways of quantifying P availability in the soil. The P available with the Resin method increased from the first to second crop in wheat (Table 2), possibly due to the positive effect of the applied gypsum recovering the P in depth, and the greater sensitivity of extractant in quantifying the availability of nutrient in the soil.

TABLE 2

Table 2. Phosphorus availability in the soil with Mehlich 1 and resin extractants after a soybean-wheat cropping system due to the P₂O₅ rates and application method (AM) (broadcast and band application), Londrina, Paraná State, Brazil.

There was a significant relationship between application methods and soybean leaf P, Ca, and S concentrations. Band application at a rate of 120 kg ha⁻¹ of P₂O₅ resulted in a higher foliar P concentration than for the other rate treatments. However, it did not differ significantly from the broadcast application at the same rate (Table 3). The both methods types, the P-M1 and P-RE were well correlated with foliar P concentration for both species (r = 0.65^* and r = 0.56^*, for soybean; r = 0.58^* and r = 0.56^*, for wheat, respectively, p ≤ 0.05), that is, the relationship between P-M1 and foliar P concentration was better than that of P-RE and foliar P concentration. Similar results were obtained by Antonangelo et al. (2019) with P rates in the same edaphoclimatic conditions. For K, there was no statistical difference; however, it was observed that P₂O₅ application using the band application method at a rate of 120 kg ha⁻¹ resulted in a 7.0% higher K concentration in soybean leaves in comparison with application using the broadcast method. Probably, there must have been an increase in root volume in depth with the P application in the band application with a consequent increase in the K concentration in the plant. The foliar Ca concentration showed a linear effect (y = 8.120 + 0.004x, R² = 0.68, p ≤ 0.05) only in the application to the broadcast, possibly due to the P source used that contained 12–14% Ca (Ribeiro et al., 1999). For leaf Mg, B, Cu, Fe, Mn, and Zn concentrations, there were no significant differences between treatments (Table 3). Rosa et al. (2015), in studies on P₂O₅ rates and application methods in soybean crop, observed that the broadcast approach was more effective at increasing the P, K, Ca, S, Mn, and Zn concentrations in the leaves than the band application method. Despite this variability, the foliar nutrient levels found in this study are within the ranges considered by Urano et al. (2006) as adequate in soybean crop. In this study, the negative P × Zn interaction was not observed in the uptake of these nutrients by plants as described by Malavolta et al. (1997) and verified by Moreira et al. (2006) in alfalfa (Medicago sativa L.) crop. Regarding the foliar nutrient concentrations of the wheat crop (Table 4), the only significant effect was of residual P₂O₅ rates on the P concentration (y = 2.080 + 0.006x, R² = 0.77, p ≤ 0.05). The leaf N, K, Ca, Mg, S, B, Cu, Fe, Mn, and Zn concentrations in wheat were not influenced by growing season and treatments, and were close to the concentrations according with Malavolta et al. (1997) and obtained by Moreira et al. (2019a) in a study on the soybean-wheat intercropping system with Cu rates cultivated under the same soil and climate conditions.

TABLE 3

Table 3. Leaf N, P, K, Ca, Mg, S, B, Cu, Fe, Mn, and Zn concentrations as influenced by P₂O₅ rates, with two application methods (AM) (broadcast or band application), on the average of two growing seasons of soybean crops, Londrina, Paraná State, Brazil.

TABLE 4

Table 4. Leaf N, P, K, Ca, Mg, S, B, Cu, Fe, Mn, and Zn concentrations as influenced by the P₂O₅ rates, with two application methods (AM) (broadcast or band application), on the average of two growing seasons of wheat crops, Londrina, Paraná State, Brazil.

Conclusions

Fertilization with a soluble phosphate source (triple superphosphate) in a soybean-wheat cropping system under NTM and subtropical conditions increased the productivity of soybeans and wheat over two growing seasons, regardless of the application method used (broadcast or band application). Given the similarity in soybean and wheat yield responses observed in this study, it is advised that the cost-benefit ratio for selecting the best method of P₂O₅ application should largely be determined by the difference in application costs and the relief of the area to be cultivated. The MER obtained for the mean of the two application methods are beyond the official recommendation for yield 3.0–4.0 Mg ha⁻¹. The soybean and wheat foliar P concentrations and the soil available P at a 0–0.1 m depth (determined using P-M1 and P-RE) were, for the most part, influenced by the P₂O₅ rates and not by the P application methods. Available P in the soil extracted by P-M1 and P-RE was highly correlated, with P-RE extracting until 5.7 times as much available P as P-M1.

Data Availability Statement

The datasets presented in this article are not readily available because no requests to access the datasets should be directed to Adônis Moreira, adonis.moreira@embrapa.br.

Author Contributions

AM: conducting the experiment, data tabulation, statistical analysis, and text writing. LM: conducting the experiment, and text writing. RP: conducting the experiment, and data tabulation. All authors contributed to the article and approved the submitted version.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

To the staff of Soil Fertility and Microbiology of Embrapa Soja for conducting the field experiments, the Santa Rita Laboratory for soil and plant analysis, and to the National Council of Technological and Scientific Development (CNPq) for providing scholarships to the first author.

References

Anderson, G. (1980). “Assessing organic phosphorus in soils,” in The role of phosphorus in agriculture, eds F. E. Khasawneh, E. C. Sample, and E. J. Kamprath (Madison, WI: American Society of Agronomy), 411–431. doi: 10.2134/1980.roleofphosphorus.c16