نوع مقاله : مقاله پژوهشی

نویسندگان

• روح اله کریمی ¹
• علیرضا آزموده اردلان ²
• سیاوش یوسفی ³

¹ استادیار گروه ژئودزی و مهندسی نقشه برداری، دانشگاه تفرش، تفرش،ایران

² استاد، گروه مهندسی نقشه برداری، قطب علمی مهندسی نقشه برداری در مقابله با سوانح طبیعی، پردیس دانشکده‌های فنی، دانشگاه تهران، ایران

³ فارغ التحصیل کارشناسی ارشد دانشکده مهندسی نقشه‌برداری و اطلاعات مکانی، پردیس دانشکده های فنی، دانشگاه تهران، تهران، ایران

چکیده

در حال حاضر بالاترین قدرت تفکیک مکانی مدل‌های ژئوپتانسیلی جهانی حدود 5 دقیقه می‌باشد، در حالی که مدل‌های توپوگرافی با قدرت تفکیک مکانی حدود 3 ثانیه و بالاتر در دسترس است. یکی از روش‌هایی که برای افزایش دقت مدل‌های ژئوپتانسیلی جهانی در تولید تابعک‌های مختلف میدان ثقل مورد استفاده قرار می‌گیرد، تلفیق این مدل‌ها با مدل‌های توپوگرافی با قدرت تفکیک مکانی بالاتر از مدل ژئوپتانسیلی است. در این مقاله هدف ارزیابی مؤلفه‌های زاویه انحراف قائم حاصل از تلفیق مدل ژئوپتانسیلی جهانی و مدل توپوگرافی با قدرت تفکیک مکانی بالا در ایران می‌باشد. تحقیق حاضر، از مدل EGM2008 با قدرت تقکیک مکانی حدود 5 دقیقه به عنوان مدل ژئوپتانسیلی جهانی، از مدل SRTM با قدرت تفکیک مکانی 3 ثانیه به عنوان مدل توپوگرافی و از مدل DTM2006 برحسب هارمونیک‌های کروی تا درجه 2190 به عنوان سطح هموار مرجع برای تولید مدل توپوگرافی باقیمانده (RTM) استفاده نموده است. روش تحقیق به این صورت است که ابتدا با استفاده از مدل جهانی، مؤلفه‌های زاویه انحراف قائم در 10 ایستگاه لاپلاس ایران محاسبه شده و سپس با استفاده از مدل توپوگرافی باقیمانده تصحیحی برای این مؤلفه‌ها بدست می‌آید. در پایان مؤلفه‌های زاویه انحراف قائم محاسبه شده توسط مدل جهانی به تنهایی و تلفیق مدل جهانی و مدل توپوگرافی باقیمانده با مؤلفه‌های زاویه انحراف قائم حاصل از مشاهدات نجومی و ژئودتیکی در 10 ایستگاه لاپلاس مقایسه می‌شوند. نتایج این مقایسه‌ها حاکی از آن است که تلفیق مدل جهانی EGM2008 و RTM باعث بهبود حدود 15% در مؤلفه شمالی-جنوبی(𝜉)  و 4/1% بهبود در مؤلفه شرقی-غربی(𝜂)در منطقه تست ایران می‌گردد.همچنین ارزیابی‌ها نشان می‌دهند که خطای نسبی در محاسبه مؤلفه𝜉 با استفاده از تلفیق مدل EGM2008 و RTM حدود 6% و در محاسبه مؤلفه 𝜂 حدود 37% است.

کلیدواژه‌ها

• ایران
• زاویه انحراف قائم
• مدل ژئوپتانسیلی جهانی
• مدل توپوگرافی باقیمانده
• SRTM
• EGM2008
• DTM2006

20.1001.1.25883860.1399.29.113.2.4

عنوان مقاله [English]

Assessment of vertical deflection derived from a combination of Global Geopotential Model and Digital Terrain Model in Iran

نویسندگان [English]

• Roohollah Karimi ¹
• Ali Reza Azmoude Ardalan ²
• Siavash Yousefi ³

¹ Assistant Professor, Department of geodesy and surveying engineering, Tafresh University, Tafresh, Iran

² Professor, Dept. of Geomatics, Center of Excellence in Geomatics Eng. and Disaster Prevention, Faculty of Eng, University of Tehran, Iran

³ M.Sc. graduate, School of Surveying and Geospatial Engineering, College of Engineering, University of Tehran, Tehran, Iran

چکیده [English]

