Next Article in Journal
Snow Cover Maps from MODIS Images at 250 m Resolution, Part 2: Validation
Next Article in Special Issue
Using Low Resolution Satellite Imagery for Yield Prediction and Yield Anomaly Detection
Previous Article in Journal
Estimating the Above-Ground Biomass in Miombo Savanna Woodlands (Mozambique, East Africa) Using L-Band Synthetic Aperture Radar Data
Previous Article in Special Issue
Investigating the Relationship between X-Band SAR Data from COSMO-SkyMed Satellite and NDVI for LAI Detection
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Testing the Temporal Ability of Landsat Imagery and Precision Agriculture Technology to Provide High Resolution Historical Estimates of Wheat Yield at the Farm Scale

School of Earth and Environmental Sciences, The University of Adelaide, Adelaide, SA 5005, Australia
*
Author to whom correspondence should be addressed.
Remote Sens. 2013, 5(4), 1549-1567; https://0-doi-org.brum.beds.ac.uk/10.3390/rs5041549
Submission received: 29 January 2013 / Revised: 7 March 2013 / Accepted: 11 March 2013 / Published: 26 March 2013
(This article belongs to the Special Issue Advances in Remote Sensing of Agriculture)

Abstract

:
The long term archiving of both Landsat imagery and wheat yield mapping datasets sensed by precision agriculture technology has the potential through the development of statistical relationships to predict high resolution estimates of wheat yield over large areas for multiple seasons. Quantifying past yield performance over different growing seasons can inform agricultural management decisions ranging from fertilizer applications at the sub-paddock scale to changes in land use at a landscape scale. However, an understanding of the magnitude of prediction errors is needed. In this study, we examine the predictive wheat yield relationships developed from Normalised Difference Vegetation Index (NDVI) acquired Landsat imagery and combine-mounted yield monitors for three Western Australian farms over different growing seasons. We further analysed their predictive capability when these relationships are used to extrapolate yield from one farm to another. Over all seasons, the best predictions were achieved with imagery acquired in September. Of the five seasons reviewed, three showed very reasonable prediction accuracies, with the low and high rainfall years providing good predictions. Medium rainfall years showed the greatest variation in prediction accuracy with marginal to poor predictions resulting from narrow ranges of measured wheat yield and NDVI values. These results demonstrate the potential benefit of fusing together two high resolution datasets to create robust wheat yield prediction models over different growing seasons, the outputs of which can be used to inform agricultural decision making.

1. Introduction

Estimating agricultural production is necessary for a wide variety of applications: early warning of food security, famine and drought [14]; supporting validation of biophysical crop models [5,6]; a substitute for yield maps to analyse yield consistency [7], and to optimise spatially explicit fertiliser application [810]. The collection of reliable and timely information of crop performance can also influence pricing policies, marketing and trading decisions [13,14] and is important for government, growers, insurance and agricultural companies [11,12] so that logistical issues can be anticipated [12].
The application of airborne or satellite remotely sensed imagery is critical for the above objectives since on-ground human observations of crops is limited. Past predictions of yield which incorporated a variety of remote sensing techniques have ranged from simple statistical relationships [11,12,1518]; more advanced relationships built on agronomic and metrological data [4,19,20], to models utilising absorbed photosynthetically active radiation [2124]. The majority of these studies have employed low resolution sensors because of their low cost, availability, extensive spatial coverage and frequent acquisition dates. This choice of imagery is particularly advantageous for broad-scale crop progression, crop canopy emergence, maturation and senescence, but sacrifices higher spatial resolution for greater temporal frequency.
Landsat imagery allows yield predictions at a higher resolution, where the pattern of measurements highlights yield performance differences due to soil type and topographic location and where large variations in yield are evident. Reasonable relationships have been observed between yields collected at the farm, field or geo-referenced hand-sampled scale and spectral vegetation indices [14,18,19,22,25]. However, spatial mismatches occur with the matching of high resolution spectral indices with broad scale farm and field estimates, as well with yield estimation at the plant level which is also labour intensive, expensive to collect [26] and usually limited to a small spatial extent. Precision agriculture and more specifically yield mapping provides an alternative method to collect yield measurements at a matching scale. Yield and spatial position are collected via a combine harvester every one to three seconds, allowing maps of yield to be produced at a high spatial resolution that is similar to the resolution of the spectral indices produced by the Landsat sensor. Several studies have reported a strong relationship between yield and Landsat and IKONOS imagery [10,3032] but most studies are from individual fields and rarely have these relationships been used to predict yields elsewhere.
Models that are able to predict yield patterns based on historic imagery have significant practical uses. Firstly, they could assist growers who have not yet adopted, or those who have recently adopted, precision agriculture technology to more readily access its benefits by extending the length of the dataset that is available for crop management decision making [34,35]. This would allow land managers to assess the spatially explicit nature of productivity, fertiliser requirements and profitability of their fields from a longer time series [36] and therefore leap-frog the time consuming process of archiving annual yield maps generally used in this assessment process. Secondly, the availability of yield estimates, at the resolution and extent of the Landsat imagery, may assist regional environmental policy decision making [37]. Yield variability over time may help to define spatially explicit scenarios such as a comparison between the returns from traditional cropping with those from an alternative, environmentally friendly land use. The major benefit of this analysis is that yield estimates are at a scale that is relevant for management [38]. Information is provided at a scale where growers make decisions about land use allocation, yet at an extent where policy decisions can be made.
The aim of this paper is to understand how wheat yield-NDVI relationships vary in time and space. The majority of past studies are based on a single paddock or farm. Here, we explicitly explore the yield-NDVI relationship across three farms over a period of five years under a range of very different rainfall conditions. The first objective of this study is to evaluate the sensitivity of the regression models created from data from different farms, years and image acquisition dates. The second objective is to evaluate if these models can be reliably extrapolated across space.

2. Study Area

The study area covers 625 km2 within the northern wheat belt of Western Australia (Figure 1), incorporating four neighbouring farms that have collected yield data for the 1998 to 2004 growing seasons.
The region is characterised by a Mediterranean climate, with cool wet winters and hot dry summers. Long term growing season average is around 320 mm, but this varies considerably within and between years. In this environment, water is the major limitation to agricultural plant productivity [39]. The landscape is predominately broad sand plains with very little relief and salty discharges situated in lower areas. Farms in the region are large, typically with over 2,000 ha of cropping area. The broad-acre cropping rotations are dominated by wheat with break crops of lupins, rapeseed (canola) and to a lesser extent barley and oats. Pastures for cattle and sheep grazing are also common, as well as small scattered stands of remnant native vegetation consisting of a mixture of evergreen shrubs and trees. Flowering and grain filling of crops occurs in spring (September) with harvest in late spring and early summer (November–December).

3. Methods

