@article { author = {فتاحی, مرتضی and حیدری, مریم}, title = {A review of luminescence dating of ancient walls and the study of the Tammisheh Wall}, journal = {Journal of the Earth and Space Physics}, volume = {42}, number = {3}, pages = {469-485}, year = {2016}, publisher = {Institute of Geophysics, University of Tehran}, issn = {2538-371X}, eissn = {2538-3906}, doi = {10.22059/jesphys.2016.58935}, abstract = {Thermoluminescence (TL) dating was established in Oxford to fulfil the demand for dating ceramics and other heated man made products (Aitken, 1985).  It determines the time elapsed since the mineral grains have been heated. Following the establishment of TL dating method, it was used to date fired materials such as pottery, flint and etc. Optically stimulated luminescence (OSL) determines the time elapsed since mineral grains have been last exposed to light just before burial (Aitken, 1989). Ancient walls are mainly built for defensive purposes.  Some walls were destroyed and erected again and some had been extended.  There is no accurate historical information about the time and purpose of the construction of the most of the ancient walls. Dating methods are key important tools in archaeological investigations.  Absolute date for ancient building including walls, facilitate accurate time of construction and interpretation about the reasons that the monument was built. The proud history of association of TL and OSL dating with archaeology is well explained in literature (e.g., Roberts, R.G. 1997; Feather j. K. 2003); Monuments and settlements were also dated by TL and OSL including three Minoan palaces in Crete (Liritzis and Thomas, 1980); stoneware from the period of the early Thai kingdom of Sukhothai (Robertson and Prescott, 1987, 1988); the remains of the strategically-located Roman city of Carnuntum in Austria (Erlach and Vana, 1988); the Terracotta Army in Xi'an, China ( Lu et al., 1988); fifteen Mesolithic to medieval habitation sites in Denmark, Sweden and Finland (Mejdahl, 1989); the foundations of the archaic fortifications on Palatino Hill in Rome itself (Bacci et al., 1990); the Spanish exploration of western Canada in the 18th century AD (Williams et al., 1991); stone lines and cairns on the Arnhem Land plateau in northern Australia (Chaloupka, 1993); Mycenean wall in Greece (Liritzis, 1994); temples and citadels in Sri Lanka (Abeyratne, 1994); the temple of Apollo in Delphi (Liritzis et al., 1997); the archaeological sites of Tall Abu Fahd and Tall Qsubi in the Middle Euphrates Valley, in Syria (Sanchez et al., 2008); the city wall system of the oasis in NW Saudi Arabia (Klasen et al., 2011); the domed Stupas in Sri Lanka (Bailief et al., 2013). Following the explanation of basis of OSL dating, this article presents the application of OSL for dating Tammisheh wall in Iran, as a case study. The Tammisheh Wall in northern Iran was around 12 km long from the Elburz Mountains into the south- of the Caspian Sea and is lined by a fort in dry land in the south and another fort that is currently under water, in the north.  A canal, conducted water along most of the wall. Tammīsheh Wall was in a good state of repair in the mid-7th century, because its inhabitants resisted against the powerful army of Sa‘īd b. al-‘Āṣ, son-inlaw of Caliph ‘Uthmān and governor of Kūfa, when it approached the Tammīsheh land corridor from the east in AD 650/651. (Sauer et al., 2014). The town of Tammīsheh continued to be occupied until at least the 13th. (Bivar and Fehérvári 1966). While Ibn Isfandiyār had attributed, the Tammīsheh to Peroz, Ibn Rusta (the early 10th-century author) suggested that the Tammīsheh Wall was built by Khusro I (AD 531–579). The Proposed dates, founder or period for construction of the wall is shown in Table 1. Until 2005, there had been no absolute date for the Tammisheh Wall.  14C is one of the most accurate methods but requires organic material. However, the wall is made of bricks that perfectly suits OSL dating. There are quartz and feldspar in the bricks. Quartz and feldspar are 2 main luminescence dosimeters. As bricks are cooked in ovens, the luminescence clock of quartz and feldspar inside the bricks will reset, due to temperature more than 300°C. The sediments contain quartz and feldspar. When sediments are exposed to day light, their luminescence clock will reset. The luminescence clock will start again as soon as bricks cool down or the sediment that was exposed to light is buried.  Therefore, we employed both OSL and carbon fourteen methods. We collected one OSL sample sediment below foundations of Tammisheh Wall (Fig 4) which unfortunately overestimated the real age. We also sampled one OSL sample and a charcoal from a preserved section of the kiln in Trench F (Fig 8). Another OSL sample and three radio carbon samples were collected from the trench “O” of Bansaran fort (Figs 14-17). All samples were transferred and dated in Oxford. The dating result of samples is presented in Table 1. The result of OSL and 14C ages proved conclusively that Tammisheh wall and Bansaran fort had been built during the Sasanian period.  }, keywords = {: Dating,Luminescence,Ancient walls,Tammisheh wall}, title_fa = {مروری برسن یابی دیوارهای باستانی با استفاده از روش درخشایی و عرضه پژوهش های صورت پذیرفته در مورد دیوارتمیشه}, abstract_fa = {با توجه به اینکه مقاله ای علمی به زبان فارسی در مورد نقش لومینسانس در سن یابی دیوارهای باستانی و همچنین دیوار تمیشه و اهمیت سن یابی آن وجود ندارد، در این نوشتار پس از مروری بر پژوهش هایی که در دنیا در باره سن یابی بناهای تاریخی به روش لومینسانس انجام شده است به دیوار تمیشه پرداخته و توضیح داده می شود که دیوار تمیشه با طول تقریبی دیوار 12 کیلومتر دیواری شمالی جنوبی است که از گوشه ی جنوب شرقی دریای خزر در نزدیکی بندر گز شروع شده و به دامنه ی رشته کوه البرز ختم می شود. در کنار دیواردر خشکی قلعه بنساران وجود داشته که حدود 700-500 متر با دیوار فاصله دارد و درغرب و جنوب دیوار تمیشه در جایی که دیوار به پایه کوه می رسد، واقع شده است. انتهای دیگر دیوار در دریا به یک جزیره ی کوچک منتهی می شده است که در آن یک قلعه ی فرضی به نام پل خشتی یا خشت سر وجود داشته است. اطلاعات تاریخی در مورد زمان ساخت دیوار تمیشه توافق ندارند. در سال 1384 (2005 میلادی) چهار نمونه لومینسانس برای سن یابی برداشت شد. نتایج سن یابی با سن های بدست آمده توسط روش رادیو کربن همخوانی داشت و دوران ساخت این دیوار را در زمان ساسانیان نشان داد.}, keywords_fa = {سن یابی,لومینسانس,دیوار های باستانی,دیوار تمیشه}, url = {https://jesphys.ut.ac.ir/article_58935.html}, eprint = {https://jesphys.ut.ac.ir/article_58935_f7ad3b6da93a750cc601d6f2f63a2fd5.pdf} } @article { author = {khaje, faeze and taraz, majid and neamati, majid and bahrampour, Ali reza}, title = {Designing possibility of a seismometer using Fiber Bragg Grating and metal bellows}, journal = {Journal of the Earth and Space Physics}, volume = {42}, number = {3}, pages = {487-498}, year = {2016}, publisher = {Institute of Geophysics, University of Tehran}, issn = {2538-371X}, eissn = {2538-3906}, doi = {10.22059/jesphys.2016.58137}, abstract = {In this study, we have theoretically investigated designing possibility of a seismometer using Fiber Bragg Grating (FBG) and Metal Bellows pairs attached on a mechanical system. This new seismometer can record earthquakes according to sensitivity of fiber optic to changes in physical parameters such as stress and strain. Firstly, in order to understand the quantity and quality of sensor performance, with a mathematically brief description, the effect of stress and strain on the FBG reflection spectrum is examined. FBG is an intrinsic sensing element which can be photo-inscribed into a silica fiber. The basic principle of this operation commonly uses FBG based sensor system. The operation is monitoring of the shift in wavelength of the Bragg signal following changing in stress, strain and also temperature. The Bragg wavelength, or resonance condition of a Grating, is                                                                                                                                                                                                                   (1) The wavelength of the optical signal reflected by the Bragg Grating depends on the FBG physical parameters ( the fiber's effective refractive index, and  is the wave length of the FBG). These parameters are changed if the Grating is subjected to mechanical deformation or temperature variation. The measured strain response at constant temperature is                                                                                                                                              (2)                                                                                                                                                               (3) The structure of the seismometer consists of an inertial mass, supported by a L-shaped aluminum cantilever beam, connected to the structure base by metal bellows and an FBG element. In case of ground acceleration, the inertial mass moves in the vertical direction, imposing a compressed or stretch of the optical fiber. This deformation induces variation on the FBG Bragg wavelength. From the system work-energy concept analysis, the seismometer undamped natural frequency can be written as:                                                                                                                                                  (4) with natural frequency about 7.8 Hz. Ideally, we can consider a seismometer as a black box whose input is ground motion (represented by a kinematic variable: displacement, velocity or acceleration) and its output is displacement Bragg wavelength. Assuming a ground acceleration ,, by using the frequency range of 2011 (MW 6.0) Rigan earthquake waveform, strain (5) and displacement of amplitudes Bragg wavelength (6) and phase change (7) for different angular frequencies and amplitudes of the ground acceleration (for different damping coefficient 0.4, 0.6, 0.8, 1) are obtained. The relationships between them are investigated and for different angular frequencies, they are plotted and compared.                                                                                                                                               (5)                                                                                                                                       (6)                                                                                                                                                                     (7) As  (undamped), the solution has an increasing amplitude as , which is called resonance. If  (under damped), the signal tend to ring at near seismometer period. For  the signal is critically damped and oscillation is minimized. For  (overdamped), no oscillations occur, but the mass returns to rest more slowly. The possibility of designing a seismometer is investigated which works in near critical damping (0.6, optimal damping coefficient) which records a part of earthquake waves with frequency close to the natural frequency of the seismometer with the maximum amplitude and without phase-shift. Also, to prove performance and ability of seismometer with the mentioned method two other earthquakes, 2010 (MW 6.4) Rigan earthquake and 2011 (MW 7.8) Saravan earthquakes were investigated. Also each of the above events recorded and processed at the same time by several seismic stations and compared with seismometer Fiber Bragg Grating. A seismometer which uses Fiber Bragg Grating has a few number of advantages such as immunity to electromagnetic interference, their capability to transmit signal over long distance without any additional amplifiers, highly accurate and digital output, lighter weight, high sensitivity, low power consumption, lower cost and a wide range of dynamics.}, keywords = {Keywords: earthquake,Stress,strain,seismometer,Fiber Bragg Grating,Metal Bellows}, title_fa = {امکان سنجی طراحی دستگاه لرزه‌نگار با به‌کارگیری توری براگ فیبر نوری و فلز ارتجاعی}, abstract_fa = {هدف از این پژوهش، طراحی لرزه‌نگاری با به‌کارگیری توری براگ فیبر نوری و فلز ارتجاعی در یک سیستم مکانیکی می‌باشد. این لرزه‌نگار، با توجه به ویژگی و حساسیتی که حسگر فیبر نوری نسبت به پارامترهایی از قبیل اثرات تنش و کرنش دارد، زمین‌لرزه را ثبت می‌کند. برای ثبت امواج زمین‌لرزه باید فرکانس طبیعی دستگاه نزدیک فرکانس ارتعاش زمین باشد که فرکانس طبیعی دستگاه بهبود بخشیده و حدود ۷/۸ هرتز تعیین شده است. این لرزه‌نگار، شتاب زمین‌لرزه را به‌ عنوان ورودی گرفته و از طریق سیستم مکانیکی به کرنش در توری براگ فیبر نوری منتقل می‌کند و خروجی آن جابه‌جایی طول موج براگ می‌باشد. به عنوان نمونه، با به‌کارگیری طیف بسامدی شتاب زمین‌لرزه منطقه ریگان استان کرمان که در سال ۲۰۱۱ میلادی روی داده، بسامدها و دامنه‌های شتاب زمین‌لرزه‌ها محاسبه شده است. سپس، کرنش و جابه‌جایی دامنه طول موج براگ و جابه‌جایی فاز به ازای ضریب میرایی‌های مختلف دستگاه و به ازای بسامدها و دامنه‌های مختلف شتاب زمین‌لرزه‌ها بررسی گردیده و برحسب بسامدهای مختلف شتاب زمین‌لرزه رسم شده‌اند. همچنین برای اثبات کارایی و قابلیت لرزه‌نگار توری براگ فیبر نوری دو زمین‌لرزه دیگر یکی در منطقه ریگان در سال ۲۰۱۰ میلادی و دیگری زمین‌لرزه سراوان در سال ۲۰۱۱ میلادی به روش بالا بررسی شده است. همچنین هر کدام از رویدادهای لرزه‌ایی توسط چند ایستگاه دیگر ثبت و پردازش و همزمان با لرزه‌نگار توری براگ فیبر نوری مقایسه شده است.}, keywords_fa = {زمین‌لرزه,تنش,تنجش,لرزه‌نگار,توری براگ فیبر نوری,فلز ارتجاعی}, url = {https://jesphys.ut.ac.ir/article_58137.html}, eprint = {https://jesphys.ut.ac.ir/article_58137_35296928397d1ca0abebc81b9e924e3c.pdf} } @article { author = {Nikrouz, Ramin and Jahanjouy, Saber and Siahkohi, Hamid Reza}, title = {Using Fourier based vectors in Matching Pursuit (MP) algorithm to interpolate seismic data and reconstruct traces on uniform grids}, journal = {Journal of the Earth and Space Physics}, volume = {42}, number = {3}, pages = {499-511}, year = {2016}, publisher = {Institute of Geophysics, University of Tehran}, issn = {2538-371X}, eissn = {2538-3906}, doi = {10.22059/jesphys.2016.58889}, abstract = {Many steps of seismic data processing sequence suppose that data sets are sampled in time and spatial dimensions uniformly. Today, this assumption is true but only in time dimension. Modern seismic exploration equipment permits seismic data sets to be sampled uniformly and densely in time dimension. However, along spatial dimension uniform and dense sampling are not possible because of operating constraints, failure of equipment, topography conditions or commercial problems. It has been proved that the results of most of seismic data processing techniques are dependent on regularity, adequate sampling and density of input data sets. The fact that we need to interpolate seismic data sets causes several new-born approaches in this field. In most of the available seismic processing software, this task is done by ‘binning’ the data. This operation is one of the error sources of seismic sections. Moreover, there are some other different computational techniques to interpolate and reconstruct seismic data on a regular grid. Some of these approaches reconstruct seismic data at the given points using physical concepts of wave propagation and solving Kirchhoff's formula. In spite of practicability of these methods, need of initial accurate information about velocity model, geological property and high computational efforts restrict the domain of operation for these methods. Nowadays various mathematical methods are provided using the design of prediction filters, mathematical transformation and some other methods use rank reduction of data matrix to interpolate seismic data. According to their utilized assumptions, computational cost, noise, sampling type, and density of input data, each of these methods have their own constraints in performance and artifacts in final results which should be recognized. In science and engineering branches, a well-known algorithm that deals with signals is Matching Pursuit (MP). Originally, MP has been introduced to time-frequency transformation and finding the frequency content of signals. This transformation represents a signal as a linear composition of vectors that are available in a complete bank of time-frequency atoms (also called Dictionary). MP is an iterative algorithm that at each iteration finds a base vector in the dictionary that best matches to the signal, then subtracts the image of signal along this vector from the signal and updates the signal. This process will be continued until the remained signal is negligible. Originally, to have a good decomposition, this dictionary should contain a vast amount and kinds of wavelets like Gabor functions that each has its own dilation, modulation and translation.Heretofore MP is used to produce a single frequency seismic attribute in geophysics. For seismic data reconstruction and interpolation purposes, sine functions are applied as base vectors. The process of interpolation by MP that uses sine functions needs to solve a Lomb-Scargle periodogram at each iteration that may need to have many computations. Due the lots of works that have been done on this subject, today multi-dimension and multi-component seismic data set can be interpolated using sine functions at MP. Other functions that can be used as MP’s base vectors are Fourier coefficients. Here, after some brief explanation about MP’s algorithm and formulations we use Fourier coefficients as the base vectors of MP, interpolate and reconstruct some synthetic and real two and three dimensional seismic data. Despite of some random noises that are due to calculation and other estimations,the traces are reproduced acceptably. The results show that amplitude and frequency contents of events are well preserved. The noticeable point is that the traces that reproduced at original sampling points are nearly identical to original traces. This property and ability to interpolate data with completely non-uniform sampling grid separates Fourier MP from many of previous interpolation methods. Cautiously picking of several base function simultaneously is proposed to reduce needed iterations and speed up the algorithm. Windowing the input data and using an antialiasing mask are proposed to achieve the assumption of sparse frequency content and linearity of events and remove aliasing effects.  }, keywords = {Matching Pursuit,Time-frequency Transform,interpolation,Reconstruction,non uniform sampling}, title_fa = {استفاده از بردارهای پایه ی فوریه در الگوریتم تعقیب تطابق برای درون یابی داده ی لرزه ای و بازسازی ردلرزه ها بر روی شبکه های منظم}, abstract_fa = {در اغلب مراحل پردازش داده‌های لرزه‌ای فرض بر یکنواخت بودن نمونه برداری زمانی و مکانی است. از اینرو بازسازی داده‌های لرزه‌ای مفقود، درونیابی گاف بین ردلرزه‌ها و مرتب نمودن ردلرزه ها بر روی یک شبکة منظم و متراکم یکی از گام‌های مهم در فرآیند پردازش داده‌های لرزه‌ای جهت بالا بردن دقّت و صحت مراحل بعدی پردازش است. اهمیت این موضوع سبب معرفی روش‌های زیادی در این راستا شده است. ایده‌ی اصلی برخی از روش‌های درونیابی و بازسازی پیدا کردن الگویی است که تغییرات داده را در حوزه‌ی اصلی داده و یا در حوزه‌ی یکی از تبدیلات نشان می‌دهد و از این الگو جهت درونیابی داده در نقاط دلخواه استفاده می‌شود. نمایش زمان-بسامد داده‌های لرزه‌ای یکی از حوزه‌هایی است که تاکنون جهت بررسی خواص ژئوفیزیکی ساختارها و مخازن مورد استفاده قرار گرفته است. در این تحقیق پس از بررسی مختصر آلگوریتم تعقیب تطابق که یکی از روش‌های موجود برای یافتن توزیع زمان-بسامد داده‌ها در علوم و مهندسی است، آلگوریتمی برای بازسازی داده‌های لرزه‌ای با استفاده از تعقیب تطابق و بردارهای پایه‌ی فوریه شرح داده شده است. نتایج بدست آمده از درونیابی داده‌های لرزه‌ای مصنوعی و واقعی و بازسازی آنها بر روی شبکه‌ی منظم حاکی از توانایی روش معرفی شده در بازسازی ساختارهای پیچیده و شیبدار و افزایش دقت نتایج با افزایش ابعاد دادة ورودی است.}, keywords_fa = {تعقیب تطابق,تبدیل زمان بسامد,درونیابی,بازسازی,نمونه برداری نامنظم}, url = {https://jesphys.ut.ac.ir/article_58889.html}, eprint = {https://jesphys.ut.ac.ir/article_58889_bf08bc7032f8206d66be41fb6f2220ff.pdf} } @article { author = {Khas Ahmadi, Shahriyar and Gholami, Ali}, title = {Seismic velocity analysis using low-rank approximation of the Kernel function and a Butterfly algorithm}, journal = {Journal of the Earth and Space Physics}, volume = {42}, number = {3}, pages = {513-522}, year = {2016}, publisher = {Institute of Geophysics, University of Tehran}, issn = {2538-371X}, eissn = {2538-3906}, doi = {10.22059/jesphys.2016.58907}, abstract = {Velocity analysis is one of the most important step in seismic data processing. It affects not only many processing steps directly and indirectly, but also is known as a primary interpretation of the data. However, it can also be assumed as one of the most time consuming processing step. The conventional velocity analysis method measures the energy amplitude along hyperbolic trajectories within a velocity interval and creates a velocity model. In this procedure, the data from time-offset domain is mapped to time-velocity or time-slowness domain. For a number of , and time, offset and velocity samples respectively,  computations is necessary to obtain a velocity model. However, in the presence of large size data and model parameters, computing the velocity spectrum using conventional method would be a time consuming task. On the other hand, in order to improve the initial velocity model obtained in the processing steps, usually velocity analysis is conducted several times during the processing of the seismic data. Hence, there should be a better way to compute the velocity model in a much less time computation. In this paper, we introduce the Butterfly algorithm for fast computation of hyperbolic Radon transform (HRT), as a kind of time variant operator, with an application in seismic velocity analysis. In seismic data processing, Radon transforms map the overlapping data in seismic gathers to another domain which they can be separated. Among different types of Radon transforms, the HRT has the most similarity to the seismic events and hence, produce the most accurate approximation in the velocity spectrum. However, its time-variant kernel prohibits its fast computation especially for large size data. Unlike time-invariant operators which use the convolution theorem in the Fourier domain to compute the velocity domain for each frequency separately and therefore efficiently, Fourier transform of time-variant operators is a function of both frequency and time and using the convolution theorem is not applicable. The Butterfly algorithm can be used as a fast solver for the Fourier Integral Operators (FIO), so reformulating the HRT integral in the Fourier domain as FIO makes it possible to use this algorithm to overcome the problem of the time-variant kernel. The basis of this solution is the existence of low-rank approximations of the kernel when it is restricted to subdomains in data and model spaces. Subdividing the model and data domain properly to smaller subdomains admits low-rank approximations of the kernel. These low-rank approximations enable us to obtain functions of only one variable, time or frequency, which approximate the kernel. This decoupling of time and frequency variables allows fast computation of the HRT integral. In order to do the subdivision properly, a pair of quad trees, one for each data and model domains, is used to restrict the domains in a level-base structure in which the size of data domain subsets are increasing while the size of model domain subsets are decreasing in each level. The Butterfly algorithm is used to compute the kernel equivalent functions in each level of these quad trees for each subdomain. Finally, at the last level, the Radon panel or velocity model is obtained. The complexity of this method for two dimensional data is  in which depends on data and model variables range. As it was demonstrated in the synthetic and the real numerical examples,  complexity results in reduction of computation time in several orders relative to the conventional method.}, keywords = {Seismic velocity analysis,hyperbolic radon transform,Butterfly algorithm,low-rank approximation}, title_fa = {تحلیل سرعت لرزه‌ای با استفاده از تقریب رتبه کم تابع هسته ای و الگوریتم پروانه‌ای}, abstract_fa = {تحلیل سرعت را می‌توان یکی‌ از اساسی‌‌ترین مراحل پردازش داده‌های لرزه‌ای دانست چرا که نه تنها بسیاری از مراحل پردازش را به صورت مستقیم و غیرمستقیم تحت تاثیر قرار می‌دهد، بلکه می‌توان آن را تفسیری اولیه‌ از داده‌ها تلقی کرد. اما همچنان می‌توان آنرا یکی‌ از زمان‌گیرترین مراحل پردازش نیز دانست. روش معمول تحلیل سرعت با اندازه‌گیری دامنه انرژی در امتداد مسیر‌های هذلولی شکل و به ازای یک بازه از مقادیر سرعت سعی‌ در بدست آوردن یک مدل سرعت دارد که در صورت بزرگ بودن ابعاد داده محاسبه طیف سرعت به وسیله این روش بسیار زمان‌گیر خواهد بود. در این مقاله به معرفی‌ الگوریتم پروانه‌ای (Butterfly Algorithm) در حل سریع تبدیل رادون هذلولی به منظور تحلیل سرعت داده‌های لرزه‌ای می‌پردازیم. اساس این الگوریتم تغییر رابطه تبدیل رادون هذلولی به صورت یک عملگر انتگرال فوریه و سپس بدست آوردن تقریب‌هایی‌ با رتبه پائین (low-rank approximation) برای کرنل مسئله خواهد بود. همگرایی این روش برای داده‌های دو بعدی برابر خواهد بود که به بازه تغییرات متغیرهای داده و مدل بستگی دارد. همانطور که در مثال‌های عددی غیرواقعی و واقعی نشان داده خواهد شد، این همگرایی باعث کاهش زمان محاسبه مدل سرعت تا چندین برابر نسبت به روش معمول می‌شود.}, keywords_fa = {تحلیل سرعت لرزه ای,تبدیل رادون هذلولی,الگوریتم پروانه‌ای,تقریب رتبه پایین}, url = {https://jesphys.ut.ac.ir/article_58907.html}, eprint = {https://jesphys.ut.ac.ir/article_58907_9d5f07c099a2532ae07de25f807236fa.pdf} } @article { author = {ابراهیم زاده اردستانی, وحید and Vatankhah, Saeed and سلطان آبادی, رضوان}, title = {Comparison between Tikhonov regularization and truncated SVD in gravity data inversion}, journal = {Journal of the Earth and Space Physics}, volume = {42}, number = {3}, pages = {523-534}, year = {2016}, publisher = {Institute of Geophysics, University of Tehran}, issn = {2538-371X}, eissn = {2538-3906}, doi = {10.22059/jesphys.2016.57794}, abstract = {In this paper the 3D inversion of gravity data using two different regularization methods, namely Tikhonov regularization and truncated singular value decomposition (TSVD), is considered. The earth under the survey area is modeled using a large number of rectangular prisms, in which the size of the prisms are kept fixed during the inversion and the values of densities of the prisms are the model parameters to be determined. A depth weighting matrix is used to counteract the natural decay of the kernel, so the inversion obtains reliable information about the source distribution with respect to depth. To generate a sharp and focused model, the minimum support (MS) constraint is used, which minimizes the total area with non zero departure of the model parameters from a given a priori model. Then, the application of iteratively reweighted least square algorithm is required to deal with non-linearity introduced by MS constraint. At each iteration of the inversion, a priori variable weighting matrix is updated using model parameters obtained at the previous iteration. We use the singular value decomposition (SVD) for computing Tikhonov solution, which also helps us to compare the results with the solution obtained by TSVD. Thus, the algorithms presented here are suitable for small to moderate size problems, where it is feasible to compute the SVD. In Tikhonov regularization method, the optimal regularization parameter at each iteration is obtained by application of the parameter-choice method. The method is based on the statistical distribution of the minimum of the Tikhonov function. For weighting of the data fidelity by a known Gaussian noise distribution on the measured data and, when the regularization term is considered to be weighted by unknown inverse covariance information on the model parameters, the minimum of the Tikhonov functional becomes a random variable that follows a  distribution. Then, a Newton root-finding algorithm can be used to find the regularization parameter. For truncated SVD regularization, the Picard plot is used to find a suitable value of truncation index. In math literature, a plot of singular values together with SVD and solution coefficients is often referred to as Picard plot. To test the algorithms, a density model which consists of a dipping dike embedded in a uniform half-space is used. The surface gravity anomaly produced by this model is contaminated with three different noise levels, and are used as input for introduced inversion algorithms. The results indicate that the algorithms are able to recover the geometry and density distribution of the original model. In general, the reconstructed model is more sparse using TSVD method as compare with Tikhonov solution. This especially happens for high noise level, where there is an important difference between two solutions. In this case, while TSVD produces a sparse model, the solution of Tikhonov regularization is not sparse. Furthermore, the number of iterations, which is required to terminate the algorithms, is more for TSVD as compare with Tikhonov method. This feature, along with automatic determination of regularization parameter, makes the implementation of the Tikhonov regularization method faster than TSVD. The inversion methods are used on real gravity data acquired over the Gotvand dam site in the south-west of Iran. Tertiary deposits of the Gachsaran formation are the dominant geological structure in this area, and it is mainly comprised of marl, gypsum, anhydrite and halite. There are several solution cavities in the area so that relative negative anomalies are distinguishable in the residual map. A window of residual map consists of 640 gridded data, which includes three negative anomalies, that is selected for modeling. The reconstructed models are shown and compare with results obtained by bore holes.}, keywords = {Gravimetry,inversion,Tikhonove Regularization,TSVD,Regularization parameter,Gotvand}, title_fa = {مقایسه بین منظم‌سازی‌ تیخونف و تجزیه مقدار تکین قطع شده در وارون‌سازی داده‌های گرانی سنجی}, abstract_fa = {در این مقاله وارون سازی سه بعدی گرانی سنجی با استفاده از دو روش منظم سازی تیخونف و تجزیه مقدار تکین قطع شده بررسی شده است. ماتریس وزن دهی عمقی و قید Minimum Support در وارون سازی مورد استفاده قرار گرفته اند. در منظم سازی تیخونف از روش برای تعیین پارامتر تنظیم استفاده شده است در حالیکه ترسیم پیکارد برای تعیین اندیس قطع در روش تجزیه مقدار تکین قطع شده به کار رفته است. داده مصنوعی تولید شده توسط یک دایک شیب دار با اضافه کردن نوفه در سه سطح مختلف به عنوان ورودی برای وارون سازی استفاده شده است. نتایج حاصل از مدل سازی دلالت بر آن دارد که الگوریتم های ارائه شده به طور کلی توانایی خوبی برای بازسازی شیب ومرزهای مدل دارند. این نتایج نشان می دهند که تُنکی مدل با استفاده از روش تجزیه مقدار تکین قطع شده بیشتر است، همزمان تعداد تکرارهای بیشتری برای توقف وارون سازی در این روش مورد نیاز است. روش تیخونف برای نوفه سطح بالا مدلی هموار بازسازی کرده است در حالیکه برای دو سطح نوفه دیگر مدل ساخته شده بسیار نزدیک به مدل اصلی می باشد. در روش تجزیه مقدار تکین قطع شده، در مقایسه با روش تیخونف، نیاز به ترسیم پیکارد سبب کاهش سرعت اجرای الگوریتم می گردد. الگوریتم های ارائه شده بر روی داده برداشت شده در ناحیه سد گتوند به کار رفته اند و نتایج حاصل از مدل سازی نشان داده شده است.}, keywords_fa = {گرانی سنجی,وارون سازی,منظم سازی تیخونف,تجزیه مقدار تکین قطع شده,پارامتر تنظیم,گتوند}, url = {https://jesphys.ut.ac.ir/article_57794.html}, eprint = {https://jesphys.ut.ac.ir/article_57794_42040ba3a4df688e1483a506b0be2012.pdf} } @article { author = {Hajian, Alireza and زمردیان, حسین}, title = {Estimation of depth and shape of subsurface cavities via Multi Adaptive Neuro-Fuzzy Interference System using Gravity data}, journal = {Journal of the Earth and Space Physics}, volume = {42}, number = {3}, pages = {535-548}, year = {2016}, publisher = {Institute of Geophysics, University of Tehran}, issn = {2538-371X}, eissn = {2538-3906}, doi = {10.22059/jesphys.2016.55631}, abstract = {In common classical methods of cavity depth estimation through microgravity data, usually when a pre-geometrical model is considered for the cavity shape, the simple geometrical models of sphere, vertical cylinder and horizontal cylinder are commonly used. It is obviously an important fact that in real conditions the shapes of the cavities are not exactly sphere, horizontal cylinder or vertical cylinder but are near or to some extent near to these simple models. The linguistic variables “near to” or “to some extent near to” are consisting of fuzzy concepts and can be described as “fuzzy” variables. The membership degree of each fuzzy variable shows how much the variable is near to the mentioned simple shapes. Using the fuzzy variables leads to better results with more accuracy because in real conditions the nature of the cavities shape is “fuzzy” so that their shape is not exactly but near to the mentioned simple shapes. Consequently, in this paper in order to help the gravity data interpreter to enhance the accuracy of his/her interpretation a neuro-fuzzy model namely Multi Adaptive Neuro-Fuzzy Interference System (MANFIS) is used. When the neural network alone is used the challenge is its black-box operation so that there is no possibility for sensitive analysis but neuro-fuzzy networks consist of the sensitive analysis via the if-then fuzzy rules achieved during the training process. In Multi Adaptive Neuro-Fuzzy Interference System, the network output before the de-fuzzification stage, shows the interpreter how much the cavity shape is near to sphere, horizontal cylinder or vertical cylinder. In this research, two Adaptive Neuro-Fuzzy Interference System (ANFIS) models were paralleled to configure a Multi Adaptive Neuro-Fuzzy Interference System (MANFIS) so that one output of the designed MANFIS is the shape factor and the other is the depth of the cavity. The inputs of the MANFIS are some of the important features selected from the gravity signal along the selected principle profiles of the residual gravity map. In order to evaluate the designed MANFIS in the presence of noise in gravity data, the method was tested for synthetic data with 5% and 10% level of noise. The results showed that the joint neural networks and fuzzy logic makes it a suitable tool to help the interpreter to improve the accuracy of shape and depth estimation of cavities. Furthermore, the method is more robust to noise which were tested for two different noise levels one with low level of noise and other with medium level of noise added to the synthetic gravity data. Despite the available classical methods or net neural methods, here without any pre-assumption about the shape of the cavity, both the shape factor and depth are estimated. In is necessary to mention that the value of the estimated shape factor implies that which of the geometrical models among sphere, vertical cylinder or horizontal cylinder are better fitted to the real shape of the subsurface cavity.  After checking and confirming the accuracy of the designed MANFIS for synthetic data, the method was tested for real data through micro-gravity data over a gravity site located in Great Bahama Free Port, west of North America. The results are very near to the available borehole and extracted data.}, keywords = {Gravity,Adaptive Neuro-fuzzy Interference System,depth estimation,Shape factor}, title_fa = {برآورد عمق وشکل حفره های زیرزمینی با استفاده از دستگاه واسط عصبی فازی تطبیقی چندگانه با داده های گرانی سنجی}, abstract_fa = {در این مقاله به منظور اکتشاف حفرات زیرزمینی با شکلهای نزدیک به کره، استوانه افقی یا عمودی ودر راستای بالابردن دقت نتایج تفسیر بی هنجاریهای گرانی ،کمک به تجربیات مفسر و مقاومت بیشتر در برابر سطوح متفاوت نوفه ، از شبکه عصبی-فازی تطبیقی چند گانه MANFIS استفاده شده است. در این پژوهش با قرار گرفتن دو سیستم عصبی فازی تطبیق پذیر به صورت موازی با یکدیگر یک شبکه عصبی-فازی تطبیق پذیر چند گانه طراحی شد که خروجی یکی فاکتور شکل حفره زیرزمینی وخروجی دیگری عمق مربوط به حفره می باشد. به منظور امتحان دقت عملکرد شبکه عصبی فازی طراحی شده در حضور نویز، روش ارائه شده ابتدا برای داده های مصنوعی با 5 درصد و10 درصد نویز مورد آزمون قرار گرفت . در مجموع نتایج نشان داد استفاده توام از شبکه های عصبی ومنطق فازی علاوه بر آن که ابزاری مفید جهت کمک به مفسردر مرحله تفسیر عمق وشکل حفره های زیرزمینی از روی داده های گرانی است، بلکه صحت تفسیر بی هنجاری های گرانی را نیزافزایش می دهد . همچنین بر خلاف روشهای موجود با رهیافت عصبی محض در اینجا بدون پیش فرض شکل درباره چشمه گرانی امکان تخمین شکل چشمه علاوه بر تخمین عمق آن وجود دارد. پس از اطمینان از صحت عملکرد شبکه عصبی- فازی طراحی شده برای داده های مصنوعی، به منظور امتحان روی داده های واقعی برای داده های گرانی سایت بند باهاما تست ومورد ارزیابی قرار گرفت که با نتایج واقعی حاصل از گمانه زنی ها و حفاری های موجود تطابق خوبی دارد}, keywords_fa = {تخمین عمق,شبکه عصبی-فازی,ANFIS,گرانی سنجی}, url = {https://jesphys.ut.ac.ir/article_55631.html}, eprint = {https://jesphys.ut.ac.ir/article_55631_be791cef569f51d840c05462f61d43c9.pdf} } @article { author = {Motavalli Anbaran, Seyed Hani and انتظارسعادت, وحید and جاماسب, علی}, title = {Isostasy of Iranian Plateau}, journal = {Journal of the Earth and Space Physics}, volume = {42}, number = {3}, pages = {549-562}, year = {2016}, publisher = {Institute of Geophysics, University of Tehran}, issn = {2538-371X}, eissn = {2538-3906}, doi = {10.22059/jesphys.2016.55835}, abstract = {Using global free air gravity and topography data, first we have calculated a Bouguer Anomaly (BA) map for Iranian plateau and then computed the residual isostasy anomaly map under the Airy-Heiskannen assumption. The term residual is used as to reflect the assumption of local isostatic compensation in contrast to a regional isostatic compensation. The value with which the gravity effect of the compensation mass (the root/anti-root in the Airy model) is calculated are chosen under careful considerations as to produce reliable results. The resulting residual isostasy map is then used to qualitatively interpret the isostatic highs and lows corresponding to crustal/lithospheric features of the Iranian plateau. The study area is a complex region as a result of its still active tectonics which is mostly driven by the continent-continent collision of the Arabian and Eurasian plates. The five most important tectonic settings in Iran are Zagros Mountains, an active belt formed as the result of the collision extending from south-west Iran along the Persian Gulf; Alborz Mountains, a young belt with an average topography of 3-5 km extending nearly in east-west direction, Makran, in south-east Iran, north of the Iran-Arabian plate boundary where an active subduction is taking place; Caspian Sea, with an oceanic crust covered with an average 15-20 km sediment layer at the Iran-Eurasia plate boundaries and Kopeh-Dagh mountains, an uplifted region as a consequence of converging continental plates. Our results indicate that the Zagros Mountains have reached an isostatic equilibrium but the scenario is slightly different for the Alborz chain. It seems that the isostatic equilibrium is not fully reached in the Alborz due to the observation of a continuous isostatic high (positive) anomaly which extends to north-west Iran, however, it may also be partly caused by a simple folding. In the southern Caspian region, there is an enormous isostatic low (negative), for the cause of which we have considered two possible reasons. First, the effect of the sediment layer on the gravity signal due to its negative density contrast. Second, we considered the deficiency in the rock mass at the base of the lithosphere due to an anti-plume or the downward flow of the lithospheric materials towards the mantle which may also explain the high depth of the southern Caspian Basin. Subduction zones are usually characterized with negative isostatic anomalies, but in the case of the subduction of the oceanic lithosphere of the Caspian under the continental crust of the Eurasia, there is no apparent negative isostatic anomaly in our map. We believe that this is probably due to the fact that the subduction is still young while in order to observe a negative effect on the residual isostasy anomaly map, the subducting slab must be in a deep position, in other words, be of older age. The subduction of the Makran, on the other hand, has caused a negative isostatic residual anomaly. This low anomaly is also partly due to the uplift of the Makran area. A high-low (positive-negative) residual isostasy anomaly pairs corresponds to suture zones. An example of which is seen for the Zagros-Bitlis suture zone which marks the continental collision of the Arabian-Eurasian plates. Our map also shows a negative residual isostatic anomaly in the Kopeh-Dagh Mountains, which we interpret as the uplift caused by the convergence of the Iranian and Eurasian plates. It must be noted that every high/low residual isostatic anomaly may not be interpreted as isostatically over/under compensated areas. on the contrary, it could be and usually is related to a geological feature of lithosphere/mantle scale.}, keywords = {Isostasy,Iranian plateau,Bouguer anomaly map,Isostatic equilibrium}, title_fa = {هم ایستایی فلات ایران}, abstract_fa = {فلات ایران یکی از مناطق پیچیده ازنظر زمین‌شناسی است که از دیرباز فعالیت‌های تکتونیکی آن موردتوجه دانشمندان علوم زمین بوده است. با استفاده از نقشه آنومالی ایزوستازی باقی‌مانده می‌توان به دید بهتری از واقعیت‌های زمین‌شناسی دست یافت و پدیده‌های زمین-شناسی و عوارض بزرگ‌مقیاس را با دقت بیشتری تحلیل کرد. برای تهیه نقشه آنومالی ایزوستازی باقی‌مانده می‌بایست ابتدا نقشه آنومالی بوگه محاسبه گردد. در این مقاله با استفاده از داده‌های ماهواره‌ای توپوگرافی و آنومالی هوای آزاد، آنومالی بوگه فلات ایران و مناطق اطراف به صورت یک‌پارچه محاسبه و در ادامه نقشه آنومالی ایزوستازی با استفاده از مدل ایری-هیسکانن تهیه گردید. درنهایت نقشه آنومالی ایزوستازی باقی‌مانده حاصل تفریق نقشه آنومالی ایزوستازی از نقشه آنومالی بوگه است.با توجه به اطلاعات به‌دست‌آمده از نقشه آنومالی ایزوستازی باقی‌مانده مشاهده می‌شود که در منطقه مکران وجود آنومالی منفی در ارتباط با وجود فرورانش و بالاآمدگی در این ناحیه است. رشته‌کوه‌های زاگرس با تقریب خوبی به تعادل ایزوستازی رسیده‌اند و جفت آنومالی مثبت منفی موجود در مرز بین دریا و خشکی در ارتباط با زون زمین‌درز بیتلیس-زاگرس است. گسل اصلی کپه‌داغ محل برخورد دو پوسته‌ی قاره‌ای اوراسیا و ایران به یکدیگر می‌باشد. در مورد رشته‌کوه‌های البرز به نظر می‌رسد که تعادل به طورکامل به دست نیامده است و در مورد خزر جفت آنومالی مثبت و منفی در نقشه آنومالی ایزوستازی باقی‌مانده نشان‌دهنده‌ی مرز صفحه‌ی ایران و صفحه‌ی اوراسیاست.}, keywords_fa = {ایزوستازی,فلات ایران,نقشه آنومالی بوگه,تعادل ایزوستازی}, url = {https://jesphys.ut.ac.ir/article_55835.html}, eprint = {https://jesphys.ut.ac.ir/article_55835_470ac5ec1f23747c99239a8a78efbfd9.pdf} } @article { author = {Izadi Yazdanabadi, Marjan and Hojat, Azadeh and Ranjbar, Hojatollah and Karimi Nasab, Saeed}, title = {Important parameters to locate new geomagnetic observatories to develop geomagnetic studies in Iran}, journal = {Journal of the Earth and Space Physics}, volume = {42}, number = {3}, pages = {563-574}, year = {2016}, publisher = {Institute of Geophysics, University of Tehran}, issn = {2538-371X}, eissn = {2538-3906}, doi = {10.22059/jesphys.2016.59866}, abstract = {Geomagnetic observatories are constructed to continuously record the earth’s geomagnetic field. The importance of such buildings has resulted in an increasing construction of geomagnetic observatories all over the world. Although geomagnetic observatories are sometimes unknown to most people and even some scientists, these very important data centers have been measuring the geomagnetic field for about 500 years. The number of standard observatories is now approaching 150, but there are still very few standard geomagnetic observatories in the Middle East. Despite the large area of Iran, Tehran geomagnetic observatory is the only observatory in the country. It was constructed by the Institute of Geophysics in 1961. Unfortunately, in the recent decade, due to the old instruments used in this observatory and considering the expansion of Tehran city, the data obtained at this geomagnetic observatory are subjected to many noise sources. Therefore, its data are not reliable anymore and Tehran observatory is removed from the list of world’s standard geomagnetic observatories. In this paper, which is aimed as an informative research note, the necessity of constructing new geomagnetic observatories in Iran is discussed from various aspects. The most important advantages of constructing a geomagnetic observatory include: providing a data center for the correction of magnetic exploration data, prediction of magnetic storms and providing a constant monitoring of the geomagnetic field variations as one of the probable earthquake precursors. Considering the various advantages of continuous recording of the earth’s geomagnetic field, it is completely necessary to construct new geomagnetic observatories in Iran. The first step to construct a standard geomagnetic observatory that can record highly reliable data, is to select an optimum area for the construction of the observatory. From the very early observations of the geomagnetic field, it was understood that the observatories must be far enough from natural fluctuations resulted from volcanic rocks or mineral deposits. It was also discovered that artificial noises can disturb the data recorded in a geomagnetic observatory. However, there is still no comprehensive report discussing the various criteria to be studied while selecting a site for the construction of these observatories. This paper presents different parameters that should be carefully investigated to select optimum sites for the construction of geomagnetic observatories. Considering the type of data that should be recorded in geomagnetic observatories, knowledge of the criteria affecting the geomagnetic data is critical. Various criteria such as the magnetic intensity, variations of the electrical conductivity of subsurface soils, artificial sources of geomagnetic changes, development of cities, geology, topography, access roads, and underground possible economical deposits can affect selection of the optimum sites for the construction of geomagnetic observatories. Keeping a proper distance away from magnetic anomalies and controlling the homogeneity of the electrical conductivity of subsurface soils fall among the most important factors to be considered at the first stage. Then, artificial sources of the geomagnetic disturbance and urban development patterns must be carefully considered. It should be noted that finding the areas that can fully satisfy all the criteria might be impossible, but the minimum requirements should be satisfied to construct an observatory. Closure or relocation of several geomagnetic observatories all over the world is an experience showing how carefully the site selection of these structures should be carried out. Having understood the necessity of constructing new geomagnetic observatories in Iran, Kerman Province is considered as one of the favorable areas for the construction of an observatory. A comprehensive research is being carried out by authors to carefully acquire and interpret all the required data to find the best sites for this purpose.}, keywords = {continuous recording,field component,geomagnetic field,geomagnetic observatory,Iran,Site Selection}, title_fa = {پارامترهای مهم در مکان یابی رصدخانه های ژئومغناطیسی جدید به منظور توسعه پژوهش های ژئومغناطیسی در ایران}, abstract_fa = {علی‌رغم وسعت زیاد کشور ایران، رصدخانه ژئومغناطیسی تهران تنها رصدخانه در کشور می‌باشد که در دهه گذشته، به‌دلیل قدیمی بودن دستگاه‌های موجود در این رصدخانه و نیز گسترش شهر و قرارگیری این مرکز در معرض منابع نوفه‌ای متعدد، متاسفانه داده‌های آن مورد اعتماد نبوده و از لیست رصدخانه‌های ژئومغناطیسی استاندارد دنیا حذف شده است. با توجه به مزایای متعدد ثبت پیوسته میدان مغناطیسی و مؤلفه‌های آن در رصدخانه‌های ژئومغناطیسی، انجام مطالعات علمی برای تعیین محل‌های بهینه جهت احداث رصدخانه‌های ژئومغناطیسی جدید در بسیاری از کشورهای خاورمیانه و به‌خصوص ایران، بسیار ضروری می‌باشد. در این مقاله، ضمن بیان ضرورت احداث رصدخانه‌های ژئومغناطیسی جدید در کشور، به معرفی مهم‌ترین معیارهای موثر در تعیین محل احداث رصدخانه-های ژئومغناطیسی پرداخته شده است. معیارهای متعددی از قبیل شدت میدان مغناطیسی، تغییرات هدایت‌ویژه الکتریکی خاک‌های زیرسطحی، منابع غیرطبیعی اختلال در میدان ژئومغناطیسی، توسعه شهرسازی، راه‌های دسترسی، زمین‌شناسی منطقه، توپوگرافی منطقه و عدم وجود ذخائر اقتصادی، در انتخاب محل مناسب جهت احداث یک رصدخانه ژئومغناطیسی تأثیرگذار هستند. ممکن است گردآوری کلیه شرایط لازم برای یک محل، به‌طور کامل میسر نباشد، اما چنانچه حداقل شرایط لازم نیز میسر نگردد، به هیچ وجه نباید اقدام به ساخت رصدخانه نمود. تعطیلی و جابه‌جایی بسیاری از رصدخانه‌ها تجربه‌ای است تا مکان رصدخانه‌های مغناطیسی جدید با دقت زیادی انتخاب شود. در ادامه مطالعات این یادداشت پژوهشی، تحقیقی منسجم توسط نویسندگان مقاله در حال انجام است که به جمع‌آوری دقیق لایه‌های اطلاعاتی مختلف جهت تعیین محل‌های مناسب برای احداث یک رصدخانه ژئومغناطیسی در استان کرمان می‌پردازد.}, keywords_fa = {رصدخانه ژئومغناطیسی,میدان مغناطیسی,مؤلفه میدان,ثبت پیوسته,تعیین محل,ایران}, url = {https://jesphys.ut.ac.ir/article_59866.html}, eprint = {https://jesphys.ut.ac.ir/article_59866_498538b66855ac31cfdd26471b6988fc.pdf} } @article { author = {Oskooi, Behrooz and Mansoori, Isa and پدرسن, لاوست and جواهری, رحمان}, title = {2-D and 3-D inversion of the magnetotelluric data to explore hydrocarbon structures in the Sehqanat oil field, SW Iran}, journal = {Journal of the Earth and Space Physics}, volume = {42}, number = {3}, pages = {575-586}, year = {2016}, publisher = {Institute of Geophysics, University of Tehran}, issn = {2538-371X}, eissn = {2538-3906}, doi = {10.22059/jesphys.2016.58919}, abstract = {Among all the geophysical techniques, the magnetotelluric method has improved considerably in recent years and is widely being used in hydrocarbon exploration especially in regions where reflection seismic has difficulties. Areas which are covered with high velocity rocks in the near surface are most popular cases. A huge high resolution magnetotelluric investigation was conducted in the Sehqanat oil field, SW of Iran, in 2013 to map geoelectrical structure of the region from surface down to several kilometers. The Sehqanat oil field is located in sedimentary Zagros zone which is encompasses more than 95 percent of Iran’s oil fields. The main geological interface which is targeted to be imaged with magnetotelluric method, due to the large resistivity contrast (based on the well logs information), is the contact between the highly conductive evaporites of Gachsaran formation and the more resistive underlying carbonates of Asmari formation. Regarding the large thickness of the high-velocity (ca. 4500 m/s) and heterogeneous Gachsaran Formation outcropping in the Sehqanat oil field and several adjoining oil fields in the study area, imaging of the underlying layers is difficult with the reflection seismic technique. On the other hand, the big contrast of the electrical resistivity between the Gachsaran Formation and the underlying layers is favourable for MT exploration. The geoelectrical contrast is well documented from the full-set log measured along the explorative Sehqanat well. The high velocity and very heterogeneous Gachsaran formation is exposed on the surface and has a varying thickness from 500 meter to more than two kilometers in the region and also covers the Asmari formation which is the main reservoir in SW oil fields of Iran, as a cap rock. Geologically, the Sehqanat oil field has been formed by a gentle and moderate-size anticline called “Sehqanat” which its structural shape, due to the low quality of reflection seismic data, is not clearly known for geologists. The Sehqanat anticline acts as a structural oil trap from aspect of the petroleum geology. In order to collect more geophysical information about the subsurface morphology of the Gachsaran-Asmari formations boundary as well as Sehqanat anticline, broadband magnetotelluric data were acquired at more than 600 stations along five parallel southwest-northeast profiles crossing the main geological trend of the study area. Transient electromagnetic data were also acquired over 400 stations along the mentioned profiles to be used for static correction of magnetotelluric data. Dimensionality and strike analysis of the MT data show 3-D effects in a considerable amount of sites and periods. Therefore in order to get a comprehensive view through the subsurface resistivity distribution of the Sehqanat oil field, two- and three-dimensional inversions were performed on the magnetotelluric data. The 2-D and more precisely 3-D resistivity models, resolved the Gachsaran-Asmari formations boundary as a transition zone from high conductivity to more resistivity range. The Sehqanat anticline has also been delineated throughout the 2-D and 3-D resistivity models as a resistive dome-shaped body corresponded to the middle parts of MT acquisition profiles. Correlation of the magnetotelluric resistivity models with the adjacent 2-D reflection seismic sections is remarkable, letting us to accomplish more reliable interpretation of subsurface geology of the survey area.  }, keywords = {Asmari Formation,Gachsaran formation,magnetotelluric,reflection seismic,Sehqanat oil field,transient electromagnetic}, title_fa = {وارون‌سازی دو و سه‌بُعدی داده‌های مگنتوتلوریک به‌منظور اکتشاف ساختارهای هیدروکربنی در میدان نفتی سه‌قنات، جنوب‌غربی ایران}, abstract_fa = {کاربرد روش مگنتوتلوریک در زمینه اکتشافات هیدروکربوری در مناطقی که روش‌ لرزه‌نگاری انعکاسی با مشکلاتی مواجه شده، قابل توجه است. در همین راستا یک مطالعه مگنتوتلوریک با حجم وسیع در میدان نفتی سه‌قنات متشکل از یک نفتگیر تاقدیسی، و در حوضه رسوبی زاگرس واقع در جنوب غربی ایران برای به نقشه در آوردن ساختار ژئوالکتریکی منطقه انجام شد. سازند کربناته آسماری، سنگ مخزن و سازند تبخیری گچساران، سنگ پوش را تشکیل می‌دهد. در میدان نفتی سه‌قنات، رخنمون سازند تبخیری گچساران که دارای خواص پلاستیک و ناهمگن می‌باشد، باعث شده است تا کیفیت داده‌های لرزه‌ای افت کرده و تفسیر مقاطع لرزه‌ای برانبارش با مشکلاتی مواجه شود. نگار‌های چاه‌پیمایی استخراجی از چاه اکتشافی سه‌قنات، به دلیل تباین مقاومت ویژه الکتریکی سازند گچساران با سازند‌های زیرین، استفاده از روش مگنتوتلوریک را به عنوان مکمل داده‌های روش لرزه‌نگاری توجیه می‌کند. به منظور انجام تصحیح جابجایی پایا نیز ، اندازه‌گیری‌های الکترومغناطیسی در حوزه زمان انجام شده است. مطالعه بعد و امتداد داده‌های مگنتوتلوریک حاکی از وجود مقادیر متنابهی از تاثیر ساختارهای سه‌بعدی در منطقه مورد مطالعه می‌باشد. با انجام وارون‌سازی دو و سه بعدی داده‌های مگنتوتلوریک، مرز سازندهای گچساران و آسماری به مانند تاقدیس نفتی سه‌قنات در مدل‌های مقاومت ویژه الکتریکی آشکار شده است. در مدل‌های مذکور، مقادیر مقاومت ویژه الکتریکی در عبور از سازند گچساران به سازند آسماری با افزایش قابل ملاحظه‌ای همراه می‌باشد. همبستگی مدل‌های مقاومت ویژه الکتریکی با مقاطع برانبارش لرزه‌ای تفسیر جامعی از تشکیلات زمین‌شناسی تحت‌الارضی میدان نفتی سه‌قنات میسر کرده و نقاط مبهم نتایج روش لرزه‌نگاری انعکاسی مرتفع گردیده است.}, keywords_fa = {سازند آسماری,سازند گچساران,مقاومت ویژه الکتریکی,مگنتوتلوریک,میدان نفتی سه‌قنات,وارون‌سازی}, url = {https://jesphys.ut.ac.ir/article_58919.html}, eprint = {https://jesphys.ut.ac.ir/article_58919_6c0361883414f83eca84d18e3c470f11.pdf} } @article { author = {Khazaei, Safa and khazai, safe}, title = {Detection of underground tunnels using electrical resistivity and refraction seismic tomography methods}, journal = {Journal of the Earth and Space Physics}, volume = {42}, number = {3}, pages = {587-606}, year = {2016}, publisher = {Institute of Geophysics, University of Tehran}, issn = {2538-371X}, eissn = {2538-3906}, doi = {10.22059/jesphys.2016.58914}, abstract = {Geophysical techniques remain the only ways to remotely and non-destructively sense the earth’s near subsurface and as such have the most promising prospect for rapid and accurate detection of underground tunnels. Today, electrical resistivity and seismic refraction geophysical methods have been greatly developed to identify structures and underground cavities. In this study, the ability of these methods to detect tunnels has been investigated using a case study. Seismic methods are sensitive to velocity and density changes of the rock, while the electrical response is dependent upon the electrical resistivity of the rock. In this paper, we present a case study using electrical resistivity tomography (ERT) data and refraction seismic tomography (SRT) data of a tunnel site. Also, to evaluate the capabilities of resistivity method as the main method to detect buried structures in this paper, geo-electrical abnormalities of a rectangular block through various simulations are examined. Electrical resistivity methods utilize direct currents or low frequency alternating currents to investigate the electrical properties of the subsurface. In the resistivity method, the source is artificially-generated electric current introduced into the ground using electrodes. The potential differences are measured at the surface and the pattern of potential differences provides information on the distribution of subsurface electrical resistivity. In near-surface refraction tomography, the travel times of seismic energy recorded at the surface by multiple source-receiver combinations are used to generate an optimized model of the distribution of seismic velocities in the subsurface. In ERT, the forward problem uses the finite-element method to compute the electric potential response of the earth due to a given input electric current. The inverse algorithm iteratively finds the best distribution of subsurface resistivity that best fits the observed data. The purpose of this research is to determine the location of underground tunnel using inverse modeling of ERT and SRT data, and evaluate the capabilities of resistivity method to determine the location of underground tunnels using geo-electrical abnormalities of a rectangular block through various simulations. Geophysical field surveys were performed at a site with a known tunnel. The tunnel is a 1 mx1.6 m concrete lined tunnel about 80m long. The data presented here was collected where the tunnel is at a depth of about 6m. Based on prior knowledge of the tunnel location the surveys are approximately perpendicular to the tunnel and were purposely centered on the approximate location of the tunnel in order maximize the geophysical sampling in the vicinity of the tunnel. The seismic refraction survey was performed using 96 geophones. A geophone spacing of 0.5m was used for a total spread length of 47.5m. The electrode layout consisted of 50 electrodes in a 1m dipole-dipole configuration for a total 49m spread length. In this study, we chose to use the dipole-dipole configuration due to its good lateral resolution. In this research, to determine the exact location of the tunnel, data obtained from ERT and SRT survey were inverted using Res2Dinv software and Rayfract software, respectively. The results of both methods show abnormalities in the tunnel under test site. The tunnel shows up in the electrical imaging as regions of high resistivity since both the concrete and air of the tunnel are higher resistivity than the conductive weathered rock. In practice, the resistive abnormally of the tunnel gets smoothed out and is larger than the actual tunnel. Therefore, in the ERT, the tunnel coincides with one of the high resistivity anomalies, but a second, shallower resistive abnormally of unknown provenance appears just to the east. Based on the results of the seismic survey, the velocity tomographic image is inadequate for tunnel detection as the smoothing inherent to the tomographic calculations results in only slight changes in velocity near the tunnel location. Instead, the ray coverage density mapping associated with ray tracing displays small regions of low coverage associated with the tunnel. In this instance the use of both methods would suggest that this second ERT abnormality is not a tunnel and illustrates how the use of both seismic refraction and ERT can be used to increase the reliability of detecting tunnels. Finally, with the simulation of resistivity data obtained from a rectangular block, the effects of various parameters such as depth of the tunnel, overburden conductivity, thickness of overburden on identifying underground tunnel, wasere investigated; this is done to clarify the status and ability of this method in detecting underground cavities. To approximate simulated data to the fact, five percent extra noise was added to the data. To this end geological models, which can be a target in the ground like an underground tunnel, were produced in Res2Dmod software. These synthetic models were provided with reverse modeling using Res2Dinv software. The results of simulation and modeling show that the electrical resistivity method is most widely geophysical method used for detecting underground tunnels.  }, keywords = {detection,Underground tunnels,Electrical resistivity,seismic refraction}, title_fa = {آشکارسازی تونل‌های زیرزمینی با استفاده از روش‌های توموگرافی مقاومت‌ویژه الکتریکی و لرزه‌نگاری شکستی}, abstract_fa = {روش‌های لرزه‌ای انکساری و مقاومت‌ویژه الکتریکی از توانمندترین روش‌های ژئوفیزیکی در آشکارسازی سازه‌ها و حفره‌های زیرزمینی از جمله تونل‌های زیرزمینی هستند. در این تحقیق بر اساس مدل‌سازی معکوس، قابلیت این روش‌ها برای آشکارسازی تونل با انجام یک مطالعه موردی بررسی شده است. در مطالعه‌ی موردی، داده‌های مربوط به هر دو روش تقریباً به طور هم‌زمان با آفست 3 متر برداشت شده است. نتایج به دست آمده از هر دو روش بی‌هنجاری‌هایی را در محل تونل مورد نظر نشان داده است. با این وجود، تخمین محل تونل با ترکیب نتایج به دست آمده از هر دو روش نتایج قابل اعتمادتری ارائه داد. همچنین در این تحقیق، بی‌هنجاری‌های ژئوالکتریکی حاصـل از یک بلوک مستطیلی با انجام شبیه‌سازی‌های متنوع مورد بررسی و تحلیل قرار گرفته است. بر اساس نتایج به دست آمده، بلوک مستطیلی مورد نظر در عمق‌های بالاتر از 140 متر، قابل آشکارسازی نیست. همچنین، هر اندازه مقاومت‌ویژه روباره کاهش و یا ضخامت روباره افزایش یابد، احتمال آشکارسازی تونل کمتر می‌شود. نتایج نشان می‌دهند که روش مقاومت‌ویژه الکتریکی از توانمندترین روش‌های ژئوفیزیکی در شناسایی تونل‌های زیرزمینی است.}, keywords_fa = {آشکارسازی,تونل‌های زیرزمینی,مقاومت‌ویژه الکتریکی,لرزه‌نگاری انکساری}, url = {https://jesphys.ut.ac.ir/article_58914.html}, eprint = {https://jesphys.ut.ac.ir/article_58914_363be8916612e77b892e1799bdcfa505.pdf} } @article { author = {Haji Aghajany, Saeid and Voosoghi, Behzad}, title = {Comparison of ERA-Interim global meteorological reanalysis model and MERIS data to reduce the tropospheric effects in InSAR displacement velocity fields}, journal = {Journal of the Earth and Space Physics}, volume = {42}, number = {3}, pages = {607-618}, year = {2016}, publisher = {Institute of Geophysics, University of Tehran}, issn = {2538-371X}, eissn = {2538-3906}, doi = {10.22059/jesphys.2016.58884}, abstract = {The main limitation of radar interferometry in measuring ground displacements is due to phase propagation delays caused by the troposphere. Tropospheric water vapor is a major limitation for high precision Interferometric Synthetic Aperture Radar (InSAR)applications. Using global meteorological reanalysis models and MERISdata are two methods that can be used for correcting the effects of the troposphere. The purpose of this study is comparison of these two methods. The MERIS instrument is located on the platform of ENVISAT satellite and measures reflected solar radiation. We adopt the use of data of this sensor for our study as MERIS provides a direct estimate of atmospheric water simultaneously with Synthetic Aperture Radar (SAR) acquisition. ERA-Interim is a global atmospheric model computed by the ECMWF based on a 4D-Var assimilation of global surface and satellite meteorological data. This reanalysis provides values of several meteorological parameters on a global 70 km grid from 1989 to the present day, at 0 am, 6 am, 12 pm and 6 pm Universal Time Coordinated) UTC (daily. The vertical stratification is described on 37 pressure levels, densely spaced at low elevation (interval of 25 hPa), with the highest level around 50 km (1 hPa). For each acquisition date, we select the ERA-Interim and MERIS outputs that are of the closest time spans to the SAR acquisition time. A Kriging interpolation in the horizontal dimensions and a spline interpolation along altitude is then applied to produce a map of the predicted delay. Total delay maps at epoch of acquisitions are then combined by pairs to produce differential delay maps corresponding to each interferogram. The use of the precise formulation of the single path delay and of the profiles of temperature, water vapor and dry air partial pressure is of importance to compute an accurate delay function.We used the two ENVISAT radar acquisitions of a region in the north west of Iran. We calculated the displacement velocity field and we corrected it (by ERA-Interim and MERIS data) and then compared the results with reported displacement velocities of GPS stations. The RMSE of two methods were 1.84 mm and 2.37 mm. The maximum difference between two methods is about 7.7 mm. This difference could be due to the presence of cloudy pixels in the MERIS data. The minimum difference between two methods is about 0.3 mm. The reason for this difference is negligible horizontal changes in the tropospheric indicators. The results show that cloudy weather and changes in the troposphere indicators, are the most important factors in the accuracy of the results of the two methods.  }, keywords = {Troposphere,InSAR,ERA-Interim,MERIS}, title_fa = {مقایسه داده های بازتحلیل ERA-Interim و حس گر MERIS در کاهش اثر لایه وردسپهر موجود در میدان های سرعت جابجایی تداخل سنجی راداری}, abstract_fa = {دو مورد از روش‌های متداول موجود برای شناسایی اختلاف فاز مربوط به وردسپهر، استفاده از داده های اتمسفری بازتحلیل شده و همچنین نقشه‌های پهنه بندی بخار آب و فشار بدست آمده از سنجنده MERIS است. هدف این مقاله مقایسه این دو روش با یکدیگر است. به منظور دستیابی به میزان اختلاف فاز از دو تصویر راداری مربوط به منطقه‌ای در شمال غرب ایران استفاده گردید. علاوه بر بدست آوردن میزان اختلاف فاز دو تصویر بدون تصحیحات از طریق روش‌های مذکور، میزان این اختلاف پس از اعمال تصحیحات وردسپهری از طریق داده های ERA-Interim و داده های MERIS نیز بدست آمد و نتایج بدست آمده با میزان جابجایی گزارش شده از ایستگاه های GPS موجود در منطقه مقایسه گردید. این نتایج حاکی از آن است که مقدار RMSE در حالت استفاده از داده‌های Interim ERA- برابر 84/1 میلیمتر و در حالت استفاده از داده‌های سنجنده MERIS برابر 37/2 میلیمتر می‌باشد. بیشینه اختلاف بین دو روش تصحیح برابر 7/7 میلیمتر می‌باشد که علت این اختلاف را می‌توان به وجود و عدم وجود پیکسل‌های ابری در آن ناحیه و در نتیجه متغیر بودن دقت تصاویرMERIS در این دو حالت نسبت داد. همچنین کمترین اختلاف بدست آمده برابر 3/0 میلیمتر است که ناشی از کم بودن تغییرات مسطحاتی عوامل جوی در ناحیه می‌باشد. نتایج نشان می‌دهد که در حالت وجود ابرهای کومه‌ای در منطقه استفاده از داده‌های ERA-Interim نتایج مطلوبتری به همراه خواهد داشت.}, keywords_fa = {بخارآب,تداخل‌سنجی راداری,وردسپهر,ERA-Interim,MERIS}, url = {https://jesphys.ut.ac.ir/article_58884.html}, eprint = {https://jesphys.ut.ac.ir/article_58884_e8086d2dfa6fee67e882dcded257d3f3.pdf} } @article { author = {عطایی, فرنوش and ایران نژاد, پرویز and مزرعه فراهانی, مجید and علیزاده چوبری, امید}, title = {Distribution of mineral dust on the global scale and its impacts on radiation fluxes as simulated by the WRF/Chem model}, journal = {Journal of the Earth and Space Physics}, volume = {42}, number = {3}, pages = {619-632}, year = {2016}, publisher = {Institute of Geophysics, University of Tehran}, issn = {2538-371X}, eissn = {2538-3906}, doi = {10.22059/jesphys.2016.55676}, abstract = {Mineral dust is produced from both natural and anthropogenic sources. Dust aerosols can be transported over long distances in the atmosphere. They reduce the incident shortwave radiation to the surface by absorbing and scattering the solar radiation; thereby leading to a cooling effect at the surface and lower tropospheric temperature. On the other hand, by absorption and re-emission of longwave radiation, they increase the net surface longwave radiation at the surface. This direct interaction of dust aerosols with shortwave and longwave radiation, known as the direct radiative impact, plays a key role in the radiation budget of the atmosphere. Although mineral dust is one of the most significant aerosols in the atmosphere, according to the Intergovernmental Panel on Climate Change (IPCC, 2007), uncertainty in its spatial distribution and radiative forcing, remains as a great challenge in climate studies. In the present study, the Weather Research and Forecasting with Chemistry (WRF-Chem) regional model is used to simulate distribution of mineral dust and its impacts on radiation fluxes on the global scale. The model was executed using 335 × 168 horizontal grid points with a horizontal spacing of 120 km, and 28 vertical levels for January and July 2011. The National Centers for Environmental Prediction (NCEP) Final Analysis (FNL) re-analysis data were used as meteorological initial conditions. The GOCART (Goddard Global Ozone Chemistry Aerosol Radiation and Transport) simple aerosol scheme was used for the simulation of dust emission and airborne dust distribution. Two experiments were conducted: the control simulation with no dust; and the interactive simulation for which dust aerosols feedback to the atmosphere. Differences between these two simulations indicate the perturbation of radiation by dust. Results indicate that the concentration of dust particles is generally much higher in the Northern Hemisphere than the Southern Hemisphere. The main sources of dust are located over the Sahara and Sahel, the Middle East, and East Asia, especially the Gobi Desert of China and Mongolia. The Eyre Basin in central Australia was identified as the most important source of dust in the Southern Hemisphere. Over the Sahara, dust emission was most intense in January, but substantially decreased in July. In contrast, in response to drier soils and higher wind speeds, sources of dust in the Middle East were more active in July than January. The Gobi Desert was also found to have much more dust activity in January than July, primarily due to stronger wind speeds during this month. On the global scale, monthly-averaged dust optical depth (DOD) was estimated to be 0.046 and 0.069 in January and July, respectively. Globally, perturbation of shortwave and longwave radiation by dust at the top of the atmosphere (TOA) was estimated to be -1.84 and 1.34 W m-2 in January, and -2.38 and 0.68 W m-2 in July, respectively. At the surface, it was estimated that perturbation of shortwave and longwave radiation to be -2.07 and 0.82 W m-2 in January, and -4.14 and 1.02 W m-2 in July, respectively. It was also found that perturbation of radiation is larger closer to the sources of dust. For instance, the perturbation of shortwave radiation exceeds -20 W m-2 over the Sahara. Globally, we identified that dust has a negative effect on the shortwave, but a positive effect on the longwave radiation at the surface. However, in snow covered regions (such as over the Tibetan Plateau, northern parts of the Scandinavia and the United States in January) deposition of dust on the surface increases the net shortwave radiation reaching the surface (due to reduction of surface albedo) and decreases net longwave radiation by increasing outgoing longwave radiation from the surface.  }, keywords = {Mineral Dust,Dust Optical Depth,Dust Emission,Perturbation of Radiation,WRF_Chem}, title_fa = {شبیه‌سازی توزیع ذرات گردوخاک در مقیاس جهانی و تاثیر آنها بر شارش‌های‌ تابشی با استفاده از مدل WRF/Chem}, abstract_fa = {در مطالعه حاضر توزیع ذرات گردوخاک و تاثیر مستقیم‌شان بر شارش‌های تابشی در مقیاس جهانی با استفاده از مدل عددی WRF-Chem، جفت شده با طرحواره گسیل گردوخاک GOCART، برای ماه‌های ژانویه و ژوئیه 2011 شبیه‌سازی شده است. دو شبیه‌سازی اجرا گردیده است: پیکربندی بدون ذرات گردوخاک و پیکربندی که شامل ذرات گردوخاک و تاثیر مستقیم‌شان بر تراز تابشی است. اختلاف نتایج دو آزمایش عمق نوری ذرات گردوخاک و پریشدگی تابش‌ طول موج کوتاه و بلند توسط ذرات گردوخاک را به‌دست می‌دهد. میانگین جهانی عمق نوری ذرات گردوخاک در μm 55/0 برای ماه‌های ژانویه و ژوئیه به‌ترتیب 046/0 و 069/0 برآورد شد که بیانگر بیشتر بودن ذرات گردوخاک معلق در جو در ماه ژوئیه نسبت به ژانویه می‌باشد. میانگین جهانی پریشیدگی تابش طول موج کوتاه توسط ذرات معلق گردوخاک و در شرایط آسمان صاف در سقف جوّ و سطح زمین در ژانویه به‌ترتیب Wm-284/1- و Wm-207/2- و در ژوئیه بهترتیب Wm-238/2- و Wm-214/4- محاسبه گردید. پریشیدگی تابش طول موج بلند مثبت، و در ژانویه Wm-234/1 در سقف جو ّو Wm-282/0 در سطح زمین و در ژوئیه Wm-286/0 در سقف جوّ و Wm-2 02/1 در سطح زمین برآورد شد. این مقادیر نشان می‌دهد که پریشیدگی تابش خالص توسط ذرات گردوخاک در سطح زمین و سقف جوّ منفی است، یعنی این ذرات باعث سرمایش سامانه زمین- جوّ می‌شوند. پریشیدگی تابشی منفی در سطح زمین ناشی از خاموشی (پراکنش و جذب) تابش خورشیدی ورودی توسط ذرات گردوخاک و پریشیدگی تابشی منفی در سقف جوّ عمدتا ناشی از افزایش طول موج کوتاه خروجی است.}, keywords_fa = {ذرات گردوخاک,عمق نوری گردوخاک,چشمه‌های گسیل گردوخاک,پریشیدگی تابشی,WRF_Chem}, url = {https://jesphys.ut.ac.ir/article_55676.html}, eprint = {https://jesphys.ut.ac.ir/article_55676_fa943c52f87cdaaac6903474ead55c4f.pdf} } @article { author = {Mostafazadeh, Raoof and Zabihi, Mohsen}, title = {Comparison of SPI and SPEI indices to meteorological drought assessment using R programming (Case study: Kurdistan Province)}, journal = {Journal of the Earth and Space Physics}, volume = {42}, number = {3}, pages = {633-643}, year = {2016}, publisher = {Institute of Geophysics, University of Tehran}, issn = {2538-371X}, eissn = {2538-3906}, doi = {10.22059/jesphys.2016.57881}, abstract = {Drought indices are commonly used to quantify and assess drought characteristics. The Standardized Precipitation Index (SPI), and recently introduced Standardized Precipitation-Evapotranspiration Index (SPEI) are considered as universal meteorological drought indices which allow comparisons of drought conditions across different climate regions. The SPI captures anomalies in precipitation, whereas the SPEI estimates anomalies in climatic water balance that incorporates temperature. The main aim of this study is to compare historical drought occurrence based on SPI and SPEI indices using R programming. The SPEI index is used because of multi-scalar nature of index and the advantage of identifying the multi-temporal nature of droughts. According to data availability, seven synoptic stations were selected for a drought analysis across Kurdistan Province. The two-parameter gamma distribution was used for calculating SPI across the study period (1995-2013) and stations within the study area. The potential evapotranspiration (PET) was computed using the Thornthwaite's equation, and then the SPEI is calculated at a monthly temporal resolution using SPEI package in R software. The SPI and SPEI values are calculated and then the statistical analysis along with significant level scatter plot was performed. The graphical plot of 3-month SPI and SPEI values were prepared to visualize the capabilities of used indices in determination of wet and dry spells over studied stations. The relationships between computed SPI and SPEI values were analyzed using correlation coefficient and p-value at each station. The results indicated that some differences in the pattern and sequence of wet and dry spells exist based on calculated indices. Also, the SPEI index identified the longer wet and dry spell conditions than SPI in almost all cases. The results of the comparative analysis indicated that the SPI and SPEI were varied between 0.19 (p<0.01), and 0.52 (P<0.01), which were statistically significant in all stations. A very low correlation between the SPI and SPEI was identified in Saghez station (correlation coefficient= 0.19), which seems to be due to evapotranspiration and moisture loss during spring/summer with the increasing temperatures that is accounted for by SPEI. The highest correlation coefficient was calculated between SPI and SPEI in the Sanadaj station (0.52%). Since, the SPEI accounts temperature in defining drought spells, therefore, it is advisable to use SPEI instead of SPI for drought assessment. According to the graphical interpretation of the results, there was large difference between the droughts depicted by the precipitation-based SPI and the temperature influenced SPEI. Also the SPEI captured the influence of temperature and depicted severe and longer duration droughts which provide support for better performance and reliability of the SPEI index. It should be noted that in terms of lack of data, evapotranspiration can be calculated by simple methods such as Thornthwaite, but considering detailed available data, the Hargreaves and Penman methods can be used to determine drought occurrences in the SPEI calculation. The calculation of drought indices should be simple and statistically reliable, in this regard, SPEI indicators in different climatic conditions and climate change issues has a significant advantage. Also, more climate variables are needed to calculate the SPEI index than the SPI index. Also, the calculated evapotranspiration value is sensitive to the used method and requires a longer data period with natural variabilities. Further research is recommended in other climatic regions which is needed for comparison of SPEI with other common drought indices to draw comprehensive conclusions.}, keywords = {Meteorological drought,Potential evapotranspiration,R software,Standardized Precipitation Index,Wet-Dry Spell Analysis}, title_fa = {تحلیل و مقایسه شاخص‌های ‏SPI‏ و ‏SPEI‏ در ارزیابی خشک سالی هواشناسی با استفاده از نرم‌افزار R (بررسی موردی: استان کردستان)}, abstract_fa = {شاخص SPEI خشکسالی بیلان آبی را با درنظر گرفتن دما و تبخیر و تعرق برآورد می‌نماید. هدف پژوهش حاضر، مقایسه الگوی وقوع خشکسالی براساس شاخص‌های SPI و SPEI به‌دلیل ماهیت چند مقیاسی بودن و قابلیت تحلیل خشکسالی در مقیاس‌های زمانی می‌باشد. در این راستا هفت ایستگاه سینوپتیک استان کردستان انتخاب و شاخص SPI بر اساس توزیع گامای دومتغیره محاسبه گردید. تبخیر و تعرق پتانسیل جهت استفاده در شاخص SPEI توسط معادله تورنوایت محاسبه و سپس مقادیر شاخص‌های SPI و SPEI در مقیاس زمانی 3 ماهه با استفاده از برنامه نویسی R به‌دست آمد. مقادیر SPI و SPEI در قالب گراف‌های توالی دوره‌ها رسم و ارتباط آن‌ها با تحلیل همبستگی مورد آزمون قرار گرفت. نتایج مقادیر شاخص‌های خشکسالی با ضریب همبستگی و سطح معنی‌داری به‌صورت ماتریس پراکندگی ابرنقاط ارائه شد. نتایج نشان داد که ضمن وجود تفاوت در الگوی وقوع دوره‌های ترسالی و خشکسالی SPI و SPEI، شاخص SPEI دوره‌های طولانی‌تر خشکی را در اکثر ایستگاه‌ها مشخص نموده است. بر اساس نتایج، ضرایب همبستگی بین 0/19 در ایستگاه سقز و 0/52 در ایستگاه سنندج و معنی‌دار از نظر آماری (در سطح 99 درصد)، متغیر هستند. به‌نظر می‌رسد همبستگی ضعیف (0/19) میان شاخص‌های SPI و SPEI در ایستگاه سقز ناشی از افزایش تبخیر و تعرق در فصول بهار و تابستان محاسبه شده توسط شاخص SPEI است. با توجه به تاثیر دما در محاسبه SPEI، تفاوت قابل توجهی بین مقادیر شاخص استاندارد مبتنی بر بارش و شاخص تبخیر و تعرق مبتنی بر دما، زمینه کارایی و صحت شاخص مذکور را فراهم می‌نماید.}, keywords_fa = {تبخیر و تعرق پتانسیل,دوره‌های خشک و مرطوب,خشکسالی هواشناسی,شاخص بارش استاندارد,نرم‌افزار R}, url = {https://jesphys.ut.ac.ir/article_57881.html}, eprint = {https://jesphys.ut.ac.ir/article_57881_2878d7e873427f99f68d02763f2f2d2a.pdf} } @article { author = {Khalilabadi, Mohammad Reza}, title = {A numerical study of internal tide generation due to interaction of barotropic tide with bottom topography in the Oman Gulf}, journal = {Journal of the Earth and Space Physics}, volume = {42}, number = {3}, pages = {645-656}, year = {2016}, publisher = {Institute of Geophysics, University of Tehran}, issn = {2538-371X}, eissn = {2538-3906}, doi = {10.22059/jesphys.2016.57903}, abstract = {A major part of tidal energy is usually dissipated by the interaction of tidal currents with bottom topography. Gulf of Oman is a marginal sea which has variable topography and its dominant tidal constituent is M2 semi-diurnal tide. In this paper, the interaction of barotropic tidal current with bottom topography is evaluated. This phenomenon causes the formation of internal tide. Internal tide is a large scale and baroclinic phenomena which causes long wave oscillations of water column. Whereas the M2 semi-diurnal constituent is dominant, therefore this constituent is the main force for formation of internal tide in the Gulf of Oman. In this paper, numerical modeling of internal tide due to M2 semi-diurnal constituent is  presented. This modeling is done using iTides model. This model is a software package which produces the internal wave field from the barotropic tide. The iTides package provides a graphical user interface (GUI) that combines all the theoretical elements necessary for producing a desired internal tide field given a set of system parameters. The model setup begins by specifying the pathway of the file to the topography, where the topography shape is specified by the horizontal coordinate assigned x, topography height h and the topographic change dh (both h and dh are functions of x). We use this dataset to define the maximum depth of the problem Ho as the maximum depth reached by the topography. Then, the density stratification must be given. The final step requires the user to specify the tide (forcing) frequency and the Coriolis frequency of the problem. After all parameters have been declared, the topography shape and stratification profile can be reviewed. This work has been done by implementation of iTides numerical model, which is a borotropic model, and forcing of an oscillating tidal current with semi-diurnal period. The model results show the formation of internal tide with a wavelength of order of O (10km) which reduces to O (1km) when reaches the shallow water. According to studied profiles of stability frequency, density stratification is quite stable in the Gulf of Oman and this Gulf is capable for formation of internal tide. Internal tide wavelength is of order of tens kilometers which reduces to a few kilometers when reaching the shallow zones. In the results, also the energy dissipation over the topography is visible. Most of internal tide energy is related to first modes. This phenomenon is mostly extended to deep zones, but for shallow zones internal tide energy is considerable between 1 to 3th internal tide modes. This fact may be due to iterated reflection of internal tidal beams from continental shelf and coastal shallow waters. The maximum energy Flux of primary modes in deep water (with depth of about 3000 meters) reaches to 20 kiloWatt per meter, whereas by decreasing the water depth, this amount of energy Flux reduces. The amount of internal tidal energy Flux in the Strait of Hormuz shallow water reaches bellow five kiloWatt per meter.    }, keywords = {topography,Barotropic,internal tide,Gulf of Oman,iTides model}, title_fa = {یک مطالعه عددی تولید امواج درونی در اثر برهمکنش مدهای فشارورد جزر و مدی با توپوگرافی کف در خلیج عمان و تنگة هرمز با استفاده از مدل iTides}, abstract_fa = {بخش عمده‌ای از انرژی جزرومدی از طریق اندرکنش جریانهای جزرومدی با توپوگرافی بستر مستهلک می‌شود. دریای عمان یک دریای حاشیه‌ای است که مولفه جزرومدی غالب در آن مولفه نیمه‌روزانه M2 است و دارای توپوگرافی متغیری است. لذا مطالعه جزرومد درونی و انرژی حاصل از آن اهمیت پیدا می‌کند. جزرومد درونی یک پدیده‌ی بزرگ مقیاس و باروتروپیک است که باعث نوسان ستون آب با طول موج بلند می‌شود. با توجه به اینکه در خلیج‌عمان، مولفه جزرومدی نیمه روزانه M2 مولفه غالب است، لذا عامل اصلی شکل‌گیری جزرومد درونی در خلیج عمان، این مولفه است. در این مقاله، اندر‌کنش یک جریان جزرومدی باروتروپیکی با توپوگرافی بستر ارزیابی می‌گردد. این کار به وسیله مدل عددی iTides که یک مدل باروتروپیکی است و اعمال یک جریان جزرومدی نوسانی با پریود نیمه روزانه انجام گرفته است. نتایج مدل‌سازی، شکل‌گیری جزرومد درونی با طول موجی از مرتبه 10 کیلومتر را نشان می‌دهد که وقتی به آب کم‌عمق می‌رسد طول موج آن به مرتبه کیلومتر کاهش می‌یابد.}, keywords_fa = {توپوگرافی,جزرومد درونی,باروتروپیک,خلیج عمان,مدل iTides}, url = {https://jesphys.ut.ac.ir/article_57903.html}, eprint = {https://jesphys.ut.ac.ir/article_57903_59dd2b9365d8616c488d17bb0d4a8e28.pdf} } @article { author = {miri, morteza and رضیئی, طیب and رحیمی, مجتبی}, title = {Evaluation and statistically comparison of TRMM and GPCC datasets with observed precipitation in Iran}, journal = {Journal of the Earth and Space Physics}, volume = {42}, number = {3}, pages = {657-672}, year = {2016}, publisher = {Institute of Geophysics, University of Tehran}, issn = {2538-371X}, eissn = {2538-3906}, doi = {10.22059/jesphys.2016.56102}, abstract = {The lack of reliable and updated precipitation datasets is the most important limitation in the study of many climatological and hydrological subjects, including climate change and temporal variability of precipitation in many data sparse areas around the globe. This is particularly valid for Iran where vast areas of central-eastern country that host the Iranian deserts, suffer from an inadequate network of rain-gage stations, required for climatological studies. The highlands of the mountainous regions of western and northern Iran have the same problem and limited representative stations are available for high elevation areas of these regions. One of solution to overcome this obstacle is to use available gridded precipitation datasets that have proved their representativeness for many different parts of the world. Among many available precipitation datasets are the Global Precipitation Climatology Center (GPCC) and theTropical Rainfall Measuring Mission (TRMM) that have been widely used in many researches, indicating their accurate estimation of precipitation values and intera-annual variation for the regions studied. The GPCC is a gage based dataset that is routinely creating through interpolation of worldwide precipitation stations combined with satellite records, whereas the TRMM is a purely remote sensed data developed by joint collaboration between NASA and the Japan Aerospace Exploration Agency (JAXA). The representativeness and performance of the GPCC and TRMM-3B43V7 precipitation datasets in estimating precipitation amounts at the locations of 46 Iranian synoptic stations distributed across the country is herein examined. Spatial resolutions of TRMM-3B43V7 and GPCC datasets used in this study are respectively 0.25 × 0.25 and 0.5 × 0.5 latitude and longitude. For each station, the closest grid point of each of the datasets to the station coordinates were chosen for statistically comparison analysis. To evaluate the performance of these datasets in comparison with the observed precipitation records at the considered locations we have used R squared, the Nash–Sutcliffe model efficiency coefficient, RMSE, Bias, B slope of the regression and the standardized RMSE indicators. The performances of the datasets were also graphically represented through scatter plots of the established regression between the observation and each of the two used datasets. The results of the statistical indicators were represented through plotting the indicators over the map of Iran to ease revealing spatial tendency of the indicators and explaining the possible geographical role in controlling the spatial variation of the indicators. The results revealed that both GPCC and TRMM-3B43V7 perform well in majority of the studied stations with strong correlation coefficients. However, it was found that the TRMM-3B43V7 underestimates precipitation in some stations located in the coastal areas of the Caspian Sea as well as in some stations along the Persian Gulf and the Oman seas, indicating that TRMM-3B43V7 is somewhat inefficient in adequately estimating precipitation in the coastal areas; which is very likely due to being unable to remove the effect of sea atmosphere interaction in stations nearby the seas. Contrarily, in some locations mostly situated in northwestern and northeastern mountainous areas of the country the TRMM-3B43V7 moderately over estimates the observed precipitation. Similarly, the GPCC well estimates precipitation in almost all stations with very high correlation coefficient and Nash–Sutcliffe model efficiency coefficient. Similar to TRMM-3B43V7, again it was found that the GPCC underestimates precipitation in most stations located along the coastal areas of the Caspian Sea. As for TRMM-3B43V7, the over-estimations of GPCC are mostly observed in northwestern Iran which is very likely due to not incorporating enough stations from high elevation areas of western Iran by the GPCC. On the whole, the results indicate that both datasets perform well in most locations of Iran and can be confidentially used in climatological and hydrological studies with or without the observation data. The results also indicate that the GPCC perform better in areas that share a denser network of stations with GPCC and vice versa. However, the very good results achieved with TRMM-3B43V7 that are completely independent from the observation indicates a promising future in having much improved remotely sensed precipitation records that well match the observed precipitation in very remote areas having no rain gages.}, keywords = {Precipitation,TRMM,GPCC,statistical indicators,Iran}, title_fa = {ارزیابی و مقایسه آماری داده‌های بارش TRMM و GPCC با داده‌های مشاهده‌ای در ایران}, abstract_fa = {پژوهش حاضر با هدف ارزیابی دقت داده‌های بارش سنجنده TRMM-3B43 و داده‌های شبکه‌بندی شده GPCC در برآورد بارش واقعی ایستگاه‌های همدیدی کشور به انجام رسیده است. برای این منظور داده‌های ماهانه بارش 46 ایستگاه همدیدی ایران با پراکنش مناسب در سطح کشور، داده‌های بارش سنجنده TRMM و داده‌های بارش GPCC برای دوره مشترک آماری 2010-1998 از تارنماهای مربوطه دریافت و استفاده شد. دقت مکانی داده‌های سنجنده TRMM و GPCCبه ترتیب 25/0×25/0 و 5/0×5/0 درجه جغرافیایی است. برای ارزیابی دقت این داده‌ها از آماره‌های ضریب تعیین(r2)، مجذور میانگین مربع خطا (Rmse)، شیب خط(Slope)، اریبی(Bias) و ضریب کارایی مدل(EF) استفاده شد. مقایسه‌های آماری انجام شده نشان داد اگرچه داده‌های TRMM در برخی مناطق مانند ایستگاههای سواحل خلیج فارس و شمال غرب ایران و بصورت موردی برای ایستگاه‌هایی مانند تهران بارش را بیشتر و یا کمتر از مقدار واقعی برآورد می‌کند، اما در مجموع برآورد بارش به وسیله TRMM در بیشتر ایستگاههای مورد مطالعه از دقت خوبی برخوردار است. ارزیابی داده‌های شبکه بندی شده GPCC نیز نتایج مشابه‌ای را بدست داد که بیانگر دقت مناسب داده‌های GPCC در سطح ایران است. بیشترین میزان ضریب همبستگی برای مناطق شمال شرق، غرب میانه و شمال غرب ایران بدست آمد که دلیل آن تراکم زیاد ایستگاههای باران سنجی در این مناطق می‌باشد که GPCC از آن برای تولید این داده‌ها بهره برده است. بررسی توزیع زمانی بارش ماهانه TRMM و GPCC در مقایسه با داده‌های مشاهده‌ای نیز نشان داد که هر دو این داده‌ها به خوبی روند تغیرات بارش ماهانه داده‌های مشاهده‌ای را شبیه سازی می‌کنند.}, keywords_fa = {بارش,GPCC,TRMM,آزمون‌های آماری,ایران}, url = {https://jesphys.ut.ac.ir/article_56102.html}, eprint = {https://jesphys.ut.ac.ir/article_56102_16cd5573b75a14aa62740b63c4925e0e.pdf} } @article { author = {Moosavi, Seyyed Shafi and فرجزاده, منوچهر and قویدل, یوسف and بیدختی, عباسعلی علی اکبر}, title = {The synoptic analysis of ozone mini-hole events over central Iran (Esfahan)}, journal = {Journal of the Earth and Space Physics}, volume = {42}, number = {3}, pages = {673-686}, year = {2016}, publisher = {Institute of Geophysics, University of Tehran}, issn = {2538-371X}, eissn = {2538-3906}, doi = {10.22059/jesphys.2016.58890}, abstract = {In this study, the daily Total Ozone Column (TOC) measured by the instruments of TOMS (2001) and OMI (2005-2011) satellites and Brewer ground station (2002-2004) is used to investigate the extreme ozone mini-holes over Esfahan. Based on previous reports on validation of the TOC data products, it is found that there is no problem with homogenization of data records, which was provided by the above measuring instruments. Firstly, it is shown that the TOC monthly mean and standard deviation over central Iran depend on the seasonal cycle with maximum values of 298 and 27 DU in winter and minimum values of 270 and 8 DU in summer, respectively. The difference between the maximum and minimum climatological monthly means is 53 DU. Regarding the absolute values of TOC, the maximum (minimum) amplitude is related to the winter season with 169 DU in Feb (summer with 39 DU in Aug). Due to the minus twice standard deviation of the monthly average which is known as the threshold chosen to identify the possible ozone mini-holes, 25 events are detected during the study period with maximum concentrations, of which 16 and 7 cases occurred in autumn and winter seasons, respectively. The most occurrences of ozone mini-hole are seen in 2005 and 2011 with 7 and 6 events, respectively. It is worthwhile to mention that the lowest levels of ozone in Arctic were also seen during the two mentioned years from 2001 to 2011. Nevertheless, no mini-holes were detected for three years 2003, 2004 and 2009. The range of ozone negative anomalies is confined from around 24% in winter (Jan) to 6% in summer (Aug). However, it was reported that ozone mini-holes in some regions have reduced the TOC up to 40% of climatology mean of mid and high latitudes over the northern hemisphere. It is found that during ozone mini-hole events, the Tropopause Height (TH) tends to move upwards (with a maximum of 5.5 km higher than monthly average on 10 March 2008) which in turn leads to decrease in the temperature and pressure of TH. Similar to its seasonal cycle, the low observed values of the tropopause temperature and pressure in summer is stronger than winter season. In general, the ranges of temperature (pressure) in the thermal tropopause during low ozone events becomes from -2.3°C (-27 hPa) in February 2006 to -15.5°C (-115 hPa) in March 2008. However, the mentioned above pattern almost explains the maximum events, the observed ozone mini-holes in January 2002 don not show similar anomalies in TH. It is more probably that low ozone events during the January of 2002 are more related to the meridional transport of air masses with climatology low ozone from the subtropical latitude which is poleward near the tropopause. Backward trajectory analysis also showed that the origins of poor ozone air masses in the spring/summer (autumn/winter) seasons are related to the eastern areas (western areas) of Iran. On 7 Jan 2002 at 16 km altitude (on 16 Oct 2011 at 22km altitude), the lower part of trajectory analysis, is more characterized by horizontal movement of poor ozone air mass from lower latitude (higher latitude). During the two extreme low ozone events over Esfahan which approximately correspond to the deepest events and eventful periods, two broad ridges are seen over coastal line of North-West Europe along with two deep troughs in the eastern-central Mediterranean Sea. The blocking ozone mini-holes over North-West Europe are related to the upward movement of geopotential height in the upper troposphere and lower stratosphere (UTLS) region which is in agreement with both the advection of poor ozone air from the sub-tropical (7 Jan 2002) and the higher latitudes (16 Oct 2011) toward the mid latitudes over central Iran.  }, keywords = {Central Iran,Total Ozone Column (TOC),Ozone mini-hole,tropopause height (TH)}, title_fa = {تحلیل همدیدی رخداد حفره‌های کوچک ازن در منطقه ایران مرکزی (اصفهان)}, abstract_fa = {این تحقیق با بهره‌گیری از داده‌های روزانه ازن پوشن‌سپهر روی اصفهان که از طریق سنجنده‌های زمینی بروئر و ماهواره TOMS و OMI از سال 2001 تا 2011 اندازه‌گیری شده، جهت شناسایی و مطالعه شدیدترین رخداد کاهش ازن (حفره کوچک ازن) بکار گرفته شد. با بکارگیری آستانه منفی دوبرابر انحراف معیار ماهانه، 25 رخداد حفره‌ کوچک ازن ‌در طول دوره مورد مطالعه نمایان شد که بیشترین تمرکز آنها با شانزده و هفت مورد به ترتیب در پاییز و زمستان اتفاق افتاده است. دامنه ناهنجاری منفی آنها نیز از عمق 24% در زمستان تا 6% در تابستان در نوسان بوده است. این بررسی نشان داد که در مواقع شکل‌گیری حفره‌های کوچک ازن، ارتفاع وردایست (TH) به تراز بالاتر جو هدایت شده و همزمان دما و فشار آن نسبت به میانگین درازمدت کاهش محسوسی می‌یابد. نقشه‌های سطوح زمین‌پتانسیل در ناحیه وردسپهر بالایی و پوشن‌سپهر پایینی (UTLS) مرتبط با دو حادثۀ کاهش شدید ازن روی فلات ایران، پشته عمیقی را روی شمالغربی اروپا که با فرودی در شرق دریای مدیترانه همراهی داشته نشان داده است. علاوه بر الگوی همدیدی مشاهده شده در ناحیه UTLS روی فلات مرکزی ایران که موجب تسهیل در وزش افقی هوای ازن کم از منشأ عرض‌های جنب‌حاره‌ (حادثه 7 ژانویه 2002) و عرض‌های بالا (حادثه 16 اکتبر 2011) شده است، می‌توان به عامل دینامیک دوم که با صعود محلی سطوح هم‌آنتروپی موجب تشدید کاهش ازن پوشن‌سپهر می‌شود همزمان برای شکل‌گیری چنین رخدادهایی مهم قلمداد نمود.}, keywords_fa = {ایران مرکزی,ستون کلی ازن,حفره کوچک ازن,ارتفاع وردایست}, url = {https://jesphys.ut.ac.ir/article_58890.html}, eprint = {https://jesphys.ut.ac.ir/article_58890_5313659df378ed237727e47977cc6a56.pdf} } @article { author = {Ghahreman, Nozar and babaeian, iman and tabatabaei, mohmadreza}, title = {Evaluation the post processed outputs of dynamic models in estimation potential evapotranspiration changes under RCP scenarios (Case Study: Mashhad plain)}, journal = {Journal of the Earth and Space Physics}, volume = {42}, number = {3}, pages = {687-696}, year = {2016}, publisher = {Institute of Geophysics, University of Tehran}, issn = {2538-371X}, eissn = {2538-3906}, doi = {10.22059/jesphys.2016.57730}, abstract = {As a direct consequence of warmer temperatures, the hydrologic cycle will undergo significant impact with accompanying changes in the rates of precipitation and evaporation. Climate change will cause changes in climate variable such as precipitation, temperature, sunshine hours, wind speed and etc. So as a result of climate variable change, the related variable such as potential evapotranspiration will change too. As the soft computing skills increased in recent decades, more number of climate models has been developed for weather and climate predictions which have significantly improved the quality and quantity of projections. This notable increase in number of climate models has enabled the scientists to estimate a wide range of main climate variables such as precipitation and temperature in fine temporal and spatial resolutions. Although the uncertainty in model outputs still remains a main challenge. Upon the release of new scenarios based on radiative forcing which are known as Representative Concentration Pathway scenarios (RCP scenarios), by Intergovernmental panel on climate change (IPCC) in fifth assessment report (AR5), a new set of 42 global climate models (GCMs) have been proposed for future climate projections. Apart from increased number of available models, three main sources of uncertainty including: measurement error, variability, and model structure, that have been explained and studied in AR5.The aim of the current study is to investigate of changes of potential evapotranspiration (ET) over Mashhad plain, Northeast of Iran in future period 2021-2070 under two RCP scenarios i.e. RCP4.5 and RCP8.5. The main synoptic station in the region is Mashhad Station located at 59◦ 38ﹶE, 36◦ 18ﹶN, with elevation of 990 m. above M.S.L. The required meteorological data including maximum and minimum temperature, sunshine hours,wind speed for period of 1991 to 2005 were obtained from Iran Meteorological Organization for ET calculation using FAO Penman-Monteith (hereafter, FAO-PM) equation. Besides, the downscaled historical data of potential evapotranspiration provided by Swedish Meteorological and Hydrological Institute (SMHI) have been retrieved for the baseline period of 1991-2005.Then these historical estimated data were compared with those estimated using FAO-PM equation. The historical ET values were post-processed using a statistical proposed method for more accuracy. By completion of this part, the accuracy of historical dataset provided by SMHI was confirmed and used for further comparisons. In the second section the ET variations for future period of 2021 to 2070 under two RCP scenarios of 4.5 and 8.5 was studied. The results showed better estimation of ET during warm months. Statistical comparisons using T-test revealed significant differences between historical and estimated values of ET in months of February, March and December. The correlation coefficient between post processed and observed data showed similar results as in T-test. Since the historical dataset of potential evapotranspiration provided by SMHI was acceptable, it was used for the analysis during future period (2021-2070) under RCP4.5 and RCP8.5 scenarios compared to baseline observed data. The result of this part showed that the highest increase of potential evapotranspiration would be for January by 15.4% and 16.4% under RCP4.5 and RCP8.5 scenarios respectively and October would experience lowest decrease by -12.5% and -10.0% decrease, respectively. In general ET increase will be more under RCP 8.5 scenario comparing to RCP 4.5}, keywords = {Potential evapotranspiration,climate change,RCP scenarios,Mashhad}, title_fa = {ارزیابی پس‌پردازش برون‌دادهای دینامیکی‌مدل‌های اقلیمی در برآورد تغییرات تبخیر ‌تعرق پتانسیل تحت سناریوهای واداشت تابشی (بررسی موردی: دشت مشهد)}, abstract_fa = {تبخیر تعرق پتانسیل پارامتری کلیدی در محاسبات نیاز آبی، مطالعات اقلیمی و مدلهای هیدرولوژیکی است. هدف از این مطالعه بررسی میزان تغییرات این کمیت بر اساس سناریوهای جدید واداشت تابشی (RCP) و مدلهای اقلیمی پس پردازش شده پروژه CORDEXدر 50 سال آتی (2021-2070)در دشت مشهد می باشد. به این منظور ، تبخیر تعرق پتانسیل با استفاده از معادله پنمن-مانتیث و داده های هواشناسی تاریخی دوره پایه (1991-2005) ایستگاه سینوپتیک مشهد برآورد شده و با خروجی دو مدل اقلیمی منتخب از مجموعه مدلهای CMIP5 مقایسه گردیده است. به منظور کم کردن خطاهای موجود در برآوردهای مدل، عمل پس‎پردازش بر روی داده‎های برآورد شده انجام شد، نتایج حاصل از پس پردازش نشان داد که مدل در ماه‎های گرم سال توانایی بهتری در برآورد این تبخیرتعرق پتانسیل نسبت به ماه‎های سرد سال دارد. همچنین مقادیر تبخیرتعرق پتانسیل در دشت مشهد طی سال‎های 2070 ـ 2021 تحت سناریوهای RCP4.5 وRCP8.5 ، تهیه و روند آن‎ها بررسی گردید.بر اساس نتایج حاصله، بیشترین افزایش مقدار تبخیر تعرق پتانسیل در ماه ژانویه تحت دو سناریوی RCP4.5 و RCP8.5 به ترتیب با 4/15% و 4/16% افزایش نسبت به دوره پایه و نیز بیشترین مقدار کاهش در ماه اکتبر با 5/12 % و 0/10 % کاهش نسبت به دوره پایه مشاهده شد.}, keywords_fa = {تبخیر تعرق پتانسیل,تغییر اقلیم,سناریوهای RCP,مشهد}, url = {https://jesphys.ut.ac.ir/article_57730.html}, eprint = {https://jesphys.ut.ac.ir/article_57730_a767647406e30dac3f0098321db46ed5.pdf} } @article { author = {masoompour samakosh, jafar and میری, مرتضی and رحیمی, مجتبی}, title = {Statistical-synoptic analysis of thunderstorm in the Southern Coast of Iran}, journal = {Journal of the Earth and Space Physics}, volume = {42}, number = {3}, pages = {697-708}, year = {2016}, publisher = {Institute of Geophysics, University of Tehran}, issn = {2538-371X}, eissn = {2538-3906}, doi = {10.22059/jesphys.2016.58912}, abstract = {Thunderstorms are regarded not only as a significant weather event but also as a key element in water and electricity cycles of the atmosphere. Generally, researchers consider the intense weather instability as a result of convection in lower levels of the atmosphere with high levels enough of humidity. Usually statistic instability, the humidity of lower levels of the atmosphere and lifting mechanisms near the ground are the main factors leading to convection. Moreover, the combination of three factors, instability, humidity and convergence in lower levels of the atmosphere plays an important role in increasing the possibility of thunderstorms. Accompanying phenomena like lightning, tornado, hail, winds, heavy precipitations (Changnon, 2001 and 1925) and hazardous atmospheric phenomena like turbulence, freezing, and wind sheering make considerable irrecoverable damages to natural and human environments, therefore recognizing the features of these phenomena have always been attracting the attention of researchers. The present study aims at recognizing statistic of thermodynamic, and synoptic features of thunderstorms of southern coasts of Iran. Referring to the archive of National Meteorological Organization, hourly data of atmospheric phenomena of 10 synoptic stations during a common twenty-year period (1995- 2014) were extracted. The data were processed in temporal scales of year, season and month. The data of upper atmosphere (radio-sound data), available in the website of Wayoming University, were applied to investigate the thermodynamic features of the occurred thunderstorms. The thermodynamic features include KI, SI, TT, LI, CAPE indices and skew- T chart in RAOB software. The days with the occurrence of thunderstorms had 5 mm or more at least in two stations that were selected to find synoptic patterns. As the samples were limited, synoptic patterns were done manually. The required maps were prepared using the data of geopotential height in 1000 - 500 hPa levels. Besides wind components u and v and sea level pressure, extracted from NCEP/ NCAR website, were mapped by GrADS software. Checking yearly frequency of thunderstorm occurrence in the southern coasts of Iran showed that the frequency of occurrence of storms in Booshehr station was more than its frequency in Hormozgan station. Moreover, the thunderstorms of Booshehr have a better chronological orderas it occurs during all common years. However, except for BandarAbbas, there is no chronological order for this phenomena. Therefore, it can be said that the occurrence of thunderstorm in the western coasts of the south of Iran has higher frequency than the central and the eastern regions, making it a potential area in this region for storm formation. The largest number of thunderstorm occurrence in seasonal scale is recorded for fall with 45% and winter with 43% respectively. Following seasonal conditions, the largest number of thunderstorm occurrence in monthly scale is recorded for cold months. In Hormozgan station, November, December, and January have more frequencies, while in Booshehr station January, February and March have more frequencies. Analyzing the applied instability indices showed that there was a slight extreme and great CAPE (more than 2500) in Bandarabas station. Besides the values of convection indices TT and KI for most of the thunderstorms suggested the possibility of convection occurrence. Instability indices LI and SI for the occurred thunderstorms reveals conditions of limited instability. Synoptic analysis shows not only the dominance of the westerly winds extending to the south of Saudi Arabia but also the location of divergent region and positive vorticity advection region in the studied region, making instability conditions raising air. The spread of the westerly winds is either due to formation of blocking system in the atmospheric middle level or their meridional blowing and cold air advection from Europe or the north of Asia to the east of Mediterranean. Statistical findings reveal that the occurrence of thunderstorms of western coasts of the Persian Gulf, have higher potential, and more frequency than the central and eastern regions. In seasonal scale, the largest number of occurrences is recorded for fall and winter respectively, while there is no substantial difference in different hours of day and night in hourly scale. As a matter of fact, they are possible to happen all the times. Synoptic analyses show that there is the dominance of two patterns of blocking systems and westerlies trough in the middle of atmosphere leading to instability and rising air in the studied region. The divergent region and positive vorticity advection region in the studied regions make instability condition and hence rising air. Based on the findings of thermodynamic indices, it can be said that convective activities and local instabilities are rarely responsible for thunderstorm occurrence in the region. Also for the occurrence of severe convective activities and relatively high instability, extreme instability and extreme severe instability is coincident with limited thunderstorm occurrences.}, keywords = {Instability indices,RAOB,Blocking System,Advection,Thunderstorm}, title_fa = {واکاوی آماری – همدیدی توفان‌های تندری سواحل جنوبی ایران}, abstract_fa = {هدف مطالعه پیشرو بررسی آماری، ترمودینامیکی و همدیدی توفان‌های تندری سواحل جنوبی کشور طی دوره 1995-2014 می‌باشد. بدین منظور از داده‌های ایستگاه‌های همدید استان‌های بوشهر و هرمزگان استفاده شد. جهت تحلیل همدیدی از داده‌های ارتفاع ژئوپتانسیلی ترازهای 500-1000 هکتوپاسکال، مؤلفۀ u، v و فشار سطح دریا استفاده شد. برای تحلیل شرایط ترمودینامیکی نیز داده‌های جو بالا توسط شاخص‌هایی نظیر CAPE، LI، TT بررسی شدند. نتایج نشان داد روند مکانی این پدیده از غرب به شرق کاهش می‌یابد و پتانسیل رخداد آن در استان بوشهر بیش از هرمزگان است. ازنظر زمانی، بیشترین رخداد این پدیده در پاییز با 45 و زمستان با 43 درصد ثبت‌شده است. در مقیاس ساعتی تفاوت قابل‌توجهی در ساعت مختلف مشاهده نمی‌شود و امکان رخداد آن در تمامی ساعت‌ها به‌ویژه ساعات صبح محلی وجود دارد. نتایج همدیدی حاکی از نفوذ بادهای غربی تا جنوب عربستان و قرارگیری ناحیه واگرایی و منطقه فرارفت تاوایی مثبت روی منطقه مطالعاتی است که شرایط ناپایداری و صعود هوا را ایجاد می‌کند. این گسترش بادهای غربی یا ناشی از شکل‌گیری سامانه‌های بندالی در تراز میانی جو و یا ناشی از وزش نصف‌النهاری آن‌ها و فرارفت هوای سرد از اروپا یا شمال آسیا روی شرق مدیترانه است. مقادیر حاصل از شاخص‌های ترمودینامیکی نشان داد که در این منطقه، همرفت و ناپایداری‌های ایجادشده تحت تأثیر شرایط محلی، عامل اصلی رخداد این پدیده نمی‌باشند و در صورت فراهم شدن شرایط مناسب در ترازهای بالایی جو به‌عنوان عامل تشدید کنند این پدیده محسوب می‌شوند.}, keywords_fa = {شاخص‌های ناپایداری,RAOB,سامانه بندالی,فرارفت,توفان تندری}, url = {https://jesphys.ut.ac.ir/article_58912.html}, eprint = {https://jesphys.ut.ac.ir/article_58912_d0aedb57e577199d90cd02a2045c17f3.pdf} }