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

    結(jié)合土地利用與克里格插值的區(qū)域土壤有機(jī)碳空間表征方法

    2014-09-21 08:04:53顧成軍史學(xué)正于東升
    水土保持研究 2014年2期
    關(guān)鍵詞:分異克里土壤有機(jī)

    顧成軍, 史學(xué)正, 于東升

    (1.土壤與農(nóng)業(yè)可持續(xù)發(fā)展國家重點(diǎn)實(shí)驗(yàn)室, 中國科學(xué)院南京土壤研究所, 南京 210008;2.中國科學(xué)院大學(xué), 北京 100049; 3.滁州學(xué)院 地理信息與旅游學(xué)院, 安徽 滁州 239000)

    結(jié)合土地利用與克里格插值的區(qū)域土壤有機(jī)碳空間表征方法

    顧成軍1,2,3, 史學(xué)正1, 于東升1

    (1.土壤與農(nóng)業(yè)可持續(xù)發(fā)展國家重點(diǎn)實(shí)驗(yàn)室, 中國科學(xué)院南京土壤研究所, 南京 210008;2.中國科學(xué)院大學(xué), 北京 100049; 3.滁州學(xué)院 地理信息與旅游學(xué)院, 安徽 滁州 239000)

    如何利用有限的土壤采樣點(diǎn)準(zhǔn)確預(yù)測土壤屬性一直是研究的熱點(diǎn)?;诤颖笔∪珖诙瓮寥榔詹閿?shù)據(jù),運(yùn)用方差分析研究了使用土地利用和克里格插值來表征區(qū)域表層土壤有機(jī)碳空間分布的可行性。結(jié)果表明:土地利用是表層土壤有機(jī)碳密度空間分異的一個(gè)影響因子,可以直接使用土地利用來反映表層土壤有機(jī)碳密度在空間上的分布和變化。但單一的土地利用對土壤有機(jī)碳密度空間分異的反映能力還不是很高,而土地利用與普通克里格插值結(jié)合以后能顯著地提高對土壤有機(jī)碳空間分異的解釋水平。因此區(qū)域空間上可以將土地利用與普通克里格插值結(jié)合來表征表層土壤有機(jī)碳密度的空間分布。

    土地利用; 土壤有機(jī)碳; 克里格插值; 空間分布

    準(zhǔn)確表征土壤有機(jī)碳的空間分布是土壤有機(jī)碳研究的基礎(chǔ)。從已有研究來看,土壤有機(jī)碳的空間分布在空間尺度較小且地形條件相對簡單的區(qū)域可以使用基于地統(tǒng)計(jì)分析的克里格插值進(jìn)行預(yù)測[1],而在空間尺度較大且環(huán)境條件復(fù)雜的區(qū)域多用土壤有機(jī)碳的相關(guān)因子來反映[2]。為進(jìn)一步提高預(yù)測精度,一些學(xué)者還通過將土壤有機(jī)碳相關(guān)因子與克里格插值結(jié)合的方法來表征土壤有機(jī)碳的空間分布[3]。

    區(qū)域尺度上不同土地利用的土壤有機(jī)碳含量對比研究很多,大多數(shù)研究結(jié)果都表明不同土地利用的土壤有機(jī)碳存在顯著差異[4-7]。因此土地利用有時(shí)候能夠反映出土壤有機(jī)碳的空間差異,是區(qū)域土壤有機(jī)碳空間分異的一個(gè)影響因子。土地利用也被一些學(xué)者作為土壤有機(jī)碳的空間分異單元用,來表征土壤有機(jī)碳的空間分布或進(jìn)行碳庫估算[2,8]。但不足的是相對于其他相關(guān)因子如土壤類型,單一的土地利用對土壤有機(jī)碳的空間預(yù)測能力還比較低。

    目前將土地利用與克里格插值結(jié)合起來表征區(qū)域土壤有機(jī)碳空間分布的研究還不多,是否可以通過普通克里格插值來提高土地利用對土壤有機(jī)碳的空間表征精度尚不確定。因此本文選擇河北省作為研究區(qū)域,基于全國第二次土壤普查數(shù)據(jù),對比土地利用與普通克里格插值結(jié)合前后對土壤有機(jī)碳空間預(yù)測精度的差異,分析使用克里格插值來提高土地利用對土壤有機(jī)碳空間預(yù)測精度的可行性,為區(qū)域土壤有機(jī)碳空間分布的準(zhǔn)確表征提供依據(jù)。

    1 材料與方法

    1.1 研究區(qū)概況

    河北省位于北緯36°03′—42°40′,東經(jīng)113°27′—119°50′,總面積18.6萬km2。該省氣候上屬溫帶大陸性季風(fēng)氣候,大部分地區(qū)四季分明。年平均氣溫0~13℃,自北向南升高。年平均降水量300~800 mm,地區(qū)分布不均勻。水熱狀況自東南向西北呈現(xiàn)有規(guī)律的變化趨勢。全省地勢西北高,東南低。有高原、山地和平原三種地貌類型,分別占該省面積的9.3%,49.5%,41.2%。植被和土地利用上顯示出不同地貌單元之間的明顯差異性。高原主要是以多年生草本植物為主的溫帶草原,部分水分較好的地區(qū)分布有森林;山地主要是溫帶落葉闊葉林和灌草叢,寬闊平坦的盆地和谷地多辟為農(nóng)田;平原主要是耕地,是該省糧棉的主產(chǎn)區(qū)。

    1.2 資料來源

    土地利用空間數(shù)據(jù)來自中國科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心提供的20世紀(jì)80年代1∶10萬河北省土地利用數(shù)據(jù)庫(http:∥www.resdc.cn/dataResource/dataResource.asp)。一級類型為耕地、林地、草地、水域、城鄉(xiāng)工礦居民用地和未利用地6大類。經(jīng)統(tǒng)計(jì),河北省土地利用主要以耕地、林地和草地為主,分別占52.71%,19.58%和18.14%,這三類土地利用共占河北省總面積的90.43%。為便于討論,本研究只選擇耕地、林地和草地三類土地利用進(jìn)行分析。由于其它三類土地利用的圖斑比較破碎,在表征河北省土壤有機(jī)碳空間分布時(shí),先用ArcGIS軟件對這三類圖斑進(jìn)行了融合處理。土壤剖面數(shù)據(jù)來自第二次全國土壤普查[9],包括《中國土種志》、《河北土種志》和河北省地級市土種志,共計(jì)典型土壤剖面851個(gè)。土壤剖面信息包括剖面地理位置、土壤類型、母質(zhì)類型、剖面描述、土層劃分、土層厚度、有機(jī)質(zhì)含量等。土壤剖面的土地利用是根據(jù)剖面的土地利用和植被覆蓋描述按照耕地、林地和草地三種土地利用類型進(jìn)行分類。

    對于一定深度D(cm)的土壤有機(jī)碳密度DSOC的計(jì)算公式如下:

    (1)

    式中:DSOC——厚度為D土層的SOC密度(kg/m2);n——土層數(shù);δi%——第i層>2 mm礫石含量;ρi——第i層土壤容重(g/cm3);Ci——第i層SOC含量(g/kg),由土壤有機(jī)質(zhì)含量乘以0.58(Bemmelen轉(zhuǎn)換系數(shù))得到[10];Ti——第i層土層厚度(cm)。由于土地利用主要對表層土壤有機(jī)碳密度影響比較大[11-12],所以本文只估算0—20 cm深的土壤有機(jī)碳密度。

    1.3 研究方法

    (1) 土地利用法。對不同土地利用的土壤有機(jī)碳密度進(jìn)行方差分析,如果有顯著差異,說明土地利用對土壤有機(jī)碳密度有顯著影響,可以使用土地利用作為分類變量直接對土壤有機(jī)碳密度的空間分異進(jìn)行預(yù)測。土地利用對土壤有機(jī)碳空間分異的解釋能力根據(jù)不同土地利用土壤有機(jī)碳密度回歸方程的決定系數(shù)R2來確定[13]。

    (2) 土地利用與克里格插值結(jié)合的方法。先將剖面土壤有機(jī)碳密度去除土地利用類型的均值,然后再用普通克里格插值對殘差進(jìn)行進(jìn)一步預(yù)測??臻g上土壤有機(jī)碳密度最終預(yù)測值為土地利用類型均值與殘差預(yù)測值之和。預(yù)測能力通過剖面點(diǎn)實(shí)測值與預(yù)測值回歸分析的決定系數(shù)R2判定。普通克里格插值預(yù)測的方法和原理可以參考相關(guān)文獻(xiàn)[14]。

    描述性統(tǒng)計(jì)、方差分析和回歸分析采用SPSS 13.0軟件,半方差函數(shù)及理論模型擬合利用地統(tǒng)計(jì)學(xué)軟件GS+7.0,克里格插值是在ArcGIS 9.2下完成。

    2 結(jié)果與分析

    根據(jù)851個(gè)土壤剖面統(tǒng)計(jì),河北省表層土壤有機(jī)碳密度最小值為0.09 kg/m2,最大值為16.50 kg/m2,平均土壤有機(jī)碳密度2.44 kg/m2,變異系數(shù)91.2%,為高度變異。

    2.1 土地利用表征土壤有機(jī)碳空間分布

    河北省主要土地利用表層平均土壤有機(jī)碳密度見表1,三種土地利用平均土壤有機(jī)碳密度以林地最大6.68 kg/m2,其次是草地2.67 kg/m2,耕地最小為1.84 kg/m2。變異系數(shù)為56.5%~85.8%,平均變異系數(shù)68.4%,均小于河北省土壤有機(jī)碳密度變異系數(shù)。表2方差分析結(jié)果表明三種土地利用的平均土壤有機(jī)碳密度間有顯著差異(P<0.001),說明土地利用對土壤有機(jī)碳密度有顯著影響,是土壤有機(jī)碳密度空間分異的一個(gè)影響因子,因此土地利用能夠反映出河北省表層土壤有機(jī)碳密度在空間上的分布和變化。

    土地利用表征的河北省土壤有機(jī)碳密度空間分布如圖1a所示。從圖1a可以看出河北省表層土壤有機(jī)碳密度空間分布的總體趨勢是西北部主要是林地和草地,土壤有機(jī)碳密度較高,而東南部主要是耕地,土壤有機(jī)碳密度較低。回歸分析表明土地利用能夠解釋剖面土壤有機(jī)碳密度空間分異的28.2%(P<0.001)。

    表1 河北省主要土地利用的土壤有機(jī)碳密度統(tǒng)計(jì)

    表2 不同土地利用土壤有機(jī)碳密度的方差分析

    圖1 不同方法表征的河北省土壤有機(jī)碳密度空間分布

    2.2土地利用與克里格插值結(jié)合表征土壤有機(jī)碳空間分布

    對剖面土壤有機(jī)碳密度去除土地利用均值后的殘差進(jìn)行K-S檢驗(yàn),結(jié)果為正態(tài)分布。殘差半方差擬合模型及參數(shù)見表3,擬合模型為指數(shù)模型,決定系數(shù)R2為0.705,模型擬合結(jié)果較好。模型的塊金方差占基臺方差的50%,屬于中等程度的空間自相關(guān),因此可以使用普通克里值插值對殘差進(jìn)行空間預(yù)測,結(jié)果如圖1b所示。從預(yù)測結(jié)果來,河北省西北部表層土壤有機(jī)碳密度的殘差空間分布值較高,而東南部殘差的分布值相對較低。說明根據(jù)土地利用來預(yù)測河北省土壤有機(jī)碳密度的空間分布時(shí),在西北部,對土壤有機(jī)碳密度的估計(jì)值總體上要低于實(shí)際值,而東南部總體上會(huì)高估土壤有機(jī)碳密度值。

    表3 土壤有機(jī)碳密度殘差半方差函數(shù)擬合模型參數(shù)

    將土地利用表征的土壤有機(jī)碳密度空間分布與普通克里格插值預(yù)測的殘差空間分布進(jìn)行疊加,如圖1c所示。與土地利用表征的河北省土壤有機(jī)碳密度空間分布圖1a相比,土地利用與普通克里格插值結(jié)合不僅能夠反映出土壤有機(jī)碳密度不同土地利用之間的差異,而且能夠反映出同一土地利用內(nèi)土壤有機(jī)碳密度的空間分異,因此要比單獨(dú)使用土地利用表征的土壤有機(jī)碳密度的空間分布效果要好?;貧w分析表明,土地利用與普通克里格插值結(jié)合能夠解釋剖面土壤有機(jī)碳密度空間分異的43.0%(P<0.001)。

    3 討 論

    在土壤有機(jī)碳的各相關(guān)因子中,以土壤類型尤其是低級土壤分類對土壤有機(jī)碳空間分異的反映能力較好,碳庫估算也多以土壤類型作為土壤有機(jī)碳的基本空間分異單元進(jìn)行估算[15-16]。由于土壤有機(jī)碳作為碳源或碳匯可能影響全球氣候,所以有必要去量化土壤有機(jī)碳的變化[17]。而土地利用是影響土壤有機(jī)碳含量高低的一個(gè)重要因子,當(dāng)土地利用發(fā)生改變以后,土壤有機(jī)碳含量會(huì)隨之發(fā)生變化[18]。目前,隨著遙感技術(shù)的發(fā)展,土地利用的空間變化能夠通過遙感數(shù)字影象快速獲取,而且土地利用的類別往往比較少。因此相對于其它相關(guān)因子例如土壤類型等,以土地利用表征土壤有機(jī)碳的空間分布或進(jìn)行制圖不僅方便快捷,還可以通過不同時(shí)期土地利用的對比來預(yù)測土壤有機(jī)碳的動(dòng)態(tài)變化。

    目前很多的研究也表明了土地利用對區(qū)域土壤有機(jī)碳空間分異有很好的解釋能力,是區(qū)域土壤有機(jī)碳空間分異的一個(gè)重要影響因子。例如Bell等[2]對英國紐卡斯?fàn)柲喜?55 km2)的研究表明,土地利用對土壤有機(jī)碳空間預(yù)測能力能達(dá)到44%。但在較大的空間區(qū)域,土地利用對土壤有機(jī)碳空間分異的解釋能力相對較差,例如張忠啟[19]對江西余江縣的研究表明土地利用只能獨(dú)立解釋土壤有機(jī)碳含量空間變異的32.3%,本研究中土地利用對河北省土壤有機(jī)碳的空間變異解釋能力僅為28.2%,同一土地利用類型內(nèi)還存在較大的空間變異,平均變異系數(shù)達(dá)到68.4%(表1)。而且與常用的土壤有機(jī)碳相關(guān)因子如土壤類型相比較,土地利用對土壤有機(jī)碳空間分異的解釋能力往往不如低級土壤分類高。例如Bell等[2]的研究表明土系對土壤有機(jī)碳的空間預(yù)測能力為48%,張忠啟[19]在江西余江縣研究表明土屬能夠解釋土壤有機(jī)碳空間分異的41.4%,都要高于相同區(qū)域土地利用對土壤有機(jī)碳空間分異的解釋能力。因此對區(qū)域土壤有機(jī)碳的空間預(yù)測,尤其是較大的空間區(qū)域,單獨(dú)使用土地利用效果還不是很好,有必要在土地利用的基礎(chǔ)上,進(jìn)一步提高土壤有機(jī)碳的空間預(yù)測精度。

    將土地利用與其它相關(guān)因子結(jié)合無疑是提高預(yù)測精度的一種方法,但這需要獲取其它相關(guān)因子的資料。而本研究將土地利用與克里格插值結(jié)合,不僅不需要獲取其它的空間數(shù)據(jù)和屬性數(shù)據(jù),還能將土壤有機(jī)碳密度空間分異的解釋程度顯著提高到43.0%,在精度上與張忠啟和Bell等研究的低級土壤分類土屬和土系對土壤有機(jī)碳空間變異的解釋能力也已比較接近[2,19]。這說明克里格插值能夠用來優(yōu)化土地利用對土壤有機(jī)碳的空間預(yù)測,是提高區(qū)域土壤有機(jī)碳空間表征精度的一個(gè)有效方法。

    4 結(jié) 論

    河北省主要的三個(gè)土地利用耕地、林地和草地表層土壤有機(jī)碳密度之間具有顯著差異(P<0.001),說明土地利用是表層土壤有機(jī)碳密度空間分異的一個(gè)影響因子。但土地利用僅能預(yù)測河北省表層土壤有機(jī)碳密度空間分異的28.2%,而土地利用與普通克里格插值結(jié)合以后能將對土壤有機(jī)碳密度空間分異的預(yù)測能力顯著提高到43.0%。因此區(qū)域空間上可以在土地利用的基礎(chǔ)上再結(jié)合普通克里格插值來優(yōu)化表層土壤有機(jī)碳密度的空間表征。

    [1] 申廣榮,葛曉燁,黃秀梅.上海崇明島表層土壤有機(jī)碳密度的空間分布特征及碳儲(chǔ)量估算[J].上海交通大學(xué)學(xué)報(bào):農(nóng)業(yè)科學(xué)版,2011,29(6):61-66.

    [2] Bell M J, Worrall F. Estimating a region′s soil organic carbon baseline: The undervalued role of land-management[J]. Geoderma,2009,152(1/2):74-84.

    [3] 趙永存,史學(xué)正,于東升,等.不同方法預(yù)測河北省土壤有機(jī)碳密度空間分布特征的研究[J].土壤學(xué)報(bào),2005,42(3):379-385.

    [4] Chaplot V, Bouahom B, Valentin C. Soil organic carbon stocks in Laos: spatial variations and controlling factors[J]. Global Change Biology,2010,16(4):1380-1393.

    [5] White II D A, Welty-Bernard A, Rasmussen C, et al. Vegetation controls on soil organic carbon dynamics in an arid, hyperthermic ecosystem[J]. Geoderma,2009,150(1):214-223.

    [6] Zhang M, Zhang X K, Liang W J, et al. Distribution of soil organic carbon fractions along the altitudinal gradient in Changbai Mountain, China[M]. Pedosphere,2011,21(5):615-620.

    [7] Chiti T, Gardin L, Perugini L, et al. Soil organic carbon stock assessment for the different cropland land uses in Italy[J]. Biology and Fertility of Soils,2012,48(1):9-17.

    [8] 徐小明.通榆縣土壤有機(jī)碳儲(chǔ)量估算及時(shí)空分布研究[D].長春:吉林大學(xué),2008.

    [9] 全國土壤普查辦公室.中國土種志(第四卷)[M].北京:中國農(nóng)業(yè)出版社,1995.

    [10] 文啟孝.土壤有機(jī)質(zhì)研究方法[M].北京:農(nóng)業(yè)出版社,1984.

    [11] Cheng H H, Kimble J M. Characterization of soil organic carbon pools[M]∥Lal R, Kimble J M, Follett R F, et al. Assessment Methods for Soil Carbon. London Lewis,2001:117-129.

    [12] Kimble J M, Grossman R B, Samson-Liebiog S E. Methodology for sampling and preparation for soil carbon determinations[M]∥Lal R, Kimble J M, Follett R F, et al. Assessment Methods for Soil Carbon. London: Lewis,2001:15-29.

    [13] 李麗霞,郜艷暉,張瑛.啞變量在統(tǒng)計(jì)分析中的應(yīng)用[J].數(shù)理醫(yī)藥學(xué)雜志,2006,19(1):51-53.

    [14] 范勝龍,黃炎和,林金石.表征土壤有機(jī)碳區(qū)域分布的優(yōu)化空間插值模型研究:以福建省龍海市為例[J].水土保持研究,2012,18(6):1-5.

    [15] 姜小三,潘劍君,李學(xué)林.江蘇表層土壤有機(jī)碳密度和儲(chǔ)量估算和空間分布分析[J].土壤通報(bào),2005,36(4):501-503.

    [16] 袁芳,趙小敏,樂麗紅,等.江西省表層土壤有機(jī)碳庫儲(chǔ)量估算與空間分布特征[J].生態(tài)環(huán)境,2008,17(1):268-272.

    [17] Stolt M H, Drohan P J, Richardson M J. Insights and approaches for mapping soil organic carbon as a dynamic soil property[J]. Soil Science Society of America Journal,2010,74(5):1685-1689.

    [18] Jeny H, Raychaudhuri S P. Effect of Climate and Cultivation on Nitrogen and Organic Matter Reserves in Indian Soils[M]. Indian Council of Agricultural Research, New Delhi,1960.

    [19] 張忠啟.樣點(diǎn)布置模式及密度對揭示土壤有機(jī)碳空間變異的影響[D].南京:中國科學(xué)院南京土壤研究所,2010.

    IndicatingMethodonSpatialDistributionofRegionalSoilOrganicCarbonbyCombiningLandUseandKriging

    GU Cheng-jun1,2,3, SHI Xue-zheng1, YU Dong-sheng1

    (1.StateKeyLaboratoryofSoilandSustainableAgriculture,InstituteofSoilScience,ChineseAcademyofSciences,Nanjing210008,China; 2.UniversityoftheChineseAcademyofSciences,Beijing100049,China; 3.GeographicInformationandTourismCollege,ChuzhouUniversity,Chuzhou239000,China)

    How to use limited soil samplings to predicate an accurate spatial distribution of soil properties is always a hot issue. Method of using land use and Kriging to indicate soil organic carbon spatial distribution was explored by variance analysis with the data of Hebei Province of the Second National Soil Survey. The results showed that land use was a factor impacting soil organic carbon spatial variation and could be used to reflect the spatial distribution of soil organic carbon. But the independent reflecting ability of land use to variation of SOC density in the surface was not very high, land use and ordinary Kriging interpolation in combination could significantly increase explanatory strength of SOC spatial variation. Therefore, Kriging interpolation could be combined with land use to indicate the spatial distribution of regional soil organic carbon in the topsoil.

    land use; soil organic carbon; Kriging; spatial distribution

    2013-08-07

    :2013-08-27

    中國科學(xué)院戰(zhàn)略性先導(dǎo)科技專項(xiàng)項(xiàng)日(XDAO5050509);國家重點(diǎn)基礎(chǔ)研究發(fā)展規(guī)劃(973計(jì)劃)項(xiàng)目(2010CB950702)

    顧成軍(1978—),男,安徽滁州人,博士研究生,研究方向:土壤資源與遙感信息。E-mail:lacjg@163.com

    史學(xué)正(1959—),男,浙江嵊縣人,博士,研究員,研究方向:土壤資源數(shù)字化管理。E-mail:xzshi@issas.ac.cn

    S153.6

    :A

    :1005-3409(2014)02-0039-04

    猜你喜歡
    分異克里土壤有機(jī)
    今晚不能去你家玩啦!
    知識窗(2023年12期)2024-01-03 01:38:55
    我可以咬一口嗎?
    知識窗(2023年2期)2023-03-05 11:28:27
    你今天真好看
    你今天真好看
    讀者(2018年24期)2018-12-04 03:01:34
    重慶市臭氧時(shí)空分異及其影響因素研究
    西雙版納橡膠林土壤有機(jī)碳分布特征研究
    秸稈還田的土壤有機(jī)碳周轉(zhuǎn)特征
    平泉縣下營坊雜巖體分異演化及其成巖成礦
    土壤有機(jī)碳轉(zhuǎn)化研究及其進(jìn)展
    AMDIS在土壤有機(jī)污染物鑒別中表征性統(tǒng)計(jì)量的探究
    成人高潮视频无遮挡免费网站| 在线观看免费午夜福利视频| 国产av不卡久久| 全区人妻精品视频| 淫妇啪啪啪对白视频| 亚洲精品久久国产高清桃花| 岛国视频午夜一区免费看| 在线十欧美十亚洲十日本专区| 国产av一区在线观看免费| 黄色片一级片一级黄色片| 午夜福利欧美成人| 国产精品亚洲一级av第二区| 成人一区二区视频在线观看| 99久久综合精品五月天人人| 中亚洲国语对白在线视频| 日本免费a在线| 一二三四社区在线视频社区8| 亚洲国产看品久久| 午夜福利在线在线| 麻豆久久精品国产亚洲av| 又紧又爽又黄一区二区| 欧美乱码精品一区二区三区| 日韩欧美在线乱码| 亚洲av成人精品一区久久| 日韩有码中文字幕| 性色av乱码一区二区三区2| 久久精品夜夜夜夜夜久久蜜豆| 国产三级在线视频| 亚洲精品一区av在线观看| 曰老女人黄片| 国产亚洲精品一区二区www| 亚洲午夜理论影院| 亚洲av成人精品一区久久| 欧美极品一区二区三区四区| 免费在线观看亚洲国产| 亚洲人成电影免费在线| 欧美日韩福利视频一区二区| 后天国语完整版免费观看| 欧美成人性av电影在线观看| avwww免费| 日韩精品青青久久久久久| 国产亚洲av高清不卡| 久久精品国产综合久久久| 国产探花在线观看一区二区| 日本黄色片子视频| 老司机午夜十八禁免费视频| 亚洲成av人片免费观看| 老鸭窝网址在线观看| 97超视频在线观看视频| 99热只有精品国产| 天堂√8在线中文| 久久精品国产亚洲av香蕉五月| 热99在线观看视频| 色哟哟哟哟哟哟| 欧美日韩精品网址| 怎么达到女性高潮| 久久久久国产一级毛片高清牌| 国产成人影院久久av| 极品教师在线免费播放| 九色成人免费人妻av| 久久这里只有精品中国| 日本精品一区二区三区蜜桃| 欧美黄色片欧美黄色片| 免费看日本二区| 欧美色视频一区免费| 色av中文字幕| 99热精品在线国产| 成人高潮视频无遮挡免费网站| 99视频精品全部免费 在线 | 99久久精品热视频| 美女免费视频网站| 在线观看66精品国产| 国产精品99久久久久久久久| 天天一区二区日本电影三级| 国产激情偷乱视频一区二区| 亚洲国产中文字幕在线视频| 欧美日韩亚洲国产一区二区在线观看| 国产又黄又爽又无遮挡在线| 九色国产91popny在线| 国产午夜精品论理片| 婷婷亚洲欧美| 麻豆成人av在线观看| 一进一出抽搐gif免费好疼| 人人妻,人人澡人人爽秒播| www日本黄色视频网| 国产成年人精品一区二区| 国产精品爽爽va在线观看网站| a在线观看视频网站| 99久久精品国产亚洲精品| 一本精品99久久精品77| 小蜜桃在线观看免费完整版高清| 757午夜福利合集在线观看| 日本黄色视频三级网站网址| 99国产极品粉嫩在线观看| 午夜a级毛片| 免费高清视频大片| 男女之事视频高清在线观看| 桃红色精品国产亚洲av| 国内精品久久久久久久电影| 亚洲国产中文字幕在线视频| 麻豆一二三区av精品| 欧美激情久久久久久爽电影| 中文亚洲av片在线观看爽| 伊人久久大香线蕉亚洲五| 日本五十路高清| 国产av在哪里看| 国产av在哪里看| 国产精品久久久久久亚洲av鲁大| 伊人久久大香线蕉亚洲五| 18禁观看日本| 男人舔奶头视频| 久久精品人妻少妇| 欧美三级亚洲精品| 日韩中文字幕欧美一区二区| 国产黄片美女视频| 超碰成人久久| 两性午夜刺激爽爽歪歪视频在线观看| netflix在线观看网站| 变态另类成人亚洲欧美熟女| 12—13女人毛片做爰片一| 国产一区二区三区视频了| 手机成人av网站| 香蕉久久夜色| 欧美绝顶高潮抽搐喷水| 香蕉久久夜色| 欧美乱色亚洲激情| 免费看光身美女| 久久精品亚洲精品国产色婷小说| 99国产精品99久久久久| 国产探花在线观看一区二区| 国产精品野战在线观看| 日韩大尺度精品在线看网址| 九色国产91popny在线| 国产精品一区二区免费欧美| 91在线观看av| 国产精品98久久久久久宅男小说| 一区二区三区高清视频在线| 国产一区在线观看成人免费| 热99在线观看视频| 亚洲精品粉嫩美女一区| 亚洲性夜色夜夜综合| 久久天躁狠狠躁夜夜2o2o| 亚洲国产欧美一区二区综合| 后天国语完整版免费观看| 三级毛片av免费| 人妻久久中文字幕网| av黄色大香蕉| 啦啦啦观看免费观看视频高清| 成人永久免费在线观看视频| 少妇熟女aⅴ在线视频| 超碰成人久久| 国产黄片美女视频| 一二三四在线观看免费中文在| 精品一区二区三区视频在线观看免费| 国产乱人视频| 成人三级黄色视频| 国产伦精品一区二区三区视频9 | 18禁黄网站禁片午夜丰满| 男插女下体视频免费在线播放| 久9热在线精品视频| 狂野欧美白嫩少妇大欣赏| 午夜激情欧美在线| 成人三级做爰电影| 亚洲在线自拍视频| 一个人看视频在线观看www免费 | 久久久久国产一级毛片高清牌| 最新在线观看一区二区三区| 精品久久久久久久久久久久久| 欧洲精品卡2卡3卡4卡5卡区| 亚洲欧美激情综合另类| 法律面前人人平等表现在哪些方面| 国产精品一区二区免费欧美| 久久久久免费精品人妻一区二区| 国产午夜精品论理片| 999久久久国产精品视频| 色吧在线观看| 午夜亚洲福利在线播放| 国产美女午夜福利| 亚洲欧洲精品一区二区精品久久久| 欧美乱色亚洲激情| 嫩草影视91久久| 日韩欧美三级三区| 亚洲国产精品成人综合色| 国产日本99.免费观看| 精品久久久久久久久久免费视频| 特级一级黄色大片| 国产爱豆传媒在线观看| 久久精品国产清高在天天线| 性色av乱码一区二区三区2| 1024手机看黄色片| 久久久久精品国产欧美久久久| 免费看光身美女| 国产av在哪里看| 国产三级黄色录像| 日本一二三区视频观看| 午夜精品久久久久久毛片777| 很黄的视频免费| 欧美国产日韩亚洲一区| 亚洲av成人精品一区久久| 亚洲av片天天在线观看| 99国产精品一区二区三区| 婷婷亚洲欧美| 一二三四在线观看免费中文在| 亚洲中文字幕日韩| 男人舔女人的私密视频| 丁香欧美五月| 91麻豆精品激情在线观看国产| 久久久久久久久久黄片| 亚洲熟妇熟女久久| 国产精品av久久久久免费| 欧美zozozo另类| 嫁个100分男人电影在线观看| 久久久国产成人免费| 中文字幕精品亚洲无线码一区| 女人被狂操c到高潮| 真人做人爱边吃奶动态| av在线蜜桃| 国产一区二区激情短视频| 一个人免费在线观看电影 | 精品无人区乱码1区二区| 欧美成人免费av一区二区三区| 最近最新中文字幕大全电影3| 国产激情欧美一区二区| 欧美乱色亚洲激情| 亚洲成人免费电影在线观看| 日韩欧美 国产精品| av天堂在线播放| 草草在线视频免费看| 亚洲熟妇熟女久久| 搡老妇女老女人老熟妇| 一本久久中文字幕| 国产美女午夜福利| 欧美日韩亚洲国产一区二区在线观看| 黑人欧美特级aaaaaa片| 亚洲av片天天在线观看| av福利片在线观看| 大型黄色视频在线免费观看| 日韩 欧美 亚洲 中文字幕| 五月玫瑰六月丁香| 久久久精品欧美日韩精品| 久久九九热精品免费| 99热这里只有精品一区 | 人妻夜夜爽99麻豆av| 色在线成人网| 伊人久久大香线蕉亚洲五| 美女大奶头视频| 日韩人妻高清精品专区| 人人妻人人澡欧美一区二区| 一进一出好大好爽视频| 欧美性猛交黑人性爽| 嫁个100分男人电影在线观看| 欧美3d第一页| 中文字幕久久专区| 欧美激情久久久久久爽电影| 两个人的视频大全免费| 99re在线观看精品视频| 亚洲精品美女久久久久99蜜臀| 日韩欧美免费精品| 69av精品久久久久久| 18禁黄网站禁片午夜丰满| 亚洲天堂国产精品一区在线| 黑人操中国人逼视频| 中文亚洲av片在线观看爽| 国产三级在线视频| 99精品欧美一区二区三区四区| 男女之事视频高清在线观看| bbb黄色大片| 精品久久久久久久久久久久久| 国产伦在线观看视频一区| 国产精品久久久久久久电影 | 欧美三级亚洲精品| 少妇人妻一区二区三区视频| a级毛片a级免费在线| 十八禁人妻一区二区| 变态另类成人亚洲欧美熟女| 99国产精品一区二区三区| 午夜视频精品福利| 久久天堂一区二区三区四区| 日日干狠狠操夜夜爽| 久久久久性生活片| 免费一级毛片在线播放高清视频| 黄色成人免费大全| 美女大奶头视频| 亚洲精品中文字幕一二三四区| 最近在线观看免费完整版| 欧美一级毛片孕妇| 成人精品一区二区免费| 国产精品久久久久久亚洲av鲁大| 天堂√8在线中文| 欧美色视频一区免费| 久久精品国产综合久久久| 少妇人妻一区二区三区视频| 国产精品女同一区二区软件 | 亚洲人与动物交配视频| 亚洲中文日韩欧美视频| 999久久久精品免费观看国产| 九九久久精品国产亚洲av麻豆 | 91在线精品国自产拍蜜月 | 国产精品 国内视频| 禁无遮挡网站| 法律面前人人平等表现在哪些方面| 免费在线观看视频国产中文字幕亚洲| av中文乱码字幕在线| 舔av片在线| 久久久久国产一级毛片高清牌| 色老头精品视频在线观看| 国内揄拍国产精品人妻在线| 99riav亚洲国产免费| 一级毛片女人18水好多| www.自偷自拍.com| 欧美日韩乱码在线| 蜜桃久久精品国产亚洲av| 九九在线视频观看精品| 深夜精品福利| 一级作爱视频免费观看| 99久久99久久久精品蜜桃| 成人三级黄色视频| 国产亚洲av嫩草精品影院| 99国产精品99久久久久| 久久人人精品亚洲av| 免费观看人在逋| 亚洲avbb在线观看| 欧美三级亚洲精品| 国产亚洲欧美在线一区二区| 亚洲人成伊人成综合网2020| 欧美绝顶高潮抽搐喷水| 婷婷六月久久综合丁香| 好男人在线观看高清免费视频| 波多野结衣巨乳人妻| 日韩欧美精品v在线| 亚洲av第一区精品v没综合| 久久伊人香网站| 国内精品美女久久久久久| 九九热线精品视视频播放| av女优亚洲男人天堂 | www.自偷自拍.com| 国产av不卡久久| 久久性视频一级片| 色尼玛亚洲综合影院| 一个人看视频在线观看www免费 | 精品不卡国产一区二区三区| 欧美在线黄色| 亚洲性夜色夜夜综合| 国产精品美女特级片免费视频播放器 | 久久久久性生活片| 成人一区二区视频在线观看| 日本黄大片高清| 亚洲激情在线av| 欧美在线一区亚洲| 91麻豆av在线| 一本综合久久免费| av中文乱码字幕在线| 成年女人永久免费观看视频| 一本一本综合久久| 久久国产精品人妻蜜桃| 一二三四在线观看免费中文在| 欧美在线黄色| 国产黄色小视频在线观看| 午夜激情欧美在线| 国产欧美日韩精品一区二区| 亚洲电影在线观看av| 久久久久久久久久黄片| 男人的好看免费观看在线视频| av黄色大香蕉| 全区人妻精品视频| 88av欧美| 欧美日韩福利视频一区二区| 国产1区2区3区精品| 久久欧美精品欧美久久欧美| 久久精品影院6| 无人区码免费观看不卡| 中文资源天堂在线| 日本与韩国留学比较| 搡老妇女老女人老熟妇| 人妻久久中文字幕网| 激情在线观看视频在线高清| 精品国产三级普通话版| 亚洲精品久久国产高清桃花| 真实男女啪啪啪动态图| 亚洲精品在线美女| 国产成人福利小说| 中文字幕精品亚洲无线码一区| 欧美乱色亚洲激情| 99re在线观看精品视频| 成人无遮挡网站| 国产av一区在线观看免费| 亚洲国产高清在线一区二区三| 午夜免费观看网址| tocl精华| 99国产精品99久久久久| 午夜免费成人在线视频| 久久精品综合一区二区三区| 亚洲人与动物交配视频| 黑人操中国人逼视频| 国产精品久久久av美女十八| 麻豆av在线久日| 麻豆一二三区av精品| 久久久国产欧美日韩av| 亚洲av片天天在线观看| 美女免费视频网站| 一个人观看的视频www高清免费观看 | 免费看十八禁软件| 亚洲人成网站高清观看| 精品电影一区二区在线| 噜噜噜噜噜久久久久久91| 首页视频小说图片口味搜索| 国产精品一区二区三区四区免费观看 | 国内毛片毛片毛片毛片毛片| 欧美成人免费av一区二区三区| 一级黄色大片毛片| 深夜精品福利| 99久久国产精品久久久| 久久久久免费精品人妻一区二区| 亚洲自拍偷在线| 欧美zozozo另类| 一本精品99久久精品77| 两个人的视频大全免费| or卡值多少钱| 99久国产av精品| 久久国产精品影院| 亚洲精品美女久久久久99蜜臀| 非洲黑人性xxxx精品又粗又长| 亚洲专区字幕在线| av欧美777| 国产精品久久视频播放| 日韩欧美国产在线观看| 欧美日韩黄片免| 床上黄色一级片| 一进一出抽搐动态| 国产亚洲精品av在线| 国产三级在线视频| 18禁美女被吸乳视频| 法律面前人人平等表现在哪些方面| 国产亚洲精品久久久久久毛片| 露出奶头的视频| 亚洲一区高清亚洲精品| 人妻夜夜爽99麻豆av| 亚洲性夜色夜夜综合| 精品99又大又爽又粗少妇毛片 | 操出白浆在线播放| 日日夜夜操网爽| 亚洲av成人不卡在线观看播放网| 一夜夜www| 老司机福利观看| 又紧又爽又黄一区二区| 欧美黑人欧美精品刺激| 在线免费观看不下载黄p国产 | 久久久水蜜桃国产精品网| 日韩中文字幕欧美一区二区| 久久久色成人| 成人精品一区二区免费| 国产精品香港三级国产av潘金莲| 国产亚洲av嫩草精品影院| 欧美绝顶高潮抽搐喷水| 成年女人永久免费观看视频| 九色国产91popny在线| 日本与韩国留学比较| 我要搜黄色片| 久久国产精品人妻蜜桃| 人人妻,人人澡人人爽秒播| 欧美黑人巨大hd| 波多野结衣巨乳人妻| 麻豆国产97在线/欧美| 亚洲欧洲精品一区二区精品久久久| 色在线成人网| 亚洲人与动物交配视频| 亚洲精品在线美女| 亚洲国产欧美一区二区综合| 国产精品亚洲美女久久久| 国产一区二区在线av高清观看| 天堂影院成人在线观看| 中文字幕人成人乱码亚洲影| 亚洲成a人片在线一区二区| 1000部很黄的大片| 又大又爽又粗| 成年版毛片免费区| 日本免费a在线| 看免费av毛片| 国产91精品成人一区二区三区| 人妻久久中文字幕网| 午夜影院日韩av| 99久久国产精品久久久| 身体一侧抽搐| 久久久久国内视频| 色老头精品视频在线观看| 国产精品av久久久久免费| 国产精品电影一区二区三区| 99热只有精品国产| 蜜桃久久精品国产亚洲av| 午夜福利18| 精品无人区乱码1区二区| 黄色 视频免费看| 一本久久中文字幕| av福利片在线观看| 亚洲aⅴ乱码一区二区在线播放| 欧美日本视频| 亚洲av电影在线进入| 麻豆一二三区av精品| 国产av不卡久久| 亚洲国产精品久久男人天堂| 亚洲美女黄片视频| 可以在线观看的亚洲视频| 一级黄色大片毛片| 免费大片18禁| 国产激情欧美一区二区| 亚洲美女视频黄频| 久久99热这里只有精品18| 国产精品 欧美亚洲| 国产高清视频在线播放一区| 久久久水蜜桃国产精品网| 欧美极品一区二区三区四区| 亚洲专区中文字幕在线| 成人欧美大片| 午夜福利欧美成人| 99精品久久久久人妻精品| 国产精品1区2区在线观看.| 国产精品av久久久久免费| 夜夜躁狠狠躁天天躁| 最新中文字幕久久久久 | 最近在线观看免费完整版| 在线十欧美十亚洲十日本专区| 国产精品九九99| 午夜久久久久精精品| 少妇人妻一区二区三区视频| 亚洲人成电影免费在线| 国产精品98久久久久久宅男小说| 成在线人永久免费视频| 国内精品久久久久精免费| 婷婷六月久久综合丁香| 久久久国产成人精品二区| 18禁观看日本| 毛片女人毛片| 变态另类成人亚洲欧美熟女| 18禁美女被吸乳视频| 精品久久久久久久人妻蜜臀av| 国产91精品成人一区二区三区| 99久国产av精品| 欧美又色又爽又黄视频| 1024手机看黄色片| 久久香蕉精品热| 日本 欧美在线| 少妇熟女aⅴ在线视频| 亚洲性夜色夜夜综合| 国产精品电影一区二区三区| 国产私拍福利视频在线观看| 99久久综合精品五月天人人| 亚洲精品在线观看二区| 午夜精品在线福利| 身体一侧抽搐| bbb黄色大片| 精品人妻1区二区| 久久精品国产99精品国产亚洲性色| 伊人久久大香线蕉亚洲五| 午夜福利视频1000在线观看| av天堂中文字幕网| cao死你这个sao货| 免费在线观看视频国产中文字幕亚洲| 亚洲狠狠婷婷综合久久图片| 天堂av国产一区二区熟女人妻| 一二三四社区在线视频社区8| 久久久精品大字幕| 亚洲av免费在线观看| 天堂网av新在线| 国产一区二区在线av高清观看| 丁香六月欧美| 久久久久久久午夜电影| 真人做人爱边吃奶动态| 视频区欧美日本亚洲| 国产1区2区3区精品| 久久久久久大精品| 窝窝影院91人妻| 亚洲精品国产精品久久久不卡| 高潮久久久久久久久久久不卡| 欧美黄色片欧美黄色片| 国产 一区 欧美 日韩| 天天躁日日操中文字幕| 国产免费男女视频| 亚洲七黄色美女视频| www日本在线高清视频| 精品免费久久久久久久清纯| 日本熟妇午夜| 亚洲成a人片在线一区二区| 久久久久久九九精品二区国产| 好男人在线观看高清免费视频| 日本黄色视频三级网站网址| 黄片大片在线免费观看| 成年女人看的毛片在线观看| 黄片大片在线免费观看| 老汉色∧v一级毛片| 可以在线观看的亚洲视频| 午夜福利在线在线| 日本黄色片子视频| 美女大奶头视频| 熟女少妇亚洲综合色aaa.| 日本 av在线| 亚洲第一电影网av| 久久精品aⅴ一区二区三区四区| 99热这里只有精品一区 | 成年女人永久免费观看视频| 一级毛片精品| 女生性感内裤真人,穿戴方法视频| 窝窝影院91人妻| 色综合站精品国产| 一本综合久久免费| 欧美在线一区亚洲| 亚洲色图av天堂| 亚洲精品久久国产高清桃花| 国产精品亚洲一级av第二区| 18禁裸乳无遮挡免费网站照片| 国产亚洲av高清不卡| 十八禁人妻一区二区| 国产99白浆流出| 在线观看免费午夜福利视频| 亚洲精品在线观看二区| 狂野欧美激情性xxxx|