Category Archives: work

地表过程建模项目完成验收

原题 中科院知识创新工程重要方向项目群“地表过程集成系统研究”各项目顺利通过验收

9月9日,中科院知识创新工程重要方向项目群“地表过程集成系统研究”各项目结题验收会议在地理资源所召开。验收组由8位专家组成:国家自然科学基金委宋长青研究员、冷疏影研究员,中国科学院遥感应用研究所施建成研究员,中国科学院寒区旱区环境与工程研究所李新研究员,中国科学院地理科学与资源研究所周成虎研究员,中国科学院水土保持研究所邵明安研究员,中国科学院新疆生态与地理研究所陈亚宁研究员,中国科学院南京土壤研究所张甘霖研究员。项目组织部门中科院资环局副局长冯仁国、国土与遥感处处长翟金良、副处长段晓男出席了会议。

首先,冯仁国在会上致辞,他介绍了项目的立项背景,强调了中国科学院对纯科学问题研究的重视及投入,简述了地表过程的集成研究的研究背景,同时感谢各项目和课题的承担单位及研究人员的辛勤工作,并对做好课题验收及后续项目验收工作提出了要求。

中科院寒区旱区环境与工程研究所南卓铜研究员、中科院遥感应用研究所陈良富研究员、中科院南京地理与湖泊研究所张奇研究员、中科院地理资源所刘慧研究员先后代表四个项目做了结题报告,向专家组汇报了课题的目标任务及完成的主要成果、应用情况、经费支出等情况,并对专家组的疑问与质询进行了逐一解答。

专家组在听取了项目的汇报,查阅相关技术报告等项目验收材料,并进行质询和讨论后,对四个项目分别给出验收意见。专家组一致认为,各项目都完成了项目任务书规定的研究内容,达到了考核指标要求,验收文档齐全,同意通过验收。

出处:北京地理所

叶城阿里区域的冻土分布图

根据科技部基础调查项目的西昆仑、改则两个调查区的野外调查数据,分别用冻土下限—高程关系和扩展的冻结数模型做了叶城—阿里区域的多年冻土分布图,因为扩展的冻结数模型是根据两个调查区进行了参数率定,从结果上看好像更合理一些,因此采用这个图。

image

图件用了250m的MODIS LST,90m STRM DEM数据,出图比例尺是1:170万(实际上可以做到100万,但因为是做到A1幅面上了,最大只能放到170万)。

我估计这应该是这区域目前最精确的冻土分布图了。图件提交给某给水部队作为该区找水打井之用。

打井也不容易

上午银川某部队水文勘测的M总和W工过来找我们。他们要给叶城—阿里区域边防驻地解决用水问题,通俗的讲就是打个水井。但因为在冻土区可能面临的问题他们没有经验,要求我们给提供一些有关冻土分布、厚度等资料以及冻土对地下水资源影响的有关专家经验。谈话里讲到冻土区打井可能遇到的一些问题,感觉很有意思并且深受启发,有机会可以在这方面深入做一些工作。

比如在多年冻土区打水井,一般需要取冻土层下承压层的地下水,那么如果打穿的话,水在压力作用下会上涌,区别于非冻土区,经过一段时间井里的水可能会结冰,那么井就废了。怎么解决这个问题? 在井下增加加温装置可能不大实用,边防驻地没有大功率发电机,而且费用不菲。一个可能的方案是水冻结前,使用水泵再抽一次水,利用地下水本身的正温防止水井冻结。然而这里还有问题,多长时间抽一次呢? 可能需要做一些实验建一些模型进行计算,得到合理的一个时间间隔。

这些实际的工程需求平时都没有注意到。埋首实验室辛苦工作的同时,如果能结合一些实际需求做一些工作也是很有意思的。

参加揭牌仪式

跟亮子、立宗本是因为一个项目工作过去的西宁,被通知12.2参加青海省生态环境遥感监测中心成立的揭牌仪式。

青海省生态环境监测中心揭牌仪式
青海马省长、环保杨厅、北京环保部领导参加的新单位揭牌仪式

swat 2005 源代码的编译

swat 2005源代码是可以从官方网站上下载到的。尝试了gfortran 编译,首先删除了*.mod,使用 gfortran –c modparm.f 重新生成 parm.mod。在main.f里注释了第一行,即include ‘modparm.f’。否则会出现重复定义的问题。最后编译,可以成功,但在运行时出现 Fortran runtime error,大致是从文件里读取end of file之类的错误。

安装了intel fortran 编译器,用的是官方的30天试用版。在visual studio 2010里建立了fortran console application项目,将swat 2005里的全部*.f 加到里面。同样重新生成 parm.mod,并修改了main.f (同上),编译,试运行,可以用。

GLDAS v1 0.25度3小时降水数据

学生下载了gldas v1 0.25x 0.25度的 3-hourly 数据,据介绍里面是有rainfall数据的。但用grads或者matlab read_grib 读取时,发现没有此项,如图:

image
(click for enlarged view)

这很奇怪,因为相关说明里分明是有rainfall一项的,如在http://hydro1.sci.gsfc.nasa.gov/dods/GLDAS_NOAH025SUBP_3H.info介绍的,

image

而且对比两者,里面的参数名也不一致。

仔细研究gldas数据对应的元数据,发现每个文件每项里都有对应的最大值最小值,试着读取并对比,发现第7项4LFTX即是rainf 一项,最大值一致。看样子是read_grib等代码是专门for WMO的,跟gldas的描述不一致。所以在读取gldas等数据时,还需要注意了。

ftp://hydro1.sci.gsfc.nasa.gov/data/s4pa/GLDAS/GLDAS_NOAH025_3H_E1.002/doc/README.GLDAS2.pdf找到parameter table的相关信息,如下图

image

注意其中的pds ids,rainfall rate是132,对应于图1的第7个记录,#=132。

大概是因为read_grib的参数表是基于ncep的某种格式,读取gldas需要配置自定义的参数表。

基金终于赶出来了

本来是要写另一个的,结果三天前换成决策支持系统的。累的呀,不过终于在截止期前也赶出来了。感谢书功、彦博、小余、小曹的大力帮忙,光我一个人那累死也忙不出来的。

这期间还搞了一个三江源生态环境监测的本子,今天315也是截止日期,因为要查新报告和合作协议,今天也是赶出来的,让小余租了个车给西宁那边给送过去了。感谢小余和老杨。

希望过了这阵子可以稍轻松一些了。