美女一级毛片免费看看,国产精品一级二级,色台湾色综合网站

其他數(shù)據(jù)論文 II 區(qū)論文（已發(fā)表） ? 版本 ZH4 Vol 7 (2) 2022

GRACE精密星間重力位差數(shù)據(jù)集（2002年4月至2016年7月）

A dataset of GRACE intersatellite geopotential differences from April, 2002 to July, 2016

鐘波，李瓊，陳劍利，羅志才

ZHONG Bo, LI Qiong, CHEN Jianli, LUO Zhicai

?>>

DOI:? 10.11922/11-6035.csd.2022.0010.zh
PID:? 21.86101.1/11-6035.csd.2022.0010.zh
:? 10.11922/sciencedb.j00001.00365
?

： 2022 - 01 - 22

： 2022 - 04 - 22

： 2022 - 03 - 10

： 2022 - 06 - 30

6829 48 0

摘要&關(guān)鍵詞

摘要：GRACE（Gravity Recovery and Climate Experiment）星間重力位差是一種沿衛(wèi)星軌道的原位觀測量，它具有明顯的地球物理含義。本研究基于改進的能量平衡方程和移去恢復(fù)技術(shù)RCR（remove-compute-restore)，利用GRACE Level-1B核心載荷數(shù)據(jù)（KBR星間距離變率、衛(wèi)星精密軌道、非保守力加速度計和衛(wèi)星姿態(tài)等數(shù)據(jù)）估計了2002年4月至2016年7月期間的精密星間重力位差數(shù)據(jù)，并通過數(shù)據(jù)質(zhì)量控制手段對估計結(jié)果進行了質(zhì)量標(biāo)記。這些數(shù)據(jù)按月進行存儲，共有155個月的數(shù)據(jù)文件，提供了GRACE雙星A和B的地心地固坐標(biāo)、扣除平均重力場（GGM05C模型）的星間重力位差和質(zhì)量標(biāo)記等數(shù)據(jù)。同時，將GRACE星間重力位差反演的時變重力場和地表質(zhì)量變化模型與官方發(fā)布模型進行比較，驗證了數(shù)據(jù)集的可靠性。本數(shù)據(jù)集進一步豐富了GRACE數(shù)據(jù)產(chǎn)品的種類，可為靜態(tài)和時變重力場建模、全球和局部地表質(zhì)量變化反演及相關(guān)地學(xué)應(yīng)用研究提供有力的數(shù)據(jù)支撐。

關(guān)鍵詞：GRACE；星間重力位差；改進的能量平衡方程；Level-1B數(shù)據(jù)；地表質(zhì)量變化

Abstract & Keywords

Abstract:?Gravity Recovery and Climate Experiment (GRACE) intersatellite geopotential differences are in-situ gravimetric observables with explicit geophysical interpretation. Based on the improved energy balance equation and remove-compute-restore (RCR) technique, we estimated the precise intersatellite geopotential differences from April, 2002 to July, 2016 from GRACE Level-1B data (including K-band ranging (KBR) range-rate, satellite orbits, accelerometer and attitude measurements), and marked the estimate results by means of data quality control method. The dataset includes 155 data files with monthly interval, and each data file includes the positions of GRACE satellites A and B in Earth-centered Earth-fixed (ECEF) coordinate system, intersatellite geopotential difference after deducting the mean gravity field model (i.e., GGM05C), and the corresponding quality flags. Furthermore, the reliability of the dataset is validated by comparing the time variable gravity field and surface mass change models retrieved by GRACE intersatellite geopotential differences with the official released models. This dataset further enriches the GRACE data products, and can provide data support for static and time-varying gravity field modeling, global and regional surface mass change inversion, as well as related studies in geoscience.

Keywords:?GRACE;?intersatellite geopotential differences;?improved energy balance equation;?Level-1B data;?surface mass variations

數(shù)據(jù)庫（集）基本信息簡介

