Building Climate-Resilient Cotton Production System for Changing Climate Scenarios Using the DSSAT Model

: Cotton production is highly vulnerable to climate change, and heat stress is a major constraint in the cotton zone of Punjab, Pakistan. Adaptation is perceived as a critical step to deal with forecasted and unexpected climatic conditions. The objective of this study was to standardize and authenticate a cotton crop model based on climate and crop husbandry data in order to develop an adaptation package for cotton crop production in the wake of climate change. For the study, the data were collected from the cotton-growing areas of Punjab, viz. Bahawalpur and Khanewal. After the calibration and validation against ﬁeld data, the Cropping System Model CSM–CROPGRO– Cotton in the shell of the decision support system for agro-technology transfer (DSSAT) was run with a future climate generated under two representative concentrations pathways (RCPs), viz. RCPs 4.5 and 8.5 with ﬁve global circulation models (GCMs). The whole study showed that a model is an artistic tool for examining the temporal variation in cotton and determining the potential impact of planting dates on crop growth, phenology, and yield. The results showed that the future climate would have drastic effects on cotton production in the project area. Reduction in seed cotton yield (SCY) was 25.7% and 32.2% under RCPs 4.5 and 8.5, respectively. The comparison of ﬁve GCMs showed that a hot/wet climate would be more damaging than other scenarios. The simulations with different production options showed that a 10% and 5% increase in nitrogen and plant population, respectively, compared to the present would be the best strategy in the future. The model further suggested that planting conducted 15 days earlier, combined with the use of water and nitrogen (fertigation), would help to improve yield with 10% less water under the future climate. Overall, the proposed adaptation package would help to recover 33% and 37% of damages in SCY due to the climate change scenarios of RCP 4.5 and 8.5, respectively. Furthermore, the proposed package would also help the farmers increase crop yield by 7.5% over baseline (current) yield.


Introduction
Agriculture is highly dependent on climate and, therefore, climate change could have major effects on crop yields and thus food and fiber supply. It is predicted that the global temperature is expected to increase by 2.9 • C to 5.5 • C by 2060, ultimately causing a considerable reduction in crop production [1]. It is anticipated that the adverse impacts of climate change on the agricultural sector will exacerbate the incidence of rural poverty. The impacts on poverty are likely to be especially severe in developing countries such as Pakistan, where the agricultural sector is an important source of livelihood for a majority of the rural farming community. The impacts of climate change for Pakistan's farmers are likely to be problematic for the food security of many rural households. Pakistan is expected to be one of the countries most affected by climate change in South Asia [2]. As (iii) develop an adaptation package for cotton crops in the wake of climate change using a crop simulation model such as DSSAT.

The Project Area
Two districts of Punjab, Bahawalpur and Khanewal, were selected due to them possessing large areas under cotton cultivation and a high agricultural vulnerability to climate change [10,11]. Field surveys were conducted by the World Wide Fund for Nature (WWF)-Pakistan in the project area of Bahawalpur and Khanewal regarding cotton crop husbandry.

Collection of Field Survey Data
The primary data were collected and assessed in this study. Before the data collection, the field facilitators (FFs) were trained to conduct the preliminary survey. For primary data, the farmers in the studied area were interviewed through well-structured and comprehensive questionnaires to gather the data of the farmers' adoptive practices for cotton cultivation. For secondary data, the Soil Survey of Pakistan (SSOP) and the Pakistan Meteorological Department (PMD) were contacted.

Experimental Site
The experiment was conducted at the agronomic research area of the University of Agriculture, Faisalabad using RCBD with three replications (n = 3) during the year 2019. Two factors were studied in the experiment: the sowing date (22 April,7 May,and 22 May) in main plots and the cotton cultivars (NIAB-878, FH-Lalazar, and BS-15) in the subplot. The size of every trial unit was 4.5 m × 8 m, keeping the row-row distance at 75 cm or the plant-plant distance at 22.5 cm. The experimental site was ploughed with the disc harrow and then with one application of the rotavator to make the soil porous, followed by two cultivations and planking. After this, the seedbed was prepared with the help of a bed-shaper. The seed was sown using the choppa method at a rate of 20 kg ha −1 . After seed sowing, a combination of pre-emergence herbicides, S-metolachlor and pendimethalin (@ 900 mL acre −1 ), were applied to overcome the infestation of weeds. The irrigation scheduling was done according to the recommendations of the Punjab Agriculture Department. After sowing, low irrigation was applied immediately to get better germination. After four days, interval irrigation was conducted to ensure good emergence. Subsequent other irrigations were applied at an interval of fifteen days throughout crop duration. Nitrogen (N), phosphorus (P), and potassium (K) fertilizers were applied at a rate of 150:125:100 kg ha −1 , respectively, whereas nitrogen application was done in three equal split doses. The phosphorus and potassium, along with one split of nitrogen, were applied at the time of sowing while the remaining nitrogen splits were side dressed at the time of square formation and flowering. For all of the treatments, all other agronomic operations such as inter-culture operations and disease and insect pest control were kept regular and uniform. Data on crop phenology, growth, radiation use efficiency, and SCY were collected.

