0.   引言

【研究意义】大气颗粒物是指悬浮在空气中微小固体和液体小滴的混合物，是雾、烟和空气尘埃的主要成分。其中按颗粒物直径可分为总悬浮颗粒物（TSP，d≤100 μm）、PM₁₀（d≤10 μm）和PM_2.5（d≤2.5 μm）^[1]。由于城市化进程加快以及能源消耗不断增加，我国大气颗粒物污染日益严重，尤以冬季为甚，危害巨大^[2-5]。植物能够阻挡、过滤和吸滞大气颗粒物^{[6, 7]}。不同的绿化树种，对于不同粒径颗粒物的滞尘能力有明显差异，在改善城市大气颗粒物污染中起重要作用^{[1, 4, 8]}。【前人研究进展】目前，国内外学者在植物滞尘能力方面进行了一定的研究，例如植物的滞尘能力受到环境因子、生理特性和叶片结构的影响^[9-12]。不同地区环境污染程度存在较大差异，植物为了适应环境会改变自身的生理特性^[13]。在大气污染严重的地区，植物叶片滞留的大量颗粒物会导致其光合作用减弱，进而使叶片中的叶绿素a、b破坏分解^[14]。同时在该地区，植物叶片内活性氧的产生和清除平衡也会受到破坏，引发膜脂过氧化，最终导致丙二醛含量升高^[14]。而且，在植物叶片周围悬浮的颗粒物很容易造成植物叶片的气孔堵塞，此时气孔密度增加有利于提高叶片单位叶面积的气孔交换能力，使其能适应污染严重的环境^[13]。因此，环境污染程度也影响着叶片表面结构的变化^{[10, 15, 16]}。植物叶表结构会影响其滞尘能力，叶表粗糙多毛的树种的滞尘能力显著高于叶表光滑的树种^[17]。此外，叶表滞留的颗粒物集中分布在中脉和凹槽附近^[18]。【本研究切入点】郑州由于冬季采暖，大气颗粒物污染严重，而部分植物存在落叶无法滞尘的情况，此时常绿树种在冬季滞尘的作用显得尤为重要^[19-21]。目前，郑州市常绿树种综合抗污染性相关研究较少。【拟解决的关键问题】以郑州市7种常见的常绿树种为研究对象，测定并计算出各树种单位叶面积滞留不同粒径颗粒物（TSP、PM₁₀和PM_2.5）的质量，比较分析不同污染程度下各树种叶片的生理指标（叶绿素a、b和丙二醛）和叶表面微形态结构，探讨不同污染程度下不同树种的抗污染能力，为筛选缓解郑州市空气粉尘污染、提升空气质量的园林绿化树种提供了一定的理论依据。

1.   材料与方法

1.1   试验地点

郑州市地处河南省中北部（34°16′～34°58′N，112°42′～114°14′E），属暖温带大陆性季风气候，四季分明，年降水量652.9 mm，且多集中在夏季，冬季降水较少，大气颗粒物污染频发^[22]。本次采样地点设置在郑州市中原区，选取3个有代表性的采样区，分别为新技术产业开发区热力公司（简称工厂）、陇海西路（简称道路）和郑州植物园（简称公园）。现场测得的空气颗粒物含量分别为146.4 μg·m⁻³、101.3 μg·m⁻³、53.8 μg·m⁻³，对应设为高、中和低3个不同污染等级（图1）。

图  1  采样区域分布

注：1、公园；2、道路；3、工厂。图5同。

Figure  1.  Schematic distribution of sampling area

Note: 1. park; 2. road; 3. factory. The same was applied in Fig.5.

下载: 全尺寸图片 幻灯片

1.2   供试树种

在3处不同污染程度采样区选择7种郑州市典型常绿树种为研究对象，分别为针叶树种：圆柏（Sabina chinensis）和雪松（Cedrus deodara），阔叶树种：女贞（Ligustrum lucidum）、石楠（Photinia serratifolia）、大叶黄杨（Buxus megistophylla）、海桐（Pittosporum tobira）和南天竹（Nandina domestica）。其中，2种针叶树种在郑州市主要公园的平均出现频率为100%，5种阔叶树种平均出现的频率为80%。每种树种随机选择3株，要求生长健壮无病害，生长状况一致。

