1,718
Views
2
CrossRef citations to date
0
Altmetric
Agronomy & Crop Ecology

Canola production and effect on soil chemical properties in response to different residue levels from three biannual crop rotations

ORCID Icon, , , , , & show all
Pages 287-296 | Received 12 May 2020, Accepted 02 Nov 2020, Published online: 06 Dec 2020

ABSTRACT

Conservation agriculture using crop rotations and residue management is increasingly accepted in conventional agricultural systems. However, the effects on crop productivity and soil properties of different rotations, levels of use, and residue management have yet to be fully understood. Our study considered the effect of three previous crops (bread wheat, durum wheat, and corn) and four levels of residue incorporation (0%, 50%, 100%, and 200%) in a split plot design with four replicates on canola production parameters and soil chemical properties at the end of this crop. Production parameters and soil chemical properties were mostly affected by the previous crop and less by the residue incorporation rates under the conditions of this experiment; however, neither of these factors affected grain yield (ranging from 4.0 to 4.8 Mg ha−1). Canola residue production was higher (9.7%) after the corn crop, regardless of the residue incorporation rate, but the harvest index was lower after this crop. Most of the soil chemical properties were also affected, revealing increased organic matter and exchangeable K when the previous crop was corn and decreased concentrations of the cations with basic reaction when the previous crop was durum wheat.

Highlights

  • Bread wheat and durum wheat residue incorporation did not affect the grain yield of the canola crop.

  • Canola residue production was positively affected by corn residue incorporation at different rates.

Introduction

The use of crop rotations with residue incorporation at different levels is an increasingly common technique in agriculture (Hiel et al., Citation2018; Karunakaran & Behera, Citation2016; Motazedian et al., Citation2019). It has benefitted soil fertility (Hirzel et al., Citation2020; Naab et al., Citation2017; Virk et al., Citation2019), decreased water and nutrient requirements in the next crop rotation (Aulakh et al., Citation2012; Selim, Citation2019), mitigated erosive effects, and increased conservation of soil resources (Blanco-Canqui & Lal, Citation2009; Selim, Citation2019; Stella et al., Citation2019). It has also lowered pathogen pressure in the following crop (Kerdraon et al., Citation2019), reduced greenhouse gas emissions associated with the burning of residues or the exclusion of residues from the soil (Liu et al., Citation2019; Malhi et al., Citation2006; Yang et al., Citation2015). However, the use of crop rotations with residue incorporation can result in significant differences in soil properties, yield and nutritional composition of the following crop (Brar et al., Citation2013; Hiel et al., Citation2018; Hirzel et al., Citation2019, Citation2020; Liu et al., Citation2019; Rong et al., Citation2018; Virk et al., Citation2017). Furthermore, the decomposition of the residue in the soil is carried out by the microbial biomass (mainly by enzymatic activities that allow the use of C, N, and other nutrients in the residue as a substrate for its metabolism). Its activity depends on some characteristics of the residue such as size, quantity, depth of incorporation, and C:N ratio (Rong et al., Citation2018; Wei et al., Citation2015; Zhao et al., Citation2016).

In a study conducted on wheat-pea rotation with wheat residue on surface, both soil organic carbon and total nitrogen in the top 10 cm soil profile was higher compared to the conventional system (Awale et al., Citation2018).

For a field experiment with three crops (winter wheat, faba bean, and corn cultivated) for six cropping seasons, using different tillage systems (reduced tillage at 10 cm depth vs. conventional tillage with ploughing at 25 cm depth), were reported effect positives on soil chemical properties (Hiel et al., Citation2018). They also used different residue management strategies (amount of crop residue returned to the soil and incorporation vs. exportation of residues), and showed that the main differences on soil chemical properties were mostly from the tillage practice and less from the restitution or removal of residues. Hirzel et al. (Citation2020) conducted experiments with two-year crop rotations (canola-corn and bean-corn) with different residue incorporation rates of canola or bean and showed that the residue incorporation rate did not affect corn grain yield. However, all canola residue levels increased soil pH and Mg, and the highest residue level reduced soil S. In addition, the soil P concentration increased proportionally with the bean residue level, and the highest bean residue level (19.0 Mg ha−1) increased soil S.

