典型西南岩溶区工业城市降水酸化缓解原因探究——以柳州市为例

删除或更新信息，请邮件至freekaoyan#163.com(#换成@)

本站小编 Free考研考试/2021-12-31

朱海燕^1,2, 辛存林¹, 刘齐³, 吴亮君², 于奭²
1. 西北师范大学地理与环境科学学院, 兰州 730070;
2. 中国地质科学院岩溶地质研究所/自然资源部、广西岩溶动力学重点实验室, 桂林 541004;
3. 柳州市环境保护监测站, 柳州 545001
收稿日期: 2018-06-25; 修回日期: 2018-08-13; 录用日期: 2018-08-13
基金项目: 地质调查项目（No.DD20160305）；国家自然科学基金（No.41402324）
作者简介: 朱海燕(1994-), 女, E-mail:1447969040@qq.com
通讯作者（责任作者）: 于奭, E-mail:yushi@karst.ac.cn

摘要: 利用柳州市环境保护监测站2013-2017年大气降水监测数据，分析近年来柳州岩溶区工业城市降水酸化缓解原因，并对离子来源进行探究.结果表明：①统计柳州近5年的降水pH值发现研究区降水酸化问题有了明显的改善，同时具有季节性差异，春冬季降水酸化相对夏秋季较重.②通过分析降水中主要离子组成及变化特征、计算中和因子（NF）和相对酸度（FA）发现，降水中主要离子浓度出现减少趋势，酸化类型逐渐向硫-硝酸混合型过渡，且99.97%的酸性离子被NH₄⁺、Ca²⁺等离子所中和，所以得出酸性物质的减少和碱性离子的中和作用可能是降水酸化改善的重要因素.③采用Hysplit模型对研究区降水酸化较严重的2016年（共60场）的降水水汽来源进行后向轨迹追踪及聚类分析，显示降水离子组成与浓度不仅与当地污染源有关，还与来自大陆内部西北和西南远距离物质输送有关.上述结果指示由于喀斯特化学风化作用、柳州市严格执行环保规定、大环境改善导致远程传输污染减弱等原因使得柳州市降水酸化情况得到极大改善.此外值得注意的是柳州市酸雨类型逐渐向硫-硝酸混合型过渡，指示随着工业污染的控制得到减弱，日益增长的机动车尾气排放或许成为研究区下一步治理重点方向.
关键词:柳州pH离子特征降雨酸化中和作用离子来源后向轨迹模拟
Mitigation of precipitation acidity in typical industrial areas of southwest Karst areas: A case study in Liuzhou
ZHU Haiyan^1,2, XIN Cunlin¹, LIU Qi³, WU Liangjun², YU Shi²
1. College of Geography and Environmental Science, Northwest Normal University, Lanzhou 730070;
2. Key Laboratory of Karst Dynamics, Ministry of Land and Resources & Guangxi; Institute of Karst Geology, Chinese Academy of Geological Sciences, Guilin 541004;
3. Environmental Monitoring Station of Liuzhou City, Liuzhou 545001
Received 25 June 2018; received in revised from 13 August 2018; accepted 13 August 2018
Abstract: According to the monitoring data of atmospheric precipitation in Liuzhou from 2013 to 2017, the effects of acid precipitation and its mitigation in the typical industrial karst areas of Southwest China were analyzed to trace the source of ions for the history of acidic precipitation. Our findings revealed that, in the past five years, the pH values increased slowly and had the obvious seasonal variation (pH values were lower in winter and spring than those in summer and fall). In this respect, the main ionic composition and variation of precipitation, the neutralization factor (NF) and the fractional acidity (AF) indicated that the ionic concentration of precipitation was still sulfur type with a reducing trend to sulfuric-nitrous mixed type and the 99.97% acid ions were neutralized by NH₄⁺ and Ca²⁺ ions, suggesting that it might be the important factor to mitigating acid precipitation. Besides that, the Hysplit model was applied to trace and cluster the precipitation with 60 times in 2016. It was found that the compositions and concentrations of ions about Liuzhou precipitation were not only controlled by local pollution sources, but also were affected by the transmission distance of pollution sources from the northwest and southwest in the interior of continents. Talked together, our results revealed that the reduced acid precipitation in Liuzhou might be as a result of the combined action of karst chemical weathering, strict implementing environmental policy, and the reduction of long-term pollutant. It should be noted that the precipitation has a reducing trend to sulfuric-nitrous mixed type, suggesting that vehicle exhaust should be focused on the future.
Keywords: LiuzhoupHionic characteristicsprecipitation acidificationneutralizationionic sourcebackward trajectory simulation
1 引言(Introduction)降水酸化已被视为人类21世纪所面临的4个主要环境问题之一, 大气降水酸化及其降水离子组成一直是国内外研究的重点(Lei et al., 1998).欧美国家从1950s开始对酸雨进行研究, 在1970s建立了酸雨观测网, 随后东亚地区也陆续建立, 我国在1982年和1989年也建立健全了全国酸雨观测网(Galloway et al., 1976; 王文兴, 1994; Ding et al., 1997; 丁国安等, 2004).国内长江以南-云贵高原以东地区是降水酸化影响的主要区域, 包括浙江、上海、江西、福建等, 目前降水酸化问题相对于20世纪90年代已呈现出很大程度地改善.降水酸化形成机制复杂, 影响因素众多, 不仅与当地污染源有关, 也与碱性物质的中和作用及中、长尺度物质传输有关.Araguás(1998)、Krupa(2002)、Akkoyunlu等(2003)、Okuda等(2004)、Al-Khashman(2005)等国内外专家关于降水酸化的分析研究主要从3方面展开：①降水的酸性程度、离子化学组成浓度的时空分布和形成机理；②污染物的来源追踪和跨国输送；③降水酸化的清除及治理过程.Báez等(2007)采用数学统计方法对墨西哥降水离子来源进行研究, 表明SO₄^2-、NO₃^-和NH₄⁺主要是人为来源.Al-Khashman等(2005)发现约旦的大气降水中Ca²⁺浓度最高, SO₄^2-来自远距离跨国输送占很大比例.Salam等(2003)发现SO₄^2-和NO₃^-是达卡降水中的主要酸性离子, NO₃^-的主要来源是汽车尾气排放.王文兴(1994)发现南方大部分酸雨区形成的重要原因是大气污染物局部排放和中长距离传输的叠加.
柳州作为广西最大的工业基地和区域经济中心, 由于其工业、城市布局、地形及周边重酸雨区的影响等, 致使酸雨一度成为困扰柳州环境问题中的顽疾, 几乎“十雨九酸”, 由于其严峻的酸雨问题及特殊的岩溶地质背景, 引起了****们对该区降水的广泛关注, 梁宏温等(2002)对柳州市1996—2001年的降水进行化学特征研究, 发现SO₄^2-是主要酸性离子, 酸雨污染类型为硫酸型；郭雅思(2017)对柳州市酸雨化学特征及人为硫源对碳酸盐岩风化的影响进行了深入的研究；董蕙青等(2005)从气象学方面研究发现广西76%酸雨受高空槽、切变线和锋面天气系统的影响；吴起鑫等(2011)和周晓得等(2017)分别对西南喀斯特农村地区和西南酸雨区降水进行了相关的化学研究.这些研究成果为探究柳州以往酸雨的化学特征与影响因素奠定了良好的研究基础, 但近年来柳州出现降水酸度改善的趋势, 降水离子组成和浓度也发生显著变化, 降水特征逐渐向硫-硝混合型转变, 与近二十余年形成了鲜明的对比(程佳惠等, 2015).本文利用2013—2017年柳州市环境保护监测站降水监测数据, 结合国内外理论, 通过对降水离子的变化特征、中和能力及来源等进行分析, 探究柳州市降水酸化的变化特征及改善原因, 并以2016年为例进行降水气团后向轨迹模拟追踪, 讨论周边区域对该地区降水的影响.通过以上研究手段旨在为进一步治理降雨酸化问题提供科学依据.
2 采样和分析方法(Sample collection and analysis methods)2.1 研究区概况柳州市(23°54′N~26°03′N, 108°32′E~110°28′E)位于广西壮族自治区中部偏东北地区(图 1), 处于柳江中游, 湘桂走廊的西南部.属于典型的亚热带季风气候, 常年被季风环流影响, 夏季高温多雨, 冬季干燥少雨, 年均气温在20.4 ℃左右, 年均降水量在1400~1500 mm之间.地区主要发育泥盆系和石炭系的灰岩和白云岩, 碳酸盐岩分布十分广泛, 岩溶地貌极为发育.研究区雨水采样点设置在柳州市环境保护监测站顶楼, 坐标为109°25′21.48″E, 24°22′11.22″N, 属于城区取样点；采样点周围较空旷, 半径50 m范围内无高大建筑物或树木.
图 1(Fig. 1)

