Multi-layer perceptron-Markov chain-based artificial neural network for modelling future land-specific carbon emission pattern and its influences on surface temperature

Reliable and accurate environmental state prediction can help in long-term sustainable planning and management. Enormous land-use/ land-cover (LULC) transformation has been increasing the carbon emissions (CEs) and land surface temperature (LST) around the world. The study aimed to (i) examine the influences of land specific CEs on LST dynamics and (ii) simulate future potential LULC, CEs and LST pattern of Khulna City Corporation. Landsat satellite images of the year 2000, 2010 and 2020 were used to derive LULC, LST and CEs pattern and change. The correlation between land-use indices (NDBI, NDVI, NDWI) and LST was examined to explore the impacts of LULC change on LST. Unplanned urbanization has increased 11.79 Km2(26.10%) buildup areas and 25,268 tons of CEs during 2000–2020. The calculated R2 value indicates the strong positive correlation between CEs and LST. To simulate the future LULC, CEs and LST pattern for the year 2030 and 2040, multi-layer perceptron-Markov chain (MLP-MC)-based artificial neural network model was utilized with the accuracy rate of 94.12%, 99% and 98.48% for LULC, LST and CEs model, respectively. The simulation shows that by 2040, buildup area will increase to 87.33%, net CEs will increase by 19.82 × 104tons, and carbon absorptions will decrease by 23. 55 × 104tons and 69.54% of the total study area's LST will be above 390C. Such predictions signify the necessity of implementing a sustainable urban development plan immediately for the sustainable, habitable and sound urban environment.


