王恩康， 池 源

(國家海洋局 第一海洋研究所， 山東 青島 266061)

廟島群島南部島群無居民海島植被凈初級生產(chǎn)力評估

王恩康， 池 源*

(國家海洋局 第一海洋研究所， 山東 青島 266061)

無居民海島具有顯著的戰(zhàn)略地位和重要的生態(tài)價(jià)值，植被凈初級生產(chǎn)力(NPP)評估對維護(hù)無居民海島生態(tài)系統(tǒng)穩(wěn)定性和活力具有重要意義.以廟島群島南部島群中螳螂島等5個(gè)無居民海島為研究對象，通過遙感手段，對海島NPP總量和密度及其空間分布特征進(jìn)行分析.結(jié)果顯示：1) 全年NPP總量為134.22 tC·a-1，平均密度為274.13 gC·m-2·a-1，并表現(xiàn)出明顯的季節(jié)差異，各島NPP平均密度由大到小依次為牛砣子島、羊砣子島、螳螂島、南砣子島和擋浪島；2) 不同地表覆蓋類型NPP平均密度由大到小依次為闊葉林、針葉林、草地、建設(shè)用地和裸地，各島NPP平均密度與裸地面積占比呈明顯的負(fù)相關(guān)，區(qū)域NPP密度與海拔、坡度、坡向和歸一化植被指數(shù)呈顯著正相關(guān)；3) 黑松和刺槐是北方無居民海島人工林的理想樹種，持續(xù)的人工林種植是提升無居民海島NPP的重要途徑.

無居民海島； 凈初級生產(chǎn)力； 空間分布； 海島人工林

根據(jù)是否屬于居民戶籍管理的住址登記地，我國海島可分為有居民海島和無居民海島，其中無居民海島數(shù)量眾多，占全國海島總數(shù)的90%以上[1].無居民海島是我國海洋國土的重要組成部分，是海洋經(jīng)濟(jì)發(fā)展的橋頭堡，是保護(hù)海洋環(huán)境的重要平臺，也是捍衛(wèi)國家主權(quán)、保障國防安全的重要支撐[2-4].同時(shí)，由于無居民海島特殊的地理位置、有限的面積大小和明顯的空間隔離，其生態(tài)系統(tǒng)具有顯著的脆弱性特征，主要表現(xiàn)在其更容易受到外界干擾的損害，且難以通過自身調(diào)節(jié)能力和自組織能力恢復(fù)[5-7].凈初級生產(chǎn)力(net primary productivity, NPP)作為碳循環(huán)的重要組分，是指植物去除呼吸消耗后剩余光合作用生產(chǎn)的有機(jī)質(zhì)[8].海島NPP不但能夠表征海島植物群落的生產(chǎn)能力[9]，也是反映海島生物棲息地適宜性的基礎(chǔ)因子，是海島生態(tài)系統(tǒng)穩(wěn)定性和活力的重要指標(biāo)[10].當(dāng)前，NPP的研究涉及草地、森林、濕地、農(nóng)田等不同生態(tài)系統(tǒng)[11-15]，在人類活動較為劇烈的有居民海島上也有相關(guān)研究[10]，但關(guān)于無居民海島NPP特征及其空間分布的研究鮮見報(bào)道.

以廟島群島南部島群中5個(gè)無居民海島為研究區(qū)，基于遙感影像和氣象數(shù)據(jù)，對無居民海島NPP總量、密度及其空間分布進(jìn)行分析，以期揭示典型無居民海島NPP特征及其主要影響因子，為維護(hù)無居民海島生態(tài)系統(tǒng)提供參考.

1研究區(qū)概況

