山西省氣溶膠光學厚度時空變化特征及氣候效應分析

王雁，郭偉，閆世明，李瑩

山西省氣象科學研究所，山西 太原 030002

由溫室氣體排放引起的全球氣候變暖已經引起廣泛關注，伴隨著溫室氣體的排放，人類活動還產生了大量的氣溶膠粒子（Boucher et al.，2001）。這些氣溶膠粒子不僅能降低大氣能見度（Cao et al.，2012）、影響大氣環境（陳永橋等，2005）、危害人類健康（Marambio et al.，2010；郭新彪等，2013），還能夠吸收和散射太陽輻射，部分氣溶膠粒子能夠形成云凝結核，增加云的光學厚度和云層反射率，從而降低了到達地面的短波輻射，導致地面氣溫下降（Ridley et al.，2014；廖禮等，2015），因此關于氣溶膠氣候效應的研究一直是氣候預測的關鍵部分。近年來，關于氣溶膠氣候效應的研究多集中在硫酸鹽氣溶膠（Zhang et al.，2012；王東東等，2014）、沙塵氣溶膠（劉沖等，2015；宿興濤等，2016）和黑炭氣溶膠（蔡子穎等，2011；黃文彥等，2015）等方面，研究方法主要是利用中尺度模式或區域氣候模式模擬氣溶膠對氣候變化的影響（黃文彥等2015；劉沖等，2015；宿興濤等，2016）。吳蓬萍等（2011）利用發展的 RIEMS對硫酸鹽氣溶膠第一間接輻射及其氣候效應進行研究，指出東亞地區硫酸鹽氣溶膠的第一間接輻射強迫為負，總體上使地表溫度降低，夏秋季平均降溫幅度在重慶及周邊地區最大可達1.5 ℃。四川盆地由于其特殊的地理和氣候，氣溶膠光學厚度（AOD）為全國之冠，基于衛星遙感結果的多年平均值在0.7左右（李成才等，2003a；Guo et al.，2011），因此，部分學者將四川盆地 80年代中后期的降溫歸因于大氣氣溶膠對低層大氣的反向散射（張天宇等，2017）。除四川盆地外，中國華北、長江流域、華南等地也是AOD的高值區（李成才等，2003b）。然而，除四川盆地外，中國其他地區卻鮮有關于氣溶膠大量排放而引起地表氣溫下降的報道。近年來，華北地區的霧霾天氣廣受關注，過去以煤炭消費為主的能源結構造成了嚴重的大氣污染，而不同地區排放的氣溶膠粒子化學組成和混合狀態均存在顯著差異，因此其在氣候變化中的作用仍難以被定量描述（石廣玉等，2008）。本研究根據邱金桓等（2001）于 2001年發展的應用地面氣象能見度和水汽壓信息確定大氣柱氣溶膠光學厚度的參數化模式，反演山西省1981—2016年 AOD，建立長序列資料，分析山西省工業顯著發展時期AOD的變化特征及其在氣候變化中的作用，為氣溶膠氣候效應研究提供參考。

1 資料和方法

1.1 研究區概況及資料選取

山西省位于黃土高原東緣，境內地形復雜，高差懸殊，分布有山地、丘陵、盆地等，以煤炭為主的工業結構造成山西省較為嚴重的大氣污染（何維燦等，2016）。大氣中的氣溶膠具有高度可變性，且其分布受地形特征影響顯著，因此利用較為密集的氣象站點能更準確地反映復雜地形下的氣溶膠光學厚度分布特征。本研究收集整理了山西省 109個地面氣象觀測站（山西省全部氣象臺站）1981—2016年逐日的水汽壓、天氣現象、每天4次能見度觀測資料，逐日平均氣溫、最高氣溫和最低氣溫資料，以及大同、太原、侯馬3站的太陽輻射資料。所有站點的氣象資料完整度高（缺測率低于0.02%），均來源于山西省氣象信息中心。20世紀80年代以來是全球增溫最為迅速的時期（沙萬英等，2002；魏鳳英等，2003），且同期也是山西省經濟快速發展、氣溫顯著變化的階段，因此本研究時間序列選擇1981—2016年。

1.2 能見度資料篩選

利用能見度估算氣溶膠光學厚度，需要排除降水和霧等自然天氣現象對能見度的影響。參考秦世廣等（2010）的方法對逐日觀測資料作以下處理：（1）每日觀測資料中選取 14:00的能見度觀測值；（2）將能見度觀測時段有降水（雪）和霧的記錄排除；（3）剔除相對濕度高于90%的資料。對數據資料進行篩選，篩選后109個站點觀測資料使用率平均為91.8%。

