稀土开采产生大量堆积尾砂,出现强烈水土流失过程,伴生严重水环境和地质灾害。采用工程和非工程措施结合手段开展稀土尾砂地植被修复是尾砂地灾害防治的常见手段。评价矿区植被恢复效果及成因可为尾砂治理措施的规划设计提供理论依据。
NDVI时程变化
一、引言
稀土开采产生大量堆积尾砂,出现强烈水土流失过程,伴生严重水环境和地质灾害。采用工程和非工程措施结合手段开展稀土尾砂地植被修复是尾砂地灾害防治的常见手段。评价矿区植被恢复效果及成因可为尾砂治理措施的规划设计提供理论依据。本文以江西省寻乌县离子型轻稀土尾砂区为研究区,收集研究区1989-2015年Landsat数据和寻乌站气象数据,实地调查研究区典型尾砂适生植被和堆积尾砂的地物波谱数据,解算典型植被的NDVI值,分析矿区NDVI的年际和年内时程变化规律。
二、基于Landsat数据的NDVI时程变化
在利用遥感技术对植被信息的研究中,NDVI是植被指数中运用最为广泛的一种,它能很好地反映了植被生长状态以及植被空间分布情况,与植被分布密度呈线性相关。许多研究表明NDVI与植被覆盖度、叶面积指数(LAI)和光合作用等植被参数有关,NDVI的时间变化曲线能有效地表现出植被的季节变化、年际变化和人为活动影响变化。因此,在监测废弃矿区植被恢复情况的过程中NDVI起到有效地作用。本文选取归一化植被指数(NDVI)来反映稀土矿尾砂地植被变化情况,其计算形式为近红外波段(NIR)与红光波段(R)的差值比上这两个波段数值之和,公式为:
表1Landsat TM/OLI NDVI 计算公式
研究区NDVI时程变化分析分为两个部分进行,第一部分为整体矿区NDVI时程变化分析,第二部分选取植被修复时间较早的“稀土尾砂试验区”双茶亭矿区分析NDVI时程变化规律。
2.1整体矿区NDVI时程变化
基于ENVI软件中BandMath 计算NDVI结果,将整体矿区NDVI值提取出来并导出成ASCII文件,方便之后分析使用。计算NDVI平均值,诊断NDVI时序趋势,揭示整体矿区NDVI时程变化规律,见图1。图1表明:1989-2015年整体矿区NDVI值在0.12-0.5间变化,平均值为0.31;NDVI最大值和最小值所在年份分别为1995年和2013年,分别为0.5和0.12,变化幅度为0.38;NDVI值相对时间的线性回归决定系数R2为0.79,取显著性水平为0.05,相应临界值R205.0 为0.08,决定系数大于相应临界值,线性回归趋势显著,NDVI值的倾向率为-0.08(10a)-1,表明NDVI值显著减小;NDVI值的Mann-kendall趋势分析检验值Z为-6.68,取显著性水平为0.05,相应临界值Z05.0 为1.96,检验值的绝对值Z大于相应临界值,Mann-kendall趋势检验显著,表明整体年NDVI值呈显著性下降趋势。
图1 整体矿区 NDVI 时程变化
NDVI时程变化可分为三个阶段,第一阶段为稳定期,时间为1989-2005年,此阶段NDVI值相对时间的线性回归决定系数R2为0.052,取显著性水平为0.05,相应临界值R205.0 为0.157,决定系数小于相应临界值,NDVI变化趋势不显著,平均NDVI为0.43,整体稀土矿区存在一定的开采,但是整体NDVI值未大幅度减小;第二阶段为剧烈开采期,时间为2006-2008年,此阶段NDVI值相对时间的线性回归决定系数R2为0.69,取显著性水平为0.05,相应临界值R205.0 为0.658,决定系数大于相应临界值,NDVI呈显著下降趋势,倾向率为-0.3(10a)-1,平均NDVI为0.29,较稳定期下降了0.14,此阶段稀土矿区进行剧烈的采矿活动,人类活动对稀土矿区影响较大,严重破坏了矿区植被,导致NDVI值迅速下降;第三阶段为植被修复期,时间为2009-2015年,此阶段NDVI值相对时间的线性回归决定系数R2为0.395,取显著性水平为0.05,相应临界值R205.0 为0.197,决定系数大于相应临界值,NDVI呈显著下降趋势,倾向率为-0.046(10a)-1,整个时间段NDVI呈现较低值,平均NDVI为0.19,较剧烈开采期下降了0.1,寻乌县2009年对稀土矿区进行植被修复治理,使研究区内的2009年和2010年的植被数量有所增加,对于整体矿区恢复效果一般。
2.2双茶亭矿区NDVI时程变化
在整体矿区NDVI影像数据中提取以种植桉树为主的双茶亭矿区NDVI,计算NDVI平均值,诊断NDVI时序趋势,揭示双茶亭矿区NDVI时程变化规律,见图2。按照整体矿区NDVI变化阶段对双茶亭矿区NDVI时程变化进行分析。
图2 双茶亭矿区 NDVI 时程变化
结果表明:1989-2015年双茶亭矿区NDVI值在0.17-0.65间变化,平均值为0.4;NDVI最大值和最小值所在年份分别为2015年和2008年,分别为0.65和0.17,变化幅度为0.48;稳定期NDVI值相对时间的线性回归决定系数R2为0.321,取显著性水平为0.05,相应临界值R205.0 为0.156,决定系数大于相应临界值,NDVI呈显著下降趋势,倾向率为-0.049(10a)-1,平均NDVI为0.44,稀土开采对双茶亭矿区内植被存在一定影响,植被呈现减少趋势;剧烈开采期NDVI值相对时间的线性回归决定系数R2为0.732,取显著性水平为0.05,相应临界值R205.0 为0.57,决定系数大于相应临界值,NDVI呈显著下降趋势,倾向率为-0.63(10a)-1,平均NDVI为0.28,较稳定期减少0.16,和整体矿区情况类似,剧烈稀土开采活动严重破坏双茶亭矿区内植被,导致植被迅速减少;植被修复期NDVI相对时间的线性回归决定系数R2为0.884,取显著性水平为0.05,相应临界值R205.0 为0.187,决定系数大于相应临界值,NDVI呈显著上升趋势,倾向率为0.25(10a)-1,平均NDVI为0.39,较剧烈开采期增加0.11,双茶亭植被修复效果较好,对比整体矿区可得,油茶配置模式植被修复治理效果较差,桉树配置模式效果较好。
三、基于地物波谱的典型地物NDVI特征值分析
基于研究区各典型地物的地物波谱特性,利用红光波段和近红外波段的平均反射率计算相应地物的NDVI,见图3。图3表明:在研究区典型地物中,尾砂相应NDVI值为0.03,芒萁为0.69,芭茅为0.53,芒草为0.52,油茶为0.69,木荷为0.8,竹子为0.79,马尾松为0.75,桉树为0.6,樟树为0.52,其中木荷NDVI值最大,尾砂NDVI值最小,尾砂存在一定的背景值。由于研究区地物波谱数采集时间为2019年1月与7月,分别属于冬季和夏季,则获取各地物波谱特性计算出的NDVI会有一定的波动范围,其中波动范围较大的为芒萁、芭茅和芒草,这三种植被属草本类植被,说明对于NDVI,草本类植被受季节变化影响较大,乔木类植被对相对较小,即草本NDVI在年内变化更为剧烈,乔木类植被NDVI年内较为稳定。
图3典型地物NDVI
四、基于小波分解的NDVI时程变化
4.1小波分析理论概述
傅里叶变换是最常用的分析信号工具,用于推导不同幅度和频率下特定正弦波函数的组合。但是傅里叶变换存在不足之处,在分析信号特征时会忽略信号的时间信息,通过变换结果无法判断某个信号发生的时间。以傅里叶变换为数学基础,理论物理学家A·Grossman建立了小波变换的理论体系。小波变换分析信号时可通过小波基函数的变换分析信号的局部特征,信号分析结果具有比傅里叶变换更良好的局部特性。分解后的信号呈现低频和高频特性两部分内容,低频信号具有较高的频率分辨率和较低的时间分辨率,高频信号具有较高的时间分辨率和较低的频率分辨率,所以有着“数学显微镜”之称。正是小波变换拥有着这些特性,使小波变换在信号处理、数值计算和流体力学等领域得到了广泛的应用。
4.2基于小波分解的NDVI多尺度时程变化
为了更好研究NDVI的时程变化规律,基于Matlab软件,本文将NDVI变化曲线作为一种信号曲线,对其进行小波分解,分析其低频信号与高频信号的特点。根据NDVI数据特点,采用db3小波函数对整体矿区NDVI和双茶亭NDVI进行分解,结果如图4和图5所示。
图4 整体矿区 NDVI 信号和小波分解部分注:(a)表示原始信号,(b)表示低频信号,(c)表示高频信号。
(1)低频信号按照NDVI变化阶段对低频信号进行分析。图4表明:稳定期整体矿区低频NDVI值处于较高水平,进入剧烈开采期后低频NDVI值迅速下降,在植被修复期保持稳定;低频信号整体趋势与整体矿区NDVI时程变化趋势一致,说明低频信号表示NDVI的年际变化特征,将年内波动信息滤除,呈现与原始信号一致的三个变化阶段。图5表明:双茶亭矿区低频NDVI值呈现微下降-剧烈下降-上升的三个阶段变化,与双茶亭矿区NDVI时程变化趋势一致。
(2)高频信号高频信号表示为NDVI的年内变化,按照NDVI变化阶段对高频信号进行分析。图4表明:稳定期整体矿区高频NDVI值呈波动变化,波动幅度较大,说明在此阶段整体矿区年内NDVI变化较为剧烈;剧烈开采期整体矿区高频NDVI值也呈波动变化,但是波动幅度较小,说明在此阶段矿区年内NDVI呈波动变化,但是由于人类活动较为剧烈,使植被数量迅速较少,导致NDVI变化幅度变小;植被修复期整体矿区高频NDVI值在此阶段初期呈水平变化,后期呈波动变化,说明研究区年内NDVI在植被修复期初期年内变化较小,这是由于剧烈开采期过后植被破坏过于严重,导致矿区内植被数量稀少,年内植被变化较小,后期由于植被修复措施的实施,植被得到一定的恢复,然后呈现出了年内波动变化。图5表明:双茶亭矿区高频NDVI值三个时期都呈波动变化,这是由于在剧烈开采后,植被措施的迅速实施,使双茶亭矿区内植被得到良好的补充,导致在剧烈开采期之后、植被恢复期初植被表现出良好的年内波动变化。
图5双茶亭矿区NDVI信号和小波分解部分注:(a)表示原始信号,(b)表示低频信号,(c)表示高频信号。
结合草本植被NDVI在年内变化更为剧烈,乔木类植被NDVI年内较为稳定的结论可得:稳定期与剧烈开采期NDVI波动较大是由于草本植被变化剧烈引起的,说明矿区植被变化的主要原因为草本植被的变化;在植被修复期初期高频NDVI值呈水平变化特征,说明此时矿区内草本类植被数量稀少,原因为剧烈稀土开采将草本植被破坏严重,几乎不存在草本类植被,在经过人工植被修复后,矿区植被得到一定的补充,草本植被开始增长,最后呈现波动特征。总体表明:矿区内对NDVI影响最大的为草本类植被,基于此情况,矿区植被恢复措施的植被配置模式可针对草本类植被进行选择,并针对草本植被进行维养。
四、讨论
本章基于LandsatTM/OLI 数据和研究区典型地物波谱数据,计算了研究区1989-2015年的NDVI值以及典型地物的NDVI值,总体分为三个部分:
(1)整体矿区NDVI时程变化和双茶亭矿区NDVI时程变化分析。整体矿区NDVI呈下降趋势,可分为三个阶段,即稳定期、剧烈开采期和植被修复期,其中整体矿区NDVI在稳定期呈波动变化,在剧烈开采期和植被修复期呈显著下降趋势,双茶亭矿区NDVI在稳定期和剧烈开采期呈显著下降趋势,在植被修复期呈显著上升趋势。油茶配置模式尾砂修复治理效果较差,桉树配置模式效果较好。
(2)典型地物NDVI值的分析。典型地物可分为三种稀土尾砂、草本植被和乔木类植被,其中尾砂和乔木类植被年内NDVI值变化幅度较小,草本植被年内NDVI值变化较大。
(3)对NDVI时程变化曲线进行小波分解。结果分解为低频信号和高频信号,且分别表示NDVI的年际变化和年内变化特征。通过分析整体矿区NDVI高频信号发现了对研究区NDVI产生影响的主要因素为草本类植被,此结论可指导后期继续进行植被修复,针对草本植被,选择更良好的植被配置模式
推荐:
地物光谱仪iSpecField-HH/NIR/WNIR
地物光谱仪是莱森光学专门用于野外遥感测量、土壤环境、矿物地质勘探等领域的最新明星产品,独有的光路设计,噪声校准技术、可以实时自动校准暗电流,采用了固定全息光栅一次性分光,测试速度快,最短积分时间最短可达20μs,操作灵活、便携方便、光谱测试速度快、光谱数据准确,广泛应用于遥感测量、农作物监测、森林研究、海洋学研究和矿物勘察等各领域。
审核编辑 黄宇
-
地物光谱仪
+关注
关注
0文章
75浏览量
3241
发布评论请先 登录
相关推荐
评论