長(zhǎng)三角海洋資源環(huán)境一體化治理策略研究
——基于海洋生態(tài)系統(tǒng)服務(wù)價(jià)值影響分析

賀義雄 ， 宋偉鳴 楊 帆

(1. 浙江海洋大學(xué) 經(jīng)濟(jì)與管理學(xué)院， 浙江 舟山 316022； 2. 南方海洋科學(xué)與工程廣東省實(shí)驗(yàn)室(珠海)， 廣東珠海 519000)

長(zhǎng)江三角洲(簡(jiǎn)稱長(zhǎng)三角)區(qū)域一體化始于1982年，至今已走過了將近40年的歷程。從1982年最初設(shè)想的上海、南京、寧波、蘇州與杭州這一范圍逐漸擴(kuò)容為當(dāng)前江浙滬皖“三省一市”全域， 一體化的發(fā)展取得了顯著成績(jī)， 從基礎(chǔ)設(shè)施和旅游開始， 已逐漸拓展至產(chǎn)業(yè)、社會(huì)、環(huán)保、文化等諸多領(lǐng)域。目前， 關(guān)于長(zhǎng)三角一體化的研究在長(zhǎng)三角一體化發(fā)展水平或趨勢(shì)的測(cè)度[1-3]、長(zhǎng)三角一體化的動(dòng)力機(jī)制與制約因素[4-7]、長(zhǎng)三角一體化的國(guó)際經(jīng)驗(yàn)借鑒[8-10]、長(zhǎng)三角一體化的發(fā)展思路[11]、長(zhǎng)三角一體化的制度建設(shè)[12-14]、長(zhǎng)三角一體化領(lǐng)域發(fā)展的措施建議[15-16]等方面積累了豐富的成果。

另一方面， 長(zhǎng)江三角洲瀕臨黃海與東海， 海洋資源豐富， 海洋經(jīng)濟(jì)發(fā)達(dá)。長(zhǎng)三角的一體化不僅只是陸域的一體化， 還應(yīng)包括海域的一體化。目前， 這方面的研究主要集中在海洋經(jīng)濟(jì)一體化發(fā)展與海洋資源環(huán)境合作治理兩個(gè)方面。其中， 海洋經(jīng)濟(jì)一體化發(fā)展相關(guān)研究成果較全面[17-20]。關(guān)于海洋資源環(huán)境治理這一長(zhǎng)三角區(qū)域海洋一體化的關(guān)鍵組成部分[21]的研究則主要探討了治理機(jī)制構(gòu)建、完善與制度安排等方面的問題[22-24]， 而對(duì)于治理領(lǐng)域的選擇， 特別是重點(diǎn)治理領(lǐng)域的確定等影響最終治理效果的其他重要因素的研究較少涉及。

海洋治理是各主體共同管理人類對(duì)其開發(fā)利用實(shí)踐活動(dòng)的過程[25]。因此對(duì)長(zhǎng)三角海洋資源環(huán)境一體化治理的實(shí)質(zhì)應(yīng)是對(duì)該區(qū)域開發(fā)利用海洋資源環(huán)境的人類實(shí)踐活動(dòng)的一體化治理。另外， 生態(tài)系統(tǒng)服務(wù)是人類維持所依存的自然環(huán)境條件[26-27]， 因此可以說作為其在海洋領(lǐng)域具體解釋的海洋生態(tài)系統(tǒng)服務(wù)[28]的價(jià)值就是利用統(tǒng)一標(biāo)準(zhǔn)對(duì)與人類息息相關(guān)的海洋資源環(huán)境情況的反映?；诖?， 本文利用2009—2018年長(zhǎng)三角各沿海城市的面板數(shù)據(jù)， 在核算各城市的海洋生態(tài)系統(tǒng)服務(wù)價(jià)值基礎(chǔ)上， 通過靜態(tài)面板與動(dòng)態(tài)面板的綜合分析， 研究人類的行為活動(dòng)對(duì)區(qū)域海洋生態(tài)系統(tǒng)服務(wù)價(jià)值的影響， 得出人類的行為對(duì)長(zhǎng)三角海洋資源環(huán)境的具體作用情況， 據(jù)此提出長(zhǎng)三角海洋資源環(huán)境一體化治理的相關(guān)政策建議。試圖回答影響長(zhǎng)三角海洋資源環(huán)境狀況的主要因素是什么、一體化治理的首選領(lǐng)域是什么等問題。文章成果可以豐富長(zhǎng)三角一體化發(fā)展的理論研究，同時(shí)為促進(jìn)區(qū)域一體化更高層次的發(fā)展[29]提供參考借鑒。

1 研究區(qū)概況、指標(biāo)與研究方法、數(shù)據(jù)來源

1.1 區(qū)域概況與研究范圍

