趙俊芳，徐 慧，孔祥娜,2，郭建平，余衛(wèi)國

?

基于MODIS和AERONET的氣溶膠地表直接輻射效應(yīng)評價*

趙俊芳1，徐 慧1，孔祥娜1,2，郭建平1，余衛(wèi)國3

（1．中國氣象科學(xué)研究院災(zāi)害天氣國家重點實驗室，北京 100081；2．四川農(nóng)業(yè)大學(xué)資源學(xué)院，成都 611130；3．安徽省農(nóng)業(yè)氣象中心，合肥 232001）

氣候變化背景下進一步開展高污染地區(qū)氣溶膠直接輻射效應(yīng)的研究，對于該地區(qū)的大氣環(huán)境監(jiān)測、氣候變化評估以及農(nóng)業(yè)生產(chǎn)布局等具有重要意義。利用全球氣溶膠監(jiān)測網(wǎng)（AERONET）多年觀測資料以及MODIS地表反照率數(shù)據(jù)，借助6S（Second Simulation of a Satellite Signal in the Solar Spectrum）輻射傳輸模式，定量評估2001年1月?2016年12月中國AERONET北京、香河和太湖3個典型高污染站點多年晴空條件下的氣溶膠直接輻射效應(yīng)。結(jié)果表明：（1）各個高污染站點氣溶膠PM2.5質(zhì)量濃度季節(jié)變化特征差異明顯，秋冬污染較為嚴(yán)重，其中高值主要集中在1月、2月、11月和12月；（2）與無氣溶膠影響相比，氣溶膠致使各站點地表直接輻射年際變化較明顯。在氣溶膠影響下，北京、香河和太湖3個站點2001?2016年地表日平均直接輻射年際變化較明顯，均呈波動增加趨勢，且香河站日平均直接輻射增加最大（621.14W·m?2），太湖站次之（743.29W·m?2），北京站最?。?88.14W·m?2）。（3）氣溶膠影響下各站點地表直接輻射明顯降低，且氣溶膠對各站點影響差異較大。2001?2016年，北京、太湖和香河站年平均地表直接輻射分別降低32.29%、24.01%和15.07%。其中，氣溶膠對北京站的地表直接輻射影響最大，香河站最小。（4）近15a來，北京、香河和太湖3個站點氣溶膠地表輻射效應(yīng)均呈現(xiàn)增加趨勢。

MODIS；AERONET；氣溶膠；PM2.5；直接輻射效應(yīng)

大氣氣溶膠是指懸浮在大氣中的固態(tài)和液態(tài)顆粒物的總稱，主要是指六大類7種氣溶膠粒子，包括沙塵氣溶膠、碳氣溶膠（黑碳和有機碳氣溶膠）、硫酸鹽氣溶膠、硝酸鹽氣溶膠、銨鹽氣溶膠和海鹽氣溶膠[1]。大氣氣溶膠的氣候效應(yīng)早在20世紀(jì)70年代已被意識到[2]，80年代末開始受到科學(xué)界的普遍關(guān)注[3]。氣溶膠通過兩條途徑影響氣候系統(tǒng)：一是通過吸收和反射太陽與地球輻射影響氣候系統(tǒng)。這種影響稱之為氣溶膠的直接輻射效應(yīng)；另一種是作為云滴中的云凝結(jié)核改變云的光學(xué)特性和生命周期，這種影響成為氣溶膠的間接輻射效應(yīng)。氣溶膠的輻射效應(yīng)不僅通過改變光合有效輻射的總量，還通過改變光譜分布、光的質(zhì)量以及在直接和散射成分之間的分配，來影響植被的光合作用，并最終影響產(chǎn)量。因此，氣候變化背景下大氣氣溶膠對大氣輻射的影響越來越受到科學(xué)家們的重視[4?6]。

