Climate Change Risk of Urban Growth and Land Use/Land Cover Conversion: An In-Depth Review of the Recent Research in Iran

: This research is the ﬁrst literature review of the past three decades’ studies on the effects of urban development and land use/land cover (LULC) change on Iran’s climate change. For this purpose, 67 articles were found, evaluated, and classiﬁed according to the spatial and temporal scale, case study, period, data type, climatic factor, methodology, and meteorological data. Moreover, the reviewed literature methodologies were classiﬁed according to the purpose, method, and data source. According to the spatial-scale results, national- and city-level studies had the lowest and highest numbers, respectively. Tehran was the most case studies because Tehran is Iran’s capital and the largest metropolitan city. In terms of the temporal scale, studies predicting future changes (urban development and climate change) included 5% of the total literature. Satellite images were the most applied data in the reviewed literature (58%). Overall, 79% of the studies used temperature-related factors to explain the climatic impacts of urban growth and LULC conversion. Spatial modeling with 52% publications was the most used method, while numerical modeling with 12% studies was the least used method. This review showed broad study gaps in applying numerical models, neighborhood scales, urban micro-scale parameters, and long-term projections forecasts due to rapid urban development in Iran compared to the rest of the world. Therefore, our synthesis will assist researchers in facilitating better design for future studies in Iran and similar countries.


Introduction
One significant consequence of the increasing human population and urbanization tendency is the global-scale challenge of uncontrolled urban expansion. The world's urban population was about 30% in 1950 and surpassed 50% in 2014. It is also estimated to reach 60% by 2030, accounting for five billion of the world's eight billion people [1]. Urban expansion influences local, regional, and global climate through two main mechanisms. First, changes in urban characteristics, such as increasing urban heat, affect the local temperature, atmospheric flows, and the intensity and frequency of lightning strikes. These changes alter the precipitation pattern toward torrential rains and storms during warm seasons and incur severe, costly damages. Moreover, temperature changes may reduce the urban environments' climate comfort and incur extra per capita expenditures due to increasing energy consumption (electricity and water) in warm months. Review studies show that urbanization reduces green space, increases impervious surfaces, and alters albedo compared to rural surfaces [2][3][4]. Second, changes in air pollutants' transport and dispersion affect the regional and global climate [5].
Increasing air pollutants' concentrations and their associated environmental, health, social, and economic consequences are among the most critical challenges facing many urban areas and their populations worldwide. A large body of research exists on the effects of urban physical expansion on urban climate and air quality in the literature. These studies have focused mainly on the characteristics of urban atmospheric turbulence in the planetary boundary layer (PBL), differences in meteorological variables between urban and rural areas (such as air temperature, wind speed, heat island, and precipitation), or observational study techniques [6][7][8][9][10][11]. With increasing technical accuracy since the late 2000s, most review articles have begun to explain the characteristics of micro-urban parameters derived from numerical climate models and climate prediction [12][13][14][15][16][17][18][19]. Studies were directed toward understanding the contribution of large urban areas to air pollution and dispersion and global warming, health impacts of air pollution, the impacts of small-and large-scale urban expansion on climate, and urban weather and climate modeling [5]. In other words, a scale shift can be noted from coarse to fine in the investigation of urban expansion and climate indicators. The focus of this research is on the micro-scale modeling of the effects of large urban areas on climatic elements, such as the spatio-temporal distribution of heat, dispersion of urban pollutants, urban airflow patterns, and deterioration of the urban heat balance due to the changing urban canopy, increasing or decreasing turbulent flows in the urban's atmospheric boundary layer that occur due to the urban surface heterogeneity and chemical transformation of pollutants during emission processes.
Despite the wealth of studies in this field, our knowledge is still inadequate in some critical parts of the world, such as Iran, which has faced unprecedented and, partly, uncontrolled urban expansion in recent decades. As a developing country, Iran has faced increasing population growth and urbanization. In 2016, the population of Iran was estimated to be 79 million, of which 59 million were urban residents. From 2011 to 2016, Iran's urban population increased by 5 million [20]. According to the World Bank, the Iran urbanization ratio increased from 26% in 1882 to 34% in 1961 and 770% in 2012 and has surpassed 73%, which is about double the world average.
Population growth, rapid socio-economic development, and the subsequent growth of urban areas and land use/land cover (LULC) change in Iran have resulted in local and regional climate change, increasing weather anomalies and other related consequences, such as sudden rainfall, severe storms, drought, and urban mean temperature increase [21]. Many Iranian urban areas are struggling with increasing air pollutant and dust concentrations. These problems have reduced the urban environments' climatic comfort and led to the consumption of excessive energy, water, and other resources and other environmentaleconomic challenges.
The global importance of mitigating climate change and creating compatible, flexible, and resilient urban environments as one of the Sustainable Development Goals highlights the importance of understanding the causes and consequences of urban-growth-induced climate change. This importance is even more significant in Iran as a rapidly urbanizing arid and semiarid country [22]. Hence, this review aimed to collate information about the effect of urban expansion and LULC conversion on climate change in the reviewed literature to identify their strengths and weaknesses and shape future research direction. Specifically, this review addressed the following goals: (1) identify the trend in urban growth research and the effect of conversion of LULC on climate change in Iran compared to other countries, (2) identify more climatic factors and data in these articles, (3) investigate methodological features and spatial and temporal scales in these articles, and (4) identify the essential scientific gaps in these articles. The review content also included the methodology and article classification, result interpretation, discussion, recognition of gaps, and the overall conclusion.

