新疆天山天池景区表土磁性特征的时空变化①

常 洁，陈学刚,2*，董 煜,2

(1 新疆师范大学地理科学与旅游学院，乌鲁木齐 830054；2 新疆师范大学新疆维吾尔自治区干旱区湖泊环境与资源重点实验室，乌鲁木齐 830054)

旅游过程中旅游景区是游客游览的核心要素。但是，随着景区内游客活动量的增加，在消耗景区内外资源的同时，景区污染物的排放也显著增加。景区内土地践踏、旅游设施建设、尾气排放、旅游垃圾等一系列的旅游活动导致的生态污染直接影响到人们的健康[1]。土壤作为旅游景区生态圈中游客直接接触的部分，在维持生态系统平衡与稳定的同时也极易受到旅游活动的干扰。而景区内表土对旅游活动的敏感性响应会通过土壤理化性质的改变表现出来[2]，所以对景区内土壤污染的研究是评价景区环境质量的重要手段[3]。土壤磁性作为土壤性质的一个重要组成部分，能够表征人为活动影响下的表土磁性特征与自然成土下的表土磁性特征的差异性[4-5]。如，在景区自助游带动下，私家车和区间车增多等一系列能源消耗与排放的旅游活动中产生的污染物通常含磁性矿物颗粒，而这些矿物颗粒主要是与交通活动如尾气排放、轮胎磨损、垃圾焚烧等有关的Pb、Cu、Zn的污染来源[6]。研究显示，污染物中磁学参数值与重金属含量具有显著正相关性，这使得磁学参数不仅成为重金属元素含量的代用指标同时也能指示人为污染[7]。由于磁学参数能够灵敏地判断样品中磁性矿物含量、矿物种类和粒径大小，且通过磁学参数分析可推断磁性颗粒受时空变化影响的过程，因而其也成为环境污染变化评估的代用指标[8-9]。磁学参数测量方法有迅速、灵敏、对环境无破坏性、费用低等优势，因而被广泛应用于土壤污染源的测定和追踪、污染评价等方面[10]。当前，国内外学者针对土壤性质的时空变化分析主要集中在土壤的养分[11-12]、重金属[13]和盐分[14]

等，而针对土壤磁性的研究主要在单一时间点的分析并仅见于农田土壤磁化率空间变异研究[15]，关于景区内土壤磁性时空变异的报道尚不多见。本文在前人研究的基础上对天池景区表土磁性随时间变化引起的空间分布变化进行研究，并绘制空间分布图，以期能得到天池景区表土磁性污染的时空变异特征，明晰表土的污染特点、污染程度、磁性分布规律等。

1 材料与方法

1.1 研究区概况

天山天池风景区位于天山东段博格达主峰北坡的新疆阜康市南，距阜康市区37 km，距首府97 km，地理坐标 88°00′ ~ 88°20′E，43°45′ ~ 43°59′N，面积248 km2，海拔1 910 m。研究区气候属典型的大陆性干旱气候，年平均气温2.55 ℃，年降水量为500 ~600 mm；区内土壤类型以灰褐色森林土、山地栗钙土、山地棕钙土为主；地质构造方面，天山天池属东天山博格达造山带，呈东–西走向夹持于准噶尔和吐鲁番–哈密两个中、新生代盆地之间，处于博格达晚古生代弧后断陷盆地与准噶尔中央地块的接合部，断层、褶皱发育，由一系列推覆构造岩片组成[16]。丰富的自然景观和人文资源，不仅为天池景区旅游活动的开展提供了良好的环境基础，且使其在新疆旅游业发展中居于重要地位，其于2007年被评为国家5A级旅游景区。

1.2 样品采集和分析方法

2015年8月和2017年5月分别对天山天池景区内表土进行样品采集，采样点的选取主要以旅游交通区、旅游活动频繁区和旅游活动稀少的未干扰区为主。采用以100 m间隔的系统采样法均匀地在景区内进行布点采样，分别采集79个(2015年)、35个(2017年)表土样品，其中2017年采样点与2015年采样点重合。每个样品利用梅花点采样法采集表土0 ~ 10 cm，每个样品混合而成1 kg左右的样品，装入双层聚乙烯塑料密封袋中，标注样品信息及编号。

