فهرست مطالب alireza vafaeinezhad

Drought is one of the most important natural disasters with devastating and harmful effects in various economic, social, and environmental fields. Due to the repetitive behavior of this phenomenon, if the appropriate solutions are not implemented, its destructive effects can remain in the region for years after its occurrence. Most natural disasters, such as floods, earthquakes, hurricanes, and landslides in the short term, can cause severe financial and human damage to society, but droughts are slow-moving and creepy in nature, and their devastating effects appear gradually and over a longer period of time. Therefore, by modeling drought, it is possible to provide plans for drought preparation and reduce the damage caused by it. In this study, computational intelligence algorithms of Multi-Layer Perceptron neural network, Generalized Regression Neural Network, Support Vector Regression with support kernel, and Support Vector regression with the proposed kernel (Support Vector) Regression New kernel has been used to model the drought using the Standardized Precipitation Index. The modeling results, in most cases, showed better performance of the proposed SVR_N model than other models. The values of RMSE and R2 were 0.093 and 0.991, respectively, and the GRNN, MLP, and SVR models performed better in modeling after SVR_N, respectively. Modeling of drought phenomenon in modeling is supported by vector regression method.

A drought crisis is a dry period of climate that can occur anywhere globally and with any climate. Although this crisis starts slowly, it can have a serious impact on health, agricultural products, the economy, energy, and the environment for a long time to come. Drought severely threatens human livelihood and health and increases the risk of various diseases. Therefore, modeling and predicting drought is one of the most important and serious issues in the scientific community. In the past, mathematical and statistical models such as simple regression, Auto-regression (AR), moving average (MA), and ARIMA were used to model the drought. In recent years, machine learning methods and computational intelligence to model and predict drought have been of great interest to scientists. Computational intelligence algorithms that have been previously considered by scientists to model drought include multilayer perceptron neural network, RBF neural network, support vector machine, fuzzy, and ANFIS methods. In this research, the purpose of modeling and predicting drought is by using three neural network algorithms, including multilayer perceptron, RBF neural network, and generalized regression neural. The drought index used in this research is the standardized precipitation index (SPI). In this research, the wavelet technique in combination with artificial neural network algorithms for modeling and predicting drought in 10 synoptic stations in Iran (Abadan, Babolsar, Bandar Abbas, Kerman, Mashhad, Rasht, Saqez, Tehran, Tabriz, and Zahedan) have been used in different climates and with suitable spatial distribution throughout Iran.

This study, initially using monthly precipitation data between 1961 and 2017, SPI drought index in time scales of 3, 6, 12, 18, 24, and 48 months through programming in soft environment MATLAB software implemented. The results of this step were validated using the available scientific software MDM and Drinc. Then, prediction models were designed using the Markov chain. In this study, a total of six computational intelligence models, including three single models of multilayer perceptron neural network (MLP), radial basis function neural network (RBF), and generalized regression neural network (GRNN), and three hybrids wavelet models with these three models (WMLP-WRBF-WGRNN) have been used to model and predict the SPI index in 10 stations of this research. In implementing all these six models, the MATLAB software programming environment has been used. In this study, four types of discrete wavelets were used, including Daubechies, Symlets, Coiflets, and Biorthogonal. Due to the better performance of the Dobbies wavelet, this type of wavelet was used as a final option in the research. In the Daubechies wavelet used between levels 1 to 45, level 3 showed the best performance among different SPI time scales; therefore, the Daubechies level 3 wavelet was used in all hybrid models of this study. After training all six algorithms used, the evaluation criteria of coefficient of determination (R2) and root mean square error (RMSE) was used to measure the difference between actual and estimated values.