氣溶膠的輻射效應(yīng)是指由氣溶膠所引起的出射輻射通量的變化，通常可分為直接輻射效應(yīng)和間接輻射效應(yīng)[7]。其中，氣溶膠直接輻射效應(yīng)指大氣中氣溶膠的含量變化引起地－氣系統(tǒng)能量平衡的擾動，在數(shù)值上定義為氣溶膠含量變化時所產(chǎn)生的平均凈輻射變化[8]，也就是在有氣溶膠和無氣溶膠存在時輻射通量的變化。在氣溶膠的輻射效應(yīng)評估方面，國內(nèi)外近年來一直大力發(fā)展衛(wèi)星遙感與地基光學(xué)遙感相結(jié)合的觀測技術(shù)與反演方法[7?10]。中高分辨率的成像光譜儀（MODIS、HIRDLS、GLI）、多角度的成像光譜儀（MISR）等星載探測器的一個重要應(yīng)用方向正是探測全球（包括陸地）的氣溶膠光學(xué)特性。在地基遙感方面，近年來較重要的研究進展之一是建立了AERONET（AErosol RObotic NET work）太陽光度計全球探測網(wǎng)絡(luò)由此提取氣溶膠的光學(xué)特性，并用于對衛(wèi)星遙感氣溶膠的檢驗[11?12]。Boiyo等[12]利用東非肯尼亞3個AERONET站點MODIS Terra、MISR和OMI數(shù)據(jù)反演了氣溶膠光學(xué)厚度，并對其反演精度進行了比較。張志薇等[10]借助平面平行輻射傳輸模式評估了蘭州大學(xué)半干旱氣候與環(huán)境觀測站SACOL、香河和太湖地區(qū)大氣層頂氣溶膠日平均直接輻射強迫分別為（?7.2±7.1）、（?11.7±9.5）和（?14.9±8.3）W·m?2，地表處分別為（?21.9±12.2）、（?36.5±29.5）和（?42.2±21.0）W·m?2，大氣中分別為（14.7±9.7）、（24.8±23.0）和（27.3±16.5）W·m?2。Tang等[2]分析了1981?2010年氣溶膠光學(xué)厚度變化及其對全國輻射的影響。陳林等[9]利用MODIS和MISＲ資料，以及RAD01輻射模式評估了中國陸地地區(qū)的氣溶膠直接輻射效應(yīng)，結(jié)果表明陸地上氣溶膠對大氣有很強的加熱作用，其強吸收性導(dǎo)致地面輻射強迫為大氣層頂輻射強迫的3倍。

由于氣溶膠是由具有不同譜分布、形狀、化學(xué)組成和光學(xué)性質(zhì)的物質(zhì)構(gòu)成的，且其濃度的時空變化可達幾個數(shù)量級而且缺乏充分的觀測資料，采用相對較新的觀測資料對氣溶膠的輻射效應(yīng)進行準(zhǔn)確模擬，有助于對當(dāng)前和未來的氣候變化進行評估。從目前情況看，多數(shù)研究針對全球或者大區(qū)域尺度，對高污染地區(qū)氣溶膠的輻射效應(yīng)研究報道仍有限。因此，本研究利用2001年1月?2016年12月全球氣溶膠監(jiān)測網(wǎng)（AERONET）多年觀測資料以及MODIS地表反照率數(shù)據(jù)，借助6S輻射傳輸模式，定量評價2001年1月?2016年12月中國AERONET北京、香河和太湖3個典型高污染站點多年晴空條件下的氣溶膠輻射效應(yīng)，以期對該地區(qū)的大氣環(huán)境監(jiān)測、污染治理以及氣候變化評估等提供理論支持。

1 資料與方法

1.1 研究區(qū)概況

選擇中國地區(qū)觀測時間序列較長，且區(qū)域代表性較強的北京、香河和太湖3個AERONET站點，站點信息見表1。北京站和香河站分別位于華北平原西北邊緣和北部，屬于溫帶季風(fēng)氣候，雨熱同期，夏季高溫多雨，冬季寒冷干燥；春、秋短促，水熱條件較好[13]，站點周圍均為農(nóng)業(yè)生產(chǎn)區(qū)，地表植被覆蓋較好，主要種植玉米、小麥、大豆等作物，是中國重要的糧棉油生產(chǎn)基地。太湖站位于太湖北側(cè)湖岸，氣候溫和濕潤，土壤肥沃，站點周圍是中國重要的商品糧基地和桑蠶基地，主要種植水稻、小麥、桑、茶葉等。近年來，隨著經(jīng)濟的快速發(fā)展和污染物排放大量增加，這些站點周圍空氣質(zhì)量惡化的問題日益突出[14?15]。

表1 研究所選AERONET站點信息

1.2 資料來源與處理

1.2.1 Aerosol Robotic Network（AERONET）數(shù)據(jù)

全球氣溶膠自動觀測網(wǎng)（AErosol RObotic NETwork，AERONET）是由美國NASA和法國國家科學(xué)研究中心（Centre National de la Recherche Scientifique，CNRS）發(fā)起，由世界大量研究機構(gòu)聯(lián)合建立的地基氣溶膠監(jiān)測網(wǎng)絡(luò)，該監(jiān)測網(wǎng)采用法國CIMEL公司生產(chǎn)的CE-318型全自動太陽直接/散射輻射計，對太陽7個波段（340、380、440、670、870、940和1020nm）的直接和散射輻射進行測量。CE-318既可以自動進行太陽直射輻射觀測，也可以進行平緯圈（ALM）和主平面（PPL）天空光掃描觀測。CE-318測得的太陽直射輻射數(shù)據(jù)可以用來計算氣溶膠光學(xué)厚度AOD、大氣透過率和水汽總量等信息，而主平面和平緯圈天空掃描數(shù)據(jù)可以用來反演氣溶膠粒子的單次散射反照率、不對稱因子和散射相函數(shù)等光學(xué)參數(shù)[16]。觀測儀器每15min采集一次340、380、440、500、675、870和1020nm通道處的氣溶膠光學(xué)特性。由于受地表參數(shù)和氣溶膠前向散射影響較小[17]，AERONET反演的氣溶膠光學(xué)厚度（Aerosol Optical Depth，AOD）反演精度可達0.02～0.01[18]，因而被較多地應(yīng)用為AOD基準(zhǔn)數(shù)據(jù)。其數(shù)據(jù)分為未篩選的Level 1.0產(chǎn)品、經(jīng)云檢測的Level 1.5產(chǎn)品，以及經(jīng)過云檢測和質(zhì)量控制的Level 2.0產(chǎn)品[16]。