廟島群島南部島群位于山東半島北側(cè)，黃海和渤海交匯處，是我國候鳥遷徙的重要通道和停歇地，是長島國家級自然保護(hù)區(qū)的重要組成，也是山東省長島縣的政治、人口和經(jīng)濟(jì)中心.擋浪島、螳螂島、羊砣子島、牛砣子島和南砣子島是南部島群中面積較大的5個(gè)無居民海島(大于0.05 km2)，其余無居民海島面積較小，植被覆蓋率較低.本文以上述5個(gè)無居民海島作為研究對象(圖1).研究區(qū)所在區(qū)域?yàn)闇貛Ъ撅L(fēng)氣候區(qū)，多年平均氣溫約為12℃，多年平均降水量約為540 mm；海島均為基巖海島，土壤以棕壤土為主，土壤多礫石，土質(zhì)較差[16]；海島森林均為人工林，包括以黑松為優(yōu)勢種的針葉林和以刺槐為主的闊葉林，海島人工林外緣分布著天然草地，海島周邊多為裸地.5個(gè)無居民海島中，擋浪島和螳螂島位于北長山島和小黑山島之間，其中擋浪島上建有小型碼頭以開展生態(tài)旅游，同時(shí)擁有燈塔以及太陽能和風(fēng)能發(fā)電設(shè)施；螳螂島上建有簡易的養(yǎng)殖看護(hù)房和小型風(fēng)力發(fā)電設(shè)施.羊砣子島和牛砣子島為廟島附屬海島，位于廟島西側(cè)，島上均有養(yǎng)殖看護(hù)房；南砣子島為大黑山島附屬海島，位于大黑山島南側(cè)，島上同樣建有養(yǎng)殖看護(hù)房和小型風(fēng)力發(fā)電裝置.

2材料與方法

2.1數(shù)據(jù)來源與處理

采用WorldView-1衛(wèi)星2013年全色波段遙感影像，空間分辨率為0.45 m.通過ArcGIS10.0軟件提取海島輪廓，得到海島面積、周長等基本信息；進(jìn)而依據(jù)人機(jī)交互解譯結(jié)果，將海島地表覆蓋分為以下5類：針葉林、闊葉林、草地、裸地和建設(shè)用地(圖2).2014年7月開展無居民海島繞島勘察，對研究區(qū)地表覆蓋類型進(jìn)行驗(yàn)證，根據(jù)驗(yàn)證結(jié)果校正現(xiàn)有數(shù)據(jù).

采用Landsat8衛(wèi)星4個(gè)時(shí)期(2013年4月21日、8月11日、11月15日和2014年1月2日，分別代表春夏秋冬四季)無居民海島所在區(qū)域的無云影像，空間分辨率為30 m.利用ENVI4.7軟件對影像進(jìn)行裁切、輻射定標(biāo)、波段運(yùn)算得到歸一化植被指數(shù)(Normalized difference vegetation index， NDVI).無居民海島的高程、坡度和坡向，采用2011年公布的Aster GDEM第二版DEM數(shù)據(jù)，水平分辨率30 m，通過ArcGIS10.0進(jìn)行提取.

太陽輻射、降雨量、氣溫、相對濕度等氣象數(shù)據(jù)來自長島縣氣象站和煙臺福山氣象站.

2.2NPP估算

根據(jù)CASA模型和相關(guān)研究，采用遙感影像和氣象數(shù)據(jù)估算無居民海島NPP及其空間特征，具體方法如下[10，17]：

NPP(x，t)=APAR(x，t)×ξ(x，t)，

(1)

APAR(x，t)=PAR(x，t)×FPAR(x，t)，

(2)

ξ(x，t)=ft(t)×fw(t)×ξmax，

(3)

式中,NPP(x，t)為x點(diǎn)t月凈初級生產(chǎn)力(gC·m-2·month-1)；APAR(x，t)是指x點(diǎn)t月吸收的光合有效輻射(MJ·m-2·month-1)；ξ(x，t) 表示x點(diǎn)t月的實(shí)際光能利用率(gC·MJ-1)；PAR(x，t)為x點(diǎn)t月的光合有效輻射(MJ·m-2·month-1)；FPAR(x，t)為x點(diǎn)t月光合有效輻射吸收比例(%)；ft(t)和fw(t)分別為研究區(qū)t月的氣溫脅迫因子和水分脅迫因子(%)；ξmax為最大光能利用率(gC·MJ-1).通過ArcGIS10.0對不同季節(jié)的NPP及其空間分布特征進(jìn)行計(jì)算，進(jìn)而得到全年無居民海島NPP總量、密度及其分布特征；采用單因素方差分析法(One-way ANOVA)探討不同季節(jié)和不同海島NPP密度的差異.

2.3主要影響因子分析

