Journal of Spatial Analysis

---

Air pollution has become one of the main problems of cities. Among the sources of air pollution, vehicular traffic plays an important role. Planning for efficient management and control of the air pollution caused by vehicular traffic requires accurate information on spatio-temporal dispersion of the pollutions. This research studies 3D spatio-temporal dispersion of NOx pollution caused by vehicular traffic at Valieasr-Fatemi intersection resides in Tehran, Iran. It is selected for being crowded and having the required meteorological and pollution data sensed by the Air Quality Control Corp. of Tehran Municipality.

This study uses GRAL that is a local micro-scale air dispersion model defined based on Euleran-Lagrangian dispersion models. It investigates the level of spatio-temporal autocorrelation generated by GRAL simulations at both 2D and 3D modes and discusses how it adapts with the reality.

Adopting the GRAL air pollution dispersion model, streets are defined as the linear source of pollution of NOx caused by vehicular traffic. The traffic rate is estimated based on street areas and directions, the designed average traffic velocity, traffic volume and car passage counting at the intersection. The 3D geometry of the buildings is also added to the model. All the required data that were available for winter of 2007 are gathered and introduced into the model.

The model is executed at 9 heights vary from 1.7 m to 52.5 m. These heights are defined covering a range from an average human level height to average building height and above. These levels are considered both separately in 2D mode and integrated into a 3D mode. The formation of NOx clusters is investigated analyzing their autocorrelation using Moran Index at global and local scale.

The calculated Moran-I at global scale at each 9 levels of heights, varies from 0.7 to 0.9 that depicts the validity of the GRAL model adopted to simulate the expected autocorrelation of pollution density affected by spatial issues. The Moran-I increases at higher levels as less air turbulence happens. However the result show that the turbulence increases temporarily at about 10m to 15m which are the average building heights. At local scale, the Moran-I/Anselin shows that HH clusters dominate at lower levels, around streets central areas that are farther from the buildings, and around the intersections. At higher levels, esp. higher than buildings average height, the LL clusters dominate. However the HH clusters formed around intersections, while are shrank, are still visible at high levels. The turbulence caused by building fronts and their down wash effect is also shown in the result as no definite cluster is formed near the buildings front and back.

The autocorrelation analysis is also carried for an integrated 3D model consists of all the 9 levels of heights. Considering the weight matrix for a 20m 2D neighborhood and 1m/s dispersion of the pollution vertically, the global calculated Moran-I equals 0.229 which shows existence of a spatio-temporal autocorrelation of the results generated by GRAL. At local scale the results show that the HH clusters have higher temporal dispersion rate than LL clusters.

---

Understanding the changes in extreme precipitation over a region is very important for adaptation strategies to climate change. One of the most important topics in this field is detection and attribution of climate change. Over the past two decades, there has been an increasing interest for scientists, engineers and policy makers to study about the effects of external forcing to the climatic variables and associated natural resources and human systems and whether such effects have surpassed the influence of the climate’s natural internal variability. The definitions used in the 5th assessment report were taken from the IPCC guidance paper on detection and attribution, and were stated as follows: “Detection of change is defined as the process of demonstrating that climate or a system affected by climate has changed in some defined statistical sense without providing a reason for that change. An identified change is detected in observations if its likelihood of occurrence by chance due to internal variability alone is determined to be small.  Attribution is defined as the process of evaluating the relative contributions of multiple causal factors to a change or event with an assignment of statistical confidence”. Detection and attribution of human-induced climate change provide a formal tool to decipher the complex causes of climate change. In this study the optimal fingerprinting detection and attribution have been attempted to investigate the changes in the annual maximum of daily precipitation and the annual maximum of 5-day consecutive precipitation amount over the southwest of Iran.