图 1 柳州研究区及数据采集点位置示意 Fig. 1Research area and data point location of Liuzhou

2.2 数据来源样品的采集与保存按照国家标准GB 13580.2—1992《大气降水样品的采集与保存》规范执行.降雨大于5 mm进行采样分析(每月1~2次)；降水以及降雨量通过瑞士DIGITEL公司生产的DRA-12型自动湿沉降采样器收集和监测(仪器换新所以改动原文内容)，该仪器配置有标准雨量筒和降水自动采集器，自动监测降雨量和采集降水.降水的pH值通过德国WTW公司生产的水化学自动监测仪(pH/cond340i测试仪)进行现场测定，分辨率分别为0. 01pH单位；降雨的年、月平均pH值采用氢离子浓度雨量加权法；阳离子(K⁺、Na⁺、Ca²⁺、Mg²⁺和NH₄⁺)采用戴安ICS1500离子色谱仪分析测定, 阴离子(SO₄^2-、NO₃^-和Cl^-)用万通MIC离子色谱仪分析测定, 阴阳离子仪器测试精度为0.01 mg·L^-1, 分析误差小于5%.
2.3 后向型轨迹模拟与聚类分析运用美国大气实验室开发的HYSPLIT模型模拟柳州2016年(60场)降水的后向气流轨迹, 设定起始点为柳州市中心(109.42°E, 24.33°N), 起始高度为1500 m, 模拟时刻为降雨日09:00, 对应的UTC时间为01:00, 进行48 h后向轨迹模拟追踪和聚类分析.
3 结果和讨论(Results and discussion)3.1 酸雨变化特征分析3.1.1 酸雨年变化特征研究发现, 在自然因素下, 大气降水由于溶解了自然产生的CO₂、NO_x和SO₂, 也会出现降水酸化问题, pH值在5.0~5.6之间, 降水酸化具有自然和人为双重影响性(Galloway et al., 1993).为了描述不同酸度降水的年变化特征, 对降水酸度进行分级统计处理, 具体分级标准为：强酸性(pH＜4. 5)、中度酸性(4. 5≤pH＜5. 0), 弱酸性(5. 0≤pH＜5. 6) (程新金等, 1998).通过计算, 如图 2所示, 2013—2017年柳州市降雨量加权平均pH值始终保持在5.18~6.11之间, pH均值为5.75；强酸雨频率和酸雨频率变化趋势基本保持一致, 酸雨和强酸雨频率均值分别为28.6%、6%；2014年酸化程度最低, pH均值为6.11, 酸雨频率为14.5%, 全年没有出现强酸雨；2016年酸化情况最严重, pH均值为5.18, 酸雨频率和强酸雨频率分别为23%、6.4%.梁宏温等(2002)、韦翠云等(2008)研究表明柳州1996—2006年pH＜5.3, 酸雨频率占70%左右, 相比可得柳州地区近年来降水酸化问题总体呈现出良好的态势.
图 2(Fig. 2)

