傅杰超, 李偉芳,趙柯,王奇,胡慧

象山港海岸帶高功能景觀退化及環(huán)境效應研究

傅杰超1,2,3, 李偉芳2,3,4,*,趙柯1,王奇1,2,3,胡慧1,2,3

1. 寧波大學地理與空間信息技術系, 寧波 315211 2. 寧波大學陸海國土空間利用與治理研究中心, 寧波 315211 3. 浙江省新型重點專業(yè)智庫寧波大學東海研究院, 寧波 315211 4. 寧波大學土木與環(huán)境工程學院, 寧波 315211

以象山港為研究區(qū), 基于1995年、2005年和2015年的土地利用數(shù)據(jù), 采用空間統(tǒng)計、標準差橢圓、景觀格局指數(shù)等方法分析象山港海岸帶高功能景觀時空演變特征, 并從生態(tài)服務價值和污染物排放變化角度探討其環(huán)境效應。結果表明: (1)象山港海岸帶高功能景觀總規(guī)模呈加速減小趨勢; (2)林地、耕地景觀重心向西南轉移, 景觀破碎化加劇, 水域、灘涂景觀分別向東南、東北擴張, 破碎度總體上呈減小趨勢; (3)受高功能景觀規(guī)模和空間格局變化的影響, 象山港海岸帶高功能景觀生態(tài)服務總價值降低, 研究區(qū)污染物排放總量增加, 不利于海岸帶可持續(xù)發(fā)展。

高功能景觀; 空間格局; 環(huán)境效應; 海岸帶; 象山港

0 前言

海岸帶具有巨大的生態(tài)價值和經(jīng)濟價值[1], 因受海陸作用以及人類開發(fā)活動的共同影響[2], 海岸帶生態(tài)環(huán)境變化顯著[3], 成為我國一大生態(tài)脆弱帶和敏感區(qū)[4], 引發(fā)國內(nèi)外專家學者[5-6]的廣泛關注。景觀格局強烈影響生態(tài)過程[7], 景觀的組成及其變化會引起生態(tài)系統(tǒng)物質(zhì)和能量的變化[8], 因此對景觀組成及其空間格局演變的分析是探討生境變化的基礎。目前我國有關景觀分類系統(tǒng)[9]、景觀格局分析方法[10]的研究已經(jīng)較為成熟, 后者主要有空間統(tǒng)計分析、轉移矩陣分析和景觀格局指數(shù)分析[11]等。在景觀格局演變的生態(tài)環(huán)境效應研究方面, 既有對生態(tài)環(huán)境效應的綜合評價研究[12], 也有從大氣環(huán)境效應[13]、水環(huán)境效應[14]、生態(tài)服務效應[15]等方面進行的分維度研究。且隨著海岸帶的開發(fā)利用, 已有越來越多的學者以海岸帶為背景研究景觀格局變化及其環(huán)境效應, 如海岸帶景觀格局對生境質(zhì)量、生態(tài)安全的影響[16-17]等。從研究對象來看, 現(xiàn)有研究大多以海岸帶所有景觀類型或單一景觀類型[18]為研究客體, 很少有研究聚焦相對具有較高生態(tài)功能的景觀群體。因此, 本研究以高生態(tài)功能景觀(指受到人類活動干預較少, 其生產(chǎn)能力和自我調(diào)節(jié)能力都比較強的景觀類型, 如林地、水域等)為切口, 旨在探討高功能景觀的時空演變特征及其與低功能景觀之間的轉化關系。

象山港區(qū)域是全國海洋生態(tài)文明示范區(qū), 生態(tài)功能突出, 海岸帶景觀分布以高功能景觀為主導, 但近年來隨著城鎮(zhèn)化和工業(yè)化的發(fā)展, 特別是在海洋經(jīng)濟乃至灣區(qū)經(jīng)濟的驅(qū)動下, 高功能景觀出現(xiàn)退化趨勢, 故以象山港作為研究區(qū)研究海岸帶高功能景觀退化特征及其環(huán)境效應, 對指導藍色生態(tài)海灣建設、促進海岸帶保護與合理利用具有重要意義。

1 研究區(qū)概況與數(shù)據(jù)來源

1.1 研究區(qū)概況

