1. 引言
空气负离子是带有负电荷的气体分子和轻离子团的总称 [1],被称为“空气维生素”。空气负离子主要为负氧离子,它具有杀菌、降尘、清洁空气等功效,对人体健康有积极意义。目前,空气负氧离子水平已成为评价旅游环境空气清洁度的重要指标 [2] [3] [4],其研究为旅游环境评价、生态康养旅游规划提供了科学依据。对空气负氧离子浓度时空分布特征的研究表明 [5] [6] [7] [8] [9],负氧离子浓度空间分布不均匀,由市中心向郊区逐渐增大,随海拔高度的增加而增大;负氧离子浓度有明显的季节变化和日变化规律。我国各地气象部门积极开展了负氧离子浓度的监测和分析工作,北京、天津和杭州等地已经开展了负氧离子浓度预报业务,并取得了较好的社会效益。
然而,以往对负氧离子的研究所使用的观测资料大多为间断性的,数据周期为几个月,连续性不好、周期短、统计规律不明显,而且将研究结果应用于实际业务,建立负氧离子预测模型方面的研究也较少 [10] [11] [12] [13]。温州地处东南沿海,境内山川秀美,旅游资源丰富,但缺少对旅游生态环境的研究,负氧离子与气象因素的关系方面尚未深入分析和研究。本文利用连续观测的负氧离子资料通过对温州地区负氧离子浓度分布特征的分析,及其与气象影响因素相关性的研究,通过多元线性回归方法,建立温州市负氧离子浓度预报方程,实现负氧离子的定量化预报预测,为生态康养旅游提供科学依据和参考信息。
2. 资料与方法
2.1. 研究资料
负氧离子浓度测量采用北京华思通科技有限公司WST-5B型大气负氧离子负氧离子监测设备。温州地区共布设了10个负氧离子气象观测点,其中乐清、永嘉、温州、泰顺、文成、平阳、瑞安、苍南和洞头这9个站位于各县市区气象观测站,温州站代表温州市区,大罗山站位于温州市区大罗山风景区内,负氧离子气象观测点的分布见图1。
Figure 1. The spatial distribution of oxygen anion concentrations from March 2019 to February 2020 in Wenzhou
图1. 2019年3月~2020年2月温州地区负氧离子浓度分布
华思通负氧离子数据从2018年3月开始,考虑到设备建成初期出现运行不稳定的情况,故选取2019年3月~2020年2月整1年的负氧离子浓度观测资料,观测间隔为每小时一次。监测的离子值分为两档,即迁移速率(即离子运动速度)大于等于1.0 cm∙s−1的大离子和大于等于0.4 cm∙s−1的小离子两档。因小离子的医疗保健等作用更突出,本文选取小离子为研究对象。在2019年3月~2020年2月期间,乐清站共有8537个观测值,大罗山站共有8482个观测值,永嘉站共有8539个观测值,温州站共有8546个观测值,泰顺站共有8556个观测值,文成站共有8555个观测值,平阳站共有8214个观测值,瑞安站共有8544个观测值,苍南站共有8526个观测值,洞头站共有8537个观测值。剔除异常值之后,最终得到可利用的乐清站8426个观测值、大罗山站7176个观测值、永嘉站8478个观测值、温州站8340个观测值、泰顺站8152个观测值、文成站7674个观测值、平阳站6994个观测值、瑞安站7951个观测值、苍南站7843个观测值、洞头站7989个观测值。温度、相对湿度、气压等气象资料选用2019年3月~2020年2月每小时一次的自动站观测数据。
2.2. 研究方法
在回归分析中,如果有两个或两个以上的自变量,就称为多元回归。事实上,一种现象常常是与多个因素相联系的,由多个自变量的最优组合共同来预测或估计因变量,比只用一个自变量进行预测或估计更有效,更符合实际。因此多元线性回归比一元线性回归的实用意义更大。假设因变量Y和k个自变量
之间存在线性相关关系,Y和X之间的函数关系可以表示为:
(1)
式中,
为回归常数,
为回归系数,
为回归残差。若对因变量进行n次观测,其观测值为
,其中
,则n个观测值可以表示为:
(2)
式中,
为回归残差,假定它们相互独立且服从同一正态分布
。
多元线性回归分析就是通过回归拟合求解回归参数
的估计值,建立多元线性回归方程,并对建立的回归方程及其回归系数的显著性进行检验。多元线性回归的数据通常需要满足以下统计假设。① 正态性:对于固定的自变量取值,因变量值呈正态分布;② 线性:因变量和自变量之间线性相关;③ 同方差性:因变量的方差不随自变量的水平不同而变化;④ 独立性:因变量之间相互独立,不存在多重共线性。
本文利用多元回归分析方法,将负氧离子浓度作为因变量,雨量、气温、湿度、风速等气象因素作为自变量,采用Matlab软件进行多元线性回归拟合而建立温州市负氧离子浓度的预测模型。
3. 结果分析
3.1. 负氧离子浓度时空分布特征
3.1.1. 地域分布特征
2019年3月~2020年2月温州地区平均负氧离子分布见图1,泰顺的负氧离子浓度最高,为102个/cm3,其次为大罗山75个/cm3,这两个站均为高山站,海拔高度在500 m以上。平阳、洞头、乐清的负氧离子浓度较低,均在50个/cm3以下,其中平阳最低,为39个/cm3。永嘉、温州、瑞安、文成、苍南的负氧离子浓度较为接近,在50~60个/cm3之间。
研究表明 [3] [4] [5] [6] [7],空气负氧离子浓度与人类活动、周围绿化环境等有关,市中心附近的负氧离子浓度远小于郊区和植被覆盖率高的森林山区。温州的观测站中,泰顺和大罗山处于高海拔山区,周围山林茂密,植物数量多,环境污染少,空气中负氧离子浓度相对较高。
3.1.2. 负氧离子浓度季节和月变化
以3~5月、6~8月、9~11月和12月至翌年2月分别代表春夏秋冬四季。由表1可见,乐清、大罗山、永嘉、平阳、瑞安负氧离子浓度以冬季最高,春季最低;温州、泰顺、洞头负氧离子浓度以夏季最高,春季最低;文成负氧离子浓度以冬季最高,夏季最低;苍南负氧离子浓度以秋季最高,春季最低。
负氧离子浓度的上述季节变化可能与天气状况有关,春季温州大部分县市区大雾、阴天日数较多,空气质量较差;而冬季温州大部分县市区平均风速大,烟尘很快被稀释,冬季晴天日数也远远多于春季。曾曙才等 [14] 的研究表明,晴天负氧离子浓度要高于阴天。夏季负氧离子浓度高,与夏季雷雨天气多、空气湿度大、草木繁茂吸收气溶胶粒子多等有关。
Table 1. The seasonal distribution of oxygen anion concentrations in Wenzhou (unit: number per cubic centimeter)
表1. 温州地区负氧离子浓度季节分布(单位:个/cm3)
图2为温州地区空气负氧离子浓度的月变化,可见除了泰顺表现为多峰特征外,其余站点的月变化表现为单峰或双峰特征。乐清负氧离子浓度的月变化表现为双峰特征,峰值出现在冬季的1月、2月和夏季的8月;大罗山负氧离子浓度的月变化表现为单峰特征,峰值出现在冬季的1月;永嘉负氧离子浓度的月变化表现为单峰特征,峰值出现在冬季的2月;温州负氧离子浓度的月变化表现为双峰特征,峰值出现在冬季的2月和夏季的7月;文成、平阳、瑞安负氧离子浓度的月变化表现为单峰特征,峰值出现在冬季的1月;苍南负氧离子浓度的月变化表现为双峰特征,峰值出现在冬季的1月和秋季的9月;洞头负氧离子浓度的月变化表现为单峰特征,峰值出现在夏季的6月;泰顺负氧离子浓度的月变化表现为多峰特征,峰值出现在冬季的1月和2月,夏季的6月和8月,以及秋季的11月。
3.1.3. 负氧离子浓度日变化
图3为温州地区各站2019年3月~2020年2月平均的负氧离子浓度日变化。由图可见,乐清、永嘉、温州、平阳、瑞安和洞头负氧离子浓度的日变化相对较小;大罗山、泰顺、文成和苍南负氧离子浓度的日变化相对明显,大罗山和苍南负氧离子浓度分布总体上呈U型,上午9时之前较高,其后开始下降,14时左右降到最低,然后浓度开始逐渐增大,直到第二日上午;泰顺和文成负氧离子浓度分布总体上呈W型,早晨和傍晚处于低值区,14时左右达到峰值。
大罗山和苍南夜间负氧离子浓度较高,白天负氧离子浓度较低,最大值出现在早晨时段,最小值出现在中午时段,因此建议人们到大罗山和苍南康养旅游时在早晨出门,可以呼吸到负氧离子含量相对较高的清新空气。泰顺和文成中午时段负氧离子浓度较高,早晨和傍晚负氧离子浓度较低,因此建议人们到泰顺和文成康养旅游时在中午出门,可以呼吸到负氧离子含量相对较高的清新空气。
3.2. 负氧离子与气象因素的关系及预测模型
3.2.1. 负氧离子浓度与气象因素的关系
对温州地区各站做2019年3月~2020年2月日平均负氧离子浓度与气温、相对湿度、风速等气象因素的相关性分析,选用日雨量、日平均气温、日最高气温、日最低气温、日最大风速、日极大风速、日平均湿度等7个气象因素。
由表2可知,乐清站日雨量、最大风速和极大风速这3个气象因子与负氧离子浓度的相关性显著,均为正相关,相关系数在0.4以上;大罗山站与平均气温、最高气温、最低气温、平均湿度显著负相关(相关系数为−0.34~−0.2),与极大风速显著正相关(相关系数为0.22);永嘉站与平均气温、最高气温、最低气温显著负相关(相关系数为−0.39~−0.24),与日雨量、最大风速、极大风速显著正相关(相关系数为0.21~0.26);文成站与平均气温、最高气温、最低气温显著负相关(相关系数为−0.54~−0.43);平阳站与平均气温、最高气温、最低气温显著负相关(相关系数为−0.34~−0.31),与最大风速、极大风速显著正相关(相关系数为0.41~0.44);瑞安站与极大风速显著正相关(相关系数为0.25);温州站、泰顺站、苍南站和洞头站负氧离子浓度与上述7个气象因素相关性不明显。
Table 2. Correlation coefficient between diurnal oxygen anion concentration and meteorological elements in Wenzhou
表2. 温州地区日负氧离子浓度与气象因素的相关系数
3.2.2. 负氧离子预测模型
由3.2.1节可知,温州站、泰顺站、苍南站和洞头站负氧离子浓度变化与气象因素基本无关,因此无法通过筛选多元回归预报因子建立预测方程。乐清站选择日雨量和极大风速作为多元回归预报因子。大罗山站选择平均气温、最高气温、最低气温、极大风速、平均湿度作为多元回归预报因子。永嘉站选择日雨量、平均气温、最高气温、最低气温、极大风速作为多元回归预报因子。文成站选择平均气温、最高气温、最低气温作为多元回归预报因子。平阳站选择平均气温、最高气温、最低气温、极大风速作为多元回归预报因子。瑞安站选择极大风速作为一元回归预报因子。
本文利用入选预报因子和负氧离子浓度资料,进行多元回归计算,最后得到乐清、大罗山、永嘉、文成、平阳和瑞安日负氧离子浓度预报回归方程:
乐清回归方程:
(3)
大罗山回归方程:
(4)
永嘉回归方程:
(5)
文成回归方程:
(6)
平阳回归方程:
(7)
瑞安回归方程:
(8)
式(3)~式(8)中,y为空气负氧离子浓度;R为日雨量;T为日平均气温;
为日最高气温;
为日最低气温;V为日极大风速;RH为日平均相对湿度。
选取2020年4~6月作为负氧离子浓度预测模型的检验样本,其中剔除异常观测值之后,乐清和文成仅选取4月检验样本,瑞安选取4~5月检验样本。将式(3)~式(8)中气象要素的观测值分别代入相应的回归方程,计算出模型预报的空气负氧离子浓度,并与实际观测值进行对比(图4)。从图4可以看出,模型预报值和实际观测值的变化趋势有较好的一致性,乐清、大罗山、永嘉、文成、平阳和瑞安观测与预报值的相关系数分别为0.31、0.33、0.3、0.23、0.2和0.45,各站模型计算的预报值与观测值显著相关。除个别样本外,预报值与观测值在量级上较为接近,方程具有较好的预报能力。
需要说明的是,在个别极值样本中,预报与观测值的误差略大,这可能是因为影响负氧离子浓度的要素众多,预测模型只引入了有限的气象因子,故模型难以全面体现。另外,个别观测极值可能是观测仪器出现故障导致的,观测值本身不准确。
(a)(b) (c)(d) (e)(f)
Figure 4. Comparison between the observed and simulated oxygen anion concentration in Yueqing (a), Daluoshan (b), Yongjia (c), Wencheng (d), Pingyang (e) and Ruian (f) of Wenzhou
图4. 温州地区乐清(a)、大罗山(b)、永嘉(c)、文成(d)、平阳(e)和瑞安(f)负氧离子浓度观测值和预报值比较
温州、泰顺、苍南和洞头负氧离子浓度变化与气象因素相关性小,因此无法通过筛选多元回归预报因子建立预测方程,但可以通过温州与瑞安、泰顺与大罗山、苍南与瑞安、洞头与乐清负氧离子浓度接近的关系以及上述6个回归方程实现温州、泰顺、苍南和洞头负氧离子浓度的预测。
4. 结论与讨论
(1) 温州地区泰顺的年均负氧离子浓度最高,为102个/cm3,其次为大罗山75个/cm3,这两个站均为高山站。平阳、洞头、乐清的负氧离子浓度较低,均在50个/cm3以下,其中平阳最低,为39个/cm3。永嘉、温州、瑞安、文成、苍南的负氧离子浓度较为接近,在50~60个/cm3之间。
(2) 乐清、大罗山、永嘉、平阳、瑞安负氧离子浓度以冬季最高,春季最低;温州、泰顺、洞头负氧离子浓度以夏季最高,春季最低;文成负氧离子浓度以冬季最高,夏季最低;苍南负氧离子浓度以秋季最高,春季最低。
(3) 乐清、永嘉、温州、平阳、瑞安和洞头负氧离子浓度的日变化相对较小;大罗山和苍南夜间负氧离子浓度较高,白天负氧离子浓度较低,最大值出现在早晨时段,最小值出现在中午时段,因此建议人们到大罗山和苍南康养旅游时在早晨出门,可以呼吸到负氧离子含量相对较高的清新空气;泰顺和文成中午时段负氧离子浓度较高,早晨和傍晚负氧离子浓度较低,因此建议人们到泰顺和文成康养旅游时在中午出门,可以呼吸到负氧离子含量相对较高的清新空气。
(4) 影响负氧离子浓度的气象因子主要是雨量、气温、风速和湿度。总体上,负氧离子浓度与雨量、风速呈正相关关系,与气温、湿度呈负相关关系。
(5) 根据空气负氧离子浓度与气象条件之间的相关性,采用多元回归方法建立了温州市负氧离子预测模型,经检验模型预报效果良好。负氧离子预测模型的建立实现了温州地区负氧离子浓度的定量化预报,这为提高生态旅游气象服务能力具有重要意义。
由于观测资料问题,所建立的预测模型较粗略,有待于进一步改进。随着观测资料的完善,对负氧离子浓度将做进一步分析研究,模型预报准确率也将进一步提高。
基金项目
温州市气象局项目(2018ZD01)资助。