图 2 酸雨程度逐年变化情况 Fig. 2Yearly variation of intensity and frequency of acid rain, serious acid rain in Liuzhou

3.1.2 酸雨月变化特征如图 3所示, 柳州市降水酸化问题具有明显的季节性差异, 春冬季降水酸度及酸雨频率明显高于夏秋季, 这个差异与我国西南岩溶区存在的降水酸化时间变化一致(刘君峰等, 2006), 分析其原因可能是夏秋季受西南季风的影响, 降水雨量多且强度大, 对污染物的稀释作用明显(张红安等, 2010).近年来柳州降水月pH值的变化趋势基本与酸雨频率和强酸雨频率呈负相关, 春夏季酸雨频率与强酸雨频率持续下降, 对应的pH值则上升；秋冬季酸雨频率变化幅度十分剧烈, 变化范围在3.12%~68%之间；6、7月降水酸度最低, pH出现峰值, 为6.36、6.38, 与此对应的酸雨频率为5.9%、3.12%, 均无强酸雨发生；12月降水酸度最高, pH出现最低值, 为4.89, 酸雨频率和强酸雨频率最高, 分别为68%、63.1%, 变化的一致性较好.
图 3(Fig. 3)

图 3 酸雨程度逐月变化情况 Fig. 3Monthly variation of intensity and frequency of acid rain, serious acid rain in Liuzhou

3.2 降水离子组成分析近年来柳州降水酸化问题有了显著的改善, pH均值达到5.75, 根据研究区所处的岩溶地质背景, 故推断柳州降水pH值偏高可能有两个原因:①酸性离子的减少；②碱性离子较强的中和作用(王自发等, 2007).据表 1所示, 近5年柳州降水中主要离子浓度在年度间出现下降趋势, 与90年代相比均出现了明显下降, 进一步将该区降水与西南、华南、西北等地区相比(表 2), 发现柳州降水主要酸性离子的浓度较低, 这可能受到了柳州近年来为整改城市环境, 大力调整产业结构和能源结构的影响, 因而初步推测柳州降水酸化改善的一个重要原因是酸性离子的减少.如图 4所示, NH₄⁺和Ca²⁺是柳州降水中浓度最高的阳离子, 浓度加权均值分别为58.98、44 μeq·L^-1, 分别占总离子含量的23.3%和17.4%, 两者之和约占总阳离子含量的82.4%；SO₄^2-和NO₃^-则是降水中主要的酸性离子, 浓度加权均值分别为82.35、26.26 μeq·L^-1, 约占总阴离子含量的86%, 两者共占总离子含量的43%, 尤其SO₄^2-占很大比重, 占总离子含量的32.5%, 占总阴离子含量的67%左右.SO₄^2-、NO₃^-浓度的比值不仅可以反映酸性离子浓度的变化, 还是判别酸雨类型的重要依据(程佳惠等, 2015).如图 5所示, 柳州近5年降水的SO₄^2-与NO₃^-的比值出现波动, 但总体在下降, 柳州市大气降水污染有向煤烟型与机动车尾气混合型过渡的趋势, 因此, 应加强机动车尾气的治理工作, 逐步完善机动车排气污染监督管理体系, 逐步提高在用车总体装备水平, 改善大气环境质量.
表 1(Table 1)

