海南熱帶雨林國家公園土地利用碳排放時(shí)空演變

海南熱帶雨林國家公園土地利用碳排放時(shí)空演變目錄1引言 12研究區(qū)概況與數(shù)據(jù)源 22.1研究區(qū)概況 22.2數(shù)據(jù)來源 33研究方法 33.1土地利用動(dòng)態(tài)度 33.2碳排放系數(shù)法 43.3空間自相關(guān) 44結(jié)果與分析 54.1土地利用變化分析 54.2土地利用碳排放時(shí)空演變特征 64.2.1土地利用碳排放時(shí)間變化特征 64.2.2土地利用碳排放空間分布特征 74.2.3土地利用碳排放空間分異特征 85結(jié)論與建議 105.1結(jié)論 105.2建議 105.3不足 116參考文獻(xiàn) 11致謝 12摘要：本文將海南熱帶雨林國家公園作為研究對(duì)象，以1990、1999、2011、2021年四個(gè)時(shí)相的Landsat遙感影像數(shù)據(jù)作為基礎(chǔ)數(shù)據(jù)源，利用碳排放系數(shù)法、空間可視化、空間自相關(guān)等方法分析海南熱帶雨林國家公園土地利用碳排放時(shí)空分布格局與演化特征，為推動(dòng)區(qū)域經(jīng)濟(jì)社會(huì)綠色低碳發(fā)展提供決策參考。結(jié)果表明，（1）1990-2021年，草地、耕地、水體面積呈減少趨勢，建設(shè)用地、林地、水體、未利用地面積呈增加趨勢，建設(shè)用地和草地的變化程度較為顯著；（2）30年來，海南熱帶雨林國家公園的土地利用凈碳排放量整體呈先略微下降后快速上升的趨勢，累計(jì)增加了73916噸。碳源增長幅度大，碳匯小幅度增長，其中林地為主要碳匯，建設(shè)用地為主要碳源，碳排放量與建設(shè)用地具有較強(qiáng)的相關(guān)性；（3）海南熱帶雨林國家公園各鄉(xiāng)鎮(zhèn)土地利用碳排放空間上呈中間低四周高的分布格局，且高值區(qū)有從四周向中部不斷演化的趨勢。（4）空間關(guān)聯(lián)上，1990年海南熱帶雨林國家公園各鄉(xiāng)鎮(zhèn)凈碳排放和碳源在空間上具有不顯著性，在1999-2021年空間相關(guān)性不斷降低，接近于均衡分布。碳匯在1990-2021年呈現(xiàn)顯著的空間正相關(guān)關(guān)系。關(guān)鍵詞：海南熱帶雨林國家公園；土地利用變化；碳排放；空間關(guān)聯(lián)；時(shí)空演變1引言全球氣候變暖是當(dāng)今社會(huì)密切關(guān)注的問題之一，不少學(xué)者認(rèn)為人類活動(dòng)范圍的擴(kuò)大是導(dǎo)致全球二氧化碳濃度升高的主要原因[1]，有些學(xué)者則持其他看法，認(rèn)為土地利用是造成溫室氣體排放量增長的主要因素[2]。無可爭議的是，全球變暖與二氧化碳濃度上升密切相關(guān)。相關(guān)研究表明，區(qū)域土地利用活動(dòng)產(chǎn)生的碳排放量占人類活動(dòng)碳排放總量的1/3，成為區(qū)域碳排放的重要碳源，對(duì)于區(qū)域的碳排放格局有著重要影響[3]。2018年，IPCC發(fā)布報(bào)告強(qiáng)調(diào)了土地在氣候系統(tǒng)中起著重要作用。2021年，中國提出“碳達(dá)峰、碳中和”目標(biāo)，明確了碳減排任務(wù)。因此，深入研究土地利用過程中碳排放的時(shí)空演變具有重大現(xiàn)實(shí)意義。早在上世紀(jì)90年代國外學(xué)者就開始研究土地利用變化與碳排放之間的關(guān)系，21世紀(jì)初以來國內(nèi)學(xué)者展開了大量有關(guān)碳排放問題的研究并取得了一定成果。這些研究主要集中于區(qū)域碳補(bǔ)償分區(qū)、碳排放時(shí)空特征、區(qū)域碳排放核算方法、碳排放影響因素等方面。全國尺度上，葛全勝等探討了中國300年間土地利用對(duì)碳循環(huán)產(chǎn)生的影響[4]；賴力在構(gòu)建碳排放清單的同時(shí)考慮了陸地生態(tài)系統(tǒng)碳排放（碳吸收）和人為源碳排放[5]；SuM等測算了28個(gè)歐洲國家20年間的溫室氣體碳排放量[6]；SohlTL等為緩解美國溫室氣體排放構(gòu)建了土地利用變化模擬框架[7]。區(qū)域尺度上，張秀梅等[8]、嚴(yán)志翰等[9]、魏艷茹等[10]、舒心等[11]、夏四友[12]等分別對(duì)江蘇省、浙江省、福建省、長三角城市群、京津冀城市群等的碳排放和碳補(bǔ)償分區(qū)進(jìn)行了研究；城市尺度上，藍(lán)家程等[13]研究了重慶市的碳排放效應(yīng)；JoHK[14]以韓國為例分析了城市綠地生態(tài)系統(tǒng)的碳匯作用。總的來說，土地利用碳排放的定量估算、低碳優(yōu)化方法是碳排放研究的重要內(nèi)容，為區(qū)域可持續(xù)發(fā)展提供理論依據(jù)。然而，當(dāng)前國內(nèi)相關(guān)文獻(xiàn)在研究尺度上主要集中于國家、地方和省級(jí)宏觀層面，而市級(jí)和縣域微觀層面的研究相對(duì)缺乏；在研究區(qū)域的選擇上，多數(shù)研究傾向于關(guān)注經(jīng)濟(jì)發(fā)達(dá)和能源消費(fèi)量大的華東、華北、東北地區(qū)，而對(duì)于自然條件優(yōu)越、以天然生態(tài)系統(tǒng)為主的南方地區(qū)研究尚顯薄弱，有待進(jìn)一步深入。因此，關(guān)于海南熱帶雨林國家公園土地利用的碳排放時(shí)空演變的研究在國內(nèi)外都具有一定的重要性和研究價(jià)值。未來的研究可以進(jìn)一步加強(qiáng)學(xué)術(shù)交流和合作，積極突破研究瓶頸，為推動(dòng)土地利用的可持續(xù)發(fā)展和氣候變化的適應(yīng)與應(yīng)對(duì)提供更深入的理論和實(shí)踐支持?；诖耍狙芯恳院Ｄ蠠釒в炅謬夜珗@為研究單元，利用1990、1999、2011、2021年四期LandsatTM遙感影像為數(shù)據(jù)源，運(yùn)用碳排放系數(shù)法、空間可視化、空間自相關(guān)等方法分析海南熱帶雨林國家公園土地利用碳排放和碳吸收的時(shí)空分布格局與演化特征。針對(duì)土地利用碳排放狀況的區(qū)域差異，采取適宜的碳補(bǔ)償策略與方針，為推動(dòng)區(qū)域經(jīng)濟(jì)社會(huì)綠色低碳發(fā)展提供決策參考。2研究區(qū)概況與數(shù)據(jù)源2.1研究區(qū)概況海南熱帶雨林國家公園位于海南島中部山區(qū)，地處南北熱量和東西水分分界的交匯處(圖1)，生態(tài)系統(tǒng)類型以熱帶雨林為主。國家公園總面積4269km2，約占海南陸域總面積的13%，森林覆蓋率為95.85%，其中包括超過95%的未被破壞的原始森林及占到55%的大陸性島嶼型熱帶雨林。擁有中國分布最集中、保存最完好、連片面積最大的大陸性島嶼型熱帶雨林[15]。陽光充足且光線充沛，全年溫度介于22.5至26.0攝氏度之間，并且常年降水量穩(wěn)定在1759毫米左右。海南熱帶雨林國家公園坐落在海南島南側(cè)的五指山-黎母嶺，最高山峰五指山（1867m），黎母嶺最高峰為鸚哥嶺（1811m），屬于中山-丘陵地貌。主要河流有萬泉河和昌化江。其中熱帶雨林分布最廣，是該公園內(nèi)最主要的植被類型，形成了獨(dú)特、極具價(jià)值的森林景觀、動(dòng)物景觀、生物資源、清潔水源、清新空氣等生態(tài)產(chǎn)品。熱帶雨林可分為熱帶低地雨林、熱帶山地雨林、熱帶針葉林和高山云霧林。公園內(nèi)野生維管束植物220科1142屬3577種，其中國家一級(jí)保護(hù)植物坡壘、伯樂樹、海南蘇鐵、葫蘆蘇鐵、龍尾蘇鐵、臺(tái)灣蘇鐵等6種，國家二級(jí)保護(hù)植物桫欏、土沉香、降香黃檀、海南紫荊木、蝴蝶樹等34種，特有植物有尖峰青岡、霸王玉蘭、吊羅山蘿芙木、五指山含笑、海南菊、海南翠柏、雅加松等428種[15]。生物物種豐富，生物多樣性高。圖1研究區(qū)位置圖2.2數(shù)據(jù)來源本文的土地利用數(shù)據(jù)來源于地理空間數(shù)據(jù)云()，選取云量較少、覆蓋較好、夏半年的遙感影像，獲取海南熱帶雨林國家公園1990年、1999年、2011年、2021年四期Landsat分辨率為30m的遙感影像，利用envi軟件對(duì)遙感影像進(jìn)行大氣校正、輻射校正、幾何校正等預(yù)處理，然后通過監(jiān)督分類與目視解譯相結(jié)合的方法，將土地利用類型分為耕地、林地、草地、水體、建設(shè)用地和未利用地六個(gè)地類。最后，構(gòu)建分類混淆矩陣，同時(shí)使用Kappa系數(shù)對(duì)經(jīng)過解譯的遙感影像進(jìn)行精度檢驗(yàn)。結(jié)果表明，四期數(shù)據(jù)Kappa系數(shù)均超過0.85，在精度上滿足使用的要求。3研究方法3.1土地利用動(dòng)態(tài)度采用單一土地利用動(dòng)態(tài)度反映研究區(qū)一定時(shí)間范圍內(nèi)某種土地利用類型的動(dòng)態(tài)變化情況[16]。土地利用動(dòng)態(tài)度越小，土地利用動(dòng)態(tài)變化越小。具體的計(jì)算公式如式（1）：K=U其中，K表示某一研究時(shí)段某一種土地利用類型的動(dòng)態(tài)度；Ua為研究期初期某種地類的面積；Ub為研究末期某種地類的面積；t表示研究周期。3.2碳排放系數(shù)法采用碳排放系數(shù)法計(jì)算公園碳排放量，碳排放量計(jì)算公式如式（2）所示：Ci式中，Ci表示第i種地類產(chǎn)生的碳排放量（t）；Ai表示第i種地類的面積（hm2）；αi表示第i種地類的碳排放系數(shù)。當(dāng)C參考以往研究成果，結(jié)合研究區(qū)的實(shí)際情況，估算得到各土地利用類型的碳排放系數(shù)（表1）。表1不同土地利用類型碳排放系數(shù)（單位：t/hm2）土地利用類型土地利用碳排放系數(shù)參考來源耕地0.497彭文甫等[17]林地-0.581彭文甫等[17]草地-0.021孫雷剛等[18]水域-0.257賴力[5]建設(shè)用地36.07譚潔等[19]未利用地-0.005賴力[5]注：負(fù)數(shù)表示碳匯。3.3空間自相關(guān)空間自相關(guān)主要探索特定地理空間范圍內(nèi)一組變量與其他觀測變量之間潛在的相互依賴性，包括全局空間自相關(guān)和局部空間自相關(guān)[20]。采用全局Moran'sI指數(shù)分析空間關(guān)聯(lián)特征，而利用局部自相關(guān)進(jìn)行聚類分析,可以揭示研究區(qū)域內(nèi)要素的高值與低值集聚分布,可用于要素時(shí)空演變特征分析。全局Moran'sI計(jì)算公式如式（3）：I=n式中，I為全局Moran'sI值；n為研究對(duì)象的個(gè)數(shù)；xi和xj為目標(biāo)屬性特征在研究對(duì)象i和j上的觀測值；Wij為i和j之間的空間權(quán)重，鄰接是為1，反之為0；4結(jié)果與分析4.1土地利用變化分析利用envi對(duì)海南熱帶雨林國家公園1990年、1999年、2011年和2021年遙感影像解譯獲取土地利用數(shù)據(jù)。由圖2和表1可知，林地是海南熱帶雨林國家公園的主導(dǎo)土地利用類型，在各年度均占區(qū)域土地總面積的85%以上，草地和耕地面積次之。1990—2021年研究區(qū)林地面積總體上緩慢增加，水體面積前期增加較快后期緩慢增加，建設(shè)用地面積和未利用地增加幅度較大且在2011-2021年增長速率最高。1990-2021年建設(shè)用地面積增加了27.09km2，增長率最高，達(dá)474%，其次是水體面積，增加了65.89km2，增長率達(dá)450%;草地和耕地面積前期緩慢下降后期下降較快，其中草地面積下降371.33km2，降幅為96%，耕地面積下降62.89km2,降幅為66%。圖2海南熱帶雨林國家公園1990-2021年土地利用類型表2海南熱帶雨林國家公園土地利用類型面積及比例土地利用類型1990年1999年2011年2021年面積/km2比例/%面積/km2比例/%面積/km2比例/%面積/km2比例/%耕地94.502.1267.911.5358.801.3231.600.71林地3949.0488.784016.8090.274064.5291.414289.3796.40草地383.188.61311.867.01245.845.5311.850.27水體14.640.3348.711.0974.391.6780.531.81建設(shè)用地5.700.134.220.092.650.0632.800.74未利用地1.250.030.050.000.030.003.380.08由表3可知，1990-2021年海南熱帶雨林國家公園土地利用結(jié)構(gòu)不斷發(fā)生變化,土地利用變化表現(xiàn)為“四增兩減”，即林地、建設(shè)用地、未利用地和水體整體呈增加趨勢,草地和耕地整體呈減少態(tài)勢。其中建設(shè)用地增加的速度最快，動(dòng)態(tài)度達(dá)15.32%,表明土地利用開發(fā)的強(qiáng)度在增大,城鎮(zhèn)面積在快速擴(kuò)張,城鎮(zhèn)化進(jìn)程不斷加快。草地減少的速度最快，31年間動(dòng)態(tài)度為-3.13%，可能存在著過度放牧以及城市擴(kuò)張導(dǎo)致的草場退化問題。表3海南熱帶雨林國家公園各土地利用類型動(dòng)態(tài)度（單位:%）土地利用類型1990-1999年1999-2011年2011-2021年1990-2021年耕地-3.13-1.12-4.63-2.15林地50.28草地-2.07-1.76-9.52-3.13水體25.864.390.8314.52建設(shè)用地-2.89-3.10113.8315.32未利用地-10.70-2.781095.005.534.2土地利用碳排放時(shí)空演變特征4.2.1土地利用碳排放時(shí)間變化特征由公式（2）估算出海南熱帶雨林國家公園1990-2021年土地利用碳排放量和碳吸收量以及凈碳排放量（表3）。由表3可知，1990-2021年海南熱帶雨林國家公園的凈碳排放量均為負(fù)值，表現(xiàn)為碳吸收，土地利用的凈碳排放量受到碳匯和碳源的影響。從碳源方面來看，公園總碳源量在1990-2021年整體呈先略微下降后上升的趨勢，累計(jì)增加了94604t。碳源包括耕地和建設(shè)用地，其中耕地碳排放量呈持續(xù)遞減趨勢，但耕地的碳排放較弱，對(duì)整體碳排放量影響較小。建設(shè)用地碳排放量由1990年20575t增長至2021年118305t，增長5.75倍，其中2011-2021年增長速度最快達(dá)1238%。建設(shè)用地面積增加是引起公園總碳源量增加的主要因素。從碳匯方面來看，公園總碳匯量在1990-2021年呈增加趨勢，碳匯量從229439t增長至249212t。其中，林地對(duì)碳匯貢獻(xiàn)率最大，占碳匯總量的99%，林地為主要碳匯。從整體看，碳匯對(duì)碳的吸收量遠(yuǎn)遠(yuǎn)高于碳源的碳排放量，凈碳排放為負(fù)值，且在1990-2011年呈持續(xù)遞減的趨勢，2011-2021年為快速增長階段，與建設(shè)用地變化趨勢一致。因此，1990-2021年建設(shè)用地面積快速增加導(dǎo)致的碳排放量增長是導(dǎo)致海南熱帶雨林國家公園凈碳排放量增長的關(guān)鍵因素。表4海南熱帶雨林國家公園1990-2021年土地利用碳排放量（單位：102t）年份碳源碳匯凈碳排放總量耕地建設(shè)用地合計(jì)林地草地水域未利用地合計(jì)1990年46.97205.75252.72-2294.39-8.05-3.76-0.01-2306.2-2053.491999年33.75152.15185.90-2333.76-6.55-12.520.00-2352.83-2166.922011年29.2295.54124.76-2361.49-5.16-19.120.00-2385.77-2261.012021年15.711183.051198.76-2492.12-0.25-20.70-0.02-2513.09-1314.334.2.2土地利用碳排放空間分布特征為了能夠更加直觀清楚的看出海南熱帶雨林國家公園1990-2021年內(nèi)部區(qū)域間土地利用凈碳排放量空間分布差異，將研究區(qū)46個(gè)鄉(xiāng)鎮(zhèn)單元凈碳排放分為低值（-19500t~-16000t）、較低值（-16000t~-7200t）、中值（-7200t~3700t）、較高值（-3700t~-1200t）、高值(-1200t~2700t)五個(gè)等級(jí)，并通過ArcGIS進(jìn)行空間可視化分析，結(jié)果如圖4所示。整體來看，1990-2021年多年平均情形下公園土地利用凈碳排放量空間上呈中間低四周高的分布格局，且高值區(qū)呈現(xiàn)出由四周向中部逐漸擴(kuò)大的趨勢，公園凈碳排放量顯著變化區(qū)域與公園建設(shè)用地顯著增加區(qū)域基本保持一致。1990-2011年凈碳排放量總體穩(wěn)定，空間分布類型變動(dòng)較少，主要為長征鎮(zhèn)由高值區(qū)下降為較高值區(qū)，東河鎮(zhèn)由中值區(qū)突破為高值區(qū)，牙叉鎮(zhèn)由中值區(qū)突破為較高值區(qū)。2011-2021年凈碳排放量整體增長顯著，尤其以什運(yùn)鄉(xiāng)、本號(hào)鎮(zhèn)、板橋鎮(zhèn)和江邊鎮(zhèn)最為明顯，凈碳排放總量由-226100.62t上升至-131433.17t。2021年低中值區(qū)集中分布在樂東黎族自治縣的尖峰鎮(zhèn)、昌江黎族自治縣的王下鄉(xiāng)、白沙黎族自治縣的南開鄉(xiāng)、五指山市的通什鎮(zhèn)等，主要原因是多處于山地地區(qū)，土地利用程度低，城鎮(zhèn)發(fā)展水平低，林地面積比重大，主要承擔(dān)生態(tài)功能，碳吸收量大，碳匯能力強(qiáng)。高值區(qū)集中分布在經(jīng)濟(jì)發(fā)展水平高、能源消耗量高的瓊中黎族自治縣的和平鎮(zhèn)、什運(yùn)鄉(xiāng)和營根鎮(zhèn)、東方市的江邊鎮(zhèn)和東河鎮(zhèn)、樂東黎族自治縣的抱由鎮(zhèn)、昌江黎族自治縣的打安鎮(zhèn)、萬寧市的南橋鎮(zhèn)等，主要原因是土地利用程度高，人口密集，經(jīng)濟(jì)發(fā)展水平高，城鎮(zhèn)面積大,城鎮(zhèn)化進(jìn)程不斷加快，建設(shè)用地面積增加，產(chǎn)生的碳排放量大。圖31990-2021年海南熱帶雨林國家公園土地利用碳排放空間變化4.2.3土地利用碳排放空間分異特征為了檢驗(yàn)海南熱帶雨林國家公園各區(qū)域間土地利用碳排放是否存在空間相關(guān)性，本研究使用全局Moran'sI指數(shù)分別對(duì)公園各鄉(xiāng)鎮(zhèn)土地利用凈碳排放、碳源、碳匯的空間關(guān)聯(lián)度進(jìn)行分析，結(jié)果如表5、6、7所示。表5海南熱帶雨林國家公園各鄉(xiāng)鎮(zhèn)凈碳排放全局自相關(guān)統(tǒng)計(jì)分析年份1990年1999年2011年2021年Moran'sI指數(shù)0.1516200.0170130.0770570.030541Z值1.6809010.3630010.9159340.491388P值0.1093990.7166040.3597020.623152方差0.0117920.0116820.0117490.011530表6海南熱帶雨林國家公園各鄉(xiāng)鎮(zhèn)碳源全局自相關(guān)統(tǒng)計(jì)分析年份1990年1999年2011年2021年Moran'sI指數(shù)0.1557580.1364260.1341460.139523Z值1.6354551.4577181.4382291.483902P值0.1019540.1449180.1503690.137835方差0.0118430.0118450.0118210.011881由表5和表6可知，1990-2021年凈碳排放和碳源的Moran'sI指數(shù)均大于0，但P值和Z值并未通過顯著性檢驗(yàn)，可見海南熱帶雨林國家公園各鄉(xiāng)鎮(zhèn)凈碳排放和碳源整體的相關(guān)性并不顯著。其中，凈碳排放Moran'sI指數(shù)從1990年的0.151620降低到2021年的0.030541，整體呈下降趨勢，說明隨著時(shí)間的推移，空間集聚現(xiàn)象在減弱，各鄉(xiāng)鎮(zhèn)的凈碳排放差距在逐漸縮小，有均衡分布的傾向。表7海南熱帶雨林國家公園各鄉(xiāng)鎮(zhèn)碳匯全局自相關(guān)統(tǒng)計(jì)分析年份1990年1999年2011年2021年Moran'sI指數(shù)0.0073130.2333220.1039220.149787Z值0.3243342.9355752.3961641.617437P值0.7456850.0033290.0165680.105784方差0.0082930.0075780.0027710.01131由表7可知，1990-2021年海南熱帶雨林國家公園各鄉(xiāng)鎮(zhèn)碳匯的Moran'sI指數(shù)均為正值，且整體呈上升趨勢，其中1999年、2011年Z值大于臨界值，P值小于0.1，通過了95%的顯著性檢驗(yàn)，表明碳匯之間存在較強(qiáng)的空間正相關(guān)關(guān)系。為了進(jìn)一步探究46個(gè)鄉(xiāng)鎮(zhèn)之間碳匯在空間上的高低值集聚情況，對(duì)具有顯著空間自相關(guān)關(guān)系的海南熱帶雨林國家公園1999年和2011年碳匯進(jìn)行局部空間自相關(guān)分析，結(jié)果如圖4所示。圖4海南熱帶雨林國家公園碳匯LISA聚類由圖4可知，碳匯高—高集聚區(qū)1999年主要分布于白沙黎族自治縣的青松鄉(xiāng)，2011年轉(zhuǎn)為五指山市的暢好鄉(xiāng)，由于該區(qū)域林地占比高，且地形以山地丘陵為主，碳匯作用顯著。低—低集聚區(qū)主要分布在東方市的江邊鎮(zhèn)和東河鎮(zhèn)，由于該區(qū)域城鎮(zhèn)化與工業(yè)化水平均較高，建設(shè)用地占比較高，林草地占比相對(duì)較少，碳匯作用弱?？傮w來看，高—高集聚區(qū)分布在林地、草地等具有碳匯作用的生態(tài)用地占比較高的鄉(xiāng)鎮(zhèn)，低—低集聚區(qū)分布在林草地占比較少、經(jīng)濟(jì)發(fā)展水平較高、建設(shè)用地占比較高的鄉(xiāng)鎮(zhèn)，說明海南熱帶雨林國家公園碳匯在部分地區(qū)存在一定集聚性。5結(jié)論與建議5.1結(jié)論（1）林地是海南熱帶雨林國家公園的主導(dǎo)土地利用類型，1990-2021年，建設(shè)用地、林地、水體、未利用地面積呈增加趨勢，耕地、草地、水體面積呈減少趨勢，建設(shè)用地和草地的變化程度較為顯著。（2）海南熱帶雨林國家公園的土地利用碳排放量總體呈先略微下降后快速上升的趨勢。林地為主要碳匯，且占區(qū)域總面積的85%以上，因此碳匯對(duì)碳的吸收量遠(yuǎn)遠(yuǎn)高于碳源的碳排放量，凈碳排放表現(xiàn)為碳吸收，空間上呈中間低四周高的分布格局。建設(shè)用地為主要碳源，公園凈碳排放量顯著變化區(qū)域與公園建設(shè)用地顯著增加區(qū)域基本保持一致。（3）從海南熱帶雨林國家公園各鄉(xiāng)鎮(zhèn)凈碳排放量的空間分布來看，凈碳排放量高值區(qū)域主要分布在以瓊中黎族自治縣的和平鎮(zhèn)、什運(yùn)鄉(xiāng)和營根鎮(zhèn)、東方市的江邊鎮(zhèn)和東河鎮(zhèn)、樂東黎族自治縣的抱由鎮(zhèn)、昌江黎族自治縣的打安鎮(zhèn)、萬寧市的南橋鎮(zhèn)為主的經(jīng)濟(jì)發(fā)展水平高、能源消耗量高的區(qū)域，凈碳排放量低值區(qū)域主要集中在樂東黎族自治縣的尖峰鎮(zhèn)、昌江黎族自治縣的王下鄉(xiāng)、白沙黎族自治縣的南開鄉(xiāng)、五指山市的通什鎮(zhèn)為主城鎮(zhèn)發(fā)展水平低，林地面積比重大的區(qū)域。碳排放高值區(qū)和較高值區(qū)呈現(xiàn)由四周向中部不斷演化的趨勢。（4）1990年海南熱帶雨林國家公園各鄉(xiāng)鎮(zhèn)凈碳排放和碳源在空間上具有不顯著的特征，且在1999-2021年空間相關(guān)性不斷降低，接近于均衡分布。碳匯在1999-2011年呈現(xiàn)顯著的空間正相關(guān)關(guān)系。5.2建議（1）合理配置各類用地的比例，嚴(yán)格控制建設(shè)用地能源消耗，在不斷增加森林、草地、水域和未利用地碳匯的基礎(chǔ)上，重點(diǎn)開發(fā)水域和濕地碳匯，探索建立以國家公園為主體、自然保護(hù)區(qū)為基礎(chǔ)的特色自然保護(hù)地體系，增強(qiáng)生態(tài)系統(tǒng)固碳能力。(2)在東方市、陵水黎族自治縣、萬寧市、瓊中黎族苗族自治縣等碳排放相對(duì)較高的地區(qū)，著力構(gòu)建清潔低碳安全高效的能源體系，加快推進(jìn)節(jié)能降耗;提高綠色節(jié)能產(chǎn)業(yè)的比重，推動(dòng)工業(yè)領(lǐng)域低碳轉(zhuǎn)型升級(jí)和現(xiàn)代化發(fā)展；推動(dòng)建筑業(yè)綠色發(fā)展，并構(gòu)建綠色低碳綜合交通運(yùn)輸體系與低碳生活方式。5.3不足本研究在測算土地利用碳排放時(shí)，對(duì)碳排放系數(shù)的確定主要參考其他已有文獻(xiàn)，但由于研究區(qū)域在自然、社會(huì)經(jīng)濟(jì)、能源等方面存在一定差異，此各地區(qū)的碳排放系數(shù)應(yīng)有所差別，但由于碳排放系數(shù)的確定具有一定難度，探究出適用于海南熱帶雨林國家公園實(shí)際的碳排放系數(shù)是尚需解決的難題。6參考文獻(xiàn)[1]GhoshS,DindaS,DasChatterjeeN,etal.Spatial-explicitcarbonemission-sequestrationbalanceestimationandevaluationofemissionsusceptiblezonesinanEasternHimalayancityusingpressure-sensitivity-resilienceframework:Anapproachtowardsachievinglowcarboncities[J].JournalofCleanerProduction,2022,336:130417.[2]韓驥，周翔，象偉寧.土地利用碳排放效應(yīng)及其低碳管理研究進(jìn)展[J].生態(tài)學(xué)報(bào)，2016，36(4):1152-1161.[3]楊慶媛.土地利用變化與碳循環(huán)[J].中國土地科學(xué),2010,24(10):7-12.[4]葛全勝，戴君虎，何凡能，等.過去300年中國土地利用、土地覆被變化與碳循環(huán)研究[J].中國科學(xué)，2008，38(2):197-210.[5]賴力.中國土地利用的碳排放效應(yīng)研究［D].江蘇南京：南京大學(xué)，2010.[6]SuM,PauleitS,YinX,etal.GreenhousegasemissionaccountingforEUmemberstatesfrom1991to2012[J].AppliedEnergy,2016,184:759-768.[7]SohlTL,SleeterBM,ZhuZ,etal.Aland-useandland-cove