Along the same lines, Swanepoel et al. (Citation2019) conducted an experiment for two consecutive seasons with a canola crop and three wheat residue levels (high, intermediate, and low at 5.1–6.4, 4.3–5.3, and 1.5–1.9 Mg ha−1, respectively) and two tillage systems (tine openers and disc openers). They indicated that the residue level of the previous crop affected the initial plant population and biomass production up to 60 days after sowing without affecting the yield of the last crop, and concluded that canola can tolerate some level of residue. However, low yields were obtained, which fluctuated between 0.9 and 1.3 Mg ha−1. These authors also indicated that canola would perform best when located in a rotation sequence following a low residue-producing crop or a crop that has been grazed or harvested for biomass with relatively little residue left on the field.

Residue characteristics had been shown to effect soil microbial biomass (Rong et al., Citation2018). These authors demonstrated that the 13–24 cm straw length, 10–17 cm burial depth, and the 370–650 g m−2 amount of straw were the most effective to improve soil microbial functional diversity and enzyme activity among 23 evaluated treatments.

There are few reports short-term rotations with different residue incorporation levels on both the canola crop (Brassica napus L.) yield and soil chemical properties at the end of the crop. Therefore, the objective of the present study was to determine the effect of three crop rotations and four levels of residue on grain yield, residue and dry matter production, harvest index of canola, and soil chemical properties at the end of the crop in the third year of a long-term biannual crop rotation.

Material and methods

The experiment was conducted for three consecutive seasons from 2016–2019 at the Santa Rosa Experimental Station (36º31ʹ S; 71º54ʹ W), INIA-Quilamapu, Chillán, Chile. The soil is volcanic (Melanoxerand). There is a temperate Mediterranean climate characterized by a hot dry summer and cold wet winter. Precipitation was concentrated in winter and spring with values of 605, 563, and 730 mm for the 2016, 2017, and 2018 seasons, respectively. The mean temperature was 12.8, 13.2, and 13.4°C and annual evaporation was 1.023, 1.041, and 980 mm for the 2016, 2017, and 2018 seasons, respectively.

Experiment management

The design of this long-term experiment consists of biannual rotations with three crop combinations: canola-bread wheat, canola-durum wheat, and canola-corn in which residues of the previous crop are incorporated at rates of 0%, 50%, 100%, and 200%; the basic design has been maintained over time. Given this basic design, the present article focuses on the cultivation of canola as the third crop in this biannual rotation (2018–2019 season) ().

Figure 1. Canola plot for the biannual rotation experiment

Figure 1. Canola plot for the biannual rotation experiment

At the start of the biannual rotation experiment, lime was applied at the rate of 3000 kg ha−1 prior to sowing canola in April 2016 to correct soil acidity (). There were three previous crops before canola: 1) bread wheat cv. Pandora-INIA sown on 5 July 2017 and harvested on 20 January 2018, 2) durum wheat cv. Queule-INIA sown on 20 July 2017 and harvested on 27 January 2018, and 3) corn cv. DK-469 (Dekalb) sown on 25 October 2017 and harvested on 22 April 2018.

Table 1. Soil chemical properties at the 0–0.2 m depth before study initiation (2015–2016 season) and before different residue incorporation levels of wheat, durum wheat, and corn crops (2018–2019 season)

The experimental unit was a 40 m long and 14 m wide (560 m2) plot with 0.2, 0.2, and 0.7 m inter-row spacing for bread wheat, durum wheat, and corn, respectively. Nitrogen, P (P2O5), and K (K2O) fertilization rates were 240, 120, and 120 kg ha−1 in bread wheat and durum wheat and 350, 120, and 120 kg ha−1 in corn in accordance with soil chemical properties (). In the three crops, P and K were applied 100% at sowing, while N was applied 15%, 45%, and 40% at the sowing, tillering, and flag leaf stages in bread wheat and durum wheat and 40% at sowing and 60% at the six leaves stage in corn. Fertilizer sources were urea, triple superphosphate, and potassium chloride. In addition, Mg, S, Zn, and B were applied at rates of 30:33:4:2 kg ha−1 before sowing in all crops, based on the soil chemical analysis (), with magnesium sulfate, zinc sulfate, and calcium borate fertilizers. Once the three crops were harvested, the residues were incorporated at rates of 0%, 50%, 100%, and 200% (April 2018) in the same experimental unit, and the plot was divided into four split-plots 20 m long and 7 m wide (140 m2). Residue production was 8.0, 8.5, and 16.4 Mg ha−1 for bread wheat, durum wheat, and corn, respectively. The machinery used to grind and incorporate residues were a displaceable mulcher (Tornado 310, Maschio Gaspardo, Campodarsego, Italy) and a compact disk harrow (Rubin 9, Lemken GmbH and Co. KG, Alpen, Germany), respectively.

