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

干旱脅迫對3種冷季型草坪草光譜反射率及生理特征的影響

2016-02-05 03:01:26趙志軍單貴蓮段新慧初曉輝
草原與草坪 2016年6期

趙志軍,單貴蓮,段新慧,姜 華,任 健,陳 功,初曉輝

(1.云南省草地動物科學研究院,云南 昆明 650212; 2.云南農業大學 草業科學系,云南 昆明 650201)

干旱脅迫對3種冷季型草坪草光譜反射率及生理特征的影響

趙志軍1,單貴蓮2,段新慧2,姜 華2,任 健2,陳 功2,初曉輝2

(1.云南省草地動物科學研究院,云南 昆明 650212; 2.云南農業大學 草業科學系,云南 昆明 650201)

為探討高光譜分析技術在草坪草水分監測上的應用,以草地早熟禾、高羊茅、多年生黑麥草為試驗材料,測定了干旱脅迫下3種冷季型草坪草生理指標及光譜反射率。結果表明:隨著干旱脅迫程度的加深,供試草坪草葉片相對含水量和葉綠素含量顯著下降(P<0.05),丙二醛含量、脯氨酸含量及敏感波段的(550和760 nm)的光譜反射率顯著上升(P<0.05),3種草坪草抗旱性能的強弱順序為:高羊茅>多年生黑麥草>草地早熟禾。3種草坪草敏感波段光譜反射率(R550 nm,R760 nm)及高光譜參數(Rg,R0,GNDVI)與生理指標間呈顯著或極顯著相關,以R550 nm與生理指標間的相關系數絕對值最大。同時R550 nm構建了供試草坪草生理指標的估測模型,經檢驗,模型估測值與實測值相關性達顯著或極顯著水平。

冷季型草坪草;干旱脅迫;生理特征;高光譜反射率;高光譜監測模型

高光譜分析技術是當今遙感領域的前沿技術,具有分析速度快、精度高、結果穩定、對目標物無破壞等優點[1-2]。同時,因高光譜的波段范圍很寬,如美國ASD Field spec FR 2500光譜儀可達到350~2500 nm[3],而肉眼的視覺僅為380~780 nm,因此,通過高光譜技術可以看到肉眼看不到的變化。居于高光譜遙感技術的優點,近年來,此技術在農業許多領域(如作物營養狀況的定量分析、產量估測以及草坪色澤、外觀質量及營養狀況評價等領域)得到了廣泛的應用。

在草坪質量、色澤及營養狀況方面,前人以草地早熟禾(Poapretensis)、高羊茅(Festucaarundinacea)、多年生黑麥草(Loliumperennel)、馬蹄金(Dichondrarepens)等為材料,開展了利用高光譜技術評價草坪質量、色澤及營養狀況方面的研究[4-9],結果表明,應用高光譜遙感技術監測草坪的色澤、外觀質量及營養狀況是確實可行的。但高光譜分析技術是否能準確監測草坪草的水分狀況及抗旱生理指標,目前鮮見報道。因此,以草地早熟禾、高羊茅、多年生黑麥草為試驗材料,開展干旱脅迫下3種冷季型草坪草生理生化指標及光譜反射特征變化的研究,分析草坪草抗旱生理生化指標與光譜反射率間的相關性,篩選出與草坪草抗旱指標密切相關的波段,在此基礎上建立草坪草抗旱生理指標定量評價的高光譜模型。

1 材料和方法

1.1 供試材料

供試材料為3種常用的冷季型草坪草種,分別為高羊茅品種美洲虎(Festucaarundinaceacv.Jaquar)、多年生黑麥草品種明星(Loliumperennelcv.Star)和草地早熟禾品種肯塔基(Poapretensiscv.Kentucky)。試驗草種均由克勞沃集團提供。

1.2 試驗設計

試驗于遮雨溫室進行。2012年6月20日,將采集于云南農業大學草坪基地的土壤(0~5 cm表層土,田間持水量約65%)裝入直徑30.5 cm,高20 cm的塑料花盆,將供試材料均勻播種于花盆,覆土1 cm,正常管理。2012年8月10日進行水分脅迫試驗,供試材料分別隔5 d(T5),10 d(T10),15 d(T15)澆水1次,同時設澆水1次/d的供試材料為對照(T1),澆水均勻噴灌至飽和(底部有水珠滲出)。每處理20次重復。9月10日處理結束,開展草坪草抗旱生理指標及光譜反射率的測定。

