Crop yield forecasting plays a vital role in coping with the challenges of the impacts of climate change on agriculture. Improvements in the timeliness and accuracy of yield forecasting by incorporating near real-time remote sensing data and the use of sophisticated statistical methods can improve our capacity to respond effectively to these challenges. The objectives of this study were (i) to investigate the use of derived vegetation indices for the yield forecasting of spring wheat (Triticum aestivum L.) from the Moderate resolution Imaging Spectroradiometer (MODIS) at the ecodistrict scale across Western Canada with the Integrated Canadian Crop Yield Forecaster (ICCYF); and (ii) to compare the ICCYF-model based forecasts and their accuracy across two spatial scales-the ecodistrict and Census Agricultural Region (CAR), namely in CAR with previously reported ICCYF weak performance. Ecodistricts are areas with distinct climate, soil, landscape and ecological aspects, whereas CARs are census-based/statistically-delineated areas. Agroclimate variables combined respectively with MODIS-NDVI and MODIS-EVI indices were used as inputs for the in-season yield forecasting of spring wheat during the 2000–2010 period. Regression models were built based on a procedure of a leave-one-year-out. The results showed that both agroclimate + MODIS-NDVI and agroclimate + MODIS-EVI performed equally well predicting spring wheat yield at the ECD scale. The mean absolute error percentages (MAPE) of the models selected from both the two data sets ranged from 2% to 33% over the study period. The model efficiency index (MEI) varied between −1.1 and 0.99 and −1.8 and 0.99, respectively for the agroclimate + MODIS-NDVI and agroclimate + MODIS-EVI data sets. Moreover, significant improvement in forecasting skill (with decreasing MAPE of 40% and 5 times increasing MEI, on average) was obtained at the finer, ecodistrict spatial scale, compared to the coarser CAR scale. Forecast models need to consider the distribution of extreme values of predictor variables to improve the selection of remote sensing indices. Our findings indicate that statistical-based forecasting error could be significantly reduced by making use of MODIS-EVI and NDVI indices at different times in the crop growing season and within different sub-regions.

Cropland productivity is impacted by climate. Knowledge on spatial-temporal patterns of the impacts at the regional scale is extremely important for improving crop management under limiting climatic factors. The aim of this study was to investigate the effects of climate variability on cropland productivity in the Canadian Prairies between 2000 and 2013 based on time series of MODIS (Moderate Resolution Imaging Spectroradiometer) FAPAR (Fraction of Absorbed Photosynthetically Active Radiation) product. Key phenological metrics, including the start (SOS) and end of growing season (EOS), and the cumulative FAPAR (CFAPAR) during the growing season (between SOS and EOS), were extracted and calculated from the FAPAR time series with the Parametric Double Hyperbolic Tangent (PDHT) method. The Mann-Kendall test was employed to assess the trends of cropland productivity and climatic variables, and partial correlation analysis was conducted to explore the potential links between climate variability and cropland productivity. An assessment using crop yield statistical data showed that CFAPAR can be taken as a surrogate of cropland productivity in the Canadian Prairies. Cropland productivity showed an increasing trend in most areas of Canadian Prairies, in general, during the period from 2000 to 2013. Interannual variability in cropland productivity on the Canadian Prairies was influenced positively by rainfall variation and negatively by mean air temperature.

Satellite soil moisture is a critical variable for identifying susceptibility to hydroclimatic risks such as drought, dryness, and excess moisture. Satellite soil moisture data from the Soil Moisture Active/Passive (SMAP) mission was used to evaluate the sensitivity to hydroclimatic risk events in Canada. The SMAP soil moisture data sets in general capture relative moisture trends with the best estimates from the passive-only derived soil moisture and little difference between the data at different spatial resolutions. In general, SMAP data sets overestimated the magnitude of moisture at the wet extremes of wetting events. A soil moisture difference from average (SMDA) was calculated from SMAP and historical Soil Moisture and Ocean Salinity (SMOS) data showed a relatively good delineation of hydroclimatic risk events, although caution must be taken due to the large variability in the data within risk categories. Satellite soil moisture data sets are more sensitive to short term water shortages than longer term water deficits. This was not improved by adding “memory” to satellite soil moisture indices to improve the sensitivity of the data to drought, and there is a large variability in satellite soil moisture values with the same drought severity rating.