موسسه ژئوفیزیک دانشگاه تهرانفیزیک زمین و فضا2538-371X35120090421Evaluation of Iran Strong Motion Network (ISMN) recording and its effects on the improvement of earthquake location in Zanjan and its neighboring regionsارزیابی کیفیت ثبت دادههای شتابنگاری سازمان تحقیقات ساختمان و مسکن و تأثیر آن در بهبود مکانیابی زلزلهها در محدوده استان زنجان1167996410.22059/jesphys.2009.79964FAعبدالرضاقدسدانشیار، گروه علوم زمین، دانشگاه علوم پایه زنجان، ایرانروحالهعسکریمربی، گروه علوم زمین، دانشگاه علوم پایه، زنجان ایرانJournal Article20210221By merging the data from the University of Tehran’s Iran Seismic Telemetry Network (ISTN) with those from Iran Strong Motion Network (ISMN), we investigate the improvement of earthquake location accuracy for the events which happened in the Zanjan province and its neighboring regions (the region bounded by 35º-37.5º N latitudes and 46.5º-50.5º E longitudes and for the period 1996-2006). Due to insufficient distribution of seismic stations in the study region, most of the events have poor
azimuthal coverage and location accuracy. The events in the Zanjan province are largely recorded by Tehran and Tabriz seismic sub-networks of ISTN in the east and west of the province, respectively. ISMN has a very good coverage within Zanjan and its neighboring provinces, and thus has a promising potential to improve the location accuracy of the events within the Zanjan province.
In this study, we assess the improvement of location accuracy by merging the ISMN data with the catalog of relocated events for Zanjan and its neighboring provinces (Askari and Ghods, 2007). The catalog is consists of 304 events for the period 1996-2006, with local magnitudes larger than 3.1, and RMS of less than 0.7 s. The catalog complete for magnitudes larger than 3.5. ISMN is an offline network and does not use GPS for precise timing of its waveform records. We first associate ISMN waveforms based on their time tag but later check if the phase readings are consistent with the phase readings from ISTN. Due to the lack of precise timing of ISMN waveforms, only relative phase arrival of Sg-Pg could be participated in the location procedure. We were able to associate 403 ISMN records with 76 events of the catalog.
We find that the ISMN data cannot significantly improve earthquake location in the Zanjan region. Out of 76 events with the ISMN data, azimuthal coverage and epicentral accuracy of 7 events could be improved significantly. This is primarily related to the fact that most of the regional faults are more or less aligned with the direction of the spread of the weak-motion seismic networks in east and west of the province. We found that ISMN could detect all events with magnitude above 4 within the study region. According to our results, the ISMN stations have been maintained properly, and we could not detect significant data loss. However, we could detect several problems in archiving the ISMN data. We have also found that the ISMN instruments do not have enough resolution for accurate recording of seismic amplitudes. This implies that picking of the first arrivals on ISMN waveforms may have errors in the range of 0.2-0.4 seconds.در این تحقیق هدف این است که با تلفیق دادههای شبکه شتابنگاری با دادههای شبکههای لرزهنگاری کشوری وابسته به مؤسسه ژئوفیزیک دانشگاه تهران، امکان بهبود مکانیابی زلزلههای استان زنجان و نواحی مجاور (چهارگوش محدود به عرضهای جغرافیایی 35 تا 5ر37 درجه شمالی و طولهای جفرافیایی 5ر46 تا 5ر50 درجه شرقی) مورد ارزیابی قرار گیرد. نتایج این تحقیق روشن ساخت که دادههای شبکه شتابنگاری در منطقه بررسی دارای پتانسیل نسبتاً ضعیفی برای بهبود مکانیابی زلزلهها است؛ زیرا اکثر گسلهای منطقه مورد بررسی در راستای بیشترین پوشش آزیموتی شبکههای لرزهنگاری قرار گرفتهاند و اضافه کردن دادههای شتابنگاری، تأثیر اندکی در بهبود پوشش آزیموتی زلزلهها دارد. با توجه به نتایج این تحقیق، عملکرد شبکه در نگهداری دستگاههای شبکه و استخراج و نگهداری دادهها در سطح قابل قبولی است. هر چند که اشکالات مهمی در نحوة عرضه دادههای رقومی، قابل مشاهده است. همچنین با بررسی شکل موجهای رکوردهای شتابنگاری مشاهده شد که دامنه امواج ثبت شده دارای دقت تفکیک کافی نیست و نمیتوان زمان رسیدهای امواج لرزهای را با دقت بیشتر از 2ر0 ثانیه تعیین کرد.https://jesphys.ut.ac.ir/article_79964_3a6ca4e7acf595a62c96b80ffe537f3e.pdfموسسه ژئوفیزیک دانشگاه تهرانفیزیک زمین و فضا2538-371X35120090421Paleostress analysis around the Lar Dam (Central Alborz), to recognize the structures involved in water escapeتحلیل تنش دیرین در اطراف سد لار (البرز مرکزی)؛ به منظور شناخت ساختارهای مؤثر در فرار آب17307996510.22059/jesphys.2009.79965FAصفیهامیدیاندانش آموختة کارشناسی ارشد پترولوژی، دانشکدة زمینشناسی، پردیس علوم، دانشگاه تهران، ایرانمحسنالیاسیاستادیار دانشکدة زمینشناسی، پردیس علوم، دانشگاه تهران، ایرانجمشیدحسنزادهدانشیار دانشکدة زمینشناسی، پردیس علوم، دانشگاه تهران، ایرانمژگانزارعینژادرئیس گروه اطلاعات مرجع (GIS)، سازمان زمینشناسی و اکتشافات معدنی کشور، تهران، ایرانJournal Article20210221Lar Dam is located 85 km to the north-east of Tehran. It supplies a fraction of agricultural water in Mazandaran and the drinking water in Tehran. This dam was built in 1980 and 2 water escape ways around it have been explored; one of them is under the right shoulder
of the dam and the other is located north of the dam toward the Haraz road (Ab-e-Ask region). For this reason, new researches are designed for shifting it to west of the same zone (Gozal Darreh). But our investigations show that structural factor (Faulting) is effective for water escape and this agent has expanded to the whole of this region and it has covered new places as well.
Geologically, the studied area is located on calcic formations related to the second and third geological periods (Jurassic and Cretaceous). These formations have a W-NW trend parallel to the general trend in the Central Alborz. In the studied region, folds, thrusts and reverse faults have the same trend is due to continuous pressure from the Arabian plate on the Iranian plate.
Originally, the aim of our research was to obtain the direction of effective stresses that are responsible for forming and expanding the structural factor for water escape. So we considered some structures such as faults, joints and fractures which are useful for reconstructing the tectonic events (the diverse directions of stress with respect to relative time) in this region. "Inversion Method" is the base of tectonic software and is designed with respect to some of these structures. The important point which is also the target of using the slickenside is to exploit the data and return step by step to reach the initial stress conditions.
In order to obtain the direction of the main stress axis in the time of affecting stress using the inversion analysis, a considerable amount of structural data was compiled. All slickenside was categorized into 12 groups, according to Yamaji software (2005) that is known to MIM (Multiple Inversion Method). Four final parameters form software are which are the results of solving reduced stress tensor. , the shape of stress field, is a quantitative parameter. So the examination of the trend of changes in the shape of stress field will be reconstructible.
Field work resulting in the structural data and the computer analysis performed afterwards shows intensive changes in the rate of stress field shape all around the dam. The path of these changes of the rate of the stress field shape shows the existence of a new fault beneath the dam which is probably a natural channel for water escape. These abrupt changes are from a prolate stress field shape (Φ=5/0-1) to an oblate stress field shape (Φ=0-5/0) with a linear WNW trend. This path is parallel to the linear trend of the pitch of the sink holes located at the basement of the dam. The first channel of the water escape is beneath of the right shoulder. Based upon detailed studies, the proposed position for constructing a new dam, at the western end of the Lar Dam (Gozal Darreh), is still located on the continuation of the new fault we discovered. In conclusion water would escape through this fault anyway.محدوده مورد تحقیق از لحاظ موقعیت زمینشناسی در سازندهای اکثراً آهکی دوران دوم (ژوراسیک و کرتاسه) با روند تقریباً<br /> W-NW در جنوب غرب آتشفشان دماوند قرار دارد. در گسترة طرح همسو با روند کلی البرز مرکزی، چینها و گسلهای رانده و معکوس با راستای W-NW (ناشی از جهت فشارش صفحة عربی به صفحة ایران) عامل شکلگیری ساختارهای غالب در منطقه است. بررسیهای میدانی، برداشت اطلاعات ساختاری و تحلیل رایانهای نتایج حاصل از آن، تغییرات شدید میزان گرادیان تنش در اطراف دریاچه را مشخص میکند. مسیر این تغییرات در جهت تکیهگاه راست، همراه با اطلاعات صحرایی به روند گسلی اشاره دارد که عامل احتمالی فرار آب از دریاچه است، به نحوی که حدفاصل این تغییرات ناگهانی از شکل میدان تنش دوکی شکل (1- 5ر0) به کلوچهای شکل (5ر0-0) به صورت یک مسیر خطی بارز با روند W-NW است. این مسیر گسلی با روند قرارگیری خطی حفرههای فروکش و اولین مکان فرار آب از زیر تکیهگاه راست همراستا است. این بررسی به شرایط مشابه حاکم بر نواحی اطراف سد، از جمله محل احداث سد جدید (گزلدره) نیز اشاره دارد.https://jesphys.ut.ac.ir/article_79965_7bb712c4006f3fa2b371c5b548082edd.pdfموسسه ژئوفیزیک دانشگاه تهرانفیزیک زمین و فضا2538-371X35120090421Determination of Uranium anomalies in Barandagh region by using airborne radiometry dataتعیین ناهنجاریهای اورانیم در منطقه برندق با استفاده از دادههای رادیومتری هوابرد31447996910.22059/jesphys.2009.79969FAاحمدرضالکزاییدانشآموخته کارشناسی ارشد ژئوفیزیک، مؤسسة ژئوفیزیک دانشگاه تهران، ایرانمجیدنبیبیدهندیدانشیار، گروه فیزیک زمین، مؤسسة ژئوفیزیک دانشگاه تهران، ایران0000-0002-9555-8327افشارضیاءظریفیعضو هیأت علمی دانشگاه آزاد اسلامی واحد لاهیجان، ایرانفرخشادیگانیکارشناس ارشد، سازمان انرژی اتمی ایرانمحمدکاظمحفیظیدانشیار، گروه فیزیک زمین، مؤسسة ژئوفیزیک دانشگاه تهران، ایرانJournal Article20210221Some new statistical techniques are gaining favor and momentum in the separation of background concentration values from anomalous values in determining the economics of extraction of Uranium deposits. In exploration of minerals and feasibility of
exploration procedures, old conventional methods are replaced by new ones which have roots in natural distribution patterns. One of the methods is usage of fractal geometry in the separation of various statistical populations used in these studies such as different background values, threshold limits and anomalous values.
In this paper, in the first step, separation of anomaly values has been performed by means of classical statistics. Then the tables of frequency distribution of Uranium, Thorium and Potassium have been classified, and the frequency distribution histograms have been plotted. The statistical parameters of these three elements have then been estimated. Then separation of anomaly values has been performed based on dispersion around the average. In the second step, separation of anomaly values has been performed by using fractal method based on concentration-area curves.
In this work a comparison of classical statistical method has been made with fractal techniques of separation and grouping of the various values. The data used in this study were the airborne acquired geophysical data of the area which has been based on gamma ray emission of radioactive nuclides present in natural earth’s environment. At the first stage of the study the statistical parameters such as mean, mode, median dispersion of the mean, standard deviation, skewness, kurtosis of the data were plotted and the relative values were calculated. At the second stage, using fractal geometry techniques, concentration-area mathematical model of fractal curves were drawn after the interpolation of X, Y and Z digitized data had been constructed. On the basis of the concentration-area model the fractal dimensions were calculated and separation of various statistical populations was made on the basis of tangent values drawn to the fractal curves. The trends of variations in various statistical populations representing the interpreted concentration values were made using the above two procedures, and the advantages and disadvantages of the methods are described. Finally, based on both classical statistics and fractal methods, anomaly maps are plotted in which the anomaly values are separated from background values for all three radioactive elements of Uranium, Thorium and Potassium.در اولین مراحل اکتشاف اورانیم، مهمترین قسمت کار اکتشافی، استفاده از دادههای رادیومتری هوابرد در تعیین بیهنجاریها است. در این مقاله ابتدا با روش آمار کلاسیک و با استفاده از محاسبه پارامترهای آماری روی دادههای برداشت شده ژئوفیزیک هوایی در منطقه برندق، جدایش جوامع بیهنجاری صورت گرفته است. سپس جدولهای توزیع فراوانی عناصر اورانیم، توریم و پتاسیم و هیستوگرامهای توزیع فراوانی این عناصر ترسیم شده است. پارامترهای آماری این عناصر محاسبه شده و در نهایت جدایش جوامع بیهنجاری براساس پراکندگی حول میانگین صورتگرفته است. در روش دوم بر اساس هندسه فراکتالی و با استفاده از نمودارهای تمام لگاریتمی عیار- مساحتِ بهدست آمده از دادههای رقومی و نقشههای همشدت رادیومتری، جدایش پلهای محیطهای متفاوت (زمینه، حدّ آستانهای، بیهنجاری) صورت گرفته است و در آخر نقشههای مربوط به مناطق بیهنجاری و معرفی اندیسهای معدنیِ قابل بررسی عناصر اورانیم، توریم و پتاسیم برای ادامه کار اکتشافی با استفاده از هر دو روش آمار کلاسیک و فراکتالی مورد بررسی قرار میگیرد.https://jesphys.ut.ac.ir/article_79969_d5307b73a823b594214ae84b7a415bea.pdfموسسه ژئوفیزیک دانشگاه تهرانفیزیک زمین و فضا2538-371X35120090421Short period fluctuations of seismicity around Tehran inferred from "a" and "b" valuesبررسی تغییرات کوتاه دوره لرزهخیزی گستره تهران با استفاده از پارامترهای a وb45577997010.22059/jesphys.2009.79970FAمحمداشتری جعفریمربی، گروه فیزیک زمین، مؤسسة ژئوفیزیک دانشگاه تهران، ایرانJournal Article20210221Earthquake size distribution follows a power law whose slope is known as the b-value and its constant is named as the a-value. The b-value fluctuations have been theoretically studied in laboratories and practically investigated in several seismotectonics zones e.g. volcanic areas, continental rifts and mines which also present different stress regimes. The b-value can explain the relative density of large and small events which has found many applications in seismic hazard studies, spatio-temporal prediction and earthquake physics. On the other hand the a-value is concerned with regional seismicity level, so studying these parameters can be of great help in an area just like Tehran where there is a high concentration of people and social-economical activities. To begin this study we extracted the events from the Tehran Digital Seismic Network database. Processing followed by removing time-dependent quakes under the examination of the Poissionian assumption and later by computing the magnitude of completeness using the goodness of fit method. Then a-value and b-value changes were mapped in time and space. The b-value temporal changes are not significant during the period of data which may be under the control of local effects but reduction by depth exists. Both values show a change of around 51.5E. Meanwhile the b-value map shows a reduction toward regions with high density of thrust and strike-slip faults.پراکندگی بزرگی زلزلهها از رابطه توانی پیروی میکند. شیب این رابطه مقدار-b و ثابت آن مقدار-a نامیده میشوند. تغییرهای مقدار-b بهصورت نظری در آزمایشگاهها و بهشکل عملی در زونهای گوناگون لرزهزمینساختی مانند ناحیههای آتشفشانی، ریفتهای قارهایی و معدنها که بیانگر رژیمهای متفاوت تنش نیز هستند بررسی شده است. مقدار-b نمایشگر پراکندگی نسبی رویداد زمینلرزههای کوچک و بزرگ نسبت به یکدیگر است و کاربردهای فراوانی در زمینههایی مانند خطر لرزهایی، پیشبینی زمانی- مکانی و فیزیک زمینلرزهها یافته است. مقدار-a نیز با سطح لرزهخیزی ناحیهایی در ارتباط است. پس بررسی این دو پارامتر در گستره تهران که محل تمرکز بخش بزرگی از جمعیت و کوششهای اقتصادی و اجتماعی کشور است میتواند به شناخت سیمای لرزهزمینساختی آن کمک فراوانی کند. برای بررسی نوسانهای مقدار-b از دادههای شبکه لرزهنگاری رقمی تهران سود جستهایم. در این راستا پس از حذف رویداد های وابسته به زمان این شبکه با فرض پواسونی بودن آنها بزرگی کاملی و تغییرهای آن به روش نیکویی برازش محاسبه شد. سپس مقدار-b در حوزه زمان- مکان مورد توجه قرار گرفت. در دوره فعالیت این شبکه نوسانهای مقدار-b چشمگیر نبوده است و میتوان آنرا به تغییرهای محلی نسبت داد. در حوزه مکان مانند بعضی گسترههای قارهای این پارامتر با ژرفا کاهش یافته است. نقشه نوسانهای سطحی آن نیز بیانگر کاهش در محل تجمع بیشتر گسلهای رورانده و امتدادلغز <br /> است. مقدار-a محاسبه شده نیز گویای میزان لرزهخیزی منطقهایی است. هر دو مقدار نوسانهایی را حول طول 5ر51 شرقی نشان میدهند.https://jesphys.ut.ac.ir/article_79970_a3b7fdb160862ebc978758da767142e1.pdfموسسه ژئوفیزیک دانشگاه تهرانفیزیک زمین و فضا2538-371X35120090421Modeling the co-seismic deformation field of a fault and determining the sensitivity of the geometrical and physical parameters of the model to this deformation fieldمدلسازی میدان جابهجایی هملرزه یک گسل و تعیین حساسیت پارامترهای هندسی و فیزیکی مدل به میدان جابهجایی آن59737997110.22059/jesphys.2009.79971FAسینانوریدانشآموخته کارشناسی ارشد نقشهبرداری، دانشکده مهندسی نقشهبرداری، دانشگاه صنعتی خواجه نصیرالدین طوسی، تهران، ایرانبهزادوثوقیاستادیار، دانشکده مهندسی نقشهبرداری، دانشگاه صنعتی خواجه نصیرالدین طوسی، تهران، ایران0000-0002-5667-0447امیر مسعودابوالقاسماستادیار، دانشکده مهندسی نقشهبرداری، دانشگاه صنعتی خواجه نصیرالدین طوسی، تهران، ایرانJournal Article20210221The study of faulting and its resultant earth surface deformations is an essential research field in Iran. The major part of this large country is located in very active seismic zones. Any study in this field can enable us to mitigate the risk of earthquake hazards.
3D modeling of displacement and surface deformation caused by earthquake faulting based on a homogeneous, isotropic, elastic half space model is the main aim of this paper. The paper focuses on the modeling of a 3D co-seismic deformation field caused by stress accumulation and its release along seismogenic faults based on a homogeneous elastic half-space model. The most commonly used analytic models of fault deformations have been based on the dislocation solutions of Okada (1985, 1992). This dislocation model is used to investigate surface deformations which are generated by strike-slip and dip-slip faulting. A method of sensitivity analysis is applied to determine sensitivity of the model and its resultant displacement field with respect to the change of parameters of the model.
During the last decades, powerful new models have been developed and deployed with encouraging results for improving knowledge of fault system behavior and its consequent earthquake hazards. Fault displacement models based on elastic dislocation theory have been used to calculate displacements and strains due to co-seismic slip events. In the elastic dislocation theory, faults are considered as displacement discontinuities or dislocations in an otherwise continuous elastic medium. In this approach, faults are represented as surfaces across which there is a discontinuity in the elastic displacement field.
The elastic dislocation theory is conceptually valid for modeling co-seismic deformations. The elastic dislocation formulation of Okada is used in our models, which expresses the displacement field U(x, y, z) at any given point as a function of fault parameters (slip, dip, strike, length, and width) and the elastic constants within the continuum, for rectangular fault panels with horizontal upper and lower edges. The Okada formulation is mathematically robust and tractable, and these attributes make it suitable for rapid, iterative, forward numerical modeling.
In the first section of the paper, the relationship between surface deformation and dislocation theory will be summarized using representation formula. The dislocation theory can be described as that part of the theory of elasticity dealing with surfaces across which the displacement field is discontinuous, the suggestion seems reasonable. As commonly done in mathematical physics, it is necessary for simplicity to make some assumptions. Here the curvature of the earth, its gravity, temperature, magnetism and non-homogeneity are neglected and a semi-infinite medium which is homogeneous and isotropic is considered. For this modeling, the fault parameters that must be considered, are dislocation amount, length, width, depth and dip angle for fault plane.
This model can calculate displacements at every depth and the free surface specially. Here, the first displacement field is calculated for a simulated fault and sensitivity analysis is carried out for model parameters. The dislocation model provides us with surface deformation fields generated by strike-slip and dip-slip faulting and the vector maps of horizontal and vertical displacement fields can be represented.
In the next step, the sensitivity analysis is done to determine the sensitivity of the model and its deformation behavior with respect to any fault parameters. Then the results of the analysis for both cases of strike and dip-slip faults are compared. The analysis shows that the model has maximum sensitivity to dislocation parameter and minimum sensitivity to lame coefficients.
The numerical results of the analysis show that when the amount of dislocation increases the range and area of the surface deformation are greater. The horizontal displacements are more sensitive to the change of the dislocation amount in comparison with the vertical displacements. The results of the analysis result are summarized in the following table.
Result of sensitivity analysis
No.
Parameters
1
Dislocation (U)
2
Depth of fault (c)
3
Dip Angle (δ)
4
Width of fault (w)
5
Length of fault (L)
6
Lame coefficient (λ , µ)
The model can be applied for simulating co-seismic deformation fields of faulting to prepare a hazard map for the investigated fault in case of any consequent earthquake due to fault motions and for use in any further planning. This knowledge will translate into tangible societal benefits by providing the basis for more effective hazard assessments and mitigation efforts.شناخت حرکات گسل و نتایج حاصل از این حرکات نظیر تغییر شکل سطحی، در کشور ما که دارای گسلهای فعال و مناطق زلزلهخیزی است ضرورت و اهمیت فراوانی دارد. این تحقیق روی مدلسازی سهبعدی تغییر شکلهای هملرزه ایجاد شده در اثر حرکت گسل در یک نیمفضای کشسان هموژن است. این مدلسازی اغلب براساس نظریة جداشدگی صورت میپذیرد و بیشتر مدلهای تحلیلی تغییر شکل گسل نیز بر پایه همین نظریه که محقق ژاپنی اکادا در طول دهه هشتاد آن را بسط داد و در 1992 ارائه کرد، پایهریزی شدهاند. بر پایه این مدل، جابهجایی حاصل از وقوع حرکت در یک گسل امتدادلغز و یا شیبلغز (زلزله) را میتوانیم با استفاده از دادههای در دسترس مربوط به مشخصات گسل مورد نظر بهدست آوریم و نتایج حاصل را تحت عنوان تغییر شکل حاصل از زلزله ارائه کنیم.
این مدل قابلیت بهدست آوردن تغییر شکل و جابهجایی در هر عمقی نسبت به سطح آزاد (سطح زمین) را دارد که در اینجا هدف اصلی بهدست آوردن این مقادیر برای سطح زمین است. یکی از نتایج حاصل از این مدلسازی را میتوان تعیین میزان حساسیت مدل و تغییر شکلهای مسطحاتی به پارامترهای متفاوت یک گسل بیان و این ویژگیها را براساس میزان تأثیر آنها در خروجی مدل دستهبندی کرد. براساس نتایج حاصل از این آنالیز مشاهده میشود که بیشترین میزان حساسیت به تغییر در کمیت جداشدگی و کمترین آنها در ضرایب لامه هستند که جابهجایی به تغییرات اعمال شده در این ضرایب کمترین میزان حساسیت را نشان میدهد. با تعیین میزان حساسیت و مقدار تغییری که در جابهجایی در اثر تغییر در آن پارامتر ایجاد میشود میتوان به تعیین پارامترهای یک گسل نا شناخته در منطقهای که عامل ایجاد زلزله شده است پرداخت. این مدل قابلیت بهکارگیری برای هر نوع گسل و با هر خصوصیات هندسی و در هر منطقهای را دارد و میزان جابهجاییهای هملرزهای برای نقاط واقع در حوزه حرکت گسل یا یک زلزله میتوان بر پایه آن بهدست آورد.https://jesphys.ut.ac.ir/article_79971_b73a9d0f917d7272981046a65c79147c.pdfموسسه ژئوفیزیک دانشگاه تهرانفیزیک زمین و فضا2538-371X35120090421Focal mechanism analysis using synthetic seismogramsبررسی سازوکار زمینلرزهها با استفاده از لرزهنگاشتهای مصنوعی75887997210.22059/jesphys.2009.79972FAمحمدرضاحاتمیاستادیار، گروه فیزیک زمین، مؤسسة ژئوفیزیک دانشگاه تهران، ایرانظاهرحسینشمالیاستادیار، گروه فیزیک زمین، مؤسسة ژئوفیزیک دانشگاه تهران، ایرانغلامجوان دلوییاستادیار، پژوهشگاه بینالمللی زلزلهشناسی و مهندسی زلزله، تهران، ایرانJournal Article20210221According to the representation theorem, the elastic displacement due to a point source is given by the following equation:
(1)
In this equation, is the component of the displacement, is the source time function which indicates how the energy is released during the earthquake process, , is the Green's function that describes the propagation path effects between the source located at and the station at . are the moment tensor components . In this equation all components of the moment tensor are assumed to have the same time functionality. Computation of the Green's function is the most important step for producing synthetic seismograms.
1. Simulation of an Earthquake and its Linear Inversion for a Completely Shear Source (Pure Double-Couple)
1-1 Earthquake Simulation: For the simulation purpose, a source with definite rake, dip and strike was assumed. The synthetic seismograms were then calculated using the wavenumber integration method based on a 12km source depth.
1-2 Linear Inversion of the Seismograms for Determination of the Earthquake source parameters under Pure Double-Couple Conditions: The moment tensor is an overall indication of the earthquake source i.e. the volume variations and a variety of shear source in different directions are included in the moment tensor. Consequently, an earthquake due to a pure double-couple can be considered as a special case of the moment tensor.
For big earthquakes, i.e. Mw>=6.0, the source time function also should be included in the inversion. In such cases, the unknown model parameters are the 6 components of the moment tensor and the source time function components.
In this section, the inversion was carried out under pure double-couple conditions. In other words it was assumed that the sum of moment tensor diagonal components was zero and the eigenvalues of the moment tensor were assumed to be 1, 0, and.
Different parameters may affect the solution in the linear inversion method e.g. the frequency band, the number of components involved and the way the stations are distributed. The effects of the most important parameters are discussed below.
A. The Effect of Frequency Band on the Inversion Solution: The linear inversion method was applied for different frequencies. In most cases, the earthquake specifications obtained were in complete agreement with the simulated earthquake which indicates that in the tested frequency bands, all the frequencies were below the corner frequency of the earthquake source. The solution obtained was very similar to the original mechanism. The fault plane resulted was perpendicular to the original plane that indicates the major fault plane and minor fault plane.
B. The Effect of the Number of Components Involved: The linear inversion method was applied for three different cases. First, the vertical components were used. Second, only the radial components were applied. Third, the inversion was varied solely by tangential components. In the first and second cases, the source coordinates were obtained precisely. But for the third one, the solution was completely wrong even though the depth was determined precisely.
C. The Effect of the Station Distribution Pattern: The effect of the station distribution was examined based on the configuration of seismographic stations in different quadrants. It was concluded that the solution can be obtained even though the data are from one quadrant but the more quadrant are involved the less error there will be.
2. Simulation of an Earthquake and its Linear Inversion for a Non-Pure Double-Couple (without volume variation)
2.1 Earthquake Simulation: In this mode, a five layered crustal model was used for simulation of synthetic seismograms in eight stations distributed in the four quadrants with respect to the source.
2.2 Linear Inversion of the Seismograms for Determination of the Earthquake source parameters under Non-Pure Double-Couple Conditions: In contrast to the linear inversion in 1.2, in this case, the inversion was used in conditions which were more similar to the reality. In the natural mode, the elastic waves traverse layers in the earth about which we don’t have enough information. Therefore, in earthquake mechanism determinations, models are used which are much more simplistic than the real one.
Like the other case, each of the parameters affecting the solution in the linear inversion method will be examined briefly
A. The Effect of Frequency Band on the Inversion Solution: The linear inversion method was applied for different frequencies. In every case, a parabolic function with the <br /> corner frequency of 0.2Hz was used. In all the cases, the earthquake specifications were restored and, as in the other case, the depth was determined within 2 or 3 km of the original depth.
B. The Effect of the Number of Components Involved: Again the linear inversion method was applied for three different cases. And the results indicated that given the mechanism, the epicenter distances and the depth of the earthquake, the vertical components are of greater amplitudes, and therefore stabilize the solution to the inversion method.
C. The Effect of the Station Distribution Pattern: Again the effect of the station distribution was examined based on the configuration of seismographic stations in different quadrants and accordingly the solution could be obtained even though the data were from one quadrant but the more quadrants there were involved the less error there was.
Conclusion: Two different types of sources (pure and non-pure double-couple) were <br /> used to produce synthetic seismograms based on the wavenumber integration method <br /> for a given velocity model. In both cases, the source model was obtained precisely depending on the conditions. Generally, it can be concluded that if the velocity <br /> model is not precise its effect can be seen on CLVD as well as the depth. In addition, if the stations are distributed in at least two quadrants, more precise solutions will be gained.در این پژوهش مفاهیم اصلی سازوکار چشمه زمینلرزه و عوامل مؤثر بر آن با استفاده از تولید لرزهنگاشت مصنوعی و برگردان خطی آنها در حوزه زمان بررسی میشود. در ابتدا از لرزهنگاشتهای مصنوعی تولید شده ناشی از چشمههای دو زوج نیروی خالص که نمایانگر اکثر سازوکارهای زلزلههای ایجاد شده با عوامل تکتونیکی هستند استفاده شد. همچنین از چشمههای زمینلرزه دو زوج نیرو ناخالص (Double-Couple+compensated linear vector dipole, CLVD) که تحت تأثیر عوامل تکتونیکی و غیر تکتونیکی (نظیر آتشفشان) هستند نیز بهره گرفته شد. لرزهنگاشت مصنوعی براساس روش انتگرالگیری روی عدد موج برای مدل سرعت پوسته فرضی محاسبه شده است. درتعیین سازوکار به روش برگردان، اثر محدوه بسامدی، غیر دقیق بودن مدل سرعت، تعداد مؤلفهها و نحوه پراکندگی ایستگاهها مورد بررسی قرار گرفته است.https://jesphys.ut.ac.ir/article_79972_ec8643187d580cddd7c59f97ad0fed7f.pdfموسسه ژئوفیزیک دانشگاه تهرانفیزیک زمین و فضا2538-371X35120090421Numerical simulation and experimental investigation of the thermal diffusivity of frozen soil under different moisture contents conditionsشبیهسازی عددی و بررسی آزمایشگاهی ضریب پخشیدگی گرمایی خاک یخزده در شرایط رطوبتی متفاوت89997997310.22059/jesphys.2009.79973FAیونسخوشخودانشآموخته کارشناسی ارشد هواشناسی کشاورزی، دانشکده مهندسی خاک و آب، پردیس کشاورزی و منابع طبیعی دانشگاه تهرانعلیخلیلیاستاد بازنشسته، دانشکده مهندسی خاک و آب، پردیس کشاورزی و منابع طبیعی دانشگاه تهرانحسنرحیمیاستاد بازنشسته، دانشکده مهندسی خاک و آب، پردیس کشاورزی و منابع طبیعی دانشگاه تهرانپرویزایراننژاداستادیار، گروه فیزیک فضا، مؤسسة ژئوفیزیک دانشگاه تهران، ایرانJournal Article20210221Soil thermal diffusivity is considered as the most important thermal characteristic of the
soil which indicates the gradient of its warming due to a unit change in its temperature. Several methods are available to determine soil thermal diffusivity from observed temperature variations. Most of these methods are based on solutions of the one-dimensional conduction heat equation with constant diffusivity and thus apply to uniform soils only.
In the absence of local heat sources or sinks, the equation describing conductive heat transfer in a one-dimensional isotropic medium is:
(1)
where is temperature, is time, is the soil depth, and () is the thermal diffusivity of the soil , equal to and being the volumetric heat capacity . Several methods have been developed for estimating the soil diffusivity using equation (1). Horton et al. (1983) have tested six methods and concluded that Harmonic Equation and Numerical Method provide the most accurate results among all. The finite difference is considered as the most applicable method for numerical solution of the heat conduction equation in soils. For approximation of partial derivatives using finite differences, different algorithms may be used. In the present research, the Crank-Nicolson method which has a high degree of accuracy was employed. Using the above method, equation (1) can be discretised as:
(2)
with
and .
where and indicate the depth node and time step, respectively.
In the present work, for numerical solution of the equation, time intervals of 1 second (Δt = 1s) and spatial intervals of 1cm (Δz = 1cm) have been employed. With application of the Crank-Nicolson method, a set of simultaneous equations will be produced for each time interval. This set of simultaneous equations can be solved using different methods. In the present research the Tri-Diagonal Matrix Algorithm (TDMA) method has been employed. When the initial and boundary conditions are known and soil temperatures at different points have been measured, the soil thermal diffusivity can be determined using a trial and error technique. The approach is based on solving equation (2) iteratively by changing , and determining the -value based on which the calculated values of temperature best match observations. In the present work the criterion used for choosing is minimizing the Root Mean Square Error (RMSE) of the calculated (C<sub>i</sub>) against observed (Mi) temperature:
The above procedure was conducted to evaluate the thermal diffusivity (α) of a silty soil with mass moisture contents of 5, 10, 15 and 20 per cent. A chamber with dimensions of 500×500×800 mm was made and its walls and bottom were carefully insulated using layers of plasto-foam sheets with a thickness of 100mm, to minimize the exchange of heat with the surrounding environment. The chamber was filled with the soil of given texture and moisture. To eliminate evaporation from the top surface of the soil, it was covered by a plastic sheet. Temperature was measured using seven thermometers installed at the depths of 50, 110, 170, 250, 350, and 500 mm, as well as at the top surface of the soils. The sensors were connected to a computer, where soil temperatures were recorded at 1 min intervals. A frost condition in soil was simulated with the use of a cooling system located at the top of the soil that was able to produce temperatures as low as -20 <sup>o</sup>C.
Thermal diffusivity was estimated for two different thermal conditions in the soil profile: one with temperatures lower than -2 <sup>o</sup>C throughout the soil at all times, and the other having temperatures lower than zero degrees centigrade at some and higher than zero at other depths.
According to the results, the model used in this study led to a low RMSE (between 0.41 to 0.71 <sup>o</sup>C) and reasonable predictions of soil temperature for the first case (i.e. temperatures lower than -2 <sup>o</sup>C at all times and depths). The results showed that the values of α increased with increasing moisture content up to a critical point and then decreased. The maximum value of α occurred at 15 percent moisture content.
The model failed to estimate the soil temperature profile within an acceptable range of error in the second case. The range of RMSE values between the simulated and measured temperature in this case was found between 1.58 to 2.76 <sup>o</sup>C. This failure was attributed to the fact that the assumptions made in solving the heat conduction equation, namely homogeneity of soil and lack of sources and sinks of heat within the soil, were not fulfilled for the second case.در این تحقیق ضریب پخشیدگی گرمایی خاک سیلت رسی در رطوبتهای وزنی 5، 10، 15، و 20 درصد در حالت وقوع پدیده یخبندان در خاک مورد بررسی قرار گرفته است. برای برآورد ضریب پخشیدگی گرمایی ابتدا معادله پخش گرما در خاک به روش کرانک- نیکسلون برای مقادیر متفاوت فرضی در محدوده تغییرات آن حل شد. سپس با استفاده از معیار جذر متوسط مربعات خطا (RMSE) بین مقادیر دمای خاک محاسباتی و مشاهداتی در زمانها و اعماق متفاوت، مقدار بهینه به روش سعی و خطا برای رطوبتهای پیشگفته استخراج شد. برای جمعآوری دادههای مورد نیاز در حل معادله پخش گرما در خاک، از 7 حسگر گرمایی در اعماق گوناگون خاک استفاده شد. برای دستیابی به شرایط یخبندان در خاک از دستگاه تبرید که قادر به تولید سرما تا دمای 20- درجه سلسیوس بود استفاده به عمل آمد. استخراج برای دو مقطع مجزا از دادههای مشاهداتی به انجام رسید. مقطع اول که نمایانگر شرایط وقوع یخبندان کامل خاک بود از بین دادههایی انتخاب شد که در آن دمای خاک در کلیه اعماق خاک و در زمانهای متفاوت مقادیری منفی و کمتر از 2- درجه سلسیوس را داشت و مقطع دوم که نمایانگر شرایط نبود یخبندان کامل در خاک بود طوری انتخاب شد که در آن دمای خاک در بخشی از مقطع انتخاب شده مثبت و در بخشی دیگر منفی بود. طبق نتایج بهدست آمده، مدل بهکار گرفته شده در این تحقیق در حالت یخبندان کامل منجر به حصول مقادیر RMSE کم و بهدست آمدن نتایجی منطقی و قابل قبول شد ولی در حالت نبود یخبندان کامل مقادیر RMSE به طرز قابل توجهی زیاد شد و این امر سبب ناکارامد شدن مدل و بهدست آمدن نتایجی غیرمنطقی شد که وقوع این امر به نقض فرضیات مورد نیاز معادله و اعمال خطاهای اضافی در اثر آن در حالت اخیر نسبت داده شد.https://jesphys.ut.ac.ir/article_79973_a383249774120011cdcc467cb4139aa2.pdfموسسه ژئوفیزیک دانشگاه تهرانفیزیک زمین و فضا2538-371X35120090421The role of convection parameterization in the simulation of the winter temperature and precipitation fields over Iran using Regional Climate Model (RegCM3)نقش روشهای متفاوت پارامترسازی همرفت در شبیهسازی میدانهای دما و بارش زمستانی با مدل منطقهای- اقلیمی RegCM در منطقه ایران1011207997410.22059/jesphys.2009.79974FAپرویزایراننژاداستادیار، گروه فیزیک فضا، مؤسسة ژئوفیزیک دانشگاه تهران، ایرانفرهنگاحمدیگیویاستادیار، گروه فیزیک فضا، مؤسسة ژئوفیزیک دانشگاه تهران، ایرانروزبهپازوکیدانشآموخته کارشناسی ارشد هواشناسی، مؤسسة ژئوفیزیک دانشگاه تهران، ایرانJournal Article20210221Convection affects the climate through its role in the redistribution of energy and moisture in the atmosphere, and subsequently producing clouds and precipitation. Even with the recent improvements in computational power, the numerical models of the atmosphere still have to run in spatial resolutions that are too coarse to capture the local scale processes such as convection. For this reason, and because of the importance of convection for the surface climate, parameterization schemes have been developed to empirically upscale convection to the scale of the model grid areas.
This study aims at evaluating the impact of different types of parameterizing convection on the simulations by Version 3 of the Regional Climate Model (RegCM3: Dickinson et al., 1989; Giorgi, 1989) of precipitation and air temperature. The convection schemes currently coupled with RegCM3 are Anthes (1997), Betts (1986) and Grell (1993), which may be closed using either the Arakawa-Schubert (1974) or the Fritch and Chappell (1980) closure schemes. The simulations are conducted for the four-month period of December 1998 to March 1999 (inclusive) with 45 × 45 km grid spacing over a domain having 60 and 70 grid points along, latitude and longitude respectively and centered in Iran at 34 ˚E and 48 ˚N. The initial and boundary conditions are derived from the NCEP/NCAR reanalysis. RegCM3 was run four times, keeping all the components of the model and the initial and boundary conditions the same, by each time coupling one the convection schemes (Anthes, Betts, Arakawa, Fritch and Chappell) with the model. To minimize the impact of possibly incorrect initial conditions, we assumed one month as the model's spin-up period, and analyzed the results for the three months of January to March 1999.
The simulated monthly mean precipitation and air temperature as well as the spatial distribution of the model outputs using different schemes are intercompared and compared with observations from the Climate Research Unit (CRU) at the University of East Anglia, United kindom. The results show that the grid-scale mean monthly and mean seasonal (winter) near-surface air temperature simulated by RegCM3 coupled with different convection schemes agree very well with the corresponding observed values. The slope of the regression line of the simulated mean winter temperature against observations is very close to one and varies between 0.973 (Fritch-Chappell) and 0.997 (Kuo-Anthes), with the coefficient of determination (R2) in the range of 0.941 to 0.944, respectively. The degree of agreement between the simulated monthly mean temperature with observations for the three months is somewhat lower that of the mean seasonal with the slope of regression line varying between 0.908 and 0.939 and coefficient of determination between 0.935 and 0.938. It is concluded that RegCM3 is highly effectual in simulating air temperature, irrespective of the type of the convection scheme used. The differences between the simulated mean temperature using different schemes are very small.
On the other hand, the effectiveness of RegCM3 in simulating monthly and seasonal precipitation is very low. Differences between simulated monthly and seasonal precipitation using the four convection schemes are negligible. Although the geographical distribution pattern of precipitation is well simulated by the model, the simulated monthly and seasonal precipitation regressed against observation shows that RegCM3 generally underestimates precipitation during the winter months. The slope of the regression lines significantly differs from unity, varying between about 0.570 and about 0.715. The highest coefficient of determination found for the four schemes, during the three months and the season is smaller than 0.280.
Given the insignificant differences among the model simulations using any of the four convection schemes, the simplest form of convection parameterization with the lowest computational costs, i.e. the Kuo-Anthes scheme, proved to be most appropriate available scheme for medium-range weather predictions in Iran.همرفت از راه باز توزیع انرژی و رطوبت در جوّ و در نتیجه ایجاد ابر و بارش در اقلیم نقش ایفا میکند. حتی با پیشرفتهای اخیر در توان محاسباتی، هنوز هم مدلهای عددی جوّ با درجه تفکیکی اجرا میشوند که برای نشان دادن پدیدههای با مقیاس محلی نظیر همرفت، بسیار کوچک است. به این دلیل و همچنین به دلیل اهمیت همرفت در اقلیم سطحی طرحوارههای پارامترسازی تهیه شدهاند که بهطور تجربی اثر همرفت را در مقیاس شبکهای مدل بیان میکنند.
هدف تحقیق حاضر ارزیابی تأثیر روشهای متفاوت پارامترسازی همرفت در شبیهسازی میدانهای بارش و دمای هوا با نسخه سوم مدل اقلیمی منطقهای (RegCM3) است. شبیهسازیها با گام مکانی km45×km45 برای دوره چهار ماهه دسامبر 1998 تا مارس 1999 برای منطقهای به مرکز 34 درجه شرقی و 48 درجه شمالی، شامل 60 نقطه شبکه در راستای نصفالنهار و 70 نقطه در راستای مدار، صورت گرفته است. شرایط مرزی برگرفته از دادههای بازتحلیل شده NCEP/NCAR است. در این مقاله، فقط نتایج سه ماه ژانویه تا مارس 1999 تحلیل شدهاند. مدل RegCM چهار بار برای سه ماه زمستان با شرایط مرزی یکسان ولی هر بار با استفاده از یکی از طرحوارههای آراکاوا- شوبرت، فریچ- چپل، بتس- میلر و کو– آنتس اجرا شده است. شبیهسازیهای میانگینهای ماهانه بارش و دمای هوا و هم چنین توزیع مکانی نتایج مدل با استفاده از طرحوارههای مختلف با یکدیگر و با دادههای مشاهداتی واحد پژوهشهای اقلیمی (CRU) دانشگاه ایستآنجلیا (بریتانیا) مقایسه شده است. نتایج شبیهسازیهای دما و بارش با استفاده از طرحوارههای مختلف همرفت در دوره مورد مطالعه در ایران بسیار شبیه یکدیگر است و توان پیشبینی مدل RegCM برای میدان دما بالاتر از بارش است. مقادیر بارش مدل تفاوت بسیار زیاد از مشاهدات دارد، ولی مدل توزیع جغرافیایی بارش را به خوبی شبیهسازی میکند. با توجه به ساختار سادهتر و هزینه محاسباتی کمتر طرحواره کو- آنتس نسبت به طرحوارههای دیگر، با وجود عدم تفاوت چشمگیر در نتایج مدل، بهنظر میرسد این طرحواره در مطالعه اقلیمی و منطقهای در بازه های فصلی، لااقل در منطقه ایران بهترین گزینه باشد.https://jesphys.ut.ac.ir/article_79974_8b32f41af27b68ef1fc85122b636468c.pdfموسسه ژئوفیزیک دانشگاه تهرانفیزیک زمین و فضا2538-371X35120090421The 500 hpa atmospheric centers of action and circulation patterns over the Middle East and their relationship with precipitation in Iranمراکز فعالیت و الگوهای گردش جو زمستانه تراز 500 هکتوپاسکال روی خاورمیانه و ارتباط آنها با بارش ایران1211417997510.22059/jesphys.2009.79975FAطیبرضیئیاستادیار، مرکز تحقیقات حفاظت خاک و آبخیزداری کشور، تهران، ایرانعباسمفیدیاستادیار، گروه جغرافیای مؤسسه آموزش عالی طبرستان، چالوس، ایرانآذرزرّیناستادیار گروه جغرافیای مؤسسه آموزش عالی طبرستان، چالوس، ایرانJournal Article20210221It is well-known that regional weather and climate around the globe are strongly influenced by large-scale atmospheric circulation patterns. The centers of action corresponding to different levels of atmosphere play an essential role in controlling the climate of different climatic regions around the globe. Over the years several efforts have been made to identify the main centers of action and large-scale atmospheric circulation patterns leading to precipitation events and to study how their variability can affect the frequency and intensity of precipitation. Hence, using weather types or circulation patterns one is able to investigate and explain the physical causes for the frequency and intensity variation of precipitation over a region.
Reviewing synoptic studies in relation to Iran have suggested that in spite of many subjective circulation classifications being implemented using observational data mainly on the monthly basis; few attempts have been done objectively using reanalysis data, especially on the basis of daily data. Hence, this paper aims to identify the main centers of action and circulation patterns in relation to winter precipitation variability over Iran.
To recognize winter atmospheric circulation patterns over the Middle East, the mean daily 500 gph for December, January, February and March were retrieved and used from the National Centers for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR) reanalysis archive, covering the period from January 1965 to 2000, accounting for 4355 days during a period of 36 years. Subsequently, daily precipitation rates for some selected days were also retrieved from the NCEP/NCAR reanalysis archive in order to assess the influences of the identified circulation patterns on precipitation in Iran.
To classify the 4355 days and extract the main circulation patterns, S-mode PCA was applied to the data matrix and the resultant 9 leading PCs were retained based on the scree test. The retained PCs were then rotated using the varimax criterion. By plotting the rotated PC loadings the centers of action for 500 hpa level that control the winter climate of the Middle East were identified. To analyze the corresponding synoptic characteristics of the identified centers of action, 10 days with the highest PC score (positive phase) were selected for each PC. The composite map of the selected days and the corresponding vorticity maps for 500 and 1000 hpa levels were represented as the winter circulation patterns. Finally, by composing the precipitation rates associated to the each aforementioned selected 10 days, the relationship between identified synoptic circulation patterns over the Middle East and the winter precipitation in Iran was investigated.
The results indicate that the spatial pattern of winter precipitation in Iran, with the exception of the southern coastal areas of the Caspian Sea, is properly governed by the 500 hpa circulation patterns. The results also show that the large spread water deficit and dry periods in Iran are related to the northward displacement of the Arabian high pressure at mid troposphere level over the western part of the Middle East. Moreover, the results indicate that the deepening of the westerly wave and the increasing positive vorticity in the area between Iran and the southern part of the Red Sea, accompanied with the development and/or reinforcement of a high pressure in the area between eastern Saudi Arabia and the central part of the Red Sea is responsible for widespread precipitation occurrences over vast areas of western and southwestern parts of Iran. Investigation of the relationship between synoptic circulation patterns and the regional scale winter precipitation in Iran shows that the precipitation occurrence in the Caspian region is mostly related to the position and strength of the lower atmosphere high pressures rather than the mid troposphere circulation patterns. This is evident if we consider that in about 4 out of 9 identified circulation patterns the southern coastal areas of the Caspian Sea get remarkable precipitation due to the predominance of a high pressure system and an increase in the negative vorticity in mid troposphere over the western part of the Caspian Sea as well as the development and/or prolonging anticyclonic circulation and northerly flows over the Caspian Sea.به منظور شناسایی الگوهای گردش جو زمستانه روی خاورمیانه و تعیین میزان ارتباط آنها با بارشهای ایران از دادههای ارتفاع ژئوپتانسیل تراز 500 هکتوپاسکالِ ماههای دسامبر، ژانویه، فوریه و مارس برای تعداد 4355 روز در یک دوره آماری 36 ساله (1965-2000) استفاده شد. به منظور گروهبندی دادهها و استخراج الگوهای اصلی، تحلیل مؤلفههای اصلی با آرایه S روی ماتریس دادهها صورت گرفت و به کمک آزمون غربالی (Scree Test) 9 مؤلفه نخست انتخاب و به روش وریمکس (Varimax) چرخانده شد. با ترسیم نقشه بارگویه (Loading) مؤلفههای چرخشیافته مراکز فعالیت (Centers of Action) تراز 500 هکتوپاسکال که اقلیم زمستانه خاورمیانه را کنترل میکنند، شناسایی شد. به منظور تحلیل ویژگیهای همدیدی مرتبط با این مراکز فعالیت، برای هر مؤلفه تعداد 10 روز با بالاترین نمره استاندارد (فاز مثبت) تعیین و نقشههای ترکیبی (Composite) ارتفاع ژئوپتانسیل و تاوایی نسبی آنها برای تراز500 هکتوپاسکال و نیز نقشه های ترکیبی فشار سطح دریا و تاوایی نسبی آنها در تراز1000 هکتوپاسکال تهیه و درحکم الگوهای گردش جو زمستانه خاورمیانه ارائه شد. سرانجام ارتباط بین الگوهای همدیدی زمستانه خاورمیانه و بارش ایران زمین مورد بررسی و تحلیل قرار گرفت.
نتایج بررسیها بیانگر آن است که الگوی فضایی توزیع بارش زمستانه ایران به غیر از سواحل جنوبی دریای خزر در سایر مناطق کشور به خوبی از الگوهای گردش تراز500 هکتوپاسکال پیروی میکند. بررسیها روشن ساخت که دورههای خشک و کم بارش فراگیر ایران با تقویت و گسترش شمالسوی پرفشار عربستان در ترازهای میانی وردسپهر روی نیمه غربی خاورمیانه در ارتباط است. همچنین یافته ها نشاندهنده آن است که عمیق شدن ناوة موج غربی و افزایش تاوایی مثبت در حد فاصل غرب ایران تا نیمه شمالی دریای سرخ به همراه شکلگیری و تقویت پرفشار در حد فاصل شرق عربستان تا بخشهای میانی دریای عرب امکان وقوع بارشهای فراگیر و قابل ملاحظه را در بخشهای وسیعی از غرب و جنوب غرب ایران فراهم میآورد. بررسی ارتباط بین الگوهای همدیدی و بارشهای مقیاس منطقهای حاکی از آن است که وقوع بارش در منطقه خزری بیش از آنکه با الگوی گردش <br /> وردسپهر میانی مرتبط باشد به موقعیت و شدت مراکز پرفشار ترازهای زیرین جو وابسته است. بهطوری که در 4 الگو از 9 الگوی شناسایی شده، بهواسطه استقرار مرکز پرفشار و افزایش تاوایی منفی در ترازهای زیرین جو بر جانب غربی دریای خزر و به دنبال <br /> آن شکلگیری و تداوم گردش واچرخندی و جریانهای شمالی روی این دریا، سواحل جنوبی دریای خزر از بارش قابل ملاحظهای برخوردار شده است.https://jesphys.ut.ac.ir/article_79975_efea08a2d7e4c031b476735a08459474.pdfموسسه ژئوفیزیک دانشگاه تهرانفیزیک زمین و فضا2538-371X35120090421Orbit integration in non-inertial framesمحاسبه مدار در چارچوبهای نالَخت187997610.22059/jesphys.2009.79976FAمهدیاسحاقدانشجوی دکتری، مؤسسه سلطنتی صنعتی استکهلم، سوئدJournal Article20210221A precise orbit of a low Earth orbiting satellite helps us to compute the long wavelength structure part of the gravity field of the Earth. There are different methods and frames for orbit integration depending on the problem and satellite mission.
In this paper, the dynamic equations of the satellite motion are presented in different frames of navigation. A simple numerical study on a satellite orbit in local frames is also included. In these frames the geodetic coordinate of the satellite is directly integrated. Numerical studies confirm that the north-east-down frame is not stable for orbit integration either. However, the paper shows how to solve this problem by choosing a wander frame and its ability.مدار دقیق ماهوارههای نزدیک سطح زمین در مطالعه ساختار بلند میدان جاذبه زمین بسیار مفیدند. روشها و چارچوبهای گوناگون برای محاسبه چنین مدارهای موجود هست که با توجه به مسئله مورد نظر و ماموریت ماهواره انتخاب میشوند. در این مقاله معادلات حرکت یک ماهواره در چارچوبهای مختلف ارائه میشوند. محاسبات عددی بیانگر این مطلب هست که چارچوبهای محلی متمایل به شمال زمین مناسب برای محاسبه مدار نیستند ولی استفاده از چارچوب معلق که حالت خاصی از چارچوب محلی هست امکان حل مدار را در چنین چارچوبهائی نیز حاصل مینماید.
چارچوبهای انتگرالگیری مدار عبارتاند از : چارچوب لَخت، چارچوب بیضوئی، محلی و معلق. چارچوب لخت یکی از سادهترین چارچوبهای انتگرالگیری است و به تناوب از سوی محققین گوناگون در زمینه ناوبری و ژئودزی فضائی مورد استفاده قرار میگیرد. چارچوب بیضوی بسیار مشابه به چارچوب لخت است ولی بهجای انتگرالگیری نسبت به مختصات دکارتی ماهواره، محاسبات مستقیما روی عرض، طول و ارتفاع ژئودتیکی صورت میگیرد. چارچوبهای محلی بیشتر مناسب برای ناوبری هواپیما و روشهای هوابرد گرانیسنجیاند و به علت همگرائی نصفالنهارها در مناطق نزدیک قطب انتگرالگیری از شتاب در باند ارتفاع دچار ناپایداری میشود و چنانچه طول مدت پرواز طولانی باشد، ناپایداری از این باند روی محاسبه عرض و طول ژئودتیکی نیز تاثیر میگذارد و در نتیجه آنها را نیز ناپایدار میسازد. یکی از روشهای حل این ناپایداری استفاده از چارچوب معلق است که به چارچوب محلی شباهت دارد ولی وجود یک آزیموت معلق جهت مقابله با همگرائی نصفالنهارها ناپایداری را برطرف میکند. به عبارت دیگر چارچوب معلق انعطافپذیرتر از محلی است. مبحث ناوبری ماهوارهای کاملاً مشابه با ناوبری هوائی است، با این تفاوت که ماهواره به تناوب از نزدیکی قطبها میگذرد و درگیر با تقارب نصفالنهارها است و در هر گردش به دور زمین، به ناپایداری در باند ارتفاعی دچار میشود و چون ماموریت ماهوارهها طولانی است، قطعاً این ناپایداری بر موقعیت مسطحاتی ماهواره نیز مؤثر خواهد بود. همانطور ذکر شد یکی از روشهای حل این ناپایداری، استفاده از چارچوب معلق است. در این مقاله ما معادلات حرکت ماهواره را به این چارچوب گسترش داده و به شکل زیر سادهسازی کردهایم:
که در این رابطه شعاع انحنای نصفالنهاری و شعاع انحنای دایره قائم اولیه در نقطهای با عرض و ارتفاع ژئودتیک h و است. عرض ژئودتیک، اولین خروج از مرکز، دومین خروج از مرکز بیضوی مرجع، آزیموت معلق و , و مؤلفههای بردار سرعت، ، و مؤلفههای بردار گرانش، ، مؤلفههای بردار شتاب در چارچوب معلق است.
شکل 1 مدار یک ماهواره نزدیک سطح زمین در یک دور گردش به دور زمین را در چارچوب معلق نمایش میدهد. همانطور که در شکل مشاهده میشود هیچگونه ناپایداری در باند ارتفاعی و همچنین سرعت ارتفاعی وجود ندارد که بیانگر صحت محاسبهها و کارآمد بودن چارچوب معلق در حل مدار ماهواره است. همانطور که شکل نشان میدهد آزیموت معلق، برای جبران تقارب نصفالنهارها تا 160 درجه تغییر میکند.https://jesphys.ut.ac.ir/article_79976_bd5390a0ee10da09cc8f4de70a31de52.pdfموسسه ژئوفیزیک دانشگاه تهرانفیزیک زمین و فضا2538-371X35120090421Detection of subsurface Qanats by Artificial Neural Network via Microgravity dataاکتشاف قناتهای زیرزمینی مدفون از طریق شبکههای عصبی مصنوعی و با استفاده از دادههای میکروگرانیسنجی9157997710.22059/jesphys.2009.79977FAعلیرضاحاجیانمربی، گروه فیزیک، دانشگاه آزاد اسلامی، واحد نجفآباد، اصفهان، ایرانوحیدابراهیمزاده اردستانیدانشیار، گروه فیزیک زمین، مؤسسة ژئوفیزیک دانشگاه تهران و قطب علمی مهندسی نقشهبرداری و مقابله با سوانح طبیعی، ایرانکارولوکساستاد، گروه برق، قطب علمی کنترل و هوش مصنوعی، دانشکده برق و کامپیوتر، دانشگاه تهران، ایرانسیدمرتضیسقائیاننژاددانشیار، گروه برق، دانشکده برق و کامپیوتر، دانشگاه صنعتی اصفهان، ایرانJournal Article20210221A full automatic algorithm is designed to detect subsurface Qanats (sub terrains) via Artificial Neural Networks .We first gained the residual gravity anomaly from microgravity data and then applied it to a Multi Layer Perceptron (MLP) which was trained for the models of sphere and cylinder.
As a field example, the depth of a subsurface Qanat buried under the north entrance of the Geophysics Institute is determined through MLP (trained with noisy data).در این مقاله یک الگوریتم هوشمند جهت اکتشاف قناتهای زیرزمینی مدفون با شبکههای عصبی و با استفاده از دادههای میکروگرانیسنجی ارائه شده است.
به منظور برآورد عمق و اندازه قناتهای زیرسطحی از روی بیهنجاری (آنومالی) گرانی باقیمانده یک شبکه عصبی مصنوعی با سرپرست، از نوع پرسپترون چندلایه (MLP) طراحی شد. از آنجاکه در طراحی شبکه عصبی سرعت پردازش دادهها از اهمیت <br /> خاصی برخوردار است و تعداد ورودیهای زیاد باعث پیچیدگی غیر منطقی توپولوژی شبکه میشود، به جای اعمال همة دادههای تصحیح شده، میکروگرانی درحکم ورودی، مجموعهای مشخصههای مناسب (Features) از روی آنومالی باقیمانده دادههای میکروگرانی استخراج میشود، سپس با توجه به مدلهای کره و استوانه که نزدیکترین مدلها به قناتهای مدفون هستند، مجموعهای از دادههای آموزشی که برای آموزش شبکه عصبی طراحی شدهاند، مورد استفاده قرار میگیرند، در واقع شبکه عصبی طراحیشده پس از این آموزش قادر خواهد بود که با توجه به مشخصههای استخراج شده از روی بیهنجاری باقیمانده، عمق و شعاع قنات مدفون را بهدست آورد.
از آنجاکه قاعده کلاسیک خاصی برای انتخاب تعداد نورونها در لایه پنهان شبکه عصبی چندلایه وجود ندارد، شبکههای عصبی چندلایه گوناگونی با تعداد نورونهای متفاوت در لایه پنهان مورد آزمایش قرار گرفت و نمودارهای عملکرد شبکه در هر حالت بهدست آمد تا از روی آن بهترین مقدار تعداد نورونها در لایه پنهان حاصل شود.
پس از این مرحله ابتدا با استفاده از مجموعهای دادههای مصنوعی، شبکه عصبی طراحیشده مورد آزمون قرار گرفت. <br /> سپس خروجیهای شبکه با استفاده از دادههای مصنوعی نوفهدار برای مدلهای کره و استوانه بررسی شد که عملکرد مناسبی را نشان داد.
همچنین، عمق قنات زیرزمینی مدفون واقع در ورودی شمالی مؤسسه ژئوفیزیک درحکم نمونهای عملی با شبکه عصبی طراحی شده، بهدست آمد که با مقدار واقعی آن انطباق خوبی داشت.https://jesphys.ut.ac.ir/article_79977_7c67200d9342048b353733d50fb5b051.pdfموسسه ژئوفیزیک دانشگاه تهرانفیزیک زمین و فضا2538-371X35120090421Gravity field implied density modeling of topography, for precise determination of the geoidتعیین دقیق ژئوئید برای مدل چگالی بهدست آمده از میدان ثقل زمین17317997810.22059/jesphys.2009.79978FAمهدینجفیعلمداریدانشیار، دانشکده مهندسی نقشهبردارای، دانشگاه صنعتی خواجهنصیرالدین طوسی، تهران، ایرانمرتضیصدیقیکارشناس ارشد سازمان نقشهبردداری کشور، تهران، ایرانسیدهاشمطباطباییمدیر عملیات لرزهنگاری شرکت ملی نفت ایرانJournal Article20210221Precise determination of the geoid using the Stokes-Helmert approach requires a density distribution model within the topography. The model is used for precise evaluation of topographical indirect effects on gravity and potential applied in transforming between the real and the Helmert spaces. The range of mass density variation within the topography is between 1000 and 3100 kg.m<sup>-3</sup>. Assigning the global average value of kg.m<sup>-3</sup> at a point, instead of its real point value, may cause errors of decimeter magnitude in the geoid determination. A regional-local gravity anomaly separation technique using Bouguer gravity anomaly (BA) along with the Free air gravity Anomaly (FA) in the region of Iran are used to estimate the local topographical effect on gravity from the observed anomalies after eliminating the non-density origin long wavelength including isostatic features into the observed BA. A Global Geopotential Model (GGM) is also used to eliminate the deep sited density-origin long wavelength features from the observed anomalies as well. Then, the power spectral analysis, apparent density mapping, and forward modeling techniques are used to convert the local topographical effect on gravity into the corresponding 3-D density model (GRADEN) model in the region. The model showed thorough correlation with the superficial geological density (GEODEN) model at the surface level, provided that the reliable digitization of the model is in order. The GRADEN model minus the constant density demonstrates contributions up to a meter in mountainous areas and 7cm in the RMS scale to the geoid in the region.در تعیین دقیق ژئوئید به روش استوکس- هلمرت، به مدل توزیع چگالی برای توپوگرافی زمین نیاز است. مدل چگالی توپوگرافی برای تعیین دقیق اثرات توپوگرافی روی شتاب ثقل و پتانسیل و تبدیل بیهنجاری (آنامولی) جاذبه مشاهداتی (فضای واقعی) به بیهنجاری هلمرت (فضای هلمرت) بهکار میرود. مقدار عددی چگالی در سطح زمین بین 1 تا 1ر3 (gcm<sup>-3</sup>) متغیر است ولی در بیشتر موارد از
مقدار متوسط 67/2 گرم بر سانتیمتر مکعب برای تعیین ژئوئید مورد قرار میگیرد که موجب خطا در محاسبه ارتفاع ژئوئید در حدود دسیمتر میشود.
برای تعیین چگالی توپوگرافی، ابتدا اثرات طول موجهای بلند موجود در بیهنجاری جاذبه بوگه ناشی از بیهنجاریهای چگالی واقع در اعماق و دیگر اثرات طول موج بلند نظیر اثر ایزوستاسی و لایههای سبک با گستره وسیع احتمالی در منتل بالایی و عوارض سبک ساختارهای زمینشناسی با گستره وسیع احتمالی منطقهای را با استفاده از مفاهیم فیزیکی و در مقایسه بیهنجاری جاذبه بوگه با بیهنجاری جاذبه هوای آزاد در منطقه حذف میکنیم و سپس با استفاده از مدلهای جهانی ژئوپتانسیل، بیهنجاری باقیمانده محلی را جدا میسازیم. با استفاده از آنالیز طیفی و بهکارگیری فیلترهای متفاوت، این جداسازی را ادامه میدهیم و آنامولیهای جدا شده را منتسب به اعماق متفاوت نسبت میدهیم و مدل سهبعدی پوسته (تا عمق حدود 20 کیلومتر) را ایجاد میکنیم. در روندی تکراری و با کمک روش معکوس، مدل چگالی را تصحیح میکنیم، یعنی مقدار تفاوت جواب مدل با مشاهدات را به حد قابل قبولی کاهش میدهیم.
میزان اثر بیهنجاری چگالی در ارتفاع ژئوئید حداقل و حداکثر 73- و 103 سانتیمتر با انحراف معیار و RMS به ترتیب 7 و 8 سانتیمتر است. با توجه به اعداد بهدست آمده میتوان نتیجه گرفت که برای تعیین ژئوئید دقیق اطلاع از میزان چگالی واقعی توپوگرافی لازم است.https://jesphys.ut.ac.ir/article_79978_94ff38d75ace1e1a7b3744fe977930a8.pdf