象山港是位于我國浙江省寧波市東南部的狹長型半封閉式深水港灣, 海岸線曲折且水域?qū)掗?。地理位置處?9°24′—29°46′N, 121°25′—122°00′E之間, 東側為舟山群島, 南鄰三門灣, 北靠杭州灣。海岸線長達406 km, 其中大陸岸線占297 km。象山港跨越鄞州、北侖、象山、奉化和寧海縣(市、區(qū)), 本文以象山港沿岸的鄉(xiāng)鎮(zhèn)行政區(qū)為研究區(qū)(圖1), 總面積1682.30 km2。象山港海岸帶土地利用類型以林地為主, 占研究區(qū)總面積60%左右, 耕地大多分布在地勢比較平緩的中部地區(qū), 建設用地、水域、灘涂和養(yǎng)殖用地主要分布在沿海區(qū)域, 未利用地多靠近水域, 分布較零散。依托豐富的港口資源、水產(chǎn)資源、旅游資源等優(yōu)勢條件, 象山港經(jīng)濟蓬勃發(fā)展, 但同時也導致了高功能景觀類型逐漸向低功能景觀類型退化, 生態(tài)環(huán)境壓力加劇。

1.2 數(shù)據(jù)來源

本研究所需的象山港海岸帶土地利用數(shù)據(jù)提取自地理國情監(jiān)測云平臺(http://www.dsac.cn/)1995、2005、2015年3期浙江省土地利用現(xiàn)狀數(shù)據(jù)。結合本文的研究目的與《土地利用現(xiàn)狀分類標準》(GB/T21010—2007), 將研究區(qū)的土地利用類型合并為林地、耕地、水域、灘涂、建設用地、養(yǎng)殖用地和未利用地七大類(如圖2), 其具體劃分依據(jù)和標準參考相關研究文獻[11]。同時根據(jù)研究需要, 將林地、耕地、水域和灘涂歸為高功能景觀, 將建設用地、養(yǎng)殖用地和未利用地歸為低功能景觀。社會經(jīng)濟數(shù)據(jù)來自1995、2005、2015年《寧波統(tǒng)計年鑒》以及象山港區(qū)域空間保護和利用規(guī)劃報告、象山港海洋環(huán)境公報、寧波市環(huán)境狀況公報和相關部門的統(tǒng)計數(shù)據(jù)。

圖1 研究區(qū)位置

Figure 1 The location of the study area

圖2 1995年、2005年、2015年象山港海岸帶土地利用現(xiàn)狀圖

Figure 2 Status of land use in the coastal zone of Xiangshan Bay in 1995,2005 and 2015

2 研究方法

2.1 景觀格局指數(shù)

本文選取斑塊數(shù)量()、平均斑塊面積()、斑塊密度()、邊界密度()、形狀指數(shù)()、聚集度指數(shù)()6項景觀格局指數(shù)對象山港海岸帶高功能景觀格局進行分析, 各指標計算公式及生態(tài)含義參考相關研究文獻[19–20]。

2.2 生態(tài)服務價值估算

高功能景觀生態(tài)服務價值計算公式[21]如下:

其中,為研究區(qū)高功能景觀生態(tài)服務價值總量,為高功能景觀類型,V為高功能景觀生態(tài)服務價值系數(shù), 即第類高功能景觀單位面積生態(tài)服務價值,A為研究區(qū)第類高功能景觀的面積。

本研究對象山港海岸帶高功能景觀單位面積生態(tài)服務價值V的估算, 采用謝高地等[22]學者的評估方法, 將林地、耕地、水域、灘涂分別與其評估體系中的森林、農(nóng)田、河流/湖泊、濕地對應, 據(jù)此得到各高功能景觀單位面積生態(tài)服務價值當量, 并根據(jù)1995—2015年象山港區(qū)域糧食作物的播種面積和產(chǎn)值, 計算出單個生態(tài)服務價值當量因子的經(jīng)濟價值為1821.78元·hm–2, 繼而得到修正后的高功能景觀生態(tài)服務價值系數(shù)(表1)。

表1 象山港海岸帶高功能景觀生態(tài)服務價值系數(shù) (元·hm–2)

2.3 不同生態(tài)功能斑塊的產(chǎn)污分析

本研究對經(jīng)營性景觀地類所產(chǎn)生污染物的排放強度的估算, 采用排放系數(shù)法[23], 計算公式如下:

Q=C*?A(2)

其中,Q表示某一地類變化所產(chǎn)生的污染物排放量,C表示該地類的污染物排放系數(shù),表示由高功能景觀轉化成該種景觀的面積增量。

為了進行直觀的景觀變化產(chǎn)生的污染物排放量的比較, 本研究排放系數(shù)以地類單位面積污染物排放量為主, 計算1995—2005年和2005—2015年兩個時段由地類面積的改變引起的污染物排放量的大小。選取象山港海岸帶有污染物排放能力的耕地及低功能景觀類型的養(yǎng)殖用地和建設用地, 通過由高功能景觀類型的面積轉化, 計算出生境轉化成經(jīng)營性地類所帶來的污染物排放量的變化。

2.3.1 耕地污染物排放系數(shù)

農(nóng)業(yè)面源污染對海岸帶水質(zhì)和土壤的影響比較大, 尤其是發(fā)生暴雨徑流和農(nóng)田瀝水時, 水流進入河道和地下水, 最終匯入海洋[24]。計算耕地污染物排放系數(shù)的公式如下:

其中,表示耕地污染物排放系數(shù),表示施肥總量,表示耕地面積,表示化肥流失系數(shù), 本文參考邱君[25]學者的研究取0.5。

2.3.2 養(yǎng)殖用地污染物排放系數(shù)

養(yǎng)殖用地投放的飼料中含有大量的N、P、COD等物質(zhì), 其作為主要的污染源, 污染物排放量的大小會直接影響到象山港海岸帶的環(huán)境質(zhì)量。本文以無機氮作為養(yǎng)殖用地污染物的代表, 參考蔡惠文[26]學者對象山港養(yǎng)殖排污廢水中N含量的測算值計算污染物排放系數(shù), 并結合歷年象山港海域海水中無機氮含量變化趨勢, 最終以2004年的計算值表示1995—2015年養(yǎng)殖用地污染物排放系數(shù)。

2.3.3 建設用地生活污水排放系數(shù)

生活污水的排放量通常與人口數(shù)量有直接的聯(lián)系, 本文以生活用水量在扣除損耗[27]后建立與住宅建筑面積的關系, 計算生活污水排放系數(shù), 計算公式如下:

其中,為生活污水排放系數(shù),表示生活污水排放總量,表示住宅面積,表示生活用水總量, a表示生活用水損失率, 一般取值在0.15—0.20之間, 本文a取值為中間值0.18。

2.3.4 建設用地工業(yè)廢水/廢氣排放系數(shù)

建設用地工業(yè)廢水/廢氣排放系數(shù)計算公式如下:

其中,表示工業(yè)廢水/廢氣排放系數(shù),表示工業(yè)廢水/廢氣排放量,表示工礦區(qū)面積。

3 結果與分析

3.1 高功能景觀斑塊退化特征

3.1.1 高功能景觀規(guī)模變化

1995年至2015年期間, 象山港海岸帶高生態(tài)功能的景觀類型總體上逐漸向低生態(tài)功能的景觀類型退化。結合表2和圖3可知, 1995年高功能景觀總面積為1516.92 km2, 2015年面積為1430.37 km2, 平均每年減少4.328 km2, 減少度為5.71%。其中, 林地在高功能景觀類型中占地面積最廣, 呈持續(xù)緩慢減少的整體趨勢。1995—2005年間和2005—2015年間, 林地面積分別減少13.87 km2和13.36 km2, 兩個10年間林地面積的減少速率基本相等。耕地是象山港海岸帶景觀面積減少最快的類型, 且1995—2015年間呈現(xiàn)加速減少的趨勢, 林地面積減少速率在各類海岸帶景觀類型中僅次于耕地。水域面積持續(xù)增加, 灘涂面積呈現(xiàn)先增長后減少的趨勢, 但總體變化較不明顯。

從高功能景觀退化的方向來看(表3), 1995— 2005年高功能景觀林地向低功能景觀建設用地和養(yǎng)殖用地轉化, 且絕大部分都轉為建設用地, 2005— 2015年林地倍速向建設用地轉化, 人類開發(fā)活動對林地的影響增加。1995—2015年耕地基本轉為建設用地, 這主要是由于象山港沿海城鎮(zhèn)經(jīng)濟快速發(fā)展, 快速城市化占用了大量的耕地資源, 也有少數(shù)耕地轉為養(yǎng)殖用地和未利用地, 前者出于養(yǎng)殖業(yè)的發(fā)展, 后者主要因為耕地使用不當或保護不力使得耕地質(zhì)量下降最終導致耕地退化。水域在1995—2015年期間未有向低功能景觀轉化的趨勢。灘涂在1995—2005年以低功能景觀(養(yǎng)殖用地、未利用地)轉入為主, 在2005—2015年則以向低功能景觀(建設用地)轉出為主導形式, 可見建設用地擴張對生態(tài)用地的侵占與威脅。

表2 象山港海岸帶高功能景觀面積(1995—2015)

圖3 象山港海岸帶景觀類型變化(1995—2015)

Figure 3 Landscape type changes in the coastal zone of Xiangshan Bay(1995-2015)

3.1.2 高功能景觀空間格局變化

使用標準差橢圓空間計量方法[28], 計算象山港海岸帶各景觀類型標準差橢圓參數(shù)值并將標準差橢圓動態(tài)變化可視化(圖4), 通過不同時期景觀重心的轉移分析1995—2015年研究區(qū)內(nèi)高功能景觀分布的變化情況。受象山港地理位置及其港灣形態(tài)的影響, 研究區(qū)內(nèi)各類景觀的標準差橢圓長軸均呈現(xiàn)東北—西南走向, 短軸則呈現(xiàn)西北—東南走向。1995—2015年林地短軸減少842.53 m, 長軸增加2137.79 m, 旋轉角由1995年的53.4°偏轉為2015年的53.2°, 標準差橢圓向逆時針方向偏轉0.2°, 總體上林地景觀方向性增強, 其東北—西南走向的分布態(tài)勢更加明顯, 重心呈現(xiàn)向西南轉移的趨勢。耕地短軸減少664.73 m, 長軸增加1716.94 m, 方向性同林地一樣增強, 但標準差橢圓向順時針方向偏轉0.1°, 重心亦向西南方向轉移。水域短軸伸長, 長軸縮短, 方向性減弱, 向與分布態(tài)勢相反的方向(西北—東南方向)擴張, 橢圓圈層擺動幅度較大, 景觀重心向東南偏移。灘涂長短軸增減幅度最小, 東北—西南走向的方向性變化不明顯, 橢圓向順時針方向偏轉0.4°, 景觀重心向東北方向轉移。同時, 對比低功能景觀的標準差橢圓動態(tài)變化, 可以發(fā)現(xiàn)建設用地向西南和東南方向擴展的趨勢比較明顯, 這主要是由象山港西南部和東南部建設用地的擴張速度較快所導致的, 而建設用地在擴張過程中大量占用林地和耕地, 正是導致林地和耕地景觀重心向西南方向轉移的主要原因。

根據(jù)1995年、2005年及2015年3個時期的土地利用數(shù)據(jù), 結合Fragstats 4軟件, 統(tǒng)計出3個時期各高功能景觀的景觀格局指數(shù)(圖5)。從中可以看出, 近20年象山港海岸帶各高功能景觀類型空間格局發(fā)生了較大變化。林地、耕地斑塊數(shù)量在1995—2015年期間逐漸增加, 斑塊類型面積卻逐年減少, 說明這兩類景觀在空間分布上趨于分散, 破碎化程度加深。其次, 林地、耕地的平均斑塊面積有所減小, 大斑塊被破碎成小斑塊, 斑塊密度和形狀指數(shù)有所增加, 聚集度下降, 進一步證明了林地和耕地日益嚴重的景觀破碎性。其中, 耕地景觀邊界密度先增大后減小, 形狀指數(shù)的增加趨勢在2005—2015年有所減緩, 這是由于耕地在被建設用地占用過程中連續(xù)性下降, 斑塊形狀趨于不規(guī)則, 但后期出于對耕地的保護政策, 其破碎化進程得到了一定程度的控制。水域景觀斑塊數(shù)量、平均斑塊面積均呈遞增趨勢, 可見象山港海岸帶河流、湖泊、水庫等水域面積容量逐年擴大, 斑塊密度和邊界密度這兩項衡量破碎化程度的指標的計算值雖然有所增加, 但結合聚集度和平均斑塊面積來看, 水域景觀的破碎化水平還是呈下降態(tài)勢, 因為小斑塊數(shù)量雖然略有增加, 但大面積斑塊的比例也在上升, 斑塊趨于集聚化。灘涂景觀的各項景觀格局指數(shù)變化幅度較小, 空間格局變化較不明顯。1995—2005年灘涂總面積和平均斑塊面積都呈增加趨勢, 這是由灘涂自然淤漲所致, 但后期人類對灘涂的圍墾利用規(guī)模超過了其自然淤漲速度, 導致灘涂總面積和平均斑塊面積減小, 斑塊密度、邊界密度和形狀指數(shù)增大, 聚集度下降, 即灘涂景觀破碎度增加, 形狀趨向復雜, 分布趨向離散。

表3 象山港海岸帶高功能景觀轉換面積變化(1995—2015)(km2)

注: 表中的正值代表低功能景觀轉成高功能景觀; 負值代表高功能景觀轉成低功能景觀。

圖4 象山港海岸帶各類型景觀標準差橢圓動態(tài)變化(1995—2015)

Figure 4 Dynamic variation of the standard deviation of various types of landscapes in the coastal zone of Xiangshan Bay(1995-2015)

圖5 象山港海岸帶高功能景觀各類型景觀格局指數(shù)(1995—2015)

Figure 5 Landscape pattern index of various types of high functional landscape in Xiangshan Bay coastal zone(1995-2015)

3.2 高功能景觀退化的環(huán)境效應

3.2.1 生態(tài)服務價值變化

景觀生態(tài)系統(tǒng)具有氣體調(diào)節(jié)、氣候調(diào)節(jié)、環(huán)境凈化、土壤保持、生物多樣性維持等生態(tài)服務功能[29], 景觀類型不同其單位面積生態(tài)系統(tǒng)服務價值也不同, 根據(jù)公式(1)計算得到象山港海岸帶高功能景觀生態(tài)系統(tǒng)服務價值變化情況(圖6)。從圖中可以看出, 1995—2015年象山港海岸帶林地、耕地景觀生態(tài)系統(tǒng)服務價值持續(xù)降低, 其根本原因是大量林地、耕地景觀轉出為建設用地等低功能景觀, 這也是導致象山港海岸帶高功能景觀生態(tài)系統(tǒng)服務價值總量從61.36億元持續(xù)減少至60.73億元的主要原因, 即使生態(tài)服務價值系數(shù)較高的水域和灘涂在這20年間生態(tài)系統(tǒng)服務價值總體上有所增加, 也改變不了高功能景觀生態(tài)服務總價值持續(xù)降低的趨勢。

另一方面, 象山港海岸帶景觀格局特別是景觀破碎度的變化, 也會對局部地區(qū)生態(tài)環(huán)境及其生態(tài)價值造成影響。在景觀中, 斑塊面積的大小可以影響到物種的豐富度、數(shù)量、食物鏈及生產(chǎn)能力。礦物質(zhì)養(yǎng)分和能量的總量與斑塊面積成正比, 所以面積較大的斑塊所含有的礦物質(zhì)養(yǎng)分和能量比面積小的斑塊多, 其次所影響的物種豐富度、數(shù)量、生產(chǎn)能力等也隨著斑塊面積的增大而增加[30]。斑塊面積越大, 生態(tài)環(huán)境多樣性就越大。而面積較小的斑塊物種少, 食物鏈短, 生態(tài)環(huán)境質(zhì)量較差, 物種滅絕的概率就會偏高[31]。高功能景觀斑塊相對來說具有較好的生態(tài)能力, 適宜物種的棲息、生存和繁殖, 是各種生物生存比較理想的區(qū)域。研究區(qū)內(nèi)林地和耕地等高功能景觀的平均斑塊面積減少, 景觀破碎度增加, 尤其是當部分林地和耕地轉化為建設用地時, 建設用地會隔斷物種之間的交流, 降低景觀之間的連通性, 不利于生物多樣性的維持。

圖6 象山港海岸帶高功能景觀生態(tài)服務價值變化(1995—2015)

Figure 6 Change of ecosystem service value of high functional landscape in Xiangshan Bay coastal zone (1995-2015)

3.2.2 污染物排放變化

高功能景觀類型在退化成低功能景觀類型的過程中, 會產(chǎn)生相應的物質(zhì)交換, 不僅自身對污染物的自凈能力有所降低, 同時還會輸送污染物到近海岸區(qū)域。以沼澤為例, 當沼澤轉化為農(nóng)田時, 所產(chǎn)生的農(nóng)田面源污染負荷將會對近海岸區(qū)域產(chǎn)生一定的影響。通過公式(2)可以估算出高功能景觀在轉變過程中成為污染物排放地類對環(huán)境所排放的污染物變化(表4), 表中正值表示污染物排放量有所增加, 負值表示污染物排放量有所減少。

1995—2015年, 由于面積減少, 耕地流失到環(huán)境中的化肥量減少了1497 t, 緩解了由N、P元素引起的水體富營養(yǎng)化與硝酸鹽污染。養(yǎng)殖用地N、P、COD等污染物排放量先增加后減少, 近20年期間總計減少了11.8 t無機氮排放量, 這是依據(jù)區(qū)域保護規(guī)劃嚴格控制養(yǎng)殖數(shù)量的結果。建設用地所排放的生活污水、工業(yè)廢水和工業(yè)廢氣分別增加了18.49×106t、14.70×106t和11.89×105萬標立米, 對周邊環(huán)境影響巨大。綜合來看, 高功能景觀在向低功能景觀轉換的過程中, 所帶來的污染排放量在總體上仍呈現(xiàn)增長趨勢。污染物會直接或間接影響象山港海岸帶的大氣環(huán)境、土壤環(huán)境、水環(huán)境、海洋環(huán)境等, 降低環(huán)境質(zhì)量, 加劇海岸帶的生態(tài)風險和環(huán)境壓力。

4 結論與討論

本文重點揭示了象山港海岸帶高功能景觀的退化特征, 并從生態(tài)服務價值和污染物排放角度闡述其對生態(tài)環(huán)境的影響, 可為海岸帶景觀格局優(yōu)化及可持續(xù)發(fā)展提供決策支持。研究表明: (1)象山港海岸帶高功能景觀以林地和耕地為主, 總規(guī)模呈加速減小趨勢, 其中耕地減少速度最快。(2)1995—2015年象山港海岸帶水域景觀破碎化水平呈下降態(tài)勢, 灘涂景觀破碎化程度先減小后增加, 林地和耕地景觀由于保護不力破碎化加劇, 再加上建設用地的無序擴張, 林地和耕地趨于分散化, 景觀連通性下降。(3)象山港海岸帶高功能景觀在退化過程中, 不僅會降低高功能景觀生態(tài)系統(tǒng)服務總價值, 同時轉化后的養(yǎng)殖用地、建設用地也會排放出大量污染物, 給海岸帶生態(tài)環(huán)境帶來負面影響。

表4 象山港海岸帶高功能景觀退化排污估算(1995—2015)

注: 耕地屬于高功能景觀, 此表計算耕地向低功能景觀轉出過程中的產(chǎn)污變化。表中除建設用地的工業(yè)廢氣排放量以萬標立米計量外, 其它地類的污染物排放量均以t計量。

本研究也存在一定的局限性, 有關景觀規(guī)模和空間格局演變對生態(tài)環(huán)境的影響還有待進一步深入探討, 如何綜合當量因子法、專家評估法、價值量評價法、償付意愿法等[32]改進生態(tài)價值評估模型, 以及通過環(huán)境科學和數(shù)學等學科知識的結合選擇更優(yōu)的排污量計算模型使結果更加精確, 是今后研究需要突破的重點和難點。

[1] 楊魯永. 連云港海岸帶土地利用變化及其生態(tài)承載力研究[D]. 南京: 南京大學, 2014.

[2] 李偉芳, 陳陽, 馬仁鋒, 等. 發(fā)展?jié)摿σ暯堑暮０稁恋乩媚Ｊ健院贾轂衬习稙槔齕J]. 地理研究, 2016, 35(6): 1061–1073.