Introduction
Components of verticaldeflection, i.e., North-South component  and East-West component ,are used for accurate determination of geoid or quasigeoid. Moreover, vertical deflection components area useful source for determination of variations in subsurface density and geophysical interpretations. Generally, there are two definitions for verticaldeflection. According to Helmert definition, vertical deflection at any given pointis the angle between the actualgravity vector (actual plumb line) and a line that is normal to the reference ellipsoid(a straight line perpendicular to the surface of reference ellipsoid). Another definition of vertical deflection is proposed by Molodensky. According this definition, vertical deflection at any given point is the angle between actualgravity vector and normal gravity vector (normal plumb line). Some relations have been introduced to convert Molodensky vertical deflection to Helmert vertical deflection. Helmert vertical deflection is estimated using astrogeodetic observations (combination of astronomical and geodetic observations).
Presently, global geopotential models (GGMs) have been expanded to the degree of2190, which is equivalenttoabout 5-min spatial resolution. Vertical deflectionat any point on the Earth can be calculated using the GGM. The resulting vertical deflection is consistent with Molodensky definition.Unfortunately, accuracy of GGMs is not sufficient for estimation of verticaldeflection.In other words, since GGMs are expanded up to a limited degree due to their resolution, omission error(or truncation error) occurs in computation of the earth’s various gravity field functionals, such as the geoidal height and verticaldeflection. Combining GGM with a digital terrain model (DTM) is a method used to reduce omission error.It should be noted that DTM has a higher spatial resolution as compared to GGM.In this method, the omitted signals of GGM can be modeled using residual terrain model (RTM) derived from subtracting high resolution DTM from a reference smooth surface. The reference smooth surface is obtained from eitherapplying average operator to DTM or expanding global topography into spherical harmonics. Fortunately, DTMs with spatial resolution of 3seconds or more,and reference smooth surface based on 2190 degree spherical harmonics are publicly available.
The present study seeks to assess vertical deflectionderived from a combination of GGM and DTM in Iran. Previously, Jekeli(1999) has studied EGM96 geopotential model with the aim of computingvertical deflection in the USA. Hirt(2010) and Hirt et al. (2010a) have assessed vertical deflection in Europe and the Alps using a combination of EGM2008 and RTM models.In Iran, GO_CONS_GCF_2_TIM_R4, a GOCE-only model, and EGM2008 geopotential model have been used toobtain vertical deflection and the results have been evaluated byKiamehr and Chavoshi-Nezhad(2014).
 
Materials & Methods
To implement the present study,a EGM2008 model with a spatial resolution of about 5-min is selected asGGM and a SRTM model with 3-sec spatial resolution is considered as DTM. To obtain RTM, DTM2006 model based on2190 degree spherical harmonicsis selected as the reference smooth surface.To compute the residual topography effect, prism method was used in an ellipsoidalmulti-cylindrical equal-area map projection system. First, we compute vertical deflectionusing EGM2008 model. It is also calculated using a combination of EGM2008 model and RTM(EGM2008/RTM method). In the next step, vertical deflection derived from the first method (EGM2008 model) and the second one (combination of EGM2008 model and RTM) are compared with vertical deflectionderived from astrogeodetic observations in 10 available Laplace stations in Iran.
 
Results & Discussion
Results indicate that there is a 1.2sec difference between North-South component of vertical deflection (i.e.) obtained from EGM2008 model and astrogeodetic observations.With RTM, this will reach 1 sec, which shows a 15% improvement. Moreover, there is a5.7secdifference between East-West component of vertical deflection () obtained from EGM2008 model and astrogeodetic observations, while this value will reach 5.6sec using RTM. Improvement in East-West component () is1.4%, which is smaller than the improvement of North-South component (). Based on the computations, we found that values of  and  in the Laplace stations canreach 17sec (RMS=7sec) and 15sec (RMS=8sec), respectively. Therefore, it is concluded that the relative error ofNorth-South component ()computation using EGM2008/RTM method is about 6% and the relative error ofEast-West component ()computation is about 37%.
 
Conclusion
The present research has studied the RTM effect on the improvement of GGM used for the determination of vertical deflectionin Iran. To performthe study, EGM2008 model with around 5-min spatial resolution was selected as GGM. RTM is also derived from subtracting the DTM2006 model (based on2190 degree spherical harmonics)from the 3-sec spatial resolutionSRTM model. Numerical findings indicate that a combination of RTM and GGM can improve the results of vertical deflectioncomputation, as compared to the results obtained from GGM-only approach. The improvement in North-South component of vertical deflection () is about15%and East-West component of the vertical deflection () undergoes about 1.4% improvement. In general, EGM2008 model and its combination with RTM have been more successful in the computation of  component as compared to computationin the geographical region of Iran. There is no clear explanation for this difference, but it can be due to errors in theastronomical or geodetic observations oflongitude in Laplace stations.

کلیدواژه‌ها [English]

• Iran
• Deflection of the vertical
• Global Geopotential Model
• Residual Terrain Model
• EGM2008
• SRTM
• DTM2006