xiaolixiaoli123的个人博客分享 http://blog.sciencenet.cn/u/xiaolixiaoli123

博文

Python与Noah-MP陆面过程模型融合技术及在站点、区域模拟

已有 1051 次阅读 2023-5-29 14:37 |个人分类:生态|系统分类:科普集锦

模型简介:

陆面模型 Noah-MP(TheCommunity Noah Land Surface Model with Multi -Parameterization Options)

是最新合进WRF 的陆面过程方案。具有多种可选物理参数化方案的优势。现有的研究大多数是基于不同陆面方案 WRF 合和不同区域进行模拟分析。

陆面过程是气候系统的重要组成部分, 影响大气环流和气候变化。陆面过程模式中人类活动、 生物物理和生物化学过程的合理描述有助于理解陆面与大气之间相互作用机制

重点

了解陆表过程的主要研究内容以及陆面模型在生态水文研究中的地位和作用;理解Noah-MP模型的原理,掌握Noah-MP模型所需的系统环境与编译环境的搭建方法,熟悉linux系统操作环境及Python基本知识;掌握利用Python语言进行站点与区域数据的处理、单站和区域的模拟、模拟结果的输出和后续分析及可视化等方法。

模型环境搭建、基于Python语言的站点和区域数据处理、站点模拟、区域模拟。区域模拟将针对全球、北美、中国等不同地区的特有数据集进行实战操作。

模型案例

陆面过程的主要研究内容(陆表能量平衡、水循环、碳循环等),介绍陆面过程研究的重要性。

1685325851179.jpg

陆面过程模型的发展历程、基本原理、常用陆面过程模型等。

1685325863697.jpg

Noah-MP模型

Noah-MP模型的发展历史、模型结构及主要模块。

模型运行环境配置、下载与安装

模型以及配套软件的下载

