Satellite detection of increases in global land surface evapotranspiration during 1984–2007

Abstract

As a key component of digital earth, remotely sensed data provides the compelling evidence that the amount of water vapour transferred from the entire global surface to the atmosphere increased from 1984 to 2007. The validation results from the earlier evapotranspiration (ET) estimation algorithm based on net radiation (R _n ), Normalised Difference Vegetation Index (NDVI), air temperature and diurnal air temperature range (DTaR) showed good agreement between estimated monthly ET and ground-measured ET from 20 flux towers. Our analysis indicates that the estimated actual ET has increased on average over the entire global land surface except for Antarctica during 1984–2007. However, this increasing trend disappears after 2000 and the reason may be that the decline in net radiation and NDVI during this period depleted surface soil moisture. Moreover, the good correspondence between the precipitation trend and the change in ET in arid and semi-arid regions indicated that surface moisture linked to precipitation affects ET. The input parameters R _n , T _air, NDVI and DTaR show substantial spatio-temporal variability that is almost consistent with that of actual ET from 1984 to 2007 and contribute most significantly to the variation in actual ET.     

HJ-1卫星数据估算地表能量平衡

利用遥感数据估算和监测地表能量平衡的研究和应用一直以来都是学术探索的前沿和焦点,目前针对国际主流卫星数据(如MODIS,Landsat TM等)的算法研究有很多,且在不同的尺度上发展了一系列遥感技术和应用方法。本文以环境一号卫星(HJ-1)数据为主要数据源,研究和分析了地表能量平衡遥感估算方法。在下行短波和长波辐射,以及地表净辐射估算的基础上,重点研究了地表显热和潜热通量(蒸散)的估算方法,针对辐射温度代替空气动力学温度所引起的误差校正,对比分析了基于KB-1系数的热力学粗糙度长度,以及附加阻抗两种方法。使用2010年HJ-1B数据反演得到海河流域地表能量平衡各分量瞬时值,利用3个地面站点的测量数据对反演结果进行验证和误差分析。结果表明下行短波辐射的反演误差是地表净辐射反演误差的主要来源,地表温度的反演精度以及地表粗糙度的模拟精度是影响显热通量反演精度的主要因素。利用两种方法估算的显热通量趋势基本一致,但KB-1系数法的反演结果一般低于附加阻抗法。在下垫面较为复杂的区域,模型结果误差较大,需要更加精细的模型以刻画地表非均匀性的影响。     

Examining the impacts of urbanization on surface radiation using Landsat imagery

Metropolitan Beijing is facing many environmental problems such as haze and urban heat island due to the rapid urbanization. Surface shortwave, longwave, and net radiations are key components of the surface-atmosphere radiation budget. Since megacities are affected by the thermal radiation of complex landscape structures and atmospheric environments, quantitative and spatially explicit retrieval from remotely sensed data remains a challenge. We collected the surface radiation fluxes from seven fixed sites representing different land-use types to calibrate the local parameters for remotely sensed retrieval of net radiation. We proposed a remote sensing–based surface radiation retrieval method by embedding the underlying land covers and integrating the observational data. The improved method is feasible to accurately retrieve surface radiation and delineate spatial characteristics in metropolitan areas. The accuracy evaluation indicated that the difference between remotely sensed and in situ observed net radiation ranged within 0~± 40 W· m^?2. The root mean squared error of the estimated net surface radiation was 32.71 W· m^?2. The strongly spatial heterogeneity of surface radiation components in metropolitan Beijing was closely related to land-cover patterns from urban area to outskirts. We also found that the surface net radiation had a decreasing trend from 1984 to 2014, and the net radiation in the urban area was lower than that in the outskirts. According to the surface radiation budgets, urbanization resulted in the cooling effect in net radiation flux in the daytime, which was stemmed from low atmospheric transmittances from massive aerosol concentration and high surface albedo from light building materials.     

Estimation of daily average net radiation and its variation over West Bengal,India using MODIS products