Materials and Methods
The main framework of this paper was to review articles published exclusively on the relationship of climate variability with urban growth and LULC conversion in Iran. To do this, important national (Magiran, Scientific Information Dataset (SID), Civilica) and international (Scopus, Web of Science, and Google Scholar) databases were evaluated from 1991 to 2020. The keywords used were "urban growth," "land use land cover," and "climate change." The reference list of the identified articles was also investigated to retrieve additional pertinent papers. We restricted the document type to "research articles," "conference papers," and "reviews" published in English and Persian languages. All studies derived from the national and international search engines (n = 103) were stored in a spreadsheet, and we were checked for whether they were the same as one another. Thirty-six papers were duplicated among databases, and therefore, they were excluded. In total, 67 articles published over the past three decades were selected and categorized in terms of spatial and temporal scales, case study, period, data type, climatic factors, methodology, and meteorological data, as displayed in Figure 1.

Results
Table S1 (in the Supplementary Materials) shows the reviewed literature on the impact of urban development and land use change in Iran according to the spatial and temporal scales, case studies, period, data type, climatic factors, and methodology.

Spatial Scale
The studies' spatial scale included local (i.e., part or the entire extent of the city), metropolitan area (containing two or more adjacent urban regions or a metropolitan), regional, and national scales. Our results showed that most studies (68%) were conducted at the local scale, followed by the metropolitan area scale with 28%. However, regional-and national-scale studies were minimal, accounting for 3% and 1% of the studies, respectively. Case studies are shown in Figure 2, indicating a clear skew in selecting case studies.
Sustainability 2022, 14, x FOR PEER REVIEW 5 of 23 regional-and national-scale studies were minimal, accounting for 3% and 1% of the studies, respectively. Case studies are shown in Figure 2, indicating a clear skew in selecting case studies.

Temporal scale
The articles were categorized into four temporal scales: past, present, past-to-present, and future. Past and present studies were published in the recent and distant past, respectively. Studies conducting a temporal analysis were classified as past-to-present. As shown in Figure 3, most of the articles fell into the past-to-present temporal scale and only 5% fell into the future temporal scale. Moreover, according to Table S1, periods in these studies were divided into five categories in which studies with only 1-year data were 32%, 2-10-year data were 16%, 11-20-year data were 12%, 21-30-year data were 15%, and long-term (>30 years) data were 23%.

Temporal Scale
The articles were categorized into four temporal scales: past, present, past-to-present, and future. Past and present studies were published in the recent and distant past, respectively. Studies conducting a temporal analysis were classified as past-to-present. As shown in Figure 3, most of the articles fell into the past-to-present temporal scale and only 5% fell into the future temporal scale. Moreover, according to Table S1, periods in these studies were divided into five categories in which studies with only 1-year data were 32%, 2-10-year data were 16%, 11-20-year data were 12%, 21-30-year data were 15%, and long-term (>30 years) data were 23%.

