Fertilization strategies and liming in no-till integrated crop–livestock systems: effects on phosphorus and potassium use efficiency

ABSTRACT

In an integrated crop-livestock system (ICLS), system fertilization exploits the nutrient cycling imposed by animal grazing and increases the system efficiency. An increasingly popular approach to fertilization in southern Brazil is anticipating P and K requirements for soybeans into the pasture phase. This can increase the use efficiency of these nutrients in ICLS based on meat production in winter and soybean in summer. This study aimed to evaluate the effect of fertilization strategy, grazing and soil acidity correction on herbage and animal production, soybean yield, P and K contents in soil and plant tissue, and P and K use and economic efficiency. In 2017, a field experiment was established on an Acrisol ( Argissolo Vermelho distrófico ) double-cropped with soybean and Italian ryegrass under no-tillage. Herbage and animal production, soybean yield, available P and K contents, and P and K plant tissue status were determined. Available P and K in the soil were unaffected by grazing and fertilization strategy. Conversely, system fertilization and liming increased the P and K contents of aboveground Italian ryegrass biomass. Additionally, the available K budget in the soil was 2.7 times smaller in the integrated system with system fertilization than in the specialized system with conventional fertilization, possibly due to K fixation in non-exchangeable forms. By contrast, the available P budget in the soil was not affected by treatments and was positive with all systems. The use of ICLS increased economic return, and P and K use efficiency for protein production. System fertilization did not affect soybean yield, but it increased the total herbage production of Italian ryegrass. Despite this, sheep live weight did not increase. Using ICLS in combination with system fertilization provides an effective nutrient management strategy with a higher potential for sustainable food production when compared with conventional fertilization.

soybean yield; sheep grazing; animal production; nutrient management

INTRODUCTION

In weathered tropical and subtropical soils, phosphorus (P) and potassium (K) contents are typically low, resulting in a strong dependency on nonrenewable nutrient inputs ( Stewart et al., 2005 ; CQFS-RS/SC, 2016 ; Dhillon et al., 2019 ). Despite being one of the world’s largest food, fiber, and biofuel producers, Brazil is heavily dependent on fertilizer imports. In 2018 Brazil used approximately 35.5 million tons of fertilizers, 27.5 million (77 %) of which were imported, P mainly from Middle East countries and K from Canada ( ANDA, 2019 ). For that reason, soil fertility management needs to be better planned to increase these nutrients’ efficiency use.

Soil fertility management in tropical and subtropical areas, in general, aims to maintain adequate levels of available nutrients by building up soil P and K contents above thresholds values necessary for effective crop development ( CQFS-RS/SC, 2016 ), which usually requires using large amounts of fertilizers. A subsequent fertilization process, which maintains the soil’s nutrient status, usually requires less fertilizer input (usually the amounts needed to replenish nutrients removed for grain, fiber, or meat production in addition to losses by erosion, runoff and/or leaching) ( CQFS-RS/SC, 2016 ; Pauletti and Motta, 2019 ).

Increasing P and K use efficiency entails neutralizing soil acidity ( Scanlan et al., 2017 ). An appropriate soil pH reduces P adsorption by Fe and Al oxides, and increases P availability to plants as a result, and also promote better root growth, increasing water use efficiency, and allowing plants to uptake nutrients from deeper soil layers, increasing nutrient efficiency ( Gustafsson et al., 2012 ; Bai et al., 2017 ; Penn and Camberato, 2019 ; Alves et al., 2021 ; Bossolani et al., 2021 ). Amending soil acidity by liming also increases the cation exchange capacity (CEC) ( Huang et al., 2020 ), mainly through deprotonation of surface functional groups present in organic matter —which is especially important in sandy soils to avoid K losses by leaching. Soil P and K availability are not only affected by acidity amendments, but also by factors like the animal component in integrated crop–livestock systems (ICLS).

Globally, ICLS is used to obtain more food per unit of land ( Moraes et al., 2014 ), improve P and K cycling ( Assmann et al., 2017 ), and increase the availability of these two nutrients in the soil ( Ferreira et al., 2011 ; Deiss et al., 2016 ). Nevertheless, P and K efficiency can be further improved in ILCS, particularly if fertilizers are applied at the optimum time. Integrated crop–livestock systems typically alternates between a grain production phase with a higher nutrient exportation and a meat production phase with lower exportation. A long-term (14 years) study in Southern Brazil has shown that grain crops such as corn and soybean grown in the summer season may export up to 95 % of P and K from the soil, and sheep meat produced in winter pasture can export only 5 % ( Alves et al., 2019 ). Thus, new fertilization strategies, which exploit nutrient cycling, transfers between organic and mineral phase, and increased soil biological activity in ICLS, are being envisioned with system fertilization ( Ferreira et al., 2011 ; Deiss et al., 2016 ; Farias et al., 2020 ; Sekaran et al., 2021 ).

With conventional fertilization, P and K are supplied when the grain crops are sown as they require and export larger amounts of nutrients than the winter pasture ( CQFS-RS/SC, 2016 ). Applying P and K to the soil at this time rapidly increase their availability and uptake by soybean plants. Following soybean harvest, the soil may return to a state of decreased availability of P and K, so the winter crop may not benefit from the fertilizer initially applied. By contrast, fertilization during the winter pasture phase can increase the soil’s available P and K levels throughout the growth period. As P and K exports through meat are minimal, it is expected that the nutrient availability in the soil will be sufficient for soybean production after the grazing period. Therefore, this system-focused strategy based on periods of high and low nutrient exports may increase P and K use efficiency.

Some studies have shown that soybean yields are unlikely to respond to fertilization in soils with high P and K availability ( Boring et al., 2018 ). By contrast, anticipation fertilization can boost herbage production, especially if pasture is grazed ( Farias et al., 2020 ). In fact, grazing in pastures with well-managed fertilization increases net aboveground primary production and root production ( Souza et al., 2008 ; López-Mársico et al., 2015 ). Thus, because it boosts growth through multiple defoliation cycles, grazing can increase nutrient uptake from pasture relative to annual crops ( Ruess et al., 1983 ; Chapin and McNaughton, 1989 ).

This study aimed to assess the effects of different P and K fertilization strategies, performed at different times: at summer soybean crop sowing – conventional fertilization, or at Italian ryegrass establishment in winter – system fertilization. These fertilization strategies were combined with grazing/ungrazed pasture and lime/control. Then, we evaluate the effect of these factors on herbage and animal production, soybean yield, soil available P and K budgets, plant nutrient uptake and use, and economic efficiency after three years, in a no-tilled Acrisol in a subtropical Brazilian region.

MATERIALS AND METHODS

Experimental area and site

The field experiment started in April 2017, and it was set up at the Experimental Agronomic Station of the Federal University of Rio Grande do Sul (EEA-UFRGS) in Eldorado do Sul, Rio Grande do Sul State, Brazil (30° 05’ S, 51° 39’ W, 46 m above sea level).

The climate in the region is classified as subtropical humid (Cfa) according to Köppen’s system ( Alvares et al., 2013 ). Mean precipitation, air temperature and soil moisture during the experimental period (April 2017 to March 2020), obtained from the EEA-UFRGS mobile automatic station (Weather Watch 2000, Campbell Scientific, Inc.), are presented in figure 1 . The soil in the study site is Argissolo Vermelho distrófico ( Santos et al., 2013 ), which corresponds to Acrisol ( IUSS Working Group WRB, 2015 ); its main chemical and physical properties are summarized in table 1 .

