0 引言

地下水是重要淡水资源.地下水储量变化常规监测方法包括压力计测量与地面监测网测量.压力计测量覆盖点到局部空间尺度，对实现大尺度地下水位变化监测存在困难.地面监测网测量根据研究区给水度，依靠实测水位数据可估计地下水储量在中长时空尺度上的动态变化特征(Rodell and Famiglietti, 2001)，但这种方式需在局部测量之间插值或非采样区域外推计算，并严格依赖样点测量的精度与可靠性，尤其在监测网动态维护与水位统测方面消耗大量人力物力(Tiwari et al., 2009; 罗志才等, 2012).除了上述两种观测手段之外，地下水文模型也可用于估计地下水储量变化，但在模型刻画及数据获取方面存在很大的局限性(Tiwari et al., 2009; Famiglietti et al., 2011).

重力反演与气候实验卫星(Gravity Recovery and Climate Experiment, GRACE)计划的出现，为监测地球表层质量变化提供了新手段(Tapley et al., 2004; Wahr et al., 2004).在较短时间尺度上，地球重力场主要敏感于表层大气、海洋及陆地水之间的质量交换，而陆地时变重力场主要由陆地水质量变化引起(Tapley et al., 2004).目前GRACE卫星重力技术在陆地水储量变化研究等方面得到了广泛应用(Rodell et al., 2001; Swenson and Wahr, 2002; Swenson et al., 2003; Wahr et al., 2004; 胡小工等, 2006; 罗志才等, 2012).

陆地水储量变化信号中混叠GRACE重力场模型的高阶项噪声(Wahr et al., 2006)，以及特定阶数之间的相关误差(Swenson and Wahr, 2006).球谐系数高阶次截断与滤波方法使用泄露了研究区水质量变化的真实信号.目前已有研究针对GRACE滤波后的结果进行尺度因子标定，以减弱信号的损失程度(Landerer and Swenson, 2012).Klees等(2006)提出了依靠水文模型原始格网时间序列的多时刻尺度因子法，该方法对水文模型的依赖程度较高；冯伟等(2012)则进一步研究了GRACE提取亚马逊流域陆地水变化的单一尺度因子法，该方法可对GRACE结果误差进行有效改正.

昆明的江女士为挂专家号就被折腾得不轻。江女士的儿子左脚患上甲沟炎，为此耽误了不少学习时间。她听朋友说北京某大医院的专家治疗甲沟炎特别有效，就带儿子来到北京。

陆地水变化包括地下水，相应的时空变化可利用GRACE结合GLDAS、CPC、CLM4.0或CLM4.5等地表水文模型进行最新的定量估算(Tiwari et al., 2009; 冯伟等, 2017).对于地下水超采地区，如印度北部、加利福尼亚中谷、巴基斯坦及华北平原等地，GRACE技术可实现地下水储量的变化监测，监测结果可与常规监测手段及地下水文模型进行比对，结果表明两者一致性较好(Tiwari et al., 2009; Famiglietti et al., 2011; Iqbal et al., 2016; 冯伟等, 2017).

关中地区对地下水资源开发利用频繁，已引发局部多处地面沉降(陶虹等, 2013).

监测地下水动态变化对实现该地区水资源可持续利用具有重要意义.目前针对关中地区地下水的相关研究主要集中在地下水位变化规律及化学性质演变等方面(Yuan et al., 2011; Ma et al., 2012; 陶虹等, 2013; 贺旭波等, 2015; Duan et al., 2016).关于该地区地下水储量变化特征的研究则相对较少.

本文利用GRACE时变重力场及单一尺度因子法反演了关中地区地下水储量变化，与地下水实测资料及全球地下水文模型进行对比分析，并探讨降雨对地下水储量变化的影响.

1 数据与方法

1.1 GRACE数据后处理

本文采用美国德克萨斯大学空间研究中心(Center for Space Research, University of Texas at Austin, CSR)提供的GRACE level-2(RL05)月重力场模型数据，该数据为正则化球谐系数，扣除了固体潮、海潮、固体极潮、非潮汐大气和海洋影响，以及由于日月等其他行星引起的重力扰动(Bettadpur et al., 2007).GRACE时变重力场反演的陆地水储量变化以等效水高变化Δhw(φ,λ)表示为(Jekeli et al., 1981; Wahr et al., 1998)：

