作者:陈天舒,山东大学环境科学(大气化学)博士在读
发现网上关于用Hysplit做后推气流轨迹聚类分析的教程较少并且不详细,因此萌生写个速成指南的想法。
后推气流轨迹的聚类指的是将一段时间内的后退气流轨迹按照一定方法进行分类归纳。聚类分析可以用于分析目标地点的气团的来向构成以及占比。
我们常在文献中看到这样的图,用于分析研究地点的气团传输机制:
也会看到这样的表,分析不同来源的气团的污染物浓度:
这些分析通过Hysplit可以做到。
分为四个步骤:气象数据准备;按所需时间段批量计算后推气流轨迹;对后退气流轨迹数据进行聚类;将聚类结果与污染物浓度数据结合分析。
一、下载气象数据
要计算后退气流轨迹,最重要的是气象数据,一般常用的气象数据为GDAS1。
下载地址为:下载gdas1气象数据
文件格式为“gdas1.英文月份缩写+年份.第几周”。例如,“gdas1.apr05.w3”指的是2005年四月第三周的气象数据。
其他Hysplit支持的气象数据格式详见:
READY - Gridded Data Archivesready.arl.noaa.gov二、按所需时间段批量计算后退气流轨迹
分为两小步:
1、初始参数设置
打开Hysplit主界面->“Trajectory"->“Setup Run”
无论是跑单次的后推气流轨迹还是跑一段时间的,都需要在Setup Run的页面里进行参数的预先设