采集的表土样品经过放置在阴凉、干燥、通风、无尘土污染的室内进行自然风干后，拣出杂质、植物根系、大块石砾等，过1 mm的尼龙筛，称重并用塑料保鲜膜包紧后装入10 cm3的测磁专用样品盒内压实封口，供测试使用。样品送至中国地质大学古地磁实验检测室进行测试，样品的低频磁化率χLF(976Hz)和高频磁化率χHF(15 616 Hz)采用卡帕桥MFK1-FA多频磁化率仪测量，根据测量结果计算频率磁化率χFD，χFD(%)=(χLF–χHF)/χLF×100。表土样品的等温剩磁(IRM)和饱和等温剩磁(SIRM)利用IM-10-30脉冲磁化仪和旋转磁力仪测量，SIRM外加磁场强度为1 000 mT；并测量样品的–20、–100、–300 mT反向磁场的等温剩磁。软剩磁SOFT=(SIRM – IRM–20mT)/(2·m)(其中m代表质量)和硬剩磁HIRM=(SIRM +IRM–300mT)/(2·m)，SIRM/χLF以及比值参数通过已测磁性参数计算得到，Sratio(%)= –IRM–300mT/SIRM×100。

1.3 数据处理

样品磁性参数的统计分析利用SPSS 20.0 软件完成，磁学参数的统计图和空间分布图利用 Graph Prism 6和ArcGIS 10. 0软件完成。

2 结果

2.1 磁性矿物类型参数

研究区两次样品的磁性参数描述性统计见表1。SOFT是反映低矫顽力的亚铁磁性矿物(如磁铁矿或磁赤铁矿)含量的指标[17-19]，HIRM则用于估计样品中高矫顽力的不完全反铁磁性矿物(如赤铁矿和针铁矿)含量[20]。从表1可以看出，研究区土壤SOFT的变化范围2015年为(134.62 ~ 1 063.04)×10–5Am2/kg(平均值 283.29×10–5Am2/kg)，2017年为(108.98 ~831.57)×10–5Am2/kg(平均值 374.15×10–5Am2/kg)；HIRM的变化范围2015年为(17.86 ~ 110.58)×10–5Am2/kg(平均值 43.96×10–5Am2/kg)，2017年为(11.48 ~ 144.39)×10–5Am2/kg(平均值 60.58×10–5Am2/kg)。Sratio可用来评价低矫顽力软剩磁与高矫顽力硬剩磁矿物的相对含量[21]，其值越接近1时，表示低矫顽力的亚铁磁性矿物(如磁铁矿或磁赤铁矿)是样品中主要的磁性矿物[22]。当Sratio<0.8，则说明样品中含有较高比例的硬磁性矿物，反之，0.8

2.2 磁性矿物含量参数

χLF、SIRM与样品中磁性矿物含量呈正比。χLF大体反映环境样品中磁性矿物(包括抗磁性、顺磁性、不完全反铁磁性和亚铁磁性矿物)的含量[24]。SIRM与χLF的区别在于其不受顺磁性和抗磁性物质的影响，主要表明亚铁磁性和不完全反铁磁性矿物的贡献[25]。天池景区表土的χLF、SIRM值变化范围2015年 分 别 为(46.05 ~ 210.32)×10–8m3/kg(平 均 值95.25×10–8m3/kg)、(516.56 ~ 4 219.12)×10–5Am2/kg(平均值 1 110.86×10–5Am2/kg)；2017年分别为(43.20 ~326.54)×10–8m3/kg(平 均 值 113.04×10–8m3/kg)、(454.24 ~ 4 241.40)×10–5Am2/kg(平均值 1 574.40×10–5Am2/kg)。

2.3 磁性矿物粒度参数

如果样品中磁铁矿主导亚铁磁性矿物，SIRM/χLF可用于估计磁性矿物粒度，其值越高，表明磁性矿物颗粒越粗，其低值指示样品以多畴和超顺磁性颗粒(SP)为主，低值到中值时以假单畴颗粒主导，其高值说明稳定单畴颗粒占优[26]。天池景区表土SIRM/χLF的平均值2015年、2017年分别为11.62×103、13.83×103A/m，属于低值到中值过渡范围；同时，2017年平均值高于2015年，表明粗颗粒磁性矿物增加。χFD反映细小的SP对磁化率的贡献，值越高表明整体磁性颗粒越细；χFD值的大小也指示SP比例的多少，同时能够辨别土壤颗粒是自然形成还是人为活动所致[27]。研究区土壤χFD变化范围2015年为1.54% ~11.08%(平均值 7.59%)，2017年为1.53% ~ 13.71%(平均值 5.75%)。Dearing[28]指出，若χFD<2%，样品中基本无SP的存在；2%<χFD<11%，样品中SP含量介于10% ~ 50%；11%<χFD<14%，样品中SP含量>75%，样品中以SP颗粒为主。研究区2015年样品中χFD% 值在8% ~ 11% 的样品占比较高，而2017年样品中χFD值在2% ~ 8% 的样品占比较高，说明2017年样品中粗颗粒比例明显增多，两次样品中多数样品的磁性载体含有10% ~ 50% 不等的SP颗粒。

