A Time Series Analysis of Forest Cover and Land Surface Temperature Change Over Dudpukuria-Dhopachari Wildlife Sanctuary Using Landsat Imagery

Forest cover change is an important criterion as it affects the environmental balance whereas land surface temperature is a significant parameter within the earth climate system. Spatio-temporal change of forest cover can be detected and land surface temperature can be retrieved by applying remote sensing technology. The present study aimed to capture the impact of forest cover change on land surface temperature in Dudpukuria-Dhopachari Wildlife Sanctuary (DDWS), Bangladesh, using multi-spectral and multi-temporal satellite data. To avoid the biasness in the calculation, leaf flash time was targeted for collecting Landsat images from United States Geological Survey (USGS) Earth Explorer and, based on availability, images were collected purposively which ones had closer time period:1990 (March 5, 1990), 2000 (February 5, 2000), 2010 (February 24, 2010) and 2020 (March 23, 2020). Unsupervised classification was applied over the images Landsat 4–5 Thematic Mapper (TM), 7 Enhanced Thematic Mapper Plus (ETM+), and 8 Operational Land Imager (OLI) and Thermal Infrared Sensor (TIRS) data for detecting forest cover change. To retrieve the land surface temperature, Mono Window Algorithm (MWA) method was applied over similar images. Maximum forest degradation was observed in 2010 and the change found was 17% as compared to 1990. After 2010, the forest started to flourish. Land surface temperature dramatically changes over the time period. The highest land surface temperature in the forested area was observed in 2020 (32.2°C) and it was changed 7.7°C from that of the 1990 (24.5°C). In every 10 years, almost 2.3°C–3.0°C temperature change was detected. In the first three decades, a reverse relationship was observed between land surface temperature and forest cover; however, in the last decade, land surface temperature was found to increase with the increase of forest cover. Thus, the results of the study revealed that land surface temperature may not be relevant with the local forest cover change directly. It can be estimated from the results that local forest cover change may have limited impact on local temperature rather than global forest cover change, whereas global warming could play a vital role in changing land surface temperature locally as well as globally.

Introduction

Remote sensing is an important tool to analyse the spatial and temporal changes of land coverage (Petropoulos et al., 2015; Chew et al., 2016; Sekertekin and Bonafoni, 2020). This is a potential technology for inventory and monitoring of world’s natural resources (Goetz and Dubayah, 2011; Sharma et al., 2021). Forest is an important natural component. It maintains the environmental and ecological balance (Liu et al., 2020; Yuan et al., 2020) by preventing unwanted climate change and keeping ecosystems balance between human beings, plants, animals, and other abiotic components. Satellite imagery is a useful unbiased space borne photography as it detects the changes of land coverage in a periodic repetitive manner. Thus, forest cover changes and damages are eventually assessed by applying remote sensing technology (Zeng et al., 2020).

Remote sensing technology is an important new source of data that offers accurate geo-referencing information for harvest and analysis with comparatively low-cost but faster precise results (Lein, 2011). Among other branches, thermal remote sensing is the one that deals with Thermal Infrared (TIR) data of the Electromagnetic (EM) spectrum (Al-doski et al., 2016; Messina and Modica, 2020). The ground emits radiation and the thermal remote sensing captures it. The land surface temperature represents the temperature of the earth’s surface and is estimated from radiance measurements (Sekertekin and Bonafoni, 2020; Wan Mohd Jaafar et al., 2020). Remotely sensed thermal infrared data analyse the large scale temporal and spatial land surface temperature (LST) from local as well as global levels.

Radiative Transfer Equation (RTE) method (Saratoon et al., 2013), Single Channel Algorithm (SCA) (Dong et al., 2017), and Mono Window Algorithm (MWA) (Qin et al., 2001; Kumari et al., 2020) are the several methods for LST retrieval using Landsat 5 Thematic Mapper (TM), 7 Enhanced Thematic Mapper Plus (ETM+), and 8 Operational Land Imager (OLI) and Thermal Infrared Sensor (TIRS) data (García-Santos et al., 2018; Hao et al., 2019). In this study, Mono Window Algorithm (MWA) was evaluated for Landsat 5, 7, and 8 imageries.

Forests in developing countries are under greater risk as the per capita land is limited here (Fagan et al., 2020; Huang et al., 2020). Agricultural expansion in the forest land and illegal harvesting from the natural forest area are the major drivers of deforestation (Islam and Sato, 2012; Kissinger et al., 2012; Hishe et al., 2021). Forest degradation in an alarming rate is a global issue now (Halme et al., 2013; Morales-Hidalgo et al., 2015; Hite and Seitz, 2021; Wagner et al., 2021).

Bangladesh is a small developing country with massive population that puts the forest coverage under greater dangers (Singh et al., 2020; Van Schendel, 2020). Forest coverage here is as an open treasure to collect resources. To protect these resources, the government of Bangladesh has declared different forest areas as protected area at different times (Chowdhury and Koike, 2010; Chowdhury, 2014; Mukul et al., 2016; Rahman and Islam, 2021). In this special consideration, there are 49 types of protected areas under the broader classification of national park, wildlife sanctuary, special biodiversity conservation area, marine protected area, vulture safe zone, botanical garden, safari park, eco-park, and aviary park (Rahman et al., 2016; BFD, 2021).

Protected area is a clearly defined marine and terrestrial conservation place, recognised, dedicated, and managed, through legal or other effective means, to achieve the long-term conservation of nature with associated ecosystem services and cultural values (Dudley et al., 2010; IUCN and UNEP-WCMC, 2014; Mitchell et al., 2018). Protected areas are very important as they help to protect species from extinction, boost carbon storage, maintain regular weather condition, and provide natural habitats for plants and animals (Mackey et al., 2008; Stolton et al., 2015; Xu et al., 2017). Thus, it is important to analyze Landsat data with land cover classification methods and LST retrieval methods over protected areas.

Several studies have been done by different researchers on land cover classification (Rodriguez-Galiano et al., 2012; Cetin and Sevik, 2020; Chowdhury et al., 2020; Kaya and Görgün, 2020; MohanRajan et al., 2020; Munthali et al., 2020; Olorunfemi et al., 2020) and LST retrieval (Kayet et al., 2016; Sekertekin and Bonafoni, 2020; Thakur et al., 2021; Rashid et al., 2021). Land cover changes and resulting land surface temperature in Bangladesh have been studied sporadically by various researchers (Dewan and Corner, 2012; Ahmed et al., 2013; Chaudhuri and Mishra, 2016; Ullah et al., 2019; Gazi et al., 2020; Kafy et al., 2020). In this continuation, the GIS and RS techniques were applied in the present study to detect the forest damages and land surface temperature changes, and to invent their interrelationship over Dudpukuria-Dhopachari Wildlife Sanctuary (DDWS) since 1990 to 2020 as no research have been done on land cover classification and LST retrieval over this protected area.

Materials and Methods

Study Location

