فهرست مطالب

ژئوفیزیک ایران - سال دهم شماره 4 (پیاپی 32، زمستان 1395)

مجله ژئوفیزیک ایران
سال دهم شماره 4 (پیاپی 32، زمستان 1395)

  • تاریخ انتشار: 1395/12/07
  • تعداد عناوین: 10
|
  • فرهنگ احمدی گیوی *، مرجانه حقیقی، عباسعلی علی اکبری بیدختی صفحات 1-11
    مطالعه اثرات دینامیکی اصطکاک بستر بر حرکت توده آب مستلزم آن است که ویژگی های لایه مرزی موسوم به لایه اکمن شناسایی شود. لایه اکمن در بستر شیب دار شامل جریان هایی به نام شارش اکمن به صورت فراشیب یا فرو شیب است که این جریان ها می توانند سبب افزایش آمیختگی و تغییر در ضخامت لایه مرزی اکمن شوند. هم چنین نیروی اصطکاک در بستر سبب کاهش سرعت جریان و به دنبال آن کاهش نیروی کوریولیس شده و در نتیجه جریان زمینگرد از بین می رود و فرایند چرخ کاهی رخ می دهد. در این پژوهش به منظور شناخت ویژگی های لایه اکمن در بستر شیب دار غرب تنگه هرمز به بررسی تاثیر شیب بستر بر توزیع پارامتر های فیزیکی آب و پدیده های حاصل از آن پرداخته شده است. برای انجام این پژوهش، با استفاده از داده های میدانیCTD مربوط به غرب تنگه هرمز در سال های 1384و 1385، موجود در موسسه ملی اقیانوس شناسی، و به کمک نرم افزارهای مرتبط، نیم رخ های تغییرات دما، شوری، چگالی و سرعت های کشندی با عمق و هم چنین مقاطع عرضی چگالی ترسیم شدند. بررسی این نیم رخ ها نشان می دهد که در بستر تنگه هرمز، جریان فراشیب است و پربندهای هم چگالی با افزایش شیب بستر تقریبا به صورت عمود بر شیب قرار می گیرند. کم ترین زمان چرخ کاهی شارش اکمن دو ساعت محاسبه شده و مربوط به ایستگاهی است که بیشترین شیب را دارد. به علاوه، هرچه شیب بستر بیشتر و چینه بندی چگالی قوی تر باشد، چرخ کاهی شارش اکمن سریع تر اتفاق می افتد. بیشینه ضریب پخش پیچکی در بین ایستگاه های مورد مطالعه حدود m2s-1 005/0 به دست آمده است و در ایستگاهی که ضریب پخش پیچکی آن از همه بزرگ تر است، عدد اکمن نیز به طور قابل ملاحظه ای از سایر ایستگاه ها بزرگ تر است که نشان دهنده جریان ناپایا و آمیختگی ناشی از چینش لایه اکمن می باشد.
    کلیدواژگان: لایه اکمن، تنگه هرمز، عدد اکمن، ضریب پخش پیچکی، بستر شیب دار، پارامتر های فیزیکی
  • پرویز ایران نژاد*، فرهنگ احمدی گیوی، نسیم نیکویی صفحات 12-27
    در این پژوهش به بررسی تغییرپذیری زمانی و مکانی و بی هنجاری دمای هوا در ارتفاع دو متر از سطح زمین در فصل زمستان در سه مقیاس شامل کل ایران، شش منطقه مشخص شده بر اساس طول و عرض جغرافیایی و در مقیاس یاخته شبکه داده های بازتحلیل NCEP/NCAR به صورت میانگین ماهانه در سال های 1950 تا 2010 پرداخته شده است. نخست داده های بازتحلیل شده دما از نظر کیفیت، بررسی و سال های پیش از 1967 به دلیل کیفیت پایین از مطالعه حذف شد. بررسی ها نشان داد که تغییرات میانگین دمای زمستان ایران سنجه مناسبی برای بیان تغییرات بین سالانه میانگین زمستانی دمای مناطق مختلف است. از 44 سال مطالعه شده، بزرگ ترین بی هنجاری های مثبت دما در زمستان های 1970، 1979، 1999 و 2010 و بزرگ ترین بی هنجاری های منفی در زمستان های 1972، 1973، 1989 و 2008 رخ داده است. بررسی میانگین الگوی بی هنجاری میدان های ارتفاع ژئوپتانسیلی، دما و باد در زمستان های بی هنجار گرم و سرد در حوزه بین 10 تا 60 درجه شمالی و 10 درجه غربی تا 80 درجه شرقی نشان داد که مراکز بیشینه بی هنجاری این کمیت ها در سال های گرم و سرد ایران کاملا عکس یکدیگرند. دو مرکز بیشینه بی هنجاری ارتفاع ژئوپتانسیلی در سال های بی هنجار سرد و گرم ایران، یکی در منطقه شمال غرب اروپا و دیگری در شمال شرق دریای کاسپین قرار دارند. بر اساس تفاوت بی هنجاری ارتفاع ژئوپتانسیلی تراز 250 هکتوپاسکال در این دو مرکز، شاخصی برای بیان بی هنجاری دمای زمستان در ایران معرفی شد. تحلیل همبستگی نشان داد که همبستگی مثبت معنی دار بین میانگین دمای زمستان ایران و هم چنین میانگین دمای زمستان مناطق مختلف کشور و شاخص معرفی شده وجود دارد. به دلیل شباهت الگوی دورپیوند به دست آمده در این پژوهش با دو دورپیوند EA/WR و NCP همبستگی بین شاخص های این دو دورپیوند و دمای زمستان ایران تحلیل شد. نتایج نشان دهنده توانایی بالاتر شاخص معرفی شده نسبت به دو شاخص دیگر برای بیان دمای زمستان ایران است.
    کلیدواژگان: ایران، اروپا، دریای کاسپین، بی هنجاری دما، بی هنجاری ارتفاع ژئوپتانسیلی، دورپیوند، EA، WR، NCP
  • حجت الله دانشمند*، پیمان محمودی صفحات 28-47
    هدف اصلی این مطالعه آشکارسازی رفتار دوره ای خشکسالی های ماهانه ایران است. برای دستیابی به این هدف از آمار بارش ماهانه 41 ایستگاه همدیدی کشور برای یک دوره آماری 31 ساله (2013–1983) و برای شناسایی فراوانی درجات مختلف خشکسالی ها از شاخص خشکسالی موثر (EDI) استفاده شد. هدف اصلی این شاخص اختصاص یک ارزش عددی به هر رویداد بارندگی در مقیاس زمانی ماهانه می باشد. در تحلیل طیفی و شناسایی دوره های غالب موجود در سری های زمانی خشکسالی ایستگاه های مورد مطالعه از تبدیل فوریه سریع (FFT) تقسیمات بسامدی با مبنای دو استفاده شد. نتایج نشان داد که دوره های غالب موجود در سری های زمانی خشکسالی های ایران بسیار متنوع بوده و از دوره های 2 تا 30 ساله را شامل می شوند. به طور کلی برای دوره های بلندمدت، می توان در کل ایران دو دوره 10 ساله و 30 ساله را به عنوان دوره های غالب در نظر گرفت. در دوره های بلند مدت، فراوانی این دو دوره در مولفه اول بیشینه است و در دو مولفه دیگر کاهش می یابد. اما درخصوص دوره های کوتاه مدت، نمی توان برای کل ایران یک دوره غالب مشخص کرد. نکته حائز اهمیت در مورد دوره های کوتاه مدت عدم وجود دوره های درون سالی، به استثنای یک ایستگاه (ایستگاه بابلسر)، در کل ایستگاه های مورد مطالعه است.
    کلیدواژگان: خشکسالی، شاخص خشکسالی موثر، تحلیل طیفی، توان، بسامد
  • فروغ کلوندی، سید خلیل متقی*، اسماعیل شبانیان صفحات 48-61
    مطالعه ساختار سرعتی هر منطقه کمک شایانی به شناخت وضعیت لرزه زمین ساختی آن می کند. تعیین پارامترهایی چون ضخامت پوسته و سنگ کره، مکان یابی دقیق زمین لرزه ها و پیش بینی تحولات آینده زمین ساختی هر منطقه در گرو داشتن مدل سرعتی موثق از منطقه است. فلات ایران گستره وسیعی بین دو صفحه عربی و اوراسیا است و دگرشکلی آن ناشی از همگرایی این دو صفحه است. در این میان، کمربند کوهزایی زاگرس یکی از فعال ترین و جوان ترین کوهزادهای موجود در مرز های همگرایی قاره ای در جهان است که شناخت ویژگی های ساختاری این منطقه می تواند به فهم نحوه تغییر شکل در مرحله اول برخورد قاره ای و کوهزایی کمک کند. در این تحقیق ساختار سنگ کره در پهنه برخوردی زاگرس شمالی، با استفاده از وارون سازی هم زمان توابع گیرنده و منحنی های پاشش امواج سطحی مورد مطالعه قرار می گیرد. منحنی های پاشش از مطالعه ای که در سال 2014 از طریق برش نگاری (توموگرافی) رحیمی و همکاران (2014) انجام دادند، برگرفته و توابع گیرنده از 161 دورلرزه با فاصله رومرکزی °30 تا °95 و بزرگای بیش از 5 محاسبه شده است. دورلرزه ها را 38 ایستگاه لرزه نگاری موقت در امتداد پروفایلی به طول 400 ~ کیلومتر در حد فاصل شهرهای ایلام تا قم ثبت کرده اند. در این مطالعه با استفاده از فرایند وارون سازی هم زمان توابع گیرنده و منحنی های پاشش، مدل های سرعتی یک بعدی موج برشی در زیر هر ایستگاه به دست می آید. سپس، از کنار هم قرار دادن این مدل های یک بعدی، یک مدل دوبعدی سرعت برای سنگ کره زیرین ایستگاه های شبکه حاصل می شود. در این مدل یک زبانه کم سرعت در درون پوسته در مدل سرعتی قابل تشخیص است که از محل گسل اصلی عهد حاضر شروع می شود و با فاصله گرفتن از این گسل در راستای شمال شرقی به عمق های بیشتر می رود. این بی هنجاری در فاصله 200 کیلومتری شمال شرق گسل در عمق 35 ~ کیلومتر قابل مشاهده است. مدل سرعتی نشان می دهد که مقدار ضخامت پوسته از ابتدای خط اندازه گیری در بخش های جنوبی زاگرس (محدوده ایلام)، 43 کیلومتر است که با حرکت به سمت شمال شرق به میزان 57 کیلومتر در زیر گسل اصلی عهد حاضر می رسد. در سنندج - سیرجان و ارومیه - دختر به تدریج بر این ضخامت افزوده شده و در مرز این دو ناحیه به بیشینه مقدار 62 کیلومتر می رسد. در زیر ایران مرکزی میزان ضخامت پوسته کاهش یافته و در انتهای خط اندازه گیری در زیر ایران مرکزی این ضخامت به 42 کیلومتر می رسد. مدل سرعتی به دست آمده، همچنین، اطلاعات خوبی از بخش سنگ کره ای گوشته ارائه می دهد. این مدل نشان می دهد که سنگ کره پرسرعت زیرین زاگرس شمالی در زیر سنندج - سیرجان و ارومیه دختر و لبه جنوبی ایران مرکزی گسترش یافته است. چنین مشاهده ای در تطابق با زیرراندگی مشاهده شده در پوسته است و می تواند شاهدی بر زیرراندگی بلوک عربی به زیر ایران مرکزی باشد.
    کلیدواژگان: توابع گیرنده، منحنی های پاشش موج ریلی، وارون سازی هم زمان، زیرراندگی، زاگرس شمالی
  • سیده صفورا صدیق مروستی، آناند نانادسیکان، عباسعلی علی اکبری بیدختی*، سرمد قادر، محمد صدیق مرتضوی صفحات 62-74
    در سال های اخیر، افزایشی در شکوفایی کشند قرمز در شمال غربی دریای عرب و دریای عمان مشاهده شده است که این سوال را در ذهن به وجود می آورد: آیا تغییرات اقلیمی باعث افزایش این روند شده است یا خیر؟ با هدف پاسخ به این سوال، در این پژوهش از مدل های مختلف سامانه زمین و داده های ماهواره ای استفاده شده است. داده های ماهواره ای دو شکوفایی را در این منطقه نشان می دهند، شکوفایی زمستانی که بیشینه آن در ماه فوریه و شکوفایی تابستانی که بیشینه آن در ماه سپتامبر است. تغییرات درون سالانه ای نیز در شکوفایی زمستانی حاصل از اثر پیچک های چرخندی که توزیع مکانی متفاوتی از یک سال به سال دیگر دارند، مشخص شده است. دو مدل با تفکیک کم (°1) با بخش زیست زمین شیمیایی تقریبا پیچیده (مدل TOPAZ) چرخه سالانه را نشان می دهند ولی قادر به نمایش پیچک ها و تغییرات درون سالانه نیستند. مدل های با تفکیک بیشتر (GFDL CM2.6) که توانایی مدل سازی پیچک ها را دارند همراه با بخش زیست زمین شیمیایی ساده تر (مدل miniBLING) تغییرات درون سالانه بزرگ تری را نشان می دهند، اما مقدار شکوفایی زمستانی را بیش از اندازه پیش بینی می کنند. این مدل اگرچه رابطه بین شکوفایی و پیچک را در بخش جنوبی به خوبی نشان می دهد، اما در بخش شمالی منطقه موفق نیست. این امر به دلیل نداشتن توانایی در مدل سازی دماشیب (ترموکلاین) و غذاشیب (نوتری کلاین) قوی و درست در آن مناطق می باشد.
    کلیدواژگان: کشند قرمز، شکوفایی، مدل سازی عددی، پیچک، مدل سامانه زمین، زیست زمین شیمیایی
  • فریده کنعانی، مصطفی فضلی*، ویکتوریا عزتیان صفحات 75-84
    ازون سطح زمین در طی فرایندهای نورشیمیایی اکسیژن هوا در حضور نور خورشید و NOx به عنوان یکی از مهم ترین عوامل موثر، تولید می شود. در این پژوهش به بررسی ازون در شهر اصفهان در آذرماه سال 1388 پرداخته شده، که بنابر مشاهدات سازمان هواشناسی، در طی این بازه دودمه نورشیمیایی بر فراز اصفهان قابل مشاهده بود. همه بررسی ها در ساعت های تابش آفتاب (15-6 محلی) انجام گرفته است. در این پژوهش با استفاده از روش تقریب حالت پایا رابطه سرعت دیفرانسیل تابعی از سه متغیر NO،NO2 و O3 برای بررسی تغییرات غلظت ازون به دست آمده است. انرژی فعال سازی، مرتبه واکنش و سازوکار تشکیل واکنش ازون در وردسپهر مطالعه و محاسبه شده اند . NO و O3 سطح زمین به عنوان شاخصی برای تشکیل دودمه نورشیمیایی درنظر گرفته می شوند، لذا تولید دودمه بر مبنای تغییرات غلظت NO وO3 مورد توجه قرار گرفت. بررسی روابط تولید دودمه برمبنای غلظت ها و بررسی براساس زمان، نشان دهنده پیروی تولید دودمه از روابط غیرخطی و درجه دوم است. بنا به معادله آهنگ تغییر زمانی ازون، d[O3]/dt و نتایج انرژی فعال سازی می توان تولید ازون را در شهر اصفهان از واکنش هایی که دارای سینتیک کلی مرتبه سوم هستند، در نظر گرفت. از طرف دیگر افزایش غلظت NOx در سطح زمین در اثر واکنش های نورکافتی (فوتولیز) افزایش غلظت ازون را سبب گردیده است.
    کلیدواژگان: ازون وردسپهری، سرعت واکنش، انرژی فعال سازی، مدل نورشیمیایی، الگوریتم تولید دودمه
  • فاطمه اعظم پور، امین روشندل کاهو *، علیرضا احمدی فرد، علیرضا حیدریان صفحات 85-96
    گسل ها یکی از مهم ترین رویدادهای زمین شناسی هستند که در اکتشافات معدنی، هیدروکربنی و مطالعات زمین شناختی از اهمیت بسیار زیادی برخوردارند. گسل ها به منزله یک تله ساختاری هیدروکربنی می باشند و از این رو شناخت و تعیین محل دقیق آنها عامل مهمی در توصیف مخازن هیدروکربنی است. در داده های لرزه ای بازتابی، گسل ها به صورت ناپیوستگی در امتداد بازتابنده ها و تغییر دامنه مشاهده می شوند. به دلیل ماهیت داده های لرزه ای، شناسایی صفحه گسل در این داده ها، فرایندی دشوار است و تعیین محل دقیق آنها به تنهایی از روی داده لرزه ای به ندرت امکان پذیر است. نشانگرهای لرزه ای، ابزاری هستند که اطلاعات نهفته در داده های لرزه ای را آشکار می نمایند و امکان تفسیر کمی را در اختیار مفسر لرزه نگاری فراهم می آورند. در این مقاله، صافی سوبل دو بعدی و سه بعدی به منزله یک نشانگر لرزه ای برای شناسایی و تعیین موقعیت گسل ها در داده لرزه ای بازتابی معرفی می شود. این صافی، یکی از ابزارهای پرکاربرد در فرایند تشخیص لبه در پردازش تصاویر است که گرادیان تصویر را محاسبه می کند و به تغییرات دامنه حساس است. به منظور بررسی عملکرد این روش در شناسایی گسل ها، نشانگر سوبل بر روی داده های لرزه ای مصنوعی و واقعی در دو حالت دو بعدی و سه بعدی مورد بررسی قرار می گیرد و نتایج حاصل با نتایج نشانگر همدوسی مقایسه می شود. بررسی نتایج به دست آمده نشان می دهد که نشانگر سوبل می تواند به مثابه یک نشانگر مناسب برای شناسایی گسل ها در داده های لرزه ای بازتابی مورد استفاده قرار بگیرد.
    کلیدواژگان: لرزه نگاری بازتابی، گسل، تشخیص لبه، نشانگر سوبل
  • مجتبی نقوی، حبیب رحیمی *، علی مرادی صفحات 97-111
    در مطالعات مرتبط با زلزله شناسی مهندسی برای برآورد خطر لرزه ای در مناطق مختلف و همچنین تعیین دقیق بزرگا، شبیه سازی جنبش نیرومند زمین و مطالعه انرژی مخرب حاصل از زمین لرزه در حوزه نزدیک تا متوسط، تخمین ضریب کیفیت امواج برشی و فشاری نقشی اساسی دارد. بدین منظور با توجه به قرار گرفتن چندین شهر بزرگ در پهنه شمال غربی فلات ایران، ضریب کیفیت امواج فشاری و برشی با استفاده از داده های هفده ایستگاه لرزه نگاری با بیش از 13000 زمین لرزه ثبت شده در شبکه مرکز لرزه نگاری کشوری (IRSC) و شبکه ملی لرزه نگاری نوار پهن ایران (INSN) برآورد شد. برای سه گروه داده به ترتیب با فاصله رو مرکزی کمتر از 100 کیلومتر، از 100 تا 200 کیلومتر و صفر تا 200 کیلومتر بررسی تغییرات جذب امواج حجمی در نه نوار بسامدی، با بسامدهای مرکزی در 3، 5، 7، 10، 14، 20، 28، 38 و 47 هرتز انجام گرفت و برای هرکدام از ایستگاه ها برای بسامدهای مختلف به طور جداگانه ضریب کیفیت تخمین زده شد. در این مطالعه از روش کدای بهنجار شده برای برآورد ضریب جذب امواج فشاری و برشی حاصل از زمین لرزه به عنوان روشی قابل اعتماد با توجه به حذف شدن اثرات چشمه و ساختگاه در مدل سازی انجام گرفته، استفاده شد. برای پهنه شمال غربی ایران مقادیر وابستگی بسامدی ضریب کیفیت امواج فشاری و برشی در همه ایستگاه ها برآورد شد به طوری که مقدار متوسط آن برای محدوده شمال غرب فلات ایران از و پیروی می کند. با توجه به مقادیر پایین و لذا جذب بالای امواج برشی و فشاری در پهنه شمال غرب ایران دامنه امواج حین عبور از زمین به شدت تضعیف می شود که این اثر جذب امواج لرزه ای، باعث کاهش خسارات ناشی از زمین لرزه ها در فواصل مناسب از گسل ها در زمان وقوع زمین لرزه خواهد شد.
    کلیدواژگان: موج های برشی و فشاری، شمال غرب فلات ایران، جذب، ضریب کیفیت
  • مهدی راستگو، حبیب رحیمی*، حسین حمزه لو حمزه لو، وحید غلامی صفحات 112-127
    در ناحیه جنوب ‏شرق ایران که شامل بخش غربی زون فرورانش فعال مکران است به دلایل مختلفی همچون لرزه‏ خیزی اندک و نبود تراکم ایستگاه ‏های لرزه ‏نگاری دائم، مطالعات اندکی در رابطه با ساختار سرعت موج برشی انجام شده است. هدف از این پژوهش، برآورد ساختار سرعت موج برشی و تغییرات عمق موهو در زیر چهار ایستگاه لرزه‏ نگاری نوارپهن ZHSF، KRBR، CHBR و BNDS واقع در محدوده جنوب‏ شرق ایران است. برای این منظور از روش وارون‏ سازی هم زمان منحنی ‏های پاشش امواج سطحی حاصل از زمین ‏لرزه های محلی ناحیه ریگان، به همراه تابع گیرنده موج P حاصل از لرزه ‏نگاشت‏ های زمین‏ لرزه ‏های دورلرز ثبت شده در چهار ایستگاه لرزه‏ نگاری ذکر شده، استفاده می‏شود. عمق موهو به دست آمده از این روش برای چهار ایستگاه ZHSF، KRBR، CHBR و BNDS به ترتیب برابر است با 4±38، 6±46، 2±26 و 5±56 کیلومتر. کم عمق بودن موهو و نیز بالاتر بودن مقادیر سرعت موج برشی به دست آمده در زیر ایستگاه CHBR با فرورانش کم عمق لایه اقیانوسی سرعت بالای صفحه عربی به زیر حاشیه جنوبی مکران مطابقت دارد. در عین حال ضخامت زیاد پوسته در محدوده ایستگاه BNDS، در تطابق با توپوگرافی مکان این ایستگاه و نیز برخورد پوسته ضخیم قاره ای صفحه عربی با زاگرس در این ناحیه می‏ باشد.
    کلیدواژگان: جنوب شرق ایران، سرعت موج برشی، عمق موهو، پاشش امواج سطحی، تابع گیرنده موج P
  • خلیل متقی*، زهرا ضرونی زاده، عبدالرضا قدس صفحات 128-141
    در این مقاله با تلفیق داده های ثبت شده توسط 8 ایستگاه دائم شبکه تبریز متعلق به شبکه لرزه نگاری کشوری موسسه ژئوفیزیک دانشگاه تهران در بازه زمانی 1996 تا 2009 و داده های ثبت شده توسط 16 ایستگاه شبکه موقت لرزه نگاری دانشگاه تحصیلات تکمیلی علوم پایه زنجان در بازه زمانی سال های 2008 تا 2011 در شمال غرب ایران، 3514 نگاشت از 943 زمین لرزه با بزرگای بین 6/1 تا 2/5 تحلیل و شکل منحنی کاهندگی با استفاده از الگوریتمRobust Lowess به صورت گرافیکی ارائه شد. شکل منحنی کاهندگی نشان می دهد که تغییر ضریب گسترش هندسی در فاصله کانونی 70 کیلومتر اتفاق می افتد و سبب می شود تا یک مدل کاهندگی دو خطی برای برازش در بسامد های کم (کمتر از 5 هرتز) استفاده شود. ضریب گسترش هندسی تا فاصله 70 کیلومتر برای بسامد های مختلف به طور متوسط 28/0±1/1- به دست آمد. ضریب گسترش هندسی در فواصل بیش از 70 کیلومتر در بسامد های کمتر از 5 هرتز به طور متوسط 27/0±44/0- به دست آمد که در مقایسه با مقادیر گزارش شده برای سایر نواحی دنیا مقدار کوچکی است و نشان از وجود مرز ناپیوستگی موهو با تباین سرعتی ضعیف تر نسبت به آن مناطق دارد. ضریب گسترش هندسی در بسامدهای 5 هرتز و بالاتر، در فواصل کمتر از 70 کیلومتر و بیش از آن به هم شبیه است و مدل کاهندگی را از حالت دوخطی به حالت یک تکه تبدیل می کند. با استفاده از ضریب کاهندگی جذب ذاتی، تابع کیفیت موج برشی به صورت معرفی می شود.
    کلیدواژگان: رابطه کاهندگی، تابع کیفیت، شمال غرب ایران
