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

    利用日光誘導(dǎo)葉綠素?zé)晒夤浪愀邷馗珊祻?fù)合脅迫下夏玉米生理生態(tài)參數(shù)

    2022-07-16 11:47:23倪鋒謝鵬飛褚榮浩李萌沙修竹彭金龍申雙和
    江蘇農(nóng)業(yè)學(xué)報 2022年3期
    關(guān)鍵詞:夏玉米

    倪鋒 謝鵬飛 褚榮浩 李萌 沙修竹 彭金龍 申雙和

    摘要: 高溫和干旱會對農(nóng)業(yè)生產(chǎn)活動及陸地生態(tài)系統(tǒng)產(chǎn)生嚴(yán)重威脅,并最終對經(jīng)濟(jì)和社會發(fā)展產(chǎn)生不利影響。本研究以夏玉米為研究對象,開展高溫干旱復(fù)合脅迫田間試驗,首先對比分析各脅迫處理下夏玉米形態(tài)和生理生態(tài)響應(yīng)特征,之后引入日光誘導(dǎo)葉綠素?zé)晒猓⊿IF),探究各脅迫處理下夏玉米冠層SIF值日變化和不同生育階段變化特征,最終揭示SIF值與生理生態(tài)參數(shù)之間的關(guān)聯(lián)性。結(jié)果表明,高溫、干旱以及高溫干旱復(fù)合脅迫均會抑制夏玉米株高和葉面積,造成產(chǎn)量下降。在生理生態(tài)參數(shù)方面,各脅迫處理下夏玉米葉片含水量( LWC )和葉綠素含量的變化規(guī)律基本一致,隨各脅迫程度的加重總體均呈現(xiàn)逐漸下降趨勢。夏玉米冠層SIF值在日尺度上呈現(xiàn)“單峰”型變化特征,在不同生育階段呈現(xiàn)逐漸下降趨勢?;?/3的冠層SIF值與對應(yīng) LWC 、葉綠素含量之間構(gòu)建的線性模型的 R? 2 值分別為0.837和0.509( P < 0.05);且基于上述線性模型和1/3的冠層SIF值估算的 LWC 和葉綠素含量估算值與實測值之間擬合的 R? 2 值分別為0.827和0.726( P < 0.05),驗證了利用SIF值估算高溫干旱復(fù)合脅迫下夏玉米生理生態(tài)參數(shù)( LWC 和葉綠素含量)的可靠性。

    關(guān)鍵詞: 夏玉米; 高溫干旱復(fù)合脅迫; 生理生態(tài)參數(shù); 日光誘導(dǎo)葉綠素?zé)晒猓⊿IF)

    中圖分類號: S166;S514?? 文獻(xiàn)標(biāo)識碼: A?? 文章編號: 1000-4440(2022)03-0587-10

    Estimation of physiological and ecological parameters of summer maize under combined stress of heat and drought using solar-induced chlorophyll fluorescence

    NI Feng 1 , XIE Peng-fei 1 , CHU Rong-hao 2 , LI Meng 1 , SHA Xiu-zhu 3 , PENG Jin-long 1 , SHEN Shuang-he 4

    (1.School of Resources and Environment, Anhui Agricultural University, Hefei 230036, China; 2.Anhui Public Meteorological Service Center/Anhui Meteorological Bureau, Hefei 230031, China; 3.Weather Modification Center of Henan Province, Zhengzhou 450003, China; 4.School of Applied Meteorology, Nanjing University of Information Science and Technology, Nanjing 210044, China)

    Abstract: Heat and drought will pose a serious threat to agricultural activities and terrestrial ecosystems, and have a negative impact on economic and social development ultimately. In this study, summer maize was taken as the research object to carry out field experiments under combined stress of heat and drought. Firstly, the morphological and ecophysiological response characteristics of summer maize under various stresses were compared and analyzed. Then solar-induced chlorophyll fluorescence (SIF) was introduced to investigate the diurnal variation of SIF values and the variation characteristics at different growth stages of summer maize under different stress treatments. Finally, the correlation between SIF values and physiological and ecological parameters was revealed. The results showed that heat, drought and combined stress of heat and drought inhibited plant height, leaf area and yield of summer maize. In terms of physiological and ecological parameters, the changes of leaf water content ( LWC ) and chlorophyll content of summer maize under different stress treatments were basically consistent, presenting a gradually downward trend on the whole. On daily scale, the canopy SIF values of summer maize showed a "single peak" variation, and gradually decreased at different growth stages. In the correlation analysis, the ?R? 2 ?values fitting between 2/3 canopy SIF values and ?LWC ?and chlorophyll content were 0.837 and 0.509, respectively( P < 0.05). The ?R? 2 ?values fitting between the measured value and the estimated value of the remaining 1/3 ?LWC ?and chlorophyll content were 0.827 and 0.726, respectively( P < 0.05), which further verified the reliability of estimating physiological and ecological parameters ( LWC ?and chlorophyll content) of summer maize under combined stress of heat and drought based? on SIF value.

    Key words: summer maize; combined stress of heat and drought; physiological and ecological parameters; solar-induced chlorophyll fluorescence(SIF)

    近年來,全球氣候變暖加劇,干旱已成為備受關(guān)注的氣候問題 [1-2] 。準(zhǔn)確及時地監(jiān)測大規(guī)模干旱事件,對研究植被變化及其對氣候變化的響應(yīng)、保證糧食安全具有重要意義 [3] 。此外,溫度脅迫是限制作物產(chǎn)量的最重要非生物脅迫因子之一。在全球氣候變化的背景下,極端高、低溫事件發(fā)生的頻率和強度均有所上升。而玉米作為全球最重要的糧食作物之一,其在生長發(fā)育過程中對環(huán)境脅迫較為敏感,其中高溫與干旱脅迫均是玉米生長過程中常見的環(huán)境脅迫因子,二者均會對玉米產(chǎn)量和品質(zhì)造成不同程度的影響。

    在有關(guān)高溫干旱脅迫研究中,方子松等 [4] 研究發(fā)現(xiàn),在干旱脅迫下,薏苡的株高和莖粗會受到不同程度的抑制。莊克章等 [5] 發(fā)現(xiàn)玉米的株高在干旱脅迫下會呈現(xiàn)明顯的下降趨勢,且脅迫越嚴(yán)重下降幅度越大。崔麗娜等 [6] 研究發(fā)現(xiàn),在高溫脅迫下,夏玉米的生育期會明顯縮短,其產(chǎn)量也會下降。在生理生態(tài)參數(shù)研究中,邵宇航等 [7] 和費立偉 [8] 發(fā)現(xiàn)小麥葉片葉綠素含量在高溫脅迫下會有所降低。郝召君等 [9] 研究指出,受高溫脅迫時,芍藥的葉綠素含量會明顯下降。

    傳統(tǒng)植被指數(shù)[如歸一化差異植被指數(shù)( NDVI )和增強植被指數(shù)( EVI )]已被廣泛用來監(jiān)測植被的生理狀態(tài)和水分脅迫對植被的影響。其雖然能夠表征植被的生理狀態(tài),但只能反映植被色素含量的變化,與植物的實際光合作用無直接關(guān)聯(lián)。因此,傳統(tǒng)植被指數(shù)對植被水分缺失情況反映較慢,不能即時反映植被對水分脅迫的動態(tài)響應(yīng)過程。葉綠素?zé)晒馐侨~綠素分子在吸收光能后,在極短的時間內(nèi)重新釋放的紅光和近紅外波段的光子 [10] 。由于獲取手段的差異性,葉綠素?zé)晒庵饕ㄖ鲃雍捅粍诱T導(dǎo)的葉綠素?zé)晒?[11] 。其中,被動誘導(dǎo)葉綠素?zé)晒馐且蕴柟鉃榧ぐl(fā)光源,在自然條件下獲取的葉綠素?zé)晒?,即日光誘導(dǎo)葉綠素?zé)晒猓⊿olar-induced chlorophyll fluorescence, SIF)。由于葉綠素?zé)晒獍l(fā)射量與光合作用效率直接相關(guān),可以提供關(guān)于光化學(xué)反應(yīng)和熱耗散量子效率的有價值的信息 [12] 。因此,SIF可以作為植被瞬時光合作用活性的理想探針,在植物環(huán)境脅迫監(jiān)測研究中應(yīng)用 [13] 。

    Dobrowski等 [14] 研究發(fā)現(xiàn),在紅邊光譜區(qū)域計算的簡單反射率指數(shù)可以追蹤水分脅迫引起的SIF變化,從而提供與植物生理過程直接相關(guān)的作物水分脅迫的快速評估。Alexander等 [15] 研究指出,在葉片和冠層尺度上,SIF信號強度的下降與水分脅迫有關(guān)。Daumard等 [16] 通過對SIF的連續(xù)地面觀測進(jìn)一步指出,SIF可以較早地檢測到高粱田發(fā)生的可逆水分脅迫,且從航空圖像SIF信號中可以檢測高粱在水分脅迫條件下發(fā)生的生理變化。Liu等 [17] 和Lee等 [18] 研究發(fā)現(xiàn),SIF和傳統(tǒng)植被指數(shù)在干旱脅迫下會同時降低,但在不同干旱脅迫程度下二者的遞減趨勢不一致。Wang等 [19] 研究指出,衛(wèi)星SIF在準(zhǔn)確、及時監(jiān)測干旱和高溫?zé)崂朔矫婢哂泻艽蟮陌l(fā)展?jié)摿?。Song等 [20] 利用衛(wèi)星SIF檢測印度恒河平原冬小麥對熱脅迫的早期響應(yīng)。

    目前大多數(shù)研究僅考慮作物受單一高溫或干旱脅迫,然而在現(xiàn)實情況下,高溫和干旱往往會相伴而生,因此有必要探究高溫干旱復(fù)合脅迫下夏玉米形態(tài)和生理生態(tài)參數(shù)響應(yīng)特征,并揭示SIF與生理生態(tài)參數(shù)之間的關(guān)聯(lián)性。本研究擬通過設(shè)計夏玉米不同等級高溫干旱復(fù)合脅迫田間試驗,在觀測數(shù)據(jù)的基礎(chǔ)上,探究高溫干旱復(fù)合脅迫下夏玉米形態(tài)(株高、葉面積和產(chǎn)量)、生理生態(tài)指標(biāo)[葉片含水量( LWC )、葉綠素含量]和冠層SIF響應(yīng)特征,對比分析不同脅迫處理下夏玉米形態(tài)指標(biāo)和生理生態(tài)指標(biāo)之間的差異,之后對比分析不同脅迫處理下冠層SIF響應(yīng)差異,并將冠層SIF與各生理生態(tài)指標(biāo)進(jìn)行回歸分析,建立回歸模型并驗證,揭示夏玉米冠層SIF在高溫干旱復(fù)合脅迫監(jiān)測應(yīng)用上的可能性和適宜性,以期為利用冠層SIF估算高溫干旱復(fù)合脅迫下夏玉米的生理生態(tài)參數(shù)提供理論依據(jù)。

    1 材料與方法

    1.1 試驗區(qū)概況

    試驗地位于安徽農(nóng)業(yè)大學(xué)農(nóng)萃園農(nóng)業(yè)氣象試驗基地(31°51′N,117°15′E),本試驗中玉米播種方式為穴播,于2020年6月30日完成播種,小區(qū)種植7行×6列,每穴內(nèi)放2~ 3粒玉米種子。之后根據(jù)出苗情況進(jìn)行間苗和補苗,施肥、除草等田間管理視玉米生長情況進(jìn)行。玉米形態(tài)觀測和數(shù)據(jù)采集于2020年7-10月夏玉米生長季進(jìn)行。本試驗采用室外大田種植方式,共設(shè)置12個小區(qū),每個小區(qū)長2.0 m,寬1.8 m。為避免受到自然降水的影響,在試驗田干旱處理小區(qū)上方搭建移動式遮雨棚;為避免受到區(qū)間和邊緣效應(yīng)的影響,小區(qū)與小區(qū)間以及小區(qū)與遮雨棚邊緣之間的間隔距離均為0.5 m。如圖1所示,每個對照組都有一個小區(qū)為重復(fù)處理,最終選取9個小區(qū)進(jìn)行觀測。

    脅迫處理在玉米拔節(jié)后7 d左右進(jìn)行,在此之前各小區(qū)充分供水,即將土壤水分灌溉至田間最大持水量( FC )的80%~ 90%水平,以保證夏玉米前期的正常生長,在拔節(jié)后7 d左右按照 FC 的百分比進(jìn)行灌溉處理(表1)。本試驗共設(shè)置三個等級的灌溉處理,分別為重度干旱(W1, 50%~ 60% ?FC )、輕度干旱(W2, 65%~ 75% ?FC )和充分供水對照組(W3, 80%~ 90% ?FC )。溫度處理共設(shè)置兩個等級,分別為H1和H2,增溫處理主要采用塑料薄膜覆蓋的方式,在有增溫處理的小區(qū)四周架設(shè)鍍鋅鋼管架,隨著夏玉米的生長不斷調(diào)整鋼管架上覆膜的位置和高度。其中,H1處理的覆膜方式為24 h四周半覆蓋,H2處理為24 h四周全覆蓋。本研究分別將W1、W1H1和W1H2,W2、W2H1和W2H2,W3、W3H1和W3H2,W1、W2和W3設(shè)為第一組、第二組、第三組和第四組,其中第一組為高溫和重度干旱復(fù)合脅迫處理,第二組為高溫和輕度干旱復(fù)合脅迫處理,第三組為不同溫度的高溫脅迫處理,第四組為不同水分的干旱脅迫處理。

    1.2 數(shù)據(jù)采集與處理

    1.2.1 夏玉米形態(tài)指標(biāo)與生理生態(tài)參數(shù) 生育期:參照《農(nóng)業(yè)氣象觀測規(guī)范(上卷)》進(jìn)行觀測。株高:每個小區(qū)選取10株有代表性的玉米植株,分別測量其高度。葉面積:在每個處理中選取有代表性的玉米植株,分別采集3片穗上葉、3片穗位葉和3片穗下葉,置于保鮮袋中密封,及時取回室內(nèi),測量玉米葉片的長和寬(最寬處),并計算單葉的葉面積。

    葉面積=長×寬×0.7(1)

    葉片含水量( LWC ):在測完葉片長、寬后,先稱取葉片鮮質(zhì)量,然后烘干。前1 h的烘干溫度設(shè)置在105 ℃左右,對葉片進(jìn)行殺青,之后烘干溫度設(shè)定在80 ℃左右,待葉片質(zhì)量不變時稱取干質(zhì)量。

    LWC=FW-DWFW(2)

    其中 FW 為葉片鮮質(zhì)量(單位:g); DW 為葉片干質(zhì)量(單位:g); LWC 為葉片含水量。

    產(chǎn)量及其構(gòu)成因素:主要包括玉米果穗長、果穗粗、禿尖比、株籽粒質(zhì)量、百粒質(zhì)量和理論產(chǎn)量 [21]

    葉綠素含量:使用SPAD葉綠素儀進(jìn)行測定,每個小區(qū)選取3株,選取每株玉米穗上葉、穗位葉、穗下葉,每個葉片均勻選取5個不同部位,分別測定葉綠素值,去掉異常數(shù)據(jù)后記錄其平均值。

    1.2.2 冠層SIF數(shù)據(jù) 利用美國Ocean Optics公司生產(chǎn)的超高分辨率光譜儀(型號為QE-Pro)進(jìn)行觀測:按照初始積分時間采集一條太陽入射光譜,然后記錄一條暗電流,隨后采集一條冠層反射光光譜,最后重新采集一條太陽入射光譜,至此一個觀測循環(huán)結(jié)束,時間為8∶00- 16∶00 ,每兩次觀測間隔1 h。通過獲得的光譜數(shù)據(jù)計算得到輻射亮度值和反射率。

    由于太陽輻射會在近紅外與紅外波段形成三個明顯的吸收波段(656.3 nm、687.0 nm與760.0 nm),這與SIF發(fā)射的波段互相重疊,因此可選取這三個波段對SIF進(jìn)行研究。Moya等 [22] 和劉良云等 [23] 研究認(rèn)為687.0 nm與760.0 nm波段處的夫瑯和費暗線較其他波段更為明顯,而760.0 nm波段處的夫瑯和費暗線在寬度與深度上均比687.0 nm波段處大一些 [22-23] 。因此,選取760.0 nm波段來研究SIF更為合適,本文也選取該波段對夏玉米冠層SIF進(jìn)行研究。

    目前SIF的獲取方法較多,常見的有FLD算法、3FLD算法、cFLD算法和SFM算法等。Damm等 [24] 和Meroni等 [25] 將FLD算法和SFM算法進(jìn)行對比,發(fā)現(xiàn)在任何噪聲干擾情況下,SFM算法的精度要高于FLD算法。以上幾種算法都需要相應(yīng)的假設(shè)作為前提條件,并不能完全符合植被所處的真實環(huán)境狀況,因此在估算熒光值過程中均會產(chǎn)生一定的誤差 [26-27] ,在使用過程中要根據(jù)具體需求來選擇最優(yōu)算法。本研究最終選取SFM算法作為SIF的提取方法。

    1.2.3 數(shù)據(jù)處理 采用Microsoft Excel 2007對原始數(shù)據(jù)進(jìn)行預(yù)處理,利用SPSS 22.0對各類數(shù)據(jù)進(jìn)行分析及回歸擬合等,采用Origin 2016完成文中各插圖的繪制。

    2 結(jié)果分析

    2.1 夏玉米形態(tài)、產(chǎn)量及其構(gòu)成因素變化特征

    2.1.1 株高變化特征 本研究選取夏玉米開花期T1(8月18日)、灌漿期T2(9月9日)和成熟期T3(9月24日)的觀測數(shù)據(jù),比較不同生育階段夏玉米株高的變化情況。從表2可以看出,隨著生育期的推進(jìn),所有處理小區(qū)夏玉米株高均呈先升后降的變化趨勢。其中,T1至T2階段株高不斷上升的主要原因是由于夏玉米處于營養(yǎng)生長階段,該時期夏玉米株高會一直上升。此階段W3、W3H1、W3H2 3個處理夏玉米株高的生長速度要明顯快于重度和輕度干旱處理。T2至T3階段夏玉米株高呈現(xiàn)下降趨勢,主要是由于該時期是夏玉米果實生長階段,莖、葉等其他器官基本停止生長并開始逐漸衰老。

    T1階段各處理夏玉米株高之間無明顯差異,而T2和T3階段差異明顯,T3階段呈現(xiàn)出W1> W1H1> W1H2、W2> W2H1> W2H2、W3> W3H1> W3H2、W3> W2> W1的變化趨勢,該現(xiàn)象表明干旱與高溫脅迫均會抑制夏玉米株高的生長,重度干旱脅迫影響要大于輕度干旱脅迫,且干旱脅迫影響大于高溫脅迫、高溫干旱復(fù)合脅迫影響大于單一脅迫。但是,T2階段株高表現(xiàn)為W2H2> W2H1> W2和W1> W2,且W2、W2H1和W2H2的標(biāo)準(zhǔn)差較大,分別為±21.0 cm、±23.1 cm和±28.2 cm,這可能是由于人工取樣造成的差異。

    2.1.2 葉面積變化特征 如表3所示,所有處理夏玉米穗上葉的葉面積從T1至T3階段均呈現(xiàn)出先增后減的變化趨勢,穗位葉和穗下葉的葉面積均呈現(xiàn)逐漸增加的變化趨勢。T2至T3階段夏玉米穗上葉葉面積出現(xiàn)下降趨勢,穗位葉和穗下葉呈現(xiàn)上升趨勢,其原因可能是由于穗上葉衰老的速度要快于穗位葉和穗下葉。在大部分脅迫處理下,夏玉米穗上葉、穗位葉、穗下葉3葉的葉面積在T1和T2階段表現(xiàn)為穗位葉> 穗上葉> 穗下葉,在T3階段則表現(xiàn)為穗位葉> 穗下葉> 穗上葉。此外,各處理葉面積之間也存在明顯差異,大致表現(xiàn)為W1> W1H1> W1H2、W2> W2H1> W2H2、W3> W3H1> W3H2、W3> W2> W1的變化趨勢。

    2.1.3 產(chǎn)量及其構(gòu)成因素變化特征 如表4所示,在不同脅迫處理下,夏玉米產(chǎn)量及其構(gòu)成因素均呈現(xiàn)出不同程度的下降趨勢,除禿尖比之外,其他參數(shù)基本呈現(xiàn)W1> W1H1> W1H2、W2> W2H1> W2H2、W3> W3H1> W3H2、W3> W2> W1的變化趨勢,禿尖比呈現(xiàn)出W1H2> W1H1> W1、W2H2> W2H1> W2、W1> W2> W3、W3H2> W3H1> W3的變化趨勢。該現(xiàn)象表明干旱與高溫脅迫均會造成果穗變細(xì)、變短,禿尖比增大,株籽粒質(zhì)量下降,最終導(dǎo)致夏玉米減產(chǎn)。其中,W1、W1H1和W1H2 3個處理的理論產(chǎn)量較對照(W3)分別下降了81.9%、82.5%和84.2%。

    2.2 夏玉米不同生育階段生理生態(tài)參數(shù)變化特征

    2.2.1 ?LWC 變化特征 ?LWC 表示葉片的相對含水量,是指征葉片水分含量的一個最為直接的生理指標(biāo) [28] 。葉片是植物進(jìn)行光合作用積累有機物的地方,其含水量變化將對葉片光合能力產(chǎn)生重要影響 [29] 。如圖2所示,各處理的夏玉米 LWC 在T1至T3階段均呈現(xiàn)出逐漸下降趨勢。其中,T1至T2階段的下降速率較快,T2至T3階段的下降速率較緩慢。各處理 LWC 在T1階段的差異較為明顯,具體表現(xiàn)為W1> W1H1> W1H2、W2> W2H1> W2H2、W3> W3H1> W3H2和W3> W2> W1。在T1和T2階段,6個干旱脅迫處理的 LWC 均明顯低于3個非干旱脅迫處理;在T3階段,各處理間 LWC 差異變小。

    2.2.2 葉綠素含量變化特征 葉綠素作為植物進(jìn)行光合作用的主要色素,其在光能的吸收及轉(zhuǎn)化過程中起著重要作用 [28] 。葉綠素含量的多少可以反映葉片的光合效率、作物生長速度、產(chǎn)量的高低以及受脅迫狀態(tài) [28] 。SPAD值與葉片葉綠素含量的相關(guān)性極顯著,并且在測量過程中既不會破壞植物葉片又方便快捷 [30] ,因此葉綠素含量可以用SPAD值進(jìn)行表示。如圖3所示,各處理夏玉米葉片的葉綠素含量隨生育期的推進(jìn)均呈現(xiàn)逐漸下降的變化趨勢。相同生育階段,各處理之間差異明顯,具體表現(xiàn)為W1> W1H1> W1H2、W2> W2H1> W2H2、W3> W3H1> W3H2和W3> W2> W1,表明干旱脅迫與高溫脅迫抑制夏玉米生長。

    2.3 夏玉米冠層SIF值變化特征

    2.3.1 日變化特征 本研究以8月28日(抽絲期)觀測數(shù)據(jù)為例,比較不同脅迫處理冠層SIF日變化情況。如圖4所示,所有處理的夏玉米冠層SIF值日變化曲線均呈現(xiàn)“單峰”型變化特征:8∶00- 12∶00 ?SIF值逐漸升高,并在12∶00 左右達(dá)到峰值。這主要是由于此時間段內(nèi)太陽光逐漸增強,夏玉米吸收的太陽輻射逐漸增加;12∶00- 16∶00 夏玉米冠層SIF值逐漸下降,這主要是由于該時間段內(nèi)太陽光逐漸減弱,夏玉米吸收的太陽輻射逐漸減少,最終導(dǎo)致冠層SIF值降低,并在16∶00 左右達(dá)到極小值。在全天大部分時間內(nèi),各處理間的SIF值呈現(xiàn)出W1> W1H1> W1H2、W2> W2H1> W2H2、W3> W3H1> W3H2、W3> W2> W1的變化趨勢。

    2.3.2 不同生育階段變化特征

    冠層SIF值均為各生育階段中午12∶00 的數(shù)值。如圖5所示,隨著生育期的推進(jìn),夏玉米冠層SIF值呈現(xiàn)出T1> T2> T3的變化趨勢,這可能是由于隨著夏玉米的生長,葉綠素含量逐漸下降、冠層結(jié)構(gòu)發(fā)生改變等原因造成。各處理間差異明顯,總體上表現(xiàn)為重度干旱脅迫的影響程度要大于輕度干旱脅迫,高溫干旱復(fù)合脅迫影響比單一脅迫影響嚴(yán)重。

    2.4 夏玉米冠層SIF值與 LWC 和葉綠素含量之間的相關(guān)性分析

    2.4.1 模型建立 ?LWC 是反應(yīng)植物水分虧缺狀況的有效指標(biāo) [28] ,其數(shù)值在高溫、干旱脅迫下呈現(xiàn)下降趨勢。因此,研究冠層SIF值與 LWC 之間的相關(guān)性能夠驗證在高溫干旱復(fù)合脅迫下利用SIF值估算夏玉米生理生態(tài)參數(shù)的可行性。本研究將T1至T3階段觀測得到的2/3夏玉米冠層SIF值與 LWC 進(jìn)行線性擬合。如圖6所示,夏玉米 LWC 隨冠層SIF值的增加而增加,冠層SIF值與 LWC 之間線性擬合的 R? 2 值為0.837(該值為調(diào)整后的值,消除了自變量之間的影響, P < 0.05),散點分布總體上比較集中,表明冠層SIF值與 LWC 之間的擬合度較好且二者之間的敏感性較高。此外,在冠層SIF值較高時, LWC 持續(xù)升高,分布相對比較分散;而當(dāng)冠層SIF值較低時, LWC 分布比較集中,其中冠層SIF值集中在1.4~ 1.6 μW/(cm 2 ·nm·sr) 時, LWC 浮動較小。上述結(jié)果表明,當(dāng)冠層SIF值達(dá)到一定數(shù)值時,可以較好地反映植被生理狀態(tài)變化。

    葉綠素含量的多少可以反映出葉片的光合效率、作物生長速度、產(chǎn)量的高低以及受脅迫狀態(tài) [28] 。因此,研究夏玉米冠層SIF值與葉綠素含量之間的相關(guān)性,可以為探究SIF值與作物光合作用之間的關(guān)聯(lián)機制提供理論依據(jù)。如圖7所示,夏玉米葉綠素含量隨冠層SIF值的增加而增加,且二者之間線性擬合的 R? 2 值為0.509(該值為調(diào)整后的值,消除了自變量之間的影響, P < 0.05)。此外,二者之間的散點分布總體上比較集中,說明冠層SIF值對葉綠素含量的敏感性較高。在冠層SIF值較高時,葉綠素含量持續(xù)升高,分布相對比較集中;而在冠層SIF值較低時,葉綠素含量分布比較分散,其中SIF值在1.0~ 1.5 μW/(cm 2 ·nm·sr) 時,葉綠素含量波動較大。上述結(jié)果表明,當(dāng)冠層SIF值不是很低時,可以反映植物生理狀態(tài)變化。

    2.4.2 模型驗證 為了驗證上述線性擬合模型的精度,本研究利用剩余1/3的 LWC 和SPAD實測值和通過線性模型得到的對應(yīng)擬合值,建立線性回歸模型進(jìn)行相關(guān)分析。如圖8所示, LWC 和SPAD實測值與擬合值之間的決定系數(shù)分別為0.827和0.726,表明SIF值與 LWC 和葉綠素含量之間建立的線性模型模擬效果較好,即利用SIF值估算高溫干旱復(fù)合脅迫下的 LWC 和葉綠素含量均具有一定的可靠性。

    3 結(jié)論與討論

    本研究基于2020年7-9月間夏玉米高溫干旱復(fù)合脅迫田間試驗,對高溫干旱復(fù)合脅迫下夏玉米生理生態(tài)響應(yīng)進(jìn)行深入分析,同時基于SIF反演算法探究了高溫干旱復(fù)合脅迫下夏玉米冠層SIF的響應(yīng)特征,最終揭示利用冠層SIF值估算高溫干旱復(fù)合脅迫下生理生態(tài)參數(shù)的可行性,主要結(jié)論如下:(1)從不同生育階段看,各處理株高均呈先升后降變化趨勢。開花期株高差異較小,灌漿期和成熟期差異較大。所有處理穗上葉葉面積先增后減,穗位葉和穗下葉葉面積呈逐漸增加趨勢。在產(chǎn)量形成方面,隨干旱脅迫和高溫脅迫程度的加重,除禿尖比外,其他參數(shù)總體呈現(xiàn)逐漸減少的變化趨勢。干旱與高溫脅迫均會抑制夏玉米株高、葉面積和產(chǎn)量,其中重度干旱脅迫影響大于輕度干旱脅迫,干旱脅迫影響大于高溫脅迫,高溫干旱復(fù)合脅迫較單一脅迫影響嚴(yán)重。(2)不同生育階段夏玉米 LWC 和葉綠素含量的變化規(guī)律基本一致,總體呈現(xiàn)逐漸下降的變化趨勢,其中 LWC 在T1至T2階段下降速度較快,在T2至T3階段下降速度較緩慢。干旱與高溫脅迫均會對 LWC 和葉綠素含量產(chǎn)生一定影響,其中重度干旱影響大于輕度干旱,高溫干旱復(fù)合脅迫影響重于單一脅迫,且隨著脅迫加劇,影響越嚴(yán)重。(3)在日尺度上,夏玉米冠層SIF值日變化曲線呈現(xiàn)“單峰”型變化特征;從不同生育階段看,冠層SIF值呈現(xiàn)T1>T2>T3的變化趨勢。干旱與高溫脅迫均會抑制夏玉米冠層SIF值,其中重度干旱脅迫影響大于輕度干旱脅迫,高溫干旱復(fù)合脅迫影響要重于單一脅迫。冠層SIF值可以及時準(zhǔn)確地監(jiān)測到高溫干旱脅迫信息,脅迫越重,冠層SIF值響應(yīng)越劇烈。(4)夏玉米 LWC 和葉綠素含量隨SIF值的增加而增加,SIF值與 LWC 和葉綠素含量之間建立的線性模型模擬效果較好,因此利用SIF值估算高溫干旱復(fù)合脅迫下的 LWC 和葉綠素含量均具有一定的可靠性。

    本研究發(fā)現(xiàn)干旱與高溫脅迫均會抑制夏玉米株高、葉面積和產(chǎn)量,其中重度干旱脅迫影響大于輕度干旱脅迫,干旱脅迫影響大于高溫脅迫,高溫干旱復(fù)合脅迫較單一脅迫影響嚴(yán)重。這一結(jié)果與劉紅芳等 [31] 的研究結(jié)論基本一致,其研究發(fā)現(xiàn)玉米的株高、葉面積、莖粗等形態(tài)指標(biāo)會受到干旱脅迫的影響。不同的是,本研究中按照梯度設(shè)計高溫干旱復(fù)合脅迫,試驗結(jié)果更加顯著。宋晉輝等 [32] 基于設(shè)計的干旱梯度試驗研究發(fā)現(xiàn),干旱脅迫會阻礙玉米株高和葉面積增長,致使產(chǎn)量與干物質(zhì)下降且干旱脅迫越嚴(yán)重,影響越大。

    干旱與高溫脅迫均會對夏玉米 LWC 和葉綠素含量產(chǎn)生影響,且隨著脅迫的加劇,影響越嚴(yán)重。這一結(jié)果與宋賀等 [33] 的研究結(jié)果相似,其研究發(fā)現(xiàn)干旱會影響植物葉綠素的生物合成,持續(xù)的干旱降低了葉片的葉綠素含量,其中在持續(xù)輕度干旱中,葉綠素含量有一個上升的趨勢,而持續(xù)的重度干旱較輕度干旱影響更大,葉綠素含量呈持續(xù)下降趨勢。SIF本身由葉綠素發(fā)射,其值對于葉綠素含量的變化更為敏感。印玉明等 [34] 研究結(jié)果表明SIF值在監(jiān)測水稻冠層和單葉葉綠素含量時效果較好,且基于SIF值的輻射傳輸模型也可以用來反演水稻冠層和單葉的葉綠素含量。此外,前人研究發(fā)現(xiàn),SIF信號強度的下降與水分脅迫有關(guān),其可以較早地檢測到作物的可逆水分脅迫,且從航空圖像SIF信號中可以檢測到作物在水分脅迫條件下的生理變化 [15-16] 。綜上,SIF值在一定程度上能夠較好地估算作物 LWC 和葉綠素含量。

    Sun等 [35] 研究發(fā)現(xiàn)SIF值對干旱的發(fā)生較為敏感,可以用于及時監(jiān)測農(nóng)業(yè)生態(tài)系統(tǒng)上發(fā)生的干旱。Chen等 [36] 研究發(fā)現(xiàn),由于降水突然下降,加上異常高溫,蒸發(fā)蒸騰作用使土壤水分迅速枯竭,導(dǎo)致干旱的迅速出現(xiàn)。此外,SIF值對植被結(jié)構(gòu)和生理生化參數(shù)的變化十分敏感,具有極佳的動態(tài)干旱監(jiān)測潛力。本試驗中冠層SIF值可以及時準(zhǔn)確地監(jiān)測到高溫干旱脅迫信息,且脅迫越重,冠層SIF響應(yīng)越劇烈。

    在今后的研究中,還需從以下角度進(jìn)一步探究高溫干旱復(fù)合脅迫下夏玉米生理生態(tài):(1)完善并繼續(xù)開展夏玉米高溫干旱復(fù)合脅迫田間試驗,強化試驗的可重復(fù)性,增強試驗結(jié)果的可信度;(2)引入SIF自動觀測系統(tǒng)。在人工測量SIF時,時間間隔長且誤差較大,因此在條件允許的情況下,可引入SIF自動觀測系統(tǒng),進(jìn)而實現(xiàn)夏玉米冠層SIF值的實時連續(xù)觀測,以便更好地提高時間精度、減少測量誤差。(3)本文沒有開展葉片尺度的SIF值觀測,在今后的研究中,可同步開展葉片和冠層尺度SIF值觀測研究,為尺度效應(yīng)相關(guān)研究提供理論支撐。

    參考文獻(xiàn):

    [1] 郭艷春. 干旱災(zāi)害風(fēng)險評估及管理對策探析[J]. 地下水, 2019, 44(4): 148-149.

    [2] 紀(jì)瑞鵬,于文穎,馮 銳,等. 作物對干旱脅迫的響應(yīng)過程與早期識別技術(shù)研究進(jìn)展[J]. 災(zāi)害學(xué), 2019, 34(2): 153-160.

    [3] 胡鵬飛. 基于多源遙感數(shù)據(jù)的農(nóng)業(yè)旱情監(jiān)測方法比較研究[D]. 蘭州:西北師范大學(xué), 2019.

    [4] 方子松,龍 華,金隆英,等. 干旱脅迫對云南薏苡新品種形態(tài)和生理指標(biāo)的影響[J]. 農(nóng)業(yè)科學(xué)與技術(shù)(英文版), 2019, 20(1):22-29.

    [5] 莊克章,吳榮華,張春艷,等. 苗期干旱及復(fù)水對玉米生長和生理特性的影響[J]. 山東農(nóng)業(yè)科學(xué), 2020, 52(350): 61-66.

    [6] 崔麗娜,董樹亭. 不同氮肥處理下高溫脅迫對夏玉米產(chǎn)量及葉片超微構(gòu)造的影響[J]. 玉米科學(xué), 2020, 28(137): 96-101.

    [7] 邵宇航,石祖梁,張 姍,等. 高溫脅迫下鎂對小麥旗葉光合特性及產(chǎn)量的影響[J]. 麥類作物學(xué)報, 2018, 38(7): 802-808.

    [8] 費立偉. 晚播對冬小麥灌漿后期高溫脅迫下光合能力和產(chǎn)量的影響[D]. 泰安:山東農(nóng)業(yè)大學(xué), 2020.

    [9] 郝召君,周春華,劉 定,等. 高溫脅迫對芍藥光合作用、葉綠素?zé)晒馓匦约俺⒔Y(jié)構(gòu)的影響[J]. 分子植物育種, 2017,15(6):2359-2367.

    [10] 關(guān)琳琳. 基于葉綠素?zé)晒獾闹脖豢偝跫壣a(chǎn)力估算[D]. 北京:中國科學(xué)院大學(xué), 2017.

    [11] 劉雷震,武建軍,周洪奎,等. 葉綠素?zé)晒饧捌湓谒置{迫監(jiān)測中的研究進(jìn)展[J]. 光譜學(xué)與光譜分析, 2017, 37(9): 2780-2787.

    [12] 程占慧,劉良云. 基于葉綠素?zé)晒獍l(fā)射光譜的光能利用率探測[J]. 農(nóng)業(yè)工程學(xué)報, 2010, 26(S2): 74-80.

    [13] 詹春暉,章釗穎,張永光. 日光誘導(dǎo)葉綠素?zé)晒廨椛鋫鬏斈P脱芯窟M(jìn)展[J]. 遙感學(xué)報, 2020, 24(8): 945-957.

    [14] DOBROWSKI S Z, PUSHNIK J C, ZARCO-TEJADA P J, et al. Simple reflectance indices track heat and water stress-induced changes in steady-state chlorophyll fluorescence at the canopy scale[J]. Remote Sensing of Environment, 2005, 97: 403-414.

    [15] ALEXANDER A, MALENOVSKY Z, OLEJNICKOVA J, et al. Meta-analysis assessing potential of steady-state chlorophyll fluorescence for remote sensing detection of plant water, temperature and nitrogen stress[J]. Remote Sensing of Environment, 2015, 168: 420-436.

    [16] DAUMARD F, CHAMPAGNE S, FOURNIER A, et al. A field platform for continuous measurement of canopy fluorescence[J]. IEEE Transactions on Geoscience and Remote Sensing, 2010, 48: 3358-3368.

    [17] LIU L, YANG X, ZHOU H, et al. Evaluating the utility of solar-induced chlorophyll fluorescence for drought monitoring by comparison with NDVI derived from wheat canopy[J]. Science of the Total Environment, 2018, 625: 1208-1217.

    [18] LEE J E, FRANKENBERG C, VANDERTOL C, et al. Forest productivity and water stress in Amazonia: observations from GOSAT chlorophyll fluorescence[J]. Proceedings of the Royal Society B(Biological Sciences), 2013, 280: 20130171.

    [19] WANG X, QIU B, LI W, et al. Impacts of drought and heatwave on the terrestrial ecosystem in China as revealed by satellite solar-induced chlorophyll fluorescence[J]. Science of the Total Environment, 2019, 693: 133627.

    [20] SONG L, LUIS G, GUAN K, et al. Satellite sun-induced chlorophyll fluorescence detects early response of winter wheat to heat stress in the Indian Indo-Gangetic Plains[J]. Global Change Biology, 2018, 24(9): 4023-4037.

    [21] 王江山. 生態(tài)與農(nóng)業(yè)氣象[M]. 北京:氣象出版社, 2005.

    [22] MOYA I, CAMENEN L, EVAIN S, et al. A new instrument for passive remote sensing[J]. Remote Sensing of Environment, 2004, 91: 186-197.

    [23] 劉良云,張永江,王紀(jì)華,等. 利用夫瑯和費暗線探測自然光條件下的植被光合作用熒光研究[J]. 遙感學(xué)報, 2006, 10(1): 130-137.

    [24] DAMM A, ERLER A, HILLEN W, et al. Modeling the impact of spectral sensor configurations on the FLD retrieval accuracy of sun-induced chlorophyll fluorescence[J]. Remote Sensing of Environment, 2011, 115: 1882-1892.

    [25] MERONI M, BUSETTO L, COLOMBO R, et al. Performance of spectral fitting methods for vegetation fluorescence quantification[J]. Remote Sensing of Environment, 2010, 114: 363-374.

    [26] 問靜怡. 基于日光誘導(dǎo)熒光的水稻銅脅迫遙感監(jiān)測機理研究[D]. 杭州:杭州師范大學(xué), 2019.

    [27] 王 冉,劉志剛,楊沛琦. 植物日光誘導(dǎo)葉綠素?zé)晒獾倪b感原理及研究進(jìn)展[J]. 地球科學(xué)進(jìn)展, 2012, 27(11): 1221-1228.

    [28] 李 萌. 夏玉米理化參數(shù)對連續(xù)水分脅迫的響應(yīng)特征及遙感監(jiān)測[D].南京:南京信息工程大學(xué), 2018.

    [29] 陳秀青,楊 琦,韓景曄, 等. 基于葉冠尺度高光譜的冬小麥葉片含水量估算[J]. 光譜學(xué)與光譜分析, 2020, 40(3): 233-239.

    [30] 王秋玲. 夏玉米生理生態(tài)與生長特性對干旱過程的響應(yīng)研究[D].北京:中國氣象科學(xué)研究院, 2015.

    [31] 劉紅芳,邸仕忠,姚啟倫. 熱帶和溫帶玉米對干旱的形態(tài)生理應(yīng)答[J]. 安徽大學(xué)學(xué)報(自然科學(xué)版), 2020, 44(3): 92-99.

    [32] 宋晉輝,甕巧云,呂愛枝,等. 拔節(jié)期干旱脅迫對青貯玉米生育與品質(zhì)的影響[J]. 中國農(nóng)業(yè)科技導(dǎo)報, 2020, 22(154): 167-173.

    [33] 宋 賀,蔣延玲,許振柱,等. 玉米光合生理參數(shù)對全生育期干旱與拔節(jié)后干旱過程的響應(yīng)[J]. 生態(tài)學(xué)報, 2019, 39(7): 2405-2415.

    [34] 印玉明,王永清,馬春晨,等. 利用日光誘導(dǎo)葉綠素?zé)晒獗O(jiān)測水稻葉片葉綠素含量[J]. 農(nóng)業(yè)工程學(xué)報, 2021, 37(12): 169-180.

    [35] SUN Y, FU R, DICKINSON R, et al. Drought onset mechanisms revealed by satellite solar-induced chlorophyll fluorescence: Insights from two contrasting extreme events[J]. Journal of Geophysical Research(Biogeosciences), 2016, 120: 2427-2440.

    [36] CHEN X, MO X, ZHANG Y, et al. Drought detection and assessment with solar-induced chlorophyll fluorescence in summer maize growth period over North China Plain[J]. Ecological Indicators, 2019, 104: 347-356.

    (責(zé)任編輯:蔣永忠)

    收稿日期:2021-08-16

    基金項目:國家自然科學(xué)基金項目(41905100);安徽省自然科學(xué)基金項目(2108085QD157、1908085QD171);國家重點研發(fā)計劃項目(2018YFD0300905);安徽農(nóng)業(yè)大學(xué)青年基金重點項目(2018zd07);安徽農(nóng)業(yè)大學(xué)引進(jìn)與穩(wěn)定人才資助項目(yj2018-57)

    作者簡介:倪 鋒(1997-),男,安徽安慶人,碩士研究生,主要從事農(nóng)業(yè)氣象研究。(E-mail)2934608403@qq.com

    通訊作者:李 萌,(E-mail)mengli@ahau.edu.cn;褚榮浩,(E-mail)ronghao_chu@163.com

    猜你喜歡
    夏玉米
    控釋摻混肥對夏玉米氮肥利用效率和氣態(tài)氮損失的影響
    巴西2020—2021年度夏玉米產(chǎn)量預(yù)計減少17.7%
    小麥?zhǔn)崭钪?如何種植夏玉米才能高產(chǎn)
    夏玉米高產(chǎn)的關(guān)鍵栽培技術(shù)措施
    長期定位施肥對夏玉米鉀素吸收及土壤鉀素動態(tài)變化的影響
    相對SPAD值用于不同品種夏玉米氮肥管理的研究
    冬小麥、夏玉米一體化壟作的養(yǎng)分利用研究
    施氮量對超高產(chǎn)夏玉米產(chǎn)量及氮素吸收利用的影響
    不同密閉材料硫包衣尿素氮素釋放特性及對夏玉米生長的影響
    天津市2011年良種補貼推介品種目錄
    亚洲,一卡二卡三卡| 国产无遮挡羞羞视频在线观看| 大香蕉久久网| videossex国产| 欧美精品一区二区大全| 熟女人妻精品中文字幕| 亚洲精品乱码久久久久久按摩| 国产精品一区www在线观看| 久久人妻熟女aⅴ| 亚洲成人一二三区av| 视频中文字幕在线观看| 日韩精品有码人妻一区| 精品99又大又爽又粗少妇毛片| 这个男人来自地球电影免费观看 | 欧美亚洲 丝袜 人妻 在线| 成人二区视频| 一级毛片电影观看| 一级,二级,三级黄色视频| 人人妻人人看人人澡| 偷拍熟女少妇极品色| 大陆偷拍与自拍| a 毛片基地| 欧美三级亚洲精品| 婷婷色综合www| 国产熟女欧美一区二区| 又爽又黄a免费视频| 美女主播在线视频| 内射极品少妇av片p| 亚洲av免费高清在线观看| 99热国产这里只有精品6| 国产91av在线免费观看| 亚洲内射少妇av| 国产亚洲午夜精品一区二区久久| 亚洲国产欧美日韩在线播放 | av又黄又爽大尺度在线免费看| 成年人免费黄色播放视频 | 国产免费视频播放在线视频| av又黄又爽大尺度在线免费看| 亚洲中文av在线| 韩国高清视频一区二区三区| 激情五月婷婷亚洲| .国产精品久久| 久久精品国产亚洲网站| 亚洲欧美日韩东京热| 最近2019中文字幕mv第一页| 激情五月婷婷亚洲| av免费观看日本| 成人二区视频| 狂野欧美激情性xxxx在线观看| 成年人免费黄色播放视频 | 成人漫画全彩无遮挡| 成人毛片60女人毛片免费| 人妻 亚洲 视频| 丝袜脚勾引网站| 精品一区在线观看国产| 国产永久视频网站| 国产乱来视频区| 黑人猛操日本美女一级片| 亚洲av日韩在线播放| 欧美丝袜亚洲另类| 一本大道久久a久久精品| 久久99精品国语久久久| 日韩视频在线欧美| 久久青草综合色| 欧美人与善性xxx| 老熟女久久久| 大香蕉97超碰在线| 久久久久视频综合| 国产熟女欧美一区二区| 久久久久精品性色| 最近最新中文字幕免费大全7| 熟女av电影| 精品久久久久久久久av| 午夜福利,免费看| 女性被躁到高潮视频| 少妇人妻久久综合中文| 人人妻人人澡人人爽人人夜夜| 熟妇人妻不卡中文字幕| 午夜福利视频精品| 亚洲真实伦在线观看| 欧美日本中文国产一区发布| 国产黄色视频一区二区在线观看| 春色校园在线视频观看| 交换朋友夫妻互换小说| 一级爰片在线观看| 亚洲成人手机| 久久久国产一区二区| 91午夜精品亚洲一区二区三区| 在线亚洲精品国产二区图片欧美 | 黄色毛片三级朝国网站 | 少妇的逼水好多| 国产精品久久久久久精品电影小说| 免费播放大片免费观看视频在线观看| 在线看a的网站| 日韩精品免费视频一区二区三区 | 久久久久久久久久成人| 午夜福利影视在线免费观看| 婷婷色综合大香蕉| 久热这里只有精品99| 少妇人妻精品综合一区二区| 97在线视频观看| av在线app专区| 精品久久久久久电影网| 亚洲精品aⅴ在线观看| 热re99久久精品国产66热6| 91精品国产国语对白视频| 国产成人精品婷婷| 久久狼人影院| 在现免费观看毛片| 亚洲天堂av无毛| 十八禁网站网址无遮挡 | 免费黄频网站在线观看国产| 大片免费播放器 马上看| 成人国产av品久久久| 国产真实伦视频高清在线观看| 国产成人一区二区在线| 韩国高清视频一区二区三区| 国产精品久久久久成人av| 七月丁香在线播放| 91午夜精品亚洲一区二区三区| 亚洲图色成人| 春色校园在线视频观看| 日韩大片免费观看网站| 成人亚洲精品一区在线观看| 国产精品嫩草影院av在线观看| 亚洲伊人久久精品综合| 国产色婷婷99| av国产久精品久网站免费入址| 日韩中字成人| 一级,二级,三级黄色视频| 国产在线免费精品| 欧美成人午夜免费资源| 夜夜骑夜夜射夜夜干| 欧美国产精品一级二级三级 | 国产精品福利在线免费观看| 国产 精品1| 日韩大片免费观看网站| 久久99精品国语久久久| 午夜视频国产福利| 国产免费一区二区三区四区乱码| 久久热精品热| 成人免费观看视频高清| 丝瓜视频免费看黄片| 国产视频内射| 六月丁香七月| 日本爱情动作片www.在线观看| 视频区图区小说| freevideosex欧美| 伊人久久精品亚洲午夜| 久久久久网色| 日韩伦理黄色片| 亚洲欧美日韩东京热| 97在线人人人人妻| 曰老女人黄片| 亚洲色图综合在线观看| 午夜免费鲁丝| 中国三级夫妇交换| 国产极品天堂在线| 亚洲欧美精品自产自拍| 午夜免费男女啪啪视频观看| 亚洲国产欧美日韩在线播放 | 啦啦啦啦在线视频资源| 2018国产大陆天天弄谢| 22中文网久久字幕| a级一级毛片免费在线观看| 三级国产精品片| 插阴视频在线观看视频| 日日啪夜夜撸| 午夜影院在线不卡| a 毛片基地| 午夜福利视频精品| 一本色道久久久久久精品综合| 国产欧美亚洲国产| 亚洲人成网站在线观看播放| 丝袜喷水一区| 国产无遮挡羞羞视频在线观看| 狂野欧美激情性bbbbbb| 欧美xxxx性猛交bbbb| 国产高清国产精品国产三级| 亚洲av二区三区四区| 在线精品无人区一区二区三| 丝袜在线中文字幕| 久久韩国三级中文字幕| 国产有黄有色有爽视频| 国产欧美日韩一区二区三区在线 | 日韩欧美 国产精品| 秋霞伦理黄片| 国产伦精品一区二区三区四那| 极品少妇高潮喷水抽搐| 最后的刺客免费高清国语| 国产欧美另类精品又又久久亚洲欧美| 国产伦精品一区二区三区四那| 国产日韩欧美视频二区| 色视频在线一区二区三区| 免费久久久久久久精品成人欧美视频 | 久久久久久久久久久久大奶| 久久久久久久久久人人人人人人| 日韩,欧美,国产一区二区三区| 秋霞伦理黄片| 亚洲av不卡在线观看| 中国三级夫妇交换| 两个人的视频大全免费| 一个人看视频在线观看www免费| 亚洲久久久国产精品| 一级二级三级毛片免费看| 国产精品国产三级国产av玫瑰| 午夜av观看不卡| 精品国产一区二区三区久久久樱花| 久久婷婷青草| 热re99久久精品国产66热6| 丰满少妇做爰视频| av天堂久久9| 国国产精品蜜臀av免费| 国产亚洲午夜精品一区二区久久| 久久久久精品久久久久真实原创| 街头女战士在线观看网站| 亚洲精品aⅴ在线观看| 久久久久国产网址| 乱码一卡2卡4卡精品| 免费观看无遮挡的男女| 五月玫瑰六月丁香| 99热国产这里只有精品6| 这个男人来自地球电影免费观看 | 国产 精品1| 国产综合精华液| 午夜久久久在线观看| 国产伦精品一区二区三区视频9| 欧美日韩精品成人综合77777| 婷婷色麻豆天堂久久| 免费人成在线观看视频色| 成人美女网站在线观看视频| av女优亚洲男人天堂| 一本色道久久久久久精品综合| 最新中文字幕久久久久| 久久99一区二区三区| 亚洲成色77777| 久久亚洲国产成人精品v| 色94色欧美一区二区| 人妻夜夜爽99麻豆av| .国产精品久久| 国产精品一区二区三区四区免费观看| 久久精品国产亚洲av涩爱| 国产在线男女| 免费人成在线观看视频色| 国产高清三级在线| 国产91av在线免费观看| 久久99热6这里只有精品| 黄色视频在线播放观看不卡| 日本猛色少妇xxxxx猛交久久| 午夜激情久久久久久久| 欧美最新免费一区二区三区| 少妇人妻精品综合一区二区| 有码 亚洲区| av福利片在线| 日韩av在线免费看完整版不卡| 国产精品国产三级国产专区5o| 国产av国产精品国产| 日韩伦理黄色片| 精品久久久噜噜| 中文乱码字字幕精品一区二区三区| 人妻人人澡人人爽人人| 国产高清有码在线观看视频| 校园人妻丝袜中文字幕| 日日摸夜夜添夜夜爱| 又大又黄又爽视频免费| 亚洲熟女精品中文字幕| 精品国产一区二区三区久久久樱花| 国产探花极品一区二区| 看免费成人av毛片| 亚洲精品成人av观看孕妇| 日日爽夜夜爽网站| 美女xxoo啪啪120秒动态图| 欧美国产精品一级二级三级 | 一本—道久久a久久精品蜜桃钙片| 国产成人免费观看mmmm| 国产一区二区在线观看av| 亚洲国产精品专区欧美| 日韩三级伦理在线观看| 免费大片18禁| av天堂久久9| 又粗又硬又长又爽又黄的视频| 这个男人来自地球电影免费观看 | 在线观看免费视频网站a站| 久久婷婷青草| 亚洲成色77777| 精品国产一区二区三区久久久樱花| 欧美日韩一区二区视频在线观看视频在线| 日本-黄色视频高清免费观看| 两个人的视频大全免费| 国产免费一级a男人的天堂| 日本vs欧美在线观看视频 | 欧美日韩视频精品一区| 国产精品一区二区性色av| 亚洲一区二区三区欧美精品| 欧美xxⅹ黑人| 欧美日韩一区二区视频在线观看视频在线| 欧美日本中文国产一区发布| 亚洲欧美日韩东京热| 成人亚洲欧美一区二区av| 一级毛片aaaaaa免费看小| 国产精品麻豆人妻色哟哟久久| 亚洲精品日本国产第一区| 18禁在线播放成人免费| 免费久久久久久久精品成人欧美视频 | 老司机亚洲免费影院| 国产在视频线精品| 国语对白做爰xxxⅹ性视频网站| 黄色视频在线播放观看不卡| 婷婷色av中文字幕| 欧美日韩亚洲高清精品| 丝袜在线中文字幕| 久久99热6这里只有精品| 一级爰片在线观看| 国产在线男女| 免费少妇av软件| 国产成人精品一,二区| 中文在线观看免费www的网站| 久久精品国产鲁丝片午夜精品| 免费少妇av软件| 熟女电影av网| 国产熟女午夜一区二区三区 | 国产精品福利在线免费观看| 最黄视频免费看| 日本vs欧美在线观看视频 | 色视频在线一区二区三区| av国产久精品久网站免费入址| 韩国av在线不卡| 亚洲欧美一区二区三区国产| 丝瓜视频免费看黄片| 插阴视频在线观看视频| 国产一区亚洲一区在线观看| 91精品伊人久久大香线蕉| 一区二区三区四区激情视频| 嫩草影院新地址| 美女xxoo啪啪120秒动态图| 成人无遮挡网站| 日韩免费高清中文字幕av| 黄色视频在线播放观看不卡| 久久久久久久大尺度免费视频| 亚洲人成网站在线播| 黄色配什么色好看| 我的老师免费观看完整版| 只有这里有精品99| 精华霜和精华液先用哪个| 黄色配什么色好看| 国产在线视频一区二区| 国产在线免费精品| av福利片在线观看| 一级毛片黄色毛片免费观看视频| 肉色欧美久久久久久久蜜桃| 91午夜精品亚洲一区二区三区| 国产 一区精品| 十分钟在线观看高清视频www | 搡女人真爽免费视频火全软件| 桃花免费在线播放| 免费av中文字幕在线| 青青草视频在线视频观看| 日本欧美国产在线视频| 日韩一区二区三区影片| 午夜精品国产一区二区电影| 国产精品麻豆人妻色哟哟久久| 亚洲性久久影院| 国产熟女欧美一区二区| 蜜臀久久99精品久久宅男| 久久久久久伊人网av| 亚洲一级一片aⅴ在线观看| 国产日韩一区二区三区精品不卡 | 熟女人妻精品中文字幕| 国产亚洲一区二区精品| 精品99又大又爽又粗少妇毛片| 成年女人在线观看亚洲视频| 欧美日本中文国产一区发布| 在线天堂最新版资源| 黄片无遮挡物在线观看| 久久久久久伊人网av| 国产亚洲欧美精品永久| 天堂8中文在线网| 我要看黄色一级片免费的| 在线亚洲精品国产二区图片欧美 | 99久久综合免费| 狂野欧美白嫩少妇大欣赏| 久久久久久久久久久丰满| 国产一区二区三区综合在线观看 | 国产探花极品一区二区| 另类亚洲欧美激情| 欧美日韩综合久久久久久| 免费人妻精品一区二区三区视频| 日韩中文字幕视频在线看片| 国国产精品蜜臀av免费| 成年人午夜在线观看视频| 亚洲欧美一区二区三区国产| 内地一区二区视频在线| 麻豆成人av视频| 大话2 男鬼变身卡| 自拍欧美九色日韩亚洲蝌蚪91 | 久久久久久久精品精品| 女性被躁到高潮视频| 乱人伦中国视频| 97超碰精品成人国产| 91aial.com中文字幕在线观看| 一区在线观看完整版| 免费人妻精品一区二区三区视频| 久久这里有精品视频免费| 日日摸夜夜添夜夜添av毛片| 十八禁网站网址无遮挡 | 亚洲av成人精品一二三区| 久久久久久久精品精品| 亚洲欧美中文字幕日韩二区| 亚洲精品乱码久久久v下载方式| av视频免费观看在线观看| 老司机影院毛片| 亚洲精品自拍成人| 街头女战士在线观看网站| 欧美日韩精品成人综合77777| 国产精品国产三级国产av玫瑰| 国产日韩一区二区三区精品不卡 | 国产男人的电影天堂91| 高清av免费在线| 偷拍熟女少妇极品色| 国产 精品1| 99热这里只有精品一区| 美女中出高潮动态图| 亚洲国产欧美日韩在线播放 | 精品一品国产午夜福利视频| 久久久久久人妻| 人妻系列 视频| 三级国产精品片| 性色avwww在线观看| 国产精品一二三区在线看| 又黄又爽又刺激的免费视频.| 国产白丝娇喘喷水9色精品| 国产黄色免费在线视频| 91午夜精品亚洲一区二区三区| 在线观看美女被高潮喷水网站| 亚洲欧美一区二区三区国产| 大又大粗又爽又黄少妇毛片口| 91精品一卡2卡3卡4卡| 国内少妇人妻偷人精品xxx网站| 国产真实伦视频高清在线观看| 免费av中文字幕在线| 欧美精品人与动牲交sv欧美| 国产真实伦视频高清在线观看| 午夜老司机福利剧场| 国精品久久久久久国模美| 晚上一个人看的免费电影| 国产av精品麻豆| 日韩强制内射视频| 亚洲av在线观看美女高潮| 日本免费在线观看一区| 哪个播放器可以免费观看大片| 看非洲黑人一级黄片| 亚洲内射少妇av| 一区二区三区免费毛片| 51国产日韩欧美| 中文字幕人妻熟人妻熟丝袜美| 黄色一级大片看看| 国产午夜精品久久久久久一区二区三区| 国产国拍精品亚洲av在线观看| 夜夜骑夜夜射夜夜干| 韩国高清视频一区二区三区| 国产男女超爽视频在线观看| 美女主播在线视频| 日本av手机在线免费观看| 亚洲精品成人av观看孕妇| 80岁老熟妇乱子伦牲交| 亚洲电影在线观看av| 成人国产av品久久久| 久久热精品热| 精品国产一区二区久久| 夫妻午夜视频| 国内精品宾馆在线| 免费黄频网站在线观看国产| 老司机影院毛片| 少妇人妻 视频| 国产免费视频播放在线视频| 如何舔出高潮| 岛国毛片在线播放| 日韩伦理黄色片| 精品一区在线观看国产| 亚洲欧美一区二区三区国产| 精品一区二区三卡| 色94色欧美一区二区| 女人精品久久久久毛片| 美女内射精品一级片tv| 国产无遮挡羞羞视频在线观看| 亚洲人成网站在线观看播放| 男女免费视频国产| 啦啦啦在线观看免费高清www| 色视频在线一区二区三区| 亚洲欧美清纯卡通| 丝瓜视频免费看黄片| 欧美高清成人免费视频www| 插逼视频在线观看| 观看av在线不卡| 蜜臀久久99精品久久宅男| 成人18禁高潮啪啪吃奶动态图 | 久久这里有精品视频免费| 国产欧美亚洲国产| av又黄又爽大尺度在线免费看| 高清欧美精品videossex| 国产一级毛片在线| 国产探花极品一区二区| 国产一区二区三区av在线| 亚洲熟女精品中文字幕| 国产精品国产三级国产av玫瑰| 色婷婷久久久亚洲欧美| 又黄又爽又刺激的免费视频.| 欧美区成人在线视频| 99热国产这里只有精品6| 国产高清有码在线观看视频| 在现免费观看毛片| 国产欧美日韩一区二区三区在线 | 精华霜和精华液先用哪个| 久久人人爽av亚洲精品天堂| 熟女人妻精品中文字幕| 国产欧美日韩精品一区二区| 日韩人妻高清精品专区| 中国国产av一级| 欧美 日韩 精品 国产| 日韩av在线免费看完整版不卡| 亚洲av电影在线观看一区二区三区| 丰满饥渴人妻一区二区三| 久久亚洲国产成人精品v| 少妇丰满av| 亚洲欧美清纯卡通| 人人妻人人看人人澡| 在线观看免费日韩欧美大片 | 日本91视频免费播放| 黄色日韩在线| 国产男人的电影天堂91| 一区二区三区乱码不卡18| 中文字幕人妻熟人妻熟丝袜美| 男女免费视频国产| 久久久久人妻精品一区果冻| 高清欧美精品videossex| 久久精品国产亚洲av天美| 黄色日韩在线| 日韩大片免费观看网站| 国产 精品1| 国产在线视频一区二区| 色婷婷av一区二区三区视频| 在现免费观看毛片| av线在线观看网站| 久久久久久久国产电影| 亚洲电影在线观看av| 欧美另类一区| 亚洲精品国产av蜜桃| 全区人妻精品视频| av福利片在线| 国产亚洲最大av| 免费黄色在线免费观看| 久热这里只有精品99| 久久女婷五月综合色啪小说| 精品久久久久久久久av| 性色av一级| 久久久久国产网址| 国产一区亚洲一区在线观看| 免费大片18禁| 两个人免费观看高清视频 | 观看av在线不卡| 在线观看免费高清a一片| 91成人精品电影| 美女大奶头黄色视频| 夫妻午夜视频| 99热全是精品| 久久久亚洲精品成人影院| 亚洲av电影在线观看一区二区三区| 一本一本综合久久| 在线观看免费视频网站a站| 欧美日韩视频精品一区| 99久国产av精品国产电影| 丝袜在线中文字幕| 在线天堂最新版资源| 激情五月婷婷亚洲| 特大巨黑吊av在线直播| 少妇被粗大猛烈的视频| 亚洲精品乱码久久久v下载方式| 观看av在线不卡| 久久国内精品自在自线图片| 欧美变态另类bdsm刘玥| 看非洲黑人一级黄片| 99久国产av精品国产电影| 十八禁高潮呻吟视频 | 在线观看免费视频网站a站| 亚洲av国产av综合av卡| 我的女老师完整版在线观看| 午夜精品国产一区二区电影| 人妻夜夜爽99麻豆av| 久久青草综合色| 亚洲三级黄色毛片| 欧美精品一区二区大全| 久久青草综合色| 精品视频人人做人人爽| 亚洲熟女精品中文字幕| 亚洲精品中文字幕在线视频 | 丝袜喷水一区| 久久免费观看电影| 一区二区三区免费毛片| 国产成人一区二区在线| 国产国拍精品亚洲av在线观看| 亚洲人成网站在线播| 亚洲精品乱码久久久久久按摩| 欧美另类一区| 三级国产精品片| 国产精品福利在线免费观看| 午夜老司机福利剧场| 中文在线观看免费www的网站| 久久久午夜欧美精品| 夜夜骑夜夜射夜夜干| 国产欧美亚洲国产| 99久久中文字幕三级久久日本| 国产乱来视频区| 全区人妻精品视频| 97超视频在线观看视频| 嘟嘟电影网在线观看|