[3] DONEY S C. The growing human footprint on coastal and open–ocean biogeochemistry[J]. Science, 2010, 328(5985): 1512–1516.

[4] 駱永明. 中國海岸帶可持續(xù)發(fā)展中的生態(tài)環(huán)境問題與海岸科學發(fā)展[J]. 中國科學院院刊, 2016, 31(10): 1133– 1142.

[5] IVAN V, ELIZABETH E , JAVIER L, et al. Tropical land–sea couplings: Role of watershed deforestation, mangrove estuary processing, and marine inputs on N fluxes in coastal Pacific Panama[J]. Science of the Total Environment, 2018, 630: 126–140.

[6] 宋曉帥, 于開寧, 吳振, 等. 萊州灣海岸帶工程地質(zhì)環(huán)境質(zhì)量分區(qū)[J]. 海洋地質(zhì)與第四紀地質(zhì), 2019, 39(2): 79–89.

[7] 劉怡娜, 孔令橋, 肖燚, 等. 長江流域景觀格局與生態(tài)系統(tǒng)水質(zhì)凈化服務的關系[J]. 生態(tài)學報, 2019, 39(3): 844– 852.

[8] TURNER M G, GARDNER R H, O'NEILL R V. Landscape ecology in theory and practice[M]. New York: Springer, 2001.