The results of this study showed that computational intelligence methods have high accuracy in modeling and predicting the SPI drought index. In the first stage, the results showed that the individual MLP, RBF, and GRNN models, if properly trained, have close results in modeling and predicting the SPI drought index. In the next step, it was observed that the wavelet technique would improve the modeling results. In using the wavelet technique in combination with three single models MLP, RBF, and GRNN, the choice of wavelet type is also more effective in modeling, so in this research, the first of the four types of discrete wavelets Daubechies, Symlet, Qoiflet, and Biorthogonal in combination with Three single models of this research were used and the results of these four types of wavelets showed the relative superiority of the Daubechies wavelet over the other three wavelets. In using the Daubechies wavelet, since this wavelet has 45 times and the choice of order was also effective in modeling, it was observed by testing the wavelet 45 times that the 3rd wavelet, in general, has higher accuracy in all time scales of SPI index, 3, 6, 12, 18, 24 and 48 months and also in all three algorithms MLP, RBF, and GRNN. Therefore, in this research, the third-order Daubechies wavelet was used in all three algorithms of this research, as well as in all time scales. The results showed that combining the wavelet technique with all three models MLP, RBF, and GRNN will improve the results. The research graphs showed that for the quarterly time scale, the values obtained from the single model prediction in MLP and RBF modeling have a somewhat one-month phase difference compared to the hybrid model, while in the GRNN model, this prediction difference is negligible. The modeling results for both single and hybrid modeling modes indicate that there is no phase difference between the single and hybrid modeling methods in time scales of 6, 12, 18, 24, and 48. For the 12- and 24-month time scales, the single GRNN model had more fluctuations and errors in SPI monthly modeling and forecasting, while the hybrid model in these two-time scales had much better behavior in monthly modeling and forecasting. Distribution diagrams of data related to observational SPI of Abadan station showed that the modeling results for single and hybrid modes in 3 and 6-month time scales are less accurate than other time scales and fit line separation, and its uncertainty is higher than others. However, in all neural network models and in all time scales, the hybrid method has shown more accuracy. The numerical results of the study indicate that in all SPIs and stations under study, the differential values of R2 are positive, which indicates higher values of R2 in the hybrid model than in single neural network modeling, which indicates an improvement in hybrid modeling compared to individual models. Also, the differential values of RMSE are negative in all studied models and stations, which indicates that the amount of RMSE in predicting hybrid models is lower than individual neural network models. In the research graphs, it can be seen that the amount of differences in RMSE and R2 indicates a greater difference in time scales 3 and 6 than the time scales 12, 18, 24, and 48, which somehow goes back to the nature of the data of these time scales. The most significant improvement in R2 and RMSE is from the 3-month low to the 48-month high, respectively.

 From the findings of this study, it can be concluded that artificial neural network algorithms are efficient methods for modeling and predicting the SPI drought index. The use of wavelets in all three models of artificial neural networks will also improve the results. It can also be concluded that for better modeling of the SPI drought index, it is necessary to select the optimal wavelet type and order. From the results of this study, it can be concluded that the wavelet technique has a greater impact on the lower time scales, i.e., 3 and 6 months, than the higher scales, i.e., 24 and 48 months.

A method to reduce the absorption of solar radiation and prevent the creation of urban heat islands is to increase shade by vegetation. A shadow creating on buildings, causes houses to cool down, reduces energy consumption and costs, increases the value of houses, and creates a proper visual effect and a sense of well-being and vitality. Although economically, the amount of savings due to shade and cooling of the air for a tree during its lifetime in different climatic regions is different and depends on the type of tree, the amount of shade during the day and in different seasons of the year, but its effect on energy savings and costs are definite.  The subject of the present study is strategic planning to increase the shade coverage of trees in urban residential areas. A simple way to create plenty of shade is to plant numerous trees around buildings. However, this method is impractical in many areas that face water shortages due to its high costs. In addition, the presence of additional shadows on the rooftop of the buildings will reduce the ability to be exposed to sunlight and the potential of using solar panels to generate electricity. So the main challenge is using a method that can provide maximum shade coverage on the facade surface and minimum shadow coverage on the rooftop with a few trees in optimal locations. The issue of locating trees with the aim of optimizing shade coverage, i.e. maximizing shade coverage on facades and opening components, and minimizing shadow coverage on the rooftop, is a Non-deterministic Polynomial hard (NP-hard) problem and has no exact solution. Therefore, the 3D Geographic Information System and the Ant Colony Optimization algorithm have been used for this purpose. Previous studies have often examined the effects of tree canopy shade on a single building. But in most cities in Iran, buildings are connected together and form a building block. So, instead of a single building, a building block is examined. In addition, in most previous studies, the effect of shade coverage of a maximum of two trees on the building has been investigated; while in this study, we examine the effect of shade coverage of 15 trees on the building block. None of the studies on optimizing the shade of trees on the facade of the building has used the meta-heuristic optimization methods and its combination with GIS. In this study, a hybrid model of GIS in a three-dimensional environment and ACO is used for maximizing the shade of trees on the facade and opening components of buildings, and minimizing the shade of trees on the rooftop.

 Two types of data are required to perform the analysis; The building block properties, for example, dimensions, position, and size of the facade, rooftop, and opening components, and the tree properties (height and position). 3D GIS and ACO algorithms have been used to model tree shade coverage optimization. 3D GIS provides abilities for storing, analyzing, and creating 3D topologies, and ACO is used to summarize real-world conditions in a mathematical problem. GIS and trigonometric rules have been used to store geographical information and spatial topology. After storing the position, composition, and description information of 2D and 3D objects by topological data, Duffie and Beckman relations (2013) is used to extract the position of the shadow. Then, according to Church and Revelle, the Maximal Covering Location Problem (MCLP) is defined. For the following 3 reasons, ACO has been used for three-dimensional optimization; 1) The complex trigonometric rules in calculating the shadow coverage on buildings, 2) There is no deterministic solution for optimization problems because of nonlinear constraints including trigonometric functions, 3) The existence of continuous space around the building block that It is possible to place a tree in any position. The details of the steps are; 1) Define the set of possible locations for the tree based on the height, diameter of the canopy, and around space of the building block, 2) Use a method to place the first tree in all possible places around the building block during hot hours on certain days of the summer and calculate the maximum shade coverage on the building block based on the weight of the building components, 3) Remove the places that may be done in the tree canopy to prevent overlapping of tree canopies, 4) Repeat steps 2 and 3 to place the next trees in the possible places around the building block until the number of trees reaches the desired number of trees to create shade. Considering the infinite possible positions, a simplification step is required to limit the number of available positions. Therefore, the constant space is reduced to possible positions for locating Ni trees with two-meter spacing in the N-S and E-W directions. Further, the possible tree positions in front of the opening components are eliminated to make daylight available, have an outlook from the building, and comment through the doors. The minimum spacing of two meters between the trees and the building is set to prevent unnecessary shading on the rooftop.