數(shù)據(jù)庫（集）名稱	GRACE精密星間重力位差數(shù)據(jù)集（2002年4月至2016年7月）
數(shù)據(jù)作者	鐘波、李瓊、陳劍利、羅志才
數(shù)據(jù)通信作者	鐘波（bzhong@whu.edu.cn）
數(shù)據(jù)時間范圍	2002年4月–2016年7月
地理區(qū)域	89°S–89°N，180°W–180°E，除兩極空白以外的全球所有區(qū)域
采樣間隔	5 s
數(shù)據(jù)量	4.39 GB
數(shù)據(jù)格式	*.txt
數(shù)據(jù)服務(wù)系統(tǒng)網(wǎng)址	http://www.doi.org/10.11922/sciencedb.j00001.00365
基金項目	國家自然科學(xué)基金（41974015、41474019）
數(shù)據(jù)庫（集）組成	本數(shù)據(jù)集包括2002年4月至2016年7月期間按月存儲的GRACE星間擾動重力位差數(shù)據(jù)，扣除因衛(wèi)星儀器故障或其他因素導(dǎo)致數(shù)據(jù)缺失的17個月份以外，共有155個月的數(shù)據(jù)文件。每個數(shù)據(jù)文件提供了該月GRACE雙星A和B的地心地固坐標(biāo)、扣除平均重力場（GGM05C模型）的星間重力位差（即擾動位差：\({T}_{AB}{=T}_{B}-{T}_{A}\)）和質(zhì)量標(biāo)記等數(shù)據(jù)。每個文件以“GRACE_GPD_YYYYMM.txt”形式命名，其中“GRACE_GPD_”表示GRACE星間重力位差，“YYYY”表示年，“MM”表示月，文件的原始格式為.txt，壓縮后的格式為.zip。

Dataset Profile

Title	A dataset of GRACE intersatellite geopotential differences from April, 2002 to July, 2016
Data corresponding author	ZHONG Bo (bzhong@whu.edu.cn)
Data authors	ZHONG Bo, LI Qiong, CHEN Jianli, LUO Zhicai
Time range	April 2002 – July 2016
Geographical scope	89°S–89°N, 180°W–180°E, global regions except the polar gaps
Sampling rate	5 s
Data volume	4.39 GB
Data format	*.txt
Data service system	http://www.doi.org/10.11922/sciencedb.j00001.00365
Sources of funding	National Natural Science Foundation of China (41974015, 41474019)
Dataset composition	The dataset includes 155 data files of GRACE intersatellite geopotential differences with monthly interval from April, 2002 to July, 2016 except the 17 months whose data are not available due to the instrument failure or other influencing factors. Each data file contains the positions of GRACE satellites A and B in Earth-centered Earth-fixed (ECEF) coordinate system, intersatellite geopotential differences after deducting the mean gravity field model GGM05C (i.e., the disturbance potential difference T_AB=T_B-T_A), and the corresponding quality flags. Each data file is named in the form of “GRACE_GPD_YYYYMM.txt”, where “GRACE_GPD_”, “YYYY” and “MM” represents the GRACE intersatellite geopotential difference (GPD), the year and the month respectively. The original format of each data file is “.txt”, and the compressed format is “.zip”.

引言

2002年3月至2017年10月，美德聯(lián)合實施的GRACE（Gravity Recovery and Climate Experiment）衛(wèi)星重力測量任務(wù)搜集了約15年的衛(wèi)星跟蹤衛(wèi)星觀測數(shù)據(jù)，其核心載荷數(shù)據(jù)包括星載GPS（Global Positioning System）接收機觀測的衛(wèi)星精密軌道數(shù)據(jù)、KBR（K-band ranging）微波測距系統(tǒng)觀測的星間精密距離變率數(shù)據(jù)、星載加速度計觀測的衛(wèi)星非保守力數(shù)據(jù)和恒星敏感系統(tǒng)觀測的衛(wèi)星姿態(tài)數(shù)據(jù)等。GRACE重力衛(wèi)星相關(guān)研究機構(gòu)CSR（Center for Space Research）、JPL（Jet Propulsion Laboratory）和GFZ（Geoforschungszentrum）等利用這些數(shù)據(jù)反演了一系列的月間隔時變重力場球諧系數(shù)模型（Spherical Harmonic Solution, SH solution）和地表質(zhì)量變化模型（Mascon solution），并廣泛應(yīng)用于陸地水循環(huán)、冰川冰蓋質(zhì)量平衡、海平面變化和海底壓力變化及地震等全球變化研究領(lǐng)域^[1-3]。