|
  • Farhang Ahmadi Givi *, Marjaneh Haghighi, Abbasali Aliakbari Bidokhti Pages 1-11
    A thorough appreciation of the dynamic effects of bottom friction on water mass flow requires an understanding of the boundary layer characteristics, known as Ekman layer. Ekman layer on a sloping bed includes upslope or downslope flows that may intensify mixing or change the thickness of the Ekman boundary layer. The bottom friction can also reduce the fluid current, thereby reducing the Coriolis force and destroying the geostrophic current, named as the process of spin down.
    This study is conducted to investigate the impact of bottom slope on the distribution of physical parameters of water and the resulting phenomena, in order to identify the characteristics of the Ekman layer on the bottom slopes of the western Strait of Hormuz. The research utilizes the CTD (conductivity, temperature, and depth) field data of the western Strait of Hormuz in 2005 gathered by the National Institute of Oceanography of Iran. The vertical profiles of temperature, salinity, density, tidal currents and the horizontal cross sections of density were made by using Excel and MATLAB software.
    The results indicate that the current on the bottom of Strait of Hormuz is of an upslope type, and the pycnoclines become almost perpendicular to the slope with the increase of slope of the bottom. The minimum Ekman flux spin down is estimated of two hr and is related to a station with the maximum slope. Greater bottom slope and stronger stratification are found to remarkably shorten the shutdown timescale of the Ekman layer. We estimated the eddy diffusivity, , between zero and 0/005 m2s-1. A positive correlation was also found between eddy diffusion coefficient and Ekman layer parameters, which is indicative of the non-stationary flow and mixing due to Ekman layer stoppage.
    The horizontal scale that characterizes the dimension of a boundary current is the Rossby radius of deformation. Since baroclinic flow involves a number of internal modes, there will be a Rossby radius of deformation for each mode. The average Rossby radius of deformation in the south of the island for the first baroclinic mode is estimated at about 5.6 km, and buoyancy frequency in the Strait of Hormuz is 0.02 ( ). Spin down time ( ) of the stations located in the south of Qeshm for the first baroclinic and barotropic modes were also calculated. Results show that the greater the depth of the boundary layer on the bottom, the larger the spin down time for the first baroclinic mode. When Ekman layer is arrested in the barotropic mode, it seems that the spin down time for the first baroclinic mode becomes very large; thus we can conclude that the spin down does not happen in the arrested Ekman layer. Typical thickness of the Ekman layer in the Strait of Hormuz varies between 0.5 and 22 meters.
    Since the most important driving force of water up the slope is the tide, the tidal level changes about the long mean for the periods of ten years, one year, one month and one day are examined in the area. The results show that in the northern Strait of Hormuz, the maximum tidal range is about 4.73 meters and the average neap and spring tides are 2.15 and 3.53 meters, respectively. The calculated tidal currents and the related profiles for all the stations indicate that the maximum tidal speed is seen at the top of the bottom mixed layer. As the bottom is approached the speed of tidal current decreases remarkably and turbulence dominates which represent the state of the surface layer.
    Keywords: Ekman layer, Strait of Hormuz, Ekman number, eddy diffusion coefficient, sloping bottom
  • Parviz Irannejad*, Farhang Ahmadi Givi, Nasim Nikouei Pages 12-27
    Air temperature is one of the principal climate variables with important direct and indirect socio-economic impacts. Among different factors affecting the climate variability of different regions, the low-frequency atmospheric phenomena have attracted the attention of climatologists because of their role in the fluctuations of the climate system in time scales on the order of years and decades. Since these phenomena affect the atmospheric conditions of regions far from their sources, they are called teleconnection patterns. Teleconnection patterns cause large-scale changes in the atmospheric circulation and affect air temperature, precipitation, storm tracks and the position and intensity of jet streams. Therefore, they are of great significance for determining regional climate anomalies.
    In this study, anomalies in the mean winter surface (2 m level) temperature of Iran are investigated for the period 1950 to 2010. Winter of each year is defined as January and February of that year and December of the previous year. The study is conducted in three different scales, including Iran as the whole, six different regions of Iran separated according to three latitude and two longitude bands and the grid scale, using the monthly mean NCEP/NCAR reanalysis data. First, the reanalysis data of temperature were quality controlled, and years before 1967 were excluded from the study because of the poor quality of the data. Statistical analysis showed that there are significant positive correlations between the mean winter temperatures of different regions and between the mean winter temperature of each region and that averaged over Iran. The largest positive temperature anomalies in the 44 winters occurred in 1970, 1979, 1999 and 2010 and the largest negative anomalies occurred in the winters of 1972, 1973, 1989 and 2008.
    Analysis of the mean winter geopotential height fields at different pressure levels in a domain covering 10–60 °N and 10°W –80 °E showed that there is a center of positive (negative) height anomaly over North-West Europe (North-East Atlantic) and a center of negative (positive) height anomaly over North-East Caspian Sea when the mean winter surface air temperature of Iran is anomalously low (high). The correlation coefficient is higher for upper atmospheric levels, and stays significant at 95% confidence level down to the 700 hPa level. In this regard, a new teleconnection index is introduced in the paper, to study the mean winter temperature anomalies in Iran, calculated as the mean winter 250hPa geopotential height anomaly over 50–60°E, 40–50°N minus that over 0–10 °E, 45–55 °N. This resembles a teleconnection pattern which is similar to that of the East Atlantic/West Russia (EA/WR: Barnston and Livezey, 1987; Krichak et al., 2012) and the North Sea–Caspian Pattern (NCP; Kutil and Benaroch, 2002). The new index has a significant correlation with the EA/WR, and particularly with the NCP index. We found a significant positive correlation between the average winter temperature in different regions of Iran and the introduced index which was much higher than those between the mean winter temperature of Iran and both EA/WR and NCP teleconnection indices. Comparison between the temperature and geopotential height anomalies at different levels reveals that the middle troposphere temperature is anomalously high (low) in winters that geopotential height is positively (negatively) anomalous. Also, the wind field shows a positive (negative) anomaly over the Caspian Sea and a negative (positive) anomaly over Europe when the mean winter air temperature in Iran is higher (lower) than normal.
    Keywords: Iran, Europe, Caspian Sea, temperature anomaly, geopotential height anomaly, teleconnection, EA, WR, NCP
  • Hojatolah Daneshmand*, Peyman Mahmoudi Pages 28-47
    Identification of periodic droughts can provide a scientific tool for predicting the occurrence of this environmental risk. One of the most important methods that can be used for detection of periodic droughts is the spectral analysis or frequency-domain analysis. The purpose of this study is to detect the periodic behavior of Iran’s monthly droughts. To this aim, use has been made of the monthly rainfall statistics of 41 synoptic stations for a period of 31 years (1983–2013) obtained from the IR of Iran Meteorological Organization (IRIMO). First, to identify considerable differences in droughts, the “effective drought index” (EDI) was employed. Assigning a numerical value to each rainfall event on a monthly time scale is the main objective of the EDI in order to compare the areas with different climates based on it. Then the spectral analysis by Fourier transform was used to identify the dominant available periods in the stations’ drought time series. In this analysis, Radix-2 Fast Fourier Transform was used. Since Fourier transform leads to difficulties in the spectral analysis of non-stationary series, first the stationary condition of the EDI monthly time series was investigated for all of the 41 stations used in this study. In those stations with a non-stationary condition, the Box–Cox procedure was applied to make them stationary.
    The results were classified into five areas. The results showed that the dominant periods in the time series of Iran’s droughts vary greatly from periods of 2 to 30 years. In addition to the dominant short-term periods, droughts in the northwest of Iran have two dominant long-term periods of 10 and 30 years. Droughts in the southern coast of the Caspian Sea have shorter periodicities. So with the exception of the two stations of Ramsar and Ardebil with the periodicities of 30 and 15 years, the other stations have periods of less than 10 years. The dominant periods in most stations in the northern districts of central Iran are short-term ones, and the longest periodicity in this area is 15 years. The southern districts of the central Iran have a more consistent spectral behavior. In this area, the dominant period with the highest power is the period of 30 years. All of the stations in this area with the exception of Mashhad station have a 30-year period in their first three components. The 10-year period, as the second dominant period, can also be seen in half of the stations in this area. The southeast area does not exhibit any specific behavior for the dominant periods. The longest periodicity in Chabahar station is a 10- year period. The two stations in Zabol and Iranshahr have shown the longest period of 30 years and the long-term period seen in Zahedan is a 15-year one. In addition to these long-term periods, a period of 2.5 to 3 years can play the role of the dominant short-term period in this area. In general, for the long-term behavior, the two periods of 10 and 30 years can be considered as being the dominant periods across Iran. Among the long-term periods, the frequency of occurrence of these two periods is maximum in the first component and reduces in the other two components. In the third component, the minimum frequency of occurrence of these two periods is seen.
    Keywords: Drought, effective drought index, spectral analysis, cycle, Iran
  • Forough Kalvandi, Khalil Motaghi *, Esmaeil Shabanian Pages 48-61
    The Zagros mountain belt, situated on the northern margin of the Arabian plate, is one of the youngest continental collision belts. This belt was formed by a collision between the Arabian plate and the Central Iranian micro-continent. In this study, we used data from 38 temporary seismological stations installed on a 400 km long profile from May to November 2003.The trend of the profile is N58°E across northern Zagros and part of the Central Iran. The stations are part of Zagros03 profile (Paul et al., 2010) operated by the International Institute of Earthquake Engineering and Seismology (IIEES) of Iran in collaboration with CNRS - Université Joseph Fourier, France. We examine the structure of the lithosphere, across the profile by analysis of P-wave receiver functions and Rayleigh wave fundamental mode phase velocity dispersion curves. Joint inversion of Rayleigh wave phase velocity dispersion and receiver functions have been used to estimate the velocity structure beneath 28 seismic stations. Receiver functions are time-series computed from the three-component body-wave seismograms and are sensitive to the earth structure near the receiver station. They are composed of P- to S-wave conversions in discontinuities under the stations. These converted waves are isolated by deconvolving the vertical component of a teleseismic P-wave record from its radial component. For each event, a 120s time-window centered at the direct P arrival is selected and used for the calculation of the receiver function. The deconvolution used is the iterative deconvolution method of Ligorria and Ammon (1999). Surface waves arise from the presence (boundary conditions) of the stress-free surface of the Earth, and in the presence of layering, they are dispersed. They provide valuable information on the absolute S-wave velocity, but they are relatively insensitive to sharp velocity contrasts. On the other hand, receiver functions are sensitive to S-wave velocity contrasts, which give rise to converted phases, but allow for a substantial trade-off between the depth and velocity above an impedance change. Combining them in a joint inversion process bridges the resolution gaps associated with each data set. We jointly inverted the stacked receiver function and surface wave dispersion data. We employ the program joint96 which is available in the software package “Computer Program in Seismology” (Herrmann and Ammon, 2003). In this study, we try to calculate the Moho depth and velocity structure in the north Zagros collision zone using the joint inversion of receiver function and surface wave dispersions. Receiver functions are calculated using teleseismic events of magnitude greater than 5.1, located between 30◦ and 95◦ epicentral distances. The fundamental mode Rayleigh-wave group velocities are extracted from thetomographic study conducted by Rahimi et al. (2014). The 1D velocity models resolved by joint inversion are juxtaposed, and a 2D velocity model is obtained. Results obtained from the 2D model reveal that the thickness of the sediments beneath the Zagros is 12 km, the Moho depth beneath this region of Zagros is 43-57 km, which increases towards Sanandaj–Sirjan zone and Urumieh–Dokhtar magmatic arc and reaches an expanse of 62 km and then decreases in the central Iran with a depth of 42 km. The velocity model confirms the presence of a crustal root and a thick high-velocity lithosphere beneath and north of the suture. These evidences imply that the Arabian plate continues to underthrust beneath the Central Iran.
    The Zagros mountain belt, situated on the northern margin of the Arabian plate, is one of the youngest continental collision belts. This belt was formed by a collision between the Arabian plate and the Central Iranian micro-continent. In this study, we used data from 38 temporary seismological stations installed on a 400 km long profile from May to November 2003.The trend of the profile is N58°E across northern Zagros and part of the Central Iran. The stations are part of Zagros03 profile (Paul et al., 2010) operated by the International Institute of Earthquake Engineering and Seismology (IIEES) of Iran in collaboration with CNRS - Université Joseph Fourier, France. We examine the structure of the lithosphere, across the profile by analysis of P-wave receiver functions and Rayleigh wave fundamental mode phase velocity dispersion curves. Joint inversion of Rayleigh wave phase velocity dispersion and receiver functions have been used to estimate the velocity structure beneath 28 seismic stations. Receiver functions are time-series computed from the three-component body-wave seismograms and are sensitive to the earth structure near the receiver station. They are composed of P- to S-wave conversions in discontinuities under the stations. These converted waves are isolated by deconvolving the vertical component of a teleseismic P-wave record from its radial component. For each event, a 120s time-window centered at the direct P arrival is selected and used for the calculation of the receiver function. The deconvolution used is the iterative deconvolution method of Ligorria and Ammon (1999). Surface waves arise from the presence (boundary conditions) of the stress-free surface of the Earth, and in the presence of layering, they are dispersed. They provide valuable information on the absolute S-wave velocity, but they are relatively insensitive to sharp velocity contrasts. On the other hand, receiver functions are sensitive to S-wave velocity contrasts, which give rise to converted phases, but allow for a substantial trade-off between the depth and velocity above an impedance change. Combining them in a joint inversion process bridges the resolution gaps associated with each data set. We jointly inverted the stacked receiver function and surface wave dispersion data. We employ the program joint96 which is available in the software package “Computer Program in Seismology” (Herrmann and Ammon, 2003). In this study, we try to calculate the Moho depth and velocity structure in the north Zagros collision zone using the joint inversion of receiver function and surface wave dispersions. Receiver functions are calculated using teleseismic events of magnitude greater than 5.1, located between 30◦ and 95◦ epicentral distances. The fundamental mode Rayleigh-wave group velocities are extracted from thetomographic study conducted by Rahimi et al. (2014). The 1D velocity models resolved by joint inversion are juxtaposed, and a 2D velocity model is obtained. Results obtained from the 2D model reveal that the thickness of the sediments beneath the Zagros is 12 km, the Moho depth beneath this region of Zagros is 43-57 km, which increases towards Sanandaj–Sirjan zone and Urumieh–Dokhtar magmatic arc and reaches an expanse of 62 km and then decreases in the central Iran with a depth of 42 km. The velocity model confirms the presence of a crustal root and a thick high-velocity lithosphere beneath and north of the suture. These evidences imply that the Arabian plate continues to underthrust beneath the Central Iran.
    Keywords: receiver functions, surface wave dispersion, joint inversion, velocity structure, north Zagroscollision
  • Seyedeh Safoora Seddigh Marvasti, Anand Gnanadesikan, Abbasali Aliakbari Bidokhti *, Sarmad Ghader, Mohammad Seddigh Mortezavi Pages 62-74
    This study evaluates the performance of Earth system models for accurately simulating the phytoplankton productivity and bloom dynamics in the Oman Sea and the northwest of Arabian Sea. Satellite data (SeaWIFS ocean color) show two climatological blooms in this region, a wintertime bloom peaking in February and a summertime bloom peaking in September. On a regional scale, interannual variability of the wintertime bloom is dominated by cyclonic eddies which vary in location from year to year. During the wintertime, while both cooling in the winter and eddies control the blooms, variability in bloom location will arise from variability in the location of eddies, and so may not be predictable. In contrast, during the Southwest Monsoon, the dominant upwelling associated with the intense environmental forcing supersedes the effects of eddies, and the activity of the cold eddies is not pronounced. We consider numerical results from five different 3-D global Earth system models, which are denoted by CORE-TOPAZ, Coupled-TOPAZ, Coupled-BLING, Coupled-miniBLING, and the Geophysical Fluid Dynamics Laboratory (GFDL) Climate Model version 2.6 (CM2.6 miniBLING). Two coarse (1° grid resolution) models with a relatively complex biogeochemistry (TOPAZ: Tracers of Ocean Productivity with Allometric Zooplankton) capture the annual cycle but fail to capture both the eddies and the interannual variability. The results showed that the models differ from the observational data in terms of interannual variability. The low-resolution models (CORE- and coupled-TOPAZ) provide an almost uniform seasonal coefficient of variation, while both the data and eddy resolving CM2.6 models show higher interannual variability and seasonal changes. The coefficients of variabilities are particularly higher during the winter and summer blooms in the observations, while the low-resolution models do not see these signals. In other words, the low-resolution models fail to attain enough variability, while the high-resolution models (i.e. CM2.6) produce too much interannual variability. Accordingly, eddies are necessary to explain the variability in the data as opposed to the low-resolution models, but that the high-resolution model does not properly capture this variability. An eddy-resolving model (GFDL CM2.6) with a simpler biogeochemistry (miniBLING) displays larger interannual variability, but overestimates the wintertime bloom and captures eddy-bloom coupling in the south but not in the north. The models fail to capture both the magnitude of the wintertime bloom and its modulation by the eddies in part because of their failure to capture the observed sharp thermocline/nutricline in this region. In the wintertime, this leads to the excessive convective supply of nutrients and too strong of a bloom. However, for a few cases, eddies with blooms at the center are tracked in the southern part of the domain. For the model to simulate the observed wintertime blooms within cyclones, it will be necessary to represent this relatively unusual nutrient structure as well as the cyclonic or cold eddies. Both the temperature and mixed layer biases in the northern part of the Arabian Sea may result from having too much water from the Persian Gulf in this region. This is a challenge in the northern Arabian Sea as it requires capturing the details of the outflow from the Persian Gulf, something that is poorly done in global models.
    Keywords: red tide, bloom, numerical modeling, eddy, Earth system models, biogeochemistry
  • Farideh Canaani, Mostafa Fazli *, Victoria Ezzatian Pages 75-84
    The concentration of ground-level ozone is the result of thousands of complex chemical reactions. Basically, an increase of ozone concentration occurs in the presence of NOx, VOCs, and the sun’s radiation. This study deals with analyzing the ground-level ozone in Isfahan city from November 22 to December 21 (a full solar month) in 2009. According to the observations made by the Isfahan Meteorological Organization, photochemical smog was visible over the city during this month. The data used in this study include NO2, NO and O3 concentrations and the meteorological variables of temperature, relative humidity and wind speed which have been measured in Isfahan in 2009. The analyses were carried out for the sunshine hours in two time periods of 9 to 12 am and 6 am to 15 pm whose main characteristics are:A) During 6–15 period: the sun rises from 6 am and by becoming closer to the dusk, i.e. about 15 pm, both the radiation intensity, and temperature decrease;
    B) During 9–12 period: the higher temperature, radiation intensity, and traffic are the effective factors in the emission of pollutants when compared to the other hours of the day.
    The days under study are classified based on maximum, minimum and average ozone concentration. In order to analyze the tropospheric ozone and smog creation, in this study, the photochemical and semi-empirical models were used. The kinetic and mechanism of a number of photochemical reactions effective in ozone formation were taken into account in order to analyze the changes in ozone concentration. Calculations were carried out by using the Excel and Matlab software programs. Making use of the steady-state approximation method and considering oxygen atom in the steady state, the reaction rates have been computed. The differential relation obtained (d[O3]/dt = k2k1[NO2]- k3[O3][NO]) is a function of three variables of NO, NO2 and O3 concentrations. The amounts of reaction rate (d[O3]/dt) and also the rate constants k2k1 and k3 were also calculated. Analysis of the experimental relation between the activation energy and the results obtained from calculations indicate that the reactions that take place in the troposphere can be considered rank 3 reactions. In the troposphere, the quantum of energy (hν), which is released in some reactions, is very strong. The activation energies obtained for all days of this study include negative values, and this confirms the fact that the energy of the photons of the sun is needed to change NO2 and O2 to O3 in the troposphere. Based on the negative activation energies obtained, we can consider the reaction NO2㭯 hν→NO as the mechanism for the tropospheric ozone production in Isfahan.
    The creation of photochemical smog, SP(t) with t denoting the time, during the mentioned days was studied based on the Jonson’s semi-empirical model. The relations obtained based on the changes in NO and O3 concentrations with respect to time show that smog creation follows a quadratic nonlinear relation.
    In general, the increase of the concentration of pollutants on the ground as a result of photolysis reactions has led to the production of ozone concentrations. The results achieved from the analysis of reaction rate, smog creation, and the resulting curves indicate that ozone concentration has not been uniformly increasing or decreasing during the studied days, but there were both the increasing and decreasing trends. In general, the photochemical reactions taken place in the atmosphere of Isfahan city have caused both production and loss of ozone. Consequently, the investigation showed that the changes of the ozone concentrations under the effects of the solar radiation followed the same pattern in autumn 2009. On the other hand, the changes in the ozone concentration on the ground level caused changes in smog creation during the studied time. With regard to the above-mentioned arguments and on the basis of the effect of nearly the same conditions, a constant process prevailed. It could be predicted that this pattern will be the same in the future years.
    Keywords: tropospheric ozone, ozone reaction rate, activation energy, photochemical model, smog production algorithm
  • Fatemeh Azampour, Amin Roshandel Kahoo*, Alireza Ahmadyfard, Alireza Heydarian Pages 85-96
    Faults are among the most important geological events for hydrocarbon and mineral exploration and geological studies. They are considered as the major hydrocarbon traps whose detection in seismic reflection data have important roles in hydrocarbon reservoir characterization and geomechanical studies of reservoirs. According to previous studies, approximately 75% of the oil fields are associated with faults. Faults can influence the efficiency of hydrocarbon reservoirs by improving the permeability of a porous medium. Faults are defined as a discontinuity along a geological layer or geological event which is the result of the failure of that layer or event against the tension exerted on it.
    Reflection seismology is one of the best geophysical methods for hydrocarbon exploration. Various methods have been introduced to identify faults in seismic reflection data. In seismic reflection data, faults can be seen as discontinuities along the reflectors. One of the conventional methods of interpretation of the faults is a manual interpretation of seismic horizons which is a highly time-consuming process. Due to the nature of the seismic data, identification of faults plane in this data is a difficult process.
    Seismic attributes are considered as useful tools that will reveal hidden information in seismic data which can help the interpreters to detect faults. In fact, a seismic attribute can be used as a filter that makes the structural and stratigraphic information more apparent from the seismic data. There are several seismic attributes for faults detection in reflection seismic data such as coherency, curvature, chaos, and variance. The seismic coherence attribute is one of the most common attributes for fault detection. This attribute is calculated by different criteria such as cross-correlation, semblance, eigenstructure, gradient structure tensor.
    In this paper, we used a new form of a seismic attribute for edge detection. It is the Sobel filter which is a widely used tool in image processing, computer vision, and edge detection problems. It is the first derivative of the image and is sensitive to amplitude changes. It is a discrete differentiation operator, computing an approximation of the gradient of the image intensity function by convolving the isotropic operators. These operators are extensible to higher sizes and dimensions. Therefore, the Sobel attributes can be used in two- and three-dimensional seismic data. It can be used in seismic data for identification of the geological events such as faults, salt dome, and buried channels.
    The efficiency of dip guided Sobel attribute is evaluated by applying it to both synthetic and real seismic data. We compared the obtained results with the results of conventional seismic coherency attributes. In this paper, we used the eigenstructure-based coherency attribute. Comparison of the results of synthetic models in both noise-free and noisy cases in two and three dimensions show that the dip guided Sobel filter can be a good alternative for coherence attributes. In comparison with the coherence attributes, dip guided Sobel attribute is a short run time process and has large stability against noise. In real seismic tests in two and three dimensions, the obtained results from two attributes show that the dip guided Sobel attribute performs better than the eigenstructure-based coherency attribute.
    Keywords: reflection seismology, Fault, edge detection, Sobel attribute
  • Mojtaba Naghavi, Habib Rahimi *, Ali Moradi Pages 97-111
    The purpose of this study is to estimate compressional and shear wave quality factors of seismic waves by using local earthquakes occurred in the NW of Iranian Plateau. In seismological engineering studies, quality factor estimation of body and shear waves plays an important role in seismic risk assessment of different areas, determining the exact magnitude of the earthquake, strong ground motion simulation and study of destructive energy of earthquake from near to intermediate region. In this study, earthquakes recorded in the Iranian Seismological Center (IRSC) and Iranian National Broadband Seismic Network (INSN) for the longitudinal band from 43°E to 53°E and the latitude band from 36°N to 40°N were used. Among the 17 stations, 14 stations belong to the IRSC and the rest belong to the INSN. Due to the presence of some big cities in the northwestern part of Iranian Plateau, quality factors of body and shear waves were estimated by using the data of 17 seismological stations including 13000 recorded earthquakes of the IRSC and INSN.
    This region of intense deformation is situated between two thrust belts of the Caucasus to the north and the Zagros Mountains to the south. The NW of Iranian Plateau is a part of Turkish–Iranian plateau and includes historical and destructive earthquakes and two volcanoes with lots of thermal springs around. The North Tabriz Fault (NTF) is one of the active faults in NW Iran that has a clear surface expression.
    Seismic quiescence and large historical earthquakes in the region in more than the two last centuries have increased the seismic risk of this region. To estimate seismic hazard in an area, a two-step process is needed. First, we must understand the nature of the earthquake sources that generate potentially hazardous ground motion. This includes knowledge of the distribution of seismic source zones, predominant fault mechanisms and return times of large events. Second, we must understand the effects of the transmitting medium (the Earth) on the seismic waves. A synthesis of the source and path effects will allow us to calculate the ground motion at a given site. Seismic attenuation is also caused by intrinsic mechanisms that convert the wave energy to heat through friction, viscosity, and thermal relaxation processes. Scattering redistributes wave energy within the medium but does not remove energy from the overall wavefield. In contrast, intrinsic attenuation mechanisms convert the wave energy to heat through friction, viscosity, and thermal relaxation processes. Energy loss caused by inelastic behavior is called inherent or internal attenuation and is determined by the inverse of the Q parameter. Large values of quality factor mean that attenuation is low and when Q is equal to zero, attenuation is very high.
    Aki (1980) used the normalized Coda for the first time in order to estimate absorption amplitude of the S waves. Since then, this method has frequently been used in seismological studies for estimation of the absorption parameters of seismic waves (see, for example, Yoshimoto, 1993; Hatzidimitriou, 1995).
    For three categories of data with epicentral distances less than 100 km, from 100 to 200 km and 0 to 200 km, attenuation variation investigation of body waves was carried out in 9 frequency bands with central frequencies of 3, 5, 7, 10, 14, 20, 28, 38 and 47 Hz and the quality factor was estimated in different frequencies for each station, separately. For the northwesern part of Iran, the frequency dependence of the body and shear wave quality factors in all stations were estimated so that their average values are quantified as Qp=55f 0.84 and
    Qs=38f 0.93 , respectively. Due to the low values of the Q parameter and thus high attenuation values of body and shear waves in North West of Iranian Plateau, the amplitude of the propagated waves are decreased severely in the interested area when these waves pass through it. The attenuation effect of seismic waves would reduce the damages caused by the earthquakes at appropriate distances from the faults at the time of probable earthquake occurence.
    Keywords: compressive, shear waves, NW of Iranian Plateau, Attenuation, quality factor
  • Mehdi Rastgoo, Habib Rahimi *, Hosein Hamzehloo, Vahid Gholami Pages 112-127
    The southeast region of Iran includes the western part of Makran as an active subduction zone on the North side of the Oman Sea. The velocity structure of this region is not well understood because of its low-level seismicity and a small number of permanent seismic stations in this region. The main purpose of this study is the estimation of shear wave velocity structure and Moho depth variations in the southeast of Iran. To this end, we apply the joint inversion process of Rayleigh and Love waves dispersion curves and the “P-wave receiver function” (PRF) around four permanent broadband seismic stations of the region named ZHSF, KRBR, CHBR, and BNDS. To find the group velocities of surface waves in the southeast of Iran, the “frequency–time analysis” (FTAN) is applied to the waveforms of 40 local earthquakes, which occurred in Rigan region and recorded by the four permanent broadband seismic stations. These local earthquakes include two main shocks of 20 December 2010 (ML=6.4) and 27 January 2011 (ML=6.2) accompanied by their foreshocks and aftershocks in Rigan region. Therefore, the group velocities of fundamental modes of Rayleigh and Love waves are calculated in the period range from 5 to 60 sec in the paths among the four seismic stations and the sources of Rigan earthquakes. Also, to calculate P-wave receiver functions around the four permanent broadband seismic stations, 485 teleseismic earthquakes with suitable signal-to-noise ratios and epicentral distances between 30° and 95° related to the region, are selected. The radial PRFs are computed by deconvolving the vertical component from the radial component based on the iterative deconvolution method (Ligorria and Ammon, 1999). After preparing these two groups of data, we can determine the shear wave velocity structure and Moho depth vicinity of each seismic station by applying the joint inversion process to the dispersion data and the PRF data related to each seismic station (using the joint96 program; Herrmann and Ammon, 2007). Based on the results obtained in this study, the average group velocity of surface waves was estimated at less than 3.5 kms-1 in the period range from 5 to 60 sec. The lowest average group velocity of surface waves was obtained in the paths between the CHBR station and the sources of the Rigan earthquakes. Also, the Moho depths beneath the ZHSF, KRBR, CHBR, and BNDS stations were estimated 38±4, 46±6, 26±2 and 56±5 km, respectively. The minimum thickness of crust beneath CHBR station as well as the higher velocity of shear wave estimated beneath this station, are consistent with the shallow subduction of a high-velocity oceanic crust of Arabian plate beneath the south side of Makran. Furthermore, the thicker crust beneath the KRBR station and the lower velocity of shear wave estimated in this area, when compared with the area encompassing the CHBR station, is due to the existence of magmatic assemblage in the vicinity of the KRBR station. These results are consistent with the crust thickening from the south to the north of Makran. The maximum thickness beneath the BNDS is due to the location of this station being in the southeast of Zagros mountain belt, where the thick continental Arabian plate collides with Zagros. This collision leads to thickening of crust in Zagros.
    Keywords: southeast of Iran, shear wave velocity, Moho depth, surface waves dispersion, P, Wave Receiver Function
  • Khalil Motaghi *, Zahra Zarunizadeh, Abdolreza Ghods Pages 128-141
    Estimation of seismic wave attenuation due to anelasticity and geometrical spreading has attracted major interests among earthquake engineering community in recent decades. The choice of ground-motion model has a significant impact on hazard estimates in an active seismic zone such as the NW Iran. Estimation of ground motion for a typical frequency range of 0.5–10 Hz is required for the proper design of earthquake resistant structures and facilities and is considered as input for engineering stochastic ground motion relationships. For seismological purposes, appropriate attenuation models make it possible to calculate more accurately the source parameters such as magnitude and seismic moment. The NW Iran has experienced very few large events during the operation of the accelerometer network of the Building and Housing Research Center (BHRC). The BHRC network has been operating since 1973 but has recorded ground acceleration for few events in the study area, because of the low seismicity rate. The availability of the abundant weak-motion waveform data from the short-period local seismograph network of the Institute of Geophysics of the University of Tehran (IGUT) provides an opportunity to derive a new and more reliable ground-motion relationship for small events to complement those of strong-motion results. In this study, we analysed 3514 records of 943 small and moderate events that were recorded by 8 permanent stations of Tabriz network (belonging to the IGUT) and 16 temporary stations of the Institute for Advanced Studies in Basic Sciences (IASBS) to prepare a dataset including week ground-motion spectral amplitudes for different magnitudes and hypocentral distances. We graphically found the distance at which the nature of geometrical spreading attenuation changes significantly using a locally weighted scatter-plot smoothing called robust LOWESS. A bilinear function with a hinge at distance of about 70 km describes the geometric spreading attenuation with distance. Geometrical spreading and intrinsic attenuation coefficients were calculated using nonlinear regression in different frequencies and an average value of was found as geometrical spreading coefficient for distance range of 10–70 km. This value is consistent with geometrical spreading in a layered Earth. The average geometrical spreading coefficient of was found for the frequency range 0.79–5 Hz and the distance range of 70–200 km. This value is smaller than the values reported for other regions in the world (e.g. .09 for Central Alborz: Motaghi and Ghods, 2012; .2 for North Iran: Motazedian, 2006; .2 for SE Canada and the NE United States: Atkinson, 2004; .1 for SE Australia: Allen et al., 2007) and indicates that the velocity contrast in the Moho discontinuity is smaller than that in the other regions. The low-velocity uppermost mantle in NW Iran was manifested by different types of tomographic results obtained for the region. The geometrical spreading coefficient does not change before and after 70 km distance for frequencies ≥ 5 Hz. Thus, the attenuation relationship in this frequency range changed from bilinear to linear function. Using anelastic attenuation coefficients calculated at different frequencies, the shear-wave quality factor, , obtained equal to for frequencies greater than 1.5 Hz. In fact, the values show a U-shaped behavior in all of the frequency ranges and the function that describes it is defined as .
    Estimation of seismic wave attenuation due to anelasticity and geometrical spreading has attracted major interests among earthquake engineering community in recent decades. The choice of ground-motion model has a significant impact on hazard estimates in an active seismic zone such as the NW Iran. Estimation of ground motion for a typical frequency range of 0.5–10 Hz is required for the proper design of earthquake resistant structures and facilities and is considered as input for engineering stochastic ground motion relationships. For seismological purposes, appropriate attenuation models make it possible to calculate more accurately the source parameters such as magnitude and seismic moment. The NW Iran has experienced very few large events during the operation of the accelerometer network of the Building and Housing Research Center (BHRC). The BHRC network has been operating since 1973 but has recorded ground acceleration for few events in the study area, because of the low seismicity rate. The availability of the abundant weak-motion waveform data from the short-period local seismograph network of the Institute of Geophysics of the University of Tehran (IGUT) provides an opportunity to derive a new and more reliable ground-motion relationship for small events to complement those of strong-motion results. In this study, we analysed 3514 records of 943 small and moderate events that were recorded by 8 permanent stations of Tabriz network (belonging to the IGUT) and 16 temporary stations of the Institute for Advanced Studies in Basic Sciences (IASBS) to prepare a dataset including week ground-motion spectral amplitudes for different magnitudes and hypocentral distances. We graphically found the distance at which the nature of geometrical spreading attenuation changes significantly using a locally weighted scatter-plot smoothing called robust LOWESS. A bilinear function with a hinge at distance of about 70 km describes the geometric spreading attenuation with distance. Geometrical spreading and intrinsic attenuation coefficients were calculated using nonlinear regression in different frequencies and an average value of was found as geometrical spreading coefficient for distance range of 10–70 km. This value is consistent with geometrical spreading in a layered Earth. The average geometrical spreading coefficient of was found for the frequency range 0.79–5 Hz and the distance range of 70–200 km. This value is smaller than the values reported for other regions in the world (e.g. .09 for Central Alborz: Motaghi and Ghods, 2012; .2 for North Iran: Motazedian, 2006; .2 for SE Canada and the NE United States: Atkinson, 2004; .1 for SE Australia: Allen et al., 2007) and indicates that the velocity contrast in the Moho discontinuity is smaller than that in the other regions. The low-velocity uppermost mantle in NW Iran was manifested by different types of tomographic results obtained for the region. The geometrical spreading coefficient does not change before and after 70 km distance for frequencies ≥ 5 Hz. Thus, the attenuation relationship in this frequency range changed from bilinear to linear function. Using anelastic attenuation coefficients calculated at different frequencies, the shear-wave quality factor, , obtained equal to for frequencies greater than 1.5 Hz. In fact, the values show a U-shaped behavior in all of the frequency ranges and the function that describes it is defined as .
    Keywords: seismic attenuation relation, quality factor, NW Iran