重慶市植被凈初級(jí)生產(chǎn)力的遙感估算及分析

簡(jiǎn)太敏

摘要[目的]遙感估算和分析2007年重慶市植被凈初級(jí)生產(chǎn)力（NPP）。[方法]基于CASA模型，結(jié)合SPOT/VEGETATION NDVI S10 數(shù)據(jù)、植被覆蓋分類圖、氣象數(shù)據(jù)等，對(duì)2007年重慶市NPP進(jìn)行遙感估算，在GIS及RS的支持下，進(jìn)行數(shù)據(jù)的前期處理，對(duì)2007年重慶市NPP分布狀況及年內(nèi)變化情況進(jìn)行分析。[結(jié)果] NPP平均值為112.343～799.391 g C/m2，研究區(qū)植被凈初級(jí)生產(chǎn)力區(qū)域變化明顯，季相差異顯著，其中高值主要出現(xiàn)在渝東南、東北區(qū)域，低值主要存在于渝西北區(qū)域；NPP值夏季（6～8月）﹥春季（3～5月）﹥秋季（9～11月）﹥冬季（12～2月），年內(nèi)闊葉林植被NPP最高，水生植被NPP生產(chǎn)量最低。[結(jié)論] 該NPP估算方法能較好地模擬重慶市的NPP及其時(shí)空格局，估算結(jié)果比較符合實(shí)際，能應(yīng)用于NPP估算項(xiàng)目中。

關(guān)鍵詞植被凈初級(jí)生產(chǎn)力；CASA模型；遙感估算；重慶市

中圖分類號(hào)S181文獻(xiàn)標(biāo)識(shí)碼A文章編號(hào)0517-6611（2016）04-080-03

Remote Sensing Estimation and Analysis of Vegetation Net Primary Productivity in Chongqing City

JIAN Taimin（Ningxia Normal University， Guyuan， Ningxia 756099）

Abstract[Objective] To estimate and analyze the vegetation net primary productivity （NPP） in Chongqing City by remote sensing. [Method] Based on CASA model， remote sensing estimation of NPP in Chongqing in 2007 was carried out by using meteorological data， SPOT/VEGETATION NDVI S10 data and vegetation classification figure. Under the support of GIS and RS， data pretreatment was adopted， and NPP distribution and change situation were analyzed in Chongqing in 2007. [Result] The NPP average value was 112.343-799.391 gC/m2； vegetation NPP changed significantly in research region； seasonal differences were significant. The high values appeared in southeast and northeast areas of Chongqing； while low values appeared in northwest Chongqing. The NPP values were in the order of summer （JuneAugust）﹥ spring （MarchMay） ﹥ autumn （SeptemberNovember） ﹥ winter （December -February）. NPP of broadleaved forest vegetation was the maximum； but NPP of aquatic vegetation was the minimum. [Conclusion] This NPP estimation method can simulate the NPP in Chongqing City and its spatial and temporal pattern. The estimation results are relatively realistic， and can be used for the NPP estimation.

Key wordsNet Primary Productivity （NPP）； CASA model； Remote sensing estimation； Chongqing

