per
انجمن ملی ژئوفیزیک ایران
مجله ژئوفیزیک ایران
2008-0336
2783-168X
2017-01-08
10
4
1
11
41686
Research Article
مطالعه لایه اکمن بهدامافتاده در بستر سواحل غرب تنگه هرمز
Study of arrested Ekman layer on the bottom of the western coastal areas of the Strait of the Hormuz
فرهنگ احمدی گیوی
ahmadig@ut.ac.ir
1
مرجانه حقیقی
2
عباسعلی علی اکبری بیدختی
bidokhti@ut.ac.ir
3
موسسه ژئوفیزیک دانشگاه تهران
واحد علوم و تحقیقات، دانشگاه آزاد اسلامی
موسسه ژئوفیزیک دانشگاه تهران
مطالعه اثرات دینامیکی اصطکاک بستر بر حرکت توده آب مستلزم آن است که ویژگیهای لایه مرزی موسوم به لایه اکمن شناسایی شود. لایه اکمن در بستر شیبدار شامل جریانهایی به نام شارش اکمن بهصورت فراشیب یا فروشیب است که این جریانها میتوانند سبب افزایش آمیختگی و تغییر در ضخامت لایه مرزی اکمن شوند. همچنین نیروی اصطکاک در بستر سبب کاهش سرعت جریان و به دنبال آن کاهش نیروی کوریولیس شده و در نتیجه جریان زمینگرد از بین میرود و فرایند چرخکاهی رخ میدهد. در این پژوهش به منظور شناخت ویژگیهای لایه اکمن در بستر شیبدار غرب تنگه هرمز به بررسی تاثیر شیب بستر بر توزیع پارامترهای فیزیکی آب و پدیدههای حاصل از آن پرداخته شده است. برای انجام این پژوهش، با استفاده از دادههای میدانیCTD مربوط به غرب تنگه هرمز در سالهای 1384و 1385، موجود در موسسه ملی اقیانوسشناسی، و بهکمک نرمافزارهای مرتبط، نیمرُخهای تغییرات دما، شوری، چگالی و سرعتهای کشندی با عمق و همچنین مقاطع عرضی چگالی ترسیم شدند. بررسی این نیمرُخها نشان میدهد که در بستر تنگه هرمز، جریان فراشیب است و پربندهای همچگالی با افزایش شیب بستر تقریبا به صورت عمود بر شیب قرار میگیرند. کمترین زمان چرخکاهی شارش اکمن دو ساعت محاسبه شده و مربوط به ایستگاهی است که بیشترین شیب را دارد. بهعلاوه، هرچه شیب بستر بیشتر و چینهبندی چگالی قویتر باشد، چرخکاهی شارش اکمن سریعتر اتفاق میافتد. بیشینه ضریب پخش پیچکی در بین ایستگاههای مورد مطالعه حدود m2s-1 005/0 بهدست آمده است و در ایستگاهی که ضریب پخش پیچکی آن از همه بزرگتر است، عدد اکمن نیز بهطور قابل ملاحظهای از سایر ایستگاهها بزرگتر است که نشاندهنده جریان ناپایا و آمیختگی ناشی از چینش لایه اکمن میباشد.
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.
https://www.ijgeophysics.ir/article_41686_46e3e9f123b34af0768c86e318197dc3.pdf
لایه اکمن
تنگه هرمز
عدد اکمن
ضریب پخش پیچکی
بستر شیبدار
پارامترهای فیزیکی
Ekman layer
Strait of Hormuz
Ekman number
eddy diffusion coefficient
sloping bottom
per
انجمن ملی ژئوفیزیک ایران
مجله ژئوفیزیک ایران
2008-0336
2783-168X
2017-01-08
10
4
12
27
41687
Research Article
مطالعه بیهنجاری دمای زمستان ایران با استفاده از دادههای بازتحلیل NCEP/NCAR
A study of winter temperature anomalies in Iran by using the NCEP/NCAR reanalysis dataset
پرویز ایران نژاد
piran@ut.ac.ir
1
فرهنگ احمدیگیوی
2
نسیم نیکویی
nasim.nikoo@yahoo.com
3
مؤسسه ژئوفیزیک دانشگاه تهران
موسسه ژئوفیزیک دانشگاه تهران
موسسه ژئوفیزیک دانشگاه تهران
در این پژوهش به بررسی تغییرپذیری زمانی و مکانی و بیهنجاری دمای هوا در ارتفاع دو متر از سطح زمین در فصل زمستان در سه مقیاس شامل کل ایران، شش منطقه مشخص شده بر اساس طول و عرض جغرافیایی و در مقیاس یاخته شبکه دادههای بازتحلیل NCEP/NCAR بهصورت میانگین ماهانه در سالهای 1950 تا 2010 پرداخته شده است. نخست دادههای بازتحلیلشده دما از نظر کیفیت، بررسی و سالهای پیش از 1967 بهدلیل کیفیت پایین از مطالعه حذف شد. بررسیها نشان داد که تغییرات میانگین دمای زمستان ایران سنجه مناسبی برای بیان تغییرات بینسالانه میانگین زمستانی دمای مناطق مختلف است. از 44 سال مطالعه شده، بزرگترین بیهنجاریهای مثبت دما در زمستانهای 1970، 1979، 1999 و 2010 و بزرگترین بیهنجاریهای منفی در زمستانهای 1972، 1973، 1989 و 2008 رخ داده است. بررسی میانگین الگوی بیهنجاری میدانهای ارتفاع ژئوپتانسیلی، دما و باد در زمستانهای بیهنجار گرم و سرد در حوزه بین 10 تا 60 درجه شمالی و 10 درجه غربی تا 80 درجه شرقی نشان داد که مراکز بیشینه بیهنجاری این کمیتها در سالهای گرم و سرد ایران کاملاً عکس یکدیگرند. دو مرکز بیشینه بیهنجاری ارتفاع ژئوپتانسیلی در سالهای بیهنجار سرد و گرم ایران، یکی در منطقه شمالغرب اروپا و دیگری در شمالشرق دریای کاسپین قرار دارند. بر اساس تفاوت بیهنجاری ارتفاع ژئوپتانسیلی تراز 250 هکتوپاسکال در این دو مرکز، شاخصی برای بیان بیهنجاری دمای زمستان در ایران معرفی شد. تحلیل همبستگی نشان داد که همبستگی مثبت معنیدار بین میانگین دمای زمستان ایران و همچنین میانگین دمای زمستان مناطق مختلف کشور و شاخص معرفی شده وجود دارد. بهدلیل شباهت الگوی دورپیوند بهدست آمده در این پژوهش با دو دورپیوند EA/WR و NCP همبستگی بین شاخصهای این دو دورپیوند و دمای زمستان ایران تحلیل شد. نتایج نشاندهنده توانایی بالاتر شاخص معرفیشده نسبت به دو شاخص دیگر برای بیان دمای زمستان ایران است.
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.
https://www.ijgeophysics.ir/article_41687_055a8759703e69951aacd3dcec2b9c78.pdf
ایران
اروپا
دریای کاسپین
بیهنجاری دما
بیهنجاری ارتفاع ژئوپتانسیلی
دورپیوند
EA/WR
NCP
Iran
Europe
Caspian Sea
temperature anomaly
geopotential height anomaly
teleconnection
EA/WR
NCP
per
انجمن ملی ژئوفیزیک ایران
مجله ژئوفیزیک ایران
2008-0336
2783-168X
2017-01-08
10
4
28
47
41688
Research Article
تحلیل طیفی خشکسالیهای ایران
A spectral analysis of Iran's droughts
حجت اله دانشمند
daneshmand@hamoon.usb.ac.ir
1
پیمان محمودی
p_mahmoudi@gep.usb.ac.ir
2
دانشگاه سیستان و بلوچستان
دانشگاه سیستان و بلوچستان
هدف اصلی این مطالعه آشکارسازی رفتار دورهای خشکسالیهای ماهانه ایران است. برای دستیابی به این هدف از آمار بارش ماهانه 41 ایستگاه همدیدی کشور برای یک دوره آماری 31 ساله (2013–1983) و برای شناسایی فراوانی درجات مختلف خشکسالیها از شاخص خشکسالی مؤثر (EDI) استفاده شد. هدف اصلی این شاخص اختصاص یک ارزش عددی به هر رویداد بارندگی در مقیاس زمانی ماهانه میباشد. در تحلیل طیفی و شناسایی دورههای غالب موجود در سریهای زمانی خشکسالی ایستگاههای مورد مطالعه از تبدیل فوریه سریع (FFT) تقسیمات بسامدی با مبنای دو استفاده شد. نتایج نشان داد که دورههای غالب موجود در سریهای زمانی خشکسالیهای ایران بسیار متنوع بوده و از دورههای 2 تا 30 ساله را شامل میشوند. بهطور کلی برای دورههای بلندمدت، میتوان در کل ایران دو دوره 10 ساله و 30 ساله را بهعنوان دورههای غالب در نظر گرفت. در دورههای بلندمدت، فراوانی این دو دوره در مولفه اول بیشینه است و در دو مولفه دیگر کاهش مییابد. اما درخصوص دورههای کوتاهمدت، نمیتوان برای کل ایران یک دوره غالب مشخص کرد. نکته حائز اهمیت در مورد دورههای کوتاه مدت عدم وجود دورههای درونسالی، به استثنای یک ایستگاه (ایستگاه بابلسر)، در کل ایستگاههای مورد مطالعه است.
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.
https://www.ijgeophysics.ir/article_41688_4d51fd170bf270eda0068f119db17ab6.pdf
خشکسالی
شاخص خشکسالی مؤثر
تحلیل طیفی
توان
بسامد
Drought
effective drought index
Spectral analysis
Cycle
Iran
per
انجمن ملی ژئوفیزیک ایران
مجله ژئوفیزیک ایران
2008-0336
2783-168X
2017-01-08
10
4
48
61
41689
Research Article
مطالعه ساختار سنگکره در منطقه برخوردی زاگرس شمالی با استفاده از وارونسازی همزمان توابع گیرنده و منحنیهای پاشش امواج سطحی
Lithosphere structure in the north Zagros collision revealed by joint inversion of P receiver function and surface wave dispersion
فروغ کلوندی
foroghka@gmail.com
1
سید خلیل متقی
kmotaghi@iasbs.ac.ir
2
اسماعیل شبانیان
shabanian@iasbs.ac.ir
3
دانشگاه تحصیلات تکمیلی علوم پایه زنجان
دانشگاه تحصیلات تکمیلی در علوم پایه زنجان
دانشگاه تحصیلات تکمیلی علوم پایه زنجان
مطالعه ساختار سرعتی هر منطقه کمک شایانی به شناخت وضعیت لرزه زمینساختی آن میکند. تعیین پارامترهایی چون ضخامت پوسته و سنگکره، مکانیابی دقیق زمینلرزهها و پیشبینی تحولات آینده زمینساختی هر منطقه در گرو داشتن مدل سرعتی موثق از منطقه است. فلات ایران گستره وسیعی بین دو صفحه عربی و اوراسیا است و دگرشکلی آن ناشی از همگرایی این دو صفحه است. در این میان، کمربند کوهزایی زاگرس یکی از فعالترین و جوانترین کوهزادهای موجود در مرزهای همگرایی قارهای در جهان است که شناخت ویژگیهای ساختاری این منطقه میتواند به فهم نحوه تغییر شکل در مرحله اول برخورد قارهای و کوهزایی کمک کند. در این تحقیق ساختار سنگکره در پهنه برخوردی زاگرس شمالی، با استفاده از وارونسازی همزمان توابع گیرنده و منحنیهای پاشش امواج سطحی مورد مطالعه قرار میگیرد. منحنیهای پاشش از مطالعهای که در سال 2014 از طریق بُرشنگاری (توموگرافی) رحیمی و همکاران (2014) انجام دادند، برگرفته و توابع گیرنده از 161 دورلرزه با فاصله رومرکزی °30 تا °95 و بزرگای بیش از 5 محاسبه شده است. دورلرزهها را 38 ایستگاه لرزهنگاری موقت در امتداد پروفایلی به طول 400 ~ کیلومتر در حد فاصل شهرهای ایلام تا قم ثبت کردهاند. در این مطالعه با استفاده از فرایند وارونسازی همزمان توابع گیرنده و منحنیهای پاشش، مدلهای سرعتی یکبُعدی موج برشی در زیر هر ایستگاه بهدست میآید. سپس، از کنار هم قرار دادن این مدلهای یکبُعدی، یک مدل دوبعدی سرعت برای سنگکره زیرین ایستگاههای شبکه حاصل میشود. در این مدل یک زبانه کمسرعت در درون پوسته در مدل سرعتی قابل تشخیص است که از محل گسل اصلی عهد حاضر شروع میشود و با فاصله گرفتن از این گسل در راستای شمال شرقی به عمقهای بیشتر میرود. این بیهنجاری در فاصله 200 کیلومتری شمال شرق گسل در عمق 35 ~ کیلومتر قابل مشاهده است. مدل سرعتی نشان میدهد که مقدار ضخامت پوسته از ابتدای خط اندازهگیری در بخشهای جنوبی زاگرس (محدوده ایلام)، 43 کیلومتر است که با حرکت به سمت شمال شرق به میزان 57 کیلومتر در زیر گسل اصلی عهد حاضر میرسد. در سنندج-سیرجان و ارومیه-دختر بهتدریج بر این ضخامت افزوده شده و در مرز این دو ناحیه به بیشینه مقدار 62 کیلومتر میرسد. در زیر ایران مرکزی میزان ضخامت پوسته کاهش یافته و در انتهای خط اندازهگیری در زیر ایران مرکزی این ضخامت به 42 کیلومتر میرسد. مدل سرعتی بهدست آمده، همچنین، اطلاعات خوبی از بخش سنگکرهای گوشته ارائه میدهد. این مدل نشان میدهد که سنگکره پُرسرعت زیرین زاگرس شمالی در زیر سنندج-سیرجان و ارومیهدختر و لبه جنوبی ایران مرکزی گسترش یافته است. چنین مشاهدهای در تطابق با زیرراندگی مشاهده شده در پوسته است و میتواند شاهدی بر زیرراندگی بلوک عربی به زیر ایران مرکزی باشد.
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
https://www.ijgeophysics.ir/article_41689_9a4a3c692d3442ae39fe62c90020c069.pdf
توابع گیرنده
منحنیهای پاشش موج ریلی
وارونسازی همزمان
زیرراندگی
زاگرس شمالی
receiver functions
surface wave dispersion
Joint inversion
velocity structure
north Zagros collision
per
انجمن ملی ژئوفیزیک ایران
مجله ژئوفیزیک ایران
2008-0336
2783-168X
2017-01-08
10
4
62
74
41690
Research Article
مدلهای جهانی برای بررسی شکوفایی گیادروایانی (فیتوپلانکتونی) در دریای عمان و شمالغربی دریای عرب
Global models for investigation of phytoplankton blooms in the Gulf of Oman and the northwest of Arabian Sea
سیده صفورا صدیق مروستی
safoora.seddigh@gmail.com
1
آناند نانادسیکان
2
عباسعلی علی اکبری بیدختی
bidokhti@ut.ac.ir
3
سرمد قادر
sghader@ut.ac.ir
4
محمد صدیق مرتضوی
mseddiq1@yahoo.com
5
دانشگاه آزاد اسلامی، واحد علوم و تحقیقات
دانشگاه جانزهاپکینز
موسسه ژئوفیزیک دانشگاه تهران، تهران
موسسه ژئوفیزیک دانشگاه تهران، تهران
پژوهشکده اکولوژی خلیج فارس و دریای عمان
در سالهای اخیر، افزایشی در شکوفایی کشند قرمز در شمال غربی دریای عرب و دریای عمان مشاهده شده است که این سؤال را در ذهن به وجود میآورد: آیا تغییرات اقلیمی باعث افزایش این روند شده است یا خیر؟ با هدف پاسخ به این سوال، در این پژوهش از مدلهای مختلف سامانه زمین و دادههای ماهوارهای استفاده شده است. دادههای ماهوارهای دو شکوفایی را در این منطقه نشان میدهند، شکوفایی زمستانی که بیشینه آن در ماه فوریه و شکوفایی تابستانی که بیشینه آن در ماه سپتامبر است. تغییرات درونسالانهای نیز در شکوفایی زمستانی حاصل از اثر پیچکهای چرخندی که توزیع مکانی متفاوتی از یک سال به سال دیگر دارند، مشخص شده است. دو مدل با تفکیک کم (°1) با بخش زیستزمینشیمیایی تقریباً پیچیده (مدل TOPAZ) چرخه سالانه را نشان میدهند ولی قادر به نمایش پیچکها و تغییرات درونسالانه نیستند. مدلهای با تفکیک بیشتر (GFDL CM2.6) که توانایی مدلسازی پیچکها را دارند همراه با بخش زیستزمینشیمیایی سادهتر (مدل miniBLING) تغییرات درونسالانه بزرگتری را نشان میدهند، اما مقدار شکوفایی زمستانی را بیش از اندازه پیشبینی میکنند. این مدل اگرچه رابطه بین شکوفایی و پیچک را در بخش جنوبی بهخوبی نشان میدهد، اما در بخش شمالی منطقه موفق نیست. این امر بهدلیل نداشتن توانایی در مدلسازی دماشیب (ترموکلاین) و غذاشیب (نوتریکلاین) قوی و درست در آن مناطق میباشد.
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.
https://www.ijgeophysics.ir/article_41690_8fbfcde4a42464f29ac3ac50321e0a91.pdf
کشند قرمز
شکوفایی
مدلسازی عددی
پیچک
مدل سامانه زمین
زیستزمینشیمیایی
red tide
bloom
Numerical modeling
eddy
Earth system models
biogeochemistry
per
انجمن ملی ژئوفیزیک ایران
مجله ژئوفیزیک ایران
2008-0336
2783-168X
2017-01-11
10
4
75
84
41735
Research Article
برآورد تولید اُزون سطح زمین و دودمِه نورشیمیایی با استفاده از غلظت NOx وO3 در اصفهان
Estimation of ground level ozone and photochemical smog production using concentrations of NOx and O3 in Isfahan
فریده کنعانی
faridehca743@hotmail.com
1
مصطفی فضلی
mfazli@semnan.ac.ir
2
ویکتوریا عزتیان
3
دپارتمان شیمی، دانشگاه سمنان
دپارتمان شیمی، دانشگاه سمنان
سازمان هواشناسی شهر اصفهان
اُزون سطح زمین در طی فرایندهای نورشیمیایی اکسیژن هوا در حضور نور خورشید و NOx بهعنوان یکی از مهمترین عوامل موثر، تولید میشود. در این پژوهش به بررسی اُزون در شهر اصفهان در آذرماه سال 1388 پرداخته شده، که بنابر مشاهدات سازمان هواشناسی، در طی این بازه دودمِه نورشیمیایی بر فراز اصفهان قابل مشاهده بود. همه بررسیها در ساعتهای تابش آفتاب (15-6 محلی) انجام گرفته است. در این پژوهش با استفاده از روش تقریب حالت پایا رابطه سرعت دیفرانسیل تابعی از سه متغیر NO،NO2 و O3 برای بررسی تغییرات غلظت اُزون بهدست آمده است. انرژی فعالسازی، مرتبه واکنش و سازوکار تشکیل واکنش اُزون در وردسپهر مطالعه و محاسبه شدهاند . NOو O3 سطح زمین بهعنوان شاخصی برای تشکیل دودمه نورشیمیایی درنظر گرفته میشوند، لذا تولید دودمِه بر مبنای تغییرات غلظت NO وO3 مورد توجه قرار گرفت. بررسی روابط تولید دودمِه برمبنای غلظتها و بررسی براساس زمان، نشاندهنده پیروی تولید دودمه از روابط غیرخطی و درجه دوم است. بنا به معادله آهنگ تغییر زمانی اُزون، d[O3]/dt و نتایج انرژی فعالسازی میتوان تولید اُزون را در شهر اصفهان از واکنشهایی که دارای سینتیک کلی مرتبه سوم هستند، در نظر گرفت. از طرف دیگر افزایش غلظت NOx در سطح زمین در اثر واکنشهای نورکافتی (فوتولیز) افزایش غلظت اُزون را سبب گردیده است.
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+O2+ hν→NO+O3 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.
https://www.ijgeophysics.ir/article_41735_5150136e6efbbdf76468b5e63a223b11.pdf
اُزون وردسپهری
سرعت واکنش
انرژی فعالسازی
مدل نورشیمیایی
الگوریتم تولید دودمِه
tropospheric ozone
ozone reaction rate
activation energy
photochemical model
smog production algorithm
per
انجمن ملی ژئوفیزیک ایران
مجله ژئوفیزیک ایران
2008-0336
2783-168X
2017-01-11
10
4
85
96
41736
Research Article
شناسایی گسلها در دادههای لرزهای بازتابی به کمک نشانگر سوبل
Faults detection in seismic reflection data using Sobel attribute
فاطمه اعظمپور
1
امین روشندل کاهو
roshandel@shahroodut.ac.ir
2
علیرضا احمدیفرد
3
علیرضا حیدریان
4
دانشکده مهندسی معدن، نفت و ژئوفیزیک، دانشگاه صنعتی شاهرود
دانشکده مهندسی معدن، نفت و ژئوفیزیک، دانشگاه صنعتی شاهرود
دانشکده مهندسی برق و رباتیک، دانشگاه صنعتی شاهرود
شرکت نفت فلات قاره ایران
گسلها یکی از مهمترین رویدادهای زمینشناسی هستند که در اکتشافات معدنی، هیدروکربنی و مطالعات زمینشناختی از اهمیت بسیار زیادی برخوردارند. گسلها بهمنزله یک تله ساختاری هیدروکربنی میباشند و از اینرو شناخت و تعیین محل دقیق آنها عامل مهمی در توصیف مخازن هیدروکربنی است. در دادههای لرزهای بازتابی، گسلها بهصورت ناپیوستگی در امتداد بازتابندهها و تغییر دامنه مشاهده میشوند. بهدلیل ماهیت دادههای لرزهای، شناسایی صفحه گسل در این دادهها، فرایندی دشوار است و تعیین محل دقیق آنها بهتنهایی از روی داده لرزهای بهندرت امکانپذیر است. نشانگرهای لرزهای، ابزاری هستند که اطلاعات نهفته در دادههای لرزهای را آشکار مینمایند و امکان تفسیر کمّی را در اختیار مفسر لرزهنگاری فراهم میآورند. در این مقاله، صافی سوبل دوبُعدی و سهبُعدی بهمنزله یک نشانگر لرزهای برای شناسایی و تعیین موقعیت گسلها در داده لرزهای بازتابی معرفی میشود. این صافی، یکی از ابزارهای پرکاربرد در فرایند تشخیص لبه در پردازش تصاویر است که گرادیان تصویر را محاسبه میکند و به تغییرات دامنه حساس است. بهمنظور بررسی عملکرد این روش در شناسایی گسلها، نشانگر سوبل بر روی دادههای لرزهای مصنوعی و واقعی در دو حالت دوبعدی و سهبعدی مورد بررسی قرار میگیرد و نتایج حاصل با نتایج نشانگر همدوسی مقایسه میشود. بررسی نتایج بهدست آمده نشان میدهد که نشانگر سوبل میتواند بهمثابه یک نشانگر مناسب برای شناسایی گسلها در دادههای لرزهای بازتابی مورد استفاده قرار بگیرد.
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.
https://www.ijgeophysics.ir/article_41736_901408851e654c44361f99d2e8aabe28.pdf
لرزهنگاری بازتابی
گسل
تشخیص لبه
نشانگر سوبل
reflection seismology
Fault
Edge detection
Sobel attribute
per
انجمن ملی ژئوفیزیک ایران
مجله ژئوفیزیک ایران
2008-0336
2783-168X
2017-01-11
10
4
97
111
41737
Research Article
برآورد ضریب کیفیت موجهای برشی و فشاری در پوسته شمال غرب ایران
Estimation of compressional and shear wave quality factor in North West of Iranian Plateau
مجتبی نقوی
mojtabanaghavi@ut.ac.ir
1
حبیب رحیمی
rahimih@ut.ac.ir
2
علی مرادی
asmoradi@ut.ac.ir
3
مؤسسه ژئوفیزیک، دانشگاه تهران
موسسه ژئوفیزیک دانشگاه تهران
مؤسسه ژئوفیزیک، دانشگاه تهران
در مطالعات مرتبط با زلزلهشناسی مهندسی برای برآورد خطر لرزهای در مناطق مختلف و همچنین تعیین دقیق بزرگا، شبیهسازی جنبش نیرومند زمین و مطالعه انرژی مخرب حاصل از زمینلرزه در حوزه نزدیک تا متوسط، تخمین ضریب کیفیت امواج برشی و فشاری نقشی اساسی دارد. بدین منظور با توجه به قرار گرفتن چندین شهر بزرگ در پهنه شمال غربی فلات ایران، ضریب کیفیت امواج فشاری و برشی با استفاده از دادههای هفده ایستگاه لرزهنگاری با بیش از 13000 زمینلرزه ثبت شده در شبکه مرکز لرزهنگاری کشوری (IRSC) و شبکه ملی لرزهنگاری نوارپهن ایران (INSN) برآورد شد. برای سه گروه داده بهترتیب با فاصله رومرکزی کمتر از 100 کیلومتر، از 100 تا 200 کیلومتر و صفر تا 200 کیلومتر بررسی تغییرات جذب امواج حجمی در نُه نوار بسامدی، با بسامدهای مرکزی در 3، 5، 7، 10، 14، 20، 28، 38 و 47 هرتز انجام گرفت و برای هرکدام از ایستگاهها برای بسامدهای مختلف بهطور جداگانه ضریب کیفیت تخمین زدهشد. در این مطالعه از روش کُدای بهنجارشده برای برآورد ضریب جذب امواج فشاری و برشی حاصل از زمینلرزه بهعنوان روشی قابل اعتماد با توجه به حذف شدن اثرات چشمه و ساختگاه در مدلسازی انجام گرفته، استفاده شد. برای پهنه شمال غربی ایران مقادیر وابستگی بسامدی ضریب کیفیت امواج فشاری و برشی در همه ایستگاهها برآورد شد بهطوریکه مقدار متوسط آن برای محدوده شمال غرب فلات ایران از و پیروی میکند. با توجه به مقادیر پایین و لذا جذب بالای امواج برشی و فشاری در پهنه شمال غرب ایران دامنه امواج حین عبور از زمین بهشدت تضعیف میشود که این اثر جذب امواج لرزهای، باعث کاهش خسارات ناشی از زمین لرزهها در فواصل مناسب از گسلها در زمان وقوع زمینلرزه خواهد شد.
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.
https://www.ijgeophysics.ir/article_41737_fc262dec83ac697b0df5e9381edaf75d.pdf
موجهای برشی و فشاری
شمال غرب فلات ایران
جذب
ضریب کیفیت
Compressive and shear waves
NW of Iranian Plateau
Attenuation
Quality factor
per
انجمن ملی ژئوفیزیک ایران
مجله ژئوفیزیک ایران
2008-0336
2783-168X
2017-01-11
10
4
112
127
41738
Research Article
برآورد ساختار سرعت موج برشی و تغییرات عمق موهو در جنوبشرق ایران، با استفاده از وارونسازی همزمان پاشش امواج سطحی و تابع گیرنده موج P
Estimation of shear wave velocity structure and Moho depth variations in the southeast of Iran, based on joint inversion of surface waves dispersion and P-wave receiver function
مهدی راستگو
rastgoo@ut.ac.ir
1
حبیب رحیمی
rahimih@ut.ac.ir
2
حسین حمزه لو حمزه لو
3
وحید غلامی
va.gholami@gmail.com
4
مؤسسه ژئوفیزیک دانشگاه تهران
مؤسسه ژئوفیزیک دانشگاه تهران
پژوهشگاه بینالمللی زلزلهشناسی و مهندسی زلزله
تهران، ایران
در ناحیه جنوبشرق ایران که شامل بخش غربی زون فرورانش فعال مکران است به دلایل مختلفی همچون لرزهخیزی اندک و نبود تراکم ایستگاههای لرزهنگاری دائم، مطالعات اندکی در رابطه با ساختار سرعت موج برشی انجام شده است. هدف از این پژوهش، برآورد ساختار سرعت موج برشی و تغییرات عمق موهو در زیر چهار ایستگاه لرزهنگاری نوارپهن ZHSF، KRBR، CHBR و BNDS واقع در محدوده جنوبشرق ایران است. برای این منظور از روش وارونسازی همزمان منحنیهای پاشش امواج سطحی حاصل از زمینلرزههای محلی ناحیه ریگان، به همراه تابع گیرنده موج P حاصل از لرزهنگاشتهای زمینلرزههای دورلرز ثبت شده در چهار ایستگاه لرزهنگاری ذکر شده، استفاده میشود. عمق موهو بهدست آمده از این روش برای چهار ایستگاه ZHSF، KRBR، CHBR و BNDS بهترتیب برابر است با 4±38، 6±46، 2±26 و 5±56 کیلومتر. کمعمق بودن موهو و نیز بالاتر بودن مقادیر سرعت موج برشی بهدست آمده در زیر ایستگاه CHBR با فرورانش کمعمق لایه اقیانوسی سرعت بالای صفحه عربی به زیر حاشیه جنوبی مکران مطابقت دارد. در عین حال ضخامت زیاد پوسته در محدوده ایستگاه BNDS، در تطابق با توپوگرافی مکان این ایستگاه و نیز برخورد پوسته ضخیم قارهای صفحه عربی با زاگرس در این ناحیه میباشد.
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.
https://www.ijgeophysics.ir/article_41738_52b46def5d85eff4b9cf05f61b6809ec.pdf
جنوبشرق ایران
سرعت موج برشی
عمق موهو
پاشش امواج سطحی
تابع گیرنده موج P
Southeast of Iran
Shear wave velocity
Moho Depth
surface waves dispersion
P-Wave Receiver Function
per
انجمن ملی ژئوفیزیک ایران
مجله ژئوفیزیک ایران
2008-0336
2783-168X
2017-01-11
10
4
128
141
41739
Research Article
محاسبه کاهندگی طیف دامنه جنبش زمین در ناحیه شمالغرب ایران
Attenuation of ground-motion spectral amplitudes in the NW Iran
خلیل متقی
kmotaghi@iasbs.ac.ir
1
زهرا ضرونیزاده
2
عبدالرضا قدس
aghods@iasbs.ac.ir
3
دانشگاه تحصیلات تکمیلی علوم پایه زنجان
دانشگاه تحصیلات تکمیلی علوم پایه زنجان
دانشگاه تحصیلات تکمیلی علوم پایه زنجان
در این مقاله با تلفیق دادههای ثبتشده توسط 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 کیلومتر و بیش از آن به هم شبیه است و مدل کاهندگی را از حالت دوخطی به حالت یکتکه تبدیل میکند. با استفاده از ضریب کاهندگی جذب ذاتی، تابع کیفیت موج برشی بهصورت معرفی میشود.
در این مقاله با تلفیق دادههای ثبتشده توسط 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 کیلومتر و بیش از آن به هم شبیه است و مدل کاهندگی را از حالت دوخطی به حالت یکتکه تبدیل میکند. با استفاده از ضریب کاهندگی جذب ذاتی، تابع کیفیت موج برشی بهصورت معرفی میشود.
واژههای کلیدی: رابطه کاهندگی، تابع کیفیت، شمالغرب ایران
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. +0.09 for Central Alborz: Motaghi and Ghods, 2012; +0.2 for North Iran: Motazedian, 2006; +0.2 for SE Canada and the NE United States: Atkinson, 2004; +0.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. +0.09 for Central Alborz: Motaghi and Ghods, 2012; +0.2 for North Iran: Motazedian, 2006; +0.2 for SE Canada and the NE United States: Atkinson, 2004; +0.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 .
https://www.ijgeophysics.ir/article_41739_cbb454c57197e5de322d151a85504361.pdf
رابطه کاهندگی
تابع کیفیت
شمالغرب ایران
seismic attenuation relation
Quality factor
NW Iran