長(zhǎng)江三角洲是長(zhǎng)江入海之前的沖積平原， 具有活躍的經(jīng)濟(jì)發(fā)展態(tài)勢(shì)、高程度的對(duì)外開放水平及強(qiáng)有力的創(chuàng)新能力。依據(jù)國(guó)家《長(zhǎng)江三角洲區(qū)域一體化發(fā)展規(guī)劃綱要》(2019)， 一體化規(guī)劃范圍包括上海市、江蘇省、浙江省、安徽省的全部區(qū)域， 共35.8萬平方公里。區(qū)域內(nèi)， 沿海城市有上海市， 江蘇省連云港市、鹽城市、南通市， 浙江省嘉興市、寧波市、舟山市、臺(tái)州市、溫州市等9個(gè)城市1地理上， 浙江省杭州市與紹興市也是沿海城市， 但因?yàn)槠錄]有管轄海域， 因此本文未將這兩座城市納入研究范圍。進(jìn)一步地， 研究范圍內(nèi)的9個(gè)城市很多實(shí)際上只是某些區(qū)(或縣)為沿海區(qū)域， 并非整個(gè)城市。但考慮到我國(guó)海洋產(chǎn)業(yè)起步較晚， 涉海統(tǒng)計(jì)方面存在的不足之處， 以縣域(或以下)為研究單元很大程度上將會(huì)受制于數(shù)據(jù)的不可獲取等無法克服的研究瓶頸問題， 因此本文的研究單元設(shè)定為地級(jí)及以上城市。。

1.2 海洋生態(tài)系統(tǒng)服務(wù)價(jià)值核算指標(biāo)與方法

目前， 包括我國(guó)學(xué)者在內(nèi)的很多學(xué)者已對(duì)生態(tài)系統(tǒng)服務(wù)進(jìn)行了系統(tǒng)的分類[30-34]， 其中引用最多的、較權(quán)威的是聯(lián)合國(guó)的千年生態(tài)系統(tǒng)評(píng)估(the millennium ecosystem assessment， MA)中的分類體系， 其將生態(tài)系統(tǒng)服務(wù)分為供給、調(diào)節(jié)、文化和支持4大類， 同時(shí)每一類中又包含有不同的小項(xiàng)[35-36]。這一分類方式被廣泛認(rèn)可， 包括涉及海洋生態(tài)系統(tǒng)服務(wù)在內(nèi)的絕大部分研究均以上述4大類為基本依據(jù)對(duì)各種服務(wù)進(jìn)行劃分， 只是對(duì)相關(guān)定義和具體分類有一些調(diào)整。如在由聯(lián)合國(guó)環(huán)境規(guī)劃署(UNEP)主持的生態(tài)系統(tǒng)與生物多樣性經(jīng)濟(jì)學(xué)(the economics of ecosystems and biodiversity， TEEB)研究中， De Groot等[37]認(rèn)為支持服務(wù)應(yīng)屬于一種生態(tài)過程， 因此其使用生境(有些文獻(xiàn)翻譯成“棲息地”)服務(wù)進(jìn)行替代，從而構(gòu)成了新的生態(tài)系統(tǒng)服務(wù)分類體系。B?hnke-Henrichs等[38]遵照TEEB的分類， 基于選取的一定指標(biāo)， 劃分了海洋生態(tài)系統(tǒng)服務(wù)的具體種類。我國(guó)學(xué)者陳尚等[39]對(duì)海洋生態(tài)系統(tǒng)服務(wù)的分類也沿用了MA的劃分方式。

鑒于此， 并參考相關(guān)研究成果， 針對(duì)研究區(qū)的范圍2當(dāng)前， 我國(guó)各沿海地區(qū)管轄的海洋區(qū)域均是根據(jù)當(dāng)?shù)睾Ｑ箝_發(fā)利用的實(shí)際需要從陸地向海逐漸擴(kuò)展確定的[40]， 而這種擴(kuò)展還尚未超出渤海、黃海、東海、南海等近海區(qū)域[41]范圍， 因此本文是以近海為研究的區(qū)域范圍界定。， 同時(shí)考慮基礎(chǔ)資料可獲取性等實(shí)際情況， 本文將長(zhǎng)三角地區(qū)海洋生態(tài)系統(tǒng)的服務(wù)劃分為供給、調(diào)節(jié)、文化和支持等4大類13小項(xiàng)， 并整合市場(chǎng)法、成果參照法、替代成本法、影子工程法、費(fèi)用支出法、效益轉(zhuǎn)移法等方法對(duì)其價(jià)值進(jìn)行核算。具體如表1所示。

1.3 海洋生態(tài)系統(tǒng)服務(wù)價(jià)值影響分析指標(biāo)與方法

1.3.1 指標(biāo)選取

