實操：利用Hysplit繪制后向軌跡圖本文主要從如何下載氣象數據、如何聚類分析、如何更換自動出的后向軌跡底圖、如何疊加各方向來源的污染物濃度圖幾個方面進行介紹。類似下圖。一、提前下載安裝離線版hysplit一、提前下載安裝離線版hysplit軟件直接在美國國家海洋和大氣管理局選擇適合自己電腦的版本，網址/documents/Tutorial/html/index.html二、氣象數據的準備（當時氣象數據的下載走了很多彎路，因為不知道怎么下載，下載多少數據夠用。一旦你數據準備沒做好，后面跑模型很容易出問題）①、ftp數據下載網址：/pub/archives/gdas1.v1/（需要注意有的時候這個網址沒法直接進去，則需要掛vpn）②、下載的數據形式如下：gdas1.英文月份縮寫+年份.第幾周gdaslgdasljar19.w1l1gdaslgdasljar19.w1l1gdasljan19m2201班/1412:55W1文件584,853KB2019/9/916:47W1文件584,853KB2019/9/91732W2文件584,853KB③、一般根據你的采樣時間及你想做的后向軌跡后推時間確定下載的文件個數比如，你的采樣時間是2018年7月1日-7月30，然后做的72小時后向軌跡模型分析，

2018年7月1日后推三天的時間是6月28,所以你下載數據的時候，起始文件應該是2018年6月28所在周的文件，即gdas1.jun18.w5，同理可以得到結束文件是gdas1.jul18.w4（為了使得運行更順利，一般可以在起始文件和結束文件前、后再多下載一個周的數據文件），所以我們需要下載的數據文件有：gdas1.jun18.w4，gdas1.jun18.w5，gdas1.jul18.w1，gdas1.jul18.w2，gdas1.jul18.w3，gdas1.jul18.w4，gdas1.jul18.w5三、按所需時間段批量計算后退氣流軌跡①、初始參數設置：打開①、初始參數設置：打開Hysplit主界面->"Trajectory"->"SetupRun”②、參數設置：/'IrajectorySetupStartingtrue(￥YM4DDHH(hlh]):1701QL□□Numbarofstarting"lo-eatiDris:2====>Setupfitarti-ng1looatirmsSaveasNumbarofstarting"lo-eatiDris:2====>Setupfitarti-ng1looatirmsSaveas（1）startingtime（采樣開始時間）：按年月日小時的順序輸入起始時間，注意空格。(2)起始位置的數量：一般填1,同時點擊右側的Setupstartinglocation,填寫坐標經緯度以及起始高度，(經緯度用小數表示，高度用100、500、1000)。本次示例為39.1°N，117.1°E,100/StartingLocation￡...Setup1StartingLocations

LatitudeLongitudeHeight(m-AGL)Location1:Location1:39.1117.110Q(3)總運行時間：按實際需求填寫，一般和你要研究的物種的大氣化學壽命相關，選擇24h或36h或48h或72h等。(4)方向：選擇后推“Back”。(5)模型層高：填寫模型的最大高度，一般按默認即可。(6)垂直運動模式：點擊“Select"，選擇3=constantdensity",按輸入的氣象數據內置的垂直運動模式即可。(7)輸出：軌跡結果輸出目錄。注意一個細節，在目錄最后要加上幾個字母作為這次計算的輸出文件的前綴，以便下一步計算的時候識別并鎖定。新建一個文件夾，文件名稱為test。(8)添加氣象數據：按所要計算的時間范圍選擇多個氣象文件。(9)可點擊“saveas”保存配置文件，以便下一次使用。點擊“save”保存配置。③、運行批量計算(1)返回Hysplit主界面->“SpecialRuns"->"Daily”選擇一天中每隔幾個小時計算一次，直到多少天后截止。這個和你后續的目的相關，如果你想將其與氣態污染物濃度相聯系，比如你的氣態污染物濃度數據是每

