馬 雯，丁建麗?，白 婷

(1. 新疆大學(xué)地理與遙感科學(xué)學(xué)院智慧城市與環(huán)境建模自治區(qū)普通高校重點(diǎn)實(shí)驗(yàn)室，新疆烏魯木齊 830017；2. 新疆大學(xué)綠洲生態(tài)教育部重點(diǎn)實(shí)驗(yàn)室，新疆烏魯木齊 830017)

0 引言

大氣氣溶膠是由懸浮在大氣中的微小顆粒（液體和固體）組成的多相體系，其粒徑在0.001～100 μm之間[1]．雖然氣溶膠在大氣中只占一小部分，但它在氣候變化中扮演著重要的角色[2]．氣溶膠通過(guò)散射和吸收太陽(yáng)輻射直接改變地球大氣系統(tǒng)的能量平衡[3?4]，還能通過(guò)改變?cè)频奈⑽锢硖匦浴勖敖邓^(guò)程來(lái)間接影響氣候[5?6]．因此，在全球和區(qū)域范圍內(nèi)已經(jīng)開展了一系列研究來(lái)估算氣溶膠光學(xué)特性，例如氣溶膠光學(xué)厚度（AOD）、Angstr¨om波長(zhǎng)指數(shù)（AE）和單次散射反照率（SSA）[7?9]．主要是通過(guò)地基和星載儀器對(duì)氣溶膠進(jìn)行長(zhǎng)期的觀測(cè)和量化，詳細(xì)了解氣溶膠的時(shí)空變化，但氣溶膠仍然是輻射傳輸和氣候變化評(píng)估中最大的不確定因素之一[10]．

近些年來(lái)，國(guó)內(nèi)外學(xué)者針對(duì)氣溶膠與降水相互作用做了很多研究，認(rèn)為氣溶膠的增加既能促進(jìn)降水，也能抑制降水[11]．Kuhlmann等[12]的研究指出，亞洲地區(qū)人類活動(dòng)排放的氣溶膠減少了該區(qū)域的降水量．張喆等[13]的研究也提到2002―2012年間新疆地區(qū)AOD的增加抑制了降水的產(chǎn)生．然而高星星等[14]則指出，氣溶膠的增加雖然抑制了塔克拉瑪干沙漠地區(qū)的夏季降水，卻促進(jìn)了中國(guó)其它區(qū)域的夏季降水，且主要影響大雨發(fā)生概率．另外在沙塵暴天氣條件下，降水不但會(huì)提前出現(xiàn)，還使累計(jì)降水增加47%[15]．由于氣溶膠對(duì)氣候的多重影響，氣溶膠的增加或減少也可能對(duì)區(qū)域植被生長(zhǎng)產(chǎn)生影響．一方面氣溶膠沉降對(duì)植物的影響，往往由于大氣沉降對(duì)植物的生態(tài)效應(yīng)很復(fù)雜且難以識(shí)別；另一方面氣溶膠通過(guò)改變輻射平衡影響植物的生長(zhǎng)，有促進(jìn)也有抑制[16]．但目前只有少數(shù)研究探討了氣溶膠對(duì)植被生長(zhǎng)方面的影響，很少研究氣溶膠與氣象參數(shù)和植被生長(zhǎng)的相互作用．

新疆作為絲綢之路經(jīng)濟(jì)帶核心區(qū)，是典型的干旱半干旱區(qū)域．區(qū)域內(nèi)有古爾班通古特沙漠和塔克拉瑪干沙漠，使其成為沙塵氣溶膠的主要來(lái)源之一[17?18]．隨著新疆首府快速的城市化和工業(yè)化，空氣質(zhì)量嚴(yán)重惡化，這可能是由氣溶膠引起的．然而針對(duì)新疆氣溶膠的研究，多集中于分析AOD的時(shí)空特征[19?20]及其與地面顆粒物質(zhì)量濃度的關(guān)系[21?23]，對(duì)于新疆地區(qū)氣溶膠變化對(duì)植被和降水的影響研究仍然缺乏．因此，本文利用2003―2019年MODIS 3 km氣溶膠反演數(shù)據(jù)在深度分析新疆地區(qū)氣溶膠光學(xué)特性時(shí)間變化的基礎(chǔ)上，進(jìn)一步對(duì)氣溶膠進(jìn)行了分類．然后，從時(shí)間尺度上探討了AOD與不同時(shí)期歸一化植被指數(shù)（NDV I）和降水量（P）的關(guān)系，以期為我國(guó)西北大范圍干旱半干旱氣候背景下的氣候變化模擬、實(shí)施污染控制和相關(guān)產(chǎn)業(yè)轉(zhuǎn)型提供科學(xué)依據(jù)．

