Impact of large-scale vegetation restoration project on summer land surface temperature on the Loess Plateau, China

A large-scale afforestation project has been carried out since 1999 in the Loess Plateau of China. However, vegetation-induced changes in land surface temperature (LST) through the changing land surface energy balance have not been well documented. Using satellite measurements, this study quantified the contribution of vegetation restoration to the changes in summer LST and analyzed the effects of different vegetation restoration patterns on LST during both daytime and nighttime. The results show that the average daytime LST decreased by 4.3°C in the vegetation restoration area while the average nighttime LST increased by 1.4°C. The contributions of the vegetation restoration project to the changes in daytime LST and nighttime LST are 58% and 60%, respectively, which are far greater than the impact of climate change. The vegetation restoration pattern of cropland (CR) converting into artificial forest (AF) has a cooling effect during daytime and a warming effect at nighttime, while the conversion of CR to grassland has an opposite effect compared with the conversion of CR to AF. Our results indicate that increasing evapotranspiration caused by the vegetation restoration on the Loess Plateau is the controlling factor of daytime LST change, while the nighttime LST change is affected by soil humidity and air humidity.

Spatio-Temporal Variability in Land Surface Temperature and Its Relationship with Vegetation Types over Ibadan, South-Western Nigeria

Ibadan has experienced a rapid urbanization over the past few decades due to heavy influx of people from different parts of the country as a result of improved economy of the region. This development induced a great change in land use and land cover over the region which has become a major environmental concern recently. This study assessed Land Surface Temperature (LST) and its spatio-temporal relationship with land cover type over Ibadan. Land use/Land cover dynamics were assessed using index maps generated from Landsat Satellite data (TM, ETM+ and OLI) of Ibadan. The corrected thermal Infrared bands of the Landsat data were used to retrieve LST. The results revealed a notable increase in built-up areas from 5.64% of the total land cover area in 1984 to 14.05% in 2014. This change has caused increase in surface temperature of Ibadan from 3.56?C to 8.54?C between 1984 and 2014 respectively. The study recorded a continuous decrease in the vegetal part of Ibadan (from 43.28% in 1984 to 14.76 in 2014) which could be attributed to anthropogenic activities as the vegetated land area lost was been converted to other form of land use. The change was found to be positively correlated to the surface temperature intensity over the region with correlation coefficient, r value of 0.9251, 0.8256 and 0.7017 in 1984, 2000 and 2014 respectively. It is recommended that Policies should be considered for planting trees, new guidelines for urban landscape design and construction.

Vegetation Cover Density and Land Surface Temperature Interrelationship Using Satellite Data, Case Study of Wadi Bisha, South KSA

The aim of this study is to identify the relationship between Vegetation Cover (VC) and the land Surface Temperature (LST), using satellite data of Wadi Bisha, south the Kingdome of Saudi Arabia (KSA). The Landsat 7 Thematic Mapper (ETM) thermal band (band 6) was used for calculating the (LST) values. The near-infrared (NIR) and red band (bands 3 and 4 respectively) were used for estimating the vegetation cover. ERDAS Imagine 9.3 and ArcGIS 10.2 were used in the current study. The results of the study show that the increase of vegetation cover (VC) coincides with decrease of (LST), while the decrease in vegetation cover is linked with increase of (LST). It was found that there was no vegetation observed in areas practiced the highest temperature of 49℃, while areas of lowest temperature of 28℃ were characterized by dense vegetation cover. Thus, a quite significant correlation is approved between the (VC) and the (LST), based on the validation of (50) locations. It was concluded that availability and continuity of Satellite remote sensing data was required for elaborating a continuous monitoring of vegetation cover conditions and mapping was recommended in Wadi Bisha. Operational monitoring is recommended to ensure the adoption of flexible land cover validation protocols.

Spatial Variation in CO_(2) Concentration Improves the Simulated Surface Air Temperature Increase in the Northern Hemisphere

