تخمین عمق موهو براساس وارون سازی داده های آنامولی جاذبه از منابع مختلف مطالعه موردی: منطقه خراسان

در این مقاله به بررسی رفتار عمق موهو با استفاده از داده‌ های آنامولی جاذبه برمبنای روش پارکراولدنبرگ پرداخته می ‌شود. فرمولی که توسط oldenburg از طریق ادغام با روش parker موسوم به روش پارکراولدنبرگ در اینجا بازنویسی شده تا به روش تکراری معکوس تبدیل فوریه آنامولی جاذبه، نتیجه حاصل شود. از آنجایی که این روش بر اساس تبدیل سریع فوریه بنا نهاده شده است، بنابراین دارای سرعت بسیار بالایی است که می‌ توان از آن برای محاسبه‌ی مدل‌هایی با تعداد بسیار بالای نقاط بدون صرف زمان زیاد برای محاسبات استفاده کرد. همچنین در صورت استفاده از میدان ثقلی با کیفیت بالا می‌ توان به نتایج خوبی در این روند دست یافت. در این پژوهش آنامولی‌های جاذبه حاصل از مدل‌ های ژئوپتانسیلی egm08, egm96 و یکی از مدل‌های ژئوپتانسیل جهانی گوسمبنا (بر اساس داده‌های ثقل‌سنجی ماهواره جهانی goce تنها حاصل شده است) وعلاوه برآن از داده‌های ثقل‌سنجی زمینی تهیه شده توسط سازمان نقشه‌برداری در منطقه خراسان استفاده شده است. بوسیله این داده‌ها یک شبکه  سلولی به منظور تولید میدان ثقل و تخمین عمق موهو ایجاد شده است. بررسی نتایج حاصل از محاسبه عمق موهو در این منطقه نشان می‌دهد که مدل عمق موهو بدست آمده از داده‌های سازمان نقشه برداری نسبت به دیگر مدل ها اختلاف زیادی دارد که به دلیل تعداد محدود نقاط مشاهدات برای رسیدن به مدل درونیابی میدان ثقل است. اما از اختلاف نتایج عمق موهو حاصل از مدل‌egm08 نسبت به مدل‌های egm96 و مدل goce مقدار rms بترتیب 1.66 و 1.07 کیلومتر در عمق موهو بدست آمده است که این بهبود دقت را می ‌توان ناشی ازکیفیت و رزولوشن مدل‌های ژئوپتانسیلی دانست. همچنین در مقایسه نتایج حاصل از مدل goce با مدل egm96 مقدار rms برابر با 0.85 کیلومتر می‌ باشد که بدلیل نزدیکی و کیفیت دو میدان مورد استفاده نسبت به هم است.

Estimation of Moho depth based on the inversion of gravity anomaly data from different sourcesCase study: Khorasan region

Extended Abstract Introduction Estimation of the Moho depth and thickness of the crust using the gravity anomaly datais one of the basic researches in the geophysics and geology sciences.   Materials & Methods Based on many geophysical studies, the threedimensional thickness determination of the density variationinterface using the gravity anomaly is a common method. One practical instanceis the modeling of the crustal discontinuity like Mohorovicic discontinuity using the gravity anomalies. To analyze the anomalies associated with these crustal discontinuities, many techniques are used. Among the common methods generallyused for estimating the Moho depth and studying thecrustal structure arethe analysis of surface and body wavesof the earthquakesrecorded at theseismological stations, the analysis of postseismic waves, the gravity data inversion method and thermal analysis. In these cases, the inversion of the filtered gravity anomalies for determining the interface geometry of the density variations is one of the main goals. Different researches have proposed different methods for calculating the interface geometry of the density variationsbased on thegravity anomaly. Many of them approximate an irregular body with several cubic prismelements withconstant density. The overall gravity field of the bodyis calculated based on the sum of the gravity field effects of the prisms. Some methods such as Oldenburg (1974) have been developed based on the rewrite of Parker’s forward method (Parker, 1973). Based on the Parker’s method,the Fourier transform of the gravity anomaly is consideredas an outcomeof thesum of theFourier transforms of the createddepth powersrelated tothe gravity anomaly. Oldenburg shows that theParker’s formula can be rewritten to determine the geometry of the density interface from thegravity anomaly data. In this method, the Parker’s formula inversion is used to calculate the gravity anomaly created by an uneven layer of materials based on the Fourier series. Oldenburg rewrote this formula to calculate the interface depth of the density with undulating geometry using thegravity anomaly based on an iterative method (ParkerOldenburg method). Therefore, the topography ofthe densityinterface is estimatedbased onan iterative inversion method, which is repeated until an acceptable solution is obtained. According to the method (Oldenburg, 1974), the process is convergedin casethe depth of the interface is greater than zero and is not removed from the topography. Moreover, the range of theinterface variations should be less than the average depth of the interface. When a specific number of iterations is performed or the difference between two successful approximations is less than a specific value, the iterative procedure ends. In general, this gravity anomaly modeled by the inversion method should be very similar to the input gravity anomaly in the first stage. This paper investigates the Moho depth behavior using gravity anomaly data based on the ParkerOldenburg method. The formula rewritten by Oldenburg through integration with the Parker’smethod called the ParkerOldenburg method is used here to obtain the results by the iterative inversion method oftheFourier transform of the gravity anomaly. Since this method is based on the Fast Fourier Transform(FFT), it has a very high speed which can be used to compute models with a very high number of points without spending too much time on computation. Good results can also be achieved by using a highquality gravity field.   Results & Discussion In this study, the gravity anomalies derived from EGM08, EGM96 geopotential models and one of the GOCEbased global geopotential models (obtained only from the global satellite gravimetry data of GOCE), as well as those derived from terrestrial gravity data provided by the National Cartography Center (NCC) have been usedin Khorasan region. A  cell grid has been createdto generate the gravity field and estimate the Moho depth. Investigation of the results obtained from theMoho depth calculation in this region shows that the Moho depth model obtained from NCC data is very different from other models due to the limited number of observation points to reach the gravity field interpolation model. The difference of theMoho depth derived from the EGM08 model and the onederived from theEGM96 and GOCE models, gave 1.66 and 1.07 km for the RMS values, respectively. This accuracy improvement can be attributed to the quality and resolution of the geopotential models. Furthermore, comparing the results of the GOCE model with the EGM96 model, the RMS value is 0.85 km which is due to the close proximity of the two models’ qualities.   Conclusion In this paper, the Moho depth model has beenobtained based on the ParkerOldenburg method using the gravity anomaly data forKhorasan region. In this method,the Fourier transform ofthe gravity anomalies accelerates themodeling for a large number of points. On the other hand, the highquality of the models for the production of anomaly, results in the production of thehighly precise geometry of the density interface to a certain extent.