(1)

A股下跌主要原因，一方面是受中美贸易摩擦的扰动，另一方面自身内部问题也比较突出。首先，2017年开始整个金融大环境就是防风险去杠杆。防风险的重要工作是整治金融乱象，2017年以来，先后限制银行理财、规范保险资金以及券商资管入市，打压游资炒作，这些都是政治金融乱象的措施。进入2018年在防风险的基础上，去杠杆成为首要任务，去杠杆会导致流动性紧缩，进而影响增量资金入市，很明显观察到，今年两市日均成交额基本维持在3000-5000亿的地量水平。

除了空间分布总体上一致之外，研究区地下水储量变化整体均呈现长期减少趋势，如图4所示，2003—2013年GRACE估计结果-0.37 cm·a-1与WGHM结果-0.35 cm·a-1十分吻合.2003—2008年期间，GRACE监测到研究区地下水整体逐年上升，这与陶虹等(2013)给出的同期水位埋深年值数据变化结果较为吻合，该年值数据经线性拟合后也表现为地下水的小幅回升.但在WGHM估计的地下水变化中难以体现，这可能由于WGHM地下水模型覆盖全球空间尺度，对区域整体部分时段的结果估算精度有待进一步改善.

+ΔSn msinm λ]·

改造前，制硫溶剂再生（小再生）单元使用的溶剂（MDEA）未并入溶剂再生（大再生）单元系统，而改造后，小再生溶剂并入大再生系统。改造前后溶剂再生单元蒸汽单耗如图5所示。改造前的溶剂再生（大再生）单元蒸汽单耗（107kg蒸汽/t溶剂）比改造后蒸汽单耗（100kg蒸汽/t溶剂）多消耗7kg蒸汽/t溶剂，所以改造后，蒸汽单耗有少量下降，这与装置负荷有一定的关系，相当于提高了大再生的负荷，降低了装置能耗。

经过去相关误差改正后的球谐系数为：

(2)

式(2)中，逐阶、次的滤波核函数Wn、Wm通过递归计算得到，即

"You have a chance to the top! Corinne told me..." Yes, he is the first Chinese reach the top, and is also the first Chinese Master Sommelier in the world.

且

关于汽车轻量化的定义，不能仅以减轻质量多少来衡量，同时需要考虑设计车身的尺寸以及功能.主要可以从以下3个方面进行论述：①满足使用功能的汽车通过减轻质量实现轻量化设计；②功能不能完全满足使用要求的汽车通过完善功能实现汽车轻量化设计，比如改进动力学性能、NVH(Noise、Vibration、Harshness)性能等；③提高和改进性能的同时也使汽车减轻质量从而实现汽车轻量化.除去减轻质量外，汽车的轻量化设计还包含优化结构和价格合理化.

当l为阶数n或次数m时即为Wn和Wm，若Wm=1时，公式(2)演化成各向同性高斯滤波，若Wn=1, Wm=1，公式(2)退化成公式(1)；r为平滑半径.

本文采用的去相关滤波方法，基于最小二乘估计对l≥15阶的球谐系数拟合第l阶的相关误差，以第l阶为滑动窗中心向边界两端各反向延拓1/2窗宽的球谐系数，使得该方法在重力场模型系数的更高次数m上得到应用.建立法方程组求解三次多项式的未知参数a0、a1、a2、a3；b0、b1、b2、b3，构建相关误差模型为：

本项目根据中国石油天然气股份有限公司长庆油田分公司第一采气厂第三净化厂现场工艺状况及相关操作参数，自主设计并搭建了1套模拟现场脱除天然气中酸性组分的MDEA溶液吸收性能评价装置。

(3)

其次，泛北部湾经济合作区的成立还会促使区域文化竞争力的提升，打造区域文化品牌。文明是一个国家、一个民族生存的方式总结，是民族的灵魂和智慧的结晶。中国和东南亚一些国家都是历史悠久、具有丰富文明成果的国家。但是相对来讲，东亚和东南亚国家的文明还不足以与欧美一些国家相抗衡，或者形成对话。泛北部湾经济合作区的建设有利于各国、各地区开展丰富多彩的文化交流，有利于各国、各地区学习到先进、优秀的他国文化，提升自身的精神面貌。这种相互学习、团结合作、取长补短的合作方式，一定会使泛北部湾经济合作区的文化整体实力得到提高，在世界上也能迅速形成一种文化品牌，增强区域文化整体竞争力。