Dudpukuria-Dhopachari Wildlife Sanctuary (DDWS) is situated along the borderline of Chittagong, Rangamati, and Bandarban districts, Bangladesh. The area extends from N 22°09′ to N 22°22′ and E 92°05′ to E 92°10′ (Shahadat et al., 2014; Hossain et al., 2017). This evergreen and semi-evergreen tropical natural forest was declared as a protected area of Bangladesh in 2010 and covers 4,717 ha of land area (Shahadat et al., 2014; BFD, 2021). This vital hilly forest ecosystem also serves as a watershed of two significant rivers of Bangladesh: Karnaphuli and Sangu (Shahadat et al., 2014). Dudpukuria-Dhopachari Wildlife Sanctuary also preserves 386 tons of CO₂ per hectare and is the harbor of 600 plant species and 400 vertebrates including 200 birds and 45 mammals. This tropical hilly forest is generally dominated by two tree species Dipterocarpus turbinatus (Garjan) and Artocarpus chama (Chapalish) (Hassan et al., 2011). The forest is also the herd of globally endangered Elephas maximus (Asian Elephant), Hoolock hoolock (Western Hoolock Gibbon), and globally vulnerable Arctictis binturong (Binturong) (Karim and Ahsan, 2016). People around this wildlife sanctuary are living in 28 villages, five of which are ethnic villages. These people work with Bangladesh Forest Department (BFD) to protect the forest from destruction in a co-managed manner. They help to expand the forest cover area by doing afforestation and assisting natural regeneration (ANR). These people also keep away intruders, illegal trespassers, and hunters by regular patrolling in the forest.

Data Collection

Landsat images were collected from the website of United States Geological Survey (USGS) Earth Explorer. In the study, forest cover changes were detected from 1990 to 2020. Landsat 4–5 Thematic Mapper (TM) data sets were downloaded and used for analysing the forest coverage of 1990 and 2010. For the years 2000 and 2020, Landsat 7 Enhanced Thematic Mapper Plus (ETM+) image and Landsat 8 Operational Land Imager (OLI) and Thermal Infrared Sensor (TIRS) images were analysed correspondingly. To avoid bias, the imageries were downloaded in the similar available time period. Mostly the immediate time around new leaf flushing was preferred for image downloading. Hence, the downloaded images were deliberately chosen and acquired in February and March based on image availability. In order to get the best results, image data were collected on March 5, 1990, February 5, 2000, February 24, 2010, and March 23, 2020.

Classification

After download, the images were pre-processed and then analysed using ArcGIS 10.3 software. In order to process and classify smoothly, a mask of the forested area was obtained from the image. To understand the severity of forest degradation, the image was broadly classified into only two land coverage, namely, natural forest coverage, and cultivated and bare land. The homesteads were too sporadic and scattered on the agricultural land or under the natural forest coverage. Thus, all homesteads did not cover at least three pixels of the image. That is why the supervised classification was avoided, as it needed at least three pixels for count. A canal-type river was found in the natural forest area and some sporadic ponds were detected throughout the villages. All the waterbodies get dry during February and March each year. Calculation of the area of waterbodies could be biased in this time period. Moreover, the present study tends to find out the forest cover change and the river situated within the forested area. In order to avoid misinterpretation and for proper detection of natural forest area changes, minor classified groups are intentionally evaded. Thus, the settlements were excluded from counting and waterbodies were also deliberately included in the broader natural forest class. In the unsupervised classification, the base map was used for proper analysis and verification. The classification process has been outlined in Figure 1.

FIGURE 1

Figure 1. Flow chart of land use and land cover classification process in Landsat 5, 7, and 8.

Retrieving Land Surface Temperature

Several methods are applicable in retrieving LST. Among all methods, four are mentionable: Mono Window Algorithm (MWA) (Qin et al., 2001; Kumari et al., 2020), Single Channel Algorithm (SCA) (Dong et al., 2017), Radiative Transfer Remote Sensing 2020, 12, 294 5 of 32 Equation (RTE) method (Saratoon et al., 2013), and Split Window Algorithm (SWA) (Rozenstein et al., 2014; Wang et al., 2019) using Landsat 5 Thematic Mapper (TM), 7 Enhanced Thematic Mapper Plus (ETM+), and 8 Operational Land Imager (OLI) and Thermal Infrared Sensor (TIRS) data (García-Santos et al., 2018; Hao et al., 2019). In this study, the MWA method (Qin et al., 2001; Guha et al., 2019; Guha and Govil, 2020; Kumari et al., 2020) was chosen as it is easily applicable to Landsat 5 TM, 7 ETM+ and 8 OLI/TIRS data. Mono-window algorithm method is different from SWA method and more suitable to apply as it retrieves LST from only one thermal band data. Moreover, compared to other methods, where several atmospheric parameters are required for LST retrieval, the mono-window algorithm requires only two: transmittance and mean atmospheric temperature parameters (Qin et al., 2001). This method has the following steps to follow. After retrieving land surface temperature from each images it was classified into five classes.

Conversion of the Digital Number to Top of Atmosphere Spectral Radiance

Top of Atmosphere (TOA) Reflectance is a unitless measurement that provides the radiation ratio reflected from the incident solar radiation on a given surface. The Landsat Thematic Mapper (TM) and Enhanced Thematic Mapper Plus (ETM+) sensors capture reflected solar energy. These data are then converted to radiance and finally, rescale this data into an 8-bit digital number (DN) with a range between 0 and 255 (Qin et al., 2001; USGS, 2021). Top of Atmosphere can be calculated by using the equations (1a) and (1b):

For Landsat 8 OLI/TIRS

$$\mathrm{L}\mathrm{\lambda }=\mathrm{ML}\times \mathrm{Qcal}+\mathrm{AL}-\mathrm{Oi}(\mathrm{1a})$$

where Lλ = Top of atmosphere spectral radiance [Watts/ (m² × sr × μm)], ML = Radiance multiplicative Band (No.), AL = Radiance Add Band (No.), Qcal = Quantized and calibrated standard product pixel values (DN), Oi = correction value for Band 10 is 0.29.

For Landsat 5 TM and 7 ETM+

$${\mathrm{L}}_{\mathrm{\lambda }}=\left(\frac{{\text{LMAX}}_{\mathrm{\lambda }}-{\text{LMIN}}_{\mathrm{\lambda }}}{\mathrm{QCALMAX}-\mathrm{QCALMIN}}\right)\times \left(\mathrm{QCAL}-\mathrm{QCALMIN}\right)+{\text{LMIN}}_{\mathrm{\lambda }}(\mathrm{1b})$$

where L_λ = Spectral Radiance, QCAL = Quantized calibrated pixel value in DN, LMAX_λ = Spectral radiance scaled to QCALMAX in [Watts/ (m² × sr × μm)], LMIN_λ = Spectral radiance scaled to QCALMIN in [Watts/ (m² × sr × μm)], QCALMIN = Minimum quantized calibrated pixel value (corresponding to LMIN_λ) in DN, QCALMAX = Maximum quantized calibrated pixel value (corresponding to LMAX_λ) in DN = 255

Conversion of Spectral Radiance to At-Satellite Brightness Temperatures