MATLAB environment is used to optimize the shade coverage of trees using the ACO algorithm. For this purpose, properties of the buildings block such as length, width, height, are modeled in a struct in MATLAB. This struct has separate matrices for the north, east, south, and west views of the building block. Another matrix is also used to model the rooftop. Each element of the mentioned matrices is equal to 10× 10 cm from the surface of the building block and has a value of zero. To model the dimensions and location of doors and windows in each facade, another struct includes separate matrices for each facade is used. In these matrices, the amount of elements in the location of doors and windows is one. The characteristics of the sun in the study area are used, including azimuth and altitude of the sun on the studied days in 15-minute intervals from 9 to 15 hours.  The shadow is created on building components, by placing the tree in any of the possible locations, and movement of the sun. The elements of the matrices equivalent to the shaded building components change from zero to one. The sum of the values of the matrix elements determines the amount of shadow created by the tree on each component of the building. The sum of the point multiplication of the door/window matrix elements in the facade matrix elements determines the amount of shadow created on the doors/windows. The objective function is defined and the ACO algorithm is used to maximize the shadow coverage of trees on the facade, doors/windows, and minimize the shadow coverage on the rooftop. The results of the ACO show that the optimal shade coverage on the buildings block, which creates the most shade on the facade and doors and windows and the least shade on the roof, depends on the number of trees and the position of the doors and windows in buildings block. In general, as the number of trees increases, the amount of shadow created on the building block components increases.

The results of the ACO showed that for buildings, in the northern hemisphere, the trees in the north of the buildings have no effect on casting shadows on the components of the building. Due to the fact that in arid and tropical regions there are restrictions on planting trees, finding a suitable position for trees plays an important role in optimizing the shade coverage. Due to the high heat transfer through the doors and windows compared to the facade and rooftop, the higher weight is considered for these components in the objective function. Finding the optimal position of the trees depends a lot on the position of the doors and windows in the building to create the most shadow on these components. For a buildings block with the number and dimensions of buildings assumed in the research and according to the dimensions and position of doors and windows, planting a tree in one of the positions K10, K16, K22, or K28 creates the most optimal shade. These positions are 2 meters from south of the buildings and in the middle of two windows. On average, this tree provides 7.48, 9.22, and 0.85% shade respectively on the facade, doors /windows, and rooftop from 9 to 15 o'clock in four days studied. In the case of planting two trees, two positions from positions K10, K16, K22, or K28 still provide the optimal shade. On average, these two trees provide 13.88%, 18.64%, and 1.69% of shade respectively on the whole facade, doors /windows, and rooftop at 9:00 AM to 3:00 PM. In the case of three trees, positions K8, K18, and K22, in the case of four trees, positions K14, K20, K26, and K32, in the case of five trees, positions K8, K14, K20, K26, and K32 create the optimal shadow. Shading coverage in the case of three trees, is 21.07, 28.54, and 2.54%, respectively on the facade, doors/windows, and rooftop, in the case of four trees, is 24.96, 35.36 and 3.39% respectively on the façade, doors/windows, and rooftop and in the case of five trees is 33.26, 44.70 and 3.95% respectively on the facade, doors/windows, and rooftop. By planting five trees, more than 88% of the south façade and more than 90% of the south façade doors/windows of the building will be covered with shade. However, due to the goal of optimizing the shadow on the building and the greater weight of the doors and windows, the ACO has optimized the position of the trees in such a way that more surfaces of the doors and windows are exposed to the shadows. Due to the fact that in the case of five trees, 90% of the southern facade is in the shade of trees, in the case of six trees, in addition to the southern facade, the eastern and western facades are also considered for planting trees. So that the positions K8, K14, K20, and K30 are chosen in the distance of 2 meters from the south and the position of H2 is chosen in the distance f 2 meters from the west, and the position of H36 is chosen in the distance of 2 meters from the east. On average, these trees provide 33.95%, 42.29%, and 3.64% shade respectively on the facade, doors/windows, and rooftop.

One of the important issues after an earthquake occurrence is to transfer residents to temporary accommodations. These points are areas of the city that have a higher safety factor against earthquake reoccurrence and have good access to management and relief centers. Due to the critical condition in the area after the earthquake, with the transfer of the township to these areas, relief procedure can be achieved with the highest speed and minimum risk. The purpose of this study is to determine the suitable locations for the construction of temporary accommodation facilities for injured people after earthquake occurrence in Karaj. To this end, firstly, the main controlling parameters in selection of the appropriate places for temporary accommodation are determined, and are classified into two groups of the main criteria of compatibility and incompatibility. The total number of the sub-criteria are twelve. Due to the spatial nature of this problem, fuzzy logic system with multi-criteria decision-making method is used to determine the most suitable points. After fuzzification of the criteria, their participation weight is determined and then, they are prioritized using fuzzy hierarchical analysis method. In the next step, desirability map of the region is prepared for each of the criteria. Eventually, by combining the impact of the studied criteria, a descriptive map of the region relative to the points of interest for temporary accommodation is generated. Based on this map, seven points are selected as suitable points for the construction of a temporary accommodation camp, and by employing the analytical hierarchy procedure with Expert Choice software, the optimal point is chosen.