Crop Modeling
The Decision Support System for Agrotechnology Transfer (DSSAT) is a software application program comprised of crop simulation models for over 42 crops (as of Version 4.7.5) as well as tools to facilitate the effective use of the models. The tools include database management programs for soil, weather, crop management and experimental data, utilities, and application programs. The crop simulation models simulate growth, development, and yield as a function of the soil-plant-atmosphere dynamics [13][14][15].
In the current study, the CSM-CROPGRO-Cotton model in the shell of DSSAT [24,25] was used to study the impact of climate change on cotton performance. Furthermore, this model was also used to optimize the production options such as the sowing date, variety selection, plant population, fertilizer rate, and irrigation requirement for developing an adaptation package to enhance climate resilience in cotton in Bahawalpur and Khanewal.
Calibration is a process of adjusting some model parameters to the local conditions. It is also necessary for getting genetic coefficients for new cultivars used in the modeling study. The model was calibrated with the experimental data that included phenology, biomass, LAI, and yield components. The model calibration was started to optimize the most important parameters in soil and genotype files under the non-limiting conditions of irrigation and nitrogen to obtain the crop coefficients potential, which improves the simulations. The soil parameters were adjusted for the best simulation of soil moisture and fertility. It included the soil fertility factor (SLPF), the soil root growth factor (SRGF) (as cotton is a deep-rooted crop), and the drained upper limit (SDUL) of each layer. The SLPF was adjusted before the genetic cultivar coefficients calibration, as it affects the overall crop growth and cottonseed yield and the simulation with a suitable SLPF would be expected to be close to the recorded observation. The simulated cottonseed yield was compared with the observation of all experimental treatments for the adjustment of the SLPF with the lowest root mean square error (RMSE) and a higher index of agreement (d-statistic). Finally, cultivar coefficients were estimated for each variety (Table 1). Cultivar coefficients were determined successively starting from the coefficients dealing with the vegetative growth and phenology of the plant followed by those described by SCY against the best-performing sowing date. To select the most suitable set of coefficients, the iterative approach proposed by Hunt et al. was employed [26].  5  5  5  5  5  5  5  5  5  5  3-9  4  EM-FL  46  45  40  42  43  44  43  44  46  43  35 to 50  38  FL-SH  11  13  12  15  16  14  15  16  12  14  05-20  12  FL-SD  24  25  28  30  26  20  21  26  22  27  10-30  15  SD-PM  48  49  48  50  52  46  45  51  43  45  30-60  42  FL-LF  68  70  72  69  70  69  71  66  65  73  35-80  75 Growth LFMAX  THRSH  65  68  63  70  71  73  69  72  70  72  40-90  , SIZLF: Maximum size of the full leaf (cm 2 ), XFRT: Maximum fraction of daily growth that is partitioned to seed + shell, SFDUR: Seed filling duration for pod cohort at standard growth conditions, PODUR: Time required for cultivar to reach final boll load under optimal conditions, THRSH The maximum ratio of (seed/(seed + shell)) at maturity.
To check the accuracy of the model simulations, it was run with data recorded against the remaining two sowing date treatments. The model was further validated with data collected during field surveys and experiments conducted by WWF-Pakistan in the project area (i.e., Bahawalpur and Khanewal), and the field information regarding crop husbandry was used as input data for the model. The genetic coefficients for other varieties found in the field were also estimated using this field survey data. Simulation performance was evaluated by calculating different statistic indices like RMSE [27] and mean percentage difference (MPD) across all locations. For individual sowing, the date error (%) between simulated and observed grain yield was calculated. The time course simulation of crop biomass and LAI was assessed by d-statistics, an aggregate overall indicator [28]. These measurements were calculated as where P i and O i are the predicted and observed values for studied variables, respectively, and n is the number of observations. A linear regression analysis between the observed and simulated grain yield and biomass at harvest was done to evaluate the performance of the model at different locations. Model performance improved as R 2 and d-index value approached unity, while RMSE, MPD, and error proceeded to zero.

