Precision agriculture (PA) recognizes and manages intra-field spatial variability with the desired outcome of increasing profitability and reduced environmental impact  . Site Specific Crop Management (SSCM) is a form of PA that focuses on managing Spatial crop yield variability through matching agricultural inputs with the site potential. SSCM subdivides a cropping field into uniformly manageable zones, based on quantitative measurement of yield limiting factors    .
Spatial variability in crop yield is the result of complex interaction of factors influencing crop growth that include soil (such as nutrients, soil water availability), topographical factors (such as elevation) and climatic factors (such as rainfall and temperature)  . In addition to their spatial variability, some of these factors such as climate have temporal variability, which causes the spatial pattern of crop yield to vary from season to season. The optimal choice of PA over uniform farm management requires a sound understanding of such temporal variability, and is well framed by Whelan et al.  as a “null hypothesis” for precision agriculture i.e. ‘‘Given the large temporal variation evident in crop yield relative to the scale of a single field, then the optimal risk aversion strategy is uniform management’’. The degree of variability, and whether PA can be technically and economically beneficial to manage the variability, are the most important issues to be considered  . Hence, SSCM requires a comprehensive understanding of both spatial and temporal variability in crop growth.
In Mediterranean environments like South Australia, yield variability of rain-fed crops is often controlled by soil water availability. The soil property Plant Available Water-holding Capacity (PAWC) explains a high degree of intra-field spatial variability    . Hence, it is suggested as a basis for management zone delineation  . However, the impact of PAWC on crop yield is highly determined by the seasonal rainfall. For example, high PAWC soils in dry seasons may have small differences in yield compared with low PAWC soils, as they rarely filled to their capacity. On the other hand in years with good opening rains, but with decreasing water availability during the growing season, large differences in yield may occur between high and low PAWC soils. In such way, soil PAWC interacts with the seasonal rainfall and controls the spatial and temporal crop yield variability.
While delineation of agricultural management zones vary in terms of the information used, it is generally based on soil and other yield determining factors  , with crop yield data often the primary source of information. However, interpreting the temporally variable pattern of multi-year yield maps and making use of them for PA purpose is a challenging task for farmers. Researchers have developed a number of different approaches to analyse multiple year yield data. The use of spatial and temporal variograms, where the semi variance is a function of spatial lag and temporal lag   , and development of spatial trend and temporal stability maps   are two examples of approaches for interpretation of multiyear yield data. However, the availability of multiple year yield data, representing variable climatic conditions, has also been identified as a limitation    . As an alternative approach, a few studies have used remote sensing technology to observe spatial crop yield variability   .
Remote sensing vegetation indices have potential to assess crop growth variability by quantifying relative growth and health condition of the crop. The Normalized Difference Vegetation Index (NDVI) is one of the most widely used indices to quantify vegetation vigour from the spectral reflectance of vegetation. NDVI has been used for many crop-monitoring applications     . In precision agriculture, NDVI has been used as a surrogate for crop yield for SSCM zone delineation    . While single images are useful for yield estimation, the inter annual growth variability resulting from soil-climate interaction can produce spurious results depending on the image date selection  . Multi-temporal NDVI, on the other hand, provides an additional temporal dimension to uncover the vegetation dynamics. Remote sensing phenology estimates phenological growth stages including the start of season and end of season from multi-temporal vegetation index data   . The derived metrics may not necessarily correspond directly to conventional, ground-based phenological events, but they provide important information about the vegetation growth dynamics that can be associated with environmental factors such as soil properties [e.g.    ].
The aim of this paper is to explore the potential of remote sensing derived phenological metrics for assessing spatio-temporal variability in cropping fields and to understand the soil-climate interactions that strongly influence crop yield. The specific objective of the project was to examine the relationship between spatial and temporal trend of phenological metrics and predefined management zones, in a South Australian cropping field where intra-field variability is strongly attributed to soil PAWC. The approach presented in this paper provides a pathway for future studies in utilizing phenological metrics for management zone delineation.
2. Materials and Methods
2.1. Study Area
The study was conducted in the rain-fed cropping region of South Australia. The region experiences a Mediterranean climate, with hot summers and wet winters with average annual rainfall of approximately 325 mm and average growing season (April-October) rainfall of approximately 241 mm  . In this region, sowing starts in late March to May, following sufficient rainfall for seeding. Following seeding, the crop germinates, grows and progressively increases in cover to reach peak greenness in September. The crops ripen and reach senescence in October and harvested in November.
The study site was a field close to Minnipa town, in Upper Eyre Peninsula (Figure 1). The field is approximately 65 ha in area, predominantly cultivated for wheat interspersed with some years of barley and pastures. This field was chosen for its soil types, which are representative of those across the wider region. It has also been studied as a focus site for Eyre Peninsula Agricultural Research Foundation (EPARF), where numbers of agricultural research and developmental trials have been undertaken  . The field is characterized by sandy loam to sandy clay loam soils, with the land zones of sandy rises reported to perform well in dry years, and shallow flats, which rarely perform well regardless of crop or pasture choice  .
Derivation of crop phenological information from remote sensing imagery requires a high temporal frequency of images. Moderate Resolution Imaging Spectroradiometer (MODIS) is one of the widely used sensors for phenological
Figure 1. Location map of the study area on Eyre Peninsula, South Australia [Source:  ].
studies. For this study, we used the 16 days composite NDVI image product (MOD13Q1), which has 250 m spatial resolution. The MOD13Q1 product is derived from daily NDVI data using the Constrained View Angle Maximum Value Composite algorithm which extracts the maximum NDVI value for each pixel within the 16 day interval to create a cloud-free composite image    . For this study, a total of 345 images between 01/01/2001 and 18/12/2015 were downloaded from NASA Land Processes Distributed Active Archive Center (LP DAAC) website.
Overlaying the farm boundary on the MODIS NDVI image, using ArcGIS 10.2.1 software  , the area intersects 17 pixels to varying degrees. In order to minimize signal contamination only the six pixels with more than 75% of the pixel area lying within the field boundary were considered (Figure 2(b)).
2.2.2. Management Zone Map
The farm at Minnipa has been a focus site for low rainfall cropping research  . Management zones in this field are well understood and have been delineated on the basis of a long record of scientific observations. The field has been subdivided into three management zones: good, medium and poor, using historic yield data, soil Electro Magnetic survey (EM38) and elevation maps  (Figure 2(a)). The good zone is characterized by red light sandy clay loam with the maximum rooting depth of 80 cm and PAWC of 108 mm. The soil type in the medium zone is red loam with a constrained maximum rooting depth of 60 cm and PAWC of 74 mm. The poor zone, in contrast, has only 40 cm maximum rooting depth with soil type of red sandy clay loam (Minnipa heavy) and PAWC of 57 mm   .
Intersecting the MODIS pixel grids with the management zone map, the pixels were assigned continuous management zone values between 1 (poor zone) and 3 (good zone), based on their relative coverage of the zones. For example,
Figure 2. Minnipa farm (a) management zone map  (b) the 6 selected pixels overlaying the management zone map.
a pixel covered 50% by good and 50% by medium zones was assigned the zone value of 2.5, (50/100 * 3 + 50/100 * 2). The resulting management zone values for the pixels range from 2.4 to 1.43. These values may not represent quantitative measure, rather they are intended to provide ordinal ranking of relative performance between the pixels. Figure 2(b) shows the selected six pixels. The management ranks for these pixels are presented in Table 1.
2.3.1. Overview of the Approach
Figure 3 summarizes the steps followed in our analysis. Firstly, the phenological
Table 1. The selected MODIS pixels with the proportion of good, medium and poor zones and their calculated zone values.
Figure 3. Conceptual workflow of the analysis.
metrics were derived from the MODIS NDVI data, their standard scores were calculated, and then the spatial trend and temporal variance of the metrics were calculated. Finally, the correlation between the management zone values and the spatial trend (temporal mean) and temporal variability (temporal variance) of each of the metrics were assessed.
2.3.2. Extraction of Phenological Metrics
The phenological metrics were extracted using “CropPhenology” package  , in the R software environment  . Figure 4 shows representation of the phenological metrics on the NDVI growth dynamic curve. The inferred physiological descriptions of the metrics are summarized in Table 2. The PhenoMetrics function of the CropPhenology package takes the time series of MODIS NDVI imagery and the boundary shape file of the study area and provides outputs of 15 phenological metrics in raster file format.
For the values of a given pixel to be comparable across time, we normalized the values relative to the mean of the field, using the standard score. The standard score value at the given pixel indicates how the pixel value deviates from the field mean in the measure of the standard deviation. Blackmore et al.  have used similar standardization to compare the relative yields among different crops. Accordingly, we calculated the pixel’s standard score for each phenological metrics, using Equation (1).
Figure 4. Illustration of NDVI dynamics and phenological metrics.
Table 2. Description of phenological metrics and their relation to yield.
where Standard scoreijk is a standard score value of metric k, at pixel i for year j, Valueijk is the value of metric k at pixel i for a year j, meanjk and sdjk is the mean and standard deviation of the metric k across all pixels of year j.
2.3.3. Spatial Trend and Temporal Variability of Phenological Metrics
Spatial trend map
The spatial trend at a pixel was calculated as the average of the pixel value over the year of interest (Equation (2)). It is adapted from Blackmoreet. Al.  . The spatial trend values of the phenological metrics from each pixel produced data for spatial trend map of that phenological metric, which shows the average value of the metric across observation years. A total of 15 spatial trend maps were created for the 15 metrics.
where Spatial trendik is the mean value of phenological metric k for the pixel i, standard Scoreijk is the standard score value for phenological metric k at of pixel i on year j, and n is the total number of observation year.
Temporal variability map
The temporal variability map was created using the temporal variance overthe observation years. It is adapted from Blackmoreet. Al.  . Temporal variability is classified into two categories: inter-year variability and relative temporal variability  . Inter-year variability is the temporal variability caused by the change in annual rainfall, which affects the overall productivity of the field. The second category described how each part of the field behaved relative to the other parts from year to year. In this analysis, we are interested in the second category, which assesses the relative variability of growth at a particular pixel, over the study period. This relative temporal variability was calculated using a temporal variance, a modified variance measured in a function of time   . The temporal variance of all the metrics were calculated for each pixel across years, using Equation (3). And the temporal variance value of the metric from each pixels produced the temporal variability map of the phenological metric. Accordingly, 15 temporal variability maps were created for the 15 metrics. The temporal variance of each metric across the season revealed the temporal variability of the growth conditions at the various phenological stages. The calculated temporal variance has low values if the pixel is stable over time and often has values close to the mean of the observation seasons. On the other hand, if the metric is temporally unstable, with high values in some years and low values in others, its temporal variance will be high.
where Temporal Varianceik is the temporal variance at pixel i of metric k, standard scoreijk is the standard score of metric k at pixel i for year j, Spatial trendi is the Spatial trend of pixel i and n is the number of observation year.
2.3.4. Relationship between Management Zone and Trends of Phenological Metrics
In this study, the monotonic relationship between the management zone values and the trends of the 15 phenological metrics were assessed. Spearman rank correlation was used to test the relative direction and strength of the relationship between the management zone pixel values and phenological metrics, in R computing environment. Spearman correlation ranks the variables and provides the correlation coefficient, rho (ρ) that indicates the strength of the correlation  .
3. Results and Discussion
3.1. Relationship between Management Zone and Spatial trend of Phenological Metrics
The Spearman rank correlations between management zone and the spatial trend of phenological metrics are shown in Table 3. In this analysis, we have undertaken multiple tests on a single dataset. This can introduce the probability of making type I errors, family wise error rate (FWER). To address this problem, we used the Bonferroni correction that adjusts the significance level (p value)  , by dividing the p value by the number of comparisons made, which results a new p value of 0.0333.
The results indicate that the spatial trend of Max Value is significantly
Table 3. Spearman correlation coefficient and statistical significance (p value) of the relationship between spatial trend of phenological metrics and the management zone.
correlated with the management zone values (ρ(13) = 0.82, p = 0.024). Although the associations are not significant, the GreenUpSlope and MaxT showed high correlation coefficients, with different directions of association. The Asymmetry, TINDVIBeforeMax and BeforeMaxT metrics also showed moderate association with the management zone values.
This result corroborates previous observations that have recognized the ability of NDVI to estimate crop yield    . Although strong associations between NDVI and yield have been reported at differing stages during the growth period, these dates are mostly around heading and flowering stages         . As crops often attain their maximum NDVI close to heading and flowering stages  , the high association of MaxValue with the management zone can be directly linked to these observations.
The results also show trends of increasing Max Time and decreasing of GreenUpSlope, moving from the poor to good management zone. As the soil PAWC strongly relates with yield potential in South Australia, the observed relationship is likely be related with water availability. Numbers of researchers reported that water stress has a strong influence on crop growth stages with evident observation of faster growth to the next growth stages in water stressed crops. Under similar climatic conditions, plants in low PAWC soils will likely to face water stress between rainfall events, as there is less water stored in the soil profile compared with the high PAWC soils. The early Maximum (low MaxT) can then be related to the hasten growth of the crop to anthesis due to water stress in the poor zone     . This is also reflected as faster rate of greenness, higher GreenUpSlope, in low PAWC soil.
Similar observations of low Max Time and high GreenUpSlope were observed in our previous study  . Unlike our current approach of averaging values across the years, the previous study considered paired comparison of pixels across number of years. The pared comparison study showed higher Max Values in low PAWC soils than the high PAWC soils. Such pattern of paired comparison observed in the previous study may average out thorough the observation years, in the current study, to reflect average higher Max Value for higher PAWC soils.
The positive moderate association of Asymmetry, TINDVIBeforeMax and BeforeMaxT metrics, indicated that the good zones have higher values of these metrics than the poor zone. Generally, this observation shows that the good zone experienced longer time before the NDVI peak, higher area under the curve between the Onset and Maximum NDVI, and higher difference between TINDVIBeforeMax and TINDVIAfetrMax than poor zone. These differences are related with the high and late NDVI peak in the good zone than the poor zone, which is explained above. Furthermore, the trend of TINDVIBeforeMax in our analysis agrees with pevious field studies that confirm the strong association of pre-anthesis growth with yield potential, as it is the time when the number of grains are determined  .
3.2. Relationship between Management Zone and Temporal Variability of Phenological Metrics
The seasons considered in this study, 2001-2015, were characterized by variable rainfall amount and seasonality. The temporal variance of the phenological metrics measures the relative variability of the metrics across the variable rainfall. The Spearman correlation analysis assessed the trend of temporal variability of phenological metrics across management zones (Table 4).
Table 4. Spearman correlation coefficients and statistical significance (p value) of the relationship between temporal variability of phenological metrics and the management zone.
The results indicate that the temporal variability of TINDVI has a significant positive correlation with the management zone, indicating an increase of temporal variability from poor to good management zones. This variability demonstrates the interaction of the spatial soil variability with the temporal rainfall seasonality. In good seasons, the spatial variability of soil PAWC becomes significant, as the high PAWC soils store more water than the low PAWC soils. However, in dry seasons there will be less difference in stored water between the contrasting soils. This results in larger spatial variability in good seasons compared with dry seasons   . Temporally, the good zone with high PAWC soil will have higher variability across variable seasons than the poor zone with low PAWC soils  .
Previous research has reported that TINDVI can be a good indicator for crop production    . However, the sensitivity of TINDVI for crop production depends on rainfall seasonality; it can be a good indicator for crop production when water is the limiting factor for crop growth  . Considering the profound correlation between TINDVI and yield in cropping fields, the high temporal variability in TINDVI in good zone than in poor zone coincide well with the yield data observation in the study area. Yield data from EPARF research shows higher difference in good zone (3.1 t/ha) between the dry year 2008 and the wet year 2010 than the poor zone (2.13 t/ha)   .
The relationship between the temporal variability of MaxV and the management zone showed weak correlation. We have observed in the previous section that the spatial trend (temporal mean) of this metrics showed a strong correlation with management zone values, when the low correlation of temporal variability with the management zone can be interpreted as both the good and poor zones are variable to a similar degree with the higher mean values from good zone than the poor zone. Thus, the MaxV metrics is more sensitive to spatial yield variability than temporal variability; on the other hand, TINDVI better reflects the temporal yield variability.
3.3. Summary of Indicative Phenological Metrics for Crop Field Management
Generally, our analysis demonstrates the potential of phenological metrics to assess the intra-field crop growth variability and to provide a comprehensive understanding of soil-climate interactions. Single vegetation index images have been successfully utilized to recognize yield variability across a crop field [e.g.  ]. However, the single image approach has been criticized for lacking information on intra-seasonal growth dynamics  . On the other hand, multi- temporal images and the derived phenological metrics uncover the intra-annual biophysical properties of the crop across the field, as driven by soil-climate interaction. This potentially provides a better understanding of both inter-annual and intra-annual variability, which is a key factor to improve management practices.
Furthermore, the method presented in this analysis show the efficacy of phenological metrics to recognize crop growth variability. Considering the increasing availability of remote sensing imagery, the spatio-temporal variability estimation using phenological metrics can provide valuable information for PA suitability assessment in areas where there is limited availability of yield data.
Currently, MODIS is the most appropriate imagery for such analysis due to its high temporal resolution. However, its coarse spatial resolution limits the application of this method to broad acre cropping. The technology of integrating the higher spatial resolution of sensors such as Landsat or Sentinel with the temporal resolution of MODIS to generate high temporal/high spatial resolution datasets  may allow future application for smaller fields. Whilst these fused spatio-temporal datasets are presently not available to the public at broad extents, once available, may be utilized at a scale typical for yield mapping in precision agriculture today. In such a way, this paper provides a basis for further analysis to utilize phenological metrics for crop management support systems.
In this study, the spatial and temporal variability of crop growth was assessed using remote sensing phenological metrics in relation to the relative yield potential of sites in the field. Phenological metrics provide information about the temporal and spatial variability of plant growth across cropping fields. Our study demonstrates the potential of satellite based phenological metrics to provide information about growth variability in cropping fields, which provide relevant information for agricultural management. The results of our analysis show that Time-integrated NDVI reflects seasonal effect on crop growth, whereas the magnitude of NDVI peak strongly reflects the soil quality and showed spatial variability of long-term site conditions across the field.
Crop yield response to PAWC is strongly influenced by the amount and seasonal distribution of rainfall, which is temporally variable. The effect is stronger when there is sufficient rainfall to fill the soil profile at the start of the season and the plants are dependent on deep stored soil-water as the season progresses. Phenological metrics provide comprehensive insight of the spatio-temporal crop growth variability across the variable seasons, which advances our understanding of soil-climate interaction across the cropping field.
The method presented here provides a pathway towards better estimation of spatio-temporal crop growth variability, which is vital for PA success. It can be used to develop future, more detailed studies to fully utilize the potential of phenological indicators for site characterization through unravelling the complex spatio-temporal soil-climate interactions.
The Authors wish to thank the Eyre Peninsula Agricultural Research Foundation particularly Linden Masters for data support and the Faculty of Sciences, University of Adelaide for research funding.
 Bramley, R.G.V. (2009) Lessons from Nearly 20 Years of Precision Agriculture Research, Development, and Adoption as a Guide to Its Appropriate Application. Crop and Pasture Science, 60, 197-217.
 Chang, J., Clay, D.E., Carlson, C.G., Clay, S.A., Malo, D.D., Berg, R., Kleinjan, J. and Wiebold, W. (2003) Different Techniques to Identify Management Zones Impact Nitrogen and Phosphorus Sampling Variability. Agronomy Journal, 95, 1550-1559.
 Mzuku, M.K., Reich, R., Inman, R., Smith, D. and MacDonald, F. (2005) Spatial Variability of Measured Soil Properties across Site-Specific Management Zones. Soil Science Society of America Journal, 69, 1572.
 Xiang, L.I., Yu-chun, P.A.N., Zhong-qiang, G.E., Chun-jiang, Z. (2007) Delineation and Scale Effect of Precision Agriculture Management Zones Using Yield Monitor Data over Four Years. Agricultural Sciences in China, 6, 180-188.
 Corwin, D.L. (2013) Site-Specific Management and Delineating Management Zones. In: Oliver, M., Bishop, T. and Marchant, B., Eds., Precision Agriculture for Sustainability and Environmental Protection. Taylor and Francis, Florence,
 Jochinke, D.C., Noonon, B.J., Wachsmann, N.G. and Norton, R.M. (2007) The Adoption of Precision Agriculture in an Australian Broadacre Cropping System—Challenges and Opportunities. Field Crops Research, 104, 68-76.
 Oliver, Y., Wong, M., Robertson, R. and Wittwer, K. (2006) PAWC Determines Spatial Variability in Grain Yield and Nitrogen Requirement by Interacting with Rainfall on Northern WA Sandplain. Proceedings 13th Australian Society of Agronomy Conference, Perth, Western Australia, 10-14 September 2006, 5.
 Oliver, Y.M., Robertson, M.J., Stone, P.J. and Whitbread, A. (2009) Improving Estimates of Water-Limited Yield of Wheat by Accounting for Soil Type and Within-Season Rainfall. Crop and Pasture Science, 60, 1137-1146.
 Rab, M., Fisher, P., Armstrong, R., Abuzar, M., Robinson, N. and Chandra, S. (2009) Advances in Precision Agriculture in South-Eastern Australia. IV. Spatial Variability in Plant-Available Water Capacity of Soil and Its Relationship with Yield in Site-Specific Management Zones. Crop and Pasture Science, 60, 885-900.
 Florin, M.J., McBratney, A.B. and Whelan, B.M. (2009) Quantification and Comparison of Wheat Yield Variation across Space and Time. European Journal of Agronomy, 30, 212-219.
 McBratney, A.X.B., Whelan, B.M. and Shatar, T.M. (2007) Variability and Uncertainty in Spatial, Temporal and Spatiotemporal Crop-Yield and Related Data. In: Lake, J.V., Bock, G.R. and Goode, J.A., Eds., Ciba Foundation Symposium 210-Precision Agriculture: Spatial and Temporal Variability of Environmental Quality, John Wiley & Sons, New York, 141-160.
 Blackmore, S., Godwin, R.J. and Fountas, S. (2003) The Analysis of Spatial and Temporal Trends in Yield Map Data over Six Years. Biosystems Engineering, 84, 455-466.
 Kaspar, T.C., Colvin, T.S., Jaynes, D.B., Karlen, D.L., James, D.E., Meek, D.W., Pulido, D. and Butler, H. (2003) Relationship between Six Years of Corn Yields and Terrain Attributes. Precision Agriculture, 4, 87-101.
 Kaspar, T.C., Pulido, D.J., Fenton, T.E., Colvin, T.S., Karlen, D.L., Jaynes, D.B. and Meek, D.W. (2004) Relationship of Corn and Soybean Yield to Soil and Terrain Properties. Agronomy Journal, 96, 700-709.
 Abuzar, M., Rampant, P. and Fisher, P. (2004) Measuring Spatial Variability of Crops and Soils at Sub-Paddock Scale Using Remote Sensing Technologies. Proceedings of Geoscience and Remote Sensing Symposium, Anchorage, 20-24 September 2004, 1633-1636.
 Anwar, M.R., O’Leary, G.J., Rab, M.A., Fisher, P.D. and Armstrong, R.D. (2009) Advances in Precision Agriculture in South-Eastern Australia. V. Effect of Seasonal Conditions on Wheat and Barley Yield Response to Applied Nitrogen across Management Zones. Crop and Pasture Science, 60, 901-911.
 Lyle, G., Lewis, M. and Ostendorf, B. (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, 1549-1567.
 Mkhabela, M.S., Bullock, P., Raj, S., Wang, S. and Yang, Y. (2011) Crop Yield Forecasting on the Canadian Prairies Using MODIS NDVI Data. Agricultural and Forest Meteorology, 151, 385-393.
 Smith, R., Adams, J., Stephens, D. and Hick, P. (1995) Forecasting Wheat Yield in a Mediterranean-Type Environment from the NOAA Satellite. Australian Journal of Agricultural Research, 46, 113-125.
 Basso, B., Ritchie, J.T., Pierce, F.J., Braga, R.P. and Jones, J.W. (2001) Spatial Validation of Crop Models for Precision Agriculture. Agricultural Systems, 68, 97-112.
 Buttafuoco, G., Castrignanò, A., Cucci, G., Rinaldi, M. and Ruggieri, S. (2015) An Approach to Delineate Management Zones in a Durum Wheat Field: Validation Using Remote Sensing and Yield Mapping. Proceedings of 10th European Conference on Precision Agriculture, Beit-Dagan, 12-16 July 2015, 241-247.
 Maynard, J.J. and Levi, M.R. (2017) Hyper-Temporal Remote Sensing for Digital Soil Mapping: Characterizing Soil-Vegetation Response to Climatic Variability. Geoderma, 285, 94-109.
 Zhang, X., Friedl, M.A., Schaaf, C.B., Strahler, A.H., Hodges, J.C.F., Gao, F., Reed, B.C. and Huete, A. (2003) Monitoring Vegetation Phenology Using MODIS. Remote Sensing of Environment, 84, 471-475.
 Reed, B.C., Schwartz, M.D. and Xiao, X. (2009) Remote Sensing Phenology: Status and the Way Forward. In: Reed, B.C., Schwartz, M.D. and Xiao, X., Eds., Phenology of Ecosystem Processes, Springer, New York, 231-246.
 Araya, S., Lyle, G., Lewis, M. and Ostendorf, B. (2016) Phenologic Metrics Derived from MODIS NDVI as Indicators for Plant Available Water-Holding Capacity. Ecological Indicators, 60, 1263-1272.
 Latta, R., Wilhelm, N. and Treloar, P. (2013) Summary of Paddock North 1 VRT Study at MAC. In: Scholz, N., Cook, A., Latta, R., Wilhelm, N., Crettenden, J., Holbery, S. and Masters, L., Eds., Eyre Peninsula Farming Systems Summary, Grains Research & Development Corporation, Barton, 112-114.
 Huete, A., Didan, K., Miura, T., Rodriguez, E.P., Gao, X. and Ferreira, L.G. (2002) Overview of the Radiometric and Biophysical Performance of the MODIS Vegetation Indices. Remote Sensing of Environment, 83, 195-213.
 Vermote, E.F., El Saleous, N.Z. and Justice, C.O. (2002) Atmospheric Correction of MODIS Data in the Visible to Middle Infrared: First Results. Remote Sensing of Environment, 83, 97-111.
 Doraiswamy, P.C., Sinclair, T.R., Hollinger, S., Akhmedov, B., Stern, A. and Prueger, J. (2005) Application of MODIS Derived Parameters for Regional Crop Yield Assessment. Remote Sensing of Environment, 97, 192-202.
 Quarmby, N.A., Milnes, M., Hindle, T.L. and Silleos, N. (1993) The Use of Multi-Temporal NDVI Measurements from AVHRR Data for Crop Yield Estimation and Prediction. International Journal of Remote Sensing, 14, 199-210.
 Peralta, N., Assefa, Y., Du, J., Barden, C. and Ciampitti, I. (2016) Mid-Season High-Resolution Satellite Imagery for Forecasting Site-Specific Corn Yield. Remote Sensing, 8, 848.
 Schut, A.G.T., Stephens, D.J., Stovold, R.G.H., Adams, M. and Craig, R.L. (2009) Improved Wheat Yield and Production Forecasting with a Moisture Stress Index, AVHRR and MODIS Data. Crop and Pasture Science, 60, 60-70.
 Shanahan, J.F., Schepers, J.S., Francis, D.D., Varvel, G.E., Wilhelm, W.W., Tringe, J.M., Schlemmer, M.R. and Major, D.J. (2001) Use of Remote-Sensing Imagery to Estimate Corn Grain Yield. Agronomy Journal, 93, 583-589.
 Tagarakis, A., Liakos, V., Fountas, S., Koundouras, S. and Gemtos, T.A. (2012) Management Zones Delineation Using Fuzzy Clustering Techniques in Grapevines. Precision Agriculture, 14, 18-39.
 Sakamoto, T., Yokozawa, M., Toritani, H., Shibayama, M., Ishitsuka, N. and Ohno, H. (2005) A Crop Phenology Detection Method Using Time-Series MODIS Data. Remote Sensing of Environment, 96, 366-374.
 Araya, S., Ostendorf, B., Lyle, G. and Lewis, M. (2013) Crop Phenology Based on MODIS Satellite Imagery as an Indicator of Plant Available Water Content. Proceedings of 20th International Congress on Modelling and Simulation, Adelaide, 1-6 December 2013, 1896-1902.
 McMaster, G.S., Wilhelm, W.W. and Frank, A.B. (2005) Developmental Sequences for Simulating Crop Phenology for Water-Limiting Conditions. Australian Journal of Agricultural Research, 56, 1277-1288.
 Poole, N. and Hunt, J. (2014) Advancing the Management of Crop Canopies—Keeping Crops Greener for Longer. In: Perry, M., Ed., A Cereal Crop Management Guide, Grains Research & Development Corporation, Barton, 1-11.
 Wong, M. and Asseng, S. (2006) Determining the Causes of Spatial and Temporal Variability of Wheat Yields at Sub-Field Scale Using a New Method of Upscaling a Crop Model. Plant and Soil, 283, 203-215.
 Hill, M.J. and Donald, G.E. (2003) Estimating Spatio-Temporal Patterns of Agricultural Productivity in Fragmented Landscapes Using AVHRR NDVI Time Series. Remote Sensing of Environment, 84, 367-384.
 Paterson, C., Latta, R., Treloar, P., Shepperd, W. and Richter, I. (2011) Can Adjusting Zones within N1 Paddock at Minnipa Improve VRT Outcomes? In: Scholz, N., Cook, A., Latta, R., Wilhelm, N., Paterson, C., Crettenden, J. and McNeill, A., Eds., Eyre Peninsula Farming Systems Summary, Grains Research & Development Corporation, Barton,102-106.
 Gao, F., Anderson, M.C., Zhang, X., Yang, Z., Alfieri, J.G., Kustas, W.P., Mueller, R., Johnson, D.M. and Prueger, J.H. (2017) Toward Mapping Crop Progress at Field Scales through Fusion of Landsat and MODIS Imagery. Remote Sensing of Environment, 188, 9-25.