考虑高程的降水空间插值方法研究——以巴河流域为例.pdf

中国农村水利水电China Rural Water and Hydropower考虑高程的降水空间插值方法研究以巴河流域为例鞠琴1，2，刘小妮1，2，连子旭1，张译尹1，吴金雨1，2，尼玛扎西3，刘晓丽1，4（1.河海大学 水灾害防御全国重点实验室，江苏 南京 210098；2.河海大学水安全与水科学协同创新中心，江苏 南京 210098；3.西藏自治区水文水资源勘测局阿里水文水资源分局，西藏 拉萨 850032；4.安徽农业大学工学院，安徽 合肥 230036）摘要：在气候变暖的背景下，构建质量可靠的降水数据，对水资源管理与水土流失预防与治理至关重要。基于站点观测数据，利用空间插值方法是获得可靠的降水数据的重要途径。为提高空间插值方法对水文站点稀疏、地形复杂的山区流域降水数据精度，以巴河流域为例，采用径流系数作为评价指标，评估GLADS、GSMaP以及再分析数据集CMFD等3种降水数据的适用性，并考虑高程变化对降水空间垂直分布的影响，提出了一种基于高程递推的反距离加权法，分析研究区不同高程降水的空间变化规律。结果表明：GLADS、GSMaP以及再分析数据集CMFD降水数据对巴河流域的年均降水量的模拟存在不同程度的低估现象，模拟精度从大到小排序为CMFD GSMaP GLADS，均不能准确的刻画巴河流域降水量随高程的变化。巴河流域水文观测站点较少，代表性有限，遥感数据反演精度偏低，基于高程递推的反距离加权法模拟的巴河流域降水空间分布再现效果较好，明显优于传统的反距离加权法。该方法得到的降水量能够更好的反映高程变化对巴河流域降水空间分布的影响，巴河流域年均降水量为 1 522 mm，径流系数为 0.92，降水量主要来源于 4 0006 000 m高海拔地区，相对贡献率最大，占比高达89.5%，4 000 m以下海拔对应的降水量占比为10.1%。研究可为水文资料短缺的高寒山区流域的水文模拟与预报提供参考。关键词：降水；反距离加权；径流系数；高程递推；巴河流域中图分类号：P426.6 文献标识码：A DOI：10.12396/znsd.222065鞠 琴，刘小妮，连子旭，等.考虑高程的降水空间插值方法研究：以巴河流域为例 J.中国农村水利水电，2023（8）：18-24.DOI：10.12396/znsd.222065.JU Q，LIU X N，LIAN Z X，et al.Research on the spatial interpolation method of precipitation considering elevation：a case study of Bahe River Basin J.China Rural Water and Hydropower，2023（8）：18-24.DOI：10.12396/znsd.222065.Research on the Spatial Interpolation Method of Precipitation Considering Elevation：A Case Study of Bahe River BasinJU Qin1，2，LIU Xiao-ni1，2，LIAN Zi-xu1，ZHANG Yi-yin1，WU Jin-yu1，2，NIMA Zha-xi3，LIU Xiao-li1，4（1.The National Key Laboratory of Water Disaster Prevention，Hohai University，Nanjing 210098，Jiangsu Province，China；2.Cooperative Innovation Center for Water Safety and Hydro Science，Hohai University，Nanjing 210098，Jiangsu Province，China；3.Ali Hydrology and water Resources Branch of the Tibet Autonomous Region Bureau of Hydrology and water Resources Survey，Lhasa 850032，Tibet，China；4.School of Engineering，Anhui Agricultural University，Hefei 230036，Anhui Province，China）Abstract：In the context of climate warming，the construction of reliable precipitation data is crucial for water resources management and soil 文章编号：1007-2284（2023）08-0018-07收稿日期：2022-10-22基金项目：国家重点研发计划资助项目（2021YFC3201104）；国家自然科学基金项目（52179013）；“一带一路”水与可持续发展科技基金项目（2021491011）。作者简介：鞠琴（1980-），女，副教授，博士，主要从事变化环境下水文循环响应研究。E-mail：。通讯作者：刘小妮（1998-），女，硕士研究生，主要从事水文水资源及气候变化方面的研究。E-mail：。水文水资源18考虑高程的降水空间插值方法研究以巴河流域为例 鞠琴 刘小妮 连子旭 等erosion prevention and control.Spatial interpolation is an important way to obtain reliable precipitation data based on site observation data.To improve the accuracy of precipitation data in alpine mountainous watersheds with sparse hydro-meteorological stations，this paper uses runoff coefficient as an evaluation index to assesse the applicability of GLADS，GSMaP，and the reanalysis dataset CMFD in the Bahe River Basin.Then，considering the influence of elevation on the precipitation variations，an improved inverse distance weighting method based on elevation recurrence is proposed to analyze the spatial variations of precipitation at different elevations.The results show that the average annual precipitation simulated by GLADS，GSMaP and reanalysis data set CMFD precipitation data is underestimated to varying degrees，and their simulation accuracy is ranked as CMFD GSMaP GLADS，which is not capable of characterizing the variations of precipitation by elevation in the Bahe River Basin.There are few hydrological observation stations in the Bahe River Basin with limited representativeness，resulting in low accuracy of remote sensing data inversion.The precipitation data interpolated by the improved inverse distance weighting method effectively reproduced the spatial distribution of precipitation variation by elevation，which is significantly better than the traditional inverse distance weighting method.The simulation of precipitation showed that this method can well reflect the influence of elevation change on the spatial distribution of precipitation.The average precipitation in the Bahe River Basin calculated by the improved inverse distance weighting method was 1 522 mm with a runoff coefficient of 0.92，and the largest relative contribution rate of precipitation in areas with the elevation of 4 0006 000 m was obviously higher than other areas，accounting for 89.5%，and that of precipitation below 4 000 m was 10.1%.This paper provides a new technical support for hydrological simulation and forecasting in alpine mountainous watersheds with scarce hydrological information.Key words：precipitation；inverse distance weighting method；runoff coefficient；elevation recursion；Bahe river basin0引 言准确、精细化的降水数据对研究流域水文特征、产汇流机理和建立流域水文模型至关重要1，2。高原山区降水时空分布差异性与水量平衡问题的研究一直是水文学研究的难点3。高原山区地形起伏较大，土壤与植被等下垫面条件复杂，使降水空间分布存在显著的差异性4，导致高原山区出现降水径流水量不平衡问题。对于水文站点稀少，资料短缺或者无观测资料的高原山区，常采用遥感降水数据，主要包括TRMM、GLADS、GSMaP、CMFD和GPCP等5，为水文、气象、灾害防治等领域的研究提供了重要的途径。目前，传统的降水空间插值方法主要有反距离加权法、克里金法、样条函数法等6，由于流域局部地形、高山与植被以及其他影响因素引起的空间差异7，导致模拟的降水空间分布结果误差很大8。因此，针对影响降水空间分布的地理因素，国内外学者开展了研究，陈裕迪等9研究表明，五种降水数据在红河流域的模拟精度在不同程度上均受到高程的影响；朱浩楠等10采用考虑地形Cokriging与PRISM算法对重庆市降水进行插值分析，发现PRISM模拟效果主要依赖于样本数量；李豪等11基于考虑了空间位置、地形等因子的混合地理加权回归克里格模型，发现其对四川省的降水数据模拟精度优于克里格插值方法。由于降水的空间复杂性，学者们认识到地形地貌因素对降水插值精度的影响12，但高原山区缺乏降水实测数据或高程影响权重较小等原因，使高程梯度变化对降水插值效果的影响有待进一步研究13。例如，刘田等14发现气象站数量稀少的高原山区不能通过对站点数据插值获取准确的高空间分辨率的降水数据；周思儒等15基于8种降水数据集对青藏高原降水的时空变化进行研究，发现CMFD数据集模拟质量最高，相对误差为13.64%，模拟的各降水数据集均在地形较为复杂或降水量较低的干旱地区有着更高的相对误差，此外稀疏的实测站点也对数据集质量影响较大；Bai等16评估了12种降水数据在青藏高原的适用性，发现不同数据源的质量存在较大差异，多源降水数据模拟效果优于单源降水数据集。以上研究无论是从遥感数据还是站点数据插值的方法模拟高原山区的降水，都很少考虑高程梯度变化对降水空间分布的影响，使得高原山区的降水模拟精度普遍不高。针对此问题，提出了一种基于高程递推的反距离加权法，以此开展对高原山区的降水模拟研究。巴河是尼洋河最大支流，发源于念青唐古拉山脉南麓冰川，流域内高海拔地区发育有大量冰川，海拔5 700 m以上高度带，终年积雪，属于典型的高寒山区。因高寒山区高海拔地区长期缺乏有效的降水观测数据，其降水变化缺乏了解，因此，以巴河流域为研究区域，评估再分析数据集 CMFD、遥感数据GLADS和GSMaP三种类型的降水数据集的适用性，利用基于高程递推的反距离加权插值方法分析巴河流域降水空间分布特征，旨在为缺资料的高原山区水文过程模拟与水资源规划提供参考。1研究区与数据1.1研究区概况巴河是尼洋河最大支流，地处雅鲁藏布江中下游，流域面积约 4 200 km2，位于巴河桥站以北，占尼洋河流域面积的23.9%17，地理位置如图1所示。其南北宽约80 km，东西长约130 km，平均海拔4 000 m以上，流域落差近3 000 m，平均比降为3%，是平原区最大值的30倍左右，地形起伏变化较大18。流域内气象水文测站十分稀少，站点分布不均，缺乏全面、同步的观测数据。同时缺乏冰雪径流和土壤冻融的观测资料，对巴河流域的水文过程的模拟较为困难。1.2数据资料采用的降水数据分别为水文站实测降水数据、GLADS、GSMaP以及再分析数据集CMFD降水数据，各降水数据基本情况见表1。尼洋河流域在山区修建了一批用于中小河流山洪预19考虑高程的降水空间插值方法研究以巴河流域为例 鞠琴 刘小妮 连子旭 等警的雨量站，巴河桥水文站以上新建14站：霞拉、错高、错高1、罗巴、黑村、拉如、冲戈、雪卡、崩嘎、根田、英达、剥为、布如、朱拉。由于巴河流域15个观测站点数据在2014年以后，因维护困难出现观测中断，为便于与实测数据进行对比分析，评估降水数据的模拟效果，本文以2014年巴河流域观测降水数据进行讨论分析。2研究方法2.1基于高程递推的反距离加权法反距离平方加权法（Inverse Distance Weight，IDW）是计算降水空间分布的常用的空间插值方法之一，其计算原理是将离散的站点数据根据距离分配不同的权重，分配至整个流域面，从而模拟出流域的降水空间分布情况19。为探讨高原山区降水的复杂性及其随高程变化的规律，在传统的反距离加权空间插值方法的基础上，基于高程递推公式，提出了一种改进的反距离加权法，以便准确分析高程变化对降水空间分布的影响。首先假设巴河流域降水与高程存在一定的相关关系，即流域中不同高程对应的降水标准值符合统一的线性关系，表达式如下：Ph=a h+b（1）式中：Ph为h处的降水标准值，mm；h为计算点高程，m；a为降水与高程线性系数，mm/m；b为截距，mm。考虑到地形、植被等影响，流域各站点的降水标准值与实际降水仍有一定的误差，引入校正系数进行校正，校正系数与计算点水平位置相关，表达式如下：Pa=Ph C（2）式中：Pa为计算点实际降水量，mm；C为降水量校正系数。若已有控制站点的时段降水资料，则计算流域各点降水模拟值表达式如下：Pe=P0 CeC0(a he+ba h0+b)（3）式中：Pe与P0分别为计算点、控制站的时段降水量，mm；he与h0分别为计算点、控制站高程，m；Ce与C0分别为计算点、控制站降水量校正系数。高程递推法在计算中小流域山区降水量时，误差相对较小，而研究较大的流域时，仅根据降水与高程的线性关系进行递推，会导致递变速率难以适用，且多个递变速率计算繁杂。同时，为了削弱地形、植被、距离等其他因素的影响带来的误差，利用反距离加权法进行校正处理，计算降水量校正系数的表达式为：Ce=n=1Nd-2n Cnn=1Nd-2n（4）式中：Cn为已有资料参证站校正系数；dn为第n个参证站到计算点距离，m；N为参证站总数。2.2模拟评价指标-径流系数王建群等20研究了尼洋河流域 1997-2002年的降水径流水量平衡问题，发现流域雨量观测站设立不合理，观测降水量小于实测径流量，其中巴河流域的巴河桥水文站多年平均径流深1 482.2 mm，观测流域平均降水744.0 mm，多年降水径流存在严重的水量不平衡现象，其采用径流系数对流域内降水径流特征进行评估21。而巴河流域是尼洋河流域的子流域，位于藏东南地区，属于青藏高原地区，地形复杂，实测数据相对较少，所以本文也采用径流系数作为评价指标来评估降水的模拟精度，进一步探究改进的反距离加权方法对不同高程下的巴河流域降水空间分布的影响。3结果分析3.1降水数据比较根据GLADS、GSMaP与再分析数据集CMFD三种降水数据计算巴河流域年均降水量与空间分布情况，并与实测数据进行对比，以此反映流域的降水空间分布特征，分别见表2、图2。由表2可知，GLADS降水数据模拟的巴河流域降水量大部分小于500 mm，年均降水量最小为462.8 mm，与站点实测值误差较大，径流系数高达 3.03，明显低估了巴河流域的降水量。GSMaP降水数据模拟的巴河流域降水量空间变化幅度不大，河图1巴河流域地理位置图Fig.1 Map of Bahe River basin location表1降水数据基本情况表Tab.1 Basic information of precipitation data降水产品GLADSGSMaPCMFD观测数据全球陆面数据同化系统全球卫星降水制图项目中国区域高时空分辨率地面气象要素驱动数据集霞拉、错高、错高1、罗巴、黑村、拉如、冲戈、雪卡、崩嘎、根田、英达、剥为、布如、朱拉、巴河桥等水文站点空间分辨率0.250.250.10.10.10.1/数据来源http：/ldas.gsfc.nasa.govhttps：/sharaku.eorc.jaxa.jp/GSMaP_crest/index.htmlhttp：/水文年鉴20考虑高程的降水空间插值方法研究以巴河流域为例 鞠琴 刘小妮 连子旭 等谷地区降水量在 750900 mm 之间，流域年均降水量为 827.5 mm，径流系数为1.69，低于站点实测值。再分析数据集CMFD模拟的流域年均降水量为888.0 mm，径流系数为1.58，略高于巴河流域实测径流系数。总体来看，再分析数据集CMFD对巴河流域降水的模拟效果较 GSMaP、GLADS降水数据均有所提高，3种降水数据对巴河流域降水的模拟精度从大到小排序为CMFD GSMaP GLADS。由图2可知，3种降水数据模拟的降水空间分布差异性显著，计算的年均降水量与径流系数模拟结果与实际降水情况不符，无法准确反映高原山区高程对降水空间差异的影响，且卫星降水与再分析资料空间分辨率粗糙，对地形剧烈变化的高原山区进行模拟时存在较大的不确定性，不能得到准确的巴河流域实际的降水空间变化特征。3.2降水与高程相关性分析为探究巴河流域降水随高程的变化规律，更好的反映巴河流域降水的垂直分布情况，分析了巴河流域站点降水与高程之间的相关关系，结果如图3所示。由图3可知，巴河流域站点降水与高程的相关系数R为0.62，两者存在一定的线性相关关系，满足高程递推关系。大部分站点的降水量与高程呈正相关关系，少数站点距离趋势线较远，利用反距离加权法将这些距离变化结合到高程递推法中，用来校正降水插值过程中存在的误差，提出了一种基于高程递推的反距离加权法，以此得到流域的降水空间分布情况。3.3改进的反距离平方加权法适用性分析针对高寒山区降水观测缺乏、遥感数据反演精度偏低的情况，采用传统的反距离加权法与改进的反距离加权法分别计算巴河流域的降水空间分布情况，如图4所示。由图4（a）可知，利用反距离加权法对实测站点数据进行空间插值分析，得到研究区年均降水量为1 048.2 mm，径流系数为1.34，低估了巴河流域的降水量。传统的反距离加权法依赖于站点的海拔高度与站点数量，而巴河流域站点观测的区域代表性有限，说明对流域降水的低估很有可能来自高海拔地区。因此，仅利用反距离加权法无法展现巴河流域降水量随高程之间的变化，不能更好的解释降水径流的关系，难以得到较为准确的降水空间分布情况。这与李大伟等22分析秦岭中部山区降水变化，发现反距离加权法不能准确的反映降水空间变化显著的高海拔山区流域实际降水情况的研究结论一致。表2降水数据对巴河流域模拟结果统计表Tab.2 Simulation results of precipitation data for the Bahe River basin降水数据GLADSGSMaPCMFD年均降水量/mm462.8827.5888.0径流系数3.031.691.58空间分辨率0.250.250.10.10.10.1图2巴河流域降水量空间分布图Fig.2 Spatial distribution of precipitation in the Bahe River basin图3巴河流域站点降水与高程关系Fig.3 Precipitation and elevation relationships at sites in the Bahe River basin21考虑高程的降水空间插值方法研究以巴河流域为例 鞠琴 刘小妮 连子旭 等由图4（b）可知，受高程变化的影响，巴河流域年均降水量随海拔的升降而明显变化，改进的反距离加权法计算的降水空间分布与实际情况更加符合，平均降水量达到1 522 mm，径流系数为0.92。由全国径流系数分布图11可知，巴河流域多年平均径流系数处于0.61.0之间，计算结果满足流域水量平衡关系，此方法得到的降水径流关系较合理。表明改进的反距离加权法能准确的模拟流域不同高程的降水空间分布情况，在巴河流域适用性较好。3.4不同高程降水量变化根据改进的反距离加权法的计算结果，进一步分析高程对巴河流域降水空间分布的影响，由巴河流域降水与高程的累积频率分布曲线（图5）可知，降水与高程的累积频率分布曲线变化较为一致，降水量的频率分布主要集中在1 2001 800 mm之间，高程集中在4 0005 500 m之间，说明高程变化与降水空间分布密切相关。对巴河流域不同高程的降水量进行统计，结果见表3。巴河流域位于4 0005 000 m高程的流域面积占比最大，降水量相对贡献率最大，降水量占比为58.8%。降水的主要贡献高程位于4 0006 000 m，降水量占比高达89.6%，表明巴河流域降水量主要来源于高海拔地区的影响。4讨 论GLADS、GSMaP与再分析数据集CMFD降水数据本身空间分辨率较低23，遥感产品对低海拔地区的模拟表现优于高海拔地区，无法有效评估高海拔山区的降雨量，尤其是4 000 m海拔以上区域24。巴河流域实测站点属于较低海拔地区，仅有的15个站点均分布在4 000 m以下地区，往往对高原山区高海拔地区的降水量存在低估现象25。传统的反距离加权法以及常规的站点空间插值方法，其模拟精度主要受流域站点空间分布密度的限制26，忽略了高程、地形等因素对降水的影响，仅利用分布不均的15个站点观测数据不足以代表整个巴河流域的降水空间分布特征。由于巴河流域缺乏长系列的降水观测资料，流域年均降水量受高程变化的影响，存在明显的水量不平衡现象。提出的改进的反距离加权法不仅能够满足巴河流域水量平衡问题，使流域降水径流关系更加合理，模拟的降水空间分布特征更能反映不同高程对应的实际降水情况。模拟的时间序列取决于控制站点降水资料序列长度，不受其他因素干扰，有着较强的适应能力，可以为高原山区分布式水文模型提供前期降水资料的处理27。基于此方法，开展后续研究将考虑冰川融雪对径流的补给以及水汽通量对降水的影响，增加降水评估指标，进一步完善降水空间变化的适用性评估。5结 论为探讨水文资料短缺、地形复杂的巴河流域的降水空间分布情况，评估再分析数据集CMFD、遥感数据GLADS和GSMaP三种类型的降水数据集的适用性，同时，考虑高程对降水空间分布的影响，对比了反距离加权法与改进的反距离加权法对流域降水空间分布的模拟能力，分析了巴河流域不同高程的降水变化，主要结论如下。图5巴河流域高程与降水累积频率曲线Fig.5 Cumulative frequency of elevation and precipitation curves in the Bahe River basin图4巴河流域降水空间分布图Fig.4 Spatial distribution of precipitation in the Bahe River basin表3巴河流域不同高程的降水量统计表Tab.3 Precipitation at different elevations in the Bahe River basin高程/m6 000平均降水量/mm1 116.81 503.11 801.32 189.4面积占比/%13.859.826.10.3降水量贡献/mm154.2899.3470.55.4降水量占比/%10.158.830.80.4栅格个数7813 3851 4781422考虑高程的降水空间插值方法研究以巴河流域为例 鞠琴 刘小妮 连子旭 等（1）GLADS、GSMaP以及再分析数据集 CMFD 三种降水数据估算的巴河流域的年均降水量均以低估为主，径流系数分别为3.03、1.69、1.58，不满足流域水量平衡，未考虑不同高程对降水空间分布差异的影响，难以对地形起伏变化较大的高原山区进行模拟，均不适用于巴河流域。（2）基于高程递推的反距离加权法计算的平均降水量为 1 522 mm，径流系数为0.92，满足巴河流域的水量平衡，模拟结果优于传统的反距离加权法，在巴河流域适用性较好。（3）改进的反距离加权法能够细致地刻画巴河流域不同高程梯度的降水情况，巴河流域降水量主要来源于高海拔地区，其中 4 0006 000 m 海拔地区降水量占比高达 89.5%，4 000 5 000 m海拔的流域面积占比最大。参考文献：1 JU Q，ZHANG R，WANG G，et al.Simulating the freezing-thawing processes based on MODIS data in the Three-River source region，China J.Frontiers in Environmental Science，2022，10（1）：1 559.2 罗上学，张美玲，聂雅梅，等.基于 CEEMDAN-LSTM 模型的郑州市月降水量预测 J.水利规划与设计，2022，220（2）：45-50.LUO S X，ZHANG M L，NIE Y M，et al.Prediction of monthly precipitation in Zhengzhou based on CEEMDAN-LSTM model J.Water Resources Planning and Design，2022，220（2）：45-50.3 何晓波，叶柏生，丁永建.青藏高原唐古拉山区降水观测误差修正分析 J.水科学进展，2009，20（3）：403-408.HE X B，YE B S，DING Y J.Analysis of error correction of precipitation observation in Tanggula Mountains of Qinghai-Tibet plateau J.Advances in Water Science，2009，20（3）：403-408.4 ZHANG G X，SU X L，AYANTOBO O，et al.Spatial interpolation of daily precipitation based on modified ADW method for gauge-scarce mountainous regions：A case study in the Shiyang River basinJ.Atmospheric Research，2021，247：105 167.5 任宾宾，王平，邱绍柱，等.基于Co-kriging的呼伦贝尔市2018年降水量空间分布特征模拟及插值精度分析 J.西南大学学报（自然科学版），2021，43（11）：162-171.REN B B，WANG P，QIU S Z，et al.Simulation of spatial distribution characteristics and interpolation accuracy of precipitation in Hulunbuir City in 2018 based on co-kriging J.Journal of Southwest University（Natural Science Edition），2021，43（11）：162-171.6 符静，易臻照，赵莎.基于GIS空间插值的降水分布模拟方法比较研究 J.中国农村水利水电，2021（1）：94-97+104.FU J，YI Z Z，ZHAO S.Comparative research on precipitation distribution simulation methods based on GIS spatial interpolation J.China Rural Water and Hydropower，2021（1）：94-97+104.7 杨家亮，崔璨，李淑雅，等.湖南省多尺度降水序列混沌特性空间分异 J.水利规划与设计，2022，226（8）：33-37.YANG J L，CUI C，LI S Y，et al.Spatial differentiation of chaotic characteristics of multi-scale precipitation series in Hunan provinceJ.Water Resources Planning and Design，2022，226（8）：33-37.8 文迁，谭国良，罗嗣林.降水分布受地形影响的分析 J.水文，1997（增1）：64-66.WEN Q，TAN G L，LUO S L.Analysis of precipitation distribution affected by terrain J.Journal of China Hydrology，1997（S1）：64-66.9 陈裕迪，王洁，马燮铫，等.五种降水数据产品在元江红河流域的适应性评估 J.水电能源科学，2022，40（3）：9-12+16.CHEN Y D，WANG J，MA X C，et al.Adaptability assessment of five precipitation data products in the Yuanjiang-Honghe River basin J.Water Resources and Power，2022，40（3）：9-12+16.10 朱浩楠，刘晓冉，李永华，等.考虑地形的空间插值算法在复杂下垫面地区气温和降水精细化插值的评估 J.气象，2020，46（5）：655-665.ZHU H N，LIU X R，LI Y H，et al.Evaluation of refined interpolation of temperature and precipitation in complex underlying surface areas by spatial interpolation algorithm considering terrain J.Meteorological Monthly，2020，46（5）：655-665.11 李豪，刘涛，徐精文.基于混合地理加权回归与克里格的区域降水量空间插值方法J.中国农业气象，2018，39（10）：674-684.LI H，LIU T，XU J W.Spatial Interpolation method of regional precipitation based on mixed geographically weighted regression and kriging J.Chinese Journal of Agrometeorology，2018，39（10）：674-684.12 李慧晴，叶爱中.基于地形加权的降水空间插值方法研究 J.武汉大学学报（工学版），2021，54（1）：28-37.LI H Q，YE A Z.Spatial interpolation method of precipitation based on topographic weighting J.Engineering Journal of Wuhan University，2021，54（1）：28-37.13 李琼，魏加华，安娟，等.基于地形因子的TRMM 3B43降水数据在黄河源区的融合校正研究 J.应用基础与工程科学学报，2018，26（6）：1 147-1 163.LI Q，WEI J H，AN J，et al.Fusion correction of TRMM 3B43 precipitation data based on terrain factors in the source region of the Yellow River J.Journal of Basic Science and Engineering，2018，26（6）：1 147-1 163.14 刘田，阳坤，秦军，等.青藏高原中、东部气象站降水资料时间序列的构建与应用 J.高原气象，2018，37（6）：1 449-1 457.LIU T，YANG K，QIN J，et al.Construction and applications of time series of monthly precipitation at weather stations in the central and Eastern Qinghai-Tibetan plateauJ.Plateau Meteorology，2018，37（6）：1 449-1 457.15 周思儒，信忠保.基于多源数据的近40年青藏高原降水量时空变化研究 J/OL.长江科学院院报：1-11.2023-03-01.http：/ S R，XIN Z B.Study on temporal and spatial variation of precipitation in Qinghai-Tibet plateau in recent 40 years based on multi-source dataJ/OL.Journal of Yangtze River Scientific Research Institute：1-11.2023-03-01.http：/ BAII L，WEN Y，SHI C，et al.Which precipitation product works best in the Qinghai-T