Parameterization of SSM Model to Analyze Wheat Growth and Yield Potential Under Pakdasht Conditions

Simulation models of crops are used for experimental and complementary research on field projects. These models are also useful for interpreting the results and examining agricultural systems under different environmental and management conditions. The aim of this study was to describe a model for wheat (SSM), guarantee wheat cultivars in a genetic discussion in the Pakdasht environment, and present the results of its evaluation. The model of phenological stages, growth, and aging of leaf area and the production and distribution of dry matter simulates water function and balance. The SSM model simulates the growth stages of the plant in response to environmental factors, heat, and the ability to access solar radiation. In order to evaluate the SSM model, field experiment data of two wheat cultivars—SW and Pishtaz—were used as factorial, based on a randomized complete block design with four replications. Subsequently, the parameters were evaluated, the model was tested in accordance with independent data, and the results indicate its acceptance for the main aspects of crops compared to the observed experiments—for example, for SW, we have 1830 GDD to 2310 GDD from pollination to treatment and extinction factor in Pishtaz is 0.71 and PLAPOW coefficient is 1.6484±.063, which can finally be used to simulate these figures.

The models of the crops have long been used in many countries to simulate the growth of crops. 1,2 because of the fact that the modeling process of crops is a process that requires a lot of time to estimate the parameters and finally run and evaluate the model. Choosing a better model with more reliable results among the various models of growth and development of crops will save a lot of workload in terms of studies.
It is necessary to analyze the factors that limit genetic factors, including genetic, climatic, soil, water, and potential yield of the plant, to study and investigate ways of increasing wheat yield. 3 Such a study and analysis can be carried out using laboratory and field experiments around the world in a few years. As is known, these experiment fields are time-consuming and costly. Today's information technology makes possible to simulate the climatic, biological, and physical phenomena at different stages of the research; these simulation models allow the application of various scientific aspects, including soil chemistry, soil physics, meteorology, crop physiology, and crop improvement by mathematical equations to predict growth and yield. 4,5 These models are good tools for simulations of crops. 6,7 One important point to consider is that these models are complementary to field experiments and are by no means an alternative to field experiments. 8 A basic model for wheat can simulate the phenology, dry matter to various plant components and yield. The SSM Model is used for wheat. 9 With the help of this model, it is possible to analyze the crops and genetic, environmental, and management constraints. 10 The purpose of this study was to evaluate the efficiency and use of SSM model under Pakdasht conditions and determine the genetic coefficients (parameterization) in SW and Pishtaz wheat cultivars.

Materials and Methods
To simulate the development of wheat yield using the SSM model in the province of Tehran, field experiment data were used in 2016-17 (research farm of Tehran University-Abureihan Campus). In this research, two Pishtaz and SW cultivars were used (as used in Tehran province). Before the experiment, depths of 0 to 30 cm were sampled and the chemical and physical properties of soil were determined. Initially, a germination test was performed for the two cultivars and, after a week, SW had 94% and Pishtaz had 85% vigor. According to the densities considered, 50, 100, 150, 200, and 400 plants per square meter of seeds were accurately counted and spaced 20 cm on the cultivation row.
In this research, traits such as the number of leaves, leaf area, number of green leaves left and fall in the main stem, received radiation ratio, phonological stages, dry weight (green and yellow leaf, stem, pod and seed), plant height, yield, and yield components were measured and calculated.

