@article { author = {متقی, سیدخلیل and قدس, عبدالرضا and سیاهکوهی, حمیدرضا}, title = {Visual Presentation of Seismic Attenuation Curves in the Tehran Region and Moho Thickness Estimation through the Curves}, journal = {Journal of the Earth and Space Physics}, volume = {36}, number = {2}, pages = {-}, year = {2010}, publisher = {Institute of Geophysics, University of Tehran}, issn = {2538-371X}, eissn = {2538-3906}, doi = {}, abstract = {The central Alborz Mountains of northern Iran (Figure 1) is a zone of intense active continental deformation and is the most heavily populated region of Iran, including the mega-city of Tehran with a population of over 10 million. The region is affected by numerous active faults, some of them having great seismic potential, and documented historical seismicity. Hence, the central Alborz is a prime candidate for seismic hazard evaluation studies. The choice of ground-motion model has a significant impact on risk estimates for a distance range of 50 to 200 km within an active seismic zone such as the central Alborz. To visualize the attenuation behavior of calculated amplitude and to find a best fit attenuation curve, we used the statistical LOcally-WEighted regression Scatterplot Smoother (LOWESS) robust curve-fitting method of Cleveland et al. (1992). The method requires no assumptions regarding the underlying form of the relationship between amplitude of acceleration spectrum and distance. The use of a theoretical function can distort the fit by forcing a number of hinge points rather than by the natural hinge points in the data trend. We only used 220 seismograms from 22 selected events with azimuthal gap less than 180 degrees. The data was bounded by azimuthal gap constrain to decrease successfully the radiation pattern effects on obtained spectral amplitudes. Two significant hinges are observed in 106 and 191 km. The first hinge is due to SmS phases received from Moho discontinuity beyond 106 km. The second hinge is related to the distance at which surface wave propagation becomes dominant. A small hinge was recognized in 55 km. We believe that it is because of a shallower discontinuity, possibly the Conrad discontinuity. To find the reasons behind the existence of the hinges, we model reflection coefficient versus distance for an Earth model having the Conrad or Moho discontinuity. We can estimate the depth of discontinuities using the calculated distances for the hinges (55 and 106 km). Using a simple model and mentioned hinges (55 and 106 km), the thickness of the Conrad and Moho are estimated as 24.5 and 46.5 km, respectively. The Moho thickness agrees well with that obtained from 1D tomography of the study region.}, keywords = {Conrad discontinuity,Moho discontinuity,Robust Lowess algorithm,Seismic attenuation curve}, title_fa = {نمایش شکل منحنی کاهندگی لرزه‌ای در ناحیه تهران و برآورد ضخامت موهو از روی آن}, abstract_fa = {تعیین شکل منحنی کاهندگی در یک ناحیه ناشی از ناکشسان بودن زمین و گسترش هندسی موج در چندین دهه اخیر همواره مورد توجه بوده است. با استفاده از دستگاه‌های سرعت‌نگاشت شبکه لرزه‌نگاری رقومی مؤسسة ژئوفیزیک دانشگاه تهران، 220 نگاشت زلزله با بزرگی بین 2ر3 تا 9ر3، رویداده در ناحیه تهران بررسی و شکل منحنی کاهندگی در این ناحیه و در بسامد‌های متفاوت با استفاده از الگوریتمRobust Lowess عرضه می‌شود. از روی شکل منحنی کاهندگی و مقایسه این منحنی‌ها در بسامد‌های گوناگون، اثرات ناپیوستگی‌های مهم سنگ‌سپهر نشان داده می‌شود و نقاط تغییر ضریب گسترش هندسی که ناشی از وجود این ناپیوستگی‌ها است در 106 و 191 کیلومتر به‌دست می‌آید. در پایان از روی محل شکستگی‌ها در منحنی‌های کاهندگی اندازه ضخامت ناپیوستگی‌های کُنراد برابر با 5ر24 کیلومتر و موهو 5ر46 کیلومتر برآورد می‌شود.}, keywords_fa = {الگوریتم Robust Lowess,منحنی کاهندگی دامنه موج لرزه‌ای,ناپیوستگی کُنراد,ناپیوستگی موهو}, url = {https://jesphys.ut.ac.ir/article_21552.html}, eprint = {https://jesphys.ut.ac.ir/article_21552_db52cb1f2a441dea5f8aa559156523d0.pdf} } @article { author = {Azadmanesh, A and Gheitanchi, Mohammad Reza}, title = {Source Parameters of The 26 April 2007 Bam Earthquake Obtained By The Near Field Accelerograms Data}, journal = {Journal of the Earth and Space Physics}, volume = {36}, number = {2}, pages = {-}, year = {2010}, publisher = {Institute of Geophysics, University of Tehran}, issn = {2538-371X}, eissn = {2538-3906}, doi = {}, abstract = {On 26 December 2003 at 05:26LT (01:26GMT), a catastrophic earthquake Mw 6.6, seismic moment (6 – 9 ×10 Nm) struck the city of Bam in the Kerman province of southeast Iran. The intense shaking in the city caused the complete collapse of nearly every building in the central parts of the city including many of the newer buildings, killing about 40000 people officially. The city lies to the east of the Nyband – Gowk – Sarvestan fault system and on which several powerful earthquakes have occurred over the past 23 years [Berberian et al., 1984; Berberian and Qorashi, 1994; Berberian et al., 2001]. There are no recorded historical earthquakes at Bam. Most of the citadel of Arg-e-Bam, one of world heritage sites inscribed by UNESCO, which was constructed by mud brick about 2000 years ago was destroyed in this earthquake. Four powerful earthquakes (M>6) since 1981 have occurred near the Gowk fault zone [Berberian et al., 1984, 2001]. The Gowk fault zone, a predominantly right–lateral strike–slip zone that extends from 50 km west of Bam northward [Walker and Jackson, 2002], has also been associated with several large historical earthquakes [Ambraseys and Melville, 1982; Berberian and Yeats, 1999]. After 26 December 2003 many aftershocks and earthquakes have occurred in Bam, southeast Iran. The last event with M=4.9 (NEIC) occurred on 26 April 2007. In this study, we calculated the ground motion parameters of the latest event using digital accelerograms of the Building and Housing Research Center (BHRC). We calculated seismic Moment (Mo) and Moment magnitude (Mw) of the 03/26/2007 Bam earthquake. For acceleration, Velocity and Displacement diagrams; we determined flat portion (K) and corner frequency (fc) by correction of base line and band pass filtering of Fourier spectrum. We calculated average seismic Moment (Mo=3.6×1016 Nm) and Mw 5; the answers conform to the results by other methods computation. We determined the location of the epicenter, our motivation for these computations was the determination of the difference of epicentral location with Maximum 36 km. In this study, we determined the focal mechanism of this earthquake; strike angle (N51?W), Dip89 ? westwards and Rake 138?, by using first P motion of the accelerograms and seismograms Data. The occurrence of a strong earthquake such as that of the 26 April 2007 indicates that the Bam region is still seismically active.}, keywords = {Accelerograms Data,Bam earthquake,Earthquake location,Focal mechanism,Moment magnitude,seismic moment,Source parameters}, title_fa = {محاسبه پارامترهای منبع زمین‌لرزه ششم فروردین 1386 بم با استفاده از داده‌های شتاب‌نگاری میدان نزدیک و تعیین سازوکارکانونی زمین‌لرزه}, abstract_fa = {به‌دنبال زمین‌لرزه مخرب 5 دی ماه 1382 بم در جنوب شرقی ایران، زمین‌لرزه‌ها و پس‌لرزه‌های زیادی در این منطقه به وقوع پیوسته است که آخرین آنها تا زمان این تحقیق، زمین‌لرزه 6 فروردین 1386 است. در این تحقیق به کمک داده‌های رقمی شتاب‌نگاری مرکز تحقیقات ساختمان و مسکن، نمودارهای شتاب، سرعت و جابه‌جایی نسبت به زمان، طیف فوریه شتاب با تصحیح خط مبنا و فیلتر مناسب ترسیم شد. پس از تصحیحات لازم پارامترهای منبع از جمله، بسامد گوشه و قسمت تخت نمودارتعیین شد. با استفاده از رابطه گشتاور لرزه‌ای، Mo و رابطه کاناموری مقدار بزرگی گشتاوری (Mw) ایستگاه‌های شتاب‌نگاری محاسبه شد. پس از میانگین‌گیری برای گشتاور لرزه‌ای مقدار 16 10 × 6ر3 Mo= نیوتن‌متر و برای بزرگای گشتاوری مقدار 5=Mw به‌دست آمد، که با نتایج محاسبات سایر روش‌ها مطابقت و همخوانی دارد. همچنین با استفاده از داده‌های رقمی شتاب‌نگاری رومرکز این زمین‌لرزه تعیین محل شد. انگیزه ما از این تحقیق، اثبات وجود تفاوت حداکثر 36 کیلومتری در تعیین محل زمین‌لرزه از سوی مراکز لرزه‌نگاری معتبر داخل و خارج است. در این تحقیق سازوکار کانونی این زمین‌لرزه به کمک اولین رسید موج p ایستگاه‌های شتاب‌نگاری و لرزه‌نگاری داخل کشور امتدادلغز راست‌گرد، با روند ?W51 Nشیب 89 درجه غربی و زاویه ریک 138درجه به‌دست آمد. وقوع زمین‌لرزه شش فروردین 1386 و خردلرزه‌های پس از آن حاکی از فعال بودن منطقه ازحیث لرزه‌خیزی است.}, keywords_fa = {بزرگی گشتاوری,پارامترهای منبع,تعیین محل زمین‌لرزه,داده‌های شتاب‌نگاری,زمین‌لرزه بم,گشتاور لرزه‌ای}, url = {https://jesphys.ut.ac.ir/article_21553.html}, eprint = {https://jesphys.ut.ac.ir/article_21553_0678f6c4d7d53fa0cd9c182bfe2e0ceb.pdf} } @article { author = {Amini, Ahmad and Fazel Bidgoli, yaser and Rokn Abadi, Hasan}, title = {Linear filter for upward- down ward continuation by using wavelet transform and its application in magnetic data processing}, journal = {Journal of the Earth and Space Physics}, volume = {36}, number = {2}, pages = {-}, year = {2010}, publisher = {Institute of Geophysics, University of Tehran}, issn = {2538-371X}, eissn = {2538-3906}, doi = {}, abstract = {Linear filters are used for a wide range of magnetic science including noise attenuation, spatial derivatives, upward and downward continuation and reduction to the pole. The majority of these filters are spatially-invariant, meaning that the filter has a constant wave number response over the whole signal. In contrast, many signals which encountered real problems such as magnetic signals, typically exhibit a spatially-varying wave number content which motivates us to design filters with spatially-varying wave number responses. This leads to better preservation of anomaly gradients in the calculated derivatives than is possible using conventional Fourier or space domain smoothing techniques. Fourier transform because of its stability, the simple physical interpretation of the transform coefficients and the diagonalisation of spatially invariant linear operators in the Fourier domain play an important role in magnetic processing. However, Fourier filters cannot be designed to adapt to local properties of the signal or to generate spatially-varying filters. One of the methods for generating spatially-varying filters is based on the continuous wavelet transform (CWT) which provides new powerful tools in magnetic data processing. The wavelet transform is ideal for analysing signals such as magnetic signals that contain short duration, transient features. Wavelet techniques can be used to provide solutions to problems that are difficult or impossible to solve using conventional global techniques such as Fourier-based methods. The wavelet transform preserves both spatial and wave number information about a signal allowing us to design a range of spatially-varying filters that act on the wavelet coefficients. This method provides robust and efficient new frameworks for designing filters that are impractical to implement using conventional space or wave number domain techniques. This method is compared with other techniques in upward-downward continuation. We demonstrate the application of spatially-varying scale filters to the problem of upward and downward continuation from a level observation surface to a new irregular-height surface. Downward continuation is the most difficult of these operations as it is highly numerically unstable and is very sensitive to high wave number noise. For comparison, conventional methods of downward continuation, such as the Taylor-series and chessboard methods are used which for stabilizing a global low-pass filter are applied to the data to attenuate any high wave number noise that may create difficulties in the continuation procedure. The wavelet implementation produces a superior result compared with conventional techniques such as Taylor-series and chessboard algorithms. In this study, the wavelet approach combined with the exponential smoothing filter produces sharper images than either the chessboard or Taylor-series methods that are clearly evident in the case study and synthetic examples. In contrast, the Taylor-series cosine roll-off filter is designed to ensure that the downward continuation is stable over the largest continuation distance. The chessboard method should theoretically be able to behave like the wavelet method by adjusting the amount of smoothing applied to each downward continuation slice. However, the sliding-rule filter does not appear to be as effective as the wavelet exponential filter. The large differences in the performance of each of these downward continuation methods highlight the significance in the choice of smoothing method. In addition, the wavelet exponential filter has the advantage of being locally adapted to the signal, which means that we do not need to oversmooth the signal when the local downward continuation distances are small but the chessboard method suffers from oversmoothing which is needed to prevent artifacts from continuing downward below the shallowest magnetic sources. One advantage of the wavelet and chessboard methods over the Taylor-series method is simplicity to automatically generate the parameters needed to design the smoothing filters for the wavelet and chessboard methods. In contrast, the Taylor-series method requires some trial-and-error intervention by the user. With a careful choice of smoothing parameters, the Taylor-series method can be designed to perform equally as well as the wavelet method. However, the choice of these parameters is often difficult when the downward continuation distances are large.}, keywords = {Linear filter,Magnetic data,Upward-downward continuation,Wavelet transform,Wave number}, title_fa = {فیلتر خطی برای ادام? فراسو و فروسو به کمک تبدیل موجک و کاربرد آن در پردازش داده‌های مغناطیسی}, abstract_fa = {تبدیل‌های موجک عرص? جدیدی برای پردازش داده‌های مغناطیسی فراهم می‌آورند. این تبدیل‌ها، اطلاعات حاصل از وابستگی فضایی عدد موج‌ها را حفظ می‌کنند و به‌این‌صورت می‌توان به کمک ضرایب موجک، فیلترهای متغیر با فضا را طراحی کرد. در این بررسی، به کمک تبدیل موجک پیوسته، فیلتر خطی یک‌بُعدی و دوبُعدی متقارن شعاعی با پاسخ‌های عدد موجی وابسته به فضا، ساخته می‌شود. یکی از کاربردهای این فیلتر استفاده به‌منزلة عملگر ادام? فراسو و فروسو است. بسیاری از فیلترها برای استفاده در حوز? فضایی یا عدد موجی، غیرِ کاربردی هستند؛ فیلتر موجک، چارچوب قوی و مؤثر جدیدی را برای بهبود این‌گونه فیلترها فراهم می‌آورد. مشکلات موجود ناشی از وجود نوفه که به‌شدت ادام? فروسو را تحت تاثیر قرار می‌دهد، با این فیلتر کاهش می‌یابد. در انتها این فیلتر با دیگر روش‌های متداول برای ادام? فراسو و فروسو مقایسه می‌شود.}, keywords_fa = {ادام? فراسو و فروسو,تبدیل موجک پیوسته,داده‌های مغناطیسی,عدد موج,فیلتر خطی}, url = {https://jesphys.ut.ac.ir/article_21554.html}, eprint = {https://jesphys.ut.ac.ir/article_21554_dc99b2fbdd1c2b8cf01e2ab79c614674.pdf} } @article { author = {Arasteh, Hoda and Javaherian, A}, title = {Imaging thin beds using attributes achieved in spectral decomposition by short time Fourier transforms}, journal = {Journal of the Earth and Space Physics}, volume = {36}, number = {2}, pages = {-}, year = {2010}, publisher = {Institute of Geophysics, University of Tehran}, issn = {2538-371X}, eissn = {2538-3906}, doi = {}, abstract = {One of the most fundamental reservoir characteristics is the thickness. The analysis of thin bed tuning on seismic reflectivity has been studied extensively by Wides (1973) and Neidel and Pogiaggliomi (1977), who discussed the limits of seismic resolution. During the past decade, the industry has developed a plethora of new attributes in studying thin beds by employing spectral decomposition (Peyton et al., 1998; Partyka et al., 1999), and attributes which are obtained from it (Marfurt and Kirlin, 2001). Spectral decomposition refers to all methods that generate frequency spectrums consisting of amplitude spectrum, phase spectrum, change of phase with frequency and power spectrum in windows with the center of each time sample of a trace. These methods are used in studying geological features, thin beds, hydrocarbon reservoirs and noise attenuation. The most important of these methods are short time Fourier transform (STFT), continuous wavelet transforms, S-transform, Wigner-Ville distribution and matching pursuit decomposition. The result of a trace spectral decomposition is a time-frequency map. In this paper, the STFT method is applied for imaging thin beds in 2D seismic sections. In this method, a short time window with constant size is multiplied by a trace for each of its time samples of it then Fourier transform is applied on it. So, there is a local frequency spectrum for each time sample. The main factor of this method to be considered is the effect of the existence of two boundaries of a thin layer inside a time window. This affects local frequency spectrum and changes it so that by choosing some quantities of local spectrum as attributes it is possible to image the thin bed. The spectral attributes that are studied are the peak frequency, peak amplitude and change of phase with frequency. Detectable layers are related to the size of the time window. In this paper, frequency spectrums of a bed response and their reflection with minimum and zero phase wavelets are considered. Then, the effect of changing of thickness on frequency spectrums is studied. For controlling the main factor of this method, STFT algorithm is applied on a synthetic seismic section. It is shown that with the peak frequency and the peak amplitude obtained from a small portion of each local amplitude spectrum, and local change of phase with frequency in a particular frequency a good conclusion can be obtained. After that the STFT algorithm is applied on a part of a real seismic section and the place of layers with thickness up to 10 ms, between 10 to 18 ms and more than 18 ms are defined. It should be noted that the main idea of detecting thin beds on seismic sections is deference of frequency spectrum of a bed and a boundary. It is necessary to use more than one attribute for achieving more accurate results. It is shown that in particular frequency portions, thin beds are detected better. The minimum thickness that can be detected by this method depends on the frequency content of the seismic wavelet and the time sampling interval. It should be considered that applying STFT in more stages and less difference of size of windows in alternative stages increases accuracy of estimated thickness ranges.}, keywords = {Local change of phase with frequency,Peak amplitude,Peak frequency,Short time Fourier transform,Spectral decomposition}, title_fa = {تصویرسازی لایه‌های نازک با استفاده از نشانگرهای به‌دست آمده از تجزیه طیفی به روش تبدیل فوریه زمان کوتاه}, abstract_fa = {تجزیه طیفی به تمامی روش‌هایی گویند که برای هر پنجره کوچک به مرکز یک نمونه زمانی از ردلرزه، طیف بسامدی (طیف دامنه، طیف فاز یا تغییرات فاز با بسامد و طیف انرژی) را می‌دهد. بنابراین حاصل کار تجزیه طیفی یک ردلرزه صفحه نمایش زمان- بسامد است. در این مقاله تجزیه طیفی به روش تبدیل فوریه زمان کوتاه در تصویرسازی لایه‌های نازک در مقاطع لرزه‌ای بررسی شده است. شیوه کار بدین صورت است که از صفحه نمایش زمان- بسامد مربوط به هر ردلرزه، نشانگرهای طیفی استخراج می‌شوند تا لایه نازک را تصویر کنند. نشانگرهای مورد بررسی در این تحقیق شامل بسامد قله، دامنه قله و تغییرات محلی فاز با بسامد هستند. الگوریتم روش تبدیل فوریه زمان کوتاه با نرم‌افزار MATLAB طراحی شده است. برای ارزیابی این روش در بررسی لایه‌های نازک، الگوریتم طراحی شده ابتدا روی یک مقطع مصنوعی و سپس روی قسمتی از یک مقطع لرزه‌ای واقعی با فاز صفر اعمال شد. نشانگرهای بسامد قله و دامنه قله در بازه بسامدی کوچک و نشانگر تغییرات محلی فاز با بسامد به ازای یک بسامد خاص نتیجه قابل قبولی عرضه کردند.}, keywords_fa = {بسامد قله,تبدیل فوریه زمان کوتاه,تجزیه طیفی,تغییرات محلی فاز با بسامد,دامنه قله}, url = {https://jesphys.ut.ac.ir/article_21555.html}, eprint = {https://jesphys.ut.ac.ir/article_21555_507e9008b7d7a42dedb3d3f5ebdbf168.pdf} } @article { author = {Sardar, hadi and Memarian, Hosein and Mahbaz, S Bijan}, title = {Porosity Estimation from Seismic Data in the Asmari Reservoir of Mansuri Oil Field: A Rock Physics Modeling Approach}, journal = {Journal of the Earth and Space Physics}, volume = {36}, number = {2}, pages = {-}, year = {2010}, publisher = {Institute of Geophysics, University of Tehran}, issn = {2538-371X}, eissn = {2538-3906}, doi = {}, abstract = {The ultimate goal of the reservoir geophysicist is to determine the reservoir properties (lithology and porosity) and conditions (pore fluid and pressure) from seismic data. This goal can be achieved by applying rock physics transforms to a volume of seismically-derived elastic properties. The basis for establishing rock physics transforms are controlled experiments where the reservoir and elastic sediment properties are measured on the same sample of reservoir under the same conditions. Such experimental data come from well logs and cores. The purpose of this paper is to apply rock physics diagnostic to well log data from the Asmari reservoir in Mansuri oil field and then used to quantify porosity from acoustic impedance inversion volumes. The results show that Nur-Dvorkin (1996) model and Raymer-Greenberg-Castagna (1997) are good for predicting porosity in sandstone and the carbonate zone of the Asmari reservoir respectively.}, keywords = {Acoustic impedance,Asmari Reservoir,Gassmann equation,Nur-Dvorkin model,Raymer-Greenberg-Castagna Model,Rock Physics}, title_fa = {برآورد تخلخل از داده‌های لرزه‌ای با استفاده از مدل‌سازی فیزیک سنگ در مخزن آسماری میدان منصوری}, abstract_fa = {هدف نهایی متخصصان ژئوفیزیک مخزن، تعیین خصوصیات مخزنی، نظیر سنگ‌شناسی و تخلخل و شرایط آن، نظیر فشار و نحوه توزیع سیال با استفاده از داده‌های لرزه‌ای است. برای دستیابی به این هدف می‌توان مدل‌های فیزیک سنگ را بر حجم خصوصیات کشسانی به‌دست آمده از داده‌های لرزه‌ای اعمال کرد. اساس مدل‌های فیزیک سنگ، تعیین خصوصیات کشسانی و مخزنی در شرایط یکسان به لحاظ زمین‌شناسی و با استفاده از آزمایش‌های کنترل شده است. داده‌های آزمایشگاهی مورد نظر غالبا از نگار و یا مغزه به‌دست می‌آیند. هدف از این مقاله، استفاده از تحقیقات فیزیک سنگ و نگار به‌منظور محاسبه کمّی تخلخل از حجم امپدانس صوتی حاصل از معکوس‌سازی داده‌های لرزه‌ای در مخزن آسماری میدان منصوری است. نتایج این بررسی نشان می‌دهد که مدل نیور- ورکین (1996) برای محاسبه تخلخل در بخش ماسه‌سنگی و مدل رایمر-گرینبرگ-کاستاگنا (1997) برای پیش‌بینی تخلخل در بخش کربناته مخزن آسماری، مناسب است و این دو مدل، مقدار تخلخل را در بخش مخزنی سازند آسماری تا %30 برآورد می‌کنند. این مقدار تطابق خوبی با مقادیر اعلام شده در گزارش‌های فیزیک ‌سنگ مخزن، نشان می‌دهد.}, keywords_fa = {امپدانس صوتی,فیزیک سنگ,مخزن آسماری,مدل رایمر-گرینبرگ-کاستاگنا,مدل نیور- ورکین,معادله گاسمان}, url = {https://jesphys.ut.ac.ir/article_21556.html}, eprint = {https://jesphys.ut.ac.ir/article_21556_5c95bdde62e41615c1a01191db0deac0.pdf} } @article { author = {Joata Bayrami, A and Hosseinzadeh Guya, Nser and mirzaei, mahmood}, title = {Determination of the electrical conductivity and temperature of the upper mantle using Sq field}, journal = {Journal of the Earth and Space Physics}, volume = {36}, number = {2}, pages = {-}, year = {2010}, publisher = {Institute of Geophysics, University of Tehran}, issn = {2538-371X}, eissn = {2538-3906}, doi = {}, abstract = {A fluctuating electric current flowing in the Earth’s atmosphere causes corresponding electric currents to flow in the conducting Earth below the source current. The depth of penetration of the induced currents is determined by the characteristics of the source currents as well as the distribution of electrically conducting materials in the Earth. At the Earth’s surface observatories measure the composite of external (source) and internal (induced) field components from the currents. The quiet daily field variations, called Sq for 'solar quiet-time', provide a natural signal source with frequencies appropriate to upper mantle conductivity studies. This paper is concerned with the quiet-day field variations, their separation into external and internal contributions, and the use of this separation to profile the electrical conductivity of the Earth's upper mantle. For the situation in which field measurements are available about a spherical surface that separates the source from the induced currents (and a current doesn’t flow across this surface), Gauss (1838) devised a special solution of the differential electromagnetic field equations that is separable in the spherical coordinates r,? and ? In Gauss’s solution, the field terms that represent radial dependence appear as two series-one with increasing powers of the sphere radius, r, and one with increasing powers of 1/r. As the value of r increases (outward from the sphere) the first series produces increased field strength, as if approaching external current sources. As the value of r decreases (toward the sphere center) the second series of 1/r terms indicate increased field strength, as if approaching internal current sources. Gauss had devised the way to separately represent the currents that were external and internal to his analysis sphere. For these external and internal series, there are individual spherical harmonic analysis (SHA) polynomial (Legendre) terms, each having two indices, degree m and order n. The full field is then represented as paired (external and internal Legendre terms) elemental parts, each satisfying the physical laws. Gauss applied the SHA method to the global field observations and verified that most of the Earth’s main field originated from internal sources. By the turn of the century, the method was also used to show that the daily, quiet-time geomagnetic field variations came mostly from sources of current external to the Earth (Schuster, 1889 and 1908); this finding led to the discovery of the ionosphere. We will analyze the quiet fields with this SHA to separate the Sq ionospheric source currents from the induced currents within the Earth. On days undisturbed by solar-terrestrial field and particle activity the geomagnetic records from a surface observatory display a smooth variation of field during the daylight hours. These variations are dominated by 24-, 12-, 8-, and 6-hr spectral components that change slowly from day to day through the seasons. Some authors prefer to select quiet days by a limiting value of the day's Ap (e.g., Ap = 10). Others take a fixed number (e.g., five) of the quietest days for a given month, whatever the values may be. We have preferred to select those days for which no Kp index exceeded a limiting level (e.g., 2+) because active records can cause unrealistic conductivity computations. The 19 North American observatories are selected in this study. There are 60 days in 1997 in which the global geomagnetic disturbance index, Kp, have all 8 daily values less than 2+. These days are taken as preliminary "quiet day" recordings. All observatories have 60-min sample records. The original recordings of field are in Universal Time (UT) as orthogonal north, east, and into-the-earth components of field as X, Y, and Z. The data for each component are Fourier analyzed for each quiet day. In the method outlined by Schmucker (1970) for profiling the Earth's substructure, formulas are developed that provide the depth (d) and conductivity ( ) of apparent layers (substitute conductors) that would produce surface-field relationships similar to the observed components. These profile values, (1) and (2) need to be determined for each n, m set of SHA coefficients using the real and imaginary parts of a complex induction transfer function, , given as: (3) Schmucker (1970) showed that the transfer function is obtained from the ratio of the field components, for a given n, m, as: (4) (5) where X, Y and Z are the northward, eastward, and into the Earth field components in gammas, R is the Earth's radius in kilometers, ? is the colatitude at the field measurement location, and is the Schmidt normalized associated Legendre polynomial. To use the above conductivity-depth formulation, it is now necessary to find the individual n, m SHA terms of the three field components at the field-measurement location. With the fitted potential function (V) given as: (6) in which the cosine (A) and sine (B) coefficients of the expansion for the external (ex) and internal (in) parts are taken to be: (7) Campbell and Anderssen (1983) showed that Equations (1) and (2) may be determined directly from the SHA coefficients with the expressions (8) (9) Electrical conductivity properties of the upper mantle for a North American sector of the Earth have been determined using the 24-, 12-, 8-, and 6-hr spectral components of the quiet-day geomagnetic variations. Spherical harmonic coefficients obtained from an analysis of the three components of the quiet daily variation (Sq) field for the solar-quiet year of 1997 were applied to a modeling procedure that was modified from Schmucher's (1970) publication. From a depth of about 100 to 650 km, the conductivity, , may be represented by with d is the depth in kilometers. Small perturbations of conductivity indicating some layering at 100 to 400 and 400 to 650 km correspond to the similar behavior of the Earth's density in these regions. From temperature-depth models we infer that the multiphase bulk properties of the expected silicates in these regions behave approximately as with T as the temperature in Kelvin.}, keywords = {Electrical conductivity,Electromagnetic response,Skin depth,Sq field,Substitute conductor}, title_fa = {تعیین رسانش الکتریکی و درجه حرارت گوشته‌ فوقانی با استفاده از تغییرات میدان}, abstract_fa = {رسانش الکتریکی گوشته فوقانی در قطاع امریکای شمالی با استفاده از مؤلفه‌های طیفی 24، 12، 8 و 6 ساعته تغییرات میدان به‌دست آمده است. ضرایب هماهنگ کروی حاصل از تحلیل سه مؤلفه تغییرات آرام روزانه میدان برای سال آرام خورشیدی 1997 در یک روش مدل‌سازی به‌کار برده شده‌اند که اصلاح شده روش اشموکر (1970) است. از عمق 100 تا 650 کیلومتر، هدایت الکتریکی ? برحسب زیمنس بر متر را می‌توان با رابطه نشان داد که در آن d، عمق برحسب کیلومتر است. از مدل‌های دما عمق می‌توان نتیجه گرفت که خواص سیلیکات‌های موجود در این نواحی منجر به رابطه تقریبی می‌شوند که T دما در مقیاس کلوین است.}, keywords_fa = {پاسخ الکترومغناطیسی,رسانای هم‌ارز,رسانش الکتریکی,عمق پوستی,میدان sq}, url = {https://jesphys.ut.ac.ir/article_21557.html}, eprint = {https://jesphys.ut.ac.ir/article_21557_91877c3b6d618484cfa5a50751741d01.pdf} } @article { author = {Mashayekhi, Robab and Irannejad, Parviz and Aliakbari-Bidokhti, A}, title = {The simulation of aerosols and its radiative forcing using the new coupled system of aerosol HAM model with the Weather Research and Forecasting (WRF) model}, journal = {Journal of the Earth and Space Physics}, volume = {36}, number = {2}, pages = {-}, year = {2010}, publisher = {Institute of Geophysics, University of Tehran}, issn = {2538-371X}, eissn = {2538-3906}, doi = {}, abstract = {A new coupled system of the aerosol HAM model and the Weather Research and Forecasting (WRF) model is introduced and tested in this paper. Unlike the two other aerosol schemes currently coupled in the WRF model (MADE/Sorgam and MOSAIC), the HAM applies a new "pseudomodal" approach for the representation of the aerosol particles. The aerosol population in this model is represented by the superposition of seven modes based on size and solubility. The particle size is assumed to be distributed log-normally. The seven modes are categorized into four geometrical size classes, ranging from the nucleation, Aitken and accumulation modes to coarse modes. The aerosols are also divided into two types of internally mixed and water soluble particles (four modes), and externally mixed and insoluble particles (three modes). This classification makes possible the prediction of the hygroscopic properties of initially insoluble aerosol compounds which controls their atmospheric lifetimes and also their interactions with clouds. The WRF-HAM model includes the various microphysical processes of condensation of sulfuric acid, nucleation and new particle formation, coagulation of aerosol particles, and the thermodynamical equilibrium of aerosols with the water vapor. The main removal processes for the aerosol particles in the coupled WRF-HAM model are gravitational sedimentation and dry deposition. The model also considers the in-cloud scavenging of aerosol particles by precipitation within the convective cumulus clouds. The main global aerosol compounds including sulfate, black carbon, particulate organic carbon (POM), sea salt and mineral dust have been considered in this study. The emission fluxes of different aerosol compounds are based on the prescribed Emission Inventory for the Aerosol Model Inter-comparison Experiment B, AEROCOM representative for the year 2000. The simulations are carried out for a 6-day simulation period from 6 to 12 May 2006 over a domain with 30-km grid spacing, covering south-western Asia, North Africa and some parts of Europe. The diurnal variation of the simulated hourly PM10 mass concentration at Tehran is qualitatively close to the hourly observations made by the Air Quality Control Company (AQCC) of the Municipality of Tehran. The model captures diurnal cycle and the magnitude of the observed PM10 concentration during most of the simulation period. Coupling WRF model with HAM aerosol scheme improves the hourly PM10 mass concentration compared to the case where WRF was coupled with MADE scheme. A negative radiative forcing and cooling of the atmosphere are found mainly over the regions of high emission of mineral dust. The absorption of shortwave radiation by black carbon causes warming effects in some regions with positive radiative forcing. The inclusion of aerosol feedback in the shortwave radiation scheme improves the simulated daily mean shortwave radiation fluxes in Tehran. The difference between the simulated and observed mean daily downward shortwave radiation, temperature and surface pressure by the HAM model is smaller than by the aerosol MADE scheme. Compared to WRF-MADE, using the coupled WRF-HAM model improves the simulation of downward shortwave radiation by up to 40 Wm . The spatial variation in the simulated mean optical depth of aerosols at 500 nm wavelength by the WRF-HAM model is qualitatively close to the measurements by MODIS instrument. Furthermore, the comparison of the simulated aerosol optical depth by HAM and MADE aerosol schemes and the observations in Solar Village site of the global AERONET Network shows that the HAM model highly outperforms MADE.}, keywords = {Aerosol HAM model,aerosols,Radiative forcing,WRF model}, title_fa = {شبیه‌سازی هواویزها و واداشت تابشی ناشی از آنها با استفاده از مدل جفت‌شده هواویز HAM و مدل میان‌مقیاس پیش‌بینی وضع هوا WRF}, abstract_fa = {در این پژوهش سیستم جفت‌شده جدیدی از طرحواره هواویز HAM با مدل میان‌مقیاس پیش‌بینی عددی وضع هوا WRF معرفی شده است. در مقایسه با طرحواره‌های هواویز به‌کاررفته در مدل WRF، طرحواره هواویز HAM از دستاورد جدیدی تحت عنوان روش شبه مدال برای توزیع اندازه ذرات هواویز استفاده می‌کند. پنج نوع اصلی هواویزهای جهانی شامل سولفات، کربن سیاه، کربن آلی، ذرات گردوغبار و نمک دریا در این مدل در نظر گرفته شده‌اند. شبیه‌سازی‌های اولیه حاصل از مدل جفت‌شده WRF-HAM برای یک دوره شبیه‌سازی شش روزه از 6 مه تا 12 مه 2006 آورده شده است و نتایج با شبیه‌سازی‌های حاصل از طرحواره هواویز موجود در نسخه دوم مدل WRF، یعنی طرحواره MADE، مقایسه شده‌اند. غلظت جرمی شبیه‌سازی شده ذرات PM10 در منطقه تهران با استفاده از مدل جفت شده جدید، بهبود قابل ملاحظه‌ای نسبت به طرحواره هواویز MADE نشان می‌دهد. وارد کردن هواویزها در شبیه‌سازی‌های مدل، منجر به واداشت تابشی منفی و ایجاد سرمایش به‌ویژه در مناطق با مقادیر بزرگ شار گسیل ذرات گرد و غبار و واداشت تابشی مثبت و تولید گرمایش در مناطق دارای هواویزهای کربن سیاه در حوزه شبیه‌سازی می‌شود. اختلاف شار تابش طول موج کوتاه پایین‌سوی شبیه‌سازی شده و مقادیر مشاهداتی در منطقه تهران با وجود هواویزها کوچک‌تر شده و این بهبود در طرحواره جدید HAM نسبت به طرحواره MADE قابل ملاحظه‌تر است. شار تابش موج کوتاه پایین‌سو تا بیش از 40 وات بر مترمربع در شبیه‌سازی‌های HAM بهبود می‌یابد. تغییرات کوچکی نیز در دمای جو و فشار سطحی با در نظر گرفتن هواویزها در مدل ایجاد می‌شوند. متوسط روزانه عمق نوری شبیه‌سازی شده با مدل، توافق خوبی با توزیع مکانی مشاهداتی از داده‌های ماهواره MODIS نشان می‌دهد. از سوی دیگر مقادیر شبیه‌سازی شده عمق نوری در طول موج 500 نانومتر با طرحواره HAM با مقادیر مشاهداتی حاصل از داده‌های جهانی AERONET در ایستگاه سولار ویلیج نیز بهبود قابل ملاحظه‌ای نسبت به عملکرد MADE نشان می‌دهد.}, keywords_fa = {طرحواره هواویز HAM,مدل میان‌مقیاس WRF و واداشت تابشی,هواویز}, url = {https://jesphys.ut.ac.ir/article_21558.html}, eprint = {https://jesphys.ut.ac.ir/article_21558_7ad182d3c2226a06540a0b3812bcb636.pdf} } @article { author = {Esmaeilpour, Sajjad and Riahi, Mohammad Ali and Emamzadeh, A and Sokooti, Reza}, title = {Seismic AVO Modeling Using In situ Samples for Reservoir Study}, journal = {Journal of the Earth and Space Physics}, volume = {36}, number = {2}, pages = {-}, year = {2010}, publisher = {Institute of Geophysics, University of Tehran}, issn = {2538-371X}, eissn = {2538-3906}, doi = {}, abstract = {Distinguishable amplitude phenomena on surface seismic data often resulted from contrasts in the elastic parameters of subsurface layers. Various techniques have been involved to analyze and highlight such phenomena for their potential use as "Direct Hydrocarbon Indicator (DHI)". More recently, other techniques have been developed based on the variation of reflection coefficient with angle of incidence, conventionally called Amplitude-Versus-Offset (AVO). During the last twenty years the significance of AVO analysis for studying seismic reflection in oil exploration has been considered more importantly. In this work, first, a seismic line from a gas field and also a well are selected to indicate the results of the application of AVO analysis for detection of hydrocarbon reservoir in this field. In this project, using well logs and information obtained by core analysis, a synthetic seismogram has been built applying Zoeppritz equation. And Using Hampson-Russell software, AVO attributes have been extracted from synthetic seismogram. Then anomalies of these attributes have been investigated and compared with the anomalies from AVO attributes which were extracted from real seismic data to characterize the reservoir. It has been seen that the extracted attributes of the synthetic seismogram confirm the anomalies from real seismic data. Finally according to the obtained result, observed anomalies can be interpreted as a Gas Cap for this reservoir. This study is useful to identify reservoir and nonreservoirs and the results of this study are considered as input for detailed reservoir studies. In particular, knowing the reservoir physical and saturating fluid properties is of great importance in making plans for developing the reservoir.}, keywords = {attribute,AVO,Compressional Wave Velocity,Offset,Pre-stack,Shear wave velocity,Synthetic seismogram}, title_fa = {مدلسازی لرزه‌ای AVO با استفاده از داده‌های مغزه برای مطالعه مخزن}, abstract_fa = {امروزه به‌کارگیری آنالیز تغییرات دامنه لرزه‌ای با دورافت (AVO) در تفسیر کیفی و به‌ویژه کمّی داده‌های لرزه‌ای از مراحل اولیه اکتشافات مخازن هیدروکربوری گرفته تا تولید آنها، جایگاه ویژه‌ای در صنایع بالادستی نفت را به خود اختصاص داده است. استخراج خواص مخزنی شامل: تخلخل، محتوای سیال، فشار، جداسازی مرز سیالات هیدروکربوری در مخزن، تعیین رخساره‌های سنگی (به‌ویژه شیلی) و تعیین‌جهت و چگالی شکستگی‌های مخزنی و همچنین تعیین خواص کشسانی سنگ مخزن شامل مدول یانگ، مدول بالک، ضریب پواسون و مانند آن از جمله کاربردهای استفاده از آنالیز تغییرات دامنه لرزه‌ای با دورافت است. آنالیز تغییرات دامنه لرزه‌ای با دورافت با به‌کارگیری نشانگرهای لرزه‌ای ویژه، در شبیه‌سازی هرچه دقیق‌تر مخزن، روش‌های تحقیقاتی متداول لرزه‌ای را تحت‌الشعاع قرار داده است. با استفاده از این آنالیز، شناسایی مستقیم هیدروکربورهای گازی و تعیین سطح جدایش سیالات مخزنی امکان‌پذیر است و بدین‌ترتیب مخاطره اکتشافات نفت و گاز به میزان قابل توجهی کاهش می‌یابد. همچنین نتایج حاصل از آنالیزAVO می‌تواند به‌منظور اجرای تحقیقات پایشگری لرزه‌ای (چهاربُعدی) و اجرای هر چه بهینه‌تر روش‌های ازدیاد برداشت در مراحل بعدی توسعه‌ای و تولیدی مورد استفاده قرار گیرد. در این پروژه در مدل‌سازی مستقیم با استفاده از روابط زوپریتس (1919)، شوی (1985) و دیگران و به‌کارگیری نمودارهای صوتی (تراکمی و برشی) و چگالی، لرزه‌نگاشت مصنوعی با دورافت برای چاه‌های موجود ساخته می‌شود. از آنجاکه اندازه‌گیری نمودار صوتی از امواج برشی بسیار پرهزینه است، در مواردی که این نمودار در اختیار نباشد، با به‌کارگیری روابط فیزیک‌سنگی نظری و تجربی که پیچیدگی‌های سنگ مخزن در آنها لحاظ شده باشد، مانند معادله بیوت- گسمن و معادله کاستر- توکسوز، برآورد می‌شود و برای ساخت لرزه‌نگاشت مصنوعی با دورافت، مورد استفاده قرار می‌گیرد. در تحقیق حاضر، این داده‌ها با کمک نتایج حاصل از آنالیز مغزه در آزمایشگاه محاسبه شده است. در ادامه با کمک نگارهای چاهی با روش معادله زوپریتس مدل لرزه‌ای قبل از استک ساخته می‌شود. سپس به استخراج نشانگرها روی مدل مصنوعی می‌پردازیم. در پایان مدل لرزه‌ای مصنوعی و بی‌هنجاری‌های شناخته شده را با بی‌هنجاری‌های AVO آشکار شده در نشانگرهای استخراج شده از داده‌های لرزه‌ای واقعی، مقایسه و بررسی می‌کنیم و از آنها در شناسایی مستقیم هیدروکربن در مخزن بهره می‌گیریم. همچنین پس از این تحقیق، برای سازند مورد بررسی و در مناطق مشابه، لاگ سرعت موج برشی با استفاده از لاگ سرعت موج فشارشی قابل برآورد خواهد بود.}, keywords_fa = {پیش از برانبارش,تغییرات دامنه لرزه‌ای با دورافت (AVO),دورافت,سرعت امواج لرزه‌ای برشی و طولی,لرزه نگاشت مصنوعی,نشانگرهای لرزه‌ای}, url = {https://jesphys.ut.ac.ir/article_21559.html}, eprint = {https://jesphys.ut.ac.ir/article_21559_c7e6b8bf451a7b31eb12fd2aa508577e.pdf} } @article { author = {Ardestani, Vahid}, title = {Precise Edge detection of gravity anomalies by Tilt angle filters}, journal = {Journal of the Earth and Space Physics}, volume = {36}, number = {2}, pages = {-}, year = {2010}, publisher = {Institute of Geophysics, University of Tehran}, issn = {2538-371X}, eissn = {2538-3906}, doi = {}, abstract = {Measurements of the tilt angle of potential fields can be a useful aid to their interpretation. There are several varieties of the methods in use, such as tilt angle, and the derivative of tilt angle. Tilt angle filters are used for exact edge detection of the rectangular bodies artificially by writing a MATLAB code. The gravity effects of these models are computed by standard algorithms. The edges of these complex rectangular bodies are determined appropriately by tilt angle. It is also applied on real gravity data. The residual gravity anomalies are used as the real data. Tilt angles are used to determine the edges of these residual anomalies. The results are compared with the results from the analytic signal method which is a classical edge detection method. This shows the better performance of the tilt angle filter method.}, keywords = {Edge detection,gravity anomalies,Tilt angle filters}, title_fa = {تعیین دقیق گوشه‌های آنومالی گرانی با استفاده از روش فیلتر زاویه تیلت}, abstract_fa = {اندازه‌گیری‌های زاویه تیلت میدان پتانسیل می‌تواند کمک بسیار مهمی در تفسیر داده‌های گرانی باشد. علاوه بر زاویه تیلت گرادیان زاویه تیلت را هم می‌توان مورد استفاده قرار داد. از این روش برای تعیین دقیق گوشه‌های بی‌هنجاری استفاده می‌شود. از این روش برای تعیین گوشه‌های اجسام مکعبی شکل و با استفاده از برنامه نوشته شده در محیط مطلب (Matlab) استفاده شده است. از این روش برای تعیین گوشه‌های بی‌هنجاری داده‌های واقعی استفاده شده است. در این ارتباط از داده‌ای بی‌هنجاری باقی گرانی به‌منزلة داده‌های واقعی استفاده شده است. نتایج استفاده از این روش با نتایج حاصل از روش سیگنال تحلیلی مقایسه می‌شود. نتایج نشان می‌دهد که روش زاویه تیلت در آشکارسازی گوشه‌های بی‌هنجاری بهتر عمل می‌کند.}, keywords_fa = {بی‌هنجاری گرانی,تعیین گوشه,فیلترهای زاویه تیلت}, url = {https://jesphys.ut.ac.ir/article_21560.html}, eprint = {https://jesphys.ut.ac.ir/article_21560_c7b642e418d95278b42bca1c0653e8b2.pdf} } @article { author = {Montahaie, Mansoureh and Brasse, H and Oskooi, Behroz}, title = {Crustal conductivity structure of a continental margin, from magnetotelluric investigations}, journal = {Journal of the Earth and Space Physics}, volume = {36}, number = {2}, pages = {-}, year = {2010}, publisher = {Institute of Geophysics, University of Tehran}, issn = {2538-371X}, eissn = {2538-3906}, doi = {}, abstract = {The deep internal structure of the crust can be determined using appropriate seismic and electromagnetic methods. The natural source magnetotelluric (MT) method is the most suitable electromagnetic technique for probing into the deep crust. A long period magnetotelluric data set obtained in the Southern Chilean Andes is investigated in this paper. Dimensionality analysis shows that the data may be regarded as with a strike 2D direction aligned to the N-S direction. Results of a joint inversion of different MT data types indicate relatively high conductive structures in the middle to deep crust beneath the volcanic arc.}, keywords = {dimensionality analysis,Impedance,Inversion modeling,magnetotelluric}, title_fa = {تعیین ساختار هدایت ویژه پوسته زمین در یک حاشیه قاره‌ای با استفاده از داده‌های مگنتوتلوریک}, abstract_fa = {ساختار عمیق داخلی پوسته زمین با استفاده از روش‌های مناسب لرزه‌ای و الکترومغناطیسی مورد بررسی قرار می‌گیرد. ازمیان روش‌های الکترومغناطیسی روش مگنتوتلوریک که از منبع طبیعی (نوسانات میدان مغناطیسی زمین) استفاده می‌کند مناسب‌ترین روش برای کاوش پوسته پایینی زمین است. عمق کاوش در این روش بسته به بسامد نوسان‌های مورد استفاده تغییر می‌کند و هر چه بسامد کمتر باشد، ساختارهای عمیق‌تر بررسی می‌شوند. در این مقاله داده‌های مگنتوتلوریک بلند دوره (بسامد کوتاه) که در منطقه حاشیه قاره‌ای شیلی در امتداد دو نیم‌رخ شرقی-غربی و در عرض‌های جغرافیائی‌ 39.3 و 38.9 درجه جنوبی برداشت شده‌اند، مورد بررسی قرار گرفته ‌است. قطعه‌های زمین‌شناسی مهم که این دو نیم‌رخ‌ از آنها می‌گذرند عبارت‌اند از: رشته کوه‌های موازی ساحلی (Coastal Cordillera) ، دره طولی(Longitudinal Valley) و کمان آتشفشانی که دقیقاً زیر آتشفشان‌های فعال ویلاریکا (Villarrica) و لایما (Llaima) واقع شده ‌است (شکل1). مجموعاً 32 ایستگاه در امتداد این نیم‌رخ‌ها با فاصله تقریبی 10 کیلومتر از هم واقع شده‌اند و داده‌های مگنتوتلوریک را در محدوده دوره‌ای 20.000-10 ثانیه ثبت کرده‌اند. مهم‌ترین مشکلی که در استفاده از این داده‌های MT دورة بلند وجود دارد، اثر واپیچش گالوانیکی ناشی از ناهمگنی‌های سه‌بُعدی کوچک محلی است که در قسمت‌های سطحی‌تر قرارگرفته و باعث می‌شوند که پاسخ‌های ناشی از ساختار عمیق دوبُعدی منطقه‌ای پنهان شود. میزان انحراف داده‌های مگنتوتلوریک از مدل دوبُعدی منطقه‌ای با پارامتر پیچش (skew) اندازه‌گیری می‌شود. این مرحله که از آن تحت عنوان "تحلیل ابعادی" داده‌ها (dimensionality analysis) یاد می‌شود به روش بار (بار، 1988) صورت گرفته ‌است. مقدار آستانه پارامتر پیچش در این روش برای تمایز بین ساختارهای هدایت ‌ویژه دوبُعدی و سه‌بُعدی 3ر0 است. شکل3 نمایش پربندی مقادیر پیچش محاسبه شده برای این دو نیم‌رخ است. این شکل نشان می‌دهد که مقادیر پیچش این داده‌ها به‌طورکلی بین 3ر0 و 1ر0 توزیع شده‌اند. بنابراین رویکرد دوبّعدی در مورد این مجموعه داده‌ها موجه است. به‌منظور تفسیر داده‌های مگنتوتلوریک با استفاده از مدل‌های دوبُعدی ضروری است که یک جهت استرایک ثابت برای کل نیم‌رخ تعیین شود. تفسیر دوبُعدی داده‌ها فقط پس از چرخش آنها در جهت این استرایک امکان‌پذیر خواهد بود. در مورد این مجموعه از داده‌ها استرایک ساختار هدایت ‌ویژه منطقه‌ای با استفاده از روش "تانسور فاز" برای هردو نیم‌رخ در جهت شمالی- جنوبی محاسبه شده ‌است (شکل4). مدل‌سازی معکوس با استفاده از الگوریتم "گرادیان‌های مزدوج غیرخطی" و با استفاده از کد رودی و مکی (Rodi and Mackie, 2001) صورت گرفته‌ است. فضای مدل با 92 سطر و 148 ستون مجزا شده و مدل اولیه از یک نیم‌فضای همگن با مقاومت‌ویژه 100 اهم‌متر و اقیانوس با مقاومت‌ویژه 33ر0 اهم‌متر که محدوده آن با توجه به داده‌های عمق‌سنجی مشخص و تثبیت شده، تشکیل شده ‌است. نتایج حاصل از مدل‌سازی معکوس داده‌های مربوط به این دو نیم‌رخ (شکل‌های 5 و6) بسیار شبیه به هم هستند و ساختارهای رسانای متمایزی را در پوسته میانی و پایینی نشان می‌دهند. این توده‌های رسانا به منزله سیالات موجود در نواحی گسلی منطقه و یا نواحی با گداختگی جزئی (واقع در کمان آتشفشانی) تفسیر می‌شوند.}, keywords_fa = {امپدانس,تحلیل بُعدپذیری,مدل‌سازی وارون,مگنتوتلوریک}, url = {https://jesphys.ut.ac.ir/article_21561.html}, eprint = {https://jesphys.ut.ac.ir/article_21561_2610d30afbc579d374b23164285a26a2.pdf} } @article { author = {Taghavi, Farahnaz}, title = {Linkage between Climate Change and Extreme Events in Iran}, journal = {Journal of the Earth and Space Physics}, volume = {36}, number = {2}, pages = {-}, year = {2010}, publisher = {Institute of Geophysics, University of Tehran}, issn = {2538-371X}, eissn = {2538-3906}, doi = {}, abstract = {In recent years, interest in climatic extreme events (droughts, floods and heat waves) has been growing not only due to their direct impact on many socio-economic sectors, but also as climate change indicators. In this study, the behavior of extreme events is investigated using Extreme Climate Index Software (ECIS) package. The trend analysis of extreme indices is applied to daily series of temperature and precipitation observations from 16 synoptic stations over Iran in the period 1951-2003. Results show that the number of very warm days (T40) has increased while the number of very cool days (ID) has decreased and so the return period (?ret) of cold extreme indices increased and the Tret of warm extreme indices decreased. Changes in total and extreme precipitation indices vary, depending upon geographic location. Symmetric warming in the tails of most indices is seen overall. On the basis of a systematic analysis of observed changes in indices of climate extremes, a significant proportion of Iran is increasingly affected by a significant change in extreme indices in the ~50-yr study period.}, keywords = {climate change,Extreme Events,Iran,Linkage,Return period}, title_fa = {ارتباط بین تغییر اقلیم و رویدادهای حدی}, abstract_fa = {در سال‌های اخیر توجه به رویدادهای حدی (خشکسالی‌ها، سیل‌ها و امواج شدید گرمایی)، نه‌‌تنها به علت پیامدهای اقتصادی، اجتماعی بلکه به‌‌منزلة نشانگرهای تغییر اقلیم، اهمیت زیادی پیدا کرده است. در این تحقیق برای یافتن ارتباط بین تغییر اقلیم و رویدادهای حدی، رفتار و فراوانی رویدادهای حدی با استفاده از شاخص‌هایی در بسته نرم‌افزاری بنام شاخص‌های اقلیمی حدی Extreme Climate Index Software, ECIS)) بررسی شده است. دربسته نرم‌افزاری ECIS روند شاخص‌ها براساس داده‌های روزانه دما و بارش 16 ایستگاه همدیدی (شاخص نواحی اقلیمی ایران) در دوره آماری2003-1951 محاسبه شده است. در نهایت نقشه‌های روند شاخص‌های حدی در ایران تهیه شد. نتایج تحقیق نشان می‌دهد که روند کاهشی شاخص‌های حدی سرد FD وروند افزایشی در شاخص‌های حدی گرم مانند T40 در همه نواحی اقلیمی ایران به‌جز در ناحیه شمال غرب وجود دارد و گرمایش متقارن در بیشتر دنباله سری‌های زمانی شاخص‌های حدی به‌جز در دو ایستگاه همدان و ارومیه وجود دارد .همچنین نتایج مبین افزایش قابل توجهی در تعداد روزها و شب‌های خیلی گرم و کاهش در دامنه دمایی حدهای سالانه است و همین امر سبب کوتاه شدن دوره بازگشت شاخص‌های حدی گرم مانندSU ، T40و بلندترشدن دوره بازگشت شاخص‌های حدی سرد FD وID شده است اما روند شاخص‌های حدی بارش تغییرات اندکی را نشان می‌دهد. این نتایج فراوانی رویدادهای حدی اقلیمی مانند خشک‌سالی و طولانی شدن امواج گرم در اکثر نواحی اقلیمی ایران را تایید می‌کند. به‌طورِکلی براساس نتایج تغییر در روند شاخص‌های حدی در ایران بارز است و نشانه‌های روشنی از تغییر اقلیم وجود دارد.}, keywords_fa = {ارتباط,ایران,تغییر اقلیم,دوره بازگشت,روند,رویدادهای ‌حدی}, url = {https://jesphys.ut.ac.ir/article_21562.html}, eprint = {https://jesphys.ut.ac.ir/article_21562_b9ac73b050fdcd8e4c4343b4474a622a.pdf} }