表 1 近5年柳州市降水中各种离子成分年均表 Table 1 The major ions concentration of Liuzhou in recent five years

μeq·L^-1
年份	Cl^-	NO₃^-	SO₄^2-	Na⁺	NH₄⁺	K⁺	Mg²⁺	Ca²⁺
2013	17.50	38.69	133.29	8.66	72.87	7.12	21.45	68.94
2014	12.13	18.69	77.03	4.31	35.86	3.24	6.26	20.26
2015	11.67	27.13	77.33	6.81	69.69	4.04	41.32	37.25
2016	15.18	28.69	68.91	9.11	59.38	2.88	7.25	37.25
2017	10.64	18.09	55.20	4.45	57.08	2.72	7.22	38.51
近5年均值	13.42	26.26	82.35	6.67	58.98	4.00	11.30	44.00
90年代均值(郭雅思, 2017)	31.50	85.00	820.20	41.70	251.70	20.80	96.70	423.50

表 2(Table 2)

表 2 柳州同其它城市降水化学成分比较 Table 2 Comparison of the major ions concentration in Liuzhou with other sites

μeq·L^-1
地区	pH	Cl^-	NO₃^-	SO₄^2-	NH₄⁺	K⁺	Na⁺	Ca²⁺	Mg²⁺
柳州	5.75	13.42	26.26	82.35	58.98	4.00	6.67	44.00	11.30
桂林(郭雅思等, 2016)	5.27	36.44	55.00	257.00	120.46	9.20	8.58	186.76	18.64
昆明(王红梅, 2015)	4.80	0.13	47.10	448.20	170.00	—	24.40	358.60	43.50
长沙(吴甫成等, 2000)	4.60	40.30	43.20	210.90	368.60	15.20	69.80	103.60	6.60
贵阳(Han et al., 2004)	3.78	18.03	27.74	215.00	76.67	9.72	18.26	139.72	13.16
南宁(郑凤琴等, 2007)	4.50	21.20	48.20	198.00	—	11.00	4.00	114.00	26.00
深圳(牛彧文等, 2008)	4.91	19.09	15.32	213.40	60.18	7.79	12.67	111.78	10.46
丽江(刘嘉麒等, 1996)	4.68	24.50	23.60	66.10	30.70	5.10	14.40	59.40	10.70
兰州(Xu et al., 2009)	6.60	38.70	128.80	489.70	229.80	13.80	13.80	425.60	36.60
河北(Wu et al., 2016)	7.70	27.90	74.40	208.00	57.20	7.30	12.30	886.00	46.50

图 4(Fig. 4)

图 4 柳州降水离子组成百分数 Fig. 4Ionic composition of precipitation of Liuzhou

图 5(Fig. 5)

图 5 近5年柳州SO₄^2-和NO₃^-浓度比 Fig. 5Ratios of SO₄^2-and NO₃^- in precipitationin of Liuzhou

3.3 离子中和分析目前对于酸雨区的鉴定, 按照年均降水pH＞5.6, 酸雨率在0~20%之间, 为非酸雨区；酸雨率在20%~30%之间, 为轻酸雨区进行划分；近年来柳州地区pH均值达到5.75, 酸雨和强酸雨率均值分别为28.6%、6%, 属于轻酸雨控制区, 将该区降水与全球内降水背景点丽江玉龙雪山相比(表 2), 发现柳州pH值和酸性离子浓度均高于它, 进一步将柳州主要酸性离子与中和离子进行相关性分析, 得到r=0.98(图 6), 因而初步推测柳州降水酸化改善的另一原因是碱性物质的中和效应.
图 6(Fig. 6)

图 6 中和离子、阴离子相关性 Fig. 6Linear regression of acidic ions vs.neutralized ions

通常, 降水中主要酸性离子的被中和状况可以通过计算相对酸度FA (fractional acidity)来评价(Balasubramanian et al., 2001), 见式(1).

(1)

式中, 离子浓度单位为μmol·L^-1, 当FA=1时, 说明酸性离子没有被中和, 计算得到FA =0.0003, 说明99.97%酸性离子发生了中和效应.Al-Momani等(1998)指出酸性离子主要被碱性颗粒物中的NH₃和CaCO₃所中和.降水中主要碱性离子的中和状况可以通过计算Possanzini(1988)提出的中和因子NF(nutralization factor)来评价, 见式(2).

(2)

