朱超凡，朱求安，2，王 樂，劉 佳，瞿莉莎

（1.河海大學(xué)水文水資源學(xué)院，南京 210000；2.國(guó)家地球系統(tǒng)科學(xué)數(shù)據(jù)中心，北京 100101）

自天然林保護(hù)工程（簡(jiǎn)稱天保工程）實(shí)施以來(lái)，工程區(qū)森林覆蓋率呈增加趨勢(shì)，其中長(zhǎng)江上游地區(qū)森林覆蓋率最高，其森林增加面積占全區(qū)森林增加面積的50.97%［1］。人口擴(kuò)張和城市化發(fā)展導(dǎo)致長(zhǎng)江中上游地區(qū)森林過度開發(fā)、亂砍濫伐、毀林開荒。森林遭受破壞會(huì)引起大面積水土流失以及環(huán)境退化，造成江湖、水庫(kù)淤積，降低了長(zhǎng)江本來(lái)的調(diào)蓄洪能力［2］。中國(guó)在各地推出天保工程、退耕還林還草工程，長(zhǎng)江流域的植被得到一定程度的恢復(fù)［3，4］。定期開展森林資源監(jiān)測(cè)有助于更好地評(píng)估天保工程對(duì)森林資源恢復(fù)的效益，以及為國(guó)家層面生態(tài)恢復(fù)決策建立科學(xué)基礎(chǔ)。

隨著衛(wèi)星遙感事業(yè)的發(fā)展，遙感制圖技術(shù)得到了跨越式發(fā)展，大批全球或區(qū)域尺度的遙感監(jiān)測(cè)產(chǎn)品紛紛出現(xiàn)［5］。這些遙感產(chǎn)品來(lái)源于不同衛(wèi)星監(jiān)測(cè)數(shù)據(jù)、使用不同分類系統(tǒng)與制圖技術(shù)，因此各遙感數(shù)據(jù)有自身的優(yōu)缺點(diǎn)［6，7］。土地利用數(shù)據(jù)對(duì)于揭示人類對(duì)環(huán)境的影響，探討環(huán)境、資源、生態(tài)的變化有重要意義。但土地覆蓋數(shù)據(jù)產(chǎn)品仍存在諸多問題，如空間分辨率低、精度低、不同產(chǎn)品間的一致性較差［8］。

天保工程歷時(shí)數(shù)年，早已對(duì)生態(tài)環(huán)境和人類生活產(chǎn)生重大影響。雖然現(xiàn)階段的土地覆蓋數(shù)據(jù)產(chǎn)品相對(duì)豐富，但不同產(chǎn)品間的森林覆蓋分布仍存在較大差異［9，10］。本研究通過對(duì)各土地覆蓋數(shù)據(jù)產(chǎn)品進(jìn)行不確定性分析及精度評(píng)價(jià)，系統(tǒng)評(píng)價(jià)天保工程實(shí)施前后長(zhǎng)江上游地區(qū)森林動(dòng)態(tài)變化，揭示天保工程下長(zhǎng)江上游地區(qū)森林變化過程。

1 研究區(qū)概況與方法

1.1 研究區(qū)概況

長(zhǎng)江發(fā)源于青藏高原，其上游流域地形錯(cuò)綜復(fù)雜，不僅受東南季風(fēng)和西南季風(fēng)影響，而且受青藏高原的影響，是氣候脆弱地區(qū)［11］。長(zhǎng)江上游為湖北省宜昌市以上河段，總體涵蓋9 個(gè)省（市、區(qū)）。長(zhǎng)江上游流域的面積達(dá)100 萬(wàn)km2，介于24°46′—36°31′N、90°47′—113°31′E。流域內(nèi)水系眾多，包括金沙江、岷江、烏江、嘉陵江等河流。

1.2 數(shù)據(jù)來(lái)源

