Simulating urban growth in a megalopolitan area using a patch‐based cellular automata

Although traditional cellular automata (CA)‐based models can effectively simulate urban land‐use changes, they typically ignore the spatial evolution of urban patches, due to their use of cell‐based simulation strategies. This research proposes a new patch‐based CA model to incorporate a spatial constraint based on the growth patterns of urban patches into the conventional CA model for reducing the uncertainty of the distribution of simulated new urban patches. In this model, the growth pattern of urban patches is first estimated using a developed indicator that is based on the local variations in existing urban patches. The urban growth is then simulated by integrating the estimated growth pattern and land suitability using a pattern‐calibrated method. In this method, the pattern of new urban patches is gradually calibrated toward the dominant growth pattern through the steps of the CA model. The proposed model is applied to simulate urban growth in the Tehran megalopolitan area during 2000–2006–2012. The results from this model were compared with two common models: cell‐based CA and logistic‐patch CA. The proposed model yields a degree of patch‐level agreement that is 23.4 and 7.5% higher than those of these pre‐existing models, respectively. This reveals that the patch‐based CA model simulates actual development patterns much better than the two other models.     

Statistical analysis of surface urban heat island intensity variations: A case study of Babol city,Iran

The urban heat island is considered as one of the most important climate change phenomena in urban areas, which can result in remarkable negative effects on flora, concentration of pollutants, air quality, energy and water consumption, human health, ecological and economic impacts, and even on global warming. The variation analysis of the surface urban heat island intensity (SUHII) is important for understanding the effect of urbanization and urban planning. The objective of this study was to present a new strategy based on the Shannon’s entropy and Pearson chi-square statistic to investigate the spatial and temporal variations of the SUHII. In this study, Landsat TM, ETM+, OLI and TIRS images, MODIS products, meteorological data, topographic and population maps of the Babol city, Iran, from 1985 to 2017, and air temperature data recorded by ground recorder devices in 2017 were used. First, Single-Channel algorithm was used to estimate land surface temperature (LST), and the maximum likelihood classifier was employed to classify Landsat images. Then, based on LST maps, surface urban heat island ratio index was employed to calculate the SUHII. Further, several statistical methods, such as the degree-of-freedom, degree-of-sprawl and degree-of-goodness, were used to analyse the SUHII variation along different geographic directions and in various time periods. Finally, correlation between various parameters such as air temperature, SUHII, population variation and degree-of-goodness index values were investigated. The results indicated that the SUHII value increased by 24% in Babol over different time periods. The correlation coefficient yielded 0.82 between the values of the difference between the mean air temperature of the urban and suburbs and the SUHII values for the geographic directions. Furthermore, the correlation coefficient between the population variation and the degree-of-goodness index values reached 0.8. The results suggested that the SUHII variation of Babol city had a high degree-of-freedom, high degree-of-sprawl and negative degree-of-goodness.     

Assessment of mesoscale convective systems using IR brightness temperature in the southwest of Iran

In this research, the spatial and temporal distribution of Mesoscale Convective Systems was assessed in the southwest of Iran using Global merged satellite IR brightness temperature (acquired from Meteosat, GOES, and GMS geostationary satellites) and synoptic station data. Event days were selected using a set of storm reports and precipitation criteria. The following criteria are used to determine the days with occurrence of convective systems: (1) at least one station reported 6-h precipitation exceeding 10 mm and (2) at least three stations reported phenomena related to convection (thunderstorm, lightning, and shower). MCSs were detected based on brightness temperature, maximum areal extent, and duration thresholds (228 K, 10,000 km², and 3 h, respectively). An MCS occurrence classification system is developed based on mean sea level, 850 and 500 hPa pressure patterns.

The results indicated that the highest frequency of MCSs occurred in December and April. Assessment of MCSs spatial frequency showed that MCS occurrence is strongly correlated with topography in April and May unlike the cold months. In other words, the role of Zagros Mountains in developing MCSs varies based on the season so that its impact increases with enhancement of mean monthly temperature. In addition, the occurrence of MCSs depends closely on the configuration of the Sudan Low in the southwest of Iran.

     

