Annual Variation of the Height of Urban Boundary Layer of Tehran

Document Type : Full length article


1 PhD Student in Synoptic Climatology, University of Tehran, Iran

2 Associate Professor of Climatology, University of Tehran, Iran

3 Assistant Professor of Climatology, University of Tehran, Iran

4 Professor of Climatology, University of Canterbury, New Zeland


Atmospheric Boundary Layer (ABL) is the lowest layer of the Troposphere affected by the land surface. Unlike the free atmosphere above, the land surface has a significant effect on the ABL and this layer, is the only part of the atmosphere in which the effects of friction and the diurnal variation of temperature can be observed. In fact, the ABL plays the role of a dealer of energy and mass between the land surface and the free atmosphere.
Many efforts have been undertaken to understand the behavior of the ABL due to its importance since almost all human activity (except aviation) forms in this layer. This, however, should be mentioned that those efforts are in a strong correlation with the development level of the country and the methods to study the ABL varies country by country. The study methods for ABL observations are divided into three main branches: a) Modeling and Simulation (using dynamic and numerical models); b) Numerical estimation (using atmospheric profiles); and c) Using High-tech devices (Remote Sensing, RADAR, and LiDAR).
Like most climatology phenomena, the ABL probably owns some regular pattern in different time scales. In this research, the variation of the ABL over the city of Tehran is investigated in an annual time scale. The reason to choose the city is the high population, high urban concentration, and constant integration with inversion and air pollution, and above all, lack of knowledge about the ABL over the city.
In order to conduct a thorough research, long term data were required from a variety of sources. Therefore, 30 years of data were gathered in daily time scale (at 00:00 GMT and 12:00 GMT equal to 03:30 and 15:30 Tehran Mean Time respectively). Total data records reached a sum of 10958 for each parameter. Data were collected from the European Center for Mid-range Weather Forcast (ECMWF) and the Atmospheric Science Department of the University of Wyoming. Data gathered from the ECMWF was in NetCDF format and the air profile gathered was in Notepad form.
To have solid results, the data should be of high certainty and reliability. To investigate the reliability of the ECMWF data, the ABL height was calculated for some days employing the Advanced Parcel Method which defines the top of the ABL as the height in which the virtual potential temperature is equal to that of the surface values. Using Kuchran’s method to estimate the volume of the test subject from a society, 372 days were randomly selected and the ABL height was calculated based on radiosonde profiles. Then, the correlation coefficient and the Root Mean Square Error (RMSE) between the two data sets were calculated and since the coefficients were significant, data was used further. As for the maps, data gathered from the ECMWF was averaged usnig the ECMWF tools. This tool uses the R (or Rstudio) software abd the RBridge to calculate integrated NetCDF data. The maps were drawn in ArcMap 10.
Results and discussion
Mean daily ABL height was located almost at 850 meters above ground level (AGL). While nighttime ABL descends to almost 80 meters AGL, the daytime ABL rises to 2300 meter AGL. The ABL height has experienced a total rise of 5 meters per day in the total time period of the research. Spatially, the lowest ABL heights were experienced in northeastern part of the city and the highest ABL were measured in south and southwestern part of the city. Also, the correlation between the mean daily data (the ECMWF uses eight data measurements to calculate the mean daily data) and the maximum and minimum averages were calculated. The results showed that the maximum values are of higher influence in mean daily data rather than the minimum.
The findings also indicate that, during all 6 periods of the study time scale, the position of the minimum and maximum boundary layer values are almost the same. High values occur at south and southwest of the city and low values occur at northeast of the city. There is also a midsection area that is usually extended from Ka valley to the east of the city. There is also a sign of some core area, especially the one located above the Pardisan Park that sometimes affects the ABL patterns.
Furthermore, the correlation between the ABL and some climatic parameters (e.g. sunshine hours, surface temperature, surface heat and moisture flux, relative humidity, air pressure, wind speed at 0, and 10 meter height AGL) was calculated for Tehran. The results indicated that there is some significant correlation between the ABL and climatic parameters. The highest correlation was seen between the ABL variation and the surface temperature. Seemingly, the closer the measured parameters to the surface are, the higher is the correlation coefficient.
The Boundary layer height fluctuates between 80 meters AGL at nights and 2300 meters in daylight. The average height has increased almost five meters per year. However, since 1988 to 2012 it has risen and after 2012 it experienced subsidence. Both in minimum and maximum, the highest boundary layer height has been measured in south and southwestern part of the city and the lowest values were measured in northeastern of the city. The iso-height lines are extended from northwest to southeast. This is probably due to the effect of Alburz Mountain range surrounding Tehran in northern edges. The ABL showed high correlation with some climatic parameters but the coefficient between the ABL height and some other parameters were insignificant. Moreover, regarding the correlation coefficient and the root mean square error results, the Advanced Parcel Method seem to be of enough reliability to calculate the boundary layer height based on radiosonde profiles.


