{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,3,20]],"date-time":"2026-03-20T15:36:59Z","timestamp":1774021019109,"version":"3.50.1"},"reference-count":52,"publisher":"MDPI AG","issue":"10","license":[{"start":{"date-parts":[[2014,10,13]],"date-time":"2014-10-13T00:00:00Z","timestamp":1413158400000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":["Remote Sensing"],"abstract":"<jats:p>Policy makers, government planners and agricultural market participants in Pakistan require accurate and timely information about wheat yield and production. Punjab Province is by far the most important wheat producing region in the country. The manual collection of field data and data processing for crop forecasting by the provincial government requires significant amounts of time before official reports can be released. Several studies have shown that wheat yield can be effectively forecast using satellite remote sensing data. In this study, we developed a methodology for estimating wheat yield and area for Punjab Province from freely available Landsat and MODIS satellite imagery approximately six weeks before harvest. Wheat yield was derived by regressing reported yield values against time series of four different peak-season MODIS-derived vegetation indices. We also tested deriving wheat area from the same MODIS time series using a regression-tree approach. Among the four evaluated indices, WDRVI provided more consistent and accurate yield forecasts compared to NDVI, EVI2 and saturation-adjusted normalized difference vegetation index (SANDVI). The lowest RMSE values at the district level for forecast versus reported yield were found when using six or more years of training data. Forecast yield for the 2007\/2008 to 2012\/2013 growing seasons were within 0.2% and 11.5% of final reported values. Absolute deviations of wheat area and production forecasts from reported values were slightly greater compared to using the previous year's or the three- or six-year moving average values, implying that 250-m MODIS data does not provide sufficient spatial resolution for providing improved wheat area and production forecasts.<\/jats:p>","DOI":"10.3390\/rs6109653","type":"journal-article","created":{"date-parts":[[2014,10,14]],"date-time":"2014-10-14T02:13:08Z","timestamp":1413252788000},"page":"9653-9675","update-policy":"https:\/\/doi.org\/10.3390\/mdpi_crossmark_policy","source":"Crossref","is-referenced-by-count":107,"title":["Wheat Yield Forecasting for Punjab Province from Vegetation Index Time Series and Historic Crop Statistics"],"prefix":"10.3390","volume":"6","author":[{"given":"Jan","family":"Dempewolf","sequence":"first","affiliation":[{"name":"Department of Geographical Sciences, University of Maryland, 2181 Samuel J. LeFrak Hall, College Park, MD 20742, USA"}]},{"given":"Bernard","family":"Adusei","sequence":"additional","affiliation":[{"name":"Radius Technology Group Inc., 804 Pershing Dr, Silver Spring, MD 20910, USA"}]},{"given":"Inbal","family":"Becker-Reshef","sequence":"additional","affiliation":[{"name":"Department of Geographical Sciences, University of Maryland, 2181 Samuel J. LeFrak Hall, College Park, MD 20742, USA"}]},{"ORCID":"https:\/\/orcid.org\/0000-0003-0042-2767","authenticated-orcid":false,"given":"Matthew","family":"Hansen","sequence":"additional","affiliation":[{"name":"Department of Geographical Sciences, University of Maryland, 2181 Samuel J. LeFrak Hall, College Park, MD 20742, USA"}]},{"given":"Peter","family":"Potapov","sequence":"additional","affiliation":[{"name":"Department of Geographical Sciences, University of Maryland, 2181 Samuel J. LeFrak Hall, College Park, MD 20742, USA"}]},{"given":"Ahmad","family":"Khan","sequence":"additional","affiliation":[{"name":"Department of Geographical Sciences, University of Maryland, 2181 Samuel J. LeFrak Hall, College Park, MD 20742, USA"}]},{"given":"Brian","family":"Barker","sequence":"additional","affiliation":[{"name":"Department of Geographical Sciences, University of Maryland, 2181 Samuel J. LeFrak Hall, College Park, MD 20742, USA"}]}],"member":"1968","published-online":{"date-parts":[[2014,10,13]]},"reference":[{"key":"ref_1","unstructured":"Health Nutrition and Population Statistics. Available online: http:\/\/databank.worldbank.org."},{"key":"ref_2","unstructured":"Briscoe, J., and Qamar, U. (2005). Pakistan\u2019s Water Economy: Running Dry, Oxford University Press."},{"key":"ref_3","unstructured":"Khan, S.B., Khattak, A., and Rehman, A. (2012). Agricultural Statistics of Pakistan 2010\u20132011, Pakistan Bureau of Statistics, Government of Pakistan."},{"key":"ref_4","unstructured":"Pakistan Economic Survey 2011\u20132012, Available online: http:\/\/www.finance.gov.pk\/survey_1112.html."},{"key":"ref_5","doi-asserted-by":"crossref","first-page":"71","DOI":"10.30541\/v47i1pp.71-87","article-title":"Wheat markets and price stabilisation in Pakistan: An analysis of policy options","volume":"47","author":"Dorosh","year":"2008","journal-title":"Pak. Dev. Rev"},{"key":"ref_6","unstructured":"Niaz, M.S. Wheat Policy\u2014A Success or Failure. DAWN. Available online: http:\/\/www.dawn.com\/news\/194423\/wheat-policy-a-success-or-failure\/."},{"key":"ref_7","unstructured":"Akhtar, I.U.H. Pakistan Needs a New Crop Forecasting System. Available online: http:\/\/www.scidev.net\/en\/new-technologies\/space-technology\/opinions\/pakistan-needs-a-new-crop-forecasting-system.html."},{"key":"ref_8","doi-asserted-by":"crossref","first-page":"271","DOI":"10.1016\/S0065-2113(08)70012-9","article-title":"The role of remote sensing in determining the distribution and yield of crops","volume":"27","author":"Bauer","year":"1975","journal-title":"Adv. Agron"},{"key":"ref_9","unstructured":"Erickson, J.D. (1984). The Role of Terrestrial Vegetation in the Global Carbon Cycle: Measurement by Remote Sensing, John Wiley."},{"key":"ref_10","unstructured":"Doraiswamy, P., Hodges, T., and Phinney, D. (1979). Crop Yield Literature Review for AgRISTARS Crops Corn, Soybeans, Wheat, Barley, Sorghum, Rice, Cotton, and Sunflowers, Lockheed Electronics Co.. Technical Report 73\u2013312."},{"key":"ref_11","first-page":"935","article-title":"The use of remote sensing within the MARS crop yield monitoring system of the European Commission","volume":"37","author":"Baruth","year":"2008","journal-title":"Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci"},{"key":"ref_12","doi-asserted-by":"crossref","first-page":"1589","DOI":"10.3390\/rs2061589","article-title":"Monitoring global croplands with coarse resolution earth observations: The Global Agriculture Monitoring (GLAM) project","volume":"2","author":"Justice","year":"2010","journal-title":"Remote Sens"},{"key":"ref_13","doi-asserted-by":"crossref","first-page":"321","DOI":"10.1016\/S0167-8809(02)00034-8","article-title":"A new crop yield forecasting model based on satellite measurements applied across the Indus Basin, Pakistan","volume":"94","author":"Bastiaanssen","year":"2003","journal-title":"Agric. Ecosyst. Environ"},{"key":"ref_14","unstructured":"Crop Situation Forecast. Available online: http:\/\/www.suparco.gov.pk\/pages\/pak-scms.asp."},{"key":"ref_15","doi-asserted-by":"crossref","first-page":"1844","DOI":"10.3390\/rs2071844","article-title":"Estimating global cropland extent with multi-year MODIS data","volume":"2","author":"Pittman","year":"2010","journal-title":"Remote Sens"},{"key":"ref_16","doi-asserted-by":"crossref","first-page":"123","DOI":"10.1007\/BF00058655","article-title":"Bagging predictors","volume":"24","author":"Breiman","year":"1996","journal-title":"Mach. Learn"},{"key":"ref_17","doi-asserted-by":"crossref","first-page":"209","DOI":"10.1016\/0034-4257(95)00142-5","article-title":"Global discrimination of land cover types from metrics derived from AVHRR pathfinder data","volume":"54","author":"DeFries","year":"1995","journal-title":"Remote Sens. Environ"},{"key":"ref_18","doi-asserted-by":"crossref","first-page":"303","DOI":"10.1016\/S0034-4257(02)00079-2","article-title":"Towards an operational MODIS continuous field of percent tree cover algorithm: Examples using AVHRR and MODIS data","volume":"83","author":"Hansen","year":"2002","journal-title":"Remote Sens. Environ"},{"key":"ref_19","doi-asserted-by":"crossref","first-page":"106","DOI":"10.1016\/j.rse.2011.08.027","article-title":"Quantifying forest cover loss in Democratic Republic of the Congo, 2000\u20132010, with Landsat ETM+ data","volume":"122","author":"Potapov","year":"2012","journal-title":"Remote Sens. Environ"},{"key":"ref_20","doi-asserted-by":"crossref","first-page":"412","DOI":"10.1016\/j.rse.2004.08.002","article-title":"Cropland distributions from temporal unmixing of MODIS data","volume":"93","author":"Lobell","year":"2004","journal-title":"Remote Sens. Environ"},{"key":"ref_21","doi-asserted-by":"crossref","first-page":"1190","DOI":"10.1016\/j.rse.2010.01.006","article-title":"The spatial distribution of crop types from MODIS data: Temporal unmixing using Independent Component Analysis","volume":"114","author":"Ozdogan","year":"2010","journal-title":"Remote Sens. Environ"},{"key":"ref_22","doi-asserted-by":"crossref","first-page":"1096","DOI":"10.1016\/j.rse.2007.07.019","article-title":"Large-area crop mapping using time-series MODIS 250 m NDVI data: An assessment for the U.S. Central Great Plains","volume":"112","author":"Wardlow","year":"2008","journal-title":"Remote Sens. Environ"},{"key":"ref_23","doi-asserted-by":"crossref","first-page":"1654","DOI":"10.2134\/agronj2007.0170","article-title":"Corn and soybean mapping in the United States using MODIS time-series data sets","volume":"99","author":"Chang","year":"2007","journal-title":"Agron. J"},{"key":"ref_24","doi-asserted-by":"crossref","first-page":"1312","DOI":"10.1016\/j.rse.2010.01.010","article-title":"A generalized regression-based model for forecasting winter wheat yields in Kansas and Ukraine using MODIS data","volume":"114","author":"Vermote","year":"2010","journal-title":"Remote Sens. Environ"},{"key":"ref_25","first-page":"19","article-title":"Physical geography of the punjab","volume":"11","author":"Gosal","year":"2004","journal-title":"J. Punjab Stud"},{"key":"ref_26","doi-asserted-by":"crossref","first-page":"1965","DOI":"10.1002\/joc.1276","article-title":"Very high resolution interpolated climate surfaces for global land areas","volume":"25","author":"Hijmans","year":"2005","journal-title":"Int. J. Climatol"},{"key":"ref_27","unstructured":"Steduto, P., Hsiao, T., Fereres, E., and Raes, D. (2012). Crop Yield Response to Water, FAO. Technical Report 66."},{"key":"ref_28","unstructured":"(2012). Punjab CRS: Base Line Survey\u2014Agriculture Information System, FAO\/UN SUPARCO. Technical Report."},{"key":"ref_29","doi-asserted-by":"crossref","first-page":"171","DOI":"10.1016\/0034-4257(81)90018-3","article-title":"Remote sensing of total dry-matter accumulation in winter wheat","volume":"11","author":"Tucker","year":"1981","journal-title":"Remote Sens. Environ"},{"key":"ref_30","doi-asserted-by":"crossref","first-page":"4169","DOI":"10.1080\/01431160110107653","article-title":"Wheat yield estimates using multi-temporal NDVI satellite imagery","volume":"23","author":"Labus","year":"2002","journal-title":"Int. J. Remote Sens"},{"key":"ref_31","doi-asserted-by":"crossref","first-page":"385","DOI":"10.1016\/j.agrformet.2010.11.012","article-title":"Crop yield forecasting on the Canadian Prairies using MODIS NDVI data","volume":"151","author":"Mkhabela","year":"2011","journal-title":"Agric. For. Meteorol"},{"key":"ref_32","doi-asserted-by":"crossref","unstructured":"Dempewolf, J., Adusei, B., Becker-Reshef, I., Barker, B., Potapov, P., Hansen, M., and Justice, C. (2013, January 21\u201326). Wheat production forecasting for Pakistan from satellite data. Melbourne, VIC, Australia.","DOI":"10.1109\/IGARSS.2013.6723517"},{"key":"ref_33","first-page":"403","article-title":"Regional yield estimation for winter wheat with MODIS-NDVI data in Shandong, China","volume":"10","author":"Ren","year":"2008","journal-title":"Int. J. Appl. Earth Obs. Geoinf"},{"key":"ref_34","doi-asserted-by":"crossref","first-page":"539","DOI":"10.3390\/rs5020539","article-title":"Remote sensing based yield estimation in a stochastic framework\u2014Case study of durum wheat in Tunisia","volume":"5","author":"Meroni","year":"2013","journal-title":"Remote Sens"},{"key":"ref_35","doi-asserted-by":"crossref","first-page":"115","DOI":"10.1016\/j.rse.2008.08.015","article-title":"Phenologically-tuned MODIS NDVI-based production anomaly estimates for Zimbabwe","volume":"113","author":"Funk","year":"2009","journal-title":"Remote Sens. Environ"},{"key":"ref_36","doi-asserted-by":"crossref","first-page":"127","DOI":"10.1016\/0034-4257(79)90013-0","article-title":"Red and photographic infrared linear combinations for monitoring vegetation","volume":"8","author":"Tucker","year":"1979","journal-title":"Remote Sens. Environ"},{"key":"ref_37","doi-asserted-by":"crossref","first-page":"1","DOI":"10.1016\/j.ecolind.2013.01.041","article-title":"NDVI saturation adjustment: A new approach for improving cropland performance estimates in the Greater Platte River Basin, USA","volume":"30","author":"Gu","year":"2013","journal-title":"Ecol. Indic"},{"key":"ref_38","doi-asserted-by":"crossref","first-page":"440","DOI":"10.1016\/S0034-4257(96)00112-5","article-title":"A comparison of vegetation indices over a global set of TM images for EOS-MODIS","volume":"59","author":"Huete","year":"1997","journal-title":"Remote Sens. Environ"},{"key":"ref_39","doi-asserted-by":"crossref","first-page":"609","DOI":"10.1016\/S0034-4257(00)00150-4","article-title":"Optical-biophysical relationships of vegetation spectra without background contamination","volume":"74","author":"Gao","year":"2000","journal-title":"Remote Sens. Environ"},{"key":"ref_40","doi-asserted-by":"crossref","first-page":"3833","DOI":"10.1016\/j.rse.2008.06.006","article-title":"Development of a two-band enhanced vegetation index without a blue band","volume":"112","author":"Jiang","year":"2008","journal-title":"Remote Sens. Environ"},{"key":"ref_41","doi-asserted-by":"crossref","first-page":"165","DOI":"10.1078\/0176-1617-01176","article-title":"Wide dynamic range vegetation index for remote quantification of biophysical characteristics of vegetation","volume":"161","author":"Gitelson","year":"2004","journal-title":"J. Plant Physiol"},{"key":"ref_42","doi-asserted-by":"crossref","unstructured":"Gitelson, A.A., Wardlow, B.D., Keydan, G.P., and Leavitt, B. (2007). An evaluation of MODIS 250-m data for green LAI estimation in crops. Geophys. Res. Lett.","DOI":"10.1029\/2007GL031620"},{"key":"ref_43","doi-asserted-by":"crossref","first-page":"15","DOI":"10.1016\/j.agrformet.2012.03.012","article-title":"An evaluation of MODIS 8- and 16-day composite products for monitoring maize green leaf area index","volume":"161","author":"Gitelson","year":"2012","journal-title":"Agric. Forest Meteorol"},{"key":"ref_44","doi-asserted-by":"crossref","first-page":"215","DOI":"10.1016\/j.rse.2012.12.017","article-title":"MODIS-based corn grain yield estimation model incorporating crop phenology information","volume":"131","author":"Sakamoto","year":"2013","journal-title":"Remote Sens. Environ"},{"key":"ref_45","doi-asserted-by":"crossref","first-page":"24","DOI":"10.2134\/agronj2013.0314","article-title":"Testing remote sensing approaches for assessing yield variability among maize fields","volume":"106","author":"Sibley","year":"2014","journal-title":"Agron. J"},{"key":"ref_46","doi-asserted-by":"crossref","first-page":"4155","DOI":"10.1080\/014311602320567955","article-title":"Improving an operational wheat yield model using phenological phase-based Normalized Difference Vegetation Index","volume":"23","author":"Boken","year":"2002","journal-title":"Int. J. Remote Sens"},{"key":"ref_47","doi-asserted-by":"crossref","first-page":"97","DOI":"10.4141\/P03-070","article-title":"Optimal time for remote sensing to relate to crop grain yield on the Canadian prairies","volume":"84","author":"Basnyat","year":"2004","journal-title":"Can. J. Plant Sci"},{"key":"ref_48","doi-asserted-by":"crossref","first-page":"389","DOI":"10.1111\/j.1751-5823.2005.tb00155.x","article-title":"Using remote sensing for agricultural statistics","volume":"73","author":"Carfagna","year":"2005","journal-title":"Int. Stat. Rev"},{"key":"ref_49","first-page":"1201","article-title":"The 2009 cropland data layer","volume":"76","author":"Johnson","year":"2010","journal-title":"Photogramm. Eng. Remote Sens"},{"key":"ref_50","doi-asserted-by":"crossref","first-page":"89","DOI":"10.5194\/adgeo-5-89-2005","article-title":"Comparison of different efficiency criteria for hydrological model assessment","volume":"5","author":"Krause","year":"2005","journal-title":"Adv. Geosci"},{"key":"ref_51","doi-asserted-by":"crossref","first-page":"2686","DOI":"10.1016\/j.rse.2011.05.026","article-title":"Crop specific green area index retrieval from MODIS data at regional scale by controlling pixel-target adequacy","volume":"115","author":"Duveiller","year":"2011","journal-title":"Remote Sens. Environ"},{"key":"ref_52","unstructured":"(2012). PAK-SCMS Bulletin\u2014Pakistan Satellite Based Crop Monitoring System, Government of Pakistan\u2014SUPARCO, UN-FAO. Monthly Bulletin Volume II, Issue 3."}],"container-title":["Remote Sensing"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/www.mdpi.com\/2072-4292\/6\/10\/9653\/pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2025,10,11]],"date-time":"2025-10-11T21:16:48Z","timestamp":1760217408000},"score":1,"resource":{"primary":{"URL":"https:\/\/www.mdpi.com\/2072-4292\/6\/10\/9653"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2014,10,13]]},"references-count":52,"journal-issue":{"issue":"10","published-online":{"date-parts":[[2014,10]]}},"alternative-id":["rs6109653"],"URL":"https:\/\/doi.org\/10.3390\/rs6109653","relation":{},"ISSN":["2072-4292"],"issn-type":[{"value":"2072-4292","type":"electronic"}],"subject":[],"published":{"date-parts":[[2014,10,13]]}}}