A simple and exclusively satellite-based instantaneous net radiation (INR) and its components have been estimated and it is based on the method proposed. The method has been used to check its suitability over West Bengal, India on clear sky days which is eliminating the need for ground-based observations. The method explicitly recognizes the need for spatially varied input parameters to estimate net radiation using remote sensing information and provides a distributed daily average net radiation map with finer spatial resolution. The satellite-derived radiative flux components were compared with the Eddy Flux Tower measurements located at Bonnie camp within mangrove forest. The bias, root-mean square error and R² between modelled and measured INR are 52.07, 49.5 W m^?2 and 0.74, respectively which are within the acceptable limit. The errors in the downward short wave and long wave as well as upward short-wave and long-wave radiative flux are 13.78, 3.54, 12.41 and 3.54%, respectively. Further, variations of all the radiative flux components were analysed over different eco-systems.     

Winter wheat biomass estimation using high temporal and spatial resolution satellite data combined with a light use efficiency model

Winter wheat biomass was estimated using HJ CCD and MODIS data, combined with a radiation use efficiency model. Results were validated with ground measurement data. Winter wheat biomass estimated with HJ CCD data correlated well with observed biomass in different experiments (coefficients of determination R² of 0.507, 0.556 and 0.499; n?=?48). In addition, R² values between MODIS estimated and observed biomass are 0.420, 0.502 and 0.633. Even if we downscaled biomass estimated using HJ CCD data to MODIS pixel size (9?×?9 HJ CCD pixels to approximate that MODIS pixel), R² values between estimated and observed biomass were still higher than those from MODIS. We conclude that estimation with remote sensing data, such as the HJ CCD data with high spatial resolution and shorter revisit cycle, can show more detail in spatial pattern and improve the application of remote sensing on a local scale. There is also potential for applying the approach to many other studies, including agricultural production estimation, crop growth monitoring and agricultural ecosystem carbon cycle studies.     

The Generalised Settlement Area: mapping the Earth surface in the vicinity of built-up areas

ABSTRACT

Geo-information on settlements from Earth Observation offers a base for objective and scalable monitoring of the evolution of cities and settlements, including their location, extent and other attributes. In this work, we deploy the best available global knowledge on the presence of human settlements and built-up structures derived from Earth Observation to advance the understanding of the human presence on Earth. We start from a concept of Generalised Settlement Area to identify the Earth surface within which any built-up structure is present. We further characterise the resulted map by using an agreement map among the state of the art of remote sensing products mapping built-up areas or other strictly related semantic abstractions as urban areas or artificial surfaces. The agreement map is formed by a grid of 1?km², where each cell is classified according to the number of EO-derived products reporting any positive occurrence of the abstractions related to the presence of built-up structures. The paper describes the characteristics of the Generalised Settlement Area, the differences in the agreement map across geographic regions of the world, and outlines the implications for potential users of the EO-derived products used in this study.     

基于葵花-8卫星大气产品的地表下行短波辐射计算

地表下行短波辐射DSSR(Downward Surface Shortwave Radiation)的准确估算在气候变化研究和地表太阳能估算等领域具有重要作用。新一代静止气象卫星葵花-8(Himawari-8)具有高达10 min的对地观测能力,为DSSR近实时估算提供了新机遇。然而,日本宇宙航空研究开发机构(JAXA)对外公开的葵花-8辐射产品中,没有将其反演的云、气溶胶产品作为DSSR的输入参数,从而没有形成一整套的DSSR估算算法流程,缺乏产品输出的一致性。大气中的云、气溶胶是DSSR的重要影响因子,本文重点考虑云、气溶胶对太阳辐射的影响,基于大气辐射传输模式RSTAR构建了DSSR查找表,开发了DSSR的快速计算方法,进而将JAXA葵花-8二级云、气溶胶产品(光学厚度,粒子有效半径等)作为快速化计算方法的输入参量,计算得到了DSSR。通过与JAXA葵花-8二级DSSR产品(JAXA DSSR)对比,发现两者具有很好的空间一致性。为了进一步评价本文的DSSR计算精度,分别选取了陆地(Yonsei)和海洋(0n_165e)的观测数据验证了2016年4、7、10和12月本文计算的DSSR和同时期的JAXA DSSR产品,验证结果显示两者的DSSR在两个观测站点均具有非常高的相关性(全天空、晴空和云天条件下的相关系数R均大于0.88)。在两个站点云天条件下的验证结果中,考虑了云相态并在冰云模型中使用了非球形冰晶粒子(六棱柱)来计算DSSR,获得了比JAXA DSSR更小的偏差。本文提出的快速化计算方法能快速准确地计算DSSR,可为计算地表辐射收支等研究提供重要数据支撑。     