在海洋資源環(huán)境開發(fā)利用過程中， 人類污染和破壞了海洋生態(tài)環(huán)境。而由于這種污染與破壞具有極端的外部性特征[52]， 這部分損失就由其他社會(huì)成員共同分?jǐn)偂慕?jīng)濟(jì)學(xué)角度分析， 這是海洋生態(tài)環(huán)境問題產(chǎn)生的根源[53]。因此， 為了達(dá)到維護(hù)海洋生態(tài)平衡、實(shí)現(xiàn)海洋可持續(xù)開發(fā)這一海洋治理的第一目標(biāo)[54]， 應(yīng)該將解決海洋生態(tài)環(huán)境的污染和破壞作為首要問題?；谌啦ǖ萚55]對(duì)海洋生態(tài)環(huán)境污染和破壞原因的分析， 參考已有對(duì)海洋資源環(huán)境影響評(píng)價(jià)的研究中關(guān)于人類行為活動(dòng)的指標(biāo)設(shè)定， 結(jié)合長(zhǎng)三角區(qū)域主要海洋產(chǎn)業(yè)情況[17]， 從能夠反映對(duì)海洋資源環(huán)境造成直接影響的人類行為活動(dòng)角度， 依據(jù)科學(xué)性、目標(biāo)導(dǎo)向性等基本原則， 考慮數(shù)據(jù)的可獲得性， 本文選取海洋生態(tài)系統(tǒng)服務(wù)價(jià)值影響的評(píng)價(jià)指標(biāo)如表2所示。

表2 海洋生態(tài)系統(tǒng)服務(wù)價(jià)值影響評(píng)價(jià)指標(biāo)Tab. 2 Evaluation indexes of the impact of the marine ecosystem service value

1.3.2 模型建構(gòu)

目前， 關(guān)于生態(tài)系統(tǒng)服務(wù)價(jià)值影響分析所采用的方法主要有敏感性指數(shù)[56-57]、Pearson分析[58-59]、灰色關(guān)聯(lián)分析[60]等?？紤]本文的數(shù)據(jù)為面板數(shù)據(jù)， 因此參考何愛平等[61]的研究成果， 設(shè)立面板靜態(tài)方程與動(dòng)態(tài)方程如下。

?

?

上述公式中， OESV表示海洋生態(tài)系統(tǒng)服務(wù)價(jià)值，OESVi，y–1為其1階滯后項(xiàng)，δ為滯后項(xiàng)系數(shù)。WD表示廢水排放量， IE表示工業(yè)廢氣排放量， WWDDO表示廢水直排入海量， SWD表示固體廢物傾倒丟棄量， MY表示海水養(yǎng)殖產(chǎn)量， MA表示海水養(yǎng)殖面積， CO表示海洋捕撈量， RA表示圍填海面積，i和y分別表示具體的城市和年份。β為系數(shù)矩陣，εiy表示隨機(jī)擾動(dòng)項(xiàng)。

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

對(duì)近海海域初級(jí)生產(chǎn)力的確定， 考慮基于初級(jí)生產(chǎn)力垂向歸納模型(VGPM)的算法反演， 同時(shí)由于具有長(zhǎng)時(shí)間序列和較強(qiáng)普適性已被大量用于相關(guān)方面的研究[62]， 因此本文采用遙感數(shù)據(jù)產(chǎn)品Ocean Productivity(www.science.oregonstate.edu/ocean.productivity/)得出。

其他指標(biāo)的數(shù)據(jù)來源于《中國(guó)統(tǒng)計(jì)年鑒》、《中國(guó)海洋經(jīng)濟(jì)統(tǒng)計(jì)公報(bào)》、《中國(guó)海洋年鑒》、《中國(guó)海洋統(tǒng)計(jì)年鑒》、《中國(guó)海洋環(huán)境公報(bào)》、《中國(guó)海洋災(zāi)害公報(bào)》、《中國(guó)海洋生態(tài)環(huán)境狀況公報(bào)》、《中國(guó)漁業(yè)年鑒》、《中國(guó)漁業(yè)統(tǒng)計(jì)年鑒》、《中國(guó)近岸海域環(huán)境質(zhì)量公報(bào)》、《上海統(tǒng)計(jì)年鑒》、《江蘇統(tǒng)計(jì)年鑒》、《浙江統(tǒng)計(jì)年鑒》、《上海市海洋環(huán)境質(zhì)量公報(bào)》、《浙江省海洋環(huán)境公報(bào)》、《浙江省漁業(yè)水域環(huán)境質(zhì)量公報(bào)》、《江蘇省近岸海域水環(huán)境質(zhì)量公報(bào)》及各沿海城市的統(tǒng)計(jì)年鑒等公開信息， 樣本區(qū)間為2009—2018年， 采樣頻率為次/a。同時(shí)， 對(duì)于海洋捕撈量， 基于研究范圍的界定， 剔除了遠(yuǎn)洋捕撈的數(shù)值。

2 結(jié)果和分析

2.1 海洋生態(tài)系統(tǒng)服務(wù)價(jià)值

2.1.1 各城市海洋生態(tài)系統(tǒng)服務(wù)的價(jià)值均呈增加狀態(tài)

由圖1可知， 研究期內(nèi)， 長(zhǎng)三角各沿海城市的海洋生態(tài)系統(tǒng)服務(wù)價(jià)值均呈上升態(tài)勢(shì)， 但不是一直增加， 均存在某一年份低于之前年份的個(gè)例。

圖1 各城市海洋生態(tài)系統(tǒng)服務(wù)價(jià)值(未消除貨幣通貨膨脹影響)Fig. 1 Marine ecosystem service value of each city (with impact of inflation)