1.3 測定指標與方法

1.3.1 抗旱生理指標的測定 處理結束后,取新鮮葉片測定葉片相對含水量、總葉綠素含量、丙二醛含量和脯氨酸含量。其中葉片相對含水量采用稱重法測定,葉綠素含量采用浸泡-比色法測定,丙二醛含量采用硫代巴比妥酸法測定,脯氨酸含量采用酸性茚三酮顯色法測定[10],3次重復。

1.3.2 反射光譜數據的測定 采用美國Ocean公司的HR2000光譜儀。在無風、晴朗光線充足天氣(2012年9月10日),于上午10∶00~14∶00進行光譜反射率的測定[11]。測定時,利用一個鐵架固定探頭,移動各處理的花盆,探頭垂直向下距離草坪5 cm采集光譜數據。測量期間及時利用標準白板進行校正。每處理40次重復(模型的建立和檢驗各用20組數據)。

1.4 數據分析

選取與草坪草抗旱生理指標密切相關的高光譜參數(表1)和敏感波段反射率,分析供試草坪草抗旱生理指標與高光譜參數和敏感波段反射率的相關性,并通過回歸分析建立模型,以確定抗旱生理指標與高光譜參數和敏感波段反射率之間的定量關系。所有試驗數據采用SPSS 19.0進行方差分析、顯著性檢驗、相關和回歸分析。采用Excel2003作圖。

表1 幾種高光譜參數

Table 1 Several hyperspectral parameters

高光譜參數縮寫描述(計算公式)綠峰反射率Rg波長510~560nm最大的波段反射率紅谷反射率Ro波長640~680nm最小的波段反射率綠一歸一化差GNDVI(R750?R550)/(R750+R550)

注:R750 nm、R550 nm分別表示750 nm、550 nm處反射率

2 結果與分析

2.1 3種草坪草反射光譜特征比較研究

正常管理水平下3種草坪草反射光譜特征為,在可見光藍光區(400~520 nm),3種草坪植物的光譜反射率較低,該波段的反射曲線具有平滑的形狀和較低的數值,且3種草坪植物間的反射率差異不大(圖1)。在可見光綠光區(520~600 nm),3種草坪植物的反射光譜曲線具有波峰的形態和中等的反射率數值,約在550 nm處出現綠光波段葉綠素反射峰(“綠峰”),原因是可見光綠光區(520~600 nm)是一個低光合效應的光譜帶,該波段植物對光的吸收較少,因此光譜反射率相對較高,且該波段光譜反射率的高低主要受植物體內葉綠素含量高低的影響[12]。顏色較深的高羊茅草坪在綠光區反射率相對較低,而顏色較淺的草地早熟禾草坪在綠光區的反射率相對較高,在“綠峰”處反射率差異尤其明顯,高羊茅為6.89%,黑麥草為7.74%,早熟禾為9.68%。在可見光紅光區(630~690 nm),3種草坪植物的反射光譜曲線具有波谷的形態和較低的反射率數值,在680 nm處出現紅光波段葉綠素吸收谷(“紅谷”),原因是可見光紅光區(630~690 nm)是葉綠素最強吸收的光譜帶,是紅光區域中具有最強光合活性的光譜帶,該波段光譜反射率主要受植物葉綠素含量、植物的種類和覆蓋度的影響[11]。3種草坪植物的形態、覆蓋度等較為一致,因此,可推斷該波段3種草坪植物光譜反射率的差異主要是由于顏色的差異引起,顏色較深的高羊茅和黑麥草在該波段的反射率相對較低,而顏色較淺的草地早熟禾在該波段的反射率相對較高,在“紅谷”處高羊茅的反射率為2.58%,黑麥草為2.87%,早熟禾為4.27%。在近紅外光區(760~900 nm),3種草坪植物的反射光譜曲線形成一個高反射的平臺,此波段反射率隨著波長的增加而急劇升高,以760 nm處基本增加至最大值。多年生黑麥草、高羊茅和草地早熟禾的反射率在760 nm處差異也較為顯著,760 nm處高羊茅的反射率為46.46%,黑麥草的為51.28%,早熟禾的為69.54%。

圖1 正常管理水平下3種草坪草的光譜反射率Fig.1 Spectrl reflectance of three turfgrass under normal management

