Monitoring Land Surface Temperature Change with Landsat Images during Dry Seasons in Bac Binh, Vietnam

Global warming-induced climate change evolved to be one of the most important research topics in Earth System Sciences, where remote sensing-based methods have shown great potential for detecting spatial temperature changes. This study utilized a time series of Landsat images to investigate the Land Surface Temperature (LST) of dry seasons between 1989 and 2019 in the Bac Binh district, Binh Thuan province, Vietnam. Our study aims to monitor LST change, and its relationship to land-cover change during the last 30 years. The results for the study area show that the share of Green Vegetation coverage has decreased rapidly for the dry season in recent years. The area covered by vegetation shrank between 1989 and 2019 by 29.44%. Our findings show that the LST increase and decrease trend is clearly related to the change of the main land-cover classes, namely Bare Land and Green Vegetation. For the same period, we find an average increase of absolute mean LST of 0.03 ◦C per year for over thirty years across all land-cover classes. For the dry season in 2005, the LST was extraordinarily high and the area with a LST exceeding 40 ◦C covered 64.10% of the total area. We expect that methodological approach and the findings can be applied to study change in LST, land-cover, and can contribute to climate change monitoring and forecasting of impacts in comparable regions.


Introduction
According to the Intergovernmental Panel on Climate Change (IPCC) 2018 "Special Report Global Warming of 1.5 • C", climate change is occurring at a faster rate in recent years [1]. The Earth is warming and climate change causes negative impacts across the globe, including Vietnam. The Bac Binh district is a region in Vietnam, which is undergoing urbanization processes, accompanied by increasing tourism and socio-economic development. These trends come along with land-cover changes and land-use conflicts, resulting in a reduction of natural vegetation and the increase land used for infrastructure, settlement and industry, as well as agriculture and aquaculture. Climate change, along with urbanization, causes severe impacts in Vietnam, like the intensification and durations of floods and droughts. These heat-related changes put ecosystems under stress and affect the wellbeing of the local population with extreme heat events [2][3][4][5][6].
Land Surface Temperature (LST) is one of the most important climate variables and is crucial for the determination of the radiative energy budget of the Earth's surface [7]. It can describe processes such as the heat flux between land surface and the atmosphere and is therefore a valuable indicator for climate change [8]. LST also complements near surface air-temperature data and spatial in nature, thus contributing to the Sustainable Development Goal 13, Climate Action [9]. Many studies have proven the suitability of medium-to high-resolution satellite imagery to retrieve area-wide LST [10][11][12][13][14]. Most of these studies use Landsat images for LST time series analysis to derive trends and understand climate change and urbanization effects [15][16][17]. Using only LST values from a single Landsat scene to represent one season or one whole year shows some limitation of the LST value representativeness due the local daily temperature fluctuation. To overcome this, we used all the available dry season Landsat images for each year to ensure that the derived average LST better represents the year in the trend of the 30-year study period.
While many studies have been conducted on various LST trends, there are currently few LST studies done in Vietnam. Thanh Hoan et al. [18] compared land-use types and the LST derived from a single image acquisition to understand the urban heat island effect of Hanoi city. Another work has been undertaken for the city of Can Tho in the Mekong Delta [19] using Landsat image to monitor the temperature difference between urban and suburban areas from 1996 to 2016. Several other studies have used MODIS satellite imagery for LST value extraction [20][21][22][23]. MODIS product come with a large scene coverage and a relative low spatial resolution and is therefore suitable for large study areas. For the proposed study with a focus on the Bac Binh district, we utilize a Landsat time series dataset with higher spatial resolution compared to MODIS in order to achieve a higher level of detail for the LST information. We monitor the extent of four land-cover types (Bare Land, Green Vegetation, Water Bodies, Built-up Area), and the respective ranges of LST during the last 30 years. Using the Landsat archive, we can extend our study period to 30 years, back to 1989, which is important for monitoring temperature-related phenomena.
There is only one meteorological station in the surroundings of the study area that provides air temperature measurements, therefore satellite images are a useful source of data for studying the LST patterns, land-cover change and effects of climate change in the area. The proposed approach exemplifies and discusses how to use Landsat image for LST monitoring in an area with a limited number of meteorological stations. The Landsat-derived LST is verified by MODIS-derived LST as well as the in-situ data from the nearby meteorological station.
The monitoring of LST aids the understanding of the local climate settings and trends, which are important for the development of mitigation and adaptation measures for sustainability in the Bac Binh district. The focus of this study is monitoring change in four land-cover types (Bare Land, Green Vegetation, Water Bodies, Built-up Area), LST change, and the interrelation between LST and land-cover during thirty years in Bac Binh. This approach can be used for climate change research, natural resource management, and it can be applied to comparable regions.