模型需要在Linux下运行,需提前预装Linux系统(推荐使用CentOS系统,下载地址为:

https://www.centos.org/download/

,安装教程可参考:

https://www.runoob.com/w3cnote/vmware-install-centos7.html

运行模型需要提前确定模型运行环境,如系统使用的fortran及C编译器类型等,为之后运行解压缩包下的./configure及Makefile做准备。

模型下载地址:https://github.com/NCAR/hrldas/tree/master/hrldas

1685326061784.jpg

模型以及配套软件的安装

从虚拟机出发,逐步讲解Noah-MP模型运行所需的linux环境的搭建、intel编译器的安装与配置、必要软件的下载与安装、模型构建与编译等内容,细化步骤、逐指令讲解。

上机实操

完成相关软件和linux系统的配置,逐行运行指令,直至hrldas模型搭建编译完成。熟悉linux系统环境,掌握终端(指令行)下进行文件操作的技能,为后续运行模型打下基础

Python编程环境配置

陆面过程模型的数据处理及可视化需要Python语言的协助,在进行陆面模型模拟前需要配置好对应的Python编译环境并掌握Python的基本语法及操作。

环境配置

主要讲解几种不同的环境配置,包括编译器和解释器及包管理器的选择与使用。同时还针对Python的基本语法知识进行讲解。

1685326841953.jpg

处理数据

格式化处理陆面模型运行能够识别的格式

1685327322865.jpg

Noah-MP模型单站运行

大气驱动数据的准备

驱动数据主要包括站点的风速、气温、相对湿度、气压、长波辐射、短波辐射以及降水数据。对于Noah-MP模型而言,原始驱动数据需制作成模型可识别的标准格式,才能够进行下一步的驱动数据编译,将编译结果带入模型进行运算(python脚本)。

1685327974219.jpg

数据时间格式转换

完成驱动数据在世界时和当地时之间的转换。

上机实操

根据示例数据中的站点原始数据,基于python脚本,进行数据的提取合并以及时间转换,生成站点模拟所需的.dat格式文件。

准备静态数据

完成驱动数据的制作后,还需在生成的.dat文件中添加静态数据。此部分数据主要包括站点属性以及模型信息,如站点的海拔,经纬度,土壤类型,初始状态变量的设定和转换系数等

1685328380094.jpg

1685328547940.jpg

驱动数据的编译运行

基于以上数据,生成指定时间步长的一系列.LDASIN_DOMAIN1文件,同时生成hrldas_setup_file.nc文件。

运行模型

根据研究区实际情况与模拟需求修改namelist.hrldas文件,./hrldas.exe即可启动模型,结果将以netCDF格式输出至指定文件夹内。

1685328580566.jpg

模型运行结果的可视化与分析

讲解可视化部分使用的python脚本结构和用法。

借助python netCDF(https://github.com/Unidata/netcdf4-python)或xarray(http://xarray.pydata.org/en/stable)等工具对模拟结果(netCDF格式)进行变量提取与可视化,以用于进一步分析。

1685328682204.jpg

上机实操

基于示例代码中数据可视化部分的python脚本,对叶面积指数、感热通量以及潜热通量等模拟结果进行提取可视化,在熟悉代码结构的基础上,也可对其他变量进行筛选和提取。

单站模拟作业1

基于完整的单站模拟流程,选择课程示例站点之外的一个站点,完成数据下载、变量提取、格式转换、数据编译、模型参数设定、模型运行、结果提取与导出、结果可视化等操作。巩固学习成果,对疏漏点和易错点进行交流,老师全程答疑解惑。

单站模拟作业2

选择不同的模型参数化方案,分别运行模型并对结果进行提取与可视化,对比与其他参数化方案所得结果的差异。

Noah-MP模型区域运行

Noah-MP模型区域运行整体流程与站点模拟相似,大致也可分为数据下载、数据变量提取、制作驱动数据、设置参数化方案、模型模拟与模拟结果可视化分析等几个步骤。但区域模拟的数据准备过程相较于站点运行更为繁琐,使用的数据源也更为多样,如GLDAS、NLDAS以及CLDAS等。同时模拟结果以二维形式存储,提取与可视化的方法也不同于单站。

区域运行将会讲解基于GLDAS和NLDAS等数据源制作驱动数据的模拟方法,并介绍区域运行结果可视化工具xarray的使用方法。课程结构主要分为数据下载、数据处理、数据编译、模型运行与结果可视化几个部分,下面以GLDAS为例,详细介绍讲解流程与配套练习内容:

准备大气驱动数据

Noah-MP模型的运行需要格式正确的驱动数据(气温、降水、气压、风速、辐射等)。因此用户需提前下载并处理好相关数据,以下为数据准备的简要步骤。

下载大气驱动数据-以GLDAS为例

区域驱动数据主要包括研究区的风速、气温、相对湿度、气压、长波辐射、短波辐射以及降水等。可从GLDAS官网下载相关数据:

https://disc.gsfc.nasa.gov/datasets/GLDAS_NOAH025_3H_2.1/summary?keywords=GLDAS

1685328709648.jpg

上机实操

讲解两种区域数据的获取方法,一种是基于downthemall的批量下载方法,操作简便;另一种是基于python脚本的数据抓取方式,自由化更高;大家可尝试不同的数据下载方法,获取区域驱动数据的源数据。

相关变量的提取与时间转换

对模型模拟所需变量进行提取,同时完成世界时和当地时之间的转换。

1685340823846.jpg

初始状态变量的提取与时间转换

1685340859629.jpg

风速的分解、降水数据的整合

1685341637556.jpg

1685341656452.jpg

上机实操

基于下载的netCDF4格式的源数据,分别编辑并运行以上python及perl脚本,生成变量分解后的一系列netCDF4格式文件,用于编译生成驱动数据。

准备研究区静态数据

制作geo_em_d0x.nc数据

基于WPS制作区域静态数据(包括研究区范围、研究区土地利用情况、植被覆盖度等信息)。此部分需安装WRF及WPS,并下载WPS_GEOG数据。详细流程可参考:

https://www2.mmm.ucar.edu/wrf/OnLineTutorial/compilation_tutorial.php

1685341690277.jpg

准备数字高程数据

hrldas (Noah-MP)自带的高程文件-GLDASGLDASp4_elevation_025d.nc4。一般情况下,将此文件作为模型的高程输入数据即可。 

编译与运行-生成模型驱动数据

区域原始的气象驱动数据和静态数据准备完毕后,需编译运行生成符合模型要求的驱动数据-*.LDASIN_DOMAIN1。

1685341709807.jpg

1685341739542.jpg

上机实操

修改好相关目录的目录结构,在指定目录下,准备好必须的程序与文件,编译运行生成.LDASIN_DOMAIN1格式的驱动数据。

运行hrldas (Noah-MP)模型

根据研究区特点及用户模拟需求,修改namelist.hrldas文件,Namelist.hrldas文件编辑完成后,即可运行可执行程序/hrldas.exe。模型模拟结果将输出至指定文件夹下,若运行成功,在输出目录下应包含指定模拟时间段内的指定时间步长的模拟结果-*.LDASOUT_DOMAIN1。

1685341859041.jpg

上机实操

练习设置namelist.hrldas,选择一套参数化方案,运行模型。

数据的分析与可视化

区域模拟结果亦为netCDF格式文件,借助xarray,pandas等工具进一步进行变量提取、可视化等工作。

1685341889249.jpg

上机实操

基于提供的python脚本,提取模拟结果中的相关变量并可视化,熟悉xarray的使用方法。

区域模拟作业1

基于完整的区域模拟流程,选择课程示例区域之外的一个区域,完成数据下载、变量提取、格式转换、数据编译、模型参数设定、模型运行、结果提取与导出、结果可视化等操作。巩固学习成果,对疏漏点和易错点进行交流,老师全程答疑解惑。

区域模拟作业2

选择不同的模型参数化方案,分别运行模型并对结果进行提取与可视化,对比与其他参数化方案所得结果的差异。




https://blog.sciencenet.cn/blog-3559456-1389809.html

上一篇:SolVES 模型与多技术融合【QGIS、SQL、ARCGIS、MAXENT、R】如何实现生态系统服务功能社会价值评估
下一篇:生态系统NPP及碳源、碳汇模拟(土地利用变化、未来气候变化、空间动态模拟)【CASA模型】
收藏 IP: 111.225.71.*| 热度|

1 郑永军

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-4-27 02:19

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部