GRACE任務(wù)是通過測定在同一低軌道上相距約220 km的兩顆衛(wèi)星之間的相對速率變化來感知地球重力場變化并推演地表質(zhì)量變化，其最直接的觀測量為KBR星間距離變率數(shù)據(jù)。但由于星間距離變率為幾何觀測量，其與地表質(zhì)量變化之間為復(fù)雜的非線性函數(shù)關(guān)系，并不能直接地用于地表質(zhì)量變化的反演與分析。目前，CSR、GFZ、JPL和GSFC（Goddard Space Flight Center）等機構(gòu)利用動力學(xué)積分法由GRACE星間距離變率數(shù)據(jù)及衛(wèi)星軌道、加速度計和衛(wèi)星姿態(tài)等輔助數(shù)據(jù)反演地表質(zhì)量變化模型的空間分辨率約為200–300 km、時間分辨率為個1個月^[4-6]。但由于動力學(xué)積分法的算法復(fù)雜、計算耗時，并且各種官方發(fā)布的地表質(zhì)量變化模型采用了固定的正則化約束方法進行求解，并不便于用戶采用不同的正則化方法對特定區(qū)域進行研究與分析^[7]。

GRACE星間重力位差（geopotential difference, GPD）是一種沿衛(wèi)星軌道的原位觀測量，它具有明顯的物理含義（與幾何觀測量的星間距離變率不同），并且與地表質(zhì)量變化之間為簡單的線性關(guān)系，因此可直接用于地表質(zhì)量變化的反演與分析^[8-11]。本研究基于改進的能量平衡方程^[11-13]，利用GRACE Level-1B（RL03）數(shù)據(jù)估計了2002.04–2016.07期間月間隔的共155個月（其余17個月數(shù)據(jù)缺失）的精密星間重力位差數(shù)據(jù)，可方便用于靜態(tài)和時變重力場建模、全球和區(qū)域地表質(zhì)量變化的反演與分析等研究。本研究提供的GRACE星間重力位差在本質(zhì)上是由核心載荷觀測數(shù)據(jù)計算導(dǎo)出的一種間接觀測量，因此可視為介于Level-1B（觀測值）和Level-2（應(yīng)用產(chǎn)品）之間的一種新的數(shù)據(jù)產(chǎn)品，可為大地測量學(xué)、地球物理學(xué)、水文學(xué)、冰川學(xué)、海洋學(xué)和地震學(xué)等相關(guān)地學(xué)研究提供基礎(chǔ)數(shù)據(jù)支撐，具有重要的科學(xué)意義和實用價值。

1 ? 數(shù)據(jù)采集和處理方法

1.1 ? 數(shù)據(jù)來源

本研究利用2002年4月至2016年7月的GRACE Level-1B（RL03）數(shù)據(jù)（ftp://isdcftp.gfz-potsdam.de/grace/Level-1B/JPL/）估計星間重力位差，這些數(shù)據(jù)具體包括KBR星間距離變率數(shù)據(jù)（KBR1B）、衛(wèi)星精密軌道數(shù)據(jù)（GNV1B）、非保守力加速度計數(shù)據(jù)（ACC1B）和衛(wèi)星姿態(tài)數(shù)據(jù)（SCA1B）等。其中，KBR1B和GNV1B數(shù)據(jù)的采樣間隔為5 s，而ACC1B和SCA1B數(shù)據(jù)的采樣間隔為1 s。為了使這些數(shù)據(jù)的時間間隔保持一致，本研究將ACC1B和SCA1B數(shù)據(jù)降采樣為5 s。同時，為了估計精密的GRACE星間重力位差數(shù)據(jù)，本研究參考了最新的GRACE Level-2（RL06）產(chǎn)品數(shù)據(jù)處理標(biāo)準(zhǔn)^[14]，最終采用的背景力模型如表1所示。其中，靜態(tài)重力場采用GGM05C模型計算，三體攝動采用JPL DE430星歷計算，海潮采用EOT11a模型計算，固體潮、大氣潮、固體極潮、海極潮和相對論效應(yīng)的計算與IERS 2010（International Earth Rotation Service Conventions 2010）一致，短期大氣和海洋非潮汐變化采用AOD1B（Atmosphere and Ocean De-aliasing Level-1B） RL06模型（ftp://isdcftp.gfz-potsdam.de/grace/Level-1B/GFZ/AOD/RL06/）計算。