表1 不同时期表土样品磁性参数统计结果Table 1 Statistics of magnetic parameters of topsoil samples in different periods

3 讨论

3.1 天池景区表土磁性特征

不同磁参数的双变量图是一种有用、有效的定性图形技术，无需昂贵的化学分析，也可以快速、简便地鉴别某些磁性矿物特征和畴态[26]。当磁性矿物由单一磁性矿物构成或磁性矿物的混合物以恒定比率控制时，χLF和SIRM呈线性关系。天池景区两次表土样品的χLF与SIRM变化趋势相似且呈显著正相关(r=0.75，0.74；图2)，反映了样品中磁性矿物主要是亚铁磁性矿物质和不完全反铁磁性矿物，同时较高的相关性说明表土样品中磁性矿物来源的多样性；χLF与SOFT呈极显著正相关(r=0.90，0.89)，与HIRM呈中低度相关性(r=0.45，0.64)，反映表土样品的磁学性质主要受低矫顽力亚铁磁性矿物所主导，且伴有一定不完全反磁性矿物；SOFT平均值远高于HIRM，也进一步证明了上述结果；Sratio值大于90%，表示低矫顽力的亚铁磁性矿物(如磁铁矿或磁赤铁矿)是样品中主要的磁性矿物；从SIRM/χLF结果看，两次样品的载磁矿物以亚铁磁性矿物(钛磁铁矿)为主。

天池景区2017年样品χLF的平均值高于2015年，且变化幅度大于2015年。与其他研究区表土样品相比，该景区表土样品χLF平均值略低于喀纳斯景区表土样品(平均值 122.30×10–8m3/kg)[3]、银川市城市表土样品(平均值 117×10–8m3/kg)[29]。考虑近两年自然因素的不变性，因此，2017年样品磁性的增强主要是由旅游活动产生的生活排污、汽车尾气、化石燃料燃烧等人为因素造成。

已有研究表明，自然成土土壤的χLF与χFD具有良好的正相关性，而污染土壤的相关性一般表现较高的χLF和较低的χFD[7]。2015年样品χFD与χLF无相关性(r=0.07)，2017年样品呈低相关性(r= –0.21)，表明景区表土样品磁性未受SP颗粒主导，2017年样品中粗颗粒明显增加，磁性增强。人为磁性矿物的磁性不同于自然产生的磁性矿物，因为它们具有较粗的磁性颗粒尺寸(多畴和假单畴)[30]。综合来看，天池景区表土样品中存在人为形成的假单畴粗颗粒和一定量的SP细颗粒。

3.2 天池景区表土磁性特征的时间变化

不同时间段表土样品表现出的磁性特征差异，在一定程度上可以追踪当地土壤污染的时间变化。以天池景区表土磁性参数随时间变化的关系作图(图3)，且以1% ~ 99%、25% ~ 75% 的概率作为时间变化下磁性参数变化分布范围，比较两个时间段的各项参数值，由图3可知，景区2017年表土样品的χLF、SIRM、SOFT和HIRM值明显高于2015年。研究区远离城市和工业区，受工矿污染排放影响很小，磁性矿物含量的增强可能与近两年景区旅游活动排放的污染物累积有关，如景区内汽车尾气排放及轮胎磨损、筑路材料等磁性矿物颗粒，这与其他学者研究结论类似[31]。比较两次样品的Sratio值和χFD，如图3所示，2017年样品的Sratio高于2015年，同时2017年样品的χFD整体又小于2015年，这表明虽然两组样品的磁性颗粒都以粗颗粒为主，但2017年样品的磁性颗粒矫顽力更低，磁性颗粒更粗。比较两次景区表土样品的SIRM/χLF平均值，可见，随时间的进展磁性矿物颗粒粒径增加。