2.2 干旱脅迫下3種草坪草反射光譜及生理指標動態變化

“過去五年,臺州港抓改革、重創新,促轉型、惠民生,其中港口吞吐量節節攀升,港口投資創歷史最好紀錄,核心港區初步形成,積極融入全省港口一體化的成效逐漸顯現。”臺州市港航管理局局長周建業表示,下階段,臺州將進一步整合各個港區資源,充分發揮自身優勢,更加積極地融入長江經濟帶,并通過海鐵聯運向浙西腹地輻射,為浙江西部打通大宗商品物流通道。

2.2.1 干旱脅迫下3種草坪草反射光譜特征動態變化 干旱脅迫下3種草坪草在550和760 nm處的光譜反射率差異相對較大,即3種草坪草對干旱脅迫的敏感波段為550 nm和760 nm(圖2~4)。

圖2 干旱脅迫下高羊茅的光譜反射率Fig.2 Reflectance of Tall fascue under drought stress

圖3 干旱脅迫下黑麥草的光譜反射率Fig.3 Reflectance of Perennial ryegrass under drought stress

圖4 干旱脅迫下草地早熟禾的光譜反射率Fig.4 Reflectance of Bluegrassunder drought stress

試驗進一步以不同脅迫程度3種草坪草在550和760 nm處的光譜反射率進行方差分析并作圖(圖5,6),與T1,T5相比,T10和T15處理下3種草坪草在550和760 nm處的反射率顯著增加(P<0.05)。

圖5 干旱脅迫下3種草坪草在550nm處的光譜反射率Fig.5 R550 of three turfgrass under drought stress

圖6 干旱脅迫下3種草坪草在760nm處的光譜反射率Fig.6 R760 of three turfgrass under drought stress

2.2.2 干旱脅迫下3種草坪草生理指標的變化 在干旱或高溫等逆境脅迫下,各種植物葉片相對含水量均下降,脅迫程度越深,植物葉片相對含水量越低。對于不同品種來說,處理條件一致時葉片相對含水量減少速度小者抗旱性應較強[13-14]。干旱脅迫下3種冷季型草坪草的葉片相對含水量均下降。與正常灌溉(T1)相比,輕度水分脅迫下(T5),3種冷季型草坪草的葉片相對含水量下降不明顯,中度(T10)與重度(T15)脅迫下,3種冷季型草坪草的葉片相對含水量顯著下降,以草地早熟禾葉片相對含水量下降較為顯著(圖7)。

圖7 干旱脅迫對3種冷季型草坪草葉片相對含水量Fig.7 Effect of drought stress on leaf ralative water content of three cool-season turfgrass

干旱脅迫下3種冷季型草坪草的葉綠素含量均顯著下降(P<0.05),同一脅迫程度下,3種草坪草葉綠素含量的高低順序為:高羊茅>多年生黑麥草>草地早熟禾(圖8)。

圖8 干旱脅迫對3種冷季型草坪草葉片葉綠素的含量Fig.8 Effect of drought stress on chlorophyll content of three cool-season turfgrass

植物器官在逆境條件下,往往發生膜脂過氧化作用,丙二醛(MDA)是其產物之一,通常將其作為脂質過氧化指標,用于表示細胞膜脂過氧化程度和植物對逆境條件反映的強弱[15-16]。干旱脅迫下3種冷季型草坪草MDA含量均呈上升的變化趨勢,且隨脅迫程度的加深,MDA含量增加幅度變大,若以中度和重度脅迫下各草種MDA含量作為衡量標準,則3種冷季型草坪草抗旱性高低順序為:高羊茅>多年生黑麥草>草地早熟禾(圖9)。

圖9 干旱脅迫下3種冷季型草坪草丙二醛的含量Fig.9 Effect of drought stress on MDA content of three cool-season turfgrass

圖10 干旱脅迫下3種冷季型草坪草葉片脯氨酸的含量Fig.10 Effect of drought stress on pro content of three cool-season turfgrass

2.3 敏感波段光譜反射率及高光譜參數與生理指標相關性分析