在托盘中心O点，托盘右侧B点，托盘上侧A点放置3个直径各不相同的圆柱，分别记为圆柱O、B、A。圆柱中轴位于非等腰直角三角形ABO的三个顶点上。以正方形托盘中心为原点如图建立坐标系，各圆的解析表达为：

4.合适的碳水化合物：碳水化合物的每日平均需要量为150g/d，提供的能量应占总能量的55%～65%，同时应避免摄入过多的精制糖，并且要注意粗细搭配。

Cm(l)-ΔCm(l)，

=Sm(l)-ΔSm(l).

1.2 单一尺度因子方法

在陆地水变化的时空尺度上，可认为研究区陆地水变化信号具有空间均一分布特点，且忽略区域外部信号变化影响(Chen et al., 2006).因此本文考虑基于区域整体及每个格网点的不同时段陆地水变化时间序列进行相同尺度的信号恢复.通过比较地表水文模型资料处理前后研究区的水文信号强度，估计一个信号振幅衰减比例因子即尺度因子，该尺度因子同样也代表了GRACE后处理所导致的信号衰减的幅度.这里采用单一尺度因子法对结果进行标定，使得经振幅恢复后的格点时间序列与真实时间序列尽可能接近.具体处理过程为：

(1) 提取全球水文模型GLDAS(NOAH)研究区格网数据，通过空间平均计算研究区陆地水变化时间序列ΔMT.

若将GRACE观测的地下水储量变化与地下水位变化进行等效对比，需将地下水位结果乘以含水层给水度(冯伟等, 2017).因此根据两者数值结果估计关中地区给水度约为0.06.

(3) 依据与原始时序的残差平方和最小，即

∑(ΔMT-scale·ΔST)2=min.

(4)

计算尺度因子scale，陆地水与地表水储量变化的振幅通过该尺度因子scale进行恢复.

1.3 水文数据

(1) GLDAS水文模型

活动结束要离开台商子弟学校了，才感叹：如此近的一个学校，十年来却没有来了解过，舍近求远往往错过很多机会，教育需要更多的开放！

本研究发现年龄＞60岁、冠心病、下肢动脉疾病是颈动脉狭窄的主要影响因素，P＜0.05；颈动脉狭窄与冠心病、下肢动脉疾病之间存在相关性，这也为脑卒中高危人群的筛查提供了依据。本研究显示重度颈动脉狭窄多发于冠心病、下肢静脉疾病的人群中，且年龄＞60岁的老年人成为高发人群。研究结果提示年龄＞60岁、冠心病和下肢静脉疾病患者成为颈动脉狭窄和脑卒中的高危人群，临床上对该群体进行颈动脉筛查后，能实现颈动脉狭窄和脑卒中的早发现和早治疗。

GLDAS全球水文模型由美国宇航局哥达航空中心(NASA)和美国国家环境预报中心(NCEP)共同建立(Rodell et al., 2004)，反映陆地表面土壤水和冰雪变化，空间分辨率为1°×1°，每月一值，包括NOAH、VIC、CLM、MOSAIC四种陆地表面模式.

丁辛醇装置是我公司产品结构调整的主体工程之一。该装置以丙烯和PSA制气装置提供的精合成气、氢气为原料，采用世界先进的低压羰基合成工艺生产丁辛醇。

(2) WGHM水文模型

同时将GLDAS格网数据球谐展开至与GRACE相同阶次，通过扣平均、组合滤波及尺度因子恢复法估算研究区地表水变化，为减小每种水文模型的不确定性，这里取GLDAS四种模型均值.对水储量变化格网经空间平均得到研究区整体变化结果.对比图1a两种时间序列可知：季节性变化在GRACE与GLDAS估算的水储量时间变化中表现明显，两者整体变化趋势较为一致，但由于GLDAS模型未考虑深层土壤水及地下水影响，因此模型振幅小于GRACE.图1b为基于最小二乘拟合计算研究区水储量的季节变化，振幅变化在±5.7 cm，GRACE估计的水储量最大亏损发生在每年的5—6月，最大盈余在每年的10月出现，与水文模型结果一致.

