IMPACT OF URBAN EVOLUTION ON LAND-USE CHANGE OF SARGODHA CITY, PAKISTAN

Similar documents
URBAN CHANGE DETECTION OF LAHORE (PAKISTAN) USING A TIME SERIES OF SATELLITE IMAGES SINCE 1972

1. Introduction. Chaithanya, V.V. 1, Binoy, B.V. 2, Vinod, T.R. 2. Publication Date: 8 April DOI: /cloud.ijarsg.

Applications of GIS and Remote Sensing for Analysis of Urban Heat Island

West meets East: Monitoring and modeling urbanization in China Land Cover-Land Use Change Program Science Team Meeting April 3, 2012

International Journal of Intellectual Advancements and Research in Engineering Computations

Landuse and Landcover change analysis in Selaiyur village, Tambaram taluk, Chennai

Overview of Remote Sensing in Natural Resources Mapping

AssessmentofUrbanHeatIslandUHIusingRemoteSensingandGIS

Urban Expansion of the City Kolkata since last 25 years using Remote Sensing

APPLICATION OF LAND CHANGE MODELER FOR PREDICTION OF FUTURE LAND USE LAND COVER A CASE STUDY OF VIJAYAWADA CITY

URBAN HEAT ISLAND IN SEOUL

Using Geographic Information Systems and Remote Sensing Technology to Analyze Land Use Change in Harbin, China from 2005 to 2015

Research Article A Quantitative Assessment of Surface Urban Heat Islands Using Satellite Multitemporal Data over Abeokuta, Nigeria

Remote Sensing and GIS Application in Change Detection Study Using Multi Temporal Satellite

Land cover/land use mapping and cha Mongolian plateau using remote sens. Title. Author(s) Bagan, Hasi; Yamagata, Yoshiki. Citation Japan.

Abstract: About the Author:

Comparison between Land Surface Temperature Retrieval Using Classification Based Emissivity and NDVI Based Emissivity

CHANGES IN VIJAYAWADA CITY BY REMOTE SENSING AND GIS

INTERNATIONAL JOURNAL OF GEOMATICS AND GEOSCIENCES Volume 2, No 1, 2011

Spatial Analysis of Natural Hazard and Climate Change Risks in Peri-Urban Expansion Areas of Dakar, Senegal

Geospatial technology for land cover analysis

PREDICTION OF FUTURE LAND USE LAND COVER CHANGES OF VIJAYAWADA CITY USING REMOTE SENSING AND GIS

Investigation of the Effect of Transportation Network on Urban Growth by Using Satellite Images and Geographic Information Systems

Monitoring Vegetation Growth of Spectrally Landsat Satellite Imagery ETM+ 7 & TM 5 for Western Region of Iraq by Using Remote Sensing Techniques.

Shalaby, A. & Gad, A.

Geospatial Information for Urban Sprawl Planning and Policies Implementation in Developing Country s NCR Region: A Study of NOIDA City, India

Urban Tree Canopy Assessment Purcellville, Virginia

INFLUENCE OF SEA SURFACE TEMPERATURE ON COASTAL URBAN AREA - CASE STUDY IN OSAKA BAY, JAPAN -

What is Climate? Understanding and predicting climatic changes are the basic goals of climatology.

IDENTIFICATION OF TRENDS IN LAND USE/LAND COVER CHANGES IN THE MOUNT CAMEROON FOREST REGION

Preparation of LULC map from GE images for GIS based Urban Hydrological Modeling

Remote Sensing and GIS Techniques for Monitoring Industrial Wastes for Baghdad City

Analysis of Changing Land Cover in Chittagong City Corporation Area (CCC) by Remote Sensing and GIS

Land Use/Land Cover Mapping in and around South Chennai Using Remote Sensing and GIS Techniques ABSTRACT

Landuse/Landcover Change Detection in Umshing- Mawkynroh of East Khasi Hills District, Meghalaya Using Spatial Information Technology

Abstract. TECHNOFAME- A Journal of Multidisciplinary Advance Research. Vol.2 No. 2, (2013) Received: Feb.2013; Accepted Oct.

Effect of land use/land cover changes on runoff in a river basin: a case study

M.C.PALIWAL. Department of Civil Engineering NATIONAL INSTITUTE OF TECHNICAL TEACHERS TRAINING & RESEARCH, BHOPAL (M.P.), INDIA

