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

    基于碳儲量視角的城鎮(zhèn)土地利用模式與生態(tài)效益分異特征

    2017-06-05 14:19:10郭旭東包慧敏呂春艷張衍毓
    中國土地科學(xué) 2017年4期
    關(guān)鍵詞:儲量城鎮(zhèn)土地利用

    郭旭東,常 青,劉 筱,包慧敏,孫 藝,呂春艷,張衍毓

    (1.中國土地勘測規(guī)劃院,國土資源部土地利用重點實驗室,北京 100035;2.中國農(nóng)業(yè)大學(xué)觀賞園藝與園林系,北京 100193)

    基于碳儲量視角的城鎮(zhèn)土地利用模式與生態(tài)效益分異特征

    郭旭東1,常 青2,劉 筱2,包慧敏2,孫 藝2,呂春艷1,張衍毓1

    (1.中國土地勘測規(guī)劃院,國土資源部土地利用重點實驗室,北京 100035;2.中國農(nóng)業(yè)大學(xué)觀賞園藝與園林系,北京 100193)

    研究目的:探析不同土地利用模式下碳儲量的空間差異特征,為優(yōu)化城鎮(zhèn)土地利用及其生態(tài)效益提供參考。研究方法:InVEST碳儲量修正模型、歸一化不透水地表指數(shù)(NDISI)和層次聚類分析。研究結(jié)果:(1)2014年北京市三山五園地區(qū)碳儲量為23.31×104t,單位土地面積平均碳儲量32.28 t/hm2;各街道碳儲量差異明顯。(2)以NDISI變化率表征土地利用開發(fā)強(qiáng)度,1992—2014年間上地街道開發(fā)強(qiáng)度最高,青龍橋街道(頤和園)開發(fā)強(qiáng)度最低。(3)該區(qū)各街道碳儲量與其區(qū)位、用地結(jié)構(gòu)及開發(fā)建設(shè)模式密切相關(guān),開發(fā)建設(shè)活動并非一定會降低城鎮(zhèn)土地碳匯能力。研究結(jié)論:土地利用模式影響城鎮(zhèn)土地碳儲能力,加強(qiáng)城鎮(zhèn)土地利用模式與生態(tài)效益研究,有助于合理制定城鎮(zhèn)土地利用政策。

    土地利用;碳儲量;生態(tài)效益;土地開發(fā)模式;三山五園

    1 引言

    近40年來的快速城市化過程中,中國土地資源開發(fā)利用取得了顯著的社會經(jīng)濟(jì)效益,但同時也出現(xiàn)了自然棲息地退化、水土流失、城市熱島等生態(tài)環(huán)境問題[1-2]。城鎮(zhèn)土地利用生態(tài)效益隨之成為中國土地資源可持續(xù)利用研究的熱點[3-7]。目前,國內(nèi)外學(xué)者通過構(gòu)建指標(biāo)模型、圍繞城市土地利用生態(tài)效益評價開展了大量研究,多數(shù)研究采用森林或綠地覆蓋率[4-7]、環(huán)境污染物排放量[4-7]、水土流失或生物豐富度[4-5]、土地利用規(guī)模與結(jié)構(gòu)[6-7]等指標(biāo),偶見采用生態(tài)系統(tǒng)服務(wù)或生態(tài)綠當(dāng)量指標(biāo)[3,7],對土地利用生態(tài)效益的直接度量指標(biāo)缺乏;研究內(nèi)容多集中于整個城市土地利用生態(tài)效益動態(tài)變化[3-5]或空間差異上[6-7],且對于土地利用規(guī)模與結(jié)構(gòu)指標(biāo)考慮居多[4-7],忽略了城鎮(zhèn)土地利用模式(特別是在相似土地利用規(guī)模與結(jié)構(gòu)下不同的開發(fā)建設(shè)模式)對生態(tài)系統(tǒng)服務(wù)能力的影響。不透水地表擴(kuò)張會導(dǎo)致城市地域植被與土壤有機(jī)碳的下降[8-9]、地表徑流的增加[10]和熱島效應(yīng)的增強(qiáng)[11-12];而改變土地利用開發(fā)模式有助于改善以上城市生態(tài)系統(tǒng)過程,提高土地利用服務(wù)能力。例如,低影響開發(fā)作為一類城市雨洪問題的最佳管理模式[13],已被國內(nèi)外學(xué)者廣泛認(rèn)可。相關(guān)研究表明,在一定降雨條件下,最大化地使植被等滲透地表對不透水地表進(jìn)行阻斷的土地利用模式,可使匯水區(qū)總下滲量提高99.3%,總徑流量下降52%,徑流系數(shù)下降52.4%[14];在不同的用地類型下, 植被可直接降低土壤溫度[15],而提升植被覆蓋度所起的降溫效果也有明顯差異,降溫效果從大到小依次為林地、裸地、建設(shè)用地、園地、草地、耕地和水體[16]。可見,分析不同土地利用模式下特定生態(tài)系統(tǒng)服務(wù)供給能力差異特征,有助于深入城鎮(zhèn)土地利用生態(tài)效益研究,提高城鎮(zhèn)土地利用綜合效益。

    鑒于此,本文以北京市三山五園地區(qū)為案例區(qū),選取生態(tài)系統(tǒng)服務(wù)評估中相對成熟的碳儲量評估方法,通過分析不同土地利用規(guī)模、結(jié)構(gòu)、開發(fā)方式和空間布局等土地利用模式下碳儲量的空間分異特征,從碳儲量視角揭示有助于提升生態(tài)系統(tǒng)服務(wù)能力的城鎮(zhèn)土地利用模式,以期為合理、高效利用城鎮(zhèn)土地提供參考。

    2 研究區(qū)與研究方法

    2.1 研究區(qū)概況

    “三山五園”地區(qū)是北京西北郊原皇家行宮苑囿的總稱,位于北京市海淀區(qū)(地理坐標(biāo)為30°31′—30°36′,東經(jīng)114°21′—114°28′E)[17],北至馬連洼街道,南抵四季青鎮(zhèn),東臨清華園街道,西達(dá)香山山脊線,總面積67.27 km2,涵蓋了燕園、萬柳街道、香山、四季青鎮(zhèn)、海淀等10個街道行政單元。該區(qū)域自然生態(tài)環(huán)境本底條件優(yōu)越、文化底蘊深厚,分布有香山靜宜園、萬壽山清漪園(頤和園)、玉泉山靜明園以及附近的圓明園、暢春園等全國重點歷史文化保護(hù)單位。近十幾年來,隨北京市近郊城市化進(jìn)程的加快,特別是周邊中關(guān)村、上地軟件園等重要經(jīng)濟(jì)產(chǎn)業(yè)發(fā)展中心的帶動,三山五園地區(qū)稻田、濕地等逐漸被建設(shè)用地取代,景觀空間格局破碎[17]。但即使如此,仰賴自身生態(tài)基底與歷史文化保護(hù)政策,三山五園地區(qū)依然是北京市西北部大型綠化隔離帶和楔形綠地的必要組成部分。那么,在歷史文化保護(hù)與經(jīng)濟(jì)發(fā)展推動的共軛作用下,以三山五園地區(qū)為案例區(qū)開展城鎮(zhèn)土地利用模式與生態(tài)效益分析,具有一定的代表性。

    2.2 城市生態(tài)系統(tǒng)碳儲量評估模型

    作為全球公認(rèn)的生態(tài)系統(tǒng)服務(wù)模型之一,InVEST模型中碳儲量模塊主要基于土地利用/覆被和對應(yīng)的碳庫參數(shù)評估土地利用碳儲量[18]。4大碳庫的碳密度值(即地上碳、地下碳、土壤碳和死亡有機(jī)碳碳庫)是計算碳儲量的基礎(chǔ),一般可通過地表調(diào)查和遙感估算兩種方法獲取[19]。城市土地利用碳儲量可分解為植被碳和土壤碳[20],基于大量實地調(diào)查獲取相應(yīng)森林、草地等土地利用類型的碳密度值[21-22]。由于城市內(nèi)分布著大量的綠化用地,如公園、居住區(qū)附屬綠地、街旁綠地等,針對這些園林綠地的植被碳密度和土壤碳密度的實地調(diào)查也取得一定研究成果[22-25],為進(jìn)一步優(yōu)化城市土地利用碳儲量評估方法提供了基礎(chǔ)。

    針對城市土地利用的多樣性,本文綜合自然生態(tài)系統(tǒng)、農(nóng)業(yè)生態(tài)系統(tǒng)、園林綠地碳儲量等研究成果修正碳庫參數(shù),提出基于InVEST模型計算城市土地利用/覆被類型碳儲量的修正模型,計算公式如下:

    式(1)中,Ctotal為研究區(qū)的碳儲量總數(shù)(t),Cai和Cbi分別為i類土地利用/覆被類型上植被的地上、地下碳密度值(t/hm2),可通過園林植物[22-25]、草地[22,26-29]、森林[26-28]和濕地[29]碳儲量研究成果獲取植被碳密度(表1);其中,林地地下碳密度值是依據(jù)IPCC(2006)推薦使用的根莖比(0.24)、利用地上碳密度值進(jìn)行換算,園地、公園與綠地地下碳密度值在根據(jù)根莖比和地上碳密度計算的基礎(chǔ)上再乘以0.8的換算系數(shù)[27,29-30]。Csi為i類土地利用/覆被類型上的土壤碳密度值(t/hm2),在本研究中采用北京市西北部典型土地利用類型土壤碳研究成果[23-24](表1);NDVIi為i類土地利用/覆被類型單元的NDVI值,NDVImax為土地利用/覆被單元的NDVI最大值,可基于2014年Landsat 8遙感影像數(shù)據(jù)獲??;Ai為第i類土地利用/覆被類型的面積。

    在本研究中,研究區(qū)土地利用與覆被類型可根據(jù)全國土地利用現(xiàn)狀分類標(biāo)準(zhǔn)(GB/T 21010-2007),在研究區(qū)土地利用變更調(diào)查數(shù)據(jù)的基礎(chǔ)上,基于2014年航空照片數(shù)據(jù)、輔以城市建設(shè)用地現(xiàn)狀圖,利用目視解譯方法劃分為6類,包括林地、園地、公園與綠地、水域、裸地和建設(shè)用地(表1)。

    此外,2014年Landsat 8遙感影像、航空照片數(shù)據(jù)和城市建設(shè)用地現(xiàn)狀圖,均根據(jù)土地利用變更調(diào)查數(shù)據(jù)進(jìn)行空間配準(zhǔn),以滿足后續(xù)分析使用。

    表1 土地利用類型的植物地上碳和土壤碳密度值Tab.1 The per-unit carbon storage of plants and soils for the study area

    2.3 城市土地利用模式的度量方法:基于歸一化不透水地表指數(shù)(NDISI)

    目前,基于多時段、中空間分辨率遙感影像數(shù)據(jù)(如Landsat),監(jiān)測土地利用變化是宏觀尺度上開展城市土地利用模式研究的主要方法[8-12]。但在中小尺度上,如城市建成區(qū),僅僅通過土地利用類型變化很難反映建筑密度、綠化狀況及空間配置格局等城鎮(zhèn)土地利用開發(fā)建設(shè)模式的差別。基于城市地表多時段遙感數(shù)據(jù)研究發(fā)現(xiàn),城市化進(jìn)程中植被、不透水表面與土壤三種組分此消彼長,基于混合光譜VIS模型(Vegetation-Impervious Surface-Soil,簡稱VIS模型)提取不透水地表指數(shù)可作為城市土地利用變化特征指標(biāo)之一[32-33]。不透水地表指數(shù)可以度量諸如建筑物頂部、人行道、公路、廣場、機(jī)場、停車場等多種類型的人工建筑物地表,是城市建設(shè)程度的重要指標(biāo)之一[34-35]。

    據(jù)此,本研究基于Landsat 5影像和Landsat 8影像分別提取1992年和2014年的歸一化差值不透水面指數(shù)(NDISI),并計算NDISI變化率,并將其作為度量城市建設(shè)強(qiáng)度的基本指標(biāo)。NDISI及其變化率計算公式如下:

    式(2)中,CRI為城市化前后NDISI值的變化率,N1、N2分別為快速城市化前期和后期的NDISI值,NDISI指數(shù)值可通過Landsat 8遙感影像計算,如式(3)、式(4)[36]。

    式(3)—式(4)中,NDISI為歸一化差值不透水面指數(shù),此值越高,代表城市建設(shè)強(qiáng)度越高;MNDWI為改進(jìn)型歸一化水體指數(shù),TIR、NIR、MIR1和Green分別為Landsat 8影像的熱紅外、近紅外、中紅外1和綠光波段的灰度值。

    在此基礎(chǔ)上,以街道為單位,分別從土地利用結(jié)構(gòu)與規(guī)模、開發(fā)方式與布局、區(qū)位條件等多個層面,提取各街道的NDISI變化率、土地利用/覆被構(gòu)成與格局以及區(qū)位等因素,綜合度量城鎮(zhèn)土地利用模式(表2)。其中,NDISI和NDVI空間聚集度表征土地利用/覆被構(gòu)成的空間格局特征,可采用地理信息系統(tǒng)空間聚類檢驗中的局部莫蘭指數(shù)(Local Moran Index),采用HH-LH區(qū)域占所在區(qū)域的面積比例來表征,計算公式如下:

    表2 城鎮(zhèn)土地利用模式表征指標(biāo)Tab.2 Index of urban land use pattern for the study area

    式(5)中,Ii為空間單元i的局部莫蘭指數(shù);wij是空間權(quán)重,表示空間單元i、j的鄰近關(guān)系;xi、xj為空間單元i、j的NDISI值或NDVI值,是NDISI值或NDVI值的平均值。

    根據(jù)研究區(qū)特征,區(qū)位條件以歷史文化保護(hù)和產(chǎn)業(yè)發(fā)展中心影響程度表征,可通過式(6)計算:

    式(6)中,Di是歷史文化保護(hù)中心或產(chǎn)業(yè)發(fā)展中心的影響值,wi是影響權(quán)重,歷史文化保護(hù)中心的wi可根據(jù)《北京歷史文化名城保護(hù)規(guī)劃》中劃定的世界文化遺產(chǎn)、國家級/市級/區(qū)級文物保護(hù)單位進(jìn)行量化賦值,產(chǎn)業(yè)發(fā)展中心的wi因缺少產(chǎn)業(yè)等級類型數(shù)據(jù)、暫以等權(quán)重處理,即∑wi= 1;bi是與歷史文化保護(hù)中心或產(chǎn)業(yè)發(fā)展中心的緩沖區(qū)距離,根據(jù)影響距離衰減規(guī)律進(jìn)行量化賦值。

    2.4 城鎮(zhèn)土地碳儲量與土地利用模式的相關(guān)分析

    聚類統(tǒng)計作為常用多元統(tǒng)計學(xué)方法之一,可根據(jù)樣本多個指標(biāo)的相似性和差異性,對樣本進(jìn)行特征性分類。本文通過SPSS統(tǒng)計軟件聚類分析城鎮(zhèn)土地碳儲量與土地利用模式指標(biāo)的相關(guān)關(guān)系,明確不同土地利用模式下碳儲能力差異,進(jìn)而揭示影響城鎮(zhèn)土地碳儲能力的土地利用關(guān)鍵因素。

    3 結(jié)果分析

    3.1 碳儲量空間分布特征

    根據(jù)InVEST模型碳儲量模塊修正模型,經(jīng)原始數(shù)據(jù)處理、計算,2014年三山五園地區(qū)土地利用總碳儲量為23.31×104t,單位土地面積平均碳儲量為32.28 t/hm2。全區(qū)土地碳儲量在空間上整體呈現(xiàn)西高—東低、中間高—南北低的梯度分布格局。

    具體來看,單位土地面積平均碳儲量高于全區(qū)平均水平的包括香山街道、萬柳街道(西區(qū))、四季青鎮(zhèn)、萬柳街道(東區(qū))和青龍橋街道(圓明園區(qū));其中,香山街道單位土地面積平均碳儲量值最高,約為47.41 t/hm2。其余8個街道的單位土地面積平均碳儲量均低于全區(qū)平均水平,其中青龍橋街道(頤和園區(qū))的單位面積土地碳儲量值約31.37 t/hm2,與平均水平較為接近;上地街道和青龍橋街道(北區(qū))單位面積土地碳儲量值約27 t/hm2,而清華園、燕園、海淀、曙光和馬連洼街道單位面積土地碳儲量值均低于26 t/hm2,尤其是馬連洼街道僅為22 t/hm2。

    3.2 城鎮(zhèn)土地利用模式與碳儲量的聚類分析

    根據(jù)單位土地面積碳儲量與其土地利用模式指標(biāo)的Pearson相關(guān)性(表3),通過層次聚類分析將13個街道單元劃分為高碳匯、中碳匯和低碳匯三類(圖1)。

    表3 單位土地面積碳儲量與土地利用模式指標(biāo)的Pearson相關(guān)性Tab.3 The pearson correlation coeffcient between per unit carbon-storage and land use pattern index

    3.2.1 高碳匯—城鎮(zhèn)土地利用模式特征 萬柳街道(西區(qū))與四季青鎮(zhèn)、青龍橋街道(圓明園)和香山街道、以及萬柳街道(東區(qū))、青龍橋街道(頤和園)6大街道,其內(nèi)單位土地面積碳儲量平均值高達(dá)41 t/hm2以上(表4),故稱之為高碳匯—城鎮(zhèn)土地利用模式類型。

    圖1 三山五園地區(qū)土地利用模式—碳儲能力聚類結(jié)果(a)及其主要特征(b)Fig.1 The characteristics of land use pattern with different-level carbon storages in the Three-hills and Five-gardens Area

    表4 第一類街道(高碳匯)土地利用模式與碳儲量指標(biāo)Tab.4 The land use pattern index of the sub-districts with high carbon-storage in the Three-hills and Five-gardens Area

    (1)此類地區(qū)歷史文化保護(hù)影響度平均值(9.08)遠(yuǎn)遠(yuǎn)高于全區(qū)平均水平(6.15),受產(chǎn)業(yè)經(jīng)濟(jì)發(fā)展中心的影響度(6.05)則明顯低于全區(qū)平均水平(10.88);這說明歷史文化保護(hù)政策不僅有利于提高歷史文化遺產(chǎn)保護(hù)與社會文化服務(wù)能力,而且有益于維持城鎮(zhèn)土地碳匯能力。

    (2)受歷史文化保護(hù)的高度庇護(hù),盡管產(chǎn)業(yè)經(jīng)濟(jì)發(fā)展影響程度差異較大,但此類區(qū)域建設(shè)用地總量均控制在65%以下,低于全區(qū)平均水平(68.3%);與之共軛,綠地與水域面積比合計均在35%以上,高于全區(qū)平均水平(31.7%);這說明高碳匯—城鎮(zhèn)土地利用模式中,綠地與建設(shè)用地面積比應(yīng)大于1∶2為佳。

    (3)NDISI空間聚集度平均僅為1.41,最高值僅為2.34(青龍橋街道—圓明園區(qū)),遠(yuǎn)遠(yuǎn)低于全區(qū)平均水平(5.04),這說明高碳匯—城鎮(zhèn)土地利用模式呈現(xiàn)以不透水地表為主的開發(fā)建設(shè)密度低、且相對緊湊的空間形態(tài)。

    (4)1992—2014年NDISI變化率平均值(5.54)低于全區(qū)平均水平(6.33),土地開發(fā)建設(shè)強(qiáng)度相對較低,但空間分異明顯。其中,萬柳街道(東、西區(qū))和四季青鎮(zhèn)的NDISI變化率較高,均在9.0左右,香山街道的NDISI變化率居中(約5.44),4者單位土地面積碳儲量值均高于40 t/hm2;而青龍橋街道(頤和園區(qū)、圓明園區(qū))NDISI變化率為-3.38、3.45,建設(shè)程度明顯低于前四者,但單位土地面積碳儲量卻低于前4者,約31—37 t/hm2。對比NDISI和NDVI聚集度可知,香山街道、萬柳街道(東、西區(qū))和青龍橋街道(頤和園區(qū)、圓明園區(qū))的NDISI空間聚集度較低、NDVI空間聚集度較高,而四季青鎮(zhèn)二者均較低。這說明在綠地面積相似的前提下,較高強(qiáng)度的開發(fā)建設(shè)活動并非會降低土地利用碳儲能力,這與其土地利用開發(fā)的空間格局密切相關(guān):緊湊的建設(shè)用地布局(以不透水地表為表征)配置足夠數(shù)量、相對聚集的綠地,有利于提高城市土地利用碳匯能力。

    3.2.2 中碳匯—城鎮(zhèn)土地利用模式特征 清華園街道、燕園街道和青龍橋街道(北區(qū))等三個街道,單位土地面積碳儲量平均值約25.33 t/hm2,與其他區(qū)域相比其碳儲能力居中,故稱之為中碳匯—城鎮(zhèn)土地利用模式。

    (1)根據(jù)1992—2014年NDISI變化率,清華園街道開發(fā)建設(shè)強(qiáng)度最高(7.71),建設(shè)用地面積高達(dá)90%以上,燕園和青龍橋街道(北區(qū))開發(fā)建設(shè)強(qiáng)度相當(dāng)且NDISI變化率較低,建設(shè)用地面積約80%左右。

    (2)對應(yīng)地理區(qū)域和用地結(jié)構(gòu)來看,清華園和燕園街道受產(chǎn)業(yè)經(jīng)濟(jì)發(fā)展影響較大(影響度>16),同時,燕園街道明顯受到歷史文化遺產(chǎn)保護(hù)政策的庇護(hù)(影響度=9.2);而產(chǎn)業(yè)經(jīng)濟(jì)發(fā)展對青龍橋街道(北區(qū))影響較?。?.4),歷史文化保護(hù)影響度居中(4.8)。仰賴歷史文化保護(hù)的庇護(hù),燕園街道和青龍橋街道(北區(qū))綠地面積相對較高(約15%),而清華園街道綠地面積不足10%。

    (3)綠地面積高的燕園街道單位土地面積碳儲量(23.47),卻低于青龍橋街道(北區(qū))和清華園街道。對比NDISI和NDVI空間聚集度可知,青龍橋街道(北區(qū))綠地面積比雖低于燕園街道,但其NDISI空間聚集度極低(0.22),但燕園街道卻高達(dá)10.43;而清華園街道綠地面積比僅為7.14%,但由于NDISI和NDVI空間聚集度相對較低(分別為5.46和0.8),其單位土地面積碳儲量(25.34)略高于燕園街道(表5)。這說明在綠地面積較低(10%左右)、建設(shè)面積較高(約80%及以上)的情況下,分散式的綠地布局有助于提高城鎮(zhèn)土地利用碳儲能力。

    3.2.3 低碳匯—城鎮(zhèn)土地利用模式特征 馬連洼街道、曙光街道、上地街道和海淀街道4個街道,其單位土地面積碳儲量平均值為24.21 t/hm2,在三山五園地區(qū)最低,故稱之為低碳匯—城鎮(zhèn)土地利用模式(表6)。

    (1)1992—2014年NDISI變化率平均約9.12%,城市開發(fā)建設(shè)強(qiáng)度明顯高于前兩類地區(qū);其中,馬連洼和上地街道NDISI變化率高達(dá)11%以上,說明馬連洼和上地街道是三山五園地區(qū)1992—2014年間開發(fā)建設(shè)強(qiáng)度最高的區(qū)域;而曙光和海淀街道開發(fā)建設(shè)強(qiáng)度相對較低(NDISI變化率約5%),與香山街道持平。

    (2)在開發(fā)建設(shè)強(qiáng)度高的馬連洼和上地街道中,上地街道由于產(chǎn)業(yè)經(jīng)濟(jì)發(fā)展水平較高(特別是中關(guān)村軟件園等高新技術(shù)園區(qū)),且綠地和水域面積比例較高(約20%左右,與燕園街道的相近),其單位土地面積碳儲量明顯較馬連洼街道高(約27.09 t/hm2);而馬連洼街道建設(shè)用地面積接近90%,綠地面積僅為6.16%,故其單位土地面積碳儲量僅為21.98 t/hm2。

    表5 第二類街道(中碳匯)土地利用模式與碳儲量指標(biāo)Tab.5 The land use pattern index of the sub-districts with medium carbon-storage in the Three-hills and Five-gardens Area

    表6 第三類街道(低碳匯)土地利用模式與碳儲量指標(biāo)Tab.6 The land use pattern index of the sub-districts with low carbon-storage in the Three-hills and Five-gardens Area

    (3)在曙光和海淀街道中,盡管1992—2014年開發(fā)建設(shè)強(qiáng)度較低,但二者NDISI空間聚集度均高于10,特別是海淀街道高達(dá)28.14,說明開發(fā)建設(shè)密度已很高;受原建設(shè)用地格局的影響,曙光街道綠地率約12.31%,而海淀街道綠地率卻不足5%(盡管受歷史文化保護(hù)庇護(hù)和產(chǎn)業(yè)經(jīng)濟(jì)發(fā)展推動),因此海淀街道單位土地面積碳儲量(23.66 t/hm2)依舊低于曙光街道(24.1 t/hm2)。

    (4)在綠地面積比<10%的海淀和馬連洼街道中,受歷史文化保護(hù)政策庇護(hù)和經(jīng)濟(jì)發(fā)展支持的雙重影響,海淀街道NDISI變化率較小,且NDVI和NDISI空間聚集度遠(yuǎn)高于馬連洼街道,其單位土地面積碳儲量高于馬連洼街道;對于綠地面積比>10%的曙光和上地街道,NDVI空間聚集度均很低,但曙光街道NDISI空間聚集度遠(yuǎn)高于上地街道,故其單位土地面積碳儲量低于上地街道。這也說明在綠地面積比較低且開發(fā)強(qiáng)度較大(建設(shè)用地面積較高)時,分散式的綠地與建筑布局模式有助于提高城鎮(zhèn)土地利用碳儲能力。

    4 結(jié)論與討論

    本文以碳儲量作為土地利用生態(tài)效益切入點,通過不透水地表變化率為基礎(chǔ)的土地利用模式指標(biāo)來揭示城鎮(zhèn)不同土地利用模式下其碳儲量的空間差異,進(jìn)而為優(yōu)化城鎮(zhèn)土地利用結(jié)構(gòu)、開發(fā)建設(shè)模式及空間配置提供參考。

    (1)2014年,三山五園地區(qū)土地利用碳儲總量為23.31×104t,單位土地面積平均碳儲量32.28 t/hm2;各街道土地利用碳儲量水平差異顯著,空間上呈現(xiàn)西高—東低、中間高—南北低的梯度分布格局。通過街道土地利用特征與碳儲量水平的聚類分析,可將全區(qū)土地利用模式劃分為高、中、低碳匯三類區(qū)域。

    (2)以不透水地表變化率表征土地開發(fā)建設(shè)強(qiáng)度,本研究發(fā)現(xiàn)城市開發(fā)建設(shè)活動并非一定降低土地利用碳匯能力,這與土地所在區(qū)位條件、建設(shè)類型以及開發(fā)建設(shè)模式等因素密切相關(guān):①在產(chǎn)業(yè)/經(jīng)濟(jì)發(fā)展推動城鄉(xiāng)建設(shè)的背景下,歷史文化遺產(chǎn)保護(hù)政策有利于控制城市建設(shè)用地擴(kuò)張,引導(dǎo)城市建設(shè)用地供給與規(guī)劃,保護(hù)和維續(xù)現(xiàn)存城鎮(zhèn)生態(tài)用地,進(jìn)而提高城鎮(zhèn)土地的碳匯能力。②在綠地面積比較高且數(shù)量相當(dāng)?shù)那闆r下,采用緊湊且分散式的城市開發(fā)建設(shè)模式有助于提高土地利用碳儲能力。③在綠地面積較低、高密度建設(shè)的土地利用方式下,提高建筑物立體綠化、增加園林植物種植密度,特別是種植固碳能力強(qiáng)的樹種及地被植物更有利于提高土地碳儲能力。

    (3)未來研究中可加強(qiáng)城市土地利用模式及其生態(tài)效益關(guān)系動態(tài)研究,特別是進(jìn)一步分析不同城鎮(zhèn)土地利用類型(如住宅用地、公共服務(wù)與管理用地等)及其開發(fā)建設(shè)模式對碳儲量等生態(tài)系統(tǒng)服務(wù)能力的影響,將有助于制定合理的城鎮(zhèn)土地利用供給結(jié)構(gòu)、開發(fā)建設(shè)模式以及生態(tài)環(huán)境建設(shè)策略,進(jìn)而提高城鎮(zhèn)土地利用的綜合效益。

    (References):

    [1] 歐陽志云,王橋,鄭華,等. 全國生態(tài)環(huán)境十年變化(2000—2010年)遙感調(diào)查評估[J] . 中國科學(xué)院院刊,2014(,4):462 - 466.

    [2] 趙國松,劉紀(jì)遠(yuǎn),匡文慧,等. 1990—2010年中國土地利用變化對生物多樣性保護(hù)重點區(qū)域的擾動[J] . 地理學(xué)報,2014,69(11):1640 - 1650.

    [3] 胡蒙蒙,張軍民,徐麗萍,等. 基于生態(tài)綠當(dāng)量的瑪納斯河流域土地利用生態(tài)效益研究[J] . 干旱區(qū)研究,2016,3(35):996 -1002.

    [4] 梁紅梅,劉衛(wèi)東,劉會平,等. 土地利用社會經(jīng)濟(jì)效益與生態(tài)環(huán)境效益的耦合關(guān)系——以深圳市和寧波市為例[J] . 2008,22(2):42 - 48.

    [5] 吳金華,李紀(jì)偉,梁晶晶. 土地利用程度與效益關(guān)系研究——以延安市為例[J] . 中國土地科學(xué),2011,2(58):54 - 60.

    [6] 葉敏婷,王仰麟,彭建,等. 深圳市土地利用效益變化及其區(qū)域分異[J] . 資源科學(xué),2008,2(03):401 - 408.

    [7] 王國剛,劉彥隨,方方. 環(huán)渤海地區(qū)土地利用效益綜合測度及空間分異[J] . 地理科學(xué)進(jìn)展,2013,3(24):649 - 656.

    [8] Yan Y., Kuang W., Zhang C., et al. Impacts of impervious surface expansion on soil organic carbon-a spatially explicit study[J] . Science Report 5, 17905, doi: 10.1038/srep1790(52015).

    [9] Edmondson J., Davies Z., McHugh N., et al. Organic carbon hidden in urban ecosystems[J] . Science Report 2, 963, doi: 10.1038/ srep0096(32012).

    [10] Sofia G., Roder G.,Dalla Fontana G., Tarolli P. Flood dynamics in urbanised landscapes: 100 years of climate and humans’ interaction[J] . Science Report 7, 40527, doi: 10.1038/srep4052(72017).

    [11] Mao W., Wang X., Cai J., et al. Multidimensional histogram-based information capacity analysis of urban heat island effect using Landsat 8 data[J] . Remote Sensing Letters,2016,1(07):925 - 934.

    [12] Weng Q., Lu D., Schubring J. Estimation of land surface temperature and vegetation abundance relationship for urban heat island studies[J] . Remote Sensing of Environment,2004,89:467 - 483.

    [13] 賈海峰,姚海蓉,唐穎,等. 城市降雨徑流控制LID BMPs 規(guī)劃方法及案例[J] . 水科學(xué)進(jìn)展,2014,2(52):260 - 267.

    [14] 班玉龍,孔繁花,尹海偉,等. 土地利用格局對SWMM模型匯流模式選擇及相應(yīng)產(chǎn)流特征的影響[J] . 生態(tài)學(xué)報,2016,36(14):4317 - 4326.

    [15] Edmondson J., Stott I., Davies Z., et al. Soil surface temperatures reveal moderation of the urban heat island effect by trees and shrubs[J] .Science Report 6, 33708, doi: 10.1038/srep3370(82016).

    [16] 張小飛,王仰麟,吳健生,等. 城市地域地表溫度與植被覆蓋定量關(guān)系分析——以深圳市為例[J] . 地理研究,2006,25(3):369 - 377.

    [17] 劉劍,胡立輝,李樹華. 北京“三山五園”地區(qū)景觀歷史性變遷分析[J] . 中國園林,2011(,2):54 - 58.

    [18] Kareiva P., Tallis H., Ricketts T. H., et al. Natural capital: theory and practice of mapping ecosystem services[J] . New York: Oxford University Press,2011.

    [19] 史琰,葛瀅,金荷仙. 城市植被碳固存研究進(jìn)展[J] . 林業(yè)科學(xué),2016,5(26):122 - 129.

    [20] 吳佩君,劉小平,黎夏,等. 基于InVEST 模型和元胞自動機(jī)的城市擴(kuò)張對陸地生態(tài)系統(tǒng)碳儲量影響評估——以廣東省為例[J] . 地理與地理信息科學(xué),2016,3(25):22 - 28.

    [21] Timilsina N., Escobedo F. J, Staudhammer C. L., et al. Analyzing the causal factors of carbon stores in a subtropical urban fores[tJ] . Ecological Complexity,2014,20:23 - 32.

    [22] 王迪生. 基于生物量計測的北京城區(qū)園林綠地凈碳儲量研究[D] . 北京:北京林業(yè)大學(xué),2010.

    [23] 羅上華,毛齊正,馬克明,等. 北京城市綠地表層土壤碳氮分布特征[J] . 生態(tài)學(xué)報,2014,3(420):6011 - 6019.

    [24] 張廷龍,孫睿,胡波,等. 北京西北部典型城市化地區(qū)不同土地利用類型土壤碳特征分析[J] . 北京師范大學(xué)學(xué)報(自然科學(xué)版),2010(,1):97 - 102.

    [25] 孫榮凱,陳智慧. 北京溫榆河流域周邊河岸帶景觀結(jié)構(gòu)與土壤碳儲量響應(yīng)[J] . 安徽農(nóng)業(yè)科學(xué),2016,4(413):207 - 210.

    [26] 張萍. 北京森林碳儲量研究[D] . 北京:北京林業(yè)大學(xué),2009.

    [27] 王光華. 北京森林植被固碳能力研究[D] . 北京:北京林業(yè)大學(xué),2012.

    [28] 楊芝歌,周彬,余新曉,等. 北京山區(qū)生物多樣性分析與碳儲量評估[J] . 水土保持通報,2012,3(23):42 - 46.

    [29] 奚小環(huán),楊忠芳,崔玉軍,等. 東北平原土壤有機(jī)碳分布與變化趨勢研究[J] . 地學(xué)前緣,2010(,3):213 - 221.

    [30] Mchale M., Burke I., Lefsky M., et al. Urban forest biomass estimates: is it important to use allometric relationships developed specifically for urban trees?[J] . Urban Ecosystems,2009,1(21):95 - 113.

    [31] Nowak D. Atmospheric carbon reduction by urban trees[J] . Journal of Environmental Management,1993,3(73):207 - 217.

    [32] Ridd M. K. Exploring a V-I-S( Vegetation-Impervious Surface-Soil) model for urban ecosystem analysis through remote sensing: comparative anatomy for cities[J] . International Journal of Remote Sensing,1995,1(612):2165 - 2185.

    [33] 劉珍環(huán),王仰麟,彭建. 不透水表面遙感監(jiān)測及其應(yīng)用研究進(jìn)展[J] . 地理科學(xué)進(jìn)展,2010,(9):1143 - 1152.

    [34] Arnold C., Gbbons C. Impervious surface: the emergence of a key environmental indicato[rJ] . Journal of the American Planning Association,1996,6(22):243 - 258.

    [35] 王婷,趙玲,武堅,等. 改進(jìn)型歸一化差值不透水表面指數(shù)探討——以蘭州市建成區(qū)為例[J] . 測繪與空間地理信息,2015,(3):182 - 185.

    [36] 徐涵秋. 城市不透水面與相關(guān)城市生態(tài)要素關(guān)系的定量分析[J] . 生態(tài)學(xué)報,2009(,5):2456 - 2462.

    (本文責(zé)編:王慶日)

    Urban Land Use Pattern and Its Ecological Benefit from the Perspective of Carbon Storage

    GUO Xu-dong1, CHANG Qing2, LIU Xiao2, BAO Hui-min2, SUN Yi2, LV Chun-yan1, ZHANG Yan-yu1
    (1. Key Laboratory of Land Use, Ministry of Land and Resources, Beijing 100035, China; 2. Department of Ornamental Horticulture and Landscape Architecture, China Agricultural University, Beijing 100193, China)

    Taken the carbon storage as the ecological benefit index, this paper aims to analyze the spatial differences of carbon storages among various land use pattern in the Three-hills and Five-gardens Area of Haidian District, Beijing City. Methods applied in this paper include the InVEST Carbon model, land use pattern index based on the Normalized Impervious Surface Index (NDISI) and hierarchical cluster analysis. The results show that the carbon stored in the Three-hills and Five-gardens Area in 2014 was 233.1 thousand tons in total, and the average carbon storage of land is 32.28 t/hm2. There are significant spatial differences of carbon storages among the sub-districts; for example, all of the average carbon storage of land in the sub-district of Xiangshan, Wanliu and Sijiqing are higher than 41 t/hm2, while lower than 28 t/hm2in the sub-district of Malianwa, Shuguang and Shangdi. Comparing the carbon storage with the landutilization intensity quantified by the change rate index of NDISI, this paper found that the increasing of land use intensity did not necessarily lead to the reduction of carbon storage; the carbon storage in a given area was dependent on its location, land use intensity and structure, as well as development modes. The compact construction might provide enough green spaces, which has the potential capability of higher carbon storage, such as Wanliu and Sijiqing. Therefore, from the aspect of carbon storage, not only land use intensity but also land use structure has highly influence on the ecological benefit of urban land use. In conclusion, it is supposed that studies about urban land use pattern and its ecological benefit need to be enhanced to optimize the comprehensive benefit of land use in the future.

    land use; carbon storage; ecological benefit; land development pattern; Three-hills and Five-gardens Area

    F301.2

    A

    1001-8158(2017)04-0061-10

    10.11994/zgtdkx.20170424.145441

    2016-12-19;

    2017-03-13

    國土資源部土地利用重點實驗室2016年開放基金(KLLU201602)。

    郭旭東(1972-)男,遼寧鐵嶺人,研究員。主要研究方向為土地利用變化和土地生態(tài)。E-mail: guoxudong@mail.clspi.org.cn

    常青(1978-),女,內(nèi)蒙古豐鎮(zhèn)人,副教授。主要研究方向城市生態(tài)與景觀規(guī)劃。E-mail: changqing@cau.edu.cn

    猜你喜歡
    儲量城鎮(zhèn)土地利用
    《礦產(chǎn)資源儲量技術(shù)標(biāo)準(zhǔn)》修訂對資源儲量報告編寫的影響
    2.5 MPa及以上城鎮(zhèn)燃?xì)夤艿琅c輸氣管道區(qū)別
    煤氣與熱力(2021年3期)2021-06-09 06:16:28
    基于三維軟件資源儲量估算對比研究
    文化邊城鎮(zhèn)遠(yuǎn)
    中華民居(2020年6期)2020-06-09 05:47:48
    土地利用生態(tài)系統(tǒng)服務(wù)研究進(jìn)展及啟示
    濱海縣土地利用挖潛方向在哪里
    城鎮(zhèn)排水系統(tǒng)量化指標(biāo)體系研究
    城鎮(zhèn)醫(yī)保支出為何跑贏消費支出
    概率統(tǒng)計法在儲量估算中的應(yīng)用
    斷塊油氣田(2014年5期)2014-03-11 15:33:45
    本月起實施頁巖氣儲量行業(yè)標(biāo)準(zhǔn)
    黄色日韩在线| 国产美女午夜福利| 亚洲熟妇熟女久久| 午夜福利欧美成人| av欧美777| 久久热精品热| 免费看a级黄色片| 91麻豆精品激情在线观看国产| 三级男女做爰猛烈吃奶摸视频| 在线播放国产精品三级| 免费观看精品视频网站| 99久久成人亚洲精品观看| 午夜免费激情av| 日本一二三区视频观看| 久久天躁狠狠躁夜夜2o2o| 国产精品精品国产色婷婷| 久久6这里有精品| 嫩草影院入口| 噜噜噜噜噜久久久久久91| 两个人视频免费观看高清| 日韩av在线大香蕉| 精品一区二区免费观看| 少妇丰满av| 亚洲精品日韩av片在线观看| 国产精品综合久久久久久久免费| 好男人电影高清在线观看| 亚洲第一区二区三区不卡| 午夜精品一区二区三区免费看| 18+在线观看网站| 日韩国内少妇激情av| 日本熟妇午夜| 特大巨黑吊av在线直播| 日韩欧美在线二视频| 直男gayav资源| 亚洲成av人片免费观看| 午夜激情欧美在线| 97超级碰碰碰精品色视频在线观看| 亚洲成人久久爱视频| 久久午夜亚洲精品久久| 毛片一级片免费看久久久久 | 国产亚洲精品久久久久久毛片| 成人国产综合亚洲| 日本熟妇午夜| 免费人成在线观看视频色| 免费av不卡在线播放| 欧美最新免费一区二区三区 | 亚州av有码| 欧美性感艳星| 色吧在线观看| 亚洲第一欧美日韩一区二区三区| 亚洲人成电影免费在线| 国内精品久久久久久久电影| 51国产日韩欧美| 成人鲁丝片一二三区免费| 精品午夜福利视频在线观看一区| 国产精品久久久久久久久免 | 男人舔女人下体高潮全视频| 午夜福利18| 国产一区二区三区视频了| 黄色一级大片看看| 好看av亚洲va欧美ⅴa在| 久久人妻av系列| 国产成人a区在线观看| 国产精品电影一区二区三区| 熟妇人妻久久中文字幕3abv| 国产高清三级在线| 桃色一区二区三区在线观看| 狂野欧美白嫩少妇大欣赏| 日韩 亚洲 欧美在线| 少妇人妻精品综合一区二区 | 老司机午夜十八禁免费视频| 很黄的视频免费| 少妇人妻精品综合一区二区 | 国模一区二区三区四区视频| 欧美成狂野欧美在线观看| 久久草成人影院| 亚洲久久久久久中文字幕| 精品无人区乱码1区二区| 国产精品久久电影中文字幕| 国产色爽女视频免费观看| 国产精品永久免费网站| 欧美极品一区二区三区四区| 国产成人欧美在线观看| 久久伊人香网站| 人人妻人人看人人澡| 我的老师免费观看完整版| 18禁黄网站禁片免费观看直播| 天堂√8在线中文| 亚洲最大成人手机在线| 亚洲成a人片在线一区二区| 一级av片app| 精品无人区乱码1区二区| 午夜精品久久久久久毛片777| 欧美成人免费av一区二区三区| 亚洲第一欧美日韩一区二区三区| 日韩大尺度精品在线看网址| 国产免费男女视频| 精品一区二区三区人妻视频| 午夜a级毛片| 亚洲精品粉嫩美女一区| 中亚洲国语对白在线视频| 桃色一区二区三区在线观看| 国产在线精品亚洲第一网站| 国产免费男女视频| 国产高潮美女av| 此物有八面人人有两片| 午夜亚洲福利在线播放| 亚洲天堂国产精品一区在线| 欧美+亚洲+日韩+国产| 嫩草影院精品99| 国产亚洲欧美98| 亚洲va日本ⅴa欧美va伊人久久| avwww免费| 嫁个100分男人电影在线观看| 亚洲av成人精品一区久久| 成人av一区二区三区在线看| 一个人免费在线观看的高清视频| 国产成人av教育| 国产色爽女视频免费观看| 亚洲人与动物交配视频| 亚洲精品粉嫩美女一区| 中文字幕熟女人妻在线| 精品一区二区免费观看| 综合色av麻豆| 嫩草影院新地址| 国产精品女同一区二区软件 | 久久久久久久久久黄片| 久久久久久大精品| 国模一区二区三区四区视频| 国产中年淑女户外野战色| 国产精品永久免费网站| 好看av亚洲va欧美ⅴa在| 国产高清视频在线播放一区| 亚洲美女视频黄频| 午夜亚洲福利在线播放| 日韩有码中文字幕| 国内久久婷婷六月综合欲色啪| 亚洲午夜理论影院| 日韩中文字幕欧美一区二区| 嫩草影院精品99| 欧美另类亚洲清纯唯美| 91字幕亚洲| 51午夜福利影视在线观看| 日本在线视频免费播放| 欧美+亚洲+日韩+国产| 在线观看美女被高潮喷水网站 | 亚洲精品在线观看二区| 欧美成人一区二区免费高清观看| 久久久久久久久中文| 久久国产乱子免费精品| 俺也久久电影网| 精品久久久久久久人妻蜜臀av| 国产午夜精品久久久久久一区二区三区 | 亚洲在线自拍视频| 神马国产精品三级电影在线观看| 97热精品久久久久久| 脱女人内裤的视频| 国产主播在线观看一区二区| 麻豆国产av国片精品| 亚洲成av人片在线播放无| 97人妻精品一区二区三区麻豆| 啪啪无遮挡十八禁网站| 天堂av国产一区二区熟女人妻| 亚洲乱码一区二区免费版| 啪啪无遮挡十八禁网站| 免费一级毛片在线播放高清视频| 国内揄拍国产精品人妻在线| 亚洲人成伊人成综合网2020| 黄色一级大片看看| 免费高清视频大片| 色av中文字幕| 日本免费一区二区三区高清不卡| 毛片一级片免费看久久久久 | 老司机深夜福利视频在线观看| 国产欧美日韩一区二区三| 简卡轻食公司| 伦理电影大哥的女人| 国产蜜桃级精品一区二区三区| 国产高清有码在线观看视频| 精品久久久久久久人妻蜜臀av| eeuss影院久久| 亚洲欧美日韩高清专用| 亚洲av五月六月丁香网| 亚洲最大成人av| 亚洲久久久久久中文字幕| 国产免费一级a男人的天堂| 日韩欧美免费精品| 日本一本二区三区精品| 在线观看一区二区三区| 久久久国产成人免费| 91狼人影院| 丰满的人妻完整版| bbb黄色大片| 日韩国内少妇激情av| 欧美黑人欧美精品刺激| 非洲黑人性xxxx精品又粗又长| 欧美精品啪啪一区二区三区| 琪琪午夜伦伦电影理论片6080| 亚洲avbb在线观看| 99国产极品粉嫩在线观看| 亚洲精品影视一区二区三区av| 欧美潮喷喷水| 黄片小视频在线播放| 久久久久性生活片| 不卡一级毛片| 综合色av麻豆| 久久久久精品国产欧美久久久| 亚洲精品亚洲一区二区| 毛片女人毛片| 天天一区二区日本电影三级| 成年女人永久免费观看视频| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 日本 欧美在线| 国产激情偷乱视频一区二区| 他把我摸到了高潮在线观看| 成年女人毛片免费观看观看9| 不卡一级毛片| 亚洲精品久久国产高清桃花| 熟妇人妻久久中文字幕3abv| 国产精品一区二区三区四区久久| 日韩欧美在线二视频| 国产三级中文精品| 久久久色成人| 欧美不卡视频在线免费观看| 色吧在线观看| 蜜桃久久精品国产亚洲av| 一进一出抽搐gif免费好疼| 亚洲国产精品成人综合色| 看免费av毛片| 搡女人真爽免费视频火全软件 | 亚洲狠狠婷婷综合久久图片| 一级av片app| 我的老师免费观看完整版| 一区二区三区四区激情视频 | 91久久精品国产一区二区成人| 99热精品在线国产| 美女大奶头视频| 十八禁人妻一区二区| 色5月婷婷丁香| 国产私拍福利视频在线观看| 97超级碰碰碰精品色视频在线观看| 成人特级黄色片久久久久久久| 俄罗斯特黄特色一大片| 欧美极品一区二区三区四区| 国产午夜精品久久久久久一区二区三区 | 亚洲国产高清在线一区二区三| 亚洲av日韩精品久久久久久密| 亚洲第一电影网av| 大型黄色视频在线免费观看| 性插视频无遮挡在线免费观看| 欧美日本视频| 国产精品日韩av在线免费观看| 尤物成人国产欧美一区二区三区| 一区福利在线观看| 久久精品国产99精品国产亚洲性色| 国产一区二区三区视频了| 免费人成视频x8x8入口观看| 国产精品亚洲美女久久久| 久久精品影院6| 国产久久久一区二区三区| 亚洲精品一区av在线观看| 免费观看的影片在线观看| 在线十欧美十亚洲十日本专区| 九色成人免费人妻av| 亚洲精品乱码久久久v下载方式| 神马国产精品三级电影在线观看| 精品人妻熟女av久视频| 久久国产精品影院| 欧美成人a在线观看| 国产午夜精品论理片| 亚洲av成人av| 亚洲人成网站高清观看| 国产精品免费一区二区三区在线| 国产久久久一区二区三区| 亚洲美女黄片视频| 美女 人体艺术 gogo| 欧美最新免费一区二区三区 | 赤兔流量卡办理| 99热只有精品国产| 午夜福利成人在线免费观看| 亚洲自拍偷在线| 国产淫片久久久久久久久 | 亚洲av日韩精品久久久久久密| 极品教师在线视频| 亚洲人成网站在线播放欧美日韩| a级毛片免费高清观看在线播放| 欧美日韩亚洲国产一区二区在线观看| 免费av不卡在线播放| 国产欧美日韩精品亚洲av| 首页视频小说图片口味搜索| 激情在线观看视频在线高清| 亚洲精品在线观看二区| 国产黄色小视频在线观看| 久久草成人影院| 国语自产精品视频在线第100页| 国产主播在线观看一区二区| 少妇的逼好多水| 18禁黄网站禁片午夜丰满| 久久人妻av系列| 激情在线观看视频在线高清| 精品乱码久久久久久99久播| av女优亚洲男人天堂| 91在线观看av| 久久精品国产自在天天线| 色精品久久人妻99蜜桃| 一本精品99久久精品77| 亚洲 国产 在线| 久久亚洲真实| 网址你懂的国产日韩在线| 一夜夜www| 99久久99久久久精品蜜桃| 床上黄色一级片| 久久国产精品人妻蜜桃| 亚洲激情在线av| 色av中文字幕| 一区福利在线观看| 99热这里只有精品一区| 亚洲精品日韩av片在线观看| 中出人妻视频一区二区| 亚洲精品成人久久久久久| 99久久精品国产亚洲精品| 大型黄色视频在线免费观看| 狂野欧美白嫩少妇大欣赏| 1000部很黄的大片| 国产亚洲av嫩草精品影院| 婷婷精品国产亚洲av| 91在线精品国自产拍蜜月| 欧美日韩国产亚洲二区| 俄罗斯特黄特色一大片| 亚洲,欧美,日韩| 露出奶头的视频| 99热精品在线国产| 色哟哟·www| 天美传媒精品一区二区| 熟妇人妻久久中文字幕3abv| 精品免费久久久久久久清纯| 色5月婷婷丁香| 嫩草影院精品99| 中文字幕av在线有码专区| 亚洲av成人av| 99久久精品一区二区三区| 精品午夜福利在线看| 男女那种视频在线观看| 亚州av有码| 久久国产精品人妻蜜桃| 亚洲av美国av| 日韩大尺度精品在线看网址| 精品久久国产蜜桃| 少妇的逼水好多| 国产精品日韩av在线免费观看| 精品国内亚洲2022精品成人| 九色成人免费人妻av| 亚洲三级黄色毛片| 麻豆国产97在线/欧美| 啦啦啦观看免费观看视频高清| 日韩欧美三级三区| 亚洲av成人精品一区久久| 日韩中字成人| 免费黄网站久久成人精品 | 国产成人a区在线观看| 欧美黄色片欧美黄色片| 人妻丰满熟妇av一区二区三区| 国产成年人精品一区二区| 又黄又爽又刺激的免费视频.| 亚洲综合色惰| 亚洲av成人精品一区久久| 美女高潮喷水抽搐中文字幕| 啦啦啦观看免费观看视频高清| 欧美精品国产亚洲| 日日夜夜操网爽| 99在线人妻在线中文字幕| 99热精品在线国产| 午夜福利免费观看在线| 久久久久免费精品人妻一区二区| 超碰av人人做人人爽久久| 国产精品亚洲美女久久久| 在线观看美女被高潮喷水网站 | 综合色av麻豆| 欧美乱色亚洲激情| 最近中文字幕高清免费大全6 | 日本黄色视频三级网站网址| 亚洲最大成人手机在线| 精品人妻偷拍中文字幕| 亚洲欧美激情综合另类| 午夜影院日韩av| 女人被狂操c到高潮| 91av网一区二区| 中出人妻视频一区二区| 午夜福利欧美成人| 午夜久久久久精精品| 亚洲av.av天堂| 国产高清三级在线| 国产一区二区在线av高清观看| 成人一区二区视频在线观看| 嫩草影院精品99| 国产伦一二天堂av在线观看| 最近中文字幕高清免费大全6 | 国产免费一级a男人的天堂| 国产淫片久久久久久久久 | 午夜精品在线福利| 中文字幕av成人在线电影| 99热6这里只有精品| av在线观看视频网站免费| 国产精品亚洲av一区麻豆| 亚洲欧美日韩无卡精品| 亚洲美女搞黄在线观看 | 深夜a级毛片| 色综合欧美亚洲国产小说| 精品国产三级普通话版| 十八禁国产超污无遮挡网站| 国产三级黄色录像| 搡老熟女国产l中国老女人| 国内毛片毛片毛片毛片毛片| 白带黄色成豆腐渣| 亚洲欧美清纯卡通| 中文字幕av成人在线电影| 国产日本99.免费观看| 最好的美女福利视频网| 欧美成人免费av一区二区三区| 久久亚洲精品不卡| 国产又黄又爽又无遮挡在线| 久久精品影院6| 少妇的逼水好多| 亚洲av一区综合| 成人国产综合亚洲| 禁无遮挡网站| 久久欧美精品欧美久久欧美| 特大巨黑吊av在线直播| 国产91精品成人一区二区三区| 国产成人a区在线观看| 亚洲内射少妇av| 午夜日韩欧美国产| 久久6这里有精品| 国产老妇女一区| 麻豆国产av国片精品| 国产高潮美女av| 在线看三级毛片| 日韩中字成人| 真实男女啪啪啪动态图| 国产精品三级大全| 久久精品国产亚洲av涩爱 | 国产高潮美女av| 亚洲欧美激情综合另类| 欧美激情在线99| 男人狂女人下面高潮的视频| 一级a爱片免费观看的视频| 亚洲国产精品久久男人天堂| 91麻豆精品激情在线观看国产| 亚洲三级黄色毛片| 国产精品女同一区二区软件 | 精品一区二区免费观看| 啦啦啦韩国在线观看视频| 嫩草影院精品99| 日韩精品中文字幕看吧| 亚洲最大成人手机在线| 亚洲真实伦在线观看| 蜜桃亚洲精品一区二区三区| 久久久成人免费电影| 国产色婷婷99| 免费大片18禁| 一级黄色大片毛片| av女优亚洲男人天堂| 最好的美女福利视频网| 国内久久婷婷六月综合欲色啪| 十八禁网站免费在线| 国产亚洲精品av在线| 91久久精品国产一区二区成人| 久久精品91蜜桃| 美女 人体艺术 gogo| 欧美在线黄色| 又黄又爽又刺激的免费视频.| a级毛片免费高清观看在线播放| 国产免费男女视频| 高清毛片免费观看视频网站| 久久精品人妻少妇| 精品久久久久久成人av| 久久精品影院6| 国产伦一二天堂av在线观看| 18禁裸乳无遮挡免费网站照片| 日本免费一区二区三区高清不卡| 欧美成人免费av一区二区三区| 久久中文看片网| 国产在线男女| 亚洲av熟女| 中文亚洲av片在线观看爽| 欧美日韩福利视频一区二区| 3wmmmm亚洲av在线观看| 国产高清三级在线| 18禁在线播放成人免费| 欧美区成人在线视频| 18禁黄网站禁片午夜丰满| 亚洲精品成人久久久久久| av在线老鸭窝| 免费看美女性在线毛片视频| 久久香蕉精品热| 少妇人妻精品综合一区二区 | 久久久久久国产a免费观看| 中文字幕久久专区| 制服丝袜大香蕉在线| 在线观看66精品国产| 国产真实乱freesex| 精华霜和精华液先用哪个| 欧美高清性xxxxhd video| 亚洲av免费在线观看| 在线播放无遮挡| 欧美黑人巨大hd| 国产免费一级a男人的天堂| 女人被狂操c到高潮| 我的女老师完整版在线观看| 欧美高清成人免费视频www| 日本 欧美在线| 一进一出抽搐gif免费好疼| 亚洲一区高清亚洲精品| 免费电影在线观看免费观看| 欧美高清成人免费视频www| 日本 欧美在线| 亚洲专区中文字幕在线| 精品久久久久久成人av| 亚洲,欧美精品.| 国产亚洲精品av在线| 成人无遮挡网站| 欧美又色又爽又黄视频| 成熟少妇高潮喷水视频| 亚洲内射少妇av| 国产毛片a区久久久久| 亚洲av免费高清在线观看| 亚洲真实伦在线观看| 一本精品99久久精品77| 欧美性感艳星| 亚洲国产高清在线一区二区三| 简卡轻食公司| 国产精品影院久久| 中出人妻视频一区二区| 成年版毛片免费区| 精品国产三级普通话版| 麻豆成人av在线观看| 亚洲精品粉嫩美女一区| 成人欧美大片| АⅤ资源中文在线天堂| 国产欧美日韩精品一区二区| 国产精品精品国产色婷婷| 国产伦在线观看视频一区| 91久久精品国产一区二区成人| 成人特级黄色片久久久久久久| 欧美日韩福利视频一区二区| 此物有八面人人有两片| 最新在线观看一区二区三区| 精品国产三级普通话版| 狂野欧美白嫩少妇大欣赏| 久久国产乱子伦精品免费另类| 日本在线视频免费播放| 国产国拍精品亚洲av在线观看| 国产在线男女| 18禁在线播放成人免费| 亚洲va日本ⅴa欧美va伊人久久| 欧美成人a在线观看| 1000部很黄的大片| 最新在线观看一区二区三区| 九九久久精品国产亚洲av麻豆| 亚洲人成网站高清观看| 淫秽高清视频在线观看| 国产精品国产高清国产av| 欧美成狂野欧美在线观看| 久久热精品热| 亚洲第一欧美日韩一区二区三区| 一a级毛片在线观看| 国产aⅴ精品一区二区三区波| 亚洲熟妇中文字幕五十中出| 欧美一级a爱片免费观看看| 最近视频中文字幕2019在线8| 免费看光身美女| 性插视频无遮挡在线免费观看| 免费在线观看影片大全网站| 成人鲁丝片一二三区免费| 久久久成人免费电影| a级一级毛片免费在线观看| 嫩草影院新地址| 亚洲成人免费电影在线观看| 狂野欧美白嫩少妇大欣赏| 国产不卡一卡二| 成人美女网站在线观看视频| 噜噜噜噜噜久久久久久91| 一本综合久久免费| 欧美成人一区二区免费高清观看| 亚洲精品在线美女| 国产日本99.免费观看| 久久精品国产亚洲av香蕉五月| 日日夜夜操网爽| 一区二区三区四区激情视频 | 婷婷亚洲欧美| 99在线视频只有这里精品首页| 久久国产乱子伦精品免费另类| 精品久久久久久成人av| 美女高潮喷水抽搐中文字幕| 国产av麻豆久久久久久久| 观看美女的网站| 欧美三级亚洲精品| 成人鲁丝片一二三区免费| 九九久久精品国产亚洲av麻豆| 啦啦啦韩国在线观看视频| 国产高清有码在线观看视频| 午夜福利在线观看免费完整高清在 | 色吧在线观看| 午夜日韩欧美国产| 欧美成人a在线观看| 一进一出抽搐gif免费好疼| 婷婷丁香在线五月| 99热只有精品国产| 亚洲成av人片在线播放无| 久久天躁狠狠躁夜夜2o2o| 国产精品综合久久久久久久免费| 如何舔出高潮| 欧美黑人欧美精品刺激| 免费av观看视频| 亚洲aⅴ乱码一区二区在线播放|