1 資料與方法

1.1 研究區(qū)概況

新疆地處中國(guó)西北部，屬于溫帶大陸性氣候，氣候干燥，全年降水量大多小于150 mm，區(qū)域內(nèi)沙漠、戈壁分布較廣．以天山山脈為界將新疆分為北疆和南疆，處于東天山的東疆一般由哈密市和吐魯番市組成．在南疆有塔克拉瑪干沙漠[24]，是中國(guó)最大的沙塵源區(qū)．產(chǎn)生的沙塵氣溶膠對(duì)區(qū)域氣候變化和人類生產(chǎn)生活有著重要影響．在北疆和東疆多為山谷和盆地，植被覆蓋率相對(duì)較高，城市區(qū)域相對(duì)集中，人類活動(dòng)密度大，氣溶膠分布復(fù)雜．考慮到新疆復(fù)雜的氣候背景條件和“三山兩盆”的獨(dú)特地形特征，本文將新疆分為三個(gè)區(qū)域進(jìn)行分析：北疆（North）；南疆（South）；東疆（East）．

1.2 數(shù)據(jù)來(lái)源

1.2.1 氣溶膠光學(xué)厚度

表 1 本文所使用的太陽(yáng)光度計(jì)站點(diǎn)信息

使用2003年1月至2019年12月期間空間分辨率為3 km（MYD043K）的MODIS/Aqua每日氣溶膠產(chǎn)品．為了驗(yàn)證MODIS AOD產(chǎn)品在新疆地區(qū)的質(zhì)量，選擇了3個(gè)地基遙感觀測(cè)的AOD數(shù)據(jù)．地面觀測(cè)點(diǎn)的位置如表1所示，精河站點(diǎn)位于精河氣象局，使用Cimel CE318光度計(jì)（CE318-T）進(jìn)行定點(diǎn)連續(xù)觀測(cè)；新疆大學(xué)站點(diǎn)位于烏魯木齊市，使用Microtops II太陽(yáng)光度計(jì)進(jìn)行定點(diǎn)連續(xù)觀測(cè)；喀什站點(diǎn)位于喀什市中國(guó)科學(xué)院遙感與數(shù)字地球研究所園區(qū)內(nèi)，使用全自動(dòng)太陽(yáng)天空輻射計(jì)（CE318-DP）進(jìn)行定點(diǎn)連續(xù)觀測(cè)．該數(shù)據(jù)已被納入中國(guó)建設(shè)的地面遙感觀測(cè)網(wǎng)絡(luò)（Sun-sky Radiometer Observation Network, SONET），SONET提供的大氣氣溶膠特性參數(shù)，可用于衛(wèi)星遙感產(chǎn)品驗(yàn)證[25?26]，該數(shù)據(jù)可從SONET網(wǎng)站（http://www.sonet.ac.cn/）獲取．這項(xiàng)驗(yàn)證工作可以彌補(bǔ)AERONET在中國(guó)西北地區(qū)站點(diǎn)分布較少的不足．

1.2.2 其它輔助數(shù)據(jù)

為了探討新疆地區(qū)大氣氣溶膠變化對(duì)植被和降水的影響，收集了與氣溶膠同一時(shí)間段內(nèi)的3種輔助數(shù)據(jù)．第一，利用MODIS C6歸一化植被指數(shù)（NDV I）3級(jí)月度數(shù)據(jù)，用于分析新疆上空AOD對(duì)植被生長(zhǎng)的影響．第二，從熱帶降水測(cè)量計(jì)劃（TRMM3B43）中獲取月度降水?dāng)?shù)據(jù)集，用于描述降水和AOD之間的關(guān)系．第三，新疆每年的GDP數(shù)據(jù)從《新疆統(tǒng)計(jì)年鑒》獲?。?/p>