[9] 馮異星, 羅格平, 周德成, 等. 近50a土地利用變化對干旱區(qū)典型流域景觀格局的影響——以新疆瑪納斯河流域為例[J]. 生態(tài)學報, 2010, 30(16): 4295–4305.

[10] 何改麗, 史小麗, 李加林, 等. 快速城市化背景下寧波北侖綠地系統(tǒng)演化及生態(tài)服務響應[J]. 寧波大學學報(理工版), 2017, 30(5): 89–96.

[11] 徐諒慧. 岸線開發(fā)影響下的浙江省海岸類型及景觀演化研究[D]. 寧波: 寧波大學, 2015.

[12] 李保杰, 顧和和, 紀亞洲. 礦區(qū)土地復墾景觀格局變化和生態(tài)效應[J]. 農(nóng)業(yè)工程學報, 2012, 28(3): 251–256.

[13] 苗世光, 王曉云, 蔣維楣, 等. 城市規(guī)劃中綠地布局對氣象環(huán)境的影響——以成都城市綠地規(guī)劃方案為例[J]. 城市規(guī)劃, 2013, 37(6): 41–46.

[14] 馮源嵩. 貴陽市核心區(qū)景觀格局演變的水環(huán)境效應研究[D]. 重慶: 西南大學, 2016.