Data Type
The data used in the articles included meteorological data, satellite imagery, demographic data, thematic maps, and their combination. The percentage contribution of each data type in the selected studies is shown in Figure 3. Satellite imagery was the most common type of data (58%; Figure 4). Forty percent of the studies used a combination of these data, including meteorological/demographic data, meteorological/satellite data, demographic/satellite data, meteorological data/thematic maps, demographic data/thematic maps, meteorological, demographic data / thematic maps, meteorological/ satellites/demographic data, as well as all data types. past-to-present, and future. Past and present studies were published in the recent and distant past, respectively. Studies conducting a temporal analysis were classified as past-to-present. As shown in Figure 3, most of the articles fell into the past-to-present temporal scale and only 5% fell into the future temporal scale. Moreover, according to Table S1, periods in these studies were divided into five categories in which studies with only 1-year data were 32%, 2-10-year data were 16%, 11-20-year data were 12%, 21-30-year data were 15%, and long-term (>30 years) data were 23%.   The data used in the articles included meteorological data, satellite imagery, demographic data, thematic maps, and their combination. The percentage contribution of each data type in the selected studies is shown in Figure 3. Satellite imagery was the most common type of data (58%; Figure 4). Forty percent of the studies used a combination of these data, including meteorological/demographic data, meteorological/satellite data, demographic/satellite data, meteorological data/thematic maps, demographic data/thematic maps, meteorological, demographic data / thematic maps, meteorological/ satellites/demographic data, as well as all data types.

Climatic Factors
Several climatic factors were found to be affected by urban expansion and LULC conversion, including air temperature (mean, maximum, and minimum values), wind direction and velocity, precipitation, LST, and UHIs. These and some other parameters were also combined to explore this association, such as temperature/wind, temperature/vegetation, temperature/air pollutants, temperature/relative humidity, surface temperature/heat island, temperature/cloudiness, temperature/precipitation, surface temperature/air pollutants/heat island, air pollutants/wind speed/surface temperature, temperature/precipitation/upward long-wave radiation, temperature/precipitation/relative humidity/cloudiness/wind, and surface temperature/surface long-wave flux/total ozone/precipitation/carbon density. Figure 5 shows the percentage contribution of each climatic factor in the reviewed literature. Overall, 79% of the studies used temperature-related factors to explain urban growth and LULC conversion ( Figure  5).

Climatic Factors
Several climatic factors were found to be affected by urban expansion and LULC conversion, including air temperature (mean, maximum, and minimum values), wind direction and velocity, precipitation, LST, and UHIs. These and some other parameters were also combined to explore this association, such as temperature/wind, temperature/vegetation, temperature/air pollutants, temperature/relative humidity, surface temperature/heat island, temperature/cloudiness, temperature/precipitation, surface temperature/air pollutants/heat island, air pollutants/wind speed/surface temperature, temperature/precipitation/upward long-wave radiation, temperature/precipitation/relative humidity/cloudiness/wind, and surface temperature/surface long-wave flux/total ozone/ precipitation/carbon density. Figure 5 shows the percentage contribution of each climatic factor in the reviewed literature. Overall, 79% of the studies used temperature-related factors to explain urban growth and LULC conversion ( Figure 5).

Methodology
The selected studies' methodology was broadly categorized into three modeling classes: spatial, statistical, and numerical. The percentage use of each methodology used in the selected studies is shown in Figure 6, indicating that spatial modeling of urban growth and LULC conversion was the most common methodology applied in 52% of the studies.

Methodology
The selected studies' methodology was broadly categorized into three modeling classes: spatial, statistical, and numerical. The percentage use of each methodology used in the selected studies is shown in Figure 6, indicating that spatial modeling of urban growth and LULC conversion was the most common methodology applied in 52% of the studies.

Spatial modeling
Spatial models are a methodology or set of analytical techniques used to explore the spatial relationship between geographic attributes. This methodology mainly includes image classification and production of LULCs or urban areas in different periods and their comparison, retrieval of LST or heat islands as climatic factors, and evaluation of the relationship of spatial changes or the intensity of heat islands, LST, and other climatic factors with different LULCs. Of note is that this investigation included a visual comparison of spatial layers and their quantitative evaluation, analysis, and statistical comparison.
A well-known example of climate change in urban development is UHI, referring to the higher temperature of urban areas than their surroundings [23]. UHIs affect the

Methodology
The selected studies' methodology was broadly categorized into three modeling classes: spatial, statistical, and numerical. The percentage use of each methodology used in the selected studies is shown in Figure 6, indicating that spatial modeling of urban growth and LULC conversion was the most common methodology applied in 52% of the studies. 3.5.1. Spatial modeling Spatial models are a methodology or set of analytical techniques used to explore the spatial relationship between geographic attributes. This methodology mainly includes image classification and production of LULCs or urban areas in different periods and their comparison, retrieval of LST or heat islands as climatic factors, and evaluation of the relationship of spatial changes or the intensity of heat islands, LST, and other climatic factors with different LULCs. Of note is that this investigation included a visual comparison of spatial layers and their quantitative evaluation, analysis, and statistical comparison.
A well-known example of climate change in urban development is UHI, referring to the higher temperature of urban areas than their surroundings [23]. UHIs affect the