Mapping Coastal Change Using LiDAR and Multispectral Imagery

Analysis of Rural-Urban Systems and Historic Climate Data in Mexico

Urban Hydrology - A Case Study On Water Supply And Sewage Network For Madurai Region, Using Remote Sensing & GIS Techniques

9/19/2018. Dr. Hoch. RGPL 103 Global Cities: Planning and Development.

7.1 INTRODUCTION 7.2 OBJECTIVE

2/25/2019. Taking the northern and southern hemispheres together, on average the world s population lives 24 degrees from the equator.

SPATIO-TEMPORAL ANALYSIS OF URBAN POPULATION GROWTH AND DISTRIBUTION IN AURANGABAD CITY

MAPPING LAND USE/ LAND COVER OF WEST GODAVARI DISTRICT USING NDVI TECHNIQUES AND GIS Anusha. B 1, Sridhar. P 2

CORRELATION BETWEEN URBAN HEAT ISLAND EFFECT AND THE THERMAL INERTIA USING ASTER DATA IN BEIJING, CHINA

Remote sensing Based Assessment of Urban Heat Island Phenomenon in Nagpur Metropolitan Area

Quick Response Report #126 Hurricane Floyd Flood Mapping Integrating Landsat 7 TM Satellite Imagery and DEM Data

ASSESSING THEMATIC MAP USING SAMPLING TECHNIQUE

World Geography. WG.1.1 Explain Earth s grid system and be able to locate places using degrees of latitude and longitude.

Cambridge International Examinations Cambridge International General Certificate of Secondary Education

Analysis on Climate Change of Guangzhou in Nearly 65 Years

Module 3 Indicator Land Consumption Rate to Population Growth Rate

Spatiotemporal Analysis of Noida, Greater Noida and Surrounding Areas (India) Using Remote Sensing and GIS Approaches

About the Author: Abstract:

Using Cellular Automaton to Simulate Urban Expansion in Changchun, China

International Journal of Scientific & Engineering Research, Volume 6, Issue 7, July ISSN

SPATIAL CHANGE ANALISYS BASED ON NDVI VALUES USING LANDSAT DATA: CASE STUDY IN TETOVO, MACEDONIA

Urban Expansion. Urban Expansion: a global phenomenon with local causes? Stephen Sheppard Williams College

DROUGHT ASSESSMENT USING SATELLITE DERIVED METEOROLOGICAL PARAMETERS AND NDVI IN POTOHAR REGION

Seasonal Variations of the Urban Heat Island Effect:

Spatio-temporal dynamics of the urban fringe landscapes

ScienceDirect. Local Climate Zone Study for Sustainable Megacities Development by Using Improved WUDAPT Methodology A Case Study in Guangzhou

Urban Sprawl Prediction and Change Detection Analysis in and around Thiruvannamalai Town Using Remote Sensing and GIS

VISUALIZATION URBAN SPATIAL GROWTH OF DESERT CITIES FROM SATELLITE IMAGERY: A PRELIMINARY STUDY

An Investigation of Urban form in Lahore and neighboring cities

Digital Change Detection Using Remotely Sensed Data for Monitoring Green Space Destruction in Tabriz

LAND USE CHANGE WITHIN THE BASIC DESIGN IN NORTH OF BAGHDAD

1. Introduction. S.S. Patil 1, Sachidananda 1, U.B. Angadi 2, and D.K. Prabhuraj 3

SUPPORTING INFORMATION. Ecological restoration and its effects on the

Range of Opportunities

Summer air temperature distribution and thermal environment in urban areas of Japan

INVESTIGATION LAND USE CHANGES IN MEGACITY ISTANBUL BETWEEN THE YEARS BY USING DIFFERENT TYPES OF SPATIAL DATA

URBAN SPRAWL AND SPATIO TEMPORAL ANALYSIS OF HISAR CITY IN HARYANA USING REMOTE SENSING & GIS TECNOLOGY

Advanced Image Analysis in Disaster Response

URBAN PATTERN ANALYSIS -MAJOR CITIES IN INDIA.

The Study of Impact of Urbanization on Urban Heat Island with Temperature Variation Analysis of MODIS Data Using Remote Sensing and GIS Technology

Space Application in Support of Land Management for SDG Implementation

APPLICATION OF REMOTE SENSING IN LAND USE CHANGE PATTERN IN DA NANG CITY, VIETNAM