進(jìn)一步地， 運(yùn)用指數(shù)調(diào)整法， 以2009年為基期，消除貨幣通貨膨脹的影響后3考慮海洋生態(tài)系統(tǒng)服務(wù)價(jià)值的計(jì)算方法， 氣候調(diào)節(jié)、氣體調(diào)節(jié)、廢棄物處理、科學(xué)研究等服務(wù)不用進(jìn)行通貨膨脹影響的消除。考慮2009—2018年人民幣對(duì)美元的匯率變化幅度不大， 參考黎鶴仙等[43]的計(jì)算方法， 可認(rèn)為基因資源供給、干擾調(diào)節(jié)、生物控制、精神文化、物種多樣性維持、營(yíng)養(yǎng)物質(zhì)循環(huán)等服務(wù)的計(jì)算方式與前述四種類別的計(jì)算方式相同， 因此， 未對(duì)其進(jìn)行通貨膨脹影響的消除。綜上所述， 本文只考慮了通貨膨脹對(duì)食品供給、原材料供給、旅游娛樂等服務(wù)的影響。， 得到的最終結(jié)果(圖2)也呈現(xiàn)出基本一致的態(tài)勢(shì)。這充分反映出在研究期內(nèi)各城市整體上對(duì)海洋資源環(huán)境的利用程度增加， 說明海洋資源環(huán)境對(duì)長(zhǎng)三角區(qū)域的重要性顯著增強(qiáng)。

圖2 各城市海洋生態(tài)系統(tǒng)服務(wù)價(jià)值(消除貨幣通貨膨脹影響)Fig. 2 Marine ecosystem service value of each city (without impact of inflation)

2.1.2 各城市相同類別海洋生態(tài)系統(tǒng)服務(wù)的價(jià)值變化情況相近

基于圖1的結(jié)果， 繪制各城市的海洋生態(tài)系統(tǒng)服務(wù)價(jià)值組成變化態(tài)勢(shì)圖(圖3)。據(jù)此， 分類別看， 各城市海洋生態(tài)系統(tǒng)的調(diào)節(jié)服務(wù)價(jià)值均呈現(xiàn)明顯的波動(dòng)，體現(xiàn)長(zhǎng)三角區(qū)域海域初級(jí)生產(chǎn)力不穩(wěn)定， 進(jìn)而反映區(qū)域整體海洋生態(tài)環(huán)境不平穩(wěn)。其他類別的價(jià)值均增加，特別是文化服務(wù)價(jià)值一直為上升態(tài)勢(shì)， 體現(xiàn)出各城市對(duì)海洋生態(tài)系統(tǒng)文化服務(wù)的重視， 特別大力發(fā)展旅游娛樂的成效及游客對(duì)長(zhǎng)三角區(qū)域的認(rèn)可與向往。

2.1.3 各城市海洋生態(tài)系統(tǒng)服務(wù)價(jià)值的比值相對(duì)穩(wěn)定

由圖1可知， 上海市、舟山市的海洋生態(tài)系統(tǒng)服務(wù)價(jià)值始終明顯高于其他城市， 在區(qū)域內(nèi)占有較大比例。連云港市、嘉興市始終位列后幾名， 這與其海域面積較小(表3)、食品供給服務(wù)價(jià)值的數(shù)值較低(圖3)等有關(guān)。但嘉興市的旅游娛樂的增長(zhǎng)幅度較大(圖3、圖4)， 使其總值增長(zhǎng)。而連云港市的數(shù)值最終增加幅度較小， 比值降低。其他城市的排名變化不明顯。

圖4 嘉興市海洋生態(tài)系統(tǒng)文化服務(wù)價(jià)值Fig. 4 Marine ecosystem service value of culture in Jiaxing

表3 各城市海域面積(單位: 103 km2)Tab. 3 Sea area of each city (unit: 103 km2)

2.1.4 各城市相同類別海洋生態(tài)系統(tǒng)服務(wù)的價(jià)值占比情況不同

由圖3可知， 除上海市、嘉興市、舟山市、臺(tái)州市外， 其他城市中調(diào)節(jié)服務(wù)的價(jià)值在總價(jià)值中占比一直為最大。對(duì)于上海市， 文化服務(wù)所占比例一直是最大， 其次才是調(diào)節(jié)服務(wù)。上海市的文化服務(wù)中， 旅游娛樂貢獻(xiàn)了主要部分(圖5)。嘉興市、舟山市的調(diào)節(jié)服務(wù)價(jià)值開始占有最大比例， 最后被文化服務(wù)價(jià)值反超。與這兩個(gè)城市不同， 臺(tái)州市的變化主要是由供給服務(wù)價(jià)值引起的。這說明對(duì)不同城市， 海洋資源環(huán)境發(fā)揮的具體作用不同， 同時(shí)還可能變化。因此各城市應(yīng)在注重海洋資源環(huán)境這一基本性原則的基礎(chǔ)上， 結(jié)合其發(fā)揮的不同作用及城市未來發(fā)展的政策與方向有所選擇地進(jìn)行重點(diǎn)關(guān)注。