Thermal band data can be converted from spectral radiance to top of atmosphere brightness temperature using the thermal constants K1 and K2 (USGS, 2021). Equation (2) can simplify the whole process:

$$BT={\displaystyle \frac{K2}{\ln \left(\frac{K1}{L_{\mathrm{\lambda }}}+1\right)}}-273.15(2)$$

where BT = Top of atmosphere brightness temperature (^oC), L_λ = TOA spectral radiance [Watts/ (m² × sr × μm)], K1 = K1 Constant Band (No.), K2 = K2 Constant Band (No.).

Normalized Different Variance Index

Normalized Different Variance Index (NDVI) measures surface reflectance and widely used to estimate vegetation growth and biomass (Bhandari et al., 2012; Wu et al., 2016). It is calculated as the ratio between TOA reflectance of a red band and a near-infrared (NIR) band (Viana et al., 2019). NDVI can easily be clarified in the equation (3):

$$NDVI={\displaystyle \frac{(\mathrm{NIR}-\mathrm{RED})}{(\mathrm{NIR}+\mathrm{RED})}}(3)$$

where RED = DN values from the Red band, NIR = DN values from Near-Infrared band.

Land Surface Emissivity

Land surface emissivity (LSE) is an important surface parameter (Li et al., 2013). It is the average emissivity of an earth surface element and is calculated from the emitted radiance measured from space and land surface temperature (Norman and Becker, 1995; Pal and Ziaul, 2017). The algorithms can be explained in detail in equations (4) and (5):

$$P{}_{\mathrm{\upsilon }}=[(NDVI-NDVI{}_{min})/(NDVI{}_{max}-NDVI{}_{min})]{}^2(4)$$

where: P_υ = Proportion of Vegetation, NDVI = DN values from NDVI Image, NDVImin = Minimum DN values from NDVI Image, NDVImax = Maximum DN values from NDVI Image.

$$\mathrm{\epsilon }=0.004\times P{}_{\mathrm{\upsilon }}+0.986(5)$$

where ε = Land Surface Emissivity, P_υ = Proportion of Vegetation, 0.004 and 0.986 correspond to correction values of the equation.

Land Surface Temperature

Land surface temperature (LST) is an important reference index used to measure the radiative skin temperature of the land derived from infrared radiation. It has impacts on energy cycles, ecological system balance, and human life (Yan et al., 2020). Procedure of retrieving LST can be elaborated in Equation (6):

$$LST=BT/[1+(\mathrm{\lambda }\times BT/c2)\times ln(E)](6)$$

where BT = Top of atmosphere brightness temperature (^oC), λ = Wavelength emitted radiance, E = Land Surface Emissivity, c2 = h × c/s = 1.4388 × 10^–2mK = 14388mK, h = Plank’s Constant = 6.626 × 10^–34J s, s = Boltzmann constant = 1.38 × 10^–23 JK, c = velocity of light = 2.998 × 10⁸m/s.

Accuracy Assessment

Accuracy assessment is one of the best ways to determine the quality of information derived from remotely sensed data, and to remove confusion about the predicted classification with the actual one (Yusuf et al., 2014; Zhang et al., 2016). In this process, land cover classification results are compared with the geospatial data. In the present study, resulted land cover classes were linked with Google Earth maps to test the mistakes. The pixels derived from the classified map image were compared with the same location in the geospatial map. Thus, accuracy assessment provides an overall accuracy of map as well as accuracy for each class in the map. The overall accuracy percentage was calculated by using following formula:

$$\mathrm{Overall}\mathrm{Accuracy}={\displaystyle \frac{TotalNumberofCorrectlyClassifiedPixels(Diagonal)}{TotalNumberofReferencePixels}}\times 100$$

A total of 118 sample sites from Google Earth were selected for accuracy assessment. In order to calculate the class-based accuracy besides overall accuracy, accuracy of individual classes was calculated for each class in each map. The user’s accuracy and producer’s accuracy were two approaches applied. Producer’s accuracy is the ratio of the number of correctly classified pixels in each category to the total number of reference pixels in that category. It takes into account the error of omission that refers the real land cover type that were erroneously omitted from the classified map. User’s accuracy is calculated by dividing the number of correctly classified pixels in each category by the total number of classified pixels in that category (Pal and Ziaul, 2017). Producer’s accuracy and user’s accuracy are calculated by using the following formula:

$$\mathrm{Producer}\mathrm{Accuracy}={\displaystyle \frac{\mathrm{Number}\mathrm{of}\mathrm{Correctly}\mathrm{Classified}\mathrm{Pixels}\mathrm{in}\mathrm{each}\mathrm{Category}}{\mathrm{Total}\mathrm{Number}\mathrm{of}\mathrm{Reference}\mathrm{Pixels}\mathrm{in}\mathrm{that}\mathrm{Category}}}\times 100$$

$$\mathrm{User}\mathrm{Accuracy}={\displaystyle \frac{\mathrm{Number}\mathrm{of}\mathrm{Correctly}\mathrm{Classified}\mathrm{Pixels}\mathrm{in}\mathrm{each}\mathrm{Category}}{\mathrm{Total}\mathrm{Number}\mathrm{of}\mathrm{Classified}\mathrm{Pixels}\mathrm{in}\mathrm{that}\mathrm{Category}}}\times 100$$

Besides overall accuracy, Kappa coefficient (K) is also used as an alternative measurement proposed by Verma et al. (2020). Kappa coefficient (K) was calculated by multiplying the total correctly classified pixels in the verification classes by the total pixels, from which subtracting the sum of the ground verification pixels in a class (column total) times the sum of the classified pixels in that class (row total) summed over all classes (sum of row total × column total), and dividing by the sum of the ground verification pixels in a class (column total) times the sum of the classified pixels in that class (row total) summed over all classes (sum of row total × column total) subtracting from the square of the total pixels (Pal and Ziaul, 2017). The value of Kappa generally lies between 0 and 1, where 0 denotes agreement due to chance only and 1 denotes complete agreement between the two data sets. If negative values occur, they are spurious (Congalton, 1991). Kappa coefficient (K) is formally expressed as a percentage (%). A true Kappa value is always less than or equal to 1. The more the value means, the more the accuracy. Value equal to 1 represents the perfect agreement. More precisely, poor or very poor agreement between the images determined when Kappa value is less than 0.4, fair agreement determined when values between 0.4 and 0.55, values from 0.55 to 0.70 represent good agreement, values from 0.70 to 0.85 represent very good agreement, and values higher than 0.85 represent excellent agreement (Kikon et al., 2016; Zhang et al., 2016; Pal and Ziaul, 2017; Hua and Ping, 2018). The calculation of Kappa statistic is as follows:

$$\mathrm{Kappa}\mathrm{Coefficient}={\displaystyle \frac{\left(\text{TP}\times \text{TCP}\right)-\sum (ColumnTotal\times RowTotal)}{{\text{TP}}^2-\sum (ColumnTotal\times RowTotal)}}\times 100$$

Where, TP = Total Pixels, and TCP = Total Corrected Pixels.

Results