[15] 鄒月. 西安市土地利用景觀格局變化的生態(tài)系統(tǒng)服務價值響應及優(yōu)化研究[D]. 西安: 陜西師范大學, 2018.

[16] 褚琳, 黃翀, 劉慶生, 等. 2000–2010年遼寧省海岸帶景觀格局與生境質(zhì)量變化研究[J]. 資源科學, 2015, 37(10): 1962–1972.

[17] 黃寧, 楊綿海, 林志蘭, 等. 廈門市海岸帶景觀格局變化及其對生態(tài)安全的影響[J]. 生態(tài)學雜志, 2012, 31(12): 3193–3202.

[18] 周云軒, 田波, 黃穎, 等. 我國海岸帶濕地生態(tài)系統(tǒng)退化成因及其對策[J]. 中國科學院院刊, 2016, 31(10): 1157– 1166.

[19] 郭曉妮. 康?？h景觀格局變化及其生態(tài)效應分析[D]. 北京: 首都師范大學, 2009.

[20] 鄒月, 周忠學. 西安市景觀格局演變對生態(tài)系統(tǒng)服務價值的影響[J]. 應用生態(tài)學報, 2017, 28(8): 2629–2639.

[21] 魏慧, 趙文武, 張驍, 等. 基于土地利用變化的區(qū)域生態(tài)系統(tǒng)服務價值評價——以山東省德州市為例[J]. 生態(tài)學報, 2017, 37(11): 3830–3839.