1.3 AOD反演計算方法

邱金桓等（2001）在 2001年發展了氣溶膠光學厚度模式，該模式比較適合中國特點，表達式為：

基于該式可計算結果為波長λ=0.55 μm下的光學厚度。在這個模式中，氣溶膠光學厚度依賴于地面能見度和水汽壓、海拔高度、波長和Junge譜參數。式中，V為海平面能見度；H1=0.886+0.0222V（km）；H2=3.77 km；Junge普參數v*=-2；Z為氣象臺站海拔高度；f為訂正系數。應用氣溶膠粒子濃度隨高度的分布函數，在海拔高度Z（km）上的能見度 VZ（km）與訂正到海平面的能見度 V（km）存在如下關系：

訂正系數f在不同地區存在2種不同的模式，對除東北以外的中國其他區域，f如下式所示：

式中，Pw為地面水汽壓（hPa）。根據此方法估算得到了 1981—2016年各氣象站氣溶膠光學厚度長時間序列。

2 結果與分析

2.1 山西省氣溶膠光學厚度的空間分布特征

圖1A所示為山西省1981—2016年AOD多年平均值的空間分布。山西省AOD空間分布呈現南高北低、中部區域高兩側山地低的特征，總體上同羅宇翔（2012）的研究結果相似。山西省北部的大同AOD為0.14，中部太原盆地AOD為0.33，南部地區侯馬AOD為0.46，區域差異顯著。山西省AOD在0.10～0.46之間，最低點出現在五臺山地區（圖A右上角深藍色區域），最高點為南部的侯馬市，全省平均值為0.23。從圖1A還可以看出，即使在山西省南部地區，部分站點的AOD也低于0.2，表明AOD整體上雖然有較為顯著的分布特征，但局部地區仍存差異。

山西省 AOD的空間分布與山西地形特征密切相關。圖1B所示為山西省地形圖，山西省東西兩側以山地為主，由北向南在中間區域分布著多個盆地。對AOD與海拔高度進行相關性分析，結果顯示，二者相關系數為-0.839，通過了0.01顯著性檢驗，表明海拔對AOD有顯著影響，海拔高度下降，大氣厚度增加，氣溶膠光學厚度增加。秦世廣等（2010）的研究顯示，整體上中國AOD呈現從東南向西北逐漸降低的趨勢，這也基本吻合中國的地形高差特征。

表1所示為山西省不同水平AOD區域數量，總體以中低值為主。低值區（0.1, 0.2]主要分布在高原山地及北部地區，中值區（0.2, 0.3]主要分布在忻定盆地、太原盆地、長治盆地等地區，高值區主要分布在臨汾運城盆地，AOD在0.4以上的地區分別是永濟、臨汾和侯馬。除了地形因素外，經濟發展和人類活動也對AOD有重要影響，山西省太原、臨汾和運城經濟總量常年位居山西省前三位，經濟的快速發展造成的污染物的較高排放是這些地區AOD偏高的根本原因。

表1 山西省不同水平AOD分布特征Table 1 Number of stations under different AOD level in Shanxi province

圖1 山西省1981—2016年AOD均值空間分布圖（A）及山西省地形圖（B）Fig. 1 AOD spatial distribution (A) from 1981 to 2016 and topographic (B) in Shanxi Province

2.2 山西省氣溶膠估算結果與已有研究比較

邱金桓等（2001）計算的 1990年太原 AOD為0.39，本研究計算的太原1990年AOD為0.28，36年平均為 0.33，本研究結果偏低的主要原因是數據篩選后對部分天氣條件下能見度資料進行了剔除。根據秦世廣等（2010）的研究，山西省AOD值在全國屬于中等水平，絕大部分地區高于0.3，中南部地區高于0.4，局部地區高于0.5，與本研究結果相比，整體上偏高。原因主要有兩方面，一是選取的站點位置不同，站點分布、密度也不同；二是時間序列不同，本研究更加關注20世紀80年代以后的變化。鄭小波等（2011）利用MODIS資料計算2000—2009年山西省AOD結果為0.33，同本文結果較為接近。總體上，本研究反演的山西省AOD結果與前人的研究都較接近，數據稍低主要原因是由于山西省地形復雜，且所選擇的氣象站點較為密集。

2.3 山西省氣溶膠光學厚度的時空變化

2.3.1 山西省氣溶膠光學厚度的空間變化特征