The forest coverage damage was analysed and detected from the Landsat images at 10 years interval since 1990 to 2020 (Figure 2). Images were classified into two broad classes: natural forest coverage, and cultivated and bare land. In 1990, the natural forest coverage was 85% which declined to 78% in 2000, and 68% in 2010 (Figures 2, 3). The cultivated and bare land increases with time. It is observed that maximum forest damages occurred during 2000–2010 which was lately covered by natural regeneration and with the help of some plantation by Bangladesh Forest Department (BFD) and United States Agency for International Development (USAID) funded projects, i.e., Integrated Protected Area Co-management (IPAC) and Climate Resilience Ecosystems and Livelihoods (CREL) projects by participating local people from 2010 to 2020.

FIGURE 2

Figure 2. Forest cover changes in DDWS since 1990 to 2020.

FIGURE 3

Figure 3. Forest land area (acre) changes in DDWS since 1990 to 2020.

The classification of the Landsat images for the years 1990, 2000, 2010, and 2020 revealed that the wildlife sanctuary area has been experiencing declining trends from 1990 to 2010 and dramatically changed after 2010. Since 2010 to 2020, the forest coverage had been increased and the cultivated and bare land had been decreased in the area (Figures 2, 3).

The analysis of forest coverage revealed a dynamic change in the forested area. Figure 2 shows the significant variation of forest cover changes over the study period. The natural forest coverage indicates a declining rate at first and then increased, while the cultivated and bare land increased previously with the decline of forested area and decreased when the forest coverage improved (Figure 3).

In 1990, the highest forest coverage (20938.72 acre) and lowest uncovered land (3722.90 acre) observed (Figure 3). In 2010, the lowest forest coverage (16883.95 acre) and highest forest degraded and converted land (7773.62 acre) found after analysis. From 1990 to 2000 an amount of 1787.05 acre forest land has been converted to other land uses. During 2000–2010, total converted forest land was 2267.72 acre. Since 2010 to 2020, the scenario is completely reversed. It is observed by analysing Figure 3 that a total of 1514.50 acres of bare land has been converted to forest land and increased the total forest coverage area from 16883.95 acre (2010) to 18398.45 acre (2020).

The overall accuracy and Kappa coefficients for land use and land cover maps for 1990, 2000, 2010, and 2020 were compared with each of the classifications (Table 1). In the present study, unsupervised classification was employed to conduct land cover class analysis. For the years 1990, 2000, 2010, and 2020, the overall accuracy lies between 92.59 and 96.67% and Kappa coefficient lies between 84.66 and 92.68% which signifies the excellent agreement of the classified images (Table 1). Concomitantly, the producer’s accuracy and user’s accuracy for all the classes from 1990 to 2020 exceeded 90% which ensures the highest accuracy at the time when classification was done (Table 1).

TABLE 1

Table 1. Accuracy assessments of land cover classes.

With the time period, the land surface temperature has changed significantly. It may be due to two reasons: change in forest coverage area and global warming. Global warming was considered as the most important factor in this study as after 2010 though the forest coverage increased, the temperature also increased at the same time rather than decreasing. For avoiding the biasness of analysing land surface temperature, the images were collected between the time period February 5 and March 23 based on availability of required images in the website. It was found that the temperature increased from 2.3 to 3.0°C in each decade (Figures 4–6) which is a great threat to any country.

FIGURE 4

Figure 4. Land Surface Temperature maps of DDWS.

FIGURE 5

Figure 5. Land surface temperature (°C) upward trend in each decade at DDWS.

FIGURE 6

Figure 6. Change in temperature (°C) with time in DDWS from 1990 to 2020.

In 1990 (March 5, 1990), the highest temperature found was 24.5°C and the lowest was 17.9°C. In 2000 (February 5, 2000), maximum temperature was 26.9°C and minimum was 18.7°C. During 2010 (February 24, 2010) and 2020 (March 23, 2020), the uppermost value of temperature was analysed 29.2°C and 32.2°C, whereas the lowermost was 21.1°C and 21.6°C, respectively (Table 2 and Figure 6). Mean value of temperature show a gradual upward trend (Table 2 and Figure 5). Overall temperature raised 7.7°C from 1990 to 2020 and changes found a rise of 2.3–3°C in each decade (Figure 7).

TABLE 2

Table 2. Variation of five temperature classes of DDWS in selected time periods.

FIGURE 7

Figure 7. Changing pattern of highest and lowest value of temperature in each decade.

The lowest temperature in 1990 (March 5, 1990) was 17.9°C and in 2020 (March 23, 2020), it was 21.6°C. In 30 years, a drastic change in temperature occurred. The highest value of temperature has changed 7.7°C valuing from 24.5°C in 1990 to 32.2°C in 2020. The lowest value change was 3.7°C esteeming from 17.9°C in 1990 to 21.6°C in 2020 (Table 2 and Figures 6, 7). From the analysis it can be concluded that only local forest coverage change is not responsible for dramatic temperature change in the study area.

Discussion

The relation between land use land cover change and land surface temperature change is complex and multidirectional as land use change affects the local, regional, and global climate. Forest cover change in the Dudpukuria-Dhopachari Wildlife Sanctuary was observed throughout the study time. Decrease in forest cover before 2010 and increase after 2010 due to the declaration of DDWS as protected area by the government changed the whole scenario of the studied forest. An improvement of deforestation scenario in DDWS after declaration was also reported by Rahman and Islam (2021). Forest cover seemed as strongest driver of climate change and variability (Li et al., 2016; Winckler et al., 2019). The present research showed an increase in temperature by roughly 2.3–3.0°C in 10 years in DDWS. This findings support the research of Aik et al. (2021) on Cameron Highlands, where a rise in temperature found 2.0–3.5°C in each decade and the increase was introduced by conversion of forest cover in to urban and urban in to agriculture. Choudhury et al. (2019) also reported similar research findings, where the land surface temperature increased from 0.6 to 4.3°C per decade. The studies by Gohain et al. (2021); Rashid et al. (2021); and Thakur et al. (2021) also found a negative relationship between LST and vegetation coverage. But, surprisingly, in the present study, during the first two decades, though the land surface temperature of DDWS increased with the decrease of forest cover, in the last decade, the land surface temperature increased, however, the forest cover increased.

Forest conversion has a very sensitive relationship with climate change. Commonly, land surface temperature decrease when forest cover increased. Hence, the present study provides a new turn to the existing thinking of the researchers. The present study concluded that local land cover change is not only responsible for climatic change or change of its variability. The impact of local land use and land cover change in the forest area is minor on the land surface temperature change; rather, it is the global warming which is suspected as the main reason behind the LST changes in this study. The results of the study coincide with the findings by Hansen et al. (2010) who researched on urban effects on global temperature change and found a small effect on global temperature. The similar kind of results observed by Parker (2010) and Hansen et al. (2001), who reported a modest effect of urban warming on near surface temperature. Global temperature change also includes the ocean data and polar region data. This means, the local level urbanization or forest cover change has a small impact on LST. Instead, globally, all type of local level land use and forest cover change along with ocean and polar region transition as a whole increase the local, regional, as well as global temperature. The global temperature then affects the LST. Global mean surface temperature (GMST) is poorly correlated with net radiation at the top of the atmosphere (Xie and Kosaka, 2017). It is estimated from the study that only sporadic increase of forest coverage will not be able to reduce local temperature but a combined timely increase of local forest coverage by each country may be able to halt the gradual increase of global temperature as well as local land surface temperature.

