999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

番茄、玉米套種膜下滴灌條件下農田地溫變化特征

2015-03-10 08:44:52龔雪文李仙岳史海濱田德龍彭遵原
生態學報 2015年2期

龔雪文,李仙岳,*,史海濱,田德龍,李 禎,彭遵原

1. 內蒙古農業大學水利與土木建筑工程學院, 呼和浩特 010018 2. 水利部牧區水利科學研究所, 呼和浩特 010020

番茄、玉米套種膜下滴灌條件下農田地溫變化特征

龔雪文1,李仙岳1,*,史海濱1,田德龍2,李 禎1,彭遵原1

1. 內蒙古農業大學水利與土木建筑工程學院, 呼和浩特 010018 2. 水利部牧區水利科學研究所, 呼和浩特 010020

土壤水熱狀況是影響作物生長的重要因素,以膜下滴灌及套種雙重效應為背景,針對番茄、玉米套種農田設置了高(番茄:180 mm,玉米:270 mm)、中(番茄:132 mm,玉米:202.5 mm)、低(番茄:84 mm,135 mm)3種不同的灌溉定額,采用定位監測法測定番茄帶膜間裸地(A)、覆膜番茄行間(B)、番茄帶與玉米帶膜間裸地(C)、覆膜玉米行間(D)共4個位置的不同深度主根區含水率和地溫動態變化,研究了地膜覆蓋、土壤水分以及作物遮陰對套種農田不同行間,不同深度土壤溫度的影響特征。結果表明:全生育期內,4個水平位置的地溫大小順序與土壤深度有關,10 cm以上的地溫大小順序為C>D>A>B,10 cm以下的地溫大小順序為D>B>A>C;從土壤溫度場的分布情況來看,地溫的最大變幅出現在C處,為8.10 ℃;最小變幅出現在B處,為4.71 ℃;含水率與地溫的關系表現為:5 cm以上地溫與含水率呈反比關系,20 cm以下地溫與含水率呈正比關系??梢姡追N覆膜滴灌農田合適的土壤水分狀況將有利于根區良好水熱的形成,從而更有利于作物的生長。研究結果將對西北地區滴灌農田管理有一定的參考價值。

套種; 膜下滴灌; 地溫; 含水率

套種不僅能提高單位面積的產量,而且能極大提高光、水、熱等資源的利用效率[1- 2],據報道,在內蒙古河套灌區套種已超過總種植面積的70%,然而由于套種農田對太陽輻射的有效攔截和利用不同[3],導致套種農田不同行間位置的地溫分布特征不同,而且覆膜后,塑料薄膜能夠有效將太陽能轉化為熱能,有效調節作物根區環境[4- 6]。因此將膜下滴灌技術與套種農田結合起來探討作物耕層溫度變化規律,對提高套種農田水熱利用效率具有重要的現實意義。近些年針對套種農田采用膜下滴灌技術的研究主要集中在作物系數、水分利用效率、水分遷移過程以及產量品質的影響[7- 10]等方面,而針對水、熱特性的研究主要以套種效應、覆膜效應、覆蓋效應、土壤水分效應等單項研究為主,如套種可有效利用太陽輻射提高產量[3];覆膜可有效提高土壤耕層溫度,促進作物生長[11];覆蓋可抑制地溫的大幅度波動[12];土壤含水量與土壤溫度呈顯著負相關[13]。另外,國內外學者還對覆膜條件或套種條件下的土壤熱效應及地溫的極值變化理論和模擬等方面做了研究[14- 19]。顯然,這些研究主要是針對或覆膜或套種或滴灌等單項技術開展的,而對套種農田的不同行間位置地溫分布特征尚未做深入研究。故本文主要針對套種及膜下滴灌耦合條件下多作物農田不同行間位置的地溫變化規律展開研究,并分析了套種模式下覆膜效應及土壤水分效應對地溫影響的綜合效應,本研究對于套種農田和膜下滴灌耦合作用下水熱理論發展具有一定的意義,并為該技術的進一步推廣奠定了基礎。

