基于能值的山東省海洋生態(tài)足跡研究

狄乾斌， 張海紅，曹可

（1.遼寧師范大學 海洋經濟與可持續(xù)發(fā)展研究中心，遼寧 大連 116029；2.國家海洋環(huán)境監(jiān)測中心，遼寧 大連 116026）

在可持續(xù)發(fā)展研究中，海洋經濟可持續(xù)發(fā)展的定量分析（韓增林等，2003；狄乾斌等，2009；邵桂蘭 等，2011；李夫星等，2013；方欣 等，2013） 成為不少學者的研究方向。生態(tài)足跡模型的提出（Rees，1992），為定量研究海洋經濟的可持續(xù)發(fā)展提供了新的理論和方法。

不少學者運用傳統(tǒng)的生態(tài)足跡模型（徐中民等，2003；魏明等，2006；楊天姿 等，2009） 或者修正的生態(tài)足跡模型（隋春花等，2006；楊海波等，2009；楊海真等，2010） 定量測度了不同層次區(qū)域的可持續(xù)發(fā)展狀況。20 世紀90年代初引入能值模型（藍盛芳等，2001，2002），國內學者（王建源等，2007；趙志強等，2008，2009；劉淼等，2008） 將傳統(tǒng)生態(tài)足跡模型與能值理論相結合，運用改進的能值-生態(tài)足跡模型（Emergybased Ecological Footprint，EEF），研究了不同地區(qū)的能值-生態(tài)承載力和能值-生態(tài)足跡。在能值-生態(tài)承載力（Emergy-based carrying capacity， ECC）的計算上，大多采用全球平均能值密度（張芳怡等，2006；隋春鋒等，2008），而能值-生態(tài)足跡則采用區(qū)域能值密度，但也有不少學者從不同角度改進了能值足跡法：金丹等（2010） 以地球生物圈循環(huán)中3 種主要能源在特定區(qū)域中的能值總量計算出區(qū)域能值密度，以同一能值密度計算區(qū)域能值生態(tài)承載力和生態(tài)足跡；趙志強等（2008） 將輸入系統(tǒng)的可更新資源定義為潛在生態(tài)承載力，將研究區(qū)產品產出、進出口以及勞務均納入生態(tài)承載力范圍；劉淼等（2008） 將研究區(qū)內所有人類可利用的可再生資源的能值作為生態(tài)承載力總能值，該能值與研究區(qū)域面積之比作為區(qū)域能值密度，并運用此能值密度計算生態(tài)承載力。這些研究從不同角度對區(qū)域經濟進行分析，為區(qū)域可持續(xù)發(fā)展提供決策依據(jù)。國內學者（楊山 等2011；趙美霞等，2011）分別將生態(tài)足跡模型運用到海洋可持續(xù)研究中，然而運用能值-生態(tài)足跡法在區(qū)域海洋經濟方面的研究還不多見。能值-生態(tài)足跡模型的提出，克服了生態(tài)足跡法偏生態(tài)性的特點（顧康康等，2009），從生態(tài)環(huán)境、社會發(fā)展和經濟發(fā)展角度進行評價，更能真實的反映區(qū)域的可持續(xù)發(fā)展狀況。本文基于能值理論與生態(tài)足跡模型在可持續(xù)發(fā)展評價上的可行性，將能值-生態(tài)足跡模型運用到海洋生態(tài)系統(tǒng)中，構建了海洋能值-生態(tài)足跡創(chuàng)新模型，并以山東海洋生態(tài)系統(tǒng)為例，定量分析了特定區(qū)域海洋經濟的可持續(xù)發(fā)展。

1 研究區(qū)域

山東省海域面積1.595×105km2，分布于渤海和黃海，以蓬萊角為界，向西屬于渤海海域，向東屬于黃海海域。海岸線全長3 345 km，北起魯冀交界處的漳衛(wèi)新河河口（381 543.77 N），南至魯蘇交界處的繡針河河口（350 445 N）。山東沿海屬暖溫帶季風氣候，冬季盛行偏北風，夏季盛行偏南風，海域年均風速為7.1～7.5 m/s，年均降水量在580 916 mm之間，日照時數(shù)23 002 900 hr。沿海潮汐主要為不正規(guī)半日潮和正規(guī)半日潮，海上風能、潮汐能、波浪能等可開發(fā)海洋能源資源儲量豐富。