Conclusion

The present study detected the forest damages and retrieved land surface temperature change in Dudpukuria-Dhopachari Wildlife Sanctuary by applying GIS & RS for 30 years since 1990 to 2020. The forest cover change revealed a significant damages of forest from 1990 to 2010 and a dramatic forest restoration since 2010 to 2020. The land surface temperature showed the increasing pattern in the studied time period. Since 1990 to 2020, a highest temperature increases found was 7.7°C and in each decade it was 2.3–3.0°C. The present research revealed that local forest cover change is not the only reason behind the land surface temperature change, rather, global warming seems to have major influence for this high and rapid rise in temperature.

Data Availability Statement

The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.

Author Contributions

The author confirms being the sole contributor of this work and has approved it for publication.

Conflict of Interest

The author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

Ahmed, B., Kamruzzaman, M. D., Zhu, X., Rahman, M., and Choi, K. (2013). Simulating land cover changes and their impacts on land surface temperature in Dhaka, Bangladesh. Remote Sens. 5, 5969–5998. doi: 10.3390/rs5115969

CrossRef Full Text | Google Scholar

Aik, D. H. J., Ismail, M. H., Muharam, F. M., and Alias, M. A. (2021). Evaluating the impacts of land use/land cover changes across topography against land surface temperature in Cameron Highlands. PLoS One 16:e0252111. doi: 10.1371/journal.pone.0252111

PubMed Abstract | CrossRef Full Text | Google Scholar

Al-doski, J., Mansor, S. B., Shafri, H. Z. B. M., and Zulhaidi, H. (2016). Thermal imaging for pests detecting a review. Int. J. Agric. For. Plant. 2, 10–30.

Google Scholar

BFD. (2021). Bangladesh Forest Department, Government of the People’s Republic of Bangladesh. Available online at: http://www.bforest.gov.bd/ (accessed March 16, 2021).

Google Scholar

Bhandari, A. K., Kumar, A., and Singh, G. K. (2012). Feature extraction using normalized difference vegetation index (NDVI): a case study of Jabalpur city. Procedia Technol. 6, 612–621. doi: 10.1016/j.protcy.2012.10.074

CrossRef Full Text | Google Scholar

Cetin, I. Z., and Sevik, H. (2020). Investigation of the relationship between bioclimatic comfort and land use by using GIS and RS techniques in Trabzon. Environ. Monit. Assess. 192, 1–14. doi: 10.1007/s10661-019-8029-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Chaudhuri, G., and Mishra, N. B. (2016). Spatio-temporal dynamics of land cover and land surface temperature in Ganges-Brahmaputra delta: a comparative analysis between India and Bangladesh. Appl. Geogr. 68, 68–83. doi: 10.1016/j.apgeog.2016.01.002

CrossRef Full Text | Google Scholar

Chew, C., Shah, R., Zuffada, C., Hajj, G., Masters, D., and Mannucci, A. J. (2016). Demonstrating soil moisture remote sensing with observations from the UK TechDemoSat−1 satellite mission. Geophys. Res. Lett. 43, 3317–3324. doi: 10.1002/2016GL068189

CrossRef Full Text | Google Scholar

Choudhury, D., Das, K., and Das, A. (2019). Assessment of land use land cover changes and its impact on variations of land surface temperature in Asansol-Durgapur development region. Egypt. J. Remote Sens. Space Sci. 22, 203–218. doi: 10.1016/j.ejrs.2018.05.004

CrossRef Full Text | Google Scholar

Chowdhury, M., Hasan, M. E., and Abdullah-Al-Mamun, M. M. (2020). Land use/land cover change assessment of Halda watershed using remote sensing and GIS. Egypt. J. Remote Sens. Space Sci. 23, 63–75. doi: 10.1016/j.ejrs.2018.11.003

CrossRef Full Text | Google Scholar

Chowdhury, M. S. H. (ed.) (2014). Forest Conservation in Protected Areas of Bangladesh: Policy and Community Development Perspectives, Vol. 20. Cham: Springer.

Google Scholar

Chowdhury, M. S. H., and Koike, M. (2010). An overview on the protected area system for forest conservation in Bangladesh. J. For. Res. 21, 111–118. doi: 10.1007/s11676-010-0019-x

CrossRef Full Text | Google Scholar

Congalton, R. G. (1991). A review of assessing the accuracy of classifications of remotely sensed data. Remote Sens. Environ. 37, 35–46. doi: 10.1016/0034-4257(91)90048-B

CrossRef Full Text | Google Scholar

Dewan, A. M., and Corner, R. J. (2012). “The impact of land use and land cover changes on land surface temperature in a rapidly urbanizing megacity,” in Proceedings of the 2012 IEEE International Geoscience and Remote Sensing Symposium, (Honolulu, HI: IEEE), 6337–6339.

Google Scholar

Dong, J., Crow, W. T., and Bindlish, R. (2017). The error structure of the SMAP single and dual channel soil moisture retrievals. Geophys. Res. Lett. 45, 758–765. doi: 10.1002/2017GL075656

PubMed Abstract | CrossRef Full Text | Google Scholar

Dudley, N., Parrish, J. D., Redford, K. H., and Stolton, S. (2010). The revised IUCN protected area management categories: the debate and ways forward. Oryx 44, 485–490. doi: 10.1017/S0030605310000566

CrossRef Full Text | Google Scholar

Fagan, M. E., Reid, J. L., Holland, M. B., Drew, J. G., and Zahawi, R. A. (2020). How feasible are global forest restoration commitments? Conserv. Lett. 13:e12700. doi: 10.1111/conl.12700

CrossRef Full Text | Google Scholar

García-Santos, V., Cuxart, J., Martínez-Villagrasa, D., Jiménez, M. A., and Simó, G. (2018). Comparison of three methods for estimating land surface temperature from landsat 8-tirs sensor data. Remote Sens. 10:1450. doi: 10.3390/rs10091450

CrossRef Full Text | Google Scholar

Gazi, M. Y., Rahman, M. Z., Uddin, M. M., and Rahman, F. A. (2020). Spatio-temporal dynamic land cover changes and their impacts on the urban thermal environment in the Chittagong metropolitan area, Bangladesh. GeoJournal 1–16. doi: 10.1007/s10708-020-10178-4

CrossRef Full Text | Google Scholar

Goetz, S., and Dubayah, R. (2011). Advances in remote sensing technology and implications for measuring and monitoring forest carbon stocks and change. Carbon Manag. 2, 231–244. doi: 10.4155/cmt.11.18

CrossRef Full Text | Google Scholar