表1 ? 背景模型信息

序號	背景模型	數(shù)據(jù)源	模型描述
1	靜態(tài)重力場	GGM05C	展開至360階次
2	三體攝動	JPL DE430	牛頓點質(zhì)量模型，包括日、月及其他行星
3	固體潮	IERS 2010	計算2，3，4階潮，包括與頻率相關(guān)項和無關(guān)項
4	海潮	EOT11a	展開至120階次，包括18個主波、238個次波
5	大氣潮	Ray and Ponte^[15]	S1和S2分潮波，展開至30階次
6	固體極潮	IERS 2010	C₂₁和S₂₁項，平均極采用線性模型
7	海極潮	Desai^[16]	展開至100階次
8	相對論效應(yīng)	IERS 2010	僅考慮日、月的影響
9	非潮汐大氣海洋去混頻	AOD1B RL06	展開至180階次，時間分辨率為3小時

1.2 ? 數(shù)據(jù)處理方法

為了精密估計GRACE星間重力位差，首先需要消除衛(wèi)星軌道數(shù)據(jù)和背景力模型誤差的影響。本研究將衛(wèi)星精密軌道數(shù)據(jù)（GNV1B）視為偽觀測，利用表1給出的背景力模型由動力學(xué)軌道積分方法重新估計純動力學(xué)軌道，并且在軌道積分過程中對非保守力加速度計數(shù)據(jù)（ACC1B）進行校準(zhǔn)。以每天一個弧段進行軌道積分計算，對于GRACE的每顆衛(wèi)星而言，估計的參數(shù)包括衛(wèi)星初始狀態(tài)參數(shù)（3個位置坐標(biāo)分量和3個速度坐標(biāo)分量）、加速度計尺度和偏差參數(shù)（每個坐標(biāo)分量估計1個尺度和1個偏差）。

基于改進的能量平衡方程，利用KBR星間距離變率數(shù)據(jù)（KBR1B）、重新估計的純動力學(xué)軌道、校準(zhǔn)后的非保守力加速度計數(shù)據(jù)及衛(wèi)星姿態(tài)數(shù)據(jù)（SCA1B）可直接計算GRACE星間重力位差。由于以上采用了已知背景力模型重新估計的純動力學(xué)軌道進行計算，估計的GRACE星間重力位差只包含了KBR星間距離變率觀測提供的時變重力場信息，有效避免了衛(wèi)星軌道和背景力模型的誤差影響^[11-12,17]。同時，為了消除KBR星間距離變率數(shù)據(jù)低頻系統(tǒng)誤差的影響，還需要采用與衛(wèi)星軌道周期有關(guān)經(jīng)驗參數(shù)^[8,10,12]對以上計算的GRACE星間重力位差進行擬合，以有效分離系統(tǒng)誤差的影響。

但是，引入的經(jīng)驗參數(shù)在分離系統(tǒng)誤差的同時，也會吸收GRACE星間重力位差包含的部分時變信號，特別是低頻的長波分量信號。為此，本研究基于移去恢復(fù)技術(shù)RCR（remove-compute-restore) ，首先從計算的GRACE星間重力位差中移去長波分量的貢獻（以CSR時變重力場模型的低階系數(shù)為參考），然后利用經(jīng)驗參數(shù)估計并分離系統(tǒng)誤差，最后再將移去的長波分量信號進行恢復(fù)。這樣便有效避免了經(jīng)驗參數(shù)對時變重力場信號的吸收，進而可得到高精度的全頻段的GRACE星間重力位差數(shù)據(jù)。圖1給出GRACE星間重力位差估計的數(shù)據(jù)處理流程，相關(guān)計算公式和算法可參見文獻[11]。