We developed empirical relationships between wheat yield measured on three farms with precision agriculture yield mapping capability and the Normalised Difference Vegetation Index (NDVI) derived from Landsat 5 TM and 7 ETM+ imagery. The strength of these relationships was tested at different times within the season and over different annual rainfall distributions. The predictive power of the models was tested on adjacent farms and using imagery at different times.

3.1. Characterising Annual Growing Seasons by Rainfall Distributions

Wheat production in this study area is extremely dependent on highly fluctuating rainfall. The growing season occurs between March and November and the end of season yield is determined by the plants’ ability to get water during specific phases of development corresponding to germination, vegetative and reproductive growth [40] (Table 1). Annual decisions to sow a wheat crop are determined by the amount of rain falling between March and April and the availability of water held in the soil. This variability in the initial sowing date therefore also leads to variable harvesting dates.
To encapsulate the variations in season break and finish, we characterized each growing season by the total rainfall between March and November. Figure 2 shows that 1999 had the highest annual rainfall while rainfall totals for the other years were fairly consistent, with 2003 having the next highest and the years 1998, 2001 and 2004 the lowest.

3.2. Wheat Phenology and Image Acquisition Date

In agricultural areas the phenological cycle of seeding, growth, maturity and harvesting of managed agricultural vegetation is repeated on an annual basis. Several authors [43,44] suggest that vegetation indices such as NDVI recorded at critical times during the growing season can help characterise the spatial variability in crop performance. Others [13,37] suggest that accurate wheat yield predictions are possible using only one image, provided it is acquired towards the middle of the growing season when most wheat crop canopies are fully developed. Selection of the optimal image acquisition date depends on two factors; firstly, the relationship of NDVI to particular wheat development stages, and secondly the relationship between NDVI at this time and the final wheat grain yield.
Laboratory-based spectral studies of wheat crop development have found that NDVI is very sensitive to leaf area index (LAI) between 0 and 2 [45]. After this stage (during the ripening process) there is an asymptotic relationship between NDVI, LAI and yield [14,46,47]. This indicates that the use of NDVI as a surrogate for biomass and yield prediction is best at crop stages where LAI values are less than 3 [45]. Similarly, measurement of spectral reflectance of wheat growth stages from a tractor mounted radiometer [48] showed similar asymptotic relationships of NDVI and LAI values during this time period. These findings have also been supported by satellite based studies which indicated good associations between NDVI and wheat grain yield 50 to 110 days after plant emergence [2,15,19,26,49,50]. This wide time period encapsulates critical phenological phases for predicting the final wheat yield, including the formation and organisation of the canopy structure after the emergence of the flag leaf [8,50], at the end of stem elongation and beginning of heading [47] and grain formation [51]. Other studies [11,14] and those undertaken in Australia [37,49,52] suggest that for cereal crops, biomass or LAI measures at the flowering or anthesis stage (September to November) are closely related to final yield. These studies show that wheat yield has a positive relationship with the biomass produced by the plant, i.e., low biomass leads to low yield, high biomass high yield.
For this study, cloud free Landsat 5 TM and 7 ETM+ data were acquired (path-row 113/81 and 112/81) between August and October (Table 2). Images were processed to USGS Level 1G and further orthocorrection to the Geocentric Datum of Australia 1994 and systematic radiometric corrections were made by the vendor, Geoscience Australia. Pixel size was re-sampled to 25 m. As no direct comparison was made between the images and the time between the acquisitions was short (usually 16 days) no further radiometric calibration was performed.

3.3. Wheat Yield Mapping Data

Wheat yield data was obtained from four different yield monitoring systems used on four commercial farms for different growing seasons (Table 3). Previous research has concentrated on identifying spatial variability of yield within a single or a limited number of fields. In this study, by contrast, we used wheat yield data from four farms, incorporating around 10 to 16 fields, each over 70 hectares in size. The total area cropped to wheat each year was over 1,500 hectares per farm. This yield data was aggregated by farm and year for the analysis.
The extent of data available was limited by when the yield monitoring technology was adopted by each grower (Table 3). For 2004, the last year of analysis, yield mapping data from Farm 3 was unavailable. However, the yield monitor used in Farm 1 was trialled on another neighbouring farm (Farm 4) and this was used as a replacement for the Farm 3 dataset. This farm had monitored wheat yield over approximately 600 ha.
The raw wheat yield data derived from the combine monitors was cleaned using algorithms that removed erroneous yield values associated with harvester dynamics and operator error. Data was then interpolated to a 25 m by 25 m grid using VESPER kriging software [53] under a defined protocol [54], filling in missing areas due to error removal and provided a grid of similar resolution to that of the Landsat imagery.
The majority of average yields were similar within and across farms in all seasons with the exception of 1998 where Farm 2 recorded nearly a tonne higher yield than the neighbouring Farm 1. This may be the result of calibration errors in the first year of harvest monitor use. Visual interpretation of the histograms of the wheat yields (broken into 0.25 kg/ha intervals) show that the datasets range between positive and negative skewness (see supplementary section for histograms in annual videos).

3.4. Comparison of Wheat Yields to Landsat NDVI

The Normalised Difference Vegetation Index (NDVI = (band4 – band3)/(band4 + band3)) was calculated for all images. Empirical calibration relationships were developed for each farm using interpolated wheat yield estimates and the associated NDVI value for all 25 m by 25 m grid cells. Maps of wheat yield, NDVI and scatter plots of their relationships for each farm are present in Supplementary video section. Second order polynomial regression models were fitted to the relationships (see Appendix 1). The models were forced through the origin to acknowledge that NDVI values close to zero are indicative of an absence of photosynthetic vegetation and hence cannot produce any yield. This approach simplified the models by reducing the number of fitted parameters. Each farm model was used to predict wheat yields from the NDVI values taken from the other farms in the corresponding image date. Predicted yields from the yield-NDVI models were then compared to the actual measured yield on each farm. Two criteria were used to assess the accuracy of the predictions of wheat yield: the Root Mean Square Error (RMSE) and the Nash-Sutcliffe Efficiency Criterion (E). The Nash Sutcliffe Efficiency Criterion (E) is a estimate of model performance that compares the variance around the 1:1 line between predicted (Pi) and observe values (Oi) with the variance of observations around the mean (Ō) (Equation (1)). It is defined as one minus the sum of the absolute squared difference between the predicted and observed values normalised by the variance of the observed values [55]. Values of E range between 1.0 (perfect fit) and −∞. Here, a negative efficiency indicates that the mean value of the observed yield would have been a better predictor than the model [56].
E = 1 i = 1 n ( O i P i ) 2 i = 1 n ( O i O ¯ ) 2

4. Results

4.1 Accuracy of the Regression Relationships and Model Extrapolations over Time