Spatial Modeling
Spatial models are a methodology or set of analytical techniques used to explore the spatial relationship between geographic attributes. This methodology mainly includes image classification and production of LULCs or urban areas in different periods and their comparison, retrieval of LST or heat islands as climatic factors, and evaluation of the relationship of spatial changes or the intensity of heat islands, LST, and other climatic factors with different LULCs. Of note is that this investigation included a visual comparison of spatial layers and their quantitative evaluation, analysis, and statistical comparison.
A well-known example of climate change in urban development is UHI, referring to the higher temperature of urban areas than their surroundings [23]. UHIs affect the temperature comfort of the urban environment and increase energy consumption by airconditioning systems and, therefore, reduce the air quality of urban areas compared to their surroundings. Hence, urban-growth-driven changes in climatic parameters also increase energy consumption that act as a significant source of heat emission into the atmosphere. The relationship between UHIs and energy consumption in densely populated cities, such as Washington, Los Angeles, and Dallas, showed that for an average temperature increase of 2 to 3 • F (1 to 1.5 • C), there would be an increase of 1.5 to 2% in the maximum load of electricity consumption [24]. There are also reports of heatwave deaths associated with UHIs [24]. Traditional UHI methods were mainly based on temperature data obtained from meteorological stations or vehicle-based surveys. However, remote sensing is a relatively new technique used as a suitable replacement for traditional methods in research on UHIs. Although the climatic data recorded by meteorological stations cover relatively long-term periods, their spatial resolution is low due to the limited number of stations with large geographical distances in between, reducing the accuracy of spatial generalization and thus the reliability of urban micro-climate studies. However, remote sensing techniques have the advantage of acquiring large-scale, high-spatial-resolution data over short temporal scales [25].
Remote sensing has provided a strong foundation for retrieving heat maps and determining the earth's surface radiative energy using large-scale and well-timed automated sensing of energy from the electromagnetic spectrum [26]. These images also have a wide diversity of applications in global research on urban land use analysis [27]. Using the images with appropriate resolutions, it is now possible to retrieve UHIs at different continental and regional scales, produce quantitative data on heterogeneous land surface characteristics, and better understand the relationship between temperature and urban and non-urban environments [28].
Deciding on the choice of thermal data and the type of sensor is a crucial step in spatial modeling because numerous different thermal sensors are now available for the study of urban thermal environments, such as coarse-resolution sensors of the Geodynamics Experimental Ocean Satellite (GEOS) with a spatial resolution of 4 km, the National Oceanic and Atmospheric Administration (NOAA), the Advanced Very-High-Resolution Radiometer (AVHRR), and the Moderate Resolution Imaging Spectroradiometer (MODIS) with a spatial resolution of 1 km and medium-resolution sensors of the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) and the Enhanced Thematic Mapper (ETM + ) with thermal band resolutions of 90 and 60 m, respectively. In the published studies, the land surface temperature was only retrieved from images of different medium-resolution satellites, including Landsat, Terra, and Aqua satellites and the Thematic Mapper (TM), ETM + , ASTER, MODIS, and the Thermal Infrared Sensor (TIRS) sensors.
The main algorithms used for retrieval of land surface temperature are (1) the singlewindow method, (2) the single-channel method, (3) the Surface Energy Balance Algorithm for Land (SEBAL), (4) the Landsat Science Office, and (5) Stephen-Boltzmann [29], of which SEBAL was found to be the most applied algorithm in the selected articles.
Classification techniques were applied to represent the details of LULCs by grouping image pixels into separate classes. There are three main types of classification [30]: (1) unsupervised classification, (2) supervised classification, and (3) object-oriented classification.
(1) Unsupervised classification: Pixels are divided into groups based on their reflective properties. These groups are called clusters. The user specifies the number of clusters and the type of bands for image classification. This type of classification is used when no training data are available.
(2) Supervised classification: The user labels candidate pixels as representatives of each LULC class to train the classification algorithm and perform the classification procedure. There are various supervised algorithms whose selection is a function of the classification purpose, the image properties, and the sample types were taken. The most common supervised algorithms are maximum likelihood and minimum distance. The maximum likelihood algorithm was also the most-applied classification algorithm in the reviewed literature, except for one study that used the K-mean method. This indicates the high accuracy of the classification results using the maximum likelihood algorithm.
(3) Object-oriented classification: This classification method is different from the other two types and is known as multi-resolution segmentation. In this type of classification, homogeneous objects are produced and classified by grouping pixels based on texture, content, and geometric shape. Several bands can be used for object-oriented classification. The user specifies sample points for each LULC after multi-resolution segmentation and performs calculations for object classification, allowing the software to classify the objects based on their similarity to the training points.
Supervised and unsupervised classifications were common in the 1970s. However, the object-oriented classification method is a newly developed and increasingly popular technique for image classification [30]. Our results showed that 90% of the studies used unsupervised and supervised classification methods individually or in combination to produce LULC maps. However, only 10% of the studies used the object-oriented approach.
In the evaluated studies, error matrices, ground observations, land use maps, and higher-spatial-resolution satellite images were used separately or in combination to evaluate the classification results. Of note is that hybrid methods, such as detailed field surveys and the construction of error matrices, provided more comprehensive assessments. Moreover, the validity of LST maps was commonly estimated using data from ground stations. Table 1 summarizes the spatial modeling articles by purpose and structure.