Phenology
All the stages in the development period were carried out using 10 plants and in accordance with the method of. 11 In order to measure dry weight in each sampling (every 10 days), after measuring the leaf area and counting the number of nodes in the main stem, green and yellow leaves, stems, seeds, and shells of pods were separately put in an oven at 70°C to reach constant weight. In order to estimate the radiation received, the radiation meter should be adjusted in accordance with the environment and then the radiation was measured and used in accordance with the method of. 12 An area of around 2 square meters of each plot was harvested once at full treatment stage, where 95% of the wheat was yellow and brown and the seed moisture content was 20%, to determine the yield and yield components. Now the model is named Triticum model was developed by. 13 This model has the ability to simulate the phonological stages, leaf growth and decay, distribution of dry matter, formation of yield, and the balance of soil water. The model also simulates the plant's response to environmental factors like phenoperiod, heat, nitrogen, and radiation on a daily basis and uses the available water-soil-air information. 14 This model does not take into account the effects of diseases, pests, and weeds on the plant. Nevertheless, these factors are not relevant to physiological properties and plant traits.
The quantitative method of predicting phonological growth in wheat was examined as a function of phonological period and temperature in this model.
The key phonological stages required to simulate the physiological processes of growth and yield in wheat are as follows (14) : Planting until germination Germination to tillering Stem to seed Seed to spike Spike to pollination Pollination to physiological treatment Physiological treatment until the harvest The cumulative thermal time required for one stage is expressed in terms of the minimum number of calendar days required for a physiological stage (BD), which is called the required biological day -a number between zero and one which defines a required day for development.
In this equation, f (T) and f (P) are temperature and phonological period functions respectively, which show the relative velocity of growth in each temperature and phonological period relative to the optimal state. The relative velocity response of wheat growth to the average daily temperature uses a" tooth-like function" TMP TBD or TMP TCD In this equation, TMP is the average daily temperature, Tmax is the maximum daily temperature, Tmin is the minimum daily temperature, F(T) is the temperature Function, TBD is the base temperature, TP1D is the final desired temperature, TP2D is the upper optimal temperature, and TCD is the ceiling temperature. In this survey, the minimum and maximum daily temperature, rainfall, and sunny hours were obtained from measured values from the meteorological station of Mehrabad Tehran including semi dried weather, 20 km from the research field.
The model uses a piecewise function to quantify the reaction of the development rate to the phonological period.
Where PP is the phonological period (based on hours per day), CPP is the critical phonological period( i.e.), the phonological period which begins to be reduced because of coincidence with short days of the growth, and PP sen is called the phonological period sensitivity factor. the leaf area The Figure 1 shows the rate of the appearance of philocrone. 15  ...(7) In this equation, PLAi is the current day leaf area and PLAi-1 is the cumulative leaf area of the day before. The coefficient of 10,000 has been added to convert GLAI from centimeters per square meter per day to square meters per square meter per day.
The daily reduction in leaf area index in the model is calculated as follows:

the leaf area index
Finally, leaf area index by increasing GLAI and Daily DLAI in the model can be calculated as follows: LAl t = LAI t-1 + GLAI -DLAI ...(9) dry Matter Production In this model, dry matter production is predicted by a simple RUE-based method.
(Photosynthetic Active Radation) PAR value received per day is obtained from the leaf area index and KPAR. KPAR value was considered constant for the entire growth season. The amount of dry matter produced per day is obtained from PAR multiplication received by RUE. RUE is moderated for the average daily temperature and CO 2 concentration. 16 To calculate the dry matter production per day, calculating the fraction of PAR input received by the crop canopy (FINT) is necessary. This fraction is simply calculated by using an exponential equation between radiation and receipt, which is calculated in accordance with Beyer-Buger-Lambert rule 17 : ...(10) In this equation, KPAR is the light extinction factor.   In the SSM model, a simple method is used for distributing dry matter. In this method, first, the reservoirs for the dry matter and the developmental stages in which the reservoirs are active should be identified. Three reservoirs are considered: (1) leaves, (2) stems, and (3) seeds From the stage of emergence to the beginning seed growth (BSG), the daily dry weight gain of leaves (GLF in grams per square meter per day) is affected by the daily production of dry matter (DDMP in grams per gram / m² / day) and allocation factor to the leaf (FLF in grams per gram). The amount of daily dry production that is not present on the leaves will be allocated to "stems" (GST; stem growth rate in grams per square meter per day):

Formation of yield
The modeling of the seed growth rate and the formation of yield on the basis of the concept of linear increase in the harvest index 9,19 is modified by. 20 The seed growth rate per day (SGR) (g / m² / day) is calculated as follows (13): SGR = DHI × (WTOP i-1 + DDMP i ) + DDMP i ×HI i-1 .... (18) Here, HI is the harvest index, DHI is the increase in daily harvest index in the linear phase of the seed growth, WTOP is the total dry weight (g / m²), and DDMP is the daily dry matter production rate. In the model, the first daily total dry matter produced (DBP) is added to the dry weight of the vegetative components, and then the seed growth rate is deducted from the dry weight of the vegetative components, i.e.:  factor and, in fact, indicates the difference in energy content between vegetative (WVEG) and reproductive tissues (WGRN). The actual value of the gradient of the increase in the harvest index (i.e. DHI) with respect to the maximum value of this gradient (PDHI), the total dry matter at the beginning seed growth (BSGDM) is calculated for the growth season in the simulation. This correction is based on the suggestion of (13). In the model, a tooth-like function was used to calculate the relative gradient of the increase in the harvest index (DHIF). To use this function, the critical quantities of dry matter are required at the beginning seed growth (WDHI). Therefore: ... (22) Multiplying DHIF by PDH I yields the DHI value for the growth season while simulating. During the seed growth, if the daily amount of dry matter production is less than the seed requirement (i.e. SGR), the shortage is obtained from the difference between the modified seed growth rate for dry matter energy (SGR / GCF) and daily dry matter production (DDMP). This is the same as TRANSL.