山東省是我國海洋經濟大省，其海洋生產總值一直居全國前列。2010年實現(xiàn)海洋生產總值7074.5 億元，占全國比例為17.9 %，三次海洋產業(yè)比重分別為6.3%、50.2%和43.5%。海洋漁業(yè)產量631 萬t（遠洋捕撈1.5×105t），海鹽產量2.273×107t，均居全國第一。海洋漁業(yè)、海洋化工業(yè)、海洋生物醫(yī)藥業(yè)、海洋電力業(yè)、濱海旅游業(yè)和海洋交通運輸業(yè)分別比上年增長16.4%、36.0%、68.9%、35.6 %、22.7 %和8.0 %。海洋資源儲藏量豐富，海洋經濟持續(xù)增長。

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

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

本文中可更新資源基礎數(shù)據(jù)來源于《山東省海洋功能區(qū)劃（2011-2020）》，海洋資源消費項目數(shù)據(jù)來源于《中國海洋統(tǒng)計年鑒2001-2011》、《中國海洋年鑒2011》、 《山東統(tǒng)計年鑒2011》 和《2010年山東省國民經濟和社會發(fā)展統(tǒng)計公報》。數(shù)據(jù)計算與制圖分別采用EXCEL 和WORD 完成。

2.2 海洋能值-生態(tài)足跡模型

2.2.1 海洋能值-生態(tài)足跡的概念

Folke 等（1991；1998） 將海洋生態(tài)足跡定義為人類為滿足自身需要的海產品以及服務所占用的海洋生態(tài)系統(tǒng)的表面積。Warren-Rhodes 等（2003） 給出了海洋生態(tài)足跡的數(shù)值表達式：MEF=SC/AS，SC 表示生產一定消耗量所需要的海洋面積，AS 表示某區(qū)域的實際海洋面積，以MEF 值的大小來判斷研究區(qū)域海洋生態(tài)的赤字或者盈余。國內學者陳成忠等（2008） 提出海洋足跡和海洋生物承載力兩個概念：海洋足跡是指人類單位水域海洋生物資源的捕撈量，海洋生物承載力是指海洋單位水域供給人類生物資源的能力，并運用非線性科學理論模擬了最大可持續(xù)海洋足跡。

依據(jù)上述海洋生態(tài)足跡的定義，結合由H.T.Odum 創(chuàng)立的能值理論，基于生態(tài)系統(tǒng)中的能量直接或間接源于太陽能，通過太陽能值轉換率將生態(tài)系統(tǒng)或生態(tài)經濟系統(tǒng)中不同類型、不能直接比較的能量換算成可進行比較的能值標準——太陽能值。本文將海洋能值-生態(tài)足跡定義為：人類為滿足自身需要的海產品以及服務所消耗的海洋生態(tài)系統(tǒng)的能量所占用的（通過能值轉換率和海域能值密度換算） 的海域表面積。海洋能值-生態(tài)承載力是海洋生態(tài)系統(tǒng)中可更新資源的能值。

表1 2000-2010年山東海域自然資源原始數(shù)據(jù)

海洋水體的流動性、海洋環(huán)境的復雜性以及海洋開發(fā)的立體性導致海洋經濟的研究比陸域經濟的研究要復雜很多。本文將能值-生態(tài)足跡模型應用到海洋經濟的研究中，主要基于以下幾點考慮：（1） 每個沿海地區(qū)都具有不同的區(qū)域特性，其海洋自然資源、海陸經濟發(fā)展水平和社會條件各不相同。從區(qū)域經濟的科學研究方面，可以將具有不同發(fā)展條件的海洋區(qū)域劃分為一個個相對獨立的經濟研究單位。（2） 基于海域自然能源資源由特定海域區(qū)位、氣候等條件決定，海域礦產資源是由海域生物地化循環(huán)長期累積的結果，因此短期內均不會發(fā)生較大變化。海洋生物資源以及海域生產力在短時間內容易受自然災害以及人為災害等的影響，但是基于海洋生態(tài)系統(tǒng)強大的自我調節(jié)功能，在一定的時間周期內，也不會發(fā)生較大變化。（3） 以研究地區(qū)所轄海域范圍海洋資源的產量和海洋能源的利用量作為生態(tài)足跡指標，范圍之外的資源消耗則為外界輸入能值。圖1 是特定海域海洋生態(tài)系統(tǒng)的自然資源的能值流量圖。

