重力场模型低阶系数时变研究.PDF

第六章 重力場模型低階係數時變研究 地球重力場是地球質量分佈的表現，為地球的一種物理屬性，表徵地球內部、 表面或外部各點所受地球重力作用的空間 。但地球質量會因為地核振蕩、板塊運 動、海平面變化、大氣、地殼的冰期後回彈等現象不斷的改變和重新分佈 (Torge ， 1989) ，根據地球重力場的分佈，可以研究地球內部結構、地球形狀等，並進一步 尋找地球變化的動力學和物理學原因。發 生全球質量位移時會影響地球引力場球 諧函數模型，低階球諧係數表現出地球質心位置與扁率的變化。 6-1 CHAMP軌道資料求解低階重力球諧係數之時間變化 前一章提到 CHAMP軌道資料是由德國 Muenchen 技術大學的 Svehla D.博士 所提供，日期為 2002 年 3月 11日 -2003 年 3月 10日共一年的軌道資料，數據時 間隔為 30秒，原始資料為地固座標系統 (Conventional Inertial System ，CTS) 。 因考慮到硬體的效能和 GEODYN軟體的功能，所以篩選每月三天為一組，計 12 組的軌道資料(表 6-1) ，分別求取5 階、 10 階低階次的重力球諧係數。關於 GEODYN 步驟 GIIS 之 CHAMP衛星定軌參數設定見表 6-2 。 地球非球形引力擾動可分爲帶諧擾動和田諧擾動。其中，二階帶諧係數 J 的 2 量級爲 10−3 ，其他高階帶諧項影響都要比二階帶諧項小得多，所以在考慮地球非球 形引力的攝動影響時，取J 和 J 這兩項就足夠。 2 3 低階帶諧係數 J 可經由下列公式 (6-1)(Kaula ，1996)轉換得到，求得之二階帶 n 諧係數 J 可表達地球形狀為越趨近圓形或越趨近扁形的變化。 2 76 J n −Cn 0 − 2n +1Cn 0 (6-1) 式中 n ：階次(degree) C20 ：球諧係數 C20 ：完全正規化(Fully Normalized)球諧係數 求出 12 組的軌道資料，選取每一組的J 、J 減去平均值，得 J 和 J的時間 2 3 2 3 變化，圖 6-1 、圖6-2 為分別利用 5階和 10階重力球諧係數所解算之結果，其中可 看出兩個階次的 J 、J時間變化圖形相似，皆有趨近季節性的 變化，於夏季月份 2 3 曲線走勢向上升，但因求解資料時間較短變化不顯著。 表 6-1 CHAMP衛星軌道資料分組表 Data Number GPS Week MJD Start-End date(year-month-day) 1 1158 52351-52353 2002-03-18 to 2002-03-20 2 1162 52382-52384 2002-04-18 to 2002-04-20 3 1167 52412-52414 2002-05-18 to 2002-05-20 4 1171 52443-52445 2002-06-18 to 2002-06-20 5 1175 52473-52475 2002-07-18 to 2002-07-20