1.3   叶片采集方法

采样时间为郑州冬季供暖后，于 2020年11 月雨后第7天开始（微风无降水的天气条件）在3处采样地区进行采集，各供试树种设置3个重复，共采集3次。每次采样采取分层取样法，将树冠分为上、中、下3个垂直部分和东西南北的4个方向，随机均匀地采集健康完整的叶片，单叶较大树种采取10～15片，较小的树种采取30～45片，小心放入样品袋中，并立即带回实验室处理。

1.4   叶面滞尘量的测定

采用改良的分级滤膜过滤法测定各树种叶片滞留 TSP、PM₁₀和PM_2.5的质量^[4]。通过便携式叶面积仪（YMJ-B，浙江托普）对植物叶面积进行测定；针叶树采用王会霞的排水法计算叶面积^[23]。每个树种叶面积的测定重复3次，分别计算各植物 TSP、PM₁₀、PM_2.5滞留量与叶片总面积的比值，得到单位叶面积各粒径颗粒物的滞留量。

1.5   生理指标测定

选取新鲜的叶片对植物其他生理指标进行测定。每项试验重复测定3 次，将叶片洗净后放在65 ℃的烘箱中烘干，除去叶脉后剪碎，叶绿素a、b值采用改良的80%丙酮溶液浸提比色分析测定^[14]；丙二醛（MAD）采用硫代巴比妥酸（TBA）显色法测定^[24]。

1.6   叶表面微结构观测

将叶片平放，切成1 cm×1 cm正方形，用超景深三维视频显微镜（Leica-DVM6A，德国）在60倍的放大倍数下观察各树种叶表面靠近中脉和远离中脉的颗粒物的分布情况以及叶表面的形态结构。随后刮去上表皮和叶肉组织，只留下表皮，在1200倍的放大倍数下观察并统计气孔的大小、数量等。

1.7   数据处理

采用SPSS 23.0软件进行单因素方差分析（ANOVA），确定不同树种各粒径颗粒物的滞尘能力以及生理指标的差异，并用最小显著差数法（LSD）和Dunnett’s T2检验数据之间的差异显著性，采用Pearson相关性分析法分析植物叶片滞尘能力与生理指标之间的关系，采用Excel 2010、Origin 2007绘图。

2.   结果与分析

2.1   不同程度污染下树种滞尘能力比较

如图2所示，各树种单位叶面积滞留TSP、PM₁₀、PM_2.5的变化范围分别为0.385～29.012 g∙m⁻²，0.211～13.033 g∙m⁻²和0.134～12.109 g∙m⁻²。不同树种滞尘量差异显著，针叶树种单位叶面积滞尘量显著高于阔叶树种，针叶树种中圆柏的滞尘量高于雪松，阔叶树种中单位叶面积滞尘量最高为女贞，其次依次为石楠、大叶黄杨和海桐，滞留量最低的为南天竹。其中石楠与大叶黄杨、海桐与南天竹之间滞尘量存在显著性差异（P＜0.05）。同一树种在不同污染程度采样区的滞尘量差异也较大，污染越严重的地区树种滞尘量越高。单位叶面积滞留TSP能力最强的圆柏和雪松在工厂的滞留量是其在清洁区公园的4.3倍和5.6倍，而TSP滞留能力较差的大叶黄杨和海桐在工厂的滞留量是其在公园滞留量的3.6倍和3.8倍。

图  2  7种常绿树种在3个采样区滞尘量的比较

注：不同小写字母表示在同一采样区下不同树种间差异显著（P＜0.05）；A、圆柏；B、雪松；C、女贞；D、石楠；E、大叶黄杨；F、海桐；G、南天竹。图3、图4和图5同。