RiskCity Training package on the Application of GIS for multi- hazard risk assessment in an urban environment.

INDIANA ACADEMIC STANDARDS FOR SOCIAL STUDIES, WORLD GEOGRAPHY. PAGE(S) WHERE TAUGHT (If submission is not a book, cite appropriate location(s))

Earth Science Lesson Plan Quarter 2, Week 6, Day 1

National Disaster Management Centre (NDMC) Republic of Maldives. Location

A "Local Climate Zone" based approach to urban planning in Colombo, Sri Lanka

Mario Flores, Graduate Student Department of Applied Mathematics, UTSA. EES 5053: Remote Sensing

Measuring Urban Sprawl A Case Study of China

Display data in a map-like format so that geographic patterns and interrelationships are visible

WELCOME TO GCSE GEOGRAPHY WHERE WILL IT TAKE US TODAY?

INTERNATIONAL JOURNAL OF GEOMATICS AND GEOSCIENCES Volume 7, No 1, 2016

Planning Road Networks in New Cities Using GIS: The Case of New Sohag, Egypt

THE REVISION OF 1:50000 TOPOGRAPHIC MAP OF ONITSHA METROPOLIS, ANAMBRA STATE, NIGERIA USING NIGERIASAT-1 IMAGERY

Simulating urban growth in South Asia: A SLEUTH application

Submitted to: Central Coalfields Limited Ranchi, Jharkhand. Ashoka & Piparwar OCPs, CCL

UKRAINE CITIES IN EUROPE AND CENTRAL ASIA METHODOLOGY. Public Disclosure Authorized. Public Disclosure Authorized. Public Disclosure Authorized

Joint International Mechanical, Electronic and Information Technology Conference (JIMET 2015)

Contents. Introduction Study area Data and Methodology Results Conclusions

Population Density and Growth. Distribution of people on Earth

Urban Expansion and Loss of Agricultural Land: A Remote Sensing Based Study of Shirpur City, Maharashtra

INTERNATIONAL JOURNAL OF GEOMATICS AND GEOSCIENCES Volume 4, No 1, 2013

Transcription:

Impact of Urban Evolution on Land-Use Change of Sargodha City in Pakistan J. Sc. & Tech. Univ. Peshawar, 2015, 39 (2), 29-36. IMPACT OF URBAN EVOLUTION ON LAND-USE CHANGE OF SARGODHA CITY, PAKISTAN SAJJAD HUSSAIN SAJJAD 1, SADAF HUSSAIN 1, SAFDAR ALI SHIRAZI 2, KHADIJA SHAKRULLAH 3, RABIA BATOOL 1, BASIT SAEED 1 1 Department of Earth Sciences, University of Sargodha, Pakistan 2 Department of Geography, University of the Punjab Lahore, Pakistan 3 Department of Geography, Forman Christian College (FCC), A Chartered University, Lahore Pakistan Abstract: Sargodha is the 11th largest city of Pakistan with a population size of 0.7 million. This comparatively small size city is expanding faster. The objective of this work is to the impact of rapid urbanization on land-use change of Sargodha city. For this purpose to check out the extension of the built up area of Sargodha city due to urbanization four satellite images for the years 1980, 1992, 2000 and 2013, Arc GIS software was used for mapping and analysis through land-use classification. The results showed the regular increase in built-up are of the city while vegetation cover and bare land is gradually reducing due to expansion of the city. Keywords: Urbanization, land-use change, built-up area, Sargodha, Arc GIS. Introduction The process of urbanization is one of the most important dimensions of economic, social and physical change in developing countries, when growth rate is higher in cities than rural areas, it is called urbanization. Indifferent countries of the world if half per cent of total population live in towns, country is considered urbanized. European countries were first, which became urbanized, each state used exceptional criteria to define built up area like population size and space. Round about 97% population lived in rural areas in 1800.In 1900 only 5.5 per cent lived in urban areas and half of the total world s population lived in urban area in 2002 (Long, 1998).The urban population of developing nations was 300 million in 1950, at present which has increased more than1.7 billion. The growth in urban population is three times faster than surroundings. It is expected that the population of urban areas of the world will be twice in the next 25 years, reaching nearly up to four billion. Because of the migration of the people from one place to another and the increased birth rate. Rapid growth of human population is the cause of environmental change, deforestation and desertification. The process of urbanization produces change in land surface and atmospheric properties area, because vegetation decreased and replace into asphalt and concert (Oke, 1987). Land-use and land-cover (LULC) change is a noticeable effect of human s alteration in the global ecological unit, and it has a major impact on the local, regional atmosphere (Weng, 2001). Changing land-use of the city and rapidly increasing urbanization has large impact on environment, in urban areas use of energy increase and reduces incoming radiations (Stanhill & Kalma, 1995). In cities non porous surfaces e.g. walls, roads absorb the incoming solar radiations and change into sensible heat, which increased temperature in urban areas. Built up area due to high building density generate the heat than surrounding air temperature that is 29