植被凈初級(jí)生產(chǎn)力（Net Primary Productivity，NPP）是指綠色植物在單位面積、單位時(shí)間內(nèi)所累計(jì)的有機(jī)物數(shù)量，是由光合作用所產(chǎn)生的有機(jī)質(zhì)總量扣除自養(yǎng)呼吸后的剩余部分[1-2]，通常以干重表示。NPP作為地表碳循環(huán)的重要組成部分，不僅直接反映植被群落在自然環(huán)境條件下的生產(chǎn)能力，表征陸地生態(tài)系統(tǒng)的質(zhì)量狀況，還是判定生態(tài)系統(tǒng)碳匯和調(diào)節(jié)生態(tài)過(guò)程的主要因子[3]，在全球變化及碳平衡中起著重要作用。自20世紀(jì)60年代國(guó)際地圈-生物圈（International GeosphereBiosphereProgramme，IGBP）計(jì)劃以來(lái)，全球及區(qū)域NPP得到了大規(guī)模研究，IGBP和全球變化與陸地生態(tài)系統(tǒng)（GCTE）和京都議定書中均把植被的NPP研究確定為核心內(nèi)容之一[1]。近年來(lái)，各國(guó)學(xué)者進(jìn)行了大量的植物NPP測(cè)定，且以測(cè)定資料為基礎(chǔ)，結(jié)合氣候環(huán)境因子建立模型，對(duì)植被NPP的區(qū)域分布進(jìn)行評(píng)估。而利用遙感技術(shù)所獲數(shù)據(jù)來(lái)實(shí)現(xiàn)對(duì)地表植被凈初級(jí)生產(chǎn)力的估算，是最近10年來(lái)NPP的模型建立和估算方法研究的熱點(diǎn)，由于遙感不受過(guò)程模型所需特定觀測(cè)條件的限制或因素的影響，因而得到了迅速發(fā)展。我國(guó)有關(guān)于NPP的研究起步較晚，直至20世紀(jì)70年代后期才進(jìn)行了一些不同植被的生產(chǎn)量研究，取得了一些成果，但是數(shù)據(jù)精度不一，且用來(lái)推算大面積NPP有一定困難，大多借用國(guó)外一些NPP模型或?qū)δＰ瓦M(jìn)行改進(jìn)研究[3]，在20世紀(jì)90年代末進(jìn)入新領(lǐng)域，發(fā)展速度較快，研究主要集中于植被NPP測(cè)算模型及其估算、分析。目前，仍欠缺對(duì)重慶市NPP遙感估研究。鑒于此，筆者采用基于光能利用率估算的CASA模型 [3-4]，對(duì)2007年重慶市NPP進(jìn)行遙感估算，旨在為該模型參數(shù)在我國(guó)區(qū)域尺度的改進(jìn)及有效應(yīng)用提供科學(xué)依據(jù)。

1研究區(qū)概況

重慶市位于我國(guó)西南地區(qū)的經(jīng)濟(jì)文化中心，地處我國(guó)東西部地區(qū)的過(guò)渡地帶，地理坐標(biāo)為105°11′～110°12′ E，28°10′～32°15′ N，與陜、湘、鄂、川、貴等省毗鄰，位于長(zhǎng)江上游地區(qū)，是西南區(qū)域工商業(yè)重鎮(zhèn)，水陸交通發(fā)達(dá)，是商貿(mào)口岸與水陸交通樞紐。重慶市總面積約8.24萬(wàn)km2，目前所轄40個(gè)區(qū)縣，橫跨川中方山丘陵、川東平行嶺谷和盆南邊緣山地三大自然地理區(qū)，地形以山地、丘陵為主，是亞熱帶季風(fēng)性濕潤(rùn)氣候的典型區(qū)域，全境的兩江（長(zhǎng)江、嘉陵江）干流為“山城”氣候增加了色彩。特殊的地形、氣候形成了雨熱同季、冬暖春早、夏熱多雨的特征，這種氣候適合農(nóng)作物的生長(zhǎng)發(fā)育。全年氣候宜人，年均溫在18 ℃左右，冬季平均氣溫在6～8 ℃，夏季平均氣溫在27～29 ℃，常年降雨量達(dá)1 000～1 400 mm。復(fù)雜的地形與濕熱氣候環(huán)境共同形成了重慶市特殊的生態(tài)環(huán)境。

2數(shù)據(jù)來(lái)源與研究方法

2.1數(shù)據(jù)來(lái)源該研究采用的數(shù)據(jù)分為氣象數(shù)據(jù)、遙感數(shù)據(jù)和植被覆蓋分類數(shù)據(jù)3部分。①氣象數(shù)據(jù)以氣溫、降水月值數(shù)據(jù)、地面輻射數(shù)據(jù)等為主；利用ArcGIS 10.0數(shù)據(jù)插值工具，根據(jù)各氣象站點(diǎn)的定位坐標(biāo)信息，對(duì)站點(diǎn)的點(diǎn)值氣象數(shù)據(jù)通過(guò)克里金插值法獲取像元大小與NDVI數(shù)據(jù)一致、投影相同的氣象要素面狀柵格圖。②遙感數(shù)據(jù)主要來(lái)源于國(guó)家自然科學(xué)基金委員會(huì)“中國(guó)西部環(huán)境與生態(tài)科學(xué)數(shù)據(jù)中心”，SPOT/VEGETATION S10按旬合成的空間分辨率為1.0 km的東亞地區(qū)植被歸一化指數(shù)產(chǎn)品，數(shù)據(jù)時(shí)間段為2007年1～12月，共36個(gè)時(shí)相；數(shù)據(jù)經(jīng)幾何校正后，應(yīng)用ENVI 4.7軟件的波段計(jì)算，求取NDVI月值數(shù)據(jù)，并通過(guò)掩膜處理獲取重慶市區(qū)域尺度數(shù)據(jù)。③植被覆蓋分類數(shù)據(jù)是通過(guò)掃描儀器掃描《重慶市地圖集》中的植被類型圖經(jīng)過(guò)數(shù)據(jù)矢量化得到的；對(duì)數(shù)據(jù)分類做合并處理，最終得到所使用的8類植被。