2.2.2 海洋生態(tài)系統(tǒng)自然資源的能值計算

海洋生態(tài)系統(tǒng)中可更新自然資源包括太陽輻射能、風能、雨水化學能、潮汐能、潮流能、波浪能、溫差能、鹽差能和地球旋轉能，其總能值用Er 表示。不可更新能源資源包括海洋石油、海洋天然氣，其總能值用En 表示。某自然資源的太陽能值（sej）=其有效能（J）×太陽能值轉換率（sejJ）。

圖1 特定海域海洋自然資源能值流量圖

2.2.3 海域能值密度的計算

因太陽輻射能、海洋風能、雨水化學能、波浪能、溫差能和鹽差能均是由太陽輻射作用產生的，潮汐能和潮流能是由太陽和月球對地球的萬有引力引起的，計算海域能值密度時只取其中最大值，則海域能值密度計算公式為：

式中：D 表示海域能值密度；Er[maxA,maxB，C]表示太陽輻射能、海洋風能、雨水化學能、波浪能、溫差能、鹽差能中的能值最大值與潮汐能和潮流能的能值最大值以及地球旋轉能之和；S 表示研究區(qū)海域面積。

2.2.4 海域能值/貨幣比率（sej/$） 的計算

因人類消耗的自然環(huán)境資源可用能值來衡量，經濟系統(tǒng)中循環(huán)流通的貨幣的實際購買力可以通過單位貨幣所能購買的能值量（能值/貨幣比率） 來衡量。因此，可以通過能值/貨幣比率將以貨幣形式表示的資源環(huán)境消費項目轉換為能值，進而計算其能值-生態(tài)足跡。根據(jù)藍盛芳（2002），由國家的能值/貨幣比率推出海域能值/貨幣比率的計算公式：

式（2） 中：該國全年能值投入總量包括可更新資源、不可更新資源及進口商品、資源的能值。式（3） 中，表示研究區(qū)海域可更新自然資源能值投入總量；表示不可更新自然資源的能值投入總量；表示輸入的外界資源能值。海域能值/貨幣比率則表示某海域單位海洋生產總值所對應的能值量。

2.2.5 海洋能值-生態(tài)承載力的計算

海洋能值-生態(tài)承載力（Marine Emergy Ecological Capacity） 的計算公式為：

式中：MEEC 表示海洋能值-生態(tài)承載力；Er 表示可利用的海洋生態(tài)系統(tǒng)可更新自然資源的能值；D表示海域能值密度。

2.2.6 海洋能值-生態(tài)足跡的計算

海洋能值-生態(tài)足跡（Marine Emergy Ecological Footprint） 的計算公式為：

式中，MEEC 表示海洋能值-生態(tài)足跡；Ec 表示消耗的海洋資源的總能值，海洋生態(tài)足跡消費項目包括海洋漁業(yè)（海洋捕撈和海水養(yǎng)殖）、海鹽、海洋原油、海洋天然氣、海洋化工、海洋生物醫(yī)藥、海洋電力、海水綜合利用、海洋工程建筑以及排放到海洋中的廢水、廢固。由于海洋化工是包括海鹽化工、海水化工、海洋石油化工等的生產活動，因此，為了避免重復計算，只取海洋資源的原始產量數(shù)據(jù)；海洋生物醫(yī)藥是由海水產品制成的，計算生態(tài)足跡時只取海洋捕撈和海水養(yǎng)殖產量；D 表示海域能值密度。

3 山東海洋生態(tài)系統(tǒng)能值-生態(tài)足跡分析

由于山東海域海洋自然資源的產出不只是供應山東省本地的資源消費，還包括國內以及國外一部分地區(qū)，但是外界消耗的山東海洋資源量沒有比較具體的統(tǒng)計數(shù)據(jù)，因此，本文通過計算海域總能值-生態(tài)承載力和總能值-生態(tài)足跡來進行海洋生態(tài)系統(tǒng)評價。

表2 是山東海域可更新資源能值?；?.26×1024sej/yr 的全球能值基準（Campbell et al，2004；Campbell，2000），根據(jù)公式（4） 計算得出山東海洋生態(tài)系統(tǒng)總能值-生態(tài)承載力為1.60 × 107hm2。根據(jù)公式（1） 計算得出山東海洋生態(tài)系統(tǒng)能值密度為1.24 ×1017sej/hm2。海域能值-生態(tài)承載力的計算只考慮了海洋生態(tài)系統(tǒng)中5 種可更新能源資源的總能值，海流能、波浪能、溫差能以及鹽差能等海洋能數(shù)據(jù)缺失，因此海洋能值-生態(tài)承載力計算結果偏小。

