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

    農(nóng)藝與氣候因素對青藏高原栽培大麥淀粉含量空間分布的影響

    2017-11-23 02:33:44鐘志明王建林馮西博侯維海王改花大次卓嘎
    生態(tài)學報 2017年21期
    關鍵詞:日較差大麥平均氣溫

    鐘志明,王建林,馮西博,付 剛,侯維海,王改花,大次卓嘎

    1 中國科學院地理科學與資源研究所,北京 100101 2 西藏農(nóng)牧學院, 林芝 860000

    農(nóng)藝與氣候因素對青藏高原栽培大麥淀粉含量空間分布的影響

    鐘志明1,王建林2,*,馮西博2,付 剛1,侯維海2,王改花2,大次卓嘎2

    1 中國科學院地理科學與資源研究所,北京 100101 2 西藏農(nóng)牧學院, 林芝 860000

    為了揭示不同環(huán)境因子對青藏高原栽培大麥籽粒淀粉含量(GSC)積累的影響程度,完善大麥GSC空間分異與環(huán)境因子的關系,明確青藏高原不同地區(qū)大麥品種GSC的環(huán)境效應,利用83個樣點的地理、氣候、土壤、農(nóng)藝因子數(shù)據(jù),研究了青藏高原栽培大麥GSC的分布特征,結(jié)果表明:(1)在地理水平方向上,青藏高原栽培大麥GSC的水平分布總體呈現(xiàn)出斑塊狀交錯分布和南高北低的格局,并形成了以西藏拉孜、隆子、堆龍德慶、曲水、尼木、定日、薩迦、達孜、扎囊、日喀則、墨竹工卡、貢嘎、瓊杰為中心的青藏高原中南部和以青海共和、貴德、海晏、同德為中心的青藏高原東北部等2個栽培大麥GSC高值區(qū);(2)在地理垂直方向上,栽培大麥GSC的變化呈現(xiàn)出“S”型分布格局,即在海拔3300.0—3600.0 m以上,隨著海拔的升高,栽培大麥GSC逐漸增加,在海拔4200.0 m與4500.0 m之間達到最高值,此后隨著海拔的升高略有下降;(3)影響栽培大麥GSC的因子從大到小的順序是穗密度﹥6月平均氣溫日較差>芒長﹥9月平均氣溫>1月平均氣溫>年日照時數(shù)﹥≥0℃積溫>5月平均氣溫﹥8月平均氣溫日較差>8月平均氣溫>6月平均氣溫>≥10℃積溫>6月平均月降水量>5月平均月降水量>7月平均相對濕度>8月平均相對濕度>7月平均氣溫。這一研究結(jié)果顯示,對栽培大麥GSC影響最大的是基因型,其次是氣候因素,土壤因素對GSC的影響不明顯。影響栽培大麥GSC的農(nóng)藝因子主要是穗密度和芒長,氣候因子主要是拔節(jié)抽穗期的氣溫日較差和籽粒灌漿成熟期的平均氣溫,日照和降水的影響相對較小。

    農(nóng)藝因素;氣候因素;大麥;淀粉含量;空間分布;青藏高原

    大麥是全球栽培的第四大谷類作物,其適應性強,利用價值廣泛,不僅廣泛用于食品和飼料,還用于啤酒釀造。淀粉是大麥中存儲的一種重要多糖,同時也是全球重要的食物和工業(yè)原料。淀粉是便宜的、可再生,能進行生物降解[1]。淀粉有直鏈淀粉和支鏈淀粉兩類,前者為無分支的螺旋結(jié)構(gòu),后者以24—30個葡萄糖殘基以α- 1,4-糖苷鍵首尾相連而成,在支鏈處為α- 1,6-糖苷鍵[2- 4]。不論是食用、飼用還是釀造啤酒,大麥籽粒淀粉含量(Grain Starch Content,GSC)均是一個十分重要的品質(zhì)指標[5]。近年來,一些研究者認為大麥GSC不僅明顯受其本身遺傳因素的影響,還與栽培條件及產(chǎn)地的環(huán)境條件密切相關[6- 14]。但是,不同的環(huán)境因子是如何影響大麥GSC的,目前尚存在很大爭議[15-22]。而且,這些研究多是在幾個地點進行試驗所得出的結(jié)論。截至目前,國內(nèi)外尚未見到在大尺度上進行栽培大麥GSC與環(huán)境因子關系的研究報道。

    青藏高原既是世界栽培大麥起源中心之一,又是中國重要的裸大麥(青稞)種植區(qū)。這里蘊藏著許多高GSC的種質(zhì)資源[5,23-24],急需進一步挖掘與鑒別利用。迄今為止,除個別學者測定分析過部分西藏栽培大麥GSC外,尚未見到有關青藏高原栽培大麥GSC與環(huán)境因子關系的系統(tǒng)性研究報道。為此,本文試圖以青藏高原為研究區(qū),從大尺度上研究栽培大麥GSC的空間分異規(guī)律,以期為揭示不同環(huán)境(地理、氣候、土壤)因子對青藏高原栽培大麥GSC積累的影響程度,完善大麥GSC空間分異與環(huán)境因子的關系,明確青藏高原不同地區(qū)大麥品種GSC的環(huán)境效應,為指導青藏高原大麥育種、栽培和加工提供理論依據(jù)。

    1 材料和方法

    1.1 研究區(qū)概況

    研究區(qū)涉及四川、甘肅、青海、西藏等4個省(區(qū))主要栽培大麥種植區(qū)。研究區(qū)跨越山地草原、高寒草甸草原、高寒灌叢草甸、山地針葉林、山地灌叢草原、山地荒漠、半荒漠等7個自然植被帶,年均日照時數(shù)1642—3393 h,年均溫-0.6—12.9℃,年均降水量150.0—890.0 mm。研究對象為當?shù)胤N植的栽培大麥(均為裸大麥)。

    1.2 樣品采集與分析

    (1)樣點布置:2014年6—9月,根據(jù)青藏高原大氣溫度和降水分布格局,于不同生態(tài)種植區(qū)(27°—38°N、79°—104°E)設置83個典型樣點。在樣點設置時,重點選擇長勢整齊、成片種植的栽培大麥田,所有樣點的具體地理位置如圖1所示。

    圖1 采樣點分布圖Fig.1 The locations of sampling sites

    (2)植物樣品采集與分析:采樣在栽培大麥成熟期進行。采樣前先用GPS定位、咨詢當?shù)剞r(nóng)民了解樣點所在農(nóng)田種植的栽培大麥品種、施肥情況等基本信息。采樣時在每個樣點設5—6個樣方,于每個樣方(面積不小于20 m2)中隨機采收栽培大麥籽粒2000 g,送回實驗室風干后,用酶水解法(GB5009.9- 85)測定栽培大麥的GSC。同時,在每個樣點,隨機選擇20個植株進行常規(guī)大田考種。

    (3)土壤樣品采集與分析:在每個樣點進行植物采樣的同時,設5—6個樣方,于每個樣方(面積不小于20.0 m2)中隨機采集耕作層(0—30.0 cm)土壤樣品3份,將其裝入土壤布袋,送回實驗室風干后,磨細過1 mm篩,測定土壤全氮、全磷、全鉀、速效氮、速效磷、速效鉀、pH值、有機質(zhì)等理化指標。其中:土壤全氮和速效氮采用半微量凱氏法,土壤全磷和速效磷測定采用鉬銻抗比色法,土壤全鉀和速效鉀測定采用原子吸收分光光度法,土壤有機質(zhì)采用重鉻酸鉀氧化—外加熱法、土壤pH值采用電位法進行測定。

    1.3 數(shù)據(jù)分析與處理

    首先將83個樣點栽培大麥GSC數(shù)據(jù)分成5類(即<40,40—45,45—50,50—60和>60,單位:%),然后利用ARCGIS 9.3軟件作其空間分布圖。同時,采用DPS軟件進行統(tǒng)計分析。運用單因素方差分析和最小顯著差異法進行數(shù)據(jù)組間的差異性比較,用逐步回歸分析法建立栽培大麥GSC與地理、氣候(注:氣候資料來源于中國國家氣象信息中心)、土壤、農(nóng)藝因子之間的關系方程,并用隨機森林回歸分析法綜合分析不同因子對青藏高原栽培大麥GSC的影響程度。

    2 結(jié)果與分析

    2.1 青藏高原栽培大麥GSC的分布特征

    2.1.1 青藏高原栽培大麥GSC的水平分布特征

    從圖2可以看出,栽培大麥GSC在青藏高原形成了兩個高值區(qū)。其中:一個是介于87.0°—92.0°E、29.0°—30.0°N之間,以西藏拉孜、隆子、堆龍德慶、曲水、尼木、定日、薩迦、達孜、扎囊、日喀則、墨竹工卡、貢嘎、瓊杰為中心的青藏高原中南部高值區(qū),這一區(qū)域栽培大麥GSC平均在(62.1635±5.6130)%左右;另一個是介于100.0°—101.5°E、35.0°—36.5°N之間,以青海共和、貴德、海晏、同德為中心的青藏高原東北部高值區(qū),這一區(qū)域栽培大麥GSC平均在(57.9762±2.0447)%左右。同時,從圖2還可以看出,青藏高原栽培大麥GSC的水平分布總體呈現(xiàn)出斑塊狀交錯分布和南高北低的趨勢。

    圖2 青藏高原栽培大麥GSC水平分布狀況Fig.2 GSC distribution conditions along the horizontal direction in Qinghai-Tibet Plateau

    2.1.2 青藏高原栽培大麥籽粒淀粉的垂直分布特征

    從表1可以看出,青藏高原栽培大麥GSC的平均值為(48.1675±9.9556)%,變異系數(shù)為19.69%。隨著海拔升高,栽培大麥GSC的變化呈現(xiàn)出“S”型分布格局。即在海拔3300.0—3600.0 m以下, 栽培大麥GSC隨著海拔的變化略有波動性變化,但是總體含量很低,平均值為(45.9538±9.2987)%。在海拔3300.0—3600.0 m以上,隨著海拔的升高,栽培大麥GSC逐漸增加,在海拔4200.0—4500.0 m之間達到最高值,此海拔區(qū)間栽培大麥籽粒淀粉平均含量為(50.4275±10.7836)%,此后隨著海拔的升高栽培大麥GSC略有下降。其中在3900.0—4500.0 m之間,栽培大麥GSC平均在(50.1431±10.7156)%左右。同時,從表1還可以看出,栽培大麥GSC變異系數(shù)隨著海拔的升高而呈現(xiàn)有波動性變化。

    表1 不同海拔梯度栽培大麥GSC的分布狀況

    2.2 青藏高原栽培大麥GSC與生態(tài)因子的關系

    2.2.1 栽培大麥GSC與地理因子的關系

    基于逐步回歸分析法,建立栽培大麥GSC與地理緯度(X1)、地理經(jīng)度(X2)、海拔(X3)等地理因子的關系方程如下:

    Y=26.7002+0.1311X1-0.1752X2+0.0040X3

    (N=83R=0.1646R1=0.0517R2=-0.0327R3=0.1264) (1)

    從式(1)可以看出,栽培大麥GSC(Y)與地理經(jīng)度(X1)、海拔高度(X3)呈正相關關系,而與地理緯度(X2)呈負相關關系。根據(jù)標準誤差檢驗,該回歸方程未通過=0.05水平的顯著性檢驗。表明,栽培大麥GSC隨著地理經(jīng)度和海拔高度的增加以及地理緯度的降低而增加。同時也表明,地理因子對栽培大麥GSC的影響未達到顯著性水平(P<0.05)。

    2.2.2 栽培大麥GSC與氣候因子的關系

    基于逐步回歸分析法,建立栽培大麥GSC與年均日照時數(shù)(X1)、年蒸發(fā)量(X2)、≥0℃積溫(X3)、≥10℃積溫(X4)、年均氣溫(X5)、年均相對濕度(X6)、年均降水量(X7)、1月平均氣溫(X8)、5月平均氣溫(X9)、6月平均氣溫(X10)、7月平均氣溫(X11)、8月平均氣溫(X12)、9月平均氣溫(X13)、5月平均氣溫日較差(X14)、6月平均氣溫日較差(X15)、7月平均氣溫日較差(X16)、8月平均氣溫日較差(X17)、9月平均氣溫日較差(X18)、5月平均相對濕度(X19)、6月平均相對濕度(X20)、7月平均相對濕度(X21)、8月平均相對濕度(X22)、9月平均相對濕度(X23)、5月平均降水量(X24)、6月平均月降水量(X25)、7月平均降水量(X26)、8月平均降水量(X27)、9月平均降水量(X28) 等氣候因子的關系方程如下:

    Y=68.4096+0.0066X1-7.3785X3-3.0721X4+1.2481X6+12.2109X8-8.0992X9+12.6286X10-10.3583X11+

    9.8997X12-9.5364X13-3.1067X15+3.2766X17-1.8666X21+0.3827X22-0.5219X24+0.8191X25

    (N=83R= 0.5968*R1=0.3449**R3=-0.3538**R4=-0.3576**R6=0.1716R8=0.3728**R9= -0.2635*R10=0.3262**R11=-0.4007**R12=0.293*R13=-0.2426*R15=-0.3747**R17=0.3898**R21=-0.374**R22=0.3652**R24=-0.245*R25=0.2890*)

    (2)

    從式(2)可以看出,影響栽培大麥GSC(Y)的主要氣候因子是年均日照時數(shù)(X1)、≥0℃積溫(X3)、≥10℃積溫(X4)、年均相對濕度(X6)、1月平均氣溫(X8)、5月平均氣溫(X9)、6月平均氣溫(X10)、7月平均氣溫(X11)、8月平均氣溫(X12)、9月平均氣溫(X13)、6月平均氣溫日較差(X15)、8月平均氣溫日較差(X17)、7月平均相對濕度(X21)、8月平均相對濕度(X22)、5月平均降水量(X24)、6月平均月降水量(X25),而年蒸發(fā)量(X2)、年均氣溫(X5)、年均降水量(X7)、5月平均氣溫日較差(X14)、7月平均氣溫日較差(X16)、9月平均氣溫日較差(X18)、5月平均相對濕度(X19)、6月平均相對濕度(X20)、9月平均相對濕度(X23)、7月平均降水量(X26)、8月平均降水量(X27)、9月平均降水量(X28)的影響則很小。

    同時,從式(2)還可以看出,栽培大麥GSC與年均日照時數(shù)(X1)、1月平均氣溫(X8)、6月平均氣溫(X10)、8月平均氣溫(X12)、8月平均氣溫日較差(X17)、8月平均相對濕度(X22)、6月平均月降水量(X25)的偏相關系數(shù)達到顯著或極顯著正相關差異水平,與≥0℃積溫(X3)、≥10℃積溫(X4)、5月平均氣溫(X9)、7月平均氣溫(X11)、9月平均氣溫(X13)、6月平均氣溫日較差(X15)、7月平均相對濕度(X21)、5月平均降水量(X24)的偏相關系數(shù)達到顯著或極顯著負相關差異水平,但是與年均相對濕度(X6)的偏相關關系未達到顯著性(P<0.05)差異。

    2.2.3 栽培大麥GSC與土壤因子的關系

    基于逐步回歸分析法,建立栽培大麥GSC與土壤有機質(zhì)(X1)、全氮(X2)、全磷(X3)、全鉀(X4)、速效氮(X5)、速效磷(X6)、速效鉀含量(X7)和土壤pH值(X8)等土壤因子的數(shù)學模型:

    Y=55.6769-7.1357X4-0.0260X7(N=83R=0.2523R4=-0.145R7=-0.1929)

    (3)

    從式(3)可以看出,影響栽培大麥GSC(Y)的主要土壤因子是土壤全鉀(X4)、速效鉀含量(X7),而土壤有機質(zhì)(X1)、全氮(X2)、全磷(X3)、速效氮(X5)、速效磷(X6)、和土壤pH值(X8)的影響則很小。同時從式(3)還可以看出,栽培大麥GSC與土壤全鉀(X4)、速效鉀含量(X7)呈負相關關系。根據(jù)標準誤差檢驗,該回歸方程未通過=0.05水平的顯著性檢驗。表明,栽培大麥GSC隨著土壤全鉀(X4)、速效鉀含量(X7)的減少而增加。同時也表明,土壤因子對栽培大麥GSC的整體影響未達到顯著性水平(P<0.05)。

    2.2.4 栽培大麥GSC與農(nóng)藝因子的關系

    基于逐步回歸分析法,建立栽培大麥GSC與株高(X1)、穗長(X2)、穗密度(X3)、穗形(X4)、穗姿(X5)、總小穗數(shù)(X6)、不孕小穗數(shù)(X7)、穗粒數(shù)(X8)、芒長(X9)、籽粒顏色(X10)、籽粒形狀(X11)、護穎寬度(X12)、千粒重(X13)等農(nóng)藝因子的數(shù)學模型:

    Y=1.7731+1.8331X3-3.8122X5+0.1279X6+1.0718X9-0.5077X10+0.1916X13

    (N=83R= 0.4717**R3= -0.3207**R5= -0.1885R6=0.1803R9= 0.3683R10= -0.1219R13=0.1683**)

    (4)

    從式(4)可以看出,影響栽培大麥GSC(Y)的主要農(nóng)藝因子是穗密度(X3)、穗姿(X5)、總小穗數(shù)(X6)、芒長(X9)、籽粒顏色(X10)、千粒重(X13),而株高(X1)、穗長(X2)、穗形(X4)、不孕小穗數(shù)(X7)、穗粒數(shù)(X8)、籽粒形狀(X11)、護穎寬度(X12)的影響則很小。同時從式(4)還可以看出,栽培大麥GSC與芒長(X9)的偏相關系數(shù)達到極顯著正相關差異水平,與穗密度(X3)的偏相關系數(shù)達到極顯著負相關差異水平,但是與穗姿(X5)、總小穗數(shù)(X6)、千粒重(X13)、籽粒顏色(X10)的偏相關關系均未達到顯著性(P<0.05)差異。根據(jù)標準誤差檢驗,該回歸方程通過了=0.01水平的顯著性檢驗。表明,栽培大麥GSC隨著芒長(X9)的增加和穗密度的減少而顯著增加。同時也表明,農(nóng)藝因子對栽培大麥GSC的整體影響達到極顯著性水平(P<0.01)。

    2.2.5 栽培大麥GSC與綜合因子的關系

    通過以上分析可以看出,栽培大麥GSC與年均日照時數(shù)、1月平均氣溫、6月平均氣溫、8月平均氣溫、8月平均氣溫日較差、8月平均相對濕度、6月平均月降水量、芒長的偏相關系數(shù)達到顯著或極顯著正相關差異水平,與≥0℃積溫、≥10℃積溫、5月平均氣溫、7月平均氣溫、9月平均氣溫、6月平均氣溫日較差、7月平均相對濕度、5月平均降水量、穗密度的偏相關系數(shù)達到顯著或極顯著負相關差異水平。但是,這些因子對栽培大麥GSC的影響程度如何尚不清楚。為此,我們運用隨機森林回歸分析法,對栽培大麥GSC影響顯著的年均日照時數(shù)、1月平均氣溫、6月平均氣溫、8月平均氣溫、8月平均氣溫日較差、8月平均相對濕度、6月平均月降水量、≥0℃積溫、≥10℃積溫、5月平均氣溫、7月平均氣溫、9月平均氣溫、6月平均氣溫日較差、7月平均相對濕度、5月平均降水量、穗密度、芒長等17個變量進行了隨機森林回歸分析。從表2可以看出,這17個變量對栽培大麥GSC影響從大到小的順序是穗密度﹥6月平均氣溫日較差>芒長﹥9月平均氣溫>1月平均氣溫>年日照時數(shù)﹥≥0℃積溫>5月平均氣溫﹥8月平均氣溫日較差>8月平均氣溫>6月平均氣溫>≥10℃積溫>6月平均月降水量>5月平均月降水量>7月平均相對濕度>8月平均相對濕度>7月平均氣溫。

    表2 不同因子對栽培大麥GSC影響的重要程度

    3 結(jié)論與討論

    在地理水平方向上,青藏高原栽培大麥GSC的水平分布總體呈現(xiàn)出斑塊狀交錯分布和南高北低的格局,并形成了以西藏拉孜、隆子、堆龍德慶、曲水、尼木、定日、薩迦、達孜、扎囊、日喀則、墨竹工卡、貢嘎、瓊杰為中心的青藏高原中南部和以青海共和、貴德、海晏、同德為中心的青藏高原東北部等2個栽培大麥GSC高值區(qū);在地理垂直方向上栽培大麥GSC的變化呈現(xiàn)出“S”型分布格局,即在海拔3300.0—3600.0 m以下, 栽培大麥GSC雖有波動性變化,但是總體含量很低,在海拔3300.0—3600.0 m以上,隨著海拔的升高,栽培大麥GSC逐漸增加,在海拔4200.0—4500.0 m達到最高值,此后隨著海拔的升高栽培大麥GSC略有下降。表明,今后要選擇高β淀粉的青藏高原栽培大麥品種,可能從這些區(qū)域的地方品種中進行選擇可能會有更大的收獲。要發(fā)展青藏高原青稞淀粉產(chǎn)業(yè),可能在這些地區(qū)進行規(guī)劃種植栽培大麥可能會有更好的收效。

    農(nóng)藝、氣候、土壤等因子對栽培大麥GSC影響從大到小的順序是穗密度﹥6月平均氣溫日較差>芒長﹥9月平均氣溫>1月平均氣溫>年日照時數(shù)﹥≥0℃積溫>5月平均氣溫﹥8月平均氣溫日較差>8月平均氣溫>6月平均氣溫>≥10℃積溫>6月平均月降水量>5月平均月降水量>7月平均相對濕度>8月平均相對濕度>7月平均氣溫。表明,對栽培大麥GSC影響最大的是基因型,其次是是氣候因素,是二者共同作用的結(jié)果,土壤因素對GSC的影響不明顯。影響栽培大麥GSC的農(nóng)藝因子主要是穗密度和芒長,氣候因子主要是拔節(jié)抽穗期的氣溫日較差和籽粒灌漿成熟期的平均氣溫,日照和降水的影響相對較小。這一研究結(jié)果支持大麥籽粒淀粉含量雖然主要受品種特性,即遺傳因素的控制,但環(huán)境因素對它也有明顯影響的研究結(jié)果[5- 14]。

    栽培大麥GSC與年均日照時數(shù)、6月平均月降水量的偏相關系數(shù)均達到顯著正相關差異水平,與9月份平均氣溫的偏相關系數(shù)均達到顯著負相關差異水平,與土壤含氮量的偏相關系數(shù)未達到顯著性差異水平。這與Erbs等對德國大麥的研究表明,隨著灌漿期溫度的升高和降水的減少,作物淀粉含量下降[8];Ahmed等對西藏大麥的研究表明,干旱條件下,大麥淀粉含量下降[9];Macleod和Duffus認為大麥籽粒發(fā)育中,英國大麥淀粉含量隨著溫度的升高而減少[25];Wei等對中國大麥的研究表明,隨著溫度的升高,淀粉含量減少[26]; Reinhardt等對德國大麥的研究研究表明,隨著溫度的升高,作物淀粉含量下降[27]的結(jié)論一致,但是與Holopainen等對英國二棱大麥(Hordeumvulgare)的研究表明,隨著日照時間的縮短,大麥淀粉含量增多,且二者差異未達到顯著水平[28];Erbs等對德國大麥的研究表明,隨著施氮量的增加,淀粉含量顯著下降[28]的結(jié)論相反。

    [1] Goldstein A, Annor G, Putaux J L, Hebelstrup K H, Blennow A, Bertoft E. Impact of full range of amylose contents on the architecture of starch granules. International Journal of Biological Macromolecules, 2016, 89: 305- 318.

    [2] Pérez S, Bertoft E. The molecular structures of starch components and Their contribution to the architecture of starch granules: A comprehensive review. Starch, 2010, 62(8): 389- 420.

    [3] Wang S J, Copeland L. Effect of acid hydrolysis on starch structure and functionality: a review. Critical Reviews in Food Science and Nutrition, 2015, 55(8): 1081- 1097.

    [4] Hermansson A M, Svegmark K S. Developments in the understanding of starch functionality. Trends in Food Science & Technology, 1996, 7(11): 345- 353.

    [5] 王建林. 西藏高原作物栽培學. 北京: 中國農(nóng)業(yè)出版社, 2012.

    [6] Griffey C, Brooks W, Kurantz M, Thomason W, Taylor F, Obert D, Moreau R, Flores R, Sohn M, Hicks K. Grain composition of Virginia winter barley and implications for use in feed, food, and biofuels production. Journal of Cereal Science, 2010, 51(1): 41- 49.

    [7] Fox G, Borgognone G M, Flinn P, Poulsen D. Genetic and environmental analysis of NIR 3. Feed quality predictions on genotypes of barley (HordeumvulgareL.). Field Crops Research, 2011, 120(3): 380- 386.

    [8] Erbs M, Manderscheid R, Jansen G, Seddig S, Pacholski A, Weigel H J. Effects of free-air CO2enrichment and nitrogen supply on grain quality parameters and elemental composition of wheat and barley grown in a crop rotation. Agriculture, Ecosystems & Environment, 2010, 136(1/2): 59- 68.

    [9] Ahmed I M, Cao F B, Han Y, Nadira U A, Zhang G P, Wu F B. Differential changes in grain ultrastructure, amylase, protein and amino acid profiles between Tibetan wild and cultivated barleys under drought and salinity alone and combined stress. Food Chemistry, 2013, 141(3): 2743- 2750.

    [10] Novo P, Vojislava M I, Jasmina S, Milan M. Effect of growing season and variety on quality of spring two-rowed barley. Genetika, 2014, 46(1): 59- 73.

    [11] Ahmed Z, Tetlow I J, Ahmed R, Morell M K, Emes M J. Protein interactions among enzymes of starch biosynthesis in high-amylose barley genotypes reveal differential roles of heteromeric enzyme complexes in the synthesis of A and B granules. Plant Science, 2015, 233: 95- 106.

    [12] Ahmadi A, Baker D A. The effect of water stress on the activities of key regulatory enzymes of the sucrose to starch pathway in wheat. Plant Growth Regulation, 2001, 35(1): 81- 91.

    [13] Ahmed I M, Dai H X, Zheng H X, Cao F B, Zhang G P, Sun D F, Wu F B. Genotypic differences in physiological characteristics in the tolerance to drought and salinity combined stress between Tibetan wild and cultivated barley. Plant Physiology and Biochemistry, 2013, 63: 49- 60.

    [14] Balla K, Rakszegi M, LI Z Y, Békés F, Bencze S, Veisz O. Quality of winter wheat in relation to heat and drought shock after Anthesis. Czech Journal of Food Science, 2011, 29(2): 117- 128.

    [15] Chen Z, Newman I, Zhou M, Mendham N, Zhang G, Shabala S. Screening plants for salt tolerance by measuring K+ux: a case study for barley. Plant, Cell & Environment, 2005, 28(10): 1230- 1246.

    [16] Brandt C, Balko C, Eichler-L?bermann B. Interactive effects of soil water content and phytin supply on phosphorus nutrition of different crops species. Landbauforschung, 2011, 61(4): 317- 325.

    [17] Clausen S K, Frenck G, Linden L G, Mikkelsen T N, Lunde C, J?rgensen R B. Effects of single and multifactor treatments with elevated temperature, CO2and ozone on oilseed rape and barley. Journal of Agronomy and Crop Science, 2011, 197(6): 442- 453.

    [18] Ferris R, Ellis R H, Wheeler T R, Hadley P. Effect of high temperature stress at anthesis on grain yield and biomass of field-grown crops of wheat. Annals of Botany, 1998, 82(5): 631- 639.

    [19] Gusmao M, Siddique K H M, Flower K, Nesbitt H, Veneklaas E J. Water deficit during the reproductive period of grass pea (LathyrussativusL.) reduced grain yield but maintained seed size. Journal of Agronomy Crop Science, 2012, 198(6): 430- 441.

    [20] del Moral G L F, Rharrabti Y R, Martos V, Royo C. Environmentally induced changes in amino acid composition in the grain of durum wheat grown under different water and temperature regimes in a Mediterranean environment. Journal of Agricultural and Food Chemistry, 2007, 55(20): 8144- 8151.

    [21] Goupy P, Hugues M B, Boivin P, Amiot M J. Antioxidant composition and activity of barley (Hordeumvulgare) and malt extracts and of isolated phenolic compounds. Journal of the Science of Food and Agriculture, 1999, 79(12): 1625- 1634.

    [22] Hejgaard J, Boisen S. High-lysine proteins in Hiproly barley breeding: Identification, nutritional significance and new screening methods. Hereditas, 2009, 93(2): 311- 320.

    [23] 馬得泉. 中國西藏大麥遺傳資源. 北京: 中國農(nóng)業(yè)出版社, 2000.

    [24] 王建林, 欒運芳, 大次卓嘎, 胡單. 西藏栽培大麥變種組成和分布規(guī)律研究. 中國農(nóng)業(yè)科學, 2006, 39(11): 2163- 2169.

    [25] Macleod L C, Duffus C M. Temperature effects on starch granules in developing barley grains. Journal of Cereal Science, 1988, 8(1): 29- 37.

    [26] Wei K, Dai F, Wu F B, Zhang G P. The Variation of β-amylase activity and protein fractions in barley grains as affected by genotypes and post-Anthesis temperatures. Journal of the Institute of Brewing, 2009, 115(3): 207- 213.

    [27] Reinhardt D, Jansen G, Seddig S, Eichler-Loebermann B. Temperature stress during flowering time affects yield and quality parameters of waxy barley. Landbauforschung Volkenrode, 2013, 63(1): 79- 83.

    [28] Holopainen U R M, Wilhelmson A, Home S, Poutanen K, Shewry P R. Day-length effects on protein localisation affect water absorption in barley (Hordeumvulgare) grains. Food and Agriculture, 2012, 92(15): 2944- 2951.

    EffectsofagronomicandclimaticfactorsonthespatialdistributionofstarchconcentrationinbarleycultivatedontheQinghai-TibetPlateau

    ZHONG Zhiming1, WANG Jianlin2, *, FENG Xibo2, FU Gang1, HOU Weihai2, WANG Gaihua2, Dacizhuoga2

    1InstituteofGeographicSciencesandNaturalResourcesResearch,ChineseAcademyofSciences,Beijing100101,China2CollegeofPlantScience,TibetInstituteofAgricultureandAnimalHusbandry,Linzhi860000,China

    Starch concentration is an important quality character of barley grain. In recent years, some researchers have realized that barley grain starch concentration (GSC) is affected by genetic factors, as well as by its cultivation and environmental conditions. However, the effect of different environmental factors on barley GSC remains unclear. In addition, previous studies have mostly been conducted at several sites, and large-scale studies have not addressed the relationship between cultivated barley GSC and environment factors. The Qinghai-Tibet Plateau is an important region for growing hull-less barley in China, and the relationship between barley GSC and environment factors on the Qinghai-Tibet Plateau have not been reported. Therefore, the present study investigated the distribution of barley GSC on the Qinghai-Tibet Plateau using geography, climate, soil, and agronomy data from 83 sampling sites, the cultivated barley GSC distribution characteristics in Qinghai-Tibet Plateau was studied in this research work. The aim of the study was to identify the relative effects of different environmental factors on GSC, to clarify the relationship between the distribution of barley GSC and environmental factors, and to get explicit GSC acclimation of barley varieties planted in different areas to environment. The results showed that (1) GSC exhibited a macrocosm pattern of spatial distribution along the horizontal direction on the Qinghai-Tibet Plateau, with staggered patches of different value and a tendency of “high in the south and low in the north.” In addition, two regions of higher GSC (south-central and north-eastern) included the Lazi, Longzi, Duilongdeiqin, Qushui, Nimu, Dingri, Sajia, Dazi, Zhanang, Shigatse, Mozhugongka, Gongga, and Qiongjie Counties in Tibet and the Gonghe, Guide, Haiyan, and Tongde Counties in the Qinghai Province. (2) Along the vertical direction, GSC exhibited an S-shaped distribution pattern. Indeed, GSC increased gradually with increasing altitude from 3300.0—3600.0 m, reached its maximum between 4200.0 m and 4500.0 m, and thereafter slightly decreased with increasing altitude. (3) The relative effects of the environmental factors on barley GSC were as follows: spike density > average diurnal temperature range in June > length of the awn > average temperature in September > average temperature in January > annual sunshine hours > accumulated temperature (≥0 ℃) > average temperature in May > average diurnal temperature range in August > average temperature in August > average temperature in June > accumulated temperature (≥10 ℃) > average precipitation in June > average precipitation in May > average relative humidity in July > average relative humidity in August > average temperature in July. Therefore, genotype had the greatest impact on the GSC of cultivated barley, followed by climatic factors, and the effects of soil factors effect were not apparent. Spike density and length of the awn were the most important agronomic factors, and diurnal temperature range during the period of elongating and heading and average temperature during the period of filling and maturation were the main climatic factors, whereas the effects of sunshine hours and precipitation were relatively small. These findings support previous reports that barley GSC is mainly controlled by the characteristics (i.e., genetic factors) of individual barley varieties, although the effects of environmental factors were also apparent.

    agronomic factor; climatic factor; barley; starch concentration; spatial distribution; Qinghai-Tibet Plateau

    國家自然科學基金 (31360300,31370458,31560362); 西藏自治區(qū)重點項目(XZXTCX- 2015-03,XBTSZWXK- 2015-01);西藏農(nóng)牧學院雪域英才工程人才發(fā)展支持計劃項目(2015XYB01);作物學科建設項目(2015ZWXKJS)

    2016- 08- 31; < class="emphasis_bold">網(wǎng)絡出版日期

    日期:2017- 07- 11

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

    10.5846/stxb201608311782

    鐘志明,王建林,馮西博,付剛,侯維海,王改花,大次卓嘎.農(nóng)藝與氣候因素對青藏高原栽培大麥淀粉含量空間分布的影響.生態(tài)學報,2017,37(21):7197- 7205.

    Zhong Z M, Wang J L, Feng X B, Fu G, Hou W H, Wang G H, Dacizhuoga.Effects of agronomic and climatic factors on the spatial distribution of starch concentration in barley cultivated on the Qinghai-Tibet Plateau.Acta Ecologica Sinica,2017,37(21):7197- 7205.

    猜你喜歡
    日較差大麥平均氣溫
    我的大麥哥哥
    1961—2018年黃岡氣溫日較差變化特征及其影響因子分析
    大麥蟲對聚苯乙烯塑料的生物降解和礦化作用
    Global interest in Chinese baijiu
    1961—2018 年我國氣溫日較差日數(shù)的時空演變特征及區(qū)域差異
    烏蘭縣近38年氣溫變化特征分析
    石河子墾區(qū)氣溫日較差的變化特征*
    從全球氣候變暖大背景看萊州市30a氣溫變化
    1981—2010年拐子湖地區(qū)氣溫變化特征及趨勢分析
    近50年來全球背景下青藏高原氣候變化特征分析
    99久国产av精品国产电影| 亚洲欧美中文字幕日韩二区| 免费日韩欧美在线观看| 老司机亚洲免费影院| 最近2019中文字幕mv第一页| 日韩 亚洲 欧美在线| 亚洲五月色婷婷综合| 麻豆精品久久久久久蜜桃| 韩国高清视频一区二区三区| 午夜福利乱码中文字幕| 亚洲精品视频女| 亚洲av.av天堂| 精品第一国产精品| 欧美丝袜亚洲另类| 岛国毛片在线播放| 草草在线视频免费看| 一区二区三区四区激情视频| 成年av动漫网址| 欧美日本中文国产一区发布| 男人操女人黄网站| 亚洲精品av麻豆狂野| 一级黄片播放器| 久久精品aⅴ一区二区三区四区 | 五月开心婷婷网| 两性夫妻黄色片 | 又黄又爽又刺激的免费视频.| 亚洲经典国产精华液单| 亚洲国产av影院在线观看| 午夜福利视频精品| 91成人精品电影| 欧美最新免费一区二区三区| 99久久中文字幕三级久久日本| 亚洲高清免费不卡视频| 亚洲国产精品国产精品| 最近中文字幕2019免费版| 国产麻豆69| 日日爽夜夜爽网站| 国产亚洲精品第一综合不卡 | 亚洲欧美中文字幕日韩二区| 宅男免费午夜| 国产 一区精品| 青青草视频在线视频观看| 麻豆乱淫一区二区| 欧美日本中文国产一区发布| 哪个播放器可以免费观看大片| 制服人妻中文乱码| 日本欧美国产在线视频| 成人午夜精彩视频在线观看| 国产精品 国内视频| 午夜久久久在线观看| 在线 av 中文字幕| 精品视频人人做人人爽| 国产极品天堂在线| 久久午夜综合久久蜜桃| 欧美亚洲日本最大视频资源| 色94色欧美一区二区| 国产极品粉嫩免费观看在线| 性色avwww在线观看| 亚洲精品一区蜜桃| 制服丝袜香蕉在线| 91久久精品国产一区二区三区| 亚洲欧美一区二区三区国产| 熟女av电影| 亚洲精品乱码久久久久久按摩| 日本猛色少妇xxxxx猛交久久| 亚洲av免费高清在线观看| 巨乳人妻的诱惑在线观看| 国产免费福利视频在线观看| 丰满乱子伦码专区| 狂野欧美激情性bbbbbb| 在线精品无人区一区二区三| 草草在线视频免费看| 国产免费一区二区三区四区乱码| 欧美成人午夜免费资源| 亚洲成人av在线免费| 亚洲av.av天堂| 高清欧美精品videossex| 国产xxxxx性猛交| 国产福利在线免费观看视频| 中文字幕最新亚洲高清| 18禁动态无遮挡网站| 黄色一级大片看看| 日本wwww免费看| 精品国产国语对白av| 国产日韩欧美亚洲二区| 国产精品久久久久久久电影| 午夜久久久在线观看| 少妇人妻精品综合一区二区| 97超碰精品成人国产| 免费高清在线观看视频在线观看| 亚洲精品久久午夜乱码| av在线app专区| 你懂的网址亚洲精品在线观看| 五月玫瑰六月丁香| 欧美成人午夜精品| 成人亚洲精品一区在线观看| 亚洲av福利一区| a级毛片黄视频| 精品人妻在线不人妻| 亚洲精品,欧美精品| 久久99热这里只频精品6学生| 国产精品蜜桃在线观看| 免费高清在线观看日韩| 九色成人免费人妻av| 少妇 在线观看| 熟女av电影| 国产熟女午夜一区二区三区| 人妻人人澡人人爽人人| 亚洲国产毛片av蜜桃av| 亚洲激情五月婷婷啪啪| 日日爽夜夜爽网站| 亚洲国产看品久久| 国产精品国产三级专区第一集| 亚洲情色 制服丝袜| 久热久热在线精品观看| 永久网站在线| 亚洲中文av在线| 91午夜精品亚洲一区二区三区| 亚洲av电影在线观看一区二区三区| 黄色 视频免费看| 亚洲成色77777| 成人国产av品久久久| 伦理电影免费视频| 热re99久久精品国产66热6| 在线天堂中文资源库| 少妇人妻久久综合中文| 亚洲精品第二区| 国产精品不卡视频一区二区| 欧美日韩一区二区视频在线观看视频在线| 黄色一级大片看看| 日本av免费视频播放| 国产男人的电影天堂91| 高清不卡的av网站| 美女视频免费永久观看网站| 欧美激情 高清一区二区三区| 国产福利在线免费观看视频| 极品人妻少妇av视频| 日韩中字成人| 国产男女内射视频| 日韩成人av中文字幕在线观看| 男男h啪啪无遮挡| 性色av一级| 久久这里有精品视频免费| 老司机影院毛片| 免费观看a级毛片全部| 中国国产av一级| av女优亚洲男人天堂| 内地一区二区视频在线| 欧美激情国产日韩精品一区| 母亲3免费完整高清在线观看 | 亚洲av国产av综合av卡| 精品国产国语对白av| 国产日韩欧美在线精品| av福利片在线| 欧美日韩av久久| 亚洲国产精品999| 日韩精品有码人妻一区| 免费黄色在线免费观看| 性高湖久久久久久久久免费观看| 2018国产大陆天天弄谢| 亚洲精品国产av蜜桃| 人人妻人人澡人人爽人人夜夜| 日韩中字成人| 自拍欧美九色日韩亚洲蝌蚪91| 国产色婷婷99| 欧美激情 高清一区二区三区| 美女脱内裤让男人舔精品视频| 国产成人精品一,二区| 国产精品不卡视频一区二区| 满18在线观看网站| 亚洲美女视频黄频| 日韩成人av中文字幕在线观看| 亚洲精品成人av观看孕妇| 精品一品国产午夜福利视频| av.在线天堂| 极品人妻少妇av视频| 在线观看免费高清a一片| 丝袜在线中文字幕| 久久av网站| 汤姆久久久久久久影院中文字幕| 多毛熟女@视频| 国产精品一区www在线观看| 午夜影院在线不卡| 寂寞人妻少妇视频99o| 久久久久久伊人网av| 久久精品夜色国产| 国产熟女欧美一区二区| 亚洲精品久久午夜乱码| 在线观看免费视频网站a站| 制服诱惑二区| 一二三四在线观看免费中文在 | 色婷婷av一区二区三区视频| 性色av一级| 99精国产麻豆久久婷婷| 日本av手机在线免费观看| 国产亚洲精品第一综合不卡 | 乱人伦中国视频| 精品福利永久在线观看| 欧美日韩国产mv在线观看视频| 精品人妻在线不人妻| 久久久久精品久久久久真实原创| 韩国av在线不卡| 色网站视频免费| 国产亚洲最大av| 韩国精品一区二区三区 | 三上悠亚av全集在线观看| 永久免费av网站大全| 亚洲精品久久午夜乱码| 午夜免费男女啪啪视频观看| 日韩三级伦理在线观看| 三上悠亚av全集在线观看| 男的添女的下面高潮视频| 国产1区2区3区精品| 人妻 亚洲 视频| 成人国产麻豆网| 国产精品一区www在线观看| 国产一区二区激情短视频 | 久久精品夜色国产| 亚洲国产精品一区二区三区在线| 精品少妇内射三级| 91精品三级在线观看| 日韩三级伦理在线观看| 亚洲av在线观看美女高潮| 欧美精品人与动牲交sv欧美| 婷婷色综合www| 日韩在线高清观看一区二区三区| 精品午夜福利在线看| 18禁动态无遮挡网站| 国产精品 国内视频| 精品国产国语对白av| 亚洲,欧美精品.| 国产一区二区在线观看av| 韩国av在线不卡| 在线天堂最新版资源| 黄色怎么调成土黄色| 国产又爽黄色视频| 99热网站在线观看| 18禁观看日本| 国产精品一国产av| 久久人人爽人人片av| 国产av精品麻豆| 精品第一国产精品| 少妇被粗大的猛进出69影院 | 狠狠精品人妻久久久久久综合| 男的添女的下面高潮视频| 精品卡一卡二卡四卡免费| 国产又色又爽无遮挡免| 又大又黄又爽视频免费| 春色校园在线视频观看| 99久久中文字幕三级久久日本| 人人妻人人添人人爽欧美一区卜| 秋霞在线观看毛片| 久久国产精品男人的天堂亚洲 | 精品国产乱码久久久久久小说| 精品少妇久久久久久888优播| 亚洲精品国产av蜜桃| 国产男女超爽视频在线观看| 日日撸夜夜添| www日本在线高清视频| 精品人妻在线不人妻| 日产精品乱码卡一卡2卡三| 青青草视频在线视频观看| 纯流量卡能插随身wifi吗| 国产男女超爽视频在线观看| 亚洲一级一片aⅴ在线观看| 美女国产高潮福利片在线看| a级片在线免费高清观看视频| 最黄视频免费看| a级毛片黄视频| 免费av不卡在线播放| 亚洲欧洲日产国产| 欧美丝袜亚洲另类| 国产高清国产精品国产三级| 亚洲 欧美一区二区三区| 欧美人与性动交α欧美软件 | 国产有黄有色有爽视频| 国产男人的电影天堂91| 亚洲欧美一区二区三区国产| 中国美白少妇内射xxxbb| 欧美日韩精品成人综合77777| av国产久精品久网站免费入址| 国产精品麻豆人妻色哟哟久久| 亚洲精品aⅴ在线观看| 中国国产av一级| 18禁观看日本| 欧美日韩一区二区视频在线观看视频在线| 免费不卡的大黄色大毛片视频在线观看| 国产老妇伦熟女老妇高清| 精品一区二区三卡| 欧美最新免费一区二区三区| 国产精品人妻久久久影院| 各种免费的搞黄视频| 肉色欧美久久久久久久蜜桃| 国产一区二区在线观看av| 丝袜美足系列| 久久精品国产a三级三级三级| 中文字幕另类日韩欧美亚洲嫩草| av又黄又爽大尺度在线免费看| 在线免费观看不下载黄p国产| 捣出白浆h1v1| 七月丁香在线播放| 欧美xxxx性猛交bbbb| 捣出白浆h1v1| 少妇的丰满在线观看| 熟女人妻精品中文字幕| 欧美+日韩+精品| 亚洲综合色惰| 最近的中文字幕免费完整| 99国产精品免费福利视频| 亚洲av.av天堂| 99热全是精品| 少妇高潮的动态图| 黄色毛片三级朝国网站| 精品酒店卫生间| 观看美女的网站| 日韩在线高清观看一区二区三区| 国产精品一区二区在线不卡| 少妇的逼水好多| 黑人高潮一二区| 男女午夜视频在线观看 | 欧美丝袜亚洲另类| 国产精品无大码| 免费黄色在线免费观看| 婷婷色麻豆天堂久久| 精品一区二区三区视频在线| 男女免费视频国产| 国产精品嫩草影院av在线观看| 久久国产精品大桥未久av| 卡戴珊不雅视频在线播放| 十八禁网站网址无遮挡| 看免费av毛片| 国产毛片在线视频| 久久人人爽av亚洲精品天堂| 久久亚洲国产成人精品v| 男女边吃奶边做爰视频| 久久毛片免费看一区二区三区| 日韩av免费高清视频| 久久久久久久亚洲中文字幕| 各种免费的搞黄视频| av在线app专区| 国产麻豆69| 国产日韩欧美视频二区| 欧美最新免费一区二区三区| av又黄又爽大尺度在线免费看| 国产成人aa在线观看| 久久久久精品人妻al黑| 色婷婷av一区二区三区视频| 久久久久精品性色| 高清视频免费观看一区二区| 一级毛片 在线播放| 韩国av在线不卡| 老司机影院毛片| av天堂久久9| videos熟女内射| 欧美精品一区二区免费开放| 亚洲精品国产av成人精品| 人人妻人人爽人人添夜夜欢视频| av免费在线看不卡| 建设人人有责人人尽责人人享有的| 欧美日韩综合久久久久久| 伦理电影免费视频| 精品人妻一区二区三区麻豆| 亚洲国产av影院在线观看| 久久人人97超碰香蕉20202| 视频区图区小说| 在线亚洲精品国产二区图片欧美| 一区二区av电影网| 精品亚洲成国产av| 日韩中文字幕视频在线看片| 久久婷婷青草| 国产高清国产精品国产三级| 亚洲精华国产精华液的使用体验| 中国国产av一级| 亚洲,欧美精品.| 久久久久人妻精品一区果冻| 91精品三级在线观看| 久久午夜综合久久蜜桃| 国产探花极品一区二区| 国产1区2区3区精品| 国产白丝娇喘喷水9色精品| 黄色毛片三级朝国网站| 久久精品人人爽人人爽视色| 最近手机中文字幕大全| 亚洲av日韩在线播放| 亚洲av免费高清在线观看| 最近最新中文字幕大全免费视频 | 亚洲欧美精品自产自拍| 色视频在线一区二区三区| 在线观看人妻少妇| 免费观看a级毛片全部| 美女国产高潮福利片在线看| 岛国毛片在线播放| 男女国产视频网站| 久久精品国产鲁丝片午夜精品| 丝袜人妻中文字幕| 一本久久精品| 欧美激情极品国产一区二区三区 | 成人午夜精彩视频在线观看| 亚洲伊人久久精品综合| 激情视频va一区二区三区| 成人毛片60女人毛片免费| 一二三四中文在线观看免费高清| 国产精品久久久久成人av| 国产又色又爽无遮挡免| 精品午夜福利在线看| 久久人人97超碰香蕉20202| 日产精品乱码卡一卡2卡三| 午夜老司机福利剧场| 80岁老熟妇乱子伦牲交| 亚洲国产av影院在线观看| 最近2019中文字幕mv第一页| 精品酒店卫生间| 国产日韩一区二区三区精品不卡| 国产精品久久久av美女十八| 国产国语露脸激情在线看| 亚洲欧美一区二区三区国产| 人体艺术视频欧美日本| 久久人人爽av亚洲精品天堂| 精品一品国产午夜福利视频| 国产成人91sexporn| 国产精品久久久久久精品古装| 国产成人精品无人区| 毛片一级片免费看久久久久| 亚洲成av片中文字幕在线观看 | 视频区图区小说| av在线app专区| 国产男人的电影天堂91| 97超碰精品成人国产| 日韩成人av中文字幕在线观看| 少妇高潮的动态图| 欧美 亚洲 国产 日韩一| 纵有疾风起免费观看全集完整版| 国产日韩欧美视频二区| 精品国产乱码久久久久久小说| 久久久久久人妻| 成人国产av品久久久| 建设人人有责人人尽责人人享有的| 成人黄色视频免费在线看| 国产爽快片一区二区三区| 亚洲美女黄色视频免费看| 丰满少妇做爰视频| av卡一久久| 精品亚洲成国产av| 香蕉国产在线看| 少妇人妻精品综合一区二区| 精品福利永久在线观看| 好男人视频免费观看在线| 不卡视频在线观看欧美| av又黄又爽大尺度在线免费看| 日韩制服丝袜自拍偷拍| 亚洲精品视频女| a 毛片基地| 天堂中文最新版在线下载| 国产xxxxx性猛交| 亚洲av男天堂| 久久婷婷青草| 亚洲av电影在线观看一区二区三区| 看免费成人av毛片| 我要看黄色一级片免费的| 精品福利永久在线观看| 国产精品久久久久久久久免| 亚洲一级一片aⅴ在线观看| 亚洲丝袜综合中文字幕| 国产成人av激情在线播放| 久久久国产精品麻豆| 高清欧美精品videossex| 美女大奶头黄色视频| 国产69精品久久久久777片| 婷婷色综合大香蕉| 亚洲国产毛片av蜜桃av| 国产日韩一区二区三区精品不卡| 性高湖久久久久久久久免费观看| 99久国产av精品国产电影| 2022亚洲国产成人精品| 亚洲图色成人| 亚洲内射少妇av| 91精品三级在线观看| 咕卡用的链子| 婷婷色麻豆天堂久久| 久久久久久久国产电影| 久久久国产一区二区| 久久av网站| 国产免费一区二区三区四区乱码| 一级片免费观看大全| 黄色 视频免费看| 国产精品国产三级国产av玫瑰| 午夜视频国产福利| 久久韩国三级中文字幕| 成人毛片a级毛片在线播放| 国产无遮挡羞羞视频在线观看| 少妇人妻 视频| 老女人水多毛片| 七月丁香在线播放| 人妻系列 视频| 乱人伦中国视频| 欧美日韩视频精品一区| 丰满饥渴人妻一区二区三| 大话2 男鬼变身卡| av卡一久久| 亚洲经典国产精华液单| 91精品国产国语对白视频| 黄色一级大片看看| 在线观看美女被高潮喷水网站| kizo精华| 91久久精品国产一区二区三区| 成年女人在线观看亚洲视频| 国产成人精品婷婷| 国产淫语在线视频| 一级片'在线观看视频| 一级爰片在线观看| 国产国拍精品亚洲av在线观看| 美女中出高潮动态图| 亚洲内射少妇av| 亚洲精品久久成人aⅴ小说| 午夜久久久在线观看| av视频免费观看在线观看| 精品一区在线观看国产| 久久这里有精品视频免费| 免费av中文字幕在线| 久久久精品免费免费高清| www.熟女人妻精品国产 | 久久这里只有精品19| 亚洲国产精品国产精品| 亚洲av男天堂| 国产精品麻豆人妻色哟哟久久| 18禁裸乳无遮挡动漫免费视频| 午夜免费观看性视频| 亚洲色图综合在线观看| 在线看a的网站| 蜜桃国产av成人99| 亚洲综合色网址| 看免费成人av毛片| a 毛片基地| 日韩不卡一区二区三区视频在线| 午夜视频国产福利| 久久精品国产综合久久久 | av线在线观看网站| 伦理电影免费视频| 九九在线视频观看精品| 激情视频va一区二区三区| 最黄视频免费看| 高清在线视频一区二区三区| 国产成人免费观看mmmm| 18禁观看日本| 日韩不卡一区二区三区视频在线| 另类精品久久| 在线观看免费日韩欧美大片| 97精品久久久久久久久久精品| 日韩一本色道免费dvd| a级毛色黄片| av福利片在线| 日日啪夜夜爽| 精品久久蜜臀av无| 免费人成在线观看视频色| 男女国产视频网站| 91精品三级在线观看| 黄色毛片三级朝国网站| 成人亚洲欧美一区二区av| 欧美人与性动交α欧美软件 | 欧美日韩av久久| 欧美日韩成人在线一区二区| 热re99久久精品国产66热6| 国产一区二区在线观看日韩| 高清在线视频一区二区三区| av又黄又爽大尺度在线免费看| 久热这里只有精品99| 国产精品一区二区在线不卡| 亚洲精品美女久久久久99蜜臀 | 蜜桃在线观看..| 最近手机中文字幕大全| 亚洲美女黄色视频免费看| 精品一区二区三卡| 热re99久久国产66热| 蜜臀久久99精品久久宅男| 午夜影院在线不卡| 久久 成人 亚洲| 这个男人来自地球电影免费观看 | 九九爱精品视频在线观看| 欧美精品av麻豆av| 熟女电影av网| 黑人猛操日本美女一级片| 亚洲成色77777| 天堂俺去俺来也www色官网| 免费不卡的大黄色大毛片视频在线观看| 国产精品人妻久久久久久| 99久久精品国产国产毛片| 最黄视频免费看| 国产又色又爽无遮挡免| 国产亚洲精品久久久com| 美女主播在线视频| av视频免费观看在线观看| 亚洲欧美清纯卡通| 波多野结衣一区麻豆| 精品人妻在线不人妻| 午夜老司机福利剧场| 午夜av观看不卡| 久久久精品免费免费高清| 亚洲综合精品二区| 哪个播放器可以免费观看大片| 日韩欧美一区视频在线观看| 人妻人人澡人人爽人人| 美女主播在线视频| 国产不卡av网站在线观看| 99久国产av精品国产电影| 亚洲欧美成人综合另类久久久| 亚洲欧美色中文字幕在线| 搡老乐熟女国产| 国产在线免费精品| 一区二区三区乱码不卡18| 满18在线观看网站| 国产乱来视频区| 日韩中字成人| 老女人水多毛片| 久久影院123| 少妇高潮的动态图|