Gohain, K. J., Mohammad, P., and Goswami, A. (2021). Assessing the impact of land use land cover changes on land surface temperature over Pune city, India. Quat. Int. 575, 259–269. doi: 10.1016/j.quaint.2020.04.052

CrossRef Full Text | Google Scholar

Guha, S., and Govil, H. (2020). An assessment on the relationship between land surface temperature and normalized difference vegetation index. Environ. Dev. Sustain. 1–20. doi: 10.1007/s10668-020-00657-6

CrossRef Full Text | Google Scholar

Guha, S., Govil, H., and Diwan, P. (2019). Analytical study of seasonal variability in land surface temperature with normalized difference vegetation index, normalized difference water index, normalized difference built-up index, and normalized multiband drought index. J. Appl. Remote Sens. 13:024518. doi: 10.1117/1.JRS.13.024518

CrossRef Full Text | Google Scholar

Halme, P., Allen, K. A., Auniņš, A., Bradshaw, R. H., Brūmelis, G., Čada, V., et al. (2013). Challenges of ecological restoration: lessons from forests in northern Europe. Biol. Conserv. 167, 248–256. doi: 10.1016/j.biocon.2013.08.029

CrossRef Full Text | Google Scholar

Hansen, J., Ruedy, R., Sato, M., Imhoff, M., Lawrence, W., Easterling, D., et al. (2001). A closer look at United States and global surface temperature change. J. Geophys. Res. 106, 23947–23963.

Google Scholar

Hansen, J., Ruedy, R., Sato, M., and Lo, K. (2010). Global surface temperature change. Rev. Geophys. 48, 1–29. doi: 10.1029/2010RG000345

CrossRef Full Text | Google Scholar

Hao, Y., Song, L., Cui, L., and Wang, H. (2019). A three-dimensional geometric features-based SCA algorithm for compound faults diagnosis. Measurement 134, 480–491. doi: 10.1016/j.measurement.2018.10.098

CrossRef Full Text | Google Scholar

Hassan, M. K., Tongchangya, P. K., Khalilullah, M. I., Nahid, M. I., and Feeroz, M. M. (2011). Palla’s squirrel in dudpukuria dhopachari wildlife sanctuary, Bangladesh. Small Mammal. Mail. 3, 36–37.

Google Scholar

Hishe, H., Giday, K., Van Orshoven, J., Muys, B., Taheri, F., Azadi, H., et al. (2021). Analysis of land use land cover dynamics and driving factors in Desa’a forest in Northern Ethiopia. Land Use Policy 101:105039. doi: 10.1016/j.landusepol.2020.105039

CrossRef Full Text | Google Scholar

Hite, K. A., and Seitz, J. L. (2021). Global Issues: An Introduction. Hoboken, NJ: John Wiley & Sons.

Google Scholar

Hossain, M. A., Hossain, M. K., Alam, M. S., and Al Mamun, M. A. (2017). Structural composition and distribution of tree species of Dudhpukuria-Dhopachori Wildlife Sanctuary, Chittagong, Bangladesh. J. Biodivers. Conserv. Bioresour. Manag. 3, 17–30. doi: 10.3329/jbcbm.v3i1.36757

CrossRef Full Text | Google Scholar

Hua, A. K., and Ping, O. W. (2018). The influence of land-use/land-cover changes on land surface temperature: a case study of Kuala Lumpur metropolitan city. Eur. J. Remote Sens. 51, 1049–1069. doi: 10.1080/22797254.2018.1542976

CrossRef Full Text | Google Scholar

Huang, Y., Raza, S. M. F., Hanif, I., Alharthi, M., Abbas, Q., and Zain-ul-Abidin, S. (2020). The role of forest resources, mineral resources, and oil extraction in economic progress of developing Asian economies. Resour. Policy 69, 101878. doi: 10.1016/j.resourpol.2020.101878

CrossRef Full Text | Google Scholar

Islam, K., and Sato, N. (2012). Deforestation, land conversion and illegal logging in Bangladesh: the case of the Sal (Shorea robusta) forests. IForest 5:171. doi: 10.3832/ifor0578-005

PubMed Abstract | CrossRef Full Text | Google Scholar

IUCN, and UNEP-WCMC. (2014). The World Database on Protected Areas (WDPA)[On-line]. Cambridge: UNEP-WCMC.

Google Scholar

Kafy, A. A., Rahman, M. S., Hasan, M. M., and Islam, M. (2020). Modelling future land use land cover changes and their impacts on land surface temperatures in Rajshahi, Bangladesh. Remote Sens. Appl. Soc. Environ. 18:100314. doi: 10.1016/j.rsase.2020.100314

CrossRef Full Text | Google Scholar

Karim, R., and Ahsan, F. (2016). Mammalian fauna and conservational issues of the Baraiyadhala National Park in Chittagong. Bangladesh. Open J. For. 6:123.

Google Scholar

Kaya, İA., and Görgün, E. K. (2020). Land use and land cover change monitoring in Bandırma (Turkey) using remote sensing and geographic information systems. Environ. Monit. Assess. 192, 1–18. doi: 10.1007/s10661-020-08411-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Kayet, N., Pathak, K., Chakrabarty, A., and Sahoo, S. (2016). Spatial impact of land use/land cover change on surface temperature distribution in Saranda forest, Jharkhand. Model. Earth Syst. Environ. 2, 1–10. doi: 10.1007/s40808-016-0159-x

CrossRef Full Text | Google Scholar

Kikon, N., Singh, P., Singh, S. K., and Vyas, A. (2016). Assessment of urban heat islands (UHI) of Noida City, India using multi-temporal satellite data. Sustain. Cities Soc. 22, 19–28.

Google Scholar

Kissinger, G. M., Herold, M., and De Sy, V. (2012). Drivers of Deforestation and Forest Degradation: A Synthesis Report for REDD+ Policymakers. Vancouver: Lexeme Consulting.

Google Scholar

Kumari, B., Tayyab, M., Ahmed, I. A., Baig, M. R. I., Khan, M. F., and Rahman, A. (2020). Longitudinal study of land surface temperature (LST) using mono-and split-window algorithms and its relationship with NDVI and NDBI over selected metro cities of India. Arab. J. Geosci. 13, 1–19. doi: 10.1007/s12517-020-06068-1

CrossRef Full Text | Google Scholar

Lein, J. K. (2011). Environmental Sensing: Analytical Techniques for Earth Observation. Cham: Springer Science.

Google Scholar

Li, Y., Zhao, M., Mildrexler, D. J., Motesharrei, S., Mu, Q., Kalnay, E., et al. (2016). Potential and actual impacts of deforestation and afforestation on land surface temperature. J. Geophys. Res. 121, 372–314. doi: 10.1002/2016JD024969

CrossRef Full Text | Google Scholar

Li, Z. L., Wu, H., Wang, N., Qiu, S., Sobrino, J. A., Wan, Z., et al. (2013). Land surface emissivity retrieval from satellite data. Int. J. Remote Sens. 34, 3084–3127. doi: 10.1080/01431161.2012.716540

CrossRef Full Text | Google Scholar