1.3 研究方法

由于地基觀測(cè)儀器不能直接獲取550 nm波長(zhǎng)的AOD數(shù)據(jù)，為了與MODIS產(chǎn)品進(jìn)行比較，使用AE算法得到550 nm波長(zhǎng)的AOD值[27?28]．為了進(jìn)行時(shí)空匹配，以每個(gè)地面站點(diǎn)為中心的25 km×25 km窗口內(nèi)的MYD043K AOD求取均值，與MODIS/Aqua過(guò)境時(shí)間±0.5 h內(nèi)的地基AOD數(shù)據(jù)的平均值進(jìn)行比較．采用以下統(tǒng)計(jì)參數(shù)進(jìn)行評(píng)估，包括相關(guān)系數(shù)（R）、平均絕對(duì)誤差（MAE）、均方根誤差（RMSE）、相對(duì)平均偏差（RMB），產(chǎn)品期望誤差（EE）為± (0.05+0.2×AOD地基)．

2 結(jié)果與分析

2.1 MODIS AOD與地基AOD對(duì)比分析

如圖1所示，通過(guò)時(shí)空、波段匹配得到符合要求的數(shù)據(jù)共有1 285組．其中：落在期望誤差范圍內(nèi)的數(shù)據(jù)占比為31.52%，有60.39%的數(shù)據(jù)低于期望誤差，表明該產(chǎn)品在新疆地區(qū)存在一定的低估，這主要與地面觀測(cè)點(diǎn)對(duì)地表反射率的高估有關(guān)[29]．但MYD043K AOD與地基觀測(cè)AOD數(shù)據(jù)整體趨勢(shì)較為一致，二者的相關(guān)系數(shù)為0.84，均方根誤差為0.37．這與孫曉雷等[30]對(duì)該產(chǎn)品在阿克蘇地區(qū)的驗(yàn)證結(jié)果相似（R=0.883 6），進(jìn)一步說(shuō)明MYD043K的AOD數(shù)據(jù)可以用于新疆地區(qū)氣溶膠的研究．

圖 1 新疆地區(qū)MODIS AOD與地基AOD的驗(yàn)證結(jié)果

2.2 氣溶膠光學(xué)特性的時(shí)間變化

2.2.1 氣溶膠光學(xué)特性和經(jīng)濟(jì)指標(biāo)的年際變化

圖2為新疆地區(qū)AOD、AE和經(jīng)濟(jì)發(fā)展指標(biāo)（GDP）的年際變化．2003―2019年AOD的均值為0.26，范圍在0.19～0.29之間，且以0.008/10 a的速率呈微弱遞減趨勢(shì)．AOD峰值出現(xiàn)在2006年和2014年，主要是由于頻繁的沙塵天氣．而在2006―2009年和2014―2017年期間，AOD呈現(xiàn)下降趨勢(shì)，這可能與同期GDP有一定的關(guān)系，以及與2013年以來(lái)新疆大部分地區(qū)實(shí)施的重點(diǎn)大氣污染防治項(xiàng)目有關(guān)．AE的年均值在0.91～1.03之間，平均值為0.98，整體上以0.027/10 a的速率呈下降趨勢(shì)．特別是2012年以后AE的下降趨勢(shì)，說(shuō)明提倡的綠色生態(tài)對(duì)該區(qū)域人為氣溶膠的減少有很大作用．

圖 2 2003―2019年新疆地區(qū)AOD、AE和GDP的年際變化

進(jìn)一步采用Morlet小波分析對(duì)2003―2019年新疆地區(qū)AOD的周期性進(jìn)行探討．利用Matlab軟件繪制小波系數(shù)的實(shí)部圖和小波方差圖（圖3）．其中：實(shí)部圖中的小波系數(shù)與AOD值之間是正相關(guān)關(guān)系，也就是說(shuō)小波系數(shù)越小，對(duì)應(yīng)的AOD值越小，反之越大；從小波系數(shù)實(shí)部圖中可以看出，AOD有5～8 a和9～12 a兩個(gè)振蕩期，其高值中心分別是3 a、6 a和10 a．再結(jié)合小波方差可以得到AOD隨時(shí)間變化的周期特征，即峰值最高的數(shù)值是主周期．因此，由圖3（b）可知，在時(shí)間尺度為10 a時(shí)，周期振蕩最強(qiáng)，為第一主周期．