[22] 謝高地, 甄霖, 魯春霞, 等. 一個基于專家知識的生態(tài)系統(tǒng)服務價值化方法[J]. 自然資源學報, 2008, 23(5): 911– 919.

[23] 萬峻. 渤海典型海岸帶景觀格局變化及其生態(tài)環(huán)境效應研究[D]. 北京: 中國環(huán)境科學研究院, 2004.

[24] 陳明劍. 海洋功能區(qū)劃中的空間關系模型及其GIS實現(xiàn)(以萊州灣為例)[D]. 青島: 中國海洋大學, 2003.

[25] 邱君. 我國化肥施用對水污染的影響及其調(diào)控措施[J]. 農(nóng)業(yè)經(jīng)濟問題, 2007(S1): 75–80.

[26] 蔡惠文. 象山港養(yǎng)殖環(huán)境容量研究[D]. 青島: 中國海洋大學, 2004.

[27] 方子云. 水資源保護工作手冊[M]. 南京: 河海大學出版社, 1988.

[28] 陳陽, 李偉芳, 任麗燕, 等. 空間統(tǒng)計視角下的農(nóng)村居民點分布變化及驅(qū)動因素分析——以鄞州區(qū)濱海平原為例[J]. 資源科學, 2014, 36(11): 2273–2281.

