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

    基于多元統(tǒng)計(jì)分析的小麥低溫凍害評(píng)價(jià)及水分效應(yīng)差異研究

    2022-04-14 01:47:14王洋洋劉萬代賀利任德超段劍釗胡新郭天財(cái)王永華馮偉
    中國農(nóng)業(yè)科學(xué) 2022年7期
    關(guān)鍵詞:抗寒性灌水生理

    王洋洋,劉萬代,賀利,任德超,段劍釗,胡新,郭天財(cái),王永華,馮偉

    基于多元統(tǒng)計(jì)分析的小麥低溫凍害評(píng)價(jià)及水分效應(yīng)差異研究

    王洋洋1,2,劉萬代1,2,賀利1,2,任德超3,段劍釗1,2,胡新3,郭天財(cái)1,2,王永華1,2,馮偉1,2

    1河南農(nóng)業(yè)大學(xué)農(nóng)學(xué)院/作物生長發(fā)育調(diào)控教育部重點(diǎn)實(shí)驗(yàn)室,鄭州 450046;2國家小麥工程技術(shù)研究中心,鄭州 450046;3商丘市農(nóng)林科學(xué)院小麥研究所,河南商丘 476000

    【】明確不同水分條件下低溫脅迫引起的小麥凍害程度,篩選凍害鑒定指標(biāo)并建立凍害定量評(píng)估模型,為小麥生產(chǎn)科學(xué)防控低溫凍害提供理論支撐。以弱春性品種偃展4110、蘭考198和半冬性品種鄭麥366、豐德存麥21為試驗(yàn)材料,在模擬凍害發(fā)生前一周進(jìn)行灌水(W)和不灌水(D)處理,于雌雄蕊分化期將盆栽小麥移到低溫模擬室進(jìn)行處理,設(shè)置的溫度為-2℃(T1)、-4℃(T2)、-6℃(T3)、-8℃(T4)和-10℃(T5)以及對(duì)照(CK為當(dāng)天大田溫度),低溫脅迫后的第2天測(cè)定小麥生理生化指標(biāo),將標(biāo)準(zhǔn)化后的各個(gè)生理指標(biāo)進(jìn)行主成分、隸屬函數(shù)、聚類分析和逐步回歸等多元統(tǒng)計(jì)分析。不同品種、水分、溫度下的各單項(xiàng)生理生化指標(biāo)之間大多存在著顯著相關(guān)性,通過主成分分析將19個(gè)生理生化指標(biāo)轉(zhuǎn)化為6個(gè)相互獨(dú)立的綜合指標(biāo),其貢獻(xiàn)度分別為55.972%、11.93%、7.168%、5.075%、4.236%和3.079%,代表了全部原始數(shù)據(jù)的87.459%的信息量,并根據(jù)隸屬函數(shù)算法求出各處理的凍害程度綜合評(píng)價(jià)值(值)。以值作為因變量,通過逐步回歸分析篩選出7個(gè)關(guān)鍵指標(biāo),分別為葉綠素a、葉片含水量、脯氨酸、v/m、可溶性蛋白、MDA和SOD,并確立了定量估算值的數(shù)學(xué)模型。同時(shí),將預(yù)測(cè)值與產(chǎn)量損失率進(jìn)行相關(guān)性分析,線性方程決定系數(shù)2= 0.898,表明該預(yù)測(cè)模型能夠很好地評(píng)價(jià)凍害程度。進(jìn)一步對(duì)預(yù)測(cè)值進(jìn)行聚類分析,可將不同凍害處理劃分為5類:未受凍(D-CK、W-CK)、輕度受凍(D-T1、W-T1)、中度受凍(D-T2、W-T2、W-T3)、重度受凍(D-T3、W-T4)以及特重受凍(D-T4、W-T5、D-T5),其產(chǎn)量損失率分別為0、0—10%、10%—30%、30%—50%及50%以上,相同的溫度和水分條件下弱春性品種的凍害程度重于半冬性品種,相同的品種和溫度條件下不灌水處理的凍害程度重于灌水處理。具體考察所篩選的生理評(píng)價(jià)指標(biāo),隨著低溫脅迫的加重,葉綠素a、葉片含水量、v/m表現(xiàn)下降趨勢(shì),而脯氨酸、可溶性蛋白和SOD活性表現(xiàn)先升后降的特征,MDA則為相反趨勢(shì)。生產(chǎn)中在晚霜凍易發(fā)地區(qū)應(yīng)選用半冬性品種,且依據(jù)天氣預(yù)報(bào)寒流來臨前加強(qiáng)灌水管理,并可在凍害發(fā)生時(shí)通過凍害評(píng)價(jià)指標(biāo)及定量模型及時(shí)準(zhǔn)確評(píng)估凍害發(fā)生程度,這有利于晚霜凍害的科學(xué)防控,為災(zāi)后產(chǎn)量恢復(fù)及決策管理提供依據(jù)。

    冬小麥;晚霜凍害;灌水;綜合評(píng)價(jià);評(píng)估模型

    0 引言

    【研究意義】隨著溫室效應(yīng)的加劇,全球氣溫不斷升高。在全球變暖的背景下,極端天氣發(fā)生的次數(shù)、強(qiáng)度也在不斷增加[1]。低溫頻發(fā)增加了冬小麥生產(chǎn)系統(tǒng)的不穩(wěn)定性,加劇了冬小麥產(chǎn)量損失風(fēng)險(xiǎn)[2]。小麥拔節(jié)后霜凍被稱為晚霜凍害,在我國小麥重要產(chǎn)區(qū)黃淮海地區(qū)頻繁發(fā)生,此時(shí)小麥抗凍性相比前期明顯變?nèi)?,產(chǎn)量損失風(fēng)險(xiǎn)加大,凍害一旦發(fā)生,產(chǎn)量損失可達(dá)15%—60%[3],嚴(yán)重影響我國糧食安全。因此,及時(shí)準(zhǔn)確評(píng)估晚霜凍害對(duì)冬小麥產(chǎn)量的影響,以便盡早采取補(bǔ)救措施,對(duì)挽回因凍害造成的產(chǎn)量損失具有重要意義。【前人研究進(jìn)展】研究表明,逆境脅迫對(duì)植物生長產(chǎn)生重要影響,主要通過影響滲透調(diào)節(jié)物質(zhì)、抗氧化物酶等生理指標(biāo)[3-5],進(jìn)而對(duì)葉綠體、光合結(jié)構(gòu)造成破壞,從而導(dǎo)致葉綠素降解,葉綠素?zé)晒馊鏔v/Fm和qP下降[6-8]。低溫脅迫是一種常見的逆境,植物對(duì)低溫的適應(yīng)與抗性建成是一個(gè)復(fù)雜的生物學(xué)過程,單一的生理指標(biāo)變化無法準(zhǔn)確評(píng)估逆境對(duì)植株的影響[9]。目前大多數(shù)研究是從生理學(xué)和形態(tài)學(xué)指標(biāo)等方面對(duì)逆境脅迫進(jìn)行綜合性評(píng)價(jià)[10-11]。鄒成林等[12]利用指標(biāo)之間的相關(guān)性及評(píng)分分級(jí)法,綜合評(píng)價(jià)玉米苗期的抗旱性,并提出植株含水量、葉綠素含量、酶活性、可溶性糖含量、MDA含量可作為抗旱性鑒定指標(biāo)。袁雨豪等[13]研究了鹽脅迫對(duì)糜子生理指標(biāo)的影響,利用主成分及聚類分析將多種生理指標(biāo)進(jìn)行綜合,進(jìn)而對(duì)耐鹽性進(jìn)行分類,篩選出高度耐鹽品種。李春喜等[14]通過主成分分析、隸屬函數(shù)法和聚類分析法相結(jié)合,篩選出不同抗旱性小麥品種,并對(duì)抗旱性強(qiáng)弱進(jìn)行分類。盧廣超等[15]對(duì)熒光參數(shù)及生理指標(biāo)進(jìn)行主成分分析,實(shí)現(xiàn)了4種常用植物的抗寒性分類。以上研究表明,前人利用生理指標(biāo)在植物逆境脅迫下的響應(yīng)變化,并采用主成分分析、隸屬函數(shù)法等分析方法綜合評(píng)估植株的抗逆性。有關(guān)作物晚霜凍害的研究也多利用多元統(tǒng)計(jì)分析方法對(duì)作物抗寒性進(jìn)行評(píng)價(jià)[16]。王樹剛等[9]發(fā)現(xiàn),SOD、POD、CAT、MDA和可溶性蛋白可以作為小麥品種抗凍性鑒定的指標(biāo),并利用主成分分析和聚類分析法對(duì)品種進(jìn)行抗寒性篩選。姜麗娜等[17]對(duì)小麥拔節(jié)期進(jìn)行低溫處理,提出了滲透調(diào)節(jié)物質(zhì)、抗氧化物酶等生理指標(biāo)可作為抗寒性鑒定指標(biāo),并利用隸屬函數(shù)法和極點(diǎn)排序法將不同品種的抗寒性進(jìn)行強(qiáng)弱劃分。【本研究切入點(diǎn)】前人相關(guān)研究多以單低溫處理、多品種材料為對(duì)象,利用主成分分析方法篩選種質(zhì)資源,以及評(píng)價(jià)不同品種的抗寒性強(qiáng)弱。同時(shí)寒流發(fā)生前灌水是預(yù)防小麥凍害、減輕危害的最常用的生產(chǎn)措施。但不同水分條件及低溫水平對(duì)作物凍害程度及產(chǎn)量損失的效應(yīng)是不同的,而有關(guān)這方面的研究較少。因此積極探索并明確不同水分、低溫及品種多種因素條件下晚霜凍害對(duì)小麥產(chǎn)量的影響程度及其評(píng)價(jià)方法,這對(duì)于提前評(píng)估因凍害引起的植株傷害及其產(chǎn)量損失,并及時(shí)采取預(yù)防或補(bǔ)救措施具有重要的參考價(jià)值?!緮M解決的關(guān)鍵問題】本研究設(shè)置了模擬凍害發(fā)生前灌水和不灌水兩種水分條件,選用具有代表性且生產(chǎn)上大面積推廣應(yīng)用的不同類型小麥品種,測(cè)定不同低溫處理下植株生理生化指標(biāo)的響應(yīng)變化,通過主成分分析法、隸屬函數(shù)法及聚類分析方法相結(jié)合綜合評(píng)估小麥凍害程度,進(jìn)而借助逐步回歸分析篩選出適宜的抗寒性評(píng)價(jià)指標(biāo),以期為災(zāi)后產(chǎn)量恢復(fù)以及應(yīng)時(shí)決策管理提供科學(xué)依據(jù)。

    1 材料與方法

    1.1 試驗(yàn)材料

    試驗(yàn)于2018—2020年度在河南省商丘市農(nóng)林科學(xué)院小麥試驗(yàn)基地的低溫模擬室進(jìn)行。試驗(yàn)材料選取河南省主推的4個(gè)不同類型的冬小麥品種,弱春性品種偃展4110(YZ4110)和蘭考198(LK198),半冬性品種鄭麥366(ZM366)和豐德存麥21(FDC21)。試驗(yàn)用盆為直徑25 cm、高度35 cm、兩頭開口的圓柱形空心管。試驗(yàn)用土取自商丘市農(nóng)林科學(xué)院小麥試驗(yàn)基地大田土壤,經(jīng)過細(xì)篩、曬干、稱重后裝入桶中,土壤肥力狀況平均為有機(jī)質(zhì)13.4 g·kg-1,水解氮74.2 mg·g-1,速效磷36.31 mg·g-1,速效鉀129.48 mg·g-1。播種前將盆埋在大田之中,使其頂部與地表基本持平,使得盆栽生長環(huán)境等同于大田環(huán)境。分別于2018年11月9日、2019年11月12日進(jìn)行播種,每盆均勻定株7株。低溫處理時(shí)期為雌雄蕊分化期,并設(shè)定兩個(gè)水分處理,分別為凍前1周不灌水(土壤含水量13%—15%,土壤表層干旱,D),凍前1周灌水1次(土壤含水量20%—22%,土壤水飽滿,W)。在兩個(gè)水分條件下,并設(shè)置了5個(gè)溫度梯度T1、T2、T3、T4、T5,分別對(duì)應(yīng)-2℃、-4℃、-6℃、-8℃、-10℃,以及對(duì)照CK(當(dāng)日大田氣溫分別為11℃/5℃、22℃/8℃)。低溫處理時(shí)間為8 h,每次處理為5個(gè)重復(fù)。低溫處理結(jié)束后,立即將盆栽搬到大田原位,使其恢復(fù)原來的生長環(huán)境,并于低溫脅迫后的第2天進(jìn)行取樣測(cè)定生理指標(biāo)。

    1.2 低溫模擬室

    低溫室的長寬高分別為7 m×5 m×2.6 m,頂部裝置有2條尼龍布形成的通風(fēng)管道,分別與室外的兩臺(tái)制冷壓縮機(jī)連接。在管道兩側(cè)均勻排布若干個(gè)直徑為5 cm的圓形通風(fēng)口,冷風(fēng)通過圓形通風(fēng)口在1.9 m以上的水平方向上進(jìn)行冷氣流動(dòng),并形成一個(gè)制冷空間。冷空氣均勻下沉,形成強(qiáng)輻射降溫。

    1.3 生理指標(biāo)測(cè)定

    低溫脅迫后第2天,分別取各處理以及對(duì)照品種完全展開的頂1葉用于生理指標(biāo)測(cè)定,每個(gè)處理設(shè)3次重復(fù)。葉片相對(duì)含水量(LWC)、莖相對(duì)含水量(SWC)采用烘干法,色素類指標(biāo)含量測(cè)定采用分光光度法,具體包括葉綠素a(Chla)、葉綠素b(Chlb)、類胡蘿卜素(Car)、葉綠素總量(Chla+b)、葉綠素a/b(Chla/b),可溶性糖含量(SSC)采用蒽酮比色法[18],可溶性蛋白(BCA)采用考馬斯亮藍(lán)比色法[18],游離脯氨酸含量(Pro)采用磺基水楊酸浸提-酸性茚三酮顯色法[19],丙二醛(MDA)含量采用TBA比色法[20],超氧化物歧化酶(SOD)采用NBT比色法[21],過氧化物酶(POD)采用愈創(chuàng)木酚比色法[21],過氧化氫酶(CAT)采用過氧化氫法測(cè)定[21]。

    1.4 葉綠素?zé)晒鈪?shù)測(cè)定

    葉綠素?zé)晒鈪?shù)測(cè)定采用德國MINI-PAM量子分析儀。測(cè)定在上午進(jìn)行,均選擇完全展開的頂1葉,每個(gè)處理測(cè)試3張葉片。葉片經(jīng)過30 min暗適應(yīng)后測(cè)定初始熒光(o)、最大熒光(m)、PSⅡ潛在活性(v/o)和PSⅡ最大光學(xué)效率(v/m);之后測(cè)定光化學(xué)淬滅系數(shù)(qP)、相對(duì)電子傳遞速率(ETR)和PSII實(shí)際光量子效率(Y(II))。

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

    利用Excel軟件對(duì)原始數(shù)據(jù)進(jìn)行預(yù)處理,SPSS19.0軟件對(duì)預(yù)處理后數(shù)據(jù)進(jìn)行主成分、方差、聚類和逐步回歸分析。因?yàn)椴煌霓r(nóng)藝指標(biāo)有不同的數(shù)據(jù)范圍和單位,采用數(shù)據(jù)標(biāo)準(zhǔn)化進(jìn)行轉(zhuǎn)化,每個(gè)生理指標(biāo)轉(zhuǎn)換如下:

    式中,′為每個(gè)指標(biāo)標(biāo)準(zhǔn)化之后的數(shù)據(jù),為第個(gè)品種指標(biāo)測(cè)定值,和′為第個(gè)指標(biāo)的標(biāo)準(zhǔn)偏差和第個(gè)指標(biāo)的樣本平均值。

    在數(shù)據(jù)預(yù)處理后,用SPSS19.0進(jìn)行主成分分析(PCA),用公式(2)求得每個(gè)處理各綜合指標(biāo)的隸屬函數(shù)值:

    式中,表示第個(gè)綜合指標(biāo);表示第個(gè)綜合指標(biāo)的最小值;表示第個(gè)綜合指標(biāo)的最大值。

    用公式(3)可求出各綜合指標(biāo)的權(quán)重:

    式中,表示第個(gè)綜合指標(biāo)在所有綜合指標(biāo)中的重要程度即權(quán)重;為各小麥品種第個(gè)綜合指標(biāo)的貢獻(xiàn)率。

    用公式(4)計(jì)算每個(gè)處理下的小麥的凍害程度的大小[22],以值表示:

    式中,值為各低溫處理的小麥通過用綜合指標(biāo)評(píng)價(jià)得到的凍害程度綜合評(píng)價(jià)值。

    將值作為因變量,對(duì)各個(gè)標(biāo)準(zhǔn)化的數(shù)據(jù)進(jìn)行逐步回歸分析,獲得回歸方程。

    2 結(jié)果

    2.1 凍害對(duì)不同品種和水分處理的幼穗形態(tài)及產(chǎn)量的影響

    小麥幼穗受凍的形態(tài)結(jié)果(以YZ4110為例)表明,隨著溫度不斷降低,幼穗受凍程度逐漸嚴(yán)重,幼穗逐漸失水萎縮變形,顏色變白,部分幼穗死亡,小花敗育加重,導(dǎo)致穗粒數(shù)顯著減少。在不灌水(D)條件下,T3溫度時(shí)出現(xiàn)幼穗部分死亡,當(dāng)溫度降至T4和T5時(shí),小麥幼穗死亡且不可恢復(fù);而在灌水處理(W)條件下,T3溫度脅迫的幼穗凍死率較低,受凍程度明顯輕于不灌水處理,當(dāng)溫度降至T4時(shí),幼穗才出現(xiàn)部分死亡,至T5時(shí),幼穗死亡(圖1)。

    由圖2可知,凍害條件下2年產(chǎn)量的變化趨勢(shì)相同,不同低溫導(dǎo)致的小麥產(chǎn)量損失情況因水分和品種而表現(xiàn)不同。當(dāng)溫度降至T1時(shí),不灌水(D)的弱春性品種(YZ4110和LK198)的產(chǎn)量降低,較CK差異達(dá)顯著水平,而半冬性品種(ZM366和FDC21)產(chǎn)量的降低未達(dá)到顯著水平;當(dāng)溫度降至T2及以下時(shí),兩種水分條件下的所有品種產(chǎn)量均顯著降低。從產(chǎn)量損失率均值來看,不灌水條件下 T1、T2、T3、T4—T5的產(chǎn)量損失率分別為3.33%—7.83%、17.78%—28.46%、33.56%—49.38%和67.27%—93.65%;灌水條件下的產(chǎn)量損失率分別為1.70%—4.13%、12.15%—18.60%、22.33%—32.75%和41.92%—76.04%。兩種水分條件下產(chǎn)量損失率在品種間均表現(xiàn)為LK198>YZ4110>ZM366>FDC21。

    D為不灌水處理;W為灌水處理。5個(gè)低溫處理分別為-2℃(T1),-4℃(T2),-6℃(T3),-8℃(T4)和-10℃(T5)。下同

    A為2019年,B為2020年。YZ4110、LK198、ZM366和FDC21分別為偃展4110、蘭考198、鄭麥366和豐德存麥21。誤差線上不同字母表示處理間差異顯著(P<0.05)。下同

    2.2 不同處理冬小麥生理生化指標(biāo)的變異性

    表1為各個(gè)指標(biāo)在不同溫度和水分條件下的變化。從變異系數(shù)可以看出,在不灌水條件下,因低溫和品種引起的各指標(biāo)變異系數(shù)范圍為4.87%—52.12%,而灌水條件下為3.74%—44.59%,且各指標(biāo)在不灌水的變異系數(shù)均大于灌水,這表明了低溫脅迫前進(jìn)行灌水,使得生理指標(biāo)的穩(wěn)定性增強(qiáng)。從平均值可以得出,灌水條件下滲透調(diào)節(jié)物質(zhì)含量及抗氧化物酶活性均提高。灌水下的滲透調(diào)節(jié)物質(zhì)(可溶性糖、可溶性蛋白、脯氨酸)含量較不灌水分別提高了26.80%、10.23%和13.43%;抗氧化物酶活性(CAT、POD、SOD)較不灌水增幅分別為14.06%、7.2%和12.81%。對(duì)細(xì)胞膜具有毒害作用的MDA含量在灌水條件下較不灌水降低了13.06%。而葉綠素指標(biāo)、植株含水量、葉綠素?zé)晒鈪?shù)等也是反映植株受凍程度的生理表征指標(biāo),其數(shù)值在灌水條件下均高于不灌水。

    表1 不同水分處理下低溫脅迫處理的植株各性狀比較

    Chla為葉綠素a,Chlb為葉綠素b,Car為類胡蘿卜素,Chla+b為葉綠素總量,Cha/b為葉綠素a/b,SWC為莖含水量,LWC為葉含水量,SSC為可溶性糖含量,Pro為脯氨酸,Y(II)為PSII實(shí)際光量子效率,ETR為相對(duì)電子傳遞速率,qP為光化學(xué)淬滅系數(shù),v/m為PSⅡ最大光學(xué)效率,v/o為PSⅡ潛在活性,SP為可溶性蛋白,CAT為過氧化氫酶,MDA為丙二醛,POD為過氧化物酶,SOD為超氧化物歧化酶。下同

    Chla was chlorophyll a, Chlb was chlorophyll b, Car was Carotenoids, Chla+b was total content of chlorophyll, Chla/b was chlorophyII a/b, SWC was stem water content, LWC was leaf water content, SSC was soluble sugar content, Pro was proline, Y(II) was PSII actual photochemical efficiency, ETR was relative electron transport rate, qP was photochemical quenching,v/mwas PSII maximum photochemical efficiency,v/owas PSII potential activity, SP was soluble protein, CAT was catalase, MDA was malondialdehyde, POD was peroxide dismutase, SOD was superoxide dismutase. The same as below

    通過Mann-Whitney U非參數(shù)檢驗(yàn)可知,除Car、Chla+b、SWC、Y (II)、ETR和POD外,其余生理指標(biāo)均在灌水和不灌水間存在顯著差異(<0.05),其中v/m達(dá)到極顯著水平(<0.01)。以上結(jié)果表明灌水增強(qiáng)了植株生理生化指標(biāo)對(duì)低溫的應(yīng)激響應(yīng)能力,減緩了低溫對(duì)植株造成的傷害程度,進(jìn)而提高抗寒性。

    2.3 不同生理生化指標(biāo)之間相關(guān)性分析

    測(cè)定指標(biāo)中,MDA與其他指標(biāo)間呈負(fù)相關(guān),且均達(dá)極顯著水平(<0.01),其他指標(biāo)間相關(guān)性均為正相關(guān)(表2)。Chlb與余下參數(shù)中一半指標(biāo)(9個(gè))的相關(guān)性未達(dá)顯著水平;Chla/b與SWC和BAC間相關(guān)性不顯著;SWC與5個(gè)指標(biāo)(SSC、BCA、CAT、POD和SOD)間相關(guān)性不顯著;除以上指標(biāo)外,其余各生理指標(biāo)間的相關(guān)性均達(dá)到顯著水平(<0.05)。各生理指標(biāo)之間表現(xiàn)出大小不等的相關(guān)性,且大多數(shù)指標(biāo)之間的相關(guān)性達(dá)到顯著性水平,這就導(dǎo)致了生理指標(biāo)之間反映植株的抗寒信息存在交叉,甚至發(fā)生很大程度的重疊。不同低溫下各個(gè)指標(biāo)在植株抗寒性應(yīng)激反應(yīng)中的功能不盡相同,直接應(yīng)用各個(gè)指標(biāo)對(duì)植株的抗寒性進(jìn)行評(píng)價(jià)是不準(zhǔn)確的,很難準(zhǔn)確反映出小麥的凍害嚴(yán)重程度,導(dǎo)致評(píng)價(jià)結(jié)果具有很大的局限性。為了彌補(bǔ)直接應(yīng)用各個(gè)生理指標(biāo)進(jìn)行抗寒性評(píng)價(jià)的不足,有必要進(jìn)一步采用主成分和聚類分析進(jìn)行綜合評(píng)價(jià)。

    2.4 主成分分析及隸屬函數(shù)分析

    將原始數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化,并使用SPSS軟件對(duì)標(biāo)準(zhǔn)化后的數(shù)據(jù)主成分分析。根據(jù)累計(jì)貢獻(xiàn)率>85%的原則,共提取了6個(gè)主成分,并將19個(gè)單項(xiàng)指標(biāo)轉(zhuǎn)化為6個(gè)綜合主成分(CI1—CI6),各主成分的特征向量及貢獻(xiàn)率如表3所示。前6個(gè)主成分的累計(jì)貢獻(xiàn)率高達(dá)87.459%,具有較強(qiáng)的代表性,包含了原始指標(biāo)所攜帶的大部分信息。其中,第一主成分中指標(biāo)系數(shù)較高的為Chla、Car、Chla/b和熒光參數(shù),第二主成分的指標(biāo)為CAT和SOD,第三主成分的指標(biāo)為Chlb和Chla+b,第四主成分的指標(biāo)為可溶性蛋白BCA和SWC??梢?,光合性狀、抗氧化物酶、滲透調(diào)節(jié)物質(zhì)是反映植物受凍情況最顯著的生理指標(biāo)。

    表2 不同生理生化指標(biāo)間的相關(guān)系數(shù)

    *和**分別表示0.05和0.01水平顯著相關(guān)

    * and ** are significant correlation at 0.05 and 0.01 probability levels, respectively

    表3 各綜合指標(biāo)特征值、貢獻(xiàn)率以及隸屬函數(shù)權(quán)重

    同時(shí)利用公式(2)對(duì)各低溫處理下的冬小麥綜合指標(biāo)進(jìn)行隸屬函數(shù)分析并獲得各成分的隸屬函數(shù)值U(X)及其對(duì)應(yīng)權(quán)重。并利用公式(4)計(jì)算各個(gè)低溫處理下的冬小麥凍害程度綜合評(píng)價(jià)值(值)(圖3)。由各品種的值范圍結(jié)果可知,同一水分條件下,半冬性品種的最大、最小及平均值均高于弱春性品種;而同一品種條件下,灌水的各項(xiàng)值也均高于不灌水。

    圖3 不同灌水條件下4個(gè)小麥品種凍害程度綜合評(píng)價(jià)值(F值)的變化

    為篩選出表征小麥凍害程度的重要鑒定指標(biāo),以上述的綜合評(píng)價(jià)值()作因變量,試驗(yàn)測(cè)定的19個(gè)生理生化指標(biāo)(X)為自變量,采用逐步回歸分析的方法,篩選出了7個(gè)關(guān)鍵指標(biāo),并建立了最優(yōu)化的評(píng)價(jià)凍害程度的定量數(shù)學(xué)模型:

    =0.605+0.038X1+0.031X7+0.022X9+0.06X13+0.022X15-0.025X17+0.038X19(2=0.990,=0.0001) 式中,X1、X7、X9、X13、X15、X17、X19分別為葉綠素a、葉片含水量、脯氨酸、v/m、可溶性蛋白、丙二醛和超氧化物歧化酶。

    2.5 不同處理的F預(yù)測(cè)值聚類分析

    采用最近距離法對(duì)預(yù)測(cè)值進(jìn)行聚類分析(圖4),可將48個(gè)處理的小麥凍害程度劃分為5類,第一類為未受凍(D-CK、W-CK),第二類為輕度受凍(D-T1、W-T1),第三類為中度受凍(D-T2、W-T2、W-T3),第四類為重度受凍(D-T3、W-T4),第五類為特重受凍(D-T4、W-T5、D-T5)。不同水分條件影響聚類結(jié)果,在T1和T2溫度下,水分之間無差異;當(dāng)溫度降至T3時(shí),D-T3處理為重度受凍,而W-T3處理則為中度受凍;溫度繼續(xù)降至T4時(shí),W處理處于重度受凍,而D處理則為特重受凍。輕度及以下受凍的產(chǎn)量損失率低于10%,中度受凍為10%—30%,重度受凍為30%—50%,特重受凍為50%以上。

    為驗(yàn)證值評(píng)估凍害程度的準(zhǔn)確性,進(jìn)一步分析產(chǎn)量損失率與相對(duì)值之間的關(guān)系(圖5)。產(chǎn)量損失率與值之間呈一元三次擬合關(guān)系(2=0.898),當(dāng)產(chǎn)量損失率低于10%,該階段處于輕度及以下受凍程度,主要表現(xiàn)為內(nèi)在生理指標(biāo)的應(yīng)激及響應(yīng)能力提高,以減少凍害的損傷程度,使得值趨近于1;當(dāng)產(chǎn)量損失率大于10%,此時(shí)凍害已達(dá)到中度程度,已超出生理應(yīng)激范圍,值下降趨勢(shì)逐漸明顯;當(dāng)脅迫程度處于重度以上時(shí),凍害對(duì)產(chǎn)量影響嚴(yán)重,各品種值降至最低,降幅平緩。這表明了基于植株早期生理生化變化所構(gòu)建的凍害綜合模型可以準(zhǔn)確地評(píng)價(jià)小麥后期的產(chǎn)量損失狀況。

    圖4 不同水分、品種、溫度處理下的F值聚類

    圖5 相對(duì)F值與產(chǎn)量損失率間關(guān)系

    2.6 低溫脅迫對(duì)不同處理冬小麥生理指標(biāo)的影響

    通過主成分分析及多元逐步分析篩選出7個(gè)關(guān)鍵指標(biāo),以此7個(gè)生理指標(biāo)為例,具體分析低溫脅迫的影響。由圖6可知,隨著凍害程度的加重,不同小麥品種的葉片含水量、葉綠素a、v/m均下降。同一低溫脅迫下,葉片含水量、葉綠素a、v/m的變化情況因品種和灌水而降幅不同。整體而言,T1和T2溫度時(shí)葉片含水量、葉綠素a、v/m降幅較小,當(dāng)溫度降至T3時(shí),降幅增大且不灌水和灌水間差異顯著,兩個(gè)水分條件下降幅最大和最小的品種均為LK198和FDC21。T3溫度相較于CK,葉片含水量在不灌水條件下降幅為4.75%—11.85%,而在灌水條件下降幅為3.49%—4.64%;葉綠素a含量在不灌水和灌水條件下降幅分別為17.48%—31.32%和8.46%—15.27%;v/m在不灌水和灌水條件下降幅分別為14.46%—24.07%和7.83%—15.95%。

    由圖7可知,在正常溫度CK條件下,脯氨酸和可溶性蛋白在灌水和不灌水條件下存在顯著差異,這表明不灌水下植株處于相對(duì)脅迫狀態(tài),使得滲透調(diào)節(jié)物質(zhì)含量提高。隨著脅迫的不斷加重,滲透調(diào)節(jié)物質(zhì)脯氨酸和可溶性蛋白在兩種水分條件下均表現(xiàn)出先升后降的趨勢(shì)。脯氨酸含量在不灌水條件下各品種均于T2溫度時(shí)達(dá)到最大值,且處理間差異達(dá)到顯著,其中LK198增幅最小,F(xiàn)DC21增幅最大,分別較對(duì)照增加了15.43%和27.1%,當(dāng)溫度降至T3時(shí),ZM366和FDC21的脯氨酸含量仍高于CK,而另外兩個(gè)品種已低于CK;灌水條件下脯氨酸含量出現(xiàn)最高值的溫度因品種而異,YZ4110和LK198在T2溫度時(shí)脯氨酸含量最高,增幅分別為39.94%和37.38%,而ZM366和FDC21則在T3溫度時(shí)最高,分別較CK增加了51.16%和57.86%,且4個(gè)品種均在T5溫度時(shí)脯氨酸含量開始低于CK。不灌水下的可溶性蛋白含量在T1(YZ4110和LK198)或T2(ZM366和FDC21)溫度時(shí)達(dá)最大值,其中YZ4110和LK198的最大增幅為12.52%和7.17%,而ZM366和FDC21在T2下可溶性蛋白含量較CK分別增加了31.9%和32.5%。而灌水條件下的各品種可溶性蛋白含量則均在T3溫度時(shí)達(dá)峰值,較CK增幅最大的為FDC21(70.31%),增幅最小為LK198(42.67%)。

    由圖8可知,在正常溫度CK條件下,兩種水分處理間SOD酶活性存在顯著性差異。SOD酶作為抗氧化物酶系統(tǒng)的第一道防線,隨著溫度的不斷降低,SOD酶活性呈現(xiàn)出先升后降的趨勢(shì)。不灌水條件下的各品種SOD酶活性在T2溫度均達(dá)最大值,4個(gè)品種的最大值較CK分別增加了25.39%、19.65%、29.83%和33.04%;灌水條件下各品種SOD酶活性均在T3溫度下達(dá)最大值,以上4個(gè)品種較CK增幅分別為55.84%、49.71%、62.6%和68.62%。作為膜脂過氧化最終產(chǎn)物的MDA,其含量隨溫度降低呈先降后升的趨勢(shì)。在不灌水條件下,4個(gè)品種均在T1溫度時(shí)最低,降幅分別為8.67%(YZ4110)、6.05%(LK198)、10%(ZM366)和11.53%(FDC21);而灌水條件下的MDA含量在T2溫度時(shí)最小,4個(gè)品種一致,具體降幅為FDC21(29.54%)>ZM366(29.13%)>YZ4110(25.65%)>LK198(17.86%)。

    圖6 不同灌水處理下凍害對(duì)葉片含水量、葉綠素a和Fv/Fm的影響

    圖7 不同灌水處理下凍害對(duì)脯氨酸和可溶性蛋白含量的影響

    圖8 不同灌水處理下凍害對(duì)SOD酶活性及MDA含量的影響

    3 討論

    3.1 灌水緩解低溫脅迫影響的生理機(jī)制

    低溫脅迫下植物會(huì)發(fā)生一系列生理生化反應(yīng)以響應(yīng)和抵御低溫造成的傷害[23],凍害程度超出植株承受能力后,抵御低溫脅迫的能力降低,如滲透調(diào)節(jié)物質(zhì)含量、抗氧化物酶活性降低,使得細(xì)胞喪失抗寒能力,從而導(dǎo)致植株死亡[24-28]。為了減緩植株受凍害的程度,前人通過灌水措施加強(qiáng)作物田間保溫,提高植株冠層溫度,有助于穩(wěn)定葉片含水量[29],在減緩葉綠素和光合機(jī)構(gòu)受損程度方面發(fā)揮一定作用[30],增強(qiáng)植株的抗寒性,為此有研究者將植株冠層溫度作為低溫脅迫指標(biāo)[31]。曹娜等[32]研究表明,晚稻抽穗揚(yáng)花期在受到低溫脅迫時(shí),通過灌水處理可以提高冠層溫度和葉片溫度,降低冷害脅迫程度。劉立軍等[29]認(rèn)為,低溫條件下,灌水處理的高羊茅葉片含水量降幅較小,不僅提高滲透調(diào)節(jié)物質(zhì)含量,而且降低膜脂過氧化程度,緩解電解質(zhì)滲透速度,從而減弱低溫脅迫對(duì)細(xì)胞膜的傷害。單旭東等[33]研究認(rèn)為,灌水處理通過提高抗氧化物酶活性,從而降低了MDA含量,緩解低溫脅迫對(duì)結(jié)縷草細(xì)胞膜的損傷。曹倩等[34]發(fā)現(xiàn)苜蓿在遭受低溫脅迫時(shí),適宜土壤水分提高了滲透調(diào)節(jié)物質(zhì)的含量,增強(qiáng)了植株的滲透調(diào)節(jié)能力,使得光合機(jī)構(gòu)和葉綠素體結(jié)構(gòu)破壞程度降低。在小麥凍前澆水研究表明,凍前灌水對(duì)遭受到低溫脅迫的冬小麥具有正效應(yīng),可有效緩解凍害造成的傷害[4]。但前人側(cè)重研究土壤表層濕度減輕低溫導(dǎo)致的產(chǎn)量損失程度,較少解釋增加土壤濕度進(jìn)而增強(qiáng)植株抗寒性的生理機(jī)制。植株含水量是對(duì)土壤水分狀況的最直接的反映,本試驗(yàn)中,灌溉土壤含水量高于未灌溉土壤含水量,進(jìn)而導(dǎo)致正常溫度CK條件下灌溉植株含水量高于未灌溉植株,但兩種水分處理間差異并不顯著,而可溶性蛋白、脯氨酸、SOD酶活性在水分處理間差異達(dá)到顯著水平,這在一定程度上顯示出一些生理生化指標(biāo)的變化盡管與含水量有關(guān),卻較含水量更能反映出在低溫脅迫時(shí)緩解凍害的作用,其內(nèi)在的機(jī)理還有待進(jìn)一步明確。在同一溫度下,灌水處理的葉片含水量高于不灌水處理,在正常溫度CK條件下兩種水分處理間沒有顯著差異,但隨著低溫脅迫的加重,兩種水分處理間差異水平越明顯,如T2以后差異達(dá)顯著水平,導(dǎo)致這種現(xiàn)象的原因可能與水分相態(tài)發(fā)生變化進(jìn)而影響組織結(jié)構(gòu)及水分運(yùn)輸有關(guān),在這方面還需深入研究明確。本研究中葉片含水量與MDA含量及熒光參數(shù)極顯著相關(guān),表明良好的水分狀況有利于減輕膜脂過氧化的傷害程度,維持較高的光合活性,進(jìn)而減輕低溫對(duì)植株的凍害。整體而言,灌溉處理下葉片含水量高于未灌水,表明適當(dāng)?shù)膬銮肮嗨幚碓鰪?qiáng)了細(xì)胞保水能力,并降低了細(xì)胞膜透性,同時(shí)抗氧化物酶類活性提高,緩解了膜脂過氧化程度,使得MDA含量低于不灌水處理;可溶性蛋白、脯氨酸等滲透調(diào)節(jié)物質(zhì)的累積量增加,進(jìn)一步降低了細(xì)胞冰點(diǎn),減緩了細(xì)胞膜受損程度,并且緩解了葉綠素的降解,光合受損程度降低,進(jìn)而增強(qiáng)了植株的抗寒性。各品種在灌水處理下的葉片含水量、葉綠素含量、葉綠素?zé)晒鈪?shù)等指標(biāo)在不同低溫脅迫下的降幅均低于不灌水處理;而通常作為抗寒性鑒定指標(biāo)SOD酶等抗氧化物酶、滲透調(diào)節(jié)物質(zhì)含量在相同低溫脅迫下,灌水高于不灌水,這表明灌水處理的小麥植株通過最大程度的生理應(yīng)激響應(yīng),以維持機(jī)體的正常代謝機(jī)能,減緩了受凍程度,并使得各品種能夠在嚴(yán)重脅迫下生存。本試驗(yàn)中,灌水處理的小麥在T3溫度下凍害程度僅為中度受凍,而不灌水則已處于重度受凍。這表明凍前灌水提高了冬小麥的抗寒性,緩解了低溫脅迫造成的凍害程度,生產(chǎn)中應(yīng)依據(jù)天氣預(yù)報(bào)在寒流到來前,尤其對(duì)干旱地塊進(jìn)行灌水,可以在很大程度上防御晚霜凍害。

    3.2 多因素作用下凍害程度的綜合評(píng)估

    不同的生理生化指標(biāo)對(duì)逆境脅迫的響應(yīng)程度不同,單項(xiàng)生理生化指標(biāo)無法準(zhǔn)確評(píng)估植株抗逆性[9]。前人研究多利用綜合評(píng)價(jià)方法鑒定植株的抗逆性[22,35-36]。姜麗娜等[17]分析了24個(gè)小麥品種受到拔節(jié)期低溫脅迫的生理響應(yīng),通過隸屬函數(shù)和綜合排序法將品種抗寒性分為5類。張軍等[37]通過分析低溫脅迫對(duì)生理指標(biāo)的影響效應(yīng),基于所測(cè)生理指標(biāo)的隸屬值評(píng)估小麥品種的抗寒性。但由于各指標(biāo)之間存在一定的相關(guān)性和生理重疊性,僅用隸屬函數(shù)法評(píng)估小麥的抗寒性具有較大的局限性[10]。而主成分分析法(PCA)可以將多個(gè)指標(biāo)轉(zhuǎn)化成新的個(gè)數(shù)較少且彼此獨(dú)立的綜合指標(biāo),同時(shí)根據(jù)各自貢獻(xiàn)率的大小明確每個(gè)綜合指標(biāo)的相對(duì)重要性。因此,前人利用主成分和隸屬函數(shù)相結(jié)合的方法評(píng)估植物的耐寒性,實(shí)現(xiàn)了綜合評(píng)價(jià)不同品種抗寒性強(qiáng)弱的目標(biāo)[16-17,38-39]。前人對(duì)小麥凍害程度的研究多為單一低溫條件下品種資源的抗寒性比較,而不同的水分條件及基因型品種對(duì)低溫的響應(yīng)程度及抗寒性存在明顯差異,這導(dǎo)致以前確立的綜合評(píng)價(jià)方法并不適用于不同品種、灌水及溫度等多因素綜合條件下的小麥凍害程度評(píng)估。為此,本試驗(yàn)研究小麥生產(chǎn)中3個(gè)關(guān)鍵要素(2個(gè)品種類型、2種灌水處理和6個(gè)溫度處理)條件下植株的受凍程度,但僅用單一的生理指標(biāo)或者簡(jiǎn)單相關(guān)分析凍害程度具有局限性,且缺乏準(zhǔn)確性。故利用主成分分析和隸屬函數(shù)法將19個(gè)生理指標(biāo)和熒光參數(shù)進(jìn)行綜合分析,最終得到凍害程度綜合評(píng)價(jià)值(值)。隨著溫度脅迫的不斷加重,不同品種在相同水分下的綜合值均表現(xiàn)為先升后降趨勢(shì),其原因是脅迫程度未超出細(xì)胞承受能力時(shí),抗氧化物酶和滲透調(diào)節(jié)物質(zhì)含量升高,以抵御低溫脅迫對(duì)細(xì)胞的傷害,因此值越大,表示抗寒性越強(qiáng);當(dāng)脅迫加重后,由于細(xì)胞膜破裂,抗氧化物酶和滲透調(diào)節(jié)物質(zhì)含量開始明顯下降,葉片含水量、葉綠素含量和葉綠素?zé)晒鈪?shù)等降低,進(jìn)而值表現(xiàn)逐漸變小。在相同的溫度和水分條件時(shí),值表現(xiàn)為半冬性品種(FDC21、ZM366)>弱春性品種(YZ4110、LK198),表明基因型在作物抗寒機(jī)制中發(fā)揮基礎(chǔ)作用;在相同的品種和溫度條件下,值表現(xiàn)為灌水>不灌水,表明灌水處理通過加強(qiáng)對(duì)生理指標(biāo)的應(yīng)激響應(yīng)能力,在生理生化層面上增強(qiáng)了植株的抗寒能力。進(jìn)一步以值為因變量,19個(gè)生理指標(biāo)為自變量,通過逐步回歸確立了凍害綜合指數(shù)的定量估算模型。基于該預(yù)測(cè)值進(jìn)行聚類分析,可劃分為5類:未受凍、輕度受凍、中度受凍、重度受凍和特重受凍,對(duì)應(yīng)的產(chǎn)量損失率分別為0、<10%、10%—30%、30%—50%及>50%。研究結(jié)果有利于準(zhǔn)確評(píng)估凍害發(fā)生后的小麥?zhǔn)軆龀潭燃爱a(chǎn)量損失狀況,為凍害后的補(bǔ)救措施制定提供科學(xué)依據(jù)。

    3.3 凍害綜合評(píng)估的指標(biāo)篩選與提前預(yù)測(cè)評(píng)估

    有關(guān)植物抗寒性的生理生化鑒定指標(biāo),前人相關(guān)研究并不完全一致,這與作物種類、品種類型、生態(tài)區(qū)以及生長環(huán)境有關(guān)[16,25,27,29]。作物抗逆性是基于基因和環(huán)境互作的多種生理指標(biāo)的綜合體現(xiàn),但由于生理指標(biāo)較多,各個(gè)指標(biāo)在逆境應(yīng)激反應(yīng)中的作用也不一樣,并不是所有的生理生化指標(biāo)均能有效準(zhǔn)確鑒定植物的抗逆性能。為此,前人多通過主成分-回歸分析相結(jié)合的方法篩選出抗逆性評(píng)價(jià)指標(biāo),先后在花生耐蔭性[40]、谷子耐鹽性[41]及棉花抗旱性評(píng)估[34]方面取得較好效果。目前對(duì)于小麥抗寒性指標(biāo)的篩選研究大多通過生理指標(biāo)與農(nóng)藝參數(shù)的相關(guān)系數(shù)及關(guān)聯(lián)性進(jìn)行,王濤等[42]以小麥返青率為目標(biāo),篩選出相對(duì)含水量是小麥抗寒性評(píng)價(jià)指標(biāo);趙瑞玲等[43]以小麥葉片枯黃面積為標(biāo)準(zhǔn),發(fā)現(xiàn)相對(duì)電導(dǎo)率可指示凍害程度;張自陽等[44]以小麥幼穗結(jié)實(shí)數(shù)為對(duì)象,明確了POD活性、SOD活性、MDA含量為小麥抗寒性鑒定指標(biāo);吳青霞等[45]以產(chǎn)量為標(biāo)準(zhǔn),篩選出脯氨酸、可溶性蛋白含量等為抗寒性指標(biāo)。前人基于主成分-隸屬分析-逐步回歸等多元統(tǒng)計(jì)方法篩選生理指標(biāo)的研究較少,缺少低溫脅迫下植株生理生化反應(yīng)的綜合評(píng)價(jià)指標(biāo)構(gòu)建,以及少數(shù)關(guān)鍵敏感性因子的降維。隨著全球氣候變暖及春季極端溫度波動(dòng)加大,晚霜凍害的發(fā)生頻率呈上升趨勢(shì),晚霜凍害已成為影響黃淮麥區(qū)小麥生產(chǎn)的主要災(zāi)害因子。晚霜凍害嚴(yán)重時(shí),會(huì)使冬小麥小穗大面積受傷甚至死亡,降低產(chǎn)量。由于晚霜凍害發(fā)生較晚,往往莖葉無異常表現(xiàn),受害部位主要是穗部,危害情況難以直觀地了解和評(píng)估,只有在抽穗以后才能夠觀察,而此時(shí)小麥成穗數(shù)已完全確定,對(duì)穗粒數(shù)所造成的影響也難以挽回。因此,凍害發(fā)生后及時(shí)依據(jù)少數(shù)關(guān)鍵幾個(gè)植株生理生化指標(biāo)的變化定量評(píng)估小麥凍害程度,這對(duì)于產(chǎn)量損失評(píng)估及補(bǔ)救措施制定具有重要意義。本研究設(shè)置不同水分條件、品種類型及低溫脅迫水平,于凍后第2天測(cè)定19個(gè)生理生化指標(biāo),通過主成分-隸屬函數(shù)-逐步回歸分析法篩選出7個(gè)關(guān)鍵性評(píng)估凍害程度的生理指標(biāo):葉綠素a、葉片含水量、脯氨酸、v/m、可溶性蛋白、MDA和SOD,以此確立的凍害綜合指數(shù)預(yù)測(cè)值與產(chǎn)量間方程決定系數(shù)達(dá)到0.898。小麥晚霜凍害發(fā)生后,植株生理生化的變化比外觀形態(tài)指標(biāo)更為顯著,在凍后第2天及時(shí)進(jìn)行生理指標(biāo)的采集,更能及時(shí)反映出植株的凍害程度,這為早期評(píng)估晚霜凍害對(duì)小麥造成的損害程度,從而為農(nóng)業(yè)生產(chǎn)中及時(shí)制定積極主動(dòng)的補(bǔ)救措施提供依據(jù)和決策參考。此外,本文篩選出的7個(gè)指標(biāo)不僅能代表眾多生理指標(biāo)的原始信息量,更達(dá)到了評(píng)價(jià)指標(biāo)降維的目的,減少了實(shí)際測(cè)量的工作量,這有利于準(zhǔn)確有效評(píng)估小麥凍害程度,可應(yīng)用性強(qiáng)。但是生產(chǎn)中晚霜凍害發(fā)生的時(shí)期及小麥幼穗發(fā)育所處的階段并不固定,抗寒性生理指標(biāo)是否存在差異,還需進(jìn)一步驗(yàn)證。小麥拔節(jié)后幼穗發(fā)育對(duì)低溫的抵抗力顯著下降,及時(shí)灌水改善土壤墑情是減輕晚霜凍害的主要措施,有關(guān)灌水增強(qiáng)小麥低溫應(yīng)激響應(yīng)的護(hù)穗機(jī)理,以及減輕凍害傷害的代謝調(diào)控網(wǎng)絡(luò)均有待進(jìn)一步明確,這將為小麥凍害防御與補(bǔ)救減損提供理論依據(jù)和技術(shù)參考。

    4 結(jié)論

    本研究明確了兩種類型小麥品種在不同水分和低溫處理下植株的生理生化指標(biāo)變化,利用主成分分析、隸屬函數(shù)法及多元逐步回歸構(gòu)建出凍害綜合評(píng)估模型,并篩選出葉綠素a、葉片含水量、脯氨酸、v/m、可溶性蛋白、MDA和SOD為小麥凍害程度的評(píng)價(jià)指標(biāo)。進(jìn)一步利用聚類分析將不同處理的值分為5類:第一類為未受凍(D-CK、W-CK),第二類為輕度受凍(D-T1、W-T1),第三類為中度受凍(D-T2、W-T2、W-T3),第四類為重度受凍(D-T3、W-T4),第五類為特重受凍(D-T4、W-T5、D-T5),其產(chǎn)量損失率分別對(duì)應(yīng)為0、<10%、10%—30%、30%—50%及50%以上。比較而言,弱春性品種的凍害程度重于半冬性品種,灌水處理低于不灌水處理??梢姡谛←溡装l(fā)生凍害地區(qū)應(yīng)選用半冬性品種,且根據(jù)天氣變化及早灌水,有利于減輕晚霜凍害,并依靠?jī)龊﹁b定指標(biāo)及模型為災(zāi)害評(píng)估及管理提供技術(shù)支持。

    [1] Kodra E, Steinhaeuser K, Ganguly A R. Persisting cold extremes under 21st-century warming scenarios. Geophysical Research Letters, 2011, 38(8): 1-5.

    [2] Zheng B Y, Chenu K, Fernanda D M, Chapman S C. Breeding for the future: What are the potential impacts of future frost and heat events on sowing and flowering time requirements for Australian bread wheat () varieties?. Global Change Biology, 2012, 18(9): 2899-2914.

    [3] 張永閣. 柘城縣小麥凍害發(fā)生原因及防治措施. 河南農(nóng)業(yè), 2016, 1: 45.

    ZHANG Y G. Reasons and protective measures of wheat frost damage in Zhecheng. Agriculture of Henan, 2016, 1: 45. (in Chinese)

    [4] Bali A S, Sidhu G P S. Abiotic stress-induced oxidative stress in wheat//HASANNZZAMAN M, NAHAR K, HOSSAIN M A. Wheat Production in Changing Environments. Singapore: Springer, 2019: 225-239

    [5] Upadhyaya H, Sahoo L, Panda S K. Molecular physiology of osmotic stress in plants//ROUT G R, DAS A B. Molecular Stress Physiology of Plants. India: Springer, 2013: 179-192.

    [6] Galiba G, Vanková R, Tari I, BáNFALVI Z, MSC P P, DOBREV P, BOLDIZSáR A, VáGúJFALVI A, Kocsy G. Hormones, NO, antioxidant and metabolites as key players in plant cold acclimation//IMAI R, YOSHIDA M, MASTUMOTO N. Plant and Microbe Adaptations to Cold in a Changing World. New York: Springer, 2013: 73-87.

    [7] 陳思思, 李春燕, 楊景, 徐雯, 朱新開, 郭文善, 封超年. 拔節(jié)期低溫凍害對(duì)揚(yáng)麥16光合特性及產(chǎn)量形成的影響. 揚(yáng)州大學(xué)學(xué)報(bào)(農(nóng)業(yè)與生命科學(xué)版), 2014, 35(3) : 59-64.

    CHEN S S, LI C Y, YANG J, XU W, ZHU X K, GUO W S, FENG C N. Effect of low temperature at jointing stage on photosynthetic characteristics and yield in wheat cultivar Yangmai 16. Journal of Yangzhou University(Agricultural and Life Science Edition), 2014, 35(3): 59-64. (in Chinese)

    [8] 劉蕾蕾, 紀(jì)洪亭, 劉兵, 馬吉鋒, 肖瀏駿, 湯亮, 曹衛(wèi)星, 朱艷. 拔節(jié)期和孕穗期低溫處理對(duì)小麥葉片光合及葉綠素?zé)晒馓匦缘挠绊? 中國農(nóng)業(yè)科學(xué), 2018, 51(23) : 4434-4448.

    LIU L L, JI H T, LIU B, MA J F, XIAO L J, TANG L, CAO W X, ZHU Y. Effects of jointing and booting low temperature treatments on photosynthetic and chlorophyll fluorescence characteristics in wheat leaf. Scientia Agricultura Sinica, 2018, 51(23): 4434-4448. (in Chinese)

    [9] 王樹剛, 王振林, 王平, 王海偉, 李府, 黃瑋, 武玉國, 尹燕枰. 不同小麥品種對(duì)低溫脅迫的反應(yīng)及抗凍性評(píng)價(jià). 生態(tài)學(xué)報(bào), 2011, 31(4): 1064-1072.

    WANG S G, WANG Z L, WANG P, WANG H W, LI F, HUANG W, WU Y G, YIN Y P. Evaluation of wheat freezing resistance based on the responses of the physiological indices to low temperature stress. Acta Ecologica Sinica, 2011, 31(4): 1064-1072. (in Chinese)

    [10] 張笑笑, 潘映紅, 任富莉, 蒲偉軍, 王道平, 李玉斌, 陸平, 李桂英, 朱莉. 基于多重表型分析的準(zhǔn)確評(píng)價(jià)高粱抗旱性方法的建立. 作物學(xué)報(bào), 2019, 45(11) : 1735-1745.

    ZHANG X X, PAN Y H, REN F L, PU W J, WANG D P, LI Y B, LU P, LI G Y, ZHU L. Establishment of an accurate evaluation method for drought resistance based on multilevel phenotype analysis in sorghum. Acta Agronomica Sinica, 2019, 45(11): 1735-1745. (in Chinese)

    [11] 秦嶺, 張艷亭, 陳二影, 楊延兵, 黎飛飛, 管延安. 谷子萌發(fā)期耐鹽種質(zhì)篩選及其保護(hù)酶系統(tǒng)對(duì)鹽脅迫的響應(yīng). 中國農(nóng)業(yè)科學(xué), 2019, 52(22) : 4027-4038.

    QIN L, ZHANG Y T, CHEN E Y, YANG Y B, LI F F, GUAN Y A. Screening for germplasms tolerant to salt at germination stage and response of protective enzymes to salt stress in foxtail millet. Scientia Agricultura Sinica, 2019, 52(22): 4027-4038. (in Chinese)

    [12] 鄒成林, 譚華, 黃開建, 吳永升, 鄭德波, 黃愛花, 莫潤秀, 韋新興.廣西玉米品種苗期生理生化指標(biāo)與其抗旱性的綜合評(píng)價(jià). 南方農(nóng)業(yè)學(xué)報(bào), 2015, 46(3) : 408-414.

    ZOU C L, TAN H, HUANG K J, WU Y S, ZHENG D B, HUANG A H, MO R X, WEI X X. Comprehensive evaluation of physiological and biochemical indexes and drought resistance of maize varieties during seedling stage in Guangxi. Journal of Southern Agriculture, 2015, 46(3): 408-414. (in Chinese)

    [13] 袁雨豪, 楊清華, 黨科, 楊璞, 高金鋒, 高小麗, 王鵬科, 陸平, 劉敏軒, 馮佰利. 糜子資源耐鹽性評(píng)價(jià)與鹽脅迫生理響應(yīng). 中國農(nóng)業(yè)科學(xué), 2019, 52(22) : 4066-4078.

    YUAN Y H, YANG Q H, DANG K, YANG P, GAO J F, GAO X L, WANG P K, LU P, LIU M X, FENG B L. Salt-tolerance evaluation and physiological response of salt stress of broomcorn millet (L.). Scientia Agricultura Sinica, 2019, 52(22): 4066-4078. (in Chinese)

    [14] 李春喜, 郭雪妮, 張黛靜, 劉雪晴, 劉安琪, 王艷杰. 輕度干旱脅迫下黃淮麥區(qū)不同基因型小麥的響應(yīng)分析. 生態(tài)環(huán)境學(xué)報(bào), 2016, 25(9): 1446-1452.

    LI C X, GUO X N, ZHANG D J, LIU X Q, LIU A Q, WANG Y J. Response of different genotypes wheat to light drought stress in the Huang-Huai Plain. Ecology and Environmental Sciences, 2016, 25(9):1446-1452. (in Chinese)

    [15] 盧廣超, 許建新, 薛立, 張柔, 吳彩瓊, 邵怡若. 低溫脅迫對(duì)4種幼苗的葉綠素?zé)晒馓匦缘挠绊? 中南林業(yè)科技大學(xué)學(xué)報(bào), 2014, 34(2): 44-49.

    LU G C, XU J X, XUE L, ZHANG R, WU C Q, SHAO Y R. Effects of low temperature stress on chlorophyll fluorescence characteristics of four types of tree species seedlings. Journal of Central South University of Forestry & Technology, 2014, 34(2): 44-49. (in Chinese)

    [16] 徐瀾, 高志強(qiáng), 安偉, 原亞琦, 李彥良. 冬麥春播小麥穗分化階段對(duì)低溫脅迫的響應(yīng)及耐寒性. 應(yīng)用生態(tài)學(xué)報(bào), 2015, 26(6): 1679-1686.

    XU L, GAO Z Q, AN W, YUAN Y Q, LI Y L. Low-temperature response and cold tolerance at spike differentiation stage of winter wheat varieties sowed in spring. Chinese Journal of Applied Ecology, 2015, 26(6): 1679-1686. (in Chinese)

    [17] 姜麗娜, 張黛靜, 宋飛, 劉佩, 樊婷婷, 余海波, 李春喜. 不同品種小麥葉片對(duì)拔節(jié)期低溫的生理響應(yīng)及抗寒性評(píng)價(jià). 生態(tài)學(xué)報(bào), 2014, 34(15) : 4251-4261.

    JIANG L N, ZHANG D J, SONG F, LIU P, FAN T T, YU H B, LI C X. Evaluation of cold resistance of different wheat varieties based on physiological responses of leaves to low temperature at the jointing stage. Acta Ecologica Sinica, 2014, 34(15): 4251-4261. (in Chinese)

    [18] 李合生. 植物生理生化實(shí)驗(yàn)原理和技術(shù). 北京: 高等教育出版社, 2000: 184-185, 195-197.

    LI H S. The experiment principle and technique on plant Physiology and biochemistry. Beijing: Higher Education Press, 2000: 184-185, 195-197. (in Chinese)

    [19] 中國科學(xué)院上海植物生理研究所, 上海植物生理學(xué)會(huì). 現(xiàn)代植物生理學(xué)實(shí)驗(yàn)指南. 北京: 科學(xué)出版社, 1999: 303.

    Shanghai Institute of Plant Physiology, Chinese Academy of Sciences, Plant Physiology Academy of Shanghai. Modern Plant Physiology Experiment of Directory. Beijing: Science Press, 1999: 303. (in Chinese)

    [20] 侯福林. 植物生理學(xué)實(shí)驗(yàn)教程. 北京: 科學(xué)出版社, 2004: 91.

    HOU F L. Plant Physiology Experimental Course. Beijing: Science Press, 2004: 91. (in Chinese)

    [21] 張志良, 瞿偉菁. 植物生理學(xué)實(shí)驗(yàn)指導(dǎo). 第3版. 北京: 高等教育出版社, 2003.

    ZHANG Z L, QU W J. The Experimental Guide for Plant Physiology. 3rd edition. Beijing: Higher Education Press, 2003. (in Chinese)

    [22] Zou J, Hu W, LI Y X, He J Q, Zhu H H, Zhou Z G. Screening of drought resistance indices and evaluation of drought resistance in cotton (L.). Journal of Integrative Agriculture, 2020, 19(2): 495-508.

    [23] 王秀田, 盧秋巍, 蒼晶, 包雨卓, 孟德義, 于晶, 徐慶華, 趙虎, 李想, 相智也, 謝波. 低溫馴化對(duì)冬小麥葉綠素?zé)晒馓匦约翱寡趸富钚缘挠绊? 植物生理學(xué)報(bào), 2016, 52(12): 1959-1969.

    Wang X T, Lu Q W, Cang J, Bao Y Z, Meng D Y, Yu J, Xu Q H, Zhao H, Li X, Xiang Z Y, Xie B. Effects of cold acclimation on chlorophyll fluorescence characteristics and antioxidant enzyme activity in winter wheat. Plant Physiology Journal, 2016, 52(12): 1959-1969. (in Chinese)

    [24] 李春燕, 陳思思, 徐雯, 李東升, 顧驍, 朱新開, 郭文善, 封超年. 苗期低溫脅迫對(duì)揚(yáng)麥16葉片抗氧化酶和滲透調(diào)節(jié)物質(zhì)的影響. 作物學(xué)報(bào), 2011, 37(12) : 2293-2298.

    LI C Y, CHEN S S, XU W, LI D S, GU X, ZHU X K, GUO W S, FENG C N. Effects of low temperature at seedling stage on antioxidant enzymes and cytoplasmic osmoticum of leaves in wheat cultivar Yangmai 16. Acta Agronomica Sinica, 2011, 37(12): 2293-2298. (in Chinese)

    [25] 趙慧, 趙一博, 郭江波, 辛翠花. 植物耐受低溫脅迫研究進(jìn)展. 種子, 2017, 36(5): 47-50.

    ZHAO H, ZHAO Y B, GUO J B, XIN C H. Research progress of plant under low temperature stress. Seed, 2017, 36(5): 47-50. (in Chinese)

    [26] 劉海卿, 孫萬倉, 劉自剛, 武軍艷, 錢武, 王志江, 郭仁迪, 馬驪, 侯獻(xiàn)飛, 劉林波. 北方寒旱區(qū)白菜型冬油菜抗寒性與抗旱性評(píng)價(jià)及其關(guān)系. 中國農(nóng)業(yè)科學(xué), 2015, 48(18) : 3743-3756.

    LIU H Q, SUN W C, LIU Z G, WU J Y, QIAN W, WANG Z J, GUO R D, MA L, HOU X F, LIU L B. Evaluation of drought resistance and cold resistance and research of their relationship at seedling stage of winter rapeseed (L.) in cold and arid regions in North China. Scientia Agricultura Sinica, 2015, 48(18): 3743-3756. (in Chinese)

    [27] Rubio M C, James E K, Clemente M R, Bucciarelli B, Fedorova M, Vance C P, Becana M. Localization of superoxide dismutases and hydrogen peroxide in legume root nodules. Molecular Plant-Microbe Interactions, 2004, 17(12): 1294-1305.

    [28] 王曉楠, 付連雙, 李卓夫, 孫艷麗, 王玉波, 劉燦, 王金偉, 陳禹興. 低溫馴化及封凍后不同抗寒性小麥品種的形態(tài)建成及生理基礎(chǔ)分析. 作物學(xué)報(bào), 2009, 35(7): 1313-1319.

    Wang X N, Fu L S, Li Z F, Sun Y L, Wang Y B, Liu C, Wang J W, Chen Y Y. Morphogenesis and physiological basis in wheat cultivars with different levels of cold-resistance during cold acclimation and freezing period. Acta Agronomica Sinica, 2009, 35(7): 1313-1319. (in Chinese)

    [29] 劉立軍, 陳為峰, 王凱, 張呈祥. 不同灌水量下限對(duì)高羊茅綠期及抗寒性生理指標(biāo)的影響. 應(yīng)用生態(tài)學(xué)報(bào), 2012, 23(1): 38-44.

    LIU L J, CHEN W F, WANG K, ZHANG C X. Effects of different irrigation minima on green period and cold-resistance physiological indices of Festuca arundinacea. Chinese Journal of Applied Ecology, 2012, 23(1): 38-44. (in Chinese)

    [30] 譚孟祥, 景元書, 薛楊, 曾文全. 水層深度對(duì)早稻幼穗分化期遭遇低溫過程時(shí)葉片生理特性的影響. 中國農(nóng)業(yè)氣象, 2015, 36(5): 553-560.

    Tan M X, Jing Y S, XUE Y, Zeng W Q. Effects of different water depth on leaf physiological characteristics of early rice during panicle primordium suffered to low temperature. Chinese Journal of Agrometeorology, 2015, 36(5): 553-560. (in Chinese)

    [31] 曹娜, 熊強(qiáng)強(qiáng), 陳小榮, 賀浩華, 朱昌蘭, 傅軍如, 才碩, 徐濤. 幼穗分化期低溫條件下灌水對(duì)早稻溫度調(diào)節(jié)、產(chǎn)量形成及相關(guān)生理特性的影響. 核農(nóng)學(xué)報(bào), 2018, 32(4): 760-770.

    Cao N, Xiong Q Q, Chen X R, He H H, Zhu C L, Fu J R, Cai S, Xu T. Effects of irrigation on temperature regulation, yield formation and related physiological characteristics in early rice under low temperature conditions during panicle primordium differentiation stage. Journal of Nuclear Agricultural Sciences, 2018, 32(4): 760-770. (in Chinese)

    [32] 曹娜, 陳小榮, 賀浩華, 朱昌蘭, 才碩, 徐濤, 謝亨旺, 劉方平. 抽穗揚(yáng)花期不同灌水處理對(duì)晚稻抵御低溫、產(chǎn)量和生理特性的影響.應(yīng)用生態(tài)學(xué)報(bào), 2017, 28(12): 3935-3944.

    CAO N, CHEN X R, HE H H, ZHU C L, CAI S, XU T, XIE H W, LIU F P. Effects of different irrigation treatments during heading and flowering stage on cold resistance, yield and physiological characteristics of late rice. Chinese Journal of Applied Ecology, 2017, 28(12): 3935-3944. (in Chinese)

    [33] 單旭東, 高婭楠, 張桐瑞, 程錦, 許立新. 保溫和澆水處理對(duì)結(jié)縷草冬季脅迫的緩解作用. 草地學(xué)報(bào), 2021, 29(1): 202-210.

    SHAN X D, GAO Y N, ZHANG T R, CHENG J, XU L X. Alleviating effects of heat preservation and watering treatments on winter stress of zoysia japonica. Acta Agrestia Sinica, 2021, 29(1): 202-210. (in Chinese)

    [34] 曹倩, 王顯國, 張尚雄, 劉博文, 金美艷, 寧亞明. 抗寒鍛煉前土壤水分含量對(duì)苜??购缘挠绊? 中國草地學(xué)報(bào), 2021, 43(3): 1-9.

    CAO Q, WANG X G, ZHANG S X, LIU B W, JIN M Y, NING Y M. The effect of soil water content before cold hardening on cold resistance of alfalfa. Chinese Journal of Grassland, 2021, 43(3): 1-9. (in Chinese)

    [35] 汪燦, 周棱波, 張國兵, 張立異, 徐燕, 高旭, 姜訥, 邵明波. 薏苡種質(zhì)資源苗期抗旱性鑒定及抗旱指標(biāo)篩選. 中國農(nóng)業(yè)科學(xué), 2017, 50(15) : 2872-2887.

    WANG C, ZHOU L B, ZHANG G B, ZHANG L Y, XU Y, GAO X, JIANG N, SHAO M B. Drought resistance identification and drought resistance indices screening of Job’s Tears (L.) germplasm resources at seedling stage. Scientia Agricultura Sinica, 2017, 50(15): 2872-2887. (in Chinese)

    [36] Naroui Rad m R, Abbasi m r, Fanaei H R. Evaluation of drought stress tolerance with use of stress tolerance indexes in sorghum collected germplasms national plant gene bank of Iran. Agronomy, 2009, 22(1): 11-18.

    [37] 張軍, 魯敏, 孫樹貴, 杜萬里, 劉洋, 武軍, 陳新宏. 拔節(jié)期低溫脅迫對(duì)小麥生理生化特性和產(chǎn)量的影響. 西北農(nóng)業(yè)學(xué)報(bào), 2014, 23(2) : 73-79.

    ZHANG J, LU M, SUN S G, DU W L, LIU Y, WU J, CHEN X H. Changes of physiological and biochemical parameters and grain yield at jointing stage of wheat under low temperature stress. Acta Agriculturae Boreali-occidentalis Sinica, 2014, 23(2): 73-79. (in Chinese)

    [38] 劉建新, 歐曉彬, 劉秀麗, 王金成. 過氧化氫緩解裸燕麥幼苗低溫脅迫的主成分和隸屬函數(shù)分析. 植物研究, 2018, 38(5): 748-756.

    LIU J X, OU X B, LIU X L, WANG J C. Principal component and subordinate function of the alleviating effects of hydrogen peroxide (H2O2) on low-temperature stress in naked oat () seedlings. Bulletin of Botanical Research, 2018, 38(5): 748-756. (in Chinese)

    [39] 邵怡若, 許建新, 薛立, 張柔, 吳彩瓊, 盧廣超. 低溫脅迫時(shí)間對(duì)4種幼苗生理生化及光合特性的影響. 生態(tài)學(xué)報(bào), 2013, 33(14): 4237-4247.

    SHAO Y R, XU J X, XUE L, ZHANG R, WU C Q, LU G C. Effects of low temperature stress on physiological-biochemical indexes and photosynthetic characteristics of seedlings of four plant species. Acta Ecologica Sinica, 2013, 33(14): 4237-4247. (in Chinese)

    [40] 胡廷會(huì), 成良強(qiáng), 王軍, 呂建偉, 饒慶琳. 不同基因型花生耐蔭性評(píng)價(jià)及其鑒定指標(biāo)的篩選. 中國農(nóng)業(yè)科學(xué), 2020, 53(6) : 1140-1153.

    HU T H, CHENG L Q, WANG J, Lü J W, RAO Q L. Evaluation of shade tolerance of peanut with different genotypes and screening of identification indexes. Scientia Agricultura Sinica, 2020, 53(6): 1140-1153. (in Chinese)

    [41] 陳二影, 王潤豐, 秦嶺, 楊延兵, 黎飛飛, 張華文, 王海蓮, 劉賓, 孔清華, 管延安. 谷子芽期耐鹽堿綜合鑒定及評(píng)價(jià). 作物學(xué)報(bào), 2020, 46(10): 1591-1604.

    CHEN E Y, WANG R F, QIN L, YANG Y B, LI F F, ZHANG H W, WANG H L, LIU B, KONG Q H, GUAN Y A. Comprehensive identification and evaluation of foxtail millet for saline-alkaline tolerance during germination. Acta Agronomica Sinica, 2020, 46(10): 1591-1604. (in Chinese)

    [42] 王濤, 王曉楠, 王明芳, 孫瑩璐, 李卓夫, 付連雙. 低溫下冬小麥水分含量的變化與抗寒性鑒定. 作物雜志, 2015(1): 61-66.

    Wang T, Wang X N, Wang M F, Sun Y L, Li Z F, Fu L S. Moisture content changes and cold resistance identification by winter wheat varieties at low temperature. Crops, 2015(1): 61-66. (in Chinese)

    [43] 趙瑞玲, 趙勇, 徐渴, 李嘉豪, 張樹華, 楊學(xué)舉. 室內(nèi)冷凍法鑒定小麥抗寒性的研究. 植物遺傳資源學(xué)報(bào), 2019, 20(2): 284-296.

    Zhao R L, Zhao Y, Xu K, Li J H, Zhang S H, Yang X J. Identification of cold resistance wheat varieties by laboratory freezing method. Journal of Plant Genetic Resources, 2019, 20(2): 284-296. (in Chinese)

    [44] 張自陽, 王智煜, 王斌, 王志偉, 朱啟迪, 霍云風(fēng), 茹振鋼, 劉明久. 春季穗分化階段低溫處理對(duì)不同小麥品種幼穗結(jié)實(shí)性及生理特性的影響. 華北農(nóng)學(xué)報(bào), 2019, 34(4): 130-139.

    Zhang Z Y, Wang Z Y, Wang B, Wang Z W, Zhu Q D, Huo Y F, Ru Z G, Liu M J. Effects of low temperature treatment at spring spike differentiation stage on young ear fruiting and physiological characteristics of different wheat varieties. Acta Agriculturae Boreali-Sinica, 2019, 34(4): 130-139. (in Chinese)

    [45] 吳青霞, 楊林, 邵慧, 冉從福, 楊子博, 余靜, 李立群, 李學(xué)軍. 藥隔期低溫脅迫對(duì)小麥生理及產(chǎn)量的影響. 麥類作物學(xué)報(bào), 2013, 33(4): 752-757.

    Wu Q X, Yang L, Shao H, Ran C F, Yang Z B, Yu J, Li L Q, Li X J. Effects of low temperature stress at anther connective formation phase on physiological characteristics and yield of wheat. Journal of Triticeae Crops, 2013, 33(4): 752-757. (in Chinese)

    Evaluation of low temperature freezing injury in winter wheat and difference analysis of water effect based on multivariate statistical analysis

    WANG YangYang1, 2, LIU WanDai1, 2, HE Li1, 2, REN DeChao3, DUAN JianZhao1, 2, HU Xin3, GUO TianCai1, 2, WANG YongHua1, 2, FENG Wei1, 2

    1College of Agronomy, Henan Agriculture University/Key Laboratory of Regulating and Controlling Crop Growth and Development, Ministry of Education, Zhengzhou 450046;2National Engineering Research Center for Wheat, Zhengzhou 450046;3Wheat Research Institute, Shangqiu Academy of Agriculture and Forestry Sciences, Shangqiu 476000, Henan

    In order to clarify the freezing injury degree of wheat under different water conditions caused by low temperature stress, the identification indexes and quantitative evaluation model of freezing injury were screened and established, which provided the theoretical support for prevention and control of freezing injury in wheat production.Weak spring cultivars of Yanzhan 4110 and Lankao 198, semi-winter cultivars of Zhengmai 366 and Fengdecunmai 21 were used as experimental materials. They were treated with irrigation (W) or no irrigation (D) one week before the freezing injury, respectively. Pot experiments were moved to a low-temperature simulation room during the female and male ear differentiation stages. The temperatures were set as -2℃ (T1), -4℃ (T2), -6℃ (T3), -8℃ (T4), -10℃ (T5) and control (CK is the field temperature on the same day). Physiological and biochemical indexes of wheat were measured on the second day after low temperature stress. The standardized physiological indexes were analyzed by multivariate statistical analysis, such as principal component, membership function, cluster analysis and step wise regression.There were significant correlations among the individual physiological and biochemical indexes under different cultivars, water contents and temperatures. Through principal component analysis, 19 physiological and biochemical indexes were transformed into 6 mutually independent comprehensive indexes, whose contribution degrees were 55.972%, 11.93%, 7.168%, 5.075%, 4.236% and 3.079%, respectively, representing 87.459% information of all original data. According to the membership function algorithm, the comprehensive evaluation value (value ) of freezing injury degree of each treatment was calculated. Takevalue as the dependent variable, the seven key indexes were selected by stepwise regression analysis, namely chlorophyll a, leaf water content, proline,v/m, soluble protein, MDA and SOD, and the mathematical model for quantitative estimation ofvalue was established. At the same time, the correlation betweenprediction value and yield loss rate was analyzed, and the linear equation determination coefficient2= 0.898, indicating that theprediction model could well evaluate the freezing injury degree.predicted value could be divided into five categories by further cluster analysis: non-freezing (D-CK, W-CK), mild frozen (D-T1, W-T1), moderate frozen (D-T2, W-T2, W-T3), severe frozen (D-T3, W-T4), and extremely severe frozen (D-T4, W-T5, D-T5). Corresponding yield loss rate were 0, 0-10%, 10%-30%, 30%-50% and more than 50%, respectively. Under the same temperature and moisture conditions, the freezing injury degree of weak spring varieties was heavier than that of semi-winter varieties, and the freezing injury degree of no irrigation treatment was heavier than that of irrigation treatment under the same varieties and temperature conditions. With the increasing of low temperature stress, chlorophyll a, leaf water content andv/mshowed a decreasing trend, the activities of proline, soluble protein and SOD increased first and then decreased, while MDA showed an opposite trend. According to the clustering results, under the same temperature and water conditions, the freezing injury degree of weak spring cultivars was more serious than that of semi-winter cultivars. Under the same variety and temperature conditions, the freezing damage degree without irrigation was worse than that under irrigation.Therefore, the semi-winter varieties should be selected in the areas prone to late frost in production, and the irrigation management should be strengthened before a cold wave according to the weather forecast. When freezing injury happened, the injury degree could be accurately assessed in timely through the evaluation index and quantitative model, which was conducive to prevention and control of late frost injury, and provides technology basis for production recovery and decision management after freezing disaster.

    winter wheat; frost damage degree; irrigation; comprehensive evaluation; estimation model

    10.3864/j.issn.0578-1752.2022.07.004

    2021-06-09;

    2021-10-08

    “十三五”國家重點(diǎn)研發(fā)計(jì)劃“糧食豐產(chǎn)增效科技創(chuàng)新”(2017YFD0300204)

    王洋洋,E-mail:wyy65wyy@163.com。通信作者賀利,E-mail:he-li19870308@163.com。通信作者馮偉,E-mail:fengwei78@126.com

    (責(zé)任編輯 楊鑫浩)

    猜你喜歡
    抗寒性灌水生理
    灌水取球
    番茄灌水掌握技巧
    冬季棚菜灌水四關(guān)鍵
    打破生理“平衡”
    基于BP神經(jīng)網(wǎng)絡(luò)的旋轉(zhuǎn)血泵生理控制
    媽媽們產(chǎn)后的生理煩惱
    Coco薇(2017年5期)2017-06-05 13:03:24
    棗樹抗寒性檢測(cè)方法的篩選
    巴梨的抗寒性試驗(yàn)情況初報(bào)
    中國果菜(2016年9期)2016-03-01 01:28:40
    馬鈴薯普通栽培種雜交后代抗寒性分析
    中國馬鈴薯(2015年5期)2016-01-09 06:11:11
    灌水秘笈
    国产精品精品国产色婷婷| 久久精品91无色码中文字幕| 最新在线观看一区二区三区| 国产色婷婷99| 我的老师免费观看完整版| 色噜噜av男人的天堂激情| 国产精品三级大全| 狂野欧美白嫩少妇大欣赏| 国产精品一区二区免费欧美| 国产精品99久久99久久久不卡| 国产男靠女视频免费网站| 99在线人妻在线中文字幕| 国产精华一区二区三区| 国产成人a区在线观看| 又黄又爽又免费观看的视频| 男女下面进入的视频免费午夜| a级毛片a级免费在线| 久久久久久大精品| 狂野欧美白嫩少妇大欣赏| 欧美色欧美亚洲另类二区| 真实男女啪啪啪动态图| 夜夜看夜夜爽夜夜摸| 草草在线视频免费看| 欧美性感艳星| 老司机午夜十八禁免费视频| 久久精品国产亚洲av涩爱 | 97碰自拍视频| 中文字幕熟女人妻在线| 少妇的丰满在线观看| 欧美高清成人免费视频www| 成人av一区二区三区在线看| 露出奶头的视频| 精品熟女少妇八av免费久了| 真实男女啪啪啪动态图| 麻豆久久精品国产亚洲av| 色吧在线观看| 日本一二三区视频观看| 国产激情欧美一区二区| 波多野结衣高清无吗| 女生性感内裤真人,穿戴方法视频| 亚洲成av人片免费观看| 窝窝影院91人妻| 三级国产精品欧美在线观看| 琪琪午夜伦伦电影理论片6080| 精品一区二区三区视频在线 | 亚洲无线在线观看| 两人在一起打扑克的视频| 少妇的逼水好多| 嫩草影视91久久| 国产一区二区激情短视频| 亚洲欧美日韩高清专用| 波多野结衣巨乳人妻| 免费在线观看影片大全网站| 琪琪午夜伦伦电影理论片6080| bbb黄色大片| 日日夜夜操网爽| 两个人的视频大全免费| 欧美性猛交╳xxx乱大交人| 非洲黑人性xxxx精品又粗又长| 国产精品野战在线观看| 人人妻人人澡欧美一区二区| 亚洲av熟女| 老司机午夜福利在线观看视频| 国产亚洲精品一区二区www| 亚洲真实伦在线观看| 中文字幕熟女人妻在线| 精品久久久久久成人av| 夜夜看夜夜爽夜夜摸| 久久久久久久久久黄片| 久久久国产成人免费| 人妻丰满熟妇av一区二区三区| 999久久久精品免费观看国产| 亚洲,欧美精品.| 十八禁人妻一区二区| 天堂影院成人在线观看| 成人特级黄色片久久久久久久| 国产精品一及| 精品欧美国产一区二区三| 国产色婷婷99| 日本与韩国留学比较| 制服丝袜大香蕉在线| 女人被狂操c到高潮| 又黄又粗又硬又大视频| 免费av毛片视频| 最近最新免费中文字幕在线| 亚洲av成人精品一区久久| 国产午夜福利久久久久久| 中文字幕人妻丝袜一区二区| 非洲黑人性xxxx精品又粗又长| 搡老熟女国产l中国老女人| av中文乱码字幕在线| 18禁在线播放成人免费| 国产精品 欧美亚洲| 一个人看的www免费观看视频| 国产精品香港三级国产av潘金莲| 嫩草影院入口| 一本精品99久久精品77| 精品人妻偷拍中文字幕| 性色av乱码一区二区三区2| 女人被狂操c到高潮| 岛国在线观看网站| 久久婷婷人人爽人人干人人爱| 一区二区三区免费毛片| 可以在线观看毛片的网站| 色哟哟哟哟哟哟| 国产一区二区在线av高清观看| 午夜老司机福利剧场| 婷婷丁香在线五月| 国产免费av片在线观看野外av| 小说图片视频综合网站| 91字幕亚洲| 亚洲精品成人久久久久久| 99精品欧美一区二区三区四区| 国产黄a三级三级三级人| 欧美不卡视频在线免费观看| 午夜福利在线观看免费完整高清在 | 国产精品乱码一区二三区的特点| 久久久久免费精品人妻一区二区| www日本在线高清视频| 欧美中文综合在线视频| 免费在线观看成人毛片| 久久久久久人人人人人| 午夜两性在线视频| 国产午夜福利久久久久久| 亚洲av二区三区四区| 欧美+日韩+精品| 黄色成人免费大全| 99久久九九国产精品国产免费| 久久香蕉精品热| 每晚都被弄得嗷嗷叫到高潮| av在线蜜桃| 好男人电影高清在线观看| 国产乱人伦免费视频| 真人做人爱边吃奶动态| 一区二区三区国产精品乱码| 亚洲专区国产一区二区| 少妇人妻一区二区三区视频| 男女床上黄色一级片免费看| 在线观看美女被高潮喷水网站 | 亚洲avbb在线观看| 日本撒尿小便嘘嘘汇集6| 国产在线精品亚洲第一网站| 亚洲国产精品合色在线| 大型黄色视频在线免费观看| 亚洲美女视频黄频| 九九久久精品国产亚洲av麻豆| 精华霜和精华液先用哪个| 无遮挡黄片免费观看| 久久欧美精品欧美久久欧美| 在线观看av片永久免费下载| 久久精品人妻少妇| 一进一出抽搐gif免费好疼| 亚洲国产高清在线一区二区三| 精华霜和精华液先用哪个| 日韩 欧美 亚洲 中文字幕| 亚洲av成人av| 免费观看人在逋| 女生性感内裤真人,穿戴方法视频| 热99在线观看视频| 国产伦精品一区二区三区视频9 | 天堂√8在线中文| 男女下面进入的视频免费午夜| 精品久久久久久久毛片微露脸| 啪啪无遮挡十八禁网站| 亚洲人成电影免费在线| 国产成人影院久久av| 18禁裸乳无遮挡免费网站照片| 嫩草影院精品99| 久久香蕉国产精品| 黄片小视频在线播放| 一a级毛片在线观看| 久久国产乱子伦精品免费另类| 国产爱豆传媒在线观看| 99视频精品全部免费 在线| 亚洲专区国产一区二区| 少妇高潮的动态图| 老司机深夜福利视频在线观看| 可以在线观看毛片的网站| 脱女人内裤的视频| 精品一区二区三区av网在线观看| 久久久久精品国产欧美久久久| 色综合欧美亚洲国产小说| 国产中年淑女户外野战色| 午夜激情欧美在线| 欧美日韩福利视频一区二区| 国产欧美日韩一区二区精品| 亚洲国产中文字幕在线视频| 亚洲av电影不卡..在线观看| 国产成人欧美在线观看| 好看av亚洲va欧美ⅴa在| 国产精品国产高清国产av| 亚洲va日本ⅴa欧美va伊人久久| 国产午夜精品论理片| 久久久久久久久中文| 亚洲av中文字字幕乱码综合| 中文字幕高清在线视频| 国产激情偷乱视频一区二区| 看黄色毛片网站| 麻豆一二三区av精品| a级一级毛片免费在线观看| 欧美一级a爱片免费观看看| 国产成人a区在线观看| 一级毛片女人18水好多| 成人一区二区视频在线观看| 成人一区二区视频在线观看| 国产精品电影一区二区三区| 亚洲人成网站在线播| 亚洲国产精品成人综合色| 欧美激情在线99| 一区二区三区免费毛片| 久久久久九九精品影院| 男人的好看免费观看在线视频| 免费看十八禁软件| 性色avwww在线观看| 禁无遮挡网站| 成人av在线播放网站| 亚洲国产色片| 欧美日韩福利视频一区二区| 很黄的视频免费| 亚洲人成伊人成综合网2020| 久久久久性生活片| 免费看日本二区| 欧美日韩亚洲国产一区二区在线观看| 精品99又大又爽又粗少妇毛片 | 日韩人妻高清精品专区| 伊人久久大香线蕉亚洲五| 久久久久亚洲av毛片大全| av中文乱码字幕在线| 在线视频色国产色| 欧美bdsm另类| 国产精品久久久人人做人人爽| 88av欧美| 综合色av麻豆| 91在线精品国自产拍蜜月 | 18+在线观看网站| 成人永久免费在线观看视频| 国产熟女xx| h日本视频在线播放| 亚洲成人免费电影在线观看| 亚洲中文字幕一区二区三区有码在线看| 成年女人看的毛片在线观看| 1000部很黄的大片| 嫩草影视91久久| 搡老岳熟女国产| 91麻豆精品激情在线观看国产| 国产成年人精品一区二区| 九色成人免费人妻av| 俄罗斯特黄特色一大片| 性色av乱码一区二区三区2| 久久天躁狠狠躁夜夜2o2o| 成人无遮挡网站| 麻豆成人午夜福利视频| 女警被强在线播放| 亚洲一区高清亚洲精品| 欧美中文综合在线视频| 亚洲在线自拍视频| 脱女人内裤的视频| 国产成+人综合+亚洲专区| 欧美日韩瑟瑟在线播放| 黄色成人免费大全| 三级国产精品欧美在线观看| 免费看a级黄色片| 在线观看一区二区三区| 国产精品一区二区三区四区久久| 亚洲国产欧美人成| 精品免费久久久久久久清纯| 欧美xxxx黑人xx丫x性爽| 麻豆国产av国片精品| 非洲黑人性xxxx精品又粗又长| 欧美极品一区二区三区四区| 一本精品99久久精品77| 免费一级毛片在线播放高清视频| 搡老熟女国产l中国老女人| 免费无遮挡裸体视频| 动漫黄色视频在线观看| 男人舔奶头视频| av欧美777| 久久久久久久久久黄片| 极品教师在线免费播放| 天天躁日日操中文字幕| 麻豆成人午夜福利视频| 国产单亲对白刺激| 国产在线精品亚洲第一网站| 国产亚洲精品综合一区在线观看| 日韩精品中文字幕看吧| 午夜两性在线视频| 久久欧美精品欧美久久欧美| 久久精品国产自在天天线| 少妇的逼好多水| 欧美成人性av电影在线观看| 免费观看的影片在线观看| 成年免费大片在线观看| 尤物成人国产欧美一区二区三区| 天堂av国产一区二区熟女人妻| www日本在线高清视频| 精品国产超薄肉色丝袜足j| 精品久久久久久,| 91麻豆精品激情在线观看国产| 久久久国产精品麻豆| 岛国在线免费视频观看| 国产精华一区二区三区| 亚洲精品456在线播放app | 国产不卡一卡二| 国产av不卡久久| 老鸭窝网址在线观看| x7x7x7水蜜桃| 久久久久精品国产欧美久久久| 99久久99久久久精品蜜桃| av天堂中文字幕网| 日韩精品青青久久久久久| www.熟女人妻精品国产| 免费一级毛片在线播放高清视频| 成人国产综合亚洲| 国产爱豆传媒在线观看| 热99在线观看视频| 天天添夜夜摸| 美女cb高潮喷水在线观看| 久久久国产成人精品二区| 国产高清三级在线| 又黄又粗又硬又大视频| 精品一区二区三区人妻视频| 日韩欧美 国产精品| 国产91精品成人一区二区三区| 久9热在线精品视频| 欧美乱码精品一区二区三区| 亚洲成人中文字幕在线播放| 色精品久久人妻99蜜桃| 亚洲中文日韩欧美视频| 草草在线视频免费看| 最后的刺客免费高清国语| 看免费av毛片| 久久精品亚洲精品国产色婷小说| 精品一区二区三区视频在线观看免费| 久久这里只有精品中国| 欧美日韩中文字幕国产精品一区二区三区| 嫩草影院入口| 国产精品亚洲av一区麻豆| 看免费av毛片| 99久久九九国产精品国产免费| 亚洲精品一区av在线观看| 一进一出好大好爽视频| 九九热线精品视视频播放| 亚洲激情在线av| 制服丝袜大香蕉在线| 嫩草影院精品99| 欧美乱妇无乱码| 国产高清videossex| 国产免费男女视频| 此物有八面人人有两片| 美女黄网站色视频| 99久久无色码亚洲精品果冻| 热99re8久久精品国产| 日韩大尺度精品在线看网址| 一个人看的www免费观看视频| 色在线成人网| 国产一区二区激情短视频| 午夜老司机福利剧场| 法律面前人人平等表现在哪些方面| 99久久99久久久精品蜜桃| 一区二区三区高清视频在线| 一进一出好大好爽视频| 国产高清有码在线观看视频| 午夜精品久久久久久毛片777| 99久久成人亚洲精品观看| 天天一区二区日本电影三级| 最好的美女福利视频网| 波多野结衣高清作品| 亚洲人与动物交配视频| 免费一级毛片在线播放高清视频| 亚洲专区中文字幕在线| 变态另类丝袜制服| 一本久久中文字幕| 国产 一区 欧美 日韩| 亚洲欧美日韩东京热| 国产午夜精品久久久久久一区二区三区 | 给我免费播放毛片高清在线观看| 又爽又黄无遮挡网站| 制服人妻中文乱码| 国产精品1区2区在线观看.| 色噜噜av男人的天堂激情| 最近视频中文字幕2019在线8| bbb黄色大片| 夜夜躁狠狠躁天天躁| 国产精品一区二区三区四区免费观看 | www国产在线视频色| 9191精品国产免费久久| 亚洲 国产 在线| 国产av不卡久久| 国产淫片久久久久久久久 | 高清毛片免费观看视频网站| 高清在线国产一区| 黄色成人免费大全| 久久婷婷人人爽人人干人人爱| 亚洲av熟女| 99久久综合精品五月天人人| 亚洲午夜理论影院| 91字幕亚洲| 国产一区二区亚洲精品在线观看| 给我免费播放毛片高清在线观看| 看免费av毛片| 啦啦啦免费观看视频1| 美女cb高潮喷水在线观看| 成人欧美大片| 精品一区二区三区人妻视频| 天堂√8在线中文| 最后的刺客免费高清国语| 亚洲国产中文字幕在线视频| 国内精品一区二区在线观看| 国产三级中文精品| 女生性感内裤真人,穿戴方法视频| 1024手机看黄色片| 午夜精品一区二区三区免费看| 久久香蕉精品热| 香蕉av资源在线| 制服丝袜大香蕉在线| 午夜激情福利司机影院| 两个人看的免费小视频| 在线观看免费视频日本深夜| av欧美777| 亚洲aⅴ乱码一区二区在线播放| 欧美丝袜亚洲另类 | 一个人看的www免费观看视频| 中文字幕人成人乱码亚洲影| 嫩草影院精品99| 日韩有码中文字幕| 国产精品一区二区免费欧美| 日韩成人在线观看一区二区三区| 国产国拍精品亚洲av在线观看 | 亚洲人成伊人成综合网2020| 免费av观看视频| 国内精品美女久久久久久| 亚洲熟妇熟女久久| 无人区码免费观看不卡| 亚洲中文日韩欧美视频| 欧美一区二区精品小视频在线| 午夜免费观看网址| 亚洲精品影视一区二区三区av| 国产精品久久电影中文字幕| 欧美精品啪啪一区二区三区| 亚洲片人在线观看| 好看av亚洲va欧美ⅴa在| 亚洲成人精品中文字幕电影| 搡老妇女老女人老熟妇| www国产在线视频色| 搞女人的毛片| 91av网一区二区| 99国产综合亚洲精品| 亚洲第一欧美日韩一区二区三区| 看黄色毛片网站| 欧美av亚洲av综合av国产av| 精品免费久久久久久久清纯| 婷婷丁香在线五月| 欧美不卡视频在线免费观看| 一进一出好大好爽视频| 免费观看的影片在线观看| 看免费av毛片| 欧美又色又爽又黄视频| 天美传媒精品一区二区| 熟女少妇亚洲综合色aaa.| 亚洲色图av天堂| 精品乱码久久久久久99久播| 亚洲国产色片| 成人精品一区二区免费| 久久久久久人人人人人| 嫩草影院精品99| 欧美中文日本在线观看视频| 免费无遮挡裸体视频| 他把我摸到了高潮在线观看| 欧美乱妇无乱码| 免费一级毛片在线播放高清视频| 又粗又爽又猛毛片免费看| 亚洲精品国产精品久久久不卡| av片东京热男人的天堂| 老司机深夜福利视频在线观看| 热99在线观看视频| 成年人黄色毛片网站| 免费一级毛片在线播放高清视频| 日本与韩国留学比较| 两个人视频免费观看高清| 国产精品 欧美亚洲| 国产三级中文精品| 嫁个100分男人电影在线观看| 成人高潮视频无遮挡免费网站| 麻豆国产97在线/欧美| 国产成人欧美在线观看| 国产精品影院久久| 一级黄片播放器| 黄色日韩在线| 亚洲av电影在线进入| 国产 一区 欧美 日韩| 91av网一区二区| 国产美女午夜福利| 国产老妇女一区| 久久久久久久精品吃奶| 色综合婷婷激情| 99精品欧美一区二区三区四区| 人妻久久中文字幕网| 国产三级中文精品| 久久性视频一级片| 手机成人av网站| or卡值多少钱| 精品久久久久久久久久免费视频| 精华霜和精华液先用哪个| 国产麻豆成人av免费视频| 日韩中文字幕欧美一区二区| 女人十人毛片免费观看3o分钟| 天堂av国产一区二区熟女人妻| 免费电影在线观看免费观看| 中文在线观看免费www的网站| 97人妻精品一区二区三区麻豆| 国产淫片久久久久久久久 | 国产精品 欧美亚洲| 久久国产乱子伦精品免费另类| 国内精品久久久久久久电影| 一区二区三区激情视频| 国产精品自产拍在线观看55亚洲| 欧美中文综合在线视频| 成人国产一区最新在线观看| 久久久久国内视频| 中文字幕av在线有码专区| 国产蜜桃级精品一区二区三区| aaaaa片日本免费| 国产99白浆流出| 丝袜美腿在线中文| 国产激情欧美一区二区| 3wmmmm亚洲av在线观看| 18禁黄网站禁片午夜丰满| 欧美绝顶高潮抽搐喷水| 国产亚洲精品av在线| 一卡2卡三卡四卡精品乱码亚洲| 亚洲成人久久性| 国产精品一区二区三区四区久久| 日韩欧美精品免费久久 | 最好的美女福利视频网| 亚洲无线在线观看| 女人十人毛片免费观看3o分钟| 久久久久精品国产欧美久久久| 亚洲成人免费电影在线观看| 成年女人永久免费观看视频| 精品无人区乱码1区二区| 欧美精品啪啪一区二区三区| 国产91精品成人一区二区三区| 俺也久久电影网| 老司机深夜福利视频在线观看| 国产麻豆成人av免费视频| 91在线观看av| eeuss影院久久| 国产成人欧美在线观看| 久久香蕉精品热| АⅤ资源中文在线天堂| 国产精品 国内视频| 欧美日本视频| 色噜噜av男人的天堂激情| 久久性视频一级片| 最近最新中文字幕大全电影3| 搞女人的毛片| 丰满的人妻完整版| 欧美乱妇无乱码| 午夜免费男女啪啪视频观看 | 国产一区二区三区视频了| 欧美在线一区亚洲| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 国产av不卡久久| 国产视频内射| 又粗又爽又猛毛片免费看| 最近视频中文字幕2019在线8| 国产69精品久久久久777片| 国产精品久久久久久久电影 | 性欧美人与动物交配| 精品国内亚洲2022精品成人| 中文字幕av成人在线电影| 一区福利在线观看| 欧美在线一区亚洲| 亚洲男人的天堂狠狠| 国产免费男女视频| 国产三级在线视频| 国产成年人精品一区二区| 又爽又黄无遮挡网站| 欧美日韩黄片免| 久久久久久久午夜电影| 欧美中文日本在线观看视频| 不卡一级毛片| 欧美日韩福利视频一区二区| netflix在线观看网站| 有码 亚洲区| 啦啦啦韩国在线观看视频| 国产不卡一卡二| 午夜激情欧美在线| 最好的美女福利视频网| 日日干狠狠操夜夜爽| 亚洲av成人精品一区久久| 国产欧美日韩精品亚洲av| 很黄的视频免费| 成人av一区二区三区在线看| x7x7x7水蜜桃| 亚洲人成网站在线播放欧美日韩| 亚洲av美国av| 成人鲁丝片一二三区免费| 精品人妻1区二区| 免费观看的影片在线观看| 国产精品乱码一区二三区的特点| 最近视频中文字幕2019在线8| 色尼玛亚洲综合影院| 亚洲成人中文字幕在线播放| 欧美成人免费av一区二区三区| 国产精品久久久久久亚洲av鲁大| 99riav亚洲国产免费| 国产亚洲精品一区二区www| 天堂√8在线中文| 亚洲av不卡在线观看| 人人妻,人人澡人人爽秒播| 桃色一区二区三区在线观看| 久久香蕉精品热| 午夜福利欧美成人| 国产精品久久久久久久久免 |