The canola (Brassica napus L.) crop ‘Imminent-SIS’ was sown on 25 May 2018 and harvested on 15 January 2019. The inter-row spacing in the canola crop was 0.7 m. Irrigation was applied to the canola crop at the flowering stage. Total weed control was carried out and disease control was not necessary. Nitrogen, P (P2O5), and K (K2O) fertilization rates were 160, 120, and 80 kg ha−1. Phosphorous and K were applied 100% at sowing, while N was applied 50% at sowing and 50% at the 60% covering stage. Fertilizer sources were urea, triple superphosphate, and potassium chloride. In addition, Mg, S, Zn, and B were applied at rates of 30:33:4:2 kg ha−1 before sowing with magnesium sulfate, zinc sulfate, and calcium borate fertilizers.

Canola yield and residue production

The plots were harvested manually at grain maturity and threshed with a stationary thresher. Plant samples were collected from a 2.1 m2 plot area and separated as grain and above ground residue. Grain and tissue samples were oven-dried at 70°C for 72 h. In addition, the harvest index was determined as the relationship between grain production and the sum of grain and above ground residue production (total dry matter production).

Soil analysis

Composite samples were collected manually from the 0–20 cm soil depth for each treatment on the same day canola was harvested. All samples were air-dried and sieved (2 mm mesh). Soil pH was determined in 1:2.5 soil:water extracts. Soil organic C was established by Walkley-Black wet digestion (Sadzawka et al., Citation2006). Soil inorganic N (NO3-N and NH4-N) was extracted with 2 M KCl and determined by colorimetry with a segmented flux spectrophotometer (autoanalyzer, Skalar Analytical BV, Breda, The Netherlands). Soil extractable P was 0.5 M NaHCO3 (Olsen P) by the molybdate-ascorbic acid method. Exchangeable Ca, Mg, K, and Na were determined by 1 M NH4OAc extraction followed by flame spectroscopy: absorption (Ca and Mg) and emission (K and Na). Soil exchangeable Al concentration was found with 1 M KCl extraction by absorption spectroscopy. Sulfur (SO42-) was determined with calcium phosphate 0.01 M and turbidimetry.

Experimental design and statistical analysis

The experimental design was a split plot in which the main plot was the previous crop (three crops) and the split plot was the residue level (four levels) with four replicates. Results were analyzed by ANOVA and mean were separated using Tukey’s test (p = 0.05) using the SAS PROC MIXED Model procedure (SAS Institute, Cary, North Carolina, USA). For the significant interactions, contrast analysis was used to compare separately the effects of the treatment. In addition, a Pearson correlation analysis was performed for the soil chemical properties under analysis.

Results

Soil chemical properties before sowing the canola crop

The post-harvest chemical properties of the soil for the bread wheat, durum wheat, and corn crops were similar (p > 0.05) (), and without limitations for the canola crop. The increase in pH and Ca concentration from the start of the experiment in 2015 () was due to lime application. Small variations in the other soil properties were observed from the start of the crop rotations to the end of the third season (); these were not attributed to the different crop rotations because there was no significant effect of the previous crop on the soil chemical properties before sowing the canola crop (p > 0.05) ().

Significance analysis

The analysis of significance () showed a greater effect of the previous crop and limited effects of the residue level and the interaction between both factors. The previous crop affected both the residue and total dry matter production and the harvest index in the canola crop (p < 0.05); it also affected most of the analyzed soil chemical properties and the effect was highly significant (p < 0.01) (). The soil Ca concentration was significantly affected (p < 0.05) by the previous crop, while the Al concentration was not (p > 0.05) (). The residue incorporation level had a significant effect on the harvest index, soil organic matter content (p < 0.05), and K content (p < 0.01) (). The previous crop × residue interaction only affected the harvest index, soil S concentration (p < 0.05), and soil Na concentration (p < 0.01).

Table 2. Significance testing of canola grain yield, residue and dry matter production, harvest index, and soil chemical properties after the canola crop harvest as affected by three previous crops and incorporation of four residue levels