results and discussion
During the growth period of wheat at the time of performing the test, the average maximum and minimum daily temperatures were 18.6 and 4.9°C respectively. Rainfall during the growing season at the time of performing the test was .96 millimeters . The radiation level at the time of performing the test was not noticeable from November to May.

estimation of ssM model parameters Phonological growth
The values of the parameters necessary to quantitate the reaction to temperature using a toothlike function, such as base temperature, desired lower temperature, desired upper temperature, and ceiling temperature-also known as cardinal temperatures-were considered to be 0, 25, 28, and 40, according to the available references (Table 1) (8). where y is the number of leaves in the main stem, x is time (accumulated growth day after planting), a is the curve contact point with a vertical axis (x = 0), b is linear increase rate in the number of leaves (leaf per unit increase in temperature unit), x0 is the time of the linear increase in the number of leaves, and a + bx0 shows the maximum number of leaves in the main stem. Also, with y = 0 in the equation y = a + bx, the degree of growth day after planting, in which the production of leaves starts in the main stem, is calculated as -a / b (18) . The fitting of Eq. (1) was done separately for cultivars and densities.
95% of cultivars showed a significant difference found between the coefficients of the equation in different densities and cultivars. Fig.  2-4 shows, for example, the relationship between the number of leaves in the main stem and accumulated thermal time. The end times of linear increases in the number of leaves for Pishtaz and SW cultivars were 1496.4 and 1462.2 time units respectively figure 1. The rate of increase of each leaf per unit of increase in the temperature unit for Pishtaz and SW cultivars was 0.0331 ± 0.0009 and 0.0289 ± 0.0014 leaves per unit of temperature respectively.
The amount of philocrone, which is referred to as the time-period between the emergence of healthy leaves in the plant, which is, in fact, unlike that of the leaf emergence rate, in Pishtaz, is equivalent to 32.21°C per day and in SW is equivalent to 34.6°C per day.

PlaPoW estimation
The power Equision. was used to describe the allometric regulation between the number of leaves in the main stem and the leaf area of the plant . 15