احمدی گیوی، ف.؛ ثابت‏قدم، س. و علی‏اکبری بیدختی، ع.ع. (1387). بررسی نوسان عمق لایة آمیختة جو شهری تهران با استفاده از مدل mm5 و عوامل مؤثر در آن، مجلة فیزیک زمین و فضا، 35(2): ۱۰۵-117.
قسامی، ط.؛ علی‏اکبری بیدختی، ع.ع.؛ صداقت‏کردار، ع. و صحرائیان، ف. (1389). بررسی تغییرات دمای قائم پتانسیل در چند دورة بحرانی آلودگی هوای تهران، علوم و تکنولوژی محیط زیست، 12(3): ۱۳-24.
لشکری، ح. و هدایت، پ. (1386). تحلیل الگوی سینوپتیکی اینورژن‏های شدید شهر تهران، پژوهش‏های جغرافیایی، 38(۱۴۷۳): ۶۵-82.
جهان‏بخش اصل، س. و روشنی، ر. (1393). بررسی شرایط الگوی سینوپتیکی حاکم بر وضعیت‏های وارونگی دمای بسیار شدید شهر تبریز، جغرافیا و برنامه‏ریزی، 18(۴۸): ۸۱-96.
کیخسروی، ق. و لشکری، ح. (1393). تحلیل رابطه بین ضخامت و ارتفاع وارونگی و شدت آلودگی هوا در شهر تهران، جغرافیا و برنامه‏ریزی، 18(۴۹): ۲۳۱-257.
یاوری، ح. و سلیقه، م. (1390). سطوح وارونگی در آلودگی‏های شهر تهران، نشریة تحقیقات کاربردی علوم جغرافیایی، 17(20): ۸۹-105.
کرم‏پور، م.؛ سلیقه، م.؛ طولابی‏نژاد، م. و زارعی چغابلکی، ز. (1395). بررسی آلودگی هوای شهر تهران به روش وارونگی بحرانی هافتر، نشریة تحلیلفضاییمخاطراتمحیطی، ۳(۱): ۵۱-64.
Ahmadi Givi, F.; Sabetghadam, S. and Ali Akbari Bidokhti, A. (2009). Investigation of Mixed Layer Depth Variations of Tehran using mm5 model, Journal of Earth and Space Physics, pp. 105-117.
Bachour, D. and Perez-Astudillo, D. (2013). Boundary Layer Height Measurements over Doha Using Lidar, Energy Procedia, 57: 1086-1091.
Beljaars, A.C.M. and Betts, A.K. (1992). Validation of the boundary layer representation in the ECMWF model, ECMWF Seminar, Proceedings: Validation of Models over Europe, Vol. II. Reading, UK, 7-11 September 1992.
Chen, Z.; Liu, W.; Zhang, Y.; He. J. and Ruan, J. (2011). Mixing layer height and meteorological measurements in Hefei China during the total solar eclipse of 22 July, 2009, Optics & Laser Technology, 43(1): 50-54.
Ghasami, T.; Ali Akbari Bidokhti, A.; Sedaghatkerdar, A. and Sahraiyan, F. (2012). Investigation of Potential Vertical Temperature in some accute air pollution events in Tehran, Science and Technology of the Environment, pp. 13-24.
Hennemuth, B. and Lammert, A. (2005). Determination of the Atmospheric Boundary Layer Height from Radiosonde and LIDAR Backscatter, Boundary-Layer Meteorology, 120: 181-200.
Holzworth, C.G. (1964). Estimates of mean maximum mixing depths in the contiguous United States, Monthly Weather Review, 92: 235-242.
Jahanbakhsh Asl, S. and Roshani, R. (2015). The Study of Synoptic Patterns Dominating on the Very Intensive Temperature Inversion in Tabriz, Journal of Geography and Planning, pp. 81-96.
Karampour, M.; Saligheh, M.; Toulabinejad, M. and Zarei Choghabaki, Z. (2016). Evaluation of air pollution in Tehran city by Hefter's critical Inversion method. Jsaeh., 3(1): 51-64.
Keykhosrowi, Gh. and Hasan, L. (2011). Analysis of the Relationship between the Thickness and Height of the Inversion and the Severity of Air Pollution in Tehran, Journal of Geography and Planning, pp. 231-257.
Lashkari, H. and Hedayat, P. (2008). Synoptic Analysis of Extreme Inversions in Tehran, Geographical Researches, 38(1): 65-82.
Mao, F.; Gong, W.; Song, S. and  Zhu, Z. (2013). Determination of the boundary layer top from lidar backscatter profiles using a Haar wavelet method over Wuhan, China, Optics & Laser Technology, 49: 343-349.
McKendry, I.G., D.van der Kamp, K.B.Strawbridge, A.Christen, Crawford, B. (2009). Simultaneous observations of boundary-layer aerosol layers with CL31 ceilometer and 1064/532 nm lidar,Atmospheric Environment, 43(36): 5847-5852.
Pal, S.; Xueref-Remy, I.; Ammoura, L.; Chazette, P.; Gibert, F.; Royer, P.; Dieudonné, E.; Dupont, J.C.; Haeffelin, M.; Lac, C.; Lopez, M.; Morille, Y. and Ravetta, F. (2012). Spatio-temporal variability of the atmospheric boundary layer depth over the Paris agglomeration: An assessment of the impact of the urban heat island Intensity, Atmospheric Environment, 63: 261-275.
Ribeiro, F.; Oliveira, A.; Soares, J.; Miranda, RM.; Barlage, M. and F.Chen (2018). Effect of sea breeze propagation on the urban boundary layer of the metropolitan region of Sao Paulo, Brazil, Atmospheric Research, 214: 174-188.
Seibert, P.; Beyrich, F.; Gryning, S. E.; Joffre, S.; Rasmussen A. and Tercier, P. (2000). Review and Intercomparison of Operational Methods for the Determination of the Mixing Height, Atmos. Environ., 34: 1001-1027.
Shukla, K.K.; Phanikumar, D.V.; Newsom, R. K.; Kumar, K.N.; Ratnam, M.V. and NarendraSingh, M.N. (2014). Estimation of the mixing layer height over a high altitude site in Central Himalayan region by using Doppler lidar,Journal of Atmospheric and Solar-Terrestrial Physics, 109: 48-53.
 Sugiyama, G. and Nasstrom, J.S. (1999). Methods for Determining the Height of the Atmospheric Boundary Layer, Report for Lawrence Livermore National Library, United States of America.
Wagner, P.; Schäfer, K. (2017). Influence of mixing layer height on air pollutant concentrations in an urban street canyon, Urban Climate, 22: 64-79.
Wang, Ch.; Shi, H.; Jin, L.; Chen, H. and  Wen, H. (2015). Measuring boundary-layer height under clear and cloudy conditions using three instruments, Particuology, 28: 15-21.
Yavari, H. and Saligheh, M. (2011). Inversion Levels in Tehran’s Pollution, Researches in Geographical Sciences, pp. 89-105.
Zéphoris, M.; Holin, H.; Lavie, F.; Cenac, N.; Cluzeau, M.; Delas, O.; Eideliman, F.; Gagneux, J.; Gander, A. and  Thibord, C. (2005). Ceilometer observations of aerosol layer structure above the Petit Lubéron during ESCOMPTE's IOP 2, Atmospheric Research, 74(1-4): 581-595.