(3) 地下水测井资料

地下水测井资料来自陕西省地下水位年鉴，为关中地区监测井观测的地下水位埋深数据，观测时隔为5天，水位单位为m.为方便计算，本文直接引用年鉴的年平均水位数据.

2 结果与分析

2.1 GRACE与水文模型结果分析

利用2003年1月至2014年12月的月重力场模型对关中地区水储量变化进行反演.GRACE因其自身较高的轨道倾角，对重力场低阶项C20 敏感程度较低，可使用卫星激光测距SLR数据进行替换(Cheng and Tapley, 2005).扣除多年平均获得月重力场相对变化量，计算中缺失月份经插值得到.结合去相关滤波P3M15与 300 km Fan滤波估计该地区0.5°×0.5°的陆地水储量变化.最后采用单一尺度因子法对数据结果进行标定，乘以按公式(4)计算得到的尺度因子1.588，作为研究区的最终估计.

陆地水储量变化ΔTWS实际上包括土壤水、地下水、地表水、冰雪及生物量的质量变化总和(Strassberg et al., 2009).根据陕西省水资源公报结果估算近年来关中地区主要河流年平均径流量长期变化趋势约为0.006 km3·a-1(<0.01 km3·a-1).对长期陆地水变化影响来讲，地表径流变化可忽略不计.加之研究区的冰雪与生物量的变化相对较小，因此，本文可认为陆地水变化主要受土壤水与地下水变化影响.

ΔTWS=ΔSM+ΔGWS，

(5)

式(5)中，Δ表示为水储量的相对变化，SM为土壤水含量，GWS为地下水储量.当ΔTWS和ΔSM分别依靠GRACE和陆地表面水文模型获得时，地下水储量变化ΔGWS可通过ΔGWS=ΔTWS-ΔSM进行估计.

WGHM(WaterGAP Global Hydrology Model)全球地下水文模型(2.2c)由德国法拉克福大学自然地理研究所(IPG, University of Frankfurt)研制，该模型联合月时序参数以及气候和地理数据建立而成，时段为2003—2013年，空间分辨率为0.5°×0.5°，每月一值(Döll et al., 2014).

图1 GRACE与GLDAS均值估计的关中地区(a)水储量变化和(b)水储量季节性变化. Fig.1 Water storage variations (a) and water storage seasonal variations (b) in the Guanzhong area from GRACE and average GLDAS

为了定量分析GRACE与GLDAS模型之间的相关程度，对两者相关系数进行计算，如表1所示，相关系数多数高于0.7，GLDAS模型取平均后与GRACE相关系数明显增大，达到0.8，说明GLDAS平均结果与GRACE反演的陆地水储量变化更为吻合.

表1 GRACE与GLDAS及GLDAS均值相关系数统计 Table 1 Correlation between GRACE, GLDAS and average GLDAS

水文模型时间序列GLDASVIC NOAHCLM MOSAICGLDAS(average)与GRACE相关系数(无量纲)0.750.760.670.700.80

在利用尺度因子法获取区域时变重力场信号时，需对尺度因子改正前后GRACE及GLDAS结果的泄露误差进行估算.研究区水储量变化时间序列通过格网点空间平均得到，采用NOAH水文模型改正后结果均方根与滤波后的GRACE均方根RMSGRACE及NOAH均方根RMSmodel比值的乘积，如公式(6)，计算改正后的GRACE泄露误差

(6)

式(6)中，ΔMT、scale·ΔST在公式(4)中已定义.GRACE基于研究区的泄露误差为5.0 cm，经单一尺度因子改正后GRACE泄露误差减小为4.6 cm.同理GLDAS模型平均结果经改正前后的泄露误差分别为3.7 cm、3.4 cm.

2.2 地下水储量变化分析

2.2.1 与地下水模型对比分析

GRACE与GLDAS结果的相关性分析主要针对水储量随时间的变化，但评估地下水与模型结果之间的一致性程度应综合时空特征的比较.

直接从陆地水变化中扣除GLDAS均值获得研究区地下水储量变化等效水高，经线性回归得到地下水长期变化趋势的空间分布，如图2a所示，同时将WGHM依照GRACE处理过程估计地下水储量变化.