本研究選取AERONET提供的2001年1月?2016年12月AOD數(shù)據(jù)。440nm波段AOD的瞬時數(shù)據(jù)及440?675nm波段的?ngstr?m指數(shù)指來源于云檢測和質(zhì)量控制的Level 2.0產(chǎn)品。各站觀測時段見表1。

1.2.2 MODIS數(shù)據(jù)

MODIS的工作波段涵蓋了0.553?14.235μm范圍內(nèi)的可見光、近紅外和紅外的36個波段[19]。目前運行的2個MODIS傳感器分別搭載在Terra和Aqua衛(wèi)星平臺上，提供了包括AOD產(chǎn)品MOD04（MODIS/Terra）和MYD04（MODIS/Aqua）2種主要氣溶膠數(shù)據(jù)產(chǎn)品。MODIS AOD在大氣環(huán)境、氣候變化等相關(guān)領(lǐng)域的研究已較成熟，其反演精度與空間分辨率較高，可直接提供550nm處的AOD值[6]。

采用質(zhì)量控制的Level 2.0產(chǎn)品MOD04和MYD04氣溶膠數(shù)據(jù)。該產(chǎn)品含有550nm通道處的AOD值，空間分辨率約為10km。

1.2.3 地表反照率數(shù)據(jù)

地表反照率表征地球表面對太陽輻射的反射能力，是影響氣溶膠輻射效應(yīng)的重要因素。地表反照率數(shù)據(jù)來源于MODIS的MCD43地表反照率產(chǎn)品。MCD43分500m、1km和0. 05°氣候網(wǎng)格等3種空間分辨率（分別編號A、B和C），提供BRDF模型參數(shù)、BRDF反演質(zhì)量、反照率和BRDF校正后的天底反射率NBAR（Nadir BRDF-Adjusted Reflectance）等4種產(chǎn)品（依次編號1?4），其中，反照率產(chǎn)品提供MODIS 1?7通道和可見光VIS（0.3?0.7μm）、近紅外NIR（0.7?5.0μm）和短波SW（0.3?5.0μm）等3個寬波段的黑空反照率BSA（Black-Sky Albedo）和白空反照率WSA（White-Sky Albedo）[20]。研究數(shù)據(jù)為空間分辨率0.05°、正弦投影的MCD43地表反照率產(chǎn)品。

1.3 基于6S輻射傳輸模式的氣溶膠直接輻射效應(yīng)模擬

采用6S輻射傳輸模式（Second Simulation of a Satellite Signal in the Solar Spectrum）[21]模擬到達地球表面向下輻射通量。6S輻射傳輸模型主要用來模擬無云（晴空）天氣條件下，在太陽反射波段（0.25～4.0μm）衛(wèi)星傳感器理論上應(yīng)接收到的輻射值，不但能較合理地處理空氣分子和氣溶膠的散射和吸收，而且還能產(chǎn)生連續(xù)光譜[22]。該模型對大氣的影響（二氧化碳、氧氣、臭氧、水汽等的氣體吸收和分子散射作用）、非均一性地表、二向性反射率以及邊界條件等都有細致的考慮和運算，且充分利用了分析表達式和預(yù)選大氣模式，從而使計算時間大為縮短，被廣泛用于模擬通過耦合表面大氣介質(zhì)傳輸?shù)男l(wèi)星傳感器觀測到的光譜、幾何和大氣情況的信號[23]。

基于AERONET的AOD小時數(shù)據(jù)獲取AOD日均值數(shù)據(jù)。分別輸入AERONET的AOD小時數(shù)據(jù)、Terra或Aqua的AOD數(shù)據(jù)，以及地表反照率數(shù)據(jù)，使用6S輻射傳輸模型進行地表的輻射通量密度和氣溶膠直接輻射效應(yīng)的模擬。在模擬實驗中，以小時為時間間隔，獲取氣溶膠直接輻射效應(yīng)的日變化，最后對小時的氣溶膠直接輻射效應(yīng)進行積分，獲取有氣溶膠和無氣溶膠存在時輻射通量變化的日均值。晴空條件下地表（SFC）氣溶膠直接輻射效應(yīng)指在晴空條件下，地面短波輻射通量在分別考慮氣溶膠存在和不考慮氣溶膠存在時兩者間的差值[9]，即

