一种基于水碳耦合方法的全球地表蒸散发估算方法技术

本发明专利技术的一种基于水碳耦合方法的全球地表蒸散发估算方法包括以下步骤：首先搜集预处理地面通量站数据、全球SIF数据、叶面积指数数据以及再分析气象数据集；再基于水碳耦合原理相关的植物生理学公式，建立地表蒸散发与SIF之间的关系式；然后，通过将站点尺度数据输入至模型中，并优化得到站点内蒸散发估算参数；随后，通过环境变量、站点蒸散发估算参数和两种不同回归模型，计算全球蒸散发参数；最后，将全球气象再分析产品、全球SIF与全球蒸散发估算参数代入到模型中，得到全球地表蒸散发。估算的全球地表蒸散发产品可为流域水文模拟和水资源规划提供重要且可操作性强的参考依据。水资源规划提供重要且可操作性强的参考依据。水资源规划提供重要且可操作性强的参考依据。

【技术实现步骤摘要】
一种基于水碳耦合方法的全球地表蒸散发估算方法


[0001]本专利技术涉及遥感蒸散发估算领域，涉及一种基于水碳耦合方法的全球地表蒸散发估算方法。

技术介绍

[0002]地表蒸散发是全球和区域水循环的重要组分，是能将水循环、能量循环与碳循环连接起来的关键要素。准确估算地表蒸散发对理解水碳循环和能量平衡至关重要。当前，大尺度蒸散发估算模型依赖气象数据和植被物候的遥感产品，但是这些蒸散发模型难以准确表示受生物过程调节的生态因素。因此，当前蒸散发模型的改进重点主要聚焦于对蒸散发过程中生物过程的参数化表征。
[0003]近年来，日光诱导叶绿素荧光(SIF)作为植被光合作用的理想探针，在估算蒸散发方面具有很大的潜力，可以刻画蒸散发的重要组成部分——植物蒸腾量。随着卫星遥测和机器学习技术的快速发展，基于遥感重构的叶绿素荧光产品具有较高的时空精度。这类重构数据的出现，补足了目前地面观测叶绿素荧光数据不足的缺陷，开始被广泛应用于全球生态学和水文学领域。
[0004]目前，叶绿素荧光与地表蒸散发之间的内在物理联系仍然不清，现有研究普遍存在以下问题：(1)对叶绿素荧光与地表蒸散发的处理经验化；(2)基于叶绿素荧光的蒸散发估算需要依靠大量植物生理学参数，且参数过多，部分参数无法在大尺度上获取，需要引入额外的假设条件，可能会引入较大的不确定性。总体而言，现有研究未能充分融合遥感数据与物理机制估计地表蒸散发。

技术实现思路