Statistical Analysis
Statistics refers to scientific methods used to collect, organize, summarize, classify, analyze, and interpret information [74] in various disciplines, especially in understanding the characteristics of climatic factors. The broad temporal pattern of climate change studies highlights the use of long-term data recorded at ground stations and the adoption of statistical methods that have the potential to analyze long-term data series. According to Figure 6 and Table 2, about 52% of the selected articles applied statistical analysis to examine the climate change risk of urban development and LULC changes. However, different statistical methods were used depending on the research objectives, such as trend analysis (detection of UHIs or other climatic factors and their temporal variability using station information), time series analysis, and comparison of climatic factors. The primary purposes of the articles using comparison and analysis of statistical relationships were to investigate the relationship between LULC changes and UHIs or other climatic factors, comparison of several LST retrieval methods, and analyzing their relationship with land surface changes and identification of the most influential factors affecting climatic variability. In many cases, spatial modeling and statistical analysis were used in combination.  Meteorological data were also used more frequently for statistical analysis or accurate assessment of UHI modeling results. According to the Iran Meteorological Organization, three meteorological stations record weather data in Iran, including terrestrial, coastal, and high-altitude stations. There are also different terrestrial stations: synoptic stations (complete monitoring equipment), agricultural research stations, and climatology stations. As shown in Table 2 and Figure 7, the selected studies used data from different terrestrial stations. Data from synoptic stations accounted for 59% of the studies. However, only 8% of the reviewed literature used data from air-pollution-monitoring stations that belong to the Iranian Department of Environment.

Statistical Analysis
Statistics refers to scientific methods used to collect, organize, summarize, classify, analyze, and interpret information [74] in various disciplines, especially in understanding the characteristics of climatic factors. The broad temporal pattern of climate change studies highlights the use of long-term data recorded at ground stations and the adoption of statistical methods that have the potential to analyze long-term data series. According to Figure 6 and Table 2, about 52% of the selected articles applied statistical analysis to examine the climate change risk of urban development and LULC changes. However, different statistical methods were used depending on the research objectives, such as trend analysis (detection of UHIs or other climatic factors and their temporal variability using station information), time series analysis, and comparison of climatic factors. The primary purposes of the articles using comparison and analysis of statistical relationships were to investigate the relationship between LULC changes and UHIs or other climatic factors, comparison of several LST retrieval methods, and analyzing their relationship with land surface changes and identification of the most influential factors affecting climatic variability. In many cases, spatial modeling and statistical analysis were used in combination.
Meteorological data were also used more frequently for statistical analysis or accurate assessment of UHI modeling results. According to the Iran Meteorological Organization, three meteorological stations record weather data in Iran, including terrestrial, coastal, and high-altitude stations. There are also different terrestrial stations: synoptic stations (complete monitoring equipment), agricultural research stations, and climatology stations. As shown in Table 2 and Figure 7, the selected studies used data from different terrestrial stations. Data from synoptic stations accounted for 59% of the studies. However, only 8% of the reviewed literature used data from air-pollution-monitoring stations that belong to the Iranian Department of Environment.