Sajjad Hussain Sajjad, Sadaf Hussain, Safdar Ali Shirazi, Khadija Shakrullah, Rabia Batool, Basit Saeed caused urban heat island (Kimura & Takahashi, 1991). Because of increasing use of air conditioner rise, global warming and air problems are common in cities (McCarthy et al., 2001). When in metropolitan cities, vegetation cover decreases and it is convert into asphalt structure, theses hard surfaces hold heat for long time than rural area. It is because of it that less vegetation cover in urban areas than rural areas is the major contributing factor for the phenomenon of urban heat island (Sailor, 1994). Above the last decade, major impacts is anthropogenic heat, which change land cover and change climate variables, the contribution to surface warming was not received more attention (Zhang et al., 2005). Urbanization is a long-term form of land-use and land-cover change and its extent of increase are linked to population growth and economic development. Extension of the built-up area rapidly changed because of socioeconomic activities which have numerous impacts on urban atmosphere e.g. water and air pollution. The growth rate of urban population of Pakistan is higher than the other regional countries of South Asia. Table 1 Population growth Sargodha city 1998 2011 Future's estimated population 2016 2018 2020 531,896 692,250 789,729 825,182 862,227 Source: Urban unit, 2011 The urban population growth rate of Pakistan is 3.5 per cent while it is 2.8 per cent for other South Asia. Sargodha is one of the largest cities that is considered 11 th largest city of Pakistan (PDS, 2007). The population of the city will keep on growing (Table 1). The growing population of the city not only causing to increase the demand of new infrastructure (roads, buildings), it also has greater effect to reduce the vegetation cover, which reduces photosynthesis and changes the surface roughness, albedo and hydrology (Bounoua et al., 2009). Part of this study explains the study area and methodology, part 3 elaborates the land-use change during different years and in part 4 discussions and conclusion is given. Land use change Materials and Methods For the land-use change detection, post-classification comparison change detection method is applied. It requires rectification and classification of each remotely sensed image. A change detection matrix is used to compare the resulting maps on a pixel-by-pixel after the classification of each image separately. The following image processing procedures were employed in this study: 1) data collection; 2) data preparation; 3) supervised image classification; 4) analysis and 5) preparation of change detection maps. These applications were carried out using ERDAS imagine and Arc Map software. For this purpose, land sat satellite images for the years given in table 2 were downloaded from USGS website (http://glovis.usgs.gov/). Study area Sargodha which is located between 72 38" N and 72 43"N and between 32 3" E and 32 7"E.It is located on flat surface with an altitude of 190 meter (Topo Contour, 2015). Sargodha (city of Eagles) is the 11 th most populated metropolitan city in Pakistan and the 5 th largest in Punjab (Iqbal et at., 2009) and 11 th largest city in Pakistan with population of 0.7 million (Urban Unit, 2013). Covering an area of 52 km 2 30

Impact of Urban Evolution on Land-Use Change of Sargodha City in Pakistan (Demographia, 2013), the city is connected with all the major cities of the country by road and rail. The city has hot and cold climate. The residential areas of the city were developed in a series of blocks (Punjab Development Statistics, 2011) which titled it as a plan city. Table 2 Detail dates of the satellite images considered for classification of land use change of Sargodha city. Serial No Satellite Images 1 16th May 1980 2 31st May 1992 3 25th Aug 2000 4 12th July 2013 Data Preparation For the case of land-use change detection it consists of geo referencing and sub-setting of images. Sub setting of an image is done to clip study area from a complete scene. Figure.1 shows the original image and subsets of study area that belong to four different years as a sample. Supervised Image classification Supervised Image Classification is used that is the process of sorting pixels into a finite number of individual classes, or categories of data, based on their data file values. In this process pattern recognition technique which is essential to find the patterns in data, supervised training technique which is closely controlled by the analyst and signature evaluation methods are considered. Fig.1. Original images and subset study areas of 1980, 1992, 2000 and 2013 31