2.2研究方法

2.2.1NPP的估算。采用1993年P(guān)otter提出的光能利用率模型（CarnegieAmesStanford Approach，CASA）[5]來(lái)估算NPP。CASA是由遙感、氣象、植被以及土壤類型數(shù)據(jù)共同驅(qū)動(dòng)的光能利用率模型，屬于過(guò)程模型范疇[6-7]。隨著計(jì)算機(jī)技術(shù)、遙感與地理信息系統(tǒng)等手段的不斷完善，CASA模型的應(yīng)用已成為大面積區(qū)域NPP和碳循環(huán)研究的主要手段之一[7]。CASA模型中的NPP主要由植物吸收的光和有效輻射（APAR）和光能利用率（ε）2個(gè)變量來(lái)確定：NPP（x，t）=APAR（x，t）×ε（x，t），式中，t為時(shí)間；x為空間位置；APAR（x，t）為像元x在t月吸收的光合有效輻射；ε（x，t）為像元 x在 t月的實(shí)際光能利用率。

2.2.2相關(guān)性分析。為對(duì)比2007年12個(gè)月NPP的月平均變化幅度及其與氣候因子之間的相互關(guān)系，筆者先對(duì)月均NPP、月總降水量和月均氣溫進(jìn)行標(biāo)準(zhǔn)化處理，即進(jìn)行了歸一化處理，將12個(gè)月的各因子最大值取1.0，最小值取0，然后求算以上每個(gè)因子的波動(dòng)規(guī)律與變化趨勢(shì)。應(yīng)用SPSS軟件做月NPP值與氣象因子的相關(guān)性分析，計(jì)算得到其相關(guān)性，并通過(guò)0.01顯著性檢驗(yàn)其線性相關(guān)。最終發(fā)現(xiàn)，NPP的變化與氣溫、降水的相關(guān)性均較高，達(dá)0.687 0以上，但氣溫的相關(guān)性略高于與降水的相關(guān)性。

3結(jié)果與分析

3.1NPP季相分布NPP平均值在季節(jié)分配上主要表現(xiàn)為夏季（6～8月）﹥春季（3～5月）﹥秋季（9～11月）﹥冬季（12～2月）。由于重慶市雨熱同期，主要在夏季，而夏季光照充足，以致估算NPP的氣候因素（如降雨、輻射、氣溫、光照等）在夏季最大、冬季最小，也因此影響到了NPP的季節(jié)分配，但是空間分布趨勢(shì)變化不大。春季高值區(qū)域分布范圍較大，且分布片區(qū)較密集；夏季東南部高值區(qū)減少，東北部中高值區(qū)有所增加，但不明顯；秋季高值區(qū)進(jìn)一步減少，東南及西南部大片為中高值區(qū)域，冬季高值區(qū)域比較破碎，中、低值區(qū)域面積較大，且連續(xù)性較好。

3.2NPP年內(nèi)分布對(duì)2007年重慶市NPP的估算結(jié)果為112.343～799.391 gC/m2。如圖1所示，NPP在空間分布上，與地形相關(guān)性較高，高值區(qū)域大部分分布于植被密度較高的山地地區(qū)，地質(zhì)區(qū)域分布于地勢(shì)較平坦的區(qū)域，具體情況如下：渝東北、東南地區(qū)NPP較高，渝西北較低，這與重慶市植被覆蓋類型有關(guān)，由植被覆蓋類型可知，渝東北部分布有大片針闊葉林，渝東南受湖南總輻射數(shù)據(jù)的影響（總輻射數(shù)據(jù)圖12007年重慶市NPP空間分布情況