Figure 1
Rainfall, mean air temperature and soil moisture over the three pasture (2017, 2018, and 2019) and cropping cycles (2017/2018, 2018/2019, and 2019/2020) in the ICLS evaluated in Southern Brazil.

Table 1
Soil chemical and physical properties before the beginning of the field experiment (July 2017) in an Acrisol under no-tillage in Eldorado do Sul (Rio Grande do Sul State, southern Brazil)

Soil properties	Soil layer (m)		Interpretation by CQFS-RS/SC (2016) for 0.00–0.10 m soil layer	Reference
	0.00–0.10	0.10–0.20
pH(H 2 O) (1)	3.9	4.0	Below reference value	Soil pH below the reference value for no-tillage (pH ≥6.0)
Available P (mg dm -3 ) (2)	94.2	45.8	Very high	Available P very high for a soil containing ≤200 g kg -1 of clay (30 mg dm -3 )
Available K (mg dm -3 ) (2)	96.6	68.5	High	Available K high for a CECpH 7.0 value of 7.6–15.0 cmol c dm -3 (90 mg dm −3 )
Ca 2+ (cmol c dm -3 ) (3)	1.1	0.9	Below reference value	Ca 2+ below the adequate threshold (4.0 cmol c dm -3 )
Mg 2+ (cmol c dm -3 ) (3)	0.5	0.4	Below reference value	Mg 2+ below the adequate threshold (1.0 cmol c dm -3 )
Al 3+ (cmol c dm -3 ) (3)	1.6	1.6	–	–
CECpH 7 (cmol c dm -3 ) (4)	12.4	11.8	–	–
Cation saturation (%) (5)	14.9	14.5	Below reference value	Cation saturation below the reference value (≥65 %)
Al saturation (%) (5)	48.5	53.1	Above reference value	Al saturation above the reference value (≤10 %).
Soil organic carbon (g kg -1 )	16.8	8.1	Medium	Medium soil organic carbon content (14.4–29.0 g kg -1 )
Clay (g kg -1 ) (6)	134	149	Class 4	Textural class 4 (≤200 g clay kg –1 )
Silt (g kg -1 ) (6)	239	234	–	–
Sand (g kg -1 ) (6)	627	618	–	–

• ⁽¹⁾ In (1:1, v/v) water suspension. ⁽²⁾ Available P and K extracted by Mehlich-1. ⁽³⁾ Al ³⁺ , Ca ²⁺ , and Mg ²⁺ extracted by KCl 1.0 mol L ^-1 (1:10; v/v). ⁽⁴⁾ Cation exchange capacity at pH 7.0 as calculated by combining H ⁺ , Al ³⁺ , Ca ²⁺ , Mg ²⁺ and K ⁺ . ⁽⁵⁾ Cation saturation = (Ca ²⁺ + Mg ²⁺ + K ⁺ )/(CEC _pH7 ) × 100, and Al saturation = Al ³⁺ /(CEC _effective. ) × 100. ⁽⁶⁾ Particle size distribution in clay <0.002 mm, silt 0.002–0.05 mm and sand >0.05 mm.

The experimental area has been managed under no-tillage and ICLS since 2003. In the winter season (April to September), Italian ryegrass ( Lolium multiflorum Lam.) pasture is grazed by sheep; the soil is cropped with soybean ( Glycine max ) in the summer (October to March).

Experimental design and soil management history

The experiment was established in 4.8 ha, split into 16 paddocks 0.23–0.41 ha in size. The experimental design was a completely randomized block with 4 replications in a 2 × 2 factorial system and split-plots. The first factor was animal grazing, that is, whether the land was grazed (integrated system) or ungrazed (specialized system). The second factor was P and K fertilization strategy, which involved supplying the soil with fertilizer in the soybean cropping phase (conventional fertilization) or winter pasture (Italian ryegrass) phase (system fertilization). The split-plot was the effect of acidity amendment (i.e., whether or not the soil was limed).

The effect of acidity neutralization (i.e., liming or no liming) was examined in all plots by excluding an area of 32 m ² (4 × 8 m) from limestone application to maintain the original acidity conditions (subplot). Table 1 summarizes the chemical properties of the soil at that point. Soil acidity was amended by applying limestone in the amount needed to raise the pH (H ₂ O) to 6.0, as recommended by the local Soil Fertility Committee ( CQFS-RS/SC, 2016 ). The soil was limed with 7.5 Mg ha ^-1 of dolomitic limestone [CaMg(CO ₃ ) ₂ ] with an effective neutralizing power of 72 %. Liming was performed in July 2017 on the soil surface without incorporation.

Grazing started in June or July of 2017, 2018, and 2019, and finished in October each year ( Table 2 ). Italian ryegrass was established by sowing viable seeds at a rate of 25 kg ha ^-1 with a centrifugal distributor in May of each year. Urea (45 % N) was applied at a rate of 150 kg ha ^-1 of N to all paddocks at the ryegrass V ₃ stage (3 totally expanded leaves). All treatments received N fertilization. Table 2 shows detailed information about the sheep. The stocking rate was adjusted by put-and-take method ( Mott and Lucas, 1952 ) and it was used to maintain the average sward canopy height (SCH) at 0.15 m, which provides the optimum plant structure for maximizing animal production ( Carvalho, 2013 ). The SCH was monitored at 7-day intervals using a sward stick ( Barthram, 1985 ) to measure 150 randomly chosen points per experimental unit monthly. Immediately, at the end of the pasture phase, residual ryegrass was desiccated with glyphosate herbicide before soybean was sown.

Table 2
Summary of the pasture and cropping phases in the three years of the experiment

Pasture phase	Season
	2017	2018	2019
Stocking period
Start of stocking period	June 13	June 5	July 5
End of stocking period	October 15	October 5	October 26
Grazing days	124	122	113
Italian ryegrass variety	BRS Ponteio	BRS Ponteio	BRS Ponteio
Sowing rate (kg ha –1 )	25	25	25
Nitrogen fertilization rate (kg ha -1 of N)	150	150	150
Animal information
Breed group	Corriedale	Corriedale	Corriedale
Age (months)	11	11	10
Initial live weight (kg)	25	30	31
Cropping phase	Cropping period
	2017/2018	2018/2019	2019/2020
Sowing	November 20	October 23	November 28
Harvest	April 27	April 23	April 13
Crop days	158	182	137
Soybean cultivar	DM 5958RSF IPRO	ND 5909	TMG 7063 IPRO
Sowing rate (seeds ha –1 )	255 000	255 000	255 000
Phosphorus fertilization rate (kg ha –1 of P)	30	30	25
Potassium fertilization rate (kg ha –1 of K)	58	58	67

The P and K fertilization rates were calculated from the amounts of P and K removed by soybean grains at an expected yield of 4.0 Mg ha ¹ . In the soybean crops of 2017/2018 and 2018/2019 were applied 30 kg ha ¹ of P and 58 kg ha ¹ of K, and in 2019/2020 crop season a rate of 25 kg ha ¹ of P and 67 kg ha ¹ of K was used ( Table 2 ). Fertilization rates were calculated according to CQFS-RS/SC (2016) . The application of P and K fertilizer was performed on the soil surface, both in the conventional fertilization treatment, coupled with soybean sowing, and in the system fertilization, at the establishment of the Italian ryegrass. As can be seen in table 1 , soil available P and K levels at the beginning of the experiment exceeded the critical thresholds. Table 2 shows the cultivars, sowing method, harvesting date, and sowing density for the 2017/2018, 2018/2019, and 2019/2020 seasons. Soybean was sown in rows 45 cm apart in all cropping seasons.