The increasing concentration of atmospheric CO_(2) since the Industrial Revolution has affected surface air temperature.However,the impact of the spatial distribution of atmospheric CO_(2) concentration on surface air temperature biases remains highly unclear.By incorporating the spatial distribution of satellite-derived atmospheric CO_(2) concentration in the Beijing Normal University Earth System Model,this study investigated the increase in surface air temperature since the Industrial Revolution in the Northern Hemisphere(NH) under historical conditions from 1976-2005.In comparison with the increase in surface temperature simulated using a uniform distribution of CO_(2),simulation with a nonuniform distribution of CO_(2)produced better agreement with the Climatic Research Unit(CRU) data in the NH under the historical condition relative to the baseline over the period 1901-30.Hemispheric June-July-August(JJA) surface air temperature increased by 1.28℃ ±0.29℃ in simulations with a uniform distribution of CO_(2),by 1.00℃±0.24℃ in simulations with a non-uniform distribution of CO_(2),and by 0.24℃ in the CRU data.The decrease in downward shortwave radiation in the non-uniform CO_(2) simulation was primarily attributable to reduced warming in Eurasia,combined with feedbacks resulting from increased leaf area index(LAI) and latent heat fluxes.These effects were more pronounced in the non-uniform CO_(2)simulation compared to the uniform CO_(2) simulation.Results indicate that consideration of the spatial distribution of CO_(2)concentration can reduce the overestimated increase in surface air temperature simulated by Earth system models.

Dynamic Drought Monitoring in Guangxi Using Revised Temperature Vegetation Dryness Index

Moderate resolution imaging spectroradiometer （MODIS） data are very suitable for vast extent, long term and dynamic drought monitoring for its high temporal resolution, high spectral resolution and moderate spatial resolution. The composite Enhanced Vegetation Index （EVI） and composite land surface temperature （Ts） obtained from MODIS data MOD11A2 and MOD13A2 were used to construct the EVI-Ts space. And Temperature Vegetation Dryness Index （TVDI） was calculated to evaluate the agriculture drought in Guangxi province, China in October of 2006. The results showed that the drought area in Guangxi was evidently increasing and continuously deteriorating from the middle of September to the middle of November. The TVDI, coming from the EVI-Ts space, could effectively indicate the spatial distribution and temporal evolution of drought, so that it could provide a strong technical support for the forecasting agricultural drought in south China.

Examining the Relationship Between Spatial Configurations of Urban Impervious Surfaces and Land Surface Temperature

The urban heat island(UHI) effect has significant effects on the quality of life and public health. Numerous studies have addressed the relationship between UHI and the increase in urban impervious surface area(ISA), but few of them have considered the impact of the spatial configuration of ISA on UHI. Land surface temperature(LST) may be affected not only by urban land cover, but also by neighboring land cover. The aim of this research was to investigate the effects of the abundance and spatial association of ISAs on LST. Taking Harbin City, China as an example, the impact of ISA spatial association on LST measurements was examined. The abundance of ISAs and the LST measurements were derived from Landsat Thematic Mapper(TM) imagery of 2000 and 2010, and the spatial association patterns of ISAs were calculated using the local Moran’s I index. The impacts of ISA abundance and spatial association on LST were examined using correlation analysis. The results suggested that LST has significant positive associations with both ISA abundance and the Moran’s I index of ISAs, indicating that both the abundance and spatial clustering of ISAs contribute to elevated values of LST. It was also found that LST is positively associated with clustering of high-ISA-percentage areas(i.e.,>50%) and negatively associated with clustering of low-ISA-percentage areas(i.e.,<25%). The results suggest that, in addition to the abundance of ISAs,their spatial association has a significant effect on UHIs.

Impact of land cover change on land surface temperature: A case study of Spiti Valley

Land surface temperature(LST) is the skin temperature of the earth surface. LST depends on the amount of sunlight received by any geographical area. Apart from sun light, LST is also affected by the land cover, which leads to change in land surface temperature. Impact of land cover change(LCC) on LST has been assessed using Landsat TM5, Landsat 8 TIRS/OLI and Digital Elevation Model(ASTER) for Spiti Valley, Himachal Pradesh, India. In the present study, Spiti valley was divided into three altitudinal zones to check the pattern of changing land cover along different altitudes and LST was calculated for all the four land cover categories extracted from remote sensing data for the years of 1990 and 2015. Matrix table was used as a technique to evaluate the land cover change between two different years. Matrix table shows that as a whole, about 2,151,647 ha(30%) area of Spiti valley experienced change in land cover in the last 25 years. The result also shows vegetation and water bodies increased by 107,560.2 ha(605.87%) and 45 ha(0.98%), respectively. Snow cover and barren land decreased by 19,016.5 ha(23.92%) and 88,589(14.14%), during the study period. A significant increase has been noticed in vegetation amongst all land cover types. Minimum, maximum and mean LST for three altitudinal zones have been calculated. The mean LST recorded was 11℃ in 1990 but it rose by 2℃ and reached to 13℃ in 2015. Changes in LST were obtained for each land cover categories. The mean temperature of different land cover types was calculated by averaging value of all pixels of a given land cover types. The mean LST of vegetation, barren land, snow cover and water body increased by 6℃, 9℃, 1℃, and 7℃, respectively. Further, relationships between LST, Normalized Difference Snow Index(NDSI), and Normalised Difference Vegetation Index(NDVI) were established using Linear Regression.