Fig.1Spatial distribution of vegetation NPP in Chongqing in 2007是應(yīng)用重慶市一個(gè)站點(diǎn)與周邊站點(diǎn)數(shù)據(jù)做插值獲取），以致NPP值也相應(yīng)升高。在行政單位上巫溪、武隆縣NPP較高，這與該2個(gè)縣的山地地形森林覆蓋率較高密切相關(guān)，潼南縣最低；主城區(qū)由于森林植被覆蓋度較低，但是人工綠化草叢、園林規(guī)劃及減噪林的建造致使主城區(qū)NPP為中低值。

3.3NPP植被分布對(duì)重慶市8類植被的NPP進(jìn)行估算，由表1可知NPP在各植被類型分布與葉面積、植被含水量等相關(guān)性較為吻合，具體體現(xiàn)在：NPP以闊葉林最高，達(dá)465.98 gC/m2，水生植被NPP最低，僅202.17 gC/m2，這是由于闊葉林分布面積較廣，植被葉面積大，所吸收光照、水分指數(shù)高，且最大光能利用率指數(shù)高；而水生植被不僅因?yàn)榻邮艿墓狻岜凰治蘸拖鳒p最終吸收的較少，且其最大光能利用率最低，分布密度低，單株數(shù)量少，也對(duì)其生產(chǎn)NPP的能力產(chǎn)生一定影響。

4結(jié)論與討論

（1）筆者在CASA模型的基礎(chǔ)上，利用SPOT/VEGETATION NDVI S10 數(shù)據(jù)和氣象資料建立了重慶市NPP遙感估算模型，并對(duì)2007年重慶市的NPP進(jìn)行了模擬分析。結(jié)果表明： ①建立的基于SPOT/VEGETATION NDVI S10 數(shù)據(jù)的NPP估算方法能較好地模擬重慶市的NPP及其時(shí)空格局，估算結(jié)果比較符合實(shí)際，能應(yīng)用于NPP估算項(xiàng)目中；②重慶市NPP年內(nèi)季節(jié)性變化明顯，區(qū)域內(nèi)部變化顯著，這說(shuō)明NPP變化受植被覆蓋率及其氣候因子變化的影響較大。

（2）筆者基于當(dāng)前氣象部門收集處理的大量使用的SPOT/VEGETATION NDVI S10 數(shù)據(jù)，無(wú)需經(jīng)過(guò)復(fù)雜的大氣校正和雙向反射率校正來(lái)實(shí)現(xiàn)動(dòng)態(tài)估算NPP及監(jiān)測(cè)其變化，對(duì)氣象業(yè)務(wù)發(fā)展應(yīng)用具有積極意義[8]。

（3）由于目前有關(guān)于重慶市NPP的研究鮮見報(bào)道，筆者采用涵蓋重慶區(qū)域的其他研究結(jié)果與該研究估算結(jié)果進(jìn)行比對(duì)驗(yàn)證。朱文泉[9]對(duì)1989～1993年我國(guó)植被分類結(jié)果進(jìn)行估算得出NPP在100～400 g/m2。柯金虎等[10]對(duì)1982～1999年長(zhǎng)江流域的NPP進(jìn)行估算，得到重慶市的NPP為100～400 g/m2。而筆者研究得出2007年NPP為112.343～799.391 gC/m2，且空間分布趨勢(shì)相似。董丹等[11]利用CASA 模型對(duì)西南喀斯特NPP進(jìn)行模擬計(jì)算，其中，重慶市空間變化趨勢(shì)與筆者研究得出的一致，但其計(jì)算結(jié)果顯示重慶屬于中等值地區(qū)，1999～2003年NPP值為285～480 gC/m2，低值略低于文獻(xiàn)[9]的計(jì)算值，高值略高，這可能由于區(qū)域、年份、氣候差異，以致對(duì)NPP估算結(jié)果有很大影響的光能利用率計(jì)算中的水分、氣溫因子不同，同時(shí)筆者劃分的植被類型中水生植被與水體最大光能利用率值較低，以上因素使得計(jì)算結(jié)果與以往研究有部分差異。綜上所述，該研究估算結(jié)果與相近區(qū)域空間相關(guān)研究結(jié)果基本一致。

參考文獻(xiàn)

[1] WALKER B H，STEFFEN W L，CANADELL J，et al.The terrestrial biosphere and global change：Implications for natural and managed ecosystems：A synthesis of GCTE and related research[Z].IGBP Book Series Number 4，1997.

[2] 鄭凌云，張佳華.草地凈第一性生產(chǎn)力估算的研究進(jìn)展[J].農(nóng)業(yè)工程學(xué)報(bào)，2007（23）：280-283.

