Urban expansion impacts on agricultural land and thermal environment in Larkana, Pakistan

Urban area expansion and the increase in the built-up area are major threats to agriculture, the natural environment, and ecology throughout the world. The population of the city of Larkana is rapidly increasing due to natural growth and uncontrolled migration from the surrounding areas, which have resulted in a haphazard increase in the built-up area over the fertile agricultural land and affected the city’s thermal environment. This research aims to evaluate the impacts of urban area expansion on agricultural land and the thermal environment of the study area by assessing the change in land use/ land cover (LULC) from 1990 to 2020 and land surface temperature (LST). For this purpose, the Random Forest (RF) algorithm was applied, and LST was calculated in Google Earth Engine (GEE) using Landsat imageries for the years 1990, 2000, 2010, and 2020. For the RF classifier algorithm, 3500 reference random points were generated for each year, which were then randomly divided into two datasets, i.e., a training sample consisting of 70% of the points and a validation dataset (30% of the points). After LULC classification, the results were validated for 1990, 2000, 2010, and 2020, and the accuracy was 88.3%, 89%, 90.01%, and 90.8%, respectively. The study results showed that the built-up area increased from 12.31 to 43.83 km2 while the barren land in the study region decreased from 56.51 km2 to 11.62 km2. The agricultural land was 66.66 km2 in 1990, expanded to 101.38 km2, and then decreased to 79.49 km2 in 2020. The results also revealed that most of the urban expansion in the last decade (2010-2020) took place on agricultural land. The urban thermal environment also showed a gradual increase in surface temperatur,e as recorded by the LST results. The LST results revealed that the maximum LST was 33.4°C in 1990, which increased to 36.1°C in 2020, and similarly, the minimum LST was 25.1°C, which also increased to 26.6°C in 2020. This study provides useful results for stakeholders to devise better policies and plans to control further haphazard urban expansion on fertile agricultural land in the study area.


