庞 勇,蒙诗栎,史锴源,余 涛,王晓慧,牛晓栋,赵 旦,刘良云,冯 敏,覃先林,王 强,佃袁勇

1 中国林业科学研究院资源信息研究所, 北京 100091 2 国家林业和草原局林业遥感与信息技术重点实验室, 北京 100091 3 北京林业大学信息学院, 北京 100083 4 中国科学院空天信息创新研究院, 北京 100094 5 中国科学院青藏高原研究所, 北京 100101 6 黑龙江工程学院, 哈尔滨 150008 7 华中农业大学园艺林学学院, 武汉 430070

天然林是自然界中群落最稳定、生物多样性最丰富的陆地生态系统。全面保护修复天然林是加快构建以天然林为主体的健康稳定森林生态系统、构筑国土生态安全底线的重要举措。2000年10月,国务院正式批准了《长江上游、黄河上中游地区天然林资源保护工程实施方案》和《东北、内蒙古等重点国有林区天然林资源保护工程实施方案》,实施范围涉及长江上游、黄河上中游、东北内蒙古以及新疆、海南等重点国有林区的17个省(区、市)的734个县和167个森工局。2011年5月,为严格保护南水北调重要水源地,实施范围在一期工程基础上增加了丹江口库区的11个县(市、区)。实施天然林资源保护工程(简称天保工程),是我国林业由木材生产为主向生态建设为主转变的重要标志。天保二期工程于2020年底结束,天然林保护即将转入长期工作阶段,对天然林资源开展定量监测与评估,有助于全面、科学地评估天保工程对森林资源的保护成效,为全面推进天然林保护修复提供决策依据和科学支撑。

20 世纪80 年代以来,我国已经形成由国家森林资源连续清查(简称“一类清查”)、森林资源规划设计调查(简称“二类调查”)、森林作业设计调查和各类专项检查核查所组成的森林资源调查监测体系,为森林资源管理和决策提供了有力的技术支撑[1- 2]。一类清查每5年一次,二类调查每10年一次,各省调查的年份不统一,均无法反映同一个时间点上(如2000年)的森林资源状况,且没有空间连续的森林分布图,这一问题直到2010年开启了全国森林资源管理“一张图”才得到解决,但“一张图”没有2010年之前的数据。随着对生态文明建设的高度重视和深入推进,森林资源监测面临前所未有的挑战,对现有监测体系的时效性、协同性、现势性以及信息采集能力提出了更高的要求[3]。与森林资源调查方法不同,遥感技术以其观测范围广、监测周期长、速度快和精度高等特点为天保工程监测评价提供了便利手段。Vina等[4]用MODIS的植被覆盖度产品和NPP产品评价了2000—2010年间的中国县级尺度森林覆盖变化,研究表明天保工程的实施、减产禁伐措施对森林资源恢复起到了明显作用。Pang等[5]进一步指出由于人工促进森林恢复实施区域的广泛分布和地块面积多变性,有必要使用中高分辨率遥感数据开展中国森林生态恢复工程的定量评估。Lu等[6]综合地面调查和遥感数据的研究结果表明, 2000—2010年间天保工程区的碳储量占六大生态工程区总碳储量的58.5%,且保持着最大的碳汇能力。Chen等[7]用2000—2019年的MODIS时间序列LAI产品,评价了全球植被变绿的趋势,显示中印领先全球的植被增加趋势,其中中国天然林保护区域呈现出明显的植被增加趋势。部分学者[8-12]在陕西省黄龙山、新疆维吾尔自治区阿尔泰山、四川省甘孜州道孚县和川中丘陵区等天保工程区开展了多期动态监测,以获取天然林资源的时空变化,掌握天保工程实施成效。