Assessment and planning for integrated river basin management using remote sensing,SWAT model and morphometric analysis (case study: Kaddam river basin,India)

Abstract

River basin assessment is crucial for water management and to address the watershed issues. So, an integrated river basin management and assessment model using morphometric assessment, remote sensing, GIS and SWAT model was envisaged and applied to Kaddam river basin, Telangana state, India. Morphometric results showed high drainage density ranging from 2.19 to 5.5?km²/km, with elongated fan shape having elongation ratio of 0.60–0.75 with sparse vegetation and high relief. Land use change assessment showed that 265.26?km² of forest land is converted into irrigated land and has increased sediment yields in watersheds. The calibration (r ²?=?0.74, NSE?=?0.84) and validation (r ²?=?0.72, NSE?=?0.84) of SWAT model showed that simulated and observed results were in agreement and in recommended ranges. The SWAT simulations were used to compute mean annual water and sediment yield from 1997 to 2012, along with morphometric results to categorize critical watersheds and conservation structures were proposed accordingly.     

Low-degree earth deformation from reprocessed GPS observations

Surface mass variations of low spherical harmonic degree are derived from residual displacements of continuously tracking global positioning system (GPS) sites. Reprocessed GPS observations of 14 years are adjusted to obtain surface load coefficients up to degree n _max = 6 together with station positions and velocities from a rigorous parameter combination. Amplitude and phase estimates of the degree-1 annual variations are partly in good agreement with previously published results, but also show interannual differences of up to 2 mm and about 30 days, respectively. The results of this paper reveal significant impacts from different GPS observation modeling approaches on estimated degree-1 coefficients. We obtain displacements of the center of figure (CF) relative to the center of mass (CM), Δr _CF–CM, that differ by about 10 mm in maximum when compared to those of the commonly used coordinate residual approach. Neglected higher-order ionospheric terms are found to induce artificial seasonal and long-term variations especially for the z-component of Δr _CF–CM. Daily degree-1 estimates are examined in the frequency domain to assess alias contributions from model deficiencies with regard to satellite orbits. Finally, we directly compare our estimated low-degree surface load coefficients with recent results that involve data from the Gravity Recovery and Climate Experiment (GRACE) satellite mission.     

A comparison of geographic datasets and field measurements to model soil carbon using random forests and stepwise regressions (British Columbia,Canada)

We used geographic datasets and field measurements to examine the mechanisms that affect soil carbon (SC) storage for 65 grazed and non-grazed pastures in southern interior grasslands of British Columbia, Canada. Stepwise linear regression (SR) modeling was compared with random forest (RF) modeling. Models produced with SR performed better than those produced using RF models (r² = 0.56–0.77 AIC = 0.16–0.30 for SR models; r² = 0.38–0.53 and AIC = 0.18–0.30 for RF models). The factors most significant when predicting SC were elevation, precipitation, and the normalized difference vegetation index (NDVI). NDVI was evaluated at two scales using: (1) the MOD 13Q1 (250 m/16-day resolution) NDVI data product from the moderate resolution imaging spectro-radiometer (MODIS) (NDVI_MODIS), and (2) a handheld multispectral radiometer (MSR, 1 m resolution) (NDVI_MSR) in order to understand the potential for increasing model accuracy by increasing the spatial resolution of the gridded geographic datasets. When NDVI_MSR data were used to predict SC, the percentage of the variance explained by the model was greater than for models that relied on NDVI_MODIS data (r² = 0.68 for SC for non-grazed systems, modeled with SR based on NDVI_MODIS data; r² = 0.77 for SC for non-grazed systems, modeled with SR based on NDVI_MSR data). The outcomes of this study provide the groundwork for effective monitoring of SC using geographic datasets to enable a carbon offset program for the ranching industry.     

Estimation of winter wheat crop growth parameters using time series Sentinel-1A SAR data