式中, X_i为待评价的离子浓度(μmol·L^-1), 计算得到NH₄⁺、Ca²⁺、Na⁺、K⁺、Mg²⁺的中和因子分别为0.31、0.23、0.035、0.021和0.059, 说明NH₄⁺、Ca²⁺是最重要的中和离子, NH₄⁺除了农业施肥、生物质燃烧、排泄物释放等, 喀斯特化学风化作用也是不可或缺的来源(Han et al., 2010), Ca²⁺主要是土壤来源(土壤风化作用和沙尘的中-长距离传送)(艾东升等, 2010).
3.4 离子来源3.4.1 离子相关性分析为了对离子来源进行初步判断, 对降水中主要离子做相关性分析, 据表 3所示, SO₄^2-、NO₃^-具有相似的来源和化学形成过程(Al-Khashman et al., 2005), 因而它们相关性显著(r=0.743).通常认为Na⁺主要来源于海洋, Na⁺和Cl^-相关性显著(r=0.819), 表明Cl^-也可能有海洋来源.NH₄⁺、Ca²⁺和SO₄^2-、NO₃^-都存在较好的相关性, 说明它们在中和作用中的贡献都很明显, 此结论和前面计算相对酸度(FA)和中和因子(NF)所得的结论一致.
表 3(Table 3)

表 3 降水中主要离子间的相关系数 Table 3 Spearman correlation matrix for precipitation of Liuzhou

	Cl^-	NO₃^-	SO₄^2-	NH₄⁺	K⁺	Na⁺	Ca²⁺	Mg²⁺
Cl^-	1
NO₃^-	0.792^**	1
SO₄^2-	0.660^**	0.743^**	1
NH₄⁺	0.568^**	0.767^**	0.650^**	1
K⁺	0.709^**	0.767^**	0.729^**	0.753^**	1
Na⁺	0.819^**	0.641^**	0.545^**	0.481^**	0.517^**	1
Ca²⁺	0.565^**	0.671^**	0.526^**	0.611^**	0.599^**	0.559^**	1
Mg²⁺	0.150	0.200	0.640^**	0.335^*	0.408^**	0.060	0.021	1
注：*表示在0.01水平(双侧)上显著相关, 表示在0.05水平(双侧)上显著相关.

3.4.2 富集因子为了进一步判断降水中主要离子来源, 故进行富集因子(EF)计算(Rahn et al., 1976), 通常认为Na⁺完全来自海洋, 将其作为海源参照物, Ca²⁺主要来源土壤, 作为陆源参照物(Keene et al., 1986; Kulshrestha et al., 2003).海洋富集因子(EF_seawater)和土壤富集因子(EF_soil)的计算公式分别如下所示.

(3)

(4)

式中, X代表待计算的离子浓度(μmol·L^-1), [X/Na⁺]_seawater、[X/Ca²⁺]_soil分别为来自海洋和土壤的离子成分组成(Taylor et al., 1964; Négrel et al., 2007), 当EF＞1时, 表示该离子相对于参考离子被富集, 否则被稀释.NO₃^-的海源很少, 可以忽略不计, 因此不进行EF_seawater计算；如表 4所示, SO₄^2-和NO₃^-的EF_seawater和EF_soil两者均远大于1, 表示相对于海水和土壤均被富集, 所以主要是人为来源；Cl^-的EF_seawater值是1.73, EF_soil值是196.77, 明显相对于土壤被富集, 所以主要是海洋来源；K⁺和Mg²⁺的EF_seawater＞1, EF_soil＜1, 相对于海洋富集, 所以主要是土壤来源.
表 4(Table 4)

表 4 柳州降水中主要离子富集系数计算结果 Table 4 Enrichment coefficientof ion composition of precipitationrelative to the sea water and soil

μmol·L^-1
项目	Na⁺	K⁺	Ca²⁺	Mg²⁺	Cl^-	SO₄^2-	NO₃^-
[X/Na⁺]_seawater	—	0.020	0.040	0.230	1.160	0.120	—
[X/Na⁺]_rainwater	—	0.600	3.300	0.850	2.012	6.173	—
EF_seawater	—	27.510	75.130	3.730	1.730	51.440	—
[X/Ca²⁺]_soil	0.569	0.504	—	0.561	0.003	0.019	0.002
[X/Ca²⁺]_rainwater	0.303	0.182	—	0.257	0.610	1.872	1.194
EF_soil	0.530	0.360	—	0.460	196.770	99.55	568.400

降水中主要离子由海源(SSF)、陆源(CF)和人为源(AF)3方面组成, 为了判断各端元对降水中主要离子的贡献量, 故进行不同源的计算, 公式如下所示.

(5)

(6)

(7)

计算结果如表 5所示, SO₄^2-和NO₃^-人为排放分别占到了97.06%和99.824%, 由此可知, 柳州市大气降水中主要离子受化石燃料燃烧以及机动车尾气排放等人类活动影响非常严重, 这与降水酸化逐渐向硫-硝混合型转变的结果相一致；Cl^-是反映海源特征的主要离子, 其海洋源比例占57.65%, 但也有42.34%属于人为源, 这说明柳州市采样点附近工业生产排放的废气对Cl^-浓度的贡献也不容小觑(Sigg et al., 1987)；K⁺除海源外绝大多数来自于生物的燃烧和土壤、岩石风化(Das et al., 2005)；Ca²⁺和Mg²⁺是反映陆源的特征性离子, 其陆源输入比例分别占98.66%、73.2%, 与柳州地区广泛的喀斯特化学风化作用密切相关.
表 5(Table 5)