2003—2013年地下水储量变化率呈现从区域西南方向到东北方向依次递减的阶梯分布，表现出地下水变化趋势的显著空间差异性，这与该地区降雨由南向北、由西至东逐渐递减的地理分布一致(冯晶和钱会, 2015).区域范围内地下水表现出不同程度的亏损，变化幅度在-1.2～-0.2 cm·a-1，在东北区域达到最大，南部区域略有减少.但在WGHM估计结果中，西部及北部区域大幅亏损信号减弱，整个东南部区域与GRACE结果相比保持稳定，可以看出整个分布格局仍然从西南至东北依次减少，振幅变化在-0.7～-0.2 cm·a-1，如图2b所示.GRACE与WGHM估计结果在整体空间变化趋势上均表现为从西南至东北方向逐渐递减，但在局部长期变化速率上仍存在较大的差异性，这种差异可能由于WGHM在进行模型构建时采用的全球气候强度数据在局部空间上的不确定性导致模型本身存在一定误差，目前WGHM模型仍无法准确估计局部地下水变化的空间分布和亏损速率(Döll et al., 2014; 冯伟等, 2017).

图2 2003—2013年(a)GRACE得到的和(b)WGHM模型估计的关中地区地下水储量长期变化速率空间分布 Fig.2 Spatial patterns of long-time rate of GWS variations in the Guanzhong area from GRACE (a) and WGHM (b) from 2003 to 2013

直接利用公式(1)反演的水储量变化存在较明显的高频噪声，需截断高阶项球谐系数以降低高频误差影响，通常采用60阶重力场模型进行反演计算(Swenson and Wahr, 2006; Longuevergne et al., 2010; Feng et al., 2013; 吴云龙等, 2015; Belda et al., 2015).对接近60阶的高阶高次项位系数采用平滑半径为300 km的扇形滤波进行处理(Han et al., 2005; Zhang et al., 2009).奇偶阶位系数间各自序列的相关误差表现为“南-北”条带，利用P3M15反向边界延拓去相关滤波进行相关误差改正，同时调整固定滑动窗口宽度为5点(詹金刚等，2011).因单一滤波有其局限性，现将两种滤波方法进行组合使用，组合滤波计算模型如下：

式(1)中，θ为地心余纬，(φ,λ)为地面点的地心纬度和经度；(ΔCn m,ΔSn m)为n阶m次规格化位系数变化；为经规格化的n阶m次缔合勒让德函数；kn为n阶负荷勒夫数；ρw≈103 kg·m-3为水的密度；ρe≈5.5×103 kg·m-3为固体地球平均密度；R为地球平均半径.

图3 关中地区地下水监测井点位分布图 Fig.3 Locations of monitored well stations of groundwater in the Guanzhong area

2.2.2 地下水实测资料的对比分析

研究区地下水位监测井位置分布如图3所示.观测数据经时空平均得到地下水位年变化，地下水位升降反映在地球重力场随时间的演变，采用2003—2008年地下水位数据对GRACE反演结果进行验证.已有相关研究表明，近10年来该地区地下水位持续下降局面因受控采限制及引水工程实施得到遏制.图4a显示了GRACE观测2003—2008年的地下水储量相对多年平均有所增加，增加速率为0.25 cm·a-1.地下水位在相同时间内也处于逐年上升状态，上升速率约为4.39 cm·a-1，并在2004年表现出相比前一年的明显抬升，与GRACE观测结果比较一致，如图4b所示.

图4 (a) 2003—2013年GRACE得到的关中地区地下水储量变化及分段趋势；(b) 2003—2013年WGHM估计的 关中地区地下水储量变化及长期趋势；2003—2008年实测水位数据估计的水位变化及长期趋势 Fig.4 (a) GWS variations and segmental trend in the Guanzhong area from GRACE from 2003 to 2013; (b) GWS variations and long-time trend in the Guanzhong area from WGHM from 2003 to 2013; Water level change and long- time trend from measured water level data from 2003 to 2008

(2) 同时将原始格网数据球谐展开至与GRACE相同阶次，按照GRACE后处理过程得到时间变化序列ΔST.

2.2.3 降雨影响分析