圖3 各城市海洋生態(tài)系統(tǒng)服務(wù)價(jià)值組成Fig. 3 Composition of the marine ecosystem service value of each city

圖5 上海市海洋生態(tài)系統(tǒng)文化服務(wù)價(jià)值Fig. 5 Marine ecosystem service value of culture in Shanghai

2.2 海洋生態(tài)系統(tǒng)服務(wù)價(jià)值影響因素

2.2.1 靜態(tài)面板回歸結(jié)果

為了分析表2中各指標(biāo)對(duì)長(zhǎng)三角區(qū)域海洋生態(tài)系統(tǒng)服務(wù)價(jià)值的影響， 本文采用最小二乘法(OLS)和面板模型等計(jì)量方法進(jìn)行定量評(píng)估。

表4 為OLS和面板模型的回歸結(jié)果， 其中一共包含5種模型。模型1調(diào)整之后的擬合優(yōu)度為0.685 000， 說明模型的整體擬合效果較好。WD在1%的水平上高度顯著， 說明廢水排放量是影響海洋生態(tài)系統(tǒng)服務(wù)價(jià)值的一個(gè)關(guān)鍵因素。IE、WWDDO的系數(shù)值分別為–0.008 000、0.023 800， 但在統(tǒng)計(jì)意義上不顯著， 說明工業(yè)廢氣排放量、廢水直排入海量并不是影響海洋生態(tài)系統(tǒng)服務(wù)價(jià)值的關(guān)鍵因素。SWD的系數(shù)值為–0.159 000， 在1%的水平上高度顯著， 說明總體而言， 固體廢物傾倒丟棄量每增加一個(gè)單位， 會(huì)使得海洋生態(tài)系統(tǒng)服務(wù)價(jià)值減少0.159 000個(gè)單位。MY、RA的系數(shù)值分別為–0.001 390、–0.077 700， 均在1%的水平上高度顯著， 說明海水養(yǎng)殖產(chǎn)量、圍填海面積的提高也會(huì)顯著降低海洋生態(tài)系統(tǒng)服務(wù)價(jià)值。MA、CO的系數(shù)值分別為0.003 470、0.000 722， 也均在1%的水平上高度顯著，說明海水養(yǎng)殖面積、海洋捕撈量的增加會(huì)顯著增加海洋生態(tài)系統(tǒng)服務(wù)價(jià)值。

表4 靜態(tài)面板回歸結(jié)果Tab. 4 Static panel regression results

續(xù)表

但是， 由于本文的數(shù)據(jù)是面板數(shù)據(jù)格式， 只采用OLS估計(jì)可能會(huì)導(dǎo)致遺漏變量問題， 因此模型2和模型3分別為控制個(gè)體因素的固定效應(yīng)模型和隨機(jī)效應(yīng)模型的回歸結(jié)果。依據(jù)調(diào)整之后的R2值， 模型2為0.462 000， 要低于模型3的0.713 000， 說明隨機(jī)效應(yīng)模型的擬合優(yōu)度要顯著高于固定效應(yīng)模型。然而， 模型2和模型3之間的豪斯曼檢驗(yàn)為118.850 0， 且在1%水平上高度顯著， 說明固定效應(yīng)模型要比隨機(jī)效應(yīng)模型更好。綜合對(duì)比之后， 本文采用固定效應(yīng)模型即模型2來進(jìn)行分析。WD的系數(shù)值為–0.003 600， 但只是在5%的水平上顯著， 說明廢水排放量雖然仍是影響海洋生態(tài)系統(tǒng)服務(wù)價(jià)值的一個(gè)因素， 但其重要性有所降低。IE的系數(shù)符號(hào)雖然與OLS模型的結(jié)果不一致， 同時(shí)數(shù)值也有一定出入， 但同樣為不顯著因素。WWDDO的系數(shù)值為0.030 100， 雖然數(shù)值等級(jí)與符號(hào)同OLS模型的結(jié)果一致， 但其為5%的水平上顯著， 說明廢水直排入海量是影響海洋生態(tài)系統(tǒng)服務(wù)價(jià)值的一個(gè)較重要因素。與此相反， SWD變?yōu)榱瞬伙@著因素， 這與各城市的固體廢物傾倒丟棄量從2011、2012年左右就減為0相一致。MY的系數(shù)雖然數(shù)值等級(jí)同OLS模型的結(jié)果一致， 但其符號(hào)變?yōu)檎颍?同時(shí)顯著性水平也降為5%。MA的顯著性雖然與OLS模型的結(jié)果保持了高度的統(tǒng)一性， 但其系數(shù)符號(hào)變?yōu)樨?fù)值。CO的系數(shù)雖然數(shù)值等級(jí)與符號(hào)同OLS模型的結(jié)果一致， 但是顯著性水平降為5%。RA的系數(shù)符號(hào)、顯著性水平雖然沒有變化， 但其數(shù)值有一定出入。

