{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,1,31]],"date-time":"2026-01-31T02:33:24Z","timestamp":1769826804294,"version":"3.49.0"},"reference-count":51,"publisher":"MDPI AG","issue":"1","license":[{"start":{"date-parts":[[2017,1,15]],"date-time":"2017-01-15T00:00:00Z","timestamp":1484438400000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"funder":[{"name":"Spanish Ministry of Economy and Competitiveness","award":["MTM2014-51992-R"],"award-info":[{"award-number":["MTM2014-51992-R"]}]},{"name":"Government of Navarre, Spain","award":["PI015, 2016"],"award-info":[{"award-number":["PI015, 2016"]}]},{"name":"Fundaci\u00f3n Caja Navarra-UNED Pamplona, Spain","award":["2016"],"award-info":[{"award-number":["2016"]}]}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":["Remote Sensing"],"abstract":"<jats:p>The normalized difference vegetation index (NDVI) is an important indicator for evaluating vegetation change, monitoring land surface fluxes or predicting crop models. Due to the great availability of images provided by different satellites in recent years, much attention has been devoted to testing trend changes with a time series of NDVI individual pixels. However, the spatial dependence inherent in these data is usually lost unless global scales are analyzed. In this paper, we propose incorporating both the spatial and the temporal dependence among pixels using a stochastic spatio-temporal model for estimating the NDVI distribution thoroughly. The stochastic model is a state-space model that uses meteorological data of the Climatic Research Unit (CRU TS3.10) as auxiliary information. The model will be estimated with the Expectation-Maximization (EM) algorithm. The result is a set of smoothed images providing an overall analysis of the NDVI distribution across space and time, where fluctuations generated by atmospheric disturbances, fire events, land-use\/cover changes or engineering problems from image capture are treated as random fluctuations. The illustration is carried out with the third generation of NDVI images, termed NDVI3g, of the Global Inventory Modeling and Mapping Studies (GIMMS) in continental Spain. This data are taken in bymonthly periods from January 2011 to December 2013, but the model can be applied to many other variables, countries or regions with different resolutions.<\/jats:p>","DOI":"10.3390\/rs9010076","type":"journal-article","created":{"date-parts":[[2017,1,16]],"date-time":"2017-01-16T09:44:02Z","timestamp":1484559842000},"page":"76","update-policy":"https:\/\/doi.org\/10.3390\/mdpi_crossmark_policy","source":"Crossref","is-referenced-by-count":21,"title":["Stochastic Spatio-Temporal Models for Analysing NDVI Distribution of GIMMS NDVI3g Images"],"prefix":"10.3390","volume":"9","author":[{"ORCID":"https:\/\/orcid.org\/0000-0002-0631-3919","authenticated-orcid":false,"given":"Ana","family":"Militino","sequence":"first","affiliation":[{"name":"Department of Statistics and Operations Research, Public University of Navarre, 31006 Pamplona, Spain"},{"name":"Institute for Advanced Materials (InaMat), Public University of Navarre, 31006 Pamplona, Spain"},{"name":"Department of Mathematics, Spanish Open University (UNED), 31006 Pamplona, Spain"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-3505-8400","authenticated-orcid":false,"given":"Maria","family":"Ugarte","sequence":"additional","affiliation":[{"name":"Department of Statistics and Operations Research, Public University of Navarre, 31006 Pamplona, Spain"},{"name":"Institute for Advanced Materials (InaMat), Public University of Navarre, 31006 Pamplona, Spain"},{"name":"Department of Mathematics, Spanish Open University (UNED), 31006 Pamplona, Spain"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-2796-9079","authenticated-orcid":false,"given":"Unai","family":"P\u00e9rez-Goya","sequence":"additional","affiliation":[{"name":"Department of Statistics and Operations Research, Public University of Navarre, 31006 Pamplona, Spain"}],"role":[{"role":"author","vocabulary":"crossref"}]}],"member":"1968","published-online":{"date-parts":[[2017,1,15]]},"reference":[{"key":"ref_1","doi-asserted-by":"crossref","first-page":"1","DOI":"10.1046\/j.1365-2486.2003.00507.x","article-title":"Northern hemisphere photosynthetic trends 1982\u20131999","volume":"9","author":"Slayback","year":"2003","journal-title":"Glob. Chang. Biol."},{"key":"ref_2","doi-asserted-by":"crossref","first-page":"4485","DOI":"10.1080\/01431160500168686","article-title":"An extended AVHRR 8-km NDVI dataset compatible with MODIS and SPOT vegetation NDVI data","volume":"26","author":"Tucker","year":"2005","journal-title":"Int. J. Remote Sens."},{"key":"ref_3","first-page":"309","article-title":"Monitoring vegetation systems in the Great Plains with ERTS","volume":"351","author":"Rouse","year":"1974","journal-title":"NASA Spec. Publ."},{"key":"ref_4","doi-asserted-by":"crossref","first-page":"4267","DOI":"10.1080\/01431161.2010.486414","article-title":"Global trends in NDVI-derived parameters obtained from GIMMS data","volume":"32","author":"Sobrino","year":"2011","journal-title":"Int. J. Remote Sens."},{"key":"ref_5","doi-asserted-by":"crossref","first-page":"2113","DOI":"10.3390\/rs5052113","article-title":"Trend change detection in NDVI time series: Effects of inter-annual variability and methodology","volume":"5","author":"Forkel","year":"2013","journal-title":"Remote Sens."},{"key":"ref_6","doi-asserted-by":"crossref","first-page":"5238","DOI":"10.3390\/rs6065238","article-title":"Comparison of Medium Spatial Resolution ENVISAT-MERIS and Terra-MODIS time series for vegetation decline analysis: A case study in Central Asia","volume":"6","author":"Dubovyk","year":"2014","journal-title":"Remote Sens."},{"key":"ref_7","doi-asserted-by":"crossref","first-page":"267","DOI":"10.3390\/rs8040267","article-title":"Operational Drought Monitoring in Kenya Using MODIS NDVI Time Series","volume":"8","author":"Klisch","year":"2016","journal-title":"Remote Sens."},{"key":"ref_8","doi-asserted-by":"crossref","first-page":"269","DOI":"10.3390\/rs8040269","article-title":"Corn Response to Climate Stress Detected with Satellite-Based NDVI Time Series","volume":"8","author":"Wang","year":"2016","journal-title":"Remote Sens."},{"key":"ref_9","doi-asserted-by":"crossref","first-page":"13233","DOI":"10.3390\/rs71013233","article-title":"Spatial and temporal patterns of global NDVI trends: Correlations with climate and human factors","volume":"7","author":"Liu","year":"2015","journal-title":"Remote Sens."},{"key":"ref_10","doi-asserted-by":"crossref","first-page":"279","DOI":"10.1016\/j.rse.2014.06.021","article-title":"Operational monitoring of daily evapotranspiration by the combination of MODIS NDVI and ground meteorological data: Application and evaluation in Central Italy","volume":"152","author":"Maselli","year":"2014","journal-title":"Remote Sens. Environ."},{"key":"ref_11","doi-asserted-by":"crossref","first-page":"232","DOI":"10.1016\/j.rse.2013.07.020","article-title":"Monitoring and modeling spatial and temporal patterns of grassland dynamics using time series MODIS NDVI with climate and stocking data","volume":"138","author":"Li","year":"2013","journal-title":"Remote Sens. Environ."},{"key":"ref_12","doi-asserted-by":"crossref","first-page":"692","DOI":"10.1016\/j.rse.2010.10.011","article-title":"Analysis of monotonic greening and browning trends from global NDVI time series","volume":"115","author":"Schaepman","year":"2011","journal-title":"Remote Sens. Environ."},{"key":"ref_13","doi-asserted-by":"crossref","first-page":"4217","DOI":"10.3390\/rs6054217","article-title":"Evaluation of spatiotemporal variations of global fractional vegetation cover based on GIMMS NDVI data from 1982 to 2011","volume":"6","author":"Wu","year":"2014","journal-title":"Remote Sens."},{"key":"ref_14","doi-asserted-by":"crossref","first-page":"5247","DOI":"10.1080\/01431161.2010.496470","article-title":"Changes in vegetation spring dates in the second half of the twentieth century","volume":"32","author":"Sobrino","year":"2011","journal-title":"Int. J. Remote Sens."},{"key":"ref_15","doi-asserted-by":"crossref","first-page":"599","DOI":"10.1111\/j.1467-9671.2011.01280.x","article-title":"A contextual mann-kendall approach for the assessment of trend significance in image time series","volume":"15","author":"Neeti","year":"2011","journal-title":"Trans. GIS"},{"key":"ref_16","doi-asserted-by":"crossref","first-page":"8906","DOI":"10.3390\/rs70708906","article-title":"A Temporal-Spatial Iteration Method to Reconstruct NDVI Time Series Datasets","volume":"7","author":"Xu","year":"2015","journal-title":"Remote Sens."},{"key":"ref_17","doi-asserted-by":"crossref","first-page":"833","DOI":"10.1016\/j.cageo.2004.05.006","article-title":"TIMESAT\u2014A program for analyzing time series of satellite sensor data","volume":"30","author":"Eklundh","year":"2004","journal-title":"Comput. Geosci."},{"key":"ref_18","unstructured":"Eklundh, L., and Jonsson, P. (2012). TIMESAT 3. 2 with Parallel Processing\u2013Software Manual, Lund University."},{"key":"ref_19","unstructured":"TIMESAT. A Software Package to Analyse Time-Series of Satellite Sensor Data. Available online: http:\/\/www.nateko.lu.se\/TIMESAT\/."},{"key":"ref_20","unstructured":"Cressie, N., and Wikle, C.K. (2015). Statistics for Spatio-Temporal Data, John Wiley & Sons."},{"key":"ref_21","doi-asserted-by":"crossref","first-page":"265","DOI":"10.1007\/s00704-011-0464-2","article-title":"Spatio-temporal prediction of daily temperatures using time series of MODIS LST images","volume":"107","author":"Hengl","year":"2012","journal-title":"Theor. Appl. Climatol."},{"key":"ref_22","doi-asserted-by":"crossref","first-page":"2108","DOI":"10.3390\/rs6032108","article-title":"Comparison of gross primary productivity derived from GIMMS NDVI3g, GIMMS, and MODIS in Southeast Asia","volume":"6","author":"Wang","year":"2014","journal-title":"Remote Sens."},{"key":"ref_23","doi-asserted-by":"crossref","first-page":"6929","DOI":"10.3390\/rs6086929","article-title":"A non-stationary 1981\u20132012 AVHRR NDVI3g time series","volume":"6","author":"Pinzon","year":"2014","journal-title":"Remote Sens."},{"key":"ref_24","doi-asserted-by":"crossref","first-page":"364","DOI":"10.3390\/rs8050364","article-title":"Impact of Climate Change on Vegetation Growth in Arid Northwest of China from 1982 to 2011","volume":"8","author":"Zhang","year":"2016","journal-title":"Remote Sens."},{"key":"ref_25","doi-asserted-by":"crossref","first-page":"10164","DOI":"10.3390\/rs70810164","article-title":"Effects of precipitation intensity and temperature on NDVI-based grass change over Northern China during the period from 1982 to 2011","volume":"7","author":"Yuan","year":"2015","journal-title":"Remote Sens."},{"key":"ref_26","unstructured":"Ecocast Monitoring, Modeling and Forecasting Ecosystem Change, Available online: http:\/\/ecocast.arc.nasa.gov\/data\/pub\/gimms\/3g.v0\/."},{"key":"ref_27","doi-asserted-by":"crossref","first-page":"3041","DOI":"10.3390\/rs6043041","article-title":"Vegetation greenness in northeastern brazil and its relation to ENSO warm events","volume":"6","author":"Erasmi","year":"2014","journal-title":"Remote Sens."},{"key":"ref_28","unstructured":"R Core Team (2016). R: A Language and Environment for Statistical Computing, R Foundation for Statistical Computing."},{"key":"ref_29","unstructured":"Detsch, F. Gimms: Download and Process GIMMS NDVI3g Data, R Package Version 0.5.1. Available online: https:\/\/cran.r-project.org\/web\/packages\/gimms\/gimms.pdf."},{"key":"ref_30","unstructured":"Hijmans, R.J. Raster: Geographic Data Analysis and Modeling, R Package Version 2.5-2. Available online: https:\/\/cran.r-project.org\/web\/packages\/raster\/raster.pdf."},{"key":"ref_31","unstructured":"CRU TS3.10 Climatic Research Unit. Available online: http:\/\/www.cru.uea.ac.uk\/data."},{"key":"ref_32","doi-asserted-by":"crossref","first-page":"623","DOI":"10.1002\/joc.3711","article-title":"Updated high-resolution grids of monthly climatic observations\u2014The CRU TS3. 10 Dataset","volume":"34","author":"Harris","year":"2014","journal-title":"Int. J.Climatol."},{"key":"ref_33","doi-asserted-by":"crossref","unstructured":"Ugarte, M.D., Militino, A.F., and Arnholt, A.T. (2015). Probability and Statistics with R, CRC Press. [2nd ed.].","DOI":"10.1201\/b18682"},{"key":"ref_34","doi-asserted-by":"crossref","unstructured":"Durbin, J., and Koopman, S.J. (2012). Time Series Analysis by State Space Methods, Oxford University Press.","DOI":"10.1093\/acprof:oso\/9780199641178.001.0001"},{"key":"ref_35","doi-asserted-by":"crossref","first-page":"1027","DOI":"10.1016\/j.envsoft.2009.02.009","article-title":"The EM algorithm in a distributed computing environment for modelling environmental space\u2013time data","volume":"24","author":"Cameletti","year":"2009","journal-title":"Environ. Model. Softw."},{"key":"ref_36","doi-asserted-by":"crossref","first-page":"199","DOI":"10.1016\/S0167-7152(00)00241-8","article-title":"Assessing the covariance function in geostatistics","volume":"52","author":"Militino","year":"2001","journal-title":"Stat. Probab. Lett."},{"key":"ref_37","unstructured":"Cameletti, M. Stem: Spatio-Temporal Models in R, R Package Version 1.0. Available online: https:\/\/cran.r-project.org\/web\/packages\/Stem\/Stem.pdf."},{"key":"ref_38","doi-asserted-by":"crossref","first-page":"209","DOI":"10.5194\/hess-9-209-2005","article-title":"Using a spatio-temporal dynamic state-space model with the EM algorithm to patch gaps in daily riverflow series","volume":"9","author":"Amisigo","year":"2005","journal-title":"Hydrol. Earth Syst. Sci. Discuss."},{"key":"ref_39","doi-asserted-by":"crossref","first-page":"1453","DOI":"10.1002\/joc.4068","article-title":"Interpolation of daily rainfall using spatiotemporal models and clustering","volume":"35","author":"Militino","year":"2015","journal-title":"Int. J. Climatol."},{"key":"ref_40","doi-asserted-by":"crossref","first-page":"1","DOI":"10.1111\/j.2517-6161.1977.tb01600.x","article-title":"Maximum likelihood from incomplete data via the EM algorithm","volume":"39","author":"Dempster","year":"1977","journal-title":"J. R. Stat. Soc. Ser. B B (Methodol.)"},{"key":"ref_41","doi-asserted-by":"crossref","first-page":"17","DOI":"10.1093\/biomet\/37.1-2.17","article-title":"Notes on continuous stochastic phenomena","volume":"37","author":"Moran","year":"1950","journal-title":"Biometrika"},{"key":"ref_42","doi-asserted-by":"crossref","first-page":"1","DOI":"10.18637\/jss.v063.i12","article-title":"Model-Based geostatistics the easy way","volume":"63","author":"Brown","year":"2015","journal-title":"J. Stat. Softw."},{"key":"ref_43","unstructured":"AEMET Agencia Estatal de Meteorolog\u00eda. Available online: http:\/\/www.aemet.es\/es\/portada."},{"key":"ref_44","doi-asserted-by":"crossref","first-page":"1417","DOI":"10.1080\/01431168608948945","article-title":"Characteristics of maximum-value composite images from temporal AVHRR data","volume":"7","author":"Holben","year":"1986","journal-title":"Int. J. Remote Sens."},{"key":"ref_45","doi-asserted-by":"crossref","first-page":"1886","DOI":"10.1016\/j.rse.2009.04.004","article-title":"Evaluation of earth observation based long term vegetation trends\u2014Intercomparing NDVI time series trend analysis consistency of Sahel from AVHRR GIMMS, Terra MODIS and SPOT VGT data","volume":"113","author":"Fensholt","year":"2009","journal-title":"Remote Sens. Environ."},{"key":"ref_46","doi-asserted-by":"crossref","first-page":"2701","DOI":"10.1016\/j.rse.2009.08.008","article-title":"Impact of orthorectification and spatial sampling on maximum NDVI composite data in mountain regions","volume":"113","author":"Fontana","year":"2009","journal-title":"Remote Sens. Environ."},{"key":"ref_47","doi-asserted-by":"crossref","first-page":"3827","DOI":"10.1080\/01431160010007033","article-title":"Spatial patterns of NDVI in response to precipitation and temperature in the central Great Plains","volume":"22","author":"Wang","year":"2001","journal-title":"Int. J. Remote Sens."},{"key":"ref_48","doi-asserted-by":"crossref","first-page":"3873","DOI":"10.1080\/01431160110119416","article-title":"Global correlation analysis for NDVI and climatic variables and NDVI trends: 1982\u20131990","volume":"23","author":"Ichii","year":"2002","journal-title":"Int. J. Remote Sens."},{"key":"ref_49","doi-asserted-by":"crossref","first-page":"277","DOI":"10.1016\/0273-1177(93)90559-T","article-title":"Global correlation of temperature, NDVI and precipitation","volume":"13","author":"Schultz","year":"1993","journal-title":"Adv. Space Res."},{"key":"ref_50","doi-asserted-by":"crossref","first-page":"2921","DOI":"10.1080\/014311698214352","article-title":"Global analysis of empirical relations between annual climate and seasonality of NDVI","volume":"19","author":"Potter","year":"1998","journal-title":"Int. J. Remote Sens."},{"key":"ref_51","doi-asserted-by":"crossref","first-page":"2057","DOI":"10.1080\/01431161003762363","article-title":"Temporal analysis of normalized difference vegetation index (NDVI) and land surface temperature (LST) parameters to detect changes in the Iberian land cover between 1981 and 2001","volume":"32","author":"Julien","year":"2011","journal-title":"Int. J. Remote Sens."}],"container-title":["Remote Sensing"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/www.mdpi.com\/2072-4292\/9\/1\/76\/pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2025,10,11]],"date-time":"2025-10-11T18:26:16Z","timestamp":1760207176000},"score":1,"resource":{"primary":{"URL":"https:\/\/www.mdpi.com\/2072-4292\/9\/1\/76"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2017,1,15]]},"references-count":51,"journal-issue":{"issue":"1","published-online":{"date-parts":[[2017,1]]}},"alternative-id":["rs9010076"],"URL":"https:\/\/doi.org\/10.3390\/rs9010076","relation":{},"ISSN":["2072-4292"],"issn-type":[{"value":"2072-4292","type":"electronic"}],"subject":[],"published":{"date-parts":[[2017,1,15]]}}}