大熊猫主食竹种有团竹(Fargesia obliqua)、冷箭竹(Bashania fangiana)、短锥玉山竹(Yushania brevipaniculata)、石棉玉山竹(Y. lineolata)等，因此大熊猫文化可视为四川竹文化的衍生[4]。四川省大熊猫栖息地作为世界自然遗产，跨越4个地级行政区，馆舍分布于7个自然保护区。其中雅安作为熊猫模式标本产地、熊猫文化发祥地以及“四川大熊猫栖息地”的核心区之一，现已拥有蜂桶寨国家自然保护区、碧峰峡大熊猫科研繁育基地等大熊猫基地，目前正在推进大熊猫国家公园体制试点工作。

降雨对关中地区地下水储量变化具有重要影响(贺旭波等, 2015)，本文进一步分析了两者之间的关系.采用中国气象数据网(http:∥data.cma.cn/data)提供的2003年1月—2014年12月的气象站月值降雨数据，得到关中地区整体降雨变化(分别以降水量实测值及距平值表示，单位mm).

图5a为2003—2014年地下水储量变化及同期降雨实测值，对原始时序扣除季节项并进行三个月滑动平均得到地下水年际变化，研究区地下水等效水高长期变化趋势为-0.44 cm·a-1，这与降雨长期时间分布较为一致，如图5b所示.

图5a显示了该地区每年7—10月存在大幅持续降雨，GRACE监测结果较降雨数据滞后，这主要因为降雨对地下水补给存在延迟所致.结合图5b可以发现，地下水储量年变化特征与降雨的年变化规律比较吻合.特别是在2003年，关中地区降雨量较多年平均高出286 mm，GRACE观测的地下水表现出快速回升；2007—2008年地下水先小幅上升而后发生下降突变，与同期降雨先小幅增多，而后大幅减少相一致；2012—2013年地下水处于亏损状态，同期降雨呈逐年减少趋势；2014年降雨偏多，地下水储量衰减趋于平稳.

3 结论

本文利用GRACE卫星重力数据可监测到关中地区地下水储量变化，基于单一尺度因子法研究关中地区地下水储量的时空变化特征，并与全球地下水文模型及实测资料进行了综合比较与分析，研究结论如下：

图5 (a) 2003—2014年GRACE得到的关中地区地下水储量变化及每月降雨量数据； (b) 2003—2014年关中地区年降雨异常 (a)中，黑色细虚线表示原始时间序列；黑色实线为扣除季节项并进行三个月滑动平均后的时间序列； 黑色粗虚线表示关中地区地下水储量变化长期趋势. Fig.5 (a) Time series of GWS variations from GRACE and the monthly rainfall data in the Guanzhong area from 2003 to 2014; (b) Annual rainfall anomalies in the Guanzhong area from 2003 to 2014 In (a), black thin dashed line indicates the original time series; black solid line is the time series after removing the seasonal term and sliding average of three months; black thick dashed line indicates the long-term trend of GWS variations in the Guanzhong area.

(1) GRACE反演的关中地区陆地水储量变化与GLDAS模型平均结果在时间序列上具有很强的相关性.

(2) GRACE观测的地下水长期亏损速率存在明显的空间差异性特征，与WGHM结果在时空分布上吻合较好；GRACE结果与地下水实测水位在局部时段内均出现正增长.

(3) 地下水储量变化与年降雨异常密切相关，在年际变化及长期趋势上均表现出较好的一致性.

(4) 本文给出了关中地区给水度约为0.06，该结果可为水文地质学等相关研究提供一定参考.

致谢 感谢德国法拉克福大学自然地理研究所Döll P教授提供的WGHM 2.2c全球地下水文模型.

References

Belda S, García-García D, Ferrándiz J M. 2015. On the decorrelation filtering of RL05 GRACE data for global applications. Geophysical Journal International, 200(1): 173-184.

Bettadpur S. 2007. UTCSR Level-2 gravity field product user handbook, GRACE 327-734, Center for Space Research. The University of Texas at Austin, Austin.

Chen J L, Wilson C R, Famiglietti J S, et al. 2006. Spatial sensitivity of the Gravity Recovery and Climate Experiment (GRACE) time-variable gravity observations. Journal of Geophysical Research: Solid Earth, 111(B6): 115-139, B08408, doi: 10.1029/2004JB003536.

Cheng M K, Tapley B D. 2005. Variations in the Earth's oblateness during the past 28 years. Journal of Geophysical Research: Solid Earth, 110(B9): B09402, doi: 10.1029/2004JB003028.