Soil and Weather Data for Model
The Lyallpur soil series was used for Faisalabad during the calibration of the model with experimental data. The data relating to the soil series of each location (Bahawalpur and Khanewal) were collected from the SSOP ( Table 2). The soil profile data of the sand, silt, clay, organic carbon (OC %), pH in water, cation exchange capacity (cmol kg −1 ), and total nitrogen % were collected and used as input for the model. The organic carbon percentage was computed by dividing the value of organic matter by 1.70. The parameters such as bulk density (gcm −3 ), permanent wilting point (%), and field capacity (%) were calculated by the model using the methods described by Rawls et al. [29]. The daily maximum temperature (Tmax), minimum temperature (Tmin), and rainfall and sunshine hours for the cropping season (2019) were also collected (Figure 1a,b). The solar radiations were calculated using the Angstrom formula (Angstrom 1924). These data were also used as the input data set in model simulations.

Statistical Downscaling and Climate Change Projections
Baseline climate data (1989-2019) for daily Tmax, Tmin, solar radiations and rainfall were obtained from the PMD. The climate change scenarios were generated for mid-century (2039-2069) under RCP 8.5. It was assumed that global annual greenhouse gas emissions continue to rise throughout the twenty-first century using the delta method approach. Five General Circulation Models (GCMs) were used for simulating the climate systems in response to climate change factors, e.g., greenhouse gasses (Table 3) and for representing the uncertainty of the prediction of the temperatures and precipitation for all locations from the Coupled Model Inter-comparison Project (CMIP5) [30]. In GCM simulations, the shape parameters of the gamma distribution for a wet event are not of sufficient quality; therefore, mean monthly changes were imposed to the baseline climate using a stretch distribution approach [31]. Parameters such as solar radiation and humidity were assumed to be unchanged during this process. Finally, cool/wet, cool/dry, hot/wet, hot/dry, and middle climate scenarios were developed for the GCMs. Rosenzweig et al. [32] described the detailed protocol. Table 3. Detail of the General Circulation Models (GCMs) used in climate change scenario development.

Climate Change Impact Assessment
A seasonal analysis tool in the shell of DSSAT was used to evaluate the impact of climate change on cotton for mid-century scenarios as predicted by the Pakistan Meteorological Department through the Geophysical Fluid Dynamics Laboratory Climate Model 3 (GFDL-CM3) model. Seasonal files were created from the management practices of optimum sowing dates for both locations. The baseline climate (1989-2019) and future scenarios generated by the GFDL-CM3 model were used to create the weather files for each location. The Carbon dioxide (CO 2 ) concentration of 571 ppm was used for RCP 8.5 under the mid-century, as proposed by IPCC (2013). The impact of climate change was calculated from the future and baseline to the mean yield using the following equation.
The accuracy and precision of the model simulations were further evaluated by running the model with data collected against the remaining treatments (7 May and 22 May). The surveyed data, including the crop husbandry practices of farmers in the project area (Bahawalpur and Khanewal), were used as the input data set, and the model was validated to assess its accuracy at the farmer field.

Climate Change Scenario Generations
Climate change projections for the region were generated using the output of the five GCMs from the latest CMIP5 family under RCP 4.5 and RCP 8.5 scenarios (CO 2 concentration at 571 ppm). A major indication of climate change in the target region complies with the global trend of positive increases in both maximum and minimum temperature. However, there are highly heterogeneous change patterns observed in the precipitation regime owing to the high inter-annual variability in the region. The integrated assessment of the target region is projected with the temperature changes to be highest under GCMs with hot/wet and hot/dry characteristics. For the cotton season, the highest of the changes were projected in a hot/wet climate with a 3.6 • C (3.5 • C) increase in minimum temperature and a 4.3 • C (3.8 • C) increase in maximum temperature (Table 4). Available observed data and generated future data were used for developing the weather files in DSSAT. The simulations were compared with the observed data to study the impact of climate change on cotton productivity in the project area.

