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

    DNDC模型在華北平原冬小麥區(qū)的率定和驗(yàn)證

    2018-07-09 12:32:04喬帥帥胡振華張寶忠
    中國(guó)農(nóng)村水利水電 2018年6期
    關(guān)鍵詞:土壤水分冬小麥生物量

    喬帥帥,胡振華,魏 征,張寶忠

    (1.山西農(nóng)業(yè)大學(xué)林學(xué)院,山西 晉中 030801;2.中國(guó)水利水電科學(xué)研究院 流域水循環(huán)模擬與調(diào)控國(guó)家重點(diǎn)實(shí)驗(yàn)室,北京 100038;3.國(guó)家節(jié)水灌溉北京工程技術(shù)研究中心,北京 100048)

    目前國(guó)內(nèi)外學(xué)者研究和開(kāi)發(fā)了各種作物模型,如AquaCrop模型、DSSAT模型、DNDC模型等。這些模型被應(yīng)用在世界各地,與當(dāng)?shù)赝寥馈夂?、耕作方式等結(jié)合,逐漸實(shí)現(xiàn)本土化,也使模型本身得到豐富和更好地利用。通過(guò)對(duì)作物生長(zhǎng)的模擬,可以預(yù)測(cè)作物生長(zhǎng)、器官建成、產(chǎn)量及品質(zhì)形成,提供優(yōu)化動(dòng)態(tài)目標(biāo),為實(shí)施調(diào)控與決策提供依據(jù)[1]。作物模型從系統(tǒng)角度描述作物及其環(huán)境之間的關(guān)系,具有系統(tǒng)性、機(jī)理性和效益性等優(yōu)點(diǎn)[2]。而在這些模型當(dāng)中,基于過(guò)程的生物地球化學(xué)模型—DNDC 模型具有模擬功能強(qiáng)大、輸入?yún)?shù)容易獲得、軟件界面友好等優(yōu)點(diǎn),在世界范圍內(nèi)被廣泛應(yīng)用[3]。

    DNDC模型[4-8](Denitrification-Decomposition model,脫氮-分解作用模型)是美國(guó)新罕布什爾大學(xué)陸地海洋空間研究中心李長(zhǎng)生教授等開(kāi)發(fā)研制并推廣起來(lái)的。DNDC模型通過(guò)模擬生態(tài)系統(tǒng)光合、呼吸、有機(jī)碳生產(chǎn)等,可以用來(lái)模擬預(yù)測(cè)作物生產(chǎn)過(guò)程、土壤水分狀況、土壤碳,氮的釋放與淋濕,也包括微量氣體的排放等。DNDC模型自1992年間問(wèn)世以來(lái),已在世界許多地區(qū)得到驗(yàn)證并開(kāi)展了應(yīng)用研究[9,10],在我國(guó)也有驗(yàn)證和相關(guān)方面的應(yīng)用研究[11-13]。近年來(lái)DNDC模型越來(lái)越多地被應(yīng)用于作物生長(zhǎng)過(guò)程中土壤水分、生物量和產(chǎn)量的模擬。陳海心[14]等,應(yīng)用DNDC模型模擬關(guān)中地區(qū)作物產(chǎn)量、模型預(yù)測(cè)精度較好,可用于對(duì)產(chǎn)量的模擬。劉寧[15]等利用DNDC模型模擬了溝壟覆膜種植對(duì)土壤水熱及冬小麥產(chǎn)量的影響。韓娟[16]利用DNDC模型模擬陜西省玉米產(chǎn)量及土壤水熱效應(yīng),模擬精度較好。在華北平原區(qū),不同灌溉制度下,冬小麥生長(zhǎng)期間的土壤水分、生物量和產(chǎn)量的結(jié)果差距較大,而DNDC模型在不同灌溉制度下的冬小麥土壤水分、生物量和產(chǎn)量的模擬與驗(yàn)證研究尚未見(jiàn)報(bào)道。

    1 材料與方法

    1.1 研究區(qū)概括

    試驗(yàn)在中國(guó)水利水電科學(xué)研究院大興區(qū)國(guó)家節(jié)水灌溉中心完成(見(jiàn)圖1)。該試驗(yàn)基地地處北京南約30 km,北緯39°37.25′,東經(jīng)116°25.51′,屬于半干旱大陸性季風(fēng)氣候,冬季寒冷少雨,夏季炎熱多雨,多年平均降雨量540 mm,降雨多集中在6-9月。該區(qū)全年大于10 ℃的有效積溫為4 730 ℃,全年日照時(shí)數(shù)為2 600 h。易形成冬春干旱,夏季多雨,旱澇交替的氣候特點(diǎn)。該區(qū)土壤類(lèi)型以沙壤土為主。

    圖1 大興試驗(yàn)站地理位置示意圖

    1.2 數(shù)據(jù)獲取

    1.2.1 管理數(shù)據(jù)

    試驗(yàn)間供試小麥為中麥175, 2014-2015年冬小麥播種日期為2014年10月11日,收獲日期為2015年6月11日。 2015-2016年冬小麥播種日期為2015年10月18日,收獲日期為2016年6月11日。實(shí)驗(yàn)統(tǒng)一設(shè)計(jì)為3個(gè)水分處理,3個(gè)重復(fù),共9個(gè)小區(qū),每個(gè)小區(qū)約為30 m2。所有小區(qū)的施肥、耕作措施與當(dāng)?shù)剞r(nóng)民習(xí)慣一致。不同年間灌溉制度如表1所示。

    表1 冬小麥灌溉制度 mm

    1.2.2 氣象數(shù)據(jù)

    在冬小麥種植期間,氣象數(shù)據(jù)通過(guò)實(shí)驗(yàn)站自動(dòng)氣象觀測(cè)站監(jiān)測(cè),主要監(jiān)測(cè)內(nèi)容包括降雨、風(fēng)速、氣溫、地溫、空氣濕度和太陽(yáng)輻射等,每隔半小時(shí)自動(dòng)采集數(shù)據(jù)一次。

    1.2.3 土壤水分?jǐn)?shù)據(jù)

    土壤水分通過(guò)TRIME-DATA PILOT(德國(guó)IMKO公司)土壤水分測(cè)量系統(tǒng)進(jìn)行測(cè)定,小區(qū)設(shè)有兩根trime 管,每隔3 d測(cè)量一次土壤體積含水率(0~50 cm),降雨和灌水后需加測(cè)一次。

    自覺(jué)接受黨委領(lǐng)導(dǎo)。各鎮(zhèn)(街道)人大自覺(jué)接受黨(工)委對(duì)“雙聯(lián)”工作的領(lǐng)導(dǎo),定期向黨(工)委請(qǐng)示匯報(bào)“雙聯(lián)”工作。各鎮(zhèn)(街道)黨(工)委將“雙聯(lián)”工作列入重要議程,將人大工作納入黨委總體工作部署,在辦公場(chǎng)所、辦公設(shè)備、人員配備、活動(dòng)經(jīng)費(fèi)等方面給予全力支持。

    1.2.4 地上部生物量和產(chǎn)量數(shù)據(jù)

    首先測(cè)定每個(gè)小區(qū)寬度和小區(qū)冬小麥行數(shù),在冬小麥返青后,以每隔10 d為間隔對(duì)冬小麥生物量進(jìn)行測(cè)定。每個(gè)小區(qū)隨機(jī)割取20 cm長(zhǎng)的冬小麥,從莖基部剪下,獲得完整冠部,統(tǒng)計(jì)20 cm株數(shù),按小區(qū)編號(hào),殺青(105 ℃)30 min,恒溫(75 ℃)24 h烘至恒重后稱(chēng)重。收割前在各小區(qū)取1 m2的典型樣本,計(jì)算每平方米株數(shù)及千粒重,在樣本中取30株,計(jì)算穗粒數(shù),并計(jì)算樣本內(nèi)小麥產(chǎn)量。

    1.3 DNDC模型率定及精度評(píng)價(jià)指標(biāo)

    1.3.1 參數(shù)敏感性分析

    利用獨(dú)立參數(shù)擾動(dòng)法,在冬小麥不同灌溉制度下,以生物量和產(chǎn)量為目標(biāo)對(duì)DNDC模型中的作物參數(shù)進(jìn)行敏感性分析,確定對(duì)生物量和產(chǎn)量敏感參數(shù)。假定各參數(shù)之間相互獨(dú)立,對(duì)各參數(shù)值進(jìn)行增加和減少10%,采用相對(duì)敏感度RS表示敏感性大小。

    1.3.2 模型參數(shù)率定

    對(duì)模型參數(shù)的敏感性進(jìn)行識(shí)別,進(jìn)而重點(diǎn)對(duì)模型的敏感參數(shù)進(jìn)行敏感性分析,是確定模型關(guān)鍵參數(shù)、控制模型效率非常有效的途徑[17]。應(yīng)用DNDC模型模擬任一點(diǎn)的生物地球化學(xué)過(guò)程時(shí),只要根據(jù)當(dāng)?shù)胤N植耕作情況輸入氣象、土壤以及作物等數(shù)據(jù)便可進(jìn)行一年至多年的模擬[14]。選擇2014年10月—2015年6月冬小麥種植期間的觀測(cè)數(shù)據(jù)對(duì)DNDC模型進(jìn)行率定,最后以2015年10月—2016年6月冬小麥觀測(cè)數(shù)據(jù)進(jìn)行驗(yàn)證,得到模擬大興區(qū)冬小麥生長(zhǎng)的DNDC模型參數(shù)。率定后作物參數(shù)如表2所示。

    表2 模型中冬小麥作物參數(shù)

    由DNDC模型模擬出來(lái)的生物量、產(chǎn)量和土壤水分與觀測(cè)值之間的擬合程度選用歸一化均方根誤差(NRMSE)和Willmott[18]的一致性指標(biāo)d。根據(jù)Gjettermanna[19]提出歸一化均方根誤差(NRMSE)小于25%表明模擬效果較好,25%~30%為可接受的范圍;d∈(-∞,1),d越接近于1,表明模型模擬結(jié)果越好。以上模型評(píng)價(jià)指標(biāo)具體計(jì)算公式為:

    (1)

    (2)

    2 結(jié)果與分析

    2.1 模型敏感性分析

    在對(duì)DNDC模型中冬小麥作物參數(shù)進(jìn)行敏感性分析后,分析結(jié)果(表3)表明,不同作物參數(shù)之間對(duì)模型結(jié)果影響不盡相同。在各項(xiàng)參數(shù)值改變中,作物需水量的整體改變對(duì)產(chǎn)量和生物量的影響較大,敏感度在0.626以上。籽粒生物量C/N比和最適溫度的減少對(duì)地上部生物量的影響不大,而籽粒生物量C/N比和最適溫度的增加會(huì)對(duì)地上部生物量產(chǎn)生較大的變化,敏感度在0.771~1.354之間。生長(zhǎng)積溫整體改變對(duì)生物量結(jié)果影響較小,對(duì)產(chǎn)量結(jié)果影響較大,敏感度在0.789~2.406之間。在對(duì)模型進(jìn)行模擬過(guò)程中需保證參數(shù)輸入精確性。

    表3 敏感性分析結(jié)果

    2.2 模型驗(yàn)證

    DNDC模型需要輸入的參數(shù)包括:氣候、土壤、農(nóng)作物和農(nóng)田管理等。本實(shí)驗(yàn)通過(guò)將模型原有參數(shù)與2015年冬小麥生育期實(shí)測(cè)數(shù)據(jù)結(jié)合,以生物量和產(chǎn)量為目標(biāo)對(duì)模型參數(shù)進(jìn)行敏感性分析,并進(jìn)行校正。由校正后的模型對(duì)2016年的冬小麥的土壤水分、地上部生物量和產(chǎn)量進(jìn)行模擬,并與實(shí)測(cè)結(jié)果進(jìn)行比較。最后通過(guò)模型評(píng)價(jià)指標(biāo)對(duì)模型精度進(jìn)行驗(yàn)證。

    2.2.1 冬小麥土壤含水率模擬

    對(duì)2016年不同灌溉制度下冬小麥土壤水分驗(yàn)證(圖2)可知,不同階段、灌水處理下下冬小麥土壤含水率變化差異較大,DNDC模型在冬小麥的返青期以后模擬效果較好,不同灌水處理下土壤含水率均在可接受的范圍內(nèi),決定系數(shù)R2分別為0.691、0.596和0.677。歸一均方根誤差NRMSE分別為13.5%、14.2%和11.1%,均小于14.0%,處在模擬較好范圍內(nèi)。一致性指數(shù)d分別為0.810、0.960和0.738。這些評(píng)價(jià)指標(biāo)表明DNDC模型可以較好的模擬不同灌溉制度下冬小麥生育期的土壤水分變化。

    2.2.2 冬小麥生物量模擬

    在對(duì)2016冬小麥的生物量模擬過(guò)程中,由圖3可知,隨時(shí)間變化,模擬值與實(shí)測(cè)值一致性較好,在生物量積累過(guò)程中模擬值均大于實(shí)際生物量值,但都處于可接受的范圍內(nèi)。統(tǒng)計(jì)結(jié)果表明,不同水分處理下,模型模擬值與實(shí)測(cè)值的決定系數(shù)R2分別為0.980、0.978和0.923,均大于0.920。歸一均方根誤差NRMSE分別為21%、18%和18%,一致性指數(shù)d分別為0.76、0.80和0.96。這些指標(biāo)值較好,表明DNDC模型對(duì)于生物量的模擬效果比較理想。

    圖2 2016年不同灌溉處理下土壤含水率模擬值與實(shí)測(cè)值比較

    圖3 2016年不同灌溉處理下生物量模擬值與實(shí)測(cè)值比較

    2.2.3 冬小麥產(chǎn)量結(jié)果模擬

    由DNDC模型產(chǎn)量模擬值與實(shí)測(cè)值(表4)可知:隨著灌水量的增加,產(chǎn)量也在不斷增加。當(dāng)灌溉定額為220 mm時(shí),截止到測(cè)產(chǎn)時(shí)產(chǎn)量結(jié)果差異不大,模擬結(jié)果與實(shí)測(cè)結(jié)果相吻合。模型模擬值與實(shí)測(cè)值相對(duì)誤差在18%以下,表明DNDC模型對(duì)于冬小麥產(chǎn)量模擬結(jié)果較好。

    表4 DNDC模型產(chǎn)量模擬值與實(shí)測(cè)值

    3 結(jié) 語(yǔ)

    本文通過(guò)模型敏感性分析、模型的校正和驗(yàn)證試驗(yàn)構(gòu)建了華北地區(qū)冬小麥的DNDC模型數(shù)據(jù)庫(kù)。利用國(guó)家節(jié)水灌溉北京工程技術(shù)研究中心大興試驗(yàn)站2015年冬小麥生長(zhǎng)數(shù)據(jù)對(duì)DNDC模型參數(shù)進(jìn)行了率定,并應(yīng)用于2016年冬小麥,與實(shí)測(cè)數(shù)據(jù)進(jìn)行驗(yàn)證,對(duì)模型精度進(jìn)行評(píng)價(jià)。具體結(jié)果如下:

    (1)通過(guò)模型敏感性分析可知,作物參數(shù)中作物需水量對(duì)地上部生物量和產(chǎn)量結(jié)果影響較大。

    (2)DNDC模型經(jīng)過(guò)驗(yàn)證后,可以較為準(zhǔn)確的模擬華北地區(qū)不同灌溉制度下冬小麥的土壤含水量、地上部生物量和產(chǎn)量的變化。

    (3)對(duì)于產(chǎn)量的模擬,當(dāng)灌水量達(dá)到一定量時(shí),DNDC模型模擬產(chǎn)量值差異變化不明顯,需要進(jìn)一步對(duì)DNDC模型進(jìn)行改進(jìn)。

    參考文獻(xiàn):

    [1] 李三愛(ài), 居 輝, 池寶亮. 作物生產(chǎn)潛力研究進(jìn)展[J]. 中國(guó)農(nóng)業(yè)氣象, 2005,26(2):106-111.

    [2] 楊 寧, 廖桂平. 作物生長(zhǎng)模擬研究進(jìn)展[J]. 作物研究, 2002,(5):255-257.

    [3] 王立剛,邱建軍,馬永良,等.應(yīng)用DNDC 模型分析施肥與翻耕方式對(duì)土壤有機(jī)碳含量的長(zhǎng)期影響[J]. 中國(guó)農(nóng)業(yè)大學(xué)學(xué)報(bào), 2004, 9(6):15-19.

    [4] Pathak H, Li C,Wassmann R.Greenhouse gas emissions from Indian rice fields: calibration and upscaling using the DNDC model[J]. Biogeosciences, 2005,2(2):113-123.

    [5] Tonitto C, David M B,Drinkwater L E, et al. Application of the DNDC model to tile-drained Illinois agroecosystems: model calibration, validation, and uncertainty analysis[J]. Nutrient Cycling in Agroecosystems, 2007,78(1):51-63.

    [6] Smith W N, Desjardins R L,Grant B, et al. Testing the DNDC model using N2O emissions at two experimental sites in Canade[J]. Canadian Journal of Soil Science, 2002,82(3):365-374.

    [7] Giltrap D, Saggar S, Li C, et al. Using the NZ-DNDC model to estimate agricultural N2O emissions in the Manawatu-Wanganui region[J]. Plant & Soil, 2008,309(1/2):191-209.

    [8] Cai Z, Sawamoto T, Li C, et al. Field validation of the DNDC model for greenhouse gas emissions in East Asian cropping systems[J]. Global Biogeochemical Cycles, 2003,17(4):1-11.

    [9] Smith P, Smith J U, Powlson D S, et al. A comparison of performance of nine soil organic matter models using datasets from seven long-term experiments[J]. Geoderma, 1997,81:153-225.

    [10] Li C, Frolking S, Crocker G J. Simulating trends in soil organic carbon in long-term experiments using the DNDC model[J]. Geoderma, 1997,81:45-60.

    [11] 李長(zhǎng)生.土壤碳儲(chǔ)量減少:中國(guó)農(nóng)業(yè)之隱患----中美農(nóng)業(yè)生態(tài)系統(tǒng)碳循環(huán)對(duì)比研究[J]. 第四紀(jì)研究, 2000,20(4):345-350.

    [12] 邱建軍,唐華俊,陳慶沐,等.區(qū)域農(nóng)業(yè)生態(tài)系統(tǒng)碳氮平衡研究以國(guó)家尺度為例[M]. 北京:中國(guó)人口出版社,2002:298-316.

    [13] 徐文彬, 洪業(yè)湯, 陳旭暉, 等.區(qū)域農(nóng)業(yè)土壤N2O釋放研究----以貴州省為例[J]. 中國(guó)科學(xué)(D 輯), 1999,29(5):450-456

    [14] 陳海心, 孫本華, 馮 浩,等. 應(yīng)用DNDC模型模擬關(guān)中地區(qū)農(nóng)田長(zhǎng)期施肥條件下土壤碳含量及作物產(chǎn)量[J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào), 2014,33(9):1 782-1 790.

    [15] 劉 寧,韓 娟,劉璐璐,等.應(yīng)用DNDC模型模擬溝壟集雨種植對(duì)陜西省玉米產(chǎn)量的影響[J]. 西北農(nóng)業(yè)學(xué)報(bào), 2016,25(5):691-701.

    [16] 任小龍, 賈志寬, 陳小莉,等. 模擬降雨量下溝壟微型集雨種植玉米的水溫效應(yīng)[J]. 中國(guó)農(nóng)業(yè)科學(xué), 2008,41(1):70-77.

    [17] 席 慶,李兆富,羅 川.基于擾動(dòng)分析方法的AnnAGNPS模型水文水質(zhì)參數(shù)敏感性分析[J].環(huán)境科學(xué), 2014,35(5):1 773-1 780.

    [18] Willmott C J, Ackleson S G, Davis R E, et al. Statistics for the evaluation and comparison of models[J]. Journal of Geophysical Research, 1985,90(C5):8 995-9 005.

    [19] Gjettermann B, Styczen M, Hansen H C B. Challenges in modeling dissolved organic matter dynamics in agricultural soil using DAISY[J]. Soil Biology & Biochemistry, 2008,(40):1 506-1 518.

    猜你喜歡
    土壤水分冬小麥生物量
    輪牧能有效促進(jìn)高寒草地生物量和穩(wěn)定性
    甘肅冬小麥田
    西藏高原土壤水分遙感監(jiān)測(cè)方法研究
    生物量高的富鋅酵母的開(kāi)發(fā)應(yīng)用
    冬小麥和春小麥
    中學(xué)生(2015年4期)2015-08-31 02:53:50
    基于SPOT-5遙感影像估算玉米成熟期地上生物量及其碳氮累積量
    不同覆蓋措施對(duì)棗園土壤水分和溫度的影響
    冬小麥——新冬18號(hào)
    冬小麥—新冬41號(hào)
    植被覆蓋區(qū)土壤水分反演研究——以北京市為例
    80岁老熟妇乱子伦牲交| 精品国产乱子伦一区二区三区| 亚洲欧美一区二区三区黑人| 欧美+亚洲+日韩+国产| 国产高清国产精品国产三级| 一个人观看的视频www高清免费观看 | 久久热在线av| 午夜精品久久久久久毛片777| 国产精品爽爽va在线观看网站 | 日韩精品青青久久久久久| 不卡av一区二区三区| xxx96com| www.精华液| 性欧美人与动物交配| 18禁裸乳无遮挡免费网站照片 | 国产有黄有色有爽视频| 欧美乱色亚洲激情| 国产在线精品亚洲第一网站| 欧美成狂野欧美在线观看| 亚洲一区二区三区欧美精品| 精品免费久久久久久久清纯| 悠悠久久av| 欧美不卡视频在线免费观看 | 在线av久久热| 中国美女看黄片| 亚洲av电影在线进入| 一区二区日韩欧美中文字幕| 国产精品综合久久久久久久免费 | 性少妇av在线| 亚洲人成电影观看| 国产又爽黄色视频| 午夜视频精品福利| 国产成人av激情在线播放| 性少妇av在线| 桃红色精品国产亚洲av| 99精国产麻豆久久婷婷| 久久天堂一区二区三区四区| 神马国产精品三级电影在线观看 | 一边摸一边抽搐一进一出视频| 亚洲精品久久午夜乱码| 久久伊人香网站| 欧美黄色淫秽网站| 99精品久久久久人妻精品| 亚洲一区高清亚洲精品| 欧美人与性动交α欧美精品济南到| 久久99一区二区三区| 久久精品人人爽人人爽视色| 午夜免费观看网址| 一a级毛片在线观看| 人人妻,人人澡人人爽秒播| 国产精品亚洲一级av第二区| 色婷婷久久久亚洲欧美| 男男h啪啪无遮挡| 黄色丝袜av网址大全| www国产在线视频色| 久9热在线精品视频| 老鸭窝网址在线观看| 露出奶头的视频| 88av欧美| 午夜a级毛片| 国产野战对白在线观看| 亚洲精品中文字幕在线视频| 亚洲色图av天堂| 日韩免费高清中文字幕av| 亚洲 欧美一区二区三区| 欧美中文综合在线视频| 性少妇av在线| 欧美人与性动交α欧美精品济南到| 69av精品久久久久久| 嫩草影院精品99| 久久久水蜜桃国产精品网| 午夜福利一区二区在线看| 午夜影院日韩av| 欧美一区二区精品小视频在线| 精品福利观看| netflix在线观看网站| 欧美最黄视频在线播放免费 | 免费日韩欧美在线观看| 一级毛片高清免费大全| 91老司机精品| 国产精品永久免费网站| 自拍欧美九色日韩亚洲蝌蚪91| 欧美日韩一级在线毛片| 大香蕉久久成人网| 亚洲精华国产精华精| 久久精品国产亚洲av高清一级| 午夜亚洲福利在线播放| 91成年电影在线观看| 亚洲av成人av| avwww免费| 国产成人系列免费观看| 日韩精品青青久久久久久| 国产99白浆流出| 动漫黄色视频在线观看| 国产精品免费视频内射| 日韩国内少妇激情av| 色播在线永久视频| 中文字幕av电影在线播放| 免费av中文字幕在线| 夜夜躁狠狠躁天天躁| 看片在线看免费视频| 欧美中文综合在线视频| 久久草成人影院| 成人免费观看视频高清| 国产精品二区激情视频| 国产欧美日韩一区二区三| 男人舔女人下体高潮全视频| 亚洲国产精品sss在线观看 | 长腿黑丝高跟| 亚洲熟妇中文字幕五十中出 | 成熟少妇高潮喷水视频| 欧美中文日本在线观看视频| 女人精品久久久久毛片| 99久久国产精品久久久| 99精品久久久久人妻精品| 欧美日韩一级在线毛片| 久久午夜亚洲精品久久| 在线观看免费日韩欧美大片| a级毛片在线看网站| 在线看a的网站| 亚洲精品中文字幕一二三四区| 久久中文字幕人妻熟女| 99精国产麻豆久久婷婷| 国产一区二区三区在线臀色熟女 | 两个人看的免费小视频| 午夜成年电影在线免费观看| 男男h啪啪无遮挡| 九色亚洲精品在线播放| 岛国在线观看网站| 欧美 亚洲 国产 日韩一| 欧美 亚洲 国产 日韩一| 国产精品乱码一区二三区的特点 | 男男h啪啪无遮挡| 黄色毛片三级朝国网站| 可以在线观看毛片的网站| 亚洲欧美日韩高清在线视频| 亚洲欧美日韩高清在线视频| 精品久久久久久久久久免费视频 | 国产乱人伦免费视频| 一二三四在线观看免费中文在| 自线自在国产av| 自线自在国产av| 女人高潮潮喷娇喘18禁视频| 日韩精品青青久久久久久| 久久人人97超碰香蕉20202| 成人av一区二区三区在线看| 12—13女人毛片做爰片一| 丝袜在线中文字幕| 悠悠久久av| 伊人久久大香线蕉亚洲五| 国产精品影院久久| 午夜影院日韩av| 成人三级黄色视频| 亚洲国产精品sss在线观看 | 亚洲av电影在线进入| 黄色视频,在线免费观看| 琪琪午夜伦伦电影理论片6080| 久久久精品欧美日韩精品| 免费女性裸体啪啪无遮挡网站| 满18在线观看网站| 久久精品91无色码中文字幕| 超色免费av| 日本一区二区免费在线视频| 精品久久久久久电影网| 91麻豆av在线| 91大片在线观看| 精品国内亚洲2022精品成人| 看片在线看免费视频| 丝袜美足系列| 亚洲自偷自拍图片 自拍| 日韩精品中文字幕看吧| 久热爱精品视频在线9| 国产精品98久久久久久宅男小说| 欧美乱码精品一区二区三区| 久久 成人 亚洲| 久久久久国产精品人妻aⅴ院| 色老头精品视频在线观看| 精品少妇一区二区三区视频日本电影| 又紧又爽又黄一区二区| 一区二区三区激情视频| 亚洲国产欧美日韩在线播放| 天天躁夜夜躁狠狠躁躁| 天天躁狠狠躁夜夜躁狠狠躁| 午夜福利欧美成人| 国产一卡二卡三卡精品| 欧美激情 高清一区二区三区| 日本精品一区二区三区蜜桃| 国产精品九九99| 亚洲av日韩精品久久久久久密| 香蕉丝袜av| cao死你这个sao货| 午夜免费激情av| 视频区图区小说| 狂野欧美激情性xxxx| 欧美日本亚洲视频在线播放| 男女床上黄色一级片免费看| 美女高潮到喷水免费观看| 天堂√8在线中文| 五月开心婷婷网| 国产精品日韩av在线免费观看 | 午夜福利,免费看| a级毛片在线看网站| 国产成人欧美在线观看| 日韩欧美在线二视频| 两个人看的免费小视频| 日韩大尺度精品在线看网址 | 成人影院久久| 91大片在线观看| 一级毛片女人18水好多| 一级毛片高清免费大全| av电影中文网址| 亚洲av五月六月丁香网| 亚洲欧美精品综合一区二区三区| 变态另类成人亚洲欧美熟女 | www日本在线高清视频| 手机成人av网站| 日韩中文字幕欧美一区二区| 每晚都被弄得嗷嗷叫到高潮| 成人国产一区最新在线观看| 男女午夜视频在线观看| 在线观看一区二区三区| 99精国产麻豆久久婷婷| 高清毛片免费观看视频网站 | www.精华液| 国产在线精品亚洲第一网站| 看片在线看免费视频| 国产伦一二天堂av在线观看| 一边摸一边抽搐一进一出视频| 国产视频一区二区在线看| 少妇粗大呻吟视频| 波多野结衣高清无吗| 巨乳人妻的诱惑在线观看| 少妇被粗大的猛进出69影院| 国产成年人精品一区二区 | xxx96com| 看免费av毛片| 脱女人内裤的视频| 精品国产超薄肉色丝袜足j| 国产又爽黄色视频| 中文字幕另类日韩欧美亚洲嫩草| 人成视频在线观看免费观看| 欧美日韩一级在线毛片| 桃色一区二区三区在线观看| 熟女少妇亚洲综合色aaa.| 亚洲情色 制服丝袜| 精品卡一卡二卡四卡免费| 欧美日韩瑟瑟在线播放| 老熟妇仑乱视频hdxx| 免费在线观看完整版高清| 国产成人精品无人区| 人妻久久中文字幕网| 好看av亚洲va欧美ⅴa在| 午夜免费成人在线视频| 久久中文字幕人妻熟女| 精品国产一区二区久久| 看免费av毛片| 水蜜桃什么品种好| 成人三级黄色视频| 亚洲国产精品一区二区三区在线| 狂野欧美激情性xxxx| 国产成+人综合+亚洲专区| 国产不卡一卡二| 久久人妻熟女aⅴ| 亚洲免费av在线视频| 久久热在线av| 亚洲午夜理论影院| av免费在线观看网站| 亚洲精品一区av在线观看| 岛国视频午夜一区免费看| 国内毛片毛片毛片毛片毛片| 日韩高清综合在线| 很黄的视频免费| 亚洲va日本ⅴa欧美va伊人久久| 亚洲国产精品sss在线观看 | 久久天躁狠狠躁夜夜2o2o| 日本三级黄在线观看| 757午夜福利合集在线观看| 国产主播在线观看一区二区| 51午夜福利影视在线观看| 在线av久久热| 亚洲欧美一区二区三区久久| 纯流量卡能插随身wifi吗| 99re在线观看精品视频| 国产蜜桃级精品一区二区三区| 欧美日韩中文字幕国产精品一区二区三区 | 国产精品久久电影中文字幕| 国产精品久久久久成人av| 亚洲精品一区av在线观看| 日韩三级视频一区二区三区| 久久天躁狠狠躁夜夜2o2o| 夜夜夜夜夜久久久久| 人妻丰满熟妇av一区二区三区| 97人妻天天添夜夜摸| 亚洲美女黄片视频| 色综合站精品国产| aaaaa片日本免费| 国产精品乱码一区二三区的特点 | 精品一区二区三区视频在线观看免费 | 国产免费av片在线观看野外av| 日本黄色日本黄色录像| 国产精品美女特级片免费视频播放器 | 一级黄色大片毛片| 精品一区二区三卡| 久久久久久久精品吃奶| 在线播放国产精品三级| 精品国产超薄肉色丝袜足j| 亚洲熟妇熟女久久| 欧美日韩视频精品一区| 级片在线观看| 黑人巨大精品欧美一区二区mp4| 久久国产乱子伦精品免费另类| 亚洲欧美一区二区三区黑人| 97超级碰碰碰精品色视频在线观看| 国产亚洲精品综合一区在线观看 | 久久精品91蜜桃| 国产精品免费一区二区三区在线| 亚洲五月天丁香| 91大片在线观看| 男人舔女人的私密视频| 国产一区二区激情短视频| 欧美日韩一级在线毛片| 91老司机精品| 麻豆成人av在线观看| 国产xxxxx性猛交| 亚洲熟妇熟女久久| 香蕉国产在线看| 伦理电影免费视频| 欧美国产精品va在线观看不卡| 亚洲 国产 在线| 国产精品一区二区三区四区久久 | 99久久国产精品久久久| 我的亚洲天堂| 黄网站色视频无遮挡免费观看| 两个人免费观看高清视频| 亚洲欧美一区二区三区黑人| 首页视频小说图片口味搜索| 免费在线观看视频国产中文字幕亚洲| 久久欧美精品欧美久久欧美| 女生性感内裤真人,穿戴方法视频| 制服诱惑二区| 精品国内亚洲2022精品成人| 精品电影一区二区在线| 亚洲av第一区精品v没综合| 亚洲精品国产一区二区精华液| 亚洲欧美日韩另类电影网站| 国产高清激情床上av| 国产成人精品久久二区二区91| 精品一区二区三区四区五区乱码| 亚洲av成人一区二区三| 搡老岳熟女国产| 午夜免费激情av| 淫妇啪啪啪对白视频| 自线自在国产av| 精品免费久久久久久久清纯| 黄色丝袜av网址大全| 在线观看www视频免费| 日本精品一区二区三区蜜桃| 18禁黄网站禁片午夜丰满| 美女福利国产在线| 最好的美女福利视频网| 99国产精品一区二区蜜桃av| 男女高潮啪啪啪动态图| 国产亚洲av高清不卡| 岛国视频午夜一区免费看| 人人妻人人澡人人看| 国产免费男女视频| 国产主播在线观看一区二区| tocl精华| 国产成人免费无遮挡视频| 国产野战对白在线观看| 1024视频免费在线观看| 成人18禁高潮啪啪吃奶动态图| 午夜两性在线视频| 国产三级黄色录像| 国产成人一区二区三区免费视频网站| 国产成+人综合+亚洲专区| 88av欧美| 久久久精品国产亚洲av高清涩受| 女人爽到高潮嗷嗷叫在线视频| 亚洲自拍偷在线| 亚洲精品中文字幕一二三四区| 法律面前人人平等表现在哪些方面| 黄色成人免费大全| 亚洲精品美女久久久久99蜜臀| 国产欧美日韩综合在线一区二区| 丝袜美足系列| 欧美在线黄色| 男人的好看免费观看在线视频 | 欧美大码av| www.自偷自拍.com| 久久人妻福利社区极品人妻图片| 狂野欧美激情性xxxx| 欧美日韩精品网址| a级片在线免费高清观看视频| 欧美乱妇无乱码| 国产精品野战在线观看 | 一区二区三区精品91| 50天的宝宝边吃奶边哭怎么回事| 每晚都被弄得嗷嗷叫到高潮| 国产免费现黄频在线看| 99久久精品国产亚洲精品| 岛国在线观看网站| 亚洲久久久国产精品| 精品久久久久久电影网| 日韩国内少妇激情av| 亚洲成人国产一区在线观看| av福利片在线| 久久精品亚洲熟妇少妇任你| 18禁国产床啪视频网站| 在线观看免费午夜福利视频| 欧美日韩国产mv在线观看视频| 80岁老熟妇乱子伦牲交| 一级作爱视频免费观看| 又黄又粗又硬又大视频| 国产成年人精品一区二区 | 热re99久久精品国产66热6| 久久人妻福利社区极品人妻图片| 免费观看精品视频网站| 男女做爰动态图高潮gif福利片 | 国产精品九九99| 久久久久精品国产欧美久久久| 欧美一区二区精品小视频在线| 日韩欧美免费精品| 亚洲一区二区三区不卡视频| 国产又色又爽无遮挡免费看| 成年人黄色毛片网站| 亚洲va日本ⅴa欧美va伊人久久| 国产精品二区激情视频| 精品久久久久久成人av| 搡老岳熟女国产| 人人妻人人澡人人看| 国产一区二区三区视频了| 啪啪无遮挡十八禁网站| 最新美女视频免费是黄的| 亚洲 国产 在线| 中国美女看黄片| 成人永久免费在线观看视频| 国产精品一区二区三区四区久久 | 女性被躁到高潮视频| 每晚都被弄得嗷嗷叫到高潮| 精品一区二区三区四区五区乱码| 日本欧美视频一区| 久久影院123| 久久久水蜜桃国产精品网| 亚洲av片天天在线观看| 国产精品av久久久久免费| 久久久久国产精品人妻aⅴ院| 午夜福利在线观看吧| 不卡一级毛片| 99久久精品国产亚洲精品| 精品一区二区三卡| xxx96com| 精品久久久久久久毛片微露脸| 国产成+人综合+亚洲专区| 一个人观看的视频www高清免费观看 | 久久人人爽av亚洲精品天堂| av有码第一页| 老司机亚洲免费影院| 在线观看午夜福利视频| 日本免费一区二区三区高清不卡 | 搡老熟女国产l中国老女人| 国产成人精品久久二区二区91| 久久香蕉精品热| 亚洲人成电影免费在线| 在线观看一区二区三区激情| 免费久久久久久久精品成人欧美视频| 欧美黑人欧美精品刺激| 91九色精品人成在线观看| 精品福利永久在线观看| 身体一侧抽搐| 欧美丝袜亚洲另类 | 国产真人三级小视频在线观看| 精品无人区乱码1区二区| 久久精品91蜜桃| 1024视频免费在线观看| 岛国在线观看网站| 国产av精品麻豆| 一区二区三区国产精品乱码| 亚洲熟妇熟女久久| 成人18禁高潮啪啪吃奶动态图| 免费少妇av软件| 人人妻人人澡人人看| 女同久久另类99精品国产91| 99精国产麻豆久久婷婷| 黑人欧美特级aaaaaa片| 欧美+亚洲+日韩+国产| 亚洲av第一区精品v没综合| 免费在线观看日本一区| 日韩精品免费视频一区二区三区| 欧美日韩福利视频一区二区| 女性被躁到高潮视频| 精品久久久久久久毛片微露脸| 国产高清videossex| 女人高潮潮喷娇喘18禁视频| 国产黄色免费在线视频| 亚洲午夜理论影院| 久久中文看片网| 久久久国产成人精品二区 | 久久久精品欧美日韩精品| 国产成人av激情在线播放| 国产一卡二卡三卡精品| 日本三级黄在线观看| 国产真人三级小视频在线观看| 国产成人av教育| 999久久久国产精品视频| 中文字幕人妻熟女乱码| 自线自在国产av| 国产亚洲欧美精品永久| 国产一区二区三区综合在线观看| 欧美性长视频在线观看| 亚洲色图 男人天堂 中文字幕| 国产亚洲精品第一综合不卡| 国产aⅴ精品一区二区三区波| 亚洲av电影在线进入| 亚洲欧美精品综合一区二区三区| 人人妻人人澡人人看| 在线观看免费视频网站a站| 国产精品电影一区二区三区| 男女下面插进去视频免费观看| 色哟哟哟哟哟哟| 国产成人欧美在线观看| 日本精品一区二区三区蜜桃| 亚洲自偷自拍图片 自拍| 91国产中文字幕| 男女午夜视频在线观看| 精品第一国产精品| 国产亚洲精品久久久久5区| 日本黄色视频三级网站网址| 三级毛片av免费| 国产亚洲精品第一综合不卡| 老司机午夜福利在线观看视频| 久久国产精品男人的天堂亚洲| 两人在一起打扑克的视频| 日韩精品青青久久久久久| 男人的好看免费观看在线视频 | 91在线观看av| 日本精品一区二区三区蜜桃| 多毛熟女@视频| 免费看a级黄色片| 后天国语完整版免费观看| 一二三四在线观看免费中文在| 91av网站免费观看| 国产成人精品久久二区二区免费| 国产精品一区二区三区四区久久 | 亚洲一区二区三区不卡视频| 琪琪午夜伦伦电影理论片6080| 人人澡人人妻人| 亚洲,欧美精品.| 国产片内射在线| 日韩 欧美 亚洲 中文字幕| 欧美丝袜亚洲另类 | 国产又色又爽无遮挡免费看| av天堂久久9| 亚洲九九香蕉| 成人特级黄色片久久久久久久| 性欧美人与动物交配| www.999成人在线观看| 757午夜福利合集在线观看| av福利片在线| 午夜福利一区二区在线看| 9191精品国产免费久久| 国产精华一区二区三区| 免费少妇av软件| 一级毛片精品| 午夜福利在线观看吧| 欧美日本亚洲视频在线播放| 国产成人av激情在线播放| 黑丝袜美女国产一区| 色婷婷久久久亚洲欧美| 色精品久久人妻99蜜桃| 国产视频一区二区在线看| 午夜免费激情av| 亚洲精品在线美女| 中文欧美无线码| 午夜两性在线视频| 在线观看免费日韩欧美大片| 亚洲成人免费电影在线观看| 亚洲 国产 在线| 一进一出抽搐动态| 18禁观看日本| 亚洲免费av在线视频| 国产xxxxx性猛交| 美女扒开内裤让男人捅视频| 久久国产乱子伦精品免费另类| 欧美日韩一级在线毛片| 不卡一级毛片| 国产精华一区二区三区| 亚洲精品成人av观看孕妇| 十八禁人妻一区二区| 在线观看免费视频日本深夜| 午夜免费激情av| 成人特级黄色片久久久久久久| 桃红色精品国产亚洲av| 亚洲成人免费电影在线观看| 男人的好看免费观看在线视频 | 精品卡一卡二卡四卡免费| 一个人免费在线观看的高清视频| 亚洲av片天天在线观看| 亚洲人成77777在线视频| 国产熟女xx| 狠狠狠狠99中文字幕| 一二三四在线观看免费中文在| 亚洲五月婷婷丁香| 日韩免费av在线播放| 亚洲一区高清亚洲精品| 婷婷六月久久综合丁香| 国产精品一区二区免费欧美| 欧美日韩亚洲综合一区二区三区_| 丝袜在线中文字幕| www.熟女人妻精品国产| 男女之事视频高清在线观看| 欧美+亚洲+日韩+国产| 久久久久久久久免费视频了| 夜夜爽天天搞| 无遮挡黄片免费观看|