NPP及其空間分布往往受到氣候、植被、地形、土壤等環(huán)境因子的影響[10，14-15].本文為小空間尺度的研究，氣溫、水汽壓、太陽總輻射等氣象氣候因子作為NPP估算的重要參數(shù)，在整個(gè)研究區(qū)具有均一性，直接影響NPP的整體特征，但對其空間分布不構(gòu)成影響，這里不參與影響因子分析；由于無居民海島可達(dá)性較差，登島受到各種條件的限制，目前尚未能夠獲取無居民海島的土壤樣品，由鄰近有居民海島的土壤因子分析結(jié)果可知，研究區(qū)土壤性質(zhì)表現(xiàn)出明顯的一致性，即含水量和養(yǎng)分均偏低[18]，這使得土壤對NPP空間分布的影響可能較小.本次研究主要分析地表覆蓋類型、植被生長狀況、地形對NPP空間分布的影響.

2.3.1地表覆蓋類型 以上文的地表覆蓋分類為基礎(chǔ)，分析不同地表覆蓋類型中NPP總量和平均密度；采用One-way ANOVA探討不同地表覆蓋類型中NPP密度的差異性.2.3.2植被生長狀況 NDVI能夠準(zhǔn)確地反映植被生長狀況[19].統(tǒng)計(jì)無居民海島的NDVI及其空間特征，建立每個(gè)像元NPP和NDVI的一一對應(yīng)關(guān)系，通過相關(guān)分析法探討NDVI與NPP的相關(guān)性.

2.3.3地形 采用海拔、坡度和坡向因子，同樣通過相關(guān)分析法探討地形與NPP的相關(guān)性.其中，坡向原始值為按0～360順時(shí)針增大，正北與正南分別為0與180，按向陽性原則，進(jìn)行標(biāo)準(zhǔn)化，計(jì)算公式如下：

(4)

式中，ASx與Ax分別為x點(diǎn)標(biāo)準(zhǔn)化坡向值與原始坡向值.

3結(jié)果

3.1NPP估算結(jié)果

3.1.1全年NPP5個(gè)無居民海島全年NPP總量為134.22tC·a-1，NPP密度在0～656.31gC·m-2·a-1之間，平均密度為274.13gC·m-2·a-1,如圖3所示.

3.1.2不同季節(jié)NPP不同季節(jié)NPP估算結(jié)果見表1和圖4.夏季NPP總量最大，占全年75%以上，其次為春、秋兩季，NPP總量分別占比14.8%和6.9%，而冬季僅占2.0%；One-wayANOVA結(jié)果顯示，不同季節(jié)NPP密度表現(xiàn)出顯著的差異(P<0.001)，平均密度由大到小依次為夏季、春季、秋季和冬季.

3.1.3 不同海島NPP 不同海島NPP估算結(jié)果見表2.各島NPP總量由大到小依次為羊砣子島、螳螂島、南砣子島、牛砣子島和擋浪島；NPP密度也表現(xiàn)出顯著的差異(P<0.01)，將各島依據(jù)平均密度由大到小排列為：牛砣子島、羊砣子島、螳螂島、南砣子島和擋浪島.

3.2主要影響因子

3.2.1地表覆蓋類型 不同地表覆蓋類型中， NPP總量最大與最小分別為針葉林(占40%以上)與建設(shè)用地(僅占0.3%).NPP密度表現(xiàn)出顯著的差異(P<0.01)，闊葉林NPP平均密度最高，其次為針葉林、草地、建設(shè)用地和裸地，如表3所示.

3.2.2NDVI和地形 NPP與NDVI、海拔、坡度和坡向均呈顯著正相關(guān)，其中與NDVI的相關(guān)系數(shù)高達(dá)0.886；同時(shí)，NDVI與海拔、坡度和坡向也均表現(xiàn)出顯著的正相關(guān)關(guān)系，如表4所示.

注：**指在0.01水平(雙側(cè))上顯著相關(guān).

4討論與結(jié)論