圖2所示為山西省1981—2016年AOD均值空間變化圖，此處直接用氣象站點顯示山西省 AOD空間變化特征。由圖可知，AOD值呈上升趨勢的地區占多數，達到 79個。從空間位置看，各個盆地上升站點較為密集，近 36年來，太原盆地、臨汾盆地和長治盆地AOD上升幅度均超過0.08，運城盆地上升了0.04，忻定盆地及其東側也有AOD值上升區，但幅度較小，僅為0.02。相對而言，AOD值下降的區域較為分散，太原盆地東側（呂梁山山脈）和西側（太行山中部地區）均有小范圍的集中區域，平均下降幅度在0.02左右。

圖2 山西省1981—2016年AOD空間變化圖Fig. 2 Spatial variation of AOD from 1981 to 2016 in Shanxi Province

圖3 山西省1981—2016年AOD變化趨勢圖Fig. 3 Variation trend of AOD from 1981 to 2016 in Shanxi Province

圖4 山西省不同季節AOD變化趨勢Fig. 4 Variation characteristics of AOD in different seasons in Shanxi Province

2.3.2 山西省氣溶膠厚度的長期變化趨勢

圖3和表2所示為1981—2016年山西省AOD的變化趨勢及其在不同海拔高度下的變化特征。總體上，山西省AOD上升趨勢顯著。AOD最低值出現在1986年（0.19），最高出現在2016年（0.27），然而這種上升趨勢并不是完全呈線性，在 2008—2013年之間出現了一個明顯的低谷，平均值為0.23，這可能與2008年北京舉辦奧運會大力減排以及隨后爆發的世界金融危機造成山西整體經濟下滑有關，2013年之后又表現為跳躍性上升。

根據山西省地形分布特征，將氣象站點按海拔高度劃分為 3個水平（海拔低于 500 m的站點、500～1000 m的站點和大于1000 m的站點）進行分析。從圖3可知，不同海拔下AOD的變化特征存在顯著差異，3個海拔AOD均呈顯著上升趨勢（表2），表明 AOD的增加具有全省域和整體性特征；海拔低于500 m的地區不僅是AOD值最大的地區，也是AOD上升最快的地區，對全省AOD增加貢獻最大。

2.3.3 山西省不同季節氣溶膠變化趨勢

圖4和表3所示為山西省不同季節AOD變化趨勢。按氣象上季節劃分標準：冬季為12—翌年2月，春季為3—5月，夏季為6—8月，秋季為9—11月。山西省AOD冬季最高，春秋季總體相差不大，夏季AOD最小。山西省冬季由于采暖需要，污染物排放量大，且降水稀少，風速較小，多盆地的地形特征造成了污染物的集中，使得冬季 AOD濃度偏高；春季風速偏大，易出現揚塵天氣，因此春季AOD波動性較大；夏季雨水充沛，濕沉降作用明顯，污染物濃度總體較低。從時間特征看，冬季、夏季和秋季AOD都呈顯著上升趨勢，其中冬季AOD上升最顯著，對全年的AOD變化貢獻最大；春季AOD變化趨勢不明顯。

2.4 山西省AOD的氣候效應分析

2.4.1 山西省1981—2016氣候變化特征

圖5所示為1981年以來山西省年平均氣溫的變化趨勢和突變特征。結果顯示，山西省 1981年以來年平均氣溫存在兩個不同的變化階段（圖5A），1981—1998年呈顯著上升趨勢，年均增幅在0.06 ℃左右，1998年以后氣溫增幅顯著放緩。運用M-K檢驗進行突變分析，結果顯示，年平均氣溫在1998年左右發生突變（圖5B），氣溫在1998年達到最高后，上升趨勢不再顯著。有研究指出，相對于1979—1998年，全球地表溫度增溫速率在1999年度之后的 10年間明顯偏低，此即全球變暖趨緩（停滯）現象（Chen et al.，2014；蘇京志等，2016）；盡管氣溫上升趨勢顯著減緩，但 1999—2016年平均氣溫顯著高于 1981—1998年，即近年來地表氣溫一直維持在較高水平。對山西省不同季節升溫幅度進行分析，結果顯示（圖略），在顯著升溫階段，冬季升溫幅度最顯著，其次是秋季和夏季，春季氣溫變化最不明顯。

表2 不同海拔高度AOD特征及變化趨勢Table 2 Characteristics and variation trend of AOD in different altitude in Shanxi Province

表3 山西省不同季節AOD特征及變化趨勢Table 3 Characteristics and variation trend of AOD in different seasons in Shanxi Province