[0005]本专利技术要解决的技术问题在于，针对现有技术中的缺陷，提供一种基于水碳耦合方法的全球地表蒸散发估算方法。
[0006]本专利技术解决其技术问题所采用的技术方案是：
[0007]一种基于水碳耦合方法的全球地表蒸散发估算方法，该方法包括以下步骤：
[0008]步骤1、搜集地面通量观测数据、气象再分析数据、MODIS地表叶面积指数LAI数据以及全球重建日光诱导叶绿素荧光数据；
[0009]步骤2、在站点尺度构建叶绿素荧光与蒸散发的函数关系；
[0010]步骤3、在各通量站中对叶绿素荧光与蒸散发之间的关系进行参数率定，求得各通量站下蒸散发的估算参数；
[0011]步骤4、构建不同植被气象要素与蒸散发估算参数的函数关系，构建蒸散发估算参数的空间分布数据库；
[0012]步骤5、将获取的叶绿素荧光、全球气象数据以及全球蒸散发估算参数代入至步骤2构建的模型中，进行地表蒸散发的估算。
[0013]进一步地，所述步骤1包括如下子步骤：
[0014]步骤1.1、基于ERA5
‑
Land小时尺度气象再分析数据，提取某一时间分辨率的露点温度T
d
、气温T
a
、饱和水汽压亏缺VPD、相对湿度RH、年均气温MAT、年均降水MAP以及净辐射R
n
，同时将上述数据的空间分辨率通过空间插值进行统一；
[0015]步骤1.2、基于现有年均降水以及年均潜在蒸散发数据计算干旱指数DI；
[0016]步骤1.3、获取陆地通量网数据，筛选出具有至少连续两年有效观测的站点，并提取与步骤1.1相同时间分辨率的站点饱和水汽压亏缺VPD、CO2浓度C
a
、气温T
a
和站点蒸发潜热ET数据；
[0017]步骤1.4、从MODIS数据集中获取地表叶面积指数LAI数据，并将空间分辨率转化为与步骤1.1中相同的空间分辨率；
[0018]步骤1.5、获取一定时间段全球重建SIF数据CSIF，并与陆地通量网数据站点进行匹配，得到站点尺度CSIF的时间序列。
[0019]进一步地，步骤2的具体方法为：
[0020]步骤2.1.根据双源地表蒸散发模型，地表蒸散发可以被分成两个部分：土壤蒸发以及植被蒸腾；由此，地表蒸散发可以被表示为：
[0021]ET＝T
r
+E
s
[0022]其中，E
s
为土壤蒸发，T
r
为植被蒸腾，ET为地表蒸散发；
[0023]步骤2.2.基于水碳耦合原理建立植被蒸腾T
r
与叶绿素荧光SIF的关系；
[0024]步骤2.3.构建VPD和T
a
对蒸散发的控制函数f
VPD
和f
Ta
来计算土壤蒸发，其中f
VPD
和f
Ta
均为二次函数；
[0025]步骤2.4.将步骤2.2和步骤2.3中计算得到的植被蒸腾和土壤蒸发相加得到地表蒸散发。
[0026]进一步地，GPP与SIF关系可用以下方程表达描述：
[0027]GPP＝αSIF+β
[0028]从Fick定律角度出发，冠层尺度GPP可以用以下方程描述，
[0029]GPP＝g
c
×
(C
a
‑
C
i
)
[0030]根据水碳耦合理论，H2O分子流出植物气孔与CO2公用通道，因此可以根据Fick定律描述蒸腾：
[0031]T
r
＝g
h
×
(q
i
‑
q
a
)
[0032]q
i
‑
q
a
可用以下公式表征：q
i
‑
q
a
＝VPD/P
a
[0033]通过以上方程，可求得植被蒸腾量：
[0034][0035]根据最优气孔导度理论，
[0036][0037]统一量纲后，蒸腾可表达为：
[0038][0039]式中：Γ为光合作用补偿点，由以下方程确定，
[0040][0041]因此，
[0042]其中，α和β为方程参数，g
c
为植被对CO2的导度，C
a
与C
i
分别为大气和植物气孔内部的CO2浓度，g
h
为植物气孔对H2O分子的导度，且存在g
c
＝g
h
/1.6；q
i
为植物叶片内H2O浓度，q
a
为大气内H2O浓度，P
a
为大气压强；，λ
cf
为边际水分利用效率，为常数；O2为叶内部O2分压。
[0043]进一步地，土壤蒸发采用以下公式计算，
[0044][0045]式中：Δ为饱和水汽压曲线的斜率，γ为干湿表常数，k
A
为消光系数，其值等于0.6，R
h
为相对湿度，LAI为叶面积指数，R
n
为地表净辐射，G为地表热通量；f
VPD
和f
Ta
为VPD和T
a
对蒸散发的控制函数，其表达式如下：
[0046]f
VPD
＝1
...

【技术保护点】