另一方面， 考慮到存在一些隨著時(shí)間變化而不對(duì)個(gè)體變化的遺漏變量因素， 我們?cè)谀Ｐ?和模型3的基礎(chǔ)之上進(jìn)一步考慮年份固定效應(yīng)， 即采用雙向固定模型來檢驗(yàn)。模型4和模型5分別為相應(yīng)的估計(jì)結(jié)果。同理， 首先依據(jù)調(diào)整之后的R2值， 模型4為0.578 000， 模型5為0.834 000， 但模型4和模型5之間的豪斯曼檢驗(yàn)為57.560 000， 且在1%水平上高度顯著， 說明固定效應(yīng)模型要比隨機(jī)效應(yīng)模型更好，這和模型2、模型3得出的結(jié)論是一致的。但是模型4的回歸結(jié)果顯示， 只有廢水直排入海量與圍填海面積為顯著相關(guān)因素， 同時(shí)大部分自變量的系數(shù)絕對(duì)值要低于模型2的結(jié)果， 此外前者的調(diào)整R2值只是略高于后者， 因此本文認(rèn)為模型2的擬合優(yōu)度更好。

綜上所述， 模型2是經(jīng)過對(duì)比之后遴選出的最適宜模型， 從它的回歸結(jié)果來看， 廢水排放量、廢水直排入海量、海水養(yǎng)殖產(chǎn)量、海洋捕撈量為對(duì)海洋生態(tài)系統(tǒng)服務(wù)價(jià)值產(chǎn)生影響的較關(guān)鍵指標(biāo)， 海水養(yǎng)殖面積、圍填海面積為對(duì)海洋生態(tài)系統(tǒng)服務(wù)價(jià)值產(chǎn)生影響的關(guān)鍵指標(biāo)。其中， 廢水排放量、海水養(yǎng)殖面積、圍填海面積對(duì)海洋生態(tài)系統(tǒng)服務(wù)價(jià)值起到抑制作用。

2.2.2 動(dòng)態(tài)面板回歸結(jié)果

理論上， 之前的海洋生態(tài)系統(tǒng)服務(wù)情況會(huì)影響當(dāng)期的海洋生態(tài)系統(tǒng)服務(wù)水平。因此， 為了深度檢驗(yàn)各因素與海洋生態(tài)系統(tǒng)服務(wù)價(jià)值之間的影響關(guān)系，本文加入了滯后項(xiàng)作為工具變量， 并選擇GMM動(dòng)態(tài)面板模型。具體而言， 本文使用海洋生態(tài)系統(tǒng)服務(wù)價(jià)值OESV的滯后一期值作為工具變量， 同時(shí)選擇差分GMM(DiffGMM)和系統(tǒng)GMM(SyseGMM)作為模型表達(dá)。表5表明， 四個(gè)模型的有效性檢驗(yàn)結(jié)果均顯示工具變量不存在過度識(shí)別問題， 同時(shí)AR檢驗(yàn)結(jié)果均顯示擾動(dòng)項(xiàng)無自相關(guān)問題， 模型效果較好。

表5 動(dòng)態(tài)面板回歸結(jié)果Tab. 5 Dynamic panel regression results

從海洋生態(tài)系統(tǒng)服務(wù)價(jià)值的視角看， 滯后一期即為正的高顯著水平， 表示上一期的海洋生態(tài)系統(tǒng)服務(wù)價(jià)值對(duì)本期具有明顯的促進(jìn)作用。該結(jié)果意味著上一期價(jià)值的情況較令人滿意， 決定了在本期相關(guān)行為將持續(xù)發(fā)生， 以維持類似的狀態(tài)。這一情況也促使了長(zhǎng)三角各城市的海洋生態(tài)系統(tǒng)服務(wù)價(jià)值均呈現(xiàn)出上升的態(tài)勢(shì)。

另一方面， 由于模型3、4是在模型1、2基礎(chǔ)之上考慮了年份固定效應(yīng)， 因此模型3、4更加合理。此外， 模型3、4之間還互為穩(wěn)健性檢驗(yàn)。