1 材料與方法

1.1 區域概況

2012年在內蒙古河套灌區磴口試驗站(40.20 E,107.1 N)開展了番茄套種玉米的試驗,該試驗田地處干旱半荒漠草原地帶,夏季多風干旱,冬季少雨寒冷,屬中溫帶大陸性季風氣候,海拔2004 m,年均氣溫7.2 ℃,年均日照時數2190 h,年均降水量198.3 mm,年均蒸發量2460.3 mm,年均風速2.74 m/s,多年平均無霜132 d。試驗田土壤質地為夾砂粘土,1 m內土壤平均容重為1.50 g/cm3,凋萎系數為8%,田持為24.1%。

1.2 試驗設計

供試作物為玉米(品種為中地77)和番茄(品種為屯河48),采用“一膜一管兩行”的滴灌布設方式,采用4行番茄套2行玉米的種植模式。在當地灌水經驗的基礎上,試驗設置了高、中、低3種不同的灌水處理,本試驗考慮了番茄與玉米需水量不同的特點,在每個小區中分別對番茄和玉米設定了不同的灌水定額(表1),采用精度為0.001 m3的水表(NSD6)控制灌水量,每個處理均設置了3個重復。

每個處理均布設了4組精度為0.1 ℃的水銀直角地溫計(RM-004),作物種植方式及地溫計的布設方式如圖1所示。

表1 全生育期番茄套種玉米的灌溉制度Table 1 Irrigation scheduling for intercropping between tomato and maize during the whole growing season

圖1 番茄、玉米套種方式及地溫計與滴灌帶布置圖Fig.1 Latout of intercropping between tomato and maize and geothermometer and drip irrigation system

1.3 試驗方法

1.3.1 氣象資料測量

離試驗區300 m的空曠地中設有自動氣象站(HOBO),用于觀測降水、氣溫、濕度、太陽輻射、大氣壓等。

1.3.2 土壤含水率測量

采用TRIME-PICO TDR和土鉆相結合的方法,以TRIME為主,每隔3 d測量1次,定期用土鉆取土烘干的方法加以校正,TDR與地溫計位置相同(見圖1),采用定位監測法測量含水率,測量深度為0—100 cm。

1.3.3 地溫觀測時間與方法

在6—9月,利用直角地溫計(RM-004)對5、10、15、20、25 cm處的土壤溫度進行定點觀測,測量時間與土壤含水率一致,每隔3 d測量1次,每天的觀測時間為8:00—20:00,2 h測1次。

1.3.4 作物生長指標測量

用直尺測量玉米和番茄的株高;葉面積的測量:玉米葉面積=葉長×葉寬×0.75,番茄采用折算系數法計算[20],以上指標均在每個生育期內測3次,取平均值。

1.3.5 試驗觀測日期

玉米4月20日播種,9月10日收獲,番茄5月5日移苗,8月10日采摘第1茬,8月25日采摘結束。試驗觀測時間從6月開始到9月結束,試驗觀測期內氣溫變化見圖2。

1.4 數據分析

本試驗數據采用Spss 18.0和Excel 2003進行統計分析,采用Excel 2003和Surfer 8.0 軟件繪制圖形。

2 結果與分析

2.1 套種農田不同位置及不同深度地溫日變化規律

在套種農田不同行間的不同位置土壤溫度日變化亦有所不同。通過對比分析6—9月份地溫數據可以看出,高水分處理T1各行間位置的地溫變化過程(圖3)均是隨著土層深度的增加地溫日波動呈減小趨勢,其中距地表5 cm處地溫日波動最劇烈,距地表25 cm處地溫日波動最平緩,距地表5、10、15、20、25 cm的日最高與最低地溫差值分別為:12.72、8.78、6.09、4.44、2.96 ℃,存在極顯著性差異(F=12.74005,P=0.007301)。

圖2 試驗觀測期內的氣溫變化Fig.2 The temperature change during the experimental period

