2. 浙江省环境监测中心, 杭州 310012
2. Environmental Monitoring Center of Zhejiang Province, Hangzhou 310012
能见度是表征城市空气质量优劣的重要指标, 而能见度主要取决于大气的消光作用.大气消光作用包括气体消光作用和颗粒物消光作用, 其中, 颗粒物的消光作用, 即大气颗粒物对太阳光的散射和吸收, 是影响能见度的主要因素(Watson, 2002).颗粒物的光学特性取决于颗粒物的粒径、形态、化学组成、混合状态及吸湿性质等多种物理化学性质(Meier et al., 2009; Ma et al., 2012; Liu et al., 2014; Cao et al., 2012).因此, 研究颗粒物的这些物理化学性质是定量研究太阳辐射强迫与颗粒物环境效应的关键.
在研究颗粒物的消光作用方面, 早期研究多集中于颗粒物质量浓度与消光系数之间的关系, 并且发现二者具有显著的正相关关系(Chan et al., 1999; Wang et al., 2006), 但由于颗粒物化学成分复杂且粒径差异较大, 用单一的颗粒物质量浓度并不能很好地表征消光系数的变化趋势.此外, 有研究采用Mie模型模拟颗粒物的混合状态(Cheng et al., 2008)以探究颗粒物的消光性质, 然而Mie模型作为理论模型也无法很好地定量表征颗粒物对能见度的影响.美国IMPROVE项目(Interagency Monitoring of Protected Visual Environments)利用多元回归分析方法, 构建了由多种化学物种质量浓度及其相应消光效率组成的消光系数计算公式(Sisler et al., 2000), 简称IMPROVE公式, 该公式在涉及颗粒物消光作用的相关研究中得到了广泛的应用(Wang et al., 2016; Zhou et al., 2016; Yu et al., 2016).值得注意的是, IMPROVE公式主要针对的是PM2.5的化学成分和光学性质, 而对于与大气消光影响更为密切的PM1并不适用.除此以外, 绝大多数研究都是将有机气溶胶作为一个整体计算其消光贡献(Wang et al., 2016; Zhou et al., 2016; Yu et al., 2016), 但事实上, 有机气溶胶物理化学性质十分复杂, 不同源排放的不同类型的有机气溶胶光学性质差异很大, 因此, 有必要探讨不同类型有机气溶胶对大气消光的影响.
本研究基于高时间分辨率的气溶胶观测仪器, 利用多元线性回归方法建立干气溶胶中化学组分与消光系数之间的关系, 特别是量化不同类型的有机气溶胶对颗粒物的消光贡献, 探讨各组分在不同情况下对消光系数的影响, 为我国的雾霾治理工作提供科学依据.
2 仪器与方法(Instruments and methods) 2.1 观测点位与观测时间本研究的观测地点位于浙江省金华市金东区环保大楼7楼实验室(29°6′6″N, 119°41′8″E), 采样头架设在楼顶距离地面1.5 m左右.金东区位于金华市东面, 浙江省中部.采样地点周边空旷, 无明显的工业污染源.观测时间为2011年10月31日—11月28日, 为期29 d, 中途由于电力供应问题, 导致两段较大的数据缺失.
2.2 观测仪器与数据应用EC9841氮氧化物分析仪在线监测二氧化氮, 时间分辨率设为1 min.气体中二氧化氮对太阳可见光波段的吸收最显著, 因此, 气体吸收系数bag(单位:10-4 m-1)由所测得的二氧化氮浓度(单位:ppm)按照公式bag=3.3[NO2]进行换算(姚婷婷等, 2010).由于大气中的气体浓度较为稳定, 因此, 其它气体分子造成的散射作用通常比较恒定, 一般取13 Mm-1(Cohen, 1975)作为大气气体散射系数bsg.
应用美国Aerodyne公司生产的高分辨飞行时间气溶胶质谱仪(High Resolution Time-of-Flight Aerosol Mass Spectrometer, AMS)在线监测亚微米颗粒物, 获得PM1中有机物、硫酸根、硝酸根、铵离子及氯离子等非难熔组分的质量浓度, 时间分辨率为5 min.AMS的工作原理及仪器操作维护详见文献(Canagaratna et al., 2007).应用美国Magee公司生产的七波段碳黑仪(Aethalometer, AE-31)在线监测880 nm处黑碳的质量浓度.需要注意的是, AE-31采用的切割头为PM2.5的切割头.AMS和AE-31进样口前端均设有除湿装置.
应用美国DMT公司生产的三波段光声黑碳光度计(Photo-Acoustic Soot Spectrometer, Three Wavelength, PASS-3)对颗粒物的吸收系数bap和散射系数bsp进行在线监测, 检测波长532 nm, 时间分辨率设为5 min. PASS-3进样口前端架设干燥管, 因而该仪器所测定的及后文所指的均为干气溶胶消光系数.
大气总消光系数bext由气体吸收系数bag、气体散射系数bsg、颗粒物吸收系数bap、颗粒物散射系数bsp四者加和得到.
3 结果与讨论(Results and discussion) 3.1 大气消光状况观测期间, 金华地区大气消光状况如图 1所示, 大气消光系数平均值为(427±262) Mm-1.在总消光系数中, 颗粒物散射平均占比84.9%, 颗粒物吸收平均占比9.9%, 气体散射平均占比3.0%, 气体吸收平均占比2.2%.颗粒物的总消光贡献达94.8%, 与杭州(Wang et al., 2016)、上海(Zhou et al., 2016)、南京(Yu et al., 2016)等地的研究结果较为一致, 说明颗粒物的消光作用是导致能见度下降的主要原因.若考虑相对湿度对颗粒物消光的影响, 则颗粒物实际消光在总消光系数中的比例将更为显著.
观测期间, PM1中各化学组分的时间序列与平均质量百分比如图 2所示.观测期间, PM1的质量浓度在11.0~264 μg·m-3范围内波动, 平均值为(78.5±40.4) μg·m-3.从图 2中可以看出, 各污染物的质量浓度变化幅度很大, 显示出明显的城市污染特征.其中, 2011年11月13—15日的污染过程比较明显, 各组分的浓度均有上升, 有机物、硫酸根、硝酸根的浓度抬升比较明显.在整个观测期间, 对平均质量浓度贡献最大的是有机物, 占到PM1总质量浓度的36.5%;其次是硫酸根, 占到22.6%;硝酸根、铵离子、黑碳、氯离子则分别占到15.8%、11.6%、9.7%、3.8%.需要注意的是, 本研究中AE-31采样时所用的切割头为PM2.5的切割头, 根据Lan等(2011)所测的BC粒径分布, 本研究PM1中BC质量浓度可能偏高约20%.
观测期间, 有机物是PM1中比重最大的组分, 因此, 有机气溶胶在大气消光作用中的影响不容忽视.为了了解不同来源的有机气溶胶(Organic Aerosol, OA)对大气消光的影响, 本研究利用正矩阵因子分析模型(Drewnick et al., 2005) (Positive Matrix Factorization, PMF)对观测期间AMS得到的有机气溶胶高分辨率质谱数据进行因子聚类分析, 具体方法见文献(Ulbrich et al., 2009).一共解析出4个因子:还原性气溶胶(Hydrocarbon-like OA, HOA)、半挥发氧化有机气溶胶(Semi-volatile Oxygenated OA, SV-OOA)、低挥发氧化有机气溶胶(Low-volatility Oxygenated OA, LV-OOA)、生物质燃烧排放的有机气溶胶(Biomass Burning OA, BBOA). HOA和BBOA属于一次有机气溶胶(Primary OA, POA), LV-OOA和SV-OOA属于二次有机气溶胶(Secondary OA, SOA) (Mohr et al., 2009; Canagaratna et al., 2004).
观测期间, 不同类别有机气溶胶的时间序列与平均质量百分比如图 3所示.LV-OOA、SV-OOA、HOA、BBOA的平均质量百分比分别为33.3%、32.9%、24.0%、9.8%. SOA在观测期间平均占到总有机气溶胶的66.2%, 占到PM1质量浓度的24.2%, 这反映了观测期间SOA是金华秋季大气亚微颗粒物中十分重要的污染成分, 对能见度的削弱可能有比较大的影响.
由于大气中亚微米颗粒物的粒径(0.1~1 μm)更接近太阳可见光波段的波长(0.3~0.8 μm), 因此, 与粒径更大的颗粒物相比, 亚微米颗粒物的消光作用更加显著.此外, 在之前的分析中可知, 颗粒物的消光系数占大气总消光系数的94.8%.基于这两点, 本研究运用PASS-3实测532 nm下颗粒物的吸收和散射系数与同时期HR-Tof-AMS和AE-31所获得的PM1化学成分进行多元回归统计讨论(所用软件为SPSS V18).本研究考虑到不同类型有机气溶胶的消光作用可能不相同, 因此, 在模型构造过程中, 有机气溶胶的源解析结果将代替有机物参与拟合.
在本研究中, 假设所有颗粒物皆处于外部混合状态, 且每种化学物种的消光作用是相对独立的(Lowenthal et al., 1995).另外, 由于矿物、粉尘、金属等对PM1质量浓度的贡献很少且消光效率较低(Watson, 2002; Tao et al., 2014), 这些因素在颗粒物消光拟合过程中将不予考虑.同时, 考虑到拟合公式的简洁性, 将不同老化程度的SV-OOA和LV-OOA合并为SOA, 并且拟合公式将不引入常数项.
3.3.1 颗粒物吸收系数的拟合通常研究认为, 黑碳是唯一对太阳光有全谱吸收的大气颗粒物.但有研究表明, 除黑碳外, 气溶胶中有颜色的有机物包括土壤腐殖质、类腐殖质物质、生物气溶胶等对太阳光都有吸收作用(Andreae et al., 2006).本课题组姚婷婷等(2010)在拟合深圳市冬季大气消光系数与细粒子化学组成之间的关系式时发现, 尽管有机物对光的吸收效率远低于黑碳的吸收效率, 但有机物的质量浓度却高出黑碳数倍, 结果表明, 观测期间黑碳对颗粒物吸收系数的平均贡献为73%, 而有机物的贡献则为21%.因此, 针对颗粒物吸收系数的拟合主要考虑黑碳和有机物对颗粒物吸收系数的贡献.
应用BC、SOA、HOA、BBOA一同参与颗粒物吸收系数的拟合, 得到关系式(1).该回归方程和回归系数均通过显著性检验(Sig.F=0.000, BC、SOA、BBOA Sig.T=0.000, HOA Sig.T=0.001<0.05), 判定系数R2=0.989, 标准化残差符合正态分布.
(1) |
根据吸收拟合公式(1)可以计算出BC、SOA、HOA和BBOA对吸收系数的平均贡献分别为89.9%、3.3%、1.0%、5.8%.由于HOA对吸收系数的贡献率仅为1.0%, 为了拟合公式的简洁, 舍去HOA, 应用BC、SOA和BBOA重新拟合, 得到新的关系式为:
(2) |
该回归方程和回归系数均通过显著性检验(Sig.F=0.000, BC、SOA、BBOA Sig.T=0.000), 判定系数R2=0.989, 标准化残差符合正态分布.
3.3.2 颗粒物散射系数的拟合大量研究数据表明, 颗粒物中的硫酸铵、硝酸铵、氯化铵及有机物等化学物种是散射太阳光的主要成分(Watson, 2002).此外, 硫酸盐、硝酸盐、铵盐等具有吸湿性的物质在相对湿度较高的情况下, 其散射作用会增强(Liu et al., 2014; 陶金花等, 2015; 崔蕾等, 2016).但由于本研究中颗粒物的散射系数是由PASS-3经过干燥管之后直接测得的干气溶胶散射系数, 其相对湿度控制在35%以下, 颗粒物的吸湿增长并不明显, 故在针对颗粒物散射的拟合时未将湿度影响考虑在内.
为了方便拟合, 假设颗粒物中的硫酸盐、硝酸盐及氯盐均以铵盐的形式存在, 因此, 硫酸铵、硝酸铵及氯化铵的质量浓度可以由SO42-、NO3-、Cl-按照以下公式进行换算:[(NH4)2SO4]=1.375[SO42-], [NH4NO3]=1.29[NO3-], [NH4Cl]=1.507[Cl-].
应用(NH4)2SO4、NH4NO3、NH4Cl、SOA、HOA、BBOA一同参与颗粒物散射系数的拟合, 得到的关系式见式(3).该回归方程和回归系数均通过显著性检验(Sig.F=0.000, Sig.T=0.000), 判定系数R2=0.991, 标准化残差符合正态分布.
(3) |
散射拟合公式(3)中HOA虽然通过显著性检验, 但系数为负, 显然不合理, 因此, 舍去HOA, 重新进行拟合, 得到关系式(4).该回归方程和回归系数均通过显著性检验(Sig.F=0.000, Sig.T=0.000), 判定系数R2=0.991, 标准化残差符合正态分布.
(4) |
将最终挑选出来的吸收拟合公式(2)和散射拟合公式(4)相加可以获得浙江金华秋季观测期间颗粒物总消光系数的拟合关系式:
(5) |
图 4展示了颗粒物总消光系数拟合值与实测值之间的比较, 判定系数R2=0.977, 说明该模型利用PM1化学组分数据可以很好地预测颗粒物的消光状况.从总拟合关系式(5)中可以看出, BBOA的消光效率最高, 为17.3 m2·g-1, 说明生物质燃烧排放对大气能见度的影响十分重要.硝酸铵和硫酸铵的消光效率次之.Hand等(2007)总结了48篇运用多元线性回归方法研究PM2.5组分与消光系数的报道, 并且分析换算了干气溶胶里的硫酸铵和硝酸铵平均消光效率:硫酸铵(2.8±0.5) m2·g-1, 硝酸铵(2.8±0.5) m2·g-1.需要注意的是, 如果用PM2.5中相同化学组分参与这样的拟合, 那么基于PM1与PM2.5质量浓度的比值, 相应拟合中各物种的消光效率将会变成式(5)中相应值的60%~70%(Zhang et al., 2013).因此, 换算之后, 本研究中硫酸铵和硝酸铵的消光效率与Hand等(2007)综合报道的结果相近.
根据拟合关系式(5)可以计算得出PM1各组分消光贡献随时间的变化情况和平均贡献比例, 结果如图 5所示.需要注意的是, 由于AE-31的采样头为PM2.5的切割头, 本研究PM1中BC质量浓度可能偏高约20%, 但因消光贡献是由物种的质量浓度和拟合系数共同决定的, 故BC的消光贡献并不会偏高.在整个观测期间, 干气溶胶中消光贡献最大的是硫酸铵, 消光贡献率为11.0%~57.0%, 平均贡献率为35.1%;其次是硝酸铵, 消光贡献率为3.8%~43.1%, 平均值为26.7%;SOA的贡献率为4.8%~30.3%, 平均值为14.3%;BBOA的贡献率为0~52.7%, 平均值为11.2%;BC的贡献率为3.5%~36.1%, 平均值为8.7%;氯化铵的贡献率为0.1%~24.4%, 平均值为4.0%.SOA和BBOA的总消光贡献占比达到25.5%, 说明有机物的二次转化与生物质燃烧排放是造成浙江金华地区能见度下降的重要因素.从图 5中可以发现, BBOA在特定时段具有最大贡献(11月27日), 是导致低能见度的首要因子.
在本研究中, 硫酸铵的消光贡献最高; 类似的研究发现, 深圳(姚婷婷等, 2010)是有机物的消光贡献最高; 西安(王浥尘等, 2014)则是硝酸铵; 北京(Zhang et al., 2014; Han et al., 2015; Xu et al., 2016; Wang et al., 2015)则是有机物或者硝酸铵.可见不同区域需要采用不同的污染控制策略来更有效地改善大气能见度.
4 结论(Conclusions)1) 观测期间, 浙江金华的大气消光系数为(427±262) Mm-1, 其中, 颗粒物消光系数占总消光系数的94.8%, 说明颗粒物消光是引起金华地区能见度下降的主要原因.
2) 观测期间, 金华地区PM1的质量浓度范围为11.0~264 μg·m-3, 平均值为(78.5±40.4) μg·m-3.对PM1质量浓度贡献最大的是有机物, 占到PM1总质量浓度的36.5%;其次是硫酸盐, 占到22.6%;硝酸盐、铵盐、黑碳、氯盐分别占到15.8%、11.6%、9.7%、3.8%.同时利用PMF模型解析出PM1中4类有机气溶胶, HOA、BBOA、LV-OOA、SV-OOA对有机气溶胶的贡献分别为24.0%、9.8%、33.3%、32.9%.
3) 基于PM1化学组分和有机气溶胶的PMF源解析结果, 利用多元线性回归方法, 建立了拟合优度很高的干气溶胶消光系数与PM1化学组分之间的定量关系模型:
4) 根据模型计算结果, 干气溶胶中消光贡献最大的是硫酸铵, 贡献率为35.1%;其次是硝酸铵, 贡献率为26.7%;SOA、BBOA、BC及氯化铵的消光贡献率分别为14.3%、11.2%、8.7%、4.0%. BBOA在特定时段具有最大贡献, 是导致低能见度的首要因子.
Andreae M O, Gelencser A. 2006. Black carbon or brown carbon? The nature of light-absorbing carbonaceous aerosols[J]. Atmospheric Chemistry and Physics, 6(10): 3131–3148.
DOI:10.5194/acp-6-3131-2006
|
Canagaratna M R, Jayne J T, Ghertner D A, et al. 2004. Chase studies of particulate emissions from in-use New York City vehicles[J]. Aerosol Science and Technology, 38(6): 555–573.
DOI:10.1080/02786820490465504
|
Canagaratna M R, Jayne J T, Jimenez J L, et al. 2007. Chemical and microphysical characterization of ambient aerosols with the aerodyne aerosol mass spectrometer[J]. Mass Spectrometry Reviews, 26(2): 185–222.
DOI:10.1002/(ISSN)1098-2787
|
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.
DOI:10.1016/j.atmosenv.2012.05.036
|
Chan Y C, Simpson R W, Mctainsh G H, et al. 1999. Source apportionment of visibility degradation problems in Brisbane (Australia) using the multiple linear regression techniques[J]. Atmospheric Environment, 33(19): 3237–3250.
DOI:10.1016/S1352-2310(99)00091-6
|
Cheng Y F, Wiedensohler A, Eichler H, et al. 2008. Aerosol optical properties and related chemical apportionment at Xinken in Pearl River Delta of China[J]. Atmospheric Environment, 42(25): 6351–6372.
DOI:10.1016/j.atmosenv.2008.02.034
|
Cohen A. 1975. Horizontal visibility and the measurement of atmospheric optical depth of lidar[J]. Applied Optics, 14(12): 2878–2882.
DOI:10.1364/AO.14.002878
|
崔蕾, 倪长健, 孙欢欢, 等. 2016. 成都颗粒物吸湿增长特征及订正方法研究[J]. 环境科学学报, 2016, 36(11): 3938–3943.
|
Drewnick F, Hings S S, Decarlo P, et al. 2005. A new time-of-flight aerosol mass spectrometer (TOF-AMS)-Instrument description and first field deployment[J]. Aerosol Science and Technology, 39(7): 637–658.
DOI:10.1080/02786820500182040
|
Han T, Xu W, Chen C, et al. 2015. Chemical apportionment of aerosol optical properties during the Asia-Pacific Economic Cooperation summit in Beijing, China[J]. Journal of Geophysical Research:Atmospheres, 120(23): 12, 212–281, 295.
|
Hand J L, Malm W C. 2007. Review of aerosol mass scattering efficiencies from ground-based measurements since 1990[J]. Journal of Geophysical Research, 112(D16).
DOI:10.1029/2017JD008484
|
Lan Z, Chen D, Li X, et al. 2011. Modal characteristics of carbonaceous aerosol size distribution in an urban atmosphere of South China[J]. Atmospheric Research, 100(1): 51–60.
DOI:10.1016/j.atmosres.2010.12.022
|
Liu H J, Zhao C S, Nekat B, et al. 2014. Aerosol hygroscopicity derived from size-segregated chemical composition and its parameterization in the North China Plain[J]. Atmospheric Chemistry and Physics, 14(5): 2525–2539.
DOI:10.5194/acp-14-2525-2014
|
Lowenthal D H, Rodgers C F, Saxena P, et al. 1995. Sensitivity of estimated light extinction coefficients to model assumptions and measurement errors[J]. Atmospheric Environment, 29(7): 751–766.
DOI:10.1016/1352-2310(94)00340-Q
|
Ma N, Zhao C S, Mueller T, et al. 2012. A new method to determine the mixing state of light absorbing carbonaceous using the measured aerosol optical properties and number size distributions[J]. Atmospheric Chemistry and Physics, 12(5): 2381–2397.
DOI:10.5194/acp-12-2381-2012
|
Meier J, Wehner B, Massling A, et al. 2009. Hygroscopic growth of urban aerosol particles in Beijing (China) during wintertime:a comparison of three experimental methods[J]. Atmospheric Chemistry and Physics, 9(18): 6865–6880.
DOI:10.5194/acp-9-6865-2009
|
Mohr C, Huffman J A, Cubison M J, et al. 2009. Characterization of primary organic aerosol emissions from meat cooking, trash burning, and motor vehicles with high-resolution aerosol mass spectrometry and comparison with ambient and chamber observations[J]. Environmental Science & Technology, 43(7): 2443–2449.
|
Sisler J F, Malm W C. 2000. Interpretation of trends of PM2.5 and reconstructed visibility from the IMPROVE network[J]. Journal of the Air & Waste Management Association, 50(5): 775–789.
|
Tao J, Zhang L, Cao J, et al. 2014. Characterization and source apportionment of aerosol light extinction in Chengdu, southwest China[J]. Atmospheric Environment, 95: 552–562.
DOI:10.1016/j.atmosenv.2014.07.017
|
陶金花, 王子峰, 徐谦, 等. 2015. 北京地区颗粒物质量消光吸湿增长模型研究[J]. 遥感学报, 2015, 19(1): 12–24.
DOI:10.11834/jrs.20153311 |
Ulbrich I M, Canagaratna M R, Zhang Q, et al. 2009. Interpretation of organic components from Positive Matrix Factorization of aerosol mass spectrometric data[J]. Atmospheric Chemistry and Physics, 9(9): 2891–2918.
DOI:10.5194/acp-9-2891-2009
|
王浥尘, 曹军骥, 张宁宁, 等. 2014. 西安大气细颗粒(PM1)化学组成及其对能见度的影响[J]. 地球科学与环境学报, 2014, 36(3): 94–101.
|
Wang J, Zhang Y, Feng Y, et al. 2016. Characterization and source apportionment of aerosol light extinction with a coupled model of CMB-IMPROVE in Hangzhou, Yangtze River Delta of China[J]. Atmospheric Research, 178.
|
Wang J, Zhang Y, Shao M, et al.2006.Quantitative relationship between visibility and mass concentration of PM2.5 in Beijing[J].Journal of Environmental Sciences, 18(3):475-481
|
Wang Q, Sun Y, Jiang Q, et al. 2015. Chemical composition of aerosol particles and light extinction apportionment before and during the heating season in Beijing, China[J]. Journal of Geophysical Research:Atmospheres, 120(24): 12708–12722.
DOI:10.1002/2015JD023871
|
Watson J G. 2002. Visibility:Science and regulation[J]. Journal of the Air & Waste Management Association, 52(6): 628–713.
|
姚婷婷, 黄晓锋, 何凌燕, 等. 2010. 深圳市冬季大气消光性质与细粒子化学组成的高时间分辨率观测和统计关系研究[J]. 中国科学:化学, 2010, 40(8): 1163–1171.
|
Xu X, Zhao W, Zhang Q, et al. 2016. Optical properties of atmospheric fine particles near Beijing during the HOPE-J3A campaign[J]. Atmospheric Chemistry and Physics, 16(10): 6421–6439.
DOI:10.5194/acp-16-6421-2016
|
Yu X, Ma J, An J, et al. 2016. Impacts of meteorological condition and aerosol chemical compositions on visibility impairment in Nanjing, China[J]. Journal of Cleaner Production, 131: 112–120.
DOI:10.1016/j.jclepro.2016.05.067
|
Zhang J K, Wang L L, Wang Y H, et al. 2016. Submicron aerosols during the Beijing Asia–Pacific Economic Cooperation conference in 2014[J]. Atmospheric Environment, 124: 224–231.
DOI:10.1016/j.atmosenv.2015.06.049
|
Zhang Y, Sun J, Zhang X, et al. 2013. Seasonal characterization of components and size distributions for submicron aerosols in Beijing[J]. Science China Earth Sciences, 56(5): 890–900.
DOI:10.1007/s11430-012-4515-z
|
Zhou M, Qiao L, Zhu S, et al. 2016. Chemical characteristics of fine particles and their impact on visibility impairment in Shanghai based on a 1-year period observation[J]. Journal of Environmental Sciences, 48: 151–160.
DOI:10.1016/j.jes.2016.01.022
|