[发明专利]大尺度区域实际蒸散发量的估算方法

权利要求书

1.大尺度区域实际蒸散发量的估算方法，其特征在于，所述实际蒸散发量是根据水分平衡公式计算得到；

公式1：AET＝P-R-ΔH；其中，

P表示研究区计算时段内的平均降水量，单位mm；

R是计算时段内的等效径流深，单位mm；

AET是计算时段内区域的实际蒸散发量，单位mm；

ΔH是计算时间段内的区域水量的蓄变量，即地下水和土壤水的变化量，对于半干旱厚包气带区域来说，ΔH可以近似看作是区域土壤水的变化量，单位mm；

所述ΔH是通过GRACE重力卫星数据反演得到；其中，由重力场模型球谐系数表示地球表面质量变化量的模型为：

公式2:

其中，Δh(θ,φ)为等效水高；α为地球半径；θ和φ分别为待定点的余纬和经度；为规格化的缔合Legendre多项式；和则是一组无量纲的球谐系数的变化量；l和m分别表示阶数和次数；ρ^avc和ρ^wat分别为地球的平均密度和水的密度；k_l为一阶勒夫数；

公式3：

公式4：

其中，W_n为高斯(Gaussian)平滑核函数，γ为平滑半径。

2.如权利要求1所述的大尺度区域实际蒸散发量的估算方法，其特征在于，所述土壤层为黄土高原的厚包气带土壤层。

3.如权利要求2所述的大尺度区域实际蒸散发量的估算方法，其特征在于，所述土壤层的层厚一般为50m-200m。

4.如权利要求2所述的大尺度区域实际蒸散发量的估算方法，其特征在于，所述黄土高原的降水数据由黄土高原的300个气象观测站点获取。

5.如权利要求2所述的大尺度区域实际蒸散发量的估算方法，其特征在于，所述径流深数据是从ERA-Interim再分析数据中的径流模拟结果中提取。

6.如权利要求5所述的大尺度区域实际蒸散发量的估算方法，其特征在于，所述ERA-Interim的径流数据经过所述黄土高原的20个典型无人类活动干扰站点的径流数据的实际验证，其验证结果为所述ERA-Interim包括的所述黄土高原的径流数据符合实际。

7.如权利要求6所述的大尺度区域实际蒸散发量的估算方法，其特征在于，所述20个典型无人类活动的水文站点分别位于延川、绥德、张家山、秦安、河津、刘家河、庆阳、张家山、正宁、龙门、大村、子长、吴堡、大宁、华亭、耀县、志丹、延安、万荣及柳林。