-
植被是陆地生态系统的主体,受地形、气候、土壤等多种环境因子的制约[1-3],同时,植被也是连接大气、土壤和水分的“纽带”,具有减少雨滴激溅,减缓地表径流,提高土壤保土固土功能,是控制土壤侵蚀与水土流失的重要因子。植被覆盖度是指单位面积内植被地上部分(包括叶、茎、枝)在地面的垂直投影面积百分比[4-5],是地表植被生长趋势的重要参量,也是区域生态退化监测最重要、最敏感的指标。获取地表植被覆盖度及其变化信息,对揭示地表植被时空变化规律,探讨植被覆盖变化的驱动因子,评价区域生态环境具有重要意义。植被变化与地理环境和人类活动密切相关,其中地形因子对植被生长具有重要影响。研究表明:岷江上游流域低植被覆盖度区域受坡度影响显著,而极高植被覆盖度区域受高程影响显著[6];汶川地震灾区,植被覆盖度受坡度和高程影响较大,受坡向影响较弱[7];赤水河流域,坡度和坡向对植被覆盖的影响都明显高于高程[3];王毅等[8]研究发现:贵州省普定县喀斯特地貌区植被覆盖度的分布与地形因子的空间关系具有不确定性,地形因素主要通过影响人类活动间接影响植被覆盖度。可见,地形因子对植被覆盖度变化的影响具有复杂性,不同地区影响植被覆盖度的关键地形因子不同。福建省福州市为东南沿海重要的经济中心,随着福州市经济的快速发展,其生态环境也经历了一个漫长而曲折的过程,城市生态环境及植被覆盖情况亟待提高。近年来,关于福州市植被覆盖度时空特征的研究颇多。例如,马雪莹[9]研究发现:福州市1994-2015年植被呈恢复趋势,但各辖区植被恢复速度差异较大。然而,李晶[10]表明:福州市1989-2014年植被呈退化趋势,且地表温度是植被覆盖度的重要影响因素。目前,针对福州市植被覆盖度时空变化特征与高程、坡度、坡向等地形因子关系的研究未见报道。因此,本研究以福州市的Landsat TM/OLI遥感影像为基础,采用改进的像元二分模型,探究福州市植被覆盖度时空变化特征,并分析植被覆盖度与地形因子的定量关系,以期为福州市城市植被恢复和生态环境建设提供依据。
HTML
-
福州市为福建省省会,位于福建省东部,闽江下游及沿海地区。福州地形属典型的河口盆地,福州市区位于盆地中央,四周被群山环绕,地势自西向东倾斜。福州地区跨中亚热带和南亚热带,气候类型为亚热带海洋性季风气候,年均气温为19.6 ℃,年均降水量为1 343.7 mm。本研究区为福州市5区,涵盖鼓楼区、台江区、仓山区、晋安区和马尾区,不包含长乐区,范围25°56′~26°25′N,119°10′~119°40′E,区域面积1 043 km2。
-
本研究所采用的遥感数据为Landsat 5 TM以及Landsat 8 OLI影像,原始影像数据源自地理空间数据云(http://www.gscloud.cn),影像轨道号均为119/042,空间分辨率为30 m。遥感影像的选择有2个方面要求:首先,云含量多少极易影响提取结果的精确度,应选取云含量较低的遥感影像反演植被覆盖度;其次,为减弱因植被的生长发育所引起的误差,精确地提取植被覆盖度,遥感影像需选取植被发育情况相近的时间段。因此,本研究所选取的福州市1995-2015年遥感影像均为秋季(9-11月)且云量低于5%的遥感影像。遥感数据投影方式均转换为基于WGS-84坐标系的通用横轴墨卡托投影,并进行几何校正、辐射定标和大气校正等预处理;然后,采用相同投影的矢量边界数据裁剪得到研究区域的遥感影像数据。遥感数据源信息如表 1。
影像 年份 传感器 波段 日期(月-日) Landsat 5 1995 TM 1~7 11-23 2002 11-10 2008 11-10 Landsat 8 2010 OLI 1~10 10-31 2013 10-23 2016 09-27 Table 1. Data sources of remote sensing for inversion of fractional vegetation coverage in Fuzhou City from 1995 to 2016
-
地形因子包括高程、坡度和坡向,30 m地面分辨率ASTER GDEM V2高程数据源自地理空间数据云(http://www.gscloud.cn),坡度及坡向数据主要基于数字高程模型(DEM)数据通过Arc GIS 10.2软件3D Analyst模块Aspect及Slope工具生成。
-
福州市植被覆盖度采用改进的像元二分模型来估算,其公式为[11-13]:
式(1)中:Fc为植被覆盖度;INDV为任意像元的归一化植被指数;INDVsoil为纯裸地的植被指数;INDVveg为纯植被覆盖的植被指数。归一化植被指数INDV通过式(2)计算[14-16]:
式(2)中:ρNIR为近红外波段的反射率;ρR为红光波段的反射率。
根据影像的INDV值序列,采用5%和95%置信度截取INDV的上下阈值,在置信区间内取最大值和最小值分别代替INDVveg和INDVsoil。利用ENVI 5.3软件估算福州市6期植被覆盖度,并依据SL190-2008《土壤侵蚀分类分级标准》[17],利用Arc GIS 10.2软件栅格重分类工具将研究区植被覆盖度分为5个等级:Ⅰ级为裸地,Ⅱ级为低覆盖度,Ⅲ级为中覆盖度,Ⅳ级为高覆盖度,Ⅴ级为极高覆盖度。依据研究区植被覆盖度总值与像元总数进行计算植被覆盖度均值:
式(3)中:n为像元总数;Fcm为第m个像元的植被覆盖度;$\overline{F_{\rm c}}$为植被覆盖度均值。
-
通过对2期植被覆盖度栅格图像数据进行逐像元做差值运算,并将结果进行标准化分级,来衡量植被动态变化的趋势、大小及空间特征。其计算公式[18]:
式(4)中:ΔFc为2期栅格图像数据的差值;Fca和Fcb为第a期和第b期栅格图像数据。基于1995和2015年的植被覆盖度栅格影像做差值运算,分析植被衰退恢复演变特征,并将植被覆盖度分成7个衰退恢复等级:严重退化(-100%,-30%]、中度退化(-30%,-15%]、轻微退化(-15%,-5%]、稳定(-5%,5%]、轻微改善(5%,15%]、中度改善(15%,30%]、极度改善(30%,100%]。
-
利用Arc GIS 10.2软件将6期植被覆盖度图转为点数据,将地形因子栅格数据提取至植被覆盖度点数据,获得各高程梯度、坡度梯度及坡向植被覆盖度均值,最后利用Excel 2016软件获得植被覆盖度受地形因子影响的动态变化特征曲线。
1.1. 研究区概况
1.2. 数据来源及处理
1.2.1. 植被覆盖度遥感数据
1.2.2. 地形数据
1.3. 研究方法
1.3.1. 基于像元二分模型的植被覆盖度计算
1.3.2. 基于差值法的植被覆盖度动态特征及植被恢复等级提取
1.3.3. 地形因子与植被覆盖度关系
-
1995-2015年,福州市植被覆盖度均值由57.93%增加至66.44%。1995-2002和2013-2015年,福州市植被覆盖度增速较快,平均每年提高0.88%和1.13%。然而,2010-2013年福州市植被覆盖度呈下降趋势,降低了-0.19%。福州市植被以极高覆盖度植被为主,面积565 km2,约占整个福州市面积的54.18%,其余各等级覆盖度植被面积差异不大(表 2)。裸地(Ⅰ级)和低覆盖度(Ⅱ级)植被集中分布在闽江、乌龙江沿岸以及鼓楼区、台江区等人口聚居区,高覆盖度(Ⅳ级)和极高覆盖度(Ⅴ级)植被则主要分布在北部和东部低山丘陵区域。
等级 Fc 占比/% 面积/km2 Ⅰ级(裸地) 0 10.30 107 Ⅱ级(低覆盖度) (0,25%] 11.78 123 Ⅲ级(中覆盖度) (25%,50%] 10.62 111 Ⅳ级(高覆盖度) (50%,75%] 13.11 137 Ⅴ级(极高覆盖度) (75%,100%] 54.18 565 Table 2. Classification of fractional vegetation coverage in Fuzhou City from 1995 to 2015 and area of fractional vegetation coverage
1995-2015年,福州市各时期不同等级植被覆盖度面积变化规律存在差异性(图 1,表 3)。Ⅰ级覆盖度植被面积呈波浪式增加,但1995-2002和2008-2010年表现为衰减特征。Ⅱ级和Ⅲ级覆盖度植被面积呈衰减特征,但2008-2013年出现增加趋势。Ⅳ级覆盖度植被面积呈持续衰减趋势,衰减幅度最大。虽然,2008-2013年Ⅴ级覆盖度植被面积表现出小幅度衰减,但1995-2008和2013-2015年均为大幅增加。整体而言,1995和2015年比较发现,福州市Ⅰ级和Ⅴ级覆盖度植被面积增加,分别增加13.33和194.86 km2,而Ⅱ级、Ⅲ级和Ⅳ级覆盖度植被面积分别减少12.38、59.09、136.71 km2。
Figure 1. Characteristics of fractional vegetation coverage at different levels in the period from 1995 to 2015 in Fuzhou City
Fc分级 1995-2002年 2002-2008年 2008-2010年 2010-2013年 2013-2015年 变化量/km2 变化率/% 变化量/km2 变化率/% 变化量/km2 变化率/% 变化量/km2 变化率/% 变化量/km2 变化率/% Ⅰ -1.40 -1.34 2.88 3.03 -4.47 -3.98 2.60 2.51 3.50 3.15 Ⅱ -0.85 -0.69 0.70 0.59 5.44 4.48 0.14 0.11 -10.98 -8.27 Ⅲ -7.05 -4.53 -1.90 -1.78 4.44 4.68 1.16 1.12 -5.37 -5.01 Ⅳ -10.21 -4.44 -6.49 -4.09 -4.63 -3.87 -0.82 -0.74 -7.30 -6.76 Ⅴ 19.51 4.54 4.81 0.85 -0.79 -0.13 -3.08 -0.52 20.15 3.45 说明:院负值表示不同等级覆盖度植被变化呈退化趋势 Table 3. Annual change of different fractional vegetation coverage during the different period from 1995 to 2015 in Fuzhou City
-
将1995与2015年植被覆盖度栅格数据利用差值法作差值运算,获得福州市区植被覆盖度衰退与恢复的面积信息,并根据植被覆盖度变化量,将植被覆盖度变化分为严重、中度和轻微退化、稳定、轻微、中度和极度改善7个等级(表 4)。
退化与改善等级分类 覆盖度/% 面积/km2 占比/% 严重退化 (-100,-30] 61 5.86 中度退化 (-30,-15] 50 4.78 轻微退化 (-15,-5] 66 6.34 稳定 (-5,5] 285 27.29 轻微改善 (5,15] 231 22.18 中度改善 (15,30] 197 18.93 极度改善 (30,100] 152 14.63 Table 4. Vegetation decline and recovery area and its proportion in Fuzhou City from 1995 to 2015
福州市植被主要以改善为主,轻微改善、中度改善、极度改善面积占比分别为22.18%、18.93%和14.63%,整个研究区域改善区域面积为581 km2,约占福州市面积的55.73%(图 2,表 4)。1995-2015年,福州市有些区域的植被覆盖度呈衰退趋势,面积为177 km2,约占研究区总面积的16.98%。从空间分布而言,福州市植被覆盖度改善区域主要位于北部和中部低山丘陵区域,而退化区域多分布于南部的闽江及乌龙江沿岸区域。
-
本研究依据研究区高程,以100 m为间隔将高程划分为7个梯度。依据GB/T 15772-2008《水土保持综合治理规划通则》[19]和水土流失调查中采用的8°作为缓坡和斜坡界线的方法,将坡度分为平坡(0°~5°)、缓坡(5°~8°)、斜坡(8°~15°)、陡坡(15°~25°)、急坡(25°~35°)和险坡(>35°)6个等级[20]。考虑到坡向对植被长势的影响,将研究区划分为平地(0°),阴坡(0°~67.5°,337.5°~360.0°),半阴坡(67.5°~112.5°,292.5°~337.5°),半阳坡(112.5°~157.5°,247.5°~292.5°),阳坡(157.5°~247.5°)5种坡向类型[20]。福州市地形因子空间分布如图 3。
Figure 3. Characteristics of elevations (A), and slope grades (B) and slope directions (C) in Fuzhou City
不同地形条件下植被覆盖度存在显著差异。随着高程值增加,植被覆盖度呈逐渐增加趋势(图 4A);随着坡度增加,植被覆盖度均值呈逐步增加后微弱衰减趋势(图 4B);坡向对植被覆盖度影响较小,除平地植被覆盖度仅为9.12%,其余坡向内植被覆盖度基本一致(图 4C)。
1995-2015年,不同地形条件下植被覆盖度时间变化规律亦存在差异。在第1高程(<300 m)地区,地表植被覆盖度由39.14%提高至45.88%,增长幅度为6.74%;第2~6高程(300~800 m)地区,地表植被覆盖度分别提高11.47%、11.20%、9.83%、10.15%和9.12%;第7高程(>800 m)地区地表植被覆盖度增加最快,增幅达13.71%(图 5A)。从1995到2015年,平坡区域,植被覆盖度仅增加2.09%,而缓坡、斜坡、陡坡、急坡和险坡地表植被覆盖度分别提高11.17%、12.30%、12.97%、13.93%和14.29%(图 5B)。1995-2015年,阴坡地表植被覆盖度增加幅度最小(7.61%),半阳坡地表植被覆盖度增加幅度最大(9.52%),阳坡和半阴坡地表植被覆盖度增加幅度分别为8.91%和8.97%(图 5C)。