圖1 ? GRACE星間重力位差估計的數(shù)據(jù)處理流程

Figure 1 Data processing flow for estimating GRACE intersatellite geopotential differences

2 ? 數(shù)據(jù)樣本描述

本數(shù)據(jù)集包括2002年4月–2016年7月期間按月存儲的GRACE星間重力位差數(shù)據(jù)，扣除因衛(wèi)星儀器故障或其他因素導(dǎo)致數(shù)據(jù)缺失的17個月份以外，共有155個月的數(shù)據(jù)文件，總數(shù)據(jù)量為4.39GB。每個文件以“GRACE_GPD_YYYYMM.txt”形式命名，其中“GRACE_GPD_”表示GRACE星間重力位差，“YYYY”表示年，“MM”表示月。同時，為了節(jié)省存儲空間，每個文件被壓縮為“GRACE_GPD_YYYYMM.zip”的形式。每個文件都包括頭文件信息，用于說明平均重力場模型GGM05C的相關(guān)參數(shù)、該月可用的數(shù)據(jù)天數(shù)、每列數(shù)據(jù)的具體含義及單位等。頭文件結(jié)束后，數(shù)據(jù)的每一行提供了GRACE雙星A和B的地心地固坐標(biāo)3分量（單位為m）、扣除平均重力場（GGM05C）的星間重力位差（即擾動位差：\({T}_{AB}{=T}_{B}-{T}_{A}\)，單位為m²/s²）和質(zhì)量標(biāo)記（0或1）等數(shù)據(jù)，該數(shù)據(jù)的采樣間隔為5 s。圖2繪制了2005年3月、6月、9月和12月的GRACE星間重力位差空間分布，由于該數(shù)據(jù)已經(jīng)扣除了平均重力場的貢獻，因此可以直接反映地球重力場變化及地表質(zhì)量變化。對于陸地區(qū)域的地表質(zhì)量變化來說，其主要為陸地水儲量變化（如非洲、南美和中國等地區(qū)），可以看出該變化具有明顯的季節(jié)性特征。

圖2 ? 2005年3月、6月、9月和12月的GRACE星間重力位差空間分布圖

Figure 2 Spatial patterns of GRACE intersatellite geopotential differences in March, June, September and December in 2005

3 ? 數(shù)據(jù)質(zhì)量控制和評估

由于衛(wèi)星載荷觀測誤差、背景力模型誤差和數(shù)值計算近似誤差等的影響，GRACE星間重力位差的估計精度約為0.0020~0.0025 m²/s²^[8,12] 。為此，本研究在估計GRACE星間重力位差的過程中，以衛(wèi)星繞地球的旋轉(zhuǎn)周期（約為5400 s）為一個弧段，統(tǒng)計每個弧段估計星間重力位差的殘差標(biāo)準(zhǔn)差\(\sigma \)。若\(\sigma \)>0.0025 m²/s²，則認為該弧段估計GRACE星間重力位差的質(zhì)量不可靠，并將整個弧段估計結(jié)果的質(zhì)量標(biāo)簽標(biāo)記為“1”，反之則認為該弧段估計結(jié)果的質(zhì)量符合要求，將其質(zhì)量標(biāo)簽標(biāo)記為“0”。另外，對于每個弧段的GRACE星間重力位差數(shù)據(jù)，進一步采用3倍標(biāo)準(zhǔn)差準(zhǔn)則對粗差或異常點進行判斷，如果星間重力位差的殘差大于0.0075 m²/s²，則將其數(shù)據(jù)的質(zhì)量標(biāo)簽標(biāo)記為“1”。

由于GRACE衛(wèi)星儀器故障或其他因素的影響，通常會出現(xiàn)部分觀測數(shù)據(jù)缺失或質(zhì)量較差無法應(yīng)用的情況。例如，2002年4月–2016年7月期間共有155個月的觀測數(shù)據(jù)可用，有17個月的數(shù)據(jù)缺失。同時，在這可用的155個月份的數(shù)據(jù)中，部分月份的數(shù)據(jù)也會出現(xiàn)有若干天數(shù)據(jù)缺失的情況。因此，對于生成的每個月的GRACE星間重力位差數(shù)據(jù)，在數(shù)據(jù)頭文件中給出了該月可用數(shù)據(jù)的天數(shù)信息。

