• <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
    空間
    亚洲精品乱码久久久久久按摩| 欧美成人a在线观看| 成人性生交大片免费视频hd| 丝袜喷水一区| 欧美最新免费一区二区三区| 男人舔女人下体高潮全视频| 一级毛片久久久久久久久女| 桃色一区二区三区在线观看| 国产午夜精品论理片| 22中文网久久字幕| 国产亚洲91精品色在线| 国产亚洲av嫩草精品影院| 亚洲图色成人| 变态另类成人亚洲欧美熟女| 男人舔女人下体高潮全视频| 69av精品久久久久久| 干丝袜人妻中文字幕| 久久久久久久久中文| 一区二区三区高清视频在线| 国产一区二区激情短视频| 亚洲va在线va天堂va国产| 亚洲久久久久久中文字幕| 国产伦一二天堂av在线观看| 欧美日韩在线观看h| 97在线视频观看| 日韩亚洲欧美综合| 久久久精品大字幕| 国产午夜福利久久久久久| 色综合亚洲欧美另类图片| 免费人成视频x8x8入口观看| 亚洲图色成人| 亚洲中文字幕日韩| 亚洲18禁久久av| 久久午夜福利片| 中文字幕制服av| 不卡视频在线观看欧美| 国产不卡一卡二| 欧美性猛交╳xxx乱大交人| 一进一出抽搐动态| kizo精华| 亚洲精品成人久久久久久| 18禁在线播放成人免费| 欧美3d第一页| 日日摸夜夜添夜夜添av毛片| 欧美日韩在线观看h| 日韩强制内射视频| 麻豆av噜噜一区二区三区| 在线国产一区二区在线| 国产亚洲5aaaaa淫片| 久久精品综合一区二区三区| 国产日本99.免费观看| 欧美变态另类bdsm刘玥| 欧美又色又爽又黄视频| 亚洲最大成人av| 可以在线观看的亚洲视频| 青春草视频在线免费观看| 国产一级毛片七仙女欲春2| 永久网站在线| 少妇的逼好多水| 欧美色视频一区免费| 亚洲精品乱码久久久v下载方式| 国产精品99久久久久久久久| 久久精品国产自在天天线| 91在线精品国自产拍蜜月| 91精品国产九色| 久久精品国产亚洲av涩爱 | 能在线免费看毛片的网站| 蜜桃久久精品国产亚洲av| 国产精品一区二区三区四区免费观看| 国产 一区精品| 嫩草影院新地址| 嘟嘟电影网在线观看| 特级一级黄色大片| 久久亚洲精品不卡| 久久国内精品自在自线图片| 99久久精品一区二区三区| 午夜福利高清视频| 国产精品女同一区二区软件| 一卡2卡三卡四卡精品乱码亚洲| 看十八女毛片水多多多| 亚洲av男天堂| 亚洲国产精品成人久久小说 | 一夜夜www| 国产v大片淫在线免费观看| 看免费成人av毛片| 亚洲欧美成人综合另类久久久 | 久久久久久久久久久丰满| 变态另类成人亚洲欧美熟女| 变态另类丝袜制服| 麻豆av噜噜一区二区三区| 国产探花极品一区二区| 99久久中文字幕三级久久日本| 男人的好看免费观看在线视频| 国产综合懂色| 国产精品一区二区性色av| 天堂√8在线中文| 亚洲欧美中文字幕日韩二区| 亚洲人成网站在线观看播放| 免费观看a级毛片全部| 女人十人毛片免费观看3o分钟| 亚洲av不卡在线观看| 黄色一级大片看看| 亚洲av男天堂| 非洲黑人性xxxx精品又粗又长| 97超碰精品成人国产| 国产午夜精品久久久久久一区二区三区| 69人妻影院| 免费黄网站久久成人精品| 在线a可以看的网站| 99热只有精品国产| 日本撒尿小便嘘嘘汇集6| 内射极品少妇av片p| 少妇丰满av| 日韩欧美精品免费久久| 国产精品野战在线观看| 91av网一区二区| 六月丁香七月| 成人一区二区视频在线观看| 久久99精品国语久久久| 性插视频无遮挡在线免费观看| 国产又黄又爽又无遮挡在线| 日本免费一区二区三区高清不卡| 国产黄片视频在线免费观看| 亚洲在线自拍视频| 自拍偷自拍亚洲精品老妇| 成人午夜高清在线视频| 精品久久久噜噜| 波多野结衣高清无吗| 18禁在线播放成人免费| 国产综合懂色| 自拍偷自拍亚洲精品老妇| 亚洲欧美日韩卡通动漫| 青春草视频在线免费观看| 午夜福利在线在线| 日日摸夜夜添夜夜爱| 国产乱人偷精品视频| 精品熟女少妇av免费看| 久久综合国产亚洲精品| av在线老鸭窝| 日本与韩国留学比较| 亚洲av熟女| 精品人妻一区二区三区麻豆| av在线观看视频网站免费| 日韩av在线大香蕉| 最近最新中文字幕大全电影3| 一区二区三区高清视频在线| 日日摸夜夜添夜夜爱| 久久久a久久爽久久v久久| 边亲边吃奶的免费视频| 两个人视频免费观看高清| 中文字幕人妻熟人妻熟丝袜美| 国产免费一级a男人的天堂| a级毛片免费高清观看在线播放| 啦啦啦啦在线视频资源| 亚洲中文字幕日韩| 淫秽高清视频在线观看| 直男gayav资源| 一区二区三区高清视频在线| 夜夜爽天天搞| 在线国产一区二区在线| 久久久成人免费电影| 久久国内精品自在自线图片| 在线播放无遮挡| 亚洲美女视频黄频| 99热6这里只有精品| 欧美不卡视频在线免费观看| 欧美一区二区国产精品久久精品| 高清在线视频一区二区三区 | 国产精品福利在线免费观看| 亚洲欧美清纯卡通| 亚洲,欧美,日韩| 日韩制服骚丝袜av| 久久精品国产亚洲网站| 国产熟女欧美一区二区| 日本爱情动作片www.在线观看| 国产又黄又爽又无遮挡在线| 久久99精品国语久久久| 久久久色成人| 婷婷亚洲欧美| 欧美性猛交黑人性爽| 国产精品久久久久久久久免| 一区二区三区高清视频在线| 久久午夜福利片| 成人午夜精彩视频在线观看| 色综合亚洲欧美另类图片| 可以在线观看的亚洲视频| 国内精品一区二区在线观看| 亚洲熟妇中文字幕五十中出| 欧美高清性xxxxhd video| 国产高潮美女av| 岛国在线免费视频观看| 欧美变态另类bdsm刘玥| 久久婷婷人人爽人人干人人爱| 精品一区二区三区视频在线| .国产精品久久| 一级毛片电影观看 | 联通29元200g的流量卡| 亚洲欧美成人综合另类久久久 | 人人妻人人看人人澡| avwww免费| 麻豆成人午夜福利视频| 久久精品影院6| 成年免费大片在线观看| 欧美一区二区精品小视频在线| 日日干狠狠操夜夜爽| 男女那种视频在线观看| 久久精品国产自在天天线| 国产精品久久久久久精品电影小说 | 真实男女啪啪啪动态图| 男人舔女人下体高潮全视频| 99热这里只有是精品在线观看| 黄色日韩在线| 日韩强制内射视频| 免费观看精品视频网站| 精品久久久久久久久亚洲| 国产老妇伦熟女老妇高清| 一个人看视频在线观看www免费| 日韩欧美三级三区| 亚洲美女搞黄在线观看| 国产真实伦视频高清在线观看| 男女视频在线观看网站免费| 亚洲,欧美,日韩| 日韩三级伦理在线观看| 亚洲人成网站在线播| 男女啪啪激烈高潮av片| 午夜老司机福利剧场| 色5月婷婷丁香| 97超视频在线观看视频| 男人狂女人下面高潮的视频| 久久精品人妻少妇| 国内精品宾馆在线| 国产精品一区二区在线观看99 | 中文欧美无线码| 国内精品宾馆在线| 精品不卡国产一区二区三区| 狠狠狠狠99中文字幕| 国产精品野战在线观看| 精品午夜福利在线看| 国产v大片淫在线免费观看| 国产高潮美女av| 性插视频无遮挡在线免费观看| 中文字幕av在线有码专区| 日韩精品有码人妻一区| 国产精品伦人一区二区| 亚洲真实伦在线观看| 最后的刺客免费高清国语| 中文资源天堂在线| 精品久久久久久久久av| 亚洲av成人精品一区久久| 亚洲国产精品sss在线观看| 简卡轻食公司| 免费观看在线日韩| 能在线免费看毛片的网站| 久久久久久久久久成人| 中国国产av一级| 久久久精品大字幕| 久久中文看片网| 女的被弄到高潮叫床怎么办| 国产乱人视频| 亚洲四区av| 岛国毛片在线播放| 久久久久久久午夜电影| 国内精品久久久久精免费| 国产高清激情床上av| 成年免费大片在线观看| 久久久久久九九精品二区国产| 在线观看午夜福利视频| 国产精品1区2区在线观看.| 婷婷精品国产亚洲av| 亚洲精品国产av成人精品| 日韩成人伦理影院| 天堂网av新在线| 日韩大尺度精品在线看网址| 五月伊人婷婷丁香| 麻豆成人午夜福利视频| 欧美最新免费一区二区三区| 久久久久久国产a免费观看| 一级黄片播放器| 国产精品99久久久久久久久| 草草在线视频免费看| 久久久久网色| 亚洲av二区三区四区| 国产国拍精品亚洲av在线观看| 久久精品夜色国产| 亚洲国产高清在线一区二区三| 九九热线精品视视频播放| 欧美高清性xxxxhd video| 日韩av不卡免费在线播放| 亚洲av熟女| 一级毛片aaaaaa免费看小| 美女大奶头视频| 日韩亚洲欧美综合| 国产精品1区2区在线观看.| 变态另类丝袜制服| 老女人水多毛片| 久久99蜜桃精品久久| 18禁裸乳无遮挡免费网站照片| 天堂av国产一区二区熟女人妻| 少妇丰满av| 51国产日韩欧美| 美女脱内裤让男人舔精品视频 | 久久精品国产亚洲av天美| 亚洲激情五月婷婷啪啪| 给我免费播放毛片高清在线观看| 色综合色国产| 免费观看人在逋| 成人一区二区视频在线观看| 99热网站在线观看| 午夜福利成人在线免费观看| 99国产极品粉嫩在线观看| 午夜福利视频1000在线观看| 天天一区二区日本电影三级| 成年女人永久免费观看视频| 国产成人a区在线观看| 国产真实伦视频高清在线观看| 99在线视频只有这里精品首页| 免费无遮挡裸体视频| 久久久精品94久久精品| 国产av一区在线观看免费| 亚洲自拍偷在线| 精品一区二区三区视频在线| 夜夜看夜夜爽夜夜摸| 一区二区三区四区激情视频 | 精品久久久久久久久av| 国产 一区 欧美 日韩| 看非洲黑人一级黄片| 国产视频内射| 自拍偷自拍亚洲精品老妇| 久久久国产成人免费| 亚洲人成网站在线播放欧美日韩| 亚洲精品影视一区二区三区av| 成人永久免费在线观看视频| 国产一区二区三区在线臀色熟女| a级一级毛片免费在线观看| 国产高清视频在线观看网站| 最近最新中文字幕大全电影3| 国产真实伦视频高清在线观看| 午夜激情欧美在线| 亚洲中文字幕日韩| 免费看美女性在线毛片视频| 成人一区二区视频在线观看| 亚洲18禁久久av| 美女 人体艺术 gogo| 亚洲色图av天堂| 日韩成人av中文字幕在线观看| 国产精品人妻久久久久久| 国产国拍精品亚洲av在线观看| 一本精品99久久精品77| 国产成人91sexporn| 成熟少妇高潮喷水视频| 五月伊人婷婷丁香| 精品99又大又爽又粗少妇毛片| 成人三级黄色视频| 久久精品久久久久久久性| 又粗又爽又猛毛片免费看| 久久久国产成人免费| 不卡一级毛片| 久久久久网色| 久久热精品热| 日韩 亚洲 欧美在线| 国内精品美女久久久久久| 性色avwww在线观看| 夫妻性生交免费视频一级片| 91久久精品国产一区二区三区| 中文字幕人妻熟人妻熟丝袜美| 精品国产三级普通话版| 免费黄网站久久成人精品| 3wmmmm亚洲av在线观看| 午夜精品国产一区二区电影 | 最近的中文字幕免费完整| 人妻夜夜爽99麻豆av| 国产精品嫩草影院av在线观看| 99视频精品全部免费 在线| 国产成人a∨麻豆精品| 日韩,欧美,国产一区二区三区 | 亚洲欧美精品综合久久99| 久久人人精品亚洲av| 三级毛片av免费| 国产高潮美女av| 日本欧美国产在线视频| 色哟哟哟哟哟哟| 男人和女人高潮做爰伦理| 欧美三级亚洲精品| 国产老妇伦熟女老妇高清| 久久精品91蜜桃| 午夜福利视频1000在线观看| 成人毛片60女人毛片免费| av女优亚洲男人天堂| 亚洲成人久久性| 国产老妇伦熟女老妇高清| 亚洲成人久久性| 麻豆成人午夜福利视频| 日韩精品青青久久久久久| 少妇熟女欧美另类| 国产精品福利在线免费观看| 婷婷色av中文字幕| 久久精品国产自在天天线| 亚洲av不卡在线观看| 欧美日韩精品成人综合77777| 99久国产av精品| 免费观看a级毛片全部| 欧美成人免费av一区二区三区| 97人妻精品一区二区三区麻豆| 女人被狂操c到高潮| 亚洲欧美日韩无卡精品| 亚洲av不卡在线观看| 久久亚洲精品不卡| 真实男女啪啪啪动态图| 亚洲熟妇中文字幕五十中出| 国产精品野战在线观看| 淫秽高清视频在线观看| 色5月婷婷丁香| 欧美一区二区亚洲| 久久鲁丝午夜福利片| 69人妻影院| 99久久九九国产精品国产免费| 亚洲国产精品成人综合色| 白带黄色成豆腐渣| 午夜老司机福利剧场| 九九在线视频观看精品| 男人和女人高潮做爰伦理| 欧美一区二区精品小视频在线| 国产人妻一区二区三区在| 亚洲精品日韩在线中文字幕 | 久久精品国产亚洲av香蕉五月| 亚洲图色成人| 97超视频在线观看视频| 成人一区二区视频在线观看| 亚洲最大成人av| 日韩成人伦理影院| 丝袜美腿在线中文| 国产老妇女一区| 亚洲欧美清纯卡通| 国产高潮美女av| 秋霞在线观看毛片| 久久国产乱子免费精品| 亚洲图色成人| 国产高清视频在线观看网站| 日韩欧美在线乱码| 精品免费久久久久久久清纯| 三级经典国产精品| 国产亚洲精品久久久com| 国产精品一区二区性色av| 国产精品久久久久久精品电影小说 | 人人妻人人看人人澡| 中文字幕熟女人妻在线| 成人永久免费在线观看视频| 嫩草影院入口| 99久久久亚洲精品蜜臀av| 久久人人爽人人片av| 在线观看一区二区三区| 国产精品久久电影中文字幕| 精品午夜福利在线看| 欧美最黄视频在线播放免费| 久久久精品欧美日韩精品| 国产伦一二天堂av在线观看| 大香蕉久久网| 一本久久中文字幕| 精品人妻偷拍中文字幕| 久久精品久久久久久久性| 日韩精品有码人妻一区| 真实男女啪啪啪动态图| 91久久精品国产一区二区三区| 搞女人的毛片| 好男人视频免费观看在线| 国产午夜福利久久久久久| 成人av在线播放网站| 在线观看av片永久免费下载| 亚洲欧美日韩高清专用| 欧美变态另类bdsm刘玥| 亚洲国产日韩欧美精品在线观看| 九九久久精品国产亚洲av麻豆| 一级av片app| 日韩人妻高清精品专区| 性色avwww在线观看| 一个人看的www免费观看视频| 亚洲国产精品成人综合色| 精品不卡国产一区二区三区| 国产精品一区www在线观看| 村上凉子中文字幕在线| 日本三级黄在线观看| 亚洲自偷自拍三级| 在线播放无遮挡| 日韩av在线大香蕉| 12—13女人毛片做爰片一| 桃色一区二区三区在线观看| 久久久久久久久大av| 国产在线男女| 男女视频在线观看网站免费| 身体一侧抽搐| 欧美丝袜亚洲另类| 免费观看a级毛片全部| 一本精品99久久精品77| 久久久精品欧美日韩精品| 免费无遮挡裸体视频| 亚洲欧美日韩高清专用| 黄片无遮挡物在线观看| 精品人妻一区二区三区麻豆| 男女边吃奶边做爰视频| www.av在线官网国产| 一本一本综合久久| 少妇人妻一区二区三区视频| 国产在线精品亚洲第一网站| 我要搜黄色片| 免费观看在线日韩| 干丝袜人妻中文字幕| 国产成人91sexporn| 欧美变态另类bdsm刘玥| 六月丁香七月| 亚洲精品乱码久久久久久按摩| 人妻久久中文字幕网| av又黄又爽大尺度在线免费看 | 国产精品野战在线观看| 91久久精品国产一区二区成人| 午夜福利在线观看免费完整高清在 | 成年版毛片免费区| 国产老妇伦熟女老妇高清| 亚洲国产精品sss在线观看| 国产精品麻豆人妻色哟哟久久 | 国产精品三级大全| 黄色一级大片看看| 成人av在线播放网站| 国产一区二区在线观看日韩| 国产精品电影一区二区三区| 欧美在线一区亚洲| 国产女主播在线喷水免费视频网站 | 舔av片在线| 一个人看的www免费观看视频| 久99久视频精品免费| eeuss影院久久| 99热精品在线国产| 少妇猛男粗大的猛烈进出视频 | 成人亚洲精品av一区二区| 国产国拍精品亚洲av在线观看| 亚洲最大成人中文| 亚洲欧美日韩高清专用| 看非洲黑人一级黄片| 免费人成在线观看视频色| 校园人妻丝袜中文字幕| 亚洲一区二区三区色噜噜| 国产精品久久久久久久久免| 少妇的逼好多水| 国产伦在线观看视频一区| 在线观看美女被高潮喷水网站| 伦精品一区二区三区| 久久久久久大精品| 波多野结衣巨乳人妻| 亚洲无线在线观看| 精品少妇黑人巨大在线播放 | 插阴视频在线观看视频| 欧美区成人在线视频| 国产淫片久久久久久久久| 不卡一级毛片| 1024手机看黄色片| 欧美一区二区精品小视频在线| 麻豆成人午夜福利视频| 三级国产精品欧美在线观看| 国产片特级美女逼逼视频| 亚洲av中文av极速乱| 久久久久久久久久久丰满| 国产精品爽爽va在线观看网站| 女同久久另类99精品国产91| 国产日本99.免费观看| 丰满人妻一区二区三区视频av| 欧美+日韩+精品| 毛片女人毛片| 91av网一区二区| 成人欧美大片| 性插视频无遮挡在线免费观看| 久久久午夜欧美精品| 中文在线观看免费www的网站| 国产精品精品国产色婷婷| 天堂中文最新版在线下载 | 又粗又爽又猛毛片免费看| 级片在线观看| 免费看美女性在线毛片视频| av在线蜜桃| 欧美极品一区二区三区四区| 久久久久久大精品| 免费观看在线日韩| 午夜福利成人在线免费观看| 久久久久久大精品| 日本黄色视频三级网站网址| 晚上一个人看的免费电影| 少妇的逼好多水| 夫妻性生交免费视频一级片| 性欧美人与动物交配| 久久人人爽人人爽人人片va| 少妇高潮的动态图| 非洲黑人性xxxx精品又粗又长| 亚洲av.av天堂| 亚洲在久久综合| 能在线免费观看的黄片| 精品一区二区免费观看| 色5月婷婷丁香| 在线免费十八禁| 亚洲美女视频黄频| 91精品一卡2卡3卡4卡| 久久久色成人| 一区福利在线观看| 国产成人午夜福利电影在线观看| 亚洲四区av| 草草在线视频免费看| 亚洲性久久影院| 1024手机看黄色片| 菩萨蛮人人尽说江南好唐韦庄 | 国产熟女欧美一区二区| 欧美成人免费av一区二区三区| 91久久精品国产一区二区三区| 欧美日韩一区二区视频在线观看视频在线 | 内射极品少妇av片p| 熟女人妻精品中文字幕| 久久中文看片网| 99久久九九国产精品国产免费| 国产精品一二三区在线看|