Variation and Attribution of Terrestrial Water Storage in the Yangtze River Basin
-
摘要:
气候变化和人类活动共同影响了长江流域陆地水储量,对中国经济社会、生态环境和粮食安全产生严重威胁。结合重力恢复与气候实验(gravity recovery and climate experiment, GRACE)卫星产品、地面水文观测、陆地水储量统计学模型和水文模型,分析了2002—2021年长江流域陆地水储量异常（terrestrial water storage anomaly,TWSA）的变化及归因。美国德克萨斯大学空间研究中心、喷气动力实验室和戈达德航天中心的GRACE mascon产品相对不确定性分别为5 mm、7.2 mm和3.9 mm；融合的GRACE陆地水储量产品与地面观测及流域水资源公报具有良好可比性,相关系数分别为0.59和0.86（p<0.01）,但不同子流域的对比情况存在一定差异；统计学模型和水文模型的TWSA分解结果较为一致,均揭示流域TWSA主要受降水影响；而人类活动因素,如水库蓄水和农业灌溉活动,则造成局部地区陆地水储量的明显变化。基于分离后的融合GRACE信号,发现长江流域历史TWSA已经因气候变率和人类活动而发生了明显变化。
Abstract:ObjectivesClimate variability and human activity have profoundly altered the terrestrial water storage in the Yangtze River basin, which has significantly threatened social economy, ecological environment, and food safety of China. However, the variations and attributions of the terrestrial water storage in the basin has not been comprehensively analyzed.
MethodsThis paper jointly used the GRACE (gravity recovery and climate experiment) solutions, field observations, statistical model, and large-scale hydrological simulation to perform the attribution of terrestrial water storage anomaly (TWSA) variations in the Yangtze River basin from 2002 to 2021.
ResultsThe relative uncertainties of GRACE mascon product from CSR (center for space research at University of Texas, Austin), JPL (jet propulsion laboratory), and GSFC (Goddard space flight center) are 5 mm, 7.2 mm, and 3.9 mm, respectively. The merged GRACE product from three institutions based on relative uncertainties present reasonable accuracy by comparing with the in-situ results and water resources bulletins, with the correlation coefficient of 0.59 and 0.86 (p<0.01), respectively. Similar partition results are reflected by the statistical and hydrological modeling results, in which the TWSA is mainly controlled by precipitation, while the human factors such as reservoir impoundment and agricultural irrigation can significantly affect the TWSA at the local scale.
ConclusionsThe historical TWSA of the Yangtze River basin has been altered extensively by climate variability and human activity based on the partitioned merged GRACE solution.
-
随着移动通信等技术的不断发展以及智能移动终端应用的广泛普及，位置服务已经成为构建智慧城市等方面的重要基础支撑[1-2]。人类80%~90%时间生活在室内，大型公共建筑的不断出现及其内部结构的日益复杂，促使人们对位置服务的需求从室外延伸至室内环境[3]。路径规划与导航是室内位置服务的重要功能[4]，它不仅能够辅助人们在室内空间的行为活动，而且对复杂室内及地下空间的安全管理与应急响应具有重要作用[5]。
多楼层叠加是室内外空间环境的重要区别，室内导航通常采用分楼层路径规划的方式，楼层间通过电梯等跨楼层通行设施实现连通[6]。因此，各楼层路径信息自动提取是实现室内实时导航的关键问题。目前，室内路径模型主要包括语义模型、栅格模型和网络模型等。其中，语义模型原理简单，但是不支持基于最短距离路径规划和规划结果可视化表达[7-8]；栅格模型易于创建，但是不适用于大范围室内环境[9-10]；网络模型将室内对象的连通关系抽象为具有几何和语义信息的节点及弧段，能够有效支持大范围室内环境路径的规划与导航，是当前应用最广泛的室内路径模型[11-12]之一。
室内路径主要表现为可通行的区域，如走廊、大厅以及房间内部无障碍的空间等。但是，由于网络模型需按照预先构建的几何路网计算最佳路径，从而造成路径规划结果容易存在曲折和迂回等情况[13]。另外，采用Delaunay三角网提取走廊中轴线[14-15]和Door-Door[16]生成房间内部路径的方法，难以处理空间内部存在障碍物等复杂情况[13]。然而，室内导航服务的应用场景多为大型公共建筑，空间结构复杂，现有方法难以实现室内路网的自动提取，需要通过矢量化的方式采集室内路径，从而降低了路径信息的获取效率。
面向室内路径的特征及其自动获取的难点，本文提出一种基于栅格空间描述室内可通行区域边界轮廓及其连通关系的通行区域模型，该模型不仅具备栅格模型易于创建的特点，而且兼具网络模型适用于大范围室内场景的优势。另外，结合室内地图数据特征[17]，设计了通行区域自动提取算法。通过具体试验表明，本文算法针对复杂室内环境能够有效实现通行区域的自动提取，相比广泛采用的网络模型，通行区域模型的路径规划结果更加符合人们在室内的路径行走特征。
1 通行区域模型的建模思路和过程
1.1 栅格模型与相邻栅格
栅格模型采用相同间隔对连续室内平面进行二维划分，并根据各栅格所代表实际空间范围内的室内对象特征判断是否可以通行，从而实现连续室内空间的离散化表示。如图 1所示的室内平面环境，图中A、B、C、D分别表示不同的室内空间，d1、d2、d3分别表示出入口，采用某栅格尺度对其进行离散化后结果如图 2所示。由于出入口是连接不同室内空间的重要对象，因此d1、d2、d3所处的栅格作为可以通行的缺口栅格。
本文采用8方向描述相邻栅格邻接关系，按顺时针方向采用整数1~8表示(见图 3)。基于栅格模型采用A*算法可以获取图 2中栅格P至栅格Q之间的最短路径。室内栅格中任意两个目标栅格之间直接通行经过的栅格集合称为途经栅格，具体表现为目标栅格中心点的连线相交的栅格集合(见图 4)。途经栅格是判断栅格之间是否可以直接通行和实现通行区域自动提取的重要基础，仅当两个通行栅格之间的途经栅格均可通行时，它们之间才能够直接通行。
1.2 通行区域模型构建方法
基于栅格模型进行路径规划需要不断遍历当前栅格的相邻栅格，从而判断下一步的行进方向。因此，当室内范围大或栅格尺度小，会导致栅格数目增多，进而增加了路径规划算法的搜索时间，难以满足室内实时路径规划与导航的应用需求。为此，本文在栅格模型基础上构建了通行区域模型，它是对关联相同室内对象的通行栅格进行聚合，从而形成内部任意栅格之间可以无障碍通行的区域，同时建立相邻通行区域之间的连通关系。图 5是图 1中室内平面环境采用通行区域模型表达的路径信息，其中，空白部分表示不可通行的障碍区域，红色虚线分别代表了通行区域R1与R3基于栅格空间的区域范围。
为了实现通行区域的自动提取，本文对传统栅格模型进行了扩展，主要包括关联对象、隶属区域和栅格类型3个方面。关联对象用于表示栅格语义信息，通过关联的室内对象描述，如图 5中栅格a关联对象为室内空间A。隶属区域用于表示栅格所属的通行区域，如果属于多个通行区域，那么所隶属通行区域之间可以通过该栅格实现连通，如图 5中栅格d的隶属区域为R5与R6。本文进一步将通行栅格划分为内部栅格、边界栅格和枢纽栅格3种类型。
1) 内部栅格是唯一隶属于某通行区域R的通行栅格，并且其奇数邻接方向的相邻栅格均为仅隶属于R的通行栅格。如图 5中栅格P与栅格Q分别为通行区域R1和R3的内部栅格。
2) 边界栅格是唯一隶属于某通行区域R，但不属于R内部栅格的通行栅格。边界栅格主要用于描述通行区域的范围，如图 5中栅格k为R5的边界栅格。
3) 枢纽栅格是同时隶属于多个通行区域的通行栅格，它主要用于建立相邻通行区域的连通关系和描述通行区域的范围。如图 5中栅格c、m、l为R4与R5的枢纽栅格。
如图 5所示，R1内栅格P可直接行进至枢纽栅格a，经过D1可以直接行进至枢纽栅格b，再经过R4可直接行进至与R5连通的枢纽栅格c、m、l或与D2连通的枢纽栅格。基于上述原理，采用A*算法得到栅格P至栅格Q的最短路径为P-a-b-c-d-n-p-q-Q。
2 室内通行区域的自动提取
2.1 室内栅格模型初始化
室内栅格模型初始化按照相应栅格尺度生成离散化栅格，同时初始化栅格的类型、关联对象、隶属区域等信息。假设栅格尺度为S，栅格模型初始化的主要流程如下。
1) 划分离散栅格。选取室内地图最小横坐标与纵坐标为原点，沿坐标轴方向采用间隔S将室内平面划分为一系列的离散栅格。栅格尺度越小，反映的室内环境特征越精细，但会导致栅格数目增加，从而影响通行区域的提取效率；栅格尺度较大，反映的室内环境特征越粗略，可能会导致走廊等狭窄室内空间连通性的缺失。
2) 判断栅格属性。如果某室内空间Rs为可通行的区域，那么位于Rs内部的栅格初始化为边界栅格，关联对象为Rs，隶属区域为空；与Rs边界相交的栅格初始化为障碍栅格，关联对象和隶属区域均为空。如果Rs为不可通行的区域，那么与Rs的边界及其内部相交的栅格均初始化为障碍栅格，关联对象和隶属区域均为空。
3) 生成缺口栅格。由于门等缺口依附于墙体等室内空间的边界，从而造成缺口中心点所处栅格在第2)步判断为障碍栅格。因此，如果缺口d中心点位于栅格G内部，那么设置栅格G为缺口栅格，关联对象为d。如图 6所示为图 1中室内平面环境的栅格模型初始化结果，其中缺口栅格G1、G2、G3的关联对象分别为出入口d1、d2、d3。如果初始化之后经过某种方式检验发现存在室内空间连通性缺失的情况，那么需要减小栅格尺度S，并重新进行初始化。
2.2 通行区域初次提取
建立列表RList存储提取到的通行区域，GList[R]存储区域R所有栅格，BGList[R]存储R边界及枢纽栅格。建立待扩展栅格列表SGList和种子栅格SG，通行区域初次提取主要流程如图 7所示，主要包括3个部分。
1) 获取有效种子栅格。通过遍历室内每行栅格或待扩展栅格列表，获取可用于提取通行区域的种子栅格SG，如果某栅格为障碍栅格或已存在隶属区域，则不能作为种子栅格。
2) 提取各个通行区域。以SG为起点新建区域R，通过遍历R边界栅格，不断将满足条件的相邻栅格加入R。其中，BTrav(Gk, R)用于判断Gk是否可以加入R，对于BGList[R]中任意栅格Gi，仅当Gk与Gi之间途经栅格均可通行且不为其他区域的内部及枢纽栅格时，BTrav(G, R)=true。向R添加Gk时，需将R加入Gk隶属区域，同时更新Gk及其相邻栅格的类型。
3) 建立区域连通关系。根据各个缺口区域包含的缺口栅格及其相邻栅格的隶属区域，建立缺口区域与相邻空间区域之间的连通关系。图 6的栅格模型初始化结果经过通行区域初次提取后得到的结果如图 8所示，包含R1~R8共8个空间区域和D1~D3共3个缺口区域。
2.3 通行区域邻域融合
由于室内空间边界不规则等原因，导致初次提取结果通常包含部分栅格数目较少的通行区域，如图 8中R7和R8所示。为降低通行区域复杂度，需要对满足融合阈值和融合条件的区域进行融合。其中，融合阈值是指通行区域的非枢纽栅格小于一定数目；融合条件是指区域融合后不会降低与其连通的通行区域之间的连通性。另外，邻域融合需要按照区域栅格数目由少至多进行。由于出入口是建立相邻室内空间连通性的重要对象，因此缺口区域不可融合。假设某空间区域R满足融合阈值且其连通区域数目为N，若N=0，可直接将R从RList删除；若N=1，仅当R的连通区域不为缺口区域时可以融合，需从其连通区域的枢纽栅格隶属区域中删除R，并将R从RList删除；若N>1，邻域融合的主要流程如图 9所示，主要包括3个部分。
1) 判断R是否满足融合条件。对于R的任意连通区域R1、R2，仅当二者直接连通，或通过某不为R的区域间接连通，或R存在栅格G满足BMerg(G, R, R1)=true且BMerg(G, R, R2)=true，则R满足融合条件。其中，BMerg(G, RA, RB)用于判断区域RA中栅格G至区域RB的可融性，对于BGList[RB]中任意栅格Gi，仅当G与Gi之间途经栅格均可通行且不为除RA与RB外其他区域的内部及枢纽栅格时，BMerg(G, RA, RB) =true。
2) 从R中选取有效扩展栅格。为了将R中栅格尽量融合至相邻区域，需要在R中选取相邻区域共有的栅格进行扩展。首先，从R所有栅格的隶属区域删除R，并更新栅格类型，比如R中内部栅格更新为隶属区域为空的边界栅格；然后，从R中选取隶属区域不为空的边界栅格作为有效扩展栅格，即可以通过该栅格将其相邻栅格加入隶属的通行区域。最后，当R中不存在有效扩展时，邻域融合结束，并将R从RList删除。
3) 将R栅格加入连通区域。通过遍历扩展栅格的相邻栅格，不断地将R中满足隶属区域为空等条件的边界栅格加入扩展栅格隶属的通行区域，同时更新加入边界栅格及其相邻栅格的类型。假设图 8中R7与R8满足融合阈值，R7不满足融合条件不能进行融合，R8存在栅格G(如n、u、x、z)满足BMerg(G, R8, R6)=true且BMerg(G, R8, R7)=true，因此可以进行融合，按照图 9流程进行邻域融合以后得到的结果如图 5所示。
3 试验与分析
选取北京西单大悦城一楼室内地图数据进行通行区域自动提取和路径规划试验，试验数据分为楼层边界、室内空间、出入口和室内路网4层(见图 10)。其中，室内空间包括店铺、天井等类型的面状要素；走廊为楼层边界范围内除去室内空间的部分。
3.1 通行区域提取试验与分析
由于人行走的基本步距约为0.75 m，因此分别采用1.5 m和0.75 m两种栅格尺度进行通行区域的自动提取试验。根据试验区的数据特征，以上两种栅格尺度初始化结果分别如图 11(a)和图 11(b)所示。通过图 11可以看出，采用以上两种栅格尺度进行初始化后均未造成室内空间连通性的缺失。
图 11(a)经过通行区域初次提取得到90个空间区域和35个缺口区域(见图 12(a))。设置图 12(a)融合阈值为9个栅格，经过邻域融合后剩余69个空间区域和35个缺口区域(见图 12(b))。
![]()
图 12 试验楼层1.5 m栅格尺度通行区域提取结果Figure 12. Result of Experiment Floor Traversable Region Extraction Based on 1.5 m Grid Scale图 11(b)经过初次提取得到128个空间区域和29个缺口区域(见图 13(a))。根据1.5 m栅格尺度融合阈值对应的实际室内面积，设置图 13(a)融合阈值为36个栅格，经过邻域融合后剩余75个空间区域和35个缺口区域(见图 13(b))。可以看出，栅格尺度越小，初次提取的通行区域数目相对较多，特别是几何特征较为复杂的室内空间需要进行邻域融合的通行区域数目较多；不同栅格尺度最终提取的通行区域数目差别较小，邻域融合后的通行区域能够较好地反映室内主要路径信息。
![]()
图 13 试验楼层0.75 m栅格尺度通行区域提取结果Figure 13. Result of Experiment Floor Traversable Region Extraction Based on 0.75 m Grid Scale3.2 室内路径规划试验与分析
选取“APPLE”中心为起点，终点为“NHIZ”，采用相同试验环境和A*算法，分别基于室内路网和通行区域数据进行路径规划，结果如图 14所示。其中，基于图 10中室内路网规划距离为134.51 m，路径规划平均耗时20.94 ms；基于图 12(b)中通行区域规划距离为125.99 m，路径规划平均耗时53.27 ms；基于图 13(b)中通行区域的规划距离为122.51 m，路径规划平均耗时78.04 ms。可以看出，通行区域模型相比网络模型能够有效改善路径规划结果存在较多曲折等情况，更加符合复杂室内环境路径行走特征；路径规划耗时相对较长，但能够满足实时路径规划应用需求。基于北京智慧图公司室内位置服务平台和图 13(b)提取的通行区域，实地进行了路径导航试验(见图 15)，验证了通行区域模型对满足大范围室内场景中实时路径规划与导航的适用性。
![]()
图 15 基于通行区域模型的室内路径实时规划及表达Figure 15. Indoor Route Real-Time Planning and Presentation Based on Traversable Region Model4 结语
室内路径信息自动获取是实现大规模精准室内导航服务的重要基础，也是当前亟需解决的关键技术之一。本文面向室内路径信息的特点，提出一种适用于室内路径规划与导航的通行区域模型，并设计了通行区域的自动提取算法，有效实现了室内路径信息的自动获取。
本文建立的通行区域模型相比当前广泛采用的网络模型，能够有效改善路径规划结果容易存在曲折等情况，不仅能够满足路径实时规划与导航的应用需求，而且更加符合人们在复杂室内环境的路径行走特征；另外，通行区域自动提取算法能够有效适应大型公共建筑的复杂室内环境。
本文分别选取了两种栅格尺度进行了通行区域自动提取和路径规划试验，并对试验结果进行了对比分析。面向不同室内空间特征的大型公共建筑，还需要对不同栅格尺度的适用性进行更为深入的分析，这将是本文下一步的研究内容。
http://ch.whu.edu.cn/cn/article/doi/10.13203/j.whugis20230017 -
![]()
图 1 长江流域地理位置及观测站点空间分布
Figure 1. Location of the Yangtze River Basin and Spatial Distribution of Observation Stations
![]()
图 2 统计学方法和WGHM水文模型在长江流域的陆地水储量分解结果对比
Figure 2. Comparison of TWSA Partitions Based on Statistical Model and WGHM in the Yangtze River Basin
![]()
图 3 2003-2016年长江流域不同陆地水储量成分变化趋势
Figure 3. Spatial Distribution of Trend Results in Different TWSA Components in the Yangtze River Basin from 2003 to 2016
![]()
图 4 2003-2016年长江流域不同气候因子及非气候因子变化趋势
Figure 4. Spatial Distribution of Trend Results in Different Climatic Factors in the Yangtze River Basin from 2003 to 2016
![]()
图 5 长江流域灌溉面积及水库总库容年际变化
Figure 5. Annual Variation of Total Irrigation Area and Reservoir Volume in the Yangtze River Basin
表 1 长江流域融合GRACE产品与观测结果及水资源公报之间的相关性
Table 1 Correlations Between Merged GRACE Product with Observed Results and Water Resource Bulletin Data in the Yangtze River Basin
流域 水文站 流域面积/km2 相关系数 观测结果 水资源公报 金沙江 石鼓 214 184 0.67 0.55 雅砻江 桐子林 127 624 0.63 0.45 岷沱江 高场 351 535 0.55 0.6 嘉陵江 北碚 156 124 0.36 0.65 乌江 武隆 83 035 0.28 0.66 汉江 皇庄 142 056 0.62 0.76 洞庭湖 城陵矶 261 111 0.37 0.76 鄱阳湖 湖口 162 225 0.59 0.85 长江 大通 1 705 383 0.59 0.86 注： 流域面积指各流域出口水文站的控制面积。
下载: 导出CSV
表 2 长江流域水资源公报统计水资源总量年际变化/（108 m3）
Table 2 Annual Variations of Total Water Resources from Water Resource Bulletins/（108 m3）
年份 金沙江(石鼓以上) 金沙江(石鼓以下) 岷沱江 嘉陵江 乌江 汉江 洞庭湖 鄱阳湖 长江 2004年 380 1 139 1 123 539 533 504 2 017 1 006 8 734 2005年 472 1 217 1 234 775 434 763 1 957 1 465 9 889 2006年 338 843 785 400 399 393 2 057 1 582 8 061 2007年 403 1 031 897 680 578 592 1 794 1 085 8 811 2008年 420 1 243 990 605 536 517 1 972 1 335 9 457 2009年 502 965 932 666 442 557 1 751 1 148 8 733 2010年 442 1 043 1 047 745 491 771 2 312 2 233 11 266 2011年 362 724 976 754 341 673 1 410 1 021 7 838 2012年 458 1 116 1 225 753 511 470 2 373 2 149 10 809 2013年 384 905 1 105 745 401 395 1 962 1 389 8 798 2014年 444 1 076 1 054 650 578 448 2 229 1 606 10 150 2015年 302 1 100 849 537 482 474 2 364 1 986 10 330 2016年 358 1 190 925 463 598 502 2 642 2 173 11 947 2017年 451 1 130 967 692 526 716 2 327 1 608 10 616 2018年 497 1 224 1 304 738 536 458 1 679 1 140 9 374 2019年 441 987 1 193 847 562 496 2 508 2 006 10 543 2020年 492 1 202 1 363 946 714 704 2 670 1 676 12 863 2021年 444 989 1 083 1 112 593 1 033 2 220 1 414 11 186
下载: 导出CSV
-
[1] Maggioni V. Extreme Hydroclimatic Events and Multivariate Hazards in a Changing Environment： A Remote Sensing Approach[M]. Amsterdam, Netherlands： Elsevier, 2019.
[2] Cosgrove W J, Loucks D P. Water Management： Current and Future Challenges and Research Directions[J]. Water Resources Research, 2015, 51（6）： 4823-4839.
[3] Reager J T, Gardner A S, Famiglietti J S, et al. A Decade of Sea Level Rise Slowed by Climate-Driven Hydrology[J]. Science, 2016, 351（6274）： 699-703.
[4] Tapley B D, Watkins M M, Flechtner F, et al. Contributions of GRACE to Understanding Climate Change[J]. Nature Climate Change, 2019, 5（5）： 358-369.
[5] Wada Y, de Graaf I E M, van Beek L P H. High-Resolution Modeling of Human and Climate Impacts on Global Water Resources[J]. Journal of Advances in Modeling Earth Systems, 2016, 8（2）： 735-763.
[6] Long D, Yang W T, Scanlon B R, et al. South-to-North Water Diversion Stabilizing Beijing’s Groundwater Levels[J]. Nature Communications, 2020, 11（1）： 3665.
[7] Zhang B Q, Tian L, Yang Y T, et al. Revegetation Does Not Decrease Water Yield in the Loess Plateau of China[J]. Geophysical Research Letters, 2022, 49（9）： e2022GL09802.
[8] 王文, 王鹏, 崔巍. 长江流域陆地水储量与多源水文数据对比分析[J]. 水科学进展, 2015, 26（6）： 759-768. Wang Wen, Wang Peng, Cui Wei. A Comparison of Terrestrial Water Storage Data and Multiple Hydrological Data in the Yangtze River Basin[J]. Advances in Water Science, 2015, 26（6）： 759-768.
[9] 聂宁, 张万昌, 陈豪, 等. 长江流域水储量变化的时间变化特征及归因分析[J]. 水科学进展, 2021, 32（3）： 396-407. Nie Ning, Zhang Wanchang, Chen Hao, et al. Temporal Variation Characteristics and Attribution Analysis of Terrestrial Water Storage Change in the Yangtze River Basin[J]. Advances in Water Science, 2021, 32（3）： 396-407.
[10] 马倩, 谢正辉, 陈锋, 等. 长江流域1982—2005年陆地水储量变化及时空分布特征[J]. 气候与环境研究, 2011, 16（4）： 429-440. Ma Qian, Xie Zhenghui, Chen Feng, et al. Temporal and Spatial Distribution of Terrestrial Water Storage in the Yangtze River Basin During 1982-2005[J]. Climatic and Environmental Research, 2011, 16（4）： 429-440.
[11] 冉艳红, 钟敏, 陈威, 等. 利用GRACE-FO重力卫星探测2019年长江中下游极端干旱[J]. 科学通报, 2021, 66（1）： 107-117. Ran Yanhong, Zhong Min, Chen Wei, et al. Monitoring the Extreme Drought in the Middle and Lower Reaches of the Yangtze River in 2019 from GRACE-FO Satellites[J]. Chinese Science Bulletin, 2021, 66（1）： 107-117.
[12] 李晓英, 叶根苗, 蔡晨凯, 等. 基于GRACE和MODIS数据的长江流域陆地水储量变化分析及预测[J]. 长江科学院院报, 2018, 35（5）： 130-135. Li Xiaoying, Ye Genmiao, Cai Chenkai, et al. Analysis and Prediction of the Anomaly of Terrestrial Water Storage in the Yangtze River Basin Based on MODIS and GRACE[J]. Journal of Yangtze River Scientific Research Institute, 2018, 35（5）： 130-135.
[13] 王微, 李伟伟, 王奉伟. 利用GRACE模型反演长江流域陆地水储量的变化[J]. 人民长江, 2020, 51（10）： 110-115. Wang Wei, Li Weiwei, Wang Fengwei. Terrestrial Water Storage Variations in Yangtze River Basin Based on Tongji-Grace 2018 Temporal Gravity Field Model[J]. Yangtze River, 2020, 51（10）： 110-115.
[14] Ministry of Water Resources of the People’s Republic of China. Code of Practice for Water Resources Protection Planning： SL 613—2013[S]. Beijing： China Water Resources and Hydropower Press, 2013.） 中华人民共和国水利部. 水资源保护规划编制规程： SL 613—2013[S]. 北京： 中国水利水电出版社, 2013.
[15] 王旭, 林征, 张志, 等. 基于GWR模型的北极滨海平原融冻湖表面温度空间分布模拟[J]. 武汉大学学报（信息科学版）, 2016, 41（7）： 918-924. Wang Xu, Lin Zheng, Zhang Zhi, et al. Modelling the Spatial Distribution of Lake Surface Water Temperature of the Thaw Lakes in Arctic Coastal Plain Using Geographically Weighted Regression Model[J]. Geomatics and Information Science of Wuhan University, 2016, 41（7）： 918-924.
[16] 袁乐先, 李斐, 张胜凯, 等. 利用ICESat/GLAS数据研究北极海冰干舷高度[J]. 武汉大学学报（信息科学版）, 2016, 41（9）： 1176-1182. Yuan Lexian, Li Fei, Zhang Shengkai, et al. A Study of Arctic Sea Ice Freeboard Heights from ICESat/GLAS[J]. Geomatics and Information Science of Wuhan University, 2016, 41（9）： 1176-1182.
[17] 尹国应, 张洪艳, 张良培. 2001—2019年长江中下游农业干旱遥感监测及植被敏感性分析[J]. 武汉大学学报（信息科学版）, 2022, 47（8）： 1245-1256. Yin Guoying, Zhang Hongyan, Zhang Liangpei. Remote Sensing Monitoring of Agricultural Drought and Vegetation Sensitivity Analysis in the Middle and Lower Reaches of the Yangtze River from 2001 to 2019[J]. Geomatics and Information Science of Wuhan University, 2022, 47（8）： 1245-1256.
[18] Scanlon B R, Zhang Z Z, Save H, et al. Global Evaluation of New GRACE Mascon Products for Hydrologic Applications[J]. Water Resources Research, 2016, 52（12）： 9412-9429.
[19] Mo S X, Zhong Y L, Forootan E, et al. Bayesian Convolutional Neural Networks for Predicting the Terrestrial Water Storage Anomalies During GRACE and GRACE-FO Gap[J]. Journal of Hydrology, 2022, 604： 127244.
[20] Weedon G P, Balsamo G, Bellouin N, et al. The WFDEI Meteorological Forcing Data Set： WATCH Forcing Data Methodology Applied to ERA-Interim Reanalysis Data[J]. Water Resources Research, 2014, 50（9）： 7505-7514.
[21] Xiong J H, Guo S L, Chen D L, et al. Past and Future Terrestrial Water Storage Changes in the Lower Mekong River Basin： The Influences of Climatic and Non-climatic Factors[J]. Journal of Hydrology, 2022, 612： 128275.
[22] Long D, Longuevergne L, Scanlon B R. Uncertainty in Evapotranspiration from Land Surface Modeling, Remote Sensing, and GRACE Satellites[J]. Water Resources Research, 2014, 50（2）： 1131-1151.
[23] Humphrey V, Gudmundsson L. GRACE-REC： A Reconstruction of Climate-Driven Water Storage Changes over the Last Century[J]. Earth System Science Data, 2019, 11（3）： 1153-1170.
[24] Müller Schmied H, Cáceres D, Eisner S, et al. The Global Water Resources and Use Model WaterGap V2.2d： Model Description and Evaluation[J]. Geo⁃scientific Model Development, 2021, 14（2）： 1037-1079.
[25] Zhang C, Dong J W, Ge Q S. IrriMap_CN： Annual Irrigation Maps Across China in 2000-2019 Based on Satellite Observations, Environmental Variables, and Machine Learning[J]. Remote Sensing of Environment, 2022, 280： 113184.
[26] 田小娟, 邹芳, 金双根. 气候变化和人类活动对长江流域水储量变化的影响研究[J]. 大地测量与地球动力学, 2019, 39（4）： 371-376. Tian Xiaojuan, Zou Fang, Jin Shuanggen. Impact of Climate Change and Human Activities on Water Storage Changes in the Yangtze River Basin[J]. Journal of Geodesy and Geodynamics, 2019, 39（4）： 371-376.
[27] 孟莹, 刘俊国, 王子丰, 等. 气候变化和人类活动对中国陆地水储量变化的影响研究[J]. 华北水利水电大学学报（自然科学版）, 2021, 42（4）： 47-57. Meng Ying, Liu Junguo, Wang Zifeng, et al. Climate Change and Human Impacts on Terrestrial Water Storage Anomalies in China[J]. Journal of North China University of Water Resources and Electric Power （Natural Science Edition）, 2021, 42（4）： 47-57.
-
期刊类型引用(1)
1. 陈晴空，陈治君，王殿常，兰峰，高千红，李姗泽. 三峡库区悬沙颗粒磷的粒径变化和滞留特征. 水利水电技术(中英文). 2025(01): 181-192 . 
百度学术
其他类型引用(1)
计量
- 文章访问数: 722
- HTML全文浏览量: 131
- PDF下载量: 137
- 被引次数: 2