表2 2010年山東省海洋生態(tài)系統(tǒng)可更新資源能值

表3 是山東海洋生態(tài)系統(tǒng)能值-生態(tài)足跡。由公式（5） 得出2010年總的能值-生態(tài)足跡為1.74×106hm2。海洋生態(tài)足跡的指標選取以海洋消費活動對海洋生態(tài)系統(tǒng)會造成一定的生態(tài)影響為主。計算海域能值/貨幣比率時采用的不可更新資源項目為海洋原油、海洋天然氣和海洋礦業(yè)，遠洋捕撈水產品（1.64×1021sej） 計入外界輸入的能值，則山東海洋生態(tài)系統(tǒng)能值/貨幣比率為1.93×1013sej/$。

表3 2010年山東海洋能值-生態(tài)足跡賬戶

圖2 是2000-2010年山東省海洋自然資源能值-生態(tài)足跡，包括海洋捕撈、海水養(yǎng)殖、海洋原油、海洋天然氣和海鹽。2000-2010年間，海水養(yǎng)殖生態(tài)足跡持續(xù)上升，而海洋捕撈生態(tài)足跡呈下降趨勢。海鹽足跡在這11年間上升幅度最大，海洋能源資源的生態(tài)足跡沒有太大變化。

圖2 2000-2010年山東省海洋自然資源能值-生態(tài)足跡

4 結論與討論

（1） 山東省海洋生態(tài)系統(tǒng)處于生態(tài)盈余狀態(tài)

由海洋能值-生態(tài)足跡模型得出，目前山東海域總能值-生態(tài)承載力（1.60×107hm2） 大于總能值-生態(tài)足跡（1.74×106hm2），處于生態(tài)盈余狀態(tài)。究其原因如下：1） 山東海域獨特的區(qū)位條件以及自然環(huán)境賦予其豐富的可更新能源資源，可利用的海洋生態(tài)系統(tǒng)能值比較大，具有較強的能值-生態(tài)承載力（山東海域能值密度為1.24 × 1017sej/hm2，O dum 等（1996） 計算的全球平均能值密度為3.104×1014sej/hm2，高出400 倍）；2） 區(qū)域海洋生態(tài)系統(tǒng)能量流中自然資源的能值儲量遠遠小于可更新能源資源的能值儲量，盡管山東省海洋漁業(yè)產量一直居高，但是由2010年海洋電力（風力發(fā)電、潮汐發(fā)電） 太陽能值7.43×1021sej 僅占潮汐能值的3.95%可以看出，海洋可更新能源資源利用率非常低，且從圖2 中可以看出，2001-2010年間，海洋原油、海洋天然氣的能值-生態(tài)足跡占自然資源能值-生態(tài)足跡的比例和增長幅度均非常小。因此，海洋能值-生態(tài)足跡遠小于能值-生態(tài)承載力，海洋生態(tài)系統(tǒng)處于生態(tài)盈余狀態(tài)。

（2） 山東省海洋能值-生態(tài)足跡中第二產業(yè)大于第一產業(yè)

表3 顯示，2010年山東省海洋工程建筑業(yè)能值-生態(tài)足跡（8.51×105hm2） 占總能值-生態(tài)足跡的49.02%，大于海洋自然資源的能值-生態(tài)足跡。由于近年來海洋交通運輸業(yè)、海洋旅游業(yè)、海洋科研教育服務業(yè)等海洋第三產業(yè)發(fā)展迅速，而海洋工程建筑正是這些海洋產業(yè)發(fā)展所需的基礎設施，海洋經濟快速增長的同時，必定有海洋工程建筑業(yè)的不斷增長。海洋工程建筑業(yè)不斷增長的同時，其占用的海洋系統(tǒng)空間面積也隨之增大，被占用海區(qū)以及其鄰近海區(qū)的生態(tài)系統(tǒng)平衡會遭到破壞，導致海洋自然資源產量下降，甚至一些海洋生物物種瀕臨滅絕。因此，海洋工程建筑帶來的是多方面的生態(tài)環(huán)境損害。

