Using Satellite Data to Optimize Wheat Yield and Quality under Climate Change

Climatic conditions during the grain-filling period are a major factor affecting wheat grain yield and quality. Wheat in many semi-arid and arid areas faces high-temperature stress during this period. Remote sensing can be used to monitor both crops and environmental temperature. The objective of this study was to develop a tool to optimize field management (cultivar and sowing time). Analysis of 155 cultivar experiments (from 10 growth seasons) representing different environmental conditions revealed the required degree-days for each Israeli spring wheat cultivar to reach heading (from emergence). We developed a Google Earth Engine (GEE) app to analyze time series of gap-filled 1 km MODIS land surface temperature (LSTcont). By changing the cultivar and/or emergence date in the GEE app, the farmer can “expose” each wheat field to different climatic conditions during the grain-filling period, thereafter enabling him to choose the best cultivar to be sown in the field with the right timing. This approach is expected to reduce the number of fields that suffer from heat stress during the grain-filling period. The app can be also used to assess the effects of different global warming scenarios and to plan adaptation strategies in other regions too.


Introduction
Climatic changes are likely to impact food production and pose serious challenges to food security. Remote sensing can help agricultural systems to adapt measures to negate the adverse impacts of climate change [1]. Wheat (Triticum aestivum L.) is one of the vital cereal crops grown worldwide, providing the primary caloric and nutritional source for a significant portion of the human population [2]. In order to ensure food security, there is a need for strategies and policies for the mitigation of food shortages. Timely and accurate estimates of global crop production are also essential. The adaptability of wheat to a wide range of environments has been favored by allelic diversity in genes regulating growth habits [3] i.e., temperature (VRN genes) and photoperiod response (PPD genes). Climate conditions during the grain-filling period are a major factor affecting wheat grain yield and quality. Wheat in many semi-arid and arid areas is facing high-temperature stress during this growth period. This challenge will probably be relevant in other growing areas in the future. Climatic variation is considered a main factor affecting gluten quality within a given variety when tested across different sites and/or seasons [4][5][6]. Growing Degree Days (GDD) is a main factor in crop models which are used to support decision-making and planning, particularly in relation to climate change issues [7][8][9][10].
Remote sensing can be used to monitor environmental temperatures. Surface skin temperature observations (T skin ), as obtained from satellite data, provide useful climatological data of high spatial resolution and global coverage that enhances the traditional ground observations of surface air temperature (T air ). MODIS land surface temperature (LST) products retrieved from thermal infrared bands have an accuracy of <1 • C in the range from −10 to 50 • C. Using MODIS data for LST retrieval improves the quality of the LST product and the diurnal feature in the product due to better temporal, spatial and angular coverage of clear-sky observations [11]. The 1 km MODIS-LST data can reveal the large spatial variability of LST and potentially increase precision in forecasting and applying management strategies [12]. However, the relationship between T skin and T air may be complex as it is affected by numerous factors, such as vegetation cover, incoming solar radiation and wind [13]. Therefore, data must be treated with adequate procedures, taking into account the temporal variation, crop cover and phenology, etc. [14,15].
In Mediterranean rain-fed conditions, such as in Israel, the amount of rainfall and its distribution during the growing season has a marked influence on spring wheat yield and quality [4,5,[16][17][18][19]. Drought or heat shock during the grain filling period results in smaller yields and test weight, due to kernel shriveling. Sowing is recommended in Israel in a window from 15 November to 15 December based on the short rainy season, an increasing water deficit and thermal stress during grain filling and a higher probability of chilling/frost in January-February during anthesis. In Israel, late-heading cultivars can be sown from November 1st, resulting in a high variation between the commercial field's emergence and phenological stage. Furthermore, grain quality standards are not consistently met in any given year due to adverse weather, which is the overriding factor controlling grain yield, grain test weights and protein concentration [4]. Hence, the sowing date in each field should be determined precisely, based on the wheat genetic and air temperature considering the spatial variability. Reaching heading can be predicted inseason by modeling [20] or by remote sensing [21,22]. However, there is a gap of knowledge for predicting heading date pre-sowing in advance. The objective of this study is to develop, calibrate and validate algorithms for the retrieval of wheat heading stage using time series of satellite data. We used the Google Earth Engine (GEE) platform to provide non-expert stakeholders (farmers) with a tool for optimizing yield by choosing the best sowing date and cultivar.

Materials and Methods
In order to develop an algorithm that provides estimation of the heading date for a specific cultivar, given the field location and the emergence date, the following tasks are needed:

1.
A table of GDD required to reach the heading stage from the emergence stage for each cultivar. This is described in Section 2.1.

2.
To calculate GDD for each field. This is described in Section 2.2.

3.
To disseminate the algorithm as a tool to all stakeholders (farmers and the scientific community). This is described in Section 2.3.

Wheat Fields
This study uses three datasets that contain the emergence and heading dates from different locations ( Figure 1): one dataset from Gilat Research Center is used for calibration (Gilat-CAL), and two other datasets are used for validation: one from commercial wheat fields at Saad (Saad-VAL), and another from 28 wheat variety test experiments across Israel (CVexp-VAL).

1.
The Gilat-CAL dataset contains data from 76 wheat variety test experiments from ten growth seasons (2008-2017) at Gilat Research Center. Each experiment included several spring bread wheat cultivars spanning the phenological range of Israeli elite germplasm. Each experiment typically included 8-15 cultivars, with a total of 22 different cultivars. The dataset comprises more than 1500 heading date records, taken by one person (DJ. Bonfil). Daily meteorological data from Gilat meteorological station was used. The meteorological station was typically within 1 km distance from the field, with a maximum distance of 2.5 km. This data set was used to calculate the GDD (with a temperature threshold of 0 • C) and the number of days from emergence to heading (E-H) for each cultivar.

Satellite and Numerical Weather Prediction Model Data
MODIS onboard Terra and Aqua polar-orbiting sun-synchronous satellites provide global coverage daily LST at 1 km spatial resolution at 10:30 (Terra) and 1:30 (Aqua) AM/PM local time. We used 2002-2020 Aqua LST data since it better represents the maximum and minimum diurnal LST. The MYD11A1 V6 product [23] was used to retrieve the seasonal behavior of LST in cloud-free conditions. Numerical weather prediction (NWP) models provide surface temperature regardless of cloud cover. We used NCEP Climate Forecast System Version 2 (CFSv2) NWP model surface air temperature 2 m above ground at a spatial resolution of 0.2° [24] to fill gaps in the LST dataset due to cloud cover [25]. Forecasts are initialized four times per day (0000, 0600, 1200 and 1800 UTC). The data is available from 1979 to the present.
Temporal Fourier Analysis (TFA) was used to derive the climatological temperatures of the satellite (MODIS LST) and model (CFSv2) temperatures. The TFA describes the seasonal cycles of temperature in terms of annual, bi-annual and tri-annual components (or "harmonics"), each described by its phase and amplitude. These Fourier harmonics may be recombined, providing a smoothed signal, which is regarded here as the climatological expected temperatures: (1) where t is the Julian date, LSTclim(t) is the climatological temperature at time t, is the mean annual temperature, is the amplitude of the i th harmonic component and n is the number of harmonic components. We used here the first three harmonics (n = 3), following Scharlemann et al. [25] and Lensky and Dayan [12].
is the phase and is the frequency ( 2 /365) of the i th harmonic component. TFA was applied on both MODIS LST and CSFv2 surface (2 m air) temperature to derive the climatological temperature and

Satellite and Numerical Weather Prediction Model Data
MODIS onboard Terra and Aqua polar-orbiting sun-synchronous satellites provide global coverage daily LST at 1 km spatial resolution at 10:30 (Terra) and 1:30 (Aqua) AM/PM local time. We used 2002-2020 Aqua LST data since it better represents the maximum and minimum diurnal LST. The MYD11A1 V6 product [23] was used to retrieve the seasonal behavior of LST in cloud-free conditions. Numerical weather prediction (NWP) models provide surface temperature regardless of cloud cover. We used NCEP Climate Forecast System Version 2 (CFSv2) NWP model surface air temperature 2 m above ground at a spatial resolution of 0.2 • [24] to fill gaps in the LST dataset due to cloud cover [25]. Forecasts are initialized four times per day (0000, 0600, 1200 and 1800 UTC). The data is available from 1979 to the present.
Temporal Fourier Analysis (TFA) was used to derive the climatological temperatures of the satellite (MODIS LST) and model (CFSv2) temperatures. The TFA describes the seasonal cycles of temperature in terms of annual, bi-annual and tri-annual components (or "harmonics"), each described by its phase and amplitude. These Fourier harmonics may be recombined, providing a smoothed signal, which is regarded here as the climatological expected temperatures: where t is the Julian date, LST clim (t) is the climatological temperature at time t, LST is the mean annual temperature, A i is the amplitude of the i th harmonic component and n is the number of harmonic components. We used here the first three harmonics (n = 3), following Scharlemann et al. [25] and Lensky and Dayan [12]. ϕ i is the phase and ω i is the frequency and CSFv2 surface (2 m air) temperature to derive the climatological temperature and the anomaly (the deviation of the actual temperature from the climatological temperature). Surface temperature at a specific time and date can be regarded as being composed of two components: (a) the long-term mean (climatology), and (b) the deviation from that mean due to the weather (anomaly). Shiff et al. [26] proposed a method based on these principles and provided a continuous (gap-filed) LST dataset (LST cont ) based on satellite (MODIS LST) and model (CFSv2) temperature data. Figure  the anomaly (the deviation of the actual temperature from the climatological tempera ture). Surface temperature at a specific time and date can be regarded as being compose of two components: (a) the long-term mean (climatology), and (b) the deviation from tha mean due to the weather (anomaly). Shiff et al. [26] proposed a method based on thes principles and provided a continuous (gap-filed) LST dataset (LSTcont) based on satellit (MODIS LST) and model (CFSv2) temperature data. Figure   The growing degree days (accumulated degree days with temperature threshold of 0 °C) are plotted for two emergence dates: November 1 and December 1. The differences between above average temperatures 2018 and below average 2017 resulted in eight days difference at heading for the very early cultivar and is marked as a double-headed arrow in mid-February.

Google Earth Engine (GEE)
Google Earth Engine (GEE) is a parallel computation service platform for advance image analysis that hosts a variety of remote sensing and geospatial datasets [27]. Bot products (MODIS LST and CFSv2) are available on the GEE platform. GEE leverages it cloud computing service for analysis capabilities of time series data, which are otherwis heavy consumers of time and computation resources. Our research uses ~20 years' time series of daily datasets that are all hosted on the GEE platform to calculate climatologica LST (LSTclim) and continuous gap-filled LST (LSTcont). GEE also helps researchers to easil disseminate their products to other users and researchers, which enabled us to provid code and a web app, for the benefit of farmers and the science community. The presen study provides a GEE web app that enables optimization of yield by providing the farme with the expected heading date for each combination of field, cultivar and emergence dat (see Supplementary Materials).

Results
A large difference in wheat heading dates, from late December until April, was foun in the wheat fields. Differences in meteorological conditions, genotypes, sowing an emergence dates explain this variability. Since all cultivars are spring wheat that do no require vernalization, the emergence to heading (E-H) period was mainly affected by tem perature and photoperiod (day length) response. Early and late emergence exposes th The growing degree days (accumulated degree days with temperature threshold of 0 • C) are plotted for two emergence dates: November 1 and December 1. The differences between above average temperatures 2018 and below average 2017 resulted in eight days difference at heading for the very early cultivar and is marked as a double-headed arrow in mid-February.

Google Earth Engine (GEE)
Google Earth Engine (GEE) is a parallel computation service platform for advanced image analysis that hosts a variety of remote sensing and geospatial datasets [27]. Both products (MODIS LST and CFSv2) are available on the GEE platform. GEE leverages its cloud computing service for analysis capabilities of time series data, which are otherwise heavy consumers of time and computation resources. Our research uses~20 years' timeseries of daily datasets that are all hosted on the GEE platform to calculate climatological LST (LST clim ) and continuous gap-filled LST (LST cont ). GEE also helps researchers to easily disseminate their products to other users and researchers, which enabled us to provide code and a web app, for the benefit of farmers and the science community. The present study provides a GEE web app that enables optimization of yield by providing the farmer with the expected heading date for each combination of field, cultivar and emergence date (see Supplementary Materials).

Results
A large difference in wheat heading dates, from late December until April, was found in the wheat fields. Differences in meteorological conditions, genotypes, sowing and emergence dates explain this variability. Since all cultivars are spring wheat that do not require vernalization, the emergence to heading (E-H) period was mainly affected by temperature and photoperiod (day length) response. Early and late emergence exposes the wheat to warmer conditions and longer days, therefore decreasing the E-H period (for all cultivars). We calculated the growing degree days (GDD, the accumulated daily average temperature above a base temperature threshold of 0 • C) for the E-H period. The GDD (E-H) decreases as the emergence date occurs later than mid-December (Figure 3), since late emergence exposes the wheat to warmer conditions and longer days, therefore decreasing the E-H period (for all cultivars). Phenological ranking between cultivars is kept across the environments; cultivars Yuval (Dariel) are the earliest (latest) with GDD E-H of 1000 (1385 • C), respectively. The calculated GDD for each cultivar can be seen in Figure 3.
The following five cultivars have been determined for the app development: Very early, Early, Intermediate, Late and Very late (1000, 1100, 1200, 1300 and 1400 • C, respectively). wheat to warmer conditions and longer days, therefore decreasing the E-H period (for all cultivars). We calculated the growing degree days (GDD, the accumulated daily average temperature above a base temperature threshold of 0 °C) for the E-H period. The GDD (E-H) decreases as the emergence date occurs later than mid-December (Figure 3), since late emergence exposes the wheat to warmer conditions and longer days, therefore decreasing the E-H period (for all cultivars). Phenological ranking between cultivars is kept across the environments; cultivars Yuval (Dariel) are the earliest (latest) with GDD E-H of 1000 (1385 °C), respectively. The calculated GDD for each cultivar can be seen in Figure  3. The following five cultivars have been determined for the app development: Very early, Early, Intermediate, Late and Very late (1000, 1100, 1200, 1300 and 1400 °C, respectively). Usually, heading of the earliest emerging experiment occurred under cooler conditions (Figure 4). In the 2017-2018 growing season stress conditions (high temperatures) prevailed during the heading and grain-filling period, where the late phenology cultivars faced very high temperatures above 40 °C. Similar behavior was observed each season, with warmer conditions associated with later crop heading. The Heading GEE web application (https://shilosh.users.earthengine.app/view/ emergence2heading, accessed on 21 May 2021) is based on pre-computed daily climatological MODIS LST (LST clim ) at a resolution of 1 km for the entire world. The interface includes a map to select the field and to see the predicted heading map, a calendar for selecting the emergence date and a drop-down list to select the cultivar ( Figure 5). The user can select the emergence date and cultivar and visualize a map of the predicted number of days to heading (E-H). The user can then click on any point on the map and get the date of the predicted heading at that point.  During the 2017-2020 growing seasons, heading date estimations have been validated by two datasets: (a) from one farmer (his own fields) and (b) from the extension service advisers who established experiments across different growth conditions. Both achieved similar correlation (R 2 0.53-0.6), high bias (15)(16)(17)(18) and SEP of 6 to 7 days (Figure 6a,b and Table 1). Since there was high variation between the prevailed air temperatures along the four validation growing seasons, each season exhibits a different relation (Figure 6a (Figure 6e,f, Figures 7 and 8e,f       GEE app accessibility enables farmers to use the GEE heading app in the growing seasons since 2019 for planning sowing, cultivar X field X date (i.e., Genotype X Environ ment) using the climatological LST (LSTclim). Heading during warmer (cooler) seasons i expected to be earlier (later) than predicted when using LSTclim. Decision support during the season should include updated GDD. Using actual LST data (LSTcont) reduced the bia from 15.2 to 3.5 days in Saad-VAL, and from 17.9 to 11.4 days in CVexp-VAL (Table 1). R increased from ~0.56 to ~0.7; SEP decreased in Saad-VAL to 4.7 days. Moreover, the rang of the difference decreased and reached mean ± 6 or 8 days for 80% of the fields, and ± 3 or 5 days for 50% of the fields, Saad-VAL and CVexp-VAL respectively (Figure 8c,d) LSTcont solves most of the trend shown for CVexp-VAL (Figure 8e,f), where the differenc between the estimated and actual heading was related to timing: underestimation for ear lier heading and overestimation for late heading. One of the CV-experiment data can b determined as outlier (all cultivars), as their DOY was recorded at about day 70, whil according to the temperature it should have occurred one month earlier. Still, the mean difference distribution in the CVexp-VAL data is higher than 0, showing underestimation of the GEE application.  GEE app accessibility enables farmers to use the GEE heading app in the growing seasons since 2019 for planning sowing, cultivar X field X date (i.e., Genotype X Environment) using the climatological LST (LST clim ). Heading during warmer (cooler) seasons is expected to be earlier (later) than predicted when using LST clim . Decision support during the season should include updated GDD. Using actual LST data (LST cont ) reduced the bias from 15.2 to 3.5 days in Saad-VAL, and from 17.9 to 11.4 days in CVexp-VAL (Table 1). R 2 increased from~0.56 to~0.7; SEP decreased in Saad-VAL to 4.7 days. Moreover, the range of the difference decreased and reached mean ± 6 or 8 days for 80% of the fields, and ± 3 or 5 days for 50% of the fields, Saad-VAL and CVexp-VAL respectively (Figure 8c,d). LST cont solves most of the trend shown for CVexp-VAL (Figure 8e,f), where the difference between the estimated and actual heading was related to timing: underestimation for earlier heading and overestimation for late heading. One of the CV-experiment data can be determined as outlier (all cultivars), as their DOY was recorded at about day 70, while according to the temperature it should have occurred one month earlier. Still, the mean difference distribution in the CVexp-VAL data is higher than 0, showing underestimation of the GEE application.
Temperatures are expected to rise due to climate change. To characterize how climate change could potentially impact the E-H period length, we evaluated future projections from the NASA Earth Exchange Global Daily Downscaled Projections (NEX-GDDP) ensemble [28], using the Coupled Model Intercomparison Project 5 (CMIP5 [29]) under the Representative Concentration Pathway 4.5 (RCP4.5) emission scenario, with the Australian Community Climate and Earth-System Simulator (ACCESS1-0). Figure 9 shows the projected changes on three cultivar types (Very early, Intermediate and Very late) for the period 2010-2099. The E-H period length on Gilat is expected to be shorten by 11 days with the Very late cultivar until the end of the century. from the NASA Earth Exchange Global Daily Downscaled Projections (NEX-GDDP) ensemble [28], using the Coupled Model Intercomparison Project 5 (CMIP5 [29]) under the Representative Concentration Pathway 4.5 (RCP4.5) emission scenario, with the Australian Community Climate and Earth-System Simulator (ACCESS1-0). Figure 9 shows the projected changes on three cultivar types (Very early, Intermediate and Very late) for the period 2010-2099. The E-H period length on Gilat is expected to be shorten by 11 days with the Very late cultivar until the end of the century.

Discussion
Crop development is affected by temperature and is therefore a key factor in many crop models. The ability to use temperature data from space opens a new era. However, most models use air temperature, measured 2 m above the land surface, while Earth observations provide land surface temperature. Our results show that LST can be used without any transformation during the winter, which is the wheat growing season in Israel. These results are in accordance with previous studies [14,30,31]. In other growth seasons, adjustments between GDD from meteorological stations and from satellites will probably be needed, especially for summer crops. In addition, LSTcont enables using near to realtime MODIS-LST data, so LST can be incorporated into in-season decision support systems and/or crop models. However, LSTcont cannot predict future temperatures, therefore LSTcont can only provide GDD up to the current time, while a farmer who wants to plan the right sowing date and cultivar for the next season is limited to the use of LSTclim only. LSTcont enables using LST even in cloudy conditions, though validation is needed for growing areas where clouds persist for long periods.
Spatial temperature variations are small in some growing areas but can be very high in others. For example, variations between Saad fields (Saad-VAL) are small, while the CVexp-VAL dataset is based on high variation in growth conditions ( Figure 1). Moreover, even for the Saad-VAL dataset, rain and/or irrigation creates high variation in the temperatures that the crops are exposed to during the different growth stages resulting in variance in emergence dates. The most important result achieved in this study is that the GEE based algorithm can forecast the heading date of each cultivar in every field. Therefore, it was used immediately for planning sowing, cultivar X field X date (i.e., G X E) by Israeli farmers.

Discussion
Crop development is affected by temperature and is therefore a key factor in many crop models. The ability to use temperature data from space opens a new era. However, most models use air temperature, measured 2 m above the land surface, while Earth observations provide land surface temperature. Our results show that LST can be used without any transformation during the winter, which is the wheat growing season in Israel. These results are in accordance with previous studies [14,30,31]. In other growth seasons, adjustments between GDD from meteorological stations and from satellites will probably be needed, especially for summer crops. In addition, LST cont enables using near to real-time MODIS-LST data, so LST can be incorporated into in-season decision support systems and/or crop models. However, LST cont cannot predict future temperatures, therefore LST cont can only provide GDD up to the current time, while a farmer who wants to plan the right sowing date and cultivar for the next season is limited to the use of LST clim only. LST cont enables using LST even in cloudy conditions, though validation is needed for growing areas where clouds persist for long periods.
Spatial temperature variations are small in some growing areas but can be very high in others. For example, variations between Saad fields (Saad-VAL) are small, while the CVexp-VAL dataset is based on high variation in growth conditions ( Figure 1). Moreover, even for the Saad-VAL dataset, rain and/or irrigation creates high variation in the temperatures that the crops are exposed to during the different growth stages resulting in variance in emergence dates. The most important result achieved in this study is that the GEE based algorithm can forecast the heading date of each cultivar in every field. Therefore, it was used immediately for planning sowing, cultivar X field X date (i.e., G X E) by Israeli farmers.
Several factors can affect the accuracy, where "date determination" is the most important. Under good germination conditions, the emergence date can be easily determined (±1 day). However, hard conditions due to water stress, bad seed bed (large clods), variation in sowing depth and others can lead to different determination of the emergence date by different persons. Determining the heading date itself is challenging as not all plants reach this stage at the same day. Moreover, under high stand conditions differentiation between the main stem spike and tillers cannot be done. Hence, the significance of the heading date for the same field is somehow fuzzy. In the case where one person determines all fields, the date, emergence and heading, determination would be more consistent.
Therefore, validation based on Saad-VAL reached better accuracy than CVexp-VAL which was monitored by more than one person.
Most inaccuracies in the Saad-VAL dataset are within 0 ± 5 days (Figure 8e), while for the CVexp-VAL most are within 5 ± 5 days (Figure 6f). The 5 day bias can mainly be attributed to the heading determination. Heading stages have been determined in the Gilat-CAL and Saad-VAL datasets when about 100% of the plants reached this stage and the emergence of the main tiller ear was complete, but not the tillers' ears, while most extension advisors determine the heading date when emergence of about 80% of the whole field's (including tillers) ears is complete. As tillers' ears usually emerge a few days later, this explains the 5 day bias between the two validation sets. Crop modeling to predict phenology within season using real-time data reached a similar accuracy of ±3-7 days [20]. Our results exhibit that degree days can be used with similar accuracy, but more importantly, GDD can be used out of season-and even pre-season-to assist decision making.
The LST data has a spatial resolution of 1 km per pixel. Most fields are smaller than 1 km, thus the estimated field LST represents a mixed pixel. The spatial resolution can be improved using fusion schemes with Landsat data to 30 m [32,33]. However, the fusion procedure requires time and computation resources, and since the 1 km data have been found sufficient, there is no need to improve the spatial resolution.
Other sources of inaccuracy can be related to the Gilat-CAL dataset. Although calibration has been done by one person, and all data were collected from one area, for a few cultivars insufficient data was available. Consequently, their GDD requirement was not accurate enough, leading to miscalculation by the GEE-heading app.
The main benefit of using the GEE-heading app would be achieved by reducing the number of fields that suffer from heat stress during the grain-filling period. In case of an early rain event where germination occurs immediately after sowing, farmers can use the GEE-heading app to find the best sowing date for each field. In addition, the GEE-heading app can be used for other purposes as well. Estimation of the heading date can help growers reduce or prevent herbicide application that can cause damage during anthesis. The same is the case for the application of liquid nitrogen as late top-dressing.
Data can be used for interpretation and better understanding of crop development. In cases where heading in vivo exhibits a significant delay, it may be an indication for other stress limiting developments. Limiting factor identification can assist the specific or the following crop. It can be achieved by comparing the same cultivar grown in different locations and growth conditions.
The app can also be used to assess the effects of different global warming scenarios and to plan adaptation strategies in other regions too. Therefore, the environmental conditions during grain filling in a specific field can be modified by setting the optimal sowing date and cultivar. Sowing super early cultivars is a very important option to prevent warm and dry conditions during the grain filling period, encouraging breeders to release super early cultivars. However, shortening the growth season decreases the yield production potential, as crops absorb less radiation. Breeders can learn that in future seasons the growth period, emergence until heading, is expected to be shorter by 11 days (Figure 9). Therefore, they should include late phenology lines in their breeding programs too, that perhaps would not decrease yield production during shorter growth periods. All Israeli spring wheat cultivars ( Figure 3) have identical PPD alleles resulting in similar responses to photoperiod. Hence, breeders may try to include other allele that will regulate development and delay the transition from vegetative apex to floral initiation or enlarge the stem elongation period.
Growth during warmer and cooler seasons is expected to be faster or slower than the seasonal average (during season). The next GEE app version should include within season updated in vivo GDD accumulation (i.e., LST cont ). In addition, heading app validation should be extended to different growth locations and conditions. Furthermore, since the precise field heading date determination is not clear and can differ between persons, validation should include data taken by many people. Another approach would be to test the heading app globally, first with spring wheat all over the world, and perhaps for winter wheat too.

Conclusions
The 1 km MODIS-LST data reveal the large spatial variability of LST. The GEE platform and time series of LST data were used to map the spatial variability for each field. The required GDD for each cultivar to reach heading from emergence was integrated into a GEE app and is available for use by farmers. By changing the cultivar and/or emergence date in the GEE app, the farmer can "expose" each wheat field to different climatic conditions during the grain filling period, thereafter enabling him to choose the best cultivar to be sown in the field. The main benefit would be achieved by reducing the number of fields that suffer from heat stress during the grain-filling period. The app can also be used to assess the effects of different global warming scenarios and to plan adaptation strategies in other regions as well. Therefore, the climate or at least the environment during the grain filling period can be modified.