【技术特征摘要】
1.一种基于水碳耦合方法的全球地表蒸散发估算方法，其特征在于，该方法包括以下步骤：步骤1、搜集地面通量观测数据、气象再分析数据、MODIS地表叶面积指数LAI数据以及全球重建日光诱导叶绿素荧光数据；步骤2、在站点尺度构建叶绿素荧光与蒸散发的函数关系；步骤3、在各通量站中对叶绿素荧光与蒸散发之间的关系进行参数率定，求得各通量站下蒸散发的估算参数；步骤4、构建不同植被气象要素与蒸散发估算参数的函数关系，构建蒸散发估算参数的空间分布数据库；步骤5、将获取的叶绿素荧光、全球气象数据以及全球蒸散发估算参数代入至步骤2构建的模型中，进行地表蒸散发的估算。2.根据权利要求1所述的一种基于水碳耦合方法的全球地表蒸散发估算方法，其特征在于，所述步骤1包括如下子步骤：步骤1.1、基于ERA5
‑
Land小时尺度气象再分析数据，提取某一时间分辨率的露点温度T
d
、气温T
a
、饱和水汽压亏缺VPD、相对湿度RH、年均气温MAT、年均降水MAP以及净辐射R
n
，同时将上述数据的空间分辨率通过空间插值进行统一；步骤1.2、基于现有年均降水以及年均潜在蒸散发数据计算干旱指数DI；步骤1.3、获取陆地通量网数据，筛选出具有至少连续两年有效观测的站点，并提取与步骤1.1相同时间分辨率的站点饱和水汽压亏缺VPD、CO2浓度C
a
、气温T
a
和站点蒸发潜热ET数据；步骤1.4、从MODIS数据集中获取地表叶面积指数LAI数据，并将空间分辨率转化为与步骤1.1中相同的空间分辨率；步骤1.5、获取一定时间段全球重建SIF数据CSIF，并与陆地通量网数据站点进行匹配，得到站点尺度CSIF的时间序列。3.根据权利要求1所述的一种基于水碳耦合方法的全球地表蒸散发估算方法，其特征在于，步骤2的具体方法为：步骤2.1.根据双源地表蒸散发模型，地表蒸散发可以被分成两个部分：土壤蒸发以及植被蒸腾；由此，地表蒸散发可以被表示为：ET＝T
r
+E
s
其中，E
s
为土壤蒸发，T
r
为植被蒸腾，ET为地表蒸散发；步骤2.2.基于水碳耦合原理建立植被蒸腾T
r
与叶绿素荧光SIF的关系；步骤2.3.构建VPD和T
a
对蒸散发的控制函数f
VPD
和f
Ta
来计算土壤蒸发，其中f
VPD
和f
Ta
均为二次函数；步骤2.4.将步骤2.2和步骤2.3中计算得到的植被蒸腾和土壤蒸发相加得到地表蒸散发。4.根据权利要求3所述的一种基于水碳耦合方法的全球地表蒸散发估算方法，其特征在于，GPP与SIF关系可用以下方程表达描述：GPP＝αSIF+β
从Fick定律角度出发，冠层尺度GPP可以用以下方程描述，GPP＝g
c
×
(C
a
‑
C
i
)根据水碳耦合理论，H2O分子流出植物气孔与CO2公用通道，因此可以根据Fick定律描述蒸腾：T
r
＝g
h
×
(q
i
‑
q
a
)q
i
‑
q
a
可用以下公式表征：q
i
‑
q
a
＝VPD/P
a
通过以上方程，可求得植被蒸腾量：根据最优气孔导度理论，统一量纲后，蒸腾可表达为：式中：Γ为光合作用补偿点，由以下方程确定，因此，其中，α和β为方程参数，g
c
为植被对CO2的导度，C
a
与C
i
分别为大气和植物气孔内部的CO2浓度，g
h
为植物气孔对H2O分子的导度，且存在g
c
＝g
h
/1.6；q
i
为植物叶片内H2O浓度，q
a
为大气内H2O浓度，P
a
为大气压强；，λ
cf
为边际水分利用效率，为常数；O2为叶内部O2分压。5.根据权利要求4所述的一种基于水碳耦合方法的全球地表蒸散发估算方法，其特征在于，土壤蒸发采用以下公式计...

【专利技术属性】
技术研发人员：张权，周恺，刘轩旗，尹家波，
申请(专利权)人：武汉大学，
类型：发明
国别省市：