Sajjad Hussain Sajjad, Sadaf Hussain, Safdar Ali Shirazi, Khadija Shakrullah, Rabia Batool, Basit Saeed The result of training was a set of signatures that defined a training sample or cluster. Each signature corresponded to a class, and was used with a decision rule to assign the pixels in the image file to a class. Based on the on the signatures by use of a classification decision rule, the pixels of the image were sorted into classes after the signatures were defined. To classify the pixels in which unknown pixels were assigned to classes using contours of probability, maximum-likelihood statistics was used. Maximum likelihood method was used to classify the pixels in which unknown pixels were assigned to classes using contours of probability around training areas using the maximum-likelihood statistic. The output file was informed of image file with having thematic raster layer. In this file the data about class name, class table, class values, their statistics and histogram was automatically noted. Data analysis Data was analyzed into the following steps: (i) (ii) Calculation of respective area of all land-use classes after the classification of all the maps into six land-use classes by using reclassification tool in ARC map. Preparation of graphs of all the landuse types of each year (given in combine in figure 5). Results and Discussion Urban land-use/land-cover change In this work Arc GIS software is use for mapping and analysis which gave us batter way of image classification into shape file to visualize in better and attractive approach than other software packages. To check out the extension of the built up area of Sargodha city due to urbanization four satellite images for the years 1980, 1992, 2000 and 2013 were obtained and classified to take out the desire results. The output of the classified maps helped to find the change pattern from 1980 to 2013 in a visual technique in which the built-up area, dense and sparse vegetation cover, bare soil and deep and shallow water bodies were found separately. Figure 2 represents the evolution in built up area resulted from urbanization in Sargodha city from 1980 to 2013. In 1980 the total built up area of the Sargodha city was only 61.86km 2 (Table 4) and it is further increased to 102.22km 2 in 1992. During 1980 to 1992, the built up area of the city increased 3.36 km 2 per years and the trend of urban expansion is observed to north and eastward. From 1992 to 2000, the urban land covered area was 122.84km 2.During this period; the urban sprawl is measured 2.57 km 2 per year. During 1992 to 2000, most of the expansion of the city is seen toward east and south of the urban center and less expansion observed to other two sides. In 2013, the urban built-up area reached to 152.07km 2. During 2000 to 2013, the expansion of the city is measured 2.24km 2 per year. During this period, most of the expansion of the city is observed toward eastern side of the city and somehow miner growth toward other sides is also measured. The Sargodha city expands faster during 1980 to 1992 as compare to other selected years. These facts shows physical areas expand frequently. Figure 3 shows the combined map to show the evolution of land-use change of study area during different years. It helped to understand the urban sprawl and its spatial trends. Table 3 given below summarizes the built-up area of Sargodha city during 1980, 1992, 2000 and 2013. Table 3 Evolution in built-up area of Sargodha city Year Built up area (Km) 1980 61.86 1992 102.22 2000 122.84 2013 152.07 32

Impact of Urban Evolution on Land-Use Change of Sargodha City in Pakistan Fig.2. Classified image of land-use/land-cover of Sargodha city in 1980, 1992, 2000 and 2013 Fig.3. Combine map of Sargodha city for 1980to 2013 to show the extension of built-up area during the study period 33

Sajjad Hussain Sajjad, Sadaf Hussain, Safdar Ali Shirazi, Khadija Shakrullah, Rabia Batool, Basit Saeed Table 4 summarizes the total landuse classes classified from landsat images. The given data in table shows the variation in land-cover changing over time starting from 1980 to 2013. It highlights that the built-up area had been increasing while vegetation cover had been decreasing (Figure 4). Due to increasing built-up area of the city, the adjacent bare soil within in and out of the city decreased. However, the area of water body did not change a lot. Figure 5 highlights the percentage of land-use change of Sargodha city during the studied period. Table 4 Land use /land-cover change in Sargodha city in 1980, 1992, 2000 and 2013 Class 1980 1992 2000 2013 1 Build-up 61.86 102.22 122.84 152.07 2 Dense Vegation 64.47 58.81 30.85 37.93 3 Sparse Vegation 148.88 152.11 164.29 149.41 4 Bare Soil 137.75 71.32 56.17 56.94 5 Deep Water 0 19.46 23.68 6.8 6 Shallow Water 6.82 16.01 21.78 16.63 Total 420 420 420 420 Fig.4. Relationship between built-up area and vegetation cover of Sargodha city for the period of 1980 to 2013. Fig.5. Percentage of land-use /land-cover change of Sargodha city in 1980, 1992, 2000 and 2013. 34