圖3 整個觀測期內T1的日平均地溫變化圖Fig.3 The daily variation law of soil temperature during the whole observation period of T1

套種農田不同行間位置的地溫受地膜覆蓋和作物遮陰的影響較大。透明塑料薄膜對土壤的保溫效應主要體現在土層15 cm以下,從表2中可以看出,B、D位置的地溫明顯高于A、C位置,隨著土層的加深,膜內與膜外平均地溫的差值變化呈拋物線形,且均差在20 cm處達到最大值,可見在套種農田中,地膜可隔絕土壤與外界的水分交換和顯熱交換,使土壤的熱通量在20 cm處達到最大。地面增溫的唯一來源是太陽輻射,套種農田不同行間位置地溫受作物遮陰的影響程度也不相同,這在5 cm以上土層表現的較為明顯,從圖3中可以看出,距地表5 cm處的地溫大小關系表現為C>D>A>B。這是由于C處位于高矮作物的行間,遮陰面積最小,日照時間最長,因此表層(5 cm處)的地溫也最高;而B位于低矮作物行間,遮陰面積較大,地表所受太陽直射時間最短,因而地溫也最低。可見,在立體種植農田中,高桿作物遮陰區的地溫要遠遠高于矮桿作物遮陰區的地溫;在15 cm以下,膜內的日平均地溫要高于膜外。

表2 15 cm以下膜內與膜外的日平均地溫Table 2 Average day soil temperature of covered and uncovered treatments below 15 cm

2.2 套種農田不同處理及不同位置地溫顯著性分析

6—9月份是番茄、玉米進行營養生長的重要階段,地溫的變化直接影響其根系競爭水分和養分的能力。本文通過該時期地溫數據進行平均,分析不同位置、不同深度的地溫變化規律。結果顯示:套種農田不同位置地溫的差異性受到地膜覆蓋、作物遮陰及土壤水分的綜合影響,從不同位置地溫平均值變化來看,T1、T2的番茄一側(A、B)平均地溫無差異,而玉米一側(C、D)的平均地溫存在極顯著差異(表3),T3的4個位置平均地溫均在0.05水平上顯著。說明作物遮陰和土壤水分是導致套種農田不同位置地溫出現差異的主要原因。從番茄一側的不同深度地溫變化來看,T1在25 cm范圍內均無差異,而T2和T3分別在15 cm處和15 cm以下出現了顯著差異性;從玉米一側的不同深度地溫變化來看,T1、T2、T3的地溫均在15 cm處出現了差異性。這是因為在覆膜滴灌條件下,不同土壤水分所形成的獨特濕潤體不同,導致土壤水分在10—15 cm之間出現了差異性,從而改變了土壤的熱容量,加之套種農田不同行間位置的遮陰率不同,熱量在土層中分配時,其量值有所改變,因而出現了以上差異[21]。綜上可知,套種農田在營養生長的共生期,作物遮陰和水分差異是影響地溫的主要因素,而地膜覆蓋的影響變小。

表3 不同處理及不同位置地溫顯著性分析/℃Table 3 Analyze the significance of soil temperature at different processing and positions

同列數值后不同大小寫字母表示差異達1%和5%顯著水平

2.3 套種農田不同土壤水分對地溫變化的影響

由于7月份套種農田兩種作物生長發育趨于穩定,本文為了分析不同高低作物套種不同位置地溫的差異,以及不同土壤水分對地溫的影響,在7月份對土壤水分、不同位置地溫進行了連續監測。結果顯示:

