Inﬂuence of Urban Sprawl on Microclimate of Abbottabad, Pakistan

: Urban centers are expanding rapidly due to quickly-increasing population which results in microclimate change due to lack of urban planning. Factors like degradation of green areas, trees, and vegetation are deﬁning a new regime of urban environment giving rise to a lack of drinking water and lowering water tables. Consequently, unplanned urban sprawl with all its varied facets is having adverse impacts on the environment. Rapid variations in some of its climatic factors in the immediate vicinity are alarming and need to redress at war footings. This paper is an endeavor to present a hypothesis that urban sprawl plays a vital role in impacting and the inﬂuencing the microclimate of the city or the area. In this research, geospatial approaches were adopted to identify urban sprawl. Rise in land surface temperature (LST) for the last 25 years (from 1990 to 2016) has been highlighted using Landsat (5 TM, 7 ETM+ and 8 TIRS) satellite images. Weather data collected from Meteorological Department is used to identify temperature rising trends. Result of the research clearly indicates that rapid urban sprawl has adverse impacts on microclimate by increasing LST. Built-up area has been changed from 178 to 477 sq. km within the years of 1990 to 2016. Vegetated area has decreased from 770 to 602 sq. km. Water areas have been decreased from 524 to 360 sq. km within the span of last 26 years. The bare land has increased from 494 sq. km to 742 sq. km. LST in winters has risen from 17 ◦ C to 23 ◦ C and at times shot up to 31 ◦ C, which is alarming. Annual mean air temperature increased by 3 ◦ C to 4 ◦ C from the 1980s to 2016. During summers, the average air temperature rose to a startling 33 ◦ C from 28 ◦ C and LST has had a steady rise from 28.4 ◦ C to 35 ◦ C. The alarming urban sprawl in relation to temperature rise warrants measures which are required to plan the urban planning, forestation, and ecological measures to mitigate the disastrous trends which may result in ﬂash ﬂoods, landslides, soil erosion and sediment transport endangering downstream reservoirs, water quality and depletion in ground water table in the valley.


Introduction
Microclimate change due to urban sprawl or any other reason is not a far-off story for Pakistan. It is an imminent threat as we are already experiencing its undesirable effects in towns, cities, and the country as a whole. The recent examples are the experiences of destructive floods since last decade, untimely rainfalls, hailstorms, extreme snowfall, and heat waves are all because of fact that nature is now reacting to human-caused damages. This is resulting in the form of loss of human lives, and economic and social losses. [1].
Urban sprawl is one of the most noteworthy issues of sustainability challenge for the development of cities which has a great direct or indirect impact on the micro climate This research addresses following research questions: (1) Why this research is important for Abbottabad? (2) Was the hypothesis proved correct in its results? (3) Will this research be beneficial for urban managements of Abbottabad?
Sequel to above, causes of urban sprawl and its effects on microclimate brought about various consequence of micro climatic change. The sprawl occurred due to population explosions hailed from far off areas of Khyber Pakhtunkhwa Province settled in Abbottabad. This research will be a useful guide for planning departments and policy makers to monitor further city expansion for conserving ecological system and mitigating future microclimatic conditions.

Study Area-Abbottabad
District Abbottabad was named after its founder first Deputy Commissioner Major James Abbot in 1853. This valley was the official headquarters of district of Hazara and later on became division [15]. Abbottabad is situated at an elevation of 1256 m above sea level in the north of Khyber Pakhtunkhwa (KP) province of Pakistan. Area wise, it is spread over 1967 square kilometers. It is surrounded by Mansehra District, Muzaffarabad District, Rawalpindi District, and Haripur District from north, east, south, and west respectively. It is located at a distance of 120 km from Rawalpindi, 217 km from Peshawar, and 84 km from Kohala Bridge. Map of study area is shown at Figure 1. Abbottabad is mostly a rugged mountainous region with mean maximum temperature of 23 • C and minimum temperature of 11 • C. Altitude of the region varies from 600 m to 2800 m and its climate varies accordingly. In winters, mostly snowfall is experienced above the height of 900 m which includes areas of Galyat, Barian, and sometimes in the city of Abbottabad. Monsoon rainfall occurs from the start of July till September and then a dry spell remains for two months afterwards. The River Kunhar, a famous stream, flows sideways from the north eastern boundary of the district and it links in the east with another major stream: The River Jhelum. A number of streams flow through district Abbottabad and most significant amongst them are Harro, Dor, and Siran. Besides the study area, two major parts of methodology are adopted i.e., Urban Sprawl (including NDWI, NDVI, and NDBI) and methods of LST extractions from multispectral images.
Land 2021, 10, x FOR PEER REVIEW 3 of 22 This research addresses following research questions: (1) Why this research is important for Abbottabad? (2) Was the hypothesis proved correct in its results? (3) Will this research be beneficial for urban managements of Abbottabad?
Sequel to above, causes of urban sprawl and its effects on microclimate brought about various consequence of micro climatic change. The sprawl occurred due to population explosions hailed from far off areas of Khyber Pakhtunkhwa Province settled in Abbottabad. This research will be a useful guide for planning departments and policy makers to monitor further city expansion for conserving ecological system and mitigating future microclimatic conditions.