綜合靜態(tài)分析與動(dòng)態(tài)分析的結(jié)果， 可以判斷廢水排放量、廢水直排入海量、海水養(yǎng)殖產(chǎn)量、海水養(yǎng)殖面積、海洋捕撈量、圍填海面積都是能夠?qū)Ｑ笊鷳B(tài)系統(tǒng)服務(wù)價(jià)值產(chǎn)生影響的顯著指標(biāo)。其中， 廢水排放量、海水養(yǎng)殖面積、圍填海面積對(duì)海洋生態(tài)系統(tǒng)服務(wù)價(jià)值起到抑制作用。同時(shí)， 工業(yè)廢氣排放量也有一定影響。究其原因， 廢水排放量主要反映的是陸地上的廢水排放情況， 而這些廢水通過江、河、湖等最終匯入海洋， 對(duì)海水的影響是大范圍且普遍的，所以其通過對(duì)供給服務(wù)、文化服務(wù)等的負(fù)向影響最終導(dǎo)致對(duì)海洋生態(tài)系統(tǒng)服務(wù)價(jià)值產(chǎn)生抑制作用， 但由于各類別價(jià)值比例的關(guān)系， 其系數(shù)的絕對(duì)值不是很高。考慮調(diào)節(jié)服務(wù)(特別是廢棄物處理、氣候調(diào)節(jié))價(jià)值的計(jì)算方法(表1)與調(diào)節(jié)服務(wù)所占的比例情況(圖3)， 廢水直排入海量必然會(huì)對(duì)海洋生態(tài)系統(tǒng)服務(wù)價(jià)值產(chǎn)生較大的正向影響。同理， 也可說明工業(yè)廢氣排放量的影響情況。由于海水養(yǎng)殖產(chǎn)量與海洋捕撈量對(duì)供給服務(wù)的數(shù)值變化造成明顯的促進(jìn)， 因此這兩個(gè)指標(biāo)對(duì)海洋生態(tài)系統(tǒng)服務(wù)價(jià)值呈現(xiàn)的是正向影響， 但由于供給服務(wù)所占的比例總體上不高， 所以其系數(shù)不大。圍填?；顒?dòng)會(huì)顯著減少海域面積， 因此其對(duì)海洋生態(tài)系統(tǒng)服務(wù)價(jià)值會(huì)呈現(xiàn)較顯著的負(fù)向影響。而對(duì)于海水養(yǎng)殖面積這一指標(biāo)， 可能由于區(qū)域內(nèi)海水養(yǎng)殖占用了大面積的海域， 限制了廢水等的直排入?；顒?dòng)， 同時(shí)養(yǎng)殖面積與養(yǎng)殖產(chǎn)量間不存在必然的正向聯(lián)系， 再結(jié)合調(diào)節(jié)服務(wù)的占比情況(圖3)，因此其對(duì)海洋生態(tài)系統(tǒng)服務(wù)價(jià)值產(chǎn)生了較大的負(fù)向影響。

另一方面， 考慮供給服務(wù)、調(diào)節(jié)服務(wù)價(jià)值的計(jì)算方法帶來的問題4考慮供給服務(wù)(特別是食品供給)、調(diào)節(jié)服務(wù)(特別是廢棄物處理、氣候調(diào)節(jié))價(jià)值的計(jì)算方法， 海洋捕撈量、廢水直排入海量、工業(yè)廢氣排放量等指標(biāo)的正向結(jié)果實(shí)際上還包含了其對(duì)海洋資源環(huán)境的不利影響。， 總體結(jié)果也反映出長(zhǎng)三角地區(qū)海洋資源環(huán)境的整體狀況仍需改進(jìn)。考慮指標(biāo)的種類，這與韓婕妤[63]的研究結(jié)論基本一致。

3 結(jié)論和政策建議

本文計(jì)算了2009—2018年長(zhǎng)三角各沿海城市海洋生態(tài)系統(tǒng)服務(wù)價(jià)值， 運(yùn)用面板模型， 分析了區(qū)域海洋生態(tài)系統(tǒng)服務(wù)價(jià)值的影響因素， 結(jié)果發(fā)現(xiàn)： 1) 各城市海洋生態(tài)系統(tǒng)服務(wù)價(jià)值均有增加， 同時(shí)相同類別的海洋生態(tài)系統(tǒng)服務(wù)的價(jià)值變化情況相近； 2) 各城市海洋生態(tài)系統(tǒng)服務(wù)價(jià)值的比值相對(duì)穩(wěn)定， 但不同城市相同類別海洋生態(tài)系統(tǒng)服務(wù)的價(jià)值占比不同； 3) 長(zhǎng)三角地區(qū)海洋資源環(huán)境的整體狀況仍需改進(jìn)； 4) 廢水直排入海、工業(yè)廢氣排放、海水養(yǎng)殖、圍填海等活動(dòng)對(duì)區(qū)域海洋生態(tài)系統(tǒng)服務(wù)價(jià)值有較大影響， 特別是廢水直排入海、海水養(yǎng)殖、圍填海等活動(dòng)影響顯著； 5) 不同人類活動(dòng)對(duì)區(qū)域海洋資源環(huán)境的影響機(jī)制不同， 同時(shí)作用方向也不相同。

基于上述結(jié)果， 本文提出以下政策建議：

第一， 依靠科技進(jìn)步提升海洋生態(tài)系統(tǒng)服務(wù)價(jià)值水平， 同時(shí)注重海洋生態(tài)系統(tǒng)文化服務(wù)的作用。對(duì)于各城市而言， 海域面積制約了海洋生態(tài)系統(tǒng)服務(wù)價(jià)值的進(jìn)一步增長(zhǎng)。因此， 要大幅度增加海洋生態(tài)系統(tǒng)服務(wù)價(jià)值， 只能依賴于技術(shù)進(jìn)步等非自然的要素進(jìn)行效率的提升與新服務(wù)種類的認(rèn)知和開發(fā)[28]。對(duì)此， 各沿海城市應(yīng)完善相關(guān)政策， 進(jìn)一步加大投入，并通過組建產(chǎn)學(xué)研用聯(lián)盟等方式， 進(jìn)行海洋資源環(huán)境開發(fā)利用方面的相關(guān)技術(shù)特別是綠色技術(shù)的突破，提升海洋科技水平。同時(shí)， 完善成果交易平臺(tái)、產(chǎn)業(yè)孵化器等成果轉(zhuǎn)移轉(zhuǎn)化體系建設(shè)， 并大力實(shí)施成果的推廣應(yīng)用等工作， 加強(qiáng)現(xiàn)有海洋科技成果的轉(zhuǎn)化。此外， 要推動(dòng)海洋科技人才發(fā)展體制機(jī)制創(chuàng)新。首先，要繼續(xù)重視海洋科技人才的培育和引進(jìn)工作。其次，完善并創(chuàng)新海洋科技人才的激勵(lì)機(jī)制， 使其在海洋科技創(chuàng)新、推動(dòng)海洋產(chǎn)業(yè)發(fā)展等方面發(fā)揮更加顯著的作用。