Assessment of soil acidity, soil available P and K, and P and K contents in plant tissue

Phosphorus and K content in soil and plant were determined at four different sampling times from 2017 to 2019. In the winter pasture phase, soil and plant samples were collected after 100 days of grazing to determine soil pH, and Ca ² , Mg ² , and Al ³ in September 2017 and 2018. In the cropping phase (January), samples were collected at the R ₁ soybean stage in the 2017/2018 and 2018/2019 cropping seasons ( Fehr and Caviness, 1977 ). The soil was sampled with an auger in the 0.00-0.20 m layer. Samples were dried in a forced-air oven at 45 °C, the larger lumps being crumbled, ground, and sieved (Ø = 2.0 mm). Soil pH was measured in aqueous suspensions (1:1 ratio, v/v). Available P and K were extracted with Mehlich-1, and Al ³ , Ca ² , and Mg ² with KCl 1.0 mol L ¹ ( Tedesco et al., 1995 ). Cation and Al saturation were calculated according to CQFS-RS/SC (2016) .

Italian ryegrass tissue was sampled after 100 days of grazing, six sub-samples from each paddock being combined into a composite sample. Each sample was obtained by clipping at ground level the blades inside a 0.25 m ² (0.5 × 0.5 m) quadrat. Soybean tissue was obtained at the R ₁ stage by cutting whole plants near ground level. Four sub-samples per paddock of 2 linear meters (2.0 × 0.45 m) were combined to obtain a composite sample. Both pasture and soybean samples were dried in a forced-air oven at 65 °C, weighed on an analytical balance, milled, and sieved (Ø = 0.5 mm). The P and K contents of plant tissue were determined after chemical digestion with H ₂ O ₂ + H ₂ SO ₄ according to Tedesco et al. (1995) .

Soybean, pasture, and animal production

Soybean yield, total herbage production, and live weight gain (LWG) per hectare were evaluated over three seasons. Soybean yield (kg ha ¹ ) was determined in five random sub-samples (2.0 × 0.45 m) from each plot (total area 4.5 m ² ). Samples were collected at physiological maturity, threshed to determine grain moisture, and yield was calculated adjusting the moisture level to 130 g kg ¹ .

Total herbage production as dry matter (kg DM ha ¹ ) was calculated as the daily herbage accumulation rate (kg DM day ¹ ) for each stocking period multiplied by the number of days of the period and that of stocking periods, the result being added to the initial herbage mass as determined one day before the start of the stocking period. Residual herbage mass at the end of the stocking cycle was estimated identically with herbage mass. The daily herbage accumulation rate was determined by using 4 grazing exclusion cages per experimental unit, the herbage mass inside each cage being clipped at ground level at 28-day intervals. On the other hand, the daily herbage accumulation rate was obtained as the difference between herbage mass in the grazing exclusion cage and pasture mass at the beginning of each stocking period divided by the number of days of the period.

Animals were weighed at the beginning and end of each grazing period (28 ± 3 days) to adjust stocking rates and monitor animal performance. The total LWG per hectare (kg ha ¹ ) was calculated as the difference between the final and initial weight of tester sheep multiplied by the number of animals per hectare and divided by the paddock area (ha).

Available soil P and K budgets

All P and K inputs via fertilizer and outputs through soybean grains and sheep meat were considered in the total budget. Phosphorus and K removal by soybean grains were calculated using the mean values adopted by CQFS-RS/SC (2016) , namely: 6.1 kg P Mg ^-1 of grain and 16.6 kg K Mg ^-1 of grain. The nutrient contents of sheep meat (LWG) were calculated according to Williams (2007) (i.e., on the assumption that the sheep removed 0.194 g kg ^-1 of P and 0.344 g kg ^-1 of K).

Soil available P and K budgets were calculated with provision for the initial and final available P and K contents in the 0.00-0.20 m soil layer, and all inputs (fertilizer) and outputs (LWG and grain biomass) ( Alves et al., 2019 ):

$SB=(FS-IS)-(IF-OGM)$ Eq. 1

in which: SB denotes soil budget; FS denotes final soil content (January 2019); IS denotes initial soil content (July 2017); IF means inputs via fertilizer; and OGM means outputs via grain and meat. All units have been converted to kg ha ^-1 .

Economic and use efficiency of P and K fertilization

Use efficiency (UE) for protein production per P and K fertilizer unit applied was calculated as 5.71 times ( Merrill and Watt, 1973 ) the N content of soybean ( CQFS-RS/SC, 2016 ). For live weight gain, carcass yield in Corriedale lambs was assumed to be 44.1 % ( Carvalho et al., 2006 ) and protein content 20.4 % ( Kremer et al., 2004 ). The P and K use efficiency were calculated according to equation 2 .

$\text{UE}={\mathrm{Prot}}_{\text{total}}/{\text{Nutr}}_{\text{applied}}$ Eq. 2

in which: Prot _total denotes the total amount of protein produced in soybean grains and sheep live weight over 3 years, and Nutr _aplied denotes the total amount of nutrient (P or K) applied via fertilizer in the same period.

Economic efficiency (EE) per fertilizer P and K unit applied was calculated from the economic return of soybean and meat production of sheep in US dollars (USD), using the average price for the previous three years ( CEPEA, 2020 ):

$\mathrm{E}\mathrm{E}={\mathrm{U}\mathrm{S}\mathrm{D}}_{\text{total}}/{\mathrm{N}\mathrm{u}\mathrm{t}\mathrm{r}}_{\text{applied}}$ Eq. 3

in which: USD _total is the total economic return from soybean grains and sheep LWG for the three-year period; and Nutr _aplied is the total amount of nutrient supplied via fertilizer in the same period.

Statistical analyses

Statistical analyses were performed with the software SAS ^® 9.4 ( SAS, 2015 ). The results were checked for normality with the Shapiro–Wilk test and variance homoscedasticity with the Levene test, both at a significance level of 5 %, prior to analysis of variance (ANOVA, p<0.05). When significant, differences between treatment means were evaluated with Tukey’s test, also at the 5 % significance level.

The effects included in the statistical model were fertilization strategy (conventional or system fertilization), grazing (specialized or integrated system), and liming (with or without). Fertilization strategy (F), grazing (G), liming (L), and the interactions F*G, F*L, G*L and F*G*L, were used as fixed effects, and block and its interactions as random effects. We use the PROC MIXED procedure and RANDON effect in SAS ^® 9.4 ( SAS, 2015 ). The models for available P and K in soil, total herbage production, LWG per area and soybean yield included the fixed effect of year (Y) and its interactions with other factors.

RESULTS

Soil acidity, and P and K contents of soil and plants

Soil acidity was affected by neither grazing nor fertilization ( Table 3 ). Liming increased soil pH (4.3 to 5.0), Ca ² (1.7 to 2.4 cmol _c dm ^-3 ), Mg ² (1.2 to 1.9 cmol _c dm ^-3 ) and cation saturation (26.1 to 37.3 %), and decreased Al ³ (1.0 to 0.7 cmol _c dm ^-3 ) and Al saturation (25.7 to 17.3 %), in the 0.00-0.20 m soil layer after 18 months ( Table 3 ).

Table 3
Acidity-related chemical properties of the soil in the 0.00-0.20 m layer 18 months after surface liming in the ICLS