Study Area-Abbottabad
District Abbottabad was named after its founder first Deputy Commissioner Major James Abbot in 1853. This valley was the official headquarters of district of Hazara and later on became division [15]. Abbottabad is situated at an elevation of 1256 m above sea level in the north of Khyber Pakhtunkhwa (KP) province of Pakistan. Area wise, it is spread over 1967 square kilometers. It is surrounded by Mansehra District, Muzaffarabad District, Rawalpindi District, and Haripur District from north, east, south, and west respectively. It is located at a distance of 120 km from Rawalpindi, 217 km from Peshawar, and 84 km from Kohala Bridge. Map of study area is shown at Figure 1. Abbottabad is mostly a rugged mountainous region with mean maximum temperature of 23 °C and minimum temperature of 11 °C. Altitude of the region varies from 600 m to 2800 m and its climate varies accordingly. In winters, mostly snowfall is experienced above the height of 900 m which includes areas of Galyat, Barian, and sometimes in the city of Abbottabad. Monsoon rainfall occurs from the start of July till September and then a dry spell remains for two months afterwards. The River Kunhar, a famous stream, flows sideways from the north eastern boundary of the district and it links in the east with another major stream: The River Jhelum. A number of streams flow through district Abbottabad and most significant amongst them are Harro, Dor, and Siran. Besides the study area, two major parts of methodology are adopted i.e., Urban Sprawl (including NDWI, NDVI, and NDBI) and methods of LST extractions from multispectral images.

Urban Sprawl
Urban sprawl is extracted from the multispectral images by using following methodology.

Data Collection and Pre-Processing
Satellite data was obtained by the website (www.earthexplorer.usgs.gov) owned by USGS. It provides data free of cost. TIR imagery of Landsat 5, Landsat 7 ETM+ and Landsat 8 (OLI) with path 150, row 36 was collected. Object-based classification was applied on the subset of study excluding TIR band using ERDAS Imagine including layer staking. ENVI-2014 Software was used for pre-processing. Band combination of 7-2-1, 4-3-2 and 1-4-5 were used to detect urban, vegetation and water zones in object-based classification in the study area from 1990-2016 with the interval of five years length. Resolutions of the images was 30 m. Complete details of the various images are depicted in Table 1, details of constant values for LST is given in Table 2 and processes of methodology till final output of maps are shown in a flow chart at Figure 2. However, numerical outputs of different indices are shown in Table 3.

Format
Geo Tiff Geo Tiff Geo Tiff ENVI-2014 Software was used for pre-processing. Band combination of 7-2-1, 4-3-2 and 1-4-5 were used to detect urban, vegetation and water zones in object-based classification in the study area from 1990-2016 with the interval of five years length. Resolutions of the images was 30 m. Complete details of the various images are depicted in Table 1, details of constant values for LST is given in Table 2 and processes of methodology till final output of maps are shown in a flow chart at Figure 2. However, numerical outputs of different indices are shown in Table 3.   , 1990, 1995, 2000, 2005, 2010, and 2016 are taken in which urban increase is distinguished besides other common indices. The range of land and water mask (LWM) varied from 0-255. The range of water was observed from 0 to 57; however, the range of shallow water was observed from 28 to 37. The water classes were extracted utilizing Equations (1) and (2) [16]. These equations were perfectly suitable for separating water class from the image and mapped into six maps as shown into the Figures 3-5.

NDVI (Normalized Difference of Vegetation Index)
NDVI is calculated from Equation (3). In NIR Band (band 4), it has been witnessed that reflection of vegetation was much higher than that of visible bands. For the extraction of vegetation index, red band was deducted from NIR and then divided by sum of NIR and red bands in all Landsat images, which represent the ratio from −1 to +1 and mapped in  in the dense green vegetation [17]. It was observed that high reflectance remained persistent in NIR band. Usually, values of NDVI remain between −1 and +1 [18]. NDBI was derived from Equation (4). Built-up area was calculated by deducting vegetation and water indices from the images. In e-Cognition software, NDBI was calculated from Landsat TM sensors Band 5 value is subtracted from Band 4. For this whole process, object-based classification was used, where each and every object of the built-up area was considered to classify and exported in shapefiles with the calculation from following equations [19]. Built-up area is derived from Equation (5).

Object Based Image Classification
To extract different classes from the six images of different years, object-based images classification was used by incorporating eCognition software. This classification depends upon geometry, structure information and spectrum character [20]. There are two basic approaches or steps which are used to classify an image, i.e., segmentation and classification of segmented image through classifier [21], these are:

a. Segmentation
Multiresolution Segmentation algorithm was used for multiscale segmentation. To produce polygons, the images were disintegrated into many segmentations by determining their layer weight, shape, scale parameter and compactness [22]. In segmentation, layer weight was set as "1" for each layer, scale parameter was "25", shape of each segment was "0.2" and compactness was "0.8".