Numerical Modeling
Climate models are essential tools for predicting the effects of climate change and understanding the complex relationships between land, atmosphere, and ocean [19,91]. Since the development of initial models known as general atmosphere circulation models, significant improvements have been made in the predictive potency of numerical climate models. The development of computational facilities and advances in the basics of this knowledge (physics and dynamics) have largely improved their predictive clarity and accuracy. These provide good opportunities to explore the earth and atmosphere processes more accurately and realistically and, therefore, enhance weather forecasting from a few weeks to a few days at regional and larger scales. The critical challenge in large urban areas is the simulation of small-scale urban features underused in regional-scale climate models. Recent advances in simulating street and urban-driven air turbulences have introduced a new generation of modeling using effective urban parameters [92,93]. More information about the effect of urbanization on climatic conditions and the effects of LULC on climatic conditions is available in Le Treut et al. [94].
Measurement and analysis of these micro-urban characteristics allow accurate assessment and simulation of the climate of urban areas. Accordingly, several urban canopy designs and related atmospheric models are divided into three general categories. Two categories are single-layer and multi-layer urban canopy designs with simple geometry and are well connected with numerical meteorological models. The third category incorporates the properties of buildings into a variety of computer fluid dynamics (CFD) models to investigate their thermal effects [4]. Of the articles reviewed, only 12% used different numerical models. One study evaluated the impact of urban growth on regional climate change using factor analysis (FA) and analytical network process (ANP) analysis. Furthermore, seven articles performed a detailed and local-scale investigation of the impacts of urban development on the behavior of climatic variables. Among the seven studies, in two studies, the impact of urban geometry on the intensity of UHIs was investigated using the numerical equation OKE. Four papers applied the ENVI-met model to simulate the behavior of wind velocity, temperature, and UHIs caused by urban geometry, topography, the density and height of buildings, traffic load, width, and direction of passages. Among the seven studies, one study [73] combined spatial modeling, statistical analysis, and the ENVI-met model to assess and reduce the effect of UHIs on urban development at a local scale. Another article predicted the effects of climate change on the design of residential buildings with the Givoni's bioclimatic (GB) chart and the Canadian Earth System Model (CESM) [95]. Meteorological information was most obtained from synoptic stations in these articles, while only two articles used NCEP information. In Table 3, these articles are categorized by purpose and structure. Table 3. Numerical modeling articles by purpose and structure.

Meteorological Station Models Purpose Reference
Synoptic OKE Investigating the effect of urban geometry on the intensity of UHIs [88] Synoptic ENVI-met Simulating the behavior of wind velocity caused by urban geometry [96] Synoptic ENVI-met Analyzing the effect of urban design on wind speed and temperature variability [97] NASA/GES, NOAA/NCEP Hybrid factor analysis (FA) and analytical network process (ANP) Evaluating the effect of urban growth on regional climate change [98] Synoptic ENVI-met Exploring reconfiguration scenarios of high-density urban neighborhoods on urban temperature [99] NCEP, IPCC Givoni's bioclimatic chart, the second-generation Canadian Earth System Model Projecting the impact of climate change on design recommendations for residential buildings [95] Synoptic ENVI-met Review and reduce the heat island effect of urban development [73] -OKE Calculating the urban heat island intensity based on urban geometry [26] Sustainability 2022, 14, 338 13 of 21