Soil property	Without liming	With liming
pH(H 2 O)	4.3 ± 0.3 B	5.0 ± 0.3 A
Cation saturation (%)	26.1 ± 4.2 B	37.3 ± 5.4 A
Al saturation (%)	25.7 ± 3.9 A	17.3 ± 2.8 B
Exchangeable Ca (cmol c dm -3 )	1.7 ± 0.4 B	2.4 ± 0.5 A
Exchangeable Mg (cmol c dm -3 )	1.2 ± 0.2 B	1.9 ± 0.2 A
Exchangeable Al (cmol c dm -3 )	1.0 ± 0.1 A	0.7 ± 0.1 B

• Different letters in each row denote significant differences as per Tukey’s test (p<0.05).

Available P and K in the 0.00-0.20 m soil layer was affected by neither grazing nor fertilization strategy at any time during the sampling period ( Figure 2a ). Soil available P was lower in the cropping phase of the 2018/2019 season (53 mg dm ^-3 ) than it was in the pasture phases of 2017 (91 mg dm ^-3 ) and 2018 (83 mg dm ^-3 ), and in the cropping phase of 2017/2018 (90 mg dm ^-3 ) ( Figure 2a ). Liming had no effect on soil available P (mean of 77 mg dm ^-3 ; Figure 2b ). Available K was lower in the 2017/2018 cropping phase (76 mg dm ^-3 ) than it was in the 2017 and 2018 pasture phases (111 and 122 mg dm ^-3 , respectively), and in the 2017/2018 cropping phase (100 mg dm ^-3 ) ( Figure 2c ). Also, it was lower with liming (102 mg dm ^-3 ) than without it (90 mg dm ^-3 ) ( Figure 3c ).

Figure 3
Phosphorus (P) and potassium (K) contents in aboveground biomass as affected by fertilization strategy (conventional or system) (a, b), animal grazing (specialized or integrated) (c, d), and liming (with or without) (e, f). Asterisks denote significant differences as per Tukey’s test (p<0.05).

Figure 2
Available phosphorus (P) in the 0.00-0.20 m soil layer as affected by sampling season (a) and liming (b), and available potassium (K) in the soil as affected by sampling season (c) and liming (d) (with or without liming). The orange dotted lines represent the reference values for available P (30 mg dm -3 ) and available K (90 mg dm -3 ) in soil containing ≤200 g dm -3 clay and having a CEC pH7 value of 7.6–15 cmol c dm -3 ( CQFS-RS/SC, 2016 ). Different letters denote significant differences as per Tukey’s test (p<0.05).

Phosphorus and K contents of aboveground pasture biomass were affected by grazing, fertilization strategy and liming ( Figure 3 ). Thus, P contents in the 2017 and 2018 pasture seasons were higher with system fertilization than with conventional fertilization (3.4 vs 2.6 g kg ^-1 in 2017 and 5.8 vs 4.2 g kg ^-1 in 2018; Figure 3a ). Similar results were obtained as regards grazing. Thus, P contents were higher with the integrated system than they were with the specialized system in both pasture seasons (viz., 3.3 vs 2.6 g kg ^-1 in 2017 and 5.6 vs 4.6 g kg ^-1 in 2018; Figure 3c ). Potassium contents were higher with system fertilization than with conventional fertilization (viz., 16.5 vs 19.7 g kg ^-1 in 2017 and 31.5 vs 24.9 g kg ^-1 in 2018; Figure 3b ); also, they were higher with the integrated system than with specialized system (20.3 vs 15.3 g kg ^-1 in 2017 and 32.3 vs 24.1 g kg ^-1 in 2018; Figure 3d ). Finally, P contents in the 2018 pasture season were higher with liming than without it (5.1 vs 4.6 g kg ^-1 ; Figure 3e ), and so were K contents (29.4 vs 26.9 g kg ^-1 ; Figure 3f ).

Soybean, pasture, and animal production

Average production of Italian ryegrass herbage in the studied period (2017–2019) was higher with the integrated system than with the specialized system (8616 vs 7795 kg DM ha ^-1 ; Figure 4a ). Ryegrass production was also greater with system fertilization than conventional fertilization (8879 vs 7657 kg DM ha ^-1 ; Figure 4a ). Even so, the sheep LWG per unit area was similar to both fertilization strategies. Also, LWG was higher in the pasture phases of 2017 (300 kg ha ^-1 ) and 2018 (325 kg ha ^-1 ) than it was in 2019 (213 kg ha ^-1 ) ( Figure 5a ).

Figure 4
Total herbage production of Italian ryegrass (average of the years 2017, 2018 and 2019) as affected by animal grazing (specialized or integrated) (a) and fertilization strategy (conventional or system) (b). Different letters denote significant differences as per Tukey’s test (p<0.05).

Figure 5
Live weight gain of sheep per unit area in each pasture season (a), and soybean yield as affected by crop season (b) and liming (with or without) (c). Different letters denote significant differences as per Tukey’s test (p<0.05).

Soybean yield exhibited no substantial differences between grazing or fertilization strategies ( Figure 5b ). Moreover, it decreased over time, from 2.90 Mg ha ^-1 in the 2017/2018 season to 2.64 Mg ha ^-1 in 2018/2019 and 2.51 Mg ha ^-1 in 2019/2020 ( Figure 5a ). Also, the average soybean yield for the three cropping seasons was higher with liming than without it (2.77 vs 2.59 Mg ha ^-1 ; Figure 5c ).

Soil P and K budgets

Available P budget in soil was positive and affected by none of the treatments ( Table 4 ). On the other hand, the available K budget was negative with all treatments, but higher with the system fertilization-integrated system combination (-25.6 kg K ha ^-1 ) than it was with the conventional fertilization-specialized system, conventional fertilization–integrated system and system fertilization-specialized system combinations (-68.8, -95.3 and -108.9 kg K ha ^-1 , respectively; Table 4 ). Liming, however, had no effect on the available P and K budgets.

Table 4
Available P and K budgets in the 0.00-0.20 m soil layer as affected by fertilization strategy (conventional or system fertilization) and animal effect (specialized system or integrated system) in the ICLS

Fertilization	Grazing	Soil budget (1)
		Stock of available P	Stock of available K
		--------------------------- kg ha –1 ---------------------------
Conventional	Specialized system	49.9 ± 10.8	–68.8 ± 13.2 B
	Integrated system	30.3 ± 8.6	–95.3 ± 12.1 BC
System	Specialized system	26.3 ± 7.8	–108.9 ± 19.1 C
	Integrated system	36.9 ± 9.4	–25.6 ± 8.2 A
Mean		35.8 ns	–74.6 *

• ⁽¹⁾ Calculated as SB = (FS – IS) – (IF – OGM). SB: soil budget; FS: final soil content; IF: inputs via fertilizer; IS: initial soil content; OGM: outputs via grain and meat. ^ns: not significant at p<0.05. ^* Significant at p<0.05. Different letters distinguish the soil K budget under the effect of the interaction between fertilization strategy (conventional or system fertilization) and grazing (specialized system or ICLS) as per Tukey’s test (p<0.05).

Economic and use efficiency of P and K fertilization

Increased protein production per P and K unit applied via fertilizer resulted in an increased use efficiency with integrated system (35.1 kg P kg ^-1 and 16.4 kg K kg ^-1 ) relative to specialized system (29.7 kg P kg ^-1 and 13.9 kg K kg ^-1 ) ( Figure 6a ). Also, the increased economic return, in dollars, from soybean production and sheep LWG per P and K unit applied via fertilizer resulted in increased economic efficiency with the integrated system (49.7 USD kg ^-1 P and 23.3 USD kg ^-1 K) relative to the specialized system (30.9 USD kg ^-1 P and 14.5 USD kg ^-1 ) ( Figure 6b ).