為了驗證GRACE星間重力位差的估計精度，將其用于反演時變重力場球諧系數(shù)模型^[12]和地表質(zhì)量變化模型^[11]，并與官方機構(gòu)發(fā)布模型進行比對。圖3為利用2005年2月和9月的GRACE星間重力位差數(shù)據(jù)反演的球諧系數(shù)模型（記為GPD SH）與官方發(fā)布的球諧系數(shù)模型（CSR RL06SH、JPL RL06SH和GFZ RL06SH）的位系數(shù)階均方根（Degree RMS）比較。圖4為利用2005年1月至2010年12月的GRACE星間重力位差數(shù)據(jù)反演的南美地區(qū)地表質(zhì)量變化（記為GPD Mascon）與官方發(fā)布的地表質(zhì)量變化模型（CSR RL06M、JPL RL06M和GSFC Mascon）計算得到的周年振幅、長期趨勢、均方根誤差（RMSE）和信噪比（SNR）的空間分布比較。圖3和圖4均顯示出GRACE星間重力位差反演的球諧系數(shù)模型和地表質(zhì)量變化模型與官方機構(gòu)發(fā)布模型的精度相當(dāng)，并且具有很好的一致性，驗證了GRACE星間重力位數(shù)據(jù)估計的可靠性。

圖3 ? 各種重力場球諧系數(shù)模型的位系數(shù)階RMS比較（2005年2月和9月）

圖4 ? 各種Mascon模型計算的周年振幅、長期趨勢、RMSE和SNR比較（2005–2010年）

Figure 4 Comparison of the annual amplitudes, long-term trends, RMSE and SNR from different Mascon models (from 2005 to 2010)

4 ? 數(shù)據(jù)價值

盡管GRACE星間重力位差是一種間接觀測量，但它比KBR測距系統(tǒng)提供的星間距離變率這種幾何觀測量更具有明顯的物理含義，可方便地應(yīng)用于時變重力場和地表質(zhì)量變化的反演與分析。由于GRACE星間重力位差是利用Level-1B核心載荷數(shù)據(jù)（KBR1B、GNV1B、ACC1B和SCA1B）估計得到，因此如何利用更為嚴(yán)密的估計模型并有效消除觀測數(shù)據(jù)誤差及背景力模型誤差的影響，是得到GRACE精密星間重力位差的關(guān)鍵環(huán)節(jié)。

本研究采用了更為嚴(yán)密的能量平衡方程進行計算，并基于移去恢復(fù)技術(shù)RCR有效避免了文獻^[8-10]中經(jīng)驗參數(shù)對長波部分時變重力場信號的吸收，得到了全頻段的GRACE精密星間重力位差數(shù)據(jù)^[11]。并且，通過較為嚴(yán)格的數(shù)據(jù)質(zhì)量控制手段，對估計的GRACE星間重力位差數(shù)據(jù)進行了質(zhì)量標(biāo)記。本數(shù)據(jù)集的建立進一步豐富了GRACE數(shù)據(jù)產(chǎn)品的種類，可為靜態(tài)和時變重力場建模、全球和局部地表質(zhì)量變化反演及相關(guān)地學(xué)應(yīng)用研究提供有力的數(shù)據(jù)支撐。

[1]

許厚澤, 陸洋, 鐘敏, 等．衛(wèi)星重力測量及其在地球物理環(huán)境變化監(jiān)測中的應(yīng)用[J]．中國科學(xué): 地球科學(xué), 2012, 42(06): 843–853. DOI:10.1007/s11783-011-0280-z. [XU H Z, LU Y, ZHONG M, et al. Satellite gravity and its application to monitoring geophysical environmental change[J]. Scientia Sinica (Terrae), 2012, 42(06): 843–853. DOI:10.1007/s11783-011-0280-z.]

亚洲欧美日韩精品久久久,夜色福利一区二区三区,视频日韩p影院永久免费,成人免费视频国产