Shahid Abbaspour Dam area is the largest part of Karun Bozorg which had been exposed to floods twice in 2019. Given the strategic importance of the basin and recent floods, the present study uses Landsat 7, Landsat 8, Google Earth Engine and NDWI algorithms at low cost and low time to isolate the water phenomenon from other phenomena for review purposes. The changes in the area of water areas in this basin have been addressed in the last 20 years. In these studies, the area of all water areas in this catchment area has been calculated and based on it, the changes in the area of water areas resulting from the floods of recent years; It is really efficient to obtained for manage water resources, crisis management and achieves sustainable development of the region. The results of this research show that in the mentioned basin, the highest volume of water was in 2019 and the reservoirs and dams were filled with water, so that with a time difference of one month, 25 square kilometers have been added to the volume of their water. According to this information, the risk of flooding covers residential and agricultural areas, and the lack of serious supervision and attention to this issue can cause irreparable damage in the near future.

Identify factors contributing to the spread of soil erosion, and zoning is an important tool to manage and control this phenomenon and to determine appropriate ways are to deal with spreading this phenomenon. In this research, using analysis network process model and GIS has been estimated zoning map of soil erosion in the catchment tea castle in West, East Azerbaijan Province. For this purpose, by examine the sources and expert opinion, effective factors such as slope, aspect, lithology, land use, normalized difference vegetation index (NDVI), annual precipitation and soil is provided of the geographic information system environment. Then, using Arc GIS software and extracted coefficient, the analysis network process prepared to zone map soil erosion and in five classes: very high, high, medium, low and very low. The results showed that 37.12 percent of the area that includes 118.56 square kilometers, is located in risk classes the high and very high. Furthermore, necessary to explain that in the southern and central parts of the basin is the amount of soil erosion high. These areas showed critical situations and acute erosion. And due to dam construction, Ghale Chai should be placed the priority programs and plans of soil conservation and watershed management.

Surface temperature is considered to be a substantial factor in urban climatology. Italso affects internal air temperature of buildings, energy exchange, and consequently the comfort of city life. An Urban heat island (UHI) is an urban area with a significantly higher air temperature than its surrounding rural areas due to urbanization. Annual average air temperature of an urban area with a populationof almost one million can be one to three degreeshigher than its surrounding rural areas. This phenomenon can affect societies by increasing costs of air conditioning, air pollution, heat-related illnesses, greenhouse gas emissions and decreasing water quality. Today, more than fifty percent of the world’s population live in cities, and thus, urbanization has become a key factor in global warming. Tehran, the capital of Iran and one of the world’smegacities, is selected as the case study area of the present research. A megacity is usually defined as a residential area with a total population of more than ten million. We encountered significant surface heat island (SHI) effect in this area due to rapid urbanization progress and the fact that twenty percent of population in Iran are currently living in Tehran.SHI has been usually monitored and measured by in situ observations acquired from thermometer networks. Recently, observing and monitoring SHIs using thermal remote sensing technology and satellite datahave become possible. Satellite thermal imageries, especially those witha higher resolution, have the advantage of providing a repeatable dense grid of temperature data over an urban area, and even distinctive temperature data for individual buildings.Previous studies of land surface temperatures (LST) and thermal remote sensing of urban and rural areas have been primarily conducted using AVHRR or MODIS imageries.