Döll P, Müller S H, Schuh C, et al. 2014. Global-scale assessment of groundwater depletion and related groundwater abstractions: Combining hydrological modeling with information from well observations and GRACE satellites. Water Resources Research, 50(7): 5698-5720.

Duan L, Wang W K, Sun Y B, et al. 2016. Iodine in groundwater of the Guanzhong Basin, China: Sources and hydrogeochemical controls on its distribution. Environmental Earth Sciences, 75: 970.

Famiglietti J S, Lo M, Ho S L, et al. 2011. Satellites measure recent rates of groundwater depletion in California′s Central Valley. Geophysical Research Letters, 38(3): L03403, doi: 10.1029/2010GL046442.

Feng J, Qian H. 2015. Analysis of precipitation characteristics in Guanzhong area. Journal of North China University of Water Resources and Electric Power (Natural Science Edition) (in Chinese), 36(3): 8-12, 48.

Feng W, Lemoine J M, Zhong M, et al. 2012. Terrestrial water storage changes in the Amazon basin measured by GRACE during 2002—2010. Chinese Journal of Geophysics (in Chinese), 55(3): 814-821, doi: 10.6038/j.issn.0001-5733.2012.03.011.

Feng W, Wang C Q, Mu D P, et al. 2017. Groundwater storage variations in the North China Plain from GRACE with spatial constraints. Chinese Journal of Geophysics (in Chinese), 60(5): 1630-1642, doi: 10.6038/cjg20170502.

Han S C, Shum C K, Jekeli C, et al. 2005. Non-isotropic filtering of GRACE temporal gravity for geophysical signal enhancement. Geophysical Journal International, 163(1): 18-25.

He X B, Tao F P, T H, et al. 2015. Analysis on groundwater level variation regulation of Guanzhong plain and its influencing factors. Ground Water (in Chinese), 37(3): 52-54.

Hu X G, Chen J L, Zhou Y H, et al. 2006. Seasonal water storage change of the Yangtze River basin detected by GRACE. Science in China Series D, 49(5): 483-491.

Iqbal N, Hossain F, Lee H, et al. 2016. Satellite gravimetric estimation of groundwater storage variations over Indus basin in Pakistan. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 9(8): 3524-3534.

Jekeli C. 1981. Alternative methods to smooth the Earth′s gravity field. Columbus, USA: Department of Geodetic Science and Surveying, Ohio State University.

Klees R, Zapreeva E A, Winsemius H C, et al. 2006. The bias in GRACE estimates of continental water storage variations. Hydrology and Earth System Sciences, 11(4): 1227-1241.

Landerer F W, Swenson S C. 2012. Accuracy of scaled GRACE terrestrial water storage estimates. Water Resources Research, 48(4): W04531, doi: 10.1029/2011WR011453.

Longuevergne L, Scanlon B R, Wilson C R. 2010. GRACE Hydrological estimates for small basins: Evaluating processing approaches on the High Plains Aquifer, USA. Water Resources Research, 46(11): W11517, doi: 10.1029/2009WR008564.

Luo Z C, Li Q, Zhong B. 2012. Water storage variations in Heihe river basin recovered from GRACE temporal gravity field. Acta Geodaetica et Cartographica Sinica (in Chinese), 41(5): 676-681.

Ma Z Y, Dou H P, Li T, et al. 2012. Geochemical study of groundwater mineralization in Guanzhong Basin, Shaanxi province, NW China. Desalination and Water Treatment, 42(1-3): 317-322.

Rodell M, Famiglietti J S. 2001. An analysis of terrestrial water storage variations in Illinois with implications for the Gravity Recovery and Climate Experiment (GRACE). Water Resources Research, 37(5): 1327-1339.

Rodell M, Houser P R, Jambor U, et al. 2004. The global land data assimilation system. Bulletin of the American Meteorological Society, 85(3): 381-394.

Strassberg G, Scanlon B R, Chambers D. 2009. Evaluation of groundwater storage monitoring with the GRACE satellite: Case study of the High Plains aquifer, central United States. Water Resources Research, 45(5): 195-211, W05410, doi: 10.1029/2008WR006892.

Swenson S, Wahr J. 2002. Methods for inferring regional surface-mass anomalies from Gravity Recovery and Climate Experiment (GRACE) measurements of time-variable gravity. Journal of Geophysical Research: Solid Earth, 107(B9): 2193, doi: 10.1029/2001JB000576.