Figure  2.  Dust deposition on 7 species of evergreen in 3 sampling areas

Note: Different lowercase letters indicate significant differences among different tree species in the same sampling area (P＜0.05); A. Sabina chinensis. B. Cedrus deodara. C. Ligustrum lucidum. D. Photinia serratifolia. E. Buxus megistophylla. F. Pittosporum tobira. G. Nandina domestica. The same for Fig.3, 4, 5.

下载: 全尺寸图片 幻灯片

2.2   不同污染程度对树种叶片生理特性的影响

7种树种在3个不同污染程度采样区的叶绿素含量差异较大，污染越严重叶绿素含量越低；而丙二醛含量则随污染程度的加深而增加（图3）。相同污染程度下7种树种的生理响应程度也存在差异，石楠和大叶黄杨的叶绿素a、叶绿素b含量显著高于其他绿化树种（P＜0.05），圆柏和雪松的丙二醛含量显著高于其他绿化树种（P＜0.05）。在污染严重的情况下，圆柏、雪松的叶绿素a含量变化较小，分别为0.06 mg∙g⁻¹和0.03 mg∙g⁻¹，叶绿素b含量变化较小的树种为圆柏、雪松和女贞，分别为0.01 mg∙g⁻¹、0.04 mg∙g⁻¹和0.04 mg∙g⁻¹，而海桐和南天竹的丙二醛含量变化较小，分别为3.82 μmol∙g⁻¹和2.71 μmol∙g⁻¹。

图  3  7种常绿树种在3处采样区的叶绿素a、b以及丙二醛含量比较

Figure  3.  Contents of chlorophyll a, chlorophyll b, and malondialdehyde of evergreens in 3 sampling areas

下载: 全尺寸图片 幻灯片

2.3   不同植物叶面结构分析

利用超景深三维视频显微镜对7种常绿树种叶片的主脉、侧脉等进行观察（图4）。针叶树种圆柏和雪松的叶表面附着大量的颗粒物，颗粒物滞留量明显高于阔叶树种（图4A、4B）；阔叶树种女贞和石楠的叶表面均有明显的主脉和侧脉结构，附着有较多的颗粒物（图4C、4D），二者在阔叶树种中滞留颗粒物最多；大叶黄杨上表面沟槽较浅或较平整，粉尘固着较分散，除中脉外有少量的颗粒物滞留（图4E），总体颗粒物滞留量一般；海桐的褶皱和沟槽中滞留着一定数量的颗粒物，且在除中脉外的区域仍有一些颗粒物滞留，整体上海桐颗粒物滞留量一般（图4F）；南天竹极窄的中脉有少量的颗粒物滞留，除中脉外叶表面其余部分光滑无毛，无颗粒物滞留，因此整体叶片粉尘滞留最少（图4G）。

图  4  7种常绿树种的叶面颗粒物分布

注：箭头为不同植物叶面颗粒分布集中区域（×60倍）。

Figure  4.  Distribution of particles fell on leaves of various species of evergreens

Note: Arrows are the concentrated areas of leaf particles of different plants (×60 times).

下载: 全尺寸图片 幻灯片

由表1和图5可知，同一植物在不同污染程度下气孔密度和气孔大小不同，除圆柏和雪松外，女贞、石楠、大叶黄杨、海桐、南天竹的气孔密度均随污染程度增加而变大，而气孔长度和气孔宽度均随污染程度增加而变小。随着污染程度增加，7种树种中，雪松和圆柏的气孔密度变化较小，石楠和海桐变化较大；气孔长度变化较小的为女贞和圆柏，变化较大的为雪松和海桐；南天竹和女贞的气孔宽度变化较小，海桐和圆柏变化较大。利用SPSS中一般线性模型，以污染程度和树种差异为固定因子，以气孔长度、气孔宽度和气孔密度为因变量进行方差分析，探究树种差异和污染程度对气孔指标的影响程度。方差分析结果表明（表2），污染程度对气孔密度、气孔长度和气孔宽度的影响均不显著（P＞0.05），而树种差异对气孔密度、气孔长度和气孔宽度均造成了极显著的影响（P＜0.01）。

