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

    2001—2009年中國碳排放與碳足跡時空格局

    2014-08-10 12:26:40吳文佳蔣金亮高全洲蔣海兵
    生態(tài)學(xué)報(bào) 2014年22期
    關(guān)鍵詞:赤字足跡空間

    吳文佳,蔣金亮,高全洲,蔣海兵

    (1. 北京大學(xué)城市與環(huán)境學(xué)院,北京大學(xué)地表過程分析與模擬教育部重點(diǎn)實(shí)驗(yàn)室,北京 100871;2. 南京大學(xué)地理與海洋科學(xué)學(xué)院,南京 210093;3. 中山大學(xué)地理科學(xué)與規(guī)劃學(xué)院, 廣東省城市化與地理環(huán)境空間模擬重點(diǎn)實(shí)驗(yàn)室,廣州 510275;4. 鹽城師范學(xué)院城市與資源環(huán)境學(xué)院,鹽城 224051)

    2001—2009年中國碳排放與碳足跡時空格局

    吳文佳1,蔣金亮2,高全洲3,蔣海兵4,*

    (1. 北京大學(xué)城市與環(huán)境學(xué)院,北京大學(xué)地表過程分析與模擬教育部重點(diǎn)實(shí)驗(yàn)室,北京 100871;2. 南京大學(xué)地理與海洋科學(xué)學(xué)院,南京 210093;3. 中山大學(xué)地理科學(xué)與規(guī)劃學(xué)院, 廣東省城市化與地理環(huán)境空間模擬重點(diǎn)實(shí)驗(yàn)室,廣州 510275;4. 鹽城師范學(xué)院城市與資源環(huán)境學(xué)院,鹽城 224051)

    碳排放引發(fā)的全球變暖給自然環(huán)境及人類社會都帶來了顯著影響,而碳足跡可以衡量自然生態(tài)系統(tǒng)對人類活動碳排放的響應(yīng)。為研究自然-社會二元系統(tǒng)碳動態(tài),基于MODIS(Moderate Resolution Imaging Spectroradiometer)數(shù)據(jù)和統(tǒng)計(jì)資料計(jì)算2001—2009年中國陸地植被凈初級生產(chǎn)力、能源消費(fèi)碳排放、碳足跡和碳赤字;在GIS(Geographic Information System)技術(shù)支持下,運(yùn)用空間自相關(guān)分析方法討論其時空格局;據(jù)此劃分生態(tài)經(jīng)濟(jì)區(qū)。結(jié)果表明:(1)2001—2009年全國植被凈初級生產(chǎn)力(Net Primary Production,NPP)平均值為3.32 Pg C/a(1 Pg = 1015g),呈西南地區(qū)>東南沿海>華中、華東地區(qū)>東北、華北地區(qū)>西北地區(qū)的空間格局;(2)2001—2009年全國能源消費(fèi)碳排放逐年增加,年均增長率16.7%,多年平均值2.53 Pg C/a,呈東部>中部>西部的空間格局;(3)2001—2009年全國碳足跡逐年增加,年均增長率14.7%,多年平均值6.98×106km2;具有正碳赤字(即碳源)的省份為山西、環(huán)渤海地區(qū)各省、長三角地區(qū)各省、廣東;相鄰省份碳赤字的相對大小由于互相影響而改變;(4)全國分為中東部、南部、北部、西部四個生態(tài)經(jīng)濟(jì)大區(qū)。研究結(jié)果直觀揭示了中國碳排放和碳足跡的時空動態(tài),為實(shí)現(xiàn)自然-社會二元系統(tǒng)的可持續(xù)發(fā)展提供科學(xué)依據(jù)。

    自然-社會二元系統(tǒng);MODIS;碳排放;碳足跡;空間自相關(guān);生態(tài)經(jīng)濟(jì)區(qū)劃

    21世紀(jì)以來,全球氣候變化問題受到廣泛關(guān)注,碳排放被認(rèn)為是引起全球變暖的主導(dǎo)因素[1- 2]。一方面,陸地生態(tài)系統(tǒng)具有碳源/匯雙重作用,對全球碳平衡貢獻(xiàn)巨大,了解不同生態(tài)系統(tǒng)內(nèi)部植被、土壤碳動態(tài)是碳循環(huán)研究的重要組成部分[3- 7];另一方面,隨著2005年《京都議定書》的生效和2009年哥本哈根全球氣候大會的召開,削減化石燃料碳排放成為全球共識,社會經(jīng)濟(jì)系統(tǒng)中的人為碳排放研究漸成熱點(diǎn)[8- 10]。然而,基于以上兩方面開展自然-社會二元系統(tǒng)分析才能更好地了解自然和人為因素綜合作用下的碳動態(tài)以及自然生態(tài)系統(tǒng)對社會經(jīng)濟(jì)活動的反饋能力[11]。

    由生態(tài)足跡引申而來的碳足跡概念能夠直觀衡量自然系統(tǒng)對人類活動碳排放的響應(yīng)[12],自提出以來已在家庭住戶[13]、大都市區(qū)[14]的碳足跡估算和模型比較[15]等方面取得了一定進(jìn)展;國內(nèi)學(xué)者對于不同土地利用類型[16]、不同產(chǎn)業(yè)[17- 18]的碳足跡以及全國尺度上的碳足跡時空格局[19- 20]也做出了探索。這些研究為開展碳動態(tài)的進(jìn)一步研究奠定了基礎(chǔ),但目前對于中國自然-社會二元系統(tǒng)的碳排放和碳足跡格局分析仍處于起步階段,主要薄弱點(diǎn)在于:一是時空格局僅限于定性描述,尚未應(yīng)用空間統(tǒng)計(jì)分析方法[19];二是土地利用變化和植被碳匯多來源于統(tǒng)計(jì)數(shù)據(jù)和經(jīng)驗(yàn)數(shù)據(jù)[16,21- 22],口徑不一且難以體現(xiàn)異質(zhì)性。本文首次采用分辨率高、來源穩(wěn)定、實(shí)時性強(qiáng)的MODIS(Moderate Resolution Imaging Spectroradiometer)數(shù)據(jù),引入空間自相關(guān)方法,評估2001—2009年中國碳排放和碳足跡時空動態(tài),并把研究結(jié)果應(yīng)用到生態(tài)經(jīng)濟(jì)區(qū)劃中,旨在構(gòu)建全國尺度上碳動態(tài)研究的空間統(tǒng)計(jì)方法,并為優(yōu)化各區(qū)域的生態(tài)效益和經(jīng)濟(jì)效益提供科學(xué)參考。

    1 數(shù)據(jù)與方法

    1.1 研究區(qū)域

    本文參考1961年國家計(jì)委設(shè)置的華北、東北、華東、中南、西南、西北六大經(jīng)濟(jì)協(xié)作區(qū)方案,加上港、澳、臺,將全國分為七大區(qū)域。其中,藏、港、澳、臺的能源數(shù)據(jù)缺失,不作討論。

    1.2 數(shù)據(jù)及預(yù)處理

    1.2.1 數(shù)據(jù)來源及精度

    (1)遙感數(shù)據(jù)

    本文土地利用和植被凈初級生產(chǎn)力(Net Primary Production,NPP)數(shù)據(jù)分別來自NASA(National Aeronautics and Space Administration,美國國家航空航天局)網(wǎng)站(https://lpdaac.usgs.gov/products/modis_products_table)提供的MODIS產(chǎn)品MCD12Q1和MOD17A3。MCD12Q1空間分辨率為500 m,采用監(jiān)督?jīng)Q策樹分類的方法處理年度的Terra和Aqua衛(wèi)星觀測所得數(shù)據(jù);本文僅需據(jù)其提取土地總面積信息。MOD17A3空間分辨率為1 km,時間分辨率為1a,由參考BIOME-BGC(Biome-Biogeochemical Cycles)模型與光能利用率模型建立的NPP估算模型處理Terra衛(wèi)星數(shù)據(jù)得到。較之傳統(tǒng)的回歸模型,BIOME-BGC模型是從植被機(jī)理出發(fā)建立的過程模型,使用了更多的參數(shù)和更精密的估算方法以提高精度[23- 24]。近年來,MOD17A3已對全球多種植被類型進(jìn)行了驗(yàn)證[25- 27],在我國也已被用來監(jiān)測青藏高原地區(qū)[28]、南水北調(diào)東線地區(qū)[29]、陜西省[30]的植被生長狀況以及廣東省水稻生物質(zhì)能[23]、東北地區(qū)植被農(nóng)田的時空特征[30- 31]等;其精度能夠滿足全國尺度上的NPP分析。

    (2)統(tǒng)計(jì)數(shù)據(jù)

    目前以化石能源為代表的傳統(tǒng)能源是產(chǎn)生碳排放的主導(dǎo)因素[16],故本文選取煤炭、原油、汽油、煤油、柴油、燃料油、天然氣、電力8類能源為社會經(jīng)濟(jì)系統(tǒng)碳排放的核算指標(biāo)。2001—2009年各省各類能源消費(fèi)數(shù)據(jù)來自《中國能源統(tǒng)計(jì)年鑒》。

    (3)輔助數(shù)據(jù)

    主要來自國家基礎(chǔ)地理信息中心提供的《1∶400萬中國省級政區(qū)圖》。

    1.2.2 數(shù)據(jù)預(yù)處理

    為構(gòu)建與研究區(qū)實(shí)際地理基礎(chǔ)一致的MODIS時間序列數(shù)據(jù)集,使用MRT(MODIS Reprojection Tool)批處理軟件和ArcGIS9.3提取2001—2009年土地利用和NPP數(shù)據(jù)集。MCD12Q1產(chǎn)品包含5個土地覆蓋分類方案不同的數(shù)據(jù)集,本文根據(jù)Land Cover Type1數(shù)據(jù)集,即IGBP(International Geosphere-Biosphere Program)分類方案計(jì)算各省土地總面積;MOD17A3產(chǎn)品包含3個數(shù)據(jù)集,本文提取其Gridded 1 km Annual Net Productivity數(shù)據(jù)集獲取各省逐年NPP數(shù)據(jù)(圖1)。

    圖1 數(shù)據(jù)預(yù)處理流程圖Fig.1 The flow diagram of data preprocessing

    1.3 計(jì)量模型

    1.3.1 自然生態(tài)系統(tǒng)碳動態(tài)

    隨著植物生長及其生物量增加,CO2通過植被光合作用進(jìn)入自然生態(tài)系統(tǒng)并儲存起來。本文以NPP作為自然生態(tài)系統(tǒng)碳動態(tài)的衡量指標(biāo)。NPP是指綠色植物在單位時間單位面積上扣除自氧呼吸后的有機(jī)干物質(zhì)生產(chǎn)量[32],反映了植物固定和轉(zhuǎn)化光合產(chǎn)物的效率,是判定生態(tài)系統(tǒng)碳匯、調(diào)節(jié)生態(tài)過程的主要因子和表示碳吸收的常用指標(biāo)[33],已在生態(tài)系統(tǒng)碳循環(huán)研究中得到廣泛應(yīng)用[34- 35]。

    1.3.2 社會經(jīng)濟(jì)系統(tǒng)碳排放

    構(gòu)建基于能源消費(fèi)的社會經(jīng)濟(jì)系統(tǒng)碳排放計(jì)量模型(式1)。

    (1)

    式中,CEit為i區(qū)域t年能源消費(fèi)碳排放總量(Pg C/a),Eijt為i區(qū)域t年第j種能源消費(fèi)總量(kg或kW·h),cvj為第j種能源折算標(biāo)準(zhǔn)煤的參考系數(shù)(表1),δj為第j種能源的碳排放系數(shù),取自相關(guān)文獻(xiàn)平均值(表2)。

    表1 各類能源折算標(biāo)準(zhǔn)煤的參考系數(shù)

    表2 碳排放轉(zhuǎn)換系數(shù)表

    1.3.3 碳足跡

    根據(jù)環(huán)球足跡網(wǎng)絡(luò)(Global Financial Network,GFN)的定義,碳足跡是通過光合作用吸收由化石燃料燃燒所產(chǎn)生的CO2所需要的生產(chǎn)性土地(植被)面積[42],用以測度人類能源消費(fèi)碳排放對生態(tài)空間的占用情況。當(dāng)區(qū)域生態(tài)容量不足以吸收由化石燃料燃燒所產(chǎn)生的CO2時產(chǎn)生碳赤字(CFit> 0),區(qū)域發(fā)展處于不可持續(xù)狀態(tài);反之為碳盈余。本文基于對自然生態(tài)系統(tǒng)和社會經(jīng)濟(jì)系統(tǒng)碳動態(tài)的分別測算,參考生態(tài)足跡的計(jì)算方法[9]構(gòu)建碳足跡計(jì)量模型(式2,式3):

    CFit=(CEit/NPPr)×109=CEit/(NPPit/Sit×109)×109=(CEit/NPPit)×Sit

    (2)

    ΔCFit=CFit-Sit

    (3)

    式中,NPPr為區(qū)域綜合NPP(g C m-2a-1),CFit為i區(qū)域t年的碳足跡(km2),CEit為i區(qū)域t年能源消費(fèi)碳排放總量(Pg C/a),Sit為i區(qū)域t年的土地總面積(km2),NPPit為i區(qū)域t年的總凈初級生產(chǎn)力(Pg C/a)。

    1.4 空間統(tǒng)計(jì)分析方法

    空間統(tǒng)計(jì)分析方法充分利用屬性數(shù)據(jù)的空間特性,彌補(bǔ)傳統(tǒng)計(jì)量模型在空間分析和數(shù)值相關(guān)分析上顧此失彼的缺陷。在GIS(Geographic Information System)環(huán)境下對碳動態(tài)進(jìn)行空間統(tǒng)計(jì)分析,可以更形象深刻地揭示其時空格局。本文運(yùn)用了空間統(tǒng)計(jì)分析方法中的空間自相關(guān)分析模型。根據(jù)地理學(xué)第一定律,任何事物都相關(guān),但相近的事物關(guān)聯(lián)更緊密??臻g自相關(guān)分析就是基于這一定律,假設(shè)研究區(qū)中所有值都是非獨(dú)立的,相互之間存在相關(guān)性。正空間自相關(guān)是指鄰近區(qū)域的屬性具有相似的趨勢和取值,負(fù)空間自相關(guān)是指鄰近區(qū)域的屬性具有相反的趨勢和取值[43]。

    1.4.1 空間權(quán)重矩陣

    空間權(quán)重矩陣用以定量界定區(qū)域單元之間的鄰接關(guān)系,本文采用簡單的二進(jìn)制鄰接矩陣(式4),在Geoda095i軟件中以queen標(biāo)準(zhǔn)創(chuàng)建鄰接性權(quán)重文件。

    (4)

    鑒于海南和廣東、福建和中國臺灣聯(lián)系密切,故對生成的權(quán)重文件進(jìn)行微調(diào),改成海南和廣東相鄰接、福建和中國臺灣相鄰接。

    1.4.2 全局空間自相關(guān)

    全局空間自相關(guān)用以判斷現(xiàn)象整體空間聚集性的存在與否及聚集程度的大小。Moran′sI是度量全局空間自相關(guān)的常用指標(biāo)之一:

    (5)

    (6)

    (7)

    (8)

    構(gòu)造Z統(tǒng)計(jì)量對空間自相關(guān)進(jìn)行顯著性檢驗(yàn)(式9):

    (9)

    式中,E(I)、Var(I)分別為I的期望和方差。一般當(dāng)|Z|>1.96時,差異顯著,即在95%的概率下存在空間自相關(guān)。在Geoda095i軟件中可以計(jì)算Moran′sI值,評估全局空間自相關(guān)性。

    1.4.3 局域空間自相關(guān)

    局域空間自相關(guān)能夠反映局部的空間差異程度,推算出高、低值的聚集范圍。LISA(Local Indicators of Spatial Association)是分析局域空間自相關(guān)的常用方法:

    (10)

    各字母含義及檢驗(yàn)統(tǒng)計(jì)量同上。

    在Geoda095i軟件中繪制Moran散點(diǎn)圖,其橫坐標(biāo)為各區(qū)域單元標(biāo)準(zhǔn)化處理后的屬性值,縱坐標(biāo)為標(biāo)準(zhǔn)化處理后的空間鄰接矩陣所決定的相鄰單元的屬性值的平均值。散點(diǎn)圖的4個象限表達(dá)了某一區(qū)域和其周圍區(qū)域的局域空間聯(lián)系,一、三象限為正空間自相關(guān),二、四象限為負(fù)空間自相關(guān)。

    2 結(jié)果與分析

    2.1 自然生態(tài)系統(tǒng)碳動態(tài)的時空格局

    2001—2009年全國NPP年平均值變動范圍在341—371 g C m-2a-1之間,多年平均值為362 g C m-2a-1;全國年總NPP值在3.13—3.40 Pg C/a之間,多年平均值約為3.32 Pg C/a,總體波動不大(圖2)。

    研究時段內(nèi)全國多年平均NPP值的空間分布呈西南地區(qū)>東南沿海>華中、華東地區(qū)>東北、華北地區(qū)>西北地區(qū)的格局(圖3);其中海南最大(905 g C m-2a-1),新疆最小(143 g C m-2a-1)。NPP自南向北隨緯度增加而遞減,淮河以北地區(qū)自東向西隨經(jīng)度減小而遞減,淮河以南地區(qū)隨經(jīng)度遞變趨勢不明顯。這種空間格局的形成主要與地帶性植被類型有關(guān),而水熱條件是植被生長的決定因素。自南向北隨溫度降低依次分布熱帶闊葉林、亞熱帶常綠闊葉落葉林、亞熱帶熱帶針葉林、溫帶落葉闊葉林,NPP依次減小。北部地區(qū)自東向西隨水分減小依次分布濕潤溫帶森林、內(nèi)蒙古草原、干草原,NPP逐漸減小。另外,城市化進(jìn)程增加了不透水表面比例,減少了自然植被覆蓋,因而氣候特征相當(dāng)?shù)牡貐^(qū)NPP也會因城市化程度不同而有所差異,比如北京、天津的NPP值小于河北。

    2.2 社會經(jīng)濟(jì)系統(tǒng)碳排放的時空格局

    2001—2009年全國年總能源消費(fèi)碳排放量呈逐年增加趨勢(圖2),由1.52 Pg C/a(2001年)增加到3.53 Pg C/a(2009年),年均增長率為16.6%,多年平均值為2.53 Pg C/a。碳排放總量前5a增速大于后5a,說明研究時段內(nèi)就全國范圍而言碳排放增加趨勢已開始變緩。

    研究時段內(nèi)全國平均單位面積能源消費(fèi)碳排放量呈現(xiàn)東部>中部>西部的空間分異(圖3);其中上海最大(1.17×104g C m-2a-1),青海最小(14 g C m-2a-1)。能源消費(fèi)碳排放與當(dāng)?shù)亟?jīng)濟(jì)發(fā)展需求密切相關(guān),以上海為中心的長三角經(jīng)濟(jì)圈和以北京、天津?yàn)橹行牡沫h(huán)渤海經(jīng)濟(jì)圈單位面積碳排放量最高;新疆、青海、甘肅等發(fā)展較慢的西北省份單位面積碳排放量最低。

    圖2 2001—2009年自然-社會二元系統(tǒng)碳收支時間動態(tài)Fig.2 Temporal patterns of the “nature-society” dualistic carbon fluxes in China from 2001 to 2009

    圖3 2001—2009年自然-社會二元系統(tǒng)平均碳收支空間格局Fig.3 Spatial patterns of the “nature-society” dualistic carbon fluxes in China during 2001—2009

    2.3 碳足跡、碳赤字的時空格局

    2.3.1 碳足跡、碳赤字的時間格局

    2001—2009年全國碳足跡、碳赤字均呈逐年增加趨勢(圖4),碳足跡由4.46×106km2(2001年)增加到9.69×106km2(2009年),年均增長率為14.7%,多年平均值為6.98×106km;碳赤字由-4.73×106km2(2001年)增加到5.01×105km2(2009年),多年平均值為-2.20×106km2,自2009年開始出現(xiàn)正碳赤字。碳足跡、碳赤字后5a增速小于前5a,說明許多地區(qū)已開始采取碳增匯或碳減排措施控制自然-社會二元系統(tǒng)的碳排放增加。

    2.3.2 碳足跡、碳赤字的空間格局

    2001—2009年全國平均碳赤字(圖5)以產(chǎn)煤大省山西省最大(5.03×105km2),其次為重化工-資本密集型的環(huán)渤海地區(qū)、高科技-知識密集型的長三角地區(qū)、輕紡-勞動密集型的珠三角地區(qū);而中西部地區(qū)以及東部的福建、江西、吉林、黑龍江等省具有碳盈余,新疆碳盈余最大(1.16×106km2)。2001年,具有正碳赤字的省份集中在山西、河北、遼寧、山東、河南、江蘇,中西部地區(qū)存在大量碳盈余;2005年,碳赤字在環(huán)渤海經(jīng)濟(jì)圈、長三角經(jīng)濟(jì)圈、珠三角經(jīng)濟(jì)圈明顯加劇;2009年,環(huán)渤海地區(qū)各省碳赤字進(jìn)一步增加,中西部地區(qū)的吉林、重慶、陜西也開始出現(xiàn)正碳赤字。研究時段內(nèi)具有碳赤字的省份逐年增多,說明越來越多的省份生態(tài)用地?zé)o法滿足消納碳排放的需求;總體來看,華東、華北、中南地區(qū)的碳赤字增加較快,說明這些地區(qū)生態(tài)用地被占用帶來的碳匯減少快于能源結(jié)構(gòu)調(diào)整帶來的碳源減少。

    圖4 2001—2009年全國碳足跡、碳赤字時間分布圖Fig.4 Temporal patterns of carbon footprint and carbon deficit in China from 2001 to 2009

    圖5 2001—2009年全國碳赤字空間分布圖Fig.5 Spatial patterns of carbon footprint and carbon deficit in China during 2001—2009

    2.3.3 碳赤字的空間自相關(guān)分析

    (1)Moran指數(shù)

    從全國范圍來看,各年份Moran指數(shù)的統(tǒng)計(jì)量相當(dāng)顯著(圖6),表明碳赤字在全國范圍內(nèi)存在顯著的正空間自相關(guān),存在相似值的空間集聚。從不同區(qū)域來看,西北、華東地區(qū)碳赤字的Moran指數(shù)為正且較顯著,其中經(jīng)濟(jì)發(fā)展迅速而生態(tài)用地緊缺的華東地區(qū)表現(xiàn)為明顯的高值集聚,經(jīng)濟(jì)發(fā)展較慢而生態(tài)用地充裕的西北地區(qū)則為低值集聚。華北地區(qū)碳赤字的Moran指數(shù)為負(fù)且較顯著,存在負(fù)空間自相關(guān),碳赤字表現(xiàn)為明顯的高低間隔分布,區(qū)內(nèi)各省經(jīng)濟(jì)效益和生態(tài)效益參差不齊。中南、東北、西南地區(qū)的Moran指數(shù)趨于0,空間自相關(guān)性很弱,碳赤字大致呈隨機(jī)分布。

    圖6 2001—2009年全國及各區(qū)域碳赤字的Moran指數(shù)時間分布圖Fig.6 Temporal patterns of Moran′s I in China and different regions during 2001—2009

    (2)Moran散點(diǎn)圖和LISA分析

    在Geoda095i軟件中繪制主要年份全國碳赤字的Moran散點(diǎn)圖(圖7)。圖中4個象限表達(dá)了某一區(qū)域和其周圍區(qū)域的局域空間聯(lián)系。第一象限HH(High,High)表示該區(qū)域和周圍屬性值均較高;第二象限LH(Low,High)表示該區(qū)域?qū)傩灾得黠@低于周圍;第三象限LL(Low,Low)表示該區(qū)域和周圍屬性值均較低;第四象限HL(High,Low)表示該區(qū)域?qū)傩灾得黠@高于周圍。觀察各象限包含省區(qū)的變化,可以看出2001—2009年碳赤字空間格局的演變過程包括如下幾類:(1)HH→LH,即碳赤字相對于周邊減小,如海南自1999年提出建設(shè)生態(tài)省戰(zhàn)略以來,生態(tài)旅游、綠色農(nóng)業(yè)等產(chǎn)業(yè)逐漸興起,生態(tài)效益和經(jīng)濟(jì)效益得到了綜合提升,碳赤字比周邊省份減小迅速[44];(2)HH→HL,即碳赤字相對于周邊增加,如廣東,近10年來珠三角地區(qū)的快速密集城市化導(dǎo)致該省自然資源壓力明顯高于周邊[45];(3)LH→LL,如湖南,開始時相對于周邊省份具有較低的碳赤字,隨著省際聯(lián)系的加強(qiáng),在生態(tài)保護(hù)方面對周邊省份起到了示范和帶動作用,使得周邊省份的碳赤字也減小,積極的生態(tài)效應(yīng)向周邊擴(kuò)散;(4)HL→LL,如重慶,開始時相對周邊省份具有較高的碳赤字,隨著整個西部地區(qū)生態(tài)用地的合理規(guī)劃,周邊各省積極的生態(tài)效應(yīng)向該省擴(kuò)散;(5)HL→HH,如山西、遼寧、吉林、陜西,開始時相對于周邊省份具有較高的碳赤字,隨著跨省經(jīng)濟(jì)協(xié)作的加強(qiáng),能源消費(fèi)碳排放在更大范圍內(nèi)迅速增加,消極的生態(tài)效應(yīng)向周邊擴(kuò)散。

    2.4 基于碳動態(tài)的中國生態(tài)經(jīng)濟(jì)區(qū)劃

    生態(tài)經(jīng)濟(jì)區(qū)劃是以自然-社會二元復(fù)合系統(tǒng)為對象,以生態(tài)-經(jīng)濟(jì)耦合理論為基礎(chǔ),根據(jù)自然生態(tài)因素和社會經(jīng)濟(jì)因素的特點(diǎn)及其內(nèi)在聯(lián)系所構(gòu)成的空間組合形式的相似性和差異性,總結(jié)生態(tài)經(jīng)濟(jì)功能的地域分異規(guī)律,劃分融合生態(tài)和經(jīng)濟(jì)要素的地域單元,并明確不同單元在區(qū)域開發(fā)中承擔(dān)的主要功能[46]。與2001年傅伯杰等提出的劃分3個大區(qū)(東部濕潤、半濕潤生態(tài)大區(qū),西北干旱、半干旱生態(tài)大區(qū),青藏高原高寒生態(tài)大區(qū))的中國生態(tài)區(qū)劃方案[47]以及2012年謝高地等提出的劃分4個生態(tài)大區(qū)(東北部濕潤半濕潤生態(tài)大區(qū),北部干旱半干旱生態(tài)大區(qū),南部濕潤生態(tài)大區(qū)和青藏高原生態(tài)大區(qū))的中國生態(tài)區(qū)劃方案[48]不同,本文不只依據(jù)氣候、地貌、植被特征,而是兼顧了自然和人為因素,根據(jù)2001—2009年全國碳赤字的空間格局(圖5)和局域空間自相關(guān)性(圖7),把碳赤字時空動態(tài)較一致的省區(qū)合并,劃分4個大區(qū)(圖8)。

    (1)中、東部地區(qū)

    位于長江以北、太行山以東,包括北京、天津、河北、上海、江蘇、浙江、安徽山東、河南、湖北。該區(qū)跨越中溫帶、暖溫帶、北亞熱帶3個氣候帶,水熱資源沒有長江以南地區(qū)優(yōu)越,NPP較小。長三角、環(huán)渤海兩大經(jīng)濟(jì)圈以及依托“中部崛起”戰(zhàn)略而迅速發(fā)展起來的中原城市群經(jīng)濟(jì)發(fā)展速度快,急劇增加的能源消費(fèi)催生了大量碳排放,碳赤字嚴(yán)重。研究時段內(nèi),區(qū)內(nèi)各省的碳赤字表現(xiàn)出明顯的高值集聚。為了在發(fā)展經(jīng)濟(jì)的同時兼顧生態(tài)效益,該區(qū)未來應(yīng)加強(qiáng)省際協(xié)作,調(diào)整能源結(jié)構(gòu),共同打造低碳發(fā)展模式;并保護(hù)好有限的生態(tài)用地,改善植被覆蓋狀況,延緩碳赤字的進(jìn)一步增大。

    圖7 2001—2009年全國碳赤字Moran散點(diǎn)圖及對應(yīng)省份分布Fig.7 Moran scatterplot of carbon deficit and provinces corresponding to it during 2001—2009散點(diǎn)圖橫坐標(biāo)為相應(yīng)年份碳赤字標(biāo)準(zhǔn)化值,縱坐標(biāo)為相鄰單元碳赤字標(biāo)準(zhǔn)化均值;象限圖四象限分別表示“高-高(HH)”、“低-高(LH)”、“低-低(LL)”、“高-低(HL)”

    (2)南部地區(qū)

    位于長江以南、雪峰山以東,包括福建、江西、湖南、廣東、海南、中國香港、中國澳門、中國臺灣。該區(qū)跨越中亞熱帶、南亞熱帶、熱帶三個氣候帶,水熱資源豐富,NPP大。研究時段內(nèi),區(qū)內(nèi)各省碳赤字的空間集聚不明顯,湖南、海南等省生態(tài)效益較高,廣東等省生態(tài)效益較低。未來應(yīng)充分利用優(yōu)良的自然條件,合理處理經(jīng)濟(jì)發(fā)展和生態(tài)保護(hù)之間的關(guān)系;并強(qiáng)化省際聯(lián)系,尤其是碳赤字較高和較低的省份之間的互補(bǔ),促進(jìn)低碳發(fā)展模式擴(kuò)散和資源合理利用。

    圖8 基于碳赤字空間自相關(guān)性的中國生態(tài)經(jīng)濟(jì)區(qū)劃Fig.8 Eco-economy regionalization in China based on spatial autocorrelation of carbon deficit

    (3)北部地區(qū)

    位于祁連山以東,包括山西、內(nèi)蒙古、遼寧、吉林、黑龍江、陜西、寧夏。該區(qū)域跨越寒溫帶、中溫帶兩個氣候帶,水熱條件較差,NPP較低。未來應(yīng)加強(qiáng)省際合作和互補(bǔ),充分利用內(nèi)蒙古、黑龍江兩省豐富的草地、林地資源,緩解山西、遼寧等經(jīng)濟(jì)發(fā)展較快省份的碳赤字,發(fā)揮不同省份的比較優(yōu)勢和帶動作用。

    (4)西部地區(qū)

    位于賀蘭山—六盤山—大巴山一線以西,包括重慶、四川、貴州、云南、西藏、甘肅、青海、新疆。氣候帶跨度較大,包括溫帶、亞熱帶、高原氣候區(qū)。除西南小部分地區(qū)外,NPP較小,但生態(tài)用地相當(dāng)豐富。目前該區(qū)經(jīng)濟(jì)發(fā)展相對落后。未來應(yīng)在保護(hù)好生態(tài)用地的同時,加快經(jīng)濟(jì)發(fā)展和區(qū)域協(xié)作,尤其注重在大區(qū)域內(nèi)對生態(tài)用地的跨省合理規(guī)劃和利用,從而在一定程度上緩解重慶、四川等發(fā)展相對較快的省份的碳赤字。另外還應(yīng)減少或消滅沙塵源,恢復(fù)與再建干旱區(qū)綠色植被與荒漠植被。

    3 結(jié)論

    (1)自然生態(tài)系統(tǒng)碳動態(tài)的時空格局:2001—2009年全國NPP多年平均值約為3.32 Pg C/a,時間上總體波動不大,空間上呈西南地區(qū)>東南沿海>華中、華東地區(qū)>東北、華北地區(qū)>西北地區(qū)的格局。

    (2)社會經(jīng)濟(jì)系統(tǒng)碳排放的時空格局:2001—2009年全國能源消費(fèi)碳排放量逐年增加,年平均增長率為16.7%,多年平均值為2.53 Pg C/a;呈東部>中部>西部的空間格局。

    (3)自然-社會二元系統(tǒng)碳足跡和碳赤字的時空格局:碳足跡是自然-社會二元系統(tǒng)中能源結(jié)構(gòu)、產(chǎn)業(yè)結(jié)構(gòu)、土地利用結(jié)構(gòu)的綜合反映。2001—2009年全國平均碳足跡為6.98×106km2;碳足跡、碳赤字均逐年增加,碳足跡年均增長率為14.7%,自2009年開始出現(xiàn)正碳赤字。從研究時段平均狀況來看,具有正碳赤字的省份為山西、環(huán)渤海地區(qū)各省、長三角地區(qū)各省、廣東;具有碳盈余的省份多集中在中西部地區(qū)。碳赤字在全國范圍內(nèi)存在顯著的正空間自相關(guān);華東、西北、華北地區(qū)分別表現(xiàn)為高值集聚、低值集聚、高低間隔分布,其他區(qū)域空間自相關(guān)性較弱。相鄰省份碳赤字的相對大小由于互相影響而改變,因此通過加強(qiáng)區(qū)域聯(lián)系能使高值聚集區(qū)的“谷點(diǎn)”向周邊擴(kuò)散積極生態(tài)效應(yīng),或使低值聚集區(qū)的“峰點(diǎn)”受周邊積極生態(tài)效應(yīng)的影響而改善生態(tài)狀況,從而優(yōu)化整體生態(tài)效益。

    (4)全國分為中東部、南部、北部、西部4個生態(tài)經(jīng)濟(jì)大區(qū)。其中,中東部地區(qū)經(jīng)濟(jì)發(fā)展快而生態(tài)用地相對不足,急需探索低碳發(fā)展模式;南部地區(qū)經(jīng)濟(jì)發(fā)展逐漸加快且自然條件優(yōu)越,加強(qiáng)省際合作有助于積極生態(tài)效應(yīng)的擴(kuò)散;北部地區(qū)發(fā)展不均衡且自然條件較差,應(yīng)加強(qiáng)省際互補(bǔ)和資源合理保護(hù)與利用;西部地區(qū)發(fā)展較慢而生態(tài)用地相對充裕,應(yīng)在穩(wěn)定生態(tài)效益的同時提高經(jīng)濟(jì)效益。較之以往的生態(tài)區(qū)劃方案,本文對生態(tài)大區(qū)的劃分綜合考慮了自然生態(tài)系統(tǒng)對社會經(jīng)濟(jì)系統(tǒng)的響應(yīng),更有利于揭示各區(qū)資源環(huán)境和人類活動之間的耦合關(guān)系,為實(shí)現(xiàn)自然-社會二元系統(tǒng)的可持續(xù)發(fā)展提供科學(xué)依據(jù)。然而,根據(jù)本文目前的數(shù)據(jù)條件僅能作出省域尺度上的大區(qū)劃分,對于中小尺度的生態(tài)地區(qū)、生態(tài)區(qū)的劃分仍是需要繼續(xù)研究的內(nèi)容。

    致謝:感謝北京大學(xué)城市與環(huán)境學(xué)院趙淑清研究員,周德成、曾振中博士,以及住房和城鄉(xiāng)建設(shè)部城鄉(xiāng)規(guī)劃管理中心胡毅博士對本研究的幫助。

    [1] IPCC. Climate change 2007: the physical science basis // Solomon S, Qin D, Manning M, Chen Z, Marquis M, Averyt K B, Tignor M, Miller H L, eds. Contribution of Working GroupⅠ to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge: Cambridge University Press, 2007.

    [2] Fang J Y, Zhu J L, Wang S P, Yue C, Shen H H. Global warming, human-induced carbon emission, and their uncertainties. Science in China (Series D): Earth Sciences, 2011, 54(10): 1458- 1468.

    [3] Dixon R K, Solomon A M, Brown S, Houghton R A, Trexier M C, Wisniewski J. Carbon pools and flux of global forest ecosystems. Science, 1994, 263(5144): 185- 190.

    [4] Fang J Y, Chen A P, Peng C H, Zhao S Q, Ci L J. Changes in forest biomass carbon storage in China between 1949 and 1998. Science, 2001, 292(5525): 2320- 2322.

    [5] Houghton R A. Counting terrestrial sources and sinks of carbon. Climatic Change, 2001, 48(4): 525- 534.

    [6] Li K R, Wang S Q, Cao M K. Vegetation and soil carbon storage in China. Science in China (Series D): Earth Sciences, 2004, 47(1): 49- 57.

    [7] Ge Q S, Dai J H, He F N, Pan Y, Wang M M. Land use changes and their relations with carbon cycles over the past 300 a in China. Science in China (Series D): Earth Sciences, 2008, 51(6): 871- 884.

    [8] Soytas U, Sari R, Ewing B T. Energy consumption, income, and carbon emissions in the United States. Ecological Economics, 2007, 62(3- 4): 482- 489.

    [9] Xie H Y, Chen X S, Lin K R, Hu A Y. The ecological footprint analysis of fossil energy and electricity. Acta Ecologica Sinica, 2008, 28(4): 1729- 1735.

    [10] Zhu Y B, Wang Z, Pang L, Wang L J, Zhou X P. Simulation on China′s economy and prediction on energy consumption and carbon emission under optimal growth path. Acta Geographica Sinica, 2009, 64(8): 935- 944.

    [11] Zhao R Q, Huang X J. Carbon cycle of urban system: characteristics, mechanism and theoretical framework. Acta Ecologica Sinica, 2013, 33(2): 358- 366.

    [12] Wackernagel M, Ress W E. Our Ecological Footprint: Reducing Human Impact on the Earth. Gabriola Island: New Society Publishers, 1996.

    [13] Druckman A, Jackson T. The carbon footprint of UK households 1990- 2004: A socio-economically disaggregated, quasi-multi-regional input-output model. Ecological Economics, 2009, 68(7): 2066- 2077.

    [14] Sovacool B K, Brown M A. Twelve metropolitan carbon footprints: A preliminary comparative global assessment. Energy Policy, 2010, 38(9): 4856- 4869.

    [15] Kenny T, Gray N F. Comparative performance of six carbon footprint models for use in Ireland. Environmental Impact Assessment Review, 2009, 29(1): 1- 6.

    [16] Zhao R Q, Huang X J. Carbon emission and carbon footprint of different land use types based on energy consumption of Jiangsu province. Geographical Research, 2010, 29(9): 1639- 1649.

    [17] Zhao R Q, Huang X J, Zhong T Y. Research on carbon emission intensity and carbon footprint of different industrial spaces in China. Acta Geographica Sinica, 2010, 65(9): 1048- 1057.

    [18] Fan J, Li P X, Liang Y T. Framework of financial consumption oriented research on carbon footprints. Advances in Earth Science, 2010, 25(1): 61- 68.

    [19] Chuai X W, Lai L, Huang X J, Zhao R Q, Wang W J, Chen Z G. Temporospatial changes of carbon footprint based on energy consumption in China. Journal of Geographical Sciences, 2012, 22(1): 110- 124.

    [20] Shi M J, Wang Y, Zhang Z Y, Zhou X. Regional carbon footprint and interregional transfer of carbon emissions in China. Acta Geographica Sinica, 2012, 67(10): 1327- 1338.

    [21] Li Y, Huang X J, Zhen F. Effects of land use patterns on carbon emission in Jiangsu province. Transactions of the Chinese Society of Agricultural Engineering, 2008, 24(S2): 102- 107.

    [22] Lan J C, Fu W L, Yuan B, Zhang T, Peng J T. Analysis of land use patterns on carbon emission and carbon footprint in Chongqing city. Journal of Soil and Water Conservation, 2012, 26(1): 146- 155.

    [23] Zhang X H, Li X, Shi X, Liu K. Estimation of paddy biomass energy in Guangdong province. Remote Sensing Informatioin, 2007(1): 26- 29, 21- 21.

    [24] Li D K, Fan J Z, Wang J. Variation characteristics of vegetation net primary productivity in Shaanxi Province based on MOD17A3. Chinese Journal of Ecology, 2011, 30(12): 2776- 2782.

    [25] Turner D P, Ritts W D, Cohen W B, Gower S T, Running S W, Zhao M S, Costa M H, Kirschbaum A A, Ham J M, Saleska S R, Ahl D E. Evaluation of MODIS NPP and GPP products across multiple biomes. Remote Sensing of Environment, 2006, 102(3- 4): 282- 292.

    [26] Turner D P, Ritts W D, Cohen W B, Gower S T, Zhao M S, Running S W, Wofsy S C, Urbanski S, Dunn A L, Munger J W. Scaling Gross Primary Production (GPP) over boreal and deciduous forest landscapes in support of MODIS GPP product validation. Remote Sensing of Environment, 2003, 88(3): 256- 270.

    [27] Gebremichael M, Barros A P. Evaluation of MODIS gross primary productivity (GPP) in tropical monsoon regions. Remote Sensing of Environment, 2006, 100(2): 150- 166.

    [28] Guo X Y, He Y, Shen Y P, Feng D. Analysis of the terrestrial NPP based on the MODIS in the source regions of Yangtze and Yellow Rivers from 2000 to 2004. Journal of Glaciology and Geocryology, 2006, 28(4): 512- 518.

    [29] He Y, Dong W J, Guo X Y, Cao L J, Feng D. Terrestrial NPP variation in the region of the South-North Water Diversion Project (east route). Advances in Climate Change Research, 2006, 2(5): 246- 249.

    [30] Zhu F, Liu Z M, Wang Z M, Song K S. Temporal-spatial characteristics and factors influencing crop NPP across Northeastern China. Resouces Science, 2010, 32(11): 2079- 2084.

    [31] Guo Z X, Wang Z M, Zhang B, Liu D W, Yang G, Song K S, Li F. Analysis of temporal-spatial characteristics and factors influencing vegetation NPP in Northeast China from 2000 to 2006. Resource Science, 2008, 30(8): 1226- 1235.

    [32] Fang J Y, Ke J H, Tang Z Y, Chen A P. Implications and estimations of four terrestrial productivity parameters. Acta Phytoecologica Sinica, 2001, 25(4): 414- 419.

    [33] Field C B, Behrenfeld M J, Randerson J T, Falkowski P. Primary production of the biosphere: integrating terrestrial and oceanic components. Science, 281(5374): 237- 240.

    [34] Venetoulis J, Talberth J. Refining the ecological footprint. Environment, Development and Sustainability, 2008, 10(4): 441- 469.

    [35] Fang K, Shen W B, Dong D M. Modification and predication of energy ecological footprint: A case study of Jilin Province. Geographical Research, 2011, 30(10): 1835- 1846.

    [36] Zhou J H. The Research on the Cointegration and Causality between Land Use Change and Carbon Budget in Changsha City [D]. Changsha: Hunan Normal University, 2011.

    [37] Wang X N. Study on Estimation Method of Carbon Emission to Energy Carbon Sources in China [D]. Beijing: Beijing Forest University, 2006.

    [38] ORNL. Estimates of CO2Emissions from Fossil Fuel Burning and Cement Manufacturing. Tennessee: Oak Ridge National Laboratory, 1990.

    [39] Wang G, Feng X. CO2emission reduction through energy integration. Chemical Industry and Engineering Progress, 2006, 25(12): 1467- 1470.

    [40] Energy Research Institute of NDRC. Study on the sustainable energy development and carbon emission scenarios analysis in China. 2003.

    [41] Gao S T, Zhang H Q, Yang L R, Wang Q L. Forecast method and countermeasures for greenhouse gases in China. Research of Environmental Sciences, 1994, 7(6): 56- 59.

    [42] Global Footprint Network. GFN Ecological Footprint: Glossary. Oakland, CA, USA, 2007 [2013-08- 13]. http://www.footprintnetwork.org/en/index.php/GFN/page/carbon_footprint/.

    [43] Meng B, Wang J F, Zhang W Z, Liu X H. Evaluation of regional disparity in China based on spatial analysis. Scientia Geographica Sinica, 2005, 25(4): 393- 400.

    [44] Ouyang Z Y, Wang R S, Zheng H, Lin S K. Ecological culture in Hainan ecological construction. China Population, Resources and Environment, 2002, 12(4): 72- 74.

    [45] Wang S G, Zhou Y Z. Integrated frame research for resources & environment in megalopolises and regional sustainable development: case study on pearl river delta. China Population, Resources and Environment, 2002, 12(3): 54- 59.

    [46] Group of Ecological-Zonation, Research Centre of Ecology. On the principle of Ecological-economic zonation (Draft). Acta Ecologica Sinica, 1983, 3(2): 102- 102.

    [47] Fu B J, Liu G H, Chen L D, Ma K M, Li J R. Scheme of ecological regionalization in China. Acta Ecologica Sinica, 2001, 21(1): 1- 6.

    [48] Xie G D, Zhang C S, Zhang L B, Su D, Cao S Y, Leng Y F, Xiao Y. China′s county-scale ecological regionalization. Journal of Natural Resources, 2012, 27(1): 154- 162.

    參考文獻(xiàn):

    [2] 方精云, 朱江玲, 王少鵬, 岳超, 沈?;? 全球變暖、碳排放及不確定性. 中國科學(xué)(D輯:地球科學(xué)), 2011, 41(10): 1385- 1395.

    [6] 李克讓, 王紹強(qiáng), 曹明奎. 中國植被和土壤碳貯量. 中國科學(xué)(D輯:地球科學(xué)), 2003, 33(1): 72- 80.

    [7] 葛全勝, 戴君虎, 何凡能, 潘嫄, 王夢麥. 過去300年中國土地利用、土地覆被變化與碳循環(huán)研究. 中國科學(xué)(D輯:地球科學(xué)), 2008, 38(2): 197- 210.

    [9] 謝鴻宇, 陳賢生, 林凱榮, 胡安焱. 基于碳循環(huán)的化石能源及電力生態(tài)足跡. 生態(tài)學(xué)報(bào), 2008, 28(4): 1729- 1735.

    [10] 朱永彬, 王錚, 龐麗, 王麗娟, 鄒秀萍. 基于經(jīng)濟(jì)模擬的中國能源消費(fèi)與碳排放高峰預(yù)測. 地理學(xué)報(bào), 2009, 64(8): 935- 944.

    [11] 趙榮欽, 黃賢金. 城市系統(tǒng)碳循環(huán):特征、機(jī)理與理論框架. 生態(tài)學(xué)報(bào), 2013, 33(2): 358- 366.

    [16] 趙榮欽, 黃賢金. 基于能源消費(fèi)的江蘇省土地利用碳排放與碳足跡. 地理研究, 2010, 29(9): 1639- 1649.

    [17] 趙榮欽, 黃賢金, 鐘太洋. 中國不同產(chǎn)業(yè)空間的碳排放強(qiáng)度與碳足跡分析. 地理學(xué)報(bào), 2010, 65(9): 1048- 1057.

    [18] 樊杰, 李平星, 梁育填. 個人終端消費(fèi)導(dǎo)向的碳足跡研究框架——支撐我國環(huán)境外交的碳排放研究新思路. 地球科學(xué)進(jìn)展, 2010, 25(1): 61- 68.

    [20] 石敏俊, 王妍, 張卓穎, 周新. 中國各省區(qū)碳足跡與碳排放空間轉(zhuǎn)移. 地理學(xué)報(bào), 2012, 67(10): 1327- 1338.

    [21] 李穎, 黃賢金, 甄峰. 江蘇省區(qū)域不同土地利用方式的碳排放效應(yīng)分析. 農(nóng)業(yè)工程學(xué)報(bào), 2008, 24(S2): 102- 107.

    [22] 藍(lán)家程, 傅瓦利, 袁波, 張婷, 彭景濤. 重慶市不同土地利用碳排放及碳足跡分析. 水土保持學(xué)報(bào), 2012, 26(1): 146- 155.

    [23] 張曉浩, 黎夏, 施迅, 劉凱. 廣東省水稻生物質(zhì)能的估算. 遙感信息, 2007(1): 26- 29, 21- 21.

    [24] 李登科, 范建忠, 王娟. 基于MOD17A3的陜西省植被NPP變化特征. 生態(tài)學(xué)雜志, 2011, 30(12): 2776- 2782.

    [28] 郭曉寅, 何勇, 沈永平, 馮丁. 基于MODIS資料的2000- 2004年江河源區(qū)陸地植被凈初級生產(chǎn)力分析. 冰川凍土, 2006, 28(4): 512- 518.

    [29] 何勇, 董文杰, 郭曉寅, 曹麗娟, 馮丁. 我國南水北調(diào)東線地區(qū)陸地植被NPP變化特征. 氣候變化研究進(jìn)展, 2006, 2(5): 246- 249.

    [30] 朱鋒, 劉志明, 王宗明, 宋開山. 東北地區(qū)農(nóng)田凈初級生產(chǎn)力時空特征及其影響因素分析. 資源科學(xué), 2010, 32(11): 2079- 2084.

    [31] 國志興, 王宗明, 張柏, 劉殿偉, 楊桄, 宋開山, 李方. 2000年—2006年東北地區(qū)植被NPP的時空特征及影響因素分析. 資源科學(xué), 2008, 30(8): 1226- 1235.

    [32] 方精云, 柯金虎, 唐志堯, 陳安平. 生物生產(chǎn)力的“4P”概念、估算及其相互關(guān)系. 植物生態(tài)學(xué)報(bào), 2001, 25(4): 414- 419.

    [35] 方愷, 沈萬斌, 董德明. 能源足跡核算的改進(jìn)與預(yù)測——以吉林省為例. 地理研究, 2011, 30(10): 1835- 1846.

    [36] 周軍輝. 長沙市土地利用變化與碳收支協(xié)整性及因果關(guān)系研究 [D]. 長沙: 湖南師范大學(xué), 2011.

    [37] 王雪娜. 我國能源類碳源排碳量估算辦法研究 [D]. 北京: 北京林業(yè)大學(xué), 2006.

    [39] 汪剛, 馮霄. 基于能量集成的CO2減排量的確定. 化工進(jìn)展, 2006, 25(12): 1467- 1470.

    [40] 國家發(fā)展和改革委員會能源研究所. 中國可持續(xù)發(fā)展能源暨碳排放情景分析. 北京: 國家發(fā)展和改革委員會能源研究所, 2003.

    [41] 高樹婷, 張慧琴, 楊禮榮, 王秋玲. 我國溫室氣體排放量估測初探. 環(huán)境科學(xué)研究, 1994, 7(6): 56- 59.

    [43] 孟斌, 王勁峰, 張文忠, 劉旭華. 基于空間分析方法的中國區(qū)域差異研究. 地理科學(xué), 2005, 25(4): 393- 400.

    [44] 歐陽志云, 王如松, 鄭華, 林順坤. 海南生態(tài)文化建設(shè)探討. 中國人口·資源與環(huán)境, 2002, 12(4): 72- 74.

    [45] 王樹功, 周永章. 大城市群(圈)資源環(huán)境一體化與區(qū)域可持續(xù)發(fā)展研究——以珠江三角洲城市群為例. 中國人口·資源與環(huán)境, 2002, 12(3): 54- 59.

    [46] 生態(tài)學(xué)研究中心生態(tài)區(qū)劃組. 生態(tài)經(jīng)濟(jì)區(qū)劃原則(討論稿). 生態(tài)學(xué)報(bào), 1983, 3(2): 102- 102.

    [47] 傅伯杰, 劉國華, 陳利頂, 馬克明, 李俊然. 中國生態(tài)區(qū)劃方案. 生態(tài)學(xué)報(bào), 2001, 21(1): 1- 6.

    [48] 謝高地, 張昌順, 張林波, 蘇德, 曹淑艷, 冷允法, 肖玉. 保持縣域邊界完整性的中國生態(tài)區(qū)劃方案. 自然資源學(xué)報(bào), 2012, 27(1): 154- 162.

    Spatiotemporal patterns of carbon emission and carbon footprint in China during 2001—2009

    WU Wenjia1, JIANG Jinliang2, GAO Quanzhou3, JIANG Haibing4, *

    1CollegeofUrbanEnvironmentalSciences,andKeyLaboratoryforEarthSurfaceProcessesoftheMinistryofEducation,PekingUniversity,Beijing100871,China2SchoolofGeographicandOceanographicSciences,NanjingUniversity,Nanjing210093,China3GeographyandPlanningSchool,andGuangdongProvincialKeyLaboratoryforUrbanizationandGeo-simulation,SunYat-senUniversity,Guangzhou510275,China4SchoolofUrbanandResourceEnvironment,YanchengTeachersUniversity,Yancheng224051,China

    Increasing carbon emission, which has been recognized as one major reason for enhanced global warming, is influenced by complex interactions between natural and anthropogenic processes. The terrestrial ecosystem has been reported as a carbon sink during recent decades; However, this carbon sink has been largely offset by the carbon emission from human activities such as land use change and fossil fuel consumption. An accurate quantification of the carbon fluxes in natural ecosystems in response to human activities is of critical importance for global change study. Using a carbon footprint model, this study characterized the spatiotemporal carbon patterns in the coupled system of nature-society. Specifically, we firstly calculated the net primary productivity (NPP), carbon emission from energy consumption, carbon footprint, and carbon deficit in China during 2001—2009 based on MODIS and statistical data; Then we applied spatial auto-correlation analysis method to quantify spatiotemporal patterns of carbon emission and carbon footprint; Furthermore, we made an eco-economic regionalization of China based on the carbon deficit variation. Four conclusions can be drawn as follows. (1) The mean NPP was 3.32 Pg C/a during 2001—2009 in China (1 Pg = 1015g), with small annual variations and obvious spatial heterogeneities, ranking as the southwestern > the southeastern coastal region > the central and the eastern > the northeastern and the northern > the northwestern parts of China. (2) The carbon emission from energy consumption in China increased from 1.52 Pg C/a (2001) to 3.53 Pg C/a (2009) with an annual rate of 16.7% and an average value of 2.53 Pg C/a. Geographically, carbon emission decreased from east to west. (3) Carbon footprint in China increased from 4.46×106km2(2001) to 9.69×106km2(2009) with an annual rate of 14.7% and an average value of 6.98×106km2. Positive carbon deficit was found in Shanxi, Bohai Economic Zone, Yangtze River Delta and Guangdong, while negative carbon deficit existed in most provinces of the central and the western China. On the national scale, carbon deficit showed a significantly positive spatial autocorrelation; On the regional scale, however, the eastern China had a high-high clustering, the western China had a low-low clustering, and the northern China had a high-low distribution; Carbon deficit of adjacent provinces influenced each other, suggesting that interprovincial and interregional cooperation could contribute to a positive diffusion effect from the “hotspot” of carbon deficit to the surroundings. (4) China can be divided into four eco-economic zones based on spatiotemporal patterns and local spatial auto-correlation of carbon deficit. The first was the central-eastern region, which was in urgent need of low-carbon development mode considering its rapid economic development and limited availability for ecological resources. The second was the southern region, which had experienced fast economic development in past decades and a provincial cooperation is needed to strengthen the positive ecological effect diffusion in future. The third was the northern region, having unbalanced eco-economic development, where were required resources complementation among different provinces. The fourth was the western region with slow but accelerated economic development, which should stabilize ecological benefits and improve the economic effectiveness. Overall, our integrated modeling framework presented a spatial statistic method for assessing carbon uptake/ release patterns and dynamics in the coupled system of nature-society over large areas, which may help integrate biophysical and socioeconomic processes for sustainable development in different regions.

    the coupled system of nature-society; MODIS; carbon emission; carbon footprint; spatial autocorrelation; ecological economic regionalization

    國家自然科學(xué)基金青年項(xiàng)目(41301108); 南京大學(xué)研究生科研創(chuàng)新基金項(xiàng)目(2013CL07); 江蘇研究生創(chuàng)新計(jì)劃(CXLX13_033)

    2013- 02- 10; 網(wǎng)絡(luò)出版日期:2014- 03- 17

    10.5846/stxb201302100259

    *通訊作者Corresponding author.E-mail: jianghb@igsnrr.ac.cn

    吳文佳,蔣金亮,高全洲,蔣海兵.2001—2009年中國碳排放與碳足跡時空格局.生態(tài)學(xué)報(bào),2014,34(22):6722- 6733.

    Wu W J, Jiang J L, Gao Q Z, Jiang H B.Spatiotemporal patterns of carbon emission and carbon footprint in China during 2001—2009.Acta Ecologica Sinica,2014,34(2):6722- 6733.

    猜你喜歡
    赤字足跡空間
    空間是什么?
    創(chuàng)享空間
    成長足跡
    足跡
    心聲歌刊(2018年6期)2018-01-24 00:56:12
    足跡
    社會與公益(2016年2期)2016-04-13 02:49:06
    春的足跡
    互信赤字
    QQ空間那點(diǎn)事
    學(xué)與玩(2009年6期)2009-07-03 09:15:02
    空間
    97超级碰碰碰精品色视频在线观看| 久久婷婷人人爽人人干人人爱| 日本免费a在线| 欧美乱码精品一区二区三区| 精品一区二区三区视频在线 | 桃色一区二区三区在线观看| 一夜夜www| 美女午夜性视频免费| 亚洲国产高清在线一区二区三| 啦啦啦韩国在线观看视频| av女优亚洲男人天堂 | 亚洲欧洲精品一区二区精品久久久| 成人无遮挡网站| 日本精品一区二区三区蜜桃| 日韩人妻高清精品专区| 人人妻人人看人人澡| 精品久久蜜臀av无| 1024手机看黄色片| 中亚洲国语对白在线视频| 亚洲欧美日韩卡通动漫| 淫秽高清视频在线观看| 99精品欧美一区二区三区四区| 97碰自拍视频| 国产精品九九99| 欧美在线黄色| 欧美在线黄色| 色噜噜av男人的天堂激情| 美女免费视频网站| 色吧在线观看| 久久人人精品亚洲av| 熟女电影av网| 在线十欧美十亚洲十日本专区| 真人做人爱边吃奶动态| 久久久久精品国产欧美久久久| x7x7x7水蜜桃| 女警被强在线播放| 欧美激情在线99| 一区福利在线观看| 精品人妻1区二区| 国产精品98久久久久久宅男小说| 99热精品在线国产| 黑人欧美特级aaaaaa片| 国产亚洲欧美98| 黄色视频,在线免费观看| 国产精品久久久久久精品电影| 久久精品夜夜夜夜夜久久蜜豆| 国产成人av激情在线播放| 亚洲色图 男人天堂 中文字幕| 日韩成人在线观看一区二区三区| 色尼玛亚洲综合影院| 在线看三级毛片| 狂野欧美白嫩少妇大欣赏| 国产精品久久视频播放| 岛国在线免费视频观看| 叶爱在线成人免费视频播放| 在线视频色国产色| 国产爱豆传媒在线观看| 亚洲一区二区三区不卡视频| 日日摸夜夜添夜夜添小说| 亚洲欧洲精品一区二区精品久久久| 少妇的丰满在线观看| 欧美一级毛片孕妇| 色综合欧美亚洲国产小说| 夜夜夜夜夜久久久久| 亚洲欧美日韩卡通动漫| 狂野欧美激情性xxxx| 久久热在线av| 久久性视频一级片| 最近最新免费中文字幕在线| 国产一区二区激情短视频| 久久久成人免费电影| 国产男靠女视频免费网站| 亚洲自偷自拍图片 自拍| 欧美中文综合在线视频| 精品午夜福利视频在线观看一区| 欧美色视频一区免费| 99国产综合亚洲精品| 久久国产乱子伦精品免费另类| 久久精品91蜜桃| 少妇熟女aⅴ在线视频| 国产激情偷乱视频一区二区| 亚洲avbb在线观看| 99精品久久久久人妻精品| 母亲3免费完整高清在线观看| 国内毛片毛片毛片毛片毛片| 国产精品免费一区二区三区在线| 午夜影院日韩av| 国产爱豆传媒在线观看| av天堂中文字幕网| 国产精品一区二区精品视频观看| 国语自产精品视频在线第100页| 黄色成人免费大全| 小说图片视频综合网站| 91九色精品人成在线观看| 熟妇人妻久久中文字幕3abv| 亚洲国产欧美网| 欧美大码av| 99热这里只有是精品50| 久久午夜综合久久蜜桃| 中文在线观看免费www的网站| 首页视频小说图片口味搜索| 好男人在线观看高清免费视频| 亚洲国产精品成人综合色| 18禁黄网站禁片免费观看直播| 男女午夜视频在线观看| 亚洲真实伦在线观看| 久久久久久人人人人人| 欧美不卡视频在线免费观看| 在线十欧美十亚洲十日本专区| 中文资源天堂在线| 国产久久久一区二区三区| 久久热在线av| 欧美色欧美亚洲另类二区| 日韩高清综合在线| 国产成年人精品一区二区| 亚洲18禁久久av| 99国产综合亚洲精品| 嫁个100分男人电影在线观看| 欧美色欧美亚洲另类二区| 午夜成年电影在线免费观看| 国内久久婷婷六月综合欲色啪| 国产黄色小视频在线观看| 99精品在免费线老司机午夜| 欧美成人性av电影在线观看| 精品一区二区三区视频在线观看免费| 国产69精品久久久久777片 | 欧美又色又爽又黄视频| 亚洲av成人一区二区三| 99久久综合精品五月天人人| 九九热线精品视视频播放| 国产乱人视频| 国产av麻豆久久久久久久| 亚洲人成电影免费在线| 亚洲精品美女久久av网站| 18禁观看日本| 亚洲专区国产一区二区| 可以在线观看毛片的网站| 成人性生交大片免费视频hd| 麻豆国产97在线/欧美| 91字幕亚洲| 搡老熟女国产l中国老女人| 久久精品亚洲精品国产色婷小说| 婷婷丁香在线五月| 久久久国产成人免费| 嫁个100分男人电影在线观看| 搡老妇女老女人老熟妇| 一进一出好大好爽视频| av中文乱码字幕在线| 动漫黄色视频在线观看| 亚洲中文日韩欧美视频| 人妻丰满熟妇av一区二区三区| 婷婷丁香在线五月| 高清毛片免费观看视频网站| 香蕉av资源在线| 1000部很黄的大片| 国产亚洲av高清不卡| 精品久久久久久久人妻蜜臀av| 又爽又黄无遮挡网站| a在线观看视频网站| 夜夜夜夜夜久久久久| 久久久久久国产a免费观看| 国产一区二区在线观看日韩 | 国产成人一区二区三区免费视频网站| 性色avwww在线观看| 午夜福利免费观看在线| 九色成人免费人妻av| 深夜精品福利| 女人高潮潮喷娇喘18禁视频| 亚洲专区中文字幕在线| 999久久久精品免费观看国产| 久久精品亚洲精品国产色婷小说| 婷婷亚洲欧美| 成人18禁在线播放| 精品欧美国产一区二区三| 好男人在线观看高清免费视频| tocl精华| 最近在线观看免费完整版| 最近视频中文字幕2019在线8| 亚洲国产精品999在线| 日韩免费av在线播放| 99国产精品一区二区三区| 在线国产一区二区在线| 国产精品永久免费网站| www.熟女人妻精品国产| 狂野欧美激情性xxxx| 亚洲国产精品成人综合色| 男女做爰动态图高潮gif福利片| 亚洲专区国产一区二区| 中文资源天堂在线| 国产成人av激情在线播放| 亚洲avbb在线观看| 国产欧美日韩一区二区精品| 成人午夜高清在线视频| 91久久精品国产一区二区成人 | 亚洲一区二区三区不卡视频| 天堂网av新在线| 国产激情偷乱视频一区二区| 日本在线视频免费播放| 亚洲国产色片| 国产精品日韩av在线免费观看| 午夜福利在线在线| 搡老熟女国产l中国老女人| 亚洲男人的天堂狠狠| 一个人观看的视频www高清免费观看 | 黄色视频,在线免费观看| 亚洲av日韩精品久久久久久密| 啪啪无遮挡十八禁网站| 可以在线观看的亚洲视频| 91麻豆av在线| 亚洲一区高清亚洲精品| 欧美色视频一区免费| 给我免费播放毛片高清在线观看| 色播亚洲综合网| 国产精品爽爽va在线观看网站| 亚洲精品美女久久av网站| 国产免费av片在线观看野外av| 一进一出抽搐gif免费好疼| 老司机在亚洲福利影院| 欧美日本视频| 老司机午夜福利在线观看视频| 非洲黑人性xxxx精品又粗又长| 久久久久九九精品影院| 亚洲国产欧美一区二区综合| 免费人成视频x8x8入口观看| 色精品久久人妻99蜜桃| 国产精品久久久久久人妻精品电影| 亚洲七黄色美女视频| 舔av片在线| 亚洲国产欧洲综合997久久,| 少妇丰满av| 一级a爱片免费观看的视频| 久久这里只有精品中国| 亚洲片人在线观看| 五月玫瑰六月丁香| 亚洲av成人av| 嫩草影院精品99| 99精品欧美一区二区三区四区| 婷婷六月久久综合丁香| 久久国产乱子伦精品免费另类| 精品久久久久久,| 亚洲精品一区av在线观看| 国产精品香港三级国产av潘金莲| 日韩 欧美 亚洲 中文字幕| 欧美高清成人免费视频www| 精品午夜福利视频在线观看一区| 欧美xxxx黑人xx丫x性爽| 一夜夜www| 毛片女人毛片| 极品教师在线免费播放| 国产成人aa在线观看| 少妇的逼水好多| 欧美+亚洲+日韩+国产| 国产一区二区三区在线臀色熟女| 两人在一起打扑克的视频| 久久久久久久久中文| 亚洲精品在线观看二区| 99国产综合亚洲精品| 国产69精品久久久久777片 | 国产精品av视频在线免费观看| 日韩欧美免费精品| 久久草成人影院| 国产成+人综合+亚洲专区| 国产野战对白在线观看| 男女视频在线观看网站免费| 亚洲欧洲精品一区二区精品久久久| 12—13女人毛片做爰片一| 久久久久久久午夜电影| 97超级碰碰碰精品色视频在线观看| 香蕉久久夜色| 亚洲国产精品sss在线观看| 亚洲欧美日韩卡通动漫| 国产精品久久电影中文字幕| 日韩精品中文字幕看吧| 少妇丰满av| 欧美成人一区二区免费高清观看 | 国产免费av片在线观看野外av| 性色avwww在线观看| 亚洲在线自拍视频| 每晚都被弄得嗷嗷叫到高潮| 亚洲精品粉嫩美女一区| 国产美女午夜福利| 丁香六月欧美| 国产精品影院久久| 一夜夜www| 精品午夜福利视频在线观看一区| 老熟妇乱子伦视频在线观看| 久久天堂一区二区三区四区| 欧美乱妇无乱码| 啦啦啦观看免费观看视频高清| 精品国产亚洲在线| 国产亚洲欧美在线一区二区| 青草久久国产| 国产精品电影一区二区三区| 久久精品人妻少妇| 亚洲精品美女久久av网站| 三级毛片av免费| 欧美绝顶高潮抽搐喷水| 久久这里只有精品中国| 久久久久国内视频| 精品久久久久久久久久免费视频| ponron亚洲| 亚洲性夜色夜夜综合| 欧美色欧美亚洲另类二区| 国产亚洲精品av在线| 精华霜和精华液先用哪个| 国产精品av久久久久免费| 美女高潮的动态| 亚洲va日本ⅴa欧美va伊人久久| 男人舔女人下体高潮全视频| 中亚洲国语对白在线视频| 成人18禁在线播放| 亚洲 欧美 日韩 在线 免费| 亚洲中文av在线| 亚洲一区高清亚洲精品| 国产精品亚洲美女久久久| 国产免费av片在线观看野外av| 黄色片一级片一级黄色片| 天天一区二区日本电影三级| 色吧在线观看| 禁无遮挡网站| 9191精品国产免费久久| 欧美av亚洲av综合av国产av| 亚洲av熟女| 久久午夜综合久久蜜桃| 97碰自拍视频| 极品教师在线免费播放| 熟女少妇亚洲综合色aaa.| 精品国产亚洲在线| 亚洲片人在线观看| 亚洲国产精品成人综合色| 中文在线观看免费www的网站| 国产极品精品免费视频能看的| 淫秽高清视频在线观看| www.精华液| 国产一级毛片七仙女欲春2| 亚洲av成人av| 欧美丝袜亚洲另类 | 国产1区2区3区精品| 精品人妻1区二区| 国产成人精品久久二区二区91| 欧美大码av| 国产男靠女视频免费网站| 中国美女看黄片| 国产av麻豆久久久久久久| 丁香欧美五月| 久久九九热精品免费| 熟妇人妻久久中文字幕3abv| 国产精品,欧美在线| 国产成人av教育| 亚洲av五月六月丁香网| 天堂动漫精品| 亚洲成a人片在线一区二区| 亚洲欧美日韩东京热| 美女高潮喷水抽搐中文字幕| 老司机午夜十八禁免费视频| 亚洲欧美精品综合一区二区三区| 亚洲熟妇熟女久久| 女人被狂操c到高潮| 一区二区三区激情视频| 国产欧美日韩精品亚洲av| 亚洲国产中文字幕在线视频| 免费在线观看影片大全网站| 日韩欧美免费精品| 国产综合懂色| 91字幕亚洲| 成人无遮挡网站| av福利片在线观看| 不卡一级毛片| 国内毛片毛片毛片毛片毛片| 国产高清videossex| 神马国产精品三级电影在线观看| 岛国在线免费视频观看| 亚洲人成网站高清观看| 久久人人精品亚洲av| 嫁个100分男人电影在线观看| 中文字幕av在线有码专区| 一a级毛片在线观看| 香蕉国产在线看| 国产精品一区二区三区四区免费观看 | 欧美成狂野欧美在线观看| 亚洲一区二区三区不卡视频| 国产精品一及| 国产精品亚洲av一区麻豆| 又紧又爽又黄一区二区| 亚洲18禁久久av| 老司机午夜福利在线观看视频| 国内精品美女久久久久久| or卡值多少钱| 色尼玛亚洲综合影院| 日韩欧美 国产精品| 极品教师在线免费播放| 精品无人区乱码1区二区| 又紧又爽又黄一区二区| 欧美性猛交╳xxx乱大交人| 制服人妻中文乱码| 国产乱人视频| 国产精品 国内视频| 国产在线精品亚洲第一网站| 男女午夜视频在线观看| 男人的好看免费观看在线视频| 亚洲欧美日韩东京热| 亚洲av日韩精品久久久久久密| 法律面前人人平等表现在哪些方面| www.自偷自拍.com| 性色avwww在线观看| 久久精品亚洲精品国产色婷小说| 天堂动漫精品| 99热这里只有精品一区 | 色综合站精品国产| 国产人伦9x9x在线观看| 久久国产乱子伦精品免费另类| av女优亚洲男人天堂 | 色播亚洲综合网| 好看av亚洲va欧美ⅴa在| 天堂√8在线中文| 精品午夜福利视频在线观看一区| 欧美日韩黄片免| 亚洲黑人精品在线| 亚洲午夜理论影院| 99精品欧美一区二区三区四区| 欧美日韩亚洲国产一区二区在线观看| 亚洲国产精品合色在线| 精品久久久久久久毛片微露脸| 欧美国产日韩亚洲一区| av中文乱码字幕在线| xxxwww97欧美| 亚洲av电影在线进入| 国产真人三级小视频在线观看| 天天添夜夜摸| 18禁黄网站禁片午夜丰满| 黄色片一级片一级黄色片| 麻豆一二三区av精品| 又粗又爽又猛毛片免费看| 亚洲成人久久性| 哪里可以看免费的av片| 99国产精品一区二区蜜桃av| 欧美色欧美亚洲另类二区| 国产午夜精品论理片| 亚洲真实伦在线观看| 精品福利观看| 日韩大尺度精品在线看网址| 特级一级黄色大片| 免费av不卡在线播放| 亚洲 欧美一区二区三区| 色视频www国产| 女警被强在线播放| 欧美国产日韩亚洲一区| 无人区码免费观看不卡| 亚洲性夜色夜夜综合| 免费观看人在逋| 国产淫片久久久久久久久 | 桃红色精品国产亚洲av| 国产激情久久老熟女| 国产精品国产高清国产av| 久久中文看片网| 国产午夜精品久久久久久| 国产欧美日韩精品亚洲av| 亚洲片人在线观看| 欧美黑人欧美精品刺激| 麻豆国产av国片精品| 日本一本二区三区精品| 手机成人av网站| 国产精品久久久久久精品电影| 欧美xxxx黑人xx丫x性爽| 国产一区二区三区在线臀色熟女| 男人舔女人下体高潮全视频| 12—13女人毛片做爰片一| 男女午夜视频在线观看| 国产伦精品一区二区三区视频9 | 亚洲无线观看免费| 一二三四在线观看免费中文在| 久久精品国产综合久久久| 高清在线国产一区| 国产精品亚洲av一区麻豆| 99在线视频只有这里精品首页| 日韩av在线大香蕉| 九九热线精品视视频播放| 久久婷婷人人爽人人干人人爱| 97碰自拍视频| 成年女人看的毛片在线观看| 成人永久免费在线观看视频| 国内精品久久久久精免费| 亚洲中文字幕日韩| 最新中文字幕久久久久 | 一本综合久久免费| 黑人巨大精品欧美一区二区mp4| 999久久久国产精品视频| 久久国产乱子伦精品免费另类| 1024香蕉在线观看| 久久中文字幕一级| 香蕉国产在线看| 老司机福利观看| 免费一级毛片在线播放高清视频| 亚洲最大成人中文| 国产精品综合久久久久久久免费| 美女免费视频网站| 国产欧美日韩精品亚洲av| 国产高清有码在线观看视频| 国产淫片久久久久久久久 | 国产精品免费一区二区三区在线| 99国产综合亚洲精品| 国产亚洲av嫩草精品影院| 久久草成人影院| 国产激情偷乱视频一区二区| av福利片在线观看| 99国产极品粉嫩在线观看| 国产高清激情床上av| 亚洲av日韩精品久久久久久密| 中亚洲国语对白在线视频| 99久久国产精品久久久| 中文字幕最新亚洲高清| 欧美绝顶高潮抽搐喷水| 欧美一区二区精品小视频在线| 草草在线视频免费看| 国产精品久久久人人做人人爽| 听说在线观看完整版免费高清| 一边摸一边抽搐一进一小说| 桃色一区二区三区在线观看| 男女那种视频在线观看| 欧美日韩瑟瑟在线播放| 人人妻,人人澡人人爽秒播| 欧美黑人欧美精品刺激| 啪啪无遮挡十八禁网站| 亚洲人与动物交配视频| 久久人妻av系列| 网址你懂的国产日韩在线| 99国产极品粉嫩在线观看| 国产欧美日韩精品亚洲av| 18美女黄网站色大片免费观看| x7x7x7水蜜桃| 中文字幕久久专区| 在线a可以看的网站| 久久久久久人人人人人| 精品乱码久久久久久99久播| 国产成人精品无人区| 欧美丝袜亚洲另类 | 老司机在亚洲福利影院| 亚洲成人精品中文字幕电影| 国产一区二区在线av高清观看| 一本久久中文字幕| www.熟女人妻精品国产| 999精品在线视频| 国产欧美日韩一区二区三| 99久久无色码亚洲精品果冻| 少妇的丰满在线观看| 看片在线看免费视频| 久久香蕉精品热| 国产精品av视频在线免费观看| 狠狠狠狠99中文字幕| 非洲黑人性xxxx精品又粗又长| 久久国产精品人妻蜜桃| 国内精品一区二区在线观看| 久久欧美精品欧美久久欧美| 国产欧美日韩精品亚洲av| 午夜福利视频1000在线观看| 99热这里只有精品一区 | 麻豆一二三区av精品| 久久久久久人人人人人| 长腿黑丝高跟| 日韩高清综合在线| 久久伊人香网站| 99国产精品一区二区三区| 亚洲成人精品中文字幕电影| 久久中文字幕一级| 国产精品亚洲av一区麻豆| 成人永久免费在线观看视频| 久久人人精品亚洲av| 男女之事视频高清在线观看| 热99re8久久精品国产| 国产精品电影一区二区三区| av在线蜜桃| 女生性感内裤真人,穿戴方法视频| 亚洲色图 男人天堂 中文字幕| 国产成人啪精品午夜网站| 757午夜福利合集在线观看| 最好的美女福利视频网| 天天添夜夜摸| 国产单亲对白刺激| 亚洲最大成人中文| 最新中文字幕久久久久 | 亚洲九九香蕉| 在线看三级毛片| 一本综合久久免费| 亚洲欧美日韩无卡精品| 欧美日韩乱码在线| 欧美一级毛片孕妇| 又黄又爽又免费观看的视频| 国产又色又爽无遮挡免费看| 欧美在线一区亚洲| 国产激情久久老熟女| 欧美大码av| 成人av一区二区三区在线看| 久久久久久国产a免费观看| svipshipincom国产片| 亚洲成人精品中文字幕电影| 久99久视频精品免费| h日本视频在线播放| 俄罗斯特黄特色一大片| 欧美最黄视频在线播放免费| 久久精品国产清高在天天线| 欧美3d第一页| 午夜福利高清视频| 日本一二三区视频观看| 欧美日韩精品网址| 国内毛片毛片毛片毛片毛片| 国产1区2区3区精品| 88av欧美| 蜜桃久久精品国产亚洲av| 亚洲av免费在线观看| 一级毛片高清免费大全| 国产成人精品久久二区二区91| 日韩中文字幕欧美一区二区| 久久久久国内视频| 小说图片视频综合网站| 禁无遮挡网站|