In the present study, Sentinel-1A Synthetic Aperture Radar analysis of time series data at C-band was carried out to estimate the winter wheat crop growth parameters. Five different date images were acquired during January 2015–April 2015 at different growth stages from tillering to ripening in Varanasi district, India. The winter wheat crop parameters, i.e. leaf area index, vegetation water content (VWC), fresh biomass (FB), dry biomass (DB) and plant height (PH) were estimated using random forest regression (RFR), support vector regression (SVR), artificial neural network regression (ANNR) and linear regression (LR) algorithms. The Ground Range Detected products of Interferometric Wide (IW) Swath were used at VV polarization. The three different subplots of 1 m² area were taken for the measurement of crop parameters at every growth stage. In total, 73 samples were taken as the training data-sets and 39 samples were taken as testing data-sets. The highest sensitivity (adj. R^2?=?0.95579) of backscattering with VWC was found using RFR algorithm, whereas the lowest sensitivity (adj. R^2?=?0.66201) was found for the PH using LR algorithm. Overall results indicate more accurate estimation of winter wheat parameters by the RFR algorithm followed by SVR, ANNR and LR algorithms.     

Land Surface Emissivity and Temperature Determination Using NOAA‐AVHRR Data

Abstract

Studies on land surface processes using remote sensing data gains importance in the context of Geosphere Biosphere Programme. Present study addresses the applicability of split‐window method, in a tropical environment for mapping of surface temperature over heterogeneous surface from satellite data. The accuracy of the method is about +2.2°K, which is reasonable value taking into account the atmospheric attenuation in tropical environment. An attempt has been made to derive emissivity from normalized difference vegetation index (NDVI) by taking into account the fraction of vegetation cover of each pixel, which is determined by satellite data. The emissivity values estimated from satellite data found to be in reasonable agreement with an estimated error of less than 1%. The results of the study indicate the potential use of NDVI as a modulating parameter in the land surface temperature estimation from satellite data.     

Carbon fluxes in ecosystems of Yellowstone National Park predicted from remote sensing data and simulation modeling

Background

A simulation model based on remote sensing data for spatial vegetation properties has been used to estimate ecosystem carbon fluxes across Yellowstone National Park (YNP). The CASA (Carnegie Ames Stanford Approach) model was applied at a regional scale to estimate seasonal and annual carbon fluxes as net primary production (NPP) and soil respiration components. Predicted net ecosystem production (NEP) flux of CO₂ is estimated from the model for carbon sinks and sources over multi-year periods that varied in climate and (wildfire) disturbance histories. Monthly Enhanced Vegetation Index (EVI) image coverages from the NASA Moderate Resolution Imaging Spectroradiometer (MODIS) instrument (from 2000 to 2006) were direct inputs to the model. New map products have been added to CASA from airborne remote sensing of coarse woody debris (CWD) in areas burned by wildfires over the past two decades.

Results

Model results indicated that relatively cooler and wetter summer growing seasons were the most favorable for annual plant production and net ecosystem carbon gains in representative landscapes of YNP. When summed across vegetation class areas, the predominance of evergreen forest and shrubland (sagebrush) cover was evident, with these two classes together accounting for 88% of the total annual NPP flux of 2.5 Tg C yr^-1 (1 Tg = 10¹² g) for the entire Yellowstone study area from 2000-2006. Most vegetation classes were estimated as net ecosystem sinks of atmospheric CO₂ on annual basis, making the entire study area a moderate net sink of about +0.13 Tg C yr^-1. This average sink value for forested lands nonetheless masks the contribution of areas burned during the 1988 wildfires, which were estimated as net sources of CO₂ to the atmosphere, totaling to a NEP flux of -0.04 Tg C yr^-1 for the entire burned area. Several areas burned in the 1988 wildfires were estimated to be among the lowest in overall yearly NPP, namely the Hellroaring Fire, Mink Fire, and Falls Fire areas.

Conclusions

Rates of recovery for burned forest areas to pre-1988 biomass levels were estimated from a unique combination of remote sensing and CASA model predictions. Ecosystem production and carbon fluxes in the Greater Yellowstone Ecosystem (GYE) result from complex interactions between climate, forest age structure, and disturbance-recovery patterns of the landscape.     

Spatial and Temporal Variability of Net Primary Productivity (NPP) over Terrestrial Biosphere of India Using NOAA-AVHRR Based GloPEM Model