Liu, T., Mao, P., Shi, L., Eisenhauer, N., Liu, S., Wang, X., et al. (2020). Forest canopy maintains the soil community composition under elevated nitrogen deposition. Soil Biol. Biochem. 143, 107733. doi: 10.1016/j.soilbio.2020.107733

CrossRef Full Text | Google Scholar

Mackey, B. G., Watson, J. E., Hope, G., and Gilmore, S. (2008). Climate change, biodiversity conservation, and the role of protected areas: an Australian perspective. Biodiversity 9, 11–18. doi: 10.1080/14888386.2008.9712902

CrossRef Full Text | Google Scholar

Messina, G., and Modica, G. (2020). Applications of UAV thermal imagery in precision agriculture: state of the art and future research outlook. Remote Sens. 12:1491. doi: 10.3390/rs12091491

CrossRef Full Text | Google Scholar

Mitchell, B. A., Stolton, S., Bezaury-Creel, J., Bingham, H. C., Cumming, T. L., Dudley, N., et al. (2018). Guidelines for Privately Protected Areas. Best Practice Protected Area Guidelines Series, (29). Gland: IUCN, doi: 10.2305/IUCN.CH.2018.PAG.29.en

PubMed Abstract | CrossRef Full Text | Google Scholar

MohanRajan, S. N., Loganathan, A., and Manoharan, P. (2020). Survey on land use/land cover (LU/LC) change analysis in remote sensing and GIS environment: techniques and challenges. Environ. Sci. Pollut. Res. 27, 29900–29926. doi: 10.1007/s11356-020-09091-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Morales-Hidalgo, D., Oswalt, S. N., and Somanathan, E. (2015). Status and trends in global primary forest, protected areas, and areas designated for conservation of biodiversity from the global forest resources assessment 2015. For. Ecol. Manag. 352, 68–77. doi: 10.1016/j.foreco.2015.06.011

CrossRef Full Text | Google Scholar

Mukul, S. A., Rashid, A. M., and Khan, N. A. (2016). Forest protected area systems and biodiversity conservation in Bangladesh. Preprints 1–16. doi: 10.20944/preprints201611.0101.v1

PubMed Abstract | CrossRef Full Text | Google Scholar

Munthali, M. G., Mustak, S., Adeola, A., Botai, J., Singh, S. K., and Davis, N. (2020). Modelling land use and land cover dynamics of Dedza district of Malawi using hybrid cellular automata and Markov model. Remote Sens. Appl. Soc. Environ. 17:100276. doi: 10.1016/j.rsase.2019.100276

CrossRef Full Text | Google Scholar

Norman, J. M., and Becker, F. (1995). Terminology in thermal infrared remote sensing of natural surfaces. Agric. For. Meteorol. 77, 153–166.

Google Scholar

Olorunfemi, I. E., Fasinmirin, J. T., Olufayo, A. A., and Komolafe, A. A. (2020). GIS and remote sensing-based analysis of the impacts of land use/land cover change (LULCC) on the environmental sustainability of Ekiti State, southwestern Nigeria. Environ. Dev. Sustain. 22, 661–692. doi: 10.1007/s10668-018-0214-z

CrossRef Full Text | Google Scholar

Pal, S., and Ziaul, S. (2017). Detection of land use and land cover change and land surface temperature in English Bazar urban centre. Egypt J. Remote Sens. Space Sci. 20, 125–145.

Google Scholar

Parker, D. E. (2010). Urban heat island effects on estimates of observed climate change. Wiley Interdiscip. Rev. Clim. Change 1, 123–133.

Google Scholar

Petropoulos, G. P., Kalivas, D. P., Griffiths, H. M., and Dimou, P. P. (2015). Remote sensing and GIS analysis for mapping spatio-temporal changes of erosion and deposition of two Mediterranean river deltas: the case of the Axios and Aliakmonas rivers, Greece. Int. J. Appl. Earth Obs. Geoinf. 35, 217–228. doi: 10.1016/j.jag.2014.08.004

CrossRef Full Text | Google Scholar

Qin, Z., Karnieli, A., and Berliner, P. (2001). A mono-window algorithm for retrieving land surface temperature from Landsat TM data and its application to the Israel-Egypt border region. Int. J. Remote Sens. 22, 3719–3746. doi: 10.1080/01431160010006971

CrossRef Full Text | Google Scholar

Rahman, M. F., and Islam, K. (2021). Effectiveness of protected areas in reducing deforestation and forest fragmentation in Bangladesh. J. Environ. Manag. 280:111711. doi: 10.1016/j.jenvman.2020.111711

PubMed Abstract | CrossRef Full Text | Google Scholar

Rahman, M. M., Mahmud, M. A. A., Shahidullah, M., Nath, T. K., and Jashimuddin, M. (2016). The competitiveness of the phytosociological attributes of the protected areas in Bangladesh with that in the other tropical countries. J. Sustain. For. 35, 431–450. doi: 10.1080/10549811.2016.1202841

CrossRef Full Text | Google Scholar

Rashid, K. J., Hoque, M. A., Esha, T. A., Rahman, M. A., and Paul, A. (2021). Spatiotemporal changes of vegetation and land surface temperature in the refugee camps and its surrounding areas of Bangladesh after the Rohingya influx from Myanmar. Environ. Dev. Sustain. 23, 3562–3577. doi: 10.1007/s10668-020-00733-x

CrossRef Full Text | Google Scholar

Rodriguez-Galiano, V. F., Ghimire, B., Rogan, J., Chica-Olmo, M., and Rigol-Sanchez, J. P. (2012). An assessment of the effectiveness of a random forest classifier for land-cover classification. ISPRS J. Photogramm. Remote Sens. 67, 93–104. doi: 10.1016/j.isprsjprs.2011.11.002

CrossRef Full Text | Google Scholar

Rozenstein, O., Qin, Z., Derimian, Y., and Karnieli, A. (2014). Derivation of land surface temperature for Landsat-8 TIRS using a split window algorithm. Sensors 14, 5768–5780. doi: 10.3390/s140405768

PubMed Abstract | CrossRef Full Text | Google Scholar

Saratoon, T., Tarvainen, T., Cox, B. T., and Arridge, S. R. (2013). A gradient-based method for quantitative photoacoustic tomography using the radiative transfer equation. Inverse Probl. 29:075006.

Google Scholar

Sekertekin, A., and Bonafoni, S. (2020). Land surface temperature retrieval from Landsat 5, 7, and 8 over rural areas: assessment of different retrieval algorithms and emissivity models and toolbox implementation. Remote Sens. 12:294. doi: 10.3390/rs12020294

CrossRef Full Text | Google Scholar

Shahadat, O., Ahmed, T., and Neogi, A. K. (2014). Confirmation record of a butterfly species Euploea radamanthus radamanthus Fabricius, 1793 (Family: Nymphalidae) from Dudpukuria-Dhopachari wildlife sanctuary, Chittagong, Bangladesh. Int. J. Fauna Biol. Stud. 1, 22–23.

Google Scholar