圖 3 AOD的小波系數(shù)實(shí)部圖和小波方差圖

根據(jù)小波分析得出的10 a主周期，分別分析了2012年以前和2012年以后兩個(gè)時(shí)期（2003―2012年和2013―2019年）新疆地區(qū)不同分區(qū)內(nèi)AOD和AE的變化規(guī)律．由圖4可知，在第一階段（2003―2012年），三個(gè)分區(qū)的AOD均呈上升趨勢(shì)（分別為0.004 1/a、0.002 7/a和0.003 3/a），而AE在南疆和東疆呈下降趨勢(shì)，在北疆呈上升趨勢(shì)．在第二階段（2013―2019年），東疆和南疆的AOD呈現(xiàn)下降趨勢(shì)，而AE則呈現(xiàn)上升趨勢(shì)．與之相反，北疆的AOD以0.002 1/a的速率增長(zhǎng)，而AE以0.000 5/a的速率下降．整體而言，2012年以前AOD以0.003 7/a 的速率呈增長(zhǎng)趨勢(shì)，AE則以0.004 9/a的速率呈下降趨勢(shì)；而2012年以后AOD則以0.002 9/a的速率呈下降趨勢(shì)，AE則以0.001 4/a的速率呈上升趨勢(shì)．

圖 4 AOD和AE在2012年前和2012年后的變化趨勢(shì)

2.2.2 氣溶膠光學(xué)特性的季節(jié)和月度變化

由圖5可知，新疆不同區(qū)域的AOD和AE表現(xiàn)出明顯的月變化，其中南疆、北疆和東疆具有相似的月變化特征．整體而言，新疆地區(qū)的AOD峰值（0.41±0.09）出現(xiàn)在4月，AE的峰值（1.28±0.12）出現(xiàn)在12月．而AOD的最小值（0.13±0.03）出現(xiàn)在12月，AE的最小值（0.71±0.06）出現(xiàn)在3月．AOD和AE在月際變化特征中存在一定的負(fù)相關(guān)關(guān)系. 這樣的變化特征主要?dú)w因于不同氣溶膠類型及氣溶膠的來(lái)源．

從季節(jié)上看，最高的AOD值出現(xiàn)在春季（0.36±0.03），其次是夏季（0.22±0.03）、冬季（0.19±0.05）和秋季（0.14±0.02）．最大的AE值出現(xiàn)在冬季（1.15±0.03），其次是秋季（1.09±0.15）、夏季（0.88±0.05）和春季（0.76±0.06）．由于特殊的氣候條件，春季溫度迅速回升，頻繁的大風(fēng)天氣增加了土壤粉塵，導(dǎo)致AOD值最大．此外，春季是中國(guó)西北地區(qū)揚(yáng)塵運(yùn)輸?shù)母叻迤冢死敻缮衬侵袊?guó)西北干旱和半干旱地區(qū)天然氣溶膠的主要來(lái)源[31]；因此，頻繁的沙塵天氣導(dǎo)致春季的AOD高于秋季和冬季．東疆和南疆的變化特征與全疆一致，而北疆AOD和AE最高值均出現(xiàn)在冬季，這可能是由于該區(qū)域冬季燃煤產(chǎn)生的細(xì)顆粒氣溶膠產(chǎn)生的影響．

圖 5 2003―2019年新疆地區(qū)AOD和AE月度和季節(jié)變化分布

2.3 氣溶膠類型的分類及其貢獻(xiàn)

圖 6 2003―2019年新疆地區(qū)主要?dú)馊苣z類型的季節(jié)和年度分類及其貢獻(xiàn)率