4.1討論4.1.1NPP結(jié)果探討 將本文計(jì)算的NPP平均密度與其他區(qū)域進(jìn)行對比(見表5)，可以發(fā)現(xiàn)，本研究區(qū)NPP平均密度略低于全國平均密度324 gC·m-2·a-1，低于廣東省、江蘇省、盤錦市等沿海地區(qū)，高于甘肅、青海、西藏.由此可得，廟島群島南部島群的無居民海島NPP平均密度略低于全國的平均水平，高于同緯度的西部地區(qū)，但低于東部大陸沿海地區(qū).

4.1.2主要影響因子探討 本研究區(qū)以剝蝕山丘為主要地貌特征，存在淡水資源缺乏、土壤貧瘠等問題，加上大風(fēng)、寒潮等自然擾動，植物生長條件較為惡劣，使得原生林木發(fā)育較差，現(xiàn)有林木均為人工林，包括針葉林和闊葉林，分別以黑松和刺槐為優(yōu)勢種.本文結(jié)果顯示針葉林和闊葉林的平均密度分別達(dá)483.40 gC·m-2·a-1與336.74 gC·m-2·a-1，表明了黑松和刺槐具有較高的凈初級生產(chǎn)力，對提升海島生態(tài)系統(tǒng)穩(wěn)定性具有重要作用，可作為北方無居民海島人工林建設(shè)的理想樹種.

對于有居民海島而言，城鄉(xiāng)建設(shè)深刻影響著海島的地表覆蓋和NPP空間特征，已有研究表明，隨著建設(shè)用地面積占比的上升，有居民海島NPP平均密度明顯下降[10].對于無居民海島而言，人類開發(fā)活動較少，建設(shè)用地面積很小，對海島NPP影響不大；然而，不同海島裸地面積占比具有差異，擋浪島、南砣子島、螳螂島裸地面積較大，而牛砣子島、羊砣子島植被茂密，裸地面積較小.圖5顯示了各島NPP平均密度與裸地面積占比(Proportion of Bare Land， POB)呈明顯的負(fù)相關(guān)，線性擬合系數(shù)達(dá)0.950.

在大區(qū)域或大地形梯度的研究中，地形因子通過改變氣候因子進(jìn)而影響植被生長[25]，但本研究區(qū)的無居民海島面積較小，地形梯度也不甚明顯.地形主要通過影響地表覆蓋類型進(jìn)而影響NPP，海拔或坡度較低的區(qū)域分布于海島邊緣，多為裸地；隨著海拔和坡度的上升，地表覆蓋類型轉(zhuǎn)變?yōu)椴莸鼗蛄值?，NPP也明顯上升.坡向經(jīng)標(biāo)準(zhǔn)化處理后，越接近南向數(shù)值越高，也容易獲得更多的太陽輻射，從而促進(jìn)植物發(fā)育[26].NDVI與NPP具有較強(qiáng)的顯著正相關(guān)關(guān)系，鑒于其獲取通道的便捷性，NDVI可作為無居民海島凈初級生產(chǎn)力長期監(jiān)測的重要指標(biāo).

受面積制約，無居民海島NPP總量較小，但由于其被人類開發(fā)利用的程度較低，部分海島林木茂盛，擁有較高的NPP平均密度，而有的海島林木覆蓋率較低，裸地面積較大，海島裸地一部分是由出露巖石構(gòu)成，同時(shí)擁有一定面積的裸土地.無居民海島開發(fā)壓力相對較小，可以成為海島NPP的重要提供基地；同時(shí)，部分無居民海島尚有較大的NPP提升空間，可通過繼續(xù)開展人工林種植來實(shí)現(xiàn)海島NPP的提升，進(jìn)而增強(qiáng)無居民海島生態(tài)系統(tǒng)穩(wěn)定性.

4.2結(jié)論

1) 廟島群島南部島群中螳螂島等5個(gè)無居民海島全年NPP總量為134.22 t C·a-1，NPP密度在0～656.31 g C·m-2·a-1之間，平均密度為274.13 g C·m-2·a-1.NPP總量和平均密度表現(xiàn)出明顯的季節(jié)差異，春季、夏季、秋季和冬季NPP總量分別占14.8%、76.4%、6.9%和2.0%.

2) 按各島NPP總量由大到小排序，依次為羊砣子島、螳螂島、南砣子島、牛砣子島和擋浪島；按NPP平均密度由大到小排序，則為牛砣子島、羊砣子島、螳螂島、南砣子島和擋浪島.