將篩選出敏感波段光譜反射率(R550 nm,R760 nm)及高光譜參數(Rg,R0,GNDVI)與生理生化指標進行相關分析(表2),結果表明:3種草坪草敏感波段的光譜反射率(R550 nm,R760 nm)、綠峰反射率(Rg)和紅谷反射率(Ro)與葉片相對含水量和葉綠素含量呈極顯著負相關關系(P<0.01),與MDA含量呈極顯著正相關關系(P<0.01),與Pro含量呈顯著正相關關系(P<0.05)。以R550 nm與抗旱生理指標的相關系數絕對值最大。3種草坪植物的GNDVI與葉片相對含水量呈顯著正相關關系(R2=0.689,P<0.05),與葉綠素含量呈極顯著正相關關系(R2=0.795,P<0.01),與MDA含量呈極顯著負相關關系(R2=-0.724,P<0.01),與Pro含量相關性不顯著。

2.4 生理指標高光譜監測模型的構建與檢驗

相關分析結果顯示,3種草坪草在550 nm處的光譜反射率與抗旱生理指標的相關性最為顯著,因此,以550 nm處的光譜反射率(R550 nm)與葉片相對含水量、葉綠素含量、MDA含量、Pro含量進行回歸分析,得出供試草坪草抗旱生理指標的高光譜監測模型(表3)。

表2 敏感波段光譜反射率及高光譜參數與生理指標的相關分析

Table 2 Correlation analysis between reflectance/high spectral parameters and physiological index

綠峰反射率Rg紅谷反射率R0GNDVIR550nmR760nm葉片含水量-0.908??-0.828??0.698?-0.908??-0.761??葉綠素含量-0.942??-0.929??0.795??-0.943??-0.901??丙二醛含量0.936??0.864??-0.724??0.938??0.811??脯氨酸含量0.748??0.667?-0.5270.747??0.586?

注:*和**分別表示顯著相關(P<0.05)和極顯著相關(P<0.01)

表3 草坪草生理指標的高光譜監測模型

Table 3 Estimation model of physiological index based on R550

抗旱指標模型方程擬合系數R2葉片含水量y=-4.1132x+126.110.825???葉綠素含量y=8.0359x-0.52510.917???丙二醛含量y=0.2569x2-0.323x+0.69410.890???脯氨酸含量y=0.0149x2.23510.640??

為了檢驗監測模型的準確性和可靠性,將校驗樣本550 nm的光譜反射率分別代入上述回歸模型,得到供試草坪草生理指標的估測值,再將模型估測值和實測值進行相關性分析。結果表明:供試草坪草葉片相對含水量、葉綠素含量、丙二醛含量、脯氨酸含量的估測值與實測值相關性達顯著水平(P<0.05)(圖11~14),表明以R550 nm建立的供試草坪草抗旱生理指標監測模型是準確可靠的。

圖11 葉片相對含水量的預測值和實測值的相關性(n=12)Fig.11 Correlation between estimated and measured value of RWC

圖12 葉綠素含量的預測值和實測值的相關性(n=12)Fig.12 Correlation between estimated and measured value of Chlorophyll contene

圖13 MDA含量的預測值和實測值的相關性(n=12)Fig.13 Correlation between estimated and measured value of MDA

圖14 Pro含量的預測值和實測值的相關性(n=12)Fig.14 Correlation between estimated and measured value of Pro

3 討論

草坪草的抗旱性是草坪草通過一定的抗旱方式在干旱脅迫下的生存能力,而草坪草的抗旱能力與草坪草葉片相對含水量,葉綠素含量,丙二醛和脯氨酸含量等指標有著緊密的聯系。結果表明,隨著干旱脅迫程度的加深,供試草坪草葉片相對含水量和葉綠素含量顯著下降,丙二醛和脯氨酸含量顯著上升(P<0.05),在較低的細胞水勢下維持一定程度的代謝水平,其表現的是耐旱機制[19]。表明在干旱脅迫下3種冷季型草坪草通過Pro等滲透調節物含量的增加來加強其對干旱的適應和抵抗能力,這與其他草坪草在干旱脅迫下的生理變化是一致的[20]。

有研究指出,草坪是覆蓋度較高,低矮、致密、整齊、均一的人工植被,土壤背景、枯葉等因素對其影響較小,研究高光譜遙感技術在草坪上的應用比天然植被、農作物具有更高的精確度[5]。并且,高光譜技術操作簡單、結果可靠,與傳統測定植物生理指標的方法相比,可節省人力、物力,避免人為因素的影響,并且可以用在大面積草坪水分狀況監測及定量灌溉上[9]。研究建立了草坪草抗旱生理指標的高光譜監測模型,經檢驗,模型具有較高的準確性,證明利用高光譜分析技術監測草坪草的水分狀況確實可行,與前人研究結論一致[5,9]。