b. Classifier
After setting up the semination parameters, classification method was adopted using KNN Classifier. Detail of hyperparameters of the classifier was N_neighbor: "default" (default parameter was used = 5), Weight was "uniform", auto algorithm, Metric was "minkowski" (where values of "P" was 2), Metric_params were "none" and N_jobs: "none". After that final classified images were established.

Overall Accuracy Assessment and Kappa Coefficient
Accuracy assessment was measured to assess quality of results of any classification which gives us evidence of reliability of classifier. Accuracy assessment was governed by validation based on "Ground Truthing" and "Referenced Points" [23] which are depicted in Table 4. Accuracy assessment and Kappa Coefficient were estimated with the help of ArcGIS software and statistics. Confusion metrices were established after object-based classification [24]. Overall an average of 160 points were earmarked on each classified image and the same points were validated on the ground by using GPS. Kappa coefficient was a discreet multivariate technique that was used for accuracy assessment [25]. These methods validate the results. Equations (6) and (7) were used for calculation of Kappa Coefficient (as mentioned below).

. Data Collection and Compilation
Data pertained to thermal infrared images (TIR) of Landsat 5 TM to Landsat 8 (OLI) was used in this study for LST. Data were procured from USGS resources from 1990 to 2016 with the interval of five years; the collection of data was carried out from 1994-2000 from Landsat 5 TM collection and from 2001 to 2013 from Landsat 7 ETM+ in the same pattern. The details of images for months of Jan 1994-Jan 2016 and Jun 1994-Jun 2016 are given in Table 2. Images of complete months of Jan and Jun were downloaded then all the images of each month were stacked into one image. Then images of complete five years were stacked and used for average of LST of five years.

Conversion of DN Values to Physical Units
The processing of images was carried out in floating points 32-bit. These images were processed in units from irradiance to absolute radiance. Finally, 32-bit image was converted in to integer form of 16-bit instead of floating form with the original level 1 format. With the help of metadata file (having different values of different bands), images had been converted from spectral radiance to radiance as famous as physical unit. This conversion is computed from the equation mentioned below:  M ρ = Re f lectance multiplicative scaling f actor f or the band( f rom the metadata).
A ρ = Re f lectance additive scaling f actor f or the band( f rom the metadata).

DN to Brightness Conversion for Kelvin Temperature
Temperature in Kelvin degree was computed by conversion of Thermal Infrared satellite (TIRS) data. This temperature was computed using satellite prescribed emissivity from its metadata file. Equation (10) explains procedure of Kelvin temperature from the satellite images after conversion into brightness temperature. Constant value of emissivity was 0.95. Constants of below equation are mentioned in Table 2.
= Thermal conversion constant f or the band ( f rom the metadata) K2 = Thermal conversion constant f or the band ( f rom the metadata)

Results
Two types of results are discussed in successive paras, i.e., Urban Sprawl (including Land use and Land Cover) and results of LST. Figure 3a,b are showing a major increase in built-up area. In 1990, urban area was 9.07%, vegetation 39.14%, water 26.64% and bare land were 25.13%. In 1995, urban area was increased from 9.07% to 10.15%. Vegetation decreased from 39.14% to 38.39%, water 26.64% to 25.94% and bare land was increased from 25.13% to 25.49%.

Results
Two types of results are discussed in successive paras, i.e., Urban Sprawl (including Land use and Land Cover) and results of LST. Figure 3a,b are showing a major increase in built-up area. In 1990, urban area was 9.07%, vegetation 39.14%, water 26.64% and bare land were 25.13%. In 1995, urban area was increased from 9.07% to 10.15%. Vegetation decreased from 39.14% to 38.39%, water 26.64% to 25.94% and bare land was increased from 25.13% to 25.49%. Figure 4a displays results of 2000, built-up area/urban area was increased up to 11.27%. Vegetation increased from 38.39% to 39.61% (city government planned to plant more than 700,000 trees [26]), water decreased from 25.94% to 24.377% and bare land was   11.27%. Vegetation increased from 38.39% to 39.61% (city government planned to plant more than 700,000 trees [26]), water decreased from 25.94% to 24.377% and bare land was also decreased from 25.49% to 24.733%. Figure 4b is portraying a glaring variation in land cover change. In 2005, built-up area/urban area was increased from 11.27% to 14.72%. Vegetation decreased from 39.61% to 33.04%, water decreased from 24.377% to 22.75% and bare land was also increased from 24.733% to 29.469%. Figure 4a displays results of 2000, built-up area/urban area was increased up to 11.27%. Vegetation increased from 38.39% to 39.61% (city government planned to plant more than 700,000 trees [26]), water decreased from 25.94% to 24.377% and bare land was also decreased from 25.49% to 24.733%. Figure 4b is portraying a glaring variation in land cover change. In 2005, built-up area/urban area was increased from 11.27% to 14.72%. Vegetation decreased from 39.61% to 33.04%, water decreased from 24.377% to 22.75% and bare land was also increased from 24.733% to 29.469%.    Where built-up area increased from 14.72% to 16.5%. It happened because of IDPs and people migrated from Balakot, Rawalakot and erstwhile federally administrated agencies to Abbottabad city and its fringe areas as it is a safe place and well settled in terms of all basic amenities of life [27]. However, Figure 5b is illustrating a gigantically great difference in land cover change from all the figures from Figure 3 to Figure 5. In Figure 5b, urban sprawl increased from 16.5% to 24.22%. This increase is more than 275 times from 1990 that shows that rapid urban growth in the study area. This increase brought about many other changes in micro climate as well as land cover. A glaring change resulted due to built-up that followed by vegetation decreased from 30.59% to 25.39%, water decreased from 17.81% to 12.69% and bare land was also increased from 35.03% to 37.68% that will be the part of urban dwelling in future. This change is due to a settled influx in the study area which brings about the cause of continuous microlevel change in the environment. Its proof is a decrease of vegetation and water sources as well as water reservoirs that is shown in Table 3 and its graphical representation is shown at Figure 6.
climate as well as land cover. A glaring change resulted due to built-up that followed by vegetation decreased from 30.59% to 25.39%, water decreased from 17.81% to 12.69% and bare land was also increased from 35.03% to 37.68% that will be the part of urban dwelling in future. This change is due to a settled influx in the study area which brings about the cause of continuous microlevel change in the environment. Its proof is a decrease of vegetation and water sources as well as water reservoirs that is shown in Table 3 and its graphical representation is shown at Figure 6.   climate as well as land cover. A glaring change resulted due to built-up that followed by vegetation decreased from 30.59% to 25.39%, water decreased from 17.81% to 12.69% and bare land was also increased from 35.03% to 37.68% that will be the part of urban dwelling in future. This change is due to a settled influx in the study area which brings about the cause of continuous microlevel change in the environment. Its proof is a decrease of vegetation and water sources as well as water reservoirs that is shown in Table 3 and its graphical representation is shown at Figure 6.