Production parameters in the canola crop

Canola grain yield was not affected by the evaluated treatments (p > 0.05) () and fluctuated between 4.0 and 4.8 Mg ha−1 (). Canola residue production varied between 9.5 and 11.7 Mg ha−1 (), and the highest mean values for the previous crop effect were obtained after the corn crop (p < 0.05) (), which increased 9.7% compared with the values for the other crop rotations. The total dry matter production of canola ranged from 13.6 to 16.2 Mg ha−1 (), and the highest mean values were also obtained as the effect of the previous crop after corn (p < 0.05) (), an increase of 7.8% compared to the other crop rotations. The harvest index in canola fluctuated between 27.9% and 30.9% (), and the highest mean values as the effect of the previous crop were obtained after bread wheat and durum wheat (p < 0.05). The effect of the interaction between the previous crop and the residue incorporation level on the harvest index () occurred when the previous crops were bread wheat and corn (). After bread wheat, the highest harvest index value in canola was achieved with 50% incorporation of the previous crop residue, surpassing the treatments with no residue incorporation or with 200% bread wheat residue (p < 0.05) (). After the corn crop, the higher harvest index values in canola were reached with the incorporation of 100% and 200% of the previous crop residue, outperforming only the control without residue incorporation (p < 0.05) ().

Table 3. Canola grain yield, residue and dry matter production, and harvest index as affected by three previous crops and incorporation of four residue levels

Soil chemical properties after canola crop harvest

On an average, the previous crop effect after corn exhibited the highest contents of both soil organic matter and K (p < 0.5) after the canola crop harvest (). The highest concentration of both available N and P were reached after corn, surpassing only the values obtained durum wheat crop (p < 0.05) (). When canola was cultivated after durum wheat, there were lower values of both soil Ca and Mg concentrations after the canola crop was harvested (p < 0.05) ().

Table 4. Soil chemical properties after the canola crop harvest as affected by three previous crops and incorporation of four residue levels

The previous crop × residue level interaction on soil Na concentration occurred when canola was cultivated after corn (). In this case, the highest Na concentration was obtained when 200% of the corn residue was incorporated, outperforming only the control without residue incorporation (p < 0.05) (). For the soil S concentration after the canola crop harvest, the previous crop × residue level interaction occurred when canola was cultivated after bread wheat and durum wheat (). After bread wheat, the highest soil S concentration was reached when 200% of the residue was incorporated, surpassing only the control without residue incorporation (p < 0.05) (). After durum wheat, the highest soil S concentration was obtained in the control without residue incorporation, outperforming only the treatment with 100% residue incorporation of the previous crop (p < 0.05) ().

Correlations between soil chemical properties determined at canola crop harvest

Following canola harvest, the soil chemical properties exhibited many correlations greater than 50% with a significant effect (p < 0.01) (). Soil pH was usually positively correlated with the exchangeable Ca and Mg concentrations in the soil, and was negatively correlated with both exchangeable Al and available S concentrations in the soil (). Available N exhibited positive correlations with exchangeable Mg, K, and Na concentrations in the soil (). Available P was positively correlated only with the exchangeable Al concentration in the soil (). In general, there was a positive correlation between the K, Ca, and Mg cations in the soil (). Exchangeable Ca and Mg concentrations in the soil were negatively correlated with both soil exchangeable Al and available S in the soil (). The exchangeable Mg concentration in the soil was positively correlated with the exchangeable Na concentration in the soil (). Finally, there was a positive correlation between soil exchangeable Al and available S ().

Table 5. Post-harvest soil chemical properties of the canola crop that showed a correlation greater than 50% and statistical significance for the harvest date

Discussion

Grain, residue, and total dry matter production in canola

In contrast to our results for canola grain yield, other authors (Govaerts et al., Citation2005; Hirzel et al., Citation2019; Pandiaraj et al., Citation2015) have described significant effects of the previous crop on the grain yield of the following crop. In a review of canola crops, Assefa et al. (Citation2018) suggest that water supply, balanced nutrition, early planting (for both winter and spring types) at a shallow depth (10–19 mm), high seeding rate (6 kg ha−1), and different rotations (canola every 3 or 4 years) are among the best management practices to increase yields. They indicate neither a greater nor lesser effect on grain yield of any crop before canola. Nevertheless, there was an effect of the previous crop on residue and total dry matter production and the canola harvest index, which concurs with findings mentioned by some authors evaluating total dry matter production of crops such as corn and durum wheat in different crop rotations (Hirzel et al., Citation2017; Retamal-Salgado et al., Citation2017). Other authors (Butterly et al., Citation2013; Hirzel et al., Citation2020; Kazemeini et al., Citation2014) have also reported the effects of the previous crop and the residue level on soil chemical properties at the end of the crop.