Discussion
The results yielded in climatic factors, such as temperature, UHIs, precipitation, wind direction, and velocity, agreed with those of other scientists [6,11,100,101]. According to previous studies, unplanned changes in LULC have a negative effect on the environment by increasing the LST, evaporation rates, relative humidity, and irregular rainfall, which can affect water quality and the water-energy balance, soil quality, biodiversity, and ecosystem services [102,103].
Moreover, findings of air quality parameters and pollutants were in line with those of Baklanov et al. [5] and Adimalla et al. [104] due to the focus on the impact of urban development on the concentration of contaminants. Increasing numbers of people in developing countries, plus unprecedented rates of urbanization, have had a devastating impact on surrounding land, as well as biotic and abiotic resources [105]. Several researchers contend that this pressure is behind the increasing LULC changes in these regions [105,106]. Moreover, as developing countries become more urbanized, air pollution poses serious health risks to urban residents. According to Marshall et al. [107], the population density in urban areas often results in increased traffic congestion and harmful gas emissions, thereby affecting residents' health. So, due to its population size, large extent, and severe air pollution problem, more accurate and detailed prediction is highly needed to be performed in Tehran [108]. According to Frank and Engelke [108], the rise in urban population density will increase transportation demand and air pollution. In addition, rapid population growth in cities leads to changes in LULC by displacing natural resources with urban built-up areas [109][110][111].
However, the subjects reviewed in this paper, including models, climatic factors, spatio-temporal scales, and the data source, were similar to those of Yan et al. [112]. Studies conducted in Iran were mostly large-scale, a finding in line with past research [7,8]. Based on literature reviews, remote sensing in conjunction with GIS technology can provide advanced, powerful synoptic tools to understand better and monitor the spatial, temporal, and spectral characteristics of LULC changes at micro-and macro-scales [113][114][115]. A trend has emerged to monitor and model spatial data to satisfy the needs of policymakers and planners for precisely and accurately achieving LULC information [116].
We found that 79% of the national studies used temperature-related factors (i.e., air temperature, LST, and UHI) to explain urban growth and LULC conversion. In accordance with our findings, most of the temperature-related factors are associated with LULC changes caused by widespread and unplanned urbanization, which decreases evapotranspiration and aggravates the soil moisture deficit [117,118]. As the earth's surface changes due to temperature changes, the LST is increasingly regarded as an extremely important factor [119,120]. In addition, we found that SEBAL is the most-applied algorithm for retrieval of the LST. The SEBAL model has been one of the most widely used over the past 20 years because of its solid physical foundation and ability to automatically correct atmospheric effects on the LST [121,122].
A limited number of studies (5%) predicted future urban development and LULC changes to investigate their contribution to future climate change, which is similar to the findings of Taha [9] and Fernando et al. [10] but contradicts Baklanov et al. [5] and Garuma [19], because the majority of studies in other parts of the world performed longterm predictions. Data on the location, spatial extent, rate, and driving factors of urban growth have always been problematic in developing countries, making it difficult to implement effective and suitable planning policies. As the dynamic spatio-temporal change process of urban development cannot be prevented, predicting countries' future growth could be a powerful tool to understand the effects of rapid urbanization [123]. During current years, scientists have developed many models that can be used to predict LULC change [123][124][125][126]. Dynamic models (e.g., cellular automata and Markov chain) seem to be superior to empirical or statistical models at predicting changes in LULC in the future [127,128]. To date, many software packages are available to model future LULC change patterns based on the past year's dataset of climatic parameters and LULC, such as the Land Change Modeler (LCM), Cellular Automata and Markov Chain (CA-MC), DINAMICA, SLEUTH, CLUE-S, and Modules of Land Use Change Evaluation (MOLUSCE) [109,129].
Despite technical advances in the quality of satellite images, image classification techniques are still the most crucial factor affecting feature extraction accuracy. In a study at the University of Arkansas, high-resolution color aerial images and high-resolution satellite imagery were used to compare various classification methods. In both cases, the object-oriented approach outperformed the other two classification methods because of the simultaneous application of spectral data and contextual information [130]. Classification techniques are applied to represent the details of LULC by grouping image pixels into separate classes. Our results showed that 90% of the national papers used unsupervised and supervised classification methods individually or in combination to produce LULC maps. Technically, supervised classification is more accurate and reliable than the unsupervised technique, but both are superior to object-oriented classification [30]. The supervised classification techniques include the following models: random forest (RF), radial basis function (RBF), multilayer perception (MLP), decision tree (DT), support vector machine (SVM), spectral angle mapper (SAM), fuzzy adaptive resonance theory-supervised predictive mapping (fuzzy ARTMAP), naive Bayes (NB), Mahalanobis distance (MD), maximum likelihood classifier (MLC), and fuzzy logic [130][131][132]. Models for unsupervised classification include the affinity propagation (AP) cluster algorithm, fuzzy c-means algorithms, K-means algorithm, and iterative self-organizing data (ISODATA) [133,134].
Our synthesis showed that in many cases, spatial modeling and statistical analysis were used combined with the temporal analysis of LULC. Statistical methods and RS and GIS techniques have been used as an efficient substitute for assessing urban expansion and LULC [135][136][137][138].
In terms of data sources, only seven studies used micro-urban parameters as input data of numerical models, which differed from the findings of past studies [5,11,13,18,19,110], showing the wide application of such data in the studies of this field. As mentioned in the methodology section, spatial modeling and statistical analysis were the most widely applied technique in the selected studies. Numerical models were also the least applied method, while the results of some studies (e.g., [12][13][14]16,17]) showed that micro-scale numerical modeling has been more interested in the relevant literature.

