Land surface temperature (LST) is a basic parameter in energy exchange between the land and the atmosphere, and is frequently used in many sciences such as climatology, hydrology, agriculture, ecology, etc. Time series of satellite LST data have usually deficient, missing, and unacceptable data caused by the presence of clouds in images, the presence of dust in the atmosphere, and sensor failure. In this study, the singular spectrum analysis (SSA) algorithm was used to resolve the problem of missing and outlier data caused by cloud cover. The region studied in the present research included an image frame of the Moderate Resolution Imaging Spectroradiometer (MODIS) with horizontal number 22 and vertical number 05 (h22v05). This image involved a large part of Iran, Turkmenistan, and the Caspian Sea. In this study, MODIS LST products (MOD11A1) were used during 2015 with approximately 1 km × 1 km spatial resolution and day/night LST data (daily temporal resolution). On average, the data have 36.37% gaps in each pixel profile with 730 day/night LST data. The results of the SSA algorithm in the reconstruction of LST images indicated a root mean square error (RMSE) of 2.95 Kelvin (K) between the original and reconstructed LST time series data in the study region. In general, the findings showed that the SSA algorithm using spatio-temporal interpolation can be effectively used to resolve the problem of missing data caused by cloud cover.

Monitoring vegetation changes over time is very important in dry areas such as Iran, given its pronounced drought-prone agricultural system. Vegetation indices derived from remotely sensed satellite imageries are successfully used to monitor vegetation changes at various scales. Atmospheric dust as well as airborne particles, particularly gases and clouds, significantly affect the reflection of energy from the surface, especially in visible, short and infrared wavelengths. This results in imageries with missing data (gaps) and outliers while vegetation change analysis requires integrated and complete time series data. This study investigated the performance of HANTS (Harmonic ANalysis of Time Series) algorithm and (M)-SSA ((Multi-channel) Singular Spectrum Analysis) algorithm in reconstruction of wide-gap of missing data. The time series of Normalized Difference Vegetation Index (NDVI) retrieved from Landsat TM in combination with 250m MODIS NDVI time image products are used to simulate and find periodic components of the NDVI time series from 1986 to 2000 and from 2000 to 2015, respectively. This paper presents the evaluation of the performance of gap filling capability of HANTS and M-SSA by filling artificially created gaps in data using Landsat and MODIS data. The results showed that the RMSEs (Root Mean Square Errors) between the original and reconstructed data in HANTS and M-SSA algorithms were 0.027 and 0.023 NDVI value, respectively. Further, RMSEs among 15 NDVI images extracted from the time series artificially and reconstructed by HANTS and M-SSA algorithms were 0.030 and 0.025 NDVI value, respectively. RMSEs of the original and reconstructed data in HANTS and M-SSA algorithms were 0.10 and 0.04 for time series 6, respectively. The findings of this study present a favorable option for solving the missing data challenge in NDVI time series.

Water scarcity is a critical abiotic stress factor for plants in arid and semi-arid regions, impacting crop development and production yield and quality. Monitoring water stress at finer scales (e.g., farm and plant), requires multispectral imagery with thermal capabilities at centimeter resolution. This study investigates drought stress in pistachio trees in a farm located in Yazd province, Iran, by using Unmanned Aerial Vehicle (UAV) images to quantify evapotranspiration and assess drought stress in individual trees. Images were captured on 10 July 2022, using a Matrix 300 UAV with a MicaSense Altum multispectral sensor. By employing the Surface Energy Balance Algorithm for Land (SEBAL), actual field evapotranspiration was accurately calculated (10 cm spatial resolution). Maps of the optimum crop coefficient (K_c) were developed from the Normalized Difference Vegetation Index (NDVI) based on standard evapotranspiration using the Food and Agriculture Organization (FAO) 56 methodology. The comparison between actual and standard evapotranspiration allowed us to identify drought-stressed trees. Results showed an average and maximum daily evaporation of 4.3 and 8.0 mm/day, respectively, in pistachio trees. The real crop coefficient (K_c) for pistachio was 0.66, contrasting with the FAO 56 standard of 1.17 due to the stress factor (K_s). A significant correlation was found between K_c and NDVI (R² = 0.67, p < 0.01). The regression model produced a crop coefficient map, valuable to support precise irrigation management and drought prevention, considering the heterogeneity at the farm scale.