A bias-free geodetic boundary value problem approach to height datum unification

A geodetic boundary value problem (GBVP) approach has been formulated which can be used for solving the problem of height datum unification. The developed technique is applied to a test area in Southwest Finland with approximate size of 1.5° × 3° and the bias of the corresponding local height datum (local geoid) with respect to the geoid is computed. For this purpose the bias-free potential difference and gravity difference observations of the test area are used and the offset (bias) of the height datum, i.e., Finnish Height Datum 2000 (N2000) fixed to Normaal Amsterdams Peil (NAP) as origin point, with respect to the geoid is computed. The results of this computation show that potential of the origin point of N2000, i.e., NAP, is (62636857.68 ± 0.5) (m²/s²) and as such is (0.191 ± 0.003) (m) under the geoid defined by W ₀ = 62636855.8 (m²/s²). As the validity test of our methodology, the test area is divided into two parts and the corresponding potential difference and gravity difference observations are introduced into our GBVP separately and the bias of height datums of the two parts are computed with respect to the geoid. Obtaining approximately the same bias values for the height datums of the two parts being part of one height datum with one origin point proves the validity of our approach. Besides, the latter test shows the capability of our methodology for patch-wise application.     

Discrimination of sand dunes and loess deposits using grain-size analysis in northeastern Iran

Identification and characterization of aeolian deposits in arid environments provide information on mechanisms of loess and sand accumulation. The objectives of this study were to (i) identify the distribution of aeolian deposits, (ii) discriminate loess and sand deposits using granulometric data, and (iii) describe the aeolian deposition in Sarakhs area, northeastern Iran. Particle size distributions of 26 surface samples were determined using a laser grain-size analyzer. Fine sand, very fine sand, and very coarse silt were dominant fractions in studied sediments, and the sum of these fractions ranged from 46.9% in loess deposits to 93.8% in sand dunes. The mean grain size (M _z ) of sand dunes ranged from 3.31 to 3.54 ?, which gradually changed to 4.09 to 5.50 ? in loess deposits. Sorting, skewness, and kurtosis ranged from 0.84 to 1.94 ?, 0.18 to 0.49, and 0.76 to 2.38, respectively. Aeolian deposits in the area resulted in the incorporation of Hariroud River system and Kopeh Dagh Mountains for aeolian particle production and accumulation. Alluvial comminution in Hariroud River is suggested the main mechanism of sand and silt production and flood plain environment the main reservoir of these particles. The mountains of Kopeh Dagh act as a barrier and play a key role for sand and loess accumulation.     

Assessing desertification by using soil indices

Desertification generally refers to land degradation in arid, semiarid, and dry semi-humid climatic zones. It involves five principal processes: vegetation degradation, water erosion, wind erosion, salinization and waterlogging, and soil crusting and compaction. The aim of this study is assessing desertification using soil criteria. For this purpose, nine indices including sodium absorption ratio (SAR), soil gypsum percentage, soil texture, the content of HCO₃ ^?, the percentage of the organic matter, electrical conductivity (EC), pH, the content of the soil sodium, and chloride were used. The soil samples were taken in the north of Zayandeh-Rood River in Isfahan province of Iran, using soil data randomly sampled in a depth of 0–20 cm. After assessing the normality of the samples using Kolmogorov-Smirnov test, indices were imported into GIS environment and interpolated with IDW and normal and discrete kriging methods for delineating soil characteristics maps based on MEDALUS model. In this model, the data were firstly changed from 100 to 200. Thus 100 and 200 are estimated as the best and worst quality, respectively. Then the final map of soil criteria has been created by geometric mean of its indicators. The results showed that the maximum area is related to the medium class of desertification and is equal to 44,746 ha. The areas of severe and very severe classes of desertification are equal to 30,949 and 351 ha, respectively. The results also revealed that the indices of the organic matter and soil gypsum percentage are the most influential indices which affect desertification phenomenon.     