This is achieved through the use of the Asian Precipitation—Highly Resolved Observational Data Integration Towards Evaluation of Water Resources Project(APHRODITE) dataset as observation, a climate model runs and the standard optimal fingerprint method. To evaluate the response of climate to external forcing and to estimate the internal variability of the climate system from pre-industrial runs, the Norwegian Climate Center’s Earth System Model- NorESM1-M was used. We used up scaling to remap both grid data of observations and simulations to a large pixel. This remapped pixel coverages the area of the southwest of Iran. The optimal finger printing method needs standardized values like probability index(PI) or anomalies as input data, since the magnitude of precipitation varied highly from one region to another. The General Extreme Value distribution (GEV) is used to convert time series of the Rx1day and Rx5day into corresponding time series of PI.  Then we calculated non-overlapping 5-year mean PI time series over the area study. In this research, we applied optimal fingerprinting method by using empirical orthogonal functions.  The implementation of optimal fingerprinting often involves projecting onto k leading EOFs in order to decrease the dimension of the data and improve the estimate of internal climate variability. A residual consistency test used to check if the estimated residuals in regression algorithm are consistent with the assumed internal climate variability. Indeed, as the covariance matrix of internal variability is assumed to be known in these statistical models, it is important to check whether the inferred residuals are consistent with it; such that they are a typical realization of such variability. If this test is passed, the overall statistical model can be considered suitable.

Results obtained for response to anthropogenic and natural forcing combined forcing (ALL) for Rx1day and Rx5day show that scaling factors are significantly greater than zero and consistent with unit. These results indicate that the simulated ALL response is consistent with Rx1day observed changes. Also, it is found that the changes in observed extreme precipitation during 1951-2005 lie outside the range that is expected from natural internal variability of climate alone and greenhouse gasses alone, based on NorESM1-M climate model. Such changes are consistent with those expected from anthropogenic forcing alone. The detection results are sensitive to EOFs. We estimate the anthropogenic and natural forcing combined attributable change in PI over 1951–2005 to be 1.64% [0.18%, 3.1%, >90% confidence interval] for RX1day and 2.5% [1%,4%] for RX5day.

---

Tehran metropolitan with its large population, daily migrant workforce and many students, needs to planning and designing watch/warning system to reduce the climatic problems for human health.for this purpose, we need to study the climate accurately and Since the factors affecting the climate of warm and cold periods in Iran are different, in this study , the meteorological variables of Tehran warm period (May to September 2002) turned into 4 components in Temporal Synoptic Index (TSI) using PCA Method and using P-Array and Varimax rotation.By the scores of components for each day, the clustering method (in ward method) were used and, the warm days of the year was divided into two cluster named favorable and oppressive airmasses. The average maximum air temperature that is more effective in mortality, was 36.13 ° C. Days with temperatures above  34 ° C, less pressure, mild winds , dryness and more sunshine resulted in more adverse weather conditions, which resulted in a 34% increasing in mortality compare with favorable weather. The total number of deaths from cardiovascular disease during the study period was 154046 that about  67%of deaths have been simultaneous with oppressive airmass.The epidemiological study of mortality also confirms the results of previous research in this area and shows that the incidence of mortality is higher in older people as well as in men. It is clear that not all mortality can be attributed to the effects of climate, but results show that change in climatic conditions will affect on mortality and also for study the effect of climatic hazards on human health, it is better that we study the effect of all variables together on humans.

---

Conclusion
The results showed that the length of the precipitation period in the study area is 8 months (out of 213 systems, 60 systems were recorded in January and 50 systems were recorded in December). June, July, August and December did not have a comprehensive rainfall system. The difference in height between the Mediterranean Cyclone  and the Arabian Anticyclone, i.e. the MAI index, was identified in 5 categories. The results showed that the higher the value of MAI index, the less rainfall systems and the higher the rainfall intensity. Most rainfall systems occur when the MAI index reaches 100 to 150 geopotential meters. In addition, the concentration of precipitation occurred in February when the MAI index reached 200 geopotential meters and above. In the continuation of the research, the difference between the Sudan Cyclone and Arabian Anticyclone, i.e. the SAI index, was also identified in 5 categories. The results showed that most of the systems occurred in the third category with a value of SAI index between 100 and 150 geopotential meters. Moreover, the results showed that the higher the value of the SAI index (, the lower the number of rainfall systems and the higher the average of rainfall systems.