刘随心,曹军骥,赵竹子,,沈振兴
(1.中国科学院地球环境研究所 中国科学院气溶胶化学与物理重点实验室,西安 710075;2.中国科学院地球环境研究所 黄土与第四纪地质国家重点实验室,西安 710075;3.中国科学院大学, 北京 100049;4.西安交通大学 环境科学与工程系, 西安 710049)
西安夏季碳气溶胶细粒子时间变化特征分析
刘随心1,2,曹军骥1,2,赵竹子1,2,3,沈振兴1,4
(1.中国科学院地球环境研究所 中国科学院气溶胶化学与物理重点实验室,西安 710075;2.中国科学院地球环境研究所 黄土与第四纪地质国家重点实验室,西安 710075;3.中国科学院大学, 北京 100049;4.西安交通大学 环境科学与工程系, 西安 710049)
为了探讨西安高分辨率的大气碳气溶胶细粒子的变化特征,在2008年夏季8月4日—7日太阳辐射强烈、高温晴朗天气进行每2小时的连续PM2.5样品采集,通过热光碳分析仪和TOC分析仪分析获得OC、EC、WSOC浓度变化序列,其日序列显示出机动车尾气排放为主导的昼夜变化特征。观测期OC/EC比值平均水平为6.41,应用EC示踪法估算出SOC每2小时平均浓度为11.64 μg·m−3,占OC浓度的63.37%,与WSOC占OC浓度的60.22%相当,两者比值为1.01,说明西安夏季高温和强烈的太阳辐射利于光化学反应的发生形成二次有机气溶胶,而且两者大部分为同一碳组分。分析SOC与O3以及OC/EC比值序列变化,发现以O3浓度和OC/EC比值来判断西安大气二次有机气溶胶的存在优于单纯应用OC/EC比值。
OC;EC;WSOC;SOC
当前大气中悬浮的气溶胶粒子成为城市环境主要污染物之一,尤其是细粒子PM2.5(空气动力学粒径小于2.5 μm的悬浮颗粒物),它不仅对人体健康、能见度具有重要危害,而且对全球和区域气候与环境系统产生显著影响(IPCC,2001;Jacobson,2001;Menon,2004) 。碳气溶胶是大气PM2.5的主要组分,早在上世纪末Chow et al(1993) 和Rogge et al(1993)认为在高污染城市碳组分已占PM2.520%~50%;碳气溶胶一般分为有机碳(Organic carbon,OC)和无机碳(Elemental carbon,EC)。OC包括一次气溶胶有机碳(Primary OC,POC)和二次气溶胶有机碳(Secondary OC,SOC),二次气溶胶有机碳主要是气态有机物经光化学转化产生的。EC来源于任何含碳物质不完全燃烧而排放出来的无定型碳质,它只存在于一次气溶胶中(Wolff,1981)。
伴随中国近几十年来的经济高速发展和耗能增加,大气环境污染已为城市发展的突出问题。考虑到大气颗粒物对于人体健康及能见度危害和区域及全球环境的影响,在城市中开展大气碳气溶胶的观测及研究显得尤为重要。目前一些主要城市已开展此方面研究,如He et al(2001)和Duan et al(2005)在北京、Feng et al(2006)在上海、Cao et al(2003,2004)在广州及珠三角城市区域(牛彧文等,2006)开展了大气碳气溶胶研究。在西安,自2002年以来,Cao et al(2005)对大气OC、EC进行了持续而系统的深入观测研究,获得了碳气溶胶的时空分布特征和演化规律,探讨了西安大气碳气溶胶污染来源和气象因子对其的影响,并率先在国内探索和发展运用估算二次碳气溶胶的方法。另外,西安是典型的北方城市,颗粒物污染严重,曾被联合国开发计划署选列为大气颗粒物污染研究的示范性城市(张小曳等,2001),对于其碳气溶胶的各种特征分析更具典型意义。但是西安碳气溶胶高分辨率的观测研究及来源分析相对较少,因此基于国内碳气溶胶研究情况和西安大气碳气溶胶已有的研究成果,本文应用中流量旋风式采样仪对西安大气气溶胶进行高分辨率连续观测,分析碳气溶胶变化特征和可能的来源及贡献。
1.1 采样地点和时间
西安(34°16'N,108°54'E)是典型的北方较大城市,人口约800万,地处关中平原中部,属暖温带半湿润季风气候,年平均降水量约604.2 mm,年平均温度约13.3℃。采样点位于城市西南高新产业园区域,PM2.5采样仪安置在中国科学院地球环境研究所楼顶,采样头距离地面10 m。
本文采样时间为2008年8月4日10:00—7日10:00,每2小时采集一个PM2.5样品。
1.2 样品采集和分析
高分辨率碳气溶胶样品采集应用113 L·min−1流量的旋风式采样仪(Bendix/Sensidyne 240,Florida,USA),每2小时采集一个样品。PM2.5采集均使用在900℃下预烧3小时(去除可能的含碳物质)的微纤维石英滤纸(Φ =9 mm,Whatman,England),预烧过的空白滤纸在采样前和采样后都保存在大约4℃恒温柜中,以防滤纸污染和采集后样品部分物质挥发。 所有样品的OC/EC分析采用中国科学院地球环境研究所的DRI Model 2001热光碳分析仪,该仪器采用IMPROVE(Interagency Monitoring of Protected Visual Environments)热光反射的实验方法(Cao et al,2003,2004,2005)。
水溶性有机碳(WSOC)的分析采用先提取再用TOC分析仪(Shimadzu TOC-5000A liquid analyzer)分析的方法(Yang et al,2003)。提取过程是:从采样滤膜上截取一定面积的滤膜放置于具塞离心管中,用10 mL移液管加入10 mL去离子水。再放置于恒温超声洗涤器中超声振荡60 min。并加入冰块,以防止超声器内水温增加,避免有机碳类物质的挥发损失。随后静置冰柜24 h,用 0.45 μm的针筒式微孔滤膜过滤器过滤掉试管中不溶的、悬浮颗粒,取4 mL溶液用于TOC分析。本次实验,每十张滤膜做一次重复样分析,重复样分析偏差均小于5%。
2.1 OC、EC、WSOC浓度水平及时间序列变化
2008年8月4日8:00—7日10:00每2 h PM2.5碳组分浓度分布范围和平均水平如表1所示,OC、EC浓度范围是11.96~32.99 μg·m−3和1.24~15.31 μg·m−3,其平均浓度分别为19.82 μg·m−3和3.80 μg·m−3;从碳组分浓度分布范围来看,EC最高和最低值相差达12.3倍,而OC最低、最高浓度相差仅2.8倍,由此可见,EC比OC的变率大,这可能是两者来源的相对变化存在较大差异。通过夏季高分辨率碳气溶胶观测发现,WSOC浓度平均占到OC 的60.22%,平均浓度为11.55 μg·m−3,说明有机碳大部分是可溶于水的组分。表1昼夜浓度显示,碳气溶胶各组分浓度在白天高于夜间,根据Cao et al(2005)、Shen et al(2008)、刘随心等(2009)对西安各个季节气溶胶排放源分析,在夏季机动车尾气排放是碳气溶胶主导性污染源,而机动车在白天排放远高于夜间,因此各个碳组分浓度在夏季白天高于夜间,并且白天浓度是夜间的1.2倍。对于水溶性碳而言,虽然WSOC白天浓度高于夜间,但可能由于白天有机物在大气中的累积或外来碳气溶胶传输以及夜间二次碳气溶胶的形成,使其夜间的百分比反而高于白天。
图1是加强观测期间36个PM2.5中OC、EC、WSOC和风速、风向的变化序列,从中可以看出,OC、EC、WSOC呈同步变化,而且每日具有相似变化特征,但是各碳组分浓度从4日开始逐渐降低。8月4日—7日,分别对应工作日的星期一至星期四,机动车流量变化相对稳定,对比风速和风向变化发现,5日8:00以后风速逐渐增大,而风向逐渐转为东北—东南风向,说明风速的增加增强了大气扩散能力,降低了气溶胶浓度水平,而东北—东南风向可能传输较少的污染物到采样点。从序列分布可以看到,外来大气输入对采样点气溶胶质量具有明显影响,如在5日4:00—6:00前后。图2为碳气溶胶组分小时变化序列,结合图1分析发现,日内碳组分浓度呈两个峰值分布特征,在早晨8:00前后碳组分浓度达到最高,形成较高的峰值;在下午18:00前后形成较小的峰值,两个峰值所对应的都是交通高峰期,在夜间碳浓度处于较低浓度水平,显示机动车对碳气溶胶的贡献白天高于夜间。但是WSOC/OC%却在每日午后和夜间出现明显的峰值,又说明大气中的有机碳还存在其他的来源贡献。
表1 夏季高分辨率观测PM2.5中OC、EC、WSOC浓度水平及WSOC/OC、OC/EC比值,风速Table 1 Average of high time-resolved OC, EC, WSOC, WSOC/OC, OC/EC ratios and wind speed in Xi'an during summer
图1 2008年8月4日至7日西安大气中高分辨率PM2.5中OC、EC、WSOC浓度、OC/EC比值,水溶性碳占OC百分比及风速,风向时间变化序列Fig.1 Time series of high time-resolved OC, EC, WSOC, OC/EC ratios, percentage of WSOC to OC and wind speed, wind direction
图2 西安碳气溶胶组分浓度及SOC百分比的日变化Fig.2 Diurnal variations of carbonaceous species, OC/EC ratios and percentage of SOC to OC
2.2 OC、EC、WSOC相关性及OC/EC比值变化特征
由于EC主要来源于燃烧过程,常被用来作为一次有机碳的示踪物Turpin et al(1990),因此应用OC与EC之间的关系和比值,可以区分碳气溶胶粒子的来源和判断二次气溶胶的生成。若各碳组分的相关性好,则表明其来自相同污染源,利用OC、EC、WSOC的相关性能在一定程度上对大气碳气溶胶的来源进行定性分析。图3为碳组分昼夜相关性分析,在加强采样期高分率PM2.5观测中OC、EC相关系数0.87,白天和夜间的相关系数分别为0.89和0.77,整个采样期和白天PM2.5OC、EC相关程度相当,而且高于夜间。OC、EC相关系数变化表明在较长观测期或整个夏季来看,OC、EC来源相对简单,主要反映了机动车尾气的贡献;而在较短观测尺度上,白天OC、EC来源相似,受机动车活动主导,相关系数高,但是夜间相对较小的相关系数则说明OC、EC的来源具有较大差异。夜间具有较低相关性则是因为在夜间机动车活动较少,向大气排放OC、尤其是来自机动车尾气的EC明显降低,而夏季太阳辐射强烈,温度高利于光化学反应形成二次有机碳,在当天午间或此前形成的二次有机碳在夜间累积和外来源有机物的传输,使得OC、EC来源发生差异,降低其相关程度, WSOC百分比水平在夜间的增高也说明这一现象。而OC与WSOC的相关性采样期和白昼都较好,分别为0.87、0.86和0.87,说明在夏季两者浓度变化受相似的污染源和化学机制影响。和OC、EC之间相关特征一样,EC与WSOC在整个采样期和白天相关程度高,相关系数分别为0.74、 0.77,而夜间明显降低为0.47,也指示了夜间因机动车活动较少,其他来源的有机物对大气中碳气溶胶贡献增加。
通常OC/EC比值超过2指示大气中存在二次有机碳气溶胶(Chow et al,1996;Cao et al,2004),而本次观测所有OC/EC比都超过该值,平均值达6.41,变化范围是2.15~12.30,夜间平均比高于白天(见表1),说明SOC在夏季是碳气溶胶不可忽略的来源。从图1可以看出,OC/EC最小、最大比值都出现在夜间,最小值出现在5日西风向主导时段4:00—6:00,OC、EC浓度最高,显示了外来源特征,为一次排放源主导;而最大比值所对应EC浓度相对稳定,OC浓度有较明显升高,而且主导为东风向,说明远源大气中OC较高或当天或此前大气中累积SOC。结合图1和图2,发现存在两种不同OC/EC峰值特征,一种是白天中午14:00左右,OC、EC浓度降低时,对应强烈太阳辐射和大气温度利于光化学反应生成二次有机碳气溶胶的峰值;另一种为在夜间大气中累积的一次、二次有机碳或外来源传输导致OC相对EC浓度升高形成峰值。
2.3 二次有机碳气溶胶(SOC)
相对于主要来自一次污染源排放的EC,OC除了一次排放外,有机气体经过光化学过程也可以形成有机碳,即二次有机碳。因目前还没有直接的分析方法能够区分气溶胶中的一次和二次有机物,本文根据 Turpin et al (1990)、Castro et al(1999)所应用OC/EC最小比值法来估算夏季高分辨观测大气中SOC并讨论其变化特征,其计算方法如下:
公式(1)中,OCsec是指估算的SOC浓度,而EC×(OC/EC)min则表示估算的一次有机碳,即POC;OCtot代表测量有机碳;(OC/EC)min表示观测到的最小OC/EC比值。
图3 碳气溶胶组分白天和夜间相关性Fig.3 Relationships between carbonaceous species in daytime and nighttime
观测期内(OC/EC)min为 2.15,出现在8月5日4:00,尽管在西风向下可能带入外来源传输到采样点增加OC、EC的浓度水平,但是选择该值估算SOC和前人选择OC/EC>2作为SOC生成的判断相一致,如Chow et al(1993,1996)和Cao et al(2004)夏天在珠三角的研究;另外因为夏季西安以机动车排放为主导,该比值接近机动车排放源OC/EC比值1.1(Chow et al,1996)。表2为通过最小值方法计算的SOC、POC浓度、百分比水平和平均的O3浓度,温度;图4是对应的浓度水平、百分比等的变化序列对比。高分辨率观测估算结果表明,SOC浓度在0~18.35 μg·m−3,平均浓度是11.64 μg·m−3,占OC最高达82.49%,平均百分比为63.37%;对应的POC浓度在2.67~32.92 μg·m−3,平均浓度8.18 μg·m−3,小于SOC水平,可见在西安SOC夏季对OC贡献显著。与已有的SOC研究相比,西安SOC浓度水平和百分比较高,如Cao et al(2004,2007)估算夏季中国北方城市(其中包括西安)、珠三角区域夏季城市平均SOC及其百分比分别为6.4 μg·m−3,40%;6.8 μg·m−3,42.6%;牛彧文等(2006)在深圳夏季获得的SOC及百分比分别为4.0 μg·m−3、56%。在西安白天SOC浓度水平稍高于夜间,因为在白天,太阳辐射强烈、温度高,易发生光化学反应,利于SOC形成,增加大气OC浓度;而夜间因为排放源减少,虽然有外来源传输和白天SOC在大气中的累积,但是相对于白天机动车排放相对较小。而SOC占OC百分比是夜间高于白天,由于夜间机动车活动减少,EC浓度也减少,相对增加了SOC浓度,这一变化与OC/EC比值昼夜间差异一致,与WSOC相比,SOC与其浓度在白天夜间分布一致,占OC百分比的昼夜变化相同。
图4 夏季高分辨碳组分浓度及SOC、WSOC百分比、OC/EC比值、O3和温度序列对比Fig.4 Concentrations of OC, EC, SOC, POC, WSOC, OC/EC ratios, percentage of SOC, WSOC to OC, and O3, temperature in summer
Saxena and Hildemann(1996)以及Huang et al(2006)研究发现SOC基本溶于水,而水溶性有机碳大部分为二次有机物,本研究WSOC/SOC比值在昼夜都接近1,说明西安在夏季SOC基本都溶于水。通常O3可以作为光化学反应和SOC形成的标志(Cabada et al,2004),当O3浓度升高相应地OC/EC比值也随之增大。本次高分辨率大气碳气溶胶观测在夏季晴朗天气状况开展,当温度最高时,O3浓度达到最大。但在采样期因温度高,太阳辐射活动强烈,早晨上班高峰期机动车增加大气中有机挥发性气体浓度水平,光化学反应在日出后已经进行,因此从图4可以看出,SOC峰值相对超前于O3峰值;而O3在午间达到峰值时,OC/EC比值相应出现峰值,说明中午光化学反应强烈,因此在西安,SOC的形成用O3和OC/EC作为其判断标准。
表2 应用OC/EC最小值估算夏季高分辨率估算PM2.5中SOC、POC浓度水平、SOC 、POC占OC百分比、WSOC/SOC比值, 及O3、温度Table 2 Average of high time-resolved SOC、POC、SOC/OC 、POC/OC%、WSOC/OC %、WSOC/SOC, O3and temperature in Xi'an during summer
通过西安2008年8月4日—7日开展的高分辨率大气碳气溶胶观测,分析获得以下结论:
(1)高分辨率碳气溶胶观测组分OC、EC、WSOC浓度变化清晰地反映了以机动车排放为主导的大气污染特征,同时高分辨率的浓度变化显示了外来排放源和气象影响。在机动车活动高峰期,各组分浓度出现明显峰值;而白昼机动车多少和气象差异使得各个碳组分浓度夜间低于白天。
(2)观测期OC、EC、WSOC每2小时浓度变化基本一致,相关程度较高,具有相似的排放来源,其中WSOC占到OC浓度的60.22%,说明夏季大部分OC为可溶性。
(3)每2小时的OC/EC比值都高于2,根据已有研究,说明西安夏季光化学反应强烈;应用EC示踪法获得观测期每2小时SOC平均浓度为11.64 μg·m−3,SOC占OC浓度平均为63.37%,指示西安夏季二次有机物污染突出,并且在夜间二次有机污染也非常严重。WSOC、SOC占OC百分比近似,两者比值约为1,说明西安SOC基本溶于水或WSOC基本是二次有机物。
(4)O3和OC/EC比值相结合能够作为SOC形成的指标。
刘随心,曹军骥,安芷生. 2009. 西安大气细粒子质量浓度变化特征及其影响因素[J]. 过程工程学报, 9(S2): 231–236. [Liu S X, Cao J J, An Z S. 2009. Characterization of Ambient Fine Particles (PM2.5) Concentration and Its Influential Factors in Xi'an [J]. The Chinese Journal of Process Engineering, 9(S2): 231–236.]
牛彧文,何凌燕,胡 敏,等. 2006.深圳冬、夏季大气细粒子及其二次组分的污染特征[J]. 中国科学(B 辑: 化学), 36(2): 173 –180. [Niu Y W, He L Y, Hu M, et al. 2006. Pollution characteristics of atmospheric fine particles and their secondary components in the atmosphere of Shenzhen in summer and in winter [J]. Science in China Series B: Chemistry, 49(5): 466 – 474.]
张小曳,曹军骥,苏广路,等. 2001.西安市的颗粒物污染控制.中国城市空气污染控制[M].北京:中国科学技术出版社, 222–293. [Zhang X Y, Cao J J, Su G L, et al. 2001. Particulate matter pollutants control in xian city. Air pollutants control in chinese cities [M]. China Science and Technology Press, 222–293.]
Cabada J C, Pandis S N, SubramanianR, et al. 2004. Estimating the secondary organic aerosol contribution to PM2.5using the EC tracer method [J]. Aerosol Science and Technology, 38(S1): 140 –155.
Cao J J, Lee S C, Chow J C, et al. 2007. Spatial and seasonal distributions of carbonaceous aerosols over China [J]. Journal of Geophysical Research, 112, D22S11,doi:10.1029/2006JD008205.
Cao J J, Lee S C, Ho K F, et al. 2003. Characteristics of carbonaceous aerosol in Pearl River Delta Region, China during 2001 winter period [J]. Atmospheric Environment, 37: 1451–1460.
Cao J J, Lee S C, Ho K F, et al. 2004. Spatial and seasonal variations of atmospheric organic carbon and elemental carbon in Pearl River Delta Region, China [J]. Atmospheric Environment , 38(27): 4447– 4456.
Cao J J, Wu F, Chow J C, et al. 2005. Characterization and source apportionment of atmospheric organic and elemental carbon during fall and winter of 2003 in Xi'an, China [J]. Atmospheric Chemistry and Physics, 5: 3127–3137.
Castro L M, Pio C A, Harrison R M, et al. 1999. Carbonaceous aerosol in urban and rural European atmospheres: estimation of secondary organic carbon concentrations [J]. Atmospheric Environment, 33: 2771–2781.
Chow J C, Watson J G, Lowenthal D H, et al. 1993. PM10and PM2.5compositions in California's San Joaquin Valley [J]. Aerosol Science and Technology, 18: 105–128.
Chow J C, Watson J G, Lu Z, et al. 1996. Descriptive analysis of PM2.5and PM10at regionally representative locations during SJVAQS/AUSPEX [J]. Atmospheric Environment, 30(12): 2079–2112.
Duan F K, He K B, Ma Y L, et al. 2005. Characteristics of carbonaceous aerosols in Beijing,China [J] . Chemosphere, 60: 355–364.
Feng J L, Hu M, Chan C K, et al. 2006. A comparative study of the organic matter in PM2.5from three Chinese megacities in three different climatic zones [J]. Atmospheric Environment, 40: 3983–3994.
He K B, Yang F M, Ma Y L, et al. 2001. Characteristics of PM2.5in Beijing, China [J]. Atmospheric Environment, 35: 4959 – 4970.
Huang X F, Yu J Z, He L Y, et al. 2006. Water-soluble organic carbon and oxalate in aerosols at a coastal urban site in China: Size distribution characteristics, sources, and formation mechanisms [J]. Journal of Geophysical Research, 111, D22212, doi:10.1029/2006JD007408.
IPCC (Intergovernmental Panel on Climate Change). 2001. Climate Change 2001: The Scientific Basis. 5. Aerosols, Their Direct and Indirect Effects [R]. London: Cambridge University Press, 289–348.
Jacobson M Z. 2001. Strong radiative heating due to the mixing state of black carbon in atmospheric aerosols [J]. Nature, 409: 695 – 697.
Menon S. 2004. Current uncertainties in assessing aerosol effects on climate [J]. Annul Review of Environment and Resources, 29: 1–30.
Rogge W F, Mazurek M A, Hildemann L M, et al.1993. Quantification of urban organic aerosols at a molecular level: identif cation, abundance and seasonal variation [J]. Atmospheric Environment, 27(8): 1309 –1330.
Saxena P, Hildemann L M. 1996. Water-soluble organics in atmospheric particles:a critical review of the literature and application of thermodynamics to identify candidate compounds [J]. Journal of Atmospheric Chemistry, 24: 57–109.
Shen Z X, Richard A, Cao J J, et al. 2008. Seasonal variations and evidence for the effectiveness of pollution controls on water-soluble inorganic species in total suspended particulates and fine particulate matter from Xi'an, China [J]. Air and Waste Management Association, 58: 1560 –1570.
Turpin B J, Cary R A, Huntzicker J J. 1990. An in-situ, timeresolved analyzed for aerosol organic and elemental carbon [J]. Aerosol Science Technology, 12: 161–171.
Turpin B J, Huntzicker J J. 1995. Identification of secondary aerosol episodes and quantitation of primary and secondary organic aerosol concentrations during SCAQS [J]. Atmospheric Environment, 29: 3527–3544.
Wolff G T. 1981. Particulate elemental carbon in the atmosphere [J]. Journal of Air Pollution Control Association, 31: 935–938.
Yang H, Li Q, Yu J Z. 2003. Comparison of two methods for the determination of water-soluble organic carbon in atmospheric particles [J]. Atmospheric Environment, 37: 865–870.
Temporal variations of carbonaceous aerosols in PM2.5during summer in Xi'an, China
LIU Sui-xin1,2, CAO Jun-ji1,2, ZHAO Zhu-zi1,2,3, SHEN Zhen-xing1,4
(1. Key Lab of Aerosol Chemistry & Physics, Institute of Earth Environment, Chinese Academy of Sciences, Xi'an 710075, China; 2. State Key Laboratory of Loess and Quaternary Geology, Institute of Earth Environment, Chinese Academy of Sciences, Xi'an 710075, China; 3. University of Chinese Academy of Sciences, Beijing 100049, China; 4. Department of Environmental Science and Engineering, Xi'an Jiaotong University, Xi'an 710049, China)
2-hr time-resolved carbonaceous particulate organic (OC), elemental carbon (EC) and watersoluble organics (WSOC) data measured by thermal/optical carbon and TOC analyzer in Xi'an during sunny days from 4 to 7 Aug 2008 in summer were investigated the temporal patterns of carbonaceous aerosols in PM2.5, the result suggests that the motor vehicle source is dominant sources for the variation of carbonaceous aerosols in day-night variation. Mean 2-hr average OC/EC ration is 6.41, the secondary organic carbon (SOC) was estimated using an EC tracer method and the average concentration is 11.64 μg·m−3, SOC contributions 63.37% of measured OC approximating that of WSOC, and the WSOC/SOC equaled to 1.01, it indicates that the high temperature and solar radiation are favorable case to SOC formed by photochemical reactions. In this study, ozone and OC/EC ratios were used as the indicator of photochemical activity to identify SOC production in Xi'an during summer.
OC; EC; WSOC; SOC
P402
:A
:1674-9901(2014)05-0311-08
10.7515/JEE201405002
2014-08-14
陕西省科技统筹创新工程项目(2012KTZB03-01,2011KTCQ03-04)
曹军骥,E-mail: cao@loess.llqg.ac.cn