• <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)變土地利用方式
    丁香六月天网| 国产色婷婷99| 亚洲欧洲日产国产| a级毛片在线看网站| 国产伦理片在线播放av一区| 日日啪夜夜爽| 9热在线视频观看99| 无限看片的www在线观看| svipshipincom国产片| 国产一区二区在线观看av| 99国产综合亚洲精品| kizo精华| 亚洲图色成人| 日韩中文字幕欧美一区二区 | 99香蕉大伊视频| 免费女性裸体啪啪无遮挡网站| 亚洲国产精品国产精品| 免费在线观看视频国产中文字幕亚洲 | 女人被躁到高潮嗷嗷叫费观| 啦啦啦视频在线资源免费观看| 国产欧美日韩综合在线一区二区| 久久女婷五月综合色啪小说| 少妇精品久久久久久久| 国产精品久久久久成人av| 多毛熟女@视频| 91国产中文字幕| 欧美激情高清一区二区三区 | 满18在线观看网站| 欧美黑人精品巨大| 亚洲少妇的诱惑av| 一本色道久久久久久精品综合| 丰满饥渴人妻一区二区三| 久久久久久久久免费视频了| 看非洲黑人一级黄片| 肉色欧美久久久久久久蜜桃| 蜜桃在线观看..| 狠狠精品人妻久久久久久综合| 我要看黄色一级片免费的| www.熟女人妻精品国产| 日韩伦理黄色片| 国产成人精品久久久久久| 亚洲久久久国产精品| 精品一区二区免费观看| 国产成人精品无人区| 国产成人精品无人区| 天天躁夜夜躁狠狠躁躁| 一本久久精品| 各种免费的搞黄视频| 亚洲一区二区三区欧美精品| 人体艺术视频欧美日本| 亚洲,欧美精品.| 三上悠亚av全集在线观看| av在线播放精品| 久久久精品94久久精品| 亚洲美女视频黄频| 日本色播在线视频| 亚洲av欧美aⅴ国产| 亚洲成人国产一区在线观看 | 亚洲欧美激情在线| 一本大道久久a久久精品| 久久热在线av| 国产精品人妻久久久影院| 又大又黄又爽视频免费| 男女免费视频国产| 亚洲人成77777在线视频| 国产在线一区二区三区精| 亚洲自偷自拍图片 自拍| 久久久久精品国产欧美久久久 | 国产亚洲午夜精品一区二区久久| 精品一区二区免费观看| 久久久精品区二区三区| 久久ye,这里只有精品| av免费观看日本| 成人午夜精彩视频在线观看| 亚洲国产精品成人久久小说| 熟女少妇亚洲综合色aaa.| 看非洲黑人一级黄片| 电影成人av| 午夜福利在线免费观看网站| 国产成人午夜福利电影在线观看| 嫩草影视91久久| 韩国av在线不卡| 天天躁日日躁夜夜躁夜夜| 一二三四中文在线观看免费高清| 另类亚洲欧美激情| 涩涩av久久男人的天堂| 韩国精品一区二区三区| 久久久久久久大尺度免费视频| 日本欧美国产在线视频| 国产免费视频播放在线视频| 亚洲免费av在线视频| 久久久精品国产亚洲av高清涩受| 在线观看免费日韩欧美大片| 大香蕉久久网| 不卡视频在线观看欧美| 久久久精品免费免费高清| 国产又爽黄色视频| 亚洲伊人色综图| 亚洲国产欧美日韩在线播放| 欧美亚洲日本最大视频资源| 亚洲精品成人av观看孕妇| 啦啦啦在线观看免费高清www| 最黄视频免费看| av网站免费在线观看视频| 国产在线免费精品| 蜜桃在线观看..| 美女午夜性视频免费| 久久精品国产a三级三级三级| 亚洲熟女毛片儿| 亚洲欧美清纯卡通| 啦啦啦视频在线资源免费观看| 亚洲精品国产色婷婷电影| 少妇的丰满在线观看| 免费人妻精品一区二区三区视频| 亚洲精华国产精华液的使用体验| 日韩视频在线欧美| 新久久久久国产一级毛片| 亚洲伊人色综图| 老司机深夜福利视频在线观看 | 99久久综合免费| 男女边吃奶边做爰视频| 国产精品亚洲av一区麻豆 | 国产片特级美女逼逼视频| 王馨瑶露胸无遮挡在线观看| 人人澡人人妻人| 2018国产大陆天天弄谢| 色婷婷av一区二区三区视频| 亚洲天堂av无毛| 久久久亚洲精品成人影院| 午夜福利在线免费观看网站| 波野结衣二区三区在线| 亚洲久久久国产精品| 亚洲图色成人| 天天躁狠狠躁夜夜躁狠狠躁| 日韩欧美精品免费久久| 最近手机中文字幕大全| 夜夜骑夜夜射夜夜干| 人妻人人澡人人爽人人| 美女国产高潮福利片在线看| 国产精品嫩草影院av在线观看| 又黄又粗又硬又大视频| 国产免费视频播放在线视频| 91精品国产国语对白视频| 啦啦啦中文免费视频观看日本| 在线观看三级黄色| 精品久久蜜臀av无| 美女扒开内裤让男人捅视频| 日韩中文字幕视频在线看片| 精品午夜福利在线看| 精品国产乱码久久久久久小说| 国产极品天堂在线| 母亲3免费完整高清在线观看| 男女免费视频国产| 欧美在线一区亚洲| 51午夜福利影视在线观看| 亚洲av中文av极速乱| 最新的欧美精品一区二区| 日韩大码丰满熟妇| 国产精品一国产av| 亚洲伊人色综图| 亚洲国产成人一精品久久久| 亚洲欧美激情在线| 飞空精品影院首页| 如何舔出高潮| 新久久久久国产一级毛片| 国产爽快片一区二区三区| 老司机靠b影院| 悠悠久久av| 欧美 日韩 精品 国产| 国产精品久久久久久精品古装| 韩国av在线不卡| 秋霞在线观看毛片| 日韩一区二区三区影片| 欧美精品亚洲一区二区| 日韩 亚洲 欧美在线| 国产精品国产三级专区第一集| 狠狠婷婷综合久久久久久88av| 叶爱在线成人免费视频播放| 成人亚洲欧美一区二区av| 成人亚洲欧美一区二区av| 大陆偷拍与自拍| 男女边摸边吃奶| 亚洲av中文av极速乱| 国产精品久久久久久人妻精品电影 | 别揉我奶头~嗯~啊~动态视频 | 一边摸一边抽搐一进一出视频| 一级a爱视频在线免费观看| 激情视频va一区二区三区| 不卡av一区二区三区| 日韩 亚洲 欧美在线| 啦啦啦啦在线视频资源| 久久久久久久国产电影| 三上悠亚av全集在线观看| 纯流量卡能插随身wifi吗| 国产一区有黄有色的免费视频| 久热这里只有精品99| 久久毛片免费看一区二区三区| avwww免费| 国产不卡av网站在线观看| kizo精华| 不卡视频在线观看欧美| 国产一区二区激情短视频 | 久久ye,这里只有精品| 黑人欧美特级aaaaaa片| 最近中文字幕高清免费大全6| 国产成人免费无遮挡视频| 好男人视频免费观看在线| 9191精品国产免费久久| 久久精品亚洲熟妇少妇任你| 国产精品亚洲av一区麻豆 | 少妇人妻精品综合一区二区| 久久久久久久国产电影| 欧美97在线视频| 丝瓜视频免费看黄片| 丝袜喷水一区| 丝袜喷水一区| 亚洲图色成人| 波多野结衣av一区二区av| h视频一区二区三区| 国产伦理片在线播放av一区| 欧美激情极品国产一区二区三区| 夫妻性生交免费视频一级片| 欧美精品一区二区大全| 十八禁人妻一区二区| 99精品久久久久人妻精品| 中文精品一卡2卡3卡4更新| 国产日韩一区二区三区精品不卡| 国产人伦9x9x在线观看| 中文精品一卡2卡3卡4更新| 中文字幕色久视频| 国产一区二区在线观看av| 丰满乱子伦码专区| 亚洲婷婷狠狠爱综合网| 99九九在线精品视频| 97精品久久久久久久久久精品| 自线自在国产av| 美女午夜性视频免费| 啦啦啦 在线观看视频| 满18在线观看网站| 99热国产这里只有精品6| 熟女少妇亚洲综合色aaa.| 国产亚洲av高清不卡| 捣出白浆h1v1| 99久久精品国产亚洲精品| 丰满少妇做爰视频| 国产xxxxx性猛交| 欧美精品人与动牲交sv欧美| 国产日韩一区二区三区精品不卡| 99热全是精品| 久久婷婷青草| 中文天堂在线官网| 久久久久人妻精品一区果冻| 韩国av在线不卡| 欧美av亚洲av综合av国产av | 亚洲伊人久久精品综合| 亚洲精品国产色婷婷电影| 国产精品久久久久成人av| 男人舔女人的私密视频| 18禁国产床啪视频网站| 亚洲中文av在线| 亚洲精品视频女| 日本午夜av视频| 亚洲欧洲日产国产| 丝袜美足系列| 少妇 在线观看| 国产亚洲精品第一综合不卡| 18禁国产床啪视频网站| 狠狠婷婷综合久久久久久88av| 免费观看a级毛片全部| 午夜激情久久久久久久| 久久韩国三级中文字幕| 精品国产超薄肉色丝袜足j| 亚洲国产精品国产精品| 少妇被粗大的猛进出69影院| 美女高潮到喷水免费观看| 成人毛片60女人毛片免费| 国产精品 欧美亚洲| 亚洲自偷自拍图片 自拍| 精品亚洲成a人片在线观看| 最近最新中文字幕大全免费视频 | 满18在线观看网站| 69精品国产乱码久久久| 777米奇影视久久| 日日爽夜夜爽网站| a 毛片基地| 黑人巨大精品欧美一区二区蜜桃| 亚洲精华国产精华液的使用体验| 最近最新中文字幕大全免费视频 | 如日韩欧美国产精品一区二区三区| 国产片内射在线| 夜夜骑夜夜射夜夜干| 亚洲图色成人| 美女高潮到喷水免费观看| 久久久亚洲精品成人影院| 亚洲av福利一区| 日日爽夜夜爽网站| 黄网站色视频无遮挡免费观看| 国产精品香港三级国产av潘金莲 | 热99久久久久精品小说推荐| 人体艺术视频欧美日本| 免费人妻精品一区二区三区视频| 免费久久久久久久精品成人欧美视频| 精品国产国语对白av| 国产精品麻豆人妻色哟哟久久| 国产一区二区三区综合在线观看| 欧美亚洲 丝袜 人妻 在线| 成人国产麻豆网| 综合色丁香网| 黄网站色视频无遮挡免费观看| 赤兔流量卡办理| 国产精品久久久久久人妻精品电影 | 一边摸一边抽搐一进一出视频| 国产成人欧美| 人妻 亚洲 视频| 大陆偷拍与自拍| 成人黄色视频免费在线看| 日本猛色少妇xxxxx猛交久久| 亚洲欧美一区二区三区久久| 热re99久久精品国产66热6| 日韩 欧美 亚洲 中文字幕| 亚洲成色77777| 国产一区二区三区av在线| av网站在线播放免费| 午夜日韩欧美国产| 熟妇人妻不卡中文字幕| 欧美变态另类bdsm刘玥| 大陆偷拍与自拍| 一级毛片 在线播放| 亚洲国产欧美网| 久久久久久人妻| 黑丝袜美女国产一区| 亚洲美女视频黄频| 久久精品久久精品一区二区三区| 97精品久久久久久久久久精品| 日本91视频免费播放| 国产在线视频一区二区| 伊人亚洲综合成人网| 狂野欧美激情性xxxx| 丝袜美足系列| 欧美日韩亚洲高清精品| 免费看不卡的av| 亚洲欧美精品自产自拍| 亚洲精品乱久久久久久| 亚洲国产最新在线播放| 男女午夜视频在线观看| 波野结衣二区三区在线| 国产精品秋霞免费鲁丝片| 亚洲国产最新在线播放| 国产亚洲一区二区精品| 亚洲精品一二三| 亚洲人成电影观看| 国产精品三级大全| 又黄又粗又硬又大视频| 少妇精品久久久久久久| 国产伦人伦偷精品视频| 亚洲av福利一区| 在现免费观看毛片| 男人舔女人的私密视频| 久久久久精品人妻al黑| 国产精品偷伦视频观看了| 中国国产av一级| 另类精品久久| 高清视频免费观看一区二区| 久久精品亚洲av国产电影网| 久久久久久久久免费视频了| 激情视频va一区二区三区| 国产精品久久久人人做人人爽| 久久精品国产亚洲av高清一级| 午夜久久久在线观看| 久久精品亚洲熟妇少妇任你| 在线观看一区二区三区激情| 亚洲第一区二区三区不卡| 免费av中文字幕在线| 黄色视频在线播放观看不卡| 考比视频在线观看| 叶爱在线成人免费视频播放| 亚洲国产av影院在线观看| 国产av国产精品国产| 丝袜在线中文字幕| 欧美日本中文国产一区发布| 国产极品粉嫩免费观看在线| 王馨瑶露胸无遮挡在线观看| 丁香六月天网| 在线看a的网站| 2018国产大陆天天弄谢| 女性被躁到高潮视频| 国产免费又黄又爽又色| 51午夜福利影视在线观看| 一本色道久久久久久精品综合| 欧美乱码精品一区二区三区| 十八禁高潮呻吟视频| 99精国产麻豆久久婷婷| av片东京热男人的天堂| 亚洲国产最新在线播放| 女人高潮潮喷娇喘18禁视频| 18禁动态无遮挡网站| 宅男免费午夜| 欧美亚洲 丝袜 人妻 在线| av天堂久久9| 韩国av在线不卡| 尾随美女入室| 国产一区二区三区av在线| 日本爱情动作片www.在线观看| 欧美成人精品欧美一级黄| 日本黄色日本黄色录像| 亚洲天堂av无毛| 嫩草影院入口| 黑人巨大精品欧美一区二区蜜桃| 亚洲av电影在线进入| 国产成人免费观看mmmm| 曰老女人黄片| 深夜精品福利| 99久久精品国产亚洲精品| 亚洲综合精品二区| 一级爰片在线观看| 国产黄色免费在线视频| 男女国产视频网站| 黄色 视频免费看| 99精国产麻豆久久婷婷| 男女边吃奶边做爰视频| 秋霞在线观看毛片| 国产精品国产av在线观看| 亚洲美女黄色视频免费看| 久久韩国三级中文字幕| av线在线观看网站| 满18在线观看网站| 国产一区二区 视频在线| 如何舔出高潮| 你懂的网址亚洲精品在线观看| 欧美中文综合在线视频| 国产无遮挡羞羞视频在线观看| 亚洲av电影在线进入| 亚洲国产欧美在线一区| 水蜜桃什么品种好| 国产精品二区激情视频| 人人妻人人澡人人看| 我要看黄色一级片免费的| 人妻 亚洲 视频| 如日韩欧美国产精品一区二区三区| 久久国产亚洲av麻豆专区| 在线观看人妻少妇| 天天躁夜夜躁狠狠久久av| av天堂久久9| 久久久久久人人人人人| 人体艺术视频欧美日本| 美女扒开内裤让男人捅视频| 青春草国产在线视频| 美女高潮到喷水免费观看| 色婷婷av一区二区三区视频| 中文乱码字字幕精品一区二区三区| 人妻 亚洲 视频| 两性夫妻黄色片| 99九九在线精品视频| 精品人妻一区二区三区麻豆| 人成视频在线观看免费观看| 青青草视频在线视频观看| 亚洲国产欧美网| 赤兔流量卡办理| 性少妇av在线| 水蜜桃什么品种好| 亚洲精品美女久久av网站| 18禁裸乳无遮挡动漫免费视频| 欧美乱码精品一区二区三区| 亚洲人成网站在线观看播放| 只有这里有精品99| 成人18禁高潮啪啪吃奶动态图| 亚洲专区中文字幕在线 | 午夜久久久在线观看| 亚洲国产欧美日韩在线播放| 成人亚洲欧美一区二区av| www.精华液| 女人被躁到高潮嗷嗷叫费观| 香蕉丝袜av| 青春草亚洲视频在线观看| 日日啪夜夜爽| 国产黄色免费在线视频| 精品卡一卡二卡四卡免费| 亚洲第一青青草原| 最新在线观看一区二区三区 | 狂野欧美激情性xxxx| 国产亚洲av高清不卡| 黄频高清免费视频| 久久综合国产亚洲精品| 国产亚洲精品第一综合不卡| 欧美亚洲 丝袜 人妻 在线| 国产av国产精品国产| 男女无遮挡免费网站观看| 99热国产这里只有精品6| 国产精品99久久99久久久不卡 | 国产男女内射视频| 久久久久久人妻| 亚洲精品久久久久久婷婷小说| 丝袜美足系列| 久久天堂一区二区三区四区| e午夜精品久久久久久久| 性色av一级| 精品免费久久久久久久清纯 | 99re6热这里在线精品视频| 一级毛片 在线播放| 国产日韩欧美视频二区| 色网站视频免费| 亚洲情色 制服丝袜| 天天影视国产精品| 看非洲黑人一级黄片| 丝袜脚勾引网站| 成人免费观看视频高清| 国产精品三级大全| www.自偷自拍.com| 欧美日韩成人在线一区二区| 午夜福利影视在线免费观看| 中文字幕人妻丝袜一区二区 | 在线观看www视频免费| 精品免费久久久久久久清纯 | 国产精品国产av在线观看| 美女主播在线视频| 欧美日韩精品网址| 国产99久久九九免费精品| 新久久久久国产一级毛片| 中文欧美无线码| 亚洲四区av| 欧美黑人欧美精品刺激| 色综合欧美亚洲国产小说| www.自偷自拍.com| 波多野结衣一区麻豆| 777久久人妻少妇嫩草av网站| 亚洲成人av在线免费| 黄色怎么调成土黄色| 国产一卡二卡三卡精品 | 七月丁香在线播放| 悠悠久久av| 麻豆av在线久日| 成年人午夜在线观看视频| 国产成人精品在线电影| 制服人妻中文乱码| h视频一区二区三区| 国产成人91sexporn| 建设人人有责人人尽责人人享有的| 久久99一区二区三区| 免费观看性生交大片5| 另类亚洲欧美激情| 大片电影免费在线观看免费| 精品亚洲乱码少妇综合久久| 哪个播放器可以免费观看大片| 99久国产av精品国产电影| 欧美黑人欧美精品刺激| 美女扒开内裤让男人捅视频| 国产野战对白在线观看| 中文字幕人妻熟女乱码| 日韩 亚洲 欧美在线| 欧美xxⅹ黑人| 午夜免费鲁丝| 尾随美女入室| 亚洲欧洲精品一区二区精品久久久 | 免费久久久久久久精品成人欧美视频| 纯流量卡能插随身wifi吗| 亚洲人成电影观看| 97人妻天天添夜夜摸| 欧美xxⅹ黑人| 久久精品国产综合久久久| 99re6热这里在线精品视频| 欧美另类一区| 国产精品香港三级国产av潘金莲 | 91国产中文字幕| 最近2019中文字幕mv第一页| 欧美日韩一区二区视频在线观看视频在线| av在线播放精品| 免费在线观看完整版高清| 男女边吃奶边做爰视频| 国产又爽黄色视频| 一级毛片 在线播放| 韩国精品一区二区三区| 成人黄色视频免费在线看| 麻豆乱淫一区二区| 伊人亚洲综合成人网| www日本在线高清视频| 老司机影院毛片| 美女福利国产在线| 欧美激情极品国产一区二区三区| 黄片播放在线免费| 晚上一个人看的免费电影| 免费观看a级毛片全部| 久久久久久免费高清国产稀缺| 亚洲精品视频女| 日韩中文字幕欧美一区二区 | 天堂中文最新版在线下载| 中文字幕av电影在线播放| 美女高潮到喷水免费观看| 91精品伊人久久大香线蕉| 精品人妻熟女毛片av久久网站| 亚洲情色 制服丝袜| 丰满迷人的少妇在线观看| 少妇被粗大的猛进出69影院| 亚洲综合色网址| 美女视频免费永久观看网站| 91国产中文字幕| 激情视频va一区二区三区| 亚洲精品,欧美精品| 建设人人有责人人尽责人人享有的| 亚洲精品日韩在线中文字幕| 少妇人妻精品综合一区二区| 久久天躁狠狠躁夜夜2o2o | 久久 成人 亚洲| 九草在线视频观看| 啦啦啦视频在线资源免费观看| 国产精品.久久久| 丝袜美足系列| 黑丝袜美女国产一区| 天天躁日日躁夜夜躁夜夜| 欧美激情高清一区二区三区 | 国产一区二区激情短视频 | 亚洲中文av在线| 国产成人a∨麻豆精品| 黑丝袜美女国产一区| 精品久久久精品久久久| 欧美激情高清一区二区三区 |