[29] 謝高地, 張彩霞, 張雷明, 等. 基于單位面積價值當量因子的生態(tài)系統(tǒng)服務價值化方法改進[J]. 自然資源學報, 2015, 30(8): 1243–1254.

[30] 傅伯杰, 陳利頂, 馬克明, 等. 景觀生態(tài)學原理及應用[M]. 北京: 中國環(huán)境科學出版社, 2005.

[31] 肖篤寧, 布仁倉, 李秀珍. 生態(tài)空間理論與景觀異質(zhì)性[J]. 生態(tài)學報, 1997, 17(5): 453–461.

[32] 束邱愷, 高永年, 劉友兆, 等. 江蘇沿海地區(qū)土地利用生態(tài)價值測算評估[J]. 地球信息科學學報, 2016, 18(6): 787–796.

Study on high functional landscape degradation and environmental effects in coastal zone of Xiangshan Bay

FU Jiechao1,2,3, LI Weifang2,3,4,*, ZHAO Ke1, WANG Qi1,2,3, HU Hui1,2,3

1. Department of Geography and Spatial Information Techniques, Ningbo University,Ningbo315211,China 2. Center for Land and Marine Spatial Utilization and Governance Research , Ningbo University, Ningbo 315211, China 3. Institute of East China Sea, Ningbo University, Ningbo 315211, China 4. Faculty of Civil and Environmental Engineering, Ningbo University,Ningbo315211,China