個小時的數據，你可以每隔1h計算一次，如果你是晝夜樣品，有晝夜氣態污染物的濃度，那你可以每隔12h運行一次，如果沒有特殊要求，可以定義為6h。然后直到多少天后則根據你的采樣階段所持續的時間,AutomatedMultipleTrajectoriesbyTime ——lXIEzecuteascripttorunmultipleiterationsthetr-ajec.tQrycalculation±ronthestarttijnsatthe.intHoatedinterval.TheCONTROLfi.lB.sli&uldhavepreviouslyconfiguredinthesetupXL-enii.Out-pLitfilename3areappendedwithtlieysax/mantL/day/hour-Startcalculatio-nsat19DI03ODandevery6hruntilRiJdays口h.rafterstartI—It|Help|ExecuteScript（2）點擊“ExecuteScript”，稍等片刻（這時候會不斷閃現對話框），直到下圖出現，表示計算完成。點擊“Exit”StartedStartedStartedStartedStartedStartedStartedStar-tedStartedStar-tedStartedcalaulaticm:oalaulatiori"StartedStartedStartedStartedStartedStartedStartedStar-tedStartedStar-tedStartedcalaulaticm:oalaulatiori"calaulaticm:oalaulation"calculatioii:oalaulatioii:aalculaticm:GmluulaitiiciXL;calaulaticm:umluulmi七ic>:n:calaulaticm:testL7O120OtitestL7012ei2testL7012818testL701290QtestL7Q1290fitestL7012912testL7D12918t?stl7013000testL7D13D0fit?s.tl7013D12testL7D13D18CailGulaitionscample七日CailGulaitionscample七日d!可以在你剛剛新建的文件夾中發現這樣的結果文件，這就是計算的后向軌跡原始文件。testl7010200?test17010212□te5t17010300test170ia312四、對后退氣流軌跡數據進行聚類（聚類只能對一個位置一個高度聚類）①、我們現在有多日多條的后推氣流軌跡，接下來要做的就是將所有軌跡按照距離進行歸類。②、返回Hysplit主界面->“SpecialRuns"->"Clustering"->"Standard"fTrajectoryClusterAnaly&is(Jan.2010) — □XStep1;InputsRvm__ID:Standai2dHourstocluster:36Timeinterval(hrs):1Trajectoryskip: 1IEndpointsfolder:同77天氣溶膜/圖片摘錄/徐文3/后 BrowseWorlkingfolder:07大氣溶峻/悶片摘錄/景文3/后 BrowseArchivefolder:E:，天氣溶胺/附片摘錄/論文3/后 BrowseProjection(for'Runclusterj?C*-AutorPolarCLaMbertCMercatorStep2;RimClusterProgzaniMakeinfileEaseName:ftestNoteIRunclusteranalysisDisplaytotalspatialvarianceViewposaxb-Lefinalnumloer?ofalustars-Criterion(%):C2Df*30Ruristep3:GetResults(repeatfordifferentnrnnJoerofclusters)NumberofClusters: 1NoteAssigntrajectorLestoalusters: pView RunDisplayMeansDisplayclustersDisplayNotUsedarchive Quit Help③、參數設置：1、配置。Hourstocluster聚類時間：需要進行聚類的時間長度。建議做多少小時的后推氣流軌跡，就聚類多少小時。(2)Timeinterval聚類時間間隔：選擇將軌跡劃分為幾段進行聚類計算。例如36小時的聚類時間+1小時的聚類時間間隔，就意味著把所有軌跡每隔1小時進行一次距離比對，共進行36小時的距離比對。(3)Trajectoryskip軌跡間隔：簡單地說就是每隔幾條納入計算范圍。例如填寫1，就是每一條都參加計算，填寫2,就是每兩條取一條進行計算。(4)Endpointsfolder端點文件夾：選取批量計算的軌跡結果所在的文件夾。(即

你新建的名稱包含前綴名test的那個文件夾workingfolder工作文件夾：輸出聚類結果的文件夾。Archivefolder結果文件夾：可以同workingfolder。2、聚類(1)首先要鎖定目標軌跡文件。這時候就需要用到我們在批量輸出軌跡時預先設置的文件前綴。例如上文提到輸出的批量軌跡的前綴是“test”，這里就填寫比蝮。點擊“MAKEINFILE”。如果成功識別鎖定我們要的那批軌跡文件，會有如下提示/Hysplit4SystemMessage — □/Hysplit4SystemMessage — □X⑵點擊“Runclusteranalysis'，會出現如下對話框依次聚類每條軌跡，完成后點擊退出即可。Pelss67outq￡75Pass68outof75Pass69out□￡75Pass70outof75Pass71outof75Pass72outo￡75Pass73outof75Pass74outo￡75Pass75outof75CoupleteCLusten—nuBLberoftrajectoriesclustered;Exit3、分簇(1)、確定簇的數量就是看兩個指標：TSV變化圖和TSV變化的變化率。(2)、通常，TSV變化率大于30%表示“不同”簇被配對了，也就是說匹配前的簇數是可能的最終結果。但是，如果30%未識別任何最終簇數，則可以使用20%的標準。(3I單擊運行按鈕以生成可能的最終簇數的文本列表(CLUSEND)上述結果一列為簇數，一列為該簇數的迭代的TSV變化相對于前一個簇數的迭代的TSV變化的變化率。（4）、簇數選擇：在前幾個聚類迭代中，TSV的變化非常大，那么對于大部分聚類來說，隨著聚類數量的減少，TSV的增加速率通常很小，并且幾乎不變。在某一時刻，TSV的變化迅速上升，表明組合的簇不是非常相似。通常會有幾個較大的上升點。在這些TSV大幅增加的點之前的聚類數量給出了可能的最終結果。所以numbersofcluster可以選擇5條。4、分簇結果確定簇數后，點擊運行即可。垂直氣層高度選擇metersAGL,完成后點擊退出即可這個時候聚類就已經完成了。我們可以點擊‘DisplayMeans”查看分類結果，也可以點擊“DisplayClusters”查看單個簇中的所有軌跡。五、結果分析（1）聚類結果會生成在working文件夾里。主要有以下幾個結果文件：①CLUSLIST_N（N是聚類的簇數）。②C1_Nmean.tdump（簇數是多少,就會有幾個這樣的文件,C1_N一直到CN_N,介紹的是具體每個簇的信息）。③Cmean1_N.tdump（每個簇的平均情況）。④clusmean1_N（簇的分布圖）。最后一步匹配濃度的時候，需要用的文件是CLUSLIST_N。（2）將聚類結果與污染物濃度數據結合分析打開CLUSLIST_N文件。文件中會列出聚類時間段里的所有軌跡，以及他們所屬的cluster。將文件中數據