Recently, most researchers use high resolution satellite imagery to monitor thermal anomalies in urban areas. The present study takes advantage of themost recentsatellite in the Landsat series (Landsat 8) to monitor SHI, and retrieve brightness temperatures and land use/cover types.Landsat 8 carries two kind of sensors: The Operational Land Imager (OLI) sensor has all former Landsat bands in addition of three new bands: a deep blue band for aerosol/coastal investigations (band 1), a shortwave infrared band for cirrus detection (band 9), and a Quality Assessment (AQ) band. The Thermal Infrared Sensor (TIRS) provides two high spatial resolution thirty-meter thermal bands (band 10 and 11). These sensors use corrected signal-to-noise ratio (SNR) radiometric performance quantized over a 12-bit dynamic range. Improved SNR performance results in a better determination of land cover type. Furthermore, Landsat 8 imageries incorporate two valuable thermal imagery bands with 10.9 µm and 12.0 µm wavelength. These two thermal bands improve estimation of SHI by incorporating split-window algorithms, and increase the probability of detectingSHI and urban climatemodification. Therefore, it is necessary to design and use new procedures to simultaneously (a) handle the two new high resolution thermal bands of Landsat 8 imageries and (b) incorporate satellite in situ measurement into precise estimation of SHI.Lately, quantitative algorithms written for urban thermal environment and their dependent factors have been studied. These include the relationship between UHI and land cover types, along with its corresponding regression model. The relation between various vegetation indices and the surface temperature was also modelled in similar works. The present paper employ a quantitative approach to detect the relationship between SHI and common land cover indices. It also seeks to select properland coverindices from indices like Normalized Difference Vegetation Index (NDVI), Enhanced Vegetation Index (EVI), Soil Adjusted Vegetation Index (SAVI), Normalized Difference Water Index (NDWI), Normalized Difference Bareness Index (NDBaI), Normalized Difference Build-up Index (NDBI), Modified Normalized Difference Water Index (MNDWI), Bare soil Index (BI), Urban Index (UI), Index based Built up Index (IBI) and Enhanced Built up and Bareness Index (EBBI). Tasseled cap transformation (TCT) which is a method used for Landsat 8 imageries, compacts spectral data into a few bands related to thecharacteristics of physical scene with minimal information loss. The three TCT components, Brightness, Greenness and Wetness, are computed and incorporated to predict SHI effect.The main objectives of this research include developing a non-linear and kernel base analysis model for urban thermal environment area using support vector regression (SVR) method, and also comparing the proposed method with linear regression model (LRM) using a linear combination of incorporated land cover indices (features). The primary aim of this paper is to establish a framework for an optimal SHI using proper land cover indices form Landsat 8 imageries. In this regard, three scenarios were developed:  a) incorporating LRM with full feature set without any feature selection; b) incorporating SVR with full feature set without any feature selection; and c) incorporating genetically selected suitable features in SVR method (GA-SVR). Findings of the present study can improve the performance of SHI estimation methods in urban areas using Landsat 8 imageries with (a) an optimal land cover indices/feature space and (b) customized genetically selected SVR parameters.

The present study selects Tehran city as its case study area. It employs a quantitative approach to explore the relationship between land surface temperature and the most common land cover indices. It also seeks to select proper (urban and vegetation) indices by incorporating supervised feature selection procedures and Landsat 8 imageries. In this regards, a genetic algorithm is applied to choose the best indices by employing kernel, support vector regression and linear regression methods. The proposed method revealed that there is a high degree of consistency between affected information and SHI dataset (RMSE=0.9324, NRMSE=0.2695 and R2=0.9315).

Landslide is one of the most important types of natural disasters,which endangers lives and financial security of many people and destroys environment and natural resources.With the present population growth and expansion of urban areas towardsteep areas and hillsides, landslide-related losses can be catastrophic. For an instance, landslides in Badakhshan Province in Afghanistan killed around 2,700 people in 2014, and a landslide in China (Shan’xiprovince)resulted in the disappearance of 64 people in 2015.Therefore, assessingthe possibility of landslides occurrence seems to becrucial. Providing zoning maps is one of the measures which makes identification of areas prone to future landslides possible. Inferences drawn from these maps can be used for land use planning, prevention of unauthorized construction activities, infrastructure development, refurbishment and restoration.

The present research selectsEast Rudbar-e Alamut (a district of Qazvin province), which is affected by landslides and instability of hillsides, as the study area. It takes advantage of Shannon entropy and information value models to develop landslide susceptibility map of the study areain GIS environment.Shannon entropy theory has been used in extensive researcheswith the aim of prioritizinginfluential factors in the probable occurrence of natural disasters such as landslide. Information value (IV) model is one of the statistical models drawn from information theory with a widespread application in the modeling of geological hazards and disaster risk assessment. Information value model aims to find a combination of significant factors anddeterminetheir impacton theoccurrence oflandslide in an area.To implement this model, relevant data and its related criteria maps were prepared. In this regard, the location of previous landslide events in the study area was determined based on the information received from Forests, Range and Watershed Management Organization. 49 landslides were identified in this way. Then, data was randomly divided into 2 categories: training data and validation data. Thus, 70% of data (35 landslides) were used to produce the models and the remaining 30% (14 landslides) were used for validation purposes. In addition to previous landslides, a collection of topographic, environmental and climatic characteristics of the study area including seven criteria of lithology, slope, distance from faults, land use, precipitation, slope-direction and elevation were selected as the most effective independent variablesto produce criteria maps with 30-meter spatial resolution. Basic information used to map these seven influential factors was obtained from Forests, Range and Watershed Management Organization, as well as the SRTM Digital Elevation Model (DEM), and used after some modifications. Considering the capability of ArcGIS in spatial data analysis, thissoftwarewas used to produce information layers and implement the models.