Adaptation Strategies for Climate-Resilient Cotton Production
Adaptation strategies were developed by modifying the crop management practices in the model to cope with possible climate change. The values showing the maximum increase in yield under the mid-century scenarios of RCP 8.5 for both locations were considered as the adaptation strategies for climate-resilient cotton production in the project area. The selected adaptation strategies were variety selection, optimization of planting density, change in planting time, optimization of the amount of nitrogen, and irrigational water ( Figure 2). The model was tested against the set of different values/practices and the best combination was selected for maximizing the crop yield under changing climates.

Soil and Weather Data
The geographical coordinates and characteristics of the dominant soil series found in the project area were given in Table 2. This data was used to simulate the cotton yield at the farmer level in the project area using the farmer's crop husbandry practices as the input data set for the model. Likewise, mean monthly changes in the mean Tmax, Tmin, and precipitation were calculated and compared with the future 30-year climate to baseline. After that, monthly changes in the standard deviation (SD) of Tmax, SD of Tmin, and rainy days (>0.1 mm) were calculated and compared with the future climate to baseline (Table 4).

Parameterization and Calibration with Experimental Data
The calibration results in Table 5 revealed that the error percentage between the observed and simulated days to anthesis for NIAB-878, FH-Lalazaa, and BS-15 was 1.67%, as the model predicted. The calibration results regarding physiological maturity with −1.08, −1.08, and −1.08 error percentages were compared to the observed days in NIAB-878, FH-Lalazaar, and BS-15, respectively. In LAI, the results indicated that error percentage differences of −1.92, −3.39, and −5.39 were calculated in NIAB-878, FH-Lalazaar, and BS-15, respectively. In NIAB-878 for TDM, a 260 kg ha −1 difference was calculated with an error percentage of 1.95, and FH-Lalazaar was calibrated with an error percent of 7.08. The difference in TDM for BS-15 was calculated with an error percentage of 5.44, whereas the average percentage difference of 4.82 was calculated in all cultivars collectively for TDM. The observed and simulated SCY was 3159.72 kg ha −1 and 3209 kg ha −1 , 2942.77 kg ha −1 and 3058 kg ha −1 , 2984.16 kg ha −1 , and 3125 kg ha −1 for NIAB-878, FH-Lalazaar, and BS-15, with an average percentage difference of 3.40. These results were reliable enough to evaluate the model against other treatments. The comparison of the time courses TDM ( Figure 3) and LAI (Figure 4) also showed that the model simulates the growth pattern of cotton in a good manner.

Evaluation of the Model with Experimental Data
The accuracy and precision of the model simulations were further evaluated by running the model with data collected against the remaining treatments (7 May and 22 May). The following parameters were compared:

Days to Anthesis
The model evaluation results indicated that the days to anthesis were in the range of good accuracy with a mean error of 9.87%. Table 6 shows the observed and simulated values of days to anthesis for two sowing dates (7 May and 22 May) in three different cotton varieties (NIAB-878, FH-Lalazaar, and BS-15). For NIAB-878, the error difference was 5.36% and 12% on both sowing dates. In the case of FH-Lalazaar, the model evaluated days to anthesis quite well. The error differences in FH-Lalazaar for the anthesis were 7.27% and 14.29% for the second and third sowings. For FH-142, the evaluation results for various sowing dates were also satisfactory. The difference in days to anthesis was calculated at 7.27% in the second and third sowings. The d-index, RMSE and mean percent difference (MPD) were calculated with the values of 0.58, 5.4, and 10.08%, respectively. The deviation in growing climatic conditions for the different cotton varieties in different sowing windows might be responsible for the different phenological (flower initiation) responses.

