CN113156542B - 一种基于闪烁法测量区域蒸散总水量的方法 - Google Patents

一种基于闪烁法测量区域蒸散总水量的方法 Download PDF

Info

Publication number
CN113156542B
CN113156542B CN202110271725.4A CN202110271725A CN113156542B CN 113156542 B CN113156542 B CN 113156542B CN 202110271725 A CN202110271725 A CN 202110271725A CN 113156542 B CN113156542 B CN 113156542B
Authority
CN
China
Prior art keywords
unit
formula
air
heat flux
vapor pressure
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202110271725.4A
Other languages
English (en)
Other versions
CN113156542A (zh
Inventor
刘聪
吴东丽
邵长亮
朱永超
杨大生
马剑哲
徐鑫
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
CMA Meteorological Observation Centre
Original Assignee
CMA Meteorological Observation Centre
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by CMA Meteorological Observation Centre filed Critical CMA Meteorological Observation Centre
Priority to CN202110271725.4A priority Critical patent/CN113156542B/zh
Publication of CN113156542A publication Critical patent/CN113156542A/zh
Application granted granted Critical
Publication of CN113156542B publication Critical patent/CN113156542B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01W—METEOROLOGY
    • G01W1/00—Meteorology
    • G01W1/02—Instruments for indicating weather conditions by measuring two or more variables, e.g. humidity, pressure, temperature, cloud cover or wind speed
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
    • G01K17/00—Measuring quantity of heat
    • G01K17/02—Calorimeters using transport of an indicating substances, e.g. evaporation calorimeters
    • G01K17/025—Calorimeters using transport of an indicating substances, e.g. evaporation calorimeters where evaporation, sublimation or condensation caused by heating or cooling, is measured

Landscapes

  • Engineering & Computer Science (AREA)
  • Environmental & Geological Engineering (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Atmospheric Sciences (AREA)
  • Biodiversity & Conservation Biology (AREA)
  • Ecology (AREA)
  • Environmental Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Combustion & Propulsion (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

本发明公开了一种基于闪烁法测量区域蒸散总水量的方法。方法是通过大口径闪烁仪的路径积分效应测得区域大气的波文比、空气密度、空气比热容、大气折射率结构参数、平均感热通量，再通过下垫面热平衡方程计算得到整体区域的平均潜热通量，最后得出区域蒸散总水量。本方法可用于对观测区域蒸散总水量值的计算。这种通过大口径闪烁仪的路径积分效应测量地表热通量，很好地解决了通量观测中单点观测设备由于下垫面不均一带来的区域总水量不能正确统计问题，可以和卫星数据进行订正，得到更大尺度的整体蒸散数据。

Description

一种基于闪烁法测量区域蒸散总水量的方法
技术领域
本发明涉及区域水资源的测量技术，尤其是涉及一种基于闪烁法测量区域蒸散总水量的方法。
背景技术
传统的地表蒸散量观测方法多是单点观测，只能提供特定位置的观测数据。由于下垫面的不均一，观测受降水、土壤、水文以及不同植被等参数变化影响严重。传统的观测方法只能通过空间插值的方法由点尺度外推为面尺度。由于下垫表面的复杂性使得该方法在进行空间尺度扩展时会产生较大误差。需要密集布点分布观测，这种情况下，价格太高，位置不同数据处理相对困难。因此大范围的观测平均值很有必要。而且得到区域蒸散量数据之后，可以和卫星数据进行订正，得到更大尺度的整体蒸散数据。
发明内容
为了解决现有技术中无法有效的对区域蒸散总水量整体计算评估的问题，本发明提供一种基于闪烁法测量区域蒸散总水量的方法。
本发明采取的技术方案是：一种基于闪烁法测量区域蒸散总水量的方法，其特征在于，所述方法是通过大口径闪烁仪的路径积分效应测得区域大气的波文比、空气密度、空气比热容、大气折射率结构参数、平均感热通量，再通过下垫面热平衡方程计算得到整体区域的平均潜热通量，最后得出区域蒸散总水量。具体步骤如下：
第一步：计算被测量区域大气的波文比、空气密度、空气比热容；
①、波文比β：
波文比β为相对湿度对闪烁影响因素的校正因数，根据气象观测数据推导计算得出，如下公式所示：
Figure BDA0002974418730000011
式(1)中，γ为干湿球常数，γ＝a0×P，a0为常数，通常取0.665×10-3，P为大气压力，单位：Pa；ΔT为不同观测层间的空气温度差，单位：K；ΔE为不同观测层间的实际水汽压差，单位：Pa；ΔT和ΔE分别有以下关系式：
ΔT＝T1-T2------------(2)
ΔE＝E1-E2------------(3)
式(2)中，T1为上层空气温度，T2为下层空气温度；
式(3)中，E1为上层水汽压，E2为下层水汽压；
上层水汽压E1通过上层空气相对湿度RH1、上层饱和水汽压Ea1和上层空气温度t1进行换算，换算公式如下：
E1＝Ea1×RH1×0.01-----------(4)
Figure BDA0002974418730000021
t1＝T1-273.15-----------------------(6)
式(4)中，E1为上层实际水汽压，单位：Pa；RH1为上层空气相对湿度，单位：％；Ea1为上层饱和水汽压，单位：Pa；
式(5)中，t1为上层空气温度，单位：℃；a1和a2为常数，a1＝6.11、a2＝17.27；
式(6)中，T1为上层空气温度，单位：K；
下层水汽压E2通过下层空气相对湿度RH2、下层饱和水汽压Ea2和下层空气温度t2进行换算，换算公式如下：
E2＝Ea2×RH2×0.01-----------(7)
Figure BDA0002974418730000022
t2＝T2-273.15-----------------------(9)
式(7)中，E2为下层实际水汽压，单位：Pa；RH2为下层空气相对湿度，单位：％；Ea2为下层饱和水汽压，单位：Pa；
式(8)中，t2为下层空气温度，单位：℃；a1和a2为常数，a1＝6.11、a2＝17.27；
式(9)中，T2为下层空气温度，单位：K。
②、空气密度ρ：
空气密度ρ的计算公式如下：
Figure BDA0002974418730000023
式(10)中，ρ为空气密度，单位：Kg/m3；RH1为上层空气相对湿度，单位：％；Ea1为上层对应的饱和水汽压，单位：Pa；T1为上层空气温度，单位：K，P为大气压力，单位：Pa；A和B为常数，A＝3.48、B＝0.378。
③、空气比热容CP：
空气比热容CP，单位J/(kg·K)，计算公式如下：
Cp＝c-dT1-------------(11)
式(11)中，T1为上层空气温度，单位：K；c和d为常数，c＝14.439、d＝0.9504×10-3。
第二步：计算温度结构参数
Figure BDA0002974418730000031
通过大口径闪烁仪测量得到的大气折射率结构参数
Figure BDA0002974418730000032
及相关的上层空气温度T1、大气压力P以及波文比β，计算温度结构参数
Figure BDA0002974418730000033
单位：K2m-2/3，计算公式如下：
Figure BDA0002974418730000034
式(12)中，β为波文比；T1为上层空气温度，单位：K；P为大气压力，单位：Pa；p0为常数，p0＝-0.78×10-6。
第三步：计算平均感热通量H
平均感热通量H计算公式如下：
Figure BDA0002974418730000035
式(13)中，H为平均感热通量，单位：W/m2；b是经验常数，b＝0.48；ρ为空气密度，单位：Kg/m3；CP为空气比热容，单位：J/kg/K；ZLAS为设备安装高度，单位：m；d为大口径闪烁仪安装位置相对地面的零值位移，单位：m；g为重力加速度，单位：m/s2；T1为上层空气温度，单位：K。
第四步：计算平均潜热通量LvE
平均潜热通量LvE计算公式如下：
LvE＝Q*-H+Gs------------(14)
式(14)中，LvE为平均潜热通量，单位：W/m2；Q*为净辐射，单位：W/m2；H为平均感热通量，单位：W/m2；Gs为土壤热通量，单位：W/m2。
第五步：数据质量控制
对计算出的平均潜热通量进行数据剔除和平滑、去噪。
第六步：计算区域蒸散总水量ET
区域蒸散总水量ET的计算公式如下：
Figure BDA0002974418730000036
式(15)中，ET为m分钟内区域蒸散总水量，单位：kg/m2；Lv为蒸散潜热，即蒸散1KG水所需的能量，单位：J/kg；至此由ET得出整片区域的蒸散总水量。
一个地区的水资源中，蒸散量是十分重要的指标之一，反映当地水资源的收支平衡。关系到生产、生活的方方面面。从而就要对蒸散量进行测量和计算。下垫表面是大气的底边界，研究大气层大气必须讨论不同下垫表面对大气层的影响。除去动量通量传输外，下垫表面对大气边界层的另两个强制因素是称之为感热通量和潜热通量通量。感热通量表示大气热传导的热量大小，潜热通量是指大气中水汽产生相变所需要的热量。也就是水汽蒸散所需要的热量。
闪烁法测量的大气湍流是基于大气光学性质起影响作用的大气密度起伏引起的大气折射率起伏所对应的光学湍流。由于大气密度起伏主要由温度起伏决定，一般大气光学湍流由大气温度场的起伏特性决定。通过大气结构参数可以反演得到区域平均感热通量。
本发明所产生的有益效果是：本方法可用于对观测区域蒸散总水量值的计算。这种通过大口径闪烁仪的路径积分效应测量地表热通量，很好地解决了通量观测中单点观测设备由于下垫面不均一带来的区域总水量不能正确统计问题，可以和卫星数据进行订正，得到更大尺度的整体蒸散数据。
附图说明
图1为本发明实施过程的流程图；
图2为本发明根据自动气象站的观测数据和大口径闪烁仪的观测数据计算得到温度结构参数
Figure BDA0002974418730000041
的数值曲线图；
图3为本发明实测一天的净辐射、土壤热通量、平均感热通量和平均潜热通量的测量时间序列及实测计算数据图；
图4是本发明实施例计算出一天时间的蒸散总水量趋势图。
具体实施方式
以下结合附图和实施例对本发明作进一步说明。
实施例-如图1所示：
实测值：
t1为上层空气温度22.4℃，t2为下层空气温度22℃；
T1为上层空气温度295.55K，T2为下层空气温度295.15K；
RH1为上层空气相对湿度62％，RH2为上层空气相对湿度65％；
P为大气压力101200pa；
Figure BDA0002974418730000042
为大气折射率结构参数1.34×10-1；
ZLAS为设备安装高度(20m)；
d为零值位移(2m)；
g为重力加速度(9.8m/s2)；
Q*为净辐射(318W/m2)；
Gs为土壤热通量(0W/m2)。
第一步，使用自动气象站的数据计算波文比，空气密度、空气比热容。
①波文比：
ΔT为不同观测层间的空气温度差(K),由公式(2)可得：
ΔT＝T1-T2＝295.55-295.15＝0.4K；
Ea1为上层饱和水汽压(pa)，由公式(5)可得：
Figure BDA0002974418730000051
式中，a1和a2为常数，a1＝6.11、a2＝17.27；
E1为上层水汽压(pa)，由公式(4)可得：
E1＝Ea1×RH1×0.01＝2262.001×62×0.01＝1402.44pa；
Ea2为下层饱和水汽压(Pa)，由公式(8)可得：
Figure BDA0002974418730000052
E2为下层水汽压(pa)，由公式(7)可得：
E2＝Ea2×RH2×0.011＝2213.596×65×0.01＝1438.838pa；
ΔE为不同观测层间的实际水汽压差(Pa)，由公式(3)可得：
ΔE＝E1-E2＝1402.44-1438.838＝-36.3973pa；
β为波文比，由公式(1)可得：
Figure BDA0002974418730000053
式中，γ为干湿球常数。γ＝a0×P，a0为常数，通常取0.665×10-3；P为大气压力(Pa)。
②空气密度ρ：
由公式(10)可得：
Figure BDA0002974418730000054
Figure BDA0002974418730000061
式中ρ为空气密度(Kg/m3)，RH1为上层空气相对湿度(％)，Ea1为上层对应的饱和水汽压(Pa)，T1为上层空气温度(K)，P为大气压力(Pa)，A和B为常数，A＝3.48、B＝0.378。
③空气比热容CP：
由公式(11)可得：
Cp＝c-dT1＝14.439-0.9504×10-3×295.55＝14.158J/(kg·K)
式中T1为上层空气温度(K)，c和d为常数，c＝14.439、d＝0.9504×10-3。
第二步，根据设备测量数据计算温度结构参数
Figure BDA0002974418730000062
通过大口径闪烁仪测量得到的大气折射率结构参数
Figure BDA0002974418730000063
及相关的上层空气温度(T1)、大气压力(P)以及第一步计算的波文比β，可以计算测量区域温度结构参数
Figure BDA0002974418730000064
由公式(12)可得：
Figure BDA0002974418730000065
上式中：β为波文比，T1为上层空气温度(K)，P为气压(Pa)，p0为常数，p0＝-0.78×10-6。
第三步，计算平均感热通量H
由公式(13)可得：
Figure BDA0002974418730000066
式中b是经验常数：b＝0.48，ρ为空气密度(Kg/m3)，CP为空气比热容(J/kg/K)，ZLAS为设备安装高度(20m)，d为零值位移(2m)，g为重力加速度(9.8m/s2)，T1为上层空气温度(K)。
第四步，计算平均潜热通量LvE
由公式(14)可得：
LvE＝Q*-H+Gs
＝318-1.288+0
＝316.7W/m2
式中Q*为净辐射(318W/m2)，H为平均感热通量(1.288W/m2)，LvE为平均潜热通量(W/m2)，Gs为土壤热通量(W/m2)。
第五步，数据质量控制
主要工作是对计算出的平均潜热通量进行数据剔除和平滑、去噪。
第六步，计算区域蒸散总水量ET
由公式(15)可得：
Figure BDA0002974418730000071
此处m取值为60，也就是计算小时内的蒸散量。其中Lv≈2.45×106J/kg，得到ET区域蒸散总水量，数据见表1。表1列举出使用本方法计算得到的2020年9月11日6：00—17：00的平均潜热通量和蒸散总水量的数据。
表1
数据时间 LvE(w/m<sup>2</sup>) ET(kg/m<sup>2</sup>)
2020-9-11_6：00 13.10618 0.005349
2020-9-11_7：00 5148.977 2.101623
2020-9-11_8：00 11993.19 4.89518
2020-9-11_9：00 12825.86 5.235044
2020-9-11_10：00 14553.94 5.940383
2020-9-11_11：00 13244.27 5.405823
2020-9-11_12：00 12964.51 5.291637
2020-9-11_13：00 11489.8 4.689716
2020-9-11_14：00 9508.184 3.880891
2020-9-11_15：00 8318.061 3.395127
2020-9-11_16：00 6680.763 2.726842
2020-9-11_17：00 967.941 0.395078
从图2可以看出：2020年9月11日6：00—18：00温度结构参数
Figure BDA0002974418730000072
(单位：K2m-2/3)变化曲线，其分别在6：00和18：00，即日出和日落时出现明显的下降。主要原因在于在日出时分附近，大气从稳定状态转变为不稳定状态；日落时分附近，大气从不稳定状态转变为稳定状态。
图3是一天的净辐射(RN)(单位：w/m2)、土壤热通量(HFP)(单位：w/m2)、平均感热通量(H)(单位：w/m2)和平均潜热通量(LvE)(单位：w/m2)测量时间序列实测计算数据。
图4表示2020年9月11日6：00—17：00的蒸散总水量(单位：kg/m2)的趋势图可以看出，在本实施例中一天的中午时间蒸散水量是最大的。

Claims (1)

1.一种基于闪烁法测量区域蒸散总水量的方法，其特征在于，所述方法是通过大口径闪烁仪的路径积分效应测得区域大气的波文比、空气密度、空气比热容、大气折射率结构参数、平均感热通量，再通过下垫面热平衡方程计算得到整体区域的平均潜热通量，最后得出区域蒸散总水量；具体步骤如下：
第一步：计算被测量区域大气的波文比、空气密度、空气比热容；
①、波文比β：
波文比β为相对湿度对闪烁影响因素的校正因数，根据气象观测数据推导计算得出，如下公式所示：
Figure FDA0003553041940000011
式(1)中，γ为干湿球常数，γ＝a0×P，a0为常数，a0＝0.665×10-3，P为大气压力，单位：Pa；ΔT为不同观测层间的空气温度差，单位：K；ΔE为不同观测层间的实际水汽压差，单位：Pa；ΔT和ΔE分别有以下关系式：
ΔT＝T1-T2-------(2)
ΔE＝E1-E2-------(3)
式(2)中，T1为上层空气温度，T2为下层空气温度；
式(3)中，E1为上层水汽压，E2为下层水汽压；
上层水汽压E1通过上层空气相对湿度RH1、上层饱和水汽压Ea1和上层空气温度t1进行换算，换算公式如下：
E1＝Ea1×RH1×0.01------(4)
Figure FDA0003553041940000012
t1＝T1-273.15--------(6)
式(4)中，E1为上层实际水汽压，单位：Pa；RH1为上层空气相对湿度，单位：％；Ea1为上层饱和水汽压，单位：Pa；
式(5)中，t1为上层空气温度，单位：℃；a1和a2为常数，a1＝6.11、a2＝17.27；
式(6)中，T1为上层空气温度，单位：K；
下层水汽压E2通过下层空气相对湿度RH2、下层饱和水汽压Ea2和下层空气温度t2进行换算，换算公式如下：
E2＝Ea2×RH2×0.01-------(7)
Figure FDA0003553041940000021
t2＝T2-273.15-----------(9)
式(7)中，E2为下层实际水汽压，单位：Pa；RH2为下层空气相对湿度，单位：％；Ea2为下层饱和水汽压，单位：Pa；
式(8)中，t2为下层空气温度，单位：℃；a1和a2为常数，a1＝6.11、a2＝17.27；
式(9)中，T2为下层空气温度，单位：K；
②、空气密度ρ：
空气密度ρ的计算公式如下：
Figure FDA0003553041940000022
式(10)中，ρ为空气密度，单位：Kg/m3；RH1为上层空气相对湿度，单位：％；Ea1为上层对应的饱和水汽压，单位：Pa；T1为上层空气温度，单位：K，P为大气压力，单位：Pa；A和B为常数，A＝3.48、B＝0.378；
③、空气比热容CP：
空气比热容CP，单位：J/(kg·K)，计算公式如下：
Cp＝c-dT1-------(11)
式(11)中，T1为上层空气温度，单位：K；c和d为常数，c＝14.439、d＝0.9504×10-3；
第二步：计算温度结构参数
Figure FDA0003553041940000031
通过大口径闪烁仪测量得到的大气折射率结构参数
Figure FDA0003553041940000032
及相关的上层空气温度T1、大气压力P
以及波文比β，计算温度结构参数
Figure FDA0003553041940000033
单位：K2m-2/3，计算公式如下：
Figure FDA0003553041940000034
式(12)中，
Figure FDA0003553041940000035
为温度结构参数，单位：K2m-2/3，β为波文比；T1为上层空气温度，单位：K；P为大气压力，单位：Pa；p0为常数，p0＝-0.78×10-6；
第三步：计算平均感热通量H
平均感热通量H计算公式如下：
Figure FDA0003553041940000036
式(13)中，H为平均感热通量，单位：W/m2；b是经验常数，b＝0.48；ρ为空气密度，单位：Kg/m3；CP为空气比热容，单位：J/kg/K；ZLAS为设备安装高度，单位：m；d为大口径闪烁仪安装位置相对地面的零值位移，单位：m；g为重力加速度，单位：m/s2；T1为上层空气温度，单位：K；
第四步：计算平均潜热通量LvE
平均潜热通量LvE计算公式如下：
LvE＝Q*-H+Gs------(14)
式(14)中，LvE为平均潜热通量，单位：W/m2；Q*为净辐射，单位：W/m2；H为平均感热通量，单位：W/m2；Gs为土壤热通量，单位：W/m2；
第五步：数据质量控制
对计算出的平均潜热通量进行数据剔除和平滑、去噪；
第六步：计算区域蒸散总水量ET
区域蒸散总水量ET的计算公式如下：
Figure FDA0003553041940000041
式(15)中，ET为m分钟内区域蒸散总水量，单位：kg/m2；Lv为蒸散潜热，即蒸散1KG水所需的能量，单位：J/kg；
至此由ET得出整片区域的蒸散总水量。
CN202110271725.4A 2021-03-12 2021-03-12 一种基于闪烁法测量区域蒸散总水量的方法 Active CN113156542B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110271725.4A CN113156542B (zh) 2021-03-12 2021-03-12 一种基于闪烁法测量区域蒸散总水量的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110271725.4A CN113156542B (zh) 2021-03-12 2021-03-12 一种基于闪烁法测量区域蒸散总水量的方法

Publications (2)

Publication Number Publication Date
CN113156542A CN113156542A (zh) 2021-07-23
CN113156542B true CN113156542B (zh) 2022-05-24

Family

ID=76887225

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110271725.4A Active CN113156542B (zh) 2021-03-12 2021-03-12 一种基于闪烁法测量区域蒸散总水量的方法

Country Status (1)

Country Link
CN (1) CN113156542B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116843193B (zh) * 2023-06-14 2024-02-20 华南理工大学 一种基于通量观测数据的人为热通量估算方法及系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104406715A (zh) * 2014-12-15 2015-03-11 重庆市勘测院 一种遥感估算地表感热/潜热通量的精度评价方法及系统
CN111366195A (zh) * 2020-03-27 2020-07-03 北京师范大学 一种地表水热通量的多尺度观测方法
CN111999251A (zh) * 2020-08-14 2020-11-27 中国水利水电科学研究院 基于热红外遥感的区域植被蒸腾和土壤蒸发反演的遥感模型方法
CN112287296A (zh) * 2020-10-13 2021-01-29 北京师范大学 一种基于双波段闪烁仪的地表水热通量测算方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8865305B2 (en) * 2010-06-16 2014-10-21 General Electric Company Core shell phosphor and method of making the same
CN111368258B (zh) * 2020-03-04 2023-02-17 中国科学院东北地理与农业生态研究所 一种湿润地区日蒸散量的估算方法
CN112362693B (zh) * 2020-10-13 2021-12-31 华中科技大学 一种基于土壤热通量推求蒸散发量的测算方法和系统

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104406715A (zh) * 2014-12-15 2015-03-11 重庆市勘测院 一种遥感估算地表感热/潜热通量的精度评价方法及系统
CN111366195A (zh) * 2020-03-27 2020-07-03 北京师范大学 一种地表水热通量的多尺度观测方法
CN111999251A (zh) * 2020-08-14 2020-11-27 中国水利水电科学研究院 基于热红外遥感的区域植被蒸腾和土壤蒸发反演的遥感模型方法
CN112287296A (zh) * 2020-10-13 2021-01-29 北京师范大学 一种基于双波段闪烁仪的地表水热通量测算方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于遥感技术的两种区域蒸散估算方法的应用比较;李星敏 等;《高原气象》;20110215;第30卷(第1期);第125-130页 *

Also Published As

Publication number Publication date
CN113156542A (zh) 2021-07-23

Similar Documents

Publication Publication Date Title
CN111273318B (zh) 一种基于抛物线的区域对流层湿延迟计算方法
Chen et al. Seasonal global water mass budget and mean sea level variations
Rohwer Evaporation from free water surfaces
Bengtsson et al. The need for a dynamical climate reanalysis
CN111368258B (zh) 一种湿润地区日蒸散量的估算方法
CN108920414B (zh) 一种利用气象资料计算局部天顶对流层湿延迟的新方法
Zeng et al. Surface skin temperature and the interplay between sensible and ground heat fluxes over arid regions
CN113156542B (zh) 一种基于闪烁法测量区域蒸散总水量的方法
CN108829975B (zh) 地表温度日内变化过程的遥感估算方法和装置
CN112362693B (zh) 一种基于土壤热通量推求蒸散发量的测算方法和系统
Lee Relative comparison of the local recalibration of the temperature-based evapotranspiration equation for the Korea Peninsula
Steinbrecht et al. Pressure and temperature differences between Vaisala RS80 and RS92 radiosonde systems
Nygård et al. Validation of eight atmospheric reanalyses in the Antarctic Peninsula region
Kochendorfer et al. Undercatch adjustments for tipping-bucket gauge measurements of solid precipitation
CN106814373A (zh) 大气加权平均温度估算及对流层延迟积分方法
Chakraborty et al. A comparative study of cloud liquid water content from radiosonde data at a tropical location
Cheng et al. An improved method for correction of air temperature measured using different radiation shields
CN102521487B (zh) 一种提高气象环境数值模拟准确率的方法
Karam et al. Climatological basin-scale Amazonian evapotranspiration estimated through a water budget analysis
CN112287296B (zh) 一种基于双波段闪烁仪的地表水热通量测算方法
Timilsina et al. Rising temperature trends across the Narayani river Basin in Central Nepal projected by CMIP6 models
Mishra et al. The impact of the time step on the intensity of ITCZ in an aquaplanet GCM
Wang et al. The statistical significance test of regional climate change caused by land use and land cover variation in West China
Keernik et al. Evaluating reanalyses performance in the Baltic Sea region by using assimilated radiosonde data
CN115184967B (zh) 一种扫描微波辐射计水汽数据的gnss校正方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant