فهرست مطالب

مجله ژئوفیزیک ایران
سال دوازدهم شماره 2 (پیاپی 39، تابستان 1397)

  • تاریخ انتشار: 1397/07/04
  • تعداد عناوین: 9
|
  • زینب ذاکری، مجید آزادی*، سرمد قادر صفحات 1-22

    در این مطالعه اثر گوارد داده های تابندگی ماهواره، داده های prepbufr که شامل مجموعه ای از داده های سطح زمین و جو بالاست و داده های GPSro با استفاده از سامانه داده گواری WRFDA در بهبود پیش بینی باد و بار گرد و خاک در مدل WRF-Chem بررسی شده است. مطالعات انجام شده روی دو مورد گرد و خاک در غرب کشور در تاریخ های 15 ژوئن 2016 و 31 اوت 2015 با منشا کشور عراق بوده است. برای هر مورد، دو آزمایش مختلف، یک آزمایش داده گواری به روش وردشی سه بعدی و با استفاده از خطای زمینه محاسبه شده برای حوزه اجرای مدل و یک آزمایش کنترلی بدون داده گواری انجام شده است. مقایسه نقشه های ماهواره با پیش بینی بار گرد و خاک مدل نشان می دهد که با انجام داده گواری مدل محل گسیل و مسیر ترابرد گرد و خاک در ساعت های اولیه پیش بینی (24 ساعت اول) را با دقت بیشتری پیش بینی می کند؛ اما در ادامه و با زیاد شدن سن پیش بینی برونداد دو آزمایش بدون داده گواری و با داده گواری بسیار شبیه هم می شوند. در ارزیابی کمی خطای سرعت باد ملاحظه می شود که میانگین قدر مطلق خطا به طور سازگار در ترازهای 850 و 700 هکتوپاسکال و 10 متری سطح زمین با انجام داده گواری تا حدود 11 درصد کاهش می یابد. به تدریج و با زیاد شدن سن پیش بینی تاثیر مثبت داده گواری کاهش می یابد.

    کلیدواژگان: داده گواری، سامانه WRFDA، داده تابندگی ماهواره، پیش بینی گرد و خاک، مدل WRF-Chem
  • خلیل کرمی ، سرمد قادر، سید وحید موسوی صفحات 23-37
    در این مطالعه،بازتاب مستقیم فعالیت های موجی از پوشن سپهر به وردسپهر یا جذب و انتشار این امواج در پوشن سپهر شناسایی می شود. به این منظور، کمیت های ضریب بازتاب پرلویتز- هارنیک، شار گرمایی، گرمایش ناگهانی پوشن سپهری، ضریب شکست امواج راسبی و احتمال انتشار بهینه امواج راسبی محاسبه، تحلیل و بررسی شده اند. برای اجرای تحقیق، از داده های بازتحلیل ECMWF (ERA-Interim) در بازه زمانی 2014-1979 و در فصل زمستان (نوامبر- مارس) در نیمکره شمالی استفاده شده است. بررسی شاخص ضریب بازتاب پرلویتز- هارنیک نشان می دهد که در بازه زمانی 34 ساله بررسی شده در 10 زمستان، پوشن سپهر حالت بازتابی و در 24 زمستان، حالت غیربازتابی (جذب یا انتشار امواج راسبی) داشته است. بررسی این شاخص نشان می دهد احتمال بازتاب مستقیم امواج راسبی از پوشن سپهر به وردسپهر در پاییز و اوایل زمستان (نوامبر- دسامبر) کمینه است. در بازه زمانی ژانویه- مارس، پوشن سپهر حالت بازتابی داشته و بیشترین احتمال جفت شدگی بین پوشن سپهر و وردسپهر از طریق بازتاب امواج راسبی در اواسط زمستان و اوایل بهار است. بررسی شاخص شار گرمایی نشان می دهد استفاده از تعریف موج به عنوان پریشیدگی از حالت میانگین مداری نسبت به تعریف موج به عنوان پریشیدگی از میانگین زمانی، مزیت دارد. در سال هایی که پوشن سپهر حالت بازتابی دارد، تنها در 33 درصد از زمستان ها احتمال رخداد پدیده گرمایش ناگهانی پوشن سپهری وجود دارد. این در حالی است که در زمستان هایی که پوشن سپهر حالت غیربازتابی (جذب یا انتشار امواج راسبی) دارد، در حدود 62 درصد از زمستان ها احتمال رخداد این پدیده وجود دارد. نتایج بررسی دو شاخص ضریب شکست امواج راسبی و احتمال انتشار بهینه امواج راسبی نشان می دهند که در حالت پوشن سپهر غیربازتابی نسبت به پوشن سپهر بازتابی، احتمال انتشار و نفوذ امواج راسبی از وردسپهر به پوشن سپهر زیرین (100-50 هکتوپاسکال) در فاصله عرض جغرافیایی 90-70 درجه شمالی بیشتر است.
    کلیدواژگان: انتشار قائم امواج راسبی، شار فعالیت موجی، ضریب شکست امواج راسبی، گرمایش ناگهانی پوشن سپهر
  • ملیحه چالاک، علیرضا محب الحجه ، محمد میرزایی صفحات 38-49
    امواج گرانی - لختی، از جمله عواملی هستند که میدان های هواشناختی را تحت تاثیر خود قرار می دهند و با انتقال تکانه، انرژی و جرم سبب ایجاد پریشیدگی در میدان های باد، دما و فشار می شوند. بنابراین شناخت چشمه های تولید این امواج و بررسی ویژگی های آنها بسیار حائز اهمیت است. با توجه به اینکه تاکنون مطالعه ای برای شناسایی امواج گرانی - لختی در یک دوره اقلیمی روی تهران صورت نگرفته است، در این پژوهش به کمک داده های دیدبانی ایستگاه مهرآباد تهران در دوره 1961 تا 2015، به شناسایی رخدادهای این پدیده در وردسپهر زیرین پرداخته شده است.
    یکی از روش های شناسایی امواج گرانی - لختی، بررسی تغییرات ناگهانی و بزرگ در فشار سطحی است. در این مطالعه، با در نظر گرفتن آستانه تغییرات ساعتی فشار به مقدار بیش از 3 هکتوپاسکال، 101 مورد به عنوان نامزد های احتمالی رخداد امواج گرانی - لختی برای یک دوره 55 ساله در ایستگاه مهر آباد تهران شناسایی شد. سپس، با بررسی داده های دیدبانی ایستگاهی و تحلیل الگوهای همدیدی مناسب برای فعالیت امواج گرانی - لختی ، تعداد 17 مورد رخداد موج گرانی - لختی تعیین شد که 13 مورد آن با همرفت و 4 مورد آن با چرخند همراه بوده اند. در ادامه، برای ارزیابی روش قبلی و تعیین ویژگی های امواج، از مدل WRF برای شبیه سازی موارد رخداد موج گرانی - لختی استفاده شد. . نتایج شبیه سازی عددی نه تنها درستی تحلیل های انجام شده را در انتخاب موارد رخداد موج گرانی - لختی تایید کرد،، بلکه ارتباط بین چشمه امواج، نحوه انتشار و همچنین تاثیرشان را را بر میدان های هوا شناختی مانند فشار، دما و باد در منطقه تهران و به خصوص در ایستگاه مهر آباد، با دقت بیشتری تعیین کرد.
    کلیدواژگان: گرانی - لختی، شبیه سازی عددی، میدان های هواشناختی
  • ناهید شادمانی، محمدعلی نصراصفهانی ، احمدرضا قاسمی صفحات 50-63
    هدف مطالعه حاضر، بررسی منابع رطوبت و مسیر حرکت آن برای دو رویداد سیل شدید در غرب و جنوب ایران است. برای این هدف از الگوریتمی استفاده شده است که در سال های اخیر برای بررسی انتقال رطوبت از نواحی حاره ای به عرض های میانی توسعه داده شده است. بارش های سیل آسای منتخب در غرب و جنوب ایران به ترتیب روزهای 6 تا 9 آبان 1394 (28 تا 31 اکتبر 2015) و 19 تا 21 آبان 1394 (10 تا 12 نوامبر 2015) رخ داده اند که خسارت های سنگینی را برجای گذاشتند. داده های مورد استفاده در این پژوهش، داده های تحلیل شدهGFS با قدرت تفکیک نیم درجه در راستاهای طول و عرض جغرافیایی بودند که برای محاسبه انتگرال قائم بزرگی شار افقی بخار آب به کار رفتند. علاوه بر این، مولفه های مداری و نصف النهاری شار نیز جهت بررسی بیشتر ترسیم و تحلیل شده اند. نتایج نشان داد که برای سیل غرب ایران، دریای عرب و دریای سرخ تقریبا سهم یکسانی در تامین منبع رطوبت بارش رخداده دارند ولی مسیر حرکت بیانگر این است که منبع اصلی و مسیر حرکت رطوبت از سوی دریای سرخ است. اما برای سیل جنوب ایران، منبع اصلی تامین رطوبت در روزهای قبل از بارش، نواحی جنوبی دریای سرخ است ولی در روزهایی که با بارش همراه بوده است، خلیج فارس نقش اساسی در تامین رطوبت آن دارد و دریای سرخ و دریای عرب در زمان بارش نقش مهمی در تامین رطوبت ندارند.
    کلیدواژگان: بارش های سیل آسا، انتگرال قائم شار بخار آب، منبع رطوبتی، دریای سرخ، سیل
  • سعید هادیلو، سعید میرزایی ، حسین هاشمی، ابراهیم سفیداری صفحات 64-77
    تحلیل رخساره های لرزه ای، فنی است برای به نقشه درآوردن خصوصیات و ویژگی های زمین شناسی با استفاده از اطلاعات لرزه ای. برای تحلیل رخساره های لرزه ای با استفاده از الگوریتم های یادگیری ماشین، نشانگرهای لرزه ای دسته بندی می شوند. این دسته بندی به شناسایی رخساره های مختلف درون مقاطع یا افق های لرزه ای می انجامد. . به دلیل ماهیت داده های لرزه ای، که همواره درجه ای از عدم قطعیت دارند، تنوع نشانگر های لرزه ای و الگوریتم های یادگیری ماشین می تواند در تحلیل رخساره های لرزه ای نتایج مختلفی ایجاد کند؛ به همین دلیل نیاز است مراحل مختلف تحلیل مانند انتخاب پارامترهای ورودی و الگوریتم های یادگیری ماشین، با دقت زیاد و با توجه به هدف تحلیل رخساره انتخاب شوند تا درنهایت، بتوان نقشه رخساره لرزه ای مناسبی تولید کرد. در این مطالعه، برای پیدا کردن مسیر بهینه تولید نقشه رخساره لرزه ای، روشی تعاملی با نظارت مفسر و با استفاده از نرم افزار بومی SeisART پیشنهاد می شود. در این روش،در روند اجرای تحلیل رخساره، پارامتر هایی مانند نوع نشانگرهای ورودی، تعداد خوشه، الگوریتم خوشه بندی و دسته بندی به صورت بهینه برای هدف موردنظر انتخاب می شوند و مفسر در یک روند بازگشتی و چرخشی می تواند نتایج تحلیل را با توجه به تغییر پارامترهای ورودی مقایسه کند، با در نظر گرفتن نتایج دانش بیشتری از رخساره های موجود به دست آورد و درنهایت، تحلیل رخساره مناسبی ارائه دهد. . ازآنجاکه در این روش، الگوریتم های مختلف با پارامترهای ورودی مختلف بر اطلاعات لرزه ای موجود اعمال می شوند و مفسر می تواند با توجه به کارایی هرکدام از پارامترها و روش ها، برای اعمال کردن یا نکردن آن تصمیم گیری کند، استفاده کردن از این الگوریتم پیشنهادی، این امکان را ایجاد می کند که رخساره های مختلف موجود در اطلاعات شناسایی شوند و ارتباط آنها با یکدیگر معین شود. روش پیشنهادی برای شناسایی رخساره های موجود در افق MSF4 از داده های لرزه ای سه بعدی F3 دریای شمال، آزموده شده است.
    کلیدواژگان: تحلیل رخساره لرزه ای، نشانگرهای لرزه ای، ماشین های یادگیری، سیستم های فازی، خوشه بندی
  • سعید سلطانی مقدم، کامران سپانلو ، میثم خیری ملومه صفحات 78-95
    زمین لرزه اخیر ملارد تهران و وقوع چهار زمین لرزه با بزرگای بیش از چهار طی یک دهه اخیر، به خوبی اهمیت مطالعه لرزه خیزی کلان شهر تهران را مشخص می کند. در این مطالعه، سعی شده با نگاهی جامع به اطلاعات لرزه ای موجود در یک دهه گذشته، لرزه خیزی شهر تهران و مناطق اطراف بررسی شود. به این منظور، از اطلاعات سه شبکه لرزه نگاری دائمی فعال در منطقه استفاده شد. با ادغام اطلاعات این سه شبکه، حذف انفجارات و داده های پرت، تعداد 2000 زمین لرزه در بانک داده اصلی قرار گرفت. با انتخاب همگن 310 رویداد به شیوه تفکیک سلولی در سه بعد، ساختار سرعتی منطقه با استفاده از الگوریتم بهینه سازی ازدحام ذرات محاسبه شد. سپس بانک داده با استفاده از برنامه غیرخطی احتمالاتی، مجددا مکان یابی شد. نتایج مکان یابی نشان می دهد 40 درصد از رویدادهای مکان یابی شده، خطای رومرکزی کمتر از 5/2 کیلومتر، خطای عمق کمتر از 5 کیلومتر و میانگین باقی مانده زمانی 24/0 ثانیه دارند. در ادامه، زمین لرزه ها براساس عدم قطعیت های محاسبه شده در سه کلاس مختلف دسته بندی شدند. برای بررسی دقیق تر هندسه گسل های فعال منطقه، چندین مقطع عمقی رسم شد. در انتها، سازوکار کانونی 38 زمین لرزه به وقوع پیوسته در تهران و محدوده اطراف آن، با روش قطبش اولین موج تراکمی، محاسبه شد. نتایج این بررسی نشان می دهد طی یک دهه گذشته، بخش های شرقی فعالیت بیشتری از خود نشان داده اند و زمین لرزه ملارد تهران در غرب، در منطقه ای با لرزه خیزی نسبتا کم، به وقوع پیوسته است. با بررسی میدان تنش، جهت تنش بیشینه، N36E، نسبت تنش برابر 12/0 و رژیم زمینساختی راستالغز فشاری به دست آمد. بر اساس زاویه لغزش در دایره مور، اکثر صفحات گسلی اصلی با روند شرقی- غربی، سازوکار چپ گرد را تایید می کنند. این نتایج، همخوانی بسیار خوبی با روندهای موجود در ناحیه شمال تهران نشان می دهند.
    کلیدواژگان: مکان یابی غیرخطی احتمالاتی، مدل سرعتی یک بعدی، سازوکار کانونی، میدان تنش
  • میثم امیری، زهرا موسوی ، خلیل متقی صفحات 96-108
    در 20 دسامبر 2010، زمین لرزه ای با بزرگای 5/6 Mw در منطقه ریگان واقع در جنوب شهر بم رخ داد. 37 روز بعد، در 27 ژانویه 2011 پسلرزه بزرگی (با بزرگای 2/6Mw) در 20 کیلومتری زمین لرزه اول اتفاق افتاد. برای مطالعه دگرشکلی پسالرزه ای پس از وقوع زمین لرزه های فوق، تصاویر ماهواره COSMO-SkyMed (از آژانس فضایی ایتالیا) در بازه زمانی 27 ژانویه 2011 و 15 جولای 2011 تهیه و تحلیل سری زمانی روی تصاویر با استفاده از روش کوتاه ترین خط مبنا (SBAS) انجام شد. نقشه میانگین سرعت به دست آمده از تحلیل سری زمانی کوتاه ترین خط مبنا، بیانگر وجود فاز پسالرزه ای با سازوکار راست بر در منطقه ریگان است. سری زمانی جا به جا یی دو طرف گسل مسبب زلزله اصلی، نشان دهنده 8 میلی متر جا به جا یی در مدت 5 ماه پس از زلزله است. پس از این مدت، فاز پسالرزه ای به اتمام می رسد، گسل قفل می شود و فاز میان لرزه ای با انباشت آرام تنش روی صفحه گسلی شروع می شود. تعیین پارامترهای چشمه فاز پسالرزه ای با استفاده از مدل اکادا (1985) ، نشان می دهد که یک توزیع لغزش با بیشینه جابه جایی 1 متر، در زیر گسیختگی زلزله اصلی (در عمق 10 کیلومتری) و همچنین در لبه های گسل متمرکز شده است. نتیجه مدل سازی نشان می دهدکه فاز پسالرزه ای مشاهده شده صرفا به دلیل وقوع پسلرزه ها در بخش های مجاور گسیختگی اصلی، روی صفحه مسبب زلزله اول است. سازوکار این دگرشکلی، پس لغزش است.
    کلیدواژگان: زلزله ریگان، فاز پسالرزه ای، تداخل سنجی راداری، نقشه سرعت متوسط، آنالیز سری زمانی
  • صدیقه مهرآور ، حبیباللهقائدامینی، سید محمد جعفر ناظم سادات صفحات 109-126
    النینو- نوسان های جنوبی (ENSO) و نوسان های مادن - جولیان (MJO) دو پدیده بزرگ مقیاس جوی - اقیانوسی در پهنه های استوایی اقیانوس آرام و هند هستند که بر نوسان های اقلیمی گستره های استوایی تا جنب حاره ای تاثیر قابل توجهی می گذارند. در این پژوهش، پس از تعریف فازهای مثبت و منفی MJO، برهم کنش این فازها با فازهای ENSO ارزیابی و تاثیر این برهم کنش بر ابرناکی و روانش رطوبت در سطح خاورمیانه و اقیانوس های آرام و هند بررسی و بازخورد این روانش رطوبتی بر رخداد بارش پاییزه 9 ایستگاه همدیدی و باران سنجی استان فارس در بازه زمانی 2013-1974 واکاوی شد. در دوران النینو، احتمال رخداد هریک از فازهای مثبت یا منفی MJO تقریبا یکسان است در حالی که در دوران لانینا، احتمال رخداد فاز مثبت MJO نزدیک به دو برابر فاز منفی است. چیرگی فاز منفی MJO در دوران النینو با افزایش کارکردهای همرفتی در گستره های استوایی غرب اقیانوس هند و تا اندازه ای در شرق اقیانوس آرام همراه است. این همزمانی، وزش بادهای نمناک شرقی را از اقیانوس هند به سوی شبه جزیره عربستان، خلیج فارس و جنوب ایران افزایش می‏دهد. در دوران النینو، بارش و احتمال رخداد آن در فاز منفی MJO در همه ایستگاه های بررسی شده به گونه معنی داری بیشتر از مقادیر مشابه در دوران فاز مثبتMJO است. در دوران لانینا، چیرگی رخداد فاز مثبت یا منفی MJO، ویژگی های بارش فارس نوسان درخور توجهی نداشت.
    با چیرگی فاز منفی MJO در دوران النینو، بیشترین و کمترین افزایش شدت بارش نسبت به فاز مثبت در این دوران، در بخش های شرقی و جنوب شرقی استان فارس به ترتیب به اندازه 20-5 و 100-80 درصد است. در بیشتر گستره های مرکزی و جنوب غربی استان، افزایش 60-40 درصدی و در بخش های شمالی و جنوب شرقی، افزایش 110-60 درصدی احتمال رخداد بارش فاز منفی MJO در دوران گرم ENSO نسبت به فاز مثبت آن در دوران گرم دیده می شود.
    کلیدواژگان: النینو، لانینا، نوسان های جنوبی، بارش پاییزه، استان فارس
  • رضا برهانی، فرهنگ احمدی گیوی * صفحات 127-146

    این پژوهش با هدف بررسی آماری - دینامیکی تاشدگی وردایست در منطقه جنوب غرب آسیا بین سال های 2000 تا 2015 انجام شده است. تاشدگی های وردایست به کمک الگوریتم ارائه شده توسط شرلاک و همکاران (2014) و با استفاده از داده های میان یابی شده مدل IFS (Cy31r2) مربوط به ECMWF شناسایی شده اند. سری زمانی انواع تاشدگی وردایست نشان می دهد که فراوانی تاشدگی در دوره مورد مطالعه روند افزایشی دارد. توزیع فصلی تاشدگی وردایست نیز بیانگر آن است که بیشترین میزان فراوانی تاشدگی در فصل تابستان و کمترین مقدار آن در پاییز وجود دارد. در فصل تابستان، تاشدگی های کم عمق و متوسط و در زمستان، تاشدگی های عمیق فراوانی نسبی بیشتری دارند. نتایج ضریب همبستگی بین میانگین ماهانه فراوانی تاشدگی و کمیت های دینامیکی، نشان دهنده تاثیر ناپایداری کژفشار در تشکیل تاشدگی های زمستانه و تاثیر عوامل ترمودینامیکی در تاشدگی های تابستانه است. برای بررسی دینامیکی تاشدگی در فصل زمستان، ماه های ژانویه 2001 و 2004 و در فصل تابستان، ماه های ژوئن 2007 و 2015 که به ترتیب کمترین و بیشترین فراوانی را در بین ماه های متناظر داشتند، انتخاب و تحلیل شدند. بررسی شار افقی فعالیت موج تراز میانی وردسپهر در ژانویه 2004، همزمان با وجود یک چشمه موج قوی در غرب مدیترانه، وجود چاهه موجی را روی اروپا آشکار می سازد. در ژانویه 2001، شار فعالیت موج ضمن تضعیف به دو شاخه تقسیم شده که یکی از شاخه ها در عرض های میانی و دیگری در جنب حاره است. مولفه استواسوی شاخه جنوبی، سبب انتقال انرژی موج به نواحی جنوب غرب آسیا شده است. قوی و ضعیف بودن چشمه ها و چاهه های فعالیت موج در این دو ماه را می توان به ترتیب به فراوانی زیاد و کم تاشدگی وردایست مربوط دانست. در فصل تابستان، فعالیت موج عرض های میانی به شدت تضعیف شده است و همزمان با تقویت موسمی هند، در عرض های جنب حاره و روی اقیانوس هند، چشمه ها و چاهه های موج ضعیف تری شکل می گیرند.

    کلیدواژگان: تاشدگی وردایست، جنوب غرب آسیا، تاوایی پتانسیلی، تحلیل دینامیکی، فعالیت موج، ناپایداری کژفشار