距地表5 cm處,土壤受外界環境的影響較大,高水分T1的平均體積含水率明顯高于T2和T3,其中T1、T2相差3.8%,T1、T3相差7.2%(圖4),T1、T2、T3之間的體積含水率呈顯著性差異(F=6.223587,P=0.020088)。而地溫則相反,高水分T1各位置的平均地溫明顯低于T2和T3,其中T1、T2的平均地溫相差0.81 ℃,T1、T3相差1.65 ℃(圖5),T1、T2、T3之間的地溫無顯著性差異(F=1.44877,P=0.282652)。3個處理的不同位置含水率均表現為B>A>D>C,地溫則為D>C>B>A的規律。不同處理的膜內位置B(或D),含水率從T1減小到T2和T3依次降低2.73%、6.04%(或2.71%、6.63%),地溫則依次上升0.95 ℃、2.02 ℃(或0.47 ℃、0.97 ℃);不同處理的膜外位置A(或C),含水率從T1減小到T2和T3依次降低3.16%、6.44%(或6.59%、9.71%),地溫則依次上升1.16 ℃、1.97 ℃(或0.64 ℃、1.65 ℃)。說明在距地表5 cm處,地溫隨含水率的降低而上升,其中膜外的上升比率是膜內的1.2倍,表明在土壤表層膜內與膜外的溫度差較小。

圖4 不同灌溉定額水平下土壤的平均含水率Fig.4 Comparison of average soil moisture content for different irrigation quota

20 cm處,土壤受外界環境的影響變小,而受水分的影響變大,高水分T1的平均含水率仍高于T2和T3,其中T1、T2相差3.61%,T1、T3相差6.9%(圖4),無顯著性差異(F=3.338497,P=0.082302)。地溫亦為T1高于T2和T3,其中T1、T2的平均地溫相差0.87 ℃,T1、T3相差1.48 ℃(圖5),T1、T2、T3之間的地溫無顯著性差異(F=1.981115,P=0.193657)。3個處理的不同位置體積含水率均表現為B>A>D>C,地溫則為D>B>C>A的規律。不同處理的膜內位置B(或D),含水率從T1減小到T2和T3依次降低2.91%、6.58%(或2.47%、5.77%),地溫則依次降低1.04 ℃、1.63 ℃(或1.03 ℃、1.52 ℃);不同處理的膜外位置A(或C),含水率從T1減小到T2和T3依次降低3.6%、6.82%(或5.47%、8.4%),地溫則依次降低1 ℃、0.9 ℃(或0.43 ℃、1.52 ℃)。說明在20 cm處,地溫隨含水率的降低而降低,其中膜外的降低比率是膜內的2.1倍,表明較深層(20 cm以下)膜內地膜波動幅度小于膜外地溫波動幅度。

2.4 套種農田不同處理及不同位置土壤溫度場的分布特征

為進一步研究套種農田不同行間位置之間地溫的傳遞效應,本文對作物生長旺期,兩種作物都發育充分條件下選擇7、8月份觀測數據進行平均后繪成套種農田地溫二維分布圖(圖6)。由于監測的是農田日地溫的平均值,從圖中可以看出,A、C膜外位置在土體內溫度變化幅度大,B、D膜內位置在土體內溫度變化幅度小,3個處理在C位置處在土壤上層地溫明顯高于其它位置,這主要是C位置兩側遮陰小,吸熱多的原故,而B、D位置盡管覆膜,但白天這兩個位置受兩側作物的遮陰影響,地溫反而沒有C位置大。在土壤下層,C位置的地溫明顯低于其他位置地溫,這是由于C位置無覆膜,保溫效果差,到了夜間,上層熱量難以向下傳遞。可見,在15 cm處C位置地溫低于B、D位置地溫,無覆膜情況下白天太陽輻射增溫效應僅對土壤上層有效,而覆膜情況下則0—25 cm都能升溫。另外對于不同水分處理溫度場的分布主要是土體內含水量不同,對于高水分處理的T1土體內平均含水量較大,因而比熱大,所以最高地溫與最低地溫的差值變小,3個處理的最高地溫與最低地溫的差值分別為6.02、7.26、8.04 ℃。

圖5 不同處理在5cm和20cm處的平均溫度日變化Fig.5 The average temperature of the diurnal variation (soil depth 5cm and 20cm)

圖6 不同處理土壤剖面的溫度分布特征Fig.6 Comparison of soil profile characteristics of temperature distribution

3 結論

