GRACE Monthly Mass Grids - Global Mascons 简介与Notebook示例¶
GRCTellus（GRCTellus）月度质量网格海洋数据集。提供了相对2004-2010年平均基线的月度引力异常值。该数据集所包含的数据是以 "等水厚度 "为单位。
GRACE Tellus（GRCTellus）月度质量数据集由JPL基于GRACE观测Level -1 数据生产而来。该数据集利用空间和时间上的先验约束，以等面积3°x3°球冠质量浓度（mascon）函数来估计全球月度重力场，以减少测量误差影响。与传统的球面谐波解决方案相比，这使得mascon场的信噪比更好。
简介
什么是 "等效水厚度"？
观测到的重力月变化是由质量月变化引起的。可以认为，质量变化集中在地球表面附近的一个非常薄的水厚度变化层中（即最多几千米厚的一层）。实际上，大部分月重力变化都是由水文水库储水量的变化、海洋、大气和陆地冰块的移动以及这些地球系统组成部分之间的质量交换引起的。它们的垂直范围是以厘米为单位的等效水厚度，远小于地球半径或以公里为单位的变化水平尺度。有些重力变化是由 "固体 "地球的质量重新分布引起的，如大地震后的重力变化，或冰川等静力调整引起的重力变化；在这些情况下，用 "等效水厚度 "来解释重力变化是不正确的，尽管仍有可能计算出这一数量（即通过消除固体地球的影响来分离与水有关的质量变化）。
在处理过程中，使用 ECMWF 大气压力场去除大气质量，因此 GRACE Tellus 地表质量网格不包含格陵兰岛和南极洲等陆地或大陆冰区的大气质量变化（当然，ECMWF 的误差除外）。海洋模式用于去除高频率（6 小时至次月级）风和压力驱动的海洋运动，否则这些运动可能会异化为月度重力解。如果模型是完美的，所得到的重力场就不能反映海洋的变化。为了在海洋上使用这些结果，GRACE Tellus 海洋底层压力场又增加了月平均海洋模型网格。这也是我们分别提供海洋和陆地网格的原因之一。其他原因包括从海洋网格中去除陆地信号、陆地网格的缩放等。
这里提供的数据是等效水厚度相对于时间平均基线的变化。基本方法见 Wahr 等人（1998 年）的解释。陆地和海洋网格采用不同的滤波器进行处理，这些滤波器可在保留真实地球物理信号的同时，以最佳方式滤除噪声。
球谐波数据版本
GRACE 是一项首创性任务，因此数据处理的修订比更成熟的卫星测量更为频繁，这并不奇怪。有三个中心是 GRACE 地面系统的一部分，负责生成 Level-2 数据（球谐波场）： CSR（德克萨斯大学/空间研究中心）、GFZ（波茨坦地质研究中心）和 JPL（喷气推进实验室）。它们的输出包括重力场和用于计算重力场的迭加场的球谐波系数。
我们在此提供的数据以这些数据中心生成的数据（GRACE 项目的 "官方 "版本）为基础，并有额外的后处理步骤，在下面链接的页面中进行了总结。其他数据中心也提供 GRACE 衍生的质量网格，可在此处找到。
目前，我们提供来自最新 GRACE 重力场的 GRACE-Tellus 数据： 来自 CSR、JPL 和 GFZ 的第 05 版。多次测试表明，RL05 比之前的 RL04 有所改进。RL05 的球面谐波被用作我们后处理步骤的输入。我们将 RL05 的 OCEAN 版本称为 "RL05.DSTvDPC1401"，将 RL05 的 LAND 版本称为 "RL05.DSTvSCS1401"。后期处理已针对 RL05 的球面谐波进行了优化。请注意，这些版本号与我们在 RL05 数据可用后不久制作的网格的版本号不同。
相应的数据链接
Information and data for OCEAN grids
Information and data for LAND grids
数据版本
除了用球谐波求解重力变化外，还有一种方法是使用质量浓度块（或 "mascons"）作为本机基函数。使用质量浓度块作为基函数的主要优点是，每个质量浓度块都有一个特定的已知地球物理位置（与球谐波系数不同，后者没有特定的定位信息）。我们可以利用这一便利特性，在数据反演过程中指定先验信息（约束条件），从内部消除重力解算中的相关误差。因此，与 RL05 无约束球面谐波解不同，有约束的 mascon 解通常不需要去条纹或平滑处理。mascon 方法还能更好地分离陆地和海洋信号。
度数 2 / 阶数 0 系数（扁平性）
度，阶）（2，0）的球谐波系数有时与卫星激光测距（SLR）的系数不一致，后者更为可靠。此处提供的网格中的 (2,0) 系数由卫星激光测距的系数代替。更多信息。
度 1 系数（Geocenter）
GRACE 无法检索度 1 的球谐波系数，这些系数与地球地心相对于地球固定参考框架的位置成比例。这里使用的是根据 Swenson 等人（2008 年）对这些系数的估计值。更多信息。
地震
大地震，例如 2004 年 12 月 24 日在印度尼西亚苏门答腊岛北部西海岸附近、2005 年 3 月 25 日在苏门答腊岛北部、2007 年 9 月 12 日在苏门答腊岛南部、2010 年 2 月 27 日在智利毛勒近海以及 2011 年 3 月 11 日在日本东北近海本州东海岸附近发生的地震，会导致地球岩石圈发生足够的位移，从而使地球重力场发生变化，全球资源观测系统能够对其进行测量。上述所有地震的震级都在 8.5 级或以上（https://earthquake.usgs.gov/earthquakes/browse/largest-world.php）。将它们引起的重力变化解释为 "等效水厚度 "的变化是不正确的。对于冰川等静力调整，我们有足够好的模型来 "校正" GRACE 数据，但与此不同的是，我们目前还没有校正此处提供的数据，以消除大地震造成的信号。用户应警惕这些地震附近的信号。使用海洋模型的 EOFs 重建 GRACE 信号所计算出的海洋网格已经过滤掉了大部分地震信号。用户可以按照 de Linage 等人（2009 年）的方法大致去除地震信号；见他们的公式 3。
去条纹和空间平滑
经过上述步骤后，数据中会出现一个误差源，其明显特征是 N-S 条纹。Swenson 和 Wahr（2006 年）观察到了与条纹相关的球谐波系数的特殊性质，并设计了一类滤波器来消除这一问题。
GRACE 卫星的飞行高度超过 400 公里。重力场随高度而减弱，短波长比长波长衰减得更厉害。因此，有必要对短波长进行平滑处理，以恢复导致 GRACE 在其高度上看到重力场的地球表面质量集。为了减少这种噪声源，这里使用了空间平均平滑器（此处为高斯）。
精度较低的月份
用户需要注意的是，当轨道接近精确重复时，月网格的误差较大。这些月份包括 2004 年 7 月至 12 月，以及 2015 年 1 月和 2 月。造成误差较大的另一个原因是某个月份缺乏数据。
代码
import aie
aie.Authenticate()
aie.Initialize()
#指定检索数据集，可设置检索的时间范围
dataset = aie.ImageCollection('NASA_GRACE_MASS_GRIDS_MASCON')
imgs = dataset.select(['lwe_thickness']);
map = aie.Map(
center=[121.02,31.17],
height=800,
zoom=1
)
vis_params = {
'bands': ['lwe_thickness'],
'min': -25,
'max': 25,
"palette":[
'#ffffff','#d0d0d0','#ababab','#ffb93d','#3977ff'
]
}
map.addLayer(
imgs,
vis_params,
'lwe_thickness',
bounds=imgs.getBounds()
)
map结果
参考文献
Chambers, D.P. and J.A. Bonin: Evaluation of Release 05 time-variable gravity coefficients over the ocean. Ocean Science 8, 859-868, 2012. www.ocean-sci.net/8/859/2012.
Chambers, D. P: Observing seasonal steric sea level variations with GRACE and satellite altimetry, J. Geophys. Res., 111 (C3), C03010, 10.1029/2005JC002914, 2006.
Cheng, M. and Tapley, B.D.: Variations in the Earth's oblateness during the past 28 years, J. Geophys Res v109, B9, 2004.
de Linage C., L. Rivera, J. Hinderer, J.-P. Boy, Y. Rogister, S. Lambotte and R. Biancale: Separation of coseismic and postseismic gravity changes for the 2004 Sumatra–Andaman earthquake from 4.6 yr of GRACE observations and modelling of the coseismic change by normal-modes summation. Geophys. J. Int. (2009) 176, 695–714, doi: 10.1111/j.1365-246X.2008.04025.x.
Swenson, S. C. and J. Wahr, Post-processing removal of correlated errors in GRACE data, Geophys. Res. Lett., 33, L08402, doi:10.1029/2005GL025285, 2006.
Swenson S.C , D. P. Chambers, and J. Wahr: Estimating geocenter variations from a combination of GRACE and ocean model output. J Geophys. Res.-Solid Earth, Vol 113, Issue: B8, Article B08410. 2008.
Wahr, J., M. Molenaar, and F. Bryan, Time-variability of the Earth's gravity field: Hydrological and oceanic effects and their possible detection using GRACE, J. Geophys. Res., 103, 32,20530,229, 1998.
3172
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