The monitoring of terrestrial carbon dynamics is important in studies related with global climate change. This paper presents results of the inter-annual variability of Net Primary Productivity (NPP) from 1981 to 2000 derived using observations from NOAA-AVHRR data using Global Production Efficiency Model (GloPEM). The GloPEM model is based on physiological principles and uses the production efficiency concept, in which the canopy absorption of photosynthetically active radiation (APAR) is used with a conversion “efficiency” to estimate Gross Primary Production (GPP). NPP derived from GloPEM model over India showed maximum NPP about 3,000 gCm⁻²year⁻¹ in west Bengal and lowest up to 500 gCm⁻²year⁻¹ in Rajasthan. The India averaged NPP varied from 1,084.7 gCm⁻²year⁻¹ to 1,390.8 gCm⁻²year⁻¹ in the corresponding years of 1983 and 1998 respectively. The regression analysis of the 20 year NPP variability showed significant increase in NPP over India (r = 0.7, F = 17.53, p < 0.001). The mean rate of increase was observed as 10.43 gCm⁻²year⁻¹. Carbon fixation ability of terrestrial ecosystem of India is increasing with rate of 34.3 TgC annually (t = 4.18, p < 0.001). The estimated net carbon fixation over Indian landmass ranged from 3.56 PgC (in 1983) to 4.57 PgC (in 1998). Grid level temporal correlation analysis showed that agricultural regions are the source of increase in terrestrial NPP of India. Parts of forest regions (Himalayan in Nepal, north east India) are relatively less influenced over the study period and showed lower or negative correlation (trend). Finding of the study would provide valuable input in understanding the global change associated with vegetation activities as a sink for atmospheric carbon dioxide.     

A practical algorithm for estimating surface soil moisture using combined optical and thermal infrared data

Surface soil moisture (SSM) is a critical variable for understanding the energy and water exchange between the land and atmosphere. A multi-linear model was recently developed to determine SSM using ellipse variables, namely, the center horizontal coordinate (x₀), center vertical coordinate (y₀), semi-major axis (a) and rotation angle (θ), derived from the elliptical relationship between diurnal cycles of land surface temperature (LST) and net surface shortwave radiation (NSSR). However, the multi-linear model has a major disadvantage. The model coefficients are calculated based on simulated data produced by a land surface model simulation that requires sufficient meteorological measurements. This study aims to determine the model coefficients directly using limited meteorological parameters rather than via the complicated simulation process, decreasing the dependence of the model coefficients on meteorological measurements. With the simulated data, a practical algorithm was developed to estimate SSM based on combined optical and thermal infrared data. The results suggest that the proposed approach can be used to determine the coefficients associated with all ellipse variables based on historical meteorological records, whereas the constant term varies daily and can only be determined using the daily maximum solar radiation in a prediction model. Simulated results from three FLUXNET sites over 30 cloud-free days revealed an average root mean square error (RMSE) of 0.042 m³/m³ when historical meteorological records were used to synchronously determine the model coefficients. In addition, estimated SSM values exhibited generally moderate accuracies (coefficient of determination R² = 0.395, RMSE = 0.061 m³/m³) compared to SSM measurements at the Yucheng Comprehensive Experimental Station.     

Quantification of LST and CO2 levels using Landsat-8 thermal bands on urban environment

Land surface temperature (LST) is an important aspect in global to regional change studies, for control of climate change and balancing of high temperature. Urbanization is one of the influencing factors increasing land surface and atmospheric temperature, by the emission of greenhouse gases (e.g. CO₂, NO and methane). In the present study, LST was derived from Landsat-8 of multitemporal data sets to analyse the spatial structure of the urban thermal environment in relation to the urban surface characteristics and land use–land cover (LULC). LST is influenced by the greenhouse gases i.e. CO₂ plays an important role in increasing the earth’s surface temperature. In order to provide the evidence of influence of CO₂ on LST, the relationship between LST, air temperature and CO₂ was analysed. Landsat-8 satellite has two thermal bands, 10 and 11. These bands were used to accurately to calculate the temperature over the study area. Results showed that the strength of correlation between ground monitoring data and satellite data was high. Based on correlation values of each month April (R² = 0.994), May (R² = 0.297) and June (R² = 0.934), observed results show that band 10 was significantly correlating with air temperature. Relationship between LST and CO₂ levels were obtained from linear regression analysis. band 11 was correlating significantly with CO₂ values in each of the months April (R² = 0.217), May (R² = 0.914) and June, (R² = 0.934), because band 11 is closer to the 15-micron band of CO₂. From the results, it was observed that band 10 can be used for calculating air temperature and band 11 can be used for estimation of greenhouse gases.     

Comparing spatiotemporal patterns in Eurasian FPAR derived from two NDVI-based methods

Abstract