圖5 山西省氣溫年際變化特征Fig. 5 Characteristics of temperature variation in 1981—2016 in Shanxi Province

2.4.2 山西省AOD對氣候變化的影響

氣溶膠能減少到達地面的太陽輻射量，從而影響地面氣溫。對大同、太原和侯馬AOD和太陽總輻射進行相關性分析，結果顯示，AOD與大同市太陽總輻射量呈負相關關系，相關系數為-0.135，未達到顯著水平，與太原市和侯馬市呈顯著負相關關系，相關系數分別為-0.525和-0.441，這表明，氣溶膠能減少到達地面的太陽輻射量，對高AOD地區太陽輻射影響更加顯著。上世紀 50年代以來，中國存在一個以四川盆地為中心的變冷區，這個變冷中心可以用AOD增加來解釋。楊琨等（2008）、張天宇等（2017）利用四川盆地的氣溫進行合成分析后指出，氣溶膠的冷卻效應可以解釋80年代到90年代中期的氣溫偏低現象，其中冬季氣溫對于AOD的響應最明顯。氣溶膠的冷卻效應在四川盆地體現為地表氣溫的下降，在中國大部分地區氣溫顯著升高的背景下，氣溶膠的冷卻效應可能體現為對氣候變暖的減緩作用。山西省氣溫變化顯著主要發生在1981—1998年之間，全省109個氣象站點，排除因遷站而導致氣溫顯著變化的五臺山站，保留其余108個站點氣溫資料，然后對108站1981—1998年的年、季節平均氣溫及日最高氣溫的增溫趨勢進行逐一分析，計算得出這 18年各站氣溫的年均增溫速率，再與同期（1981—1998年）AOD平均值進行相關分析，結果見表4。

表4 山西省AOD同1981—1988年間氣溫增速的相關關系Table 4 The correlation between AOD and air temperature increasing rate in Shanxi

在山西省顯著增溫時期，AOD與年平均氣溫增溫速率呈負相關關系，但未達到顯著水平，表明氣溶膠雖對氣候變暖有一定的減緩作用，但效果不明顯。從季節看，冬季平均氣溫增速和最高氣溫增速均與AOD呈顯著負相關關系，即AOD越大，氣溶膠對增溫的抑制作用越大，氣溫的增速越慢。

圖6 隨著AOD增大冬季平均氣溫增速、最高氣溫增速的變化特征Fig. 6 Change characteristics of average air temperature increase rate and maximum air temperature increase rate with the increase of AOD in winter

圖6 所示為1981—1998年山西省冬季平均氣溫和最高氣溫增速隨 AOD的變化特征，結果顯示，隨著AOD的增加，其增溫速率顯著下降，即山西省AOD越大的地區，冬季氣溫增速越慢。山西省冬季 AOD每增加 0.1，平均氣溫增速減小0.019 ℃·a-1，最高氣溫增速減少 0.020 ℃·a-1。一方面，冬季 AOD為全年最高；另一方面，1981—1998年間山西省冬季取暖以燃煤為主，由于當時技術所限基本無脫硫設施，因此冬季大氣氣溶膠硫酸鹽含量也較高，從而減緩了氣候變暖的進程。羅宇翔等（2012）研究顯示，中國春季AOD高值區的面積最大，其次是夏季，面積最小的是冬季，因此中國其他地區氣溶膠對增溫的抑制作用是否也發生在冬季，這有待進一步的研究。氣溶膠的環境和氣候效應是把雙刃劍（張小曳，2007），減少氣溶膠的排放有利于環境和人類健康，同時也可能降低對大氣的冷卻效應，加速氣候變暖的進程。

2.4.3 山西省AOD對氣溫日較差的影響

氣溶膠與氣溫日較差關系密切，白天大氣氣溶膠向上反射太陽短波輻射，減少到達地面的太陽輻射，阻止地表氣溫升高；夜間，大氣氣溶膠向下反射地表長波輻射，減少地表熱量散失，阻止地表氣溫降低，因此氣溶膠的存在能使氣溫日較差變小。圖7所示為山西省氣溫日較差與AOD的變化特征，二者之間的相關系數見表5。結果表明，山西省AOD與氣溫日較差呈顯著負相關關系，隨著AOD的升高，氣溫日較差有下降趨勢。就季節來看，冬季和夏季的負相關關系達到顯著水平，其中冬季負相關系數的絕對值最大。冬季氣溶膠濃度高，AOD較大，對太陽輻射影響大；夏季太原輻射最強烈，氣溶膠對太陽輻射的反射作用也最強（廖禮等，2015）。Zhang et al.（2010）應用WRF-Chem模式模擬美國大陸地區氣溶膠的氣候效應，結果表明氣溶膠分別能夠減少1月和7月的太陽輻射多達9%和16%，表明氣溶膠的冷卻效應在冬夏可能更加顯著。氣溶膠在不同季節的氣候效應的差異除了與濃度相關外，與氣溶膠本身的化學組成和混合狀態也有密切關系，這些都會影響氣溶膠的總輻射和吸收特性（張小曳，2007）。