Research Gap and Future Directions
Studies of the complex interactions between climate, air quality, and megacities represent a relatively new and essential field of research. Many international studies have been initiated during the past decade. In particular, several major projects have been recently realized in Mexico City, European cities, and urban studies and measurement networks in Asia. They quantified feedback mechanisms linking megacity air quality, local and regional climates, and global climate change. They developed improved tools for predicting highimpact weather events and air pollution levels in megacities. While significant advances have been made, new interdisciplinary research studies are needed to increase our current understanding of the interactions between emissions, air quality, and regional and global climates in megacities. These should include studies to address fundamental and applied research; bridge the spatial and temporal scales connecting local emissions, air quality, and weather with climate and global atmospheric chemistry; and, finally, be oriented toward realizing integrated urban services [5].
However, the number of studies conducted in this field is small compared to the rapid development and the increasing number of problems. Early studies dealt with the impacts of LULC conversion on climate change. Current studies have also incorporated urban development at large-scales but lack the effects of micro-scale urban changes on climate change. Therefore, it is crucial to increase awareness in this research field and provide promising directions for future studies to identify, prevent, and minimize the destructive effects of climate change and urban development. Therefore, in proportion to the rate of population growth, urban development, various economic-environmental problems, and the need to better ensure the well-being of urban residents, it is necessary to conduct targeted, comprehensive, and appropriate studies. The enormous gap in these studies is the underuse of numerical models that take small changes into account and accurately predict future climate change. Although many coarse-scale studies have been presented in the literature, limited attempts have been made to perform analysis at the neighborhood level. The spatial distribution of case study urban areas was not uniform, so 21 studies were carried out in Tehran City. Hence, it is suggested that national-scale studies be conducted using a combination of statistical, spatial, and numerical models in order to facilitate a more comprehensive understanding of the country's situation. Numerical modeling at the neighborhood level can also lead to a more accurate understanding of future climate change and contribute to the proper design of the urban fabric. Future studies recommend studying these changes in different cities, and if studies are conducted on an urban scale (such as Tehran), the study timescale, data nature, measurement factor, and measurement method should be changed. The most significant gap and rarely addressed issue in the timescale of these studies is the prediction of future urban development and LULC changes and the subsequent climate change in comparison with each other and also in comparison with world research. Therefore, it is suggested that more attention be paid to long-term predictions of these changes in order to provide a thorough basis for planning and managing these changes.
An important consideration in future studies is the accurate measurement of precipitation, wind direction and velocity, relative humidity, air pollutants, and the combination of these factors with temperature, which are less represented in these articles. If a study examines the temperature parameter, the spatial scale, temporal scale, data type, and methodology should be upgraded. Although the data used in these studies were from various sources, there is still a significant gap in using micro-scale urban data, such as the type and height of buildings, urban fabric design, urban streets, and different population compositions, that can help more detailed analysis. The use of high-spatial-resolution sensors is suggested in spatial modeling methods. It is also recommended that various methods be used in retrieving UHIs, in addition to the SEBAL algorithm. The use of the object-oriented method is also recommended in LULC classification. If pixel-based methods, such as supervised classification, were used, it would be helpful to consider the results of other algorithms, such as the neural network and support vector machine.

Conclusions
This study is the first attempt at conducting a thorough review relating to the strengths and weaknesses of published articles on the effects of urban growth and LULC change on climate change in Iran. Studies conducted in this field are in their infancy and have mostly performed coarse-scale investigations on urban development and climate change. The main results of our national synthesis are divided into the following categories: Tehran City, with 21 articles, was the most studied case area. Due to its population size, large extent, and severe air pollution problem, more accurate and detailed prediction is highly needed to be performed in Tehran.
Proper spatial distribution of studies and modeling these effects in different urban contexts are also imperative. Identification of these changes and their consequences before urban development is a must to inform land planners and policymakers regarding the effects of their decisions, a consideration that has been neglected in urban development on Iran and, thus, has created numerous environmental problems for urban residents and also in local emissions, air quality with climate, and global atmospheric chemistry.
Our synthesis showed that 12% use of numerical models, minimal application of microurban parameters, a lack of studies at the urban neighborhood scale, and future prediction of these changes in only 5% of the articles underscore further use of numerical models at local scales to predict the influence of urban micro-scale parameters on climate change. This would be helpful in the appropriate designation of the urban fabric, the comfort of citizens, and the reduction in resource consumption (energy, water) and environmentaleconomic problems. Different methods are needed to make long-term predictions of these changes at urban, regional, and national scales because it leads to a better pre-development planning process.
Due to rapid urbanization and the high levels of air pollution in most metropolitan areas, it is suggested that future research adopt a more diverse combination of data, such as urban parameters with socio-economic-environmental parameters, as well as a broader range of climatic parameters, thus gaining better insights into the changes and their consequences. Air pollution in large urban areas, dispersion of air pollutants and their effects on public health, and global climate changes are also important topics to address in future research at different spatial scales.
Finally, it should be noted that alleviating the effects of urban development, especially its climatic effects, would be possible by increasing effective research in line with world development, which requires purposeful, accurate research based on different and diverse methods.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/10 .3390/su14010338/s1, Table S1: Reviewed studies on climate change risk of urban growth and land use/land cover conversion.

Data Availability Statement:
The data that support the findings of this study are available from the first author upon reasonable request.