Prioritizing influential factors using Shannon entropy model introducesthree factors (i.e. land use, elevation and precipitation)as the most significant factorsin the occurrence of landslides in the study area. Factors of slope angle, distance from faults (almost equal to slope angle), lithology and slope-direction were in the next influential factors.Also, results of information value model indicate that looking from lithology perspective, the category of marl, calcareous sandstone, sandy limestone and minor conglomerate has an information value of 1 and thus, the highest probability of landslide occurrence. Category of basaltic volcanic rocks, along with category of well bedded green tuff and tuffaceous shale have the lowest probability of landslide occurrence with information values ​​of -2.03 and -1.70, respectively.Only two categories of theslope angle criterionhave a positive-index. The highest information value (0. 93) in this category occurs in the class of 5-12 degrees, followed by the class of 12-20 degrees. The lowest information value occurs in slopes of more than 30 degrees. Based on this observation, it can be clearly concluded that the slope angles of 5 to 20 degrees are most prone to landslides. Distance to faults criterion indicate that the category of500 to 1000-meter distance to faultshave the highest information value (1.67). Regarding land use criterion, three land uses of garden, agriculture and garden-agriculture have the highest information values ​of 2.16 and 1.59 and 1.11, respectively. Regarding precipitation, average annual rainfall of less than 400 millimeters have the highest information value (1.50). Regardingslope-direction criterion, most landslides occur in southwest, south and eastdirections.Northeast, west, and northwest directions have the lowest probability of landslide occurrence, respectively. In terms of elevation, the information value is reduced as the height increases, and the maximum information value is related to the elevations of less than 1200 meters.After assigning a weight to each criterion and related classes, the landslide risk zone map was generated based on Shannon entropy and information valuemodels. The resulting zoning map produced based on natural breaks methods dividesthe area into five classeswith very high, high, moderate, low and very low risk. Resultsof Shannon entropy modelindicate that out of 14 landslides considered as the validation data, 3, 7, 2, 1, 1 landslideshave occurred in very high, high, moderate, low and very low risk zones, respectively. Resultsof the information value modelindicatethat 8, 4, 0, 1, 1 landslideshave occurred in very high, high, moderate, low and very low risk zones, respectively.

Evaluation of results using experimental probability index indicates that with 86% experimental probability,both models of Shannon entropy and information value are effective inidentification of landslide hazard in the East Rudbar-e Alamut region. Also, considering the number of landslides in very high and high risk zones, Shannon entropy and information value modelshave an experimental probability index of 72% and 86%, respectively, which prove higher efficiency of information value model. In Shannon entropy model, total area of very high, high and moderate risk zones covers 34% and 56% of the study area,respectively. In information value model,total area of very high and high risk zones covers 20% and 29% of the study area, respectively. Based on the landslide risk zone map, high and very high risk zones are mainly located in the west of the study area

Every year, many human beings suffer from an earthquake as a near-unpredictable natural disaster and its devastating human and financial losses. Management of such crises is related both before and after the crisis. Relief and rescue is only a stage in the occurrence of disasters can be studied in advance of the crisis to provide a solution to improve the performance of relief and rescue teams during the crisis. In this study, using a spatial information system and particle swarm algorithm and simulating a presumptive earthquake, a solution is suggested for optimal management of relief and rescue teams in earthquake. In this method, an earthquake, and 32 relief workers of four operational teams in 148 housing complexes simulated to study area in Tehran. Rescuers, with the help of particle swarm algorithm in a spatial information system, were allocated relief and rescue activities, in less time, would provide relief and rescue more efficient than the empirical mode. The use of this method to optimize simulations, as well as to implement the scientific and practical structure of relief and rescue teams and activities, will be a new way to improve the quality of relief and rescue after the earthquake. The results of the proposed method of this research showed Performance improvements of about twofold.

One of the issues that most of the world's major cities face is the issue of natural disasters. The nature of the overwhelming majority of natural disasters and the need for quick and correct decision-making and implementation of operations has created knowledge of "crisis management". This knowledge refers to the set of activities that occur before, during and after the occurrence of disaster, in order to reduce the probable vulnerability caused by the occurrence of these events [5]. It is necessary to carry out all the affairs and actions necessary to achieve the goals outlined in the above definition, which requires the assumption of operational roles by operational teams [3]. Given the importance of relief and rescue at the time of natural disasters to save lives and property, the proper allocation of aid workers to activities is necessary.
In order to improve the relief and rescue operation, firstly, activities were carried out at the time of the earthquake, and comprehensive information were obtained on the post-earthquake relief and rescue mode. In order to allocate people, using optimization methods, considering the conditions of this research, is effective in improving the efficiency and effectiveness of post-earthquake relief. Hence, due to the nonlinear relations of this study, and in light of previous research, the particle swarm optimization algorithm was chosen as a suitable method for solving this problem. Moreover, also the use of a spatial information system for modeling, displaying, and updating of force information, activities, and conditions of earthquake area is suitable for optimal forces management [1].

This section examines the responsibilities of rescue workers in the earthquake crisis and important points in the earthquake relief process. Some search and rescue actors include four components of locating, evaluating, fixing, and transferring [19]. First, the location and release of individuals and the medical assessment and, if necessary, the use of primary care, emergency treatment (stabilization) and transfer to treatment centers are carried out [26]. The rescue team should have a precise program to carry out rescue operations for those in detention.
Search and Rescue Operations Management
To ensure the success of search and rescue operations in urban areas; it must be done very carefully. The relief and rescue program can be divided into five stages, respectively [26]:
Primary Identification - Data Collection (Preliminary Assessment)
Quickly assess the area (Technical Inspection)
Surface Search and Rescue in the Damaged Area (Primary Rescue)
Search and rescue by technical means (Secondary Rescue)
The systematic removal of debris (Final Collapse Lifting)