表5 山西省AOD與氣溫日較差相關關系Table 5 Correlation coefficient between AOD and temperature daily range

3 結論

山西省氣溶膠光學厚度的多年分布具有典型的地理特征，整體上呈現為由北向南增加、由中部盆地向兩側山地減小的特征，呈現較為清晰的盆地特征。AOD全省多年平均值為0.23，不同地區AOD年均值變化范圍為0.10～0.46。

山西省AOD在 1981—2016年間整體上升趨勢明顯，但并不是完全呈線性，在2008—2013年間有1個低谷。除春季外，冬季、夏季和秋季AOD均呈顯著上升趨勢，其中以冬季變化幅度最大。山西省AOD上升的地區有79個，占72%，主要分布在太原盆地、臨汾盆地、運城盆地及長治盆地等地區。AOD分布和變化與海拔密切相關，海拔越低，AOD越大，低海拔地區的 AOD上升對全省AOD增加的貢獻最大。

山西省氣溫在1981—2016年間存在1個顯著升溫過程和 1個升溫停滯過程，氣溶膠的氣候冷卻效應對山西省1981—1998年間的氣候變暖有一定的減緩作用，尤其對冬季增溫有顯著抑制作用，即山西省AOD越大的地區，冬季氣溫增速越慢，冬季AOD每增加0.1，平均氣溫增速減小0.019 ℃·a-1，最高氣溫增速減少 0.020 ℃·a-1。此外，山西省 AOD與氣溫日較差存在顯著負相關關系，隨著AOD的增加，山西省氣溫日較差有下降的趨勢；冬夏季負相關系數的絕對值大于春秋季節，負相關關系更加顯著。

圖7 山西省年、季日較差和AOD變化特征Fig. 7 Variation of year and seasonal Temperature daily range and AOD over Shanxi Province from 1981 to 2016

BOUCHER O, HAYWOOD J. 2001. On summing the components of radiative forcing of climate change [J]. Climate Dynamics, 18(3):297-302.

CAO J, WANG Q, CHOW J C, et al. 2012. Impacts of aerosol compositions on visibility impairment in Xi'an, China [J]. Atmospheric Environment,59: 559-566.

CHEN X, TUNG K K. 2014. Varying planetary heat sink led to global-warming slowdown and acceleration [J]. Science, 345(6199):897-903.

GUO J P, ZHANG X Y, WU Y R, et al. 2011. Spatio-temporal variation trends of satellite-based aerosol optical depth in China during 1980—2008 [J]. Atmospheric Environment, 45(37): 6802-6811.

LIAO L, LOU S J, FU Y, et al. 2015. Radiative forcing of aerosols and its impact on surface air temperature on the synoptic scale in eastern China [J]. Chinese Journal of Atmospheric Sciences (in Chinese),39(1): 68-82.

MARAMBIO-JONES C, HOEK E M V. 2010. A review of the antibacterial effects of silver nanomaterials and potential implications for human health and the environment [J]. Journal of Nanoparticle Research,12(5): 1531-1551.

RIDLEY D A, SOLOMON S, BARNES J E, et al. 2014. Total volcanic stratospheric aerosol optical depths and implications for global climate change [J]. Geophysical Research Letters, 41(22): 7763-7769.

ZHANG H, SHEN Z, WEI X, et al. 2012. Comparison of optical properties of nitrate and sulfate aerosol and the direct radiative forcing due to nitrate in China [J]. Atmospheric research, 113: 113-125.

ZHANG Y, WEN X Y, JANG C J. 2010. Simulating chemistry-aerosolcloud-radiation-climate feedbacks over the continental US using the online-coupled Weather Research Forecasting Model with chemistry(WRF/Chem) [J]. Atmospheric Environment, 44(29): 3568-3582.

蔡子穎, 韓素芹, 黃鶴, 等. 2011. 天津夏季黑碳氣溶膠及其吸收特性的觀測研究[J]. 中國環境科學, 31(5): 719-723.