Grain yield was normal for the area under study (Hirzel et al., Citation2019) and is within the yield range mentioned by Assefa et al. (Citation2018). This yield was much higher than values reported by Swanepoel et al. (Citation2019) in their experiment with canola using three residue levels (high, intermediate, and low at 5.1–6.4, 4.3–5.3, and 1.5–1.9 Mg ha−1, respectively) and two tillage systems in the Mediterranean climate of the Western Cape of South Africa (soil with low organic matter and high bulk density). The difference in yield and soil physicochemical properties, is largely due to the availability of irrigation, given that there is a direct relationship between canola grain yield and crop water consumption (Assefa et al., Citation2018; Swanepoel et al., Citation2019). Similar to the indicated by Swanepoel et al. (Citation2019), our result showed lack of the effect of the residue incorporation level on canola yield.

As for residue production and the contribution of these residues to the production of the total dry matter, the highest values obtained after the corn crop can be associated with a higher amount of organic C incorporated with the corn residue and the subsequent effect on soil physicochemical properties (Li et al., Citation2018). These effects can also be associated with differences in the nutritional contribution of previous crop residue or to the effect of residue composition on soil biomass, which can stimulate the increased activity of beneficial soil microorganisms with positive effects on the productivity of the following crop (Kumar et al., Citation2018; Urra et al., Citation2018; Zhang et al., Citation2018). The highest residue and total dry matter production after the corn crop negatively affected the harvest index because there were no differences in grain yield associated with the previous crop. However, the differences in the harvest index as the effect of the previous crop × residue level interaction were inconsistent when canola was cultivated after bread wheat or corn. One would have expected a higher harvest index with the lower residue rates because a greater nutrient supply associated with the incorporation of increasing residue rates could promote more vigorous plant growth (Selim, Citation2019). The increase in the harvest index in canola with the highest incorporation rates of corn residue can be associated with a rapid nutrient supply from the corn crop (0.97%, 0.14%, 0.83%, 0.21%, 0.11%, and 0.06% of N, P, K, Ca, Mg, and S, respectively [Hirzel et al., Citation2020]). The corn crop was harvested at the end of summer and incorporated under temperature and soil moisture conditions that were favorable for soil biomass activity (Rong et al., Citation2018), and for the subsequent mineralization and nutrient supply available in the next canola crop. Corn residue had a C:N ratio that was lower than wheat residue; this could have promoted mineralization and nutrient supply from the corn residue (Fiorini et al., Citation2018).

Soil chemical properties at canola crop harvest

The effect of the previous crop on increased organic matter and K concentration is associated with the residue rate, as in the case of the corn crop, due to its supply of C and K (Hirzel et al., Citation2020). However, given the greater amount of residue applied when the previous crop was corn, an increase in the other evaluated soil nutrients would have been expected at the end of the canola crop, which was only partially observed in the concentrations of available N and P and exchangeable Ca and Mg. In the case of the P, Ca, and Mg nutrients, this can be explained by their moderate (Mg) to high (P and Ca) initial concentrations in the soil and by lime application at the start of the experiment of long-term biannual rotations, which makes it difficult to achieve a significant increase associated with the supply of nutrients with residue incorporation. As for N, determining available N is not a suitable indicator of the N reserve in the active pool or total soil N; low values of less than 20 mg kg−1 are normally obtained in this type of soil at the end of a crop cycle (Hirzel et al., Citation2020). However, after corn residue incorporation it was possible to obtain a higher concentration of available N in the soil, although not significant, associated with a lower C:N ratio of corn residue compared with bread wheat residue (Fiorini et al., Citation2018). The lower soil Ca and Mg concentrations when the previous crop was durum wheat could be explained by a greater effect of soil acidification after durum wheat residue incorporation; however, soil pH at the time of the canola harvest was similar to the three previous crops.