Days to Physiological Maturity
The model evaluation results indicated that days to physiological maturity were in the range of good accuracy with a mean error of 6.29%. Table 6 shows the observed and simulated values of days to physiological maturity for two sowing dates (7 May and 22 May) in three different cotton varieties (NIAB-878, FH-Lalazaar, and BS-15). For NIAB-878, a 4.05% and 8.70% percent difference was seen after the evaluation for days to physiological maturity after sowing. In the case of FH-Lalazaar, the model evaluated days to physiological maturity after sowing quite well. The error differences in FH-Lalazaar for physiological maturity were 4.05% and 8.70% for the second and third sowings. For BS-15, the evaluation results for various sowing dates were also satisfactory. The difference in days to physiological maturity was calculated at 4.05% and 8.70% in the second and third sowings. The d-index, RMSE and mean percent difference (MPD) were calculated with the values of 0.56, 11.07, and 6.37% respectively. The deviation in growing climatic conditions in the various seedling ages in different fine cotton varieties might be responsible for the different phenological (boll opening) responses.

Leaf Area Index (LAI)
The model fairly simulated the LAI fairly well, with an average error of −6.58%. The data in Table 7 shows the RMSE value to be 0.33, while the value of the D-index among the simulated and observed values was 0.86 and the MPD value was 6.73%. The model evaluation results for NIAB-878 demonstrated that the LAI was simulated with an error difference of −9.42% and 12.97% in the second and third sowings. The error differences of 0.24% and −6.45% were measured in both sowing dates of FH-Lalazar. The evaluation results regarding BS-15 show that the crop model simulated LAI in second and third sowings with estimated error differences of −0.71% and −10.58%.

Seed Cotton Yield
After the calibration, the crop model was evaluated for SCY against two sowing dates (7 May and 22 May) and three cotton genotypes, and the simulated results were quite precise for this agronomic parameter. The model evaluation results regarding SCY for NIAB-878 demonstrated that the model over-simulated with an error difference of 4.69% and 4.26% in the second and third sowing dates, respectively. Furthermore, the mean error difference was 6.28%, collectively, in all sowing dates. The SCY simulation for FH-Lalazaar was also quite satisfactory. Error differences of 7.73% and 9.34% were measured in both of the sowing dates of FH-Lalazaar (Table 7). The evaluation results regarding BS-15 showed that the crop model over simulated SCY with an error percent of 5.65% and 10.02% for the second and third sowings, whereas the mean value of error difference was 6.82% for all the treatments. The statistical indices were calculated. The root mean square (RMSE) value was 183.53 kg ha −1 with a mean percent difference (MPD) of 6.95%, and the d-index value was 0.77 ( Table 7). The comparison of time courses TDM and LAI also showed that the model simulates the growth pattern of cotton in a good manner (Figures 5 and 6).

Total Dry Matter
The above-ground biomass of the two sowing dates (7 May and 22 May) and three cotton cultivars were evaluated, and the simulated results were quite accurate for this agronomic parameter. The model evaluation results for NIAB-878 demonstrated that aboveground biomass was simulated with an error difference of 2.65% and 5.85% in the second and third sowings. The RMSE simulation for FH-Lalazaar was also quite satisfactory. Error differences of 12.92% and 15.81% were measured in both sowing dates of FH-Lalazaar ( Table 7). The evaluation results regarding BS-15 showed that the crop model simulated above-ground biomass in the second and third sowings with estimated error differences of 11.85% and 9.93%. The calculated d-index value was 0.83, the mean percent difference (MPD) value was 9.84%, and the root means square error (RMSE) value was 963.11 kg ha −1 . The d-index value designates a significant relationship between the experiential and simulated values of this factor (Table 7).

Model Validation with Surveyed Data
The surveyed data, including the crop husbandry practices of farmers in the project areas (Bahawalpur and Khanewal), were used as the input data set, and the model was validated to assess its accuracy in the farmer's field (Figure 7). The results showed that the model accounts for the effect of crop husbandry, e.g., the differences in fertilizers, irrigation variety, and weather data for each location. The data presented in Table 8 shows that the model simulated the SCY of different cultivars well, with an average mean percent difference of −0.65% and a RMSE of 27.88 kg ha −1 for the Bahawalpur area (Table 8). Similar results were observed for the Khanewal district area (Table 9)        The results of the calibration and evaluation with an experimental data set and further validation of the model at the farmer's field showed that the CS M-CROPGRO-Cotton in the shell of DSSAT can be used for assessing the impact of climate change on cotton productivity and the optimization of crop production options for current and future climatic conditions.