天保工程区的全区域和全过程监测面临着范围面积广、空间异质性强和时间跨度大等特点,生产时空一致性强的大区域高精度森林覆盖产品具有较高难度。因此,如何准确监测我国天保工程区的森林覆盖情况并分析其时空变化成为一大挑战。当前国际上存在多套涵盖全天保工程区的全球土地覆盖遥感产品,包括:美国马里兰大学1998年发布的UMD Land Cover Classification 1998土地覆盖分类数据,该产品利用AVHRR数据,在1°,8 km和1 km三种不同空间尺度上,将全球分为14种不同土地覆盖类型[13-15]；国际卫星陆地地面气候项目(ISLSCP)数据计划II发布的IBGP全球土地覆盖数据GLCC 1992[16]和MODIS 2001[17]生产全球土地覆盖分类；欧洲航天局自2005年执行GlobCover计划,利用ENVISAT卫星的MERIS数据(300 m)生产的GlobCover 2005和GlobCover 2009两期全球陆地覆盖产品[18-19]。近年来,得益于计算机超算平台和大数据云服务的技术发展,土地覆盖产品的空间分辨率逐渐从公里级别提高到30 m甚至10 m尺度。马里兰大学Hansen等[20]发布的2000—2012年全球30 m森林覆盖和森林变化数据；Gong等[21]利用Landsat TM/ETM+获得首套全球尺度的30 m空间分辨率的土地覆盖类型数据FROM-GLC 30,并在此基础上利用Sentinel- 2影像更新生产了10 m空间分辨率的土地覆盖数据 FROM-GLC10[22]；刘涵和宫鹏等[23- 24]基于智慧遥感制图iMap架构,获取了21世纪以来全球逐年/季节的30 m土地覆盖和土地利用制图；陈军等[25-27]开放共享了2000年、2010年、2020年覆盖全球陆域的GlobeLand 30产品,产品的全球验证点精度为80.3%。吴炳方等[28-30]基于Landsat TM/ETM和HJ- 1卫星数据完成中国2010年土地覆盖数据ChinaCover 2010,并利用面向对象的向量相似性变化检测方法更新得到中国2000年土地覆盖产品ChinaCover 2000；Zhang等[31- 32]采用DataCube架构数据组织,基于全球地物图像光谱库的时序分类算法生产了全球30 m地表覆盖精细分类产品GLC_FCS 30,全球验证的产品精度为82.5%。

当前,我国天保工程资源评价的已有研究仍局限于特定区域,对全工程区系统性的研究较少,难以体现天保工程全局的实施效果及其空间分布特征。现阶段全球和全国尺度的土地覆盖产品虽相对丰富,但不同产品间的森林覆盖分布仍存在较大差异[33- 34]。针对上述问题,本文在已有公开发布的土地覆盖产品基础上,利用新的遥感数据合成算法和森林类型综合提升方法,生产了一套针对我国天保工程区森林覆盖情况的数据产品,对1997年启动前到2020年结束的23年间森林覆盖状况进行监测与评估,以系统评价天保工程实施过程中森林覆盖动态变化。

1 数据与方法

1.1 研究区概况

如图1 所示,根据天保工程分区和地理分布特征,将天保工程实施范围分为长江上游地区、黄河上中游地区、东北和内蒙古重点国有林区、新疆重点国有林区、海南重点国有林区5个分析单元。工程区内覆盖天然林面积约7330万hm2,占全国天然林总面积约69%。

图1 我国天保工程区实施范围以及验证样地点Fig.1 The coverage of natural forest protection project of China and the verification sample sites

1.2 遥感数据及预处理

研究采用美国地质调查局(USGS)Landsat 5/7/8卫星地表反射率数据,选择云覆盖低于70%的地表反射率影像作为输入数据,基于改进的最佳像元权值合成算法[35],通过Google Earth Engine(GEE)云数据平台合成了全国天保工程范围1997年和2020年两期生长季的无云遥感影像。影像包含8个通道:蓝波段(452—512 nm),绿波段(533—590 nm),红波段(636—673 nm),近红外波段(851—879 nm),短波红外1(1566—1651 nm),短波红外2(2107—2294 nm),归一化植被指数(Normalized Difference Vegetation Index,NDVI)和归一化水体指数(Normalized Difference Water Index,NDWI)。

1.3 遥感土地覆盖产品收集及预处理