Results of Overall Accuracy and Kappa Coefficient
Results of the accuracy are shown as per values of Table 4 (classification fields) the classes which have been calculated and checked are water, vegetation, urban or built-up and bare land. Producer's accuracy and User's accuracy are mentioned in the Table 4 a,f, however, overall accuracy of six classified maps is 91.2%, 81.25%, 90.86%, 93.75%, 92.5%, and 90.63%. Kappa Coefficient used which is one of the most common means of expressing classification accuracy [28]. The table shown here is known as error matrix, it compares on a category by category basis, the relationship between known referenced data (Ground Truthing) and the derived results. When K hat is ranging the values with >0.80; it is considered as strong accuracy, from 0.40 to 0.80 very good accuracy and when coefficient is <0.40 then it is believed poor accuracy [29]. Thus, accuracy of the images was proved as strong accuracy.

Results of Land Surface Temperature (LST)
LST means surface temperature on the ground or surface of the earth. This surface temperature is not easily measured and observed with naked eyes [30]. However, this surface temperature is very much measurable through various new technologies of satellites have thermal recording sensors [31]. Figure 7a showing temperature of 1990-1994 and Figure 7b depicts temperature of 1995-2000 for the months of Jan. City center is showing high temperature because of concrete surfaces. Minimum LST of Abbottabad from Jan 1995-2000 is shown that it was increased from 3 • C to 5.8 • C. In Figure 7a the temperature of Abbottabad has risen up to 20 • C. A glaring change is very much pertinent to perceptible interpretation held in the figure. In the map all three Union Councils have different temperature values from 25 • C to 26 • C. on a category by category basis, the relationship between known referenced data (Ground Truthing) and the derived results. When Khat is ranging the values with >0.80; it is considered as strong accuracy, from 0.40 to 0.80 very good accuracy and when coefficient is <0.40 then it is believed poor accuracy [29]. Thus, accuracy of the images was proved as strong accuracy.