表  1  7种常绿树种在3处采样区的叶面气孔参数

Table  1.  Stomatal parameters of evergreens in 3 sampling area

树种 Tree species	气孔形态 Stomatal morphology	采样区域 Sampling area	气孔密度 Stomatal density/(个·mm−2)	气孔长度 Stomatal length/μm	气孔宽度 Stomatal width/μm
圆柏 Sabina chinensis	长圆形 Long round	1	55.18±0.00	44.21±0.92	36.15±0.78
		2	55.18±0.00	42.37±0.30	32.60±0.39
		3	66.68±3.98	42.11±0.51	20.93±0.49
雪松 Cedrus deodara	长圆形 Long round	1	27.59±0.00	59.28±0.91	38.67±0.30
		2	27.59±0.00	44.26±0.32	37.60±0.35
		3	27.59±0.00	31.98±1.16	31.27±0.75
女贞 Ligustrum lucidum	卵圆形 Ovate	1	223.03±14.36	24.72±0.15	17.25±1.13
		2	236.82±3.98	24.41±0.22	13.72±1.07
		3	303.50±6.90	22.93±1.24	12.59±0.05
石楠 Photinia serratifolia	圆形 Round	1	257.51±15.93	23.16±0.34	22.18±0.28
		2	425.36±14.36	20.23±0.32	17.35±1.30
		3	443.75±14.36	15.20±0.43	12.76±0.38
大叶黄杨 Buxus megistophylla	圆形 Round	1	236.82±22.17	37.91±0.28	37.26±0.23
		2	239.12±7.96	26.83±0.13	25.89±1.01
		3	287.40±14.36	24.74±0.15	24.26±0.27
海桐 Pittosporum tobira	圆形 Round	1	170.14±7.96	40.55±0.45	39.62±0.20
		2	271.31±21.07	32.50±0.43	30.92±0.06
		3	298.90±21.07	25.65±1.16	24.13±0.11
南天竹 Nandina domestica	卵圆形 Ovate	1	377.07±31.86	19.91±0.31	18.66±0.31
		2	434.55±11.95	19.48±0.35	17.39±0.47
		3	448.35±20.69	15.76±0.09	14.61±0.37

下载: 导出CSV 

| 显示表格

图  5  7种常绿树种在3处采样区的叶面气孔形态观察 （×1200倍）

Figure  5.  Leaf stomata morphology of evergreens in 3 sampling areas (1200 ×)

下载: 全尺寸图片 幻灯片

表  2  树种、污染程度与气孔指标的方差分析

Table  2.  Variance analysis on tree species, pollution degree, and stoma index

固定因子 Fixed factor	因变量 Dependent variable	Ⅲ类平方和 Sum of squares of class Ⅲ	自由度 Degrees of freedom	均方 Squared value	F	显著性 Significance
污染程度 The degree of pollution	气孔长度 Stomatal length	365.339	2	182.669	1.417	0.268
	气孔宽度 Stomata width	342.450	2	171.225	2.192	0.141
	气孔密度 Stomatal density	20557.829	2	10278.915	0.457	0.640
树种差异 Tree species differences	气孔长度 Stomatal length	2052.773	6	342.129	7.571	0.001**
	气孔宽度 Stomata width	1303.944	6	217.324	4141.690	0.001**
	气孔密度 Stomatal density	386552.852	6	64425.475	23.403	0.000**
注：**表示P＜0.01。 Note: ** means P＜0.01.

下载: 导出CSV 

| 显示表格

2.4   滞尘量与生理指标的相关性及主成分分析

为了进一步分析树种滞尘能力与生理指标之间的关系，分析各树种的单位叶面积颗粒物（TSP、PM₁₀、PM_2.5）滞留量与生理指标之间的相关性。结果表明，各指标间的相关系数达到极显著相关的有7个（表3）。