On the other hand, seven steps in search and rescue operations are assumed to be considered by the savior’s people [9]:
Data collection: One of the first steps to be taken is to assess and assess the situation.
Evaluation of Damage: By looking at different angles to the buildings.
Identifying resources and accessing them: including access to facilities, equipment, and personnel.
Priority: Includes emergency diagnosis and safety assurance for the continuation of search and rescue operations. Sometimes a building should be marked in such a way that no other person enters it and waits for other forces or more facilities.
Designing a Rescue Plan: In this section, it becomes clear who and with whom the conditions will enter the building.
Guidance for search and rescue operations: Search for people under the rubble remains and caught
Evaluate progress: The situation must always be checked to assess the progress of the rescue program and to prevent any damage to the relief forces

The first attempt by Kennedy and Eberhart, simulating the social behavior of birds in 1995, presented the particle group optimization method. The components of a group follow a simple behavior. In this way, each member of the group imitates the success of their other neighbors. The purpose of such algorithms is to move members of the group to the search space and to accumulate at an optimal point (such as the source of food).

To achieve relief & rescue optimal management, close interaction is being necessary [25]. The results of this study showed that parameters such as the duration of survival under the rubble, the duration, the distance between people and the location of activities, the speed of people when moving to the goal of the relief worker is very important in fulfilling the task. With the studies and studies, finally, the relation one was designed, which is a continuous nonlinear relationship. According to the studies, the method of optimizing the congestion of particle capabilities solves these functions, and this method allocates individuals to activities in this research is optimized:(1)
In the above relationship, all parameters must follow a unit or reputation [24], “Max Injured” the most injured number among the wounded of each residential building, “Area Assigned” is the area [20], which the same ​​activity is located inside it. “Spacing” the relief distance to the operating area and the “search time” and “search speed” are respectively the duration of the work and the speed of the relief worker. If a rescuer will be sent to a region that is estimated to be several people under debris, the duration of activities will be multiplied by the number of submarines. Moreover, to achieve the final cost of an activity that requires several people, it must be summed of the costs from each who performs that activity.

The cased study is a part of the central region of ​​Tehran. The relief and rescue activities of the earthquake crisis include Searching, Light Collapse Lifting, Heavy Collapse Lifting, Primary Helping, Securing, Pointing, Securing Pilot, Air update in the rubble, reconstruction of the network of roads [6, 19]. In this research, 32 reliefworkers of four operational teams [22], and at the beginning of the operation, they are deployed at the nearest crisis management center to the study area. Figure 1 shows the first study area and the initial position of the relief workers in the study area.
Figure 2: Study area and the first location of rescuers
The following shows building and human damages data showing the initial phase of earthquake simulation, which includes 22 out of 148 damaged sites, and the descriptive information of relief workers in a hypothetical earthquake, in which 14 relief workers out of 32 relief workers, as well as the third, are shown their activities:
Fig. 2. building and human damages
Fig. 3. Descriptive information of relief workers
Regarding the parameters stated in the method of implementation (i.e.; the descriptive information of the rescuers, the activities and initial damages of the earthquake), the proposed algorithm of this research, is evaluated and calculated by using relations discussed for all the rescuers in all the housing complexes. And eventually, the allocated of relief workers to the activities was obtained. An example of the optimal mode of relief and rescue teams is showing in the figure below.
Fig. 4. Optimization of the Relief & Rescue Team
In the study area of ​​the image above, the “Rescuers 34” relate to relief workers assigned to Light Collapse Lifting activities; “Rescuers32”, relief workers, and Pointing; “Rescuers31”, rescuers assigned to Searching activities. As well as “Rescuers33” for rescue workers who are engaged in Securing Pilot and relief workers “Rescuers 37”, engaged in Primary Helping activities. The allocation of people is carried out according to the priority, and the residential areas that have more damage are in the priority of the relief effort.
In evaluating the efficiency of the proposed algorithm, the positive effect of the initial population selection method shown in the results obtained from the implementation of the proposed algorithm. Finally, a 2.2 fold improvement in the results obtained from the state that was not used by this algorithm. In the table below, the calculation of the cost function in the two modes of implementation of the proposed algorithm and its non-implementation is set, which represents the calculating the cost of the allocation in the two situations for the entire operational team.
Table 1. Comparison of the results of the proposed algorithm and its validation
Used model
Cost calculated for the entire operational team
Without using the proposed algorithm
0.564
Using the proposed algorithm
0.252