3.1 番茄套種玉米能有效利用太陽輻射提高作物根區溫度,其不同行間位置的地溫變化主要受地膜覆蓋、作物遮陰和土壤水分含量的影響。在套種農田中,地膜的保溫作用主要體現在15 cm以下,而作物遮陰和土壤水分含量則對整個耕層都有影響??梢姡侠淼奶追N模式和適宜的土壤水分含量有利于作物根區形成良好的生長環境,從而促進作物的生長。

3.2 套種農田不同深度土層土壤含水率與地溫的關系不同。本試驗通過分析番茄套種玉米表層(5 cm)及深層(20 cm)的地溫與含水率的關系得出:5 cm以上的地溫與含水率呈反比關系,20 cm以下地溫與含水率呈正比關系,可見,在套種農田中,地溫與含水率在不同土層呈現的關系也不相同,而并不像単種農田一樣,兩者呈負相關關系。

3.3 從套種農田不同行間位置地溫的傳遞規律可以看出,高矮作物行間位置(C處)在白天有明顯的增溫效應,且溫度可有效傳遞到作物兩側,而覆膜可有效防止土壤溫度的散失,因此高矮作物套種應選擇恰當的行距,從而有效提高作物根區的溫度。

[1] Willey R W. Resource use in intercropping systems. Agricultural Water Management, 1990, 17(1): 215- 231.

[2] 王仰仁, 李明思, 康紹忠. 立體種植條件下作物需水規律研究. 水利學報, 2003, 7: 90- 95.

[3] Awal M A, Koshi H, Ikeda T. Radiation interception and use by maize/peanut intercrop canopy. Agricultural and Forest Meteorology, 2006, 139(1): 74- 83.

[4] Li F M, Guo A H, Wei H. Effects of clear plastic film mulch on yield of spring wheat. Field Crops Research, 1999, 63(1): 79- 86.

[5] Choudhary V K, Bhambri M C, Pandey N, Sharma H G. Effect of drip irrigation and mulches on physiological parameters, soil temperature, picking patterns and yield in capsicum (CapsicumannuumL.). Archives of Agronomy and Soil Science, 2012, 58(3): 277- 292.

[6] 潘渝, 郭謹, 李毅, 張江輝, 丁新利. 地膜覆蓋條件下的土壤增溫特性. 水土保持研究, 2002, 9(2): 130- 134.

[7] Gao Y, Duan A W, Sun J S, Li F S, Liu Z G, Liu H, Liu Z D. Crop coefficient and water-use efficiency of winter wheat/spring maize strip intercropping. Field Crops Research, 2009, 111(1): 65- 73.

[8] Sampathkumar T, Pandian B J, Mahimairaja S. Soil moisture distribution and root characters as influenced by deficit irrigation through drip system in cotton-maize cropping sequence. Agricultural Water Management, 2012, 103: 43- 53.

[9] Zhang L, van der Wert W, Zhang S, Li B, Spiertz J H J Spiertz. Growth, yield and quality of wheat and cotton in relay strip intercropping systems. Field Crops Research, 2007, 103(3): 178- 188.

[10] Makumba W, Janssen B, Oenema O, Akinnifesi F K, Mweta D, Kwesiga F. The long-term effects of a gliricidia-maize intercropping system in Southern Malawi, on gliricidia and maize yields, and soil properties. Agriculture, Ecosystems and Environment, 2007, 118(1): 237- 243.

[11] 張治, 田富強, 鐘瑞森, 胡和平. 新疆膜下滴灌棉田生育期地溫變化規律. 農業工程學報, 2011, 27(1): 44- 51.

[12] 孫建, 劉苗, 李立軍, 劉景輝, Acharya S N. 不同耕作方式對內蒙古旱作農田土壤水熱狀況的影響. 生態學報, 2010, 30(6): 1539- 1547.

[13] 李全起, 陳雨海, 于舜章, 吳巍, 周勛波, 董慶裕, 余松烈. 灌溉與秸稈覆蓋條件下冬小麥農田小氣候特征. 作物學報, 2006, 32(2): 306- 309

