تلفيق مشاهدات GPS و ارتفاع سنجي ماهواره اي به منظور مدلسازي محلي يونسفر در ايران

Combination of GPS and Satellite Altimetry Observations for Local Ionosphere Modeling Over Iran

از مشاهدات غني و متراكم ايستگاه ­هاي دائمي GPS كه اغلب محدود به خشكي­ ها هستند، مي­توان براي مطالعۀ يونسفر و توليد نقشه‍ اي يونسفر استفاده نمود. اين نقشه­ ها نمايشگر مقادير محتواي مجموع الكتروني در راستاي قائم (VTEC) بوده و كاربرد خاص آنها در موقعيت­ يابي دقيق با استفاده از مشاهدات گيرنده ­هاي تك فركانس است. نقشه­ هاي جهاني يونسفر (GIM) كه توسط مراكز مختلف IGS توليد مي­شوند، همگي نمايشگر مقادير VTEC به ازاي هر دو ساعت، با رزولوشن مكاني °5/2 در راستاي عرض جغرافيايي و °5 در راستاي طول جغرافيايي و يا به صورت تابعي از ضرايب هارمونيك كروي تا درجه و مرتبۀ 15 و با رزولوشن زماني دو ساعت هستند. كمبود ايستگاههاي دائمي GPS در درياها منجر به افت صحت و دقت اين مدلها در آن مناطق شده و استفاده از مشاهدات دو فركانس ارتفاع سنجي ماهواره­اي به همراه مشاهدات GPS مي­تواند به عنوان راهكاري براي بهبود اين مدلها در اين مناطق محسوب گردد. در اين مقاله به منظور توليد نقشۀ محلي يونسفر ايران (LIM)، مشاهدات 21 ايستگاه دائمي GPS در روز 107 از سال 2014 (متناظر با بيشينه فعاليت خورشيدي) و توابع پايۀ بي-اسپلاين دو بعدي به عنوان توابع پايۀ مدلسازي استفاده شده­اند. مقادير مختلفي براي انتخاب سطح بهينۀ توابع پايۀ بي-اسپلاين در دو بعد مورد آزمون قرار گرفت و در نهايت، مدلسازي با سطح 1 و 2 به ترتيب در راستاي طول جغرافيائي و عرض جغرافيائي انجام شد. نتايج بدست­آمده، اختلافات عمده ­اي را نسبت به GIM نشان مي­دهند. در ادامه، مشاهدات GPS به همراه مشاهدات ارتفاع سنجي ماهواره­اي مربوط به ماهوارۀ Jason-2 در فرآيند مدلسازي محلي يونسفر در نظر گرفته شدند و براي وزن دهي نسبي بين مشاهدات دو دسته، از تكنيك برآورد مولفه­ هاي واريانس كمترين مربعات (LS-VCE) استفاده گرديد. در اين قسمت علاوه بر ضرايب مدل يونسفر محلي و باياسهاي كد تفاضلي (DCB) گيرنده ­ها، باياس ديگري بين مشاهدات GPS و ارتفاع سنجي ماهواره­اي نيز برآرود شده و با نتايج مرحلۀ قبلي مقايسه شدند. نتايج مقايسه نشان مي­دهند كه در صورت استفاده از مشاهدات ارتفاع سنجي ماهواره­اي به همراه مشاهدات GPS، دقت مدلسازي در سطح درياها بهبود خواهند يافت.

Ionosphere is the upper part of the atmosphere that extends from 80 to 1200 km above the Earth’s surface. The existing free electrones and ions in the ionosphere layer affect the signal propogation speed such as satellite positioning and satellite altimetry signals. Regardless of the fact that Dual frequency measurments can remove ionospheric delay effect, dual frequency observations of the permanent GPS stations can also be utilized to produce the ionosphere maps including the vertical total electron content (VTEC) values. For instance, International GNSS service (IGS) sub-centers produce daily global ionosphere maps (GIMs) using the GNSS data. The spatial resolution of GIMs in the latitude and longitude directions is 2.5 degree and 5.0 degree, respectively, and their temporal resolution is 2 hours. One of the IGS sub-centers, namely CODE produces the GIMs based on the spherical harmonic basis functions up to the degree and order 15. The aim of this research is to develop a local inosphere model based on the B-spline basis functions using the combined GPS and satellite altimetry observations over Iran. Accordingly, the potentiality of the B-spline basis functions for local inosphere modeling was studied at first. For this purpose, a local ionosphere model (LIM) was produced based on observation data from 16 Iranian permanent GPS stations and 5 IGS ones and B-spline basis functions. My assumptions in this modeling are as follows: first, the ionosphere is a thin shell that is located on 450 km above the Earth’s surface, second, the smoothed code station observations obtained by Bernese 5.0 software is considered as observation vector. Third, the weight matrix elements are proportional to the satellite elevation angle. Forth, the differential code biases (DCBs) for all satellites which are obtained from IGS precise products, are considered as known parameters in the equations. and the last assumption was that a simple cosine mapping function was used to convert the slant total electron content (STEC) to the VTEC. As a result, the comparison between the LIM and the GIM showed that the B-spline basis functions were more efficient than the spherical harmonic ones for local ionosphere modeling. Following the first result, a new LIM, which is based on the B-spline basis functions, was produced by integration of permanent GPS station and Jason-2 satellite altimetry observations. The GPS and satellite altimetry observations were chosen from day 107 of year 2014, according to the latest maximum solar activity. The weight matrix of the GPS and satellite altimetry observations were determined based on the least-square varince component estimation (LS-VCE) method. The results showed that the local inosphere model derived from combination of the GPS and satelite altimetry observations were more accurate than the local inosphere model derived from the GPS observations only, this is due to the fact that the dual-frequency radar altimetry data are the main source of the ionospheric observations at sea, where there is no GPS permanent station, and can be used to improve the GIMs and LIMs. Finally As by-products, the DCB values for the permanent GPS recievers and the bias term between the GPS and satellite altimetry observations were determined.