Estimation of Land Surface Temperature from Landsat-8 OLI Thermal Infrared Satellite Data. A Comparative Analysis of Two Cities in Ghana

This study employs Landsat-8 Operational Land Imager (OLI) thermal infrared satellite data to compare land surface temperature of two cities in Ghana: Accra and Kumasi. These cities have human populations above 2 million and the corresponding anthropogenic impact on their environments significantly. Images were acquired with minimum cloud cover (<10%) from both dry and rainy seasons between December to August. Image preprocessing and rectification using ArcGIS 10.8 software were used. The shapefiles of Accra and Kumasi were used to extract from the full scenes to subset the study area. Thermal band data numbers were converted to Top of Atmospheric Spectral Radiance using radiance rescaling factors. To determine the density of green on a patch of land, normalized difference vegetation index (NDVI) was calculated by using red and near-infrared bands i.e. Band 4 and Band 5. Land surface emissivity (LSE) was also calculated to determine the efficiency of transmitting thermal energy across the surface into the atmosphere. Results of the study show variation of temperatures between different locations in two urban areas. The study found Accra to have experienced higher and lower dry season and wet season temperatures, respectively. The temperature ranges corresponding to the dry and wet seasons were found to be 21.0985oC to 46.1314oC, and, 18.3437oC to 30.9693oC respectively. Results of Kumasi also show a higher range of temperatures from 32.6986oC to 19.1077oC during the dry season. In the wet season, temperatures ranged from 26.4142oC to -0.898728oC. Among the reasons for the cities of Accra and Kumasi recorded higher than corresponding rural areas’ values can be attributed to the urban heat islands’ phenomenon.

Land Characterization Analysis of Surface Temperature of Semi-Arid Mountainous City Abha, Saudi Arabia Using Remote Sensing and GIS

This knowledge of land surface temperature and its spatial variations within a city environment is of prime importance to the study of urban climate and human-environment interactions. Few studies have examined the influence of land use and terrain on the surface temperature effects of semi-arid mountainous urban areas. This study investigates the urban environment characterization and its effects on surface temperature using remote sensing. The methodologies adapted for this study are geometric and radiometric corrections of satellite data, extraction of land use/land cover and digital elevation model, estimation of vegetation density using Normalized Difference Vegetation Index (NDVI), and estimation of surface temperature and emissivity using temperature emissivity separation (TES) algorithm. Finally geospatial model and statistical techniques are used for assessing the overall impact of urban environmental characterization on urban climate of semi-arid region of Abha, Kingdom of Saudi Arabia. Herein, results reveal that the spatial distribution of surface temperature was affected by land use/land cover (LULC) and topography. The high dense built-up and commercial/industrial areas display higher surface temperature in comparison with surrounding lands. There is gradual decrease of LULC classes’ surface temperature with the increase in altitude. The cooling effect towards the surrounding urban built-up area is found increasing at the hill located vegetated area, the downward slope and valley terrain inside the recreational park. Therefore the spatial variation in surface temperature also reflected the effects of topography on LULC classes. Suitable mountainous land use utilization would help to expand the cooling effect. In the future, the outcomes of this study could be used to build environmentally sustainable urban planning suitable to semi-arid regions and to create practices that consider the local weather environment in urban planning.

An Analysis of Spatio-Temporal Trends of Land Surface Temperature in the Dhaka Metropolitan Area by Applying Landsat Images