陳永橋, 張逸, 張曉山. 2005. 北京城鄉結合部氣溶膠中水溶性離子粒徑分布和季節變化[J]. 生態學報, 25(12): 3231-3236.

郭新彪, 魏紅英. 2013. 大氣 PM2.5對健康影響的研究進展[J]. 科學通報, 58(13): 1171-1177.

何維燦, 趙尚民, 程維明. 2016. 山西省不同地貌形態類型區土地覆被變化的GIS分析[J]. 地球信息科學學報, 18 (2): 210-219.

黃文彥, 沈新勇, 王勇, 等. 2015. 亞洲地區碳氣溶膠的時空特征及其直接氣候效應[J]. 大氣科學學報, 38(4): 448-457.

李成才, 毛節泰, 劉啟漢, 等. 2003b. 利用MODIS研究中國東部地區氣溶膠光學厚度的分布和季節變化特征[J]. 科學通報, 48(19):2094-2100.

李成才, 毛節泰, 劉啟漢. 2003a. 利用MODIS遙感資料分析四川盆地氣溶膠光學厚度時空分布特征[J]. 應用氣象學報, 14(1): 1-7.

廖禮, 漏嗣佳, 符瑜, 等. 2015. 中國東部氣溶膠在天氣尺度上的輻射強迫和對地面氣溫的影響[J]. 大氣科學, 39(1): 68-82.

劉沖, 趙天良, 熊潔, 等. 2015. 1991—2010年全球沙塵氣溶膠排放量氣候特征及其大氣環流影響因子[J]. 中國沙漠, 35(4): 959-970.

羅宇翔, 陳娟, 鄭小波. 2012. 近10年中國大陸MODlS遙感氣溶膠光學厚度特征[J]. 生態環境學報, 21(5): 876-883.

馬井會, 張華, 鄭有飛, 等. 2007. 沙塵氣溶膠光學厚度的全球分布及分析[J]. 氣候與環境研究, 12(2): 156-164.

秦世廣, 石廣玉, 陳林, 等. 2010. 利用地面水平能見度估算并分析中國地區氣溶膠光學厚度長期變化特征[J]. 大氣科學, 34(2): 449-456.

邱金桓, 林耀榮. 2001. 關于中國大氣氣溶膠光學厚度的一個參數化模式[J]. 氣象學報, 59(3): 368-372.

沙萬英, 邵雪梅, 黃玫. 2002. 20世紀80年代以來中國的氣候變暖及其對自然區域界線的影響[J]. 中國科學D輯, 32(4): 317-326.

石廣玉, 王標, 張華, 等. 2008. 大氣氣溶膠的輻射與氣候效應[J]. 大氣科學, 32(4): 826-840.

蘇京志, 溫敏, 丁一匯, 等. 2016. 全球變暖趨緩研究進展[J]. 大氣科學,40(6): 1143-1153.

王東東, 朱彬, 江志紅, 等. 2014. 硫酸鹽氣溶膠直接輻射效應對東亞副熱帶季風進程的影響[J]. 大氣科學, 38(5): 897-908.

魏鳳英, 曹鴻興, 王麗. 2003. 20世紀80～90年代我國氣候增暖進程的統計事實[J]. 應用氣象學報, 14(1): 79-86.

吳蓬萍, 韓志偉. 2011. 東亞地區硫酸鹽氣溶膠間接輻射和氣候效應的數值模擬研究[J]. 大氣科學, 35(3): 547-559.

宿興濤, 王宏, 許麗人, 等. 2016. 沙塵氣溶膠直接氣候效應對東亞冬季風影響的模擬研究[J]. 大氣科學, 40(3): 551-562.

楊琨, 孫照渤, 倪東鴻. 2008. 1999—2003年我國氣溶膠光學厚度的變化特征[J]. 南京氣象學院學報, 31(1): 92-96.

張天宇, 王勇, 程炳巖. 2017. 1961—2013年四川盆地氣溶膠光學厚度的長期變化及與氣溫的聯系[J]. 環境科學學報, 37(3): 793-802.

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

趙秀娟, 陳長和, 袁鐵, 等. 2005. 蘭州冬季大氣氣溶膠光學厚度及其與能見度的關系[J]. 高原氣象, 24(4): 617-622.

鄭小波, 周成霞, 羅宇翔, 等. 2011. 中國各省區近 10年遙感氣溶膠光學厚度和變化[J]. 生態環境學報, 20(4): 595-599.