Due to the facts that the problem is considered to be grouped of the subject of this research, the effectiveness of each person's activity on the other people's activities, and the group and the category of operations, as well as the structure of the particle swarm algorithm, which allows for more repetition in less time, the proposed algorithm of this study is identified as an appropriate solution to the post-earthquake relief and rescue problem.
The structure of the particle swarm algorithm is continuous; because of the discrete structure of the present, it is implemented discretely by applying changes to the structure of this algorithm. As previously stated, the context of individuals, their specializations, the activities, and the damaged sites have the same priorities as those that were implemented in the algorithm.
Using the proposed algorithm of this research and applying the changes expressed in it, in order to optimize and implement the scientific and practical structure of relief and rescue operation activities and teams, is a novel and effective way to improve the quality of relief and rescue after it will be an earthquake. Finally, as shown in Table 1 in the findings, the proposed algorithm implementation in this study improved the 2.2% of the results from the allocation of relief workers to a state that was not used by the proposed algorithm of this study.
For future researches, the optimization methods such as simulated annealing, ant colony, genetics, and game theory are suggested.

One of the most important steps in designing seismic rehabilitation programs for urban areas in metropolises is determining the region’s risk level. There is no doubt of cast that the implementation of a seamless seismic rehabilitation program can secure urban areas against seismic damages. However, these plans face some difficulties in most cases. For instance, their execution costs are very high, or, they are really time-consuming in some cases. It is obvious that a coherent rehabilitation strategy will be effective if it can distinguish executive priorities by identifying areas with high level of risk, and focusing primarily on improving the general condition of these areas. In this paper, the seismic vulnerability of the 7th region of Tehran is studied. For this purpose, the main parameters that affect the vulnerability of the region are identified first. Then, considering the uncertainties in these parameters, their impact on the vulnerability of the region is determined. ArcGIS software is used not only to model the impact of these parameters, but also for data training. The results of this study show the vulnerability level of different parts of Tehran's 7th region. In fact, this study provides a comprehensive and meaningful overview of the current status of the studied region, which can helps the engineers to design suitable seismic rehabilitation programs. The obtained results are presented as a map which reveals the vulnerability level of the region. Based on such a map, neighborhoods or areas with high seismic hazard are identified. Considering these classification, the seismic rehabilitation techniques can be offered purposefully.

The ubiquity of mobile devices, such as smart phones and tablets, has contributed to the development of pervasive systems, including navigation and health systems. The main characteristicsof pervasive systems are the necessity of dynamic reconfiguration and proper adaption to the continuous changes in different contexts. The existence of dynamic capabilities has been considered in the design and implements of a context aware system, including context acquisition, context understanding and computing, decision making, and context presentation.Context acquisition: This domain of research focuses on using personal sensing devices which measure various parameters by means of portable devices and save them on the external/internal database for further processing. The aim of researches is collecting, sharing, and/or reusing data in other applications or through a web interface.Context understanding and computing: The most works are in the field of context monitoring, data management, understanding or computing. The ability to automated context reasoning about various types of contexts and their properties are considered using various context models and algorithms. Most applications are customized for a specific case such as air pollution, tourist, navigation, and health-care. Context presentation: This category of research has commonly focused on context-aware application adaptation. The adaption happens between the real world, the map and user’s location and orientation. A number of studies have been carried out in the field of tourist guides or navigation adapting the presentation style to the changing requirements of the user.Most studies in ubiquitous health care have only been carried out in a small number of areas and using external portable sensors and developing applications on mobile phones. A major problem with this kind of applications is collecting and sharing data, monitoring, or reasoning without having an active role in decision making in different environmental conditions. Using external tools such as portable devices is costly and limits using the systems.
This paper has focused on design and implementation of a context aware ubiquitous system which has been customized for severe environmental conditions (in particular, air pollution). Air pollution is a spatio-temporal phenomenon and it causes changes in health conditions and it increases mortality. Eclipse Kepler software, java, PHP programming language and MySQL and SQLit database and also Google Maps API was used in this research. The proposed system design approach is based on distributed architecture in the portion of data collection and processing. Data collecting is done by means of software and hardware sensors. The context aware system is able to automatically identify the user’s context and represent required data and information after computing and reasoning. Contexts based on their impact on the decision-making process can be divided into two categories: passive and active contents.We used an active context in the research such as time, location, traffic, direction, air pollution. Collecting required data is done automatically with high speed and accuracy, and data plays an active role in decision making. In the system architecture, servers were embedded to enter data automatically and only data relating to health conditions is entered manually. Processing environment was divided into two parts, in case of abounding calculations, processing is transferred to the server so that only light processing is performed on the client. At every stage of the process, the user interface provided outputs in the form of recommendations and notifications. The system represents user-friendly environment. Context information can be posted on the process server and retrieved from the history. The proposed system can become an important tool to enable patients to be aware of air pollution conditions, not only to be applied in managing and monitoring their health information, but also in decision making, finding the best solution in severe environment, sharing data and communicating with family and doctor. The application represents suitable solution for solving the shortest path problem according to spatial-temporal and traffic condition. In fact, the path with the lowest level of air pollution is chosen as the best path.The system indirectly encourages greater use of the ubiquitous health system and motivates patients to acquire an active role in their health management and helps them to improve their health condition. The information collected and posted on the server can be reused in professional station and it presents useful information to health experts. We broadly concern about patients’ privacy in the design of the system.