3) 就地表覆蓋類型而言， NPP總量最大和最小的分別是針葉林和建設(shè)用地；將地表覆蓋類型按NPP平均密度由大到小排序，則依次為闊葉林、針葉林、草地、建設(shè)用地和裸地.各島NPP平均密度與裸地面積占比呈明顯的負(fù)相關(guān)；區(qū)域NPP平均密度與海拔、坡度、坡向和NDVI均呈顯著正相關(guān).

4) 黑松和刺槐是北方無居民海島人工林建設(shè)的理想樹種，持續(xù)的人工林種植是提升無居民海島NPP的重要途徑.

[1] 周學(xué)鋒. 基于建設(shè)海洋強(qiáng)國的無居民海島管理研究[J]. 經(jīng)濟(jì)地理， 2014， 34(1): 28-34．

[2] 胡增祥， 徐文君， 高月芬. 我國無居民海島保護(hù)與利用對策[J]. 海洋開發(fā)與管理， 2004(6): 26-29．

[3] 陳秋明. 基于生態(tài)—經(jīng)濟(jì)的無居民海島開發(fā)適宜性研究[D]. 廈門：廈門大學(xué)， 2009．

[4] 姬厚德， 羅美雪， 楊順良， 等. 無居民海島保護(hù)和利用規(guī)劃中開發(fā)空間的確定方法[J]. 海洋通報(bào)， 2016， 35(1):16-20．

[5] 池 源， 石洪華， 郭 振， 等. 海島生態(tài)脆弱性的內(nèi)涵、特征及成因探析[J]. 海洋學(xué)報(bào)， 2015， 37(12): 93-105．

[6] QIE L， LEE T M， SODHI N S， et al. Dung beetle assemblages on tropical land-bridge islands: small island effect and vulnerable species[J]. Journal of Biogeography， 2011， 38(4):792-804．

[7] KATOVAI E， BURLEY A L， MAYFIELD M M. Understory plant species and functional diversity in the degraded wet tropical forests of Kolombangara Island， Solomon Islands[J]. Biological Conservation， 2012， 145(1): 214-224．

[8] LIETH H， WHITTAKER R H. Primary productivity of the biosphere[M]. New York: Springer Verlag， 1975．

[9] FIELD C B， BEHRENFELD M J， RANDERSON J T， et al. Primary production of the biosphere: integrating terrestrial and oceanic components[J]. Science， 1998， 281(5374):237-240．

[10] 池 源， 石洪華， 王曉麗， 等. 廟島群島南五島生態(tài)系統(tǒng)凈初級生產(chǎn)力空間分布及其影響因子[J]. 生態(tài)學(xué)報(bào)， 2015， 35(24): 8094-8106．

[11] PARUELO J M， EPSTEI H E ， LAUENROTH W K， et al. A NPP estimates from NDVI for the central grassland region of the United States[J]. Ecology， 1997，78: 953-958．

[12] RUNNING S W， THORNTON P E， NEMANI R， et al. Global terrestrial gross and net primary productivity from the earth observing System[A]//SALA O E， JACKSON R B， MOONEY H A et al. Methods in Ecosystem Science. New York: Springer， 2000: 44-57．

[13] 朱文泉， 潘耀忠， 張錦水. 中國陸地植被凈初級生產(chǎn)力遙感估算[J]. 植物生態(tài)學(xué)報(bào)， 2007， 31(3):413-424．

[14] 龍慧靈， 李曉兵， 王 宏， 等. 內(nèi)蒙古草原區(qū)植被凈初級生產(chǎn)力及其與氣候的關(guān)系[J]. 生態(tài)學(xué)報(bào)， 2010， 30(5): 1367-1378．

[15] 王莉雯， 衛(wèi)亞星. 盤錦濕地凈初級生產(chǎn)力時(shí)空分布特征[J]. 生態(tài)學(xué)報(bào)， 2012， 32(19):6006-6015．

[16] 石洪華， 王曉麗， 王 嬡， 等. 北長山島森林喬木層碳儲量及其影響因子[J]. 生態(tài)學(xué)報(bào)， 2013， 33(19): 6363-6372．