ΔFSFC=F(a)SFC?FSFC'（1）

FSFC= FSFC↓?FSFC↑ （2）

式中，ΔFSFC表示地表的氣溶膠直接輻效應(yīng)；F(a)SFC表示地表有氣溶膠時的凈輻射通量；FSFC表示地表沒有氣溶膠時的凈輻射通量；FSFC↓、FSFC↑分別表示地表向下和向上方向的輻射通量。大氣氣溶膠直接輻射效應(yīng)值可正可負，取決于氣溶膠粒子自身反射和對太陽輻射吸收的相對能力以及地面反照率等因素。

日均氣溶膠直接輻射效應(yīng)為1d內(nèi)的瞬時輻射效應(yīng)的平均值[15]，即

1.4 氣溶膠PM2.5質(zhì)量濃度

大氣氣溶膠是指懸浮在大氣中的固態(tài)及液態(tài)微粒。PM2.5又稱細顆粒物、細粒、細顆粒，是指大氣中懸浮著的、空氣動力學(xué)當(dāng)量直徑小于或等于2.5μm的細顆粒物，濃度值越大說明空氣污染越嚴(yán)重。PM2.5通過對可見光的散射和吸收，可以削弱大氣能見度，影響人類正常生活和生產(chǎn)。PM2.5質(zhì)量濃度實測數(shù)據(jù)來自國家氣象信息中心。首先對其不可抗拒因素導(dǎo)致缺失的監(jiān)測值進行剔除，隨后對各站點多年細顆粒物質(zhì)量濃度數(shù)據(jù)進行統(tǒng)計分析，得到細顆粒物質(zhì)量濃度的日均值，用于PM2.5質(zhì)量濃度季節(jié)變化分析。

2 結(jié)果與分析

2.1 典型站點氣溶膠PM2.5質(zhì)量濃度季節(jié)變化特征分析

由圖1可以看出，各個典型站點氣溶膠PM2.5質(zhì)量濃度季節(jié)變化特征相似，秋冬污染較為嚴(yán)重。2001?2016年期間，北京站點各月PM2.5質(zhì)量濃度的變化范圍在79.06～113.17μg·m?3，平均92.96μg·m?3，其中高值主要集中在1月、2月、11月和12月。香河站各月PM2.5質(zhì)量濃度的變化范圍在62.37～161.32μg·m?3，平均95.31μg·m?3，其中高值主要集中在1月、2月、10月、11月和12月。太湖站各月PM2.5質(zhì)量濃度的變化范圍在27.33～93.33μg·m?3，平均54.50μg·m?3，其中高值主要集中在1月、2月、11月和12月。PM2.5質(zhì)量濃度變化有持續(xù)升高和降低的過程，反映了顆粒物的堆積和擴散的過程。11月開始北方地區(qū)進入冬季采暖期，燃煤量增大，造成了大氣中污染濃度值較高，PM2.5質(zhì)量濃度出現(xiàn)第一次峰值。而12月是冷空氣頻繁的月份，氣溫出現(xiàn)短暫回暖，但連續(xù)大風(fēng)降溫致使PM2.5質(zhì)量濃度出現(xiàn)低谷。之后一段時間內(nèi)，因冬季氣候干燥寒冷，空氣自凈能力差，逆溫強度大，使得顆粒物稀釋擴散能力降低，形成PM2.5質(zhì)量濃度整體上升趨勢，翌年2月達到第2次峰值，加之春節(jié)期間的煙花爆竹燃放使得污染加重，各地PM2.5質(zhì)量濃度均有所上升。3月中旬開始北方地區(qū)結(jié)束供暖，而且春季大風(fēng)天數(shù)增多，各地污染有所緩解。

圖1 北京、香河和太湖站點2001?2016年平均月PM2.5質(zhì)量濃度的變化特征

2.2 典型站點地表直接輻射的年際變化特征

2001?2016年，基于6S輻射傳輸模式模擬的中國AERONET三個典型高污染站點地表直接輻射的年際變化特征如圖2所示。由圖可見，無氣溶膠影響下，北京、香河和太湖3個站點2001?2016年地表日平均直接輻射年際變化不太明顯?？傮w來說，太湖站日平均直接輻射最大，北京站最小。其中，近15a來北京站地表日平均直接輻射在713.90～736.93W·m?2，平均720.68W·m?2（圖2a）；香河站為726.85～739.62W·m?2，平均731.38W·m?2（圖2b）；太湖站在963.09～996.83W·m?2，平均978.01W·m?2（圖2c）。

在氣溶膠影響下，北京、香河和太湖3個站點2001?2016年地表日平均直接輻射年際變化較明顯，均呈波動增加趨勢。近15a來北京、香河和太湖3個站地表日平均直接輻射增加趨勢分別為2.673、16.93和8.314W·m?2·a?1，且香河站日平均直接輻射增加最大，北京站最小。受氣溶膠的影響，近15a來北京站地表日平均直接輻射從423.57W·m?2增至588.13W·m?2，平均488.14W·m?2；香河站從465.87W·m?2增至712.66W·m?2，平均621.14W·m?2；太湖站從689.12W·m?2增至819.41W·m?2，平均743.29W·m?2。