Paleoenvironment of geomorphic surfaces of an alluvial fan in the eastern Isfahan,Iran, in the light of micromorphology and clay mineralogy

Alluvial fans are one of the most important landforms in geomorphological and paloenvironmental studies. The objective of this study was the application of clay mineral assemblages and micromorphological properties of the studied paleosols in the geomorphic surfaces of an alluvial fan in the eastern Isfahan as proxies for paleoenvironmental and paleoclimatic changes. Micromorphology, X-ray diffraction, and scanning electron microscopy approaches were used to study the representative pedons. The results indicated that the illuviation process in the calcareous soils of the arid regions of the eastern Isfahan was probably in response to Quaternary moist conditions. There was no significant difference between clay coating properties of the studied relict and buried paleosols. Clay mineralogical study suggested that kaolinite and illite were inherited from the parent materials, while smectite and palygorskite were formed in the soil environment. Paleoargillic horizon was characterized by smectite and calcic (especially the calcrete) horizons were dominated by palygorskite. Palygorskite was accumulated by both neoformation and illuviation processes. High clay content, high intensity of smectite peak, and activity of the illuviation process in paleoargillic horizon demonstrated the seasonality of climate (rainfall) even in the moist periods of Quaternary in Central Iran. Clay mineralogical assemblages suggested a trend of increasing environmental aridity in the study area. This study, therefore, highlighted the role of clay mineralogical investigations in arid lands’ geomorphological and paleoenvironmental researches.     

Applying an optimized proxy-based workflow for fast history matching

History matching is still one of the main challenging parts of reservoir study especially in giant brown oil fields with lots of wells. In these cases, history matching with conventional manual technique needs many runs and takes months to get a match. In this work, an innovative approach was suggested for fast history matching in a real brown field. The workflow was employed based on an optimized proxy model for history matching of a field consisting of 14 active wells with multiple responses (which are production rate and pressure data) in the south part of Iran. The main important features of the proposed algorithm were defining a proxy model which is response surface method in which 21 model parameters were incorporated based on cubic centered face method. The proxy model was then optimized by one of the most famous algorithms which is genetic algorithm. Proxy model was successfully performed using 256 samples leading into p- value of 0.531 and R ² of 0.91 dataset. As a result, the proposed workflow and algorithm showed good and acceptable results for history matching of studied real model.     

An alternative direct method towards mean dynamic topography computations

An alternative “direct method” to “mean dynamic topography” (MDT) computations using satellite altimetry-derived “mean sea surface” (MSS) and “global geopotential model” (GGM), without direct application of the geoid, is devised. The developed approach, which is based on derivation of an equipotential surface of the gravity field of the Earth that fits to global MSS in least squares sense, is formulated via a constrained optimization problem. The validity of our method is numerically tested by computing a global MDT model based on DNSC08 MSS model and EGM2008 GGM as input data.     

On correct application of one-step inversion of gravity data

One of the main problems on the numerical solution of integral equations is the resolution of input data. Among the integral equations used in geodesy we have the “onestep inversion” based on the first derivative of the Poisson integral, which transforms gravity values on the Earth’s surface to the gravity potential on the reference ellipsoid. In this study, it is shown that the required spatial resolution of the input gravity data on the Earth’s surface for correct one-step inversion depends on the height of the computational region, the fact that if overlooked can cause totally wrong results. Consequently the following two major questions are posed: (i) How could one know whether the spatial resolution of the input gravity data for correct one-step inversion is sufficient? (ii) What should be done if the spatial resolution is not sufficient? By studying the behaviour of the integral kernel, an algorithm is presented which enables an appropriate answer to the former question. In order to address the latter question, a method is proposed to modify the integral kernel which overcomes the adverse effect of insufficient spatial resolution of the input gravity data. Our answers, which possess the novelty of the study, are numerically verified by means of real and simulated gravity data. The numerical results approve the efficiency of the proposed method in solving the problem of insufficient spatial resolution of the input gravity data for correct one-step inversion.