Study Area
The Bac Binh district, Binh Thuan province is located in the south central coastal region of Vietnam ( Figure 1), with an area of 1825 km 2 and a population of more than 129,374 in 2019 [24]. Cho Lau and Luong Son are two small towns of this rural district. Its topography is hilly with an average altitude of 200 m and slopes below 8 • to the west and plain sandy coastal area to the east. Bac Binh's climate is characterized by two seasons: the dry season (from November to April) and the rainy season (from May to October). The average total precipitation in the dry season in the south-central coastal region of Vietnam is very low (below 50 mm for 6 months of dry season). January and February are the driest months that may have 0 mm rainfall [25]. In recent years, the study area has undergone rapid socio-economic development [26,27]. Population growth, urban expansion, and the changes in land-cover and the increases in temperature have been creating negative effects on the environment and affecting people's health and quality of life.

Data
This study uses Landsat imagery acquired from the United States Geological Survey (USGS) between 1989 and 2019 to derive the LST and land-cover for the Bac Binh area [28]. The dataset is comprised of three generations of Landsat sensors: Landsat TM, Landsat ETM+, and Landsat OLI TIRS. The Landsat Level 1T products are ortho-corrected [29][30][31] and have a spatial resolution of 30 m for the multispectral bands. Thermal infrared bands used in the study are Band 6 for Landsat TM and Landsat ETM+ (spatial resolution of 120 m), and Band 10 for Landsat OLI TIRS (spatial resolution of 100 m). Multiple Landsat images are used for the derivation of land-cover and LST during dry season in the study area in each year of observation ( Table 1). The selected images are cloud-free scenes, with the scanning time between around 10:30 a.m. to 11:30 a.m. local time. In addition, MODIS-LST (Day-time MOD11A2) products from 2001 to 2019 were acquired and used for comparison and evaluation of the Landsat-derived LST results [32][33][34]. The MOD11A2 LST product has a temporal resolution of eight days and a spatial resolution of 1 km.

Methodology
Prior to the extraction of LST from the thermal infrared spectral band of the Landsat imagery, we applied the Fast Line-of-sight Atmospheric Analysis of Spectral Hypercubes (FLAASH) algorithm for the atmospheric correction of all Landsat images [35][36][37]. In this study, the LST calculation using Landsat images is based on the radiative transfer equation method, with the emissivity value derived through the NDVI. The process of extracting LST is illustrated in Figure 2. LST is calculated following Equation (1) [38][39][40][41]: where LST is the land surface temperature in K. The LST value in K is converted in degrees Celsius (by subtracting 273.15) for the results and discussion sections of this study; T B is the brightness temperature in K; λ is wavelength of the emitted radiance (For Landsat TM, λ = 11.457 µm. For Landsat ETM+, λ = 11.269 µm. For Landsat OLI TIRS, λ = 10.904 µm); ρ = (h*c)/σ = 1.438 * 10 −2 mK; σ is the Boltzmann constant (1.38 × 10 −23 J/K); h is the Planck's constant (6.26 × 10 −34 JS); c is the velocity of light (2.998 × 10 8 m/s); ε is the land surface emissivity (LSE). To calculate the LST from Equation (1), it is necessary to derive T B and ε.

Extract the Brightness Temperature Value (T B )
The brightness temperature, T B , can be calculating using Equations (2) and (3): where L λ is the spectral radiance and K 1 and K 2 are constants. For Landsat TM, K 1 = 607. 76 To extract T B value, it is necessary to convert the digital number value (DN) of the image to L λ . To process multi-temporal Landsat image for the study area, the conversion for Landsat TM, Landsat ETM+, Landsat OLI TIRS images is carried out using Equations (3)-(5), respectively [42][43][44]: For Landsat TM the value is calculated using band 6 applying Equation [45]: where G rescale = 0.0551584 (Wm 2 sr 1 µm 1 )/DN, and B rescale = 1.2378 (Wm 2 sr 1 µm 1 )/DN. For Landsat ETM+ the value is calculated using band 6 applying Equation [46]: where Q cal is the quantized calibrated pixel value in (DN); Q calmin is the minimum quantized calibrate pixel value corresponding to L λmin and Q calmax is the maximum quantized calibrate pixel value corresponding to L λmax ; L λmin is the minimum spectral radiance scale to Q calmin value in (Wm 2 sr 1 µm 1 ), and L λmax is the maximum spectral radiance scale to Q calmax value in (Wm 2 sr 1 µm 1 ). For Landsat OLI TIRS the value is calculated using band 10 applying Equation [47]: where M L and A L are provided conversion factors, and Q cal is the quantized calibrated pixel value.