4 結論

隨著干旱脅迫程度的增加,3種冷季型草坪草葉片相對含水量和葉綠素含量顯著下降,丙二醛、脯氨酸含量和敏感波段光譜反射率(R550 nm,R760 nm)顯著上升(P<0.05)。3種草坪草抗旱性能由強到弱的順序為:高羊茅>多年生黑麥草>草地早熟禾。3種草坪草敏感波段光譜反射率(R550 nm,R760 nm)及高光譜參數(Rg,R0,GNDVI)與抗旱生理指標間呈顯著或極顯著相關關系,以R550 nm與抗旱生理指標的相關系數絕對值最大。用以R550 nm構建的監測模型來估測供試草坪的葉片相對含水量、葉綠素含量、丙二醛含量和脯氨酸含量,其估測值與實測值具有較好相關性。

[1] 張良培,張立富.高光譜遙感[M].武漢:武漢大學出版社,2005.

[2] 王海建,洪添勝,代芬.基于高光譜圖像技術的沙梨無損檢測[M].廣州:廣東農業科學出版社,2013.

[3] 騰安國,高峰,夏新成,等.高光譜技術在農業中的應用研究進展[M].南京:南京大學出版社,2009.

[4] 張文,陳功,龍光強,等.利用高光譜參數評價單播草坪色澤的研究[J].草原與草坪,2007(5):6-10.

[5] 張文,張建利,陳功.以高光譜植被指數研究草坪色澤[J].草地學報,2008,16(5):530-535.

[6] 龍光強,陳功,郭濱.干旱脅迫下馬蹄金草坪反射光譜研究[J].云南農業大學學報,2008,23(4):468-473.

[7] 錢育蓉,楊峰,李建龍.利用高光譜數據快速估算高羊茅牧草光合色素的研究[J].草業學報,2009,18(4):94-102.

[8] 楊峰,錢育蓉,李建龍.利用高光譜技術無損傷探測高羊茅草坪草營養狀況研究[J].中國草地學報,2009,31(1):86-91.

[9] 許岳飛,金晶煒,孫小玲.基于反射光譜技術評價草坪質量模型的研究[J].草業學報,2009,18(4):256-259.

[10] 李合生,孫群,趙世杰,等.植物生理生化實驗原理和技術[M].北京:高等教育出版社,1999.

[11] 張達,鄭玉權.高光譜遙感的發展與應用[M].上海:上海出版社,2013.

[12] 張文,張建利,陳功.葉面施肥對草坪高光譜反射率及色澤的影響[J].草業學報,2009,18(2):115-121.

[13] 初曉輝,岳信龍,任健,等.水分脅迫下4種冷季型草坪抗旱性比較研究[J].草原與草坪,2012,32(4):15-19.

[14] 王太亮,韓兆勝, 龐妙甜,等.NaCl脅迫對狗牙根生理生化特性的影響[J].草原與草坪,2015,35(1):63-67.

[15] 劉佳,郁繼華,徐秉良,等.干旱氣候條件下水分脅迫對辣椒葉片生理特性的影響[J].核農學報,2012, 26(8):1197-1203.

[16] 陳偉,張苗苗,宋陽陽,等.鎘脅迫對4種草坪草抗逆生理特性的影響[J].草原與草坪,2015,35(2):65-69.

[17] 吳艷華,陳雅君,夏忠強.草坪草抗旱生理機制的研究進展[M].哈爾濱:哈爾濱出版社,2005.

[18] 陳偉,張苗苗,宋陽陽,等.干旱脅迫對伊犁絹蒿生理生化指標的影響[J].草原與草坪,2015,34(6):82-85.

[19] 李吉躍.植物耐旱性及其機理[J].北京林業大學學報,1991,13(3):92-99.

[20] 李顯利,米福貴,閆立軍.草地早熟禾不同品種抗旱性的評價分析[J].草原與草坪,2010,30(3):43-46.

Study on spectral reflectance and physiological characteristics of three cool-season turfgrass under drought stress

ZHAO Zhi-jun1,SHAN Gui-lian2,DUAN Xin-hui2,JIANG Hua2,REN Jian2,CHEN Gong2,CHU Xiao-hui2