Results of Land Surface Temperature (LST)
LST means surface temperature on the ground or surface of the earth. This surface temperature is not easily measured and observed with naked eyes [30]. However, this surface temperature is very much measurable through various new technologies of satellites have thermal recording sensors [31]. Figure 7a showing temperature of 1990-1994 and Figure 7b depicts temperature of 1995-2000 for the months of Jan. City center is showing high temperature because of concrete surfaces. Minimum LST of Abbottabad from Jan 1995-2000 is shown that it was increased from 3 °C to 5.8 °C. In Figure 7a   From Jan 2011-2016, a drastic change in temperature was witnessed because heavy snowfall was experienced in 2016 after almost 30 years as shown in Figure 9. This snowfall phenomenon is also directing us towards the existence of extreme climate changes from micro to macro level. At minimal level temperature decreased up to −4 • C to −5 • C. At maximum level the same decrease was witnessed more than 10 • C temperature within six years' time span. Complete belt of urban land cover including Abbottabad city, Mandian, Nawansher, Jhugiyan, Shimla hill, and Abbottabad villages experienced extreme cold temperature ranging up to −9.2 • C, after a heated spell which was being experienced since last 15 years. Figure 10 is showing all the LST maps in a single glance which is depicting drift of increase and decrease temperature for the month of Jan 1990 to Jan 2016. Small set of boundaries in the maps are showing Union Council boundaries of District Abbottabad. The change in temperature at micro level within the Union Council boundaries of District Abbottabad has been increased accumulatively from 1 • C to 7 • C and the same decreased in the month of Jan2016 up to 10 • C in upper limit and −4 • C to −5 • C in lower limit. It is quite eminent that temperature of city of 25 years has a rising trend having a direct impact on air temperature as well. From Jan 2011-2016, a drastic change in temperature was witnessed because heavy snowfall was experienced in 2016 after almost 30 years as shown in Figure 9. This snowfall phenomenon is also directing us towards the existence of extreme climate changes from micro to macro level. At minimal level temperature decreased up to −4 °C to −5 °C. At maximum level the same decrease was witnessed more than 10 °C temperature within six years' time span. Complete belt of urban land cover including Abbottabad city, Mandian, Nawansher, Jhugiyan, Shimla hill, and Abbottabad villages experienced extreme cold temperature ranging up to −9.2 °C, after a heated spell which was being experienced since last 15 years. Figure 10 is showing all the LST maps in a single glance which is depicting drift of increase and decrease temperature for the month of Jan 1990 to Jan 2016. Small set of boundaries in the maps are showing Union Council boundaries of District Abbottabad. The change in temperature at micro level within the Union Council boundaries of District Abbottabad has been increased accumulatively from 1 °C to 7 °C and the same decreased in the month of Jan2016 up to 10 °C in upper limit and −4 °C to −5 °C in lower limit. It is quite eminent that temperature of city of 25 years has a rising trend having a direct impact on air temperature as well.  From Jan 2011-2016, a drastic change in temperature was witnessed because heavy snowfall was experienced in 2016 after almost 30 years as shown in Figure 9. This snowfall phenomenon is also directing us towards the existence of extreme climate changes from micro to macro level. At minimal level temperature decreased up to −4 °C to −5 °C. At maximum level the same decrease was witnessed more than 10 °C temperature within six years' time span. Complete belt of urban land cover including Abbottabad city, Mandian, Nawansher, Jhugiyan, Shimla hill, and Abbottabad villages experienced extreme cold temperature ranging up to −9.2 °C, after a heated spell which was being experienced since last 15 years. Figure 10 is showing all the LST maps in a single glance which is depicting drift of increase and decrease temperature for the month of Jan 1990 to Jan 2016. Small set of boundaries in the maps are showing Union Council boundaries of District Abbottabad. The change in temperature at micro level within the Union Council boundaries of District Abbottabad has been increased accumulatively from 1 °C to 7 °C and the same decreased in the month of Jan2016 up to 10 °C in upper limit and −4 °C to −5 °C in lower limit. It is quite eminent that temperature of city of 25 years has a rising trend having a direct impact on air temperature as well.  From Jun 1990-1994, the temperature of Abbottabad city and its surrounding areas varied from 9.2 • C to 31 • C. Quality of satellite imagery taken for Jun 1990-1994 were comparatively good but cloud cover was 10%. However, cloud cover was removed by ENVY software using atmospheric correction, even then complete clouds could not be eliminated completely, because of cloud thickness and their location near to the surface. Change in temperature was observed in Figure 11b from 14 • C to 38 • C in Abbottabad city and surroundings. In Figure 11a,b, the temperature of Abbottabad has risen up to 9 • C. Average temperature of these months from 1990 to 2013 remained at 32.3 • C. While maximum temperature of the city center rose up to 41 • C which is abnormal for a hill station city like Abbottabad. From Jun 1990-1994, the temperature of Abbottabad city and its surrounding areas varied from 9.2 °C to 31 °C. Quality of satellite imagery taken for Jun 1990-1994 were comparatively good but cloud cover was 10%. However, cloud cover was removed by ENVY software using atmospheric correction, even then complete clouds could not be eliminated completely, because of cloud thickness and their location near to the surface. Change in temperature was observed in Figure 11b from 14 °C to 38 °C in Abbottabad city and surroundings. In Figure 11a,b, the temperature of Abbottabad has risen up to 9 °C. Average temperature of these months from 1990 to 2013 remained at 32.3 °C. While maximum temperature of the city center rose up to 41 °C which is abnormal for a hill station city like Abbottabad. Figure 12a-c indicate that average temperature of Jun from 2001 to 2016 remained 33.4 °C which shows 2 °C increased within last five years' time period. Temperature change within the city and surrounding areas observed from rise of 1 °C to 4 °C. Minimum temperature has been derived from the satellite images is 19 °C and maximum is 38 °C. Figure 13 shows all the LST maps which is depicting increase and decrease temperature for the month of Jun 1990 to Jun 2016. Temperature has been increased accumulatively from 1 °C to 8 °C in the city as well as the district and the same increased in the month of Jun 2011-2016 up to 8 °C in upper limits of the district and 12 °C in lower limit of the district. As per derived LST of Abbottabad city is within the bounds of 20 °C to 27 °C whereas average temperature of these months remained 33 °C as per the record of PMD (Pakistan Metrological Department) Kakul station as shown in the graphs ( Figures  14 and 15).     Figure 13 shows all the LST maps which is depicting increase and decrease temperature for the month of Jun 1990 to Jun 2016. Temperature has been increased accumulatively from 1 • C to 8 • C in the city as well as the district and the same increased in the month of Jun 2011-2016 up to 8 • C in upper limits of the district and 12 • C in lower limit of the district. As per derived LST of Abbottabad city is within the bounds of 20 • C to 27 • C whereas average temperature of these months remained 33 • C as per the record of PMD (Pakistan Metrological Department) Kakul station as shown in the graphs (Figures 14 and 15).

