Jianxiu Shen

Large-scale rainfed cropping systems (broadacre agriculture) face intensifying climate and resource stresses that undermine yield stability and farm livelihoods. Remote sensing (RS) offers critical tools for improving resilience by monitoring crop performance—productivity, phenology, and environmental stress—across large areas and timeframes. This review aims to synthesize methodological advances over the past two decades in applying RS for broadacre crop monitoring and to identify key challenges and integration opportunities. Peer-reviewed studies across diverse crops and regions were systematically examined to evaluate the strengths, limitations, and emerging trends across the three RS application themes. The review finds that (1) RS enables spatially explicit yield estimation from regional to paddock scales, with vegetation indices (VIs) and phenology-adjusted metrics closely correlated with yield. (2) Time-series analyses of RS data effectively capture phenological transitions critical for forecasting, supported by advances in curve fitting, sensor fusion, and machine learning. (3) Thermal and multispectral indices support the early detection of abiotic (drought, heat, salinity) and biotic (pests, disease) stresses, though specificity remains limited. Across themes, methodological silos and sensor integration barriers hinder holistic application. Emerging approaches, such as multi-sensor/scale fusion, RS–crop model data assimilation, and operational and big data integration, provide promising pathways toward resilience-focused decision support. Future research should define quantifiable resilience metrics, cross-theme predictive integration, and accessible tools to guide climate adaptation.

Large-scale rainfed cropping systems (broadacre agriculture) face intensifying climate and resource stresses that undermine yield stability and farm livelihoods. Remote sensing (RS) offers critical tools for improving resilience by monitoring crop performance—productivity, phenology, and environmental stress—across large areas and timeframes. This review aims to synthesize methodological advances over the past two decades in applying RS for broadacre crop monitoring and to identify key challenges and integration opportunities. Peer-reviewed studies across diverse crops and regions were systematically examined to evaluate the strengths, limitations, and emerging trends across the three RS application themes. The review finds that (1) RS enables spatially explicit yield estimation from regional to paddock scales, with vegetation indices (VIs) and phenology-adjusted metrics closely correlated with yield. (2) Time-series analyses of RS data effectively capture phenological transitions critical for forecasting, supported by advances in curve fitting, sensor fusion, and machine learning. (3) Thermal and multispectral indices support early detection of abiotic (drought, heat, salinity) and biotic (pests, disease) stresses, though specificity remains limited. Across themes, methodological silos and sensor integration barriers hinder holistic application. Emerging approaches—such as multi-sensor/scale fusion, RS–crop model data assimilation, and operational and big data integration—provide promising pathways toward resilience-focused decision support. Future research should define quantifiable resilience metrics and cross-theme predictive integration to guide climate adaptation.

Long-term maps of within-field crop yield can help farmers understand how yield varies in time and space and optimise crop management. This study investigates the use of Landsat NDVI sequences for estimating wheat yields in fields in Western Australia (WA). By fitting statistical crop growth curves, identifying the timing and intensity of phenological events, the best single integrated NDVI metric in any year was used to estimate yield. The hypotheses were that: (1) yield estimation could be improved by incorporating additional information about sowing date or break of season in statistical curve fitting for phenology detection; (2) the integrated NDVI metrics derived from phenology detection can estimate yield with greater accuracy than the observed NDVI values at one or two time points only. We tested the hypotheses using one field (~235 ha) in the WA grain belt for training and another field (~143 ha) for testing. Integrated NDVI metrics were obtained using: (1) traditional curve fitting (SPD); (2) curve fitting that incorporates sowing date information (+SD); and (3) curve fitting that incorporates rainfall-based break of season information (+BOS). Yield estimation accuracy using integrated NDVI metrics was further compared to the results using a scalable crop yield mapper (SCYM) model. We found that: (1) relationships between integrated NDVI metrics using the three curve fitting models and yield varied from year to year; (2) overall, +SD marginally improved yield estimation (r = 0.81, RMSE = 0.56 tonnes/ha compared to r = 0.80, RMSE = 0.61 tonnes/ha using SPD), but +BOS did not show obvious improvement (r = 0.80, RMSE = 0.60 tonnes/ha); (3) use of integrated NDVI metrics was more accurate than SCYM (r = 0.70, RMSE = 0.62 tonnes/ha) on average and had higher spatial and yearly consistency with actual yield than using SCYM model. We conclude that sequences of Landsat NDVI have the potential for estimation of wheat yield variation in fields in WA but they need to be combined with additional sources of data to distinguish different relationships between integrated NDVI metrics and yield in different years and locations.

Seasonal climate is the main driver of crop growth and yield in broadacre grain cropping systems. With a 40-year record of 30 m resolution images and 16-day revisits, the Landsat satellite series is ideal for producing long-term records of remotely sensed phenology to build understanding of how climate affects crop growth. However, the time-series of Landsat images exhibits gaps caused by cloud cover, which is common in wet periods when crops reach maximum growth. We propose a novel spatial–temporal approach to gap-filling that avoids data fusion. Crop growth curve estimation is used to perform temporal smoothing and incorporation of spatial weights allows spatial smoothing. We tested our approach using Landsat NDVI data acquired for an 8000 ha study area in Western Australia using a train/test approach where 157 available Landsat-7 images between 2013 and 2019 were used to train the model, and 95 at least 80% cloud-free Landsat-8 images from the same period were used to test its performance. We found that compared to nonspatial estimation, use of spatial weights in growth curve estimation improved correlation between observed and predicted NDVI by 75%, MAE by 31% and RMSE by 75%. For cropland, the correlation is improved by 58%, the MAE by 36% and the RMSE by 76%. We conclude that spatially weighted estimation of crop growth curves can be used to fill spatial and temporal gaps in Landsat NDVI for the purpose of within-field monitoring. Our approach is also applicable to other data sources and vegetation indices.

Wheat is the most important grain crop in Australia, which plays a significant role in world grain-trading market. However, climate warming, water shortage, as well as more frequent extreme weather events (e.g., heatwaves, droughts and floods), under pressure of food demand, would pose great risks to all aspects of wheat production worldwide, especially in Australia with high climate variability. This study aggregated multi-source observational data by using meteorological statistics, in-situ investigation data and the MODIS Enhanced Vegetation Index (EVI) product to explore and examine the correlation between climate variability and spatial-temporal patterns of wheat phenology metrics and productivity. The results from tests over 370 wheat trial sites showed: 1) narrower and earlier sowing and harvesting windows occurred in a drought year (2006) compared with a normal year (2005). Differences in sowing and harvesting window lengths were 9 and 5 days, respectively; 2) different weather patterns in each agro-climatic zone were followed by different remotely sensed crop EVI seasonality profiles. Crop growth was least affected by climate variability in agro-climate region E2, which is located in the south part of study area. This study reveals new information on cropland-climate relationships across the wheat belt in NSW in a changing climate.