Swenson S, Wahr J, Milly P C D. 2003. Estimated accuracies of regional water storage variations inferred from the Gravity Recovery and Climate Experiment (GRACE). Water Resources Research, 39(8): 1223, doi:10.1029/2002WR001808.

Swenson S, Wahr J. 2006. Post-processing removal of correlated errors in GRACE data. Geophysical Research Letters, 33(8): L08402, doi: 10.1029/2005GL025285.

Tao H, Tao F P, Liu W B. 2013. Characteristics and influencing factors of groundwater dynamics in Guanzhong urban agglomeration during the last 50 years. Hydrogeology & Engineering Geology (in Chinese), 40(6): 37-42, 61.

Tapley B D, Bettadpur S, Ries J C, et al. 2004. GRACE measurements of mass variability in the Earth system. Science, 305(5683): 503-505.

Tiwari V M, Wahr J, Swenson S. 2009. Dwindling groundwater resources in northern India, from satellite gravity observations. Geophysical Research Letters, 36(18): 184-201, L18401, doi: 10.1029/2009GL039401.

Wahr J, Molenaar M, Bryan F. 1998. Time variability of the Earth′s gravity field: Hydrological and oceanic effects and their possible detection using GRACE. Journal of Geophysical Research: Solid Earth, 103(B12): 30205-30230.

Wahr J, Swenson S, Zlotnicki V, et al. 2004. Time-variable gravity from GRACE: First results. Geophysical Research Letters, 31(11): 293-317, L11501, doi: 10.1029/2004GL019779.

Wahr J, Swenson S, Velicogna I. 2006. Accuracy of GRACE mass estimates. Geophysical Research Letters, 33(6): L06401, doi: 10.1029/2005GL025305.

Wu Y L, Li H, Zou Z B, et al. 2015. Investigation of water storage variation in the Heihe River using the Forward-Modeling method.Chinese Journal of Geophysics (in Chinese), 58(10): 3507-3516, doi: 10.6038/cjg20151007.

Zhan J G, Wang Y, Hao X G. 2011. Improved methods for removal of correlated errors in GRACE data. Acta Geodaetica et Cartographica Sinica (in Chinese), 40(4): 442-446, 453.

Zhang Z Z, Chao B F, Yang L, et al. 2009. An effective filtering for GRACE time-variable gravity: Fan filter. Geophysical Research Letters, 36(17): 1397-1413, L17311, doi: 10.1029/2009GL039459.

附中文参考文献

冯晶, 钱会. 2015. 关中地区降水特征分析. 华北水利水电大学学报(自然科学版), 36(3): 8-12, 48.

冯伟, Lemoine J M, 钟敏等. 2012. 利用重力卫星GRACE监测亚马逊流域2002—2010年的陆地水变化. 地球物理学报, 55(3): 814-821,doi: 10.6038/j.issn.0001-5733.2012.03.011.

冯伟, 王长青, 穆大鹏等. 2017. 基于GRACE的空间约束方法监测华北平原地下水储量变化. 地球物理学报, 60(5): 1630-1642, doi: 10.6038/cjg20170502.

贺旭波, 陶福平, 陶虹等. 2015. 关中平原地下水水位变化规律与影响因素分析. 地下水, 37(3): 52-54.

胡小工, 陈剑利, 周永宏等. 2006. 利用GRACE空间重力测量监测长江流域水储量的季节性变化. 中国科学 D辑: 地球科学, 36(3): 225-232.

罗志才, 李琼, 钟波. 2012. 利用GRACE时变重力场反演黑河流域水储量变化. 测绘学报, 41(5): 676-681.

陶虹, 陶福平, 刘文波. 2013. 关中城市群50年地下水动态变化及影响因素研究. 水文地质工程地质, 40(6): 37-42, 61.

吴云龙, 李辉, 邹正波等. 2015. 基于Forward-Modeling方法的黑河流域水储量变化特征研究. 地球物理学报, 58(10): 3507-3516, doi: 10.6038/cjg20151007.

詹金刚, 王勇, 郝晓光. 2011. GRACE时变重力位系数误差的改进去相关算法. 测绘学报, 40(4): 442-446, 453.