Bioscience Journal | 2022 | vol. 38, e38073 | ISSN 1981-3163 1 Natã Balssan MOURA1 , Ivan Ricardo CARVALHO2 , José Antônio Gonzalez DA SILVA2 , Gerusa Massuquini CONCEIÇÃO2 , Leonir Terezinha UHDE2 , Jordana SCHIAVO2 , Bruno BERNARDO3 , Nathalia Dalla Corte BERNARDI4 , Tiago Silveira DA SILVA3 , Francine LAUTENCHLEGER5 , Iandeyara Nazaroff DA ROSA3 1 Graduate Program in Agricultural Engineering, Federal University of Santa Maria, Santa Maria, Rio Grande do Sul, Brazil. 2 Department of Agrarian Studies, Regional University of the Northwest of the State of Rio Grande do Sul, Ijuí, Rio Grande do S ul, Brazil. 3 Agronomy Course, Regional University of the Northwest of the State of Rio Grande do Sul, Ijuí, Rio Grande do Sul, Brazil. 4 Postgraduate Program in Plant Physiology, Federal University of Pelotas, Pelotas, Rio Grande do Sul, Brazil. 5 State University of the Midwest, Guarapuava, Paraná, Brazil. Corresponding author: Ivan Ricardo Carvalho carvalho.irc@gmail.com How to cite: MOURA, N.B., et al. Crop succession and its reflections on soybean performance. Bioscience Journal. 2022, 38, e38073. https://doi.org/10.14393/BJ-v38n0a2022-56872 Abstract The fluctuations in yield and consequently in production occurred due to climatic adversities in the main producing states of Brazil. Farming has changed over time, and past scenarios have shown high exploitation of natural resources focusing on soil tillage and conventional seeding methods. This study aimed to determine the yield performance of soybean grown under 10 consolidated crop succession systems. The experiment was conducted during the 2018/2019 crop season, before the research project entitled “Sustainable production systems with better use of biological and natural resources, with treatments arranged in a randomized block design and four replications”. The treatments consisted of the followin g predecessor crops: Avena sativa, Avena strigosa, Triticum aestivum, Secale cereale, Brassica napus, Raphanus sativus, Avena strigosa + Raphanus sativus + Vicia sativa, Fallow, Avena strigosa + Lolium multiflorum, and Triticum aestivum – Fodder. Soybean was subsequently sown across winter crops. The succession that showed superior yield was Avena strigosa + Lolium multiflorum. This attribute was established by associating taller plants with the maximization of the number of grains per p od, hundred- grain mass, grain mass, and plant dry mass; in contrast, there was a lower emphasis on plant residue. The determining attributes for soybean yield were plant stand, plant height, the number of pods per plant, and total grain mass, with contrasts among groups composed of the succession of Avena sativa, Avena strigosa, Triticum aestivum, Secale cereale, and Brassica napus, distanced from Raphanus sativus, Avena strigosa + Raphanus sativus + Vicia sativa, Fallow, Avena strigosa + Lolium multiflorum, and Triticum aestivum - Fodder. Keywords: Crop Production, Rotation, and Succession Systems. Integrated Systems. 1. Introduction Soybean is one of the main commodities grown and traded in the world. The main soybean-producing countries are the United States of America, Brazil, and Argentina, corresponding to more than 70% of the world production of this oilseed. It is extensively used in human and animal food through the production of CROP SUCCESSION AND ITS REFLECTIONS ON SOYBEAN PERFORMANCE https://orcid.org/0000-0002-3920-1171 https://orcid.org/0000-0001-7947-4900 https://orcid.org/0000-0002-9335-2421 https://orcid.org/0000-0003-2174-305X https://orcid.org/0000-0002-7500-4714 https://orcid.org/0000-0002-3363-6166 https://orcid.org/0000-0002-1154-6864 https://orcid.org/0000-0001-7795-6631 https://orcid.org/0000-0001-8283-3118 https://orcid.org/0000-0003-0219-6062 https://orcid.org/0000-0002-5945-4328 Bioscience Journal | 2022 | vol. 38, e38073 | https://doi.org/10.14393/BJ-v38n0a2022-56872 2 Crop succession and its reflections on soybean performance oil, bran, animal feed, and biodiesel. Thus, the demand for this grain is increasing, which requires efficient management practices such as predecessor crops that maximize crop yield. The soybean crop in 2018/2019 was based on a cultivated area of more than 35 million hectares, with production exceeding 115 million tons of grains and a 3.6% reduction relative to previous harvests (CONAB 2019). These fluctuations in yield and consequently in production occurred due to climatic adversities in the main producing states of Brazil. Farming has changed over time, and past scenarios have shown the exploitation of natural resources focusing on soil tillage and conventional seeding methods. However, with the advent of no-tillage, the plant residue of the predecessor species was preserved, causing a higher nutrient accumulation and cycling in the agricultural layer of the soil, increasing organic matter, and improving the physical, chemical, and biological properties of the soil. In the dynamics of crop succession from one production cycle to another, there is a period of absence of crops in the production area. Unoccupied areas are favorable environments for the growth and development of invasive plants, inoculum of diseases, surface and wind erosion, and reduced nutrient cycling. Based on this agricultural scenario commonly observed in the producing regions of Brazil, it is interesting to use cover plants, which are intermediate species to crops of agronomic interest that will allow cycling, recycling, and surface protection (Schnitzler 2017). Currently, many areas in the southern region of Brazil express the dynamics of crop succession unrelated to the principles of crop rotation, which is sufficient to strategically position service crops, promoting agronomic and environmental benefits in the production system in which soybean is consolidated nationally. Using this rationale dynamic will promote increases in the magnitude of plant waste and organic matter in this soil, nutrient release, thermal insulation, higher accumulation of available water, deep biopores, and the minimization of favorable conditions to the seed bank of invasive plants working as physical impediments, thus reducing the effects of pathogens and nematodes intrinsic to soil conditions. Based on the lack of technical-scientific studies addressing this topic and gathering information applicable to the crop level, this study aimed to determine the yield performance of soybean grown under 10 consolidated crop succession systems. 2. Material and Methods The experiment was conducted during the 2018/2019 crop season at the Regional Institute for Rural Development - IRDeR, in Augusto Pestana, RS, Brazil, located at latitude coordinates of 28º 26 ’00.99” and longitude of 54º 00 ’30.96”. The soil is classified as a typical dystrophic Red Latosol (Oxisol) (Santos et al. 2006) and the climate is characterized by Köeppen as humid subtropical. The LG 60150 RR2 IPRO™ cultivar was used with a relative maturity group of 5.0 and an indeterminate growth habit. The exper iment was conducted before the research project entitled “Sustainable production systems with better use of biological and natural resources, with treatments arranged in a randomized block design and four replications”. The experimental units consisted of 150 m² and a useful area of 146 m². The treatments consisted of the following predecessor crops: Avena sativa, Avena strigosa, Triticum aestivum, Secale cereale, Brassica napus, Raphanus sativus, Avena strigosa + Raphanus sativus + Vicia sativa, Fallow, Avena strigosa + Lolium multiflorum, and Triticum aestivum – Fodder. Soybean was subsequently sown across winter crops. For the cultivation of winter species, the recommendation was followed, as the biomass produced by different winter crops was evaluated with a square meter of the area collected when the species started their flowering period. Later, the samples were directed to the forced-air circulation oven for 72 hours at 70ºC. Soybeans were sown on October 31, 2018, with 45-cm spacings, fertilization based on the soil analysis recommendation, and control of invasive plants, pest insects, and diseases according to crop requirements. The variables analyzed include plant stand (plants per hectare), the number of pods per plant (units), the number of grains per pod (units), hundred-grain mass (grams), grain yield (kg ha-1), and harvest index (dimensionless). The data obtained were submitted to the model of assumptions, normality, and homogeneity of the residual variances and additivity of the statistical model. Subsequently, an analysis of variance at a 5% probability was used to identify variations attributed to the effects of crop succession. Based on the Bioscience Journal | 2022 | vol. 38, e38073 | https://doi.org/10.14393/BJ-v38n0a2022-56872 3 MOURA, N.B., et al. characters that differed statistically, the multiple means comparison test by Duncan at a 5% probability was used. After the need for identifying the trend of association among the characters, a linear correlation was used for all variables with significance based on the t-test at a 5% probability. Likewise, grain yield was listed as the main character, and multiple linear regression was applied with the stepwise method to select the independent variables and determine the selection target in each succession management. Considering the various characters involved in this study, multivariate tools such as canonical correlations were applied to define the interrelationships between group I, corresponding to the primary components of productivity (grain yield, hundred-grain mass, and harvest index); and group II, corresponding to the secondary components of productivity (the number of grains per plant pod, the number of pods per plant, and mass of vegetable waste from the succession). The average Euclidean algorithm was applied to obtain the distance matrix, serving as the basis for the UPGMA cluster responsible for showing the dendrogram and the main components. 3. Results The climatic data during the experiment (Figure 1) reveal that the rainfall volume accumulated 920 mm, and the minimum and maximum temperatures varied between 7.6ºC and 35.8ºC. Figure 1. Experimental climate approach - Regional Institute for Rural Development. The analysis of variance (Table 1) revealed that using different crop successions significantly affects all the variables studied at a 5% probability, which indicates a different behavior of soybean compared to other predecessor crops from the succession management (Table 2). The linear correlations (Figure 2) indicate that taller plants tend to reduce the magnitude of pods with two grains and the harvest index, which is confirmed by the field conditions whose minimum temperatures compromised these attributes. The number of grains per pod was inversely proportional to pods with only one or two grains formed. However, for pods containing three and four grains, the hundred - grain mass and biological dry mass determined this character. Plants with a higher magnitude of pods are determined by the positive manifestation of pods with two, three, and four grains, as well as the increase in plant dry and biological masses. The reduction in the final plant population per unit area explicitly affected the proportion of pods with only one and two grains. Soybean shows high compensatory capacity due to the smaller plant population and problems from the initial establishment, among others. In contrast, this plasticity is due to the increased number of reproductive nodes, flowers, pods, and grains per plant. The increase in plant population is directly related to the grain and dry mass of the plant, which is why cultivars are specific in this regard. In this scenario, the grain yield of soybeans was related to taller plants, the number of grains per pod, hundred-grain mass, grain mass, and plant dry mass, but there was a lower emphasis on vegetable residue. 0 10 20 30 40 50 60 70 80 90 100 -5 0 5 10 15 20 25 30 35 1 5 9 1 3 1 7 2 1 2 5 2 9 3 3 3 7 4 1 4 5 4 9 5 3 5 7 6 1 6 5 6 9 7 3 7 7 8 1 8 5 8 9 9 3 9 7 1 0 1 1 0 5 1 0 9 1 1 3 1 1 7 1 2 1 1 2 5 1 2 9 1 3 3 1 3 7 1 4 1 Crop cycle Rainfall Precipitation Minimum Temperature Maximum Temperature T e m p e ra tu re ( °C ) R a in fa ll ( m m ) Bioscience Journal | 2022 | vol. 38, e38073 | https://doi.org/10.14393/BJ-v38n0a2022-56872 4 Crop succession and its reflections on soybean performance Table 1. Mean square for plant stand (STD, plants/ha-1), plant height (PH), the number of pods of 1, 2, 3, and 4 grains (NP1, NP2, NP3, NP4), the number of pods per plant (NPP), the number of grains per pod per plant (NGPP), total grain mass (TGM), hundred-grain mass (HGM), plant dry mass (PDM), biological dry mass (BDM), harvest index (HI), and grain yield (YI) of soybean grown under different crop successions - UNIJUÍ, 2019. MEAN SQUARE SV D F STD PH NP 1 NP 2 NP3 NP4 NPP NG PP TGM HG M PDM BDM HI YI MVI TREATM ENT 9 10,681. 34* 1,276. 41* 7.0 6* 82. 5* 816.9 8* 2.1* 1,365. 71* 0.1 1* 322.3 9* 20.8 9* 633.9 3* 87.5 2* 0.01 2* 1,110,583 .13* 13,68 0* BLOCK 3 11,213. 82 485.01 10. 51 18. 86 30.67 2.99 55.49 0.2 0 94.22 8.93 244.3 9 15.7 8 0.00 1 2,418,667 .7 1,250 CV (%) - 16.94 7.39 119 .8 37. 72 33.00 150. 65 27.17 7.9 29.31 11.8 2 30.4 33.3 3 7.70 11.75 17.97 Mean - 295.285 7 72.79 1.5 4 15. 05 52.69 0.57 70.49 2.3 7 23.74 14.7 6 41.84 17.7 9 0.56 5368.198 7.07 Table 2. Mean comparisons for the variables of plant stand (STD, plants/ha-1), plant height (PH), the number of pods of 1, 2, 3, and 4 grains (NP1, NP2, NP3, NP4), the number of pods per plant (NPP), the number of grains per pod per plant (NGPP), total grain mass (TGM), hundred-grain mass (HGM), plant dry mass (PDM), biological dry mass (BDM), harvest index (HI), and grain yield (YI) of soybean grown under different crop successions - UNIJUÍ, 2019. TREATMENT STD PH NP 1 NP2 NP3 NP4 NPP NGP P TGM HG M PDM BDM HI YI MVI Avena sativa 272,2 22 c 67.1 5 d 1.8 5 ab 14.9 4 abc 48.1 5 b 0.60 abcd 67.3 5 b 2.29 c 20.7 3 d 13.3 6 e 37.07 cde 16,3 40 b 0.5 5 bc 4,607. 51 d 6,90 0 bcd Avena strigosa 271,1 11 c 74.0 5 c 2.9 0 a 17.5 5 ab 53.4 5 ab 0.15 d 74.0 5 ab 2.27 c 22.1 6 cd 13.7 5 de 41.53 bcde 19,3 75 ab 0.5 3 c 5,205. 63 abcd 8,50 0 ab Triticum aestivum 323,9 88 a 78.5 5 b 1.6 0 b 13.1 0 c 45.2 5 b 0.40 bcd 61.3 0 b 2.30 c 18.8 6 d 13.4 5 e 34.60 e 15,7 40 b 0.5 4 c 4,990. 55 bcd 9,50 0 ab Secale cereale 328,8 99 a 73.6 0 c 1.0 5 b 13.1 0 c 46.7 3 b 0.35 bcd 60.8 9 b 2.31 bc 18.6 3 d 14.1 5 cde 36,18 5de 16,1 90b 0.5 5 bc 4,760. 37 cd 11,1 00 a Brassica napus 310,0 00 ab 61.6 0 e 1.0 0 b 13.5 5 bc 52.5 0 ab 0.30 cd 67.3 5 b 2.35 abc 23.1 0 bcd 15.6 3 ab 38,81 0 cde 15,7 10b 0.5 9 a 4,986. 95 bcd 7,40 0 bc Raphanus sativus 286,6 77 bc 65.0 0 de 1.5 5 b 16.0 5 abc 53.8 5 ab 0.80 abc 72.4 0 ab 2.36 abc 27.2 6 ab 15.8 9 ab 44.89 abcd 17,6 30b 0.6 0 a 5,623. 72 abcd 4,80 0 cde Avena strigosa + Raphanus sativus + Vicia sativa 275,5 66 bc 63.2 1 e 1.9 5 ab 18.6 0 a 62.3 0 a 0.90 abc 83.7 5 a 2.44 ab 30.5 1 a 16.0 4 a 50.42 a 19,9 05ab 0.5 9 a 5,840. 13 ab 8,70 0 ab Fallow 283,3 33bc 79.1 5 b 1.0 0 b 15.9 5 abc 54.6 5 ab 0.30 cd 74.5 5 ab 2.44 ab 26.6 0 abc 15.3 0 abc 45.93 abc 17,4 11 b 0.5 7 ba 5,735. 17 abc 2,60 0 e Avena strigosa + Lolium multiflorum 280,0 00 bc 81.3 5 ab 1.0 5 b 15.0 5 abc 63.5 0 a 1.10 a 81.6 5 a 2.46 a 27.5 9 a 14.7 8 bcd 49.73 ab 22,1 45 a 0.5 5 bc 6,188. 43 a 4,00 0 de Triticum aestivum – Fodder 323,9 88 a 84.3 6 a 1.4 2 b 12.3 6 c 43.9 4 b 0.94 ab 60.6 8 b 2.47 a 20.5 4 d 15.3 6ab 39.09 cde 17,4 68b 0.5 5 bc 5,743. 49 abc 7,20 0 bc Bioscience Journal | 2022 | vol. 38, e38073 | https://doi.org/10.14393/BJ-v38n0a2022-56872 5 MOURA, N.B., et al. Figure 2. Linear correlations for the attributes measured in soybeans. * significant associations at a 5% probability by the t-test. ** significant associations at a 1% probability by the t-test. Canonical correlations were extracted for two groups of productivity components, with group 1 representing the primary components and group 2 the secondary components. The maximum likelihood indicated interclass significance for the groups, which proved the close relationship between the groups (r = 0.75) of characters (Table 3). Figure 3. Dendrogram obtained with the matrix of multivariate distances. The dendrogram (Figure 3) shows the treatments that have affinities regardless of the variables. Two major groups stand out: the grouped treatments in red and blue colors. Bioscience Journal | 2022 | vol. 38, e38073 | https://doi.org/10.14393/BJ-v38n0a2022-56872 6 Crop succession and its reflections on soybean performance Table 3. Predictive models to determine yield specifically by the succession crop. Predictive model Predecessor crop Explanatory variables of the model Y Intercep t STD PH NP1 NP2 NP4 NPP NGPP TGM PDM Avena sativa 5,194.1 6 - - - 1,713.8 1 - - - - - - Avena strigosa 5,606.9 4 - - - -16.84 - - - - - Triticum aestivum - 4,807.4 9 9.5 7 88.2 4 - - - - - - - Secale cereale -10,418 - - - - - 174.6 6 - 218.8 3 - Brassica napus 6,195.5 7 - 27.2 3 - - 291.76 - - - - - Raphanus sativus 2,701.4 2 - - - - - - - - 69.8 7 Avena strigosa + Raphanus sativus + Vicia sativa 6,657.4 9 - - -408.68 - - - - - - Fallow 8,553.4 2 - - - - - - - 2,367.7 5 116.3 6 - Avena strigosa + Lolium multiflorum -542.05 - - - - 465.10 428.5 1 - - - Triticum aestivum - Fodder 5,724.1 9 - - - 92.37 - 996.71 - - - - Table 4. Canonical correlations for the primary (Group 1) and secondary (Groups 2) components of soybean. Components Canonical pair Group 1 YI 0.45 HI 0.58 HGM 0.72 Group 2 NGPP 0.61 NPP 0.16 MVI -0.38 *R 0.75733 **P <0.0193 * inter-class correlation; ** significance due to the maximum odds ratio. The biplot analysis of the main components represents the variable/treatment ratio (Figure 4). 4. Discussion For the adequate development of the soybean crop, the volume of rainfall distributed throughout its cycle becomes crucial, as well as adequate temperatures. Studies by Zanon et al. (2016) highlight that soybean needs approximately 800 mm to express their productive potential. During the development cycle, the minimum temperature was a limiting factor for crop development because 40 days after sowing, the soybean was affected by temperatures of 7.6ºC concomitantly with the phenological stages V2 to V3, whose recommended optimal temperatures are close to 31ºC (Zanon et al. 2018). Studies by Wolschick et al. (2016) define that predecessor crops such as black oats, vetch, and forage turnip consortium provide a high production of dry mass and vegetable waste, regulate soil temperature and humidity, and minimize the incidence of invasive plants. For Santos et al. (2014), strategically positioned predecessor crops enhance grain yield, the number of pods per plant, the number of grains per pod, and soybean plant height. Bioscience Journal | 2022 | vol. 38, e38073 | https://doi.org/10.14393/BJ-v38n0a2022-56872 7 MOURA, N.B., et al. Figure 4. Main components obtained to define multivariate associations among the effects of crop succession of the variables measured in soybean. Figure 5. Technical presentation of the performance of soybeans grown in different crop successions. Bioscience Journal | 2022 | vol. 38, e38073 | https://doi.org/10.14393/BJ-v38n0a2022-56872 8 Crop succession and its reflections on soybean performance Hence, the secondary characters were determinants for the primary soybean yield, with high yield and harvest indexes obtained when the number of grains per pod and pods per plant of the secondary components increased and there was a lower emphasis on the residue of succession. This behavior is due to the dynamics for establishing the network of productivity components. In this dynamic, increasing production units in the plant and grains per production unit causes an expected higher yield and harvest rate. High production of plant residue in the predecessor crop may be somewhat beneficial because when obtaining a high concentration of vegetable waste from only one species, this material will extract a high amount of nutrients that soybean may need if the base and top dressing fertilizations are ineffective. The soybean yield is determined by several aspects that can be directly or indirectly affected, as shown in Table 3, by the explanatory variables of the model, which can contribute positively or negatively to the soybean yield in different treatments. For the predecessor crops of Avena sativa and Avena strigosa + Raphanus sativus + Vicia sativa, the determining variable for grain yield was NP1, which presented a negative relationship, providing a lower number of these pods, and other components were maximized, perhaps with an increase in NP2, NP3, and NP4 grains. In the treatments with Avena strigosa and Triticum aestivum - Fodder, the NP2 variable affected yield negatively because the reduction in this component allows other components to define grain yield. The opposite occurred for the second crop, as the determinant variable for grain yield was the increase in NP2. For Triticum aestivum, the variables with the highest positive contribution to grain yield were the highest STD and PH, which indirectly influence the total number of pods per plant. The predecessor crop of Secale cereale and the NPP and TGM variables showed higher participation in soybean yield, with an increased number of pods per plant and total grain mass, which corresponds directly to a higher yield. For the treatment with Brassica napus, the yield was positively influenced by an increase in PH and a reduction in NP2. Plant height affects yield because the crop emits more reproductive nodes, resulting in a larger production area. Considering the reduction of NP2, a higher number of pods is maximized with one, three, and four grains, resulting in a higher number of total pods, which directly compose soybean yield. For the Raphanus sativus crop, the PDM variable positively affected grain yield, meaning that heavier plants tend to be more productive. As this is an indirect yield characteristic, it is mainly influenced by grain weight. For the predecessor Fallow treatment, the variables that stood out the most for yield were NGPP, which is negatively influenced, and TGM, positively influenced, that is, with the reduction of NGPP, the plant can move its photoassimilates to optimize grain weight, thus increasing TGM and promoting higher yield. For Avena strigosa + Lolium multiflorum, the NP2 variables, which are positively influenced, were discarded, and NP4 negatively influenced yield, that is, this treatment showed an increase in the first variable, which was decisive for yield, and a reduction in the other variable. This characteristic is little expressed due to abiotic factors. In the group of treatments in red, stratification occurred between similar treatments. The first one was Brassica napus and Triticum aestivum due to the similarity among the variables. Brassica napus contributes to nutrient recycling and root development, as it manages and develops well in compacted soils, contributing to developing the successor crop, which in this case was soybean. The treatment with Triticum aestivum contributes to the system due to the soil cover, which regulates soil temperature and humidity because of the low decomposition. These characteristics occurred particularly for a certain similarity between the variables, especially grain yield. Another similar treatment was Avena sativa and Secale cereale. Using grasses in crop succession or rotation systems contributes to reducing surface erosion, as well as the benefits aforementioned for Triticum aestivum. In the group of treatments in blue, Fallow and Triticum aestivum – Fodder were similar. An aspect that stands out is the low presence of soil cover in the area, considering that the animals grazed in the dual- purpose wheat, resulting in less coverage. Hence, taller plants were developed, which provided a higher number of pods and grain mass and caused a similar yield. When analyzing the two large groups, they were grouped as follows: the red ones for treatmen ts with Avena sativa, Avena strigosa, Triticum aestivum, Secale cereale, and Brassica napus; and the blue ones with Raphanus sativus, Avena strigosa + Raphanus sativus + Vicia sativa, Fallow, Avena strigosa + Lolium Bioscience Journal | 2022 | vol. 38, e38073 | https://doi.org/10.14393/BJ-v38n0a2022-56872 9 MOURA, N.B., et al. multiflorum, and Triticum aestivum - Fodder. In the first group, the grouped treatments belonged to grass families, and another group had intercropped cultures with a high capacity for nutrient recycling and the presence of animals, which also helps this process. The treatments with Avena strigosa, Triticum aestivum, and Triticum aestivum - Fodder presented a closer relationship with the PH variable, highlighting the last treatment. Table 2 shows that the dual-purpose wheat cover presented the highest mean of 84.36 cm due to the low cover production in the area where the animals remained with the referred treatment, with a winter matter production of 7,200 kg/ha -1. The dry matter production values for the respective treatments refer to the final residue after the last animal grazing. This may have favored the rapid seedling emergence, growth, and establishment. Seedling establishment is one of the essential parameters for plant survival (Taiz 2017), from which plants start photosynthesis, not depending on the nutritional reserves of cotyledons. This favors the accumulation of primary biomass that contributes to an initial rapid growth (Floss 2011). Therefore, the faster the emergence, growth, and establishment of a plant, the sooner it starts photosynthesis, increasing growth, development, and the number of reproductive nodes, which promotes a higher number of pods and grain mass, and consequently higher yield. Figure 3 shows the soybean height in the referred successions. The treatments with Avena strigosa + Lolium multiflorum and Fallow showed a closer relationship with the variables of YI, BDM, NGPP, NP4, PDM, and NP3. Highlighting the YI variable, the highest mean was found in the treatment with Oats + Ryegrass with 6,188.43 kg/ha-1. This treatment also presented higher TGM, HGM, and NPP, and the latter is the most important component to promote increases in grain yield. This is due to a potential high variation, which guarantees the phenotypic plasticity of soybeans. The number of pods produced depends on the number of flowers and pods attached (Mundstock and Thomas 2005). It should be noted that the Oats + Ryegrass treatment showed the highest number of pods with four grains compared to the other treatments. The presence of animals positively affects the area (Assman 2013). In this sense, Silva et al. (2012) highlights that yield increases when analyzed in isolated points of cattle manure, with productions of up to 4.00 Mg/ha-1 compared to areas without manure, which produced 3.33 Mg/ha-1. The same authors conclude that the presence of animals increases nutrient availability, such as P and K, favoring the highest number of pods per plant and grain yield. The rainfall graph (Figure 2) shows a good rainfall distribution in the reproductive period of the crop, causing high grain production regardless of crop succession. It should also be noted that plant yield relates to photosynthesis, photorespiration, and respiration. Therefore, yield increases require higher photosynthetic and lower respiratory rates. In regions of milder temperatures at night, plants tend to have a lower respiratory rate, thus accumulating more biomass due to the net photosynthesis balance, producing a higher yield potential (Floss 2011). When observing the maximum and minimum temperatures during the crop cycle, on the 61st day of development, the temperatures varied between 17ºC and 35ºC, and the crop was approximately at the R1 to R2 stages. Another period with lower temperatures occurred on the 85th day when temperatures varied between 14ºC and 30ºC, and the crop was in the R4 development stage. Thus, the maximum temperatures did not restrict the rate of liquid photosynthesis, with no significant increase in the photorespiration rate. As the experiment has already been performed for four years in this crop succession system, changes in soil structure may occur, increasing the infiltration rate and water availability for the crop, and maintaining soil moisture due to the coverage of previous crops, which also reduces the evapotranspiration rate. The milder temperatures previously mentioned also contribute to this lower rate. Stone et al. (2006), when studying the evapotranspiration rate in common bean crops, reported that it varies according to the amount of dry matter provided by cover crops. Moreover, the difference in evapotranspiration is associated with crop development, in which the growing of beans at stages V2 to V4 on Brachiaria straw provides a 28% reduction in evapotranspiration and a 34% reduction in the R4 stage due to the crop covering the soil. The treatments with Raphanus sativus and Avena strigosa + Raphanus sativus + Vicia sativa showed a closer relationship with the variables of NPP, TGM, HGM NP2, and HI. In the first variable, the treatment that stood out the most was Avena strigosa + Raphanus sativus + Vicia sativa, with 83.75 pods (Table 2). This probably occurred because this treatment had the lowest stand compared to the others, with 275 thousand Bioscience Journal | 2022 | vol. 38, e38073 | https://doi.org/10.14393/BJ-v38n0a2022-56872 10 Crop succession and its reflections on soybean performance plants/ha. The soybean crop has a high ability to compensate for lower plant densities, mainly forming a higher number of pods per individual, but it occurs up to a limit (Procopio et al. 2013). For the TGM and HGM variables, as previously mentioned, the development of soybeans on Raphanus sativus and Avena strigosa + Raphanus sativus + Vicia sativa presented problems with the development of plant stand and height, thus reducing the number of reproductive nodes and pods but compensating for the higher grain filling rate. This can be explained by the better source/drain ratio that causes a higher grain weight, as heavier grains also result in higher plant weight. The harvest index (HI) measures plant efficiency, showing the proportion of photoassimilates directed to the grains. The higher the value, the higher the use and targeting of these photoassimilates to the grain. The treatments showed similarities for this variable because it has a smaller height, resulting in a distribution of photoassimilates more directed to the grains and not to the total dry mass of the plant. Studies have shown that the HI of a crop is highly influenced by seeding density, water and nutrient availability, and temperature (Duarte et al. 2013). The MVI variable regarded the treatments with Avena sativa, Secale cereale, and Brassica napus. Grasses are more prominent for this variable due to the amount of dry matter produced, favoring the production system such as soil cover for a longer time and maintaining humidity and soil temperature. 5. Conclusions The succession that showed superior yield was Avena strigosa + Lolium multiflorum. This attribute was established by associating taller plants with the maximization of the number of grains per pod, hundred - grain mass, grain mass, and plant dry mass; in contrast, there was a lower emphasis on plant residue. The determining attributes for soybean yield were plant stand, plant height, the number of pods per plant, and total grain mass. There were contrasts between groups composed of the succession of Avena sativa, Avena strigosa, Triticum aestivum, Secale cereale, and Brassica napus, distanced from Raphanus sativus, Avena strigosa + Raphanus sativus + Vicia sativa, Fallow, Avena strigosa + Lolium multiflorum, and Triticum aestivum - Fodder. Authors' Contributions: MOURA, N.B.: Conception, data acquisition, data interpretation, article writing; CARVALHO, I.R.: analysis and interpretation of data, writing and review of the article; DA SILVA, J.A.G.: data interpretation, writing and review of the a rticle; CONCEIÇÃO, G. M.: interpretation of data, writing and review of the article; UHDE, L. T.: data interpretation, writing and review of the article; SCHIAVO, J.: data acquisition, data interpretation, article writing; BERNARDO, B.: data acquisition, data interpretation, article writing; BERN ARDI, N. D. C.: data acquisition, data interpretation, article writing; DA SILVA, T. S.: data acquisition, data interpretation, article writing; LAUTENCHLEGER, F.: data acquisition, data interpretation, article writing DA ROSA, I. N.: data acquisition, data interpretation, article writing. All authors have read and approved the final version of the manuscript. Conflicts of Interest: The authors declare no conflicts of interest. Ethics Approval: Not applicable. Acknowledgments: The authors would like to thank the funding for the realization of this study provided by UNIJUÍ, CAPES, CNPq and FAPERGS. References ASSMANN, J.M. Estoque de carbono e nitrogênio no solo e ciclagem de nutrientes em sistema de integração soja -bovinos de corte em plantio direto de longa duração. 2013. 151 f. Thesis (PhD in Soil Science) – Universidade Federal do Rio Grande do Sul, Porto Alegre, 2013. CONAB, Companhia Nacional de Abastecimento. Acompanhamento da Safra Brasileira: 2019. Available from: http://www.conab.gov.br DUARTE, E.A.A., MELO FILHO, P.D.A. and SANTOS, R.C. Características agronômicas e índice de colheita de diferentes genótipos de amendoim submetidos a estresse hídrico. Agronomic characteristics and harvest index of different peanut genotypes submitted to water stress. Revista Brasileira de Engenharia Agrícola e Ambiental, 2013, 17(8), 843-847. https://doi.org/10.1590/S141543662013000800007 FLOSS, E. L. Fisiologia de Plantas Cultivadas. 5. Ed. Passo Fundo; UPF, 2011. MUNDSTOK, C.M. and THOMAS, A.L. Soja: fatores que afetam o crescimento e o rendimento de grãos. 2005. 31 f. Porto Alegre: Universidade Federal do Rio Grande do Sul, 2005. http://www.conab.gov.br/ https://doi.org/10.1590/S141543662013000800007 Bioscience Journal | 2022 | vol. 38, e38073 | https://doi.org/10.14393/BJ-v38n0a2022-56872 11 MOURA, N.B., et al. PROCÓPIO, S.O., et al. Plantio cruzado na cultura da soja utilizando uma cultivar de hábito de crescimento indeterminado. Amazonian Journal of Agriculturaland Environmental Sciences. 2013, 56(4), 319-325. http://dx.doi.org/10.4322/rca.2013.048 SANTOS, H.P., et al. Rendimento de grãos e características agronômicas de soja em função de sistemas de rotação de culturas. Bragantia. 2014, 73(3), 263–273. https://doi.org/10.1590/1678-4499.0136 SANTOS, H.G., et al. Sistema brasileiro de classificação de solos. 2. ed. Rio de Janeiro: Embrapa Solos, 2006. 306. SCHNITZLER, F. Desempenho da cultura da soja sob diferentes plantas de coberturas do solo. Undergraduate Thesis (Graduation i n Agronomy), Departamento de Estudos Agrários, Universidade Regional do Noroeste do Estado do Rio Grande do Sul, Ijuí, 2017. STONE, L.F., et al. Evapotranspiration of irrigated common bean under no-tillage on different cover crop mulches. Pesquisa Agropecuaria Brasileira. 2006, 41(4), 577–582. TAIZ, L. Fisiologia e Desenvolvimento Vegetal. 6. Ed. Porto Alegre: Artmed, 2017. WOLSCHICK, H.N., et al. Cobertura do solo, produção de biomassa e acúmulo de nutrientes por plantas de cobertura. Revista de Ciências Agroveterinárias. 2016, 15(2), 134–143. ZANON, A.J., STRECK, N. A. and GRASSINI, P. Climate and management factors influence soybean yield potential in a subtropical environment. Agronomy Journal. 2016, 108, 1447-1454. https://doi.org/10.2134/agronj2015.0535 ZANON, A.J., et al. Ecofisiologia da Soja Visando Altas Produtividades. 1. ed. Santa Maria. 2011. Received: 21 August 2021 | Accepted: 22 April 2022 | Published: 9 September 2022 This is an Open Access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. http://dx.doi.org/10.4322/rca.2013.048 https://doi.org/10.1590/1678-4499.0136 https://doi.org/10.2134/agronj2015.0535