Discussion
Urban Sprawl is expanding rapidly not only in Pakistan but also in other countries of South Asia as well as Asian continent as a whole. Framework of Economic Growth (FEG) of Planning Commission of Pakistan 2011 has emphasized on the need of change of planning department for solving urban sprawl issues, to convert "lethargic" cities to innovative urban societies and hub of new commercial innovations [32]. Deterioration of urban planning is not new in Pakistan [33]. It is also available in many cities of other countries of South Asia as well other world [34]. This study indicates that urban sprawl adversely affects the microclimate, and as a result the poor community may be more vulnerable to this change. Thus, increasing of temperature may give rise to glaciers melting and glaciers related hazards may increase near such cities located in the northern areas of Pakistan. As a result, decreasing agricultural and receding livestock production, degradation of ecosystem, and biodiversity loss may hamper the urban environment.
Microclimate is determined by many factors as when sun rays strike the surface of the earth, it expands in the noon or day time and contracts in evening or at night. At noon the temperature touches its peaks and at mid night temperature decrease at minimum level. This diurnal change of temperature can only alter the weather at short span of time whereas its permanent change throughout the year and continue this time span change the climate permanently. In this scenario, Griffith Tylor's environmental determinism school of thought confirms such change of phenomenon [35]. If it happens continuously from gradual to maximized intensity in terms of time and space, it refers to global warming. Global warming occurs due to various causes and happenings including rise of pollutants in the air in the form of industrial waste, CFC gases, carbon addition from coal burning, emission of carbon-monoxide from vehicles, dust from biological weathering, and numerous other pollutant materials in the air; which absorb the heat and trap the sun radiation resultantly this heat blankets the Earth's surface that bring about the rise in surface temperature. Dataset of Abbottabad station was in mean monthly, which was converted into mean annual, and obtained only one value from twelve months of data, and then plotted. Resultantly outputs we get are explained in the succeeding paragraphs and figures. In the graphs from Figures 14 and 15, it is clearly evident that temperature in 1953 was above 24 °C, decreased up to 22 °C till 1957, then abrupt change is witnessed in the graph. Temperature fluctuation is less varied between 1953 to 1970 which is from 22.6 °C to 24.4 °C less one year which is 22 °C in 1957. After 1970, temperature decreased from