Introduction
Environmental change is currently one of the main causes of concern around the world [1]. With the expansion of population and spontaneous extension of cities, land use patterns and biological systems have changed, prompting the arrangement of metropolitan situated natural difficulties around the globe [2]. Land use and land cover (LULC) type has been transforming enormously because of the many main impetuses. Thus, greenhouse gas (GHG) emissions, carbon emissions, climate change, environmental change, ecological change, and the random condition have been expanding which are making the climate of any region inadmissible for human home [3][4][5]. This is directly and indirectly accelerating global temperature. Currently, spatial GHG concentrations have risen from a CO 2 equivalent of 280 to 450 ppm since the Industrial Revolution, but the proposed limit is 350 ppm [3]. This massive carbon emission and GHG emission have increased the earth's surface temperature with a rate of 0.2 0 C per decade in the last 30 years [4] while the average temperature increasing rate 0.0186 0 C/year has been observed in Bangladesh [5]. IPCC has projected a temperature rise of 1.1-6.4 °C by the end of the twenty-first century and of 1 °C to 1.5 °C by 2050 in Bangladesh [6]. Global warming, surface temperature rise, climate change, ecological imbalance, ice melting and the rise of sea level, etc., have been increasing the significance of research on carbon emissions around the world so that in Bangladesh. Moreover, the increase in carbon emissions (CEs) has detrimental effects on living health and sustainable development.
LULC is considered the most perilous element of the environment which not only affects ecosystems and the environment but also climate and living beings. Human activities such as urbanization and construction activities may have economic benefits but influencing LULC change and affecting sustainability [7], increasing Land Surface Temperature (LST) [8] and also accelerating CEs [9]. The impact of massive CEs could be largely irreversible in the environment for 1000 years after CEs stops. Due to population growth, modernization and for better jobs, living and lifestyle people are moving towards the city and as a result urban area is expanding. Researchers have shown that LULC transformation is one of the main reasons for increasing CEs and the increase in carbon density in the atmosphere. The calculated net CEs due to LULC change are 12.5% of the total anthropogenic carbon emissions during 1990 to 2010 and 33% of total emissions over the last 150 years [10]. The net CEs from LULC dynamics has been calculated mostly in the carbon budget, which was accounted 1.4 (range: 0.4-2.3) PgC/yr, 1.6 (0.5-2.7) PgC/ yr and 1.1 ± 0.7 PgC/yr in the year 1980s, 1990s and from 2000 to 2009, respectively [13,14]. IPCC reported the highest increase in global average temperatures in the last century because of the increase in anthropogenic concentrations of GHG that contributing to the warming of Earth's surface [1].
Topic related to LULC change such as simulation of future potential LULC patterns and their consequences on the environment and ecosystem has recently attracted interest from a wide range of literature [15][16][17][18][19]. In previous studies, researchers have explained the impact of LULC [5,20,21] CEs [1,9,[22][23][24]. No previous study has discussed land use-based CEs, their consequences and effects on the environment in the context of Bangladesh. Moreover, the impact of spatiotemporal CEs on LST dynamics and their relationship is not discussed in any studies. Researchers have illustrated the impact of LULC change in different urban and rural areas on LST in many previous studies [25][26][27][28][29][30][31]. Similar studies have been done in some cities in Bangladesh including Dhaka [32][33][34], Chattogram [11,35], Rajshahi [36] and Cumilla [7]. Few researchers have highlighted the LULC change in Khulna city [37][38][39] but the spatiotemporal change of LST pattern and impact of LULC on LST of the growing city Khulna is not illustrated.
Many scholars simulated the future LULC scenarios for different study areas. To simulate the future LULC transformation, several models, such as Markov chain, Cellular Automata [17], CLUE [19], Agent-based model [18], Multi-Layer Perceptron-Markov Chain (MLP-MC) [36], have been developed. Subedi et al. [40] and Arsanjani et al. [41] explained MLP-MC model as the most effective model for simulating future spatiotemporal changes with highprecision LULC transformation. Limited simulation studies have been found based on LULC and LST but yet no study has been conducted by modeling future CEs pattern. So far yet, no study has been carried out on Khulna City to predict the future LULC, CEs pattern and their impacts on surface temperature.
This study used artificial neural network (ANN)-based MLP-MC model to simulate the future potential LULC, CEs and LST pattern of Khulna city accurately. The present study filled all the gaps and visualized the LULC, CEs and LST change pattern, examined the impact of spatiotemporal LULC change on LST and influences of CEs pattern on LST change. The overall study signifies the importance of the implementation of proper planning regulations and will be helpful to the engineers, policymakers, urban planners and city's responsible authorities to take necessary steps for mitigating the environmental issues in this growing city.

Theory for how LST depends on Carbon emission
Climate model projections show a basic emerging relationship for our future climate: surface warming rises almost linearly with cumulative CO 2 emissions since pre-industrial times [42]. How surface warming increases with cumulative CO 2 emissions is illustrated in equation [4,43]. The global temperature response to increased atmospheric CE is quantified by measurements such as transient climate response and climate sensitivity equilibrium. The relationship between carbon emission and climate temperature was established by Matthews et al. [43] which is shown in Eq. 1.
Where CTR = Carbon-temperature response. ∆T/∆C A = Temperature change per unit atmospheric carbon increase.
∆C A /E T = Airborne fraction of cumulative carbon emissions.
The sensitivity of surface warming depends on carbon and radiative forcing. When the radioactive forcing is affected only by the CE change in the atmosphere, the mean surface temperature ∆T is measured by Eq. 2. If other non-CO 2 radiative contributions also affect the radiative forcing, the mean surface temperature ∆T is measured by Eq. 3 [42]. For more details about these terms, the readers are requested to follow [42,44].
Where ∆L = Change in the total carbon emission.
∆T/∆R = The ratio of the surface warming and the radiative forcing. ∆R/∆L = Radiative forcing and net carbon emissions ratio. ∆R/∆R CO2 = Ratio of radiative forcing and the radiative forcing from atmospheric CO 2 . ∆R CO2 /∆L = Ratio of the radiative forcing from atmospheric CO 2 and change in total carbon emission.
With rising atmospheric CO 2, the radiative heat flux at the sea surface increases logarithmically. IPCC established Eq. 4 in 1990 for measuring the surface temperature (∆T surface ) for CE [45]. Considering CO 2o as the unperturbed concentration of CO 2, where ∆T surface2 × CO 2 is the surface temperature rise for atmospheric CO 2 doubling and ranges from 2 to 4.5 K, with an average of 3 K, from a variety of climate models [45]. The details of the relationship between surface temperature and CE are described in the literature [42,43,45].

Study area
Khulna, the 3 rd largest divisional headquarters and one of the four major cities in Bangladesh are situated in the south-west corner of the country [46]. Sundarbans, the largest mangrove forest (140,000 ha) and one of the UNESCO World Heritage Sites in the world, is located in this division. Due to this, the division is environmentally sound and the city is considered as a green city [47].
The urban population and the amount of urban area of Khulna City Corporation (KCC) have been increasing at a high rate for the last few decades [39,47]. This has been upsetting the environmental and ecological balance in the KCC area [48]. KCC area is 45.154 km 2 , is located in between 24°45′ and 24°54′ north latitudes and in between 89°28′ and 89°35′ east longitudes [61]. About 1.56 million people live in KCC [46]. The city is bounded by the Bhairab river on the north side, the Rupsa River on the middle and the Pasur river on the south side which flows along the east side of the city, the Mayur on the north side and the Hatia River on the south side which flows along the west side of the city (Fig. 1). The average temperature of Khulna is 26.37 °C and is increasing by a rate of 0.005 °C/year. The annual total rainfall is 1630 mm [49]. Figure 2 shows the increasing trend of maximum and minimum average temperature in Khulna city. In the last 20 years, the highest temperature was observed 40.7 °C in April, 2016 and the average minimum atmospheric temperature was observed 22.13 °C, whereas it was 21.87 °C in the last 69 years. Both annual total rainfall and monthly total rainfall have increased in Khulna. After 2000, the annual average total rainfall measured 1875.134 mm.

Dataset preparation
Three Multi-spectral Landsat satellite images were collected from the United States Geological Survey (USGS) for the year 2000, 2010 and 2020 to extract the LULC, CEs and the LST data. The Landsat program consists of a series of satellite missions under joint supervision by NASA and USGS. Landsat satellites have the best land resolution and spectral bands to monitor land use effectively and to record LULC changes as a result of urbanization, climate change, drought, bio-mash changes, wildfire and many other natural or human changes [33,34,36]. These three satellite images were taken in the same month (April) to avoid the seasonal effects during data analysis. Due to less possibility of rainfall in Khulna during April, the accuracy or acceptability of the LST result is higher than in other months. The maximum cloud coverage was set to less than 10% during collecting the images. The Landsat satellite images were of 30 m resolution and dated 17/04/2000, 11/04/2010 and 06/04/2020 have been collected. For slope and elevation, the data were collected from google earth and imported into the GIS environment and utilized in the MLP-MC model.

Land use classification and land cover mapping
For a better quality of the images, the radiometric and atmospheric correction has been conducted and enhanced by applying majority filter techniques in ERDAS Imagine 2014 version. The generation of composite band combinations such as natural color composite, true-color composite, false-color composite is done to classify the LULC type in the study area [36]. Blue, green, red and nearinfrared bands were used for Landsat 5 TM images and Landsat 8 OLI images to find true color during data processing in ERDAS Imagine 2014 and ArcGIS 10.6 version. The land cover type has been classified into five categories (Table 1) based on the proper understanding of LULC  dynamics and applied maximum likelihood supervised classification method [50,51] to identify the land cover pattern in the study area for the selected year. Bands were used 4-3-2 (Color Infrared) and 3-2-1 (True Color) for Landsat 5 images, 4-3-2 (True Color) and 5-4-3 (Color Infrared) for Landsat 8 images. Then the land cover type change direction during the study period 2000 to 2020 has been analyzed in ArcMap 10.6 version. To check the accuracy of the classification, more acceptable and best quantitative image classification accuracy measurement techniques named Kappa statistics and confusing matrix were calculated [52,53]. Using Eqs. 5-8, Kappa coefficient, user accuracy, producer accuracy and overall accuracy of each LULC were calculated. Around 150 training sites have been randomly chosen for each image to ensure that each LULC type is covered by all spectral groups. Stratified random sampling method has capability to reduce biasness of accuracy assessment by taking equal number of sample point from each feature classes. BU = Build up area, VG = Vegetation land, VL = Vacant land, AL = Agricultural land, WB = Waterbody.

Carbon emission estimation
The research focused on estimating land-specific carbon emissions, and emission change during the year 2000 Num. of correctly classified pixels in each category Total numb of reference pixels in each category (row total) × 100% Total num of corrected classified pixels (diagonal) Total num of reference pixels × 100% to 2020. Equation 9 was used to estimate carbon emissions from each LULC type [11]. δi 's positive values indicate carbon emissions while negative values indicate carbon absorption [11]. Different LULC types were identified in this study through a supervised image classification method while the value of the carbon emission coefficient for each LULC type has been proposed in previous studies. The carbon emission coefficient for different land cover types in Table 2 has been collected from previous studies [4,54,55]. Though there are   [54][55][56][57]. Net carbon emissions were also calculated by subtracting absorptions from emissions.

Derivation of LST
The land surface temperature has been estimated for the year 2000-2020 by using Landsat Thermal bands. Landsat Here P v is the vegetation proportion calculated following Eq. (14) To obtain the LST values in Celsius (°C), 273.15 was extracted from the initial values (K).
In Landsat 5 images, band 6 is the thermal band which is used for the derivation of LST. Following steps were followed to calculate the LST from Landsat 5 images [36,58].

Land use indices analysis
To identify the influence of different LULC changes on LST, three land-use indices namely Normalized Difference Vegetation Index (NDVI), Normalized Difference Water Index (NDWI) and Normalized Difference Buildup Index (NDBI) are used. These all indices were mostly used in many previous researches. The readers are referred to Guha et al. [59] for the detailed methodology of NDBI, for NDVI and NDWI indices Grigoraș & Urițescu 2019; Kafy et al. 2020 [36,58] for NDWI index. The author used these indices for regression analysis with LST to identify the impact of different LULC types on LST.

Simulations of future scenario
To predict the future LULC scenario, Land Change Modeler (LCM) of TerrSet Geospatial Monitoring and Modeling Software and QGIS was used in MLP-MC neural network method. Because of its training guidelines, it is also called a "back-propagation" network. The LCM tool is used to evaluate losses and gains, detect categories, net shifts [36]. The MLP-MC model contains both static images contain digital numbers. The following steps were followed to estimate the LST [58].
Here L λ = TOA Spectral Radiance (W/ (m 2 × sr × μm)). M L = Radiance multiplicative scaling factor for the band. A L = Radiance additive scaling factor for the band. Q CAL = Quantized calibrated pixel value in Digital Numbers (DN).
In the second step the TOA spectral radiance (L λ ) values are converted into At-Satellite Brightness Temperature (TB) Here TB = At-Satellite Brightness Temperature, Kelvin (K). K 1 , K 2 = Thermal conversion constants for the band.
The TOA Brightness Temperature converted to LST (In Kelvin) using the formula 12.
and dynamic aspects of high-precision LULC transformation and is commonly used by many previous studies for LULC prediction [60,61]. Several factors influence the change of land cover types. In this study, 13 factors (Table 3) were considered. For the prediction of the LULC scenario, the LULC maps of the year 2000 and 2010 were used as independent variables in QGIS. The simulated data for 2020 and classified LULC data for 2020 were compared. The model provided the best outcomes in the study area for which, MLP-MC simulation for 2030 and 2040 was carried out and also determined the possible changes.
The urban disturbance and the distance of locations from the road were calculated from the GIS data, collected from the KCC authority. The data used in the other factors in this study were derived from landuse classification. These factors in Table 3 influence to change LULC pattern [36,62]. The CA-MC transition matrix data were utilized as influential factors (trend of one LULC type to another LULC type) for simulating LULC change. The LULC transition matrix for the year 2000 to 2020, the LULC were simulated for every 10 consecutive years. The readers are referred to follow the literature [36,[60][61][62] for the details of the MLP-MC-based ANN process and influencing factors in Table 3. To check the accuracy of the simulated model, percentage error, accuracy assessment and RMSE were calculated. For model validation, a statistical similarity was performed by comparing the predicted and observed results of LULC of the year 2020. Using this MLP-MC model, the carbon emissions and LST also simulated for the year 2030 and 2040. The higher value of accuracy and lower value of percentage error, the RMSE value indicates the best-fitted prediction model [10,19]. Figure 3 represents the methodological flow that was followed to fulfill the aim of this study.

Accuracy assessment
The accuracy assessment was carried out by calculating the overall accuracy, producer's accuracy, user's accuracy and Kappa Statistic value using Eqs. 5-8 (

Land cover mapping
Using the supervised image classification method, the LULC data were extracted (Table 5) and prepared map for the year 2000, 2010 and 2020 (Fig. 4). Table 5 shows that, in 2000 vegetation land covered almost 30.42% of the total area, then periodically agricultural land 11.46%, waterbody 5.08% and vacant area 0.13%. Though buildup areas were the maximum coverage of LULC type (52.19%) in 2000, unplanned urban expansion and massive construction lead to an increased buildup area to 67.97% in 2010 and 79.01% in 2020. Such an increase in the buildup area has put enormous pressure on other land cover types.
In 2020, agricultural land declined to 3.91%, vacant land 0.04% and waterbodies 4.92% in the study area. This declination of environment-friendly LULC areas such as vegetation, waterbodies and agricultural land has resulted in several environmental problems at KCC and accelerated the carbon emissions and LST. Figure 4 shows that the southern part of the KCC has the highest concentration of buildup area and the lowest presence of other land cover types. As the growth center of the KCC is located in the southern corner and maximum industry and factories are located in this area, urban expansion mostly took place in this place. That's why the landscape in this portion has experienced the most LULC transformation in the last two decades. Several factors, such as massive urban growth, unplanned urban population growth, and rural-urban migration tendency led to such changes.

LULC transformation direction analysis
The LULC change direction and transition matrix were identified using the CA Markov Model in the GIS environment. LULC change direction analysis (Fig. 5, 6) shows that vegetation area decreased by 18.3%, waterbodies 0.2%, vacant land 0.1% and agricultural land 7.55% during the year 2000 to 2020 and these contributed to the increase in buildup area by 26.1%. Land-use change direction shows the increase in agricultural land from 11.46% to 14.94%, vacant land from 0.13% to 1.37%, buildup area 52.91% to 67.97% and declination of vegetation area from 30.42% to 11.95% and waterbody area from 5.08% to 3.77% of the total study area during 2000 to 2010. In the next decades,  the agricultural land cover decreased by 4.98 sq.km (11%) and buildup area increased by 4.98 sq.km (11%). Though agricultural and vacant land increased during 2000 to 2010 but decreased in the next decades by 11% and 1.33%, respectively. Figure 5 shows that due to urban expansion during the study period, different LULC types have been transformed into buildup areas which were quantified by the CA Markov model and presented in Table 6. The LULC change matrix in Table 5 shows the significant increase in the buildup area from 23.88 Km 2 to 35.67km 2 (52.91% to 79.01%), where agricultural land was 3.6796 km 2 , vegetation land 9.3197 km 2 . Only a few portions remained unaltered.

Spatiotemporal carbon emissions estimations
Using LULC statistics in Table 5, total CEs, carbon absorptions and net carbon emissions were estimated for the year 2000, 2010 and 2020 through utilizing Eq. 9 and presented in Table 7. Figure 7 represents the spatiotemporal carbon emissions and absorptions pattern in the study area. Negative figures refer to carbon sink or absorptions, while the positive figures refer to carbon emission sources ( Table 7). The buildup area and vegetation area were the largest contributors to carbon emissions and absorptions, respectively (Table 7). Buildup areas found to contribute more than 85% of total CEs and vegetation area accounted for the absorptions of more than 75% of total carbon absorptions over the last two decades. In 2000, total CEs estimated 9430.6 tons from

LST change analysis
Spatiotemporal LST distribution was extracted from Landsat Thermal bands for all the study year. The LST value was divided into five ranges (with 4 0 C interval) (19 °C-25 °C, 25-29 °C, 29 °C-33 0 C, 33 °C-37 0 C and ≥ 37 °C) based on the extracted lower and higher LST values in different years (Fig. 8) and calculated the percentage of areas within each range (Table 8). Table 8 shows that about 32.33% (14.60 Km 2 ) and 67.47% (30.46 Km 2 ) area was in the range of ≤ 25 °C and between 25 °C and 29 °C in 2000, while only 15.04% (6.79Km 2 ) land cover was found in between 25 °C and 29 °C in 2010. In 2020, this percentage reduced to 0%. During 2010, the LST of most of the areas (77.20%) was in the range of 29 °C-33 °C. During 2000, the entire KCC area had the LST below 30 °C while in 2020, the minimum LST was 28.76 0 C with most of the area (53.70%) having LSTs in the range of 33 °C-37 °C and LST of 46.08% of total areas was in the range of 29 °C-33 °C. Rapid LULC transformation leads to increase buildup areas that increase impermeable layers which trigger LST and climate change.

Association between LULC types and LST
The average, minimum and maximum LST on different LULC types were extracted and presented in   Figure 9 shows the increasing trend of mean LST of all LULC types in the study area during 2000-2020.
The local atmospheric conditions also influence LST change. An increase in atmospheric temperatures leads to a change in LST. The study used three land cover indices namely NDVI, NDWI and NDBI, and determined the correlation between LST and these indices to examine the influence of different LULC changes on LST change.   Figure 10(A-F) shows the negative correlation between LST and NDVI, LST and NDWI which implies that the declination of waterbody and vegetation influenced the increase in LST in the study area. Figure 10(G-I) shows the positive correlation between LST and NDBI which indicates that the increase in buildup area increased the LST during the study period. The study found the lower value of NDVI and NDWI with the higher value of LST. Both NDVI and NDWI values were decreased during the study period that is the density of vegetation, the presence of moisture in the soil, waterbody areas were decreased. On the other hand, the NDBI value increased gradually and led to higher LST values. The higher correlation-coefficient value of the regression analysis between LST and NDVI (0.9583, 0.9693, 0.9494), NDWI (0.9515, 0.919, 0.9658) and NDBI (0.979, 0.9828, 0.9366) indicates that there is a strong significant influence of LULC change in the change of LST. The relation between LST and two indices (NDVI, NDWI) indicates that the higher the vegetation and waterbody density, the lower the surface temperature. Dense forests that mean higher NDVI areas cause high evapotranspiration. The relatively high R 2 values between LST and NDBI (Fig. 10) than other indices prove that increase in the buildup area has more effect on the LST growth of KCC.

Association between carbon emissions and LST
GHG emissions influence global warming. Surface temperature increases with the increase in carbon emissions [43]. This study assessed the impacts of carbon emissions on the LST change. In this communication, we have examined the correlations between CEs and LST in the GIS environment (Fig. 11). Surface temperatures showed an increasing trend with the increase in carbon emission. The lowest LST observed 19.724 °C in carbon sinks and 27.6 °C at the place with the CEs value 5.2 × 10 4 tons/ year. The correlation value was found highest (0.9979) for the year 2010. The mean, minimum and maximum LST (Fig. 8, 9) and the net CEs also increased mostly (Table 7)  The relationship between CEs and LST implies that the increasing trend of CEs in the study area is influenced to increase LST. The increase in urban population led to the increase in urban activities which accelerated LULC transformation and carbon emissions. The reduction of carbon sink areas on large scale has reduced the capacity of carbon absorptions in Khulna and increased LST.

Simulation of LULC for the year 2030 and 2040
The MLP-MC model is used to predict the future potential LULC pattern for 2030 and 2040 (Fig. 13) by using the LULC pattern of 2000-2020. The MLP-MC is a twofold process that produces potential transition maps to buildup areas for different LULC groups. Firstly, the MLP-MC model was used for simulating LULC trends in 2020 using the LULC maps of the years 2000 and 2010 to obtain a reliable and acceptable prediction result. Table 10 shows that, with a large number of iterations, all the potential LULC types showed an overall accuracy of 94% and maximum spatiotemporal match with predicted and observed LULC classes (Fig. 12), with an R 2 value of more than 0.85 and all the land cover types with more than 97% except vacant land cover (73.72%).
The predicted LULC maps for the year 2030 in Fig. 13(A) and 2040 in Fig. 13(B) showed the significant increase in buildup areas and the reduction of other land covers (Table 11). The CA Markov LULC transformation matrix for 2020-2040 is presented in Table 12. By 2030, the buildup land cover area will increase to 35.91 Km 2 (79.55%) and by 2040, it will increase to 39.42 sq.km (87.33%). Compared to 2020, 23.38% of waterbodies and 68.48% of vegetative areas will likely be transformed to buildup areas by 2040. Agricultural land, vacant land, waterbody and vegetation lands will cover 3.72%, 1.37%, 3.81% and 3.77% of the total KCC area, respectively, by 2040 in KCC (Table 11). Between the year 2030 and 2040, the vegetation land cover will decrease from 11.95% to 3.81%. The predicted land cover area demonstrated the difference between the land cover area between the year 2000 and 2040 are, agricultural land −3.597 Km 2 (7.75% decrease), buildup area + 15.54 Km 2 (34.42% increase), vacant land + 0.56 Km 2 (1.24% increase), vegetation −12.01 Km 2 (26.61% decrease) and waterbodies −0.589 Km 2 (1.30% decrease). within the next few years, KCC will face significant environmental degradation. Proper land-use planning, urban green areas and natural resource management will help to make KCC a sustainable town in the future.

Simulation of carbon emissions pattern for 2030-2040
During 2000-2020, net CEs had increased at a significant level, which is shown in Fig. 7, Table 7. Therefore, it is essential to simulate future carbon emissions and absorptions (CEA) pattern. The CEA pattern during 2000-2020 is used in the MLP-MC model to predict the future potential spatiotemporal CEA patterns for the year 2030 and 2040 which are illustrated in Fig. 14. The MLP-MC model was used for simulating the CEA pattern in 2020 using the CEA change during 2000-2010 to obtain a reliable and acceptable prediction result. The predicted 2020 CEA data were compared with the actual data for validation (Table 13). The 98.48% accuracy of the predicted data indicates the good accuracy of the CEA prediction model. The RMSE value is relatively less and the R 2 values indicate the good fit of the model [36]. The modeled CEA in the KCC area for 2030 and 2040 shows the increase in net carbon emissions. Table 14 shows that, within 2030, the estimated total emissions will be 10.06 × 10 4 tons and absorptions will be 1.67 × 10 4 tons/ year. Build up area will contribute to 97.08% (9.77 × 10 4 tons) and agricultural land to 2.98% (0.29 × 10 4 tons) of total CEs in the year 2030. The CEs from the buildup area will increase to 10.73 × 10 4 tons and from agricultural land, it will increase to 0.32 × 10 4 tons by 2040. By 2040, LULC will lose its absorption capacity and the total absorptions will be reduced to 8100 tons/year which is 23. 55 × 10 4 tons less than the amount of carbon absorbed by the different LULC types in 2020. An enormous increase in the buildup   Since the CEs simulation was performed based on the previous (2000-2020) trend, the dominations of higher CEs from buildup areas affected the simulation significantly. The greeneries or vegetative areas contribute the most to carbon absorptions [32], and this study found the massive declination of vegetative land cover in KCC in the future. This will affect the carbon absorption capacity of the study area and accelerate carbon emissions. One of the effective ways to increase carbon sequestration capacity and reduce CEs by increasing green space through plantations and preserving existing vegetative areas and waterbodies.

Simulation of surface temperature for the year 2030 and 2040
During 2000 to 2020, a remarkable increase has been observed in LST. Therefore, LST was simulated for 2030 to 2040 using the MLP-MC model which are illustrated in Fig. 15, Tables 15 and 16. Calculating R 2 , differences, accuracy and RMSE, the accuracy of the simulation was validated by the observed and predicted LST values for 2020 (  (Table 16).
Simulated LST for the year 2030 and 2040 shows the variation of mean, maximum and minimum LST at a   The simulated LST variations showed how the temperature of the KCC area will increase in the future. The increase in LST is influenced by increased CEs which pose a danger to the ecosystem and risk to human health. It will also degrade the standard of the urban environment by contributing to the depletion of the ozone layer and formulation of acid rain, which will affect the groundwater, infrastructure, crop and harms human health as well. To mitigate this CEs effect, a more sustainable and effective way is to increase urban green space or urban trees and recovery of waterbodies for ensuring ecological and environmental sustainability.

Conclusion
For the past two decades, the major cities of Bangladesh especially Khulna have undergone an enormous expansion of urban area resulting in the declination of vegetation, waterbodies and agricultural land cover areas. This study explored the use of artificial neural network (ANN) to simulate and forecast future potential LULC, CEA and LST pattern from a sequence of past three years  The simulation model stated the increase in buildup areas to 87.33% by 2040, followed by 3.81% vegetation, 3.77% waterbodies and 3.72% agricultural land. This will lead to increase net CEs of 19.82 × 10 4 tons/year and reduced carbon absorptions capacity by 23.55 × 10 4 tons/ year by the year 2040. The simulated LST represents that the LST range will be observed from 35.11 °C to 44.43 °C and in 2040, about 59.96% (27.07Km2) area will likely experience the LST range of 39 °C to 43 °C. The mean, minimum and maximum LST will increase by 4.798 °C, 6.0041 °C and 8.2795 °C, respectively, by the year 2040. This phenomenon is a clear indication of environmental unsustainability in the future due to unplanned LULC transformation in buildup areas by destroying vegetation areas and waterbodies. To slowdown this environmental degradation, the attachment of the concept of decentralization, urban greeneries and conservation of natural resources in future eco-friendly urban development plans and its proper implementation can be the potential solutions. The inclusion of this study's output will be helpful to the policymakers and responsible authorities to make Khulna city environmentally sound, inclusive, and sustainable. In this study, only land specific carbon emissions were considered for estimating CEs among the various sources of CEs. Future research may focus on how to increase the carbon sequestration capacity of existing LULC types.

Compliance with Ethical Standrads
Conflict of interest On behalf of all authors, the corresponding author states that there is no conflict of interest.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creat iveco mmons .org/licen ses/by/4.0/.