本研究收集了來(lái)自各個(gè)國(guó)家和機(jī)構(gòu)的10 種土地覆蓋數(shù)據(jù)產(chǎn)品，包括中國(guó)逐年土地利用/土地覆蓋數(shù)據(jù)產(chǎn)品（China’s annual land use/land cover datasets，CLUD-A）［12］、歐盟委員會(huì)聯(lián)合研究中心開發(fā)的數(shù)據(jù)產(chǎn)品（GLC2000）［13］、歐洲航天局生產(chǎn)的全球土地覆蓋數(shù)據(jù)產(chǎn)品（GlobCover）［14］、清華大學(xué)地球系統(tǒng)科學(xué)系宮鵬團(tuán)隊(duì)生產(chǎn)的全球首個(gè)5 km 分辨率地表數(shù)據(jù)產(chǎn)品（GLASS-GLC）、馬里蘭大學(xué)建立的全球土地覆蓋數(shù)據(jù)產(chǎn)品（UMd）［15］、美國(guó)地質(zhì)調(diào)查局為國(guó)際地圈-生物計(jì)劃建立的全球土地覆蓋數(shù)據(jù)產(chǎn)品（IGBPDISCover）［16］、歐洲航天局發(fā)布的300 m 土地利用數(shù)據(jù)產(chǎn)品（European space agency climate change initiative-land cover，ESA CCI-LC）、中國(guó)科學(xué)院空天信息創(chuàng)新研究院發(fā)布的全球30m 精細(xì)地表覆蓋動(dòng)態(tài)監(jiān)測(cè)數(shù)據(jù)產(chǎn)品（GLC_FCS30）［17］、國(guó)家基礎(chǔ)地理信息中心和多家科研所研制的數(shù)據(jù)產(chǎn)品（GlobeLand30，GL30）［18，19］、1980—2015 年中國(guó)遙感土地利用/覆蓋數(shù)據(jù)產(chǎn)品（CLUDs）。

本研究以第六次全國(guó)森林資源清查結(jié)果為參考數(shù)據(jù)，評(píng)價(jià)各數(shù)據(jù)產(chǎn)品在2000 年反映的數(shù)據(jù)產(chǎn)品精度。全國(guó)第六次森林資源清查于1999—2003 年開展，優(yōu)化完善了各省份連續(xù)清查體系，并且在清查期內(nèi)，采用3S（遙感、全球定位、地理信息系統(tǒng)）等新技術(shù)，大幅提高了清查的準(zhǔn)確度并且增加了清查的調(diào)查內(nèi)容［20］。

1.3 研究方法

1.3.1 不同土地覆蓋數(shù)據(jù)產(chǎn)品森林分布的一致性分析 將收集到的遙感數(shù)據(jù)進(jìn)行分辨率的統(tǒng)一（1 km），根據(jù)數(shù)據(jù)所在的年份（1995、2000、2005、2010 年）計(jì)算不同土地覆蓋數(shù)據(jù)產(chǎn)品中森林空間分布的一致性（C），計(jì)算公式如下：

式中，Vp為土地覆蓋數(shù)據(jù)產(chǎn)品疊加后在該柵格上分類為森林的土地覆蓋數(shù)據(jù)產(chǎn)品數(shù)；TND為該階段數(shù)據(jù)產(chǎn)品數(shù)量；本研究將各階段的森林空間分布一致性（C）分為4 個(gè)層次，包括低一致性（25%~50%）、中一致性（50%~75%）、高一致性（75%~100%）、完全一致性（C=100%）。

1.3.2 不同土地覆蓋數(shù)據(jù)產(chǎn)品精度驗(yàn)證 本研究以全國(guó)第六次森林資源清查數(shù)據(jù)為依據(jù)，利用中國(guó)省級(jí)行政邊界矢量圖，從2000 年的土地覆蓋數(shù)據(jù)產(chǎn)品中分別提取長(zhǎng)江上游9 個(gè)?。ㄊ?、區(qū)）的森林面積，然后在各區(qū)域內(nèi)比較土地覆蓋數(shù)據(jù)產(chǎn)品森林面積與清查的森林面積，定量評(píng)價(jià)土地覆蓋數(shù)據(jù)產(chǎn)品在各區(qū)域的森林面積精度差異。

1.3.3 森林變化趨勢(shì) 本研究采用一元線性回歸分析來(lái)計(jì)算每個(gè)柵格的變化趨勢(shì)。計(jì)算公式如下：

式中，K為斜率；n為年數(shù)；xi表示第i年的變量值。其中，K＞0，表明該變量呈上升趨勢(shì)，K＜0，表明該變量呈下降趨勢(shì)。

2 結(jié)果與分析

2.1 不同土地覆蓋數(shù)據(jù)產(chǎn)品森林面積對(duì)比及一致性分析

