OTHERS_CITABLE بررسی لکه‌های زیستگاهی خرس سیاه بلوچی (Ursus thibetanus gedrosianus) با استفاده از سنجه های سیمای سرزمین (مطالعه موردی: مناطق بحر آسمان و زریاب استان کرمان) تحلیل زیستگاه با استفاده از سنجه‌های سیمای سرزمین می‌تواند در راستای مدیریت بهتر زیستگاه کارآمد باشد. خرس سیاه بلوچی به‌عنوان زیرگونه درخطر انقراض در مناطق بحر آسمان و زریاب استان کرمان پراکندگی دارد. هدف از این مطالعه مدل‌سازی پراکنش گونه و ارزیابی کیفیت لکه‌های زیستگاهی با استفاده از سنجه‌های سیمای سرزمین است. مدل‌سازی پراکنش با استفاده از رویکرد اجماع (Ensemble)، حاصل از تلفیق مدل‌های فقط حضور و حضور/ شبه عدم حضور در نرم‌افزار ModEco انجام گرفت. روی مدل اجماع از حد آستانه (True Skill Statistics, TSS) استفاده شد و لکه‌های زیستگاهی با استفاده از سنجه‌های سیمای سرزمین تحلیل شدند. نتایج این مطالعه نشان داد که گستره پراکنش، فراتر از مرز مناطق حفاظت‌شده موردمطالعه است و متغیرهای مرتبط با رطوبت بیشترین تأثیر را بر روی حضور گونه داشتند. بر پایه تحلیل سنجه‌ها، لکه‌های زیستگاهی در پناهگاه حیات‌وحش زریاب از پیوستگی بالاتر و حاشیه کمتری نسبت به منطقه حفاظت‌شده بحر آسمان برخوردار هستند. در منطقه حفاظت‌شده بحر آسمان تعداد لکه‌های زیستگاهی بیشتر بوده و عدم پیوستگی لکه‌ها منجر به شکل‌گیری حاشیه شده است. نتایج این مطالعه می‌تواند در راستای مدیریت لکه‌های زیستگاهی این‌گونه درخطر انقراض مورد استفاده قرار گیرد. http://ijae.iut.ac.ir/article-1-1118-fa.pdf 2022-09-01 1 17 10.47176/ijae.11.2.10803 خرس سیاه مدل اجماع لکه‌های زیستگاهی سیمای سرزمین استان کرمان Investigating the Habitat Patches of the Baluchistan Black Bear (Ursus thibetanus gedrosianus), Using Landscape Metrics (Case Study: Bahr Asman and Zaryab Areas, Kerman Province) Habitat analysis using landscape metrics can be efficient in better management of habitat. As a critically endangered subspecies, the Baluchistan black bear is scattered in the Bahr Asman and Zaryab areas in Kerman province. The purpose of this study was to model the distribution of the sub-species and evaluate the quality of its habitat patches, using landscape metrics. Distribution modeling was conducted using an ensemble approach, resulting from the combination of presence-only and presence/pseudo-absence data in ModEco software. True Skill Statistic (TSS) threshold was applied to the ensemble model and the habitat patches were analyzed, using landscape metrics. The results showed that the distribution rage is beyond the border of the studied protected areas and the variables related to humidity had the greatest effect on the presence of the species. Based on the metrics analysis, habitat patches in the Zaryab wildlife refugee have higher connectivity and less margin than Bahr Asman protected area. In Bahr Asman protected area, the number of habitat patches is greater than the Zaryab wildlife refugee and the lack of continuity among the patches has led to the edge formation. The results of this study can be used for the habitat management of this critically endangered subspecies   http://ijae.iut.ac.ir/article-1-1118-en.pdf 2022-09-01 1 17 10.47176/ijae.11.2.10803 Black bear Ensemble model Habitat patches Landscape Kerman province H. Ruholamininejad hoseen.roholamini@gmail.com 1 Ardakan University AUTHOR M Morovati mymorovati@ardakan.ac.ir 2 Ardakan University AUTHOR P. Karami peymankarami1988@gmail.com 3 Malayer University AUTHOR
OTHERS_CITABLE پیش‌بینی تغییرات تولید خالص اولیه در مناطق مختلف رویشی ایران در دوره زمانی 2000 تا 2016 با استفاده از مدل‌های سری زمانی پوشش گیاهی از اجزاء مهم اکوسیستم‌های خاکی است که به صورت فصلی تغییر می‌کند و پارامترسازی دقیق پویایی پوشش-گیاهی با توسعه الگوهای دوره‌ای آن می‌تواند درک ما را از تعاملات پوشش گیاهی-اقلیمی تقویت کند. پژوهش حاضر با هدف بررسی و مدلسازی تغییرات تولید خالص اولیه در طول زمان در برخی از مناطق رویشی کشور شامل مناطق خزری، بلوچی، نیمه بیابانی، استپی معتدل، نیمه استپی گرم و جنگل‌های خشک و مقایسه رفتار تصادفی این مناطق با یکدیگر انجام گرفت. در پژوهش حاضر جهت بررسی تغییرات تولید خالص اولیه از تولیدات تصاویر سنجنده مودیس با کد MOD17A2 استفاده گردید. مدلسازی با استفاده از مدل سری زمانی (Seasonal Auto Regressive Intergrated Moving Average , SARIMA)  انجام گردید. بررسی توابع (Autocorrelation Function, AFC) و (Partial Autocorrelation Function, PACF) در مناطق مورد مطالعه نشان داد که مدل‌های سری‌ زمانی این مناطق، ایستا با ویژگی فصلی بودن در دوره های 12 ماهه بودند. پوشش گیاهی در منطقه خزری پایدارتر بود که نشان دهنده یک محیط پایدار با کمترین انحراف در تغییرات آب، نور و مواد غذایی می باشد. همچنین مشخص گردید بیشتر مناطق رویشی ایران را می‌توان با SARIMA مدلسازی و تغییرات آن را تا حد قابل اطمینانی پیش‌بینی کرد. مدل‌های برآورد شده برای مناطق خزری 0/83 = (Mean Relative Absolute Erorr, MRAE )،R2  =0/87 و 0/12 = (Root-Mean-Square Erorr, RMSE) و نیمه بیابانی با (0/12=RMSE 0/95=R2, 0/048=MARE) نسبت به سایر مناطق مدل های مناسب تری بودند.        http://ijae.iut.ac.ir/article-1-1119-fa.pdf 2022-09-01 19 35 10.47176/ijae.11.2.4332 توابع خود همبستگی و خودهمبستگی جزئی مدلسازی باکس و جنکینز تغییرات فصلی مدل فصلی خودهمبسته‌ی میانگین-متحرک Prediction of Net Primary Production Changes in Different Phytogeographical Regions of Iran from 2000 to 2016, Using Time Series Models Vegetation cover is an important component of terrestrial ecosystems that changes seasonally. Accurate parameterization of vegetation cover dynamics through developing indicators of periodic patterns can assist our understanding of vegetation-climate interactions. The current study was conducted to investigate and model vegetation changes in some phytogeographical regions of Iran including, Khazari, Baluchi, semi-desert, temperate steppe, warm semi-steppe and arid forest and to compare their stochastic behavior. To study the vegetation changes the net primary production (NPP) was used, based on the products of Moderate Resolution Imaging Spectroradiometer (MODIS) sensor (MOD17A2 series). Seasonal Auto Regressive Integrated Moving Average (SARIMA) time series model was used for modeling NPP. The Autocorrelation Function (ACF) and Partial Autocorrelation Function (PACF) of time series showed that these areas were static with seasonality in 12-month periods. It also showed that the vegetation in Khazari region was more stable, which indicates a stable environmental condition with the least deviation in water, light and nutrients. We also found that most of the vegetative regions of Iran can be modeled with SARIMA and its changes can be reliably predicted. Estimated models for Khazari (Root-Mean-Square Error, (RMSE) = 0.12, R2 = 0.87, Mean Relative Absolute Error (MARE) = 0.083) and semi-desert (RMSE = 0.12, R2 = 0.95, MARE = 0.048) were more suitable models than other regions.   http://ijae.iut.ac.ir/article-1-1119-en.pdf 2022-09-01 19 35 10.47176/ijae.11.2.4332 The autocorrelation function Partial autocorrelation function Box-Jenkins model Seasonality Seasonal autoregressive integrated moving average (SARIMA) F. Sayedzadeh F.sayedzadeH@na.iut.ac.ir 1 Isfahan University of technology AUTHOR S. Soltani ssoltani@iut.ac.ir 2 Isfahan University of technology AUTHOR R. Modarres reza.modarres@iut.ac.ir 3 Isfahan University of technology AUTHOR
OTHERS_CITABLE مدلسازی آشیان بوم‌شناختی افعی‌های کوهستانی تبار راده‌ای (Raddei) در ایران، قفقاز و شرق ترکیه افعی‌های کوهستانی جنس مونتی‌ویپرا (Montivipera) به‌طور کلی، و گونه‌های تبار Raddei به‌طور خاص، نمونه‌های ویژه‌ای از بوم‌زادی جدید (Neo-endemism) در فلات ایران، آناتولی و قفقاز هستند. با توجه به چالش‌های پیش روی حفاظت از گونه‌های مذکور، تعیین زیستگاه‌های مطلوب آنها، به منظور اتخاذ رویکردهای حفاظتی مناسب، ضروری است. به این منظور آشیان بوم شناختی هر گونه بر اساس چهار الگوریتم مدل تعمیم یافته خطی، مدل تعمیم یافته ارتقایی، جنگل تصادفی و بیشینه بی‌نظمی مدلسازی و در قالب یک مدل اجماعی بررسی شد. همچنین، با استفاده از یک رویکرد تجزیه به مؤلفه‌های اصلی، جدایی آشیان بوم‌شناختی گونه‌ها مورد بررسی قرار گرفت. نتایج نشان داد که مدل پراکنش گونه‌ها با (Area Under the Curve, AUC) و (True Skills Statistics, TSS) بالاتر از 0/9 عملکرد پیش‌بینی بالایی دارند. گونه‌های تبار Raddei الگوهای متفاوتی از اشغال آشیان بوم‌شناختی را نشان داده و بیشترین تمایز در گرادیان متغیرهای تغییرات فصلی درجه حرارت، بارش سالانه و تنوع ناهمواری‌ها دیده شد. دو گونه افعی لطیفی و افعی کوهرنگی، با وجود فاصله جغرافیایی زیاد، بیشترین میزان مشابهت و هم‌پوشی آشیان بوم‌شناختی را نشان دادند. با توجه به پراکنش و انعطاف‌پذیری تکاملی محدود افعی‌های کوهستانی در مناطق کوهستانی البرز و زاگرس، طرح‌ریزی اقدامات مدیریتی جهت کاهش عوامل تهدید کننده بقای طولانی مدت این گونه‌ها امری ضروری است. http://ijae.iut.ac.ir/article-1-1126-fa.pdf 2022-09-01 37 49 10.47176/ijae.11.2.8531 مدل‌سازی مطلوبیت زیستگاه نگهداشت آشیان بوم‌شناختی انعطاف‌پذیری تکاملی طرح‌ریزی حفاظت Ecological Niche Modeling of Mountain Vipers from the Raddei Clade in Iran, Caucasus and Eastern Turkey Mountain vipers of the genus Montivipera, generally, and the species of the Raddei clade, specifically, are interesting examples of species neo-endemism in Iran, Anatolia, and the Caucasus. Given the critical conservation status of these species, it is necessary to identify their suitable habitats for prioritizing conservation measures. We modeled ecological niche of each species based on four species distribution models, including generalized linear models, generalized boosting models, random forest, and maximum entropy and combined them in an ensemble model. Also, using a new principal component analysis (PCA-env), the ecological niche divergence of the species was investigated. Results indicated that all models with AUC and TSS > 0.9 had excellent predictive performance. The species of the Raddei clade showed different patterns of ecological niche occupation and the greatest differentiation was seen in the gradient of temperature seasonality, annual precipitation, and topographical ruggedness. M. latifii and M. kuhrangica, despite the great geographical distance, revealed the highest degree of niche overlap and niche similarity. Due to the restricted distribution and limited evolutionary adaptability of mountain vipers in the mountainous regions of Alborz and Zagros, it is imperative to plan conservation measures to reduce the factors that threaten the long-term survival of these species.   http://ijae.iut.ac.ir/article-1-1126-en.pdf 2022-09-01 37 49 10.47176/ijae.11.2.8531 Habitat suitability modeling Niche conservatism Evolutionary flexibility Conservation planning M. Ahmadi mahmadi@iut.ac.ir 1 Isfahan University of Technology AUTHOR M. R. Hemami mrhemami@iut.ac.ir 2 Isfahan University of Technology AUTHOR M. Kaboli mkaboli@ut.ac.ir 3 University of Tehran AUTHOR M. Malekian mmalekian@iut.ac.ir 4 Isfahan University of Technology AUTHOR
OTHERS_CITABLE ارزیابی تغییرات کاربری اراضی در مناطق با ریسک بالای تعارض خرس قهوه‌ای در استان فارس ارزیابی تغییرات صورت گرفته در زیستگاه گامی مهم در اولویت بندی مکانی برنامه های کاهش تعارض انسان و حیات وحش است. پژوهش حاضر با هدف ارزیابی تغییرات کاربری اراضی در محدوده مناطق با ریسک بالای تعارض خرس قهوه ای در استان فارس اجرا شد. در گام نخست، با استفاده از نقاط وقوع تعارض، متغیرهای پیش بینی کننده و رویکرد اجماع مدل های کوچک، نقشه مناطق با ریسک بالای تعارض تهیه شد. در گام دوم، با استفاده از روش‌ های سنجش از دور، روند تغییرات کاربری اراضی در بازه زمانی 30 ساله در محدوده مناطق با ریسک بالا ارزیابی شد. نتایج مدل سازی ریسک نشان داد که فاصله از مناطق روستایی، کریدورهای مهاجرتی و تراکم لکه های جنگلی مهم ترین متغیرها در احتمال وقوع ریسک تعارض خرس است. به ترتیب 3/75 و 6/91 درصد از منطقه در طبقات با خطر بالا و متوسط قرار گرفت. ارزیابی تغییرات کاربری نشان داد در یک بازه زمانی 30 ساله سطح باغات و اراضی کشاورزی از 12167 به 52662 هکتار افزایش یافته است. چنین تغییری می تواند احتمال تعارض خرس با جوامع بومی را افزایش دهد. بر اساس نتایج این پژوهش طرح ریزی برنامه ریزی های بین سازمانی به منظور جلوگیری از روند تخریب عرصه های طبیعی جهت کاهش تعارضات انسان و حیات وحش امری ضروری است. http://ijae.iut.ac.ir/article-1-1137-fa.pdf 2022-09-01 51 64 10.47176/ijae.11.2.833 تعارض خرس قهوه‌ای-انسان تغییرات کاربری اراضی مدل‌سازی ریسک نقاط داغ تعارض Assessing Land Use Changes in Areas with High Risk of Human-Brown Bear Conflict in Fars Province Assessing habitat changes is an important step in the spatial prioritization of management efforts, aimed at reducing conflicts. We assessed landuse/cover change in areas with high risk of human-bear conflict in Fars province. In the first step, we predicted the conflict hotspots, using bear damage incidents, a suit of predictors, and the Ensembles of Small Models (ESMs) approach. In the second step, we assessed the trend of landuse/cove changes in a 30-years period in the areas with medium to high risk of conflict, using remote sensing techniques. Results of conflict risk modeling showed that proportion of suitable habitats, distance to village, density of forest patches, and corridor bottlenecks were the main predictors, contributing to bear damaging risk. A total of 3.75 and 6.91% of the landscape were identified as the areas with high and medium risk, respectively. Assessment of landuse/cove changes showed that in a period of 30 years, the extent of croplands and orchards has increased from 12,167 to 52,662 hectares. Such a substanitial landuse/cover changes can increase the risk of bear damages. The obtained results emphasize that inter-organizational planning is an emergency effort in mitigating human-bear conflicts. http://ijae.iut.ac.ir/article-1-1137-en.pdf 2022-09-01 51 64 10.47176/ijae.11.2.833 Human-brown bear conflict Landuse changes Risk modeling Conflict hotspots R. Khosravi r-khosravi@shirazu.ac.ir 1 Shiraz University AUTHOR H. R. Pourghasemi hm_porghasemi@yahoo.com 2 Shiraz University AUTHOR Y. Moveseghi y.movaseghi76@gmail.com 3 Shiraz University AUTHOR
OTHERS_CITABLE تغییرات مکانی-زمانی اختلال پوشش گیاهی طبیعی در حوضه اهل ایمان، استان اردبیل تحقیق حاضر با هدف ارزیابی تغییرات مکانی-زمانی وضعیت اختلال در پوشش گیاهی طبیعی حوضه اهل ایمان واقع در استان اردبیل برنامه‌ریزی شده است. بدین منظور، ابتدا نقشه کاربری اراضی سه سال 2000، 2010 و 2020 از تصاویر ماهواره لندست استخراج شد. سپس هفت سنجه سیمای سرزمین (تراکم لکه، تراکم حاشیه، شاخص تکه‌شدگی، فاصله اقلیدسی نزدیک‌ترین همسایه، شاخص سرایت، نسبت محیط به مساحت و غنای لکه)، شاخص پوشش گیاهی تفاضلی نرمال شده (Normalized Difference Vegetation Index, NDVI) و تراکم جاده محاسبه شد. در نهایت شاخص اختلال (Disturbance Index, DI) از طریق مجموع حاصل‌ضرب مقادیر معیارها در وزن آن‌ها در 11 زیرحوضه محاسبه شد. شاخص اختلال (DI) در پنج طبقه خیلی کم (87–0)، کم (163–88)، متوسط (239–164)، زیاد (315-240) و خیلی زیاد (316<) طبقه‌بندی شد. بر اساس نتایج مشخص شد که از نظر کلی حوضه در سال‌های 2000، 2010 و 2020 به‌ترتیب دارای DI برابر با 77/177، 17/95 و 07/135 بوده که ‌به‌ترتیب نشان‌دهنده وضعیت اختلال متوسط، کم و کم است. همان‌گونه که مشاهده می‌شود بین سه سال مورد بررسی، از نظر وضعیت اختلال تفاوت معنی‌داری وجود ندارد. هم‌چنین، به‌طور کلی در بخش‌های شرق و شمال حوضه وضعیت بدتری از اختلال نسبت به بخش‌های جنوبی و مرکزی مشاهده شد.  http://ijae.iut.ac.ir/article-1-1106-fa.pdf 2022-09-01 65 82 10.47176/ijae.11.2.13801 پیوستگی سیمای سرزمین تخریب بوم‌شناختی رویکرد شاخص‌محور مدیریت سرزمین Spatio-Temporal Changes of Natural Vegetation Disturbance in the Ahle Iman Watershed, Ardabil Province The present study aimed to assess the spatio-temporal changes in the natural vegetation cover of Ahl Iman watershed, Ardabil province. For this purpose, land use maps of the three years (2000, 2010, and 2020) were extracted from Landsat satellite images. Then, seven landscape metrics (patch density, edge density, patch richness, splitting index, contagion index, Euclidean nearest neighbor distance, and mean perimeter–area ratio), Normalized Differential Vegetation Index (NDVI), and road density were calculated. Finally, the Disturbance Index (DI) was calculated by the sum of the values of the criteria multiplied by their weight in the 11 sub-watersheds. The disturbance index (DI) was classified into very low (0-87), low (88-163), medium (164-239), high (240-315), and very high (>316) categories. Results showed DI of 177.77, 95.17, and 135.07 for the three studied years (2000, 2010, and 2020), representing a moderate, low, and low disturbance respectively. There was no significant difference between the three years in terms of disturbance. Further, the eastern and northern parts had higher disturbance indices compared to the southern and central parts of the watershed. http://ijae.iut.ac.ir/article-1-1106-en.pdf 2022-09-01 65 82 10.47176/ijae.11.2.13801 Ecological degradation Index-based approach Land management Landscape integrity A. Hemmatzadeh abolfazlvet@gmail.com 1 University of Mohaghegh Ardabili AUTHOR A. Esmali Ouri esmaliouri@uma.ac.ir 2 University of Mohaghegh Ardabili AUTHOR R. Mostafazadeh raoofmostafazadeh@uma.ac.ir 3 University of Mohaghegh Ardabili AUTHOR M. Golshan golshan.mohammad@yahoo.com 4 University of Agricultural Sciences and Natural Resources Sari AUTHOR Z. Hazbavi z.hazbavi@uma.ac.ir 5 University of Mohaghegh Ardabili AUTHOR N. Alaei nazila.alaie96@gmail.com 6 Urmia University AUTHOR
OTHERS_CITABLE ارزیابی سلامت اکولوژیک رسوبات سطحی سواحل مرجانی خارک و خارکو (خلیج فارس، ایران) به منظور ارزیابی ریسک اکولوژیک زیستگاه مرجانی جزایر خارک و خارکو در خلیج فارس (استان بوشهر)، در شهریور 1399 نمونه-های رسوب سطحی از هفت ایستگاه جمع آوری شد. این جزایر، به دلیل داشتن آبسنگ های مرجانی دارای ارزش اکولوژیک فوق العاده ای است. میزان عناصر بالقوه سمی، دانه بندی، مواد آلی کل، فسفرکل و نیتروژن کل در رسوبات، به ترتیب با استفاده از طیف‌سنجی جرمی پلاسمای جفت‌شده القایی، الک، روش سوزاندن در کوره، اسپکتروفتومتر و کجلدال اندازه¬گیری شد. میانگین غلظت آلومینیوم و آهن (درصد) و میانگین غلظت نیکل، سرب، روی، وانادیوم، فسفر کل و نیتروژن کل (میلی گرم در کیلوگرم) در رسوبات به ترتیب 0/53±0/76، 55/35±0/0، 19±35، 1/5±2/1، 10±22، 25±40، 7/3±0/0 و 6/2±14/7 ثبت گردید. بر اساس میزان فاکتور غنی شدگی نیکل و روی (به ترتیب در حد متوسط تا خیلی شدید)، منشاء آنها در اطراف خارک می تواند فعالیت انسانی از جمله صنایع نفتی باشد. شاخص بار آلودگی (0/25- 0/06) کلیه ایستگاه ها را بدون آلودگی نشان داد. میزان نیکل در ایستگاه 1 و 7 به ترتیب از استانداردهای"محدوده اثرات متوسط" و "سطح اثرات احتمالی" بیشتر بود که نشان دهنده اثرات زیستی متوسط و احتمالی این فلز بر روی موجودات بستر است. به طور کلی کیفیت اکولوژیک رسوبات سطحی اطراف خارکو نسبت به خارک، بهتر ارزیابی گردید. http://ijae.iut.ac.ir/article-1-1141-fa.pdf 2022-09-01 83 101 10.47176/ijae.11.2.14221 رسوب جزایر خارک و خارکو عناصر بالقوه سمی فسفر کل نیتروژن کل Ecological Health Assessment of the Surface Sediments of the Coral Reefs of Khark and Kharko Islands (Persian Gulf, Iran) To assess the ecological risk of the coral reef habitats of Khark and Kharko islands in the Persian Gulf, (Bushehr province), the surface sediment samples were collected from seven stations, in September 2019. The islands have a great ecological value due to the presence of coral reefs. The amounts of potentially toxic elements, sediment texture, total organic matter, total phosphorus and total nitrogen in the sediments were measured by inductively-coupled plasma mass spectrometry, sieve analysis, furnace burning method, spectrophotometer, and Kjeldahl, respectively. The mean concentrations of Al, Fe (%) and Ni, Pb, Zn, V, TP and TN (mg/kg) in the sediments were recorded 0.76±0.53, 0.55±0.35, 35±19, 2.1±1.5, 22±10, 40±25, 0.7±0.3 and 14.7±6.2, respectively. Based on the amount of Ni and Zn enrichment factors (moderate to very sever, respectively), their source around Khark Island could be related to human activities, such as oil industries. Pollution load index (0.06-0.25) showed all stations without pollution. The amount of Ni in stations 1 and 7 were higher than the "range of moderate effect" and "level of possible effects", which indicates the possible biological effects of this element on the benthic organisms. In general, the ecological quality of surface sediments around Kharko was better than Khark Island.   http://ijae.iut.ac.ir/article-1-1141-en.pdf 2022-09-01 83 101 10.47176/ijae.11.2.14221 Sediment Khark and Kharko Isalnds Potentially toxic elements Total phosphorus Total nitrogen Gh. Ghorbanzadeh Zafarani ghorbanzadeh110@yahoo.com 1 Research Center for Environment and Sustainable Development (RCESD), DOE AUTHOR F. Hosseini Tayefeh farhadtayefeh@gmail.com 2 Research Center for Environment and Sustainable Development (RCESD), DOE AUTHOR A. Ahmadi Givi amin.ahmadi.givi@ut.ac.ir 3 University of Tehran AUTHOR M. A. Tollab gh2012eco@gmail.com 4 Department of Environment (DOE), Bushehr, Iran AUTHOR A. Saber gh.ghorbanzadeh@rcesd.ac.ir 5 Department of Environment (DOE), Bushehr, Iran AUTHOR