Discussion
Urban Sprawl is expanding rapidly not only in Pakistan but also in other countries of South Asia as well as Asian continent as a whole. Framework of Economic Growth (FEG) of Planning Commission of Pakistan 2011 has emphasized on the need of change of planning department for solving urban sprawl issues, to convert "lethargic" cities to innovative urban societies and hub of new commercial innovations [32]. Deterioration of urban planning is not new in Pakistan [33]. It is also available in many cities of other countries of South Asia as well other world [34]. This study indicates that urban sprawl adversely affects the microclimate, and as a result the poor community may be more vulnerable to this change. Thus, increasing of temperature may give rise to glaciers melting and glaciers related hazards may increase near such cities located in the northern areas of Pakistan. As a result, decreasing agricultural and receding livestock production, degradation of ecosystem, and biodiversity loss may hamper the urban environment.
Microclimate is determined by many factors as when sun rays strike the surface of the earth, it expands in the noon or day time and contracts in evening or at night. At noon the temperature touches its peaks and at mid night temperature decrease at minimum level. This diurnal change of temperature can only alter the weather at short span of time whereas its permanent change throughout the year and continue this time span change the climate permanently. In this scenario, Griffith Tylor's environmental determinism school of thought confirms such change of phenomenon [35]. If it happens continuously from gradual to maximized intensity in terms of time and space, it refers to global warming. Global warming occurs due to various causes and happenings including rise of pollutants in the air in the form of industrial waste, CFC gases, carbon addition from coal burning, emission of carbon-monoxide from vehicles, dust from biological weathering, and numerous other pollutant materials in the air; which absorb the heat and trap the sun radiation resultantly this heat blankets the Earth's surface that bring about the rise in surface temperature. Dataset of Abbottabad station was in mean monthly, which was converted into mean annual, and obtained only one value from twelve months of data, and then plotted. Resultantly outputs we get are explained in the succeeding paragraphs and figures. In the graphs from Figures 14 and 15, it is clearly evident that temperature in 1953 was above 24 • C, decreased up to 22 • C till 1957, then abrupt change is witnessed in the graph. Temperature fluctuation is less varied between 1953 to 1970 which is from 22.6 • C to 24.4 • C less one year which is 22 • C in 1957. After 1970, temperature decreased from 24.4 • C to 22 • C and remained stable in between 21.5 • C to 23.5 • C from 1971 to 1998 but after this temperature abruptly and drastically increased till 2000 onwards because of change in land cover, i.e., urban sprawl increased resultantly microclimate is changed in the context of standard deviation, the below figure shows that temperature deviated from 5.6 • C to 8.1 • C. This change is the major change in the whole mean temperature.
Land 2021, 10, x FOR PEER REVIEW 17 of 22 after this temperature abruptly and drastically increased till 2000 onwards because of change in land cover, i.e., urban sprawl increased resultantly microclimate is changed in the context of standard deviation, the below figure shows that temperature deviated from 5.6 °C to 8.1 °C. This change is the major change in the whole mean temperature.  In Figure 16, a comparison of land surface temperature calculated and air temperature collected from Kakul weather station for the month of Jun is made from the year 1990 to 2016. It can be seen that as LST is increasing, so does the air temperature. LST has risen from 31 °C in 1990 to 38 °C in 2016. In the same manner the air temperature also rose to 2 °C in 26 years which is startling, even as per the international standards. We can also observe that LST is comparatively higher than that of air temperature, reason being ground and especially hard constructed surface retain more heat than that of atmosphere. Urban sprawl rise in the past 16 years without any planning has become the main reason for this alarming rise. In Figure 17, a comparison of land surface temperature calculated and air temperature collected from Kakul weather station for the month of Jan is made from the year 1990 to 2016. LST has risen from 18 °C in 1990 to 21 °C in 2016. In the same manner the air temperature also rose by 1.7 °C in 26 years. We can also observe that LST is comparatively higher than that of air temperature, reason being ground and especially hard constructed surfaces like roads and buildings retain more heat than that of atmosphere. after this temperature abruptly and drastically increased till 2000 onwards because of change in land cover, i.e., urban sprawl increased resultantly microclimate is changed in the context of standard deviation, the below figure shows that temperature deviated from 5.6 °C to 8.1 °C. This change is the major change in the whole mean temperature.  In Figure 16, a comparison of land surface temperature calculated and air temperature collected from Kakul weather station for the month of Jun is made from the year 1990 to 2016. It can be seen that as LST is increasing, so does the air temperature. LST has risen from 31 °C in 1990 to 38 °C in 2016. In the same manner the air temperature also rose to 2 °C in 26 years which is startling, even as per the international standards. We can also observe that LST is comparatively higher than that of air temperature, reason being ground and especially hard constructed surface retain more heat than that of atmosphere. Urban sprawl rise in the past 16 years without any planning has become the main reason for this alarming rise. In Figure 17, a comparison of land surface temperature calculated and air temperature collected from Kakul weather station for the month of Jan is made from the year 1990 to 2016. LST has risen from 18 °C in 1990 to 21 °C in 2016. In the same manner the air temperature also rose by 1.7 °C in 26 years. We can also observe that LST is comparatively higher than that of air temperature, reason being ground and especially hard constructed surfaces like roads and buildings retain more heat than that of atmosphere. In Figure 16, a comparison of land surface temperature calculated and air temperature collected from Kakul weather station for the month of Jun is made from the year 1990 to 2016. It can be seen that as LST is increasing, so does the air temperature. LST has risen from 31 • C in 1990 to 38 • C in 2016. In the same manner the air temperature also rose to 2 • C in 26 years which is startling, even as per the international standards. We can also observe that LST is comparatively higher than that of air temperature, reason being ground and especially hard constructed surface retain more heat than that of atmosphere. Urban sprawl rise in the past 16 years without any planning has become the main reason for this alarming rise. In Figure 17, a comparison of land surface temperature calculated and air temperature collected from Kakul weather station for the month of Jan is made from the year 1990 to 2016. LST has risen from 18 • C in 1990 to 21 • C in 2016. In the same manner the air temperature also rose by 1.7 • C in 26 years. We can also observe that LST is comparatively higher than that of air temperature, reason being ground and especially hard constructed surfaces like roads and buildings retain more heat than that of atmosphere. Urban sprawl in the past 16 years without any planning has become the main reason for this alarming rise. Figure 18 gives clear comparison that microclimate change of temperature due to urban sprawl.
Land 2021, 10, x FOR PEER REVIEW 18 of 22 Urban sprawl in the past 16 years without any planning has become the main reason for this alarming rise. Figure 18 gives clear comparison that microclimate change of temperature due to urban sprawl.   Urban sprawl in the past 16 years without any planning has become the main reason for this alarming rise. Figure 18 gives clear comparison that microclimate change of temperature due to urban sprawl.   Besides Abbottabad, it is very evident in many cities of South and South East Asian countries like Hong Kong SAR, Shenzhen, Seoul, and Singapore have been devasted in urban planning as well as urban green environment due to rapid their growth and economic infrastructures [36].
This research will enlighten the researchers how urban sprawl is happening in one of the emerging cities of Pakistan. It advances our understanding regarding diurnal to yearly variations in temperature in a mountainous urban area that is located in the way of famous China-Pakistan economic corridor (CPEC). The study highlights the changing of surface temperature in an elevated terrain and reports how temperature between urban and non-urban lands affects not only microclimate but also energy transfer and overall habitability of urban areas [37]. Such unsustainable developments effect the natural process of energy transfer and heat up the urban materials to affect general environmental and micro climatic patterns. So, in this case, this study suggests the developments on environmental and social concerns.   Today, we are expanding cities without focusing on deadly consequences for nature and its cycles such as water distribution or availability, vegetation effects on its neighborhoods, ecosystems depletion, and gradual disruption of temperature balance. These planning practices can be witnessed in major south Asian countries like India and Pakistan, resulting in drastic influences on mean temperature distributions [38]. The decrease of vegetation areas has further aggravated the situation; often planners neglect smaller green spaces like linear strips along the roads or tree clusters planted nearby parks or recreational lands. Such micro level urban developments can provide local cooling benefits that affect the microclimate of the surrounding localities through vegetation and evapotranspiration processes. These vegetation and water scarcity have worse implications for modern planners to design the next master plans including vegetation areas and strips in down and uphill parts of such cities. These enhancements tend to produce effective cooling effects that can diminish adverse impacts of urban heat island (UHI) [39].
Studies like this research have been done previously which highlight the relationship between urban sprawl and water depletion in such a water potential area [40]. The planners can investigate such changes and find the major causes behind these variations. This study will also provide a practical guide to initiate advanced researches to relate urban sprawl with microclimate as well as environmental developments in the world.
The combination or integration of contemporary geographic information (GI) and remote sensing (RS) tools with urban planning in microclimatic assessment perspectives is a complex task but has promising tomorrow. These sorts of integrations offer climatic evaluation as a fragment of overall urban planning and no additional endeavors are needed that prove their cost-effectiveness and time effectiveness because planners can also assess the salient influences of urban processes down to microclimatic levels thus revealing overall urban morphology [41]. Such findings answer policy questions such as: Where temperature is higher and how can it be lessened via planting trees in clusters or strips? What impacts these changes caused uphill, especially in water-associated scarcity context? How can we improve such microclimatic adversities? Further to more advanced levels, a real-time environmental monitoring and assessment platform may be established through internet of things (IoT), locational transformation (LT), and GI intelligence technologies (GIIT), to address real-time socio-environmental and microclimatic problems. This level belongs to fourth industrial revolution (industry 4.0) which is far from developing countries of South Asia and especially Pakistan; we still tremendously lack geographic enlightenment at organizational level and major geographical corridor is yet to be established where smart urban planning may be integrated with such modern big datasets. This study leads us towards such paths after revealing major urban induced problem effecting physical and social environments. The use of such platforms (with IoT, LT, and GIIT) can unleash great urban planning possibilities with public and organizational participation for sustainable development [42]. Another perspective for our cities can be advised in two directions as suggested by many researchers that we have to change our manufacturing and consumption habits to achieve effective sustainable growth. Due to these habits, we are building more unsustainable land uses and causing deforestation in many areas of the world (especially developing countries) [43]. Natural environments are shrinking due to increasing residential as well as commercial demands. Thus, urban systems are complex and always evolving to create new opportunities and problems, so we need to transform our society according to modern dynamics by educating our citizens and planners to establish new effective ways of development.
In this study we have utilized Landsat data due to its free availability. Further to this, international researchers may adopt following to learn more techniques for establishing urban sprawl and LST relationships: a.
Usage of high-resolution remote sensing data which may give more accurately analysis. b.
Access hydropower protentional in the region which can be access using in-suit data of perennial streams which is abundantly available. However, in future research, we will use ground water storage (GWS) monitoring system such as GRACE (Gravity Recovery and Climate Experiment) mission to access the ground water variation to access the impact of climate variability on GWS. c.
Access data through machine learning techniques by developing application programming interfaces (APIs).