表 5 柳州降水中主要离子来源所占比例 Table 5 Source contributions for different ionic constituents in precipitation

离子	海洋来源	地壳来源	人为来源
SO₄^2-	1.94%	1.00%	97.06%
NO₃^-	—	0.18%	99.82%
Cl^-	57.65%	0.01%	42.34%
K⁺	3.64%	96.36%
Ca²⁺	1.33%	98.66%
Mg²⁺	26.78%	73.20%

3.5 离子后向气流轨迹判源不同源的降水气团会影响降水中离子的组成和浓度, 在近5年中2016年的pH最低, 故对2016年的60场降水进行后向轨迹追踪模拟和聚类分析, 结果如图 7所示, 气团来源主要分为西北、西南、南、东南、东北5个路径, 将这5个路径降水的pH和离子浓度进行雨量加权平均, 如表 6所示, W方向的气团主要影响春冬季降水, 以高空气流为主, 由于途经云南、贵州等区经济发展和能源结构相对落后的区域, 化石燃料的燃烧、汽车尾气的排放、农业施肥以及生物质的燃烧等活动使降水中SO₄^2-、NO₃^-、NH₄⁺、Ca²⁺、Mg²⁺平均浓度较高, 这部分降水主要受人为源影响.SW方向的气团主要受西南季风影响, 该气团多影响柳州的夏季降水, 约占降水总事件的45%, 因途经南宁、来宾等地区不仅是典型的岩溶酸雨区, 而且是经济发展良好工农业城市, 使得NH₄⁺、NO₃^-和K⁺浓度达到几个路径的最高值, 这部分降水主要受陆源和人为源双重影响.SE方向的气团是从西太平洋吹来的暖湿气流, 以低空气流为主, 此类降水多发生在夏秋季, 由于途经降水酸化污染较重的台湾、珠三角地区, 这也可能是降水中SO₄^2-、NO₃^-浓度较高的原因之一, 另外降水中Na⁺和Cl^-的浓度也是降水中最高的, 表现出较强的海洋性, 这部分降水主要受海源和人为源的影响.
图 7(Fig. 7)

图 7 后向轨迹聚类图 Fig. 7Clustering graph of the air mass trajectory

表 6(Table 6)

表 6 聚类轨迹主要离子分布统计 Table 6 Statistical of the main ions of clustering

μeq·L^-1
来源	pH	NH₄⁺	Ca²⁺	K⁺	Mg²⁺	Cl^-	SO₄^2-	NO₃^-
SW(45%)	5.50	92.00	33.00	3.33	7.50	12.96	75.00	38.55
S(20%)	5.29	41.33	16.50	1.28	4.17	8.17	46.88	15.16
NE(13%)	6.32	90.33	33.00	2.05	5.00	9.58	61.88	15.65
NW(12%)	5.45	84.67	45.00	2.82	6.67	9.86	93.13	38.06
SE(10%)	5.60	80.67	52.50	2.56	7.50	17.18	91.46	35.00

4 结论(Conclusions)1) 近5年柳州市降水pH值范围为5.18~6.11, 加权平均值为5.75；酸雨频率均值为28.6%, 主要以轻酸性降雨为主, 与90年代相比主要离子浓度均有了明显下降, 且降水酸化问题具有明显的季节性差异, 春冬季降水酸化污染相对夏秋季较重.将近年来柳州降水与其它地方相比, 其离子浓度低于西南、华南、北方等城市且高于丽江玉龙雪山, 而pH高于西南、华南和玉龙雪山, 低于北方, 因此推测柳州降水酸化改善的原因是酸性物质的减少和碱性物质的中和作用.
2) 降水的离子组成以SO₄^2-、NO₃^-、Ca²⁺和NH₄⁺为主, SO₄^2-、NO₃^-为柳州降水的主要酸性离子, 约占总阴离子含量的86%, 主要是人为来源.从柳州市降水主要离子SO₄^2-/NO₃^-的比值变化分析发现, 柳州市降水酸化正在向硫-硝混合型过渡, 说明城市里面机动车污染开始占据一定比例, 控制其排放是接下来改善降水酸化的一个重要方向.NH₄⁺和Ca²⁺是柳州降水的主要中和离子, 约占总阳离子含量的82.4%, 主要是土壤来源.酸性离子与中和离子相关系数达到0.98, 中和作用显著, 其中99.7%的酸性离子被中和掉, 其中NH₄⁺和Ca²⁺的中和作用最为明显, 中和系数分别为0.37、0.23, 主要受到喀斯特化学风化作用和农业活动的影响.
3) 对柳州地区降水气团进行后向轨迹模拟追踪及聚类分析后发现, 春冬季降水多受西北、西南方向气团影响；夏秋季降水多受来自西南、东南、东北方向的气团影响；经过海洋的气团所带来的降水的离子成分受海洋影响明显；来自西北、西南方向的降水的离子总浓度最高, 表明柳州市降水离子组成与浓度不仅与本地污染源有关, 还与来自大陆内部西北和西南远距离物质输送有关.

参考文献

