CN110110410A - 一种基于wrf模式与sebs模型耦合的蒸散量估算方法 - Google Patents
一种基于wrf模式与sebs模型耦合的蒸散量估算方法 Download PDFInfo
- Publication number
- CN110110410A CN110110410A CN201910337712.5A CN201910337712A CN110110410A CN 110110410 A CN110110410 A CN 110110410A CN 201910337712 A CN201910337712 A CN 201910337712A CN 110110410 A CN110110410 A CN 110110410A
- Authority
- CN
- China
- Prior art keywords
- meteorological
- evapotranspiration
- wrf
- slope
- sebs
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Abstract
本发明公开了一种基于WRF模式与SEBS模型耦合的蒸散量估算方法,包括如下步骤:步骤1)利用最新的土地覆盖数据更新WRF模式自带的下垫面信息;步骤2)利用WRF模式模拟区域气候参量,获取高分辨率气象场,提取风速、气温、相对湿度、气压数据形成气象数据集;步骤3)基于DEM提取坡向(ASPECT)、坡度(SLOPE)信息,并将DEM、ASPECT、SLOPE三个波段按顺序叠加,形成地形数据集;步骤4)利用MODIS反射率产品、地表比辐射率产品及地表温度产品,按顺序叠加形成基础数据集;步骤5)利用步骤2)、3)、4)得到的气象数据集、地形数据集和基础数据集输入SEBS模型,计算得到蒸散量。
Description
技术领域
本发明是一项陆表信息反演领域的技术,针对当前蒸散量反演模型(SEBS)输入的气象场数据分辨率较低导致蒸散量反演效果不佳的问题,本发明利用WRF模式获取高时空分辨率气象场数据作为输入,能够明显地改善基于SEBS模型的区域蒸散量估算效果,是一项具有应用价值的空间分析技术。
背景技术
蒸散量是水循环的一个相当重要的参量,用于表征陆地表面水分向大气的输送,由土壤、水面蒸发和植被蒸腾所组成,在调节局部或区域气候方面具有重要作用。就全球尺度而言,蒸散量占全球总降水的一半以上,意味着由大气输送到地表的水分将有大部分不会被存储下来,因此对于蒸散量的精确估算及时空分布的研究将对区域水资源管理和宏观调控政策的制定显得意义非凡。
传统的蒸散量估算多是基于气象观测站点的单点计算,虽然能提供相对准确的蒸散量,但由于大尺度上,下垫面和水热传输具有非均质性,使得单点计算蒸散量无法满足区域研究的需求。遥感技术因其对区域尺度信息捕获的高效准确,在非均质下垫面的区域蒸散量估算上表现出了明显的优越性。对蒸散量的定量估算一直以来都是水循环科学领域的热点话题,学者们也一直致力于算法模型的改进以求能够获得高分辨率和高精度可以兼得的更优模型,但大气环境异常复杂,与陆地的相互作用关系更是需要极多的物理参数来表达,而其中的一些参数并不容易确定。
截至目前,国际上已经形成了诸多基于地表能量平衡的蒸散量算法模型,应用最为普遍的是基于地表能量平衡理论的SEBS(Surface Energy Balance System)模型。在以往为SEBS模型提供气象数据的案例中,学者们通常利用气象站点数据进行数学插值,以得到区域尺度的网格气象数据,也有利用粗分辨率的同化分析资料通过数学插值为SEBS模型提供气象输入的研究。虽然取得了很多成果,但对这些模型的应用都有一个共同的问题,即由于气象站点往往分布稀疏,同化资料的分辨率又太粗,且数据时间间隔一般较大通常为几个小时间隔,纯粹的数学插值很难确保在得到高分辨率气象场的同时也能具有可靠的精确度,因此不能为模型提供较为精确的网格气象驱动数据,而导致计算出的蒸散量高分辨率和高精度往往不可兼得。
气象模式WRF(weather research and forecasting)能够按照固定时间间隔输出模拟数据,时间分辨率可达到小时级,空间分辨率可达到公里至百米级,其理论基础是在大气物理环境中应用动力学、热力学理论基于数学有限差分思想在充分考虑非均质下垫面的前提下模拟获得连续的气象场,显然要比不考虑下垫面的非均质性而进行简单的数学插值更精确,从而能为SEBS模型提供更为精确的气象数据输入。基于此,本发明尝试利用WRF模式模拟得到高时空分辨率的栅格气象数据代替常规地面站点数据或粗分辨率同化资料的简单数学插值,为SEBS模型提供相对可靠的气象数据输入,得到分辨率较高且精度更为可靠的蒸散量估算结果。
发明内容
本发明提出一种基于WRF模式与SEBS模型耦合的蒸散量估算方法,通过将WRF模式与SEBS模型进行耦合以改进SEBS模型所需气象信息,不再是传统的气象站点或粗分辨率的同化资料进行简单的数学插值,WRF模式能够按照固定时间间隔输出模拟数据,大大提高了气象数据的时间分辨率,其理论基础是在大气物理环境中应用动力学、热力学理论基于数学有限差分的思想获得连续的气象场,在提高空间分辨率的同时也考虑了下垫面的非均质性。总体流程如图1,该方法包括如下步骤:
步骤1)利用最新的土地覆盖数据更新WRF模式自带的下垫面信息;步骤2)利用WRF模式模拟区域气候参量,设置嵌套,输出每小时间隔的高时间分辨率和1Km高空间分辨率的气象场,提取风速、气温、相对湿度、气压数据形成SEBS所需气象数据集;步骤3)基于DEM(Digital Elevation Model)生成坡向(ASPECT)、坡度(SLOPE)信息,并将DEM、ASPECT、SLOPE三个波段按顺序叠加,形成SEBS所需地形数据集;步骤4)获取MODIS(ModerateResolution Imaging Spectroradiometer)反射率产品、地表比辐射率产品及地表温度产品,按顺序叠加形成SEBS所需基础数据集;步骤5)利用步骤2)、3)、4)所得数据集输入SEBS模型,计算得到蒸散量。
进一步地,所述步骤2)中设置嵌套是为了逐步提升气象场的空间分辨率,范围从大到小,逐步聚焦研究区,分辨率从粗到细。
附图说明
图1为基于WRF模式与SEBS模型耦合的蒸散量估算方法流程图;
图2为WRF模式嵌套网格设置示意图;
图3为利用WRF模式与SEBS模型耦合估算的蒸散量(ETa-WRF);
图4为利用气象站点数据作为输入估算的蒸散量(ETa-STA)。
具体实施方式
下面结合附图对本发明“一种基于WRF模式与SEBS模型耦合的蒸散量估算方法”作进一步说明。
本发明试图将气象模式WRF与SEBS模型建立耦合,通过在大气物理环境中基于数理理论以粗分辨率的同化分析资料为驱动,WRF模式输出高分辨率的气象场,替代传统点源或粗分辨率气象数据的直接数学插值,为SEBS模型提供相对可靠的气象输入,以获得更为可靠的蒸散量分布结果,为区域水资源管理和宏观调控政策的制定提供参考。
1)更新WRF模式下垫面
利用2015年的CCI-LC(Climate Change Initiative Land Cover)土地覆盖数据更新WRF模式下垫面,此套数据分辨率300m,全球地表覆被类型共23大类,将其按照WRF模式自带的全球24类覆被类型进行重编码,输出ASCII码后经过预处理成为WRF模式可识别的数据。
2)制作气象数据集
对研究区设置三重嵌套,即D01、D02、D03(图2),其中D01的中心经纬度为37.735°N、82.914°E,三重嵌套的网格数分别为61×51、116×51、141×111,对应网格距分别为25Km、5Km、1Km。初始及边界条件由1°×1°NCEP(National Centers for EnvironmentalPrediction)的全球分析资料提供。模拟时间为2015年7月20日-27日,设置每小时输出一次结果,取D03模拟结果进行后续处理,输出8天平均的风速、气温、相对湿度及气压,形成SEBS模型所需气象数据集。
3)制作地形数据集
基于ASTER GDEMV2 30m分辨率的DEM数据提取坡向(ASPECT)、坡度(SLOPE),并将DEM、ASPECT、SLOPE重采样至1Km,按顺序进行波段叠加后形成地形数据集。
4)制作基础数据集
下载与气象数据时间相应的MODIS 8天合成反射率产品MOD09A1(500m)、地表温度及地表比辐射率产品MOD11A2(1Km),将7个反射率波段、地表比辐射率、地表温度9个波段按顺序叠加,重采样至1Km,形成SEBS模型所需基础数据集。
5)计算蒸散量
将准备好的地形数据集、基础数据集、陆表物理参数信息及气象数据集输入SEBS模型,运行得到2015年7月20日-7月27日的8天平均蒸散量(记为ETa-WRF)(图3)。
为了比较估算效果,利用气象站点数据直接插值得到气象数据集作为输入,采用SEBS模型估算相同区域、相同时间的蒸散量(记为ETa-STA)(图4)。经过对比发现,ETa-STA存在大量的蒸散量估算错误,尤其在南部裸地区域,与实际地表类型对应的蒸散量不符,而本发明得到的ETa-WRF噪声明显比ETa-STA少,而且对比度、轮廓都更清晰。
Claims (3)
1.一种基于WRF模式与SEBS模型耦合的蒸散量估算方法,该方法包括如下步骤:步骤1)利用最新的土地覆盖数据更新WRF模式自带的下垫面信息;步骤2)利用气象模式WRF模拟区域气候参量,获取高分辨率气象场,提取风速、气温、相对湿度、气压数据形成气象数据集;步骤3)基于DEM提取坡向(ASPECT)、坡度(SLOPE)信息,并将DEM、ASPECT、SLOPE三个波段按顺序叠加,形成地形数据集;步骤4)利用MODIS反射率产品、地表比辐射率产品及地表温度产品,按顺序叠加形成基础数据集;步骤5)利用步骤2)、3)、4)得到的气象数据集、地形数据集和基础数据集输入SEBS模型,计算得到蒸散量。
2.如权利要求1所述的估算方法,其特征在于,所述步骤2)是通过设置3层嵌套来逐步提升气象场的空间分辨率。
3.如权利要求1所述的估算方法,其特征在于,所述步骤5)是利用气象模式WRF模拟的气象场充分考虑了下垫面的非均质性,相比利用气象站点数据或者粗分辨率的同化分析资料通过数学插值形成的气象场精度高,将该气象场作为SEBS模型的输入,从而提高蒸散量的估算精度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910337712.5A CN110110410A (zh) | 2019-04-25 | 2019-04-25 | 一种基于wrf模式与sebs模型耦合的蒸散量估算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910337712.5A CN110110410A (zh) | 2019-04-25 | 2019-04-25 | 一种基于wrf模式与sebs模型耦合的蒸散量估算方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110110410A true CN110110410A (zh) | 2019-08-09 |
Family
ID=67486676
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910337712.5A Pending CN110110410A (zh) | 2019-04-25 | 2019-04-25 | 一种基于wrf模式与sebs模型耦合的蒸散量估算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110110410A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110674571A (zh) * | 2019-09-18 | 2020-01-10 | 国网湖南省电力有限公司 | 一种输电线路垭口风速降尺度计算方法及其系统 |
CN111523381A (zh) * | 2020-03-13 | 2020-08-11 | 上海眼控科技股份有限公司 | 数值天气预报中土地利用信息的更新方法及设备 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102034003A (zh) * | 2010-12-16 | 2011-04-27 | 南京大学 | 基于蓄水容量曲线和topmodel的流域水文模型的设计方法 |
CN103364532A (zh) * | 2012-03-28 | 2013-10-23 | 中国环境科学研究院 | 一种定量评估区域水源涵养量的实用方法 |
US20180315100A1 (en) * | 2017-04-27 | 2018-11-01 | Lindsay Corporation | Computer-implemented method and computer program for designing and cost-estimating irrigation systems |
CN108776745A (zh) * | 2018-06-27 | 2018-11-09 | 广东电网有限责任公司电力科学研究院 | 一种基于WRF和Fluent耦合的复杂地形风场模拟方法及装置 |
-
2019
- 2019-04-25 CN CN201910337712.5A patent/CN110110410A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102034003A (zh) * | 2010-12-16 | 2011-04-27 | 南京大学 | 基于蓄水容量曲线和topmodel的流域水文模型的设计方法 |
CN103364532A (zh) * | 2012-03-28 | 2013-10-23 | 中国环境科学研究院 | 一种定量评估区域水源涵养量的实用方法 |
US20180315100A1 (en) * | 2017-04-27 | 2018-11-01 | Lindsay Corporation | Computer-implemented method and computer program for designing and cost-estimating irrigation systems |
CN108776745A (zh) * | 2018-06-27 | 2018-11-09 | 广东电网有限责任公司电力科学研究院 | 一种基于WRF和Fluent耦合的复杂地形风场模拟方法及装置 |
Non-Patent Citations (3)
Title |
---|
吴雪娇 等: "基于涡动相关仪验证的SEBS模型对黑河中游地表蒸散发的估算研究", 《冰川冻土》 * |
王凯霖: "柴达木盆地平原区土壤湿度的遥感反演及对蒸散发的影响", 《中国优秀硕士学位论文全文数据库 农业科技辑》 * |
赵德刚 等: "面向土地覆盖分类的MODIS 影像融合研究", 《遥感技术与应用》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110674571A (zh) * | 2019-09-18 | 2020-01-10 | 国网湖南省电力有限公司 | 一种输电线路垭口风速降尺度计算方法及其系统 |
CN110674571B (zh) * | 2019-09-18 | 2022-06-24 | 国网湖南省电力有限公司 | 一种输电线路垭口风速降尺度计算方法及其系统 |
CN111523381A (zh) * | 2020-03-13 | 2020-08-11 | 上海眼控科技股份有限公司 | 数值天气预报中土地利用信息的更新方法及设备 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105912836B (zh) | 一种纯遥感数据驱动的流域水循环模拟方法 | |
Lutz et al. | Comparison of climate change signals in CMIP3 and CMIP5 multi-model ensembles and implications for Central Asian glaciers | |
Mora et al. | Climate changes of hydrometeorological and hydrological extremes in the Paute basin, Ecuadorean Andes | |
CN106446564A (zh) | 一种植被净初级生产力遥感估算方法 | |
CN110276304A (zh) | 基于降尺度的高分辨率植被生产力遥感估算方法 | |
CN101614651A (zh) | 一种土壤水分监测的数据同化方法 | |
Yang et al. | Biases and improvements in three dynamical downscaling climate simulations over China | |
CN102162850A (zh) | 基于模型的小麦产量遥感监测预报方法 | |
Zhang et al. | Dynamical downscaling of the climate for the Hawaiian Islands. Part I: Present day | |
Fan et al. | Medium-range forecasting of daily reference evapotranspiration across China using numerical weather prediction outputs downscaled by extreme gradient boosting | |
CN106021868A (zh) | 一种基于多规则算法的遥感数据降尺度方法 | |
CN106019408A (zh) | 一种基于多源遥感数据的高分辨率卫星遥感估算方法 | |
CN103793511A (zh) | 一种提高风速预报精度的方法 | |
CN105321120A (zh) | 同化蒸散发和叶面积指数的区域土壤水分监测方法 | |
CN104143043A (zh) | 一种多功能气候数据模型及其应用 | |
CN102682335A (zh) | 精确确定区域对流层延迟的神经网络方法 | |
CN110110410A (zh) | 一种基于wrf模式与sebs模型耦合的蒸散量估算方法 | |
CN110599360A (zh) | 一种干旱区农作物生长季蒸散发高分辨率遥感估算方法 | |
Zhai et al. | Climate change in China affects runoff and terrestrial ecosystem water retention more than changes in leaf area index and land use/cover over the period 1982–2015 | |
Ye et al. | Flood forecasting based on TIGGE precipitation ensemble forecast | |
Nasonova et al. | Projecting changes in Russian northern river runoff due to possible climate change during the 21st century: a case study of the Northern Dvina, Taz and Indigirka rivers | |
Zhang et al. | Downscaling hourly air temperature of WRF simulations over complex topography: A case study of Chongli District in Hebei Province, China | |
Shin et al. | An overview of KMA’s operational NWP data assimilation systems | |
CN109657988A (zh) | 基于hasm和欧氏距离算法的烟叶品质分区方法 | |
Dirmeyer et al. | The second Global Soil Wetness Project (GSWP-2): Multi-model analysis and implications for our perception of the land surface |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190809 |