• <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
    灌水秘笈
    亚洲国产高清在线一区二区三| av免费观看日本| 亚洲av免费在线观看| 久久99热这里只频精品6学生 | 岛国毛片在线播放| 国产 一区 欧美 日韩| 尾随美女入室| 久久99精品国语久久久| 九九热线精品视视频播放| 久久久久网色| 联通29元200g的流量卡| 国产高潮美女av| 最近最新中文字幕免费大全7| 超碰av人人做人人爽久久| 国产成人aa在线观看| 亚洲三级黄色毛片| 亚洲中文字幕一区二区三区有码在线看| www.色视频.com| 视频中文字幕在线观看| 99久国产av精品| 午夜视频国产福利| 男女啪啪激烈高潮av片| 国产精品爽爽va在线观看网站| 水蜜桃什么品种好| 91精品一卡2卡3卡4卡| 麻豆久久精品国产亚洲av| 国产成人免费观看mmmm| 亚洲四区av| 老师上课跳d突然被开到最大视频| 色综合色国产| 久久久成人免费电影| 欧美zozozo另类| 一级av片app| 特级一级黄色大片| 免费观看a级毛片全部| 免费看美女性在线毛片视频| 久久精品国产99精品国产亚洲性色| 欧美日本视频| 久久国产乱子免费精品| 99久久无色码亚洲精品果冻| 亚洲图色成人| 国产白丝娇喘喷水9色精品| 精品久久久久久久久亚洲| 青青草视频在线视频观看| 99热这里只有是精品在线观看| 中文字幕av在线有码专区| 两个人的视频大全免费| 久久久久久大精品| 婷婷色av中文字幕| 成年版毛片免费区| 我的女老师完整版在线观看| 大香蕉97超碰在线| 美女被艹到高潮喷水动态| 亚洲av日韩在线播放| 长腿黑丝高跟| 免费在线观看成人毛片| 黄色配什么色好看| 亚洲综合色惰| 中文乱码字字幕精品一区二区三区 | 简卡轻食公司| 亚洲欧美成人精品一区二区| 亚洲成人精品中文字幕电影| 岛国在线免费视频观看| 中文在线观看免费www的网站| 波野结衣二区三区在线| 国产欧美另类精品又又久久亚洲欧美| 成人性生交大片免费视频hd| 欧美日本视频| 搡女人真爽免费视频火全软件| 日日摸夜夜添夜夜爱| 老司机福利观看| 黄色配什么色好看| 国产精品国产三级国产av玫瑰| 国语对白做爰xxxⅹ性视频网站| 午夜久久久久精精品| 欧美日韩综合久久久久久| 国产精品久久久久久精品电影| 六月丁香七月| 乱码一卡2卡4卡精品| 日韩欧美在线乱码| 亚洲av二区三区四区| 在线a可以看的网站| 天天躁夜夜躁狠狠久久av| 97人妻精品一区二区三区麻豆| 2021天堂中文幕一二区在线观| 精品人妻偷拍中文字幕| 纵有疾风起免费观看全集完整版 | 日韩在线高清观看一区二区三区| 国产精品乱码一区二三区的特点| 国语自产精品视频在线第100页| 永久免费av网站大全| 亚洲成人精品中文字幕电影| 97超视频在线观看视频| 亚洲av一区综合| 午夜精品在线福利| 国产精品久久久久久久久免| 国产精品电影一区二区三区| 国产极品精品免费视频能看的| 五月玫瑰六月丁香| 欧美丝袜亚洲另类| 赤兔流量卡办理| av在线观看视频网站免费| 日本免费a在线| 国产v大片淫在线免费观看| 亚洲色图av天堂| 菩萨蛮人人尽说江南好唐韦庄 | 伦理电影大哥的女人| 亚洲美女视频黄频| 国语对白做爰xxxⅹ性视频网站| 久久久久久久亚洲中文字幕| 国产亚洲5aaaaa淫片| 一夜夜www| 一卡2卡三卡四卡精品乱码亚洲| 国产白丝娇喘喷水9色精品| 精品一区二区免费观看| 精品人妻熟女av久视频| 久久久久久久久久成人| 男人舔女人下体高潮全视频| 在线免费十八禁| 赤兔流量卡办理| 国产真实伦视频高清在线观看| 一区二区三区四区激情视频| 成人综合一区亚洲| 汤姆久久久久久久影院中文字幕 | 一区二区三区乱码不卡18| 国产成人精品一,二区| 午夜福利高清视频| 免费搜索国产男女视频| 国产欧美日韩精品一区二区| 亚洲精品成人久久久久久| 亚洲五月天丁香| 亚洲在线观看片| 午夜日本视频在线| 国产成人a∨麻豆精品| 看黄色毛片网站| 亚洲国产精品sss在线观看| 91精品伊人久久大香线蕉| 久久久a久久爽久久v久久| 男的添女的下面高潮视频| 欧美一区二区亚洲| 国产 一区 欧美 日韩| 国产精品综合久久久久久久免费| 91在线精品国自产拍蜜月| 久久精品国产鲁丝片午夜精品| 亚洲欧美精品自产自拍| 校园人妻丝袜中文字幕| 亚洲国产精品国产精品| 蜜臀久久99精品久久宅男| 少妇猛男粗大的猛烈进出视频 | 国产精品久久电影中文字幕| 久久热精品热| av卡一久久| 联通29元200g的流量卡| 婷婷色综合大香蕉| 国产淫片久久久久久久久| 日韩成人伦理影院| 久久久成人免费电影| 在线免费十八禁| 高清视频免费观看一区二区 | 18禁裸乳无遮挡免费网站照片| 亚洲丝袜综合中文字幕| 久久精品久久久久久噜噜老黄 | 人人妻人人澡欧美一区二区| 18禁裸乳无遮挡免费网站照片| 国产伦精品一区二区三区四那| 国产 一区精品| 亚洲五月天丁香| 欧美潮喷喷水| 97人妻精品一区二区三区麻豆| 亚洲va在线va天堂va国产| 亚洲av中文字字幕乱码综合| 亚洲av男天堂| 婷婷色av中文字幕| 欧美一区二区精品小视频在线| 久久99热这里只频精品6学生 | 男人舔女人下体高潮全视频| 国产人妻一区二区三区在| 亚洲性久久影院| 亚洲综合色惰| 国产v大片淫在线免费观看| 久久精品91蜜桃| 狂野欧美激情性xxxx在线观看| 老司机福利观看| 国产一级毛片七仙女欲春2| 听说在线观看完整版免费高清| 别揉我奶头 嗯啊视频| av在线蜜桃| 18禁在线播放成人免费| 一级爰片在线观看| 日韩精品有码人妻一区| 久久99热6这里只有精品| 小说图片视频综合网站| 免费观看人在逋| 最近最新中文字幕大全电影3| 久久久午夜欧美精品| 好男人视频免费观看在线| 成人欧美大片| www.色视频.com| 99热6这里只有精品| 免费av不卡在线播放| 有码 亚洲区| 一卡2卡三卡四卡精品乱码亚洲| 日本黄大片高清| 日韩成人av中文字幕在线观看| 久久久a久久爽久久v久久| 人人妻人人澡欧美一区二区| 久久久国产成人免费| 日本熟妇午夜| 亚洲av成人av| 别揉我奶头 嗯啊视频| 国产三级中文精品| 精品欧美国产一区二区三| 国产探花极品一区二区| 国产成人精品婷婷| av线在线观看网站| 欧美日韩在线观看h| 2021少妇久久久久久久久久久| 伦精品一区二区三区| 免费不卡的大黄色大毛片视频在线观看 | 久久久国产成人精品二区| 国产av码专区亚洲av| 中文字幕熟女人妻在线| 青春草视频在线免费观看| av免费观看日本| 国产精品久久久久久精品电影小说 | 久久精品综合一区二区三区| av黄色大香蕉| 最新中文字幕久久久久| 狂野欧美白嫩少妇大欣赏| 国产黄片美女视频| 69人妻影院| 国产淫语在线视频| 亚洲精品乱久久久久久| 色尼玛亚洲综合影院| 国产三级中文精品| 在线免费十八禁| h日本视频在线播放| 成人特级av手机在线观看| 国产av一区在线观看免费| 人人妻人人澡欧美一区二区| 国产黄色视频一区二区在线观看 | 国产成人aa在线观看| 国产精品99久久久久久久久| 国产女主播在线喷水免费视频网站 | av在线播放精品| 建设人人有责人人尽责人人享有的 | 国产在线男女| 2022亚洲国产成人精品| 69人妻影院| 尤物成人国产欧美一区二区三区| 亚洲成色77777| 丰满人妻一区二区三区视频av| 国产精品熟女久久久久浪| 成人欧美大片| 欧美又色又爽又黄视频| 中文字幕熟女人妻在线| 欧美另类亚洲清纯唯美| 亚洲中文字幕日韩| 国产高清国产精品国产三级 | 亚洲性久久影院| 黄片无遮挡物在线观看| 校园人妻丝袜中文字幕| 一卡2卡三卡四卡精品乱码亚洲| 亚洲在久久综合| 久久精品人妻少妇| 七月丁香在线播放| 午夜福利在线观看吧| 99久久成人亚洲精品观看| 亚洲精品,欧美精品| 99久国产av精品国产电影| 大香蕉久久网| 中国国产av一级| 在线观看美女被高潮喷水网站| 99热这里只有是精品50| 国产成人91sexporn| 久久久国产成人精品二区| 22中文网久久字幕| 亚洲国产色片| .国产精品久久| 校园人妻丝袜中文字幕| 国产伦一二天堂av在线观看| 中文字幕熟女人妻在线| 午夜免费男女啪啪视频观看| 国产一区二区在线av高清观看| 中文在线观看免费www的网站| 哪个播放器可以免费观看大片| 午夜激情福利司机影院| 日韩av不卡免费在线播放| 国产精华一区二区三区| 国产熟女欧美一区二区| 精品久久久久久久末码| 九草在线视频观看| 免费观看性生交大片5| 看黄色毛片网站| 国产成人精品婷婷| 成人综合一区亚洲| 亚洲熟妇中文字幕五十中出| 国产色爽女视频免费观看| 搡女人真爽免费视频火全软件| 啦啦啦啦在线视频资源| 久久99蜜桃精品久久| 国产精品美女特级片免费视频播放器| 纵有疾风起免费观看全集完整版 | 日韩 亚洲 欧美在线| 亚洲中文字幕日韩| 最近手机中文字幕大全| 国产精品久久久久久久电影| 97热精品久久久久久| 免费看光身美女| 我要看日韩黄色一级片| 看黄色毛片网站| АⅤ资源中文在线天堂| 国产成人精品一,二区| 免费av毛片视频| 中文字幕制服av| 亚洲欧美精品专区久久| 99热这里只有是精品在线观看| ponron亚洲| 午夜精品在线福利| 国产免费一级a男人的天堂| 一区二区三区免费毛片| 最近中文字幕2019免费版| 69av精品久久久久久| 欧美人与善性xxx| 欧美性猛交黑人性爽| 成人一区二区视频在线观看| 色哟哟·www| 乱码一卡2卡4卡精品| 日韩制服骚丝袜av| 国产乱人偷精品视频| 日韩精品有码人妻一区| 欧美高清性xxxxhd video| 日本午夜av视频| 成人鲁丝片一二三区免费| 2022亚洲国产成人精品| 国产精品一二三区在线看| 久久精品人妻少妇| 非洲黑人性xxxx精品又粗又长| 久久久久久九九精品二区国产| 国产精品嫩草影院av在线观看| 在线a可以看的网站| 嫩草影院入口| 夫妻性生交免费视频一级片| 国语对白做爰xxxⅹ性视频网站| 亚洲成色77777| 丝袜美腿在线中文| .国产精品久久| 亚洲高清免费不卡视频| 成人亚洲欧美一区二区av| 好男人视频免费观看在线| h日本视频在线播放| 91久久精品国产一区二区成人| 特级一级黄色大片| 欧美zozozo另类| 少妇熟女aⅴ在线视频| 日本av手机在线免费观看| 中国美白少妇内射xxxbb| 丰满乱子伦码专区| 亚洲欧美成人精品一区二区| 亚洲av一区综合| 麻豆成人av视频| 最近视频中文字幕2019在线8| kizo精华| 搡老妇女老女人老熟妇| 女人久久www免费人成看片 | 国产爱豆传媒在线观看| 在线免费观看不下载黄p国产| 日韩欧美 国产精品| 一本一本综合久久| 永久网站在线| av女优亚洲男人天堂| 男的添女的下面高潮视频| 能在线免费观看的黄片| 日本五十路高清| 国产中年淑女户外野战色| 男人舔奶头视频| 三级毛片av免费| 91午夜精品亚洲一区二区三区| 国产精品美女特级片免费视频播放器| 2022亚洲国产成人精品| 啦啦啦啦在线视频资源| 欧美bdsm另类| 亚洲欧美精品自产自拍| av女优亚洲男人天堂| 午夜福利高清视频| 99热这里只有是精品在线观看| 国产精品福利在线免费观看| 亚洲av不卡在线观看| eeuss影院久久| 欧美成人午夜免费资源| 亚洲av.av天堂| 国产高清三级在线| 亚洲国产欧美人成| 国产成人一区二区在线| 久久6这里有精品| 国语自产精品视频在线第100页| 18禁在线无遮挡免费观看视频| 在线a可以看的网站| 午夜免费激情av| 亚洲电影在线观看av| 中文乱码字字幕精品一区二区三区 | 久久欧美精品欧美久久欧美| 又爽又黄无遮挡网站| 精品久久久久久久久亚洲| 国产精品久久电影中文字幕| 美女黄网站色视频| 久久这里有精品视频免费| 91久久精品国产一区二区三区| 在线播放无遮挡| 精品免费久久久久久久清纯| 波野结衣二区三区在线| 老司机影院成人| 国内少妇人妻偷人精品xxx网站| 能在线免费观看的黄片| 日韩精品青青久久久久久| 麻豆乱淫一区二区| 成人三级黄色视频| 久久精品国产鲁丝片午夜精品| 国产精品,欧美在线| 中文字幕久久专区| 亚洲av电影在线观看一区二区三区 | 啦啦啦啦在线视频资源| 亚洲欧美成人综合另类久久久 | 亚洲最大成人av| 日日摸夜夜添夜夜添av毛片| 91在线精品国自产拍蜜月| 网址你懂的国产日韩在线| 日韩欧美精品免费久久| 成人毛片a级毛片在线播放| 国产中年淑女户外野战色| 欧美日韩在线观看h| 99热这里只有精品一区| 国产一区亚洲一区在线观看| 人人妻人人看人人澡| 天美传媒精品一区二区| 少妇熟女aⅴ在线视频| 18禁在线无遮挡免费观看视频| 全区人妻精品视频| 亚洲成人久久爱视频| 国产精品福利在线免费观看| 高清毛片免费看| 久久久久久久久大av| 日韩视频在线欧美| 18禁在线播放成人免费| 国产精品一区二区在线观看99 | 国产精品国产三级国产专区5o | 日韩视频在线欧美| 一级二级三级毛片免费看| 欧美激情久久久久久爽电影| 国产精华一区二区三区| 久久久午夜欧美精品| 91精品国产九色| 国产成人91sexporn| 青春草视频在线免费观看| 亚洲最大成人av| 亚洲国产精品合色在线| 日本一二三区视频观看| 久久精品人妻少妇| 成人亚洲精品av一区二区| 国产私拍福利视频在线观看| 成人欧美大片| 嫩草影院新地址| 精品不卡国产一区二区三区| 在线天堂最新版资源| 一本一本综合久久| 国产精品福利在线免费观看| 亚洲欧美精品综合久久99| 国产私拍福利视频在线观看| 欧美变态另类bdsm刘玥| 国产高清有码在线观看视频| 你懂的网址亚洲精品在线观看 | 看非洲黑人一级黄片| 国产综合懂色| 噜噜噜噜噜久久久久久91| 亚洲欧美日韩东京热| 欧美成人一区二区免费高清观看| 日本免费一区二区三区高清不卡| 少妇熟女欧美另类| 久久精品夜夜夜夜夜久久蜜豆| 亚洲av成人精品一区久久| 91久久精品国产一区二区三区| 少妇熟女aⅴ在线视频| 最近手机中文字幕大全| 国产精品国产三级专区第一集| 亚洲av电影不卡..在线观看| 欧美精品国产亚洲| 久久久久免费精品人妻一区二区| av国产久精品久网站免费入址| 91久久精品电影网| 18禁动态无遮挡网站| 亚州av有码| 亚洲成av人片在线播放无| 国产精品久久久久久av不卡| 国产精品综合久久久久久久免费| 极品教师在线视频| 国产精品人妻久久久影院| 久久久欧美国产精品| 亚洲av熟女| 国内精品美女久久久久久| 国产精品一二三区在线看| 菩萨蛮人人尽说江南好唐韦庄 | 日韩三级伦理在线观看| 91狼人影院| videossex国产| 日韩一本色道免费dvd| 特大巨黑吊av在线直播| 长腿黑丝高跟| 免费看光身美女| 大香蕉97超碰在线| 国产精品无大码| 天堂影院成人在线观看| 中文字幕熟女人妻在线| 成人亚洲欧美一区二区av| 最近的中文字幕免费完整| 日韩精品有码人妻一区| 日韩一区二区视频免费看| 女人十人毛片免费观看3o分钟| 国产免费福利视频在线观看| 国产单亲对白刺激| 国产精品人妻久久久影院| 国产伦在线观看视频一区| 一级毛片我不卡| 91久久精品电影网| 成年av动漫网址| 精品久久久久久久久av| 欧美不卡视频在线免费观看| 国产老妇伦熟女老妇高清| 亚洲国产精品sss在线观看| 又黄又爽又刺激的免费视频.| 美女xxoo啪啪120秒动态图| 在线天堂最新版资源| 国产视频内射| 爱豆传媒免费全集在线观看| 欧美+日韩+精品| 亚洲欧美清纯卡通| 精品人妻视频免费看| 成人美女网站在线观看视频| 最近中文字幕高清免费大全6| 大话2 男鬼变身卡| 亚洲av不卡在线观看| 亚洲内射少妇av| 岛国毛片在线播放| 午夜爱爱视频在线播放| 亚洲美女视频黄频| 18+在线观看网站| 97热精品久久久久久| 日本五十路高清| 最近最新中文字幕大全电影3| 欧美一区二区亚洲| 国产免费福利视频在线观看| 女的被弄到高潮叫床怎么办| 国产国拍精品亚洲av在线观看| 黄色一级大片看看| 水蜜桃什么品种好| 菩萨蛮人人尽说江南好唐韦庄 | 国产成人福利小说| h日本视频在线播放| 秋霞在线观看毛片| 国产私拍福利视频在线观看| 欧美激情国产日韩精品一区| 好男人在线观看高清免费视频| 日本午夜av视频| 国产欧美日韩精品一区二区| 亚洲婷婷狠狠爱综合网| 只有这里有精品99| 成人美女网站在线观看视频| 99国产精品一区二区蜜桃av| 亚洲精品亚洲一区二区| 中文在线观看免费www的网站| 少妇裸体淫交视频免费看高清| 寂寞人妻少妇视频99o| 晚上一个人看的免费电影| 亚州av有码| 九草在线视频观看| 狂野欧美白嫩少妇大欣赏| 日日撸夜夜添| 熟妇人妻久久中文字幕3abv| 国产精品久久久久久精品电影小说 | 麻豆精品久久久久久蜜桃| 美女黄网站色视频| 久久久久国产网址| 国产黄色视频一区二区在线观看 | 久久久久久久久中文| 久久久久九九精品影院| av国产免费在线观看| 免费人成在线观看视频色| 亚洲高清免费不卡视频| 国产精品福利在线免费观看| 亚洲三级黄色毛片| 美女内射精品一级片tv| 亚洲乱码一区二区免费版| 欧美97在线视频| 嫩草影院入口| 亚洲乱码一区二区免费版| 联通29元200g的流量卡| 国产精品久久视频播放| 国产一区二区三区av在线| 成人欧美大片| 一区二区三区乱码不卡18| 日本午夜av视频| 国产午夜精品一二区理论片| 禁无遮挡网站| 日韩精品有码人妻一区| 尾随美女入室| 国产片特级美女逼逼视频| 国产精品久久久久久精品电影| 搞女人的毛片| 人妻少妇偷人精品九色| 色吧在线观看| 搞女人的毛片| 三级毛片av免费| 亚洲自偷自拍三级| h日本视频在线播放| 亚洲一区高清亚洲精品| 国产av不卡久久|