Introduction
The dry biomass of crops is an important source of renewable energy, mainly thermal and electrical. The Campos Gerais region, in Paraná State, stands out for its high agricultural production, such as grain production in the summer, mainly maize and soybean crops.
Maize and soybean crops produce large amounts of residual biomass, which can be used both for a no-tillage system and energy production through direct straw combustion (Xu et al., 2020).
Even though there are doubts about removing soil cover for energy generation (Moebius-Clune et al., 2008), the partial removal of biomass has irrelevant effects concerning the maintenance of organic matter content and soil conservation.
The maximum crop potential productive is reached if some factors follow the crop requirements. In addition to the nutrients and pests management, aspects such as temperature, precipitation, and solar radiation, must meet the crop requirements. Incident solar radiation (Rs) is the main energy source of in the environment, fundamental for the photosynthesis and plant development, and is vital in several areas (Chukwujindu, 2017). Understanding the relationships between solar radiation, temperature, and water in the soil are essential to obtain great quality and crop productivity (Palencia et al., 2013).
Agricultural planning is the best option to achieve the genetic potential of crops for maximum biomass production. The use of mathematical models and functions is an excellent alternative to estimate plant production, enabling the assessment of plant growth dynamics and final plant growth. Therefore, it is essential to study agricultural systems using models and data that characterize how a given system behaves in specific circumstances. The models allow the use of computational tools that help to identify the best options and the associated risks (Jones et al., 2017).
Given the context presented, the objective of this study was to calibrate and validate recommended models to estimate dry biomass of maize and soybean crops in the Campos Gerais region, located in the subtropical zone in southern Brazil, with incident solar radiation (Rs), growing degree days (GDD), days after planting (DAP) and precipitation (P) as variables.
Material and Methods
The study area is located in the Campos Gerais region, Paraná State, in two localities: i) Ponta Grossa, with an Oxisol, sandy clay texture, gently undulating relief and humid subtropical climate (Cfb); and ii) Castro, with an Oxisol, clay texture, 9.6% slope and humid subtropical climate (Cfb) (Alvares et al., 2013). Ponta Grossa and Castro have an average temperature of 17.5 °C and 16.9 °C and precipitation of 1495 mm year1 and 1553 mm yeaH, respectively.
The dry biomass and climatic data for soybean and maize crops were provided from experiments carried out in experimental stations from ABC Foundation Agricultural Research and Development, at 2006/2007 and 2007/2008 harvests. For each crop, four experimental areas were installed.
The cultivars used, sowing dates, and locations of each harvest are shown in Table 1. In the experimental areas, four samples of whole plants were collected for aerial biomass determination in the following phenological growth stages: four fully expanded leaf (V4), eight fully expanded leaf (V8), full flowering (FF), and physiological maturity (PM). For all crops, plants were collected approximately 0.01 m from the soil surface.
At the last stage, the plants were collected and divided into grains and aerial biomass. After separation, the plants were weighed in the field to quantify the fresh biomass. All the collected and weighed material was sent to the laboratory for drying in an oven at 70 °C until constant weight, for dry biomass quantification, for the first two phenological stages (V4 and V8). For the final two stages (FF and PM), due to the volume of material collected, a sample of whole plants weighed fresh was separated and sent to the laboratory for drying. The dry biomass per hectare was calculated with the collection area's total fresh biomass data and the sample moisture values.
The parameters of Streibig (1988), Gompertz and Philos (1825), Ratkowsky (1983), Brody (1945), Bertalanffy (1957), and García y García (2002) models were adjusted by the Least Squares Method and using a spreadsheet Solver. The calibration was performed with three harvests from different locations and years. Subsequently, each model was validated with harvests not used in the calibration. The models consider the input variables: days after planting (DAP), accumulated growing degree days (GDD), accumulated incident solar radiation (Rs), and accumulated precipitation (P).
Non-linear logistic model, proposed by Streibig (1988):
Non-linear model, proposed by Gompertz and Philos (1825):
Non-linear model, proposed by Ratkowsky (1983):
Non-linear model, proposed by Brody (1945).
Non-linear model, proposed by Bertalanffy (1957):
Sigmoidal model, proposed by García y García (2002):
Where: y(xi) - total dry biomass as a function of variable x (kg ha−1); xi - an independent variable value, evaluated at each i-days after planting (DAP; day), accumulated degree days until each i-day (GDDi; °C), incident solar radiation accumulated until the each i-day (Rsi; MJ m-2 day-1) and accumulated precipitation up to each i-day (Pi; mm day-1); "a," "b" and "c" - parameters of the equation (a: kg ha-1; b: dimensionless; c: Eq. (1) and Eq. (6) have the same unit as the independent variable, and Eq. (2) to Eq. (5) have the inverse of the same unit as the independent variable).
The accumulated degree days were calculated with the sum of:
Where: ADDi - accumulated degree days at each i-day (°C); DDi - degree days at each i-day.
The degree days were calculated using the method proposed by Arnold (1959), which considers only the lower basal temperature (Tbasal) when calculating the degree-days:
Where: DDi - degree days at each i-day (°C); Tmaxi - maximum temperature recorded at each i-day (°C); Tmini - minimum temperature recorded at each i-day (°C); Tbasal - lower basal temperature of the crop (°C).
The minimum temperatures required for the development of maize and soybean crops were considered to be 10 °C (Sánchez et al., 2014) and 13 °C (INMET, 2009), respectively.
The associations between observed (in the field; kg ha-1) and estimated (with models; kg ha-1) biomass values were verified in simple linear regression analyzes. For comparison, absolute (AE) and relative (RE) errors, Pearson correlation coefficient (r), indexes of agreement "d" and performance "c" were calculated (Souza, 2018). The "c" index values (c = r . d) indicate the following performance of the tested model: "Excellent" ("c" > 0,85); "Very good" (0,76 ≤ "c" ≤ 0,85); "Good" (0,66 ≤ "c" ≤ 0,75); "Medium" (0,61 ≤ "c" ≤ 0,65);
"Tolerable" (0,51≤ "c" ≤ 0,60); "Bad" (0,41 ≤ "c"≤ 0,50); and, "Terrible" ("c" ≤ 0,40).
Results and Discussion
For soybean crops, the harvest that had the most porlonged period of days after planting, incident solar radiation, and precipitation to reach the physiological maturation stage were not the ones that obtained the highest final dry biomass productivity (Table 2). The same did not occur with the accumulated degree days (GDA).
At the V4 phenological stage, the highest biomass productivity occurred in the "Ponta Grossa 2006/2007" harvest for soybean, and the harvest presented: fewer days after planting to change the stage (29 days); lower accumulated incident solar radiation (608.40 MJ m-2 day-1); less accumulated precipitation (96.20 mm); and, the second-lowest accumulated degree days (292.75 °C). At the V8 soybean stage, there were the same characteristics observed at the V4. Except for the accumulated degree days (GDD), the exact characteristics of the variables were observed in the full flowering stage (FF), in which the "Ponta Grossa 2006/2007" harvest showed the lowest variable value. At this same stage and harvest, the value of dry biomass was the second-lowest concerning other crops. In the physiological maturation stage (PM), the highest values of all variables observed in the soybean harvest (Ponta Grossa 2007/2008) occurred and, consequently, the highest productivity of dry biomass. In the "Ponta Grossa 2006/2007" soybean harvest, the second-highest dry biomass yield was obtained, and the physiological maturation stage (PM) occurred in a shorter period. The same stage and harvest showed the second-highest value of accumulated degree days (1294.30 °C) between harvests.
Except for "Castro 2007/2008" harvest, the accumulated degree days (GDD) to reach physiological maturation were higher than those found in the literature. Wazilewski et al. (2011) observed 1130 °C of accumulated degree days to reach the physiological maturation stage in a research carried out with soybeans in Cascavel city, Paraná State. Souza et al. (2010) observed a complete phenological cycle with 1200 °C of accumulated degree days evaluating an early soybean cultivar in northeastern Pará State.
At the V8 stage (Table 2) in the "Castro 2006/2007 and 2007/2008" soybean harvests, the days after planting were the same, and the accumulated incident solar radiation in the period was similar. However, the accumulated degree days were higher in the "Castro 2006/2007" harvest than the precipitation, which was higher in the "Castro 2007/2008" harvest. The "Castro 2006/2007" considering the dry biomass productivity at the V8 stage harvest presented 311.27 kg ha-1 more than the "Castro 2007/2008" harvest (Table 2). The result indicated that higher temperatures in certain phenological stages are more influential in the plants' biomass accumulation than the incident solar radiation.
The "Castro 2007/2008a" and "Ponta Grossa 2007/2008" maize harvests were the ones that had the longest time to reach physiological maturation, being 174 and 175 days after planting, respectively. Even with a difference of only one day, it was observed that the accumulated degree days, accumulated incident solar radiation, and accumulated precipitation were relatively higher in the "Ponta Grossa 2007/2008" harvest. The final dry biomass yield between harvests was also higher, with both surpassing the others (Table 2). At the full flowering stage, the "Castro 2007/2008b" and "Ponta Grossa 2007/2008" harvests presented the same days after planting (85 days). However, in the "Ponta Grossa 2007/2008" harvest, there was an increase of 78.6 °C of accumulated degree days, accumulated incident solar radiation of 169.6 MJ m-2 day-1, and a precipitation decrease of 88.7 mm compared to "Castro 2007/2008b" harvest. The "Castro 2007/2008b" harvest presented a reduction of 1164.74 kg ha-1 in biomass productivity concerning the "Ponta Grossa 2007/2008" harvest, showing the influence of variables on dry biomass productivity.
Even adjusted, the "a" parameters values were similar to dry biomass yield in the final stage of each crop season (Table 3). Each model parameters values were close between harvests, except for "a" parameter in the Brody model. The results obtained are consistent since most of the studied functions are exponential models, in which the "a" parameter represents the maximum value obtained by variable (mass, size, volume, among others) in the considered cycle.
For some variables, the "b" parameter values were close in all soybean crops (Table 3). In the Streibig model, "b" parameters varied between -6.32 and -4.79 for days after planting (DAP), accumulated degree days (GDD), and accumulated incident solar radiation (Rs). For the accumulated precipitation (P), there was a higher variation of "b" (-9.06 to -3.50). The same occurred in Gompertz and Philos, Brody, Bertalanffy, and García y García models, in which the accumulated precipitation (P) in the "Castro 2007/2008" soybean harvest showed higher values in the first collection stage. The "b" parameter values are linked to the shape of the model curves.
(1) Parameters of the equation - "a": kg ha-1; "b": dimensionless; "c": Streibig and García y García models have the same unit as the independent variable, and the others models have the inverse of the same unit as the independent variable. r - Pearson correlation coefficient (dimensionless).
The "c" parameters of the Streibig and García y García models changed between the variables but were close between the soybean harvests (Table 3). The scale of each variable influences the "c" parameter.The Gompertz and Philos, Ratkowsky, Brody, and Bertalanffy models resulted in very minimal "c" parameters values, usually below 0.1.
The Pearson correlation coefficient showed excellent results in all calibrated models, with significant variation (Table 3). Silva et al. (2014) tested Streibig model obtained excellent adjustments (r ≥0.99) between dry weed biomass and the accumulated degree days (GDD). Understanding the production of dry biomass is essential to infer about the physiological processes of plant growth.
The tested soybean cultivar has as characteristic 129 days of the average total cycle. The shortest duration of the total cycle occurred in the "Ponta Grossa 2006/2007" harvest (127 days). However, it was the second among the four crops studied that had the highest dry biomass values at the end of the cycle. The result showed that the plants responded to other factors that are not necessarily related to the time that the crop remained in the field. As observed in the four analyzed harvests, the soybean crop cycle decreased as sowing was done later, as reported by Dallacort et al. (2006) and Kumagai and Takahashi (2020).
Several studies consider that logistic models are the ones that best adjust to the growth and dry matter accumulation by plants (Maia et al., 2009). Many are used to evaluate plant growth and development among the adjusted models, using time as a variable. However, there is little information in the literature about the use of models with the variables as accumulated degree days (GDD), accumulated incident solar radiation (Rs), and accumulated precipitation (P).
The models' adjustment with the accumulated degree days (GDD) variable was adequate due to the thermal requirement being achieved in all the studied harvests. Grein et al. (2015) described for soybean, with 14 °C of basal temperature, that the phenological cycle is completed with 1075 °C GDD. The "Castro 2007/2008" harvest had the lowest accumulated degree days (GDD; 1104.30 °C), showing that the total accumulation of temperature in the cycle was not a limiting factor. However, the temperature influenced the crop cycle in the phenological stages.
As the rainfall variable is not well defined over the crop cycle and the experiment was not irrigated, the lack or excess of rain in a given season may be the cause that most influenced the models' adjustment. During the cycle of all evaluated soybean harvests, there was no water deficiency for long periods. The maximum interval without precipitation was 10 days in the "Castro 2006/2007" harvest; between 550 and 800 mm of water during the soybean cycle are necessary to obtain appropriate yields (EMBRAPA, 2013). In Table 2, it was observed that the "Ponta Grossa 2006/2007" harvest (490.80 mm) was the
only one with accumulated precipitation less than the recommended during the crop cycle. The highest accumulated precipitation occurred in the "Ponta Grossa 2007/2008" harvest (735.2 mm). Under irrigated conditions, Gomes et al. (2014) proposed a model able to predict soybean yield satisfactorily, showing that crops without water deficit resulted in less variability in productivity. Silva et al. (2018) calibrating and validating models to estimate soybean growth had limited productivity due to water stress.
The adjusted parameters for all models in the "Castro 2006/2007" maize harvest differed from the others (Table 3). The dry biomass between the full flowering and physiological maturation stages declined instead of increasing (Table 2). The unidentified problem may have occurred due to some variation or experimental error between the dry matter collection and determination. Therefore, significant adjustments for the "Castro 2006/2007" harvest would not be possible since the tested models consider that the plants' dry biomass only tends to increase and not to decrease.
As well as when adjusting the models' parameters for soybeans, the "a" parameters values of most models were close to the final dry biomass of maize crop (Table 3). In sigmoidal models, the "a" parameters indicate the maximum point reached by the dependent variable. In the Brody and García y García models, there was a high variation between the lowest and highest values. In the Brody model, the lowest and highest values were 3365.5 and 657640.5 kg ha-1, and in the García y García model, the lowest and highest values were 10000.0 and 50000.0 kg ha-1, respectively. In both models, the values occurred in the accumulated precipitation variable. The result indicated that high parameter variations make the models less efficient since the productivity estimate tends to become inaccurate. If precipitation is not well distributed throughout the crop cycle, long dry days may occur, causing a reduction in productivity. In addition, there was a requirement to condition parameters for some models and variables, restricting the maximum limits, since, in the unrestricted calibration process, the "a" parameters values became very high and inconsistent. The maximum limits imposed as a restriction were values slightly higher than the accumulated dry biomass in the final stage of each maize harvest.
The "b" parameter of the Streibig model remained the values close, except for the "Castro 2006/2007" maize harvest (Table 3). In the Gompertz model, the same occurred, but the parameter varied in the adjustment of the accumulated precipitation variable, with the lowest value being 19.16 and the highest 673.06. In the Ratkowsky model, the "b" values were higher than the others, presenting higher amplitude, with the lowest and highest value being 100.0 and 1543.7, respectively. There was also observed a restriction between the lowest and highest value of "b." However, the highest values still stood out in the "Castro 2006/2007" harvest.
In the Streibig and García y García models, the "c" parameters values changed between variables but were close between harvests. The scale of each variable is influenced by "c" values, and in the models Gompertz and Philos, Ratkowsky, Brody and Bertalanffy were very small, all below 0.1.
Campos et al. (2012) obtained a great fit between the days after planting (DAP) variable and total dry biomass accumulation in three types of weed (r ≥ 0.98). Except for the "Castro 2006/2007" maize harvest, the Pearson correlation coefficient obtained in the associations was narrow (Table 3). The lowest value was observed with the Ratkowsky model for maize crop, for the days after planting variable (DAP; r = 0.62). The decay of dry biomass in the final crop stage made it difficult to obtain suitable adjustments.
The maize cultivars with an early growth cycle, such as the one used in the present experiment, have an average cycle between 120 and 130 days. The smallest cycle occurred in the "Castro 2006/2007" harvest (146 days; Table 2), according to the technical characteristics of the cultivar. For each maize phenological stage, caloric units are required for their development. As a result, high temperatures can accelerate the plant cycle due to the daily accumulation degree (Maldaner et al., 2014). The results of the models' adjustment to the daily accumulated degree (GDD) for maize crops, were not satisfactory for the "Castro 2006/2007" harvest. For early cycle cultivars, it takes between 831 and 890 accumulated degree days to achieve pollination. In Table 2, it was observed in the full flowering stage that the lowest value of accumulated degree days was 878.75 °C, not being a factor that limited pollination.
The maize dry biomass production is directly related to solar radiation, converting it into chemical energy in photosynthesis. This process has importance in the plant physiological processes. In years with lower solar incidence and lower temperatures, there is a reduction in crop productivity (Maldaner et al., 2014).
For the maize crop to complete the cycle, between 500 and 800 mm of rain are needed (EMBRAPA, 2004). At the end of the crop cycle, the accumulated precipitation values were within this range (Table 2). During the cycle of all evaluated maize harvest, there was no water deficiency for extended time periods. The maximum interval without precipitation occurrence was 9 days in the "Castro 2006/2007" harvest.
The validation analysis was performed with the ranges of the parameters' values obtained in the calibration for each model and variable.For this purpose, harvests not used in the calibration process were considered. Even with suitable adjustments in the calibration, it was possible to verify that some harvest showed high absolute and relative errors. Table 4 shows the errors, coefficients, indexes, and performance obtained with the tested models in the associations between the observed and estimated dry biomass for the four soybean and maize phenological stages in the "Castro 2007/2008" and "Castro 2007/2008b" harvests, respectively.
For soybean, except for the models' performance for precipitation variable (Table 4), which was between "Terrible" and "Very good," the other variables showed promising performance between "Good" and "Excellent," with a predominance of "Excellent." In the individual adjustment, it was observed in the Brody and Bertalanffy models that the dry biomass values at the beginning of the crop cycle were negative for all tested variables (Figure 1). Therefore, the use of the Brody and Bertalanffy models was restricted to an specific interval.
Among all the models and variables analyzed for maize crop (Table 4), the Brody model obtained a "Terrible" performance for the precipitation variable, and the Ratkowsky model obtained a "Bad" performance for the days after planting variable. The other models performed between "Good" and "Very good." As occurred for soybean crop, in the individual adjustment of the models, it was also observed in the Brody and Bertalanffy models that the dry biomass values at the beginning of the crop cycle were negative for all the variables analyzed (Figure 1). The use of the Brody and Bertalanffy models should be restricted to an specific interval.
Conclusions
Except for the accumulated growing degree days (GDD) variable, the soybean harvests that had the highest days after planting (DAP), incident solar radiation (Rs), and precipitation (P) to reach the physiological maturation stage were not the ones that obtained the highest final dry biomass productivity. The soybean crop cycle decreased as sowing was performed at later dates in Campos Gerais. The association of dry biomass accumulation of the crops and variables tested was satisfactory in the calibration.
With the exception of the models' performance for precipitation variable (between "Terrible" and "Excellent") for soybean and maize, the other variables showed promising performances between "Bad" and "Excellent," with a predominance of "Excellent."
The occurrence of negative values of estimated dry biomass at the beginning of the harvests with Brody and Bertalanffy models indicated a restriction on their use. Even with suitable adjustments in the calibration for soybean and maize, the validation indicated high absolute and relative errors in the harvests.