Pla = Placon. Msnn plapow
This Equision. was fitted for wheat cultivars. The confidence level of 95% showed that no significant difference was found between wheat cultivars in terms of PLAPOW value. The value of determination coefficient between 0.27 and 0.81 in different densities indicates that a good correlation is found between the leaf area per plant and the leaf number in different densities (Table 3). Figure 2 shows the linear reduction of PLAPOW coefficient relative to the reduction in density in both wheat cultivars. That , PLAPOW is defined as a function of the plant density. In the model for PLAPOW, at a density of 100 plants per square meter in Pishtaz and SW, the values were 1.46 and 1.58 respectively. 21 Also used a power equation (y = x b ) . For describing the leaf area in versus the number of nodes . Contemporary both are increasing in number in the main stem. Using this power equation, were able to compare the leaf area of the plant as a function of the plastochrone index (the interval between the sequential emergences of mung bean leaves. 18 K values are calculated in this study for Pishtaz and SW cultivars .71 and .9 respectively. Figure 3 shows the relationship between the radiation ratio and leaf area index in wheat cultivars (Pishtaz and SW) "R² is the coefficient of determination".
As shown in Fig. 4, the values †of the high determination coefficient indicate an appropriate relationship between accumulated dry matter and accumulated received radiation. 20 RUE value for Pishtaz was 2.36 g / mJ and 2.9 g / mJ for SW. According to a 95% confidence level, no significant difference was found between the two cultivars in terms of the efficiency of radiation use. RUE value for all cultivars was calculated to be 1.37 g / mJ (Fig. 5). Figure 5 shows the relationship between accumulated dry matter and accumulated received PAR in wheat cultivars "Pishtaz and SW and all cultivars".
The base temperature for dry matter production, the desired lower temperature for dry matter production, the upper limit of the desired temperature for dry matter production, and the ceiling temperature for the production of dry matter, based on available sources, are 10, 20, 30, and 40°C respectively (8).

estimation of dry Matter distribution Parameters
The distribution of dry matter means the distribution of photosynthetic materials in different plant organs. 22 To obtain the distribution parameters of dry matter before the end of the leaf growth on wheat cultivars, the data of the leaf dry matter were fitted against the total dry matter (9).
where y is the total dry weight, x is dry matter of the leaf or stem, x0 is the cycle point between two dry matter distribution stages (WTOPL), b1 is the distribution coefficient during the first stage (FLIF1A), and b2 is the distribution coefficient during the second stage (FLF2B). Figure 5 shows the fitting of two-part regression for obtaining FLF1A, FEF1B, and WTOPL coefficients in the cultivars studied for wheat in this study. As can be seen, at lower levels of the total dry matter (WTOPL), a greater part of the dry matter is allocated to the leaves (the first stage), but at higher levels of the total dry matter, lower dry matter is allocated to the leaves (the second stage). The allocation coefficients for the leaf in the first (FLIF1A) and the second (FLIF1B) parts in Pishtaz were .96 and .23 respectively, and, in SW, were .95 and .12 respectively. The cycle point from the first to the second part in Pishtaz and SW was 76.48 and 93.24 grams per square meter respectively. Therefore, when the accumulated dry matter content (WTOP) reaches 76.48 g / m² in Pishtaz and 93.24 g / m² in SW, FLF changes from FLF1A to FLF1B. With regard to the confidence level, a significant difference was found between the densities and cultivars in terms of the value of FLF1A WTPOL coefficient between two cultivars. 15 The leaves have a higher concentration of nitrogen than the stems, and the maximum rate of nitrogen accumulation in Wheat has an upper limit (ceiling).
The values used in the model for FLF1A, FLF2B, and WTOPL were determined for Pishtaz and SW cultivars, based on the data obtained. non-independent evaluation of the model Conventionally, the model is evaluated by comparing the outputs simulated by the model with real-world data gathered. 23 After estimating the plant parameters of SSM model in the previous section, the experiments performed to estimate parameters (Table 4) were simulated, and the output of the model was compared with the observed points (non-independent evaluation). The begining of this section is to determine the extent to which the model predictions are consistent with the measurements. If the predictions based on predefined criteria are in agreement with the measurements, the accuracy of the parameters' estimation is confirmed, and then the model can be independently evaluated. For non-independent evaluation of the model, days to pollination and treatment, leaf area in pollination, number of leaves per main stem, dry weight in pollination, treatment and observed yield were compared with the simulation values by the model. Also, to test the results of the model, the evaluation indices, root mean square error (RMSE), coefficient of variation (CV), correlation coefficient (R), and the rate of deviation of the predicted results from 1:1 and ±18% lines were used. 24 The results of the phonological evaluation of the model showed that the root mean square error (RMSE) for the day to pollination was 61.51 days. The simulated and observed coefficient of variation (CV) was equal to 21.18, which indicates the high accuracy of the model for the prediction of phonological stages. Also, the high correlation between the simulated and observed phenology shows that the model has been successful in predicting the phenology (Fig. 6).
The model was evaluated by using the cultivars Mahdavi, Narin, Pishgam, Zagros, Pishtaz, and Sirvan which were tested in previouse years in the same condition and the results of leaf area (m² / m²) LAI, fraction of available water FASW, accumulated crown dry weight (g / m²) WTOP, accumulated seed dry weight (grams per square meter) WGRN, seed growth rate (grams per square meter per day) SGR, product growth rate (grams per square meter per day) DBP, daily transpiration (mm / day) TR were shown graphically. The cultivar Pishtaz-common to the Pakdasht area-is described as an example.
The minimum temperature is increasing by 0.03 slope per year; b: The maximum temperature is increasing by .05 per year; c: The average temperature is increasing by 0.04 per year; d: The rainfall is reducing by -0.8 annually; e: The maximum leaf area index is 0.007; f: The plant growth per year is 5.9; g: Cumulative dry weight of the grains is increasing by 2.5 per year. the model application Parameterization of Tmin, Tmax, Tmean, WVEG, WGRN, Rain, Maxlai, and DTBSG-DTHAR by SSM model for the analysis of wheat growth and yield potential under Tehran conditions from 1983 to 2015. These parameters show us that how the model works in comparison with another situation. The minimum temperature is increasing by 0.03 slope per year. The maximum temperature is increasing by .05 per year. The average temperature is increasing by 0.04 per year. The rainfall is reducing by -0.8 annually. The maximum leaf area index is 0.007. The plant growth per year is 5.9. Cumulative dry weight of the grains is increasing by 7.1 per yea.

conclusion
In this research, the wheat model (SSM) was used. This model has open source and, because it has a simple structure, can also be used for educational purposes. Since the model uses Excel for input and output, it is easy to work with. This model can be used in many cases, including research applications; crop management and education are examples of this. In this research, the required parameters of SSM model in the phonologic development, production and distribution of dry matter, leaf area growth and decay, yield using field experiment as well as available resources in other researches for wheat, (two cultivars of Pishtaz and SW) were prepared and used in the model. By using them, parameters were estimated, simulation was performed, and the output of the model was compared with the observed points (non-independent evaluation). Non-independent evaluation results, leaf area index, and number of nodes per plant, biomass, and yield showed that the model has an acceptable predictive power for the plant.
According to Tehran's meteorological data, which was used in the model from 1983 to 2015, the results are based on this principle that with the increase in maximum temperatureminimum temperature -average temperature and annual precipitation, the number of days required to reach each stage of the plant phenology has been reduced, which is reasonable.