由圖1 可知，不同土地覆蓋數(shù)據(jù)產(chǎn)品在森林面積及其變化趨勢(shì)上存在明顯差異。GLASS-GLC 數(shù)據(jù)產(chǎn)品的森林面積最大，為5.5×105km2，整體呈上升趨勢(shì)。GLC_FCS30 數(shù)據(jù)產(chǎn)品的森林面積為3.6×105~3.9×105km2，整體呈先上升后平穩(wěn)的趨勢(shì)。Globe-Land30 數(shù)據(jù)產(chǎn)品的森林面積變化范圍為3.5×105~4.0×105km2，2000—2010 年呈上升趨勢(shì)，2010—2020年呈下降趨勢(shì)。CLUDs 數(shù)據(jù)產(chǎn)品的森林面積基本保持在3.3×105km2，整體變化幅度較小。ESA CCILC 數(shù)據(jù)產(chǎn)品森林面積變化范圍為3.0×105~3.2×105km2，整體呈緩慢上升趨勢(shì)。CLUD-A 數(shù)據(jù)產(chǎn)品森林面積基本保持在2.0×105km2，整體變化幅度偏小。UMd 數(shù)據(jù)產(chǎn)品的森林面積為3.2×105km2，同一期間IGBPDISCover 數(shù)據(jù)產(chǎn)品的森林面積僅為0.8×105km2。GLC2000 數(shù)據(jù)產(chǎn)品的森林面積為2.5×105km2。GlobCover 數(shù)據(jù)產(chǎn)品森林面積為2.2×105~2.6×105km2。

圖1 長(zhǎng)江上游地區(qū)不同土地覆蓋數(shù)據(jù)產(chǎn)品的森林面積

根據(jù)1995、2000、2005、2010 年不同土地覆蓋數(shù)據(jù)產(chǎn)品的森林一致性比較結(jié)果（圖2）可知，高一致性和完全一致性的區(qū)域集中在長(zhǎng)江上游南方區(qū)域和四川盆地外圍，而對(duì)于四川盆地，從不同土地覆蓋數(shù)據(jù)產(chǎn)品中提取的森林面積存在較大差異。

由表1 可知，土地覆蓋數(shù)據(jù)產(chǎn)品森林面積完全一致的區(qū)域占整個(gè)長(zhǎng)江上游地區(qū)森林總面積的0.43%~4.17%，1995 年完全一致的區(qū)域占整個(gè)長(zhǎng)江上游地區(qū)森林總面積的4.17%，2010 年完全一致的區(qū)域占整個(gè)長(zhǎng)江上游地區(qū)森林總面積的0.43%。土地覆蓋數(shù)據(jù)產(chǎn)品森林高一致的區(qū)域占整個(gè)長(zhǎng)江上游地區(qū)森林總面積的24.93%~28.78%，2005 年高一致的區(qū)域占整個(gè)長(zhǎng)江上游地區(qū)森林總面積的28.78%，2010 年高一致的區(qū)域占整個(gè)長(zhǎng)江上游地區(qū)森林總面積的24.93%。2000 年和2010 年因?yàn)榇嬖诖罅康耐恋馗采w數(shù)據(jù)產(chǎn)品，導(dǎo)致土地覆蓋數(shù)據(jù)產(chǎn)品的一致性相對(duì)較低。

表1 長(zhǎng)江上游地區(qū)土地覆蓋數(shù)據(jù)產(chǎn)品森林一致性（單位：%）

2.2 不同土地覆蓋數(shù)據(jù)產(chǎn)品森林面積精度評(píng)價(jià)

表2 是2000 年的土地覆蓋數(shù)據(jù)產(chǎn)品與全國(guó)第六次森林資源清查在長(zhǎng)江上游地區(qū)森林面積的比對(duì)結(jié)果。不同的土地覆蓋數(shù)據(jù)產(chǎn)品與全國(guó)第六次森林資源清查在森林面積上存在一定偏差。除CLUD-A數(shù)據(jù)產(chǎn)品外，其他土地覆蓋數(shù)據(jù)產(chǎn)品的數(shù)據(jù)與全國(guó)第六次森林資源清查數(shù)據(jù)有較大偏差。偏差最大的土地覆蓋數(shù)據(jù)產(chǎn)品為GLASS-GLC，偏差為23.00%，其次為GLC_FCS30、GL30 數(shù)據(jù)產(chǎn)品，偏差分別為15.00%、11.00%。綜上，森林面積精度評(píng)價(jià)的土地覆蓋數(shù)據(jù)產(chǎn)品中，CLUD-A 數(shù)據(jù)產(chǎn)品精度最高。

表2 不同土地覆蓋數(shù)據(jù)產(chǎn)品森林面積精度評(píng)價(jià)（單位：%）

2.3 森林面積變化趨勢(shì)分析