艾东升, 郑祥民, 周立旻, 等. 2010. 近2年上海市夏季降水地球化学特征研究[J]. 环境科学, 2010, 31(9): 2002–2009.

Al-Momani I F, Aygun S, Tuncel G. 1998. Wet deposition of major ions and trace elements in the eastern Mediterranean basin[J]. Journal of Geophysical Research:Atmospheres, 103(D7): 8287–8299.DOI:10.1029/97JD03130

Al-Khashman O A. 2005. Ionic composition of wet precipitation in thePetra Region, Jordan[J]. Atmospheric Research, 78(1/2): 1–12.

Akkoyunlu B O, Tayanc M. 2003. Analyses of wet and bulk deposition in four different regions of Istanbul, Turkey[J]. Atmospheric Environment, 37(25): 3571–3579.DOI:10.1016/S1352-2310(03)00349-2

Araguás-Araguás L. 1998. Stable isotope composition of precipitation over southeast Asia[J]. Journal of Geophysical Research Atmospheres, 103(D22): 28721–28742.DOI:10.1029/98JD02582

程佳惠, 李金娟, 孙广权, 等. 2015. 典型酸雨城市降水、降尘中阴阳离子分布特征及其相关性[J]. 环境科学学报, 2015, 35(6): 1676–1682.

Báez A, Belmont R, García R, et al. 2007. Chemical composition of rainwater collected at a southwest site of Mexico City, Mexico[J]. Atmospheric Research, 86(1): 61–75.DOI:10.1016/j.atmosres.2007.03.005

Balasubramanian R, Victor T, Chun N. 2001. Chemical and statistical analysis of precipitation in singapore[J]. Water, Air, and Soil Pollution, 130(1/4): 451–456.DOI:10.1023/A:1013801805621

程新金, 黄美元. 1998. 降水化学特性的一种分类分析方法[J]. 气候与环境研究, 1998, 3(1): 82–88.DOI:10.3878/j.issn.1006-9585.1998.01.08

丁国安, 徐晓斌, 王淑凤, 等. 2004. 中国气象局酸雨网基本资料数据集及初步分析[J]. 应用气象学报, 2004, 15(z1): 85–94.DOI:10.3969/j.issn.1001-7313.2004.z1.012

董蕙青, 黄海洪, 高安宁, 等. 2005. 影响广西酸雨的大气环流特征分析[J]. 气象, 2005, 31(10): 51–55.DOI:10.7519/j.issn.1000-0526.2005.10.013

Das R, Das S N, Misra V N. 2005. Chemical composition of rainwater and dustfall at Bhubaneswar in the east coast of India[J]. Atmospheric Environment, 39(32): 5908–5916.DOI:10.1016/j.atmosenv.2005.06.030

Ding G A, Xu X B, Fang X M, et al. 1997. Current status and future of acid rain in China[J]. Chinese Science Bulletin, 42(24): 2076–2076.DOI:10.1007/BF02882950

郭雅思, 于奭, 黎泳珊, 等. 2016. 桂林市酸雨变化特征及来源分析[J]. 环境科学, 2016, 37(8): 2897–2905.

郭雅思.2017.柳州市酸雨化学特征分析及人为硫源对碳酸盐岩风化的影响[D].重庆: 西南大学

Galloway J N, Likens G E, Edgerton E S. 1976. Acid precipitation in the Northeastern United States:pH and acidity[J]. Science (New York, N.Y.), 194(4266): 722–724.DOI:10.1126/science.194.4266.722

Galloway J N, Savoie D L, Keene W C, et al. 1993. The temporal and spatial variability of scavenging ratios for NSS sulfate, nitrate, methanesulfonate and sodium in the Atmosphere over the North Atalantic Ocean[J]. Atmospheric Environment, 27(2): 235–250.DOI:10.1016/0960-1686(93)90354-2

Han G L, Liu C Q. 2004. Strontium isotope and chemical compositions of rainwaters from Guiyang, Guizhou Province, China[J]. Geochimica et Cosmochimica Acta, 68(11): A461.

Han G L, Tang Y, Wu Q X. 2010. Chemical and strontium isotope characterization of rainwater in karst virgin forest, Southwest China[J]. Atmospheric Environment, 44(2): 174–181.DOI:10.1016/j.atmosenv.2009.10.019

Keene W C, Pszenny A A P, Galloway J N. 1986. Sea-salt corrections and interpretation of constituent ratios in marine precipitation[J]. Journal of Geophysical Research:Atmospheres, 91(D6): 6647–6658.DOI:10.1029/JD091iD06p06647

Krupa S V. 2002. Sampling and physico-chemical analysis of precipitation:a review[J]. Environmental Pollution, 120(3): 565–594.DOI:10.1016/S0269-7491(02)00165-3

Kulshrestha U C, Kulshrestha M J, Sekar R. 2003. Chemical characteristics of rainwater at an urban site of south-central India[J]. Atmospheric Environment, 37(21): 3019–3026.DOI:10.1016/S1352-2310(03)00266-8