Land surface temperature (LST) is a basic determinant of the global thermal behavior of the Earth surface. LST is a vital consideration for the appraisal of gradual thermal change for urban areas to examine the strength of the thermal intensity of the surface of urban heat island (SUHI) and to see how hot the surface of the Earth would be in a particular location. In this respect, the most developed urban city like Dhaka Metropolitan Area (DMA), Bangladesh is considered for estimation of LST, and Normalized Difference Vegetation Index (NDVI) changes trend in more developed and growing developing areas. The focus of this study is to find out the critical hotspot zones for further instantaneous analysis between these two types of areas. The trends of long-term spatial and temporal LST and NDVI are estimated applying Landsat images-Landsat 5-TM and Landsat OLI_TIRS-8 for the period of 1988 to 2018 for DMA and for developed and growing developing areas during the summer season like for the month of March. The supervised classification was used to estimate the land cover categories and to generate the LST trends maps of the different percentiles of LSTs over time using the emissivity and effective at sensor brightness temperature. The study found the change in land cover patterns by different LST groups based on 50th, 75th, and 90th percentile where the maximum LST for the whole DMA went up by 2.48°C, 1.01°C, and 3.76°C for the months of March, April, and May, respectively for the period of 1988 to 2018. The highest difference in LST was found for the most recently developed area. The moderate change of LST increased in the built-up areas where LST was found more sensitive to climate change than the growing developed areas. The vegetation coverage area decreased by 6.74% in the growing, developing areas compared to the developed areas from 1988 to 2018. The findings of the study might be helpful for urban planners and researchers to take up appropriate measures to mitigate the thermal effect on urban environment.

Effects of land use and cover change on surface wind speed in China

The surface wind speed(SWS)is affected by both large-scale circulation and land use and cover change(LUCC).In China,most studies have considered the effect of large-scale circulation rather than LUCC on SWS.In this study,we evaluated the effects of LUCC on the SWS decrease during 1979-2015 over China using the observation minus reanalysis(OMR)method.There were two key findings:(1)Observed wind speed declined significantly at a rate of 0.0112 m/(s·a),whereas ERA-Interim,which can only capture the inter-annual variation of observed data,indicated a gentle downward trend.The effects of LUCC on SWS were distinct and caused a decrease of 0.0124 m/(s·a)in SWS;(2)Due to variations in the characteristics of land use types across different regions,the influence of LUCC on SWS also varied.The observed wind speed showed a rapid decline over cultivated land in Northwest China,as well as a decrease in China’s northeastern and eastern plain regions due to the urbanization.However,in the Tibetan Plateau,the impact of LUCC on wind speed was only slight and can thus be ignored.

Estimating Evapotranspiration Using Improved Fractional Vegetation Cover and Land Surface Temperature Space

Vegetation index-land surface temperature （VI-T s ） space has been widely used to estimate evapotranspiration and soil moisture. The limitation of this method is the uncertainty of the observed dry edge, which is usually fitted by scatter plots. Here, a method was used to locate true dry and wet edges based on energy balance formulation, and a simple method to estimate surface energy flux is proposed based on the improved Fractional vegetation cover-Land surface temperature （F v -T s ） space. Seventeen days of MODIS products were selected to estimate evapotranspiration and the estimated sensible heat flux （H） is compared with Large Aperture Scintillometer （LAS） data at a site in Zhengzhou, resulting in a RMSE of 44.06 W m^-2 , bias of 36.99 W m^-2 and R^2 of 0.71. The H scatter plots of estimation versus observation show clearly that most points are around the 1：1 line. Overall, the located true and wet edges are more accurate than the observed true edge. Our results can also be applied to improve the estimation of soil moisture.

Analysis of Spatio-temporal Land Surface Temperature and Normalized Difference Vegetation Index Changes in the Andassa Watershed, Blue Nile Basin, Ethiopia

Analysis of the nexus between vegetation dynamics and climatic parameters like surface temperature is essential in environmental and ecological studies and for monitoring of the natural resources. This study explored the spatio-temporal distribution of land surface temperature(LST) and Normalized Difference Vegetation Index(NDVI) and the relationship between them in the Andassa watershed from 1986 to 2016 periods using Landsat data. Monthly average air temperature data of three meteorological sites were used for validating the results. The findings of the study showed that the LST of the Andassa watershed has increased during the study periods. Overall, average LST has been rising with an increasing rate of 0.081?C per year. Other results of this study also showed that there has been a dynamic change in vegetation cover of the watershed in all seasons. There was also a negative correlation between LST and NDVI in all the studied years. From this study we can understand that there has been degradation of vegetation and intensification of LST from 1986 to 2016.