|
  • Zeinab Zakeri, majid azadi, Sarmad Ghader Pages 1-22

    Dust storms are significant phenomenon in south west Asian countries like Iraq, Syria and Iran. Mineral dust is generated by wind erosion over arid and semiarid land surfaces and is transported locally and over vast distances, causing adverse environmental and weather problems. Recent draughts over dust sources in Iraq and Syria have remarkably increased dust events in the area particularly over west of Iran. Real time prediction of dust storms especially quantitative forecasting of dust concentration has become highly desirable to alleviate its damaging consequences. In this study the impact of the assimilation of satellite radiance, prepbufr and GPSro data in the wind speed and dust load forecasts of WRF-Chem model using WRFDA system are investigated. Prepbufr data are a collection of surface and upper air observations and GPSro data are GPS radio occultation data. These data are operationally collected by the National Centers for Environmental Prediction (NCEP). Data assimilation is applied to two dust events starting from Iraq and Syria borders on August 31st 2015 and June 15th 2016. For each case, two experiments are conducted. An experiment assimilating above mentioned data with three dimensional variational (3D-Var) intermittent assimilation method and a control simulation with no assimilation. The assimilation cycles in the intermittent method consist of three subsequent analyses at 00, 06 and 12 UTC. After the last assimilation cycle, the model is integrated for 48 hours in the future. In variational data assimilation a key element to get a qualified analysis is the accurate specification of error statistics for the background forecast. For the calculation of background error, the model with the same specification for the experiments is run for the whole January 2014 at 0000 and 1200 UTC and the 12- and 24-h forecasts are used to calculate the background error using the National Meteorological Center (NMC) method with CV5 option. The horizontal resolution of the domain is 21 km with 142×130 grid points covering Iran and western neighboring countries. The model has 41 vertical levels with the model top at 25 hPa. Initial and boundary conditions are taken from NCEP Global Forecast System (GFS) model with the horizontal resolution of 0.5º×0.5º.Results show that the agreement between spatial distribution of dust load prediction of the model and Meteosat-10 satellite RGB images is improvedusing data assimilation especially in first forecast hours. Quantitative comparison of 10 m, 850 hPa and 700 hPa model wind speed with surface observation data and ERA-Interim ECMWF reanalysis data show up to 11% improvement in RMSE especially in first forecast hour times. The positive impact of data assimilation is decreased as the forecast length increases

    Keywords: data assimilation, WRFDA, satellite radiance data, dust forecast, WRF-Chem
  • Khalil Karami, Sarmad Ghader, Seyed Vahid Mousavi Pages 23-37
    We detect the downward planetary wave reflection from the stratosphere back to the troposphere in the Northern hemisphere extended winter season (November-March), using ECMWF (ERA-Interim) reanalysis data (1979-2014). In the previous studies, the wave activity is defined as departure from a long-term time mean. However, we demonstrate some of the shortcomings of the above-mentioned definition. For instance, negative values of the heat flux at the lower stratosphere does not necessarily show the downward wave propagation, but a lower upward wave propagation. Perlwitz-Harnik index of reflection is used to categorize the stratospheric wind regimes into the two distinct states: reflective and non-reflective. Negative and positive values of this index indicate a reflective and non-reflective (either absorptive or propagative) stratospheric states, respectively. Our results show that the negative values of this index during early winter (November-December) suggests that the downward wave coupling is less likely to happen in early winter and most of the direct downward wave coupling occurs during mid-winter and early spring (January-March), which is in agreement with previous studies. Furthermore, 10 stratospheric winter states (out of 34) are reflective, while the remaining states (24 winters) are non-reflective winters. Winters whose their stratosphere experiences a major Sudden Stratospheric Warming (SSW) event is identified to detect the absorptive states of the stratosphere. Our analysis suggests that only in the %33 of the winters with reflective stratosphere a SSW can occur. In the %62 and %38 of the winters with non-reflective stratosphere, a SSW event can occur or Rossby waves can propagate upward freely to the upper stratosphere, respectively. Analysis of the heat flux at the lower stratosphere (50 hPa) using two different definitions (wave as a deviation from the time mean and as a deviation from the zonal mean) provides some useful information. If the waves are defined as a deviation from the zonal mean (), during the reflective years, this quantity is negative (indicating a downward reflection of wave activity from the stratosphere to the troposphere). Similarly, during non-reflective years, this quantity has positive values (suggesting either upward wave propagation or absorption by the mean flow). While the above-mentioned definition is in harmony with our expectation, the definition of the wave as a deviation from the long-term time mean () results in an oscillation curve around zero line, without any useful information either about the upward wave propagation or a downward wave reflection. In order to understand the role of mean flow in influencing the upward wave propagation, we calculate the Rossby wave refractive index (or vertical wavenumber alternatively). Our analysis show some of the problematic features (for instance, a very noisy structure) of this index in understanding the differences of reflective and non-reflective stratospheric states. This is most probably due to the overlapping of very large or very small values of the refractive index which cancel each other and results in a noisy structure. To overcome this problem, we use a modified diagnostic tool (compared to the refractive index), probability of favorable propagation condition for Rossby waves. This index has some clear advantages over the traditional refractive index. Analysis of this index shows that during the non-reflective stratospheric states, Rossby waves are more likely to propagate upward (with higher values of the probability), in comparison to the reflective stratospheric regimes which is a superior demonstration of the influences of the basic mean flow on the upward wave propagation over the traditional refractive index.
    Keywords: upward Rossby wave propagation, wave flux, Rossby wave refractive index, sudden stratospheric warming
  • Malihe Chalak, Alireza Mohebalhojeh, Mohammad Mirzaei Pages 38-49
    Inertia–gravity waves (IGWs) play an important role in transporting energy and momentum, in making and shaping turbulence and mixing, and in influencing the mean circulation and thermal structure of the atmosphere. They cause perturbations in the main dynamical fields, such as pressure, temperature and wind. Therefore, knowledge of their sources and analysis of their characteristics are very important. Given that there has been no study to identify and simulate IGWs that may have occurred over Tehran during the last few decades, the current research aims to investigate IGWs in the lower troposphere based on the long-term data of Mehr-Abad meteorological station during the period from 1960 to 2015.
    One method for detecting IGWs is checking the large amplitude hourly surface pressure changes. At first, hourly surface pressure observations for Tehran in a 55-yr period have been used to find the distribution of large hourly pressure changes defined as falls or rises in excess of 3 hPa per hour. After making careful checks, 101 cases were identified as potential candidates of IGWs occurrences. In the second part of the paper, we constructed composite maps of mean sea level pressure, 1000-500 hPa thickness, 500 hPa geopotential height, and 850 hPa, 500 hPa and 200 hPa wind. The cyclone-related IGW composite showed that the IGW events clustered poleward of the estimated position of the surface warm front, northeastward of the surface cyclone center, downstream and within the poleward exit region of a jet streak that was upstream of a ridge. By analyzing the synoptic-scale environments in which IGWs evolve, the probable candidates for IGWs were reduced to 17, among which 13 events were associated with convection and 4 events were related to cyclone development.
    In order to evaluate our analysis and to obtain the IGWs characteristics, the Weather Research and Forecasting (WRF) model was used. We determined the position of wave packets and their propagation during the time evolution by drawing the horizontal velocity divergence. The horizontal and vertical wavelength, intrinsic frequency, the period, horizontal and vertical velocity and the group velocity were determined by drawing the cross sections of horizontal velocity divergence and making use of the dispersion relation for hydrostatic waves, as well as the relations for the horizontal and vertical components of group velocity and the intrinsic phase speed. Then we drew time series of the horizontal velocity divergence and sea level pressure, and evaluated the correlation between these two quantities. By checking the horizontal velocity divergence at 800 hPa which is near the surface of the earth, it was observed that negative/positive divergence is associated with pressure rise/fall. On the other hand, the strength or weakness of wave’s amplitude was investigated by using the time series of horizontal divergence at different pressure levels. The numerical simulation results not only confirmed the validity of the method employed to analyze and identify IGWs, but also determined the connection between the IGWs sources, their propagation as well as their effects on meteorological fields such as pressure, temperature and wind over Tehran, especially at Mehr-Abad meteorological station
    Keywords: inertia–gravity waves, numerical simulation, meteorological fields
  • Nahid Shadmani, Mohammad Ali Nasr Esfahani, AhmadReza Ghasemi Pages 50-63
    The aim of the present study is detecting the moisture sources of the heavy floods in the western and southern parts of Iran using hydrometeorological analysis. For this purpose, two heavy floods happened during 28th to 31th October and 10th to 12th November 2015 in the western and southern parts of Iran were selected. The Integrated Vapor Transport (IVT) algorithm was used to study the transfer path of the tropical moisture to the mid-latitudes region. This algorithm is also useful for detecting the moisture source of the heavy rainfall and floods. In mid-latitude locations, most of the heavy precipitation and flood events are related to intense vertically integrated horizontal water vapor transport. The Integrated Vapor Transport algorithm is in the base of accumulated atmospheric moisture transport from 1000 hPa to 300 hPa levels which is computed using the Global Forecast System (GFS) analyzed data with 0.5°0.5° resolution for all days of October to November from 2007 to 2015. Then threshold values of the IVT were computed for the selected western and southern flood events separately as median value of the maximum IVT along a meridional line at 45E between 20N to 45N for the western flood case and a diagonal line between 42E, 23N and 52E, 33N for the southern flood event one. To detecting the role of humidity associated with the Persian Gulf in the southern flood event, an additional diagonal line was also considered in the northern parts of the Persian Gulf between 56.5E, 25N and 47E, 34.5N. To detect the source and path of the moisture for the selected floods, the IVT algoritm (as mentioned in the main text) is used. In addition to the IVT, the zonal and meridional components of integrated vertical moisture flux vector were computed and analyzed to determine the direction of the moisture flux in all grid points. The results showed that the IVT algorithm is a useful tool for finding the atmospheric moisture sources of the floods and the path of the moisture for the studied area. The algorithm indicated that the Red Sea and the Indian Ocean participate equally in moisture supply for the studied flood in the western parts of Iran. But for the peak rainfalls, the results showed that the center of the Red Sea is the main moisture source. For the selected heavy rainfall and flood in the southern parts of Iran, the IVT analysis showed that in days before the flood, the end edge of the moisture path is located over south of the Red Sea but in day with heavy rainfall and flood, the moisture path is limited to the Persian Gulf and surrounding area. So, it seems that for the selected southern flood, the main source of the moisture is the Persian Gulf and the other moisture sources have lesser role in the moisture supply.
    Keywords: Integrated Vapor Transport, flood, Red Sea, humidity source
  • Saeed Hadiloo, Saeid Mirzaei, Hosein Hashemi, Ebrahim Sefidari Pages 64-77
    The analysis of seismic facies is a technique for mapping geological features and properties using seismic data. To analyze seismic facies, seismic attributes are categorized and classified using machine learning algorithms to identify different seismic facies. Seismic facies analysis due to the nature of seismic data, which always has a degree of uncertainty, can produce different results with even small changes in input parameters of the analyzing method. For this reason, it is necessary to select the different stages of analysis, including the selection of input parameters and algorithm of machine learning, with high accuracy with regard to the objective of the seismic facies analysis. In this study, an interactive method with the supervision of interpreter is proposed for producing seismic facies map, using the optimal selection of the input parameters and the the proper selection of clustering and classification algorithms. In this method, the interpreter in a recursive and rotational process can compare the results of the analysis and generate thr optimal results by changing the input parameters. The method presented in this article is implemented in SeisART software. SeisART has a complete environment for data initialization (importing seismic data and well data). A user-friendly interactive environment allows the user to implement several methods and monitor the corresponding result in 2D and 3D.
    SeisART software makes the possibility of the interpreter contribution in the whole stages of seismic facies analysis procedure. The interpreter can select the input attributes and chose the proper methods of pattern recognition to reach the best possible result. In the software, various evaluation utilities have been provided in each stage of seismic facies analysis. These utilities allow the interpreter to monitor the results of each method quantitatively and qualitatively. In the unsupervised system, clustering quality factors are used. The interpreter calculates the validation indices for different methods of clustering and identifies the proper method which has been more successful in discovering the natural grouping of patterns in the data set. Afterward, if there is structural geology information about the horizon of interest, the interpreter can decide on the clustering result with more accuracy. In the supervised system, the most proper method is feasible using minimization of training data and validation data errors. In this case, the interpreter can use geological knowledge and well data information to verify obtained results. In this method, the interpreter can obtain different results by changing the input parameters. Comparing these results, and taking into account the path leading to this result, the interpreter gains more knowledge of existing facies.
    This method has been applied to the MSF4 horizons of the 3D seismic data of the North Sea F3 and has been shown which method is more efficient for different purposes
    Keywords: software, seismic facies analysis, seismic attributes, machine learning, fuzzy system, clustering
  • Saeed SoltaniMoghadam, Kamran Sepanloo, Myesam Kheyri Moloumeh Pages 78-95
    An earthquake of local magnitude 5.2 which occurred in Malard recently, was accompanied by four earthquakes with local magnitudes largerthan 4.0 These events can be considered as a manifestation of the state of Tehran, a mega-city with more than 13 million inhabitants, in last decade. For this reason, we tried to do a comprehensive seismic study on Tehran and surrounded area using all available data. In order to do this, we used three independent datasets belonged to Iranian Seismological Center (IRSC), International Institute of Earthquake Engineering and Seismology (IIEES) and Tehran Disaster and Mitigation Organization (TDMMO). We merged these three datasets to get a uniform one, using 3D cells gridding technique and then explosions or outliers were removed. The final dataset included 2000 events. Using 310 selected events based on horizontal and depth errors, azimuthal gap, RMS and minimum number of recording stations, we calculated 1D velocity model by utilizing Particle Swarm Optimization (PSO) method. The calculated model consists of three layers. The thickness and P velocity of layers are 4.0 km and 5.4 km/s, 6.0 km and 5.8 km/s and 5.0 km and 6.05 km/s for the first, second and third layer, respectively. The final layer is a half-space with a P velocity of 6.5 km/s. The computed Vp/Vs is 1.71. Then all events were relocated using our new velocity model utilizing fully non-linear probabilistic method to get as much as possible accurate locations. The results show that 40% of all relocated events have uncertainties less than 2.5 km and 5.0 km in horizontal and vertical direction, respectively. The final calculated mean RMS is ~ 0.24 s. In order to define the geometry of the active faults, three different subsets of events were selected based on their location uncertainties, azimuthal gap, RMS and minimum number of recording stations. This helped us to use well-located events for better defining of the fault traces on map view and in depth cross-sections. We plotted four cross-sections perpendicular to the strike of the main faults in the region. The focal mechanism solutions for 38 selected events were also computed based on P-wave first polarity method. The final results show that the eastern part of the study region is more active than the western part, at least in the last decade, and surprisingly, the Malard earthquake occurred in a region without any major activity from three months before the main shock. Stress field study also reveals that the maximum stress axis is N36E and the main seismotectonic regime is left lateral structures. These results are very consistent with the main trends of the North Tehran fault
    Keywords: probabilistic non-linear location, one-dimensional velocity model, focal mechanism, stress field, Tehran region, active faults
  • Meysam Amiri, Zahra Mousavi, Khalil Motaghi Pages 96-108
    On 20 December 2010, an earthquake with Mw 6.5 occurred in Rigan, a small town in the desert south of Bam city. The earthquake epicenter was in a low population area so, luckily, it caused only few casualties. Five days later 76 aftershocks reported by Iranian Seismological Center (ISC). On 27 January 2011, another earthquake (Mw 6.2) stroke an area at ~ 20 km southwest of the first earthquake. Bam earthquake Mw 6.6 occurred in 2003 with 40,000 victims is one of the deadliest earthquakes in Iran which is located in shear zones at southeast Iran. Considering the active faults distribution of the region and aftershocks of the 2010 Rigan earthquake encouraged us to better investigate and model the post-seismic deformation related to the 2010 earthquake. Post-seismic syudying provides information about rheology of the surrounding region and improves our knowledge about the strain release after the earthquake. In this study, COSMO-SkyMed (from Italian Space Agency, ASI) images spanning the temporal interval between 27 January 2011 and 15 July 2011 are used to investigate the post-seismic deformation following both earthquakes. We applied the Small Baseline Subset (SBAS) algorithm for images to obtain the post-seismic mean velocity map and the relative deformation time series. 109 interferograms, post-seismic mean velocity map and the relative deformation time series obtained Mean velocity map shows that displacements of post-seismic phase are right lateral strike slip same as co-seismic mechanism. Time series analysis reveals a clear post-seismic signal exponentially increasing with time until reaching the rate of more than 8 mm/year which indicates the end of post-seismic phase and following inter-seismic phase, starts with steady stress accumulation. Later, we modeled the post seismic signal considering a dislocation on a finite fault in an elastic and homogeneous half-space that are the assumptions for the Okada (1985) model. Post-seismic results modeled by adopting a two-step approach: (1) a non-linear inversion performed to constrain the fault geometry parameters and considering a uniform slip, then (2) a linear inversion performed to retrieve the slip distribution on the fault plane previously obtained. The fault plane is split into 1×1 km patches along strike and down-dip. Determining fault parameters and slip distribution by Okada model, indicates that the slip is concentrated in downdip of the coseismic depth with 1.2 m slip and also at the edges of the coseismic asperity. This slip distribution indicates that “afterslip” is the mechanism for post-seismic deformation of the Rigan earthquake.
    Keywords: Rigan earthquake, post-seismic phase, InSAR, mean deformation velocity, time series analysis
  • Sedighe Mehravar, Habib Alah Ghaedamini, Mohammad Jafar Nazemosadat Pages 109-126

    The Madden-Julian Oscillation (MJO) and El Niño Southern Oscillation (ENSO) are known as the primary modes of large-scale climate variability in the equatorial parts of the Indian and Pacific Oceans, respectively. The present study has made an effort to analyze interactions of these two oscillations during boreal autumn and the feedback of these interactions on the occurrence of dry or wet periods in Fars province. Cold and warm phases of ENSO (La Niña and El Niño, respectively) are then classified according to the size of the SOI. Then positive and negative phases of the MJO were defined according to the algebraic sign of the real-time multivariate MJO series 1 (RMM1). Daily values of October-December precipitation data of 9 synoptic stations of Fars province were obtained for the 1974-2013 period.
    Precipitation amount and the occurrence probability of precipitation were investigated for La Niña, El Niño MJO positive phase, MJO negative phase and other four combinations of these phases. Since the impacts of ENSO or MJO phases on precipitation variability had been studied previously, the present study mainly motivated to investigate precipitation characteristics for the episodes that El Nino or La Nina events were concurrent with the MJO positive phase (El-P or La-P, respectively) or MJO negative phase (El-N and La-N, respectively). The characteristics of atmospheric circulation and vapor transport over the ocean waters and the Middle East are then investigated to justify the obtained results.
    Some parameters including precipitation, intensity and occurrence probability of a rainy day were computed for all phases, separately. The ratio between these parameters during opposite phases was computed and spatially mapped to evaluate the effects of the interactions.
    The results indicated that the frequency of the positive and negative MJO phases is significantly associated with the ENSO condition. Differences between the frequency of the positive and negative MJO phases during the El Niño events were found to be insignificant. Contrariwise, during the El Niño events, the frequency and occurrence probability of the positive MJO phase are about twice than corresponding statistics during the MJO negative phase. It was found that ENSO and MJO have a significant influence on autumnal precipitation in Fars province. During El Niño or negative MJO phase, the mean occurrence probability and intensity of precipitation were significantly greater than the corresponding values during La Nina or positive MJO in Fars. For almost entire parts and for the periods of El Nino, precipitation has significantly enhanced or suppressed during the negative and positive MJO phase, respectively. On the other hand, such differences are not significant when the opposite phases of the MJO have concurred with the La Nina events. Compare to El-P, mean precipitation during El-N has increased by about 300% in the southeast and about 70% in the eastern districts. When the El-N phase has prevailed, the autumnal precipitation intensity has increased about 40-110% and 4-40% over the southern and northern area of Fars, respectively
    Keywords: El Ni?o, La Ni?a, Southern Oscillation, autumn rainfall, Fars province
  • Reza Borhani, Farhang Ahmadi, Givi Pages 127-146

    In general, the southwest Asia is one of the regions with positive anomalous values of tropopause folding frequency compared to the annual average of the northern hemisphere. The frequency of folding in this region in warm seasons is higher than in cold seasons. This research was carried out with the aim of statistical and dynamical analysis of the atmospheric processes associated with the tropopause folds in the southwest Asia during 2000–2015. Identification of tropopause folding is based on the algorithm developed by Sprenger et al. (2003) and Gray (2003) and refined by Škerlak et al. (2014). The detected folds are divided into three categories as shallow, medium, and deep based on their vertical extensions.
    The time series analysis of all types of tropopause folds shows that the frequency of folding events has an increasing trend during the period of study. The most frequent folding type is as shallow or medium in the summer season but as deep in the winter. The geographical distribution of the correlation coefficients between the monthly mean folding frequencies and some relevant dynamical quantities indicates that baroclinic instability mechanism plays the main role in the occurrences of tropopause foldings in the winter, while the effects of thermodynamic factors are dominant in the summer.Dynamical study of tropopause folds in both winter and summer seasons was conducted using January 2001 and 2004 as well as June 2007 and 2015 data sets. Results show that the winter tropopause foldings are associated with the formation of intense baroclinic waves in mid-levels of troposphere, strengthening of jet streams in upper levels and subsequently the formation of surface cyclones. Also, together with the seasonal displacement of the jet, the positions of tropopause folds are moving about 10 to 15 degrees latitudinally.The analysis of horizontal wave activity flux in January 2004, reveals the presence of a strong wave source (divergence of wave activity flux) in the Western Mediterranean and, at the same time, a wave sink (convergence of flux) over Europe. In January 2001, the wave activity flux was weakened and divided into two branches, one located in middle latitudes and the other in subtropical regions that transmitted wave energy to the southwest Asia. Furthermore, the strong equatorward wave propagation in this month indicates the anticyclonic wave breaking. In the two June months, as the baroclinic waves were weakened, the intensities of wave activity flux as well as the convergence and divergence centers in the southwest Asia were decreased compared to the January cases. Comparing the above results, it can be deduced that in winter, intense baroclinic wave packets in middle latitudes cause the strength of the subtropical jet, and consequently intensification of wave breaking which are associated with the occurrences of deep tropopause folds in the southwest Asia. In summer, the weakening of baroclinic activities leads to the reduction of deep folding frequency and the folds are formed mainly as shallow type at high levels.

    Keywords: tropopause folding, southwest Asia, potential vorticity, dynamical analysis, wave activity, baroclinic instability