• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于LUCC的四川省主體功能區(qū)碳排放與生態(tài)補(bǔ)償研究*

    2019-01-03 03:14:44潘洪義
    關(guān)鍵詞:功能區(qū)四川省土地利用

    徐 婕, 潘洪義, 黃 佩

    ?

    基于LUCC的四川省主體功能區(qū)碳排放與生態(tài)補(bǔ)償研究*

    徐 婕, 潘洪義**, 黃 佩

    (四川師范大學(xué)西南土地資源評價(jià)與監(jiān)測教育部重點(diǎn)實(shí)驗(yàn)室/四川師范大學(xué)地理與資源科學(xué)學(xué)院 成都 610068)

    土地利用/土地覆被變化(LUCC)影響下的碳排放充分體現(xiàn)人類活動對生態(tài)環(huán)境的擾動程度, 由此切入的碳收支與補(bǔ)償研究對區(qū)域低碳經(jīng)濟(jì)與平衡發(fā)展具有重要意義。本文以2005年和2015年四川省土地利用數(shù)據(jù)為基礎(chǔ), 參考IPCC假定構(gòu)建碳排放模型, 借助GIS平臺運(yùn)用數(shù)學(xué)模型對四川省5個(gè)主體功能區(qū)碳排放及碳補(bǔ)償率進(jìn)行測度與分析, 采用造林成本法和碳稅率法對生態(tài)補(bǔ)償標(biāo)準(zhǔn)進(jìn)行測算。結(jié)果表明: 1)研究期內(nèi)四川省各主體功能區(qū)土地利用與主體功能區(qū)規(guī)劃相符并按其功能定位優(yōu)化發(fā)展。重點(diǎn)開發(fā)區(qū)域耕地占絕對優(yōu)勢, 建設(shè)用地面積的總量和增量最大; 重點(diǎn)生態(tài)功能區(qū)林地、草地占比最大。2)研究期內(nèi)各主體功能區(qū)的總碳排放量大幅度增加。重點(diǎn)開發(fā)區(qū)域碳源量、總碳排放量及其增量均居首位; 重點(diǎn)生態(tài)功能區(qū)碳吸收的能力減弱。3)各主體功能區(qū)建設(shè)用地碳排放強(qiáng)度遠(yuǎn)大于總碳排放強(qiáng)度, 研究期內(nèi)各主體功能區(qū)總碳排放強(qiáng)度均增加, 建設(shè)用地碳排放強(qiáng)度則有增有減。國家層面重點(diǎn)開發(fā)區(qū)域總碳排放強(qiáng)度最大, 建設(shè)用地碳排放強(qiáng)度很小。4)生態(tài)發(fā)展區(qū)碳補(bǔ)償率高, 經(jīng)濟(jì)發(fā)達(dá)區(qū)低。國家層面限制開發(fā)區(qū)域(重點(diǎn)生態(tài)功能區(qū))碳吸收補(bǔ)償系數(shù)最高。5)碳源區(qū)應(yīng)對作為碳匯區(qū)的甘孜藏族自治州按固碳價(jià)格進(jìn)行生態(tài)補(bǔ)償。碳排放量與各主體功能區(qū)功能定位高度相關(guān), 生態(tài)補(bǔ)償標(biāo)準(zhǔn)大致從成都平原向四周遞減, 限制開發(fā)區(qū)域生態(tài)補(bǔ)償標(biāo)準(zhǔn)偏低?;谔寂欧诺纳鷳B(tài)補(bǔ)償標(biāo)準(zhǔn)建立橫向財(cái)政轉(zhuǎn)移支付政策, 有利于四川省節(jié)能減排和平衡發(fā)展。

    土地利用/土地覆被變化; 主體功能區(qū); 碳排放; 碳源區(qū); 碳匯區(qū); 碳補(bǔ)償率; 生態(tài)補(bǔ)償

    土地利用/土地覆被變化(LUCC)影響著陸地生態(tài)系統(tǒng)的分布和結(jié)構(gòu), 改變著碳循環(huán)和碳儲量, 是控制碳積累速率的主導(dǎo)因子[1]。1850—1998年間, LUCC影響下的碳排放量約占同期人類活動影響的CO2排放量的1/3[2], 對生態(tài)系統(tǒng)的物質(zhì)循環(huán)和能量流動擾動較大, 生態(tài)系統(tǒng)的結(jié)構(gòu)、過程、功能和生態(tài)系統(tǒng)碳變化隨之改變[3], 已引起國際社會研究的廣泛關(guān)注。中國承諾到2030年單位GDP的碳排放要比2005年下降60%~65%[4], 十九大報(bào)告中也相應(yīng)提出建立健全綠色低碳循環(huán)發(fā)展經(jīng)濟(jì)體系, 均表明了中國發(fā)展低碳經(jīng)濟(jì)的態(tài)度與決心。隨著低碳經(jīng)濟(jì)走向當(dāng)前社會研究的熱點(diǎn), 發(fā)展低碳經(jīng)濟(jì)、減少碳排放也成為新時(shí)期經(jīng)濟(jì)社會發(fā)展的必然。

    國內(nèi)外學(xué)者對LUCC的碳排放從多角度展開了研究。其中國外起步早, 研究尺度普遍較大, 研究成果集中于1990—2010年, 為國內(nèi)眾多研究提供了數(shù)據(jù)與方法參考。1977年Bolin[5]研究發(fā)現(xiàn)全球林地面積變化較大引起碳排放量的變化, 此后許多學(xué)者就LUCC引起碳排放展開了一系列研究。Houghton[6]對1990年全球土地利用變化進(jìn)行相關(guān)研究, 得出耕地的減少導(dǎo)致增加3×109t的碳排放; Richards等[7]受其研究成果的啟發(fā), 比較研究南亞和東南亞的林業(yè)活動和化石燃料的碳排放, 發(fā)現(xiàn)前者是后者的2倍; 隨著技術(shù)進(jìn)步研究不斷深入, Guo等[8]采用meta分析方法研究了LUCC產(chǎn)生的土壤碳庫變化; Liu等[9]通過集成生物圈模擬框架, 采用低分辨率的LUCC數(shù)據(jù)分析加利福利亞森林、灌叢和草原的碳變化。而后研究逐漸引入人文經(jīng)濟(jì), 如Ravindranath等[10]通過碳排放效益研究了印度碳價(jià)格激勵(lì)和減排潛力之間的關(guān)系。國內(nèi)研究更為細(xì)化, 眾多學(xué)者著眼于全球、國家、區(qū)域等尺度上碳排放量的測算: 劉學(xué)榮等[11]分析了東北地區(qū)2000—2014年不同土地利用方式的碳排放效應(yīng); 也有學(xué)者從經(jīng)濟(jì)和社會發(fā)展角度探討碳排放效益及其社會經(jīng)濟(jì)應(yīng)用, 如王剛等[12]結(jié)合土地利用遙感圖像和社會經(jīng)濟(jì)發(fā)展現(xiàn)狀, 分析了2014年成都市縣域單元碳收支空間分布特征, 并通過建構(gòu)回歸模型探究其與經(jīng)濟(jì)發(fā)展的協(xié)同關(guān)系; 劉慧靈等[13]以低碳經(jīng)濟(jì)為導(dǎo)向研究福州市土地利用結(jié)構(gòu)優(yōu)化方案。隨著RS與GIS的深入開發(fā), 研究逐漸偏向于特定生態(tài)系統(tǒng)或單一要素的復(fù)雜機(jī)理研究: 崔盼盼等[14]運(yùn)用LMDI-Ⅰ加法數(shù)量分解模型分析2002—2012年中國大陸30個(gè)省的城鎮(zhèn)居民隱含碳排放及其驅(qū)動機(jī)制; 黎孔清等[15]采用回歸分析法、STIRPAT和GM(1,1)模型解析了湖南省農(nóng)地投入碳排放增長機(jī)理。而碳排放無疑存在著外部性, 生態(tài)保護(hù)補(bǔ)償則是一種外部性內(nèi)部化的手段[16]。碳排放的外部性為PES(payments for ecosystem services, 即生態(tài)系統(tǒng)服務(wù)付費(fèi))邏輯基礎(chǔ)中生態(tài)服務(wù)系統(tǒng)的一部分, 論證了碳源區(qū)對碳匯區(qū)進(jìn)行生態(tài)補(bǔ)償?shù)暮侠硇浴?1世紀(jì)開始, 隨著國內(nèi)生態(tài)補(bǔ)償概念的明確提出, 國內(nèi)學(xué)者也從碳排放存在的外部性切入論證生態(tài)保護(hù)補(bǔ)償?shù)膽?yīng)用, 如王雅敬等[17]采用碳排放評價(jià)模型估算長株潭經(jīng)濟(jì)圈碳排放及其生態(tài)補(bǔ)償。

    四川省是中國西南內(nèi)陸的經(jīng)濟(jì)和人口大省, 也是西部與東部地區(qū)重要的生態(tài)屏障。2010年四川省基本完成災(zāi)后重建, 2013年四川省人民政府印發(fā)《四川省主體功能區(qū)規(guī)劃》。2005—2015年四川省城市化、工業(yè)化的迅猛發(fā)展, 土地利用/土地覆被變化顯著, 省內(nèi)區(qū)域經(jīng)濟(jì)發(fā)展不平衡的同時(shí), 生態(tài)環(huán)境問題日益突出。本文從LUCC切入, 構(gòu)建碳排放模型, 基于主體功能區(qū)對2005年、2015年四川省碳排放及碳補(bǔ)償率進(jìn)行估算, 排除地震災(zāi)害特殊情況的波動性, 圍繞四川省主體功能區(qū)的規(guī)劃普遍性展開, 體現(xiàn)10年高速發(fā)展期四川省各主體功能區(qū)規(guī)劃影響下土地利用碳排放角度的生態(tài)環(huán)境狀態(tài)變化; 采用造林成本法和碳稅率法測算21個(gè)市(州)的生態(tài)補(bǔ)償標(biāo)準(zhǔn), 響應(yīng)國家主體功能區(qū)規(guī)劃戰(zhàn)略, 以基于碳排放的生態(tài)補(bǔ)償標(biāo)準(zhǔn)做參考, 以建立橫向財(cái)政轉(zhuǎn)移支付政策做驅(qū)動, 為四川省節(jié)能減排和平衡發(fā)展提供建議, 促進(jìn)生產(chǎn)系統(tǒng)和生活系統(tǒng)的循環(huán)鏈接。

    1 材料與方法

    1.1 研究區(qū)概況

    四川省(97°21¢~108°31¢E, 26°03¢~34°19¢N)位于中國西南內(nèi)陸, 是青藏高原地區(qū)向長江中下游平原的過渡地帶。中國大陸地勢第一、二級階梯分界線穿過四川省中部, 四川省整體地勢呈西高東低。川西為高原藏區(qū)和大小涼山彝區(qū), 地形以高原和山地為主, 海拔多在3 000 m以上, 經(jīng)濟(jì)落后, 是貧困縣集中地區(qū), 森林分布廣, 人口密度小, 是重要的生態(tài)屏障。川東以丘陵和盆地地形為主, 海拔高度500~2 000 m, 主要為經(jīng)濟(jì)發(fā)達(dá)的四川盆地、成都平原地區(qū), 城市化水平高, 產(chǎn)業(yè)發(fā)達(dá), 能源消耗量大, 人口密度大, 是經(jīng)濟(jì)的重點(diǎn)開發(fā)區(qū)域, 如成都市、德陽市、綿陽市等。四川總面積約48.6×104km2, 轄21個(gè)市(州), 183個(gè)縣(市、區(qū)), 截至2015年底, 四川省常住人口8 204萬人, GDP達(dá)30 053.1億元。

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

    1)主體功能分區(qū)是依據(jù)《四川省主體功能區(qū)規(guī)劃》, 將四川省劃分為省級層面重點(diǎn)開發(fā)區(qū)域、國家層面重點(diǎn)開發(fā)區(qū)域、國家層面限制開發(fā)區(qū)域(農(nóng)產(chǎn)品主產(chǎn)區(qū))、省級層面限制開發(fā)區(qū)域(重點(diǎn)生態(tài)功能區(qū))、國家層面限制開發(fā)區(qū)域(重點(diǎn)生態(tài)功能區(qū))。2013年調(diào)整的3個(gè)行政區(qū)主體功能歸為原所屬區(qū)劃所對應(yīng)的主體功能定位(圖1)。

    2)2005年與2015年四川省LUCC數(shù)據(jù)來源于中國資源與環(huán)境數(shù)據(jù)中心, 分辨率為1 km。依據(jù)《土地利用現(xiàn)狀分類》(GB/T21010—2017), 參考中國科學(xué)院資源環(huán)境數(shù)據(jù)庫中全國1∶100 000土地利用分類系統(tǒng), 結(jié)合四川省土地利用實(shí)際和研究需要, 構(gòu)建新的研究區(qū)土地利用分類系統(tǒng), 將四川省土地利用類型劃分為6個(gè)一級類和28個(gè)二級類(表1)。

    3)社會經(jīng)濟(jì)數(shù)據(jù)來源于2006年和2016年四川省統(tǒng)計(jì)年鑒、四川省各市(州)統(tǒng)計(jì)年鑒以及統(tǒng)計(jì)網(wǎng)站, 其中183個(gè)縣(市、區(qū))的能源消費(fèi)總量由各縣(市、區(qū))的GDP及其GDP能耗系數(shù)計(jì)算得出。

    1.3 研究方法

    1.3.1 碳排放量測算

    耕地的農(nóng)作物進(jìn)行光合作用吸收CO2, 在短期內(nèi)絕大多數(shù)農(nóng)作物又通過呼吸作用釋放CO2, 耕地的碳匯效果不明顯[18], 未利用地在研究區(qū)表現(xiàn)為微弱的碳吸收[19]。故本研究區(qū)主要表現(xiàn)為以林地、草地、水域、未利用地為主的碳吸收(碳匯)和以耕地、建設(shè)用地為主的碳排放(碳源)。碳排放測算模型可表達(dá)為:

    式中:為總碳排放量,為土地利用方式產(chǎn)生的碳排放量,T為第種土地利用方式(耕地、林地、草地、水域、未利用地)的土地面積,δ為對應(yīng)土地利用方式的碳排放(吸收)系數(shù),b為建設(shè)用地碳排放量。

    建設(shè)用地的碳排放(b)通過其利用過程中能源消耗的碳排放系數(shù)間接估算[20], 世界能源需求的80%~85%來源于化石燃料, 80%~90%的碳排放來源于化石燃料能源消耗[21]。聯(lián)系研究區(qū)實(shí)際化石能源消耗情況, 建設(shè)用地碳排放估算模型可表達(dá)為:

    式中:b為建設(shè)用地碳排放量,e為能源消耗總量,e為標(biāo)準(zhǔn)煤換算系數(shù),e為能源消耗碳排放系數(shù),c、p、n分別是煤品燃料、油品燃料和天然氣的碳排放系數(shù)。

    各土地利用類型和能源的碳吸收(排放)系數(shù)根據(jù)已有研究聯(lián)系四川省實(shí)際取最具參考價(jià)值多方經(jīng)驗(yàn)數(shù)據(jù)的平均值(表2)。

    圖1 四川省主體功能分區(qū)

    表1 四川省土地分類系統(tǒng)

    表2 不同土地利用及化石燃料的碳源(匯)碳排放(吸收)系數(shù)

    1.3.2 碳排放強(qiáng)度與碳吸收補(bǔ)償系數(shù)計(jì)算

    碳排放強(qiáng)度即單位面積的碳排放, 表示土地利用的碳排放狀況[30]。碳吸收補(bǔ)償系數(shù)表示碳匯吸收能力, 能表達(dá)某區(qū)域碳補(bǔ)償率的大小[31]。計(jì)算方法如下:

    式中: Ce為第種土地利用類型的碳排放強(qiáng)度, Ct為對應(yīng)土地利用類型的碳排放量,S為對應(yīng)土地利用類型的面積,D為區(qū)域碳吸收補(bǔ)償系數(shù), Cabs為對應(yīng)區(qū)域碳匯用地吸收量, Cemi為對應(yīng)區(qū)域碳源用地排放量。

    1.3.3 碳排放的生態(tài)補(bǔ)償計(jì)算

    碳排放生態(tài)補(bǔ)償標(biāo)準(zhǔn)的定量化主要依據(jù)碳匯價(jià)格[32]。常用的計(jì)算方法是造林成本法和碳稅率法。結(jié)合中國種植森林的實(shí)際情況, 取多種較準(zhǔn)確的中國造林成本法單位價(jià)格的平均值獲得較合理的碳匯價(jià)格, 即272.65元?t-1(C)[33], 以此為研究區(qū)碳排放生態(tài)補(bǔ)償標(biāo)準(zhǔn)的下限; 上限標(biāo)準(zhǔn)則采用國際常用的瑞典碳稅率150美元?t-1(C)[34], 以2015年的平均匯率換算, 即934.52元?t-1(C)。

    2 結(jié)果與分析

    2.1 四川省主體功能區(qū)土地利用變化

    根據(jù)2005年和2015年四川省的土地利用/覆被數(shù)據(jù), 對5個(gè)主體功能區(qū)的6個(gè)一級地類面積進(jìn)行分區(qū)統(tǒng)計(jì), 得到四川省各主體功能區(qū)土地利用情況(表3)。

    表3 2005年和2015年四川省主體功能區(qū)土地利用

    A為省級層面重點(diǎn)開發(fā)區(qū)域; B為國家層面重點(diǎn)開發(fā)區(qū)域; C為國家層面限制開發(fā)區(qū)域(農(nóng)產(chǎn)品主產(chǎn)區(qū)); D為省級層面限制開發(fā)區(qū)域(重點(diǎn)生態(tài)功能區(qū)); E為國家層面限制開發(fā)區(qū)域(重點(diǎn)生態(tài)功能區(qū))。A: key development areas at provincial level; B: key development areas at the national level; C: restricted development areas at the national level (the main agricultural production areas); D: restricted development areas at the provincial level (key ecological function areas); E: restricted development areas at the national level (key ecological functional areas).

    由表3可知, 兩個(gè)重點(diǎn)開發(fā)區(qū)域研究期內(nèi)土地利用結(jié)構(gòu)相似, 表現(xiàn)為耕地占絕對優(yōu)勢, 達(dá)到各自主體功能區(qū)面積的50%以上, 其次是林地。這是由于它們囊括的89個(gè)縣級行政區(qū)主要分布在成都平原、川南和川東北。成都平原地形平坦開闊, 土壤肥沃, 農(nóng)業(yè)發(fā)達(dá); 川南丘陵在保持原有林地的同時(shí)也發(fā)展梯田農(nóng)業(yè); 川東北臨秦巴山區(qū), 森林廣布。此外, 重點(diǎn)開發(fā)區(qū)域在研究期內(nèi)建設(shè)用地面積的總量和增量均遠(yuǎn)大于其他主體功能區(qū), 國家層面重點(diǎn)開發(fā)區(qū)域又大于省級層面。原因是重點(diǎn)開發(fā)區(qū)域包括成都、綿陽、德陽、攀枝花等四川省大部分經(jīng)濟(jì)發(fā)達(dá)地區(qū), 自身城市化水平高, 再加上國家和四川省的重點(diǎn)開發(fā)支持; 且國家層面比省級層面開發(fā)支持力度大, 這一趨勢符合四川省主體功能區(qū)定位。

    國家層面限制開發(fā)區(qū)域(農(nóng)產(chǎn)品主產(chǎn)區(qū))分布在四川省盆地邊緣和淺丘、河流流域的5大農(nóng)產(chǎn)區(qū), 研究期內(nèi)也表現(xiàn)為耕地占絕對優(yōu)勢、林地次之。省級和國家層面限制開發(fā)區(qū)域(重點(diǎn)生態(tài)功能區(qū))林地、草地約占各自主體功能區(qū)面積的70%~80%, 為川西高原藏區(qū)和大小涼山彝區(qū), 多高原山區(qū), 是四川省和長江中下游地區(qū)的重要生態(tài)屏障, 因交通不便等導(dǎo)致經(jīng)濟(jì)落后, 開發(fā)成本高, 貧困縣集中, 規(guī)劃定位為生態(tài)功能區(qū), 限制開發(fā), 保護(hù)生態(tài)。

    綜上所述, 10年來, 在各主體功能區(qū)建設(shè)用地面積增加、耕地面積減少的大背景下, 重點(diǎn)開發(fā)區(qū)域建設(shè)用地面積增量居于首位, 占用耕地、林地面積最多; 農(nóng)產(chǎn)品主產(chǎn)區(qū)耕地面積持有量保持力度最大, 基本持平; 重點(diǎn)生態(tài)功能區(qū)林地、草地總量增大, 且遠(yuǎn)大于其他主體功能區(qū)。由此可看, 四川省主體功能區(qū)規(guī)劃符合四川省各主體功能區(qū)土地利用情況, 10年來, 各主體功能區(qū)土地利用均按其功能定位優(yōu)化發(fā)展。其中, 重點(diǎn)生態(tài)功能區(qū)碳匯用地占比大, 其余均為碳源用地占比大。

    2.2 四川省主體功能區(qū)的土地利用碳排放

    2.2.1 四川省主體功能區(qū)的碳排放量

    根據(jù)碳排放模型, 計(jì)算得到2005年和2015年四川省5個(gè)主體功能區(qū)的土地利用碳排放量(表4), 為綜合分析四川省各主體功能區(qū)的碳排放情況, 對計(jì)算結(jié)果進(jìn)行空間表達(dá)(圖2)。

    表4 2005年和2015年四川省主體功能區(qū)的土地利用碳排放

    A為省級層面重點(diǎn)開發(fā)區(qū)域; B為國家層面重點(diǎn)開發(fā)區(qū)域; C為國家層面限制開發(fā)區(qū)域(農(nóng)產(chǎn)品主產(chǎn)區(qū)); D為省級層面限制開發(fā)區(qū)域(重點(diǎn)生態(tài)功能區(qū)); E為國家層面限制開發(fā)區(qū)域(重點(diǎn)生態(tài)功能區(qū))。A: key development areas at provincial level; B: key development areas at the national level; C: restricted development areas at the national level (the main agricultural production areas); D: restricted development areas at the provincial level (key ecological function areas); E: restricted development areas at the national level (key ecological functional areas).

    圖2 2005年和2015年四川省主體功能區(qū)碳排放格局

    A為省級層面重點(diǎn)開發(fā)區(qū)域; B為國家層面重點(diǎn)開發(fā)區(qū)域; C為國家層面限制開發(fā)區(qū)域(農(nóng)產(chǎn)品主產(chǎn)區(qū)); D為省級層面限制開發(fā)區(qū)域(重點(diǎn)生態(tài)功能區(qū)); E為國家層面限制開發(fā)區(qū)域(重點(diǎn)生態(tài)功能區(qū))。A: key development areas at provincial level; B: key development areas at the national level; C: restricted development areas at the national level (the main agricultural production areas); D: restricted development areas at the provincial level (key ecological function areas); E: restricted development areas at the national level (key ecological functional areas).

    由表4可知, 各主體功能區(qū)均表現(xiàn)為建設(shè)用地是主要碳源, 其次是耕地; 林地的碳匯效應(yīng)最為顯著, 草地次之。10年來, 各主體功能區(qū)的總碳排放量大幅度增加, 這是由于碳源量的大幅度增加, 尤其是建設(shè)用地的碳排放量增加所致。表明城市化和工業(yè)化中能量消耗量的增加是引起四川省土地利用總碳排放增加的重要原因。

    重點(diǎn)開發(fā)區(qū)域碳源量、總碳排放量及其增量均居首位, 其中國家層面重點(diǎn)開發(fā)區(qū)域總碳排放量大于但增量小于省級層面重點(diǎn)開發(fā)區(qū)域, 國家層面限制開發(fā)區(qū)域(農(nóng)產(chǎn)品主產(chǎn)區(qū))次之。三者總的碳排放量高達(dá)全省的90%。但就空間格局來看, 碳源碳匯比例變化不明顯(圖2)。重點(diǎn)開發(fā)區(qū)域在國家和四川省的支持下, 大力發(fā)展二、三產(chǎn)業(yè), 進(jìn)行城市化建設(shè), 發(fā)展經(jīng)濟(jì), 能耗量巨大, 作為主要碳源的建設(shè)用地占用具有碳匯效應(yīng)的林地、草地。省級層面重點(diǎn)開發(fā)區(qū)域開發(fā)潛力更大, 10年來碳排放量增長尤甚。

    重點(diǎn)生態(tài)功能區(qū)碳匯量最大, 達(dá)669.596×104t, 占四川省碳匯量的63%, 對抑制區(qū)域的總碳排放量有著重要作用。研究期內(nèi)碳源和碳匯比例發(fā)生了明顯變化, 碳匯所占比重減少20%~30%。國家層面限制開發(fā)區(qū)域(重點(diǎn)生態(tài)功能區(qū))甚至從2005年的碳匯區(qū)轉(zhuǎn)變?yōu)?015年的碳源區(qū)。說明生態(tài)功能區(qū)在保持生態(tài)林地、草地的同時(shí)更加重視城鎮(zhèn)建設(shè)和經(jīng)濟(jì)發(fā)展, 建設(shè)用地的碳排放增速快于林地、草地的碳吸收增速, 對生態(tài)功能區(qū)的主體功能定位不夠重視。由此可知, 建設(shè)用地在各主體功能區(qū)碳減排中具有重要意義, 故減排應(yīng)將重點(diǎn)放在保持林地, 降低建設(shè)用地的碳排放。

    2.2.2 四川省及其主體功能區(qū)的碳排放強(qiáng)度

    為反映四川省及其5個(gè)主體功能區(qū)土地利用碳排放能力的強(qiáng)弱, 引入碳排放強(qiáng)度指標(biāo)。土地利用的碳排放中建設(shè)用地是最主要的碳源, 單位面積建設(shè)用地碳排放量的大小能更好地反映一個(gè)區(qū)域第二、第三產(chǎn)業(yè)碳排放的水平和區(qū)域差距[35], 故此計(jì)算四川省及其5個(gè)主體功能區(qū)的建設(shè)用地碳排放強(qiáng)度(表5)。

    表5 2005年和2015年四川省及其主體功能區(qū)碳排放強(qiáng)度

    A為省級層面重點(diǎn)開發(fā)區(qū)域; B為國家層面重點(diǎn)開發(fā)區(qū)域; C為國家層面限制開發(fā)區(qū)域(農(nóng)產(chǎn)品主產(chǎn)區(qū)); D為省級層面限制開發(fā)區(qū)域(重點(diǎn)生態(tài)功能區(qū)); E為國家層面限制開發(fā)區(qū)域(重點(diǎn)生態(tài)功能區(qū))。A: key development areas at provincial level; B: key development areas at the national level; C: restricted development areas at the national level (the main agricultural production areas); D: restricted development areas at the provincial level (key ecological function areas); E: restricted development areas at the national level (key ecological functional areas).

    由表5可知, 各主體功能區(qū)建設(shè)用地碳排放強(qiáng)度遠(yuǎn)大于總碳排放強(qiáng)度, 但10年來, 各主體功能區(qū)總碳排放強(qiáng)度均增加, 建設(shè)用地碳排放強(qiáng)度則有增有減。

    國家層面重點(diǎn)開發(fā)區(qū)域總碳排放強(qiáng)度最大, 達(dá)0.124 5萬t?km-2, 遠(yuǎn)高于四川省0.024 1萬t?km-2的平均水平, 分別是省級層面重點(diǎn)開發(fā)區(qū)域、國家層面限制開發(fā)區(qū)域(農(nóng)產(chǎn)品主產(chǎn)區(qū))的2倍、4倍(表5)。但其建設(shè)用地碳排放強(qiáng)度卻很小, 低于四川省平均水平, 僅高于國家層面限制開發(fā)區(qū)域(重點(diǎn)生態(tài)功能區(qū))。這說明國家層面重點(diǎn)開發(fā)區(qū)域經(jīng)濟(jì)發(fā)展趨于成熟, 處于城市轉(zhuǎn)型和產(chǎn)業(yè)結(jié)構(gòu)的優(yōu)化升級新階段, 區(qū)域經(jīng)濟(jì)發(fā)展更為注重技術(shù)引進(jìn)和自主研發(fā), 積極創(chuàng)新以求建立低能耗、高效率、低排放的科學(xué)發(fā)展模式。

    省級層面重點(diǎn)開發(fā)區(qū)域、國家層面限制開發(fā)區(qū)域(農(nóng)產(chǎn)品主產(chǎn)區(qū))建設(shè)用地碳排放強(qiáng)度有明顯減少, 而國家層面重點(diǎn)開發(fā)區(qū)域、國家層面限制開發(fā)區(qū)域(重點(diǎn)生態(tài)功能區(qū))明顯增加。在總碳排放量隨著經(jīng)濟(jì)增長、產(chǎn)業(yè)和城市發(fā)展、能源消耗而增加的大局勢下, 卻可以通過內(nèi)部產(chǎn)業(yè)結(jié)構(gòu)調(diào)整和轉(zhuǎn)型升級, 發(fā)展循環(huán)經(jīng)濟(jì), 降低能耗和資源利用效率來提升碳排放水平, 控制建設(shè)用地碳排放強(qiáng)度, 這也將拉大區(qū)域差距。

    2.3 碳吸收補(bǔ)償系數(shù)與碳排放的生態(tài)補(bǔ)償測算

    2.3.1 碳吸收補(bǔ)償系數(shù)

    四川省主體功能區(qū)的碳排放有著巨大差異, 碳排放的外部性和PES邏輯基礎(chǔ)論證了建立補(bǔ)償機(jī)制的合理性。碳吸收補(bǔ)償系數(shù)表示某區(qū)域碳匯吸收能力和碳補(bǔ)償率的大小。根據(jù)2015年四川省5個(gè)主體功能區(qū)的碳匯和碳源量計(jì)算各主體功能區(qū)的碳吸收補(bǔ)償系數(shù), 并對其進(jìn)行空間表達(dá)(圖3)。

    由圖3可知, 四川省碳吸收補(bǔ)償系數(shù)值大致呈現(xiàn)西高東低的特征, 生態(tài)發(fā)展區(qū)碳補(bǔ)償率高, 經(jīng)濟(jì)發(fā)達(dá)區(qū)低。國家層面限制開發(fā)區(qū)域(重點(diǎn)生態(tài)功能區(qū))碳吸收補(bǔ)償系數(shù)最高, 接近1%, 說明其林地、草地碳匯吸收高, 抵消碳源量大, 碳補(bǔ)償率高。省級層面限制開發(fā)區(qū)域(重點(diǎn)生態(tài)功能區(qū))、國家層面限制開發(fā)區(qū)域(農(nóng)產(chǎn)品主產(chǎn)區(qū))、省級層面重點(diǎn)開發(fā)區(qū)域、國家層面重點(diǎn)開發(fā)區(qū)域依次遞減。重點(diǎn)開發(fā)區(qū)域值最低, 這是由開發(fā)強(qiáng)度高、經(jīng)濟(jì)總量大、建設(shè)用地多所致。國家層面限制開發(fā)區(qū)域(農(nóng)產(chǎn)品主產(chǎn)區(qū))值大小則與耕地?cái)?shù)量與農(nóng)產(chǎn)品生產(chǎn)有一定程度的關(guān)系。碳吸收補(bǔ)償系數(shù)高低值差距為四川省各主體功能區(qū)節(jié)能減排和平衡發(fā)展提供參考, 督促優(yōu)化單位GDP能耗、節(jié)能減排, 發(fā)展低碳經(jīng)濟(jì); 也一定程度上反映了各區(qū)域生態(tài)保護(hù)補(bǔ)償?shù)牧Χ取?/p>

    圖3 2015年四川省主體功能區(qū)碳吸收補(bǔ)償系數(shù)(D)

    A為省級層面重點(diǎn)開發(fā)區(qū)域; B為國家層面重點(diǎn)開發(fā)區(qū)域; C為國家層面限制開發(fā)區(qū)域(農(nóng)產(chǎn)品主產(chǎn)區(qū)); D為省級層面限制開發(fā)區(qū)域(重點(diǎn)生態(tài)功能區(qū)); E為國家層面限制開發(fā)區(qū)域(重點(diǎn)生態(tài)功能區(qū))。A: key development areas at provincial level; B: key development areas at the national level; C: restricted development areas at the national level (the main agricultural production areas); D: restricted development areas at the provincial level (key ecological function areas); E: restricted development areas at the national level (key ecological functional areas).

    2.3.2 碳排放的生態(tài)補(bǔ)償測算

    生態(tài)補(bǔ)償?shù)膶?shí)行需要通過財(cái)政轉(zhuǎn)移支付的形式實(shí)現(xiàn), 各主體功能區(qū)目前沒有相關(guān)機(jī)構(gòu)有權(quán)執(zhí)行, 故根據(jù)2015年四川省21個(gè)市(州)的碳排放量和不同固碳價(jià)格估算四川省各市(州)生態(tài)補(bǔ)償標(biāo)準(zhǔn)(區(qū)間表示), 并按市級行政區(qū)劃大部分面積所屬主體功能區(qū),同等條件下優(yōu)先考慮市級行政單位政府所在的主體功能區(qū)進(jìn)行歸類(表6)。

    表6 2015年四川省各市(州)生態(tài)補(bǔ)償標(biāo)準(zhǔn)區(qū)間

    A為省級層面重點(diǎn)開發(fā)區(qū)域; B為國家層面重點(diǎn)開發(fā)區(qū)域; C為國家層面限制開發(fā)區(qū)域(農(nóng)產(chǎn)品主產(chǎn)區(qū)); D為省級層面限制開發(fā)區(qū)域(重點(diǎn)生態(tài)功能區(qū)); E為國家層面限制開發(fā)區(qū)域(重點(diǎn)生態(tài)功能區(qū))。21個(gè)市(州)生態(tài)補(bǔ)償標(biāo)準(zhǔn)均以開區(qū)間表示, 單位為108元, 正值為應(yīng)提供的生態(tài)補(bǔ)償, 負(fù)值為應(yīng)得到的生態(tài)補(bǔ)償。A: key development areas at provincial level; B: key development areas at the national level; C: restricted development areas at the national level (the main agricultural production areas); D: restricted development areas at the provincial level (key ecological function areas); E: restricted development areas at the national level (key ecological functional areas). The ecological compensation standards of 21 cities (states) are all expressed in the open interval, the unit is 108¥ RMB, the positive value is the ecological compensation that should be provided, and the negative value is the ecological compensation that should be obtained.

    由表6可知, 2015年甘孜藏族自治州屬碳匯區(qū), 其他區(qū)域?qū)偬荚磪^(qū)。就碳排放角度而言, 碳源區(qū)應(yīng)對碳匯區(qū)進(jìn)行生態(tài)補(bǔ)償。成都市、樂山市、攀枝花市等碳源區(qū)應(yīng)向周圍碳匯區(qū)按生態(tài)補(bǔ)償范圍提供生態(tài)補(bǔ)償; 甘孜藏族自治州作為碳匯區(qū)應(yīng)接收的生態(tài)補(bǔ)償標(biāo)準(zhǔn)范圍是5.82億~20.56億元。區(qū)域碳排放與建設(shè)用地能耗量相關(guān)性大, 建設(shè)用地能耗量又關(guān)系到區(qū)域的經(jīng)濟(jì)發(fā)展水平, 故以單位GDP為參考要求區(qū)域執(zhí)行碳減排較合理, 如果不能完成減排目標(biāo), 則可通過資金補(bǔ)償對碳匯區(qū)執(zhí)行生態(tài)補(bǔ)償??梢越⒁哉疄橹鲗?dǎo), 市場為導(dǎo)向, 第三方生態(tài)環(huán)境保護(hù)組織作監(jiān)督, 通過橫向財(cái)政轉(zhuǎn)移支付等手段實(shí)現(xiàn)。

    3 結(jié)論與討論

    3.1 結(jié)論

    土地利用變化對碳排放量有重要影響, 通過對四川省5個(gè)主體功能區(qū)土地利用的碳排放和生態(tài)補(bǔ)償進(jìn)行分析, 得出以下結(jié)論:

    1)研究期四川省各主體功能區(qū)土地利用實(shí)際與主體功能區(qū)規(guī)劃相符并均按其功能定位優(yōu)化發(fā)展。在各主體功能區(qū)建設(shè)用地面積增加、耕地面積減少的大背景下, 重點(diǎn)開發(fā)區(qū)域建設(shè)用地面積的總量和增量均遠(yuǎn)大于其他主體功能區(qū); 農(nóng)產(chǎn)品主產(chǎn)區(qū)耕地面積持有量保持力度最大; 重點(diǎn)生態(tài)功能區(qū)林地、草地約占各自主體功能區(qū)面積的70%~80%。

    2)城市化和工業(yè)化是引起四川省土地利用總碳排放量增加的重要原因。研究期內(nèi)各主體功能區(qū)均表現(xiàn)為建設(shè)用地是主要碳源, 林地是主要碳匯。10年來, 各主體功能區(qū)的總碳排放量大幅度增加。重點(diǎn)開發(fā)區(qū)域碳源量、總碳排放量及其增量均居首位。重點(diǎn)生態(tài)功能區(qū)碳匯量最大, 占四川省碳匯量的63%, 但研究期內(nèi), 碳吸收的能力逐漸減弱。

    3)各主體功能區(qū)建設(shè)用地碳排放強(qiáng)度遠(yuǎn)大于總碳排放強(qiáng)度, 但10年來, 各主體功能區(qū)總碳排放強(qiáng)度均增加, 建設(shè)用地碳排放強(qiáng)度則有增有減。國家層面重點(diǎn)開發(fā)區(qū)域總碳排放強(qiáng)度最大, 建設(shè)用地碳排放強(qiáng)度很小。省級層面重點(diǎn)開發(fā)區(qū)域、國家層面限制開發(fā)區(qū)域(農(nóng)產(chǎn)品主產(chǎn)區(qū))建設(shè)用地碳排放強(qiáng)度有明顯減少, 而國家層面重點(diǎn)開發(fā)區(qū)域、國家層面限制開發(fā)區(qū)域(重點(diǎn)生態(tài)功能區(qū))明顯增加。

    4)生態(tài)發(fā)展區(qū)碳補(bǔ)償率高, 經(jīng)濟(jì)發(fā)達(dá)區(qū)低。國家層面限制開發(fā)區(qū)域(重點(diǎn)生態(tài)功能區(qū))碳吸收補(bǔ)償系數(shù)最高, 省級層面限制開發(fā)區(qū)域(重點(diǎn)生態(tài)功能區(qū))、國家層面限制開發(fā)區(qū)域(農(nóng)產(chǎn)品主產(chǎn)區(qū))、省級層面重點(diǎn)開發(fā)區(qū)域、國家層面重點(diǎn)開發(fā)區(qū)域依次遞減。

    5)碳源區(qū)應(yīng)對碳匯區(qū)進(jìn)行生態(tài)補(bǔ)償。成都市、樂山市、攀枝花市等碳源區(qū)應(yīng)向周圍碳匯區(qū)按生態(tài)補(bǔ)償范圍提供生態(tài)補(bǔ)償; 甘孜藏族自治州作為碳匯區(qū)應(yīng)接收的生態(tài)補(bǔ)償標(biāo)準(zhǔn)是5.82億~20.56億元, 可通過橫向財(cái)政轉(zhuǎn)移支付等手段實(shí)現(xiàn)。

    3.2 討論

    土地利用格局的變化影響著四大碳庫, 改變著區(qū)域固碳功能[36]。本文從LUCC切入, 以28個(gè)二級地類為數(shù)據(jù)源, 與較為常見的直接提取6個(gè)基本地類相比提高了結(jié)果的準(zhǔn)確性。再根據(jù)研究區(qū)實(shí)際在現(xiàn)有土地利用分類的基礎(chǔ)上建立新的土地分類系統(tǒng),自下而上進(jìn)行歸總。并且充分考慮到縣級層面的地區(qū)差異性, 通過構(gòu)建碳排放模型先測算四川省183個(gè)縣(市、區(qū))土地利用的碳排放, 再從5個(gè)主體功能區(qū)視角進(jìn)行定量分析。因此, 從數(shù)據(jù)處理階段就極大程度上提升了研究結(jié)果的可信度。

    此外, 本文突破了傳統(tǒng)的行政區(qū)界線, 選擇從區(qū)域特定的性質(zhì)定位出發(fā), 尊重自然規(guī)律, 響應(yīng)《四川省主體功能區(qū)規(guī)劃》, 劃分5個(gè)主體功能區(qū), 將規(guī)劃中各主體功能區(qū)定位及其生態(tài)熱點(diǎn)問題相結(jié)合, 在此基礎(chǔ)上建立生態(tài)補(bǔ)償標(biāo)準(zhǔn), 為促進(jìn)研究區(qū)的充分均衡發(fā)展提供參考。目前, 將土地利用碳排放與主體功能分區(qū)相結(jié)合的研究較為少見。文中還引入了碳吸收補(bǔ)償系數(shù)指標(biāo), 研究發(fā)現(xiàn)該指標(biāo)與生態(tài)補(bǔ)償標(biāo)準(zhǔn)具有高度正相關(guān)性, 兩者呈現(xiàn)的結(jié)果與各主體功能區(qū)碳收支結(jié)論也相符, 這為以后碳補(bǔ)償研究提供了借鑒。研究顯示的各主體功能區(qū)的主要碳源、碳匯、總碳排放量趨勢與四川省經(jīng)濟(jì)發(fā)展現(xiàn)狀和眾多研究結(jié)論一致。但值得注意的是, 碳排放強(qiáng)度指標(biāo)中, 各主體功能區(qū)總碳排放強(qiáng)度均增加, 建設(shè)用地碳排放強(qiáng)度卻有增有減。特別是國家層面重點(diǎn)開發(fā)區(qū)域, 總碳排放強(qiáng)度最大, 建設(shè)用地碳排放強(qiáng)度卻僅高于國家層面限制開發(fā)區(qū)域(重點(diǎn)生態(tài)功能區(qū)), 可以看出資金和技術(shù)對提升區(qū)域第二、第三產(chǎn)業(yè)碳排放水平有重要作用, 進(jìn)而對區(qū)域節(jié)能減排具有重要意義。

    但是, 在對土地利用碳排放的核算上, 各類碳排放系數(shù)的取值來源于對鄰近研究時(shí)限該區(qū)域研究成果的簡單修訂。今后應(yīng)根據(jù)研究區(qū)域?qū)嶋H的經(jīng)濟(jì)發(fā)展水平和能源結(jié)構(gòu)等進(jìn)行修正, 從而更加精確地核算研究區(qū)域的碳收支情況。再者, 合理的補(bǔ)償標(biāo)準(zhǔn)也應(yīng)在參考固碳價(jià)格的基礎(chǔ)上考慮到區(qū)域差異。

    [1] CANADELL J G. Land use effects on terrestrial carbon sources and sinks[J]. Science in China Series C: Life Sciences, 2002, 45(S1): 1–9

    [2] HOUGHTON R A. The annual net flux of carbon to the atmosphere from changes in land use 1850-1990[J]. Tellus B: Chemical and Physical Meteorology, 1999, 51(2): 298–313

    [3] 李蓓蓓, 方修琦, 葉瑜, 等. 中國東北地區(qū)過去300年耕地開墾導(dǎo)致的碳收支[J]. 中國科學(xué): 地球科學(xué), 2014, 44(9): 1987–1996LI B B, FANG X Q, YE Y, et al. Carbon emissions induced by cropland expansion in Northeast China during the past 300 years[J]. Science China: Earth Sciences, 2014, 44(9): 1987–1996

    [4] 吳昊玥, 何艷秋, 陳柔. 中國農(nóng)業(yè)碳排放績效評價(jià)及隨機(jī)性收斂研究——基于SBM-Undesirable模型與面板單位根檢驗(yàn)[J]. 中國生態(tài)農(nóng)業(yè)學(xué)報(bào), 2017, 25(9): 1381–1391 WU H Y, HE Y Q, CHEN R. Assessment of agricultural carbon emission performance and stochastic convergence in China using SBM-Undesirable model and panel unit root test[J]. Chinese Journal of Eco-Agriculture, 2017, 25(9): 1381–1391

    [5] BOLIN B. Changes of land biota and their importance for the carbon cycle[J]. Science, 1977, 196(4290): 613–615

    [6] HOUGHTON R A. Changes in terrestrial carbon over the last 135 years[M]//HEIMANN M. The Global Carbon Cycle. Berlin, Heidelberg: Springer, 1993: 139–157

    [7] RICHARDS J F, FLINT E P, DANIELS R C. Historic Land Use and Carbon Estimates for South and Southeast Asia: 1880-1980[R]. Oak Ridge, TN: Oak Ridge National Laboratory, 1994

    [8] GUO L B, GIFFORD R M. Soil carbon stocks and land use change: A meta analysis[J]. Global Change Biology, 2002, 8(4): 345–360

    [9] LIU J X, VOGELMANN J E, ZHU Z L, et al. Estimating California ecosystem carbon change using process model and land cover disturbance data: 1951-2000[J]. Ecological Modelling, 2011, 222(14): 2333–2341

    [10] RAVINDRANATH N H, MURTHY I K, CHATURVEDI R K, et al. Carbon forestry economic mitigation potential in India, by land classification[J]. Mitigation and Adaptation Strategies for Global Change, 2007, 12(6): 1027–1050

    [11] 劉學(xué)榮, 楊琳, 王穎, 等. 基于土地利用變化的東北地區(qū)碳排放效應(yīng)研究[J]. 水土保持通報(bào), 2017, 37(2): 107–114 LIU X R, YANG L, WANG Y, et al. Assessment of carbon emissions based on land use change in Northeast China[J]. Bulletin of Soil and Water Conservation, 2017, 37(2): 107–114

    [12] 王剛, 張華兵, 薛菲, 等. 成都市縣域土地利用碳收支與經(jīng)濟(jì)發(fā)展關(guān)系研究[J]. 自然資源學(xué)報(bào), 2017, 32(7): 1170–1182 WANG G, ZHANG H B, XUE F, et al. Relations between land use carbon budget and economic development at county level in Chengdu City[J]. Journal of Natural Resources, 2017, 32(7): 1170–1182

    [13] 劉慧靈, 伍世代, 韋素瓊, 等. 基于低碳經(jīng)濟(jì)導(dǎo)向的土地利用結(jié)構(gòu)優(yōu)化研究——以福建省福州市為例[J]. 水土保持通報(bào), 2017, 37(6): 202–208 LIU H L, WU S D, WEI S Q, et al. Optimization research of low-carbon economy oriented land use structure — A case study in Fuzhou City of Fujian Province[J]. Bulletin of Soil and Water Conservation, 2017, 37(6): 202–208

    [14] 崔盼盼, 張艷平, 張麗君, 等. 中國省域隱含碳排放及其驅(qū)動機(jī)理時(shí)空演變分析[J]. 自然資源學(xué)報(bào), 2018, 33(5): 879–892 CUI P P, ZHANG Y P, ZHANG L J, et al. Analysis on the spatial and temporal evolution of indirect carbon emissions and its driving mechanism in China[J]. Journal of Natural Resources, 2018, 33(5): 879–892

    [15] 黎孔清, 陳儉軍, 馬豆豆. 基于STIRPAT和GM(1, 1)模型的湖南省農(nóng)地投入碳排放增長機(jī)理及趨勢預(yù)測[J]. 長江流域資源與環(huán)境, 2018, 27(2): 345–352 LI K Q, CHEN J J, MA D D. Growth mechanism and trend forecast of carbon emission from farmland inputs in Hunan Province based on Stirpat and GM (1, 1) model[J]. Resources and Environment in the Yangtze Basin, 2018, 27(2): 345–352

    [16] 柳荻, 胡振通, 靳樂山. 生態(tài)保護(hù)補(bǔ)償?shù)姆治隹蚣苎芯烤C述[J]. 生態(tài)學(xué)報(bào), 2018, 38(2): 380–392 LIU D, HU Z T, JIN L S. Review on analytical framework of eco-compensation[J]. Acta Ecologica Sinica, 2018, 38(2): 380–392

    [17] 王雅敬, 謝炳庚, 李曉青, 等. 長株潭經(jīng)濟(jì)圈土地利用變化的碳排放量及生態(tài)補(bǔ)償研究[J]. 生態(tài)科學(xué), 2017, 36(4): 160–165 WANG Y J, XIE B G, LI X Q, et al. Research on benefit carbon emission and countervail ecological system of Chang-Zhu-Tan land exploit[J]. Ecological Science, 2017, 36(4): 160–165

    [18] GREEN C. Potential scale-related problems in estimating the costs of CO2mitigation policies[J]. Climatic Change, 2000, 44(3): 331–349

    [19] 彭文甫, 周介銘, 徐新良, 等. 基于土地利用變化的四川省碳排放與碳足跡效應(yīng)及時(shí)空格局[J]. 生態(tài)學(xué)報(bào), 2016, 36(22): 7244–7259 PENG W F, ZHOU J M, XU X L, et al. Effect of land use changes on the temporal and spatial patterns of carbon emissions and carbon footprints in the Sichuan Province of Western China, from 1990 to 2010[J]. Acta Ecologica Sinica, 2016, 36(22): 7244–7259

    [20] 方精云, 郭兆迪, 樸世龍, 等. 1981—2000年中國陸地植被碳匯的估算[J]. 中國科學(xué) D輯: 地球科學(xué), 2007, 37(6): 804–812 FANG J Y, GUO Z D, PIAO S L, et al. Terrestrial vegetation carbon sinks in China, 1981-2000[J]. Science in ChinaSeries D: Earth Sciences, 2007, 37(6): 804–812

    [21] 王秋賢, 王登杰, 顏曉妹, 等. 基于脫鉤理論的煙臺市碳排放效應(yīng)分析[J]. 水土保持通報(bào), 2015, 35(2): 313–318WANG Q X, WANG D J, YAN X M, et al. An analysis on carbon emission effects in Yantai City based on decoupling theory[J]. Bulletin of Soil and Water Conservation, 2015, 35(2): 313–318

    [22] 李穎, 黃賢金, 甄峰. 江蘇省區(qū)域不同土地利用方式的碳排放效應(yīng)分析[J]. 農(nóng)業(yè)工程學(xué)報(bào), 2008, 24(S2): 102–107 LI Y, HUANG X J, ZHEN F. Effects of land use patterns on carbon emission in Jiangsu Province[J]. Transactions of the CSAE, 2008, 24(S2): 102–107

    [23] 蘇雅麗, 張艷芳. 陜西省土地利用變化的碳排放效益研究[J]. 水土保持學(xué)報(bào), 2011, 25(1): 152–156SU Y L, ZHANG Y F. Study on effects of carbon emission by land use patterns of Shaanxi Province[J]. Journal of Soil and Water Conservation, 2011, 25(1): 152–156

    [24] 何勇. 中國氣候、陸地生態(tài)系統(tǒng)碳循環(huán)研究[M]. 北京: 氣象出版社, 2006HE Y. Research on Carbon Cycle in Climate and Terrestrial Ecosystems in China[M]. Beijing: Meteorological Press, 2006

    [25] 賴力, 黃賢金, 劉偉良, 等. 基于投入產(chǎn)出技術(shù)的區(qū)域生態(tài)足跡調(diào)整分析——以2002年江蘇為例[J]. 生態(tài)學(xué)報(bào), 2006, 26(4): 1285–1292LAI L, HUANG X J, LIU W L, et al. Adjustment for regional ecological footprint based on input-output technique: A case study of Jiangsu Province in 2002[J]. Acta Ecologica Sinica, 2006, 26(4): 1285–1292

    [26] 段曉男, 王效科, 逯非, 等. 中國濕地生態(tài)系統(tǒng)固碳現(xiàn)狀和潛力[J]. 生態(tài)學(xué)報(bào), 2008, 28(2): 463–469DUAN X N, WANG X K, LU F, et al. Carbon sequestration and its potential by wetland ecosystems in China[J]. Acta Ecologica Sinica, 2008, 28(2): 463–469

    [27] 賴力. 中國土地利用的碳排放效應(yīng)研究[D]. 南京: 南京大學(xué), 2010 LAI L. Carbon emission effect of land use in China[D]. Nanjing: Nanjing University, 2010

    [28] 徐國泉, 劉則淵, 姜照華. 中國碳排放的因素分解模型及實(shí)證分析: 1995—2004[J]. 中國人口·資源與環(huán)境, 2006, 16(6): 158–161XU G Q, LIU Z Y, JIANG Z H. Decomposition model and empirical study of carbon emissions for China, 1995–2004[J]. China Population, Resources and Environment, 2006, 16(6): 158–161

    [29] 張秀梅, 李升峰, 黃賢金, 等. 江蘇省1996年至2007年碳排放效應(yīng)及時(shí)空格局分析[J]. 資源科學(xué), 2010, 32(4): 768–775ZHANG X M, LI S F, HUANG X J, et al. Effects of carbon emissions and their spatio-temporal patterns in Jiangsu Province from 1996 to 2007[J]. Resources Science, 2010, 32(4): 768–775

    [30] TAO Y G, HUANG Z F. Review of accounting for carbon dioxide emissions from tourism at different spatial scales[J]. Acta Ecologica Sinica, 2014, 34(5): 246–254

    [31] CAI Z C, KANG G D, TSURUTA H, et al. Estimate of CH4emissions from year-round flooded rice fields during rice growing season in China[J]. Pedosphere, 2005, 15(1): 66–71

    [32] MA C B, STERN D I. China’s changing energy intensity trend: a decomposition analysis[J]. Energy Economics, 2008, 30(3): 1037–1053

    [33] 謝守紅, 王利霞, 邵珠龍. 國內(nèi)外碳排放研究綜述[J]. 干旱區(qū)地理, 2014, 37(4): 720–730XIE S H, WANG L X, SHAO Z L. Review on carbon emissions researches at home and abroad[J]. Arid Land Geography, 2014, 37(4): 720–730

    [34] 孫賢斌, 傅先蘭, 倪建華, 等. 安徽省會經(jīng)濟(jì)圈碳排放強(qiáng)度與生態(tài)補(bǔ)償研究[J]. 地域研究與開發(fā), 2012, 31(1): 135–138SUN X B, FU X L, NI J H, et al. Intensity of carbon emissions and ecological compensation in the city ecology circle in Anhui Province[J]. Areal Research and Development, 2012, 31(1): 135–138

    [35] 王亮. 鹽城市土地利用變化碳排放效益及影響因素分解研究[J]. 生態(tài)科學(xué), 2015, 34(5): 122–128 WANG L. Research for decomposition results of impact factors about carbon emission by land use patterns in Yancheng City[J]. Ecological Science, 2015, 34(5): 122–128

    [36] 劉建華, 許皞, 王耀, 等. 基于土地利用格局變化的生態(tài)風(fēng)險(xiǎn)與固碳功能評價(jià)——以河北省黃驊市為例[J]. 中國生態(tài)農(nóng)業(yè)學(xué)報(bào), 2018, 26(8): 1217–1226 LIU J H, XU H, WANG Y, et al. Evaluation of ecological risk and carbon fixation from land use change: A case study of Huanghua City, Hebei Province[J]. Chinese Journal of Eco-Agriculture, 2018, 26(8): 1217–1226

    Carbon emission and ecological compensation of main functional areas in Sichuan Provincebased on LUCC*

    XU Jie, PAN Hongyi**, HUANG Pei

    (Key Laboratory of Land Resources Evaluation and Monitoring in Southwest, Sichuan Normal University / College of Geography and Resources Science, Sichuan Normal University, Chengdu 610068, China)

    Carbon emissions under land use/cover change (LUCC) reflect the disturbance degree of environment by human activities. Therefore, carbon budget and compensation research was of great significance to regional low-carbon economy and balanced development. Based on land use data for Sichuan Province in 2005 and 2015, this paper built a new local land use classification system and carbon emission model to quantitatively analyze carbon emission and carbon offset rate of land use in 183 counties in Sichuan Province. Taking into account regional differences at county level and breaking traditional administrative regional boundaries, the study based on the five main functional areas of Sichuan Province to measure ecological compensation standards with afforestation cost and carbon tax rate methods. The results showed that: 1) during the study period, land use in each main functional area in Sichuan Province was consistent with the main functional area planning and was optimized according to the functional orientation. Cultivated lands in the key development area had absolute advantage, while the total amount and increase in construction land area was the largest. Otherwise, cultivated land holding capacity in the main agricultural production area remained the strongest. Woodland and grassland accounted for 70%–80% in the key ecological function areas. 2) During the study period, total carbon emission in each main functional area increased significantly, and the construction land was the main carbon source, and forest land was the main carbon sink. Total carbon emissions and increments in the key development area ranked first, with carbon sinks in the key ecological functional area accounting for 63% of Sichuan Province, although the ability to absorb carbon gradually weakened. 3) Carbon emission intensity of construction land of each main functional area was much greater than total carbon emission intensity. Total carbon emission intensity of each main functional area increased during the investigated period, but carbon emission intensity of construction land showed increased or decreased. Total carbon emission intensity of the key development area at national level was largest, while that of construction land was smallest. 4) The rate of carbon compensation in the ecological development areas was high, with economic development area having the reverse trend. The compensation coefficient of carbon absorption was highest for the restricted development area at national level (key ecological function area). 5) Carbon source area provided ecological compensation for carbon sink area based on the scope of ecological compensation. Ecological compensation standard for Ganzi Tibetan Autonomous Prefecture, which was a carbon sink area, was 5.82 billion to 20.56 billion Yuan. It indicated that carbon emission was correlated with functional location of each functional area. Ecological compensation standard generally decreased from Chengdu Plain to the surrounding areas, and ecological compensation standard of the restricted development areas were relatively low. On the basis of the ecological compensation standard for carbon emissions, a horizontal fiscal transfer payment policy was established as reference basis for energy conservation, emission reduction and balanced development of Sichuan Province.

    Land use/cover change; Main functional area; Carbon emission; Carbon sources area; Carbon sink area; Carbon offset rate; Ecological compensation

    , E-mail: panhongyi80@163.com

    May 22, 2018;

    Aug. 14, 2018

    F301.24

    A

    2096-6237(2019)01-0142-11

    10.13930/j.cnki.cjea.180493

    徐婕, 潘洪義, 黃佩. 基于LUCC的四川省主體功能區(qū)碳排放與生態(tài)補(bǔ)償研究[J]. 中國生態(tài)農(nóng)業(yè)學(xué)報(bào)(中英文), 2019, 27(1): 142-152

    XU J, PAN H Y, HUANG P. Carbon emission and ecological compensation of main functional areas in Sichuan Province based on LUCC[J]. Chinese Journal of Eco-Agriculture, 2019, 27(1): 142-152

    * 國家自然科學(xué)基金項(xiàng)目(41371125)、四川省教育廳項(xiàng)目(16ZB0061)、四川省哲學(xué)社會科學(xué)規(guī)劃項(xiàng)目(SC18B095)和教育部規(guī)劃基金項(xiàng)目(18XJA630005)資助

    潘洪義, 主要從事土地利用與評價(jià)研究。E-mail: panhongyi80@163.com

    徐婕, 主要從事土地利用與生態(tài)補(bǔ)償研究。E-mail: xxujie@foxmail.com

    2018-05-22

    2018-08-14

    * This study was supported by the National Natural Science Foundation of China (41371125), the Project of Department of Education of Sichuan Province (16ZB0061), Sichuan Philosophy and Social Science Planning Project (SC18B095) and the Planning Fund Project of Ministry of Education of China (18XJA630005).

    猜你喜歡
    功能區(qū)四川省土地利用
    前三季度四川省五大支柱產(chǎn)業(yè)保持平穩(wěn)較快增長
    四川冶金(2019年5期)2019-12-23 09:04:34
    四川省土木建筑學(xué)會
    四川建筑(2019年6期)2019-07-21 02:59:28
    燒結(jié)混合機(jī)各功能區(qū)的研究與優(yōu)化
    山東冶金(2019年3期)2019-07-10 00:53:54
    四川省 多舉措增強(qiáng)2500萬 農(nóng)民工獲得感
    衡水市不同功能區(qū)土壤營養(yǎng)元素的比較
    土地利用生態(tài)系統(tǒng)服務(wù)研究進(jìn)展及啟示
    濱??h土地利用挖潛方向在哪里
    水功能區(qū)納污能力計(jì)算及污染物總量控制
    主體功能區(qū)規(guī)劃實(shí)施需要解決好四個(gè)問題
    論低碳經(jīng)濟(jì)與轉(zhuǎn)變土地利用方式
    精品久久久久久成人av| 国产伦在线观看视频一区| 人妻制服诱惑在线中文字幕| 免费高清在线观看视频在线观看| 少妇的逼好多水| 久久97久久精品| 麻豆av噜噜一区二区三区| 中文字幕久久专区| 精品一区二区三区人妻视频| 国产成人福利小说| 99久国产av精品国产电影| 国产久久久一区二区三区| 精品少妇黑人巨大在线播放| 街头女战士在线观看网站| 亚洲欧洲日产国产| 51国产日韩欧美| 久久久久久久久久久丰满| 熟妇人妻久久中文字幕3abv| 五月伊人婷婷丁香| 色播亚洲综合网| 国产单亲对白刺激| 国产毛片a区久久久久| 91aial.com中文字幕在线观看| a级毛片免费高清观看在线播放| 亚洲精品日本国产第一区| 别揉我奶头 嗯啊视频| 国产精品99久久久久久久久| 看非洲黑人一级黄片| 综合色丁香网| 欧美成人精品欧美一级黄| 国产伦理片在线播放av一区| 免费观看a级毛片全部| 国产精品女同一区二区软件| 久久草成人影院| 91精品一卡2卡3卡4卡| 我的老师免费观看完整版| 日韩制服骚丝袜av| 能在线免费看毛片的网站| 久久精品国产亚洲av天美| www.av在线官网国产| 久久国产乱子免费精品| 午夜福利视频精品| 十八禁国产超污无遮挡网站| 边亲边吃奶的免费视频| 搡老乐熟女国产| 亚洲精品第二区| 国产激情偷乱视频一区二区| 蜜臀久久99精品久久宅男| 国产亚洲91精品色在线| 美女cb高潮喷水在线观看| 中文字幕人妻熟人妻熟丝袜美| 波野结衣二区三区在线| 亚洲一级一片aⅴ在线观看| 丰满少妇做爰视频| 天堂网av新在线| 天堂中文最新版在线下载 | 国产三级在线视频| 视频中文字幕在线观看| 午夜精品国产一区二区电影 | 一级毛片黄色毛片免费观看视频| 天堂√8在线中文| 毛片女人毛片| 午夜视频国产福利| 国产av在哪里看| 国产精品精品国产色婷婷| 三级国产精品片| 久久久久免费精品人妻一区二区| 国产精品无大码| 国产黄色小视频在线观看| 熟女人妻精品中文字幕| 成人亚洲精品av一区二区| 久久久久久久久久成人| 精品熟女少妇av免费看| 七月丁香在线播放| 麻豆成人午夜福利视频| 国产亚洲精品av在线| 国产av在哪里看| 综合色丁香网| 又爽又黄a免费视频| 在线a可以看的网站| 国产成人91sexporn| 一级毛片我不卡| 国产综合懂色| 国内少妇人妻偷人精品xxx网站| 欧美另类一区| 建设人人有责人人尽责人人享有的 | 丰满乱子伦码专区| 又爽又黄无遮挡网站| 不卡视频在线观看欧美| 亚洲av.av天堂| 免费观看的影片在线观看| 国产免费又黄又爽又色| 日日啪夜夜撸| 国产视频首页在线观看| 特大巨黑吊av在线直播| 99视频精品全部免费 在线| 最近最新中文字幕免费大全7| 欧美xxxx黑人xx丫x性爽| 亚洲美女搞黄在线观看| 精品不卡国产一区二区三区| 久久久久久久午夜电影| 91精品国产九色| 免费大片黄手机在线观看| 久久久久精品久久久久真实原创| 国内精品宾馆在线| 国产欧美日韩精品一区二区| 九九在线视频观看精品| 精品熟女少妇av免费看| 一级毛片黄色毛片免费观看视频| 国产永久视频网站| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 久久精品国产亚洲av涩爱| 波野结衣二区三区在线| 色5月婷婷丁香| 国内精品一区二区在线观看| 亚洲av日韩在线播放| 26uuu在线亚洲综合色| 久久久久久国产a免费观看| 毛片一级片免费看久久久久| 干丝袜人妻中文字幕| 波多野结衣巨乳人妻| 色网站视频免费| 亚洲真实伦在线观看| 日韩三级伦理在线观看| 男女视频在线观看网站免费| 秋霞在线观看毛片| 久久精品熟女亚洲av麻豆精品 | 久久久久久伊人网av| 亚洲va在线va天堂va国产| 成人综合一区亚洲| 2021天堂中文幕一二区在线观| 日本av手机在线免费观看| 午夜福利成人在线免费观看| 1000部很黄的大片| 淫秽高清视频在线观看| 91精品国产九色| 亚洲成人av在线免费| 黑人高潮一二区| 乱系列少妇在线播放| 在线播放无遮挡| 日本色播在线视频| 男人和女人高潮做爰伦理| 欧美变态另类bdsm刘玥| 国产精品av视频在线免费观看| 国产亚洲av嫩草精品影院| 午夜视频国产福利| 菩萨蛮人人尽说江南好唐韦庄| 一本一本综合久久| 亚洲欧美日韩东京热| 日日撸夜夜添| 天堂俺去俺来也www色官网 | 一二三四中文在线观看免费高清| 男女那种视频在线观看| 亚洲欧美精品专区久久| 亚洲精品国产av成人精品| 成年女人在线观看亚洲视频 | 一夜夜www| av线在线观看网站| 久久热精品热| 性插视频无遮挡在线免费观看| 日本一本二区三区精品| 日韩欧美精品v在线| 亚洲婷婷狠狠爱综合网| 亚洲精品一二三| 日韩中字成人| 午夜爱爱视频在线播放| 国产成人精品福利久久| 亚洲国产欧美人成| 色尼玛亚洲综合影院| 国产精品综合久久久久久久免费| 午夜亚洲福利在线播放| 国产色婷婷99| 欧美bdsm另类| 51国产日韩欧美| 九草在线视频观看| 激情五月婷婷亚洲| 亚洲人成网站在线播| 国产成人精品久久久久久| 最近视频中文字幕2019在线8| 精品99又大又爽又粗少妇毛片| 天堂中文最新版在线下载 | 性插视频无遮挡在线免费观看| 日韩电影二区| 日本wwww免费看| 男人爽女人下面视频在线观看| 搞女人的毛片| 午夜激情久久久久久久| 久热久热在线精品观看| 免费观看在线日韩| 少妇猛男粗大的猛烈进出视频 | 大话2 男鬼变身卡| 深爱激情五月婷婷| 夜夜爽夜夜爽视频| 亚洲精品日韩av片在线观看| 国产老妇伦熟女老妇高清| 欧美激情国产日韩精品一区| 99久久中文字幕三级久久日本| 亚洲精品成人av观看孕妇| 亚洲av二区三区四区| 麻豆国产97在线/欧美| 91精品国产九色| 欧美三级亚洲精品| 成人午夜高清在线视频| 国产在线男女| 国产精品av视频在线免费观看| 国产美女午夜福利| 国产黄片美女视频| 在线观看美女被高潮喷水网站| 尤物成人国产欧美一区二区三区| 久久国产乱子免费精品| 午夜激情欧美在线| 欧美区成人在线视频| 嫩草影院精品99| 日日啪夜夜爽| eeuss影院久久| 在现免费观看毛片| 久久久久久久久久久免费av| 国国产精品蜜臀av免费| 天堂影院成人在线观看| 哪个播放器可以免费观看大片| 国产亚洲精品久久久com| 欧美激情国产日韩精品一区| 一级毛片电影观看| 国产精品麻豆人妻色哟哟久久 | 亚洲天堂国产精品一区在线| eeuss影院久久| 国产精品一区www在线观看| 黄色配什么色好看| 中国国产av一级| 尤物成人国产欧美一区二区三区| av在线观看视频网站免费| 婷婷色综合www| 国产伦一二天堂av在线观看| 国产一区二区亚洲精品在线观看| 乱码一卡2卡4卡精品| 亚洲成人精品中文字幕电影| 在现免费观看毛片| 精品不卡国产一区二区三区| 青青草视频在线视频观看| 国产成人精品一,二区| 五月伊人婷婷丁香| 欧美xxxx性猛交bbbb| 亚洲成色77777| 青春草亚洲视频在线观看| 亚洲精品乱码久久久v下载方式| 亚洲av国产av综合av卡| 亚洲国产精品成人久久小说| 亚洲成人av在线免费| 一级毛片久久久久久久久女| 联通29元200g的流量卡| 99久久中文字幕三级久久日本| 婷婷色av中文字幕| 神马国产精品三级电影在线观看| 国产亚洲av嫩草精品影院| 久久鲁丝午夜福利片| 高清av免费在线| 国产91av在线免费观看| 日本与韩国留学比较| 99久国产av精品国产电影| 亚洲精华国产精华液的使用体验| 91av网一区二区| 高清av免费在线| 亚洲国产欧美在线一区| 99热6这里只有精品| 三级经典国产精品| 夫妻午夜视频| 亚洲av不卡在线观看| av在线蜜桃| 久久精品久久精品一区二区三区| 国产成人a区在线观看| 欧美日韩在线观看h| 97热精品久久久久久| 夜夜爽夜夜爽视频| 婷婷色麻豆天堂久久| 日韩三级伦理在线观看| 国产亚洲精品久久久com| 国产亚洲av片在线观看秒播厂 | 免费观看a级毛片全部| 性插视频无遮挡在线免费观看| 最近最新中文字幕大全电影3| 亚洲欧美一区二区三区黑人 | 精品一区在线观看国产| a级毛色黄片| 大香蕉久久网| 国产成年人精品一区二区| 久99久视频精品免费| 久久久久免费精品人妻一区二区| 草草在线视频免费看| 99久久精品国产国产毛片| 日韩一区二区视频免费看| 日产精品乱码卡一卡2卡三| 18禁在线播放成人免费| 青春草国产在线视频| 亚洲国产日韩欧美精品在线观看| 国产黄频视频在线观看| 夫妻午夜视频| 免费观看a级毛片全部| 日本av手机在线免费观看| 韩国高清视频一区二区三区| 国产伦在线观看视频一区| 亚洲美女视频黄频| 毛片女人毛片| 国产成人午夜福利电影在线观看| 欧美日韩精品成人综合77777| 午夜老司机福利剧场| 啦啦啦韩国在线观看视频| 中文欧美无线码| 91精品伊人久久大香线蕉| 日韩视频在线欧美| 少妇人妻一区二区三区视频| 亚洲精品影视一区二区三区av| 成人综合一区亚洲| 国产黄频视频在线观看| 久久午夜福利片| 欧美潮喷喷水| 国产午夜精品久久久久久一区二区三区| 能在线免费看毛片的网站| 久久这里有精品视频免费| 亚洲精品成人久久久久久| 午夜免费激情av| kizo精华| 最近视频中文字幕2019在线8| 亚洲三级黄色毛片| 国产亚洲91精品色在线| 精品国内亚洲2022精品成人| 精品久久久久久久久av| 国产欧美日韩精品一区二区| 久久久久久国产a免费观看| 中文字幕久久专区| 有码 亚洲区| 日韩视频在线欧美| 精品久久久久久久久av| 日韩视频在线欧美| 国产一区亚洲一区在线观看| 国产高清三级在线| 亚州av有码| 中文字幕人妻熟人妻熟丝袜美| 亚洲av成人精品一二三区| 欧美日韩亚洲高清精品| 日韩伦理黄色片| 如何舔出高潮| 日韩精品有码人妻一区| 熟妇人妻不卡中文字幕| 亚洲怡红院男人天堂| 99热这里只有精品一区| 国产精品av视频在线免费观看| 亚洲av电影不卡..在线观看| 五月伊人婷婷丁香| 国产黄色视频一区二区在线观看| 久久久精品免费免费高清| 国产乱人偷精品视频| 亚洲欧美精品自产自拍| 久久精品久久精品一区二区三区| 人妻系列 视频| 久久久久国产网址| 亚洲18禁久久av| 国产中年淑女户外野战色| 人妻制服诱惑在线中文字幕| 久久久成人免费电影| 亚洲欧美一区二区三区国产| 大香蕉久久网| 有码 亚洲区| 99热这里只有精品一区| 中文字幕av成人在线电影| 性插视频无遮挡在线免费观看| 久久久精品免费免费高清| 婷婷六月久久综合丁香| 亚洲精品第二区| 人体艺术视频欧美日本| 一级二级三级毛片免费看| 乱码一卡2卡4卡精品| 亚洲乱码一区二区免费版| 18禁裸乳无遮挡免费网站照片| 少妇的逼水好多| 丰满人妻一区二区三区视频av| 欧美丝袜亚洲另类| 免费观看在线日韩| 日日摸夜夜添夜夜爱| 午夜免费男女啪啪视频观看| 国产 亚洲一区二区三区 | 亚洲av免费高清在线观看| 国产午夜福利久久久久久| 2021少妇久久久久久久久久久| 国产熟女欧美一区二区| 高清欧美精品videossex| 天堂中文最新版在线下载 | 免费在线观看成人毛片| 国产一区二区三区av在线| 国产亚洲一区二区精品| 日本午夜av视频| 欧美极品一区二区三区四区| 草草在线视频免费看| 日日啪夜夜爽| 亚洲国产欧美在线一区| 波多野结衣巨乳人妻| 舔av片在线| 久久久久免费精品人妻一区二区| 国产高清有码在线观看视频| 黄色日韩在线| 久久国内精品自在自线图片| 天天躁夜夜躁狠狠久久av| 秋霞在线观看毛片| 婷婷色综合大香蕉| 又爽又黄a免费视频| 在线免费观看的www视频| 日本一二三区视频观看| 晚上一个人看的免费电影| 亚洲电影在线观看av| 色尼玛亚洲综合影院| 精品久久久久久久久亚洲| 国产亚洲一区二区精品| 精品人妻偷拍中文字幕| 插阴视频在线观看视频| 非洲黑人性xxxx精品又粗又长| 久久久久久久久久成人| 欧美精品国产亚洲| 亚洲一区高清亚洲精品| 日本黄大片高清| 一区二区三区免费毛片| 伊人久久精品亚洲午夜| 91久久精品国产一区二区三区| 国产淫语在线视频| 热99在线观看视频| 国产亚洲一区二区精品| 日韩一本色道免费dvd| 国产精品一区www在线观看| 成人毛片60女人毛片免费| 美女内射精品一级片tv| 内射极品少妇av片p| 2018国产大陆天天弄谢| 超碰97精品在线观看| 国产精品福利在线免费观看| 天堂俺去俺来也www色官网 | 永久网站在线| av卡一久久| 久久久欧美国产精品| 国产欧美另类精品又又久久亚洲欧美| 欧美 日韩 精品 国产| 久久久午夜欧美精品| 国产片特级美女逼逼视频| 插逼视频在线观看| 久久韩国三级中文字幕| 成人高潮视频无遮挡免费网站| 亚洲精品国产成人久久av| 午夜精品一区二区三区免费看| 搡女人真爽免费视频火全软件| 嘟嘟电影网在线观看| 久久久午夜欧美精品| 久久精品夜色国产| 1000部很黄的大片| 成人无遮挡网站| 久久久久久久久久成人| 国国产精品蜜臀av免费| 日本三级黄在线观看| 日韩制服骚丝袜av| 国产午夜精品论理片| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产免费视频播放在线视频 | 91在线精品国自产拍蜜月| 成人一区二区视频在线观看| 免费电影在线观看免费观看| 精品99又大又爽又粗少妇毛片| 欧美日韩在线观看h| 国产视频内射| 国产极品天堂在线| 狂野欧美白嫩少妇大欣赏| 91久久精品电影网| 国产精品熟女久久久久浪| 免费观看在线日韩| 亚洲综合精品二区| 国产老妇女一区| 一区二区三区免费毛片| 成年版毛片免费区| 欧美激情在线99| 免费大片黄手机在线观看| 天美传媒精品一区二区| 成人欧美大片| 日本wwww免费看| 久久草成人影院| 中文欧美无线码| 日韩人妻高清精品专区| 亚洲四区av| 精品熟女少妇av免费看| 欧美成人a在线观看| 91精品一卡2卡3卡4卡| 国产精品日韩av在线免费观看| 日本一本二区三区精品| 美女黄网站色视频| 天堂影院成人在线观看| 99热网站在线观看| 久久久久久久久久久丰满| 成人国产麻豆网| 亚洲久久久久久中文字幕| 国产一区二区在线观看日韩| 国产精品蜜桃在线观看| 精华霜和精华液先用哪个| 久热久热在线精品观看| 日本色播在线视频| 91精品国产九色| 五月伊人婷婷丁香| 嘟嘟电影网在线观看| 精华霜和精华液先用哪个| 蜜桃亚洲精品一区二区三区| 人人妻人人澡欧美一区二区| 少妇熟女欧美另类| 国产亚洲5aaaaa淫片| 青春草国产在线视频| 亚洲乱码一区二区免费版| 久久草成人影院| 黄片wwwwww| 精品久久久久久成人av| 亚洲伊人久久精品综合| 亚洲av在线观看美女高潮| 精品人妻偷拍中文字幕| 亚洲欧洲国产日韩| 成人亚洲精品av一区二区| 国产 一区精品| 乱系列少妇在线播放| 亚洲精品乱码久久久久久按摩| 久久精品久久久久久久性| 日韩精品青青久久久久久| 建设人人有责人人尽责人人享有的 | 永久免费av网站大全| 人妻少妇偷人精品九色| 午夜福利高清视频| 三级毛片av免费| 国产极品天堂在线| 国产单亲对白刺激| av免费观看日本| 欧美日韩亚洲高清精品| 亚洲成人精品中文字幕电影| 国产视频内射| 美女大奶头视频| 直男gayav资源| 亚州av有码| 欧美高清成人免费视频www| 国产黄片视频在线免费观看| 人人妻人人澡人人爽人人夜夜 | av国产久精品久网站免费入址| 美女脱内裤让男人舔精品视频| 日日撸夜夜添| 国产日韩欧美在线精品| 精品午夜福利在线看| 夜夜爽夜夜爽视频| 亚洲国产欧美人成| 欧美+日韩+精品| 亚洲av免费高清在线观看| 日本色播在线视频| 免费人成在线观看视频色| 日韩av在线大香蕉| 中文字幕av成人在线电影| av线在线观看网站| 我要看日韩黄色一级片| 美女脱内裤让男人舔精品视频| 天堂俺去俺来也www色官网 | 色尼玛亚洲综合影院| 淫秽高清视频在线观看| 九九在线视频观看精品| 干丝袜人妻中文字幕| 777米奇影视久久| 欧美激情久久久久久爽电影| 国产黄色视频一区二区在线观看| 亚洲国产精品专区欧美| 久久久久免费精品人妻一区二区| 99久久精品国产国产毛片| 亚洲精品日韩av片在线观看| 啦啦啦韩国在线观看视频| 国产成人a∨麻豆精品| 亚洲精品亚洲一区二区| videos熟女内射| 观看美女的网站| 麻豆久久精品国产亚洲av| 国产91av在线免费观看| 91久久精品国产一区二区三区| 九九在线视频观看精品| 免费看日本二区| 亚洲av一区综合| 国产精品久久久久久久电影| 又爽又黄无遮挡网站| 国产亚洲91精品色在线| 伦理电影大哥的女人| 97精品久久久久久久久久精品| 亚洲怡红院男人天堂| 免费看不卡的av| 亚洲欧美清纯卡通| 深夜a级毛片| 久久精品熟女亚洲av麻豆精品 | 精品少妇黑人巨大在线播放| 国产片特级美女逼逼视频| 久久久久精品久久久久真实原创| 中文字幕av在线有码专区| 欧美成人一区二区免费高清观看| 免费看a级黄色片| 国产在线一区二区三区精| 亚洲av福利一区| 免费看不卡的av| 日韩欧美 国产精品| 免费看av在线观看网站| 亚洲综合色惰| 亚洲欧美一区二区三区国产| 亚洲av国产av综合av卡| 韩国高清视频一区二区三区| 久久精品国产自在天天线| 国产色婷婷99| 欧美一区二区亚洲| h日本视频在线播放| 亚洲精品日本国产第一区| 91午夜精品亚洲一区二区三区| 午夜精品国产一区二区电影 | 男的添女的下面高潮视频| 国产乱来视频区| 久久这里只有精品中国| 丝袜喷水一区| 国产高清国产精品国产三级 | 少妇熟女欧美另类|