氣溶膠根據(jù)顆粒物的大小和吸收特性的差異分為不同的類型. 目前，已有許多學(xué)者根據(jù)AOD和AE之間關(guān)系對(duì)不同地區(qū)的氣溶膠進(jìn)行分類[32?36]．為探討中國(guó)新疆區(qū)域污染物的主要來(lái)源和成分，我們基于現(xiàn)有數(shù)據(jù)源和地理位置，主要參考了中國(guó)長(zhǎng)三角地區(qū)[37?38]、中亞五國(guó)[33]及伊朗[34]的分類閾值，實(shí)現(xiàn)對(duì)氣溶膠的初步分類．基于上述研究，本文對(duì)AOD和AE的選擇閾值如下：（1）AOD<0.2時(shí)為大陸清潔氣溶膠（CC），表示新疆地區(qū)氣溶膠的背景條件；（2）AOD>0.3和AE>1.0時(shí)為生物質(zhì)燃燒/城市-工業(yè)型的氣溶膠（BUI）；（3）AOD>0.5和AE<0.7時(shí)為沙塵氣溶膠（DD）；（4）不屬于上述任何一組的情況，被定義為混合型氣溶膠（MX）．這樣的氣溶膠分類閾值設(shè)置在Bibi等[39]的研究中已經(jīng)進(jìn)行了精度驗(yàn)證，證明了該閾值設(shè)置能夠有效地區(qū)分氣溶膠的類型．

圖6（a～e）為新疆地區(qū)不同氣溶膠類型的年度和季節(jié)分布，圖6（f）中的柱狀圖為不同氣溶膠類型的貢獻(xiàn)百分比．所有類型的氣溶膠在4個(gè)季節(jié)都有不同程度的貢獻(xiàn)．春季，DD氣溶膠類型最為豐富（24%），明顯高于其它季節(jié)，這主要是由于春季為全年中沙塵天氣的頻發(fā)期，大粒徑的沙塵粒子占比相對(duì)較大．夏季，BUI氣溶膠類型的比例相對(duì)較高，這可能與夏季的氣候條件有關(guān)，其中氣溶膠的二次轉(zhuǎn)化具有一定作用．秋季，CC氣溶膠類型的占比在四季中最大，達(dá)到59%，說(shuō)明該季節(jié)的空氣最為清潔．冬季，CC氣溶膠類型的占比相對(duì)較大，這與趙仕偉等[40]對(duì)西北地區(qū)的研究有一定的差異，可能是由于該區(qū)域冬季有很多的空值，數(shù)據(jù)處理時(shí)用最小值進(jìn)行填充導(dǎo)致的．總的來(lái)說(shuō)，在整個(gè)研究期，CC氣溶膠類型的貢獻(xiàn)率最大（47%），說(shuō)明新疆地區(qū)全年氣溶膠以小粒徑的顆粒物為主；其次為MX氣溶膠類型，貢獻(xiàn)率達(dá)到33%．

2.4 氣溶膠變化對(duì)降水和植被的影響分析

2.4.1 AOD對(duì)降水的影響

為探討氣溶膠對(duì)降水的影響，比較了新疆地區(qū)每月的AOD和降水．由于TRMM降水?dāng)?shù)據(jù)在新疆有一定的偏差，必須在數(shù)據(jù)使用前進(jìn)行校正[41]．參考已有的研究，采用加法修正模型對(duì)TRMM的原始降水?dāng)?shù)據(jù)進(jìn)行了校正[42]．因此，使用校正后的TRMM數(shù)據(jù)進(jìn)行以下分析．新疆降水的分布特點(diǎn)是“北多南少”，降水主要集中在阿爾泰山、天山和昆侖山等山區(qū)．

由圖7（a）可知，AOD與降水之間有一定的滯后性．為了解新疆降水和AOD之間的完整關(guān)系，必須考慮每月的時(shí)滯效應(yīng)．我們計(jì)算了AOD與當(dāng)月、前1個(gè)月、前2個(gè)月和前3個(gè)月降水之間的相關(guān)性．由于新疆的降水與NDV I之間存在顯著的正相關(guān)關(guān)系（R=0.719，P <0.01），因此采用偏相關(guān)分析方法來(lái)分析AOD和降水、NDV I之間的相關(guān)性．結(jié)果表明：AOD與降水的關(guān)系較為復(fù)雜（表2）．整個(gè)新疆地區(qū)AOD與降水的偏相關(guān)系數(shù)的最大值出現(xiàn)在當(dāng)前月，說(shuō)明AOD對(duì)降水沒(méi)有時(shí)滯影響．分區(qū)來(lái)看，南疆的AOD和降水也表現(xiàn)出同樣的特點(diǎn)．東疆的AOD與降水的偏相關(guān)系數(shù)顯示AOD對(duì)降水有1個(gè)月的滯后，而北疆則有2個(gè)月的滯后．