由表2 可知，研究區(qū)內(nèi)CLUD-A 數(shù)據(jù)產(chǎn)品時(shí)間序列最完整，在收集的多種土地覆蓋數(shù)據(jù)產(chǎn)品中森林面積精度最高，因此本研究選用CLUD-A 數(shù)據(jù)產(chǎn)品，并從空間尺度上評(píng)估長(zhǎng)江上游地區(qū)34 年間森林變化趨勢(shì)。圖3 是CLUD-A 數(shù)據(jù)產(chǎn)品在4 個(gè)時(shí)間段的森林變化趨勢(shì)，在1982—1990 年長(zhǎng)江上游地區(qū)森林面積呈增加趨勢(shì)，并且集中在四川省境內(nèi)；1991—2000 年長(zhǎng)江上游地區(qū)森林面積呈減少趨勢(shì)，主要變化區(qū)域仍在四川省境內(nèi)；2001—2010 年，長(zhǎng)江上游地區(qū)森林面積整體呈增加趨勢(shì)，并且森林面積變化區(qū)域的重心開始逐步向南部?jī)A斜；2011—2015 年，森林面積呈先增加后減少趨勢(shì)。

圖3 長(zhǎng)江上游地區(qū)各階段森林變化趨勢(shì)

3 小結(jié)與討論

本研究結(jié)果表明，目前雖然有眾多土地覆蓋數(shù)據(jù)產(chǎn)品，但是不同土地覆蓋數(shù)據(jù)產(chǎn)品之間存在較大的精度差異。相關(guān)研究也證明了這一點(diǎn)，牛振國(guó)等［21］發(fā)現(xiàn)Globe Cover2009 數(shù)據(jù)產(chǎn)品在中國(guó)范圍數(shù)據(jù)總體精度只有32.46%。Giri 等［22］指出GLC2000和MODIS 數(shù)據(jù)產(chǎn)品之間不同類型數(shù)據(jù)面積相差較大，空間一致性也存在差異。目前，雖然遙感技術(shù)發(fā)展迅速，各數(shù)據(jù)產(chǎn)品總體精度較高，但是不同遙感數(shù)據(jù)之間仍存在差異。造成這種差異的原因有很多，吳文斌等［23］認(rèn)為由于不同時(shí)間生產(chǎn)的遙感數(shù)據(jù)在選擇數(shù)據(jù)源和制圖技術(shù)上存在分歧，從而導(dǎo)致不同遙感數(shù)據(jù)在地表覆蓋高度混合的區(qū)域有更多誤差和不確定性。Gong 等［24］認(rèn)為由于分類系統(tǒng)、分類技術(shù)、影像獲取時(shí)間和空間分辨率等原因，致使不同遙感數(shù)據(jù)之間存在精度差異。

自從天保工程實(shí)施以來(lái)，長(zhǎng)江上游地區(qū)的森林逐漸得到恢復(fù)［25］。天保工程在長(zhǎng)江上游地區(qū)主要實(shí)施的區(qū)域?yàn)槲髂系貐^(qū)，主要包括云南省、四川省、貴州省和重慶市。在天保工程實(shí)施下西南地區(qū)森林得到恢復(fù)。國(guó)政等［26］指出西南地區(qū)的森林覆蓋率由1998 年的30.49%上升至2008 年的40.2%。本研究結(jié)果也表明，在天保工程實(shí)施下長(zhǎng)江上游地區(qū)的森林覆蓋率得到了恢復(fù)，與國(guó)政等［26］的研究結(jié)果一致。

本研究指出長(zhǎng)江上游地區(qū)的森林覆蓋率可以分為4 個(gè)階段，第一階段是1982—1990 年，森林面積呈增加趨勢(shì)；第二個(gè)階段是1991—2000 年，森林面積呈減少趨勢(shì)；第三個(gè)階段是2001—2010 年，森林面積呈增加趨勢(shì)；第四個(gè)階段是2011—2015 年，森林面積呈減少趨勢(shì)。張亮等［27］指出長(zhǎng)江流域植被覆蓋度在1980—2000 年呈微弱增加，在該階段，長(zhǎng)江上游地區(qū)的諸多子流域內(nèi)植被覆蓋度都呈先增后減的趨勢(shì)。劉揚(yáng)揚(yáng)等［28］利用NDVI 數(shù)據(jù)對(duì)2001—2010年長(zhǎng)江上游地區(qū)植被時(shí)空變化進(jìn)行探討，結(jié)果表明，在該階段長(zhǎng)江上游地區(qū)植被呈增加趨勢(shì)。劉祖英等［29］認(rèn)為2010—2015 年是第一輪退耕還林工程末期，退耕的生態(tài)補(bǔ)貼不能滿足退耕農(nóng)戶的生活需要，部分地區(qū)出現(xiàn)復(fù)耕現(xiàn)象，從而導(dǎo)致部分地區(qū)林地退化。