Figure 6
Phosphorus and potassium use efficiency for protein production (a) and economic return (b) as affected by grazing. Different letters denote significant differences as per Tukey’s test (p<0.05).

DISCUSSION

Effect of soil acidity amendment by liming

Liming has well-known benefits on plant nutrition, like greater nutrients absorption such as N, P, K, Ca, Mg, and S, which leads to higher crop yields ( Goulding, 2016 ; Holland et al., 2018 , 2019 ; Bossolani et al., 2021 ). In this study, it increased the P and K contents of Italian ryegrass aboveground biomass in the 2018 pasture season ( Figures 3a and 3b ). Amending soil acidity by liming boosts root development and shoot growth ( Fageria and Nascente, 2014 ), and increased root growth facilitates the exploration of deeper soil layers, facilitating P and K uptake by plants.

Although the applied lime rate was calculated to raise the pH to 6.0, the soil pH in the 0.00-0.20 m layer only reached 5.0 after 18 months. The limited effect of surface liming in the first few centimeters of soil is recurrent in no-tillage areas ( Rheinheimer et al., 2018 ; Miotto et al., 2019 ), and proportional to the reaction time after application. Therefore, this result is explained by the short-term effect combined with the low solubility of lime. Moreover, the faster reactivity of the finer particles in the first centimeters raises the soil pH and Ca ² content, resulting in a slow dissolution of the coarser particles of lime, hampering soil acidity correction ( Scott et al., 1992 ). For these reasons, it is recommended that in areas with chemical restriction in-depth (>20-30 % Al saturation) associated with available P content below the critical level, liming with incorporation should be always carried out ( CQFS-RS/SC, 2016 ). As the initial levels of available P in the soil were above the critical level in the experimental area, liming was carried out on the soil surface. Despite not having increased the soil pH to the target value in the 0.00-0.20 m layer (pH 6.0), and although in our study Al ³ was not completely neutralized by liming this soil layer, it is still possible to verify the beneficial effect of liming on soybean productivity ( Figure 5c ) ( Holland et al., 2019) .

Despite that soybean yield was higher in the liming treatment, the average yield of the treatments in the three harvests (2.7 Mg ha ^-1 ) was below the expected 4 Mg ha ^-1 used to calculate the required P and K rates. Therefore, the residual acidity (pH values below 5.5 and Al saturation >10 %) in the liming treatment may be the limiting factor to achieve the target soybean yield ( CQFS-RS/SC, 2016 ). In addition, another factor that has been limiting crop yields in the experimental area is water restriction, as droughts are frequently documented ( Alves et al., 2020) . Some studies under ICLS previously showed that liming does not affect soybean yield in moderately acidic soils (pH 4.8, cation saturation 56 % and Al saturation 15 % in the 0.00-0.20 m deep layer) ( Martins et al., 2014a ). These conditions, however, are contrasting from those of our soil, which had an initial pH of 4.0, cation saturation of 15 %, and Al saturation of 50 % in the 0.00-0.20 m layer, all of which led to an increased soybean yield.

Effect of sheep grazing and fertilization strategy

Higher herbage production obtained with the system fertilization results from the higher P and K contents in the aboveground biomass of Italian ryegrass ( Figures 3a and 3d ). By contrast, P and K content in soybean biomass did not differ between fertilization strategies, resulting in similar grain yields. This result testifies to the high potential of system fertilization for improving nutrient availability and plant nutrition, largely due to P and K fertilizer being applied in the pasture phase. The fact that P and K supply was maintained throughout the pasture phase was a result of the heavy cycling of nutrients caused by grazing and by the sheep returning most of the nutrients ingested through dung and urine ( Ferreira et al., 2011 ; Silva et al., 2014 ; Deiss et al., 2016 ; Assmann et al., 2017 ; Alves et al., 2019 ).

Increased total production of Italian ryegrass is consistent with the increased herbage production previously found by Farias et al. (2020) . Nitrogen, P, and K fertilization usually increase total herbage production by providing a greater supply of fodder to sheep, thereby increasing meat production ( Ihtisham et al., 2018 ). Although the increase in total herbage production resulted in no substantial increase in LWG per area here ( Figure 5a ), gains are expected to become apparent in the long term. However, using pasture residues as inputs is crucial to regulate the soil C stock to maintain soil fertility and nutrient availability and have a positive effect on soybean yield in the long term ( Alves et al., 2020 ).

The fact that soybean yield failed to respond to earlier P and K fertilization may have resulted from the curve of nutrient uptake by grain crops comprising a single cycle. In addition, the soil P and K content were above the critical level, thus the likelihood of crop response is much lower ( Oliveira Junior et al., 2016 ; CQFS-RS/SC, 2016 ). This result is interesting as it is one of the first field studies demonstrating that system fertilization can be an effective strategy to enhance system production. Unlike grain crops, pasture is continuously stimulated to grow and uptake nutrients from the soil due to defoliation ( Gastal and Lemaire, 2015 ). Grazing can increase photosynthetic and growth rates in plants ( Gifford and Marshall, 1973 ), thereby increasing the requirements for nutrient uptake over several cycles of growth stimulation during grazing ( Ruess et al., 1983 ; Chapin and McNaughton, 1989 ). This is consistent with the increased ryegrass herbage production observed when P and K were supplied in the pasture phase. This result, however, should only be expected when soil P and K levels exceed a critical threshold ( Table 1 ). Our results cannot be extrapolated to soil conditions where P and K levels are below the critical level. Thus, further studies should be conducted using soils with available P and K contents below the critical level, as well as observing if other nutrients, such as N, could be managed in system fertilization. We emphasize that for a good functioning of the production system and for the system fertilization to be efficient, it is necessary to take into account some prerequisites, such as proper fertilization, soil acidity neutralization, use of the no-till system, and the adoption of ICLS ( Anghinoni and Vezzani, 2021 ).

Once the amount of K exported by grains and meat was lower than the K inputs via fertilizers, the negative soil available K budget obtained with all treatments was possibly a result of excess K accumulating in non-exchangeable and structural forms in the soil - a frequent occurrence in soils containing 2:1 clay minerals ( Watson et al., 2002 ; Berry et al., 2003 ) as found in previous studies in the same experimental area ( Alves et al., 2019 ). In fact, 2:1 clay minerals can easily fix K in their interlayer spaces ( Ernani et al., 2007 ). Also, applied K can partly migrate to soil layers below 0.20 m. In any case, the soil K budget was 2.7 times greater with the system fertilization–integrated system combination than with the conventional fertilization–specialized system combination. This result testifies the importance of animal grazing when nutrients are applied via fertilizer at an earlier time. In fact, grazing boosts root production ( López-Mársico et al., 2015 ), thereby expanding nutrient uptake volume. In addition, continuous growth stimulation by grazing increases the requirements for nutrient uptake ( Ruess et al., 1983 ; Chapin and McNaughton, 1989 ), thus hindering potential losses of K through leaching and runoff. By contrast, the soil available P budget was positive and similar to all treatments ( Table 3 ). This was largely the result of the fertilization history in the experimental field ( Alves et al., 2019 ), leading to saturation of the most P-eager sites and reducing the ability of the soil to immobilize P, converting it into less available forms as a result.