圖2 有氣溶膠和無氣溶膠影響下北京、香河和太湖站日平均直接輻射的年際變化

2.3 氣溶膠對地表直接輻射的影響

基于6S輻射傳輸模式模擬結(jié)果顯示，氣溶膠對3個典型高污染站點地表直接輻射的影響不同，對北京站的影響最大，香河站最?。▓D3）。2001?2016年，北京、太湖和香河站平均地表直接輻射分別降低32.29%、24.01%和15.07%。氣溶膠的影響導(dǎo)致近15a北京站點地表直接輻射降低18.90%～40.80%，平均234.53W·m?2；香河站點降低2.86%～35.99%，平均110.24W·m?2；太湖站點降低17.80%～29.35%，平均232.72W·m?2。

圖3 典型站點由氣溶膠影響導(dǎo)致的日平均直接輻射降低率

2001?2016年，中國AERONET北京、香河和太湖3個典型站點氣溶膠地表直接輻射效應(yīng)分別平均為?207.45、?80.94和?216.71W·m?2。地表直接輻射效應(yīng)大的負值（絕對值）說明氣溶膠阻擋太陽光到達地表，對地表產(chǎn)生了顯著的降溫，即氣溶膠在地表有顯著的降溫作用[15]。近15a來，北京、香河和太湖3個站點氣溶膠地表直接輻射效應(yīng)呈現(xiàn)增加趨勢，對地表的降溫作用逐漸減弱，且北京和太湖站氣溶膠對地表的降溫作用遠遠大于香河站。其中，北京站氣溶膠地表直接輻射效應(yīng)從2001年的?214.24W·m?2增至2016年的?124.31W·m?2。太湖站從2005年的?219.23W·m?2增至2016年的?160.33W·m?2。香河站從2001年的?178.78W·m?2增至2016年的?18.29W·m?2。

3 結(jié)論與討論

3.1 結(jié)論

（1）北京、香河和太湖3個高污染站點氣溶膠PM2.5質(zhì)量濃度季節(jié)變化特征差異明顯，秋冬污染較為嚴(yán)重，其中高值主要集中在1月、2月、11月和12月。2001?2016年期間，北京站點各月PM2.5質(zhì)量濃度平均為92.96μg·m?3，香河站各月PM2.5質(zhì)量濃度平均為95.31μg·m?3，太湖站各月PM2.5質(zhì)量濃度平均為54.50μg·m?3。各個站點秋冬污染較為嚴(yán)重，高值主要集中在1月、2月、11月和12月。

（2）氣溶膠對各站點地表直接輻射年際變化的影響比較明顯。在氣溶膠影響下，北京、香河和太湖3個站點2001?2016年地表日平均直接輻射年際變化較明顯，均呈波動增加趨勢，增加趨勢分別為2.673、16.93和8.314W·m?2·a?1，且香河站日平均直接輻射增加最大（621.14W·m?2），北京站最小（488.14W·m?2）。然而，在無氣溶膠影響下，北京、香河和太湖3個站點2001?2016年地表日平均直接輻射年際變化不太明顯?？傮w來說，在無氣溶膠影響下太湖站日平均直接輻射最大（978.01W·m?2），香河站次之（731.38W·m?2），北京站最小（720.68W·m?2）。

（3）氣溶膠影響下各個站點地表直接輻射明顯降低，且對各站點影響差異較大。2001?2016年，北京、太湖和香河站平均地表直接輻射分別降低32.29%、24.01%和15.07%，即分別降低234.53、232.72和110.24W·m?2。其中，氣溶膠對北京站的地表直接輻射影響最大，太湖次之，香河站最小。

（4）近15a來，北京、香河和太湖3個站點氣溶膠地表輻射效應(yīng)均呈現(xiàn)增加趨勢，說明氣溶膠對地表的降溫作用逐漸減弱，且在北京和太湖站氣溶膠對地表的降溫作用遠遠大于香河站。其中，北京站氣溶膠地表直接輻射效應(yīng)從2001年的?214.24W·m?2增至2016年的?124.31W·m?2，太湖站從2005年的?219.23W·m?2增至2016年的?160.33W·m?2，香河站從2001年的?178.78 W·m?2增至2016年的?18.29W·m?2。

3.2 討論