Another reason could be a greater extraction of Ca and Mg by the canola crop when it was sown after residue incorporation of durum wheat compared with residue incorporation of bread wheat and corn. Nutrient concentrations in the canola plant in each crop rotation were not determined. Hirzel et al. (Citation2020) have reported effects of the previous crop on changes in nutrient concentration in the grain and residue of an indicator crop (bread wheat) in the same study area. A lower concentration of base cation reactions in the soil (Ca, Mg, K, and Na) are usually generated when canola is cultivated after durum wheat; this effect could be explained by soil acidification and cation exchange reactions with elements such as (H+) and Al (Al2+ and Al3+) (Fageria & Nascente, Citation2014). However, soil pH was statistically equal after each previous crop. The effects of each type of residue on the soil microbial biomass and the reactions of this biomass with basic reaction nutrients could be studied because the type of incorporated residue generates changes in the soil microbial community (Kerdraon et al., Citation2019). The lower soil S concentration when canola was cultivated after corn can also be related to the S requirements of the microbial biomass associated with the decomposition of corn residue (Kerdraon et al., Citation2019), which was not evaluated in this experiment. The effects of the previous crop × residue level interaction on the concentrations of soil Na and S are directly related to the residue rate used when canola was cultivated after corn (Na) or after bread wheat (S), and there was an erratic effect in the S concentration when canola was cultivated after durum wheat. A higher residue rate generates a greater supply of all the nutrients both for the following crop and to increase soil bioavailable reserves (Hirzel et al., Citation2020; Selim, Citation2019).

Soil chemical property relationship

Regarding the correlations between soil chemical properties, it was expected that there would be positive correlations between pH and the elements Ca and Mg (Havlin et al., Citation1999; Kunito et al., Citation2016; De Vargas et al., Citation2019) and between base cation reactions (Ca, Mg, K, and Na) (Fageria et al., Citation2014; Fageria & Nascente, Citation2014; De Vargas et al., Citation2019). This would occur when fertilization is adjusted to the soil chemical properties and there are negative correlations between pH and Al concentrations or between base cation reactions and Al (Fageria & Nascente, Citation2014; Kunito et al., Citation2016; De Vargas et al., Citation2019). However, the negative correlation between pH and S obtained under this experimental condition cannot be generalized to all soils because the highest or lowest concentration of soil available S depends on several soil formation factors (Havlin et al., Citation1999). For volcanic soils, such as the one used in the present experiment and for other soils, a negative correlation between the exchangeable Al and P concentrations would have been expected (Rampim et al., Citation2013); however, there is more evidence of a negative correlation between available P and extractable Al in the soil (Agbenin, Citation2003; Havlin et al., Citation1999). The negative correlations found between the S concentration and base cation reactions (Ca, Mg, K) are associated with the negative correlation between pH and the S concentration; in this soil, this can be explained by the S requirement of the microbial biomass in the mineralization processes of the nutrients supplied with residue incorporation (Havlin et al., Citation1999; Kerdraon et al., Citation2019). Likewise, the positive correlation between exchangeable Al and available S can be caused by the increased bioavailable concentration of both nutrients when soil pH decreases (Fageria & Nascente, Citation2014; Havlin et al., Citation1999; Kunito et al., Citation2016). The positive correlations between the available N concentration and the K, Mg, and Na concentrations can be caused by the proportional increase of these nutrients associated with the residue rate in use and the nutrient content of these residues (Hirzel et al., Citation2020).

Conclusions

Under the conditions of the present experiment, the bread wheat, durum wheat, and corn crops and their residue incorporation rates did not affect the grain yield of the canola crop. Canola residue production was positively affected by corn residue incorporation at different rates. The canola harvest index was negatively affected by the corn crop, but the values obtained presented an increase directly proportional to the corn residue level incorporate.

The use of different crops and their residue incorporation levels before the canola crop affected most of the soil chemical properties. We have highlighted the increase in organic matter and exchangeable K when the previous crop was corn and a decrease in the concentration of base cation reactions when the previous crop was durum wheat.

Finally, the production parameters and soil chemical properties in the canola crop were mostly affected by the previous crop and to a lesser extent by the residue incorporation rates.

Disclosure statement

No potential conflict of interest was reported by the authors.

Additional information

Funding

This work was supported by the Instituto de Investigaciones Agropecuarias Chile [INIA].

References