VIC水文模型入门攻略（中）VIC的数据和输入文件准备.pdfVIP

VIC⽔⽂模型⼊门攻略 （中）VIC的数据和输⼊⽂件准备 更正 ：上篇末尾关于运⾏汇 程序的部分，在 “STEHE.stnloc”⽂件中需要修改的部分应改为NONE⽽⾮none。希望没有太多⼈因此造成 困扰… VIC的输⼊⽂件都有啥 其实VIC及其汇 程序真正的交互⽅式都是通过读⼊⽂件完成的。它们都要求有⼀个⽂件存储各种设定参数，包括运⾏的⽇期时间、运⾏模 式以及所需的各种数据⽂件的位置。 上⼀篇中启动VIC的命令 ./vicNl -g /home/ .../global_ param.ST EHE 实际上就是将VIC的全局参数⽂件的路径传递给VIC，然后VIC读⼊这些全局 参数并执⾏。 VIC所需的必备输⼊⽂件有 ： 全局参数⽂件 也就是stehekin⽰例中的 global_ param.ST EHE 。 ⼟壤参数⽂件 最核⼼的参数⽂件，保存每个⽹格的⼀些基本参数和各种⼟壤参数，包括⽹格编号，经纬坐标等，具体参考⽰例中 的 soil_ param.ST EHE 。VIC就是从这⾥获得最基本的地理信息的。 植被参数⽂件 存储每个⽹格的⼟地覆盖情况，⽰例如 veg_ param.ST EHE 。 植被库⽂件 存储每种⼟地覆盖类型的⼀些参数，⽐如LAI，结构阻抗，零平⾯位移等。⽰例中为 veglib.LDAS 拿不到相关数据的情况下可 以直接拿这个⽂件⽤。 （难道就是从LDAS那⾥拿的？） ⽓象驱动数据⽂件 即meteorological forcing files （许多资料翻译作⽓象强迫、⽓象迫使，读之蛋疼，故我中⼆翻译为此）。存储每 个⽹格的历时⽓象数据，⼀个⽹格⼀个⽂件，格式为 forci g_纬度_ 经度 。 与上边全是⽂本⽂件形式的⽂件不同，这个可以为⽂本⽂件和⼆进制⽂件两种不同的类型。 其他⼀些可选的数据输⼊⽂件有⾼程带⽂件、雪线⽂件、湖泊/湿地⽂件以及⼀些基本数据的时间序列⽂件等。具体看。 VIC的数据准备 要⽣成以上相应的⽂件，你需要先弄到以下的数据 ： 研究区域内的多个⽓象站的观测数据。 包括每个⽓象站的坐标以及各 ⾃的历时观测数据，精确到每 ⽇，⾄少包含降⾬量，⽇平均⽓温或者 ⽇最⾼/最低⽓温以及风速。通常这 些数据记录都不是完全的，有的数据项因为各种原因会在某些 ⽇⼦缺数据，甚⾄⼲脆却掉某个连续时间段的记录。于是就有必要⽤上数 据处理⼯具。 率定⽤的数据。 ⼀些已经实际测得的VIC可输出的数据，⽐如某地的⼟壤含⽔量数据或者 域出⼝的 量数据。 研究区域的地理信息数据。 ⾸先研究区域范围的经纬度坐标表⽰以及DEM数据是必须的，最好还要有 域边界的shape⽂件，并且这些数据都能在某平台中打开和 处理，⽐如ArcGIS。若实在没有 域边界的图形可以⼲脆直接开辟⼀块矩形区域来研究或者⽤⾼精度DEM提取盆域好了。 DEM可以⽤SRTM计划的，对于VIC的尺度90m的分辨率应该够⽤了。。 有必要提醒下，VIC和汇 程序要求的坐标都是经纬度坐标，位置、⽹格尺⼨什么的都要⽤经纬度表⽰和计算，因此在ArcGIS⾥⾯处理 数据的时候数据框得设置为地理坐标系⽽⾮投影坐标系，获得的地理数据也要进⾏相应变换。 ⼟壤数据。 来源有很多，若⼿头没有的话可以使⽤HWSD，1km分辨率的⼟壤类型数据库，有全球和全国两个范围，可以从IIASA⽹站或者寒区旱 区国家中⼼⽹站下载，CSDN站内也有资源。 不过直接从其得到的资料只有⼟壤类型的分类，砾、砂和粘粒的⽐例和⼀些理化性质等，所以还需要通过⼟壤转换函数换算为VIC需要 的⽔⼒参数。若⼿头没有资源可以参照SPAW这个软件压缩包 （本站也有资源）内附带的那篇论⽂⾥边给出的公式，或者在已发的VIC 相关论⽂⾥找 （记得有些论⽂带了⼟壤类型对应⽔⼒参数表） ⼟地覆盖数据 （植被数据）。 ⼿头没数据可以⽤马⾥兰⼤学发布的全球 1km分辨率⼟壤覆盖类型分布图，。点⾥边的Product Search - AVHRR, Global Land Cover Product便是。 各种⼟地覆盖类型的参数数据。也就是各种⼟地覆盖类型的⼀些动⼒学参数以及不同⽉份LAI之类。其实这个通常不⽤特意准备，可以 直接使⽤stehekin这个⽰例 ⾃带的植被库⽂件veglib.LDAS，其包含了常绿针叶林、常绿阔叶林等 11种