Extract of LSE Value
The LSE can be extracted from remote sensing image based on land-cover or vegetation index. The use of NDVI is more advantageous as it can calculate the emission at the pixel level. The emission and spectral reflectance can be measured to find the empirical relationship between LSE and NDVI using Equation (6) [48][49][50]: However, Equation (6) is only effective for a homogeneous land-cover area, for example, an area of healthy vegetation only or bare land-cover only. In this study, project area is classified into four land-cover types: (1) Water Bodies, (2) Bare Land, (3) Green Vegetation, and (4) Built-up Area. We used Equation (7) to calculate the respective LSE [51]: where ε v is the LSE for the land-cover class Green Vegetation area, ε s is the LSE for Bare Land, and p v is the share of Green Vegetation per pixel with values from 0 to 1. p v is calculated by relating the NDVI values of Bare Land and Green Vegetation with the ratio shown in Equation (8) [34,[52][53][54]: where NDVI b and NDVI v are NDVI for Bare Land and Green Vegetation, respectively. The NDVI is derived applying Equation (9): where NIR and RED are pixel values in near infrared and red bands on the image, respectively. For Landsat TM: band 3 is red, and band 4 is near infrared. For Landsat ETM+: band 3 is red, and band 4 is near infrared. For Landsat OLI TIRS: band 4 is red, and band 5 is near infrared. For NDVI values < NDVI b , p v is set to 0. For NDVI values > NDVI v p v is set to 1 [13,55].

Extract of LST from MODIS Images and Comparison with Landsat-Derived LST
Google Earth Engine (GEE) was used to access LST from the MOD11A2 products during the dry seasons between the years 2001 and 2019. For the conversion from a DN value to degrees Celsius value, the following equation was applied: where LST is the Land Surface Temperature in degrees Celsius and 0.02 is the scale factor of the MODIS LST product [56]. We calculated the mean LST value of all pixels in the study area of each selected MODIS scene. Again, the mean LST values of all the scenes of each dry season is combined into dry season mean LST value that represents each year in the chain of thirty years in our study. This MODIS-derived LST results are used to evaluate the Landsat LST results.

Land-Cover Change
Four land-cover types were derived for the Bac Binh study area, (1) Bare Land (including dry vegetation on cropland, shrubland or grassland), (2) Green Vegetation (photosynthetic active vegetation), (3) Water Bodies and (4) Built-up Area. In this study, we applied a supervised classification with the maximum likelihood algorithm to classify Landsat image into the four land-cover types for each year from 1989 to 2019 [57][58][59]. The training areas were selected by visual interpretation of the land-cover classes, selected individually on each scene. For the land-cover classification, accuracy assessment was done by a visual check at three time periods using high resolution SPOT-5 (for 2000, 2010) and SPOT-6 (for 2019) satellite images. The classification accuracy at the randomly distributed check points for the three time-spots is shown in Table 2 with the user's accuracy (UA), producer's accuracy (PA), and the overall accuracy (OA), kappa coefficient [60]. The same classification approach has been applied to all Landsat images listed in Table 1. Therefore, we assume that the accuracies of Table 2 are representative for land-cover classification from 1989 to 2019.
The land-cover classes Green Vegetation and Bare Land dominate the study area with 47% and 52% in 1989 respectively. In Figure 3a (Figures 3a and 4). The average Green Vegetation area of the first decade (1989 to 1998) covered 45% of the total project area and declined to cover only 28% on average in the period of 2010 to 2019. The class Bare Land accounted for 55% coverage in the first ten years (1989 to 1998) of the total area, whereas the in the last ten years (2010 to 2019) the class represents the dominating land-cover with a share of 71.38% of the total study area.  The areas categorized as Water Bodies and Built-up Area cover less than 1% of the project area in total, however these classes have increased 4 to 5 times over the thirty years of the observed period.
This change in land-cover has occurred at a faster rate over the last ten years than it has previously and suggests that the natural forest area and other vegetated areas have been significantly reduced. At the same time, along with processes associated with socio-economic developments, the Water Bodies and Built-up Areas has increased.

LST Results and Trends
The LST was derived from the Landsat images for the entire thirty-year period from 1989 to 2019. Figure 5 illustrates the spatial distribution of the LST mean of all scenes of each dry season for selected years.
where ∆ is an absolute average increase, decrease the value of the mean LST, Y i is mean LST value, and n is the total number of observation periods equivalent to 30. We determine ∆ according to Equation (11) for the Landsat image scenes over thirty years from 1989 to 2019. The calculation of ∆ corresponds to the three LST trends: (a) the absolute average of the minimum LST values increased by 0.1 • C per year; (b) the absolute average maximum LST values increased by 0.03 • C per year; and (c) the absolute average LST increased by 0.03 • C per year over the thirty years. The result shows that minimum LST tends to have higher increase than the maximum LST and average LST. This can be explained by the close correlation of LST with the decreasing area of Green Vegetation class in the study area. During the past thirty years, the Green Vegetation area has been reduced significantly, from 1989 accounting for nearly half of the project area, to it only accounting for 18.37% in 2019. In addition, Green Vegetation is also the class with the lowest LST among four land-cover types calculated in the area.

Comparison of Landsat LST, MODIS LST and Air Temperature (AT) In-Situ Measurements
To verify the Landsat-derived LST values, MODIS LST products (MOD11A2) were derived and averaged for each the dry season of each observation [32,61,62]. The MOD11A2 product is available from 2001 to 2019. The average Landsat LST data and the respective average MODIS products are shown in Table 3:  Table 4 shows the range of LST values between minimum and maximum of the MODIS images (from 32.22 • C to 38.75 • C) is smaller than the Landsat images (from 31.53 • C to 39.63 • C). The standard deviation value of MODIS (1.85) is also smaller than the Landsat (2.29). In this study, we used MOD11A2 LST to compare with the Landsat LST images, and RMSE (root mean square error) is calculated by the following equation: where LS is mean LST value of the Landsat image, MOD is mean LST value of the MOD11A2 image, i is the number of the years, and n is the total number of periods. In addition, we calculated the relationship of LST between Landsat and MOD11A2 images by the r coefficient (Pearson correlation coefficient) to evaluate the reliability of the results. Equation (13) for calculating the coefficient r: where LS is mean LST value of the Landsat image, MOD is mean LST value of the MOD11A2 image, LS and MOD are average of Landsat and MOD11A2 mean values, i is the number of the years, and n is the total number of period. The RMSE (root mean square error) value obtained was 0.79 • C and the r (Pearson correlation coefficient) value was 0.94. With the relatively small RMSE value of 0.79 • C, we found that the LST value obtained from the Landsat image is highly reliable for the project area. In addition, the correlation coefficient r is close to 1, which also shows a very high positive correlation between the Landsat LST values and the MODIS LST values. According to the results, we are convinced that the Landsat images can be used to retrieve reliable LST for the project area.
To verify LST data derived from Landsat image, historic AT data from the nearby weather station (Phan Thiet weather station at 108 • 06 E 10 • 56 N, about 45 km to the south from the study area) were used. They were collected daily throughout the last thirty years of study period [63]. This dataset is accessible from the website of National Climatic Data Center, U.S. Department of Commerce (https://www7.ncdc.noaa.gov/CDO/cdo). Figure 9 shows, side by side, the dry season average AT recoded at the weather station and the LST from Landsat for Bac Binh district. The difference between these two sets of temperature mean value are +7.02 • C. The average absolute temperature collected at the meteorological station (AT) in the dry season also increased by 0.03 • C per year in the period 1989-2019 at the same increase rate as LST derived from the Landsat image series.

Correlation between LST and Land-Cover
The land-cover classification results from Landsat image over time show the predominance of the two main land-cover types: Green Vegetation and Bare Land in the project area. The total coverage of these two classes is more than 99% for all the thirty years. The land-cover change in the project area is mainly between these two classes Green Vegetation and Bare Land, from Green Vegetation to Bare Land. The Water Bodies and Built-up Area cover always only less than 1% and the change of these land-cover classes has limited influence in the area.
In Figure 10a, we see the steady increase of the share of Bare Land and Built-up classes and the area with LST values >35 • C. Figure 10b shows a consistent coverage decrease of Water Bodies and Green Vegetation areas and a simultaneous decrease of the area with LST values <30 • C during the observed thirty years. Applying Equation (13) to investigate the correlation coefficient between share of the study area with LST < 30 • C and the share of the area of Green Vegetation and Water Bodies, the value r = 0.50 was obtained. For the area with LST > 35 • C and the area of Bare Land and Built-up Areas classes, r = 0.94 was obtained. Both suggests a high positive correlation between the sets of data. This result shows that the area with LST > 35 • C vs. Bare Land and Built-up Areas classes has a higher correlation than the area with LST < 30 • C vs. Green Vegetation and Water Bodies cover classes.

Discussion
The land-cover classification accuracy has been validated for the years 2000, 2010 and 2019 using high-resolution SPOT images. The overall accuracy is in the range from 70% to 80% with a Kappa coefficient of 0.63 and 0.73. These values are in range with other studies [64,65], indicating an acceptable quality for the further analysis. The classification results show that classes Bare Land and Green Vegetation dominate the study area and that the land-cover has changed greatly over thirty years. The share of Bare Land area has increased from 52.0% to 80.9% (from 949 km 2 to 1,475 km 2 ), while the share of Green Vegetation cover decreased from 47.8% to 18.37% (from 872 to 335 km 2 ) in the 1989-2019 study period. Decreasing precipitation, groundwater depletion and soil degradation are identified as drivers of this land-cover change, leading to the decrease of forested areas and even to desertification processes of the Bare Land covered areas [66,67].
Another type of land use change in the area has also been addressed in local government reports, according to which, the expansion and development of urban residential areas, tourist areas and mineral resource exploitation are continuously intensified [68]. The research result shows clear evidence of these changes. The Built-up Area has increased by more than 500% during thirty years of observation ( Figure 3b). Along with that, we also found that the area covered by Water Bodies increased four-fold as result of expanding reservoirs for storing water for daily life and for irrigation. An example is the Ca Giay Lake, one of the biggest reservoirs in the area. The dam for the reservoir was constructed from 1996 to 1999 [69] that finally contributed to the increase of Water Bodies coverage from 0.07% to 0.36% (increased by 5.3 km 2 ) 1989-2019 period. The increase of Water Bodies area from 1999 marks the first operation of the dam after the construction period (Figure 3b).
For the monitoring of LST change, this study used Landsat images, which have a higher spatial resolution than MODIS images, which was used in similar previous research on regions in Vietnam [70]. To ensure the LST value representativeness, our study used multiple Landsat scenes for each dry season (2 to 4 scenes) to obtain LST values instead of using only one scene for the whole season as it has been done in other studies [71]. The average LST of 37.07 • C for the thirty years period for the Bare Land area is always higher than the average LST of the other land-cover classes. Moreover, the average LST for Green Vegetation area is 28.70 • C and always lower than the LST of other land-cover classes. The difference between the average LST values for the two largest land-cover classes Bare Land and Green Vegetation is relative high (8.37 • C). This result is consistent with the finding from previous research on the Bare Land and Green Vegetation LST [72]. The general increase in LST over the last 30 years are in line with the 2018 IPCC report. [1]. During these thirty years, the average dry season LST shows an annual increase of 0.03 • C. The area experiencing an LST < 25 • C decreased on average by 0.41% per year and the area affected by a LST > 40 • C increased in average by 0.25% per year. There has been significant human impact on the natural land-cover area in this area, with groundwater resources and forest areas decreasing as a result to human exploitation [73][74][75][76]. We assume that the increase in LST during the dry season in Bac Binh are closely linked to the observed land-cover changes, especially the conversion from Green Vegetation to Bare Land.
To evaluate the reliability of the study result, we compared the LST value derived from Landsat image with LST extracted from MODIS MOD11A2 products (Tables 3 and 4) and AT recorded from the nearby weather monitoring station (Figure 9). The MODIS LST for each year is an average value from all images of dry season. The Landsat LST and MODIS LST comparison shows an average difference of −0.04 • C. The difference between The Landsat LST and MODIS LST in our research are consistent with findings in other studies [77,78]. The comparison shows a solid correlation between the Landsat LST and MODIS LST value series with RMSE = 0.79 • C and high positive correlation r = 0.94 that indicates the reliability of the Landsat-derived LST.
When comparing the LST results derived from Landsat images and recorded AT value from the nearby weather monitoring station, we found the same temperature increase trend of 0.03 • C per year in both datasets. However, there is constant absolute temperature difference where AT value is about 7 • C lower than LST. In other published research results [79,80], it has been observed that at the same location in clear cloudless weather, the LST measured on the ground is usually higher than AT measured in 2 m height. In our research, the nearby weather monitoring station is about 45 km outside the study area and located near the coastline where the temperature in this region is expected to be lower than in our study area. The comparison of LST and AT shows a consistent temperature increase for both sets of measurements. Coincident with the results from other studies [79,80], the LST and AT values in our study area are highly correlated during the thirty years from 1989 to 2019.
The AT measurements of the Phan Thiet weather station, as well as regional climate databases (http://berkeleyearth.lbl.gov/regions/vietnam) confirm the low LST value for 2011, as well as high values for example for the years 1998 and 2015. Yet, the extremely high average LST value with 39.63 • C for 2005, which was confirmed by the MODIS MOD11A2 product for the respective year, can neither be explained by the AT measurement of the 2005 dry-period nor by regional climate statistics.
The thirty-year LST dataset allowed detailed monitoring and analyses the relation between land-cover change and the LST fluctuation in the project area. In contrast to other studies in Vietnam and worldwide [20,81,82] that investigated only at the difference in LST between different regions or different land-cover types, our research also compares the increase and decrease of trends of LST levels with the trends of land-cover change. The results show a close relationship between the change in coverage of Green Vegetation and Water Bodies versus the area with LST < 30 • C, as well as between Built-up Area and Bare Land area versus the area with LST > 35 • C. The total area of Bare Land, Built-up Area and the area with LST > 35 • C have steadily increased, while the total area of Green Vegetation, Water Bodies area and the area with LST < 30 • C are decreasing. Analyses came up with correlation value between area with LST > 35 • C and Bare Land and Built-up Areas classes, r = 0.94, that shows close correlation between the increase of the area with LST > 35 • C and the increase of Bare Land and Built-up Area. This leads to the conclusion that the land-cover change in the project area contributes significantly to the LST increase.

Conclusions
The results proved the suitability of monitoring LST using Landsat multi-temporal satellite images. The study shows that, during thirty years of social and economic development, the land-cover in the study area has been changed significantly in the Bac Binh district, Vietnam. The share of Bare Land has increased from 52.0% to 80.9% (from 949 km 2 to 1,475 km 2 ), while the share of Green Vegetation cover decreased from 47.8% to 18.37% (from 872 to 335 km 2 ) in the 1989-2019 study period. The area covered by Water Bodies and Built-up Area has increased more than four-and five-fold respectively.
We have shown that areas affected by high LST (>35 • C) have increased over the past thirty years, during which, the absolute average LST has increased on average by 0.03 • C per year. In addition, we used LST value obtained from MODIS and AT from a weather monitoring station for comparison and evaluation of the Landsat-derived LST. The results confirm the Landsat derived LST trends. MODIS LST and AT datasets also showed a 0.03 • C increase in temperature per year as the LST derived from the Landsat time-series.
The research result shows that the area with LST < 30 • C and the area of Green Vegetation and Water Bodies land-cover classes are similarly decreasing. At the same time, the opposite trend is seen with the relationship between the area with LST > 35 • C and the area of Bare Land and Built-up Area. This indicates a clear correlation between the LST trends and the change in land-cover types of the region over the 30-year study period. This is the first study on the change of LST for a thirty-year period in the Bac Binh area. The high spatial resolution and the results from the LST time series provide a detailed and reliable dataset, compared to other studies that investigate correlating trends of LST and land-cover. The study uses multiple image scenes per dry season to ensure LST representativeness and validates the results successfully with MODIS LST products and AT measurements of a nearby weather station.
In order to achieve a more comprehensive assessment and analysis on the development of LST in this study area or comparable sites, it will be necessary to perform the analysis on the entire year, including the rainy season. Yet, the results for the dry-season provide valuable information, which might contribute to a better understanding of the climate change effects and to develop adaptation strategies for water management and spatial planners to better cope with climate change and warming of the Bac Binh study area.