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

趙俊芳，徐 慧，孔祥娜,2，郭建平，余衛國

?

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

趙俊芳1，徐 慧1，孔祥娜1,2，郭建平1，余衛國3

（1．中國氣象科學研究院災害天氣國家重點實驗室，北京 100081；2．四川農業大學資源學院，成都 611130；3．安徽省農業氣象中心，合肥 232001）

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

MODIS；AERONET；氣溶膠；PM2.5；直接輻射效應

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

氣溶膠的輻射效應是指由氣溶膠所引起的出射輻射通量的變化，通常可分為直接輻射效應和間接輻射效應[7]。其中，氣溶膠直接輻射效應指大氣中氣溶膠的含量變化引起地－氣系統能量平衡的擾動，在數值上定義為氣溶膠含量變化時所產生的平均凈輻射變化[8]，也就是在有氣溶膠和無氣溶膠存在時輻射通量的變化。在氣溶膠的輻射效應評估方面，國內外近年來一直大力發展衛星遙感與地基光學遙感相結合的觀測技術與反演方法[7?10]。中高分辨率的成像光譜儀（MODIS、HIRDLS、GLI）、多角度的成像光譜儀（MISR）等星載探測器的一個重要應用方向正是探測全球（包括陸地）的氣溶膠光學特性。在地基遙感方面，近年來較重要的研究進展之一是建立了AERONET（AErosol RObotic NET work）太陽光度計全球探測網絡由此提取氣溶膠的光學特性，并用于對衛星遙感氣溶膠的檢驗[11?12]。Boiyo等[12]利用東非肯尼亞3個AERONET站點MODIS Terra、MISR和OMI數據反演了氣溶膠光學厚度，并對其反演精度進行了比較。張志薇等[10]借助平面平行輻射傳輸模式評估了蘭州大學半干旱氣候與環境觀測站SACOL、香河和太湖地區大氣層頂氣溶膠日平均直接輻射強迫分別為（?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年氣溶膠光學厚度變化及其對全國輻射的影響。陳林等[9]利用MODIS和MISＲ資料，以及RAD01輻射模式評估了中國陸地地區的氣溶膠直接輻射效應，結果表明陸地上氣溶膠對大氣有很強的加熱作用，其強吸收性導致地面輻射強迫為大氣層頂輻射強迫的3倍。

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

1 資料與方法

1.1 研究區概況

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

表1 研究所選AERONET站點信息

1.2 資料來源與處理

1.2.1 Aerosol Robotic Network（AERONET）數據

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

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

1.2.2 MODIS數據

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

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

1.2.3 地表反照率數據

地表反照率表征地球表面對太陽輻射的反射能力，是影響氣溶膠輻射效應的重要因素。地表反照率數據來源于MODIS的MCD43地表反照率產品。MCD43分500m、1km和0. 05°氣候網格等3種空間分辨率（分別編號A、B和C），提供BRDF模型參數、BRDF反演質量、反照率和BRDF校正后的天底反射率NBAR（Nadir BRDF-Adjusted Reflectance）等4種產品（依次編號1?4），其中，反照率產品提供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]。研究數據為空間分辨率0.05°、正弦投影的MCD43地表反照率產品。

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

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

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

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

FSFC= FSFC↓?FSFC↑ （2）

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

日均氣溶膠直接輻射效應為1d內的瞬時輻射效應的平均值[15]，即

1.4 氣溶膠PM2.5質量濃度

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

2 結果與分析

2.1 典型站點氣溶膠PM2.5質量濃度季節變化特征分析

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

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

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

2001?2016年，基于6S輻射傳輸模式模擬的中國AERONET三個典型高污染站點地表直接輻射的年際變化特征如圖2所示。由圖可見，無氣溶膠影響下，北京、香河和太湖3個站點2001?2016年地表日平均直接輻射年際變化不太明顯。總體來說，太湖站日平均直接輻射最大，北京站最小。其中，近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輻射傳輸模式模擬結果顯示，氣溶膠對3個典型高污染站點地表直接輻射的影響不同，對北京站的影響最大，香河站最小（圖3）。2001?2016年，北京、太湖和香河站平均地表直接輻射分別降低32.29%、24.01%和15.07%。氣溶膠的影響導致近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 典型站點由氣溶膠影響導致的日平均直接輻射降低率

2001?2016年，中國AERONET北京、香河和太湖3個典型站點氣溶膠地表直接輻射效應分別平均為?207.45、?80.94和?216.71W·m?2。地表直接輻射效應大的負值（絕對值）說明氣溶膠阻擋太陽光到達地表，對地表產生了顯著的降溫，即氣溶膠在地表有顯著的降溫作用[15]。近15a來，北京、香河和太湖3個站點氣溶膠地表直接輻射效應呈現增加趨勢，對地表的降溫作用逐漸減弱，且北京和太湖站氣溶膠對地表的降溫作用遠遠大于香河站。其中，北京站氣溶膠地表直接輻射效應從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 結論與討論

3.1 結論

（1）北京、香河和太湖3個高污染站點氣溶膠PM2.5質量濃度季節變化特征差異明顯，秋冬污染較為嚴重，其中高值主要集中在1月、2月、11月和12月。2001?2016年期間，北京站點各月PM2.5質量濃度平均為92.96μg·m?3，香河站各月PM2.5質量濃度平均為95.31μg·m?3，太湖站各月PM2.5質量濃度平均為54.50μg·m?3。各個站點秋冬污染較為嚴重，高值主要集中在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年地表日平均直接輻射年際變化不太明顯。總體來說，在無氣溶膠影響下太湖站日平均直接輻射最大（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個站點氣溶膠地表輻射效應均呈現增加趨勢，說明氣溶膠對地表的降溫作用逐漸減弱，且在北京和太湖站氣溶膠對地表的降溫作用遠遠大于香河站。其中，北京站氣溶膠地表直接輻射效應從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 討論

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

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

[1] 張小曳.中國大氣氣溶膠及其氣候效應的研究[J].地球科學進展,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模型評估北方八省春玉米生產對氣候變化的響應[J].中國農業氣象,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] 任義芳,王春乙,趙艷霞.氣溶膠輻射效應對作物及生態系統的影響研究[J].中國農業氣象,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數據的長三角地區氣溶膠直接輻射效應研[J].生態與農村環境學報,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] 鄧學良,何冬燕,潘德爐,等.衛星遙感中國海域氣溶膠直接輻射強迫[J].氣象學報,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] 陳林,石廣玉,王標,等.基于衛星觀測資料的氣溶膠直接輻射強迫研究[A].人類發展的永恒主題論文集[C].北京:第27屆中國氣象學會年會應對氣候變化分會,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] 張志薇,王宏斌,張鐳,等.中國地區3個AERONET站點氣溶膠直接輻射強迫分析[J].中國科學院大學學報,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模型成都平原氣溶膠光學厚度反演應用[J].長江大學學報（自然科學版),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] 毛節泰,李成才.氣溶膠輻射特性的觀測研究[J].氣象學報,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].應用氣象學報,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] 胡子梅.基于衛星遙感和地面監測數據的上海市氣溶膠污染研究[D].上海:華東師范大學,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

中國氣象科學研究院基本科研業務費專項（2017R001；2017Z004）；國家重點研發計劃（2017YFA0603004）

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

10.3969/j.issn.1000-6362.2018.11.001

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