本研究收集并整理了公开发布的3套共9个时期土地覆盖产品(表1),包括GlobalLand 30[25-27](http://www.globallandcover.com/)的2000、2010和2020年产品,GLC_FCS 30[31- 32](http://data.casearth.cn/sdo/detail/5fbc7904819aec1ea2dd7061)的2000、2015和2020年产品,以及ChinaCover[28-30]的2000、2010和2015年产品。GlobalLand 30产品分为耕地、森林、水体、人造地表等10个一级类； GLC_FCS 30产品分为30个精细类别,如森林分为常绿/落叶/开放落叶/密闭落叶阔叶、常绿/落叶/开放落叶/密闭落叶针叶以及混交林等；ChinaCover的6个一级类型为联合国政府间气候变化框架土地覆被类型,包括林地、草地、耕地、湿地、人工表面和其他,以及38个基于碳收支的土地覆被分类系统二级类型。由于不同土地覆盖产品的分类系统存在差异,且同一土地覆被类型的定义方式也有所不同。为了更好地利用各分类产品,经综合分析各产品的类别定义、评估不同产品的数据一致性同时兼顾天保工程土地覆盖监测具体需求,将各产品统一分类体系,分为10个一级类别,即耕地、森林、草地、灌木地、湿地、水体、苔原、人造地表、裸地、冰雪,其中,森林类型的定义为乔木覆盖且郁闭度不低于30%的土地,具体归类详见表2。

表1 土地覆盖产品收集情况

表2 不同土地覆盖产品统一分类体系

1.4 土地覆盖产品生产方法

本研究从数据处理至产品生产全链条采取全自动运算方式(图2):采用大数据云平台预处理数据、自动选取类型样本、机器学习分类、规则约束分类后处理以及专家知识判读更新等,生产了我国天保工程区1997年和2020年两期森林覆盖分类产品。在数据源构建方面,对合成的无云遥感影像以2.1°× 2.1°空间范围(瓦片)进行划分,相邻瓦片之间互为重叠0.1°,全工程区共划分为149个数据瓦片,构建出天保工程遥感数据集。

图2 天保工程区土地覆盖产品生产流程Fig.2 The workflow of land cover products in natural forest protection project

分类样本训练集决定了机器学习分类模型的分类精度,为构建可信分类样本库,在收集的3套共9个分类产品的基础上求取2000—2020年间同类土地覆盖类型的交集区域,自动生成各覆盖类型的分类样本库。分类时按照逐瓦片范围从对应区域样本库内随机抽样一定数目样本(每类别样本数目不超过1万个)来训练和验证随机森林分类模型,从而得到土地覆盖分类结果。采用规则约束的方法对土地覆盖分类结果进行分类后处理,将分类结果与目标年份接近的其他土地覆盖产品来综合决策森林的部分,即在森林并集区域里通过众数投票方式来最终确定的森林结果,森林并集之外的区域则采取其他产品投票决定。在分类后处理的基础上,经过人工目视判别得到2020年土地覆盖结果。

此外,利用了连续变化检测算法(Continuous Change Detection and Classification,CCDC)[36- 37]生成1997至2020年的长时间序列土地覆盖变化检测结果。最后,以2020年的土地覆盖结果为基准,结合CCDC变化检测结果对1997年的土地覆盖结果进行一致性分析与更新,最终得到天保工程区1997年和2020年土地覆盖产品。

1.5 土地覆盖变化分析方法

对2020和1997年的两期土地覆盖产品,在10个一级类层面使用专题图数学差值的方法生成类型变化图,通过对变化图中数值的解析和分析,得到各类型互相转换的统计矩阵。进一步对森林类型的转入和转出进行分析,得到森林类型的净变化与交互变化。

1.6 验证数据

为了从定量层面对天保工程土地覆盖产品进行准确评估,采用中国森林生态系统定位研究网络(Chinese Forest Ecosystem Research Network,CFERN)样地数据、森林资源规划设计调查数据和地面调查数据分别从像元尺度和小班尺度来验证2020年天保土地覆盖产品的分类精度。其中,(1)收集了92个CFERN森林生态站从2017到2019年调查的共1806块固定样地,固定样地形状为方形,样地面积0.04—0.06 hm2。收集的1806块固定样地中,位于天保工程区范围内的样地有1018块(图1),这些样地在2017和2019年调查的地物类型与2020年基本无变化,大部分森林生态站观测地物为森林样地(985块),少数为荒漠草原样地、灌木林样地、林内观测室和水量观测平衡场；(2)2020年东北、内蒙古地区共计8895个野外实地踏查数据,踏查范围涵盖黑龙江、吉林、辽宁、内蒙古东部以及河北北部地区,调查内容有样地位置定位、地类信息描述、拍照记录等调查,记录的地物类型涵盖耕地、林地、草地、湿地、水体、荒漠、人工表面、裸地共8类。调查的8895个野外样地中有963个位于天保工程区范围内,其中耕地、森林、湿地和草地类型占大多数(图1)；(3)2019年湖北省森林资源规划设计调查小班数据,通过将二类调查小班矢量数据转换为栅格数据,结合小班地类属性字段对栅格数据进行地类合并,共选择163875条小班记录构建地表覆盖类型检验样本。

2 结果

2.1 精度验证

2.1.1中国森林生态系统定位研究网络验证结果

本研究采用混淆矩阵计算精度结果,统计显示,天保工程土地覆盖遥感数据在CFERN的森林生态站点处精度较高,其中森林类别的用户精度为98.30%,生产精度达100%,各类型总体分类精度为98.04%,存在少量草地、灌木、耕地、人工表面误分为森林的情况。

2.1.2野外实地踏查点验证结果

将2020年东北踏查数据与土地覆盖产品进行精度统计,表3分类混淆矩阵结果显示,森林类别的分类精度均高于90%,其中用户精度为95.60%,生产精度为96.40%,各类型总体分类精度为89.82%,少部分草地、灌木和湿地被误分为森林。

表3 2020年东北地区野外实地踏查验证点分类混淆矩阵

2.1.3森林资源规划设计调查小班验证结果

图3展示了栅格化的2019年湖北省森林二类调查数据与其他土地覆盖产品的结果。以2019年湖北省森林二类调查数据为验证,与天保工程土地覆盖产品的混淆矩阵结果显示(表4),森林类别的用户精度为87.60%,生产精度为92.80%,总体分类精度为84.00%。由于一些新造未成林地、新退耕还林林地与部分耕地的光谱差异不太明显,造成严重的两者相互错分现象。灌木、裸地、草地的误分率与漏分率也较高。部分水田在森林二类调查数据中的地类为水体,而在分类影像中被分为了耕地,也是造成分类精度不高的原因。

图3 2019年湖北省森林资源规划调查数据和其他三种分类产品 Fig.3 Forest resource inventory and planning data of Hubei Province in 2019 and the other three land cover products

表4 2019年湖北省森林资源规划调查小班区划验证结果分类混淆矩阵

2.2 天保工程区森林覆盖空间分布格局

基于2020年土地覆盖产品,利用空间统计和数据分析方法,从东北及内蒙古重点国有林区、黄河上中游地区、长江上游地区、新疆和海南重点国有林区共五个区域研究天保工程区的森林分布特征,揭示天保工程区森林覆盖变化的区域分异规律。从整体看(图4),天保工程区的森林覆盖呈现东北、中部、西南丰富密集,西北稀疏量少的特点,多数森林位于在崇山峻岭之中。(1)东北及内蒙古等重点国有林区森林分布集中,主要分布在大兴安岭、小兴安岭和长白山脉地区,森林类型以落叶针叶林、落叶针阔混交林为主。区域森林覆盖程度高,森林覆盖率为69.11%,大兴安岭林区呈现森林与草地、湿地交错分布的格局,小兴安岭和长白山林区的森林与农田交错分布。(2)黄河上中游地区森林集中于我国400 mm等降水线以东的半湿润区,主要分布在祁连山脉、六盘山、贺兰山、黄龙山、吕梁山、太行山脉和秦岭地区,以常绿阔叶林、落叶阔叶林和针阔混交林为主。该区域日照充足,降水较少,森林覆盖率为15.15%,林区多呈现森林与草地、灌木交错分布的格局。(3)长江上游地区森林集中分布于山地区域,主要分布在大巴山脉、横断山脉一带,以常绿针叶林和常绿阔叶林为主。该地区温度及降水等气候条件优越,森林覆盖较高,森林覆盖率为52.54%,林区的森林与草地、农田交错分布。(4)新疆重点国有林区的森林主要分布在新疆天山山脉和阿尔泰山地区,以山地针叶林、落叶阔叶林为主,区域光照丰富,干燥少雨,受气温和降水影响,森林覆盖较低,森林覆盖率为5.82%,森林多与草地交织共存。(5)海南重点国有林区的自然条件优越,水、热、光资源丰富,森林主要分布在海南霸王岭和尖峰岭一带,以常绿阔叶林、热带山地/低地雨林为主。该区域森林覆盖率很高,森林覆盖率为74.80%。

图4 天保工程区2020年土地覆盖产品Fig.4 Land cover map of natural forest protection project in 2020

2.3 天保工程区森林覆盖变化时空分析

1997至2020年间,我国天保工程区森林覆盖整体呈增加趋势。从时间上看,截止2020年天保二期工程结束,全工程区森林面积净增加了1.59%(表5),森林净增长速率(净增长面积占1997年总森林面积的比例)为5.27%。其中,东北、内蒙古重点国有林区的森林覆盖在20多年间增加不明显,保持相对稳定态势,森林面积净增加了0.26%,森林净增长速率为0.38%；长江上游林区、黄河上中游林区的森林覆盖均表现出明显增长,森林面积分别净增加了2.74%和1.19%,森林净增长速率分别为5.49%和8.55%；新疆重点国有林区表现出较强的增长态势,森林面积净增加了1.07%,森林增长速率为22.50%； 海南重点国有林区森林覆盖增加最快,森林面积净增加23.11%,森林增长速率达44.72%。

表5 1997至2020年天保工程区森林覆盖率变化

从各区域的森林增长面积来看(图5),长江上游地区森林增长量最高,森林净增长占全工程区森林增长面积的50.97%,其次为黄河上中游地区以及新疆重点国有林区,分别占森林总增长面积的29.00%和12.35%,海南重点国有林区森林增长相对较少,占森林总增长面积的5.53%,东北及内蒙古等重点国有林区森林变化较稳定,森林增加不显著,森林净增长占森林总增长面积的2.15%。

图5 1997年至2020年天保工程区各区森林面积增长百分比 Fig.5 Percentage of forest area increasing in each reagion of natural forest protection project from 1997 to 2020

从空间分布上看(图6),森林明显增加的区域主要位于大兴安岭南部、小兴安岭和长白山脉的东部、山西西部、陕西中部和秦岭地区、宁夏南部、甘肃东南部、四川西部、重庆北部、贵州西部、云南北部以及新疆天山南部和阿尔泰山南部林区,和海南省西南部。森林减少区域主要分布在大兴安岭和横断山脉一带。

图6 1997年至2020年天保工程区森林变化空间分布图Fig.6 The spatial distribution of forest changes of natural forest protection project between 1997 and 2020

3 讨论

3.1 无云影像合成

由于光学卫星传感器的大气效应、视野限制和极易受云雾影响等原因,很难直接获得完整覆盖大区域的无云图像,影像无云合成方法主要用于解决云的影响、气溶胶污染、观测视角效应和数据利用等问题。天保工程区纬度跨越大且涉及多个不同气候带区域,这给无云影像合成带来一定挑战,本研究以地理位置和生态分区将全国天保工程区划分为若干区域,对每个区域分别采用不同参数合成植被生长季的无云遥感影像,合成出的无云无缝高质量影像为后续土地覆盖制图提供了重要且完整的图像信息基础。本文优化了最佳像元选取准则,解决了合成影像中由于多时相数据造成的“斑块”问题,提高了合成影像的辐射一致性。对于云、雪覆盖较多的海南、新疆等少数区域,把影像选取范围扩展到目标年份前后两年的数据。

3.2 已有土地覆盖产品对新分类产品的影响分析

三种共9期土地覆盖产品采用的分类策略以及分类系统不同,这种差异和不确定性在做分类系统转换前已经存在,是各产品在相同区域一致性低的主要因素之一。如在我国南部以及西南部的亚热带和热带地区,GLC_FCS 30产品相较其他产品而言识别的灌木较多森林较少,对于东北地区,ChinaCover产品的湿地类别与其他产品表现出较大差异,在西北干旱地区,各产品间的森林、草地和灌木类型相互差别都较大。同时,同系列土地覆盖产品间的不同时期也出现部分内容有较大分歧。本研究的分类训练样本来自于各产品类别相同的区域,分类产品间的不一致性问题则会造成局部地区某些类别的分类样本过少从而直接影响分类结果的稳定。虽然通过增加临近区域样本的迁移复用能略微提高分类精度,但有限样本的稳定分类依然是下一步需解决的问题。本研究生产的土地覆盖产品是一套专门针对我国天保工程区森林覆盖监测的数据产品,在公开发布的土地覆盖产品基础上构建可信分类样本库,采用森林确信分布区域(森林交集)和森林最大可能分布区域(森林并集)的约束规则来综合决策森林的部分,能更好地避免森林的误分和漏分。如图2所示的, 天保工程区土地覆盖产品生产流程中对2020年的基准产品还专门设计了目视检查的环节,通过有经验的判读人员将分类结果和高分辨率影像、合成的分类用遥感影像进行判读和比对,对自动分类不正确的斑块进行人工修改,从而确保基准产品的准确性。但由于时间和投入的因素,本文的结果仅对东北国有林区的一部分进行了检查和修改,其他地方如果也进行严格修改,有望进一步提高分类产品的精度。

3.3 干旱和半干旱区森林识别

我国生态环境脆弱的干旱和半干旱区广泛分布着的灌木林发挥了显著的生态防护作用。干旱/半干旱区的灌木林遥感识别一直是一个难点,由于灌木林生长比较低矮,且光谱曲线易和草地发生混淆,分类精度普遍比较低,基于几套全球土地覆盖产品也难以找到共同的区域,因此对该区域的森林覆盖存在低估情况。联合国粮农组织Bastin等[38]应用Google Earth平台上大量的高时空分辨率影像对210000个0.5 hm2样地进行图像解译,估计全球干旱区森林覆盖。研究表明,2015年干旱区森林比先前中低空间分辨率遥感影像估计的面积更大,总面积接近热带雨林或者北方森林的面积。对于干旱区森林尤其是灌木林的分类,有赖于森林资源基础数据、调查数据以及高空间分辨率影像来提高信息提取精度,这也是今后研究中需要提高的地方。

4 结论

本文在已公开发布的土地覆盖产品基础上,研制了我国天保工程区1997年和2020年的森林覆盖产品,分析了天保工程实施的23年间森林覆盖的空间分布及变化。

(1)利用中国森林生态站样地数据、地面调查数据以及森林资源二类调查数据验证了2020年的土地覆盖产品,统计结果表明天保工程土地覆盖产品精度较高,森林类型的用户精度分别为98.30%、95.60%、87.60%,生产精度分别为100%、96.40%、92.80%,总体分类精度为98.04%、89.82%、84.01%。

(2)自1997年开展天保工程截止2020年结束,全工程区森林覆盖呈增加趋势,森林覆盖率增加了1.59%,其中长江上游地区森林增长量最高,占全工程区森林增长面积的50.97%,其次为黄河上中游地区和新疆重点国有林区,分别占全工程区森林增长面积的29.00%和12.35%,海南重点国有林区、东北及内蒙古等重点国有林区的森林增长相对缓慢,占全工程区森林增长面积的5.53%和2.15%。

(3)从空间分布上看,森林明显增加的区域主要位于大兴安岭南部、小兴安岭和长白山脉的东部、山西西部、陕西中部和秦岭地区、宁夏南部、甘肃东南部、四川西部、重庆北部、贵州西部、云南北部、新疆天山南部和阿尔泰山南部林区,以及海南省西南部。森林减少区域主要分布在大兴安岭和横断山脉一带。

(4)整体上看,天保工程实施有效保护和恢复了我国天保工程区内的森林资源,森林覆盖范围稳定增加,天保工程实施效果显著。

致谢：中国森林生态系统定位研究网络样地数据由国家林业和草原科学数据中心提供,2020年东北、内蒙古野外实地踏查数据由中国科学院东北地理与农业生态研究所提供,特此致谢。