圖 7 2003―2019年NDV I、降水和AOD的月均值分布

表 2 新疆地區(qū)各月AOD和NDV I與降水的偏相關(guān)關(guān)系

2.4.2 AOD對(duì)NDV I的影響

新疆西部和西北部地區(qū)，如天山南北坡、阿爾泰山和塔里木盆地周圍的綠洲，植被覆蓋度較高，具有高的NDV I值．而在塔克拉瑪干沙漠、庫(kù)姆塔格沙漠和羅布泊周圍，植被稀少，是NDV I的低值區(qū)．2003―2019年，新疆地區(qū)月平均NDV I與AOD的相關(guān)系數(shù)為0.26（P <0.01）．該結(jié)果與He等[43]的研究一致，在全國(guó)范圍內(nèi)，AOD和NDV I之間沒(méi)有明顯的關(guān)系．然而，在陜西[44]和長(zhǎng)江流域[45]，發(fā)現(xiàn)AOD與NDV I之間有明顯的負(fù)相關(guān)關(guān)系．這意味著植被有能力防止顆粒物從土壤中被吹入空氣．但由于新疆大部分地區(qū)被稀疏的植被覆蓋，這可能大大降低了植被的吸收能力，導(dǎo)致AOD和NDV I之間的關(guān)系很弱．

由圖7（b）可知，AOD和NDV I之間存在一定的滯后．用與降水量相同的方法分析AOD與NDV I的關(guān)系，發(fā)現(xiàn)AOD與NDV I呈負(fù)相關(guān)關(guān)系．整個(gè)新疆地區(qū)AOD與NDV I的最大偏相關(guān)系數(shù)出現(xiàn)在前3個(gè)月，說(shuō)明新疆地區(qū)AOD對(duì)NDV I的響應(yīng)存在3個(gè)月的滯后期．分區(qū)來(lái)看，南疆和東疆也存在同樣的滯后期．然而，北疆的AOD與NDV I的最大偏相關(guān)系數(shù)是在前1個(gè)月，表明北疆的AOD與NDV I之間存在1個(gè)月的滯后．

3 結(jié)論

利用MODIS氣溶膠產(chǎn)品數(shù)據(jù)，進(jìn)行了AOD和AE的時(shí)間特征分析，并確定了氣溶膠污染的主要類型，探討了2003―2019年新疆地區(qū)氣溶膠對(duì)降水和植被的影響，得出以下結(jié)論：

（1）MYD043K AOD與地基AOD的相關(guān)系數(shù)為0.84，表明該產(chǎn)品適用于新疆地區(qū)．研究期間，AOD和AE呈不同程度的下降趨勢(shì)．從季節(jié)尺度看，AOD在春季最高、秋季最低，而AE與之相反，說(shuō)明AOD和AE存在一定的負(fù)相關(guān)關(guān)系．

（2）根據(jù)AOD和AE二者的關(guān)系，嘗試進(jìn)行氣溶膠類型的劃分．沙塵氣溶膠類型在春季的貢獻(xiàn)最大，這主要與頻發(fā)的沙塵天氣有關(guān)；總體來(lái)看，大陸清潔型氣溶膠在全年貢獻(xiàn)最大，混合型氣溶膠次之．

（3）從月尺度分析，AOD與降水和NDV I有一定的相關(guān)性．其中：AOD與NDV I呈負(fù)相關(guān)，而AOD與降水的關(guān)系比較復(fù)雜．通過(guò)AOD對(duì)NDV I和降水的時(shí)滯效應(yīng)分析，發(fā)現(xiàn)AOD對(duì)NDV I的時(shí)滯為1～3個(gè)月，AOD對(duì)降水量的時(shí)滯為0～2個(gè)月．