由圖2 可知，2000-2010年間，海洋原油、海洋天然氣的能值-生態(tài)足跡均有小幅上升，海鹽的能值-生態(tài)足跡增長幅度最大。根據(jù)《中國海洋統(tǒng)計年鑒2001-2011》，2005年之前山東省海鹽產量占全國海鹽產量比重不到50 %，2005-2010年，由56.9%增長到了69.2%。山東省海鹽產業(yè)滿足全國大部分地區(qū)海鹽需求，能值-生態(tài)足跡較大。表3中可以看出，排放到海水中的工業(yè)廢水（直排入海） 的能值-生態(tài)足跡為3 260 hm2。所以，工業(yè)廢水在一定程度上給海洋生態(tài)系統(tǒng)造成了危害，消耗了海洋生態(tài)系統(tǒng)中的部分能量。總體而言，海洋第二產業(yè)能值消耗量大于第一產業(yè)，能值-生態(tài)足跡較大。

（3） 山東省海洋水產業(yè)能值-生態(tài)足跡總體上有小幅上升，內部產業(yè)結構變化較大

圖2 顯示，2000-2010年間，海洋捕撈能值-生態(tài)足跡下降了15.4%，海水養(yǎng)殖能值-生態(tài)足跡增加了29.9%，且從2001年起，海水養(yǎng)殖業(yè)超越海洋捕撈業(yè)成為山東省海洋漁業(yè)的支柱產業(yè)。究其原因如下：1） 由于人類海洋消費活動增加導致海洋過度捕撈、海水污染等，海洋生物鏈遭到破壞，生物物種數(shù)量減少，海洋漁業(yè)資源日趨衰退，海洋捕撈量也會隨之減少。據(jù)山東省海洋與漁業(yè)廳統(tǒng)計資料顯示，山東省近海47 個主要捕撈品種，已有27 種過度利用，6 種嚴重衰退。2） 山東省致力于發(fā)展現(xiàn)代海水養(yǎng)殖技術，從1950年山東省海水養(yǎng)殖研究所的成立至2009年12月現(xiàn)代海水養(yǎng)殖產業(yè)技術創(chuàng)新戰(zhàn)略聯(lián)盟成立，從基礎理論到應用技術研究，海水養(yǎng)殖技術水平的進步以及產業(yè)規(guī)模不斷擴大，海水養(yǎng)殖產量逐年增加，使得海洋水產品產量保持穩(wěn)定。

（4） 本研究將能值-生態(tài)足跡模型運用到海洋生態(tài)系統(tǒng)中，以期能夠為海洋經濟可持續(xù)發(fā)展研究提供有效的理論方法。但是，海洋生態(tài)系統(tǒng)本身的復雜性、海洋統(tǒng)計數(shù)據(jù)的缺失以及數(shù)據(jù)本身的誤差是研究海洋生態(tài)足跡的障礙，因此，以能值-生態(tài)足跡模型來評價區(qū)域海洋經濟可持續(xù)發(fā)展狀況仍需要進一步的完善。但是總體而言，海洋能值-生態(tài)足跡模型在不斷完善的基礎上將會是一種較可靠的可持續(xù)發(fā)展評價工具，能夠為特定區(qū)域、國家以及全球海洋經濟可持續(xù)發(fā)展方向提供科學合理的依據(jù)。

Campbell D E,2000.A revised solar transformity for tidal energy received by the earth and dissipated globally: implications for emergy analysis.Center for Environmental Policy:255-264.

Campbell D E, Meisch M, Demoss T, 2004. Keeping the books for environmental systems: an energy analysis of west Virginia.Environmental Monitoring and Assessment, (94) :217-230.

Folke C, Hammer M, Jansson A M, 1991. Life -support value of ecosystems: A case study of the Baltic Sea region.Ecological Economics,3 (2) :123-137.

Folke C,Kaustsky N,Berg H,et al,1998.The ecological footprint concept for sustainable seafood production:A review Ecological Application,8 (Suppl) :63-71.

Odum H T,1996.Environmental accounting energy and decision making.New York: John&Sons:320-370.

Rees W E, 1992. Ecological footprints and appropriated carrying capacity: what urban economics leaves out. Environment and Urbanization,4 (2) :121-130.

Warren -Rhodes K, Sadovy Y, Cesaer H, 2003. Marine ecosystem appropriation in the Indo-Pacific: A case study of the live reef fish food trade.Ambio,32 (7) :481-488.