3.3 天池景区表土磁性特征的空间变化

研究区2017年样品的χLF、χFD、SIRM、SOFT和HIRM呈现出较高的CV值，均大于0.4，除了SOFT，其他变化程度都高于2015年样品，两次样品Sratio的CV值最低，反映出2017年表土样品的磁性参数变化幅度大，空间分布不均匀。研究区两次表土样品磁性矿物含量空间变化明显，在磁性矿物含量的高值区可能存在土壤污染，为比较两个时期景区表土样品磁性矿物含量的空间变化规律，圈定受游客旅游活动影响的可能污染范围和程度，对两次样品的χLF、SIRM、SOFT和HIRM值，采用反距离插值法进行空间插值，绘制出两次样品4种磁性参数的空间分布图。从图4可以看出，无论是2015年还是2017年样品的χLF值，高值都主要出现在天池北部的西王母祖庙、定海神针几个核心景点区域，并随着样品离核心景点区域距离的增加，χLF值呈逐渐下降的趋势；在西北部停车场附近及天池湖南部周边，两次采样点的χLF值均出现了相似的变化趋势。2017年的χLF值相对于2015年增长趋势明显增大，范围变广，尤其是在天池北部的核心区域、天池湖三面以及天池中部的民族文化村，2017年高值点是2015年高值点的1.5倍，且污染面积明显扩大，这可能是两年内天池景区频繁的旅游活动造成。可见，旅游活动对景区表土样品的磁性矿物含量的贡献很大。两次采样的SIRM、SOFT空间变化与χLF变化趋势基本相似，亚铁磁性矿物含量也呈现出沿环湖区高值并向外围逐渐递减的分布趋势。2017年的SIRM高值点是2015年高值点的1.8倍，2017年的SOFT高值点是2015年高值点的1.3倍，与χLF相似，天池景区2017年西南部及环湖区域的SIRM与SOFT值明显高于2015年。两次采样HIRM呈现从环湖区域值高并随着远离核心景点区域向外围减少的趋势，湖的北部值高，说明不完全反铁磁性矿物颗粒较多。对比两次采样的HIRM，2017年的高值点小于2015年的高值点，但整体上2017年大于2015年；同时出现了湖的北部核心景点区域及环湖的值明显增大。

综上所述，4种磁性参数的空间分布状况均表现出不连续分布的高值区域，两次采样点均显示湖的北部游艇停靠区、停车场及周边、核心景点区域值较高，2017年样品的磁性参数高值区与2015年高值区相似且在2015年高值区基础上范围显著扩大。可能原因是研究区北部分布着3个核心景点(定海神针的游客密集区、西王母神庙和飞龙潭)，景区内周围分布着道路、乘船码头。该区域表土频繁有尾气的排放及游客践踏、旅游垃圾等人为活动影响。环湖的西北部、西部和西南部土壤磁性参数富集程度仅次于环湖北部区域，可能是由于环湖周围停靠在码头边的游艇尾气累积及游客增多导致交通区间车增加带来的影响。湖的西南部和西北部不连续分布的高值点的增高是由于停车场周围汽车的增多导致尾气排放量增大导致，西部区域富集增大是由于靠近福寿观和民俗文化区受到了旅游餐饮和住宿区燃煤飞灰的影响。

4 结论

1)天山天池景区表土样品中磁性矿物以低矫顽力亚铁磁性矿物(磁铁矿和磁赤铁矿)为主，同时伴有赤铁矿和针铁矿组成的不完全反铁磁性矿物；样品磁性矿物含量适中，低于喀纳斯景区和银川市区；两次样品中磁性矿物主要由假单畴粗颗粒，同时含有相对较低的SP颗粒混合构成。

2)天山天池景区表土样品中磁性矿物特征随时间变化显著，2017年磁性矿物含量、粗颗粒比例高于2015年；除自然因素外，由于景区远离城市和工业区，两年的旅游活动产生的含有磁性矿物的污染物不断的累积，因而造成两次样品磁性矿物含量和粒度的差异。

3)天山天池景区表土样品的磁性参数空间分布变化明显，两期的磁性参数空间分布图呈现出在受人为活动影响较大的停车场及交通要道等旅游活动密集区的磁性矿物含量明显高于其他地区的格局。研究发现，监测景区表土磁性参数的时空变化有助于圈定可能污染范围的变化和判别污染物来源。