Climate change affects plant dynamics and functioning of terrestrial ecosystems. This study aims to investigate temporal changes in global vegetation coverage and biomes during the past three decades. We compared historic annual NDVI time series (1982, 1983, 1984 and 1985) with recent ones (2015, 2016, 2017 and 2018), captured from NOAA-AVHRR satellite observations. To correct the NDVI time series for missing data and outliers, we applied the Harmonic Analysis of Time Series (HANTS) algorithm. The NDVI time series were decomposed in their significant amplitude and phase given their periodic fluctuation, except for ever green vegetation. Our findings show that the average NDVI values in most biomes have increased significantly (F-value<0.01) by 0.05 ndvi units over during the past three decades, except in tundra, and deserts and xeric shrublands. The highest rates of change in the harmonic components were observed in the northern hemisphere, mainly above 30° latitude. Worldwide, the mean annual phase reduced by 9° corresponding to a 9 days shift in the beginning of the growing season. Annual phases in the recent time series reduced significantly as compared to the historic time series in the five major global biomes: by 14.1, 14.8, 10.6, 9.5, and 22.8 days in boreal forests/taiga; Mediterranean forests, woodlands, and scrubs; temperate conifer forests; temperate grasslands, savannas, and shrublands; and deserts, and xeric shrublands, respectively. In tropical and subtropical biomes, however, changes in the annual phase of vegetation coverage were not statistically significant. The decrease in the level of phases and acceleration of growth and changes in plant phenology indicate the increase in temperature and climate changes of the planet.

Investigating vegetation changes, especially plant phenology, can yield valuable information about global warming and climate change. Time series satellite observations and remote sensing methods offer a great source of information on distinctions and changing aspects of vegetation. The current study aimed to determine the trend and rate of changes in some phenological components of Iran’s vegetation. In this regard, the current study employed the daily NDVI (Normalized Difference Vegetation Index) product of the AVHRR sensor with a spatial resolution of 0.05° × 0.05°, named AVH13C1. Then, using the HANTS algorithm, images of amplitude zero, annual amplitude, and annual phase were prepared annually from 1982 to 2019. Using TIMESAT software, the starting, end, and length of time of growing season were calculated for each pixel time series to prepare annual maps. The Mann–Kendall statistical test was used to investigate the significance of changes during the study period. On average in the entire area of Iran, the annual phase was declining with a trend of −0.6° per year, and the time for the start and end of the season was declining by −0.3 and −0.65 days per year, respectively. Major changes were noticed in the northeast, west, and northwest regions of Iran, where the annual phase declined with a trend of −0.9° per year. Since the annual growth cycle of the plant (equivalent to 356 days) was in the form of a sinusoidal signal, and the angular changes in the sine wave were between zero and 360°, each degree of change was equivalent to 1.01 days per year. Therefore, the reduction in the annual phase by −0.9 degrees almost means a change in the time (due to the earlier negative start phase) of the start of the annual growth signal by −0.9 days per year. The time of the start and end of the growing season declined by −0.6 and −1.33 days per year, respectively. The reduction in annual phase and differences in time of the starting season from 1982 to 2019 indicate the acceleration and earlier initiation of various phenological processes in the area.

فنولوژی به مطالعة زمان وقوع رخدادهای تکرارپذیر زندگی گیاهان در رابطه با عوامل زنده و غیرزنده می‌پردازد. یکی از روش‌های بررسی تغییرات روند رشد گیاه استفاده از سری‌های زمانی تصاویر ماهواره‌ای و روش‌های سنجش از دور است. هدف از پژوهش حاضر تعیین معنای فیزیکی هریک از مؤلفه‌های سری فوریه، درمورد انواع پوشش گیاهی و فرایندهای فنولوژی گیاهان، با استفاده از تجزیة سری‌های زمانی شاخص پوشش گیاهی (NDVI) حاصل از سنجش از دور است. هدف دیگر نیز تعیین روند تغییرات مؤلفه‌های سری فوریه، در زمینة تغییرات ارتفاع از سطح دریا و درجة حرارت فصل بهار است. بدین‌منظور، در مطالعة حاضر، از محصول هشت‌روزة شاخص‌ NDVI سنجندة MODIS، با نام MOD09Q1، با قدرت تفکیک مکانی 250 متر استفاده شد. ابتدا، با استفاده از الگوریتم HANTS و سری‌های زمانی‌ یک‌سالة NDVI سنجندة MODIS در سال‌های 2013، 2014 و 2015، مؤلفه‌های سری‌های فوریه یا همان تصاویر دامنه و فاز تولید شد. سپس تغییرات و ارتباط هریک از این مؤلفه‌ها به‌نسبت تغییرات ارتفاع و دما، در منطقة شیرکوه استان یزد، بررسی شد. براساس نتایج، با کاهش 1 درجة سانتی‌گراد در میانگین دمای فصل بهار، که با افزایش ارتفاع از سطح دریا رخ می‌دهد، تأخیر هفت‌روزه‌ای در...