We developed a total of 34 regression models for the three farms over five years using a total of 13 cloud free image dates. These models were used to extrapolate yield on adjacent farms for a total of 64 extrapolation tests (Table 2). Overall, models provided a good fit with an average RMSE of 0.59 t/ha and all were below 0.85 t/ha. Most regression analyses show highly significant relationships between predicted and observed yield. R2 values vary substantially with half the regressions models having a regression coefficient below 0.3 and half above, respectively (Appendix 1). Categorising the models by day of year showed that RMSE was lowest and E was highest around the 250th day of the year (Figure 3). Average RMSE was 0.58 t/ha (Standard Deviation (SD = 0.11) and E was 0.3 (SD = 0.16) in the time period between day 245 and 255. This range would be the optimal time period for selection of imagery for the development of the yield-NDVI relationships.
Extrapolation, as expected, showed a reduced performance with an average RMSE of 0.71 t/ha and a maximum of 2.0 t/ha. Efficiency for the extrapolation models varies around −6.29 (not shown in Figure 4) and 0.49. Extrapolation showed mixed results with only 11 models out of the 64 having an efficiency of above 30%. Both RMSE and E however, also have an optimum time period, once again around the 250th day of the year (Figure 4). Average RMSE was 0.64 t/ha (SD = 0.12) and E was 0.13 (SD = 0.24) in the time period between day 245 and 255.

4.2. Yearly Yield-NDVI Relationships

Model performance for different farms and years was evaluated to determine if these models can be extrapolated in space and time. Videos of the distributions of NDVI over time, the final yield maps and a scatter plot of their relationships are available in the supplementary sections by year. Model and yield prediction efficiency criteria, Root Mean Squared Error (RMSE) and the Nash-Sutcliffe Efficiency Criteria (E) for each year and each image are summarized in the Tables 4, 5, 6 and 7. These criteria highlight the magnitude of error and prediction accuracy when the models are extrapolated. A description of the growing season, NDVI progression and a brief interpretation of the results are also included.

5. Discussion

5.1. Error, Model Accuracy and Uncertainties

The results show a reasonable relationship between NDVI and yield at different growth stages suggesting that the approach is a promising undertaking. An average RMSE of 0.58 t/ha for the models developed in early to mid September (245–255) seems reasonable and the extrapolation result of 0.64 t/ha is not substantially worse. A comparison of RMSE measured within this study compares well with those reported in the literature (Table 9). It is noticeable that yield predictions from the AVHRR or MODIS sensors produces lower prediction errors when compared to the ground-based or high resolution sensors. This indicates that some of the variability is filtered out by smoothing caused by the larger pixel size and the mismatch of data scales used in their validation.
Overall the pattern of yield predicted from NDVI shows moderate prediction accuracy, however, there is still a substantial proportion of variability remaining unexplained and a large number of models that fail to predict yield reliably when applied to other farms. This is mostly due to a bias, where yield is consistently over or under predicted due to variability in the relationships between the yield mapped and NDVI datasets. Although, we have employed an extensive data cleaning process, unbiased noise may in part be due to errors in the yield monitoring (i.e., different calibrations of the flow rates by the yield monitors on different farms). The measurement of NDVI is less error prone than yield mapping but the snapshot approach has its limitations [33] because it provides a static representation of the current variation in plant growth. This variation consists of plants in fields with different sowing dates and represents plant phenological responses within fields to differences in soil-water interactions under water stress [64]. The comparison of the snapshot to final yield measured around three months later (December) also is hazardous since the sensitivity of wheat to adverse meteorological conditions at critical growth stages, principally at flowering and grain filling is high. Conditions may occur where above ground biomass is high which is represented by a high NDVI value in the image, but the actual yield may not be commensurately large [65] due to pest infiltration, disease or if lodging occurs late in the season. While a careful selection of the image acquisition can improve the extrapolation performance of the models, such processes are very difficult to account for in a single image and the improvement of the models is inconsistent.
Given these situations, our results showed that reliable relationships are found in three out of the five seasons and mid-September being the optimal time for the image acquisition. However, annual variations in rainfall and soil-water interactions can influence these relationships and the selection of image dates. The assumption that low yields are related to low NDVI values and that high yields are related to high NDVI values at a time where the majority of the wheat plants are at maturity can be said to hold. This research provides further evidence to Lobell’s [13] suggestion that accurate wheat yield predictions are possible using only one image, provided the image is acquired towards the middle of the growing season when most wheat crop canopies are fully developed.

5.2. Implications in Using the Extrapolated Estimates for Management

Models tended to overestimate yield at low NDVI and underestimate yields at high NDVI values. This has several implications for the use of the datasets in both farm and catchment management situations. For precision agriculture management, depending on the management strategy undertaken, this degree of prediction error may not be as limiting. Growers who prefer a zone management strategy might only require accuracy of the mean yield for a zone [59]. Those who adopt variable rate applications will require greater accuracy in the estimation of yield. Further research is needed into how the selection of the model with different overestimations and underestimations of wheat yield will affect the applicability to economic returns and also precision agriculture management decisions. For catchment management, the overestimation of yield within lower ranges therefore poses empirical limitations when trying to spatially target areas where land use alternatives may be economically comparable to traditional wheat cropping. This study has not taken into account other crops types that are present in wheat farming systems as break crops. These crops may not have the same yield-NDVI relationships due to plant physiology differences. The use of crop rotations may also limit predictions, especially where wheat is not the major crop in the agricultural landscape. The major methodological constraint to our approach is the availability of data. The methodology relies on access to historical yield mapping data and the availability of cloud free images at anthesis within the agricultural region. These issues may pose limitations to this application of this methodology elsewhere.

6. Conclusions

The long term archive of Landsat imagery combined with a growing archive of wheat yield maps shows substantial potential for the prediction of historical wheat yield over large areas at resolution that is needed for farm and environmental management. We have shown the range of prediction accuracies for empirical relationships to predict wheat yield over a variety of seasons using wheat yield data collected by combine-mounted yield monitors and Landsat NDVI. Furthermore, we have tested the sensitivity of predictions to imaging date within the growing season. Overall, imagery acquired in mid-September showed stronger relationships (average R2 of 0.32 and average RMSE of 0.58 t/ha) and prediction accuracies (average RMSE of 0.64 t/ha) when compared to those derived in late August, late September and early October. Of the five seasons reviewed, three gave moderate to good prediction accuracies, while marginal yield predictive capacity was obtained with the 2003 model. Accuracy of the 1998 model was the poorest with prediction capability highlighted by efficiency criteria showing that the average farm wheat yield estimates provided better yield prediction.
The snapshot approach used in this study has the potential to produce high resolution estimates of wheat yield. This approach assumes that biomass measured within the flowering window of a wheat crop through a remote sensing proxy (NDVI) has a positive relationship with wheat yield. These overall results show that fairly robust wheat yield prediction models can be created over differing rainfall seasons. These models can then be applied to satellite imagery to derive spatially varying estimates of wheat yield and economic performance that are both broad in extent and high in resolution. The physical dimensions of this information will help inform decisions by both growers and policy makers on agronomic management, business viability and regional environmental objectives under the threat of a changing climate.

Acknowledgments

The authors would like to thank the CSIRO and Geoscience Australia for access to the remote sensing and yield data used in this paper. Greg Lyle is a PhD candidate supported by a University of Adelaide Faculty of Sciences divisional scholarship.

References

  1. Hutchinson, C.F. Uses of satellite data for famine early warning in sub-Saharan Africa. Int. J. Remote Sens 1991, 12, 1405–1421. [Google Scholar]
  2. Quarmby, N.A.; Milnes, M.; Hindle, T.L.; Silleos, N. The use of multi-temporal NDVI measurements from AVHRR data for crop yield estimation and prediction. Int. J. Remote Sens 1993, 14, 199–210. [Google Scholar]
  3. Liu, W.T.; Kogan, F.N. Monitoring regional drought using the vegetation condition index. Int. J. Remote Sens 1996, 17, 2761–2782. [Google Scholar]
  4. Prasad, A.K.; Singh, R.P.; Tare, V.; Kafatos, M. Use of vegetation index and meteorological parameters for the prediction of crop yield in India. Int. J. Remote Sens 2007, 28, 5207–5235. [Google Scholar]
  5. Moulin, S.; Bondeau, A.; Delecolle, R. Combining agricultural crop models and satellite observations: From field to regional scales. Int. J. Remote Sens 1998, 19, 1021–1036. [Google Scholar]
  6. Basso, B.; Ritchie, J.T.; Pierce, F.J.; Braga, R.P.; Jones, J.W. Spatial validation of crop models for precision agriculture. Agr. Syst 2001, 68, 97–112. [Google Scholar]
  7. Pinter, P.J.; Hatfield, J.L.; Schepers, J.S.; Barnes, E.M.; Moran, M.S.; Daughtry, C.S.T.; Upchurch, D.R. Remote sensing for crop management. Photogramm. Eng. Remote Sensing 2003, 69, 647–664. [Google Scholar]
  8. Raun, W.R.; Solie, J.B.; Johnson, G.V.; Stone, M.L.; Lukina, E.V.; Thomason, W.E.; Schepers, J.S. In-season prediction of potential grain yield in winter wheat using canopy reflectance. Agron. J 2001, 93, 131–138. [Google Scholar]
  9. Flowers, M.; Weisz, R.; Heiniger, R.; Tarleton, B.; Meijer, A. Field validation of a remote sensing technique for early nitrogen application decisions in wheat. Agron. J 2003, 95, 167–176. [Google Scholar]
  10. Reyniers, M.; Vrindts, E. Measuring wheat nitrogen status from space and ground-based platform. Int. J. Remote Sens 2006, 27, 549–567. [Google Scholar]
  11. Vicente-Serrano, S.M.; Cuadrat-Prats, J.M.; Romo, A. Early prediction of crop production using drought indices at different time-scales and remote sensing data: Application in the Ebro Valley (north-east Spain). Int. J. Remote Sens 2006, 27, 511–518. [Google Scholar]
  12. Wall, L.; Larocque, D.; Leger, P.M. The early explanatory power of NDVI in crop yield modelling. Int. J. Remote Sens 2008, 29, 2211–2225. [Google Scholar]
  13. Lobell, D.B.; Asner, G.P.; Ortiz-Monasterio, J.I.; Benning, T.L. Remote sensing of regional crop production in the Yaqui Valley, Mexico: Estimates and uncertainties. Agr. Ecosyst. Environ 2003, 94, 205–220. [Google Scholar]
  14. Liu, L.; Wang, J.; Bao, Y.; Huang, W.; Ma, Z.; Zhao, C. Predicting winter wheat condition, grain yield and protein content using multi-temporal EnviSat-ASAR and Landsat TM satellite images. Int. J. Remote Sens 2006, 27, 737–753. [Google Scholar]
  15. Benedetti, R.; Rossini, P. On the use of NDVI profiles as a tool for agricultural statistics: The case study of wheat yield estimate and forecast in Emilia Romagna. Remote Sens. Environ 1993, 45, 311–326. [Google Scholar]
  16. Sharma, T.; Sudha, K.S.; Ravi, N.; Navalgund, R.R.; Tomar, K.P.; Chakravarty, N.V.K.; Das, D.K. Procedures for wheat yield prediction using Landsat MSS and IRS-1A data. Int. J. Remote Sens 1993, 14, 2509–2518. [Google Scholar]
  17. Sridhar, V.N.; Dadhwal, V.K.; Chaudhari, K.N.; Sharma, R.; Bairagi, G.D.; Sharma, A.K. Wheat production forecasting for a predominantly unirrigated region in Madhya Pradesh, India. Int. J. Remote Sens 1994, 15, 1307–1316. [Google Scholar]
  18. Hamar, D.; Ferencz, C.; Lichtenberger, J.; Tarcsai, G.; Ferencz-Arkos, I. Yield estimation for corn and wheat in the Hungarian Great Plain using Landsat MSS data. Int. J. Remote Sens 1996, 17, 1689–1699. [Google Scholar]
  19. Rudorff, B.F.T.; Batista, G.T. Wheat yield estimation at the farm level using TM Landsat and agro-meteorological data. Int. J. Remote Sens 1991, 12, 2477–2484. [Google Scholar]
  20. Manjunath, K.R.; Potdar, M.B.; Purohit, N.L. Large area operational wheat yield model development and validation based on spectral and meteorological data. Int. J. Remote Sens 2002, 23, 3023–3038. [Google Scholar]
  21. Bastiaanssen, W.G.M.; Ali, S. A new crop yield forecasting model on satellite measurements applied across Indus Basin, Pakistan. Agr. Ecosyst. Environ 2003, 94, 321–340. [Google Scholar]
  22. Lobell, D.B.; Asner, G.P. Comparison of Earth Observing-1 ALI and Landsat ETM+ for crop identification and yield prediction in Mexico. IEEE Trans. Geosci. Remote Sens 2003, 41, 1277–1282. [Google Scholar]
  23. Patel, N.R.; Bhattacharjee, B.; Mohammed, A.J.; Tanupriya, B.; Saha, S.K. Remote sensing of regional yield assessment of wheat in Haryana, India. Int. J. Remote Sens 2006, 27, 4071–4090. [Google Scholar]
  24. Duchemin, B.; Maisongrande, P.; Boulet, G.; Benhadj, I. A simple algorithm for yield estimates: Evaluation for semi-arid irrigated winter wheat monitored with green leaf area index. Environ. Model. Softw 2008, 23, 876–892. [Google Scholar]
  25. Ferencz, C.; Bognar, P.; Lichtenberger, J.; Hamar, D.; Tarcsai, G.; Timar, G.; Molnar, G.; Pasztor, S.Z.; Steinbach, P.; Szekely, B.; et al. Crop yield estimation by satellite remote sensing. Int. J. Remote Sens 2004, 25, 4113–4149. [Google Scholar]
  26. Doraiswamy, P.C.; Hatfield, J.L.; Jackson, T.J.; Akhmedov, B.; Prueger, J.; Stern, A. Crop condition and yield simulations using Landsat and MODIS. Remote Sens. Environ 2004, 92, 548–559. [Google Scholar]
  27. Reitz, P.; Kutzbach, H.D. Investigations on a particular yield mapping system for combine harvesters. Comput. Electron. Agr 1996, 14, 137–150. [Google Scholar]
  28. Arslan, S.; Colvin, T.S. An evaluation of the response of yield monitors and combines to varying yields. Precis. Agr 2002, 3, 107–122. [Google Scholar]
  29. Lark, R.M.; Stafford, J.V.; Bolam, H.C. Limitations on the spatial resolution of yield mapping for combinable crops. J. Agr. Eng. Res 1997, 66, 183–193. [Google Scholar]
  30. Thenkabail, P.S. Biophysical and yield information for precision farming from near-real-time and historical Landsat TM images. Int. J. Remote Sens 2003, 24, 2879–2904. [Google Scholar]
  31. Dobermann, A.; Ping, J.L. Geostatistical integration of yield monitor data and remote sensing improves yield maps. Agron. J 2004, 96, 285–297. [Google Scholar]
  32. Enclona, E.A.; Thenkabail, P.S.; Celis, D.; Diekmann, J. Within-field wheat yield prediction from IKONOS data: A new matrix approach. Int. J. Remote Sens 2004, 25, 377–388. [Google Scholar]
  33. Reeves, M.C.; Zhao, M.; Running, S.W. Usefulness and limits of MODIS GPP for estimating wheat yield. Int. J. Remote Sens 2005, 26, 1403–1421. [Google Scholar]
  34. Bullock, D.S.; Bullock, D.G. From agronomic research to farm management guidelines: A primer on the economics of information and precision technology. Precis. Agr 2000, 2, 71–101. [Google Scholar]
  35. Robertson, M.J.; Lyle, G.; Bowden, J.W. Within-field variability of wheat yield and economic implications for spatially variable nutrient management. Field Crop. Res 2008, 105, 211–220. [Google Scholar]
  36. Massey, R.E.; Myers, D.B.; Kitchen, N.R.; Sudduth, K.A. Profitability maps as an input for site-specific management decision making. Agron. J 2008, 100, 52–59. [Google Scholar]
  37. Lyle, G.; Ostendorf, B. A high resolution broad scale spatial indicator of grain growing profitability for natural resource planning. Ecol. Indic 2011, 11, 209–218. [Google Scholar]
  38. Ostendorf, B. Overview: Spatial information and indicators for sustainable management of natural resources. Ecol. Indic 2011, 11, 97–102. [Google Scholar]
  39. Turner, N.C.; Asseng, S. Productivity, sustainability, and rainfall use efficiency in Australian rain-fed Mediterranean agricultural systems. Aust. J. Agr. Res 2005, 56, 1123–1136. [Google Scholar]
  40. Sadras, V.O.; Rodriguez, D. The limit to wheat water-use efficiency in eastern Australia. II. Influences of rainfall patterns. Aust. J. Agr. Res 2007, 58, 657–669. [Google Scholar]
  41. Passioura, J. Increasing crop productivity when water is scarce—From breeding to field management. Agr. Water Manag 2006, 80, 176–196. [Google Scholar]
  42. Anderson, W.K.; Garlinge, J.R. The Wheat Book—Principles and Practice; Department of Agriculture and Food: South Perth, WA, Australia, 2000; p. 322. [Google Scholar]
  43. Alexandridis, T.K.; Gitas, I.Z.; Silleos, N.G. An estimation of the optimum temporal resolution for monitoring vegetation condition on a nationwide scale using MODIS/Terra data. Int. J. Remote Sens 2008, 29, 3589–3607. [Google Scholar]
  44. Carlson, T.N.; Ripley, D.A. On the relation between NDVI, fractional vegetation cover, and leaf area index. Remote Sens. Environ. 1997, 62, 241–252. [Google Scholar]
  45. Aparicio, N.; Villegas, D.; Araus, J.L.; Casadesus, L.; Royo, C. Spectral vegetation indices as nondestructive tools for determining durum wheat yield. Agron. J 2000, 92, 83–91. [Google Scholar]
  46. Aase, J.K.; Siddoway, F.H. Assessing winter wheat dry matter production via spectral reflectance measurements. Remote Sens. Environ 1981, 11, 267–277. [Google Scholar]
  47. Aparicio, N.; Villegas, D.; Araus, J.L.; Casadesus, J.; Royo, C. Relationship between growth traits and spectral vegetation indices in durum wheat. Crop Sci 2002, 42, 1547–1555. [Google Scholar]
  48. Scotford, I.M.; Miller, P.C.H. Combination of spectral reflectance and ultrasonic sensing to monitor the growth of winter wheat. Biosys. Eng 2004, 87, 27–38. [Google Scholar]
  49. Smith, R.C.G.; Adams, J.; Stephens, D.J.; Hick, P.T. Forecasting wheat yield in a Mediterranean-type environment from the NOAA Satellite. Aus. J. Agr. Res 1995, 46, 113–125. [Google Scholar]
  50. Beeri, O.; Peled, A. Spectral indices for precise agriculture monitoring. Int. J. Remote Sens 2006, 27, 2039–2047. [Google Scholar]
  51. Boissard, P.; Pointel, J.G. Reflectance, green leaf area index and ear hydric status of wheat from anthesis until maturity. Int. J. Remote Sens 1993, 14, 2713–2729. [Google Scholar]
  52. Dawbin, K.W.; Evans, J.C.; Duggin, M.J.; Leggett, E.K. Classification of wheat areas and prediction of yields in north-western New South Wales by repetitive Landsat data. Aus. J. Agr. Res 1980, 31, 449–53. [Google Scholar]
  53. Minasny, B.; McBratney, A.B.; Whelan, B.M. VESPER Version1.62—Variogram Estimation and Spatial Prediction plus Error; Australian Centre for Precision Agriculture: Sydney, NSW, Australia, 2005; Available online: http://sydney.edu.au/agriculture/pal/software/vesper.shtml (accessed on 25 October 2012).
  54. Taylor, J.A.; McBratney, A.B.; Whelan, B.M. Establishing management classes for broadacre agricultural production. Agron. J 2007, 99, 1366–1376. [Google Scholar]
  55. Nash, J.E.; Sutcliffe, J.V. River flow forecasting through conceptual models. Part 1-A discussion of principles. J. Hydrol 1970, 10, 282–290. [Google Scholar]
  56. Krause, P.; Boyle, D.P.; Bäse, F. Comparison of different efficiency criteria for hydrological model assessment. Adv. Geosci 2005, 5, 89–97. [Google Scholar]
  57. Ludwig, F.; Asseng, S. Climate change impacts on wheat production in a Mediterranean environment in Western Australia. Agr. Syst 2006, 90, 159–179. [Google Scholar]
  58. Luo, Q.Y.; Bellotti, W.; Williams, M.; Bryan, B. Potential impact of climate change on wheat yield in South Australia. Agr. Forest Meteorol 2005, 132, 273–285. [Google Scholar]
  59. Balaghi, R.; Tychon, B.; Eerens, H.; Jlibene, M. Empirical regression models using NDVI, rainfall and temperature data for the early prediction of wheat grain yields in Morocco. Int. J. Appl. Earth Obs 2008, 10, 438–452. [Google Scholar]
  60. Reyniers, M.; Vrindts, E. Measuring wheat nitrogen status from space and ground-based platform. Int. J. Remote Sens 2006, 27, 549–567. [Google Scholar]
  61. Ren, J.; Chen, Z.; Zhou, Q.; Tang, H. Regional yield estimation for winter wheat with MODIS-NDVI data in Shandong, China. Int. J. Appl. Earth Obs 2008, 10, 403–413. [Google Scholar]
  62. Reyniers, M.; Vrindts, E.; De Baerdemaeker, J. Comparison of an aerial-based system and an on the ground continuous measuring device to predict yield of winter wheat. Eur. J. Agron 2006, 24, 87–94. [Google Scholar]
  63. Reyniers, M.; Vrindts, E. Measuring wheat nitrogen status from space and ground-based platform. Int. J. Remote Sens 2006, 27, 549–567. [Google Scholar]
  64. McMaster, G.S.; Wilhelm, W.W.; Frank, A.B. Developmental sequences for simulating crop phenology for water-limiting conditions. Aus. J. Agr. Res 2005, 56, 1277–1288. [Google Scholar]
  65. Gomez-Macpherson, H.; Richards, R.A. Effect of sowing time on yield and agronomic characteristics of wheat in south-eastern Australia. Aus. J. Agr. Res 1995, 46, 1381–1399. [Google Scholar]

Appendix

Appendix 1. Wheat yield and NDIV regression equations and strength (R2) by date for the farm models.
Appendix 1. Wheat yield and NDIV regression equations and strength (R2) by date for the farm models.
Farm ModelDateEquationGoodness of Fit (R2)
Farm 126/08/98y = −3.0628x2 + 6.0975x0.16
Farm 326/08/98y = −4.2896x2 + 8.0742x0.05
Farm 111/09/98y = −1.1425x2 + 4.6645x0.11
Farm 311/09/98y = 0.6284x2 + 4.4938x0.08
Farm 121/08/99y = 1.8846x2 + 3.8119x0.45
Farm 221/08/99y = −3.0846x2 + 6.2372x0.28
Farm 321/08/99y = −3.5731x2 + 7.1389x0.27
Farm 106/09/99y = 3.2457x2 + 2.9875x0.48
Farm 206/09/99y = −3.9921x2 + 6.9399x0.37
Farm 306/09/99y = 1.1553x2 + 4.8376x0.34
Farm 129/09/99y = −0.129x2 + 5.7751x0.45
Farm 229/09/99y = −9.3139x2 + 10.22x0.12
Farm 329/09/99y = −4.3096x2 + 8.1084x0.28
Farm 115/10/99y = −16.144x2 + 14.513x0.27
Farm 215/10/99y = −25.798x2 + 18.077x−0.11
Farm 315/10/99y = −22.382x2 + 16.51x0.11
Farm 126/08/01y = 3.483x2 + 2.6994x0.38
Farm 226/08/01y = 1.2611x2 + 3.6751x0.31
Farm 326/08/01y = −0.3573x2 + 5.1796x0.11
Farm 111/09/01y = −4.4754x2 + 9.4465x0.50
Farm 211/09/01y = −7.7556x2 + 10.426x0.44
Farm 311/09/01y = −9.3026x2 + 11.922x0.37
Farm 116/09/03y = −0.0958x2 + 4.0322x0.10
Farm 216/09/03y = 0.2875x2 + 4.5415x0.48
Farm 316/09/03y = −2.2006x2 + 5.6679x0.08
Farm 102/10/03y = −0.7643x2 + 5.132x0.24
Farm 202/10/03y = −3.3642x2 + 8.5108x0.56
Farm 302/10/03y = −10.129x2 + 10.603x−0.17
Farm 111/09/04y = 2.9902x2 + 2.452x0.52
Farm 211/09/04y = 3.2813x2 + 2.9058x0.40
Farm 4 (substitute for Farm 3)11/09/04y = −1.2235x2 + 5.1942x0.18
Farm 127/09/04y = 2.7853x2 + 3.6798x0.55
Farm 227/09/04y = 1.5456x2 + 6.2311x0.64
Farm 4 (substitute for Farm 3)27/09/04y = 1.166x2 + 5.1651x0.29
Figure 1. The four farm study area in Western Australia.
Figure 1. The four farm study area in Western Australia.
Remotesensing 05 01549f1
Figure 2. Monthly and long term average monthly rainfall for the five years of analysis.
Figure 2. Monthly and long term average monthly rainfall for the five years of analysis.
Remotesensing 05 01549f2
Figure 3. Root Mean Square Error (RMSE) and the Nash-Sutcliffe Efficiency Criteria (E) values by day of the year for the regression models.
Figure 3. Root Mean Square Error (RMSE) and the Nash-Sutcliffe Efficiency Criteria (E) values by day of the year for the regression models.
Remotesensing 05 01549f3
Figure 4. Root Mean Square Error (RMSE) and the Nash-Sutcliffe Efficiency Criteria (E) values by day of the year for the extrapolation of the regression models.
Figure 4. Root Mean Square Error (RMSE) and the Nash-Sutcliffe Efficiency Criteria (E) values by day of the year for the extrapolation of the regression models.
Remotesensing 05 01549f4
Table 1. The role of monthly rainfall in the growth of wheat [3942].
Table 1. The role of monthly rainfall in the growth of wheat [3942].
MonthsRole of Rainfall in Crop Development Stages
March to May1. Pre-sowing soil water accumulates; aids seed germination.
June2. Post sowing soil water facilitates seed germination, emergence, formation of leaves and start of tillering.
July to August3. Rainfall used by plant to enlarge leaf area; complete closure of the crop canopy achieved.
August4. Leaf development completed; rain feeds most rapid plant growth phase from start of stem elongation and ear emergence.
September5. Optimal flowering period for wheat crops; adequate rainfall critical. Early-flowering plants may have insufficient biomass to set and fill seeds; flowers prone to frost damage. Late-flowering plants may have insufficient soil water for grain filling.
October6. Rainfall contributes to wheat grain filling.
Table 2. Landsat Thematic Mapper images acquired for the five years of analysis.
Table 2. Landsat Thematic Mapper images acquired for the five years of analysis.
August (Sensor)September (Sensor)October (Sensor)
26 Aug. 1998(5)11 Sept. 1998 (5)
21 Aug. 1999 (7)6, 29 Sept. 1999 (7)15 Oct. 1999(7)
11 Sept. 2001(7)
26 Aug. 2001 (7)16 Sept. 2003(5)13 Oct. 2001
11, 27 Sept. 2004(5)2 Oct. 2003(5)
Table 3. Availability of yield mapping data by farm and average farm wheat yield (t/ha), the number of image dates, number of regression models and extrapolation tests.
Table 3. Availability of yield mapping data by farm and average farm wheat yield (t/ha), the number of image dates, number of regression models and extrapolation tests.
Farm19981999200120032004
Farm 12.162.122.212.552.20
Farm 22.332.952.672.45
Farm 33.122.072.432.38
Farm 4 (Substitute for Farm 3)2.59
# image dates24322
# Regressions models412666
# Extrapolation tests424121212
Table 4. Model and yield prediction efficiency criteria for 1998 Landsat imagery-Root Mean Square Error (RMSE) and the Nash-Sutcliffe Efficiency Criteria (E). Values in bold italics represent efficiency criteria for the calibration models. Descriptions of 1998 Growing season rainfall, NDVI progression and interpretations of results.
Table 4. Model and yield prediction efficiency criteria for 1998 Landsat imagery-Root Mean Square Error (RMSE) and the Nash-Sutcliffe Efficiency Criteria (E). Values in bold italics represent efficiency criteria for the calibration models. Descriptions of 1998 Growing season rainfall, NDVI progression and interpretations of results.
ModelDateFarm 1Farm 2Farm 3

RMSEERMSEERMSEEGrowing Season Rainfall (GSR)NDVI Progression
Farm 126/08/980.660.16n/an/a1.11−0.62321 mm; 25% in May; 74% between April and July; 20% in August (45 mm) and September (36 mm); little rainfall in October.NDVI increasing between dates indicating increase in biomass; scattering reduced between August and mid-September.
Farm 326/08/980.93−0.67n/an/a0.850.05

Farm 111/09/980.670.11n/an/a1.04−0.43
Farm 311/09/980.81−0.27n/an/a0.840.08
Interpretation: Wet conditions in May through to August suggest cropping area was not in water deficit therefore adequate soil water for crop development. Clustering of mid-September NDVI suggest flowering occurred later with potential reduction in wheat yields due to higher temperatures and water deficit later in the season. Acquisition of imagery in late August or mid-September showed poor regression and extrapolation relationships and were not optimal for model development. Imagery acquired further into the season could have provided a better relationship if it were available.
Table 5. Model and yield prediction efficiency criteria for 1999 Landsat imagery-Root Mean Square Error (RMSE) and the Nash-Sutcliffe Efficiency Criteria (E). Values in bold italics represent efficiency criteria for the calibration models. Descriptions of 1999 Growing season rainfall, NDVI progression and interpretations of results.
Table 5. Model and yield prediction efficiency criteria for 1999 Landsat imagery-Root Mean Square Error (RMSE) and the Nash-Sutcliffe Efficiency Criteria (E). Values in bold italics represent efficiency criteria for the calibration models. Descriptions of 1999 Growing season rainfall, NDVI progression and interpretations of results.
ModelDateFarm 1Farm 2Farm 3

RMSEERMSEERMSEEGrowing Season Rainfall (GSR)NDVI Progression
Farm 121/08/990.640.450.630.190.75−0.03Extremely high rainfall year with 628 mm; high distribution of rainfall over the months with 105 mm falling in March, 169 mm falling in May and 96 mm in June; September and October rainfall was >35 mm for each month.NDVI increased and distribution narrowed from August to September; comparison of September images showed slight reductions in the magnitudes and an increase in the dispersal; reductions of NDVI values in October.
Farm 221/08/990.550.390.600.280.700.11
Farm 321/08/990.640.170.670.110.640.51

Farm 106/09/990.530.480.690.040.720.05
Farm 206/09/990.630.280.560.370.640.27
Farm 306/09/990.660.200.620.240.610.34

Farm 129/09/990.530.450.77−0.210.690.14
Farm 229/09/990.670.120.660.120.650.22
Farm 329/09/990.590.300.690.040.630.28

Farm 115/10/990.610.290.86−0.492.00−6.29
Farm 215/10/990.680.090.730.110.720.06
Farm 315/10/990.620.220.74−0.150.700.11
Interpretation: Rainfall was not a constraining factor. Similar patterns of positive yield-NDVI distributions across all image dates and regression models. Their strength declined over time. Validation of the models developed in August and September showed good and consistent predictive capacity over different farm models and farm data. Models created in October provided poorer prediction accuracy.
Table 6. Model and yield prediction efficiency criteria for 2001 Landsat imagery-Root Mean Square Error (RMSE) and the Nash-Sutcliffe Efficiency Criteria (E). Values in bold italics represent efficiency criteria for the calibration models. Descriptions of 2001 Growing season rainfall, NDVI progression and interpretations of results.
Table 6. Model and yield prediction efficiency criteria for 2001 Landsat imagery-Root Mean Square Error (RMSE) and the Nash-Sutcliffe Efficiency Criteria (E). Values in bold italics represent efficiency criteria for the calibration models. Descriptions of 2001 Growing season rainfall, NDVI progression and interpretations of results.
ModelDateFarm 1Farm 2Farm 3

RMSEERMSEERMSEEGrowing Season Rainfall (GSR)NDVI Progression
Farm 126/08/010.560.380.580.230.690.07Lowest rainfall year with 300 mm; 30 mm in March; no rainfall in April and <10 mm in June; 158 mm fell between July, August and September; around 20 mm fell in October.Narrow range of August NDVI values; declining distribution of NDVI in September indicating the August image was taken at crop maturity; earlier than in other years of analysis.
Farm 226/08/010.570.350.550.310.75−0.09
Farm 326/08/010.620.230.640.070.680.11

Farm 111/09/010.500.350.510.470.620.24
Farm 211/09/010.510.490.500.440.650.16
Farm 311/09/010.580.340.580.230.570.37
Interpretation: Limited rainfall will impact on establishment of the crop and initial growth due to prolonged period of water stress. Poor to moderate yield-NDVI relationships with similar prediction capability. October efficiency criteria were not included due to their low results but visualization of the imagery in conjunction with September’s image may provide an additional benefit. For example, NDVI pattern within fields of Farm1 show differences in soil constraints or soil-water interactions given the gradual terrain of the fields. Comparison of the NDVI distribution within specific fields highlights areas which yield better and poorly under water stress. Interestingly, areas of low NDVI do not necessarily correspond to low yield, for example, around the boundary of two fields which have moderate end of season yield values. These patterns are similar on all three farms. This information can be potentially beneficial in future farm management and land use planning because it is expected that Mediterranean farming areas will become drier under climate change [57,58].
Table 7. Model and yield prediction efficiency criteria for 2003 Landsat imagery-Root Mean Square Error (RMSE) and the Nash-Sutcliffe Efficiency Criteria (E). Values in bold italics represent efficiency criteria for the calibration models. Descriptions of 2003 Growing season rainfall, NDVI progression and interpretations of results.
Table 7. Model and yield prediction efficiency criteria for 2003 Landsat imagery-Root Mean Square Error (RMSE) and the Nash-Sutcliffe Efficiency Criteria (E). Values in bold italics represent efficiency criteria for the calibration models. Descriptions of 2003 Growing season rainfall, NDVI progression and interpretations of results.
ModelDateFarm 1Farm 2Farm 3

RMSEERMSEERMSEEGrowing Season Rainfall (GSR)NDVI Progression
Farm 116/09/030.640.100.570.0970.69−0.009372 mm; March and April received 20% of GSR; >45 mm fell in each of the months May, June, July and September; 100 mm fell in August; May to August rainfall represented 80% of GSR. .September NDVI distributions similar for Farms 1 and 3; similar to 1998 results; NDVI more distributed for Farm 2; broader October distributions with high but declining NDVI values.
Farm 216/09/030.80−0.430.430.480.73−0.13
Farm 316/09/030.660.030.460.410.660.08

Farm 102/10/030.580.260.90−1.300.92−0.79
Farm 202/10/031.2−2.220.390.561.11−1.62
Farm 302/10/030.71−0.130.460.410.74−0.17
Interpretation: Significant rainfall fell between May and September. The October image for Farms 1 and 2 provided the best models for yield prediction while the relationship developed in September was the best for Farm 3. Farm model differences exist with yield-NDVI relationship for Farm 2 showing a high positive relationship in comparison to the other models. Consequently, this model performed poorly in September and October when the relationship was validated across the other farm data. In September, the models developed for Farm 3 performed reasonably well on Farms 1 and 2. October imagery shows areas of high NDVI values and corresponding high yield values are still across Farms 1 and 3. This may indicate an extension of crop growth as a result of higher August/September rainfall, markedly different sowing dates or a mislabeling of the crop type (not wheat). In some areas high NDVI values that correspond to low yield values are also seen. This could reflect crop management (e.g., weeds or pests) or weather (e.g., high temperatures related to the later flowering date) impacting on yield later in the season.
Table 8. Model and yield prediction efficiency criteria for 2004 Landsat imagery-Root Mean Square Error (RMSE) and the Nash-Sutcliffe Efficiency Criteria (E). Values in bold italics represent efficiency criteria for the calibration models. Descriptions of 2004 Growing season rainfall, NDVI progression and interpretations of results.
Table 8. Model and yield prediction efficiency criteria for 2004 Landsat imagery-Root Mean Square Error (RMSE) and the Nash-Sutcliffe Efficiency Criteria (E). Values in bold italics represent efficiency criteria for the calibration models. Descriptions of 2004 Growing season rainfall, NDVI progression and interpretations of results.
ModelDateFarm 1Farm 2Farm 4

RMSEERMSEERMSEEGrowing season rainfall (GSR)NDVI progression
Farm 111/09/040.420.470.650.200.650.09327 mm; no rainfall fell in March and only 10 mm in April; nearly 50% fell in May and June; 90% falling before September; less than 10 mm fell in October.Reductions in magnitudes and dispersal of NDVI values between the two September images.
Farm 211/09/040.520.180.560.400.680.01
Farm 411/09/040.480.310.600.330.620.18

Farm 127/09/040.390.550.540.460.670.05
Farm 227/09/040.97−1.840.440.640.81−0.40
Farm 427/09/040.510.220.630.260.580.29
Interpretation: Yield-NDVI relationships follow a positive trend across all three farms and their strength increased with acquisition date. Validation of the models showed that the image acquired in mid-September was the best yield predictor for that year while Farm 3 was the best predictor across all acquisition dates. The Farm 1 yield and NDVI map has an interesting outlier with higher NDVI values for one field. This may be the result of a very late sowing date or mislabeling of the crop type (not wheat).
Table 9. Sensor and prediction error associated with wheat yield prediction.
Table 9. Sensor and prediction error associated with wheat yield prediction.
SensorPrediction Error (t/ha)Author
AVHRR0.14–0.2Smith [48]
AVHRR0.08–0.76Balaghi et al.[59]
AwiFS0.25–0.35Patel et al.[23]
Landsat0.37–0.44Rudorff and Batista [19]
IKONOS0.53Reyniers and Vrindts [60]
MODIS0.23Ren et al.[61]
Multi-spectral radiometer:
Hand held0.5–0.9Duchemin et al.[24]
Ground based platform0.47Reyniers et al.[62]
Aerial based platform0.51Reyniers and Vrindts [63]

Share and Cite

MDPI and ACS Style

Lyle, G.; Lewis, M.; Ostendorf, B. Testing the Temporal Ability of Landsat Imagery and Precision Agriculture Technology to Provide High Resolution Historical Estimates of Wheat Yield at the Farm Scale. Remote Sens. 2013, 5, 1549-1567. https://0-doi-org.brum.beds.ac.uk/10.3390/rs5041549

AMA Style

Lyle G, Lewis M, Ostendorf B. Testing the Temporal Ability of Landsat Imagery and Precision Agriculture Technology to Provide High Resolution Historical Estimates of Wheat Yield at the Farm Scale. Remote Sensing. 2013; 5(4):1549-1567. https://0-doi-org.brum.beds.ac.uk/10.3390/rs5041549

Chicago/Turabian Style

Lyle, Greg, Megan Lewis, and Bertram Ostendorf. 2013. "Testing the Temporal Ability of Landsat Imagery and Precision Agriculture Technology to Provide High Resolution Historical Estimates of Wheat Yield at the Farm Scale" Remote Sensing 5, no. 4: 1549-1567. https://0-doi-org.brum.beds.ac.uk/10.3390/rs5041549

Article Metrics

Back to TopTop