Sharma, L. K., Gupta, R., and Pandey, P. C. (2021). “Future aspects and potential of the remote sensing technology to meet the natural resource needs”, in Advances in Remote Sensing for Natural Resource Monitoring, eds P. C. Pandey and L. K. Sharma (Hoboken, NJ: John Wiley), 445–464. doi: 10.1002/9781119616016.ch22

CrossRef Full Text | Google Scholar

Singh, S., Tanvir Hassan, S. M., Hassan, M., and Bharti, N. (2020). Urbanisation and water insecurity in the Hindu Kush Himalaya: insights from Bangladesh, India, Nepal and Pakistan. Water Policy 22, 9–32. doi: 10.2166/wp.2019.215

CrossRef Full Text | Google Scholar

Stolton, S., Dudley, N., Avcıoğlu Çokçalışkan, B., Hunter, D., Ivanić, K. Z., Kanga, E., et al. (2015). “Values and benefits of protected areas,” in Protected Area Governance and Management, eds G. L. Worboys, M. Lockwood, A. Kothari, S. Feary, and I. Pulsford (Canberra, ACT: Australian National University Press), 145–168.

Google Scholar

Thakur, S., Maity, D., Mondal, I., Basumatary, G., Ghosh, P. B., Das, P., et al. (2021). Assessment of changes in land use, land cover, and land surface temperature in the mangrove forest of Sundarbans, northeast coast of India. Environ. Dev. Sustain. 23, 1917–1943. doi: 10.1007/s10668-020-00656-7

CrossRef Full Text | Google Scholar

Ullah, S., Ahmad, K., Sajjad, R. U., Abbasi, A. M., Nazeer, A., and Tahir, A. A. (2019). Analysis and simulation of land cover changes and their impacts on land surface temperature in a lower Himalayan region. J. Environ. Manag. 245, 348–357. doi: 10.1016/j.jenvman.2019.05.063

PubMed Abstract | CrossRef Full Text | Google Scholar

USGS. (2021). United States Geological Survey. Avavilable online at: https://www.usgs.gov/core-science-systems/nli/landsat/using-usgs-landsat-level-1-data-product (accessed May 16, 2021)

Google Scholar

Van Schendel, W. (2020). A History of Bangladesh. Cambridge: Cambridge University Press.

Google Scholar

Verma, P., Raghubanshi, A., Srivastava, P. K., and Raghubanshi, A. S. (2020). Appraisal of kappa-based metrics and disagreement indices of accuracy assessment for parametric and nonparametric techniques used in LULC classification and change detection. Model. Earth Syst. Environ. 6, 1045–1059. doi: 10.1007/s40808-020-00740-x

CrossRef Full Text | Google Scholar

Viana, C. M., Oliveira, S., Oliveira, S. C., and Rocha, J. (2019). “Land use/land cover change detection and urban sprawl analysis,” in Spatial Modeling in GIS and R for earth and Environmental Sciences, eds H. R. Pourghasemi and C. Gokceoglu (Amsterdam: Elsevier), 621–651. doi: 10.1016/B978-0-12-815226-3.00029-6

CrossRef Full Text | Google Scholar

Wagner, D. L., Fox, R., Salcido, D. M., and Dyer, L. A. (2021). A window to the world of global insect declines: moth biodiversity trends are complex and heterogeneous. Proc. Natl. Acad. Sci. U. S. A. 118:e2002549117. doi: 10.1073/pnas.2002549117

PubMed Abstract | CrossRef Full Text | Google Scholar

Wan Mohd Jaafar, W. S., Abdul Maulud, K. N., Muhmad Kamarulzaman, A. M., Raihan, A., Md Sah, S., Ahmad, A., et al. (2020). The influence of deforestation on land surface temperature—a case study of Perak and Kedah, Malaysia. Forests 11:670. doi: 10.3390/f11060670

CrossRef Full Text | Google Scholar

Wang, M., He, G., Zhang, Z., Wang, G., Wang, Z., Yin, R., et al. (2019). A radiance-based split-window algorithm for land surface temperature retrieval: theory and application to MODIS data. Int. J. Appl. Earth Obs. Geoinf. 76, 204–217. doi: 10.1016/j.jag.2018.11.015

CrossRef Full Text | Google Scholar

Winckler, J., Reick, C. H., Luyssaert, S., Cescatti, A., Stoy, P. C., Lejeune, Q., et al. (2019). Different response of surface temperature and air temperature to deforestation in climate models. Earth Syst. Dyn. 10, 473–484. doi: 10.5194/esd-2018-66

CrossRef Full Text | Google Scholar

Wu, Y., Li, W., Wang, Q., Liu, Q., Yang, D., Xing, M., et al. (2016). Landslide susceptibility assessment using frequency ratio, statistical index and certainty factor models for the Gangu County, China. Arab. J. Geosci. 9:84. doi: 10.1007/s12517-015-2112-0

CrossRef Full Text | Google Scholar

Xie, S. P., and Kosaka, Y. (2017). What caused the global surface warming hiatus of 1998–2013? Curr. Clim. Change Rep. 3, 128–140. doi: 10.1007/s40641-017-0063-0

CrossRef Full Text | Google Scholar

Xu, W., Xiao, Y., Zhang, J., Yang, W., Zhang, L., Hull, V., et al. (2017). Strengthening protected areas for biodiversity and ecosystem services in China. Proc. Natl. Acad. Sci. U. S. A. 114, 1601–1606. doi: 10.1073/pnas.1620503114

PubMed Abstract | CrossRef Full Text | Google Scholar

Yan, Y., Mao, K., Shi, J., Piao, S., Shen, X., Dozier, J., et al. (2020). Driving forces of land surface temperature anomalous changes in North America in 2002–2018. Sci. Rep. 10, 1–13. doi: 10.1038/s41598-020-63701-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Yuan, Z., Ali, A., Ruiz−Benito, T., Jucker, T., Mori, A. S., Wang, S., et al. (2020). Above−and below−ground biodiversity jointly regulate temperate forest multifunctionality along a local−scale environmental gradient. J. Ecol. 108, 2012–2024. doi: 10.1111/1365-2745.13378

CrossRef Full Text | Google Scholar

Yusuf, Y. A., Pradhan, B., and Idrees, M. O. (2014). Spatio-temporal assessment of urban heat island effects in Kuala Lumpur metropolitan city using landsat images. J. Indian Soc. Remote Sens. 42, 829–837. doi: 10.1007/s12524-013-0342-8

CrossRef Full Text | Google Scholar

Zeng, L., Wardlow, B. D., Xiang, D., Hu, S., and Li, D. (2020). A review of vegetation phenological metrics extraction using time-series, multispectral satellite data. Remote Sens. Environ. 237, 111511. doi: 10.1016/j.rse.2019.111511

CrossRef Full Text | Google Scholar

Zhang, F., Tiyip, T., Kung, H., Johnson, V. C., Maimaitiyiming, M., Zhou, M., et al. (2016). Dynamics of land surface temperature (LST) in response to land use and land cover (LULC) changes in the Weigan and Kuqa river oasis, Xinjiang, China. Arab. J. Geosci. 9, 1–14.

Google Scholar