梁宏温, 温远光, 何斌, 等. 2002. 柳州市中心酸雨化学特征与成因分析[J]. 广西农业生物科学, 2002, 21(4): 261–265.

刘嘉麒. 1996. 降水背景值与酸雨定义研究[J]. 中国环境监测, 1996(5): 7–11.

刘君峰, 宋之光, 许涛. 2006. 广州地区雨水化学组成与雨水酸度主控因子研究[J]. 环境科学, 2006, 27(10): 1998–2002.DOI:10.3321/j.issn:0250-3301.2006.10.013

Lei X E, Han Z W, Zhang M G, et al. 1998. Physical, chemical, biological processes and mathematical model on air pollution[M]. Beijing: China Meteorological Press.

Négrel P, Guerrot C, Millot R. 2007. Chemical and strontium isotopecharacterization of rainwater in France:influence of sources and hydrogeochemical implications[J]. Isotopes in Environmental and Health Studies, 43(3): 179–196.DOI:10.1080/10256010701550773

牛彧文, 何凌燕, 胡敏. 2008. 深圳大气降水的化学组成特征[J]. 环境科学, 2008, 29(4): 1014–1019.DOI:10.3321/j.issn:0250-3301.2008.04.029

Okuda T, Iwase T, Ueda H, et al. 2004. Long-term trend of chemical constituents in precipitation in Tokyo metropolitan area, Japan, from 1990 to 2002[J]. Science of the Total Environment, 339(1): 127–141.

Possanzini M, Buttini P, Di Palo V. 1988. Characterization of a rural area in terms of dry and wet deposition[J]. Science of the Total Environment, 74: 111–120.DOI:10.1016/0048-9697(88)90132-5

Rahn K A. 1976. Silicon and aluminum in atmospheric aerosols:Crust-air fractionation[J]. Atmospheric Environment, 10(8): 597–601.DOI:10.1016/0004-6981(76)90044-5

Salam A, Bauer H, Kassin K, et al. 2003. Aerosol chemical characteristics of a mega-city in Southeast Asia (Dhaka-Bangladesh)[J]. Atmospheric Environment, 37(18): 2517–2528.DOI:10.1016/S1352-2310(03)00135-3

Sigg L, Stumm W, Zobrist J. 1987. The chemistry of fog:factors of regulating its composition[J]. Chimia International Journal for Chemistry, 41(41): 159–165.

Taylor S R. 1964. Abundance of chemical elements in the continental crust:a new table[J]. Geochimica Et Cosmochimica Acta, 28(8): 1273–1285.DOI:10.1016/0016-7037(64)90129-2

王红梅. 2005. 浅析昆明市酸雨污染状况及变化[J]. 云南环境科学, 2005, 24(3): 51–54.DOI:10.3969/j.issn.1673-9655.2005.03.017

王文兴. 1994. 中国酸雨成因研究[J]. 中国环境科学, 1994, 14(5): 323–329.DOI:10.3321/j.issn:1000-6923.1994.05.008

王自发, 高超, 谢付莹. 2007. 中国酸雨模式研究回顾与所面临的挑战[J]. 自然杂志, 2007, 29(2): 78–82.DOI:10.3969/j.issn.0253-9608.2007.02.004

韦翠云, 蔡丽. 2008. 柳州酸雨分布特征分析[J]. 气象研究与应用, 2008, 29(Z1): 22.DOI:10.3969/j.issn.1673-8411.2008.z1.013

吴甫成, 吴君维, 王晓燕, 等. 2000. 湖南酸雨污染特征[J]. 环境科学学报, 2000, 20(6): 807–809.DOI:10.3321/j.issn:0253-2468.2000.06.031

吴起鑫, 韩贵琳, 陶发祥, 等. 2011. 西南喀斯特农村降水化学研究:以贵州普定为例[J]. 环境科学, 2011, 32(1): 26–32.

Wu Y, Xu Z F, Liu W J. 2016. Chemical compositions of precipitation at three non-urban sites of Hebei Province, North China:Influence of terrestrial sources on ionic composition[J]. Atmospheric Research, 181: 115–123.DOI:10.1016/j.atmosres.2016.06.009

Xu Z F, Li Y S, Tang Y. 2009. Chemical and strontium isotopecharacterization of rainwater at an urban site in Loess Plateau, northwest China[J]. Atmospheric Research, 94(3): 481–490.DOI:10.1016/j.atmosres.2009.07.005

张红安, 汤洁, 于晓岚, 等. 2010. 侯马市酸雨长期变化趋势分析[J]. 环境科学学报, 2010, 30(5): 1069–1078.

郑凤琴, 孙崇智, 谢宏斌. 2007. 1991-2004年不同环流背景下南宁市降水化学成分变化特征的研究[J]. 热带气象学报, 2007, 23(6): 664–668.DOI:10.3969/j.issn.1004-4965.2007.06.021

周晓得, 徐志方, 刘文景, 等. 2017. 中国西南酸雨区降水化学特征研究进展[J]. 环境科学, 2017, 38(10): 4438–4446.