A long-term, consistent Fraction of Absorbed Photosynthetically Active Radiation (FPAR) product is necessary to study the spatial and temporal patterns of vegetation dynamics associated with climatic changes and human activities. In this study, Eurasia was selected as the study area. The relationship between FPAR and simple infrared/red ratio relationship (SR FPAR), and that between Moderate Resolution Imaging Spectroradiometer (MODIS) FPAR and a Normalised Difference Vegetation Index (NDVI) look-up table (LUT FPAR) were employed to estimate FPAR from 1982 to 2006 by different land cover types, focusing on the comparisons of spatiotemporal FPAR patterns between the two FPAR datasets. The results showed high agreement between MODIS standard FPAR and estimated FPAR in seasonal dynamics with peak values in July. The LUT FPAR was close to MODIS standard FPAR and larger than SR FPAR. The SR and LUT FPAR showed the same spatial distribution and inter-annual variation patterns and were primarily determined by land cover types. An overall increasing trend in FPAR was observed from 1982 to 2006, with reductions from 1991 to 1994 and 2000 to 2002. The inter-annual dynamics in evergreen broadleaf forests showed a decreasing trend over 25 years, while non-forest vegetation FPAR values had slow, stable growth in inter-annual variation.     

Mixed Integer-Real Valued Adjustment (IRA) Problems: GPS Initial Cycle Ambiguity Resolution by Means of the LLL Algorithm

In order to achieve to GPS solutions of first-order accuracy and integrity, carrier phase observations as well as pseudorange observations have to be adjusted with respect to a linear/linearized model. Here the problem of mixed integer-real valued parameter adjustment (IRA) is met. Indeed, integer cycle ambiguity unknowns have to be estimated and tested. At first we review the three concepts to deal with IRA: (i) DDD or triple difference observations are produced by a properly chosen difference operator and choice of basis, namely being free of integer-valued unknowns (ii) The real-valued unknown parameters are eliminated by a Gauss elimination step while the remaining integer-valued unknown parameters (initial cycle ambiguities) are determined by Quadratic Programming and (iii) a RA substitute model is firstly implemented (real-valued estimates of initial cycle ambiguities) and secondly a minimum distance map is designed which operates on the real-valued approximation of integers with respect to the integer data in a lattice. This is the place where the integer Gram-Schmidt orthogonalization by means of the LLL algorithm (modified LLL algorithm) is applied being illustrated by four examples. In particular, we prove that in general it is impossible to transform an oblique base of a lattice to an orthogonal base by Gram-Schmidt orthogonalization where its matrix enties are integer. The volume preserving Gram-Schmidt orthogonalization operator constraint to integer entries produces “almost orthogonal” bases which, in turn, can be used to produce the integer-valued unknown parameters (initial cycle ambiguities) from the LLL algorithm (modified LLL algorithm). Systematic errors generated by “almost orthogonal” lattice bases are quantified by A. K. Lenstra et al. (1982) as well as M. Pohst (1987). The solution point of Integer Least Squares generated by the LLL algorithm is = (L')⁻¹[L'◯] ∈ ℤ ^m where L is the lower triangular Gram-Schmidt matrix rounded to nearest integers, [L], and = [L'◯] are the nearest integers of L'◯, ◯ being the real valued approximation of z ∈ ℤ ^m , the m-dimensional lattice space Λ. Indeed due to “almost orthogonality” of the integer Gram-Schmidt procedure, the solution point is only suboptimal, only close to “least squares.” ? 2000 John Wiley & Sons, Inc.     

Multitemporal Landsat-MODIS fusion for cropland drought monitoring in El Salvador

Abstract

This study aims to develop an approach to characterize cropland drought conditions in El Salvador, Central America. The data were processed for 2016–2017 through three main steps: (1) reconstructing MODIS land-surface temperature (LST), (2) Landsat-MODIS data fusion and (3) drought delineation using the temperature vegetation dryness index (TVDI). The results of LST reconstruction using the random forests (RF) indicated the median RMSE value of 0.5?°C. The fusion results achieved from the STARFM compared with the reference Landsat data revealed close agreement with the correlation coefficient (r) values higher than 0.84. The TVDI results verified with that from the reference Landsat data indicated r values of 0.85 and 0.75 for 2016 and 2017, respectively. The larger very dry area was observed for the 2016 primera season due to prolonged droughts. Approximately 11.5% and 10.7% of croplands were, respectively, associated with very dry moisture condition in the 2016 and 2017 primera seasons.