氣溶膠在大氣中所占的質(zhì)量比例非常小，對輻射傳遞和氣候的影響卻非常顯著。決定氣溶膠直接輻射強迫大小的特性是它的散射系數(shù)、吸收系數(shù)、散射相函數(shù)，以及它們隨波長的變化[24]。因此要得到這些參數(shù)隨時間和空間變化的有代表性的資料，遙感觀測是取得系統(tǒng)數(shù)據(jù)的最有力手段。本研究利用全球氣溶膠監(jiān)測網(wǎng)（AERONET）多年觀測資料以及MODIS地表反照率數(shù)據(jù)，借助6S輻射傳輸模式，定量評價氣溶膠地表直接輻射效應(yīng)。結(jié)果表明，2001?2016年，在氣溶膠影響下北京、太湖和香河站平均地表直接輻射分別降低32.29%、24.01%和15.07%，近15a來各站點地表直接輻射效應(yīng)呈現(xiàn)增加趨勢，氣溶膠對地表的降溫作用逐漸減弱。以上結(jié)果與前人研究結(jié)果定性相一致[7?9]。由于目前遙感所得數(shù)據(jù)資料是以氣溶膠光學(xué)厚度為主，今后隨著新的衛(wèi)星遙感探測器的啟用，多波段、多角度、多時次、多偏振方向遙感信息的使用，完全有可能獲取全球范圍氣溶膠的各種輻射特性參數(shù)，相關(guān)研究將進一步深入。

本研究對研究區(qū)氣溶膠地表直接輻射強迫的評估還存在一定的不確定性：（1）目前對氣溶膠時空分布、自身的物理化學(xué)性質(zhì)和下墊面的光學(xué)性質(zhì)等缺乏徹底的了解?？臻g分布的不均勻性主要是由于氣溶膠在大氣中存在周期較短造成的，還取決于氣溶膠的產(chǎn)生機制和干、濕沉降過程[1]。環(huán)境相對濕度對氣溶膠輻射強迫效應(yīng)影響較大，為了更好地探討氣溶膠在環(huán)境條件下的輻射特性和氣候強迫效應(yīng)，需要更好地了解氣溶膠各組分光學(xué)性質(zhì)隨相對濕度的變化規(guī)律。（2）實際觀測的大氣氣溶膠（收集到膜上）通常包含許多種化學(xué)成分，然而單個氣溶膠有的可以只包含一種化學(xué)成分，即外部混合，或者很多種化學(xué)成分均勻地包含于一個單一氣溶膠中，即內(nèi)部混合，或者是部分外部混合，部分內(nèi)部混合，同時化學(xué)物種的分?jǐn)偙壤龝S不同粒子尺度而變化，也可以因不同地點、不同時間、或不同的來源而變化[25]。這些作用目前都沒能在氣溶膠數(shù)值模擬中得到很好的處理。隨著研究的不斷深入、資料的不斷完善和技術(shù)手段的不斷提高，今后將逐步降低研究的不確定性。（3）由于當(dāng)前大氣污染物主要地面監(jiān)測站分布集中在大中城市，數(shù)量少，檢測手段和精度也各不相同，不能滿足大面積、多種類的大氣污染監(jiān)測。結(jié)合衛(wèi)星遙感監(jiān)測手段與地面監(jiān)測站互為補充將成為未來大氣污染監(jiān)測預(yù)警的重要方式，但對氣溶膠污染的類型、光學(xué)厚度的估算等仍存在諸多不足之處[26]。因此，準(zhǔn)確監(jiān)測預(yù)測氣溶膠污染物的動態(tài)變化趨勢和幅度仍是當(dāng)前氣溶膠污染對農(nóng)作物影響研究中面臨的重要問題。

[1] 張小曳.中國大氣氣溶膠及其氣候效應(yīng)的研究[J].地球科學(xué)進展,2007,22(1):12?16.

Zhang X Y.Aerosol over China and their climate effect[J].Advances in Earth Science,2007,22(1):12?16.(in Chinese)

[2] Tang W J,Yang K,Qin J,et al.A revisit to decadal change of aerosol optical depth and its impact on global radiation over China[J].Atmospheric Environment,2017,150:106?115.

[3] IPCC.Climate change impacts, adaptation, and vulnerability(Part A):global and sectoral aspects.Contribution of working group II to the Fifth assessment report of the intergovernmentalpanel on climate change[M].Cambridge:Cambridge University Press,2014.

[4] Subba T,Gogoi M M,Pathak B,et al.Assessment of 1D and 3D model simulated radiation flux based on surface measurements and estimation of aerosol forcing and their climatological aspects[J].Atmospheric Research,2018,204:110?127.

[5] 趙俊芳,李寧,候英雨,等.基于APSIM模型評估北方八省春玉米生產(chǎn)對氣候變化的響應(yīng)[J].中國農(nóng)業(yè)氣象,2018,39(2):108?118.

Zhao J F,Li N,Hou Y Y,et al.Evaluation of response of spring maize production to climate change in the eight provinces of Northern China based on APSIM model[J].Chinese Journal of Agrometeorology,2018,39(2):108?118.(in Chinese)

[6] 任義芳,王春乙,趙艷霞.氣溶膠輻射效應(yīng)對作物及生態(tài)系統(tǒng)的影響研究[J].中國農(nóng)業(yè)氣象,2010,31(4):533?540.

