• <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ì)量的探究
    在线观看免费视频日本深夜| 日韩人妻精品一区2区三区| 国产av一区二区精品久久| 男女午夜视频在线观看| ponron亚洲| 欧美不卡视频在线免费观看 | 成年女人毛片免费观看观看9| 国产一区二区三区综合在线观看| 国产精品久久久av美女十八| 国产亚洲av高清不卡| 欧美性长视频在线观看| 久久久久久久久免费视频了| 在线播放国产精品三级| 久9热在线精品视频| 又大又爽又粗| 亚洲一码二码三码区别大吗| 51午夜福利影视在线观看| 国产又色又爽无遮挡免费看| 午夜久久久在线观看| 深夜精品福利| 国产精品亚洲一级av第二区| 男人操女人黄网站| 一区二区日韩欧美中文字幕| 在线十欧美十亚洲十日本专区| 午夜精品在线福利| 最好的美女福利视频网| 夜夜躁狠狠躁天天躁| 亚洲一区二区三区欧美精品| 18禁观看日本| 精品日产1卡2卡| 天堂√8在线中文| 国产精品亚洲一级av第二区| 亚洲欧美一区二区三区久久| 女警被强在线播放| 日本撒尿小便嘘嘘汇集6| 免费在线观看黄色视频的| 国产成人啪精品午夜网站| 91国产中文字幕| 精品福利永久在线观看| 女警被强在线播放| 亚洲欧美精品综合久久99| 欧美最黄视频在线播放免费 | 久久香蕉精品热| 国产亚洲精品久久久久久毛片| 一二三四在线观看免费中文在| 12—13女人毛片做爰片一| 成人av一区二区三区在线看| 久久欧美精品欧美久久欧美| 欧美亚洲日本最大视频资源| 午夜精品国产一区二区电影| 日韩精品青青久久久久久| 波多野结衣高清无吗| 淫妇啪啪啪对白视频| 色综合欧美亚洲国产小说| 国产精品偷伦视频观看了| 午夜福利在线观看吧| 日本免费a在线| 久久精品91无色码中文字幕| 黄色视频,在线免费观看| 精品国产一区二区三区四区第35| 亚洲av成人不卡在线观看播放网| 久久国产精品影院| 搡老熟女国产l中国老女人| 9191精品国产免费久久| 99国产精品一区二区三区| 国产精品国产高清国产av| 激情视频va一区二区三区| 欧美av亚洲av综合av国产av| 大型av网站在线播放| 自拍欧美九色日韩亚洲蝌蚪91| 日韩免费高清中文字幕av| 亚洲美女黄片视频| 久久99一区二区三区| 老司机午夜福利在线观看视频| 欧美乱色亚洲激情| 日韩欧美一区视频在线观看| 国产精华一区二区三区| 亚洲精品久久午夜乱码| 国产av在哪里看| 91成年电影在线观看| 一区在线观看完整版| 可以免费在线观看a视频的电影网站| 国产亚洲欧美精品永久| 高清av免费在线| 日日爽夜夜爽网站| 757午夜福利合集在线观看| 夜夜看夜夜爽夜夜摸 | 涩涩av久久男人的天堂| 色综合站精品国产| 99在线视频只有这里精品首页| 极品人妻少妇av视频| 国产精品 欧美亚洲| 两个人看的免费小视频| 国产精品久久电影中文字幕| 日本一区二区免费在线视频| 欧美黄色淫秽网站| 制服人妻中文乱码| 一级毛片女人18水好多| 91精品国产国语对白视频| 欧美乱色亚洲激情| 夜夜看夜夜爽夜夜摸 | 一个人免费在线观看的高清视频| 久久人妻av系列| 亚洲午夜理论影院| 999精品在线视频| 人人妻,人人澡人人爽秒播| 国产xxxxx性猛交| 婷婷六月久久综合丁香| 日本黄色日本黄色录像| 18禁美女被吸乳视频| 国产色视频综合| 97超级碰碰碰精品色视频在线观看| 午夜久久久在线观看| 亚洲色图综合在线观看| 热re99久久精品国产66热6| 国产午夜精品久久久久久| 免费在线观看亚洲国产| 国产精品一区二区精品视频观看| 国产一卡二卡三卡精品| 黄色视频不卡| 亚洲av片天天在线观看| a在线观看视频网站| 每晚都被弄得嗷嗷叫到高潮| 少妇被粗大的猛进出69影院| 桃红色精品国产亚洲av| 成人18禁高潮啪啪吃奶动态图| 成人18禁在线播放| 老熟妇乱子伦视频在线观看| 人人妻人人爽人人添夜夜欢视频| 一级,二级,三级黄色视频| 日本三级黄在线观看| 热99国产精品久久久久久7| 亚洲精品一区av在线观看| 国产黄a三级三级三级人| 精品熟女少妇八av免费久了| 在线观看舔阴道视频| 婷婷精品国产亚洲av在线| 五月开心婷婷网| 国产成人一区二区三区免费视频网站| 欧美性长视频在线观看| 亚洲成人精品中文字幕电影 | 亚洲在线自拍视频| 美女高潮喷水抽搐中文字幕| 曰老女人黄片| 精品午夜福利视频在线观看一区| 狂野欧美激情性xxxx| 国产成+人综合+亚洲专区| 波多野结衣av一区二区av| 亚洲精品中文字幕一二三四区| 精品午夜福利视频在线观看一区| 久久久久久免费高清国产稀缺| 99精品欧美一区二区三区四区| 国产欧美日韩一区二区精品| 成人永久免费在线观看视频| 丰满的人妻完整版| 99热国产这里只有精品6| 久久国产亚洲av麻豆专区| 乱人伦中国视频| 国产有黄有色有爽视频| 国产高清videossex| 侵犯人妻中文字幕一二三四区| 亚洲精品一二三| 青草久久国产| av电影中文网址| 久久人人精品亚洲av| 天堂中文最新版在线下载| 99久久综合精品五月天人人| 亚洲精品在线美女| 亚洲狠狠婷婷综合久久图片| 欧美另类亚洲清纯唯美| 999久久久精品免费观看国产| 十八禁网站免费在线| 在线观看日韩欧美| 又黄又粗又硬又大视频| 丰满的人妻完整版| 99久久国产精品久久久| www日本在线高清视频| 在线观看一区二区三区| 99riav亚洲国产免费| 亚洲情色 制服丝袜| 女人被躁到高潮嗷嗷叫费观| 91av网站免费观看| 最新美女视频免费是黄的| 国产无遮挡羞羞视频在线观看| 国产一区二区三区综合在线观看| 女人精品久久久久毛片| 亚洲熟女毛片儿| 长腿黑丝高跟| 狂野欧美激情性xxxx| 777久久人妻少妇嫩草av网站| 亚洲视频免费观看视频| 热99re8久久精品国产| 9191精品国产免费久久| 国产又色又爽无遮挡免费看| 久久午夜亚洲精品久久| 美国免费a级毛片| 国产又色又爽无遮挡免费看| 99香蕉大伊视频| 亚洲av成人一区二区三| 久久精品国产亚洲av香蕉五月| 18禁裸乳无遮挡免费网站照片 | 午夜免费激情av| 人妻丰满熟妇av一区二区三区| 国产精品影院久久| 这个男人来自地球电影免费观看| 色综合婷婷激情| 99re在线观看精品视频| 人人妻人人澡人人看| 亚洲va日本ⅴa欧美va伊人久久| 婷婷六月久久综合丁香| 日本vs欧美在线观看视频| 亚洲一区中文字幕在线| 美女高潮到喷水免费观看| 亚洲免费av在线视频| 欧美日韩精品网址| 淫秽高清视频在线观看| av天堂久久9| 男人舔女人的私密视频| 亚洲一区二区三区欧美精品| 亚洲欧美一区二区三区久久| 免费在线观看完整版高清| bbb黄色大片| 免费在线观看视频国产中文字幕亚洲| 人人妻人人爽人人添夜夜欢视频| 色综合欧美亚洲国产小说| 夜夜看夜夜爽夜夜摸 | 精品一区二区三区av网在线观看| 亚洲国产精品sss在线观看 | 亚洲精品成人av观看孕妇| 亚洲少妇的诱惑av| 麻豆成人av在线观看| 桃红色精品国产亚洲av| 一进一出好大好爽视频| 成人国产一区最新在线观看| 亚洲精品一卡2卡三卡4卡5卡| 欧美日韩精品网址| 亚洲精品一区av在线观看| 亚洲精品久久成人aⅴ小说| 99国产精品免费福利视频| av有码第一页| 国产成人精品无人区| 久久精品亚洲av国产电影网| 成人国语在线视频| 12—13女人毛片做爰片一| 日本a在线网址| 人人澡人人妻人| 亚洲美女黄片视频| 国产极品粉嫩免费观看在线| 亚洲欧美激情在线| 久久99一区二区三区| 亚洲精品久久午夜乱码| 在线视频色国产色| 欧美日本亚洲视频在线播放| 老司机午夜十八禁免费视频| 亚洲成av片中文字幕在线观看| 日日爽夜夜爽网站| 欧美乱码精品一区二区三区| 757午夜福利合集在线观看| 亚洲成人久久性| 成人国产一区最新在线观看| 国产99白浆流出| 亚洲欧美激情在线| 两性午夜刺激爽爽歪歪视频在线观看 | 男女做爰动态图高潮gif福利片 | 男人舔女人的私密视频| 亚洲精品美女久久久久99蜜臀| 999精品在线视频| 男女做爰动态图高潮gif福利片 | 天堂影院成人在线观看| 巨乳人妻的诱惑在线观看| 男人操女人黄网站| 亚洲欧美日韩另类电影网站| 国产成人精品无人区| 午夜日韩欧美国产| 免费高清视频大片| 成人18禁高潮啪啪吃奶动态图| 天天躁狠狠躁夜夜躁狠狠躁| 久久久精品国产亚洲av高清涩受| 国产无遮挡羞羞视频在线观看| 久久久久久久精品吃奶| 欧美 亚洲 国产 日韩一| www.999成人在线观看| 校园春色视频在线观看| 国产在线观看jvid| 色婷婷av一区二区三区视频| 嫩草影视91久久| 91国产中文字幕| 午夜激情av网站| www.999成人在线观看| 日韩 欧美 亚洲 中文字幕| 色婷婷av一区二区三区视频| 一级作爱视频免费观看| 男人的好看免费观看在线视频 | 大陆偷拍与自拍| 国产有黄有色有爽视频| 久久精品国产亚洲av高清一级| 69av精品久久久久久| 午夜精品久久久久久毛片777| 女人爽到高潮嗷嗷叫在线视频| 大码成人一级视频| 色播在线永久视频| 在线观看66精品国产| 亚洲三区欧美一区| 精品久久久久久久久久免费视频 | 欧美日韩av久久| 国产精品一区二区免费欧美| 热re99久久精品国产66热6| 久久久国产欧美日韩av| 中文字幕另类日韩欧美亚洲嫩草| netflix在线观看网站| 美女高潮喷水抽搐中文字幕| 最新美女视频免费是黄的| 成人三级黄色视频| 两人在一起打扑克的视频| 一区二区三区激情视频| 中亚洲国语对白在线视频| 国产伦一二天堂av在线观看| 国产精品av久久久久免费| 欧美最黄视频在线播放免费 | 99久久综合精品五月天人人| 欧美日韩国产mv在线观看视频| 看黄色毛片网站| 精品无人区乱码1区二区| 国产精品电影一区二区三区| 美国免费a级毛片| 午夜精品国产一区二区电影| 久久精品国产清高在天天线| videosex国产| 99riav亚洲国产免费| 国产麻豆69| 黄色毛片三级朝国网站| 国产午夜精品久久久久久| 91大片在线观看| 午夜福利一区二区在线看| 两人在一起打扑克的视频| 亚洲 国产 在线| 久久九九热精品免费| 亚洲欧美日韩另类电影网站| www.www免费av| 亚洲一区中文字幕在线| cao死你这个sao货| 欧美性长视频在线观看| 在线看a的网站| 成人精品一区二区免费| 99在线人妻在线中文字幕| 如日韩欧美国产精品一区二区三区| 男人操女人黄网站| 午夜福利,免费看| 丁香欧美五月| 久久中文字幕一级| 极品人妻少妇av视频| 亚洲精品国产色婷婷电影| 免费日韩欧美在线观看| 国产精品免费一区二区三区在线| 两人在一起打扑克的视频| 12—13女人毛片做爰片一| 如日韩欧美国产精品一区二区三区| 亚洲三区欧美一区| 亚洲人成电影观看| 夜夜看夜夜爽夜夜摸 | 国产精品一区二区在线不卡| 99热只有精品国产| 国产高清videossex| 少妇被粗大的猛进出69影院| 最近最新中文字幕大全电影3 | 国产精品免费视频内射| 亚洲中文av在线| 夜夜夜夜夜久久久久| 国产精品二区激情视频| 亚洲avbb在线观看| 淫妇啪啪啪对白视频| 日本黄色视频三级网站网址| 五月开心婷婷网| 亚洲免费av在线视频| 99久久99久久久精品蜜桃| 亚洲精品久久成人aⅴ小说| 99国产精品一区二区蜜桃av| av网站在线播放免费| 免费日韩欧美在线观看| 黄色视频不卡| 夜夜爽天天搞| 可以免费在线观看a视频的电影网站| 激情在线观看视频在线高清| 国产精品久久久av美女十八| 国产精华一区二区三区| 一级片'在线观看视频| 激情视频va一区二区三区| av在线天堂中文字幕 | 成人av一区二区三区在线看| 男女下面插进去视频免费观看| 久99久视频精品免费| 激情视频va一区二区三区| 丰满的人妻完整版| 久久久国产成人精品二区 | 国产又色又爽无遮挡免费看| 午夜福利,免费看| 欧美中文综合在线视频| 国产无遮挡羞羞视频在线观看| 91av网站免费观看| 国产日韩一区二区三区精品不卡| 天天躁夜夜躁狠狠躁躁| 久久人妻福利社区极品人妻图片| 人妻久久中文字幕网| 97人妻天天添夜夜摸| 国产欧美日韩精品亚洲av| 一进一出抽搐动态| 一二三四在线观看免费中文在| a级片在线免费高清观看视频| 国产精品 欧美亚洲| 另类亚洲欧美激情| 国产精品一区二区在线不卡| 啦啦啦 在线观看视频| 免费一级毛片在线播放高清视频 | 亚洲五月天丁香| 午夜两性在线视频| 久久亚洲真实| 国产精品电影一区二区三区| 久久精品91蜜桃| 免费日韩欧美在线观看| 亚洲人成77777在线视频| 免费观看精品视频网站| 国产精品久久久人人做人人爽| 操美女的视频在线观看| 97碰自拍视频| 国产精品秋霞免费鲁丝片| 亚洲三区欧美一区| 男女下面插进去视频免费观看| 亚洲成人久久性| 在线观看免费视频网站a站| 国产成人av教育| 国产精品99久久99久久久不卡| 亚洲av日韩精品久久久久久密| 黄色a级毛片大全视频| 黄色视频不卡| 亚洲专区字幕在线| 国产精品久久久久久人妻精品电影| 在线观看66精品国产| 欧美精品一区二区免费开放| 18禁观看日本| 首页视频小说图片口味搜索| 国产精品香港三级国产av潘金莲| 制服人妻中文乱码| 日韩三级视频一区二区三区| 日韩免费高清中文字幕av| 久久久久久久精品吃奶| 亚洲精品在线美女| 日韩一卡2卡3卡4卡2021年| 国产又爽黄色视频| 在线播放国产精品三级| 黄色视频不卡| 欧美 亚洲 国产 日韩一| 一级,二级,三级黄色视频| 亚洲一区二区三区色噜噜 | 一区二区三区国产精品乱码| 黄频高清免费视频| av免费在线观看网站| 欧美乱妇无乱码| 纯流量卡能插随身wifi吗| 日韩一卡2卡3卡4卡2021年| 99re在线观看精品视频| 亚洲情色 制服丝袜| 视频区欧美日本亚洲| 热99国产精品久久久久久7| 久久久久国内视频| 欧美日韩亚洲高清精品| 女人精品久久久久毛片| 一进一出好大好爽视频| xxx96com| 国产av一区二区精品久久| 国产精品久久视频播放| 日本vs欧美在线观看视频| 欧美日韩精品网址| 亚洲第一青青草原| 亚洲第一欧美日韩一区二区三区| 国产精品亚洲一级av第二区| 国产精品一区二区三区四区久久 | 国产97色在线日韩免费| 黄色怎么调成土黄色| 国产精品国产av在线观看| 亚洲免费av在线视频| 日本五十路高清| 在线观看免费午夜福利视频| 精品一区二区三区视频在线观看免费 | 99国产精品免费福利视频| 久久午夜综合久久蜜桃| 男人舔女人下体高潮全视频| 欧美丝袜亚洲另类 | 国产一区二区激情短视频| 脱女人内裤的视频| 深夜精品福利| 色尼玛亚洲综合影院| 日本欧美视频一区| 人人妻人人澡人人看| 999精品在线视频| 18禁裸乳无遮挡免费网站照片 | 久久精品国产亚洲av高清一级| 国产成人免费无遮挡视频| 久久中文看片网| 伦理电影免费视频| 国产xxxxx性猛交| 精品无人区乱码1区二区| 一进一出抽搐gif免费好疼 | 免费在线观看影片大全网站| 99riav亚洲国产免费| 欧美成人午夜精品| 女警被强在线播放| 伊人久久大香线蕉亚洲五| 一进一出抽搐动态| 一区在线观看完整版| 免费日韩欧美在线观看| 国产精品av久久久久免费| 99久久99久久久精品蜜桃| av在线天堂中文字幕 | 久久久久久免费高清国产稀缺| 男人操女人黄网站| 亚洲三区欧美一区| 国产精华一区二区三区| 丰满饥渴人妻一区二区三| 国产精品98久久久久久宅男小说| 国产一区二区三区在线臀色熟女 | 亚洲精品中文字幕在线视频| 男女下面插进去视频免费观看| 国产精品久久久久成人av| 久久中文字幕一级| 我的亚洲天堂| 国产精品久久久久久人妻精品电影| 亚洲在线自拍视频| 丰满的人妻完整版| 国产日韩一区二区三区精品不卡| 亚洲视频免费观看视频| 国产av一区在线观看免费| 一级a爱片免费观看的视频| 大码成人一级视频| 日韩欧美一区视频在线观看| av国产精品久久久久影院| 久久中文字幕人妻熟女| 一级黄色大片毛片| 97超级碰碰碰精品色视频在线观看| 两个人免费观看高清视频| 精品高清国产在线一区| 久久精品成人免费网站| 韩国精品一区二区三区| 欧美中文综合在线视频| 久久久久久大精品| 欧美日韩av久久| 人人妻人人澡人人看| 国产免费av片在线观看野外av| 亚洲国产欧美网| 99riav亚洲国产免费| 午夜影院日韩av| 精品久久久久久久久久免费视频 | 美女扒开内裤让男人捅视频| 女性生殖器流出的白浆| 久久国产精品影院| 精品一区二区三区四区五区乱码| 久热这里只有精品99| 精品卡一卡二卡四卡免费| 18禁裸乳无遮挡免费网站照片 | 欧美国产精品va在线观看不卡| 91国产中文字幕| 天天影视国产精品| 看免费av毛片| 一级片免费观看大全| 91字幕亚洲| 级片在线观看| 少妇裸体淫交视频免费看高清 | 18禁裸乳无遮挡免费网站照片 | 色老头精品视频在线观看| 久久久国产成人免费| 99re在线观看精品视频| 中文字幕高清在线视频| 麻豆av在线久日| 亚洲美女黄片视频| 精品国内亚洲2022精品成人| 69av精品久久久久久| 黑人巨大精品欧美一区二区蜜桃| 搡老乐熟女国产| 免费女性裸体啪啪无遮挡网站| 999久久久精品免费观看国产| 老司机福利观看| 久久精品人人爽人人爽视色| 中文亚洲av片在线观看爽| 亚洲免费av在线视频| 长腿黑丝高跟| 日韩精品中文字幕看吧| 日本wwww免费看| 免费在线观看完整版高清| 欧美日韩瑟瑟在线播放| 精品国产一区二区久久| 岛国视频午夜一区免费看| 亚洲一区中文字幕在线| x7x7x7水蜜桃| a级片在线免费高清观看视频| 热99国产精品久久久久久7| 18禁黄网站禁片午夜丰满| 黄色女人牲交| 中文字幕精品免费在线观看视频| 后天国语完整版免费观看| 一a级毛片在线观看| 91在线观看av| 国产色视频综合| av欧美777| 亚洲午夜理论影院| 在线观看免费日韩欧美大片| 国产不卡一卡二| 在线观看舔阴道视频| 国产欧美日韩一区二区精品| 久久精品国产亚洲av香蕉五月| 国产伦一二天堂av在线观看| av福利片在线| 国产蜜桃级精品一区二区三区| 老司机深夜福利视频在线观看| bbb黄色大片| 99精品欧美一区二区三区四区| 欧美精品亚洲一区二区|