[14] Mahrer Y, Naot O, Rawitz E, Katan J. Temperature and moisture regimes in soils mulched with transparent polyethylene. Soil Science Society of America Journal, 1984, 48(2): 362- 367.

[15] 王建東, 龔時宏, 許迪, 馬曉鵬, 于穎多. 地表滴灌條件下水熱耦合遷移數值模擬與驗證. 農業工程學報, 2010, 26(12): 66- 71.

[16] 吳叢林, 黃介生, 沈榮開. 地膜覆蓋條件下SPAC系統水熱耦合運移模型的研究. 水利學報, 2000, (11): 89- 96.

[17] Hunt H W, Fountain A G, Doran P T, Basagic H. A dynamic physical model for soil temperature and water in Taylor Valley, Antarctica. Antarctic Science, 2010, 22(4): 419- 434.

[18] 李毅, 王文焰, 門旗, 鐘新才, 謝香文. 寬地膜覆蓋條件下土壤溫度場特征. 農業工程學報, 2001, 17(3): 32- 36.

[19] 李毅, 紹明安. 新疆農田作物覆膜地溫極值的時空變化. 應用生態學報, 2004, 15(11): 2039- 2044.

[20] 劉浩, 孫景生, 段愛旺, 孫磊, 梁媛媛. 基于AutoCAD軟件確定番茄與青椒葉面積的簡易方法. 中國農學通報, 2009, 25(5): 287- 293.

[21] 李慧星, 夏自強, 馬廣慧. 含水量變化對土壤溫度和水分交換的影響研究. 河海大學學報(自然科學版), 2007, 35(2): 172- 175.

The interplanting between tomato and maize of soil temperature dynamics under mulched drip irrigation

GONG Xuewen1, LI Xianyue1,*, SHI Haibin1, TIAN Delong2, LI Zhen1, PENG Zunyuan1

1WaterConservancyandCivilEngineeringCollegeofInnerMongoliaAgriculturalUniversity,Huhhot010018,China2InstituteofWaterResourcesforPastoralArea,MinistryofWaterResources,Huhhot010020,China

It is very significant for maintaining the optimal soil moisture and soil temperature in the crop field, because that is the important foundation of developing an optimal irrigation system. The environment of reasonable soil moisture and temperature in the field will not only enhance the water use efficiency but also benefit crop growth and increase crop yield, especially for interplanting field with complex crop environments. In the experiments, the interplanting of tomato and corn under mulched drip irrigation were tested in Hetao Irrigation District in 2012. We set up three treatments which are the high irrigation quota (tomato: 180 mm, corn: 270 mm), medium irrigation quota (tomato: 132 mm, corn: 202.5 mm) and low irrigation quota (tomato: 84 mm, corn: 135 mm). The moisture content and soil temperature of different depths of root zone were measured in the middle of bare soil between tomatoes belt (A), the middle of mulched tomato lines (B), the middle of bare soil between tomatoes and corn belt (C), the middle of mulched corn lines (D) and crop growth and yield were observed in different growing seasons. The effects of plastic film, soil moisture and crop shade on soil temperature at different rows and depths of inter-planting were studied by the experiment. The results showed that it had apparent effect of mulch on heat preservation under 15 cm soil depths, but the soil temperatures of whole root zone were influenced by crop shade and soil moisture; The soil temperatures of four horizontal positions were significant variation for different the soil depths (5 cm, 10 cm, 15 cm, 20 cm, 25 cm), and the proper sequence of average soil temperature above 10 cm soil depth in different positions was: C>D>A>B, the one under 10 cm soil depth was: D>B>A>C; The study on distribution of the soil temperature in the interplanting field showed that the maximum variational amplitude of soil temperature appeared at the C position, and it would reach 8.10 ℃; the minimum variational amplitude of soil temperature appeared in the B position, and it was only 4.71 ℃. There was definitely some inverse correlation between the soil temperature and soil moisture above 5 cm soil depth, however direct proportion was observed below 20 cm soil depth. The study also found that there was strong temperature effect in the middle of lines between high and low crops (such as C position) during daytime, and the temperature can influence the crop of two sides, so mulched cultivation has obvious effects for restraining the loss of heat. It is thus clear that the optimal soil water and soil temperature are beneficial to form better soil moisture and temperature environment and high crop yield in the intercropping field under mulched drip irrigation. The study will be very meaningful for drip irrigation farm management in the Northwest Region.