表  3  叶片滞尘量与生理指标之间的相关性

Table  3.  Correlation between dust retention and physiological indices of leaves

	X1	X2	X3	X4	X5	X6
X1	1	0.917**	0.886**	−0.692	−0.739	0.966**
X2		1	0.997**	−0.531	−0.508	0.907**
X3			1	−0.478	−0.443	0.878**
X4				1	0.918**	−0.709
X5					1	−0.727
X6						1
注：**表示在0.01水平上极显著相关；X1：单位叶面积TSP滞留量； X2：单位叶面积PM10滞留量； X3：单位叶面积PM2.5滞留量； X4 ：叶绿素a含量；X5：叶绿素b含量；X6：丙二醛含量。 Note: **extremely significant correlation at P＜ 0.01; X1: TSP retention per unit leaf area; X2: PM10 retention per unit leaf area; X3: PM2.5 retention per unit leaf area; X4: chlorophyll a content; X5: chlorophyll b content; X6: malondialdehyde content.

下载: 导出CSV 

| 显示表格

将滞尘、生理指标转化为6个主成分进行分析（表4）。由滞尘量在第1主成分上具有较高载荷，说明第1主成分主要反映了植物滞尘能力。第2主成分中叶绿素a、b含量具有较高载荷，说明第2主成分主要反映了树种的生理响应。前2个主成分的累积方差贡献率为96.65%，反映出滞尘能力和生理指标95%以上的信息；因此，可以提取前2个主成分作为植物滞尘能力和生理响应评价的综合指标。

表  4  滞尘量与生理指标的主成分贡献率和因子载荷矩阵

Table  4.  Contribution rate and factor loading matrix of principal components on dust retention and physiological indicators

主成分 Principal component	成分矩阵 Component matrix						特征值 Eigen value	贡献率 Contribution/%	累计贡献率 Cumulative/%
	X1	X2	X3	X4	X5	X6
1	0.98	0.92	0.89	−0.79	−0.80	0.97	4.80	79.92	79.92
2	0.09	0.38	0.44	0.56	0.58	0.08	1.00	16.73	96.65

下载: 导出CSV 

| 显示表格

2.5   树种抗污染性综合评价

为了更好评价供试植物的综合抗污染能力，根据提取的前2个主成分的贡献率得到树种的综合得分公式：

式中：F₁表示树种在第1主成分的得分；F₂表示树种在第2主成分的得分；F表示树种的综合得分。

基于以上公式对植物的滞尘、生理以及综合抗污染能力进行评价和排序（表5），得出滞尘能力排前2名的树种为雪松和圆柏，说明雪松和圆柏的滞尘能力较强，生理响应前2名的树种为石楠和女贞，综合得分前2名为雪松和圆柏。

表  5  树种在第1、2主成分的得分以及综合得分

Table  5.  Overall score and scores of first and second principal components on varied species of plants

树种 Tree species	F1	F2	F	排名 Ranking
雪松 Cedrus deodara	3.19	0.01	2.55	1
圆柏 Sabina chinensis	3.08	−0.31	2.41	2
女贞 Ligustrum lucidum	−0.25	0.73	−0.08	3
石楠 Photinia serratifolia	−1.50	1.28	−0.98	4
大叶黄杨 Buxus megistophylla	−1.41	0.64	−1.02	5
海桐 Pittosporum tobira	−1.32	−0.63	−1.16	6
南天竹 Nandina domestica	−1.79	−1.72	−1.72	7
注：F1：树种在第1主成分（滞尘能力）的得分，随数值增加树种滞尘能力得分高；F2：树种在第2主成分（生理响应）的得分，随数值增加树种生理响应得分高；F：树种的综合得分；表中的排名是按照F综合得分进行的排名。 Note: F1: The score of tree species in the first principal component (Dust detentions) , the dust retention ability score of tree species was higher with increasing value; F2: The score of tree species in the second principal component (Physiological response) , the physiological response score of tree species was higher with increasing value; F: Composite score of tree species; The rankings in the table are based on the F composite score.

