Determination of Upland Rice Cultivar Coefficient Specific Parameters for DSSAT (Version 4.7)-CERES-Rice Crop Simulation Model and Evaluation of the Crop Model under Different Temperature Treatments conditions ()
1. Introduction
Crop simulation models have been utilized globally as an effective or planned research for decision support tools in crop productivity or resources management. Several crop models used from long time to assist crop management practices with exploring physiological processes under different environments [1]. The Decision Support Systems for Agrotechnology Transfer (DSSAT) is a popular package software comprises crop simulation models (CSM) for over 42 different crops including rice. The DSSAT software consists of: 1) data base management system for soil, weather, genetic coefficients, and management inputs; 2) crop simulation models(CSM); 3) series of utility programs; 4) series of weather generation programs; 5) strategy evaluation program to evaluate options including choice of variety, planting date, plant population density, row spacing, soil type, irrigation, fertilizer application, initial conditions on yields, water stress in the vegetative or reproductive stages of development, and net returns. DSSAT-CSM models simulate growth, development, and yield of crops as a function of the soil-plant-atmosphere-management dynamics [2].
Crop modelling can also be useful to understand the scientist, researchers define research priorities. Using a model to estimate the importance and the effect of certain parameters, a researcher can observe which factors should be more studied in future research, thus increasing the understanding of the system. Crop Simulation Models (CSM) are tools of systems that define procedures of crop growth and development as an act of weather, soil conditions, crop management and help in solving problems related to crop production [3]. CSM was used to simulate grain yield, biomass, and water balance in rice crop [4]. For dynamic crop simulation model accurate estimation of crop cultivar coefficients is the main point into use (for research as well as decision making) for improvement, identification and consequently narrowing gaps in our knowledge over crops and biophysical aspects for improved agricultural productivity.
Rice (Oryza sativa L.) rice is one of the most important crops and it represents a staple food for over half of the world’s population, with a global production of more than 700 million tons per year and a harvested area reaching 165 million ha [5]. In 2050 world’s population growth will be 10 billion and the demand for rice will grow faster than other crops [6]. There are already many challenges to attaining higher productivity of rice. Climate change e.g., high temperature and its consequences are the major challenges which affect the production of rice crop.
Upland rice suffers severely from irregular environmental factors, e.g., air temperature, drought, and precipitation [7]. Temperature is attributed by its impact on crop yield, due to expansion under excess heat stress conditions that greatly influences the growth duration and yield of the rice plant [8]. The temperature rise is one of issues of climate change that has the effect of rice production in Thailand, especially to the development and growing of upland rice plants reported by [9]. In the growing season the mean temperature, temperature sum, ranges, distribution pattern and diurnal changes or a combination of these, highly correlated with grain yields had a significant issue [10]. Due to global warming and climatic risk, the current rice productions in Thailand are in danger. To fulfil the increase rice demand of ever-growing population pressures, an estimation of likely impact is vital for planning strategies.
For estimating cultivar coefficients, numerous methods have been recognized. However, these methods need key knowledge regarding a specific crop cultivar such as planting dates, flowering dates, physiological maturity dates and final grain yield, which in most cases are not available. Genetic coefficient calculator (GENECALC) which is a sub module in the Decision Support System for Agrotechnology Transfer (DSSAT v4.7) was used to determine cultivar coefficients for new peanut lines in Thailand from standard varietal trials reported by [11]. Some researchers used generalized possibility uncertainty estimation (GLUE) method to estimate maize cultivar coefficients [12]. Hence, DSSAT v4.7 has GLUE module for estimating crop cultivar coefficients [13]. All the above-mentioned approaches to estimate crop cultivar coefficients for use in dynamic crop models need some degree of information on a particular cultivar. Therefore, in situations where there is scarcity of data from typical variety trials or other devoted experiments, repeated field experimentations would be the only option. For dynamic crop simulation model accurate estimation of crop cultivar coefficients is the main point into use (for research as well as decision making) for improvement, identification and consequently narrowing gaps in our knowledge over crops and biophysical aspects for improved agricultural productivity. Therefore, objectives of current study are: 1) to determine upland rice crop growth and development indices under optimum temperature conditions; 2) to estimate upland rice cultivar parameters and calibrate DSSAT-CERES-Rice model using the same management; 3) to evaluate DSSAT-CERES-Rice model for simulating growth and yield under the rainfed upland conditions.
2. Material and Methods
2.1. Field Experiment for Model Calibration
For model calibration, five popular upland rice genotypes namely, Dawk Pa-yawm, Mai Tahk, Bow Leb Nahng, Dawk Kha 50 and Dawk Kahm were grown on 8th July 2018 to calibrate the CSM-CERES-Rice model by using a temperature tolerant cultivar (Basmati 385). The experiment was conducted as a split plot using Complete Block Design (CRD) as main plots with three replications. The main plots had five upland rice genotypes and sub plots had three different temperature treatments (28˚C, 30˚C and 32˚C). Temperature treatments were set up by using the data logger (UA–003–09 HOBO Pendant, for Temp/Light) from starting of the experiment. Each plot had 5 rows, with 10 plants each. Standard agronomic practices were followed for treatment of fertilizers, weed and insect control. The crop management data (i.e., phenological data) required for the simulation of the crop model include planting date, (50% and 100%) germination date, flowering and maturity, tillers number, panicles number, leaf area index, grain yield, biomass and 1000 grain weight were recorded harvest.
2.2. Description of the DSSAT-CERES-Rice Model
CERES-Rice present in DSSAT v4.7which is an advanced physiologically based model, was used to calibrate and evaluate the crop simulation model. Genetic co-efficients for the five upland rice genotypes were used to calibrate CERES-Rice model. Soil analysis was carried out before started of the experiments to analyze soil fertility and to carryout proper fertilizer management. Weather parameters including maximum and minimum temperatures, rainfall with air intensity were recorded from Kho Hong Agro meteorological station. Daily solar radiation (MJm−2 day−1) was calculated by using weatherman tools in DSSAT v4.7. Dawk Kahm genotype was used as border crop to avoid the varietal errors. To run the model the following five input files were created:
2.2.1. Daily Weather Model
Maximum and minimum air temperatures, precipitation, rainfall, and solar radiation (derived from sunshine hour data) were collected from the weather station of Kho Hong Agro meteorological office, Hat Yai, Thailand.
2.2.2. Soil Data
Soil data were collected at the depth of soil characteristics at 0 - 20, 20 - 40, 40 - 60, and 60 - 80 cm before planting. Soil classes, organic carbon (%), sand, silt, clay (%), soil texture, soil pH in water, field capacity (%), organic carbon (%), cation exchange capacity, total nitrogen, potassium and phosphorus, potential root distribution and depth were taken (Table 1).
2.2.3. Management Practices
Plant density, planting date, irrigation, weeding, plant row spacing, sowing depth, amount and types of fertilizers, insecticide application was done whenever necessary.
2.2.4. Plant Profile Data
Sowing date, emergence date, flowering date, physiological maturity date, panicle initiation (when 50% and 100% of the crop had reached those stages), planting density, plant height, tops weight (grain weight), harvest index and grain yield per genotype, i.e., grain yield per area of production.
Table 1. Soil Physical, chemical and morphological analysis for model calibration.
2.2.5. Genetic Coefficients File
Table 2 showed the genetic coefficients that were determined in the CERES-Rice model with the parameters namely P1 (Time period or basic vegetative phase), P2O (Critical photoperiod), P2R (Photoperiodism coefficients), P5 (Grain filling duration coefficient), G1 (Spikelet number coefficient), G2 (Single grain weight), G3 (Tillering coefficients) and G4 (Temperature tolerance coefficient). The genetic coefficients or cultivar coefficient values as obtained through GLUE runs were replaced by calculated values against tested genotypes and saved it in RICER047.CUL file in DSSAT model. Model calibration was done several times for P2O, P2R, G3 and G4 and subsequently suitable values for the coefficients were selected.
2.3. Experiments for Model Evaluation
For model evaluation, selected five upland rice genotypes namely Dawk Pa-yawm, Mai Tahk, Bow Leb Nahng, Dawk Kha 50 and Dawk Kahm and three temperature treatments were laid out in a Completely Randomized Design (CRD) under split plot design structure with three replications in the 2017/2018 and 2018/2019 growing seasons at the research area of Plant Science Department, faculty of Natural Resources, Prince of Songkla University Hat Yai campus, Thailand. Soil samples were collected fifteen days before sowing for important chemical and physical characterization (Table 3). We must do two seasons research experimental work for model evaluation. Here planting was done on 24th
Table 2. Genetic cultivar coefficients for the DSSAT-CERES-Rice model.
Table 3. Chemical properties for model evaluation experiments.
July for the 2017/2018 and 2018/2019 seasons. Different fertilizers whenever necessary. Other management practices were carried out accordingly. The number of days to anthesis and physiological maturity, and grain filling data were collected. Moreover, grain yield and total plant biomass was measured at physiological maturity.
2.4. Statistical Analysis
The analysis of variance (ANOVA) to evaluate cultivars growth and development and the effect of different temperatures was done. Mean separation was done by Least Square Difference (LSD) for split plot design to see the varietal differences. Test of significance between the 2017/2018 and 2018/2019 experiments and simulated and measured quantities was performed by R/agricolae program [14]. This research focused to simulate the effect of different temperatures on yield performed for phenology, grain yield, tops weight, leaf area index, harvest index, and tillers number.
2.5. Model Calibration
Calibration is the process of adjusting some model parameters to local environmental conditions and obtains genetic coefficients for new cultivar used in modeling study [15]. According to [16] model calibration was done by the Root Mean Square Error (RMSE).
(1)
2.6. Model Evaluation
Model performance was evaluated by comparing the simulated versus observed values from upland rice experiment under rainfed conditions where an agreement index or d stat index [16] was used. The RMSEn gives the level of error associated with each evaluation between the observed and simulated values.
(2)
(3)
where Si = simulated value, Oi = observed value, n = number of observations,
= overall mean of observed values,
and
.
3. Results and Discussion
3.1. Weather Conditions
From Figure 1 and Figure 2 it was showed that in both 2017 and 2018 years, weather conditions had no visible variations. The maximum and minimum
Figure 1. Mean daily maximum, minimum temperature, rainfall, solar radiation, sunshine determined in 2017/2018 year.
Figure 2. Mean daily maximum, minimum temperature, rainfall, solar radiation, sunshine determined in 2018/2019 year.
temperatures ranged between 25˚C to 38˚C and 18˚C to 22˚C. Highest maximum temperature (35˚C) and highest minimum temperature (26˚C) shown at ambient temperature 28˚C and 32˚C for 2018 Solar radiation ranged from 10.2 to 20 MJ∙m−2∙d−1, rainfall was reported between 20 - 980 mm and the highest rainfall was shown in the month of November to mid-December. After that January month was second most rainfall month near about 850 mm rainfall, respectively. Extreme high temperature beyond the average temperature hampered during experiment might have negative influence on upland rice yield, crop growth and development.
3.2. Soil Condition
Different soil layers of the study site with surroundings contain sandy clay loam texture and were highly acidic. Special management e.g., irrigation due to higher clay contents and drainage system should be adopted for heavy rainfall. As the soil was acidic, phosphorus fixation was required to improve soil pH, phosphorus fertilizers were added.
3.3. Analysis of Variance Result of Yield and Yield Attributes
Although ANOVA results (Table 4) for yield contributing traits using least significant differences (LSD) test (P < 0.05 and P < 0.01) revealed that some traits showed highly significant differences with temperatures and genotypes. Grain yield, tops weight, harvest index, flowering date, and maturity date was significantly affected by temperatures. Non-significant difference for tillers number, and leaf area index occurred possibly due to the optimum input of temperatures at the early stage. Whereas all the phenological traits had shown highly significant
Table 4. Mean squares from analysis of variance for 5 upland rice genotypes.
Table 5. F values for selected variables between upland rice cultivars for 2017/2018 and 2018/2019 seasons.
variation among different genotypes. Grain yield, flowering date and maturity date had interaction effects between temperature and genotype.
3.4. Harvest Index, Grain Yield and Tops Weight
Table 5 showed that both grain yield at harvesting period and tops weight at physiological maturity varied significantly (P < 0.05) during 2017/2018 and 2018/2019 seasons, except harvest index at 50% flowering did not fluctuate. So, there was significant variation among genotypes in the 2018 and 2019 seasons. Also, there was inter seasonal variation in the tested variables for tops weight and grain yield. The similarity in the parameters for the 2018 and 2019 growing seasons suggested that at the growing condition period rainfed condition, nutrients, fertilizers, and other managements were affected the yield contributing traits. Though, sowing date for both seasons was the same 24th July, and weather elements were comparatively similar in both years, finally the similarity could not be expected.
3.5. CERES-Rice Model Calibration
Experiment for estimation cultivar coefficient was used for model calibration. Data of observed and simulated days to anthesis, physiological maturity, grain yield, by product and tops weight were collected for all the genotypes (Table 6(a) and Table 6(b)) indicating that cultivar specific parameters within the model were reasonably adjusted. Similarly, there was a good relationship among observed and simulated variables (d state values were respectively 0.97, 0.98, 1.0,
(a) (b)
Table 6. (a) Simulated and observed values of days to anthesis, days to maturity and grain yield for upland rice genotypes. (b) Simulated and observed values of by product yield and tops weight for upland rice genotypes.
1.0 and 0.99 for all the yield attributes. Dawk Pa-yawm showed high RMSE with respect to both number of anthesis days and physiological maturity than other genotypes. Mai Tahk had higher RMSE with respect to grain yield and tops weight also indicated that was high yielding genotypes and the yield is associated with the growth duration since it took longer than others to attain anthesis and physiological maturity. Besides, Bow Leb Nahng had lower RMSE with anthesis days, maturity days, grain yield, by product and tops weight indicated lower yielding genotype.
3.6. Cultivar Specific Parameters
Table 7 showed that Dawk Kahm genotype required few thermal units to complete juvenile stage (P1) compared to Mai Tahk genotype. This allowed Mai Tahk required extra time to accumulate photosynthesis before flowering and hence higher yielding genotype compared to other cultivars. Dawk Kahm was formerly bred for temperature tolerance condition indicated that few temperature units of short duration were just required to attain end of juvenile stage. Dawk Kha 50 required few temperature units from anthesis to physiological maturity (P5), contrasting Mai Tahk with highest thermal time requirement (G4). Spikelet number coefficient (G1) and single grain weight (G2) also high in Mai Tahk cultivar and lower in Dawk kahm cultivar. This resembled to differences in grain size between the two cultivars. The rate of grain development was high in Mai Tahk compared to others since this could be a temperature tolerant cultivar for rainfed environment with upland condition. Temperature tolerant coefficient (G5) was for the cultivars varied from 0.78˚C-1.09˚C for Dawk kahm to Mai Tahk genotypes.
3.7. Model Evaluation
The cultivar specific parameters found from experiments were used to evaluate CSM-CERES-Rice for simulating different temperature treatments under rainfed upland condition. The model simulated well the average number of anthesis to maturity days with high degree of agreement as indicated by the agreement index (d-statistics) (Table 8). This is an indication that the model calibration and resulting cultivar specific parameters were reasonably estimated. Generally, there were significant different (p < 0.05) between observed and simulated data at all temperature treatments and in all variables. Especially simulated yield decreased as temperature increased in both model simulation and experimental observations. This suggests that the CERES-Rice model is sensitive to climatic variables such as temperature. Furthermore t-test showed significant difference (P < 0.05) between simulated and observed yields at all temperature treatments (Table 9).
For model evaluation, a comparison was made for the five upland rice genotypes between simulated and observed grain yield at 28˚C, 30˚C and 32˚C temperature treatments (Table 7). In this study all the genotypes had shown a close agreement between observed and simulated values of grain yield data. Results showed that Mai Tahk and Bow Leb Nahng were the best cultivars due to
Table 7. Cultivar coefficients for the five upland rice genotypes.
Table 8. Observed and simulated values of yield attributes as affected by three temperature treatments under rainfed conditions.
Table 9. Observed vs. simulated grain yields at different temperature treatments.
increasing simulated grain yield and best performed with the moderate temperature at 30˚C (Mai Tahk 5500 kg∙ha−1 and Bow Leb Nahng 5000 kg∙ha−1). Whereas at 32˚C temperature treatment, a large variation for all the genotypes decreased in yield with the maximum temperature increased. Highest reduction in simulated yield was recorded for Dawk Kha 50 cultivar (3650 kg∙ha−1 at 28˚C, 4220 kg∙ha−1 at 30˚C and 3055 kg∙ha−1 at 32˚C) compared to others while minimum temperature increased at high temperature (30˚C) stage.
4. Conclusion
Unavailability of experimental data related upland rice genotypes have hampered DSSAT-CSM use in many areas and research such as biophysical resource application, economical, cost effective tactical or strategic decision-making processes. Here, simulation experiments done using these cultivar parameters predicted reasonably well the crop growth and yield under varying temperature treatments. It showed that highest simulated grain yield bearing genotype was Mai Tahk 5500 kg/ha with 30˚C and the simulated yield of upland rice ranges (3055 - 5500) kg/ha with temperature ranging 28˚C - 32˚C. The CSM-CERES-Rice model was well validated and showed reliability of simulations with different temperatures for phenology and yield attributes of upland rice genotypes. Our study assessed the DSSAT–CERES–Rice model in a CRD experiment including two factors; one of genotype (5 Thai upland rice genotypes) and another one of different temperature treatments (3 levels) to identify the variations of grain yield in southern Thailand under changing climatic (weather, soil) conditions. Study results revealed that increase in temperature, solely lead to decrease in grain yield of upland rice. Thus, CERES-Rice model could be safely used as a tool to assess different agronomic and climatic change parameters.
Acknowledgments
The authors gratefully acknowledged the financial support from the budget (NAT610046S) of Prince of Songkla University Higher Education Research promotion and Thailand’s Education Hub for ASEAN countries (TEH-AC) project office of the Higher Education Commission. It was partially supported by the Graduate School, and thankful the Field Laboratory Centre of Plant Science Department, Faculty of Natural Resources, Prince of Songkla University, Hat Yai campus, Thailand, for kindly providing the upland rice seeds.