陳成忠，林振山，王暉，2008.人類最大可持續(xù)海洋足跡的模擬.生態(tài)學報，28（2）：656-660.

陳春鋒，王宏燕，肖篤寧，等，2008.基于傳統(tǒng)生態(tài)足跡方法和能值生態(tài)足跡方法的黑龍江省可持續(xù)發(fā)展狀態(tài)比較.應用生態(tài)學報，19（11）：2544-2549.

狄乾斌，韓增林，孫迎，2009.海洋經濟可持續(xù)發(fā)展能力評價及其在遼寧省的應用.資源科學，31（2）：288-294.

方欣，張亦飛，毛翰宣，等，2013.浙江省海洋經濟可持續(xù)發(fā)展能力演進研究.海洋通報，32（2）：121-126.

顧康康，劉景雙，汪洋，2009.生態(tài)足跡模型修正研究--以遼中地區(qū)礦業(yè)城市為例.中國科學院研究生院學報，26（5）：633-638.

韓增林，劉桂春，2003.海洋經濟可持續(xù)發(fā)展的定量分析.地域研究與開發(fā)，22（3）：1-4.

金丹，卞正富，2010.基于能值的生態(tài)足跡模型及其在資源型城市的應用.生態(tài)學報，30（7）：1725-1733.

藍盛芳，欽佩，2001. 生態(tài)系統(tǒng)的能值分析.應用生態(tài)學報，12，（1）：129-131.

藍盛芳，欽佩，陸宏芳，2002.生態(tài)經濟系統(tǒng)能值分析.北京：化學工業(yè)出版社，12：1-418.

李夫星，鄭穎娟，張玉，等，2013.環(huán)渤海四省市海洋可持續(xù)發(fā)展能力比較評價.海洋通報，32（3）：338-344.

劉淼，胡遠滿，常禹，等，2008.基于能值理論的生態(tài)足跡方法改進.自然資源學報，23（3）：447-457.

邵桂蘭，韓菲，李晨，2011.基于主成分分析的海洋經濟可持續(xù)發(fā)展能力測算——以山東省2000-2008年數(shù)據(jù)為例.中國海洋大學學報（社會科學版），（6）：18-22.

隋春花，陸宏芳，鄭鳳英，2006.基于能值分析的廣東省生態(tài)經濟系統(tǒng)綜合研究.應用生態(tài)學報，17（11）：2147-2152.

王建源，陳艷春，李曼華，2007，等.基于能值分析的山東省生態(tài)足跡.生態(tài)學雜志，26（9）：1505-1510.

魏明，王洪軍，柳新偉，2006.青島市生態(tài)足跡動態(tài)研究.生態(tài)環(huán)境，15（4）：761-764.

徐中民，張志強，程國棟，等，2003.中國1999年生態(tài)足跡計算與發(fā)展能力分析.應用生態(tài)學報，14（2）：280-285.

楊海波，王宗敏，趙紅領，等，2009.基于改進模型的東營市生態(tài)足跡動態(tài)分析.應用生態(tài)學報，20（7）：1753-1758.

楊海真，李愛梅，葉田，2010.基于修正的生態(tài)足跡區(qū)域可持續(xù)發(fā)展評價.同濟大學（自然科學版），38（8）：1188-1193.

楊山，王玉婷，2011.基于生態(tài)足跡修正模型的江蘇省海洋經濟可持續(xù)發(fā)展分析.應用生態(tài)學報，22（3）：748-754.

楊天姿，婁安剛，2009.基于生態(tài)足跡模型的威海市可持續(xù)發(fā)展定量研究.安徽農業(yè)科學，37（11）：5139-5141.

張芳怡，濮勵杰，張健，2006.基于能值分析理論的生態(tài)足跡模型及應用--以江蘇省為例.自然資源學報，21（4）：654-660.

趙美霞，余克服，張喬民，等，2011.應用海洋生態(tài)足跡評價三亞珊瑚礁資源的可持續(xù)利用.熱帶海洋學報，30（2）：74-80.

趙志強，高江波，李雙成，等，2009.基于能值改進生態(tài)足跡模型的廣東省1987-2006年生態(tài)經濟系統(tǒng)分析.北京大學學報（自然科學版），（1）：76-82.

趙志強，李雙成，高陽，2008.基于能值改進的開放系統(tǒng)生態(tài)足跡模型及其應用--以深圳市為例.生態(tài)學報，28（5）：2220-2231.