Integrated system provides greater economic returns than a specialized system ( Franzluebbers, 2007 ; Sulc and Tracy, 2007 ; Oliveira et al., 2014 ). The increased protein production (use efficiency) and economic return (economic efficiency) with the integrated system arise from incorporating animals as a new source of protein and economic return into the production system ( Oliveira et al., 2014 ; Martins et al., 2014b ; Costa et al., 2014 ). The increased protein production of integrated systems is only possible with nutrient cycling by animals; in fact, most of the nutrients ingested by grazing are returned to the soil and made available to crops ( Sanderson et al., 2013 ; Alves et al., 2019 ). Although system fertilization in ICLS increased total forage production, this did not result in higher animal production, thus not resulting in higher production and economic efficiency in these first three years of the experiment. However, it is expected that the higher total production of Italian ryegrass in the long term will benefit the system as a whole and contribute to higher system efficiency.

CONCLUSION

In this short period of time evaluated, system fertilization did not result in greater efficiency in the use of P and K, although it increased the P and K contents of pasture and total herbage production. Integrated system increases nutrient use efficiency, leading to a less negative available K budget in the soil. Neither soybean yield nor sheep live weight gain per unit was influenced by fertilization strategy, but soybean yield was increased by liming. Integrated system increased the use and economic efficiency of P and K fertilization by increasing food production per fertilizer unit. There is evidence that system fertilization provides an effective choice for promoting better use of nutrients by integrated crop-livestock systems, but future studies should continue this approach to evaluate the long-term effect of this fertilizer strategy on the efficiency of use of P and K.

REFERENCES

• Alvares CA, Stape JL, Sentelhas PC, Gonçalves JLM, Sparovek G. Köppen’s climate classification map for Brazil. Meteorol Z. 2013;22:711-28. https://doi.org/10.1127/0941-2948/2013/0507
  » https://doi.org/10.1127/0941-2948/2013/0507
• Alves LA, Ambrosini VG, Denardin LGO, Flores JPM, Martins AP, Filippi D, Bremm C, Carvalho PCF, Farias GD, Ciampitti IA, Tiecher T. Biological N₂fixation by soybeans grown with or without liming on acid soils in a no-till integrated crop-livestock system. Soil Till Res. 2021;209:104923. https://doi.org/10.1016/j.still.2020.104923
  » https://doi.org/10.1016/j.still.2020.104923
• Alves LA, Denardin LGO, Martins AP, Anghinoni I, Carvalho PCF, Tiecher T. Soil acidification and P, K, Ca, and Mg budget as affected by sheep grazing and crop rotation in a long-term integrated crop–livestock system in southern Brazil. Geoderma. 2019;351:197-208. https://doi.org/10.1016/j.geoderma.2019.04.036
  » https://doi.org/10.1016/j.geoderma.2019.04.036
• Alves LA, Denardin LGO, Martins AP, Bayer C, Veloso MG, Bremm C, Carvalho PCF, Machado DR, Tiecher T. The effect of crop rotation and sheep grazing management on plant production and soil C and N stocks in a long-term integrated crop–livestock system in Southern Brazil. Soil Till Res. 2020;203:104678. https://doi.org/10.1016/j.still.2020.104678
  » https://doi.org/10.1016/j.still.2020.104678
• Anghinoni I, Vezzani FM. Systemic Soil Fertility as product of system self-organization resulting from management. Rev Bras Cienc Solo. 2021;45:e0210090. https://doi.org/10.36783/18069657rbcs20210090
  » https://doi.org/10.36783/18069657rbcs20210090
• Assmann JM, Martins AP, Anghinoni I, Denardin LGO, Nichel GH, Costa SEVGA, Franzluebbers AJ. Phosphorus and potassium cycling in a long-term no-till integrated soybean–beef cattle production system under different grazing intensities in subtropics. Nutr Cycl Agroecosys. 2017;108:21-33. https://doi.org/10.1007/s10705-016-9818-6
  » https://doi.org/10.1007/s10705-016-9818-6
• Associação Nacional para Difusão de Adubos - ANDA. Anuário estatístico do setor de fertilizantes. São Paulo, ANDA; 2019 [cited 2020 Oct 27]. Available from: https://anda.org.br/wp-content/uploads/2019/08/Principais_Indicadores_2019.pdf .
  » https://anda.org.br/wp-content/uploads/2019/08/Principais_Indicadores_2019.pdf
• Bai J, Ye X, Jia J, Zhang G, Zhao Q, Cui B, Liu X. Phosphorus sorption–desorption and effects of temperature, pH, and salinity on phosphorus sorption in marsh soils from coastal wetlands with different flooding conditions. Chemosphere. 2017;188:677-88. https://doi.org/10.1016/j.chemosphere.2017.08.117
  » https://doi.org/10.1016/j.chemosphere.2017.08.117
• Barthram GT. Experimental Techniques: The HFRO Sward Stick. In: The Hill Farming Research Organization. Penicuik: HFRO; 1985. (Biennial Report, 29-30).
• Berry PM, Stockdale EA, Sylvester‐Bradley R, Philipps L, Smith KA, Lord EI, Watson CA, Fortune S. N, P and K budgets for crop rotations on nine organic farms in the UK. Soil Use Manage. 2003;19:112-8. https://doi.org/10.1111/j.1475-2743.2003.tb00289.x
  » https://doi.org/10.1111/j.1475-2743.2003.tb00289.x
• Boring TJ, Thelen KD, Board JE, De Bruin JL, Lee CD, Naeve SL, Ross WJ, Kent WA, Ries LL. Phosphorus and potassium fertilizer application strategies in corn–soybean rotations. Agronomy. 2018;8:195. https://doi.org/10.3390/agronomy8090195
  » https://doi.org/10.3390/agronomy8090195
• Bossolani JW, Crusciol CAC, Portugal JR, Moretti LG, Garcia A, Rodrigues VA, Fonseca MC, Bernart L, Vilela RG, Mendonça LP, Reis AR. Long-term liming improves soil fertility and soybean root growth, reflecting improvements in leaf gas exchange and grain yield. Eur J Agron. 2021;128:126308. https://doi.org/10.1016/j.eja.2021.126308
  » https://doi.org/10.1016/j.eja.2021.126308
• Carvalho PCF. Harry Stobbs Memorial Lecture: Can grazing behavior support innovations in grassland management? Trop Grasslands - Forrajes Tropicales. 2013;1:137-55. https://doi.org/10.17138/tgft(1)137-155
  » https://doi.org/10.17138/tgft(1)137-155
• Carvalho PCF, Oliveira JOR, Pontes LDS, Silveira EO, Poli CHEC, Rübensam JM, Santos RJ. Características de carcaça de cordeiros em pastagem de azevém manejada em diferentes alturas. Pesq Agropec Bras. 2006;41:1193-8. https://doi.org/10.1590/S0100-204X2006000700017
  » https://doi.org/10.1590/S0100-204X2006000700017
• Centro de Estudos Avançados em Economia Aplicada - CEPEA. Preços agropecuários. Piracicaba: CEPEA; 2020 [cited 2020 Mar 27]. Available from: https://www.cepea.esalq.usp.br/br .
  » https://www.cepea.esalq.usp.br/br
• Chapin FS, McNaughton SJ. Lack of compensatory growth under phosphorus deficiency in grazing-adapted grasses from the Serengeti Plains. Oecologia. 1989;79:551-7. https://doi.org/10.1007/BF00378674
  » https://doi.org/10.1007/BF00378674