Applying the Global Disturbance Index for Detecting Vegetation Changes in Lao Tropical Forests

Land cover change is a major challenge for many developing countries. Spatiotemporal information on this change is essential for monitoring global terrestrial ecosystem carbon, climate and biosphere exchange, and land use management. A combination of LST and the EVI indices in the global disturbance index (DI) has been proven to be useful for detecting and monitoring of changes in land covers at continental scales. However, this model has not been adequately applied or assessed in tropical regions. We aimed to demonstrate and evaluate the DI algorithm used to detect spatial change in land covers in Lao tropical forests. We used the land surface temperature and enhanced vegetation index of the Moderate Resolution Imaging Spectroradiometer time-series products from 2006-2012. We used two dates Google EarthTM images in 2006 and 2012 as ground truth data for accuracy assessment of the model. This research demonstrated that the DI was capable of detecting vegetation changes during seven-year periods with high overall accuracy;however, it showed low accuracy in detecting vegetation decrease.

Intra-annual Vegetation Change Characteristics in the NDVI-Ts Space:Application to Farming-Pastoral Zone in North China

Land use/land cover change (LUCC) mapping and analysis using multi-temporal normalize difference vegetation index (NDVI) data have been well documented. Recent empirical studies have documented that many new methods with high accuracy of retrieved land surface temperature ( Ts) have been developed. Thus, the combination of land surface temperature and NDVI has the greatest potential to improve the surface vegetation dynamic monitoring. In this study, the following objectives are pursued to: (1) introduce the practical method to produce the Ts, NDVI and Ts/NDVI based on remotely sensed data; (2) investigate the different retrieved result of vegetation cover information from NDVI, Ts and Ts/NDVI data sets, and analyze the intra-annual time trajectories of different vegetation cover categories in the NDVI- Ts space for farming-pastoral zone in North China, and (3) quantitative analysis the difference in using NDVI, Ts and Ts/NDVI data sets to express information based on the indices (information entropy and averaged information grads), and evaluate the relative role of Ts/NDVI data set in the discrimination of different vegetation cover categories through comparison to traditional NDVI data set.

Semi-analytical analysis of the response of the air temperature over the land surface to the global vegetation distribution

Response of the air temperature over the land surface to the global vegetation distribution is investigated, using a three-dimensional governing equation to simulate the steady, large-scale, limited amplitude perturbation of the free, inviscid and adiabatic atmosphere. The adoption of the static equation leads to a temperature governing equation in the terrain following coordinate. With the prescribed temperature as the upper boundary condition and the radiation balance as the lower boundary condition, the semi-analytical solution of the global circulation temperature can be calculated. In this article, only the air temperature (at 2 m height) over the land surface is analyzed, and the result suggests that this model can simulate the air temperature pattern over the land surface reasonably. A better simulation occurs when a simple feedback of the albedo on the temperature is included. Two sensitivity experiments are analyzed through this model. One suggests that the air temperature over the land surface descends obviously when the land surface is covered with ice all over, while another suggests that the air temperature rises a little when the land surface is covered with forest except the ice-covered area. This model appears to be a good tool to study the response of the air temperature to the vegetation distribution. Limitations of the model are also discussed.

Feedback and contribution of vegetation, air temperature and precipitation to land surface temperature in the Yangtze River Basin considering statistical analysis

Land surface temperature(LST),especially day-night LST difference(LSTd-LSTn),is a key variable for the stability of terrestrial ecosystems,affected by vegetation and climate change.Quantifying the contribution and feedback of vegetation and climate to LST changes is critical to developing mitigation strategies.Based on LST,Normalized vegetation index(NDVI),land use(LU),air temperature(AT)and precipitation(Pre)from 2003 to 2021,partial correlation was used to analyze the response of LST to vegetation and climate.The feedback and contribution of both to LST were further quantifed by using spatial linear relationships and partial derivatives analysis.The results showed that both interannual LST(LSTy)and LSTd-LSTn responded negatively to vegetation,and vegetation had a negative feedback effect in areas with significantly altered.Vegetation was also a major contributor to the decline of LSTd-LSTn.With the advantage of positive partial correlation area of 94.99%,AT became the main driving factor and contributor to LSTy change trend.Pre contributed negatively to both LSTy and LSTd-LSTn,with contributions of-0.004℃/y and-0.022℃/y,respectively.AT played a decisive role in LST warming of YRB,which was partially mitigated by vegetation and Pre.The present research contributed'to,the,detection,of LST changes and improved understanding of the driving mechanism.