[3] 李世華，牛錚，李壁成.植被凈第一性生產(chǎn)力遙感過(guò)程模型研究[J].水土保持研究，2005，12（3）：126-128.

[4] 劉勇洪，權(quán)維俊，高燕虎.華北植被的凈初級(jí)生產(chǎn)力研究及其時(shí)空格局分析[J].自然資源學(xué)報(bào)，2010，25（4）：265-272.

[5] 李貴才.基于MODIS數(shù)據(jù)和光能利用率模型的中國(guó)陸地凈初級(jí)生產(chǎn)力估算研究[D].北京：中國(guó)科學(xué)院遙感應(yīng)用研究所，2004.

[6] POTTERC S，RANDERSON J T，F(xiàn)IELD C B，et al.Terrestrial ecosystem production：A processmodel based on global satellite and surface data [J].Global biogeochemical cycles，1993，7：811-841.

[7] 齊曄，HALL C A S.全球變化研究中的生物圈模型（Ⅰ）：初級(jí)生產(chǎn)力模擬[M]//李博.現(xiàn)代生態(tài)學(xué)講座.北京：科學(xué)出版社，1995：129-142.

[8] 王彧，黃耀，張穩(wěn)，等.中國(guó)農(nóng)業(yè)植被凈初級(jí)生產(chǎn)力模擬（Ⅱ）：模型的驗(yàn)證與凈初級(jí)生產(chǎn)力估算[J].自然資源學(xué)報(bào)，2006，21（6）：916-924.

[9] 朱文泉.中國(guó)陸地生態(tài)系統(tǒng)植被凈初級(jí)生產(chǎn)力遙感估算及其與氣候變化關(guān)系的研究[D].北京：北京師范大學(xué)，2005.

[10] 柯金虎，樸世龍，方精云.長(zhǎng)江流域植被凈第一性生產(chǎn)力及其時(shí)空格局研究[J].植物生態(tài)學(xué)報(bào)，2003，27（6）：764-770.

[11] 董丹，倪健.利用CASA 模型模擬西南喀斯特植被凈第一性生產(chǎn)力[J].生態(tài)學(xué)報(bào)，2011，31（7）：1855-1866.安徽農(nóng)業(yè)科學(xué)2016年

（上接第73頁(yè)）

[3] 楊蘇才，南忠仁，曾靜靜.土壤重金屬污染現(xiàn)狀與治理途徑研究進(jìn)展[J].安徽農(nóng)業(yè)科學(xué)，2006，34（3）：549-552.

[4] 張麗，郁建鋒，杜蘭芳.鉛脅迫對(duì)豌豆種子的毒害效應(yīng)[J].安徽農(nóng)業(yè)科學(xué)，2007，35（34）：10997-11000.

[5] 馮淑利，邵云，劉洋，等.Zn2+脅迫對(duì)小麥幼苗生理活性影響的研究[J].農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào)，2007，26（1）：140-145.

[6] 杜蘭芳，王立新，許俊，等.鑭對(duì)汞脅迫下豌豆生長(zhǎng)發(fā)育效應(yīng)的影響[J].科技通報(bào)，2007，23（5）：671-675.

[7] 張憲政.植物葉綠索含量測(cè)定丙酮乙醇混合法[J].遼寧學(xué)報(bào)，1986（3）：26-28.

[8] 薛剛，劉文紅，黃海濤，等.2，4二硝基苯肼比色法測(cè)定嬰幼兒配方食品和乳粉及其制品中總抗壞血酸的含量[J].食品研究與開發(fā)，2006，21（3）：48-50.

[9] 鄭光美.實(shí)用種子生理學(xué)[M].北京：中國(guó)農(nóng)業(yè)出版社，1997.

[10] 金璉，葉亞新，許明霞，等.鑭對(duì)鉛脅迫下小麥生理指標(biāo)的影響[J].安徽農(nóng)業(yè)科學(xué)，2005，33（12）：2260-2262.

[11] 杜蘭芳，郁建鋒，屠云霞.鑭對(duì)鉛脅迫下豌豆幼苗生長(zhǎng)發(fā)育的影響[J].北方園藝，2007（12）：14-18.

[12] 郭伯生.稀土農(nóng)用前景廣闊[J].中小企業(yè)科技，2006（1）：14-15.安徽農(nóng)業(yè)科學(xué)，Journal of Anhui Agri. Sci.2016，44（4）：83-86