• Comissão de Química e Fertilidade do Solo - CQFS-RS/SC. Manual de calagem e adubação para os Estados do Rio Grande do Sul e de Santa Catarina. 11. ed. Porto Alegre: Sociedade Brasileira de Ciência do Solo - Núcleo Regional Sul; 2016.
• Costa SEVGA, Souza ED, Anghinoni I, Carvalho PCF, Martins AP, Kunrath TR, Cecagno D, Balerini F. Impact of an integrated no-till crop–livestock system on phosphorus distribution, availability, and stock. Agr Ecosyst Environ. 2014;190:43-51. https://doi.org/10.1016/j.agee.2013.12.001
  » https://doi.org/10.1016/j.agee.2013.12.001
• Deiss L, Moraes A, Dieckow J, Franzluebbers AJ, Gatiboni LC, Sassaki LG, Carvalho PCF. Soil phosphorus compounds in integrated crop–livestock systems of subtropical Brazil. Geoderma. 2016;274:88-96. https://doi.org/10.1016/j.geoderma.2016.03.028
  » https://doi.org/10.1016/j.geoderma.2016.03.028
• Dhillon JS, Eickhoff EM, Mullen RW, Raun WR. World potassium use efficiency in cereal crops. Agron J. 2019;111:889-96. https://10.2134/agronj2018.07.0462
  » https://10.2134/agronj2018.07.0462
• Ernani PR, Almeida JA, Santos FC. Potássio. In: Novais RF, Alvarez V VH, Barros NF, Fontes RLF, Cantarutti RB, Neves JCL, editors. Fertilidade do solo. Viçosa, MG: Sociedade Brasileira de Ciência do Solo; 2007. p. 551-594.
• Fageria NK, Nascente AS. Management of soil acidity of South American soils for sustainable crop production. Adv Agron. 2014;128:221-75. https://doi.org/10.1016/B978-0-12-802139-2.00006-8
  » https://doi.org/10.1016/B978-0-12-802139-2.00006-8
• Farias GD, Dubeux JCB, Savian JV, Duarte LP, Martins AP, Tiecher T, Alves LA, Carvalho PFC, Bremm C. Integrated crop–livestock system with system fertilization approach improves food production and resource-use efficiency in agricultural lands. Agron Sustain Dev. 2020;40:39. https://doi.org/10.1007/s13593-020-00643-2
  » https://doi.org/10.1007/s13593-020-00643-2
• Fehr WR, Caviness CE. Stages of soybean development. Ames, Iowa: Cooperative Extension Service; 1977. (Special Report, 80). https://lib.dr.iastate.edu/specialreports/87 .
  » https://lib.dr.iastate.edu/specialreports/87
• Ferreira EVO, Anghinoni I, Andrighetti MH, Martins AP, Carvalho PCF. Ciclagem e balanço de potássio e produtividade de soja na integração lavoura-pecuária sob semeadura direta. Rev Bras Cienc Solo. 2011;35:161-9. https://doi.org/10.1590/S0100-06832011000100015
  » https://doi.org/10.1590/S0100-06832011000100015
• Franzluebbers AJ. Integrated crop–livestock systems in the southeastern USA. Agron J. 2007;99:361-72. https://doi.org/10.2134/agronj2006.0076
  » https://doi.org/10.2134/agronj2006.0076
• Gastal F, Lemaire G. Defoliation, shoot plasticity, sward structure and herbage utilization in pasture: Review of the underlying ecophysiological processes. Agriculture. 2015;5:1146-71. https://doi.org/10.3390/agriculture5041146
  » https://doi.org/10.3390/agriculture5041146
• Gifford RM, Marshall C. Photosynthesis and assimilate distribution in Lolium multiflorum Lam. following differential tiller defoliation. Aust J Biol Sci. 1973;26:517-26. https://doi.org/10.1071/BI9730517
  » https://doi.org/10.1071/BI9730517
• Goulding KWT. Soil acidification and the importance of liming agricultural soils with particular reference to the United Kingdom. Soil Use Manag. 2016;32:390-9. https://doi.org/10.1111/sum.12270
  » https://doi.org/10.1111/sum.12270
• Gustafsson JP, Mwamila LB, Kergoat K. The pH dependence of phosphate sorption and desorption in Swedish agricultural soils. Geoderma. 2012;189:304-11. https://doi.org/10.1016/j.geoderma.2012.05.014
  » https://doi.org/10.1016/j.geoderma.2012.05.014
• Holland JE, Bennett AE, Newton AC, White PJ, McKenzie BM, George TS, Pakeman RJ, Bailey JS, Fornara DA, Hayes RC. Liming impacts on soils, crops, and biodiversity in the UK: A review. Sci Total Environ. 2018;610:316-32. https://doi.org/10.1016/j.scitotenv.2017.08.020
  » https://doi.org/10.1016/j.scitotenv.2017.08.020
• Holland JE, White PJ, Glendining MJ, Goulding KWT, McGrath SP. Yield responses of arable crops to liming – An evaluation of relationships between yields and soil pH from a long-term liming experiment. Eur J Agron. 2019;105:176-88. https://doi.org/10.1016/j.eja.2019.02.016
  » https://doi.org/10.1016/j.eja.2019.02.016
• Huang Y, Sheng H, Zhou P, Zhang Y. Remediation of Cd-contaminated acidic paddy fields with four-year consecutive liming. Ecotox Environ Safe. 2020;188:109903. https://doi.org/10.1016/j.ecoenv.2019.109903
  » https://doi.org/10.1016/j.ecoenv.2019.109903
• Ihtisham M, Fahad S, Luo T, Larkin RM, Yin S, Chen. Optimization of nitrogen, phosphorus, and potassium fertilization rates for overseeded perennial ryegrass turf on dormant bermudagrass in a transitional climate. Front Plant Sci. 2018;9:487. https://doi.org/10.3389/fpls.2018.00487
  » https://doi.org/10.3389/fpls.2018.00487
• IUSS Working Group WRB. World reference base for soil resources 2014, update 2015: International soil classification system for naming soils and creating legends for soil maps. Rome: Food and Agriculture Organization of the United Nations; 2015. (World Soil Resources Reports, 106).
• Kremer R, Barbato G, Castro L, Rista L, Rosés L, Herrera V, Neirotti V. Effect of sire breed, year, sex, and weight on carcass characteristics of lambs. Small Ruminant Res. 2004;53:117-24. https://doi.org/10.1016/j.smallrumres.2003.09.002
  » https://doi.org/10.1016/j.smallrumres.2003.09.002
• López-Mársico L, Altesor A, Oyarzabal M, Baldassini P, Paruelo JM. Grazing increases belowground biomass and net primary production in a temperate grassland. Plant Soil. 2015;392:155-62. https://doi.org/10.1007/s11104-015-2452-2
  » https://doi.org/10.1007/s11104-015-2452-2
• Martins AP, Anghinoni I, Costa SEVA, Carlos FS, Nichel GH, Silva RAP, Carvalho PCF. Amelioration of soil acidity and soybean yield after surface lime reapplication to a long-term no-till integrated crop–livestock system under varying grazing intensities. Soil Till Res. 2014a;144:141-9. https://doi.org/10.1016/j.still.2014.07.019
  » https://doi.org/10.1016/j.still.2014.07.019
• Martins AP, Costa SEVA, Anghinoni I, Kunrath TR, Balerini F, Cecagno D, Carvalho PCDF. Soil acidification and basic cation use efficiency in an integrated no-till crop–livestock system under different grazing intensities. Agr Ecosyst Environ. 2014b;195:18-28. https://doi.org/10.1016/j.agee.2014.05.012
  » https://doi.org/10.1016/j.agee.2014.05.012