Impact of Mid Century (2040-2069) Climate Change Scenarios
The current crop area and production in the project are presented in Table 10. The crop husbandry parameter used in the X-build file is given in Table 11. The comparison of baseline and future simulations showed that climate change would reduce cotton yield ( Figure 12). After the calibration and validation of experimental and survey data, the model was run with thirty years of historical climate (1989-2019) data for simulations for the development of baseline and generated data GCM (GFDL-CM3) for the years 2040-2069 under RCP 4.5 and 8.5 for future climate change scenarios (Tables 12 and 13).

Selection of Suitable Variety
A comparison of cotton verities showed that Mubarak and Debal perform better than other varieties at the Khanewal location, while Mubarak and CIM-313 perform better than others at Bahawalpur. However, Mubarak performs better at both location and produce SCY 2907 kg ha −1 . The poorest performance (1244 kg ha −1 ) was observed in Cyto-179 ( Figure 13).

Optimization of Planting Time
The same behavior was shown by different varieties against planting dates. So the better performing variety Mubarak was used in further simulations for the development of the adaptation package. The model was run with seven planting dates (−45, −30, −15, 0, +15, +30, and +45) with an interval of 15 days. The current sowing date was 15 May (0). The model was run with 30 years of climatic data (2040-2069) generated by GCM (GFDL-CM3) for mid-century scenarios under RCP 4.5 and 8.5. The results showed that planting 15 and 30 days earlier would be the best strategy for both locations under RCP 4.5 and 8.5 (Figures 14 and 15).

Optimization of Planting Density
The model was run with three planting densities (46,500, 62,000, and 93,000 plants per hectare) developed by varying plant-to-plant distances (15, 22.5, and 30 cm) in the field. The results showed that under future climate scenarios, 62000 plants per hectare will maximize the SCY during 2040-2069 under both RCP 4.5 and 8.5 in the project area ( Figure 16). A further increase in plant population per unit area will increase the plant-to-plant competition for nutrition, space, light, air, etc. and decrease the yield.

Optimization of Amount of Nitrogen
Nitrogen (N) is an essential nutrient for better crop growth, and our soils are deficient in organic matter and nitrogen. During the early crop season, the cotton zone faces high temperatures and less rainfall. These arid environmental conditions reduce the organic matter and nitrogen in these soils. Keeping in view the importance of nitrogen, the model was run with ten nitrogen levels (25,50,75,100,125,150,175,200,225, and 250 kg ha −1 ) to discover the optimum rate for better cotton production. The results showed that nitrogen application had a linear relationship with SCY until 150 kg N ha −1 . After that addition, the crop yield is not effected (Figure 17). It can be concluded that 150 kg N ha −1 will be the best strategy for maximizing crop yield under future scenarios at both locations.

Adaptation Package for Climate-Resilient Cotton Production
Agricultural production systems are complex, interlinked, and dependent on various factors. Crop production is the climate-prone sector of the economy. The anticipation of and adaptations to climate change are important. There are certain planned and unplanned adaptations regarding climate vulnerability in agricultural systems that maintain the balance in the ecosystem and minimize economic losses. The policies regarding high development must have a synergy effect with climate change for the better adaptive capacity of the nation. To minimize the climate losses, there can be adaptation strategies on farm/regional and national levels. After doing the above sensitivity analysis against different production options, a biophysically viable adaptation package was developed for the climate-resilient cotton production in the area (Table 14). To make the adaptation package more reliable, the model was run with 30 years of climatic data (2040-2069). The results showed that with the adoption of this package, there will be 11% benefits over baseline and 33% recovery of climate change damages under RCP 4.5 projections (Table 15). The model was run with 30 years of future climate data for 2040-2069 under RCP 8.5 for the evaluation adaptation package. The results showed that the package will help the farmers increase crop yield 7.5% over baseline yield and help recover 37% of the damages due to climate change RCP 8.5 projections (Table 16).