下载: 导出CSV 

| 显示表格

3.   讨论

综合单位叶面积滞尘量和叶表面观察，发现，7种常绿树种中，针叶树的单位叶面积滞尘能力显著强于阔叶树种，这可能与针叶树本身存在蜡质微结构有关，这与史军娜等^[25]、 赵松婷等^[26]的研究结果一致。针叶树种中圆柏的单位叶面积TSP滞留能力强于雪松，且两者之间存在显著差异，圆柏的鳞形叶有密集的脊状突起，突起之间沟槽宽度与深度较明显，对粒径较大的颗粒物(TSP)滞留能力更显著^{[27, 28]}。雪松分泌的油脂对粒径较小的颗粒物(PM₁₀和PM_2.5)滞留能力更显著^[25-27]。结合单位叶面积滞尘量，阔叶树种中，女贞的滞尘能力最强，大叶黄杨、石楠和海桐叶表面较光滑，滞尘能力一般，南天竹滞尘能力最差，虽然南天竹气孔密度最大，但由于叶片表面气孔较小且平整，滞尘能力最弱^{[19, 29]}。

本研究中各树种叶片叶绿素a、b含量在污染严重的地区较低，并与滞尘量呈显著负相关（P＜0.01）。这可能是由于叶片受到大气污染后，叶片中的叶绿素a、b受到破坏而分解。裘璐函等^[30]的研究也表明随滞尘时间的增加, 大部分树种叶片叶绿素相对含量呈下降趋势。本研究中丙二醛含量与滞尘量呈显著正相关（P＜0.01），在污染严重的地区相对较高，这可能是由于污染物破坏了植物叶片内活性氧产生和清除的平衡，导致丙二醛含量升高^{[24, 31]}；此外， 俎丽红等^[32]的研究发现，树种的丙二醛含量升高还可能与大气颗粒物中重金属含量较高、供试树种年龄较大有关。污染严重的工厂的植物，气孔密度最大，长度和宽度均最小，滞尘能力较强。王会霞等^[13]、李艳梅等^[33]发现不同程度的颗粒物污染会导致植物气孔大小和密度发生变化。这可能是由于植物在遭受颗粒物污染时其气孔密度在一定阈值内增加有利于提高叶片单位叶面积的气孔交换能力，使其能适应污染严重的环境^[34]。而气孔密度与滞尘能力在整体阈值上影响趋势仍需进一步研究。不同树种在遭受颗粒物污染时生理指标响应程度不同，在严重污染的情况下，圆柏和雪松等树种的叶绿素a、b含量、气孔密度和气孔长度等生理指标的变化较小，说明二者对于颗粒物污染的抗性较强，这可能与针叶树种能分泌脂类物质提升自身的抑菌和抗污染能力有关^[27,32]。

通过主成分分析法得出树种综合抗污染能力的得分依次为雪松、圆柏、女贞、石楠、大叶黄杨、海桐和南天竹。可见，雪松和圆柏的综合抗污染能力较强，是因为二者滞尘能力强，同时在遭受污染时生理指标变化小即对颗粒物污染的抗性较强，因此二者拥有较强的综合抗污染能力。

4.   结论

（1）树种滞尘能力与叶表面结构（沟槽、中脉、毛状体和气孔等）密切相关，叶表面粗糙的树种滞尘能力更强。不同程度污染下树种的滞尘量不同，污染越严重的地区树种滞尘量越高。

（2）大部分树种的叶绿素a、b含量在污染严重地区相对较低，丙二醛含量在污染严重地区相对较高。树种的滞尘量与生理指标存在显著相关性，不同树种在遭受颗粒物污染时生理指标响应程度不同。

（3）综上所述，今后在郑州治理大气粉尘污染、进行城市绿化时，还需将滞尘能力、不同植物叶表结构进行综合考虑，应优先选用雪松、圆柏这种综合抗污染能力较强的树种从而提高城市环境质量。