Conclusions
Due to enormous microclimate change in Abbottabad, there was a dire need to know about the reason of this microclimate change. It was revealed that it was happening due to drastic increase of urban sprawl, and importance of this research was magnified. On the basis of our results, our research hypothesis was supported. Moreover, this research will be an eye-opener for the city management authorities to monitor urban sprawl and mitigate microclimate change.
Microclimatic change has been witnessed due urban sprawl, deforestation, housing colonies and rural/urban settlements. Urban sprawl is increased 178 to 477 sq. km within the years 1990 to 2016. Vegetated area has decreased from 770 to 602 sq. km. Surface and inland water confined areas had also decreased from 524 to 360 sq. km within the last 26 years. The bare land has increased from 494 sq. km to 742 sq. km which is under urban settlements near the city. Increase in temperature as well as climatic drift at micro level in the urban built-up area of Abbottabad has been observed and extreme climatic temperatures are being experienced due to unplanned urban sprawl. LST in winters has risen from 17 • C to 23 • C and at times shoots up to 31 • C which is alarming. Air temperature has also increased from the 1980s to 2016 at a steady rate with the annual mean temperature rising by 3 • C to 4 • C. During summers average air temperature rose up to 33 • C from 28 • C and LST has a steadily risen from 28.4 • C to 35 • C. Time series analysis showed that overall increase of temperature was gradual even witnessed at micro level.