Ren Y F,Wang C Y,Zhao Y X.Review on impact of atmospheric aerosol radiation effect on crops and ecological system[J].Chinese Journal of Agrometeorology,2010,31(4):533?540.(in Chinese)

[7] 周杰,陳健,張海龍,等.基于MODIS數(shù)據(jù)的長三角地區(qū)氣溶膠直接輻射效應(yīng)研[J].生態(tài)與農(nóng)村環(huán)境學(xué)報,2018,34(1):54?63．

Zhou J,Chen J,Zhang H L,et al.Study of aerosol direct radiation effects in the Yangtze river delta based on MODIS data[J].Journal of Ecology and Rural Environment,2018,34(1):54?63.(in Chinese)

[8] 鄧學(xué)良,何冬燕,潘德爐,等.衛(wèi)星遙感中國海域氣溶膠直接輻射強迫[J].氣象學(xué)報,2010,68(5):666?679.

Deng X L,He D Y,Pan D L,et al.Aerosol direct forcing estimated from satellite data over the China sea[J].Acta Meteorologica Sinica,2010,68(5):666?679.(in Chinese)

[9] 陳林,石廣玉,王標(biāo),等.基于衛(wèi)星觀測資料的氣溶膠直接輻射強迫研究[A].人類發(fā)展的永恒主題論文集[C].北京:第27屆中國氣象學(xué)會年會應(yīng)對氣候變化分會,2010:10?25.

Chen L,Shi G Y,Wang B,et al.Study on direct radiation forcing of aerosol based on satellite observation data[A].The eternal theme of human development[C].Beijing:The 27th Annual Meeting of the Chinese Meteorological Society to Deal With Climate Change,2010:10?25.(in Chinese)

[10] 張志薇,王宏斌,張鐳,等.中國地區(qū)3個AERONET站點氣溶膠直接輻射強迫分析[J].中國科學(xué)院大學(xué)學(xué)報,2014,31(3):297?305.

Zhang Z W,Wang H B,Zhang L,et al.Aerosol direct radiative forcing at three AERONET sites in China[J].Journal of University of Chinese Academy of Sciences,2014,31(3):297?305.(in Chinese)

[11] Holben B N,Eck T F,Slutsker I,et al.A federated instrument network and data archive for aerosol characterization[J].Remote Sensing of Environment,1998,66:1?16.

[12] Boiyo R,Kumar K R,Zhao T L.Statistical intercomparison and validation of multisensory aerosol optical depth retrievals over three AERONET sites in Kenya,East Africa[J].Atmospheric Research,2017,197:277?288.

[13] 余予,陳洪濱,夏祥鰲,等.青藏高原納木錯站地表反照率觀測與MODIS資料的對比分析[J].高原氣象,2010,29( 2):260?267.

Yu Y,Chen H B,Xia X A,et al.Comparison of surface albedo measurement with MODIS product at Namco Station of Tibetan Plateau[J].Plateau Meteorology,2010,29(2):260?267.(in Chinese)

[14] Guo J P,Zhang X Y,Wu Y R,et al.Spatio-temporal variation trends of satellite-based aerosol optical depth inChina during 1980?2008[J].Atmospheric Environment,2011,45(37):6802?6811.

[15] Xu H,Guo J P,Ceamanos X,et al.On the influence of the diurnal variations of aerosol content to estimate direct aerosol radiative forcing using MODIS data[J].Atmospheric Environment,2016,141:186?196.

[16] Dubovik O,Smirnov A,Holben B N,et al.Accuracy assessments of aerosol optical properties retrieved from Aerosol Robotic Network(AERONET) sun and sky radiance measurements[J].Journal of Geophysical Research:Atmospheres,2000,105(D8):9791?9806．

[17] Ge J M,Su J,Fu Q,et al.Dust aerosol forward scattering effects on ground-based aerosol optical depth retrievals[J].Journal of Quantitative Spectroscopy and Radiative Transfer,2011,112(2):310?319.

[18] Xia X A,Chen H B,Wang P C.Validation of MODIS aerosol retrievals and evaluation of potential cloud contamination in East Asia[J].Journal of Environmental Sciences,2004,16(5):832?837.

[19] King M D,Kaufman Y J,Menzel W P,et al.Remote sensing of cloud, aerosol, and water vapor properties from the Moderate Resolution Imaging Spectrometer(MODIS)[J].IEEE Transactions on Geoscience and Remote Sensing,1992,30(1):2?27.

[20] 陳愛軍,周嬋,卞林根,等.藏北高原GlobAlbedo地表反照率的精度分析[J].高原氣象,2016,35(4):887?894.

Chen A J,Zhou C,Bian L G,et al.Assessment on the accuracy of Glob Albedo over the Northern Qinghai-Xizang Plateau[J].2016,35(4):887?894.(in Chinese)