[17] POTTER C S， RANDERSON J T， FIELD C B， et al. Terrestrial ecosystem production: a process model based on global satellite and surface data [J]. Global Biogeochemical Cycles， 1993， 7(4):811-841．

[18] CHI Y， SHI H H， WANG X L， et al. Impact factors identification of spatial heterogeneity of herbaceous plant diversity on five southern islands of Miaodao Archipelago in North China[J]. Chinese Journal of Oceanology and Limnology， 2016， 34(5): 937-951．

[19] 池 源， 郭 振， 石洪華， 等. 基于森林健康的北長山島人工林生物量[J]. 中國環(huán)境科學(xué)， 2016， 36(8): 2522-2535．

[20] 羅 艷， 王春林. 基于MODIS NDVI的廣東省陸地生態(tài)系統(tǒng)凈初級生產(chǎn)力估算[J]. 生態(tài)環(huán)境學(xué)報(bào)，2009， 18(4): 1467-1471．

[21] 王駟鷂， 劉振波. 江蘇省植被凈初級生產(chǎn)力時(shí)空分布格局研究[J]. 南京信息工程大學(xué)學(xué)報(bào)(自然科學(xué)版)， 2012， 4(4):321-325．

[22] 劉春雨， 董曉峰， 劉英英， 等. 甘肅省凈初級生產(chǎn)力時(shí)空變化特征[J]. 中國人口·資源與環(huán)境， 2014， 24(1):163-170．

[23] 衛(wèi)亞星， 王莉雯， 石迎春， 等. 青海省草地資源凈初級生產(chǎn)力遙感監(jiān)測[J]. 地理科學(xué)， 2012， 32(5):621-627．

[24] 張鐿鋰， 祁 威， 周才平， 等. 青藏高原高寒草地凈初級生產(chǎn)力(NPP)時(shí)空分異[J]. 地理學(xué)報(bào)， 2013， 68(9):1197-1211．

[25] CURRIE J D， PAQUIN V .Large-scale biogeographical patterns of species richness of trees[J]. Nature， 1987， 329: 326-327．

[26] 江 源， 黃曉霞， 劉全儒， 等. 五臺山高山、亞高山草甸植物多樣性格局分析[J]. 北京師范大學(xué)學(xué)報(bào)(自然科學(xué)版)， 2009， 45(1):91-95.

Evaluation of vegetation net primary productivity on uninhabited islands in southern Miaodao Archipelago in China

WANG Enkang， CHI Yuan

(The First Institute of Oceanography， State Oceanic Administration， Qingdao， Shandong 266061)

There are important strategic functions and ecological values in uninhabited islands， and the evaluation on vegetation net primary productivity (NPP) is of great significances for maintaining the stability and vitality island ecosystem. Five uninhabited islands in southern Miaodao Archipelago， typical islands in North China， were used as the study area. The NPP and spatial distribution of these islands were analyzed based on the method of Remote Sensing. The results indicated that the total amount of NPP of the study area reached 134.22 tC·a-1， with the mean density of 274.13 gC·m-2·a-1， and the NPP differed greatly in different seasons; the mean densities of NPP of the five islands showed the descending order as Niutuozi Island， Yangtuozi Island， Tanglang Island， Nantuozi Island and Danglang Island. The land cover types included broad-leaf forest， coniferous forest， grassland， construction land and bare land in the descending order of the mean density of NPP， and the mean densities of NPP on different islands were negatively associated with the proportion of bare land; the mean densities of NPP in the study area were positively associated with altitude， slope， aspect and normalized difference vegetation index (NDVI). Pinus thunbergii and Robinia pseudoacacia are ideal tree species for island plantation in North China， and continuous plantation is one of important measures to improve the NPP of uninhabited islands．

uninhabited islands; net primary productivity; spatial distribution; island plantation

2016-10-08.

國家海洋局海洋公益性行業(yè)科研專項(xiàng)經(jīng)費(fèi)項(xiàng)目(201505012); 中央級公益性科研院所基本科研業(yè)務(wù)費(fèi)專項(xiàng)資金項(xiàng)目(2015G13)．

1000-1190(2017)02-0195-08

P76; X87

A

*通訊聯(lián)系人. E-mail: chiyuan1988@163.com.