intercropping; mulched drip irrigation; soil temperature; soil water

國家自然科學基金(51109105, 51349003); 國家科技支撐(2011BAD29B03); 內蒙古高等學?!扒嗄昕萍加⒉胖С钟媱潯? 博士后基金(2011M500547, 2012T50250)

2013- 03- 24;

日期:2014- 03- 25

10.5846/stxb201303240506

*通訊作者Corresponding author.E-mail: lixianyue80@126.com

龔雪文,李仙岳,史海濱,田德龍,李禎,彭遵原.番茄、玉米套種膜下滴灌條件下農田地溫變化特征.生態學報,2015,35(2):489- 496.

Gong X W, Li X Y, Shi H B, Tian D L, Li Z, Peng Z Y. The interplanting between tomato and maize of soil temperature dynamics under mulched drip irrigation.Acta Ecologica Sinica,2015,35(2):489- 496.

主站蜘蛛池模板: 亚洲男人天堂2020| 国产成人福利在线视老湿机| 国产在线欧美| 美女被操91视频| 久久无码av一区二区三区| 日韩午夜福利在线观看| av在线5g无码天天| 亚洲精品无码日韩国产不卡| 热re99久久精品国99热| 日本一本正道综合久久dvd| 99草精品视频| 2021国产乱人伦在线播放 | 情侣午夜国产在线一区无码| 91精品国产自产在线观看| 久久这里只有精品66| 野花国产精品入口| 又爽又大又光又色的午夜视频| 99精品国产自在现线观看| 国模视频一区二区| 免费在线成人网| 亚洲国模精品一区| 国产成本人片免费a∨短片| 国产女人在线| 国产精品人人做人人爽人人添| 一级黄色片网| 亚洲乱码精品久久久久..| 国产区网址| 久久久久九九精品影院| 亚洲 日韩 激情 无码 中出| 国产高清国内精品福利| 色视频久久| 欧美v在线| 免费一级毛片在线播放傲雪网| 青草视频网站在线观看| 国产精女同一区二区三区久| 欧美福利在线播放| 啪啪永久免费av| 亚洲精品爱草草视频在线| 亚洲日本韩在线观看| 免费a级毛片视频| 久久精品国产精品国产一区| 欧美中日韩在线| 国产人前露出系列视频| 手机在线看片不卡中文字幕| 国产一级毛片在线| 免费一级大毛片a一观看不卡 | 67194在线午夜亚洲| 久久综合丝袜长腿丝袜| 青草精品视频| 午夜性刺激在线观看免费| 丁香亚洲综合五月天婷婷| 色婷婷色丁香| 亚洲系列中文字幕一区二区| 色欲色欲久久综合网| 97国产在线视频| 女人毛片a级大学毛片免费| 精品无码一区二区三区电影| 久久久久国产精品免费免费不卡| 色综合天天综合中文网| 九色91在线视频| 精品三级在线| 国产精品一区二区不卡的视频| 精品国产乱码久久久久久一区二区| 高h视频在线| 在线国产91| 久久成人国产精品免费软件| 午夜a视频| 狠狠ⅴ日韩v欧美v天堂| 欧美精品亚洲日韩a| 精品视频福利| 久久精品免费看一| 亚洲天堂色色人体| 国产成人精品一区二区秒拍1o| 精品国产一区91在线| 国产在线观看精品| 国产精品视频猛进猛出| 精品无码人妻一区二区| 亚洲视频免费在线看| 中文字幕欧美成人免费| 色综合手机在线| 国产美女精品人人做人人爽| 国产人前露出系列视频|