[21] Vermote E F,TanréD,Deuzé J L,et al.Second simulation of the satellite signal in the solar spectrum,6S:an overview[J].Geoscience Remote Sensing,1997,35:675?686.

[22] 董智慧,苗放,葉成名,等.基于6S模型成都平原氣溶膠光學(xué)厚度反演應(yīng)用[J].長江大學(xué)學(xué)報（自然科學(xué)版),2011,8(9):44?46.

Dong Z H,Miao F,Ye C M,et al.Retrieval of aerosol optical depth in Chengdu Plain based on 6S model[J].Journal of Yangtze University (Nat Sci Edit),2011,8(9):44?46.(in Chinese)

[23] Kotchenova S Y,Vermote E F,Matarrese R,et al.Validation of avector version of the 6S radiative transfer code for atmospheric correction ofsatellite data(Part I):path radiance[J].Applied Optics,2006,45:6762?6774.

[24] 毛節(jié)泰,李成才.氣溶膠輻射特性的觀測研究[J].氣象學(xué)報,2005,63(5):622?635.

Mao J T,Li C C.Observation study of aerosol radiative properties over China[J].Acta Meteorology Sinica,2005,63(5):622?635.(in Chinese)

[25] 顏鵬,李維亮,秦瑜.近年來大氣氣溶膠模式研究綜述[J].應(yīng)用氣象學(xué)報,2004,15(5):629?640.

Yan P,Li W L,Qin Y.A review of the current development of atmospheric aerosol models[J].Journal of Applied Meteorology Science,2004,15(5):629?640.(in Chinese)

[26] 胡子梅.基于衛(wèi)星遙感和地面監(jiān)測數(shù)據(jù)的上海市氣溶膠污染研究[D].上海:華東師范大學(xué),2014.

Hu Z M.Study on aerosol pollution in Shanghai based on satellite remote sensing and ground monitoring data[D].Shanghai:East China Normal University,2014.(in Chinese)

Estimating Surface Direct Radiation Effect of Aerosol Based on MODIS and AERONET Data

ZHAO Jun-fang1, XU Hui1, KONG Xiang-na1,2, GUO Jian-ping1, XU Jing-wen2, YU Wei-guo3

(1. State Key Laboratory of Severe Weather, Chinese Academy of Meteorological Sciences, Beijing 100081, China; 2. Resources College, Sichuan Agricultural University, Chengdu 611130; 3. Anhui Agrometerological Center, Hefei 232001)

Under climate change, the research on the direct radiation effect of aerosol in high polluted areas is of great significance for atmospheric environment monitoring, climate change assessment and agricultural production distribution in these areas. In this paper, the direct radiation effects of aerosol from January 2001 to December 2016 were quantitatively evaluated in Beijing, Xianghe and Taihu three typical global aerosol monitoring network (AERONET) sites of China based on the observation data of AERONET and the surface albedo data of MODIS. 6S (Second Simulation of a Satellite Signal in the Solar) radiative transmission model was used. The results showed: (1) the PM2.5mass concentration in three typical high pollution sites showed obvious seasonal variations, and the pollution was more serious in the autumn and in the winter than other season. In especial, the high values of PM2.5mass concentration were mainly concentrated in January, February, November and December. (2) Compared with no aerosol, the interannual change of surface direct radiation caused by aerosol was obvious in each site. Affected by aerosol, the interannual change in daily average surface direct radiation in Beijing, Xianghe and Taihu stations was obviously increasing from 2001 to 2016. In especial, the maximum daily direct radiation was found in the Xianghe station (621.14W·m?2), followed by the Taihu station (743.29W·m?2) and the Beijing station (488.14W·m?2). (3) Affected by aerosol, the surface direct radiation in each station was reduced. Moreover, the effects of aerosol on the surface direct radiation in different stations were different. During the time of 2001?2016, the average surface direct radiation in Beijing, Taihu and Xianghe stations decreased by 32.29%, 24.01% and 15.07%, respectively. Among them, the aerosol had the greatest impact on the direct surface radiation of Beijing Station, followed by Taihu station and Xianghe station. (4) Over the recent 15 years, the aerosol surface radiation effect in Beijing, Xianghe and Taihu stations all increased.

MODIS; AERONET; Aerosol; PM2.5; Direct radiation effect

2018?04?13

中國氣象科學(xué)研究院基本科研業(yè)務(wù)費專項（2017R001；2017Z004）；國家重點研發(fā)計劃（2017YFA0603004）

趙俊芳（1977?），女，博士，研究員，研究方向為氣候變化、農(nóng)業(yè)氣象和碳循環(huán)等。E-mail：zhaojfcams@163.com

10.3969/j.issn.1000-6362.2018.11.001

趙俊芳,徐慧,孔祥娜,等.基于MODIS和AERONET的氣溶膠地表直接輻射效應(yīng)評價[J].中國農(nóng)業(yè)氣象,2018,39(11):693-701