另一方面， 海洋生態(tài)系統(tǒng)文化服務(wù)的作用發(fā)揮受海域面積大小的影響較輕， 因此應(yīng)進(jìn)一步注重該服務(wù)， 特別是其中的旅游娛樂的發(fā)展。這對(duì)于連云港市、嘉興市等海域面積較小、海洋生態(tài)系統(tǒng)服務(wù)價(jià)值在區(qū)域中所占比例不高的城市尤為重要。只有這樣， 才可以增強(qiáng)在長(zhǎng)三角海洋一體化發(fā)展中的地位，更好地參與區(qū)域共同發(fā)展。

第二， 精準(zhǔn)做好區(qū)域海洋資源環(huán)境的整體改善工作。廢水排放、廢水直排入海、工業(yè)廢氣排放、海洋捕撈、圍填海等是影響長(zhǎng)三角區(qū)域海洋資源環(huán)境狀況的主要因素， 特別是廢水直排入海與圍填海的影響尤其顯著。因此， 區(qū)域?qū)Ｑ筚Y源環(huán)境的共保共治等一體化治理應(yīng)著重關(guān)注這些領(lǐng)域， 并加強(qiáng)海洋生態(tài)環(huán)境的穩(wěn)定性。同時(shí)， 考慮海洋的流動(dòng)性、開發(fā)性等特點(diǎn)， 在現(xiàn)行體制下， 海洋資源環(huán)境問題的有效解決往往需要跨部門、跨區(qū)域共同合作， 因此各城市內(nèi)部要構(gòu)建良好的涉海部門溝通協(xié)調(diào)機(jī)制(如海洋綜合管理聯(lián)席會(huì)制度)， 以形成一個(gè)自上而下、長(zhǎng)期有效的制度安排； 此外， 還要注重陸海統(tǒng)籌， 加強(qiáng)陸、海不同管理部門間的溝通協(xié)作； 另一方面， 要在充分考慮各城市正當(dāng)利益訴求的前提下， 建立完善的橫向溝通協(xié)調(diào)與合作機(jī)制并深化實(shí)施， 以實(shí)現(xiàn)區(qū)域內(nèi)部的良好協(xié)同效應(yīng)。

第三， 實(shí)施差異化治理措施。已有研究表明， 一體化的實(shí)現(xiàn)需要考慮內(nèi)部各區(qū)域的差異性[64-65]， 否則簡(jiǎn)單的同一政策措施會(huì)導(dǎo)致超載、同質(zhì)化和自然資源環(huán)境的損失等問題的出現(xiàn)[66]， 從而使得一體化的效果大打折扣。因此， 不同城市， 在區(qū)域海洋資源環(huán)境整體保護(hù)、改善的基本前提下， 應(yīng)結(jié)合城市未來發(fā)展、產(chǎn)業(yè)結(jié)構(gòu)、海洋生態(tài)系統(tǒng)服務(wù)價(jià)值的重要類別判定等， 進(jìn)一步選擇重點(diǎn)治理領(lǐng)域與對(duì)象， 同時(shí)根據(jù)各地的治理水平、能力等具體情況， 采用合適的手段， 以更好地達(dá)到區(qū)域海洋一體化目標(biāo)。

本文存在一定的研究局限性： 1) 資料所限， 本文只采用10年的數(shù)據(jù)分析， 時(shí)間相對(duì)較短； 2) 本文對(duì)海洋生態(tài)系統(tǒng)服務(wù)價(jià)值計(jì)算的模型公式較為簡(jiǎn)單，考慮的因素較少(如對(duì)于非使用價(jià)值基本沒有考慮)，同時(shí)考慮同一種生態(tài)系統(tǒng)服務(wù)在不同的社會(huì)和經(jīng)濟(jì)環(huán)境中會(huì)表現(xiàn)出不同的價(jià)值， 所以最終計(jì)算出的結(jié)果不能反映其全部真實(shí)價(jià)值水平； 3) 對(duì)海洋生態(tài)系統(tǒng)服務(wù)價(jià)值的影響分析只考慮了人類行為對(duì)各城市自身的影響， 沒有考慮對(duì)其他城市(特別是相鄰城市)的交互影響， 結(jié)果難以反映全部情況。這些都是以后研究中需要進(jìn)一步解決的問題。