Introduction
Globally, the agriculture sector is facing several threats like water scarcity, rising temperatures, drought, and urban expansion. Urbanization is a major cause of the conversion of agricultural land into the built-up area due to rapid economic growth, new industrial zones, and migration from rural areas toward towns and cities around the world. The global population is continuously increasing, from 2.5 billion people in 1950 (Al-Thawwad, 2008;Chandio and Shirazi, 2022) to 8 billion today, according to United Nations reports. The unprecedented global population growth is the result of improvements in the health sector, including medication, personal hygiene, and nutrition, which resulted in an increase in human life expectancy. Global data on urbanization also show an increase in population, particularly in the last two decades (Chandio and Shirazi, 2022). Asia is home to 4.7 billion (60%) of the world's population. In North America, 82% of the population, in the Caribbean region 81%, in Europe 74% and in Oceania 68% are settled in urban areas due to the availability of basic human facilities (UN-DESA, 2018;Chandio and Shirazi, 2022). Developing countries have been facing rapid population growth and urban sprawl for the past few decades, and the trend is expected to continue in the future, which is also the case in Pakistan (Ashraf et al., 2022). High population growth and a lack of basic amenities in rural areas have caused rural-urban migration and land use/ land cover change (LULC) from one form to another, such as from agricultural to the built-up area Yousafzai et al., 2022). Globally, it is estimated that by 2030, 60% of people will be living in urbanized localities (Ashraf et al., 2022;Ashraf et al., 2023). Urbanization is not only the increase in population but also the expansion of residential and commercial areas, boosting the demand for basic amenities like transportation, education, water supply, and sanitation, and ultimately converting agricultural land and areas covered by natural vegetation into a built-up area (Aliyu and Amadu, 2017;Yousafzai et al., 2022). The increase in the built-up area results in a number of environmental issues, such as affecting local and regional temperatures and aggravating the problem of food security (Sun et al., 2018;Farid et al., 2022).
Urban expansion triggers the expansion of the built-up area, recreational parks, and public infrastructure, encroaching on the vast agricultural and open areas surrounding a city. The increase in the built-up environment due to these human interruptions raises the urban temperature, decreases soil infiltration due to the increase in impervious areas, and thus leads to water runoff, resulting in disasters like urban heat waves and urban flooding (Rahman et al., 2019;Zou et al., 2021). Surface energy balance and urban flooding have risen because of urbanization, which has led to naturally high permeable surfaces being replaced by artificial, low-permeable surfaces (Rizwan et al., 2008;Wang et al., 2019;Zou et al., 2021). Among other impacts of urban expansion, it has taken over fertile agricultural land and intensified the issue of food security (Khan, 2019). The target of zero hunger in the Sustainable Development Goals has become a challenge to achieve by the year 2030 (Gandharum et al., 2022).
Remote sensing imagery is an important tool for monitoring changes in LULC and urban expansion (Mehdi et al., 2021). In particular, their high temporal resolution helps to monitor the rate of change and identify the possible causative factors (Gul et al., 2022;Tariq et al., 2022). The long-term temporal availability of Landsat images makes it easier to identify past LULC and detect temporal changes (Butt et al., 2015). No one can deny the role of remote sensing in mapping and monitoring LULC changes over short timescales with greater accuracy at a low or no cost (Gul et al., 2022). Working with the Google Earth Engine (GEE) since 2010, the use of remote sensing images has become more attractive as this online platform provides global scale analysis capabilities and access to multiple catalogs of satellite images, i.e., Landsat, Sentinel, MODIS, and so on (Gandharum et al., 2022). Machine learning algorithms such as the Random Forest algorithm (RF) can be used for the classification of satellite images using GEE, and these algorithms provide high accuracy of more than 90% (Gandharum et al., 2022;E. Nyland et al., 2018;Thanh Noi and Kappas, 2017).
Like other developing countries in the world, Pakistan is also facing the problem of rapid population growth and urban expansion, and it is reported that this trend will continue in the future (Gul et al., 2022). The proportion of the urban population in Pakistan was 17.8% in 1951, which increased to 36% in 2017(GOP, 2017. This abrupt increase in urban population has led to the expansion of urban areas, which has had negative consequences such as the degradation of drinking water quality, the lowering of groundwater levels, food security issues, poor water supply, and poor sanitation, often creating other related problems such as urban flooding and urban heat island phenomena. Sindh is the most urbanized province as compared to other provinces of Pakistan (Ul Din and Mak, 2021). Due to the high level of urbanization, the land resources of Sindh province have been degraded, including the saline and waterlogged fertile soil, and its freshwater resources depleted, while the infrastructure facilities are in shambles (Panhwar, 2020). Larkana is the fourth most populous city in Sindh province. This study focuses on urban expansion in the city of Larkana and its impact on agricultural land. The impact of urban expansion on the thermal environment of the city is also evaluated. To date, no such study has been conducted in Larkana. This research analyzes the spatial and temporal changes in LULC, linking the changes in LULC with the changes in agricultural land and the urban thermal environment between 1990 and 2020.

Study area
Larkana, situated on the right bank of the Indus, is the fourth most populous city in Sindh province after Karachi, Hyderabad, and Sukkur and the fifteenth largest in Pakistan in terms of population. Larkana is one of Pakistan's historic and intermediate cities. In 1951, its population was only 34,000, which increased to 490,508 in 2017 (GOP, 2017). The city extends from 27°48' to 27°61' north latitude and 68°16' to 68°24' east longitude ( Figure 1). Geographically, the district has an even topography and is rich in agriculture, and it is known as the "Paris of Pakistan" (Britannica, 2013). During the British era, Larkana was known as the "Garden of Sindh", but lost this status due to rapid urbanization. Larkana's history may be traced back to the Moen-Jo Daro Indus Valley Civilization, which was founded around 2500 BC. Generally, the population has migrated three times from rural areas toward Larkana in the last three decades. First, when Zulfiqar Ali Bhutto was Prime Minister and declared the city the Paris of Pakistan, he established the Chandka Medical College Frontiers in Environmental Science frontiersin.org 02 in 1973. Secondly, when Benazir Bhutto was Prime Minister, she redesigned the city, establishing new educational institutions. In 1988, Larkana was given the status of a division and gained more significance as the administrative and service hub of northern Sindh. Thirdly, during the superfloods of 2010, when a large number of people migrated from the different areas of Kamber-Shahdadkot, Dadu, and Shikarpur districts to settle near the Larkana city area. On all these three occasions, human migration was attracted, which changed the territorial boundaries of Larkana: significant changes in land use and land cover were observed.

Materials and methods
In the current research, the Landsat satellite images from the years 1990, 2000, 2010, and 2020 from the United States Geological Survey (USGS) have been used. The Landsat imageries for the years 1990, 2000, and 2010 were sourced from the Landsat 5 Thematic Mapper (TM) satellite, whereas the image for 2020 was from the Landsat 8 Operational Land Image (OLI) satellite (Table 1).

Image classification using a random forest algorithm
The surface reflectance products of Landsat images from each year were used for land use/ land cover classification in Google Earth Engine (GEE). The cloud mask function was applied to get cloud-free pixels. The normalized difference built-up index (NDBI), the normalized difference water index (NDWI), and the normalized difference vegetation index (NDVI) were calculated for each Landsat image. The NDBI was calculated using the NIR and SWIR bands (for Landsat 5, the formula is (Band 5Band 4) / (Band 5 + Band 4); for Landsat 8, the formula is (Band 6 -Band 5) / (Band 6 + Band 5). Gao. (1996)  To classify the Landsat images for the years 1990, 2000, 2010, and 2020, the random reference points were generated for four major land use categories, i.e., barren, built-up, water bodies, and agricultural land. Nearly 3500 reference points were generated for each year, of which 70% were used as training samples to create land use maps using the RF algorithm, and the remaining 30% was used for validation purposes. After obtaining LULC maps for each year, an accuracy assessment was performed for each LULC map using validation reference points, and an error matrix was employed to produce overall accuracy for each LULC map.
For the assessment of change in LULC, the classification results were used to observe the trend of land use change in Larkana and analyzed in ArcGIS 10.8. The LULC maps were categorized into four classes: builtup areas, water bodies, agricultural land, and barren land. The built-up area includes residential, commercial, and other urban areas like roads, etc., while the agricultural land includes natural vegetation, cropland, and orchards. The results were then further assessed to calculate the temporal land use change using the following equation: where t1 is the initial time, t2 is the final time, A1 is the land use class cover at the initial time, and A2 is the class cover at the final time. Furthermore, the land surface temperature (LST) for each temporal image was calculated using Google Earth Engine. The  LST data were then downloaded, and the map was prepared in ArcGIS 10.8. The LST minimum and maximum were acquired from each temporal image for further interpretation. The LST was also analyzed with respect to each land use class in this study. The LST was extracted from the thermal bands of Landsat satellite imagery using the following procedure:

Digital number (DN) to spectral radiance conversion
The data received by thermal sensors were converted to radiance using the following equation: where L λ is the TOA spectral radiance (Watts/ (m2 * srad * μm); and ML is the band-specific multiplicative rescaling factor; A L is the band-specific additive rescaling factor; and Q is the quantized and calibrated standard product pixel values (DN).

Conversion to at-satellite brightness temperature
Band data can then be converted from spectral radiance to topof-atmosphere brightness temperature using the thermal constants provided in the metadata file by Eq. 3.
Tb K2 ∕ ln k1 lλ + 1 (3) In the above equation, T or Tb is the top-of-atmosphere brightness temperature (K); L λ is the TOA spectral radiance (Watts/ (m2 * srad * μm); K 1 is the band-specific thermal conversion constant available in the metadata; and K 2 is the band-specific thermal conversion constant available in the metadata file.

Land surface temperature
The LST was retrieved using a single-channel algorithm, using a single thermal band of the imagery. For this purpose, band 6 was used for Landsat 5 while band 10 was used for Landsat 8 LST calculations. The annual mean LST has been derived using Google Earth Engine. There is a 15 minutes time lag between Landsat 8 and Landsat 5 satellite equatorial crossing time (Mumtaz et al., 2023) but the mean annual value of LST is not that much affected by this time lag. The emissivity-corrected LST was derived using the Arties and Carnahan (1982) procedure.
In Eq. 4, BT = at satellite temperature, where w is the wavelength of the emitted radiance; p = h * c / s (1.438 * 10 −2 m K); h is the Planck constant (6.626 * 10 −34 J s); s is the Boltzmann constant; c is the speed of light (2.998 × 10 8 m/s); e is the emissivity.
LSE, or Land Surface Emissivity, can be defined as the ratio of the energy radiated from the surface of a material to that radiated by a black body at the same temperature and wavelength, under the same viewing conditions. It is calculated based on Eq. 5.
where Pv means vegetation proportion, derived using Eq. 6: The normalized difference vegetation index (NDVI) in this study was calculated using red and near-infrared bands using the following equation: Conversion of LST from Kelvin to Celsius The retrieved LST is in Kelvin and was converted to Celsius using Eq. 8.
All these formulae were set in GEE for both Landsat 5 and Landsat 8 imagery, and the mean annual LST was derived for each study period. The maximum, minimum, and mean LST for each LULC class were calculated using zonal statistics in ArcGIS.

Validation
The results of LULC maps derived through the RF algorithm in GEE were validated using the 30% sample points. For this purpose, the kappa coefficient technique was used in this study. The 30% validation dataset points were used to compare the original land use type to the classified image. For this purpose, the value of the classified image assigned to each point was extracted using ArcGIS, and then the original image and Google Earth imagery were used as reference data. An error matrix table was generated to quantify and summarize the total accuracy of the classified images. The results showed good performance of the RF algorithm for LULC classification. The overall accuracy in 1990 was 88.3%, 89% in 2000, 90.01% in 2010, and 90.8% in 2020.

Land use/ land cover analysis
The multi-temporal satellite-based LULC mapping for the years 1990, 2000, 2010, and 2020 was categorized into barren land, built-up area, water bodies, and agricultural land. Agricultural land includes orchards and cropland. For each respective year, the Landsat satellite images were converted to color composite with bands 4,3, and 2 RGB for Landsat 5 and 5,4, and 3 for Landsat 8 satellite, using the standard false color composite band combination (Figures 2, 3). The land use for the multi-temporal years was generated from the Landsat 5 and Landsat 8 images using the RF algorithm classification method in GEE. The multi-temporal satellite-based LULC mapping results exhibited a high increase in the built-up area, which was 12.31 km 2 in 1990 and rose to 19.87 km 2 in 2000, 23.38 km 2 in 2010, and 43.83 km 2 in 2020 (Figures 2, 3, 4). In 1990, the builtup area was only 8.62%, which grew to 13.99% in 2000, 16.68% in 2010, and 30.70% of the total study area in the year 2020 ( Figure 5B). During the period 1990-2010, it took two decades to double the built-up area, while during the period 2010-2020, it doubled again in just one decade. This surge in the built-up area was the result of a number of factors, including the special Frontiers in Environmental Science frontiersin.org attention given to the city of Larkana during the People's Party government in 1988-1990, 1993-1996-2013. During 1988-1990, Larkana was upgraded to a division, a number of educational and health institutions were built, and housing colonies like OPF and Sindh Green Housing Society were initiated, which boosted the expansion of the city. During the 1993-1996 government, the city was given further attention, and the Board of Intermediate and Secondary Education Larkana was established in 1995. The city continued to grow to the southwest and northeast with the construction of new roads, a bus terminal, new hospitals, some new educational institutions, planned housing settlements, and small industrial areas. In the LULC map of 1990, most of the barren land is actually the land that has been cleared for different development projects in the city area. The expansion of the built-up area in 1990-2010 was almost all around the city center with a few patches in the surrounding areas, but in 2020 the city expanded along the main canal and to the west. The other factors contributing to the increased built-up  (Figures 2, 3), and the same was confirmed by the local inhabitants, who said that there was a massive amount of waterlogged and saline land up until 1990, which was later converted to agricultural land. In 1990 and 2000, barren land was the second largest LULC category, while from 1990 to 2020, agricultural land remained dominant in the study area, which substantially decreased in the last two decades, i.e., between 2010 and 2020 (Figures 4, 5). Agricultural land accounted for 45.65% of all land in 1990, gradually increasing to 71.02 % by 2010 and decreasing again to 55% by 2020 ( Figure 5A). From 1990 to 2010, barren land was the second largest LULC category in terms of area occupied, which Some new housing projects were also launched and developed on the western side of the city after 2010. The barren land in the study region decreased from 39.55% in 1990 to 8.14% in 2020. This is due to the conversion of the built-up area into agricultural land. The fourth category of LULC is water bodies in Larkana city. From the analysis, an increase in water bodies was observed in the study area. This increase is due to the intensification of inland fish farming ( Figure 4D). The results revealed that the agricultural land in the study area increased from 1990 until 2010 (Figures 4, 5A). In 1990, the total area covered by agricultural land amounted to 46.65% and it increased to 65.48% in 2010, to 71.02% in 2010, and it decreased again to 55.69% in 2020. According to the field observation and discussion with the local elders, the increase in agricultural land in Larkana and its surrounding areas is the result of clearing and converting the waterlogged and saline areas into agricultural land. Most of the surrounding areas in the northern and southern parts of the city neighboring were once waterlogged and saline, and were later converted to agricultural land with support from the local community and the government. After 2010, the agricultural land decreased and the major reason is the continuous increase in the built-up area; the current study also found that the majority of agricultural land was used to expand the urban area. Therefore, we can observe that the built-up area has an incrementing trend. The relative change in the built-up area between 1990 and 2000 was 5.4%, which increased to 14% between 2010 and 2020 (Figures 4, 5B).
The overall LULC area in 1990 and 2020 and their respective changes in percentage have been shown in Table 2 and Figure 7. The results indicate that agricultural land increased by 19.24% during the study period, while the increase in the built-up area was 256.00%. The barren land in the study area decreased by almost 79.44% during the study period (Figure 7). The increase in the built-up area has grabbed the barren land and agricultural land, in the study area. Most of the changes in the built-up area were observed between 2010 and 2020 and the main reason was a mass migration from the surrounding areas after the 2010 floods. A slight increase in water bodies was observed from the study results, amounting to 5.23% ( Table 2). The increase in water bodies is due to an increase in fish farming in the study region.
In this study, the change from one land use category to another was calculated (Figure 8). For this purpose, the postclassification change detection technique was applied to quantify the change in each land use category between 1990 and 2020. The change matrix technique was adopted for this purpose to analyze the LULC changes. The LULC conversion from one category to another is represented by a two-way cross-matrix, and its corresponding value was calculated from 1990 to 2020. The largest LULC change, measuring 28.71 km 2 , went from barren land to agricultural land, followed by barren land to a built-up area (18.13 km 2 ). Similarly, a 12.63 km 2 area of agricultural land was converted into a built-up area during the last three decades (Figure 8). During the same period, 3.15 km 2 of agricultural land was converted into barren land.

LULC and land surface temperature
In this study, we also used the Landsat thermal band to calculate the land surface temperature (LST) of the study area. The LST map was arranged into five different classes ( Figure 9A). It can be observed that in 1990, the temperature of the study area was generally under 28°C. There were only a few patches of temperature class ≥ 32.1°C because most of Larkana was covered with agricultural/open land and the built-up area was smaller (Figures 6, 8A). As the agricultural land decreased and the builtup area increased, the surface temperature started to rise and continued to do so until 2020. In the year 2020, more than 80% of the area within the city boundary experienced high surface temperatures because impervious surfaces absorb more heat, which is the main reason for high temperatures in cities . It can be observed that in 2000 and 2010, the patches of high surface temperatures grew, but in 2020, the high surface temperature expanded over the city border zone and the surrounding temperature also remained high ( Figure 9A).
The Figure 9B revealed that the research area's mean minimum and mean maximum temperatures grew steadily over the course of the investigation. In 1990, the minimum temperature was 25.1°C,

Frontiers in Environmental Science
frontiersin.org while the maximum temperature was 33.4°C. By 2020, those numbers eventually grew to 26.6°C and 36.1°C, respectively ( Figure 9B). Therefore, we further evaluated and compared the built-up area with the mean temperature of the urban area, and we found that as the urban area increased, so did the mean surface temperature of the urban land ( Figure 8C). The mean temperature in 1990, when urban land was only 8.62% of the total area, was 29.25°C; however, as the urban area expanded, so did Larkana and the surrounding areas' mean surface temperature, reaching a maximum of 36.1°C, and 31.35°C when the urban area increased to 30.70% in 2020 ( Figure 9C). We have extended our analysis and tried to find out the difference in the mean LST of the built-up area with barren land, agricultural land, and water bodies. It was observed that barren land has a higher surface temperature than the built-up area ( Figure 10). However, the surface temperature of built-up land was higher than that of agricultural land and water bodies, and it increased consistently with time ( Figure 10). In Table 3, the minimum, maximum, and mean LST and their respective standard deviations were compared for each LULC category in each year. The results of LST show a continuous rise in minimum, maximum, and mean LST for each LULC category from 1990 to 2020.
The lowest recorded mean average temperature in Larkana was 15°C in January, with a maximum temperature of 22.7°C and a minimum temperature of 7.6°C in the same month ( Figure 10).  It gradually started to increase and reached its highest level in May and June. The highest minimum, maximum, and mean temperatures in the study area were observed during May and June, while they started to decrease in July and August. The monsoon rains occur in the study area during July and August, which decrease the temperature to an optimal level. The annual range of mean temperature in the study area is almost 21 o C ( Figure 11). The change in LST during the study period in different land use categories was calculated as shown in Table 4. The results showed an increase in LST in all land use categories from 1990 to 2020. The highest growth in mean and maximum LST was observed in the built-up area category, while the lowest increase in mean LST was observed in the water bodies category, and the lowest increment in maximum LST was observed in the barren land category (Table 4). Similarly, the lowest increase in minimum LST was observed in the water bodies category, and the highest was in the built-up area category. The highest increase in minimum LST in the barren land category was observed between 2000 and 2010, as was the highest increase in mean LST. Similar results were revealed for other land use categories, and the highest minimum and mean LST were observed between 2000 and 2010. The main reason behind this high LST might be the highest temperature in the year 2010 during May and June in Sindh, which caused a severe heat wave in different cities of Pakistan and also created an unusually low-pressure area that   Frontiers in Environmental Science frontiersin.org 10 attracted the monsoonal winds and became the cause of massive flooding (Rahman and Khan, 2013). In 2010, the highest temperature ever recorded in South Asia was 53.5 o C on 26 May at Mohenju-Daro, near the city of Larkana (Rahman and Khan, 2013). Agricultural land also showed the highest LST during 2000-2010, indicating that crops were also under high stress during this period. Similar results were observed for water bodies.
The annual mean air temperature obtained from the Larkana meteorological station for each particular year was compared with the annual mean land surface temperature derived from Landsat imagery. Both the annual mean air temperature and the mean annual LST show a gradual increase in temperature in the study area ( Figure 12). In 1990, the mean annual air temperature was 26.87°C, which increased to 28.05°C in 2020, while the mean annual LST was 29.25°C in 1990, which rose to 31.35°C in 2020 ( Figure 12).

FIGURE 11
Monthly mean, maximum, and minimum temperatures in Larkana.
Frontiers in Environmental Science frontiersin.org

Discussion
The current study analyzed the impacts of urban expansion on agricultural land and the urban thermal environment in Larkana. Larkana is the fourth largest city in Sindh after Karachi, Hyderbad, and Sukkur. For LULC classification, the Random Forest algorithm was applied in the Google Earth Engine online platform for the years 1990, 2000, 2010, and 2020, while the land surface temperature was also calculated in GEE using a single channel algorithm. In this study, the most pronounced change was the abrupt increase in the built-up area and the drastic decline in the open barren land area, and to some extent, a decrease in agricultural land. In the last three decades, such abrupt changes have been more likely due to population growth in Larkana. According to government records, in 1951, the Larkana city population was merely 34,000 persons, which increased to 490,508 persons as of the 2017 census (GOP, 2017). Such abrupt changes in the built-up area created a number of environmental issues and replaced other LULC categories like agricultural and open spaces (barren land) (Aburas et al., 2021;Yousafzai et al., 2022). No such study has been conducted in Larkana or even other intermediate cities in the Sindh province to date; therefore, the results of this study have been discussed in light of other studies conducted worldwide and in other cities of Pakistan. From the analysis, it was observed that an increase in the built-up area caused a decrease in the proportion of barren land or open areas in addition to agricultural land, and the same results are consistent with other studies conducted in different parts of the world and different cities of Pakistan, such as the Nile delta in Egypt (Radwan et al., 2019), the Tarai region in Nepal (Rimal et al., 2018), Kansas City, USA (Zubair et al., 2019), Torres Vedras, Portugal (Gomes et al., 2019), Kolkata, India (Sahana et al., 2018), Karachi (Arshad et al., 2020), Peshawar Yousafzai et al., 2022), Faisalabad (Tariq and Shu, 2020), Islamabad (Sohail et al., 2019), Hyderabad (Peerzado et al., 2019), Quetta (Khan et al., 2020), and Lahore (Fahad et al., 2021). The change in the built-up area was more pronounced between 2010 and 2020, and according to the local community, this abrupt increase in the built-up area is due to large migration after the 2010 floods  in the Kamber-Shahdadkot, Dadu, and Shikarpur districts. This migration not only increased the built-up area but also caused a sharp decline in the agricultural land between 2010 and 2020 (Figures 3-5). The increase in urban population and urban expansion in Larkana and its surroundings was further discussed with local elders, as no satisfactory literature was available from which we could determine the possible causes.
According to the elders of the area, Larkana is the home of versatile politicians like the late Zulfiqar Ali Bhutto and Benazir Bhutto, and during their time of government they improved the city's infrastructure in many areas, including health, education, and housing. As a result, people from the districts moved in search of better education, health, quality of life, and jobs. This overwhelming migration resulted in an unplanned expansion of the built-up area over agricultural land and encroached on open areas. The inner-city areas became congested due to the low availability of land for new development, resulting in other environmental and social issues. The LST for each category of LULC was also evaluated in this study, and it was observed that barren land and built-up area have the highest LST and similar results have been observed in other six urban centers in Pakistan (Dilawar et al., 2021). The mean, maximum, and minimum LST of all categories of LULC recorded a consistent increase, and a similar increase was also observed in other studies (Ibrahim and Rasul, 2017;Dilawar et al., 2021;Moazzam et al., 2022). The results revealed that between 1990 and 2020, the minimum LST increased by 1.52 o C, the maximum LST by 2.17 o C and the mean LST increased by 1.4 o C in the barren land category, while in the built-up area the minimum LST increase was 1.67°C, the maximum LST was 2.91 o C, and the mean LST increase was 2.22 o C. Similarly, during the study period, the increase in minimum LST in agricultural land was 1.51°C, in maximum LST it was 2.18 o C and in mean LST was 1.59°C. A comparable increase was observed in the land surface temperature of water bodies during the study period. The land surface temperature of the inner city and its surroundings rose due to increase in the impervious surfaces.

Conclusion
This study assessed the urban expansion impacts on agricultural land and the urban thermal environment in Larkana by performing a Random Forest algorithm for land use/ land cover classification in Google Earth Engine and calculating land surface temperature using Landsat images. The Random Forest algorithm is a useful and highly accurate technique for assessing temporal LULC change. The results revealed that the urban area expanded over barren land and agricultural land in Larkana. In the last three decades , the built-up area of the city increased from 12.31 km 2 to 43.83 km 2 , which is more than a threefold increase. This surge in the built-up area took over the surrounding open land, which was either barren land or agricultural land. This haphazard expansion in the urban built environment not only displaced agricultural land but also caused different environmental hazards like increasing land surface temperature and affecting the urban thermal environment. Barren land in the study area decreased from 56.51 km 2 to 11.62 km 2 , representing a decrease of 79.44% over three decades. The water bodies increased in the study region due to the increasing importance of inland fish farming. The LULC classification results were validated using 30% sample points data, and for all selected years, the accuracy was close to 90%, which indicates the importance and accuracy of the RF algorithm for temporal LULC classification. The LST maps were prepared for 1990, 2000, 2010, and 2020 using GEE in order to check the status of the urban thermal environment, and the results indicated a gradual increase in surface temperature in all land use categories. The cause of this gradual increase in surface temperature is the increase in built-up area over barren, agricultural land, and open green vegetation. This study's results endorse prior planning and policy implementation to control the haphazard urban expansion over fertile agricultural land.

Data availability statement
The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.