Wheat (Triticum aestivum L.) is a primary staple crop worldwide, with a projected global production of ~740 million metric tons in 2017  . According to the 2017 estimates of the United States Department of Agriculture (USDA), the U.S. ranks fourth in wheat production by country with a projected production of 49.64 million metric tons. The southern Great Plains of the U.S. (Kansas, Oklahoma, and Texas) accounts for ~30% of the U.S. wheat production, producing 18.8 million metric tons of winter wheat per year from an area of 7.5 million hectares  . Achieving sustainable and equitable food security and profitability of crop production in future relies on better understanding of climate and changes in greenhouse gas concentration  . In North America, wheat production is expected to decrease or cease in the southern latitudes and geographic shifts in production in the northern latitudes are more likely to occur  .
Atmospheric CO2 concentrations recorded at Mauna Loa, Hawaii show an increasing trend for CO2. Initially measured at 313 ppm, now at the time of writing this paper the atmospheric concentration has crossed 408 ppm (https://www.esrl.noaa.gov/gmd/ccgg/trends/). Increased emissions of anthropogenic greenhouse gases over the last few decades have exacerbated the situation causing increase in temperature and more frequent and severe weather events over larger spatial domain  . Climate variability and extremes have multifold consequences and damage the economy as well as natural systems, and could result in pronounced deleterious impacts on food security in less developed regions of the world.
General Circulation Models (GCMs) encompass state of the art knowledge and utilize the most advanced mathematical models currently available representing physical processes of the planet. They are used to understand how the global climate dynamics responds to increasing greenhouse gas concentrations. The GCMs describe the climate on a three-dimensional grid, with a resolution anywhere from 250 - 600 km. However, there is a mismatch between the temporal and spatial footprint of the GCM data and data required to run process-based crop models  . Therefore, in order to reliably assess climate change impact at finer resolution, downscaling is required on relevant meteorological variables. Several downscaling approaches are available, which can be broadly classified into statistical downscaling and dynamic downscaling. Dynamic downscaling method nests regional climate model (RCM) into the GCM to represent a given boundary forcing but are computationally very expensive. Statistical downscaling methods are seen as alternative to dynamic downscaling methods for they use empirical relationships established between large-scale and fine-scale variables using historical data  .
Several uncertainties limit the predictability of future climate. Therefore, for a better understanding of climate trajectory in future, different sets of Global Circulation Models and Representative Concentration Pathways are explored to obtain a wider range of conditions while assessing the climatic conditions. Inter model differences of GCMs are due to the way certain physical processes, responses, and feedbacks are modelled. The Representative Concentration Pathways (RCPs) are four different emission scenarios that provide time-dependent projections of atmospheric greenhouse gas concentrations. Each pathway represents a range of radiative forcing and emission characteristics, for e.g., RCP2.6 comprises a stringent mitigation scenario, RCP4.5 and RCP6.0 comprise intermediate scenarios, and RCP 8.5 represents a very high greenhouse gas emission scenario (Figure 1). Current and past anthropogenic greenhouse gas emissions along with natural climatic variability are likely to cause an increase in global surface temperature by 3˚C - 6˚C for RCP6.0 and RCP8.5 by 2100  .
Several studies have shown that rising temperatures affect crop production negatively, due to shortened growing season and decreased photosynthate accumulation in plants   . In addition, a meta-analytic summary  indicated CO2 fertilization effect to be lower than expected on C3 crops and little to no effect on C4 crops. In case of wheat, exposure to higher temperature extremes has been linked to accelerated growth, shorter growing season, and reduced yield  . Increase in ambient temperature singly or in combination of drought during reproductive and grain-filling phase in wheat, decrease grain yield and the protein content   . Since 1980, warming in Bangladesh and India is considered to have reduced wheat yields by approximately 20% of their average trend  . In the Southern Great Plains of the U.S., Oklahoma and Texas have lower and more variable wheat yields than in the rest of the U.S.  . Agricultural production in this region has always been challenged by climate. The situation gets exacerbated due to frequent drought as evapotranspiration exceeds the amount of rainfall. In addition, over the years the ratio of harvest-to-planted ratios in this region is also declining as the wheat fields get abandoned after being used as pasture because part of the crop used as pasture turn out to be more profitable than growing the crop for grain  . Since wheat is vital for the southern Great Plains of the U.S. for its impact on the local economy and food
Figure 1. Atmospheric CO2 concentration for RCP6.0 and RCP8.5, source: (Meinshausen et al., 2011).
security, it is imperative that we understand the effect of changing climate on wheat yield.
Crop simulation models may prove to be a better tool in linking bio-physical drivers of crop production with climate variability. Crop modeling offers a big advantage as various treatments can be investigated in greater detail without actually having to run the experiment on the field, which can be both time and cost prohibitive. Process-based crop simulation models, e.g.     are the suitable tools to examine impacts of climate change because of their ability to simulate and mechanistically describe the detailed influences of a wide array of weather variables, plant genotypes, environmental and edaphic factors, management practices, etc. on crop physiological and growth processes. Therefore, these crop simulation models can be very useful in understanding complex and non-linear effects of climate change on crop growth and physiological response.
Previous studies on climate change impacts on wheat production in Oklahoma   have been carried out using the IPCC Third Assessment Report Special Report on Emissions Scenarios data of 2000. No literature has been published yet for climate change impact assessment on wheat yield performance using the downscaled AR5 data for this region. Therefore, the objectives of this study were i) to simulate wheat yield under different future climate anomalies at a finer scale over wheat growing areas in Oklahoma using the latest CMIP5 GCM climate projection, and ii) to investigate how wheat yield responds to different climate scenarios for different emission pathways across landscapes of Oklahoma.
2. Materials and Methods
2.1. Study Area
The study area is Mesonet locations in Oklahoma (Figure 2). Wheat yield was simulated with Decision Support System for Agrotechnology Transfer (DSSAT)-Crop Environment Resource Synthesis (CERES)-Wheat v4.6  water and nitrogen limitations (Ywf), water limitation (Ya), and no limitations (Yp). The sowing date of October 16 (DOY 289) was used for entire simulation. A very popular hard red winter wheat cultivar Duster, developed by the Oklahoma State University Wheat Improvement Team  was used in the simulations. Duster performs well in both grain-only and dual-purpose (grazing and grain or hay production). A total of 77 active Mesonet sites were chosen with least missing weather data. The time slices of 1980-2014 and 2040-2060 were selected to represent the baseline and future climate respectively.
2.2. Crop Simulation Model
For wheat yield simulation, we used the well-established DSSAT-CERES-Wheat model    . DSSAT is a mechanistic simulation model that uses daily weather data, geographical coordinates, sowing date, fertilizer applications and soil properties to predict the growth, development, physiology and yield of crop
Figure 2. Map of study area and geographical locations of the Mesonet sites within Oklahoma used in winter wheat yield simulation.
of our interest. The DSSAT CERES-Wheat is a preferred crop model in simulating winter wheat growth at field scale in semi-arid condition  . This process-based model is able to account for the effects of climate, crop genetics, soil, and agronomic practices on crop phenology and growth.
Same management options and inputs used for baseline simulations were used for future yield projections. Simulations for winter wheat were carried out using the sequential simulation option and adopting the continuous winter wheat-summer fallow field management option. Wheat was planted in October and harvested at physiological maturity, and was immediately followed by summer fallow. The fallow period was ended a day before the planting of the crop.
2.3. Soil Data
The DSSAT CERES-Wheat requires data on soil physical and chemical properties such as soil color, horizon depth, wilting point, field capacity, bulk density, organic carbon content etc. Information on soil properties for the 77 locations was acquired from the gridded Soil Survey Geographic Database (gSSURGO). For each Mesonet station location, dominant soil type for the county was chosen as representative of the soil type. Information on several soil physical and chemical properties that serve as inputs for CERES-wheat were extracted from the gSSURGO data. The data included upper and lower depths for each horizon (cm); sand, silt, and clay percentage; bulk density (0.33 bar); pH in water. Soil profile data for each horizon was created. The Quad UI tool (http://tools.agmip.org/quadui.php) from the Agricultural Model Intercomparison and Improvement Project (AgMIP), a desktop application for windows system, was used to convert gSSURGO soil database for each site to DSSAT-ready soil data format.
2.4. Climate Data
2.4.1. Baseline Climate Data
Robust crop yield prediction can be performed when the crop simulation is driven by long term (>30 years) daily weather data. Crop modeling system to analyze crop yields and examine future climatic impacts on a regional level requires detailed weather data on a daily basis on a higher temporal and/or spatial resolution. For baseline weather data, two data sources were used, Daymet  and the Oklahoma Mesonet  . The Daymet weather data were downloaded in the form of tarred netCDF tiles for the state of Oklahoma. Values for weather variables, minimum and maximum temperature, rainfall, and incident solar radiation, and day-length from 1980-2014 were extracted at the Mesonet site locations. Information on day length was acquired to convert the incident shortwave radiation flux density (Wm−2). Since the establishment of the Oklahoma Mesonet was initiated during the mid-90s, the Mesonet does not meet the requirements for creating a complete baseline data (>30 years). Therefore, the missing data for each of the Mesonet sites was acquired from its closest National Weather Service (NWS) Cooperative Observer Program (Coop) weather station. As the NWS Coop stations do not measure solar radiation, solar radiation was estimated using RadEst 3.0  . Climate summary files were generated using a python script that produces climate summaries for the 77 sites for both Daymet and Mesonet data. Weather data were processed and adapted to the requirements of the weather input files (*.WTH) for the DSSAT model.
2.4.2. Generating Future Daily Weather Data
Regional Climate models have limited applicability because of the computational time involved for properly assigning boundary conditions to a finer-scale weather phenomena  . One way to address the issue of inherent uncertainties within the GCMs is to create multi-model ensembles as they overcome the limitations of any single model  . Multi-model ensemble averages are found to outperform individual models when tested for the historical period  . In addition, they also help distinguish the variability between the models and the uncertainties of future climate  . Outputs of 4 GCMs, their ensemble, and ensemble of 17 GCMs were linked to DSSAT to account for the uncertainty of climate projections. Detailed information on forcing data can be found on the CMIP5 website (http://cmip-pcmdi.llnl.gov/cmip5/availability.html).
For future climate data, we used data from four state-of-the-art GCMs, HadGEM2-ES, MRI-CGCM3, MRIOC5, and CSIRO-Mk3.6.0 that were included in the Coupled Model Intercomparison Project Phase 5 (CMIP5)  (Table 1). The four models were specifically chosen as they had the finest spatial resolution. In addition to the four different GCMs, an ensemble of the four GCMs along with an ensemble of 17 different GCMs were also included to generate climate trajectories from the 5th Assessment Report  . Two emission scenarios reported in the Special Report on Emissions scenarios by IPCC were used for downscaling the data for the Mesonet locations. A moderately high (RCP6.0) and high emission scenarios (RCP8.5) were chosen to create weather outputs for 2040 to 2060 for each of the target 77 sites. Choosing multiple RCPs and GCMs allows us to account for uncertainties related to future emissions trajectory as well as the simulated atmospheric response  . Unlike current weather data, the projected data should be able to represent a range of future possible scenarios complying with a range of observed and projected statistical characteristics of variable of interest.
Future daily weather data was generated with Marksim v.2  . The Marksim v.2 is a weather generator that produces daily weather variables by mimicking the observed weather variables and their statistical properties. It is a third-order Markov rainfall generator and the program has been calibrated with more than 10,000 stations worldwide, clustered into 702 climate clusters. Current version of Marksim requires the input data from the most recent IPCC AR5 models. Additional parameters such as latitude, longitude, and elevation are required for generating daily-time-step weather data for a given location. MarkSim has been applied to various climate impact studies in agriculture field. Marksim generated future weather data was to estimate the number of Helicoverpa armigera Hub for eight pigeon pea (Cajanus cajan L.) growing locations in India for 2020, 2050, and 2080  . Marksim generated data were also used in studying future
Table 1. Four state-of-the-art general circulation models (GCMs) included in the Coupled Model Intercomparison Project Phase 5 (CMIP5)  used in our study.
availability of recreational cultural ecosystem services with changing climate  . Marksim was used to generate 21 independent climate years with 30 replications to achieve a maximum inter-annual weather variability and gain better insight into future weather patterns. Creation of DSSAT weather data file format “*.WTG” using Marksim was automated and the program was executed with different seed values (starting value used by a random number generation routine to crease random numbers) for each of the 30 replications of weather file. No bias correction with respect to the baseline climatology were applied to the downscaled GCM data.
2.5. Model Inputs
A total of six future climate scenarios (four GCMs, their ensemble, and the ensemble of all available GCMs in the CMIP5 projection data) were used for RCP6.0 and RCP8.5. Future atmospheric CO2 concentration used in the simulations were obtained from the freely available RCP Database v2.0.5 (http://tntcat.iiasa.ac.at:8787/RcpDb/dsd?Action=htmlpage&page=download).
Time for emergence, anthesis, physiological maturity and wheat yield were calculated for each year of the baseline and future scenarios. Cropping system in Oklahoma and most part of the Southern Plains is done under water limited conditions. The region has faced several severe droughts, including the recent drought of 2012. Therefore, simulations were designed by considering both non-limiting and limiting growing conditions. Three scenarios of wheat growing conditions were simulated, i) Actual (Ywf)―with water and nitrogen limiting conditions, ii) Attainable (Ya)―nitrogen limiting condition, and iii) Potential (Yp)―non-limiting conditions.
2.6. Computational Requirements
There were total of 1080 individual runs for 1617 locations-years, which included 30 replications. Running simulation on a personal computer would require a significant amount of time. Therefore, the simulations were executed on Linux cluster at the High Performance Computing Center, Oklahoma State University. Replicates were averaged for each year for a location. Performing such a large simulation required 72 nodes and a total of 96-hour wall-clock time was needed to complete the simulations.
2.7. Model Evaluation
Model performance was evaluated by comparing simulated and the USDA-NASS reported winter wheat yields for the entire 35-year period. Wheat yield was investigated across the state at a fine resolution using the DSSAT-CERES Wheat model and the outputs of 4 GCMs, their ensemble, and ensemble of 17 GCMs under RCP6.0 and RCP8.5 scenarios. We used two different baseline datasets, one derived from the Oklahoma Mesonet and the other from DAYMET dataset. Considering data source, GCMs, forcing scenarios, and yield scenarios, simulation were carried out for 72 factors for 77 sites for 20 years. Simulation were replicated 30 times with the weather data generated using different seed number for each location. For analysis purpose, simulation output for the replications were averaged and yield across the sites were averaged for each year and plotted.
3. Results and Discussion
The crux of this manuscript is to assess how the DSSAT-CERES-Wheat crop model predicts wheat yield and yield changes when the model is driven by different downscaled GCMs and their ensembles across two different RCPs. Fewer studies have assessed wheat yield performance on such a finer scale. Hence, the result and discussion section will mostly focus on the differences between the projected yield changes compared with the baseline scenario yield rather than crop phenology and physiology.
3.1. Projected Climate Change
Both baseline data looked similar in distribution (Table 2). However, the distribution of Daymet recorded SRAD for baseline show higher values compared to Mesonet. Distribution of downscaled annual average data show higher variability in terms of data distribution compared to the observed baseline (Figure 3). In
Table 2. Baseline (1980-2014) average annual summary for rainfall, solar radiation, maximum and minimum temperature for Daymet and Mesonet Data at 77 Mesonet sites in Oklahoma.
Figure 3. Distribution of annual averages of maximum temperature, minimum temperature, solar radiation, and total rainfall downscaled using Marksim for the period 2040-2060 under two RCP forcings of 6.0 and 8.5 along with baseline (1980-2014) weather data for two different weather data sources: Mesonet and Daymet. Data are averaged across select 77 Mesonet locations.
addition, median values of Daymet data appear higher than Mesonet data for most of the time. Projected average annual weather data for RCP6.0 and RCP8.5 show increase in temperature (maximum and minimum), and decrease in solar radiation and total rainfall amount (Figure 3). Percent increase or decrease of projected four weather variables for their corresponding GCMs and RCPs compared to the baseline data of Mesonet and Daymet are summarized in Table 3. Percent changes of average annual rainfall, solar radiation, maximum temperature, and minimum temperature relative to the baseline climate. On average, for RCP6.0 forcing, across the different GCMs, Daymet generated projection show 17.72% ± 1.24% decline in average total rainfall amount while Mesonet generated projection show 26.7% ± 1.24% decline in average total rainfall amount relative to their baseline values. Likewise, for RCP8.5 forcing, average total rainfall is expected to decline by 16.84% ± 1.77% and 25.70 ± 1.42 using Daymet and Mesonet data respectively.
In case of SRAD, the values are projected to decline at a range of 3.15% - 7.15% in relation to the baseline data, with a higher reduction for RCP6.0 forcing. Maximum temperature is projected to increase up to 10.13% for RCP 6.0 and 14.44% for RCP8.5. On an average, rise in maximum temperature was projected to be approximately 2.11˚C RCP6.0 and 3.04˚C for RCP8.5. However, a greater increase in minimum temperature is projected, ranging from 19.59% - 36.61%. Increase of an average of 2.0˚C was computed for the RCP6.0 forcing and 2.88˚C forRCP8.5 forcing. Higher percentages are observed for RCP8.5 forcing. Coefficient of variation (CV, %) for all the estimates average ranged somewhere between 0.82 - 5.34. CV were lowest for SRAD, followed by RAIN, TMAX, and TMIN. In addition, CV were lower for RCP6.0 compared to RCP8.5 for the four weather variables.
3.2. Future Wheat Yield Simulation
For each site, historical yield was averaged across the 1984-2014 (Figure 4) and future yields were averaged across 30 different replications for each site for each year and plotted (Figure 5). Compared to the baseline yield, downscaled GCMs using Daymet data, under RCP6.0 forcing showed average increase in wheat yield by 9.65%, 2.46%, for Ywf and Ya scenario respectively. However, average decrease of 7.45% in yield was projected for Yp scenarios. Under RCP8.5 forcing, increase in yield was observed across three yield scenarios, ranging between 4.26% - 19.68% (Table 4).
There is a consistency in the projected winter wheat yield between the GCMs. Across all the yield scenarios and under the forcing scenarios of RCP6.0 and RCP8.5, MRIOC5 consistently has the lowest yield. The yields converge for Ywf scenario, however, the difference between the GCMs and between the data sources increase for Ya and Yp scenarios under both RCP forcings. However, we can’t identify which downscaled GCM is performing better. One way to further improve the credibility of this study is to generate downscaled data for each of
Table 3. Percent changes of average annual rainfall, solar radiation, maximum temperature, and minimum temperature relative to the baseline climate.
Figure 4. Box-and-whisker plots of DSSAT simulated winter wheat yield during the 30-year (1984-2014) at 77 different Mesonet locations in Oklahoma for three yield scenarios, i) Actual (Ywf)―water and nitrogen limiting, ii) Attainable (Ya)―nitrogen limiting, and iii) Potential (Yp)―non-limiting conditions using two different weather data sources, Daymet (shaded blue) and Mesonet. From top to bottom the horizontal lines in the box plots show the 95th, 75th, 50th, 25th, and 5th percentiles.
Figure 5. Projected winter wheat yield of 2040-2060 relative to 1980-2014, averaged from 77 locations across the Oklahoma Mesonet locations using different GCMs, their ensembles under the forcing scenarios of RCP6.0 and RCP8.5.
the RCPs for baseline time period, compute the annual mean and compare against the known baseline mean weather data  . In addition, yield can be also compared across different downscaled GCMs in order to identify model that
Table 4. Percent change of simulated winter grain yield for the period 2040-2060 relative to the baseline climate (1980-2014) using different General Circulation Models (GCMs) and Representative Concentration Pathways (RCPs) 6.0 and 8.5 for three different growing scenarios: actual (Ywf), attainable (Ya), and potential (Yp).
explains the weather pattern better in this region.
Comparison of seasonal cumulative rainfall and crop evapotranspiration averages across all the site, years showed that most of the instances, evapotranspiration values were always less than rainfall across all the downscaled GCM weather data, radiative forcings, and baseline data sources (Table 5). Average seasonal rainfall for simulations performed with Daymet baseline weather data for Ywf yield scenario and RCP6.0 was 569.39 ± 16.75 mm for which the estimated evapotranspiration value was 375.98 ± 2.34 mm. For the same simulation scenario under RCP8.5 forcing, the seasonal average rainfall value was 566.98 ± 22.98 mm with an estimated crop evapotranspiration value of 373.66 ± 4.36 mm. These values are similar to the baseline seasonal rainfall (561 ± 137 mm) and evapotranspiration (375 ± 46.58 mm) for both Daymet and Mesonet data (values not shown). Comparison of number of days from planting to harvest showed that the duration decreased under projected climate data (Table 6). This could result from increase in temperature stress. As the crop senses stress, it is more likely to shorten its life cycle, causing an undesirable increase in wheat grain protein content   .
Results from our simulation are similar to those reported by  . Using the downscaled data from the Regional Climate Models CanESM2 and CANRCM4 and DSSAT-CERES-Wheat model, they report an increase in spring wheat yield between 10% and 20% averaged across the Canadian Prairies. However, yield in
Table 5. Average growing season rainfall and average crop evapotranspiration generated from DSSAT-CERES-Wheat simulation performed by using six different Global Circulation Models (GCMs) under two different forcings of RCP using two different baseline data sources for actual (Yp) and attainable (Yp) growing scenarios. Values are averaged across 77 different Mesonet locations.
Table 6. Median number of days between sowing and harvest for winter wheat, over the baseline period 1980-2014 and future climate projections 2040-2060 under three different growing scenarios, actual (Ywf), attainable (Ya), and potential (Yp).
the northern latitudes are projected to increase. DeNitrofication-DeComposition (DNDC) model used to simulate spring wheat yield under projected climatic conditions found that yield increased across IPCC SRES climate change scenarios  . Previous study in wheat productivity in Central Oklahoma also reported some small changes in wheat yield ±5%, with a 5% decrease in B2a scenario and 5% increase in GGal scenarios. Increase in yield could be due to the direct physiological effects of the increased atmospheric conditions.
The significance of this study is that it helps to evaluate the potential impacts of future climate change on wheat enterprise on a local scale on different temporal scales. There are no tools currently available for preparing future weather input files from observed weather data on a local level. We undertook a spatially explicit simulation to assess how winter wheat productivity might be impacted by mid-21st century climate change. Simulations were carried out for 77 target Mesonet location across the State of Oklahoma. Some of the existing limitations in the model is that the crop model currently does not simulate the effect of high temperatures on pollen viability; the effect of CO2 on other processes besides photosynthesis and transpiration; and the interaction between temperature, drought, and increased CO2 concentration  . Likewise, uncertainties in our study reside in the prediction of the climate models used, location of site used for simulation, and simulation output of the DSSAT model. Management practices, fertilizer inputs, and cultivar were set constant and were not changed from our baseline simulations. Therefore, the results may not hold under similar climatic conditions. We also need to stress that these simulations do not take into account of pests and disease potentials, and loss due to natural calamities such as flooding. The idea of using scenarios is not to predict future but to better understand the uncertainties for judicious decision making for a range of future conditions. In our study, we chose only two scenarios to work with because of computational and storage resources and time limitations. However, we should also keep in mind that the greater proportion of uncertainty in climate change impact assessment has been credited to variations among crop models rather than variations among the downscaled GCMs  . Our objective was not to compare the differences between the crop models, but to examine the yield differences across future climate projection based on the downscaled GCMs and their ensembles.
Although simulations suggest yield increase, further study is required to examine how soil quality and additional challenges of pest, irrigation, soil degradation, etc. may impact sustained crop production on a regional level. Also, study must verify if lodging may become issue for Ya and Yp scenarios. We only examined the annual mean data across the sites, however, the crop is sensitive to fluctuation in weather variation depending on the growth stage. Water availability for crop production in future is also uncertain, therefore, crop production with less water can be problematic if the temperature is also hotter. Since Marksim is an example of a stochastic, statistical downscaling method, with low input data requirement to quantify the relationship between crop growth and weather, it may be less adaptable to different climatic conditions both spatially and temporally  . And process-based models such as the DSSAT CERES-Wheat include more non-linear processes, so that the spatial aggregation of inputs over heterogeneous landscape is less recommended   .
This paper describes a methodology for rapid synthesis of GCM-based, spatially explicit, high resolution future weather data inputs for the DSSAT crop model, for cropland area across wheat growing regions of Oklahoma on a seamless temporal scale. We examined how climate change may impact regional food security using a “bottom-up” approach. Most previous studies have focused on how climate change influence crop yield at regional and/or national levels. This study is a marked advance in the compilation of high resolution dataset, as such data offers tremendous opportunity to quantify future yield status for other crops within the region. With advancement in climate models in term of their fidelity in simulating natural processes and generating reliable output at higher spatio-temporal scales, we would be able to better assess the magnitude of the impact of climate changes on crop production and its societal impact. Based on our experience in implementing simulations for multi decadal, fine-scale study, we argue that a good understanding of scientific computing along with crop science and physiology is extremely helpful from crop modeling perspective. Focus on replicability, transparency, and code efficiency are seldom discussed in the literature. We believe that addressing these aspects is also of significance importance for advancement in crop modeling. Therefore, we conclude that by minimizing the yield gap with improved crop husbandry, increased input efficiency for water, nutrients, and pesticides, along with improved wheat cultivars, and sustainable production practice, wheat production in this region can be sustained for changing climate to meet the nutritional demand of people in years to come.
We acknowledge at support from Oklahoma Agricultural Experiment Station, Hatch project #2939; Oklahoma Wheat Commission, the National Science Foundation Grant No. OIA-1301789 and the OSU High Performance Computing Center through NSF Grant No. OCI-1126330 for providing the computational resources required for this study.