Discussion
CROPGRO-Cotton in the shell of DSSAT has been tested by researchers for the growth and yield simulation of crops sown under different climatic conditions with different crop management practices [13,25,[33][34][35]. These studies concluded that CROPGRO-Cotton is capable of estimating climatic impacts on cotton crops. The estimation of climatic impact and possible adaptations for Pakistani cotton with the model are vitally important. The CSM-CROPGRO-Cotton model was well parameterized using the data of field trials conducted at the University of Agriculture, Faisalabad and the surveyed data of farmers' fields at the Bahawalpur and Khanewal districts of Punjab-Pakistan. The Model performed well during the calibration for all studied parameters presented in Table 5. The lowest percent error showed a low variation between the observed and predicted values, which indicates the accuracy of the model during calibration. The slight over-prediction in days to flowering, TDM accumulation, and SCY and the under-prediction in days to physiological maturity and LAI might be due to the poor prediction of seedling emergence at lower soil temperature. The model computed soil temperature from the air temperature for the upper soil layer, while in deeper layers it was computed from the annual air temperature [10]. The model simulated fairly well the dynamics in SCY and TDM for most of the planting dates. Similar performance of the model was also reported by Wajid et al. [36] for SCY with lower MPD values of 5.30% and 4.38% during the calibration and evaluation, respectively, while he was studying the impact of the planting date on cotton performance. Furthermore, an analysis of seasonal simulations exposed that late planting had increased the risk of obtaining very low SCY. This risk was generated by climatic conditions such as lesser solar radiation and variations in temperature, precipitation, and photothermal units when planting was too late and too early (Tables 8 and 9).
The model validation with surveyed data including the crop husbandry practices of farmers in Bahawalpur and Khanewal showed that the model accounts for the effect of crop husbandry, e.g., differences in fertilizers, irrigation variety, and weather data for each location (Figure 6). The comparison of values showed that the model simulated the SCY of different cultivars well, with an average mean percent difference of −0.65%, and 4.9 % at Bahawalpur and Khanewal (Tables 8 and 9). The model simulation was further assessed with the farmer data set collected during 2019. Significant results were recorded with high values of R2 and D-index and low values of error and RMSE for both the Bahawalpur and Khanewal districts (Figures 7 and 8) for the year 2019. These results of calibration and evaluation with an experimental data set and the further validation of the model at farmers' fields showed that the CSM-CROPGRO-Cotton in the shell of DSSAT can be used for assessing the impact of climate change on cotton productivity and for the optimization of crop production options for current and future climatic conditions. Overall, a comparison of the simulated and observed data showed that the model simulated well under the non-limiting condition as compared to under stressed conditions.
To make the adaptation package more reliable, the model was run with 30 years of climatic data (2040-2069). Future climate data for 2040-2069 under RCP 4.5 and 8.5 will be damaging for cotton production due to high temperature and lesser rainfall, as predicted under RCP 4.5 and 4.8 (Tables 12 and 13), which will be a cause of lesser production under the arid environments of Bahawalpur and Khanewal with present crop husbandry practices. However, the results showed that RCP 8.5 will be more damaging than such types, results that were also reported by San José et al. and Abadie et al. [37,38].
The model was run with different production options for developing the adaptations to minimize the adverse effects of the changing climate. The results showed that with the adoption of this package, there will be 11% benefits over baseline and 33% recovery of climate change damages under RCP 4.5 (Table 15). The results displayed that the package will help the farmers increase crop yield 7.5% over baseline yield and help recover 37% of the damages due to climate change, under RCP 8.5 projections (Table 16). An early sowing date (15 days earlier than present) will be the best strategy along with other proposed options for the future, as late planting pushes the crop into a lower radiation time of the year. The higher temperatures at early vegetative growth phases led to earlier flowering, boll set, boll opening and the lower growth of canopy (LAI) these results are in line with Wajid et al. [36], who applied CSM-CROPGRO-Cotton for the optimization of different cotton management practices under the semi-arid climatic conditions in central Punjab.
These benefits of the adaptation package will be more if plant breeders help to evolve more heat and drought-tolerant varieties than the present-day varieties such as Mubarak. The breed is suggesting an emphasis on fast-growing varieties with lesser LAI and no stomata on leaves to avoid excessive transpiration losses.

Conclusions
The current study revealed that that the CSM-CROPGRO-Cotton model in the shell of DSSAT V. 4.7 is an artistic tool for the simulation of crop growth, phenology, and yield. A close agreement between the simulated and observed data is proof of these qualities. Furthermore, a well-calibrated CSMCROPGRO-Cotton and validated model could be used for climate change impact assessment and the optimization of site-specific production technology with the use of minimum resources. The proposed adaptation package could be proven effective in quantifying the climate change impact on SCY in the project area and in developing possible adaptation strategies for sustainable cotton production that can mitigate the harsh effects of climate change.