• Merrill AL, Watt BK. Energy value of foods: basis and derivation. Washington, DC: United States Departament of Agriculture; 1973. (Agriculture Handbook, 74).
• Miotto A, Tiecher T, Kaminski J, Brunetto G, De Conti L, Tiecher TL, Martins AP, Rheinheimer DS. Soil acidity and aluminum speciation affected by liming in the conversion of a natural pasture from the Brazilian Campos biome into no-tillage system for grain production. Arch Agron Soil Sci. 2019;66:138-51. https://doi.org/10.1080/03650340.2019.1605164
  » https://doi.org/10.1080/03650340.2019.1605164
• Moraes A, Carvalho PCF, Anghinoni I, Lustosa SBC, Costa SEVGA, Kunrath TR. Integrated crop–livestock systems in the Brazilian subtropics. Eur J Agron. 2014;57:4-9. https://10.1016/j.eja.2013.10.004
  » https://10.1016/j.eja.2013.10.004
• Mott GO, Lucas HL. The design, conduct and interpretation of grazing trials on cultivated and improved pastures. In: Proceeding of the VI International Grassland Congress; 1952; Pensylvania. Pensylvania: State College; 1952. p. 1380-95.
• Oliveira CAO, Bremm C, Anghinoni I, Moraes A, Kunrath TR, Carvalho PCF. Comparison of an integrated crop–livestock system with soybean only: Economic and production responses in southern Brazil. Renew Agr Food Syst. 2014;29:230-8. https://doi.org/10.1017/S1742170513000410
  » https://doi.org/10.1017/S1742170513000410
• Oliveira Junior AD, Castro CD, Pereira LR, Domingos CDS. Estádios fenológicos e marcha de absorção de nutrientes da soja. Londrina: Embrapa Soja; 2016.
• Pauletti V, Motta ACV. Manual de adubação e calagem para o estado do Paraná. 2nd ed. Curitiba: Sociedade Brasileira de Ciência do Solo, Núcleo Estadual Paraná; 2019.
• Penn CJ, Camberato JJ. A critical review on soil chemical processes that control how soil pH affects phosphorus availability to plants. Agriculture. 2019;9:120. https://doi.org/10.3390/agriculture9060120
  » https://doi.org/10.3390/agriculture9060120
• Rheinheimer DS, Tiecher T, Gonzatto R, Santanna MA, Brunetto G, Silva LS. Long-term effect of surface and incorporated liming in the conversion of natural grassland to no-till system for grain production in a highly acidic sandy–loam Ultisol from south Brazilian Campos. Soil Till Res. 2018;180:222-31. https://doi.org/10.1016/j.still.2018.03.014
  » https://doi.org/10.1016/j.still.2018.03.014
• Ruess RW, McNaughton SJ, Coughenour MB. The effects of clipping, nitrogen source and nitrogen concentration on the growth responses and nitrogen uptake of an East African sedge. Oecologia. 1983;59:253-61. https://doi.org/10.1007/BF00378845
  » https://doi.org/10.1007/BF00378845
• Sanderson MA, Archer D, Hendrickson J, Kronberg S, Liebig M, Nichols K, Schmer M, Tanaka D, Aguilar J. Diversification and ecosystem services for conservation agriculture: Outcomes from pastures and integrated crop–livestock systems. Renew Agr Food Syst. 2013;28:129-44. https://doi.org/10.1017/S1742170512000312
  » https://doi.org/10.1017/S1742170512000312
• Santos HG, Jacomine PKT, Anjos LHC, Oliveira VA, Lumbreras JF, Coelho MR, Almeida JÁ, Cunha TJF, Oliveira JB. 3. ed. Sistema brasileiro de classificação de solos. Brasília, DF: Embrapa; 2013.
• Scanlan CA, Brennan RF, D’Antuono MF, Sarre GA. The interaction between soil pH and phosphorus for wheat yield and the impact of lime-induced changes to soil aluminium and potassium. Soil Res. 2017;55:341-53. https://doi.org/10.1007/s13280-017-0970-2
  » https://doi.org/10.1007/s13280-017-0970-2
• Scott BJ, Conyers MK, Fisher R, Lill W. Particle size determines the efficiency of calcitic limestone in amending acidic soil. Aust J Agric Res. 1992;43:1175-85. https://doi.org/10.1071/AR9921175
  » https://doi.org/10.1071/AR9921175
• Sekaran U, Kumar S, Gonzalez-Hernandez JL. Integration of crop and livestock enhanced soil biochemical properties and microbial community structure. Geoderma. 2021;381:114686. https://doi.org/10.1016/j.geoderma.2020.114686
  » https://doi.org/10.1016/j.geoderma.2020.114686
• Silva FD, Amado TJC, Bredemeier C, Bremm C, Anghinoni I, Carvalho PCF. Pasture grazing intensity and presence or absence of cattle dung input and its relationships to soybean nutrition and yield in integrated crop–livestock systems under no-till. Eur J Agron. 2014;57:84-91. https://doi.org/10.1016/j.eja.2013.10.009
  » https://doi.org/10.1016/j.eja.2013.10.009
• Souza ED, Costa SEVGD, Lima CVSD, Anghinoni I, Meurer EJ, Carvalho PCF. Carbono orgânico e fósforo microbiano em sistema de integração agricultura–pecuária submetido a diferentes intensidades de pastejo em plantio direto. Rev Bras Cienc Solo. 2008;32:1273-82. https://doi.org/10.1590/S0100-06832008000300035
  » https://doi.org/10.1590/S0100-06832008000300035
• Statistical Analysis System - SAS. Base SAS 9.4 procedures guide. Cary: SAS Institute Inc.; 2015 [cited 2020 Apr 2020]. Available from: https://www.sas.com/en_us/software/sas9.html
  » https://www.sas.com/en_us/software/sas9.html
• Stewart WM, Dibb DW, Johnston AE, Smyth TJ. The contribution of commercial fertilizer nutrients to food production. Agron J. 2005;97:1-6. https://doi.org/10.2134/agronj2005.0001
  » https://doi.org/10.2134/agronj2005.0001
• Sulc RM, Tracy BF. Integrated crop–livestock systems in the US Corn Belt. Agron J. 2007;99:335-45. https://doi.org/10.2134/agronj2006.0086
  » https://doi.org/10.2134/agronj2006.0086
• Tedesco MJ, Gianello C, Bissani CA, Bohnen H, Volkweiss SJ. Análises de solo, plantas e outros materiais. 2. ed. Porto Alegre: Universidade Federal do Rio Grande do Sul; 1995. (Boletim técnico, 5).
• Watson CA, Bengtsson H, Ebbesvik M, Løes AK, Myrbeck A, Salomon E, Schroder J, Stockdale EA. A review of farm‐scale nutrient budgets for organic farms as a tool for management of soil fertility. Soil Use Manage. 2002;18:264-73. https://doi.org/10.1111/j.1475-2743.2002.tb00268.x
  » https://doi.org/10.1111/j.1475-2743.2002.tb00268.x
• Williams P. Nutritional composition of red meat. Nutr Diet. 2007;64:113-9. https://doi.org/10.1111/j.1747-0080.2007.00197.x
  » https://doi.org/10.1111/j.1747-0080.2007.00197.x

Edited by

• Editors: José Miguel Reichert 0000-0001-9943-2898 and Edicarlos Damacena Souza 0000-0003-3719-8615.

Publication Dates

• Publication in this collection
  06 June 2022
• Date of issue
  2022

History

• Received
  26 Aug 2021
• Accepted
  01 Apr 2022