Impact of Urban Evolution on Land-Use Change of Sargodha City in Pakistan Conclusion In this study the evolution of urbanization and land-use change of the city was investigated. Urban extension examined through six classified images for the years of 1980, 1992, 2000 and 2013. The study highlighted that how the Sargodha city expanded during last four decades. The study explained that Sargodha city is rapidly urbanizing and urban sprawl toward south and east of the city is frequent and regular trends of spatial expansion of the city are observed. From 1980 to 2013 the built-up area of the city increased from 61.86 km 2 to 152.07 km 2 that is 90.21 km 2 raise during last 33 years. The most extension of the city during the selected periods observed from 1980 to 1992. During this period the dense vegetation cover of the study area reduced from 64.47 km to 37.93 km 2 and sparse vegetation area reduce from 148.88 km 2 to 149.41 km 2. During recent period the Sargodha city becomes rapidly urbanized and the natural land converts into built-up area. So this reduction of natural vegetation land leads to disturb the energy balance within the city area. This reduction of vegetation covers also cause to increase the urban temperature from its surrounding rural areas which is known as urban heat island. References Awan, M.S. & Iqbal, N. (2009). Determinants of Urban Poverty: The Case of Medium Sized City in Pakistan. PIDE Working Papers 2010:60, pp. 1-22. Bounoua, L., Safia, A., Masek, J., Peters- Lidard, C & Imhoff, M. (2009). Impact of Urban Growth on Surface Climate: A Case Study in Oran, Algeria. Journal of Applied Meteorology & Climatology, 48(2): 217-231. Demographia (2013). Demographia World Urban Areas (World Agglomerations). 9th Annual Edition March 2013. pp. 35. Kimura, F., & Takahashi, S. (1991). The effects of land-use and anthropogenic heating on the Surface temperature in the Tokyo metropolitan area: A numerical experiment. Atmospheric Environment. Part B. Urban Atmosphere, 25(2): 155-164. Long, R. (1998). Urbanization Sociology. Encyclopedia Britannica Article. McCarthy, J. K., Canziani, OF., Leary, NA., Dokken, DJ & White, KS. (2001). Intergovernmental Panel on Climate Change 2001. Impacts, Adaptation and Vulnerability. Cambridge, UK: Cambridge University Press. Oke, T. R. (1987). Boundary layer climates. 2 nd ed. London, Routledge, 435. Palomo Del Barrio, E. (1998). Analysis of the green roofs cooling potential in buildings.energy and buildings, 27(2): 179-193. Punjab Development Statistics, (PDS). (2007). Bureau of Statistics Government of Punjab. Punjab Development Statistics. (PDS) (2011). Bureau of Statistics, Government of the Punjab. Sailor, D. J. (1994). Sensitivity of coastal meteorology and air quality to urban surface characteristics: American Meteorological Society, Boston, MA (United States). Stanhill, G., & Kalma, J. D. (1995).Solar dimming and urban heating at Hong Kong. International Journal of Climatology, 15: 933-941. 35

Sajjad Hussain Sajjad, Sadaf Hussain, Safdar Ali Shirazi, Khadija Shakrullah, Rabia Batool, Basit Saeed Topo Contour (2015) Elevation of Sargodha, Pakistan Elevation Map, Topo Contour.http://www.floodmap.net/Elevat ion/elevationmap/?gi=1166000 Extracted on 25.02.2015. Urban Unit (2011).Punjab Cities Improvement Investment Program.Sargodha City Profile.The Urban Unit P &D department, Punjab 2011; pp. 1-9. Urban Unit (2013).Punjab cities Improvement investment program.sargodha City Profile.The Urban Unit P& D department of Punjab. Weng, Q. (2001). A remote sensing-gis evaluation of urban expansion and its impact on surface temperature in the Zhujiang Delta, China. International Journal of Remote Sensing, 22 (10): 1999 2014. Zhang, J., Dong, W & Fu, C. (2005).Impact of land surface degradation in northern China and southern Mongolia on regional climate. Chinese Science Bulletin, 50(1): 75-81. 36