(1.AcademeofGrasslandandAnimalScience,Kunming650212,China; 2.DepartmentofGrasslandscience,YunnanAgriculturalUniversity,Kunming650201,China)

To test the application of spectral analysis technology for water monitoring of turfgrass,the spectral reflectance and physiological indices of three turf grass under drought stress were measured,and the correlationship between spectral reflectance and physiological indices was analyzed.The results showed the leaf relative water and chlorophyll content decreased,and the MDA,proline content,spectral reflectance at the wave of 550 nm and 760 nm of three cool-season turfgrass increased significantly with the deepening of drought stress (P<0.05),and the sequence of drought tolerance of three cool-season turfgrass was tall fescue>perennial ryegrass>bluegrass.The correlationship between spectral reflectance (R550,R760),hyperspectral parameters (Rg,R0,GNDVI) and physiological indices was significant (P<0.05 or 0.01),and the coefficients were highest between spectral reflectance R550and physiological indices.The monitoring estimation models of physiological indices were constructed based on spectral reflectance (R550).Upon testing,the correlationship between estimated and measured value was significant.

cool season turfgrass;reflectance;spectral monitoring model

2016-05-23;

2016-06-12

國家自然科學基金(31501999);云南省自然科學基金面上項目(2016FB047,2010CD061);云南省教育廳科學研究基金重點項目(2011Z033、2015Z1050);云南省教育廳科學研究基金項目(2011Y0477、2015Y204)資助

趙志軍(1979-),男,云南大理人,助理研究員,主要從事草坪草抗旱機理方面的研究。 E-mail:523241490@qq.com 初曉輝為通訊作者。

S 688.4

A

1009-5500(2016)06-0023-07

主站蜘蛛池模板: 手机成人午夜在线视频| 亚洲成综合人影院在院播放| 一区二区理伦视频| 操操操综合网| a亚洲视频| 激情爆乳一区二区| 91精品视频在线播放| 青青青视频91在线 | av午夜福利一片免费看| av免费在线观看美女叉开腿| 国产 在线视频无码| 久久中文字幕2021精品| 97视频在线观看免费视频| 青草91视频免费观看| 依依成人精品无v国产| 亚洲一级毛片| 成人午夜免费观看| 538精品在线观看| 亚洲第一成年网| 国产精品手机在线播放| 国产成人久久777777| 国产va视频| 91福利在线观看视频| 高清无码手机在线观看 | 亚洲熟女中文字幕男人总站| 视频二区亚洲精品| 黄色在线不卡| a级毛片一区二区免费视频| 亚洲欧美日韩高清综合678| 美女视频黄又黄又免费高清| 色窝窝免费一区二区三区| 日本不卡视频在线| 欧美一级高清免费a| 国产精品久久久久久久久久98| 国产欧美网站| 一级成人a毛片免费播放| 国产精品九九视频| 久久久久九九精品影院| 无码AV高清毛片中国一级毛片| 午夜三级在线| 四虎永久在线精品影院| 狠狠做深爱婷婷久久一区| 国产精品粉嫩| 强乱中文字幕在线播放不卡| 成色7777精品在线| 无码国产偷倩在线播放老年人| 国产美女免费| 午夜日b视频| 午夜视频日本| 国产Av无码精品色午夜| 99热线精品大全在线观看| 狠狠色丁婷婷综合久久| 无码网站免费观看| 亚洲av综合网| 美女一级毛片无遮挡内谢| 成人在线第一页| 伊人久久久久久久| 精品第一国产综合精品Aⅴ| 她的性爱视频| 国产精品hd在线播放| 人与鲁专区| 秋霞午夜国产精品成人片| 国产欧美日韩18| 欧美精品导航| 欧美性猛交xxxx乱大交极品| 国产凹凸视频在线观看| 高清亚洲欧美在线看| 美女毛片在线| 国产成人一区| 999国产精品永久免费视频精品久久 | 看国产毛片| 国产成人精品亚洲日本对白优播| 91亚洲免费视频| 久久永久免费人妻精品| 亚洲天堂.com| 亚洲天堂日韩av电影| 国产男女免费视频| 色综合久久久久8天国| 亚洲欧洲自拍拍偷午夜色| 国产成人精品午夜视频'| 另类专区亚洲| 在线精品亚洲一区二区古装|