Soil temperature estimation at different depths,using remotely-sensed data

Soil temperatures at different depths down the soil profile are important agro-meteorological indicators which are necessary for ecological modeling and precision agricultural activities. In this paper, using time series of soil temperature(ST) measured at different depths(0, 5, 10, 20, and 40 cm) at agro-meteorological stations in northern China as reference data, ST was estimated from land surface temperature(LST) and normalized difference vegetation index(NDVI) derived from AQUA/TERRA MODIS data, and solar declination(Ds) in univariate and multivariate linear regression models. Results showed that when daytime LST is used as predictor, the coefficient of determination(R^2) values decrease from the 0 cm layer to the 40 cm layer. Additionally, with the use of nighttime LST as predictor, the R^2 values were relatively higher at 5, 10 and 15 cm depths than those at 0, 20 and 40 cm depths. It is further observed that the multiple linear regression models for soil temperature estimation outperform the univariate linear regression models based on the root mean squared errors(RMSEs) and R^2. These results have demonstrated the potential of MODIS data in tandem with the Ds parameter for soil temperature estimation at the upper layers of the soil profile where plant roots grow in. To the best of our knowledge, this is the first attempt at the synergistic use of LST, NDVI and Ds for soil temperature estimation at different depths of the upper layers of the soil profile, representing a significant contribution to soil remote sensing.

Changes in Global Vegetation Distribution and Carbon Fluxes in Response to Global Warming:Simulated Results from IAP-DGVM in CAS-ESM2

Terrestrial ecosystems are an important part of Earth systems,and they are undergoing remarkable changes in response to global warming.This study investigates the response of the terrestrial vegetation distribution and carbon fluxes to global warming by using the new dynamic global vegetation model in the second version of the Chinese Academy of Sciences(CAS)Earth System Model(CAS-ESM2).We conducted two sets of simulations,a present-day simulation and a future simulation,which were forced by the present-day climate during 1981-2000 and the future climate during 2081-2100,respectively,as derived from RCP8.5 outputs in CMIP5.CO_(2)concentration is kept constant in all simulations to isolate CO_(2)-fertilization effects.The results show an overall increase in vegetation coverage in response to global warming,which is the net result of the greening in the mid-high latitudes and the browning in the tropics.The results also show an enhancement in carbon fluxes in response to global warming,including gross primary productivity,net primary productivity,and autotrophic respiration.We found that the changes in vegetation coverage were significantly correlated with changes in surface air temperature,reflecting the dominant role of temperature,while the changes in carbon fluxes were caused by the combined effects of leaf area index,temperature,and precipitation.This study applies the CAS-ESM2 to investigate the response of terrestrial ecosystems to climate warming.Even though the interpretation of the results is limited by isolating CO_(2)-fertilization effects,this application is still beneficial for adding to our understanding of vegetation processes and to further improve upon model parameterizations.

Analysis of Climate Variability and Relation to Vegetation in Garamba National Park from 1990-2020

Climatic factors impact vegetation. Our study was to examine and analyze the climate variability and relationship to vegetation in Garamba National Park of the Democratic Republic of the Congo over the past 30 years (1990 to 2020), then to relate the climatic variables. Mann Kendall’s non parametric test, ANOVA, and p-value tests are used to analyze existing trends and relationships between vegetation cover, climatic factors, land surface temperature (LST) and normalized difference in temperature Vegetation index (NDVI), Enhanced vegetation index (EVI) in Garamba national park which is of particular importance for the network of protected areas of the Democratic republic of Congo because its position at the northern limit of the savanna-forest mosaics gives it a unique biodiversity. The southern part of the park is dominated by grassy shrub savannas. The results showed that: 1) In Garamba, the monthly correlation coefficient of Kendall and Pearsan between temperature and precipitation are negative respectively 0.763 and −0.876 (p-value < 0.00001). 2) Annually during the three decades in Garamba, the correlation between precipitation and NDVI is significant 0.416 (Kendall) and 0.496 (Pearsan);the same between precipitation and EVI 0.291 (Kendall) and 0.496 (Pearsan) while LST and precipitation are negatively correlated (p-value < 0.00001).