Taking Xiangshan Bay as the research area, based on the land use data of Xiangshan Bay coastal zone in 1995, 2005 and 2015, using spatial statistics, standard deviation ellipse, landscape pattern indexes and other methods, we analyzed the temporal and spatial evolution characteristics of the high functional landscape of Xiangshan Bay coastal zone, and then explored its environmental effects from the perspective of ecosystem service value and pollution emission. The results showed that: 1) The total scale of high functional landscape in the coastal zone of Xiangshan Bay showed a trend of accelerated reduction. 2) The centers of gravity of forest land and cultivated land shifted to the southwest, and the landscape fragmentations of forest land and cultivated land were intensified. The waters and tidal flats expanded toward the southeast and northeast respectively, and the degrees of fragmentation of them were generally decreasing. 3) Affected by changes in the scale and spatial pattern of high functional landscapes, the total value of ecosystem services of high functional landscape in Xiangshan Bay coastal zone was reduced, and the total discharge amount of pollutants in study area was increased, which was not conducive to the sustainable development of coastal zone.

high functional landscape; spatial pattern; environmental effect; coastal zone; Xiangshan Bay

10.14108/j.cnki.1008-8873.2020.02.020

F301.24; X171.1

A

1008-8873(2020)02-166-09

2019-07-22;

2019-10-22

NSFC-浙江兩化融合聯(lián)合基金資助項目(U1609203)

傅杰超(1995—), 女, 浙江紹興人,碩士研究生, 主要從事區(qū)域研究與濱海城鎮(zhèn)發(fā)展研究, E-mail: 532480471@qq.com

李偉芳(1964—), 男, 浙江寧波人,教授, 主要從事土地資源管理和土地利用規(guī)劃研究, E-mail: liweifang@nbu.edu.cn

傅杰超, 李偉芳, 趙柯, 等. 象山港海岸帶高功能景觀退化及環(huán)境效應研究[J]. 生態(tài)科學, 2020, 39(2): 166–174.

FU Jiechao, LI Weifang, ZHAO Ke, et al. Study on high functional landscape degradation and environmental effects in coastal zone of Xiangshan Bay[J]. Ecological Science, 2020, 39(2): 166–174.