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

    秸稈覆蓋和滴灌制度對冬小麥光合特性和產(chǎn)量的影響

    2017-07-18 11:48:17張彥群王建東龔時宏孫繼文
    農(nóng)業(yè)工程學報 2017年12期
    關(guān)鍵詞:羧化旗葉冬小麥

    張彥群,王建東,龔時宏,許 迪,孫繼文,2

    ?

    秸稈覆蓋和滴灌制度對冬小麥光合特性和產(chǎn)量的影響

    張彥群1,王建東1※,龔時宏1,許 迪1,孫繼文1,2

    (1. 中國水利水電科學研究院水利研究所,北京 100048;2. 山東理工大學資源與環(huán)境工程學院,淄博 255049)

    確定不同滴灌和秸稈覆蓋模式下的光合參數(shù)是理解其節(jié)水增產(chǎn)內(nèi)因的基礎(chǔ)。該研究通過4 a連續(xù)試驗測定了覆蓋和滴灌結(jié)合的6種處理(T1~T6分別表示低水覆蓋、低水不覆蓋、中水覆蓋、中水不覆蓋、高水覆蓋和高水不覆蓋處理)冬小麥旗葉的關(guān)鍵光合參數(shù),并分析了光合參數(shù)處理間差異來源及其與作物產(chǎn)量的關(guān)系。結(jié)果表明:提高滴灌水量和覆蓋處理能顯著提高光合能力和最大羧化速率(< 0.05),而二者交互作用對光合參數(shù)影響不顯著(> 0.1)。多重比較結(jié)果顯示T5的光合參數(shù)最大(< 0.05);除2016年第一次測定外,T6和T3的相應(yīng)值略低但與T5差異不顯著(> 0.1);而T2的相應(yīng)值則均顯著低于T5(< 0.05),即高水處理下,覆蓋與否對光合參數(shù)的影響不大,中水結(jié)合覆蓋處理可使光合參數(shù)不顯著降低,而低水處理下不覆蓋則會使光合參數(shù)顯著降低。處理間光合能力和最大羧化速率的差異能用葉片氮含量來解釋,而所有年份的產(chǎn)量又均與光合能力和最大羧化速率分別線性相關(guān)(2> 0.80)。因此,通過測定葉片氮含量,可以推算光合能力和最大羧化速率,研究結(jié)果可為產(chǎn)量估算提供參考。

    作物;灌水;光合作用;冬小麥;秸稈覆蓋;葉片氮含量

    0 引 言

    華北平原是中國糧食主產(chǎn)區(qū),該區(qū)小麥和玉米產(chǎn)量占全國糧食總產(chǎn)量的1/4[1]。近年來,由于該地區(qū)水資源短缺和地下水超采嚴重,糧食作物節(jié)水灌溉工程和農(nóng)藝措施逐步推進[2]。地表滴灌將水分施于作物根區(qū)附近,具有降低土壤蒸發(fā),提高作物水分利用效率的優(yōu)勢[3]。秸稈覆蓋是改善農(nóng)田水土環(huán)境、降低土壤蒸發(fā)、提高作物產(chǎn)量采用的常見農(nóng)藝措施[4-5]。秸稈覆蓋和地表滴灌結(jié)合模式下,節(jié)水增產(chǎn)優(yōu)勢發(fā)揮的生理基礎(chǔ)仍需要明確。目前相關(guān)研究更多的是從相對表觀的生物學指標,如作物產(chǎn)量、株高等方面來分析其優(yōu)勢[6-7]。然而,作物產(chǎn)量的形成涉及到諸多生理過程,光合作用是產(chǎn)量形成的基礎(chǔ),從光合生理方面來解釋產(chǎn)量差異來源,有利于試驗結(jié)果的穩(wěn)定性和可重復性,并為機理模型的構(gòu)建奠定基礎(chǔ)[8]。此外,該地區(qū)現(xiàn)狀研究中缺少對覆蓋免耕滴灌條件下內(nèi)在光合參數(shù)的系統(tǒng)測定,制約了其作物模型精度的提高[9-10]。

    光合生理調(diào)控研究在冬小麥光合的研究中還不夠系統(tǒng)。目前,對冬小麥旗葉光合作用的研究大多局限于描述作物某一生育期或不同生育期中某一時刻的光合速率或光合速率日變化[11-12],對光合速率變化的內(nèi)因解釋較少,且每次測定所處的光、溫、濕環(huán)境各不相同,所得到的光合參數(shù)受到天氣條件制約,不同種植區(qū)域、作物品種和肥料水平之間可比性較差[13]。因此,有必要開展控制條件測定,量化旗葉的光合特征參數(shù)而非瞬時速率,并確定不同處理光合參數(shù)的差異來源,從生理水平解釋滴灌制度和覆蓋處理引起產(chǎn)量及水分利用效率差異的內(nèi)因[14-16]。小麥產(chǎn)量形成過程中,葉片從下到上逐漸衰老,冠層頂葉片,即旗葉的光合產(chǎn)物成為灌漿期同化物的主要來源[17]。覆蓋及灌溉模式可能會影響旗葉光合生理,進而導致產(chǎn)量的差別,光合生理參數(shù)又與旗葉生化特性有關(guān),葉片氮含量是重要的旗葉生化特性之一[18]。氮素是葉綠素的重要組成成分,葉綠素是光合作用進行的載體,光合速率與葉片氮含量常存在顯著相關(guān)性[19-20]。研究小麥旗葉光合特性及其影響因子對覆蓋和滴灌處理的響應(yīng),對調(diào)控旗葉光合氣體交換,穩(wěn)定提高產(chǎn)量及水分利用效率具有積極的促進作用。

    本研究選取華北平原冬小麥為研究對象,設(shè)置了不同滴灌和覆蓋處理,2013- 2016年連續(xù)4個生長季在冬小麥抽穗期后分別進行了兩次旗葉光合光響應(yīng)和CO2響應(yīng)曲線測定,獲得了表觀光量子傳遞效率、光合能力、最大羧化速率、最大電子傳遞速率等關(guān)鍵參數(shù),比較了不同處理條件下光合參數(shù)的差異,并同時測定旗葉氮含量,以期解釋水分和覆蓋處理引起光合參數(shù)差異性,進而揭示產(chǎn)量響應(yīng)的光合生理基礎(chǔ),并為該地區(qū)小麥作物模型參數(shù)確定提供參考。

    1 材料與方法

    1.1 試驗地概況

    試驗在位于北京市大興區(qū)(39°39′N,116°15′E)的中國水利水電科學研究院農(nóng)業(yè)節(jié)水灌溉試驗站進行。試驗站氣候?qū)俚湫偷陌敫珊荡箨懶约撅L氣候,多年平均降雨量540 mm,小麥生長季(10月-來年3月)的降雨量常不足100 mm。試驗田0~100 cm土層深度的土壤質(zhì)地為壤土,0~100 cm土層的平均田間持水率、土壤容重和有機質(zhì)質(zhì)量分數(shù)分別為30.58%、1.58 g/cm3和3.25 g/kg。田間試驗于2012-2013、2013-2014、2014-2015和 2015-2016年連續(xù)4個生長季開展。開展試驗前(2012年10月)土壤的硝態(tài)氮質(zhì)量分數(shù)為5.43 mg/kg。

    1.2 試驗設(shè)計

    本試驗在覆蓋和不覆蓋2種處理下分別設(shè)置3種滴灌制度,共6個處理,每個處理設(shè)4個重復,共24個小區(qū),小區(qū)面積8×10 m2,隨機排列。冬小麥于每年10月中旬播種,行距25 cm,播前根據(jù)土壤墑情,灌水50~75 mm,并施入復合肥(有效成分:N-P2O5-K2O,比例:15%-15%-15%)556 kg/hm2作為底肥,折合純N、P、K量分別為83.4、36.4和69.2 kg/hm2。播種當年11月中旬冬灌,冬灌采用地面灌,灌水量為50 mm。次年4月初返青,開始灌水處理,灌水方式為地表滴灌,一帶4行,滴頭間距為30 cm,滴頭流量為1.35 L/h。

    灌溉制度參考王建東等[21]在該地的試驗,為保證作物正常生長,所有處理返青期第一次灌水均灌至田持(30.58%),此后至成熟期前,根據(jù)計劃濕潤層的土壤含水率上下限控制來確定灌水量和灌水時間,低、中、高3種灌溉制度分別以田持的55%~75%、65%~85%和78~98%作為灌水上下限控制因素,成熟期不灌水。覆蓋處理采用上季留存粉碎后的玉米秸稈,覆蓋量為6 000 kg/hm2。除2012-2013年試驗季于2013年5月4日開始覆蓋外,其余試驗季均為播種后即用秸稈覆蓋。每個處理返青期第一次灌水時,采用文丘里施肥器施肥,施氮量為 110 kg/hm2,施用肥料為尿素。各處理編號、生育期灌水次數(shù)、灌水量及降雨量見表1,處理1-6(T1~T6)分別表示低水覆蓋、低水不覆蓋、中水覆蓋、中水不覆蓋、高水覆蓋和高水不覆蓋處理。

    1.3 旗葉光合參數(shù)及氮含量測定

    冬小麥旗葉展開后,選擇典型晴天用Li-cor 6400便攜式光合儀進行葉片光合—光響應(yīng)曲線和光合—CO2響應(yīng)曲線的測定,進而確定葉片相關(guān)光合參數(shù)。2013-2016年每年進行2次測定。為方便各年份間不同測定日期間的比較,所有測定日期均標準化為距返青期首次灌水的天數(shù)(Dafi),標準化后,2013年的測定時間為Dafi19-20和Dafi40-41、2014年的測定時間為Dafi29-30和Dafi47- 48、2015年的測定時間為Dafi35-36和Dafi54-55、2016年的測定時間為Dafi32-35和Dafi57-58。光響應(yīng)曲線測定時葉室溫度設(shè)定為26~28 ℃,葉室CO2摩爾分數(shù)設(shè)定為400mol/mol,光強設(shè)定為2200、2000、1600、1300、1000、750、500、200、100、50、20和0mol/(m2×s),從強到弱進行測定。CO2響應(yīng)曲線測定時葉室溫度設(shè)定為26~28 ℃,葉室光強設(shè)定為2000mol/(m2×s),CO2濃度設(shè)定為400、300、200、120、80、50、0、400、400、600、800、1200、1600、2000mol/mol,按次序測定。

    表1 各試驗季處理編號、返青至收獲期灌水次數(shù)、灌水量及降雨量統(tǒng)計表 Table 1 Lists of treatment coding, irrigation times, irrigation amounts and precipitation during the reviving and harvest stages for all experiment seasons

    光合—光響應(yīng)曲線測定結(jié)果用非直角雙曲線方程(式1)[22]擬合獲得相關(guān)參數(shù),進行處理間比較分析。其中,光強0~150mol/(m2×s)的數(shù)據(jù)采用直線擬合,直線的斜率為表觀光量子效率(,無量綱),直線與縱坐標的交點為呼吸速率(d,mol/(m2×s)),回歸直線與橫坐標的交點為光補償點(PARc,mol/(m2×s))。確定上述參數(shù)后,采用固定曲率(=0.9)的非直角雙曲線進行擬合,得到最大光合速率,即光合能力(max,mol/(m2×s))。非直角雙曲線模型的表達式為:

    式中為植物的總光合速率(mol/(m2×s)),n為植物的凈光合速率(mol/(m2×s)),PAR為光合有效輻射(mol/ (m2×s)),為曲線的曲率(無量綱),曲率越大,曲線的彎曲程度越大。

    光合—CO2曲線測定結(jié)果采用Sharkey等完善了的Farquhar和von Caemmerer提出的光合模型擬合[23]。該模型指出n的不同階段分別受2個過程限制:Rubisco酶的羧化效率(c,mol/(m2×s))和電子傳遞速率 (J,mol/ (m2×s))。各個過程的表達式見式(2)~式(6):

    (3)

    (4)

    (6)

    式中c和o分別為Rubisco酶的羧化和加氧速率,單位均為mol/(m2×s);*為無暗呼吸時的CO2補償點(mol/mol);d為呼吸速率(mol/(m2×s));cmax和max分別為最大羧化速率和最大電子傳遞速率,單位均為mol/(m2×s);i為細胞間CO2摩爾分數(shù)(mol/mol);c(Pa)和o(kPa)分別為CO2和O2的米氏常數(shù);為氧氣分壓(kPa);是電子傳遞速率(mol/(m2×s)),與、PAR及max有關(guān);通過擬合,可以得到光合參數(shù)cmax和max。

    光合測定完成后將旗葉取下,掃描葉面積,105 ℃殺青,65 ℃下烘干至恒量,研磨,過100目篩,采用硫酸-水楊酸消煮-凱式定氮儀法確定葉片氮質(zhì)量分數(shù)(mass,%)。收獲時各小區(qū)單收單打計產(chǎn),小區(qū)小麥產(chǎn)量以實收產(chǎn)來計算,籽粒曬干至質(zhì)量無變化時測其產(chǎn)量。

    1.4 數(shù)據(jù)分析

    光合—光響應(yīng)曲線和光合—CO2響應(yīng)曲線的擬合均采用SPSS18.0軟件的自定義非線性擬合進行,參數(shù)確定后首先采用SPSS18.0軟件的雙因素方差分析確定覆蓋和滴灌處理及其交互作用對光合參數(shù)的影響,在確定不存在二者交互作用的前提下,采用單因素方差分析及dacan多重比較方法進行所有處理的差異顯著性分析及處理間的均值的兩兩比較。在分析光合參數(shù)與葉片氮含量關(guān)系及產(chǎn)量與光合參數(shù)的關(guān)系時,首先采用SPSS18.0軟件的協(xié)方差分析確定不同年份之間回歸直線斜率和截距是否存在顯著差異,不存在顯著差異的情況下,再采用Sigma-Plot12.0軟件將所有年份統(tǒng)一回歸分析,確定回歸方程和統(tǒng)計參數(shù)。文中所有作圖均采用SigmaPlot12.0軟件進行。

    2 結(jié)果與分析

    2.1 光合—光響應(yīng)曲線參數(shù)差異

    不同滴灌水量和覆蓋處理不同程度影響了表觀光量子效率和光合能力max(表2)。表2中打亂年份順序,按照測定日距返青后首次灌水的天數(shù)(Dafi)排列,依次反映旗葉展開至成熟期間不同生育階段光合參數(shù)的處理間差異情況(表3采用同樣的排列順序)。其中,2013年每個處理的光曲線只測定了一條,無法進行均值比較,2014-2016年的每次測定均至少測定兩條光曲線,可以進行均值比較及方差分析。

    滴灌水量對值的顯著影響僅表現(xiàn)在抽穗灌漿后期的兩次測定(Dafi47-48和54-55,表2)中,而并未發(fā)現(xiàn)覆蓋處理對值的顯著影響,也未發(fā)現(xiàn)滴灌水量和覆蓋處理的交互作用。6個處理間的比較結(jié)果也顯示,僅抽穗灌漿后期最后一次測定的值處理間差異顯著(= 0.038),倒數(shù)第二次測定的值處理間差異接近顯著(= 0.078)。

    表2 滴灌水量和覆蓋處理及其交互作用對表觀光量子效率α和光合能力Amax影響的顯著性水平 Table 2 Results of two-way analysis of variance of irrigation amount, mulching and interaction effects of above two treatments on apparent quantum efficiency (α) and photosynthetic capacity (Amax)

    注:**表示值小于0.05,*表示值介于0.05~0.1之間,下同。

    Notes: ** stands for differences among treatments are significant at= 0.05, * stands for differences among treatments are significant at 0.05<< 0.1. The same below.

    圖1顯示了值的處理間多重比較結(jié)果。2014年,Dafi47-48測定中,T1的值與T5差異接近顯著,T5的值較高,而T2、T6的值則處于上述兩處理之間,且與之均無顯著差異。2015年,Dafi54-55測定中,滴灌水量之間值差異顯著,低水的T1、T2的值顯著低于高水處理T5、T6,且顯著低于中水覆蓋處理T3,T4的值介于處理1、2和T3、T5、T6之間,且與之均無顯著差異。T1~T6的值4個生長季平均分別為:0.043、0.041、0.048、0.048、0.050和0.047。陸佩玲和于強研究指出,大田條件下小麥的值一般在0.05~0.07之間[24],本研究只有2015-2016年度高水和覆蓋處理的值在0.06左右,而其余年份,尤其是低水不覆蓋處理的值較低,表明覆蓋和提高滴灌水量有利于提高。

    滴灌水量除對2016年Dafi32-35測定中max的影響接近顯著(=0.077)外,其余測定日中,滴灌水量對max均影響顯著(<0.05,在0.001-0.019之間,表2)。覆蓋處理除對2014年Dafi29-30的max值影響不顯著和對2015年Dafi35-36的max值影響接近顯著(=0.094)外,其余測定日中,覆蓋對max均影響顯著(<0.05,在0.001~0.032之間)。然而,研究并未發(fā)現(xiàn)不同滴灌水量和覆蓋處理的交互作用。6個處理間的比較結(jié)果顯示,所有測定日的max值處理間差異均顯著(<0.05,在0.005~ 0.049之間)。

    圖2顯示了max值的處理間多重比較結(jié)果。2014年,Dafi29-30測定中,T2的max值顯著低于T5和T6,T1介于T2和T5、T6之間,與之均無顯著差別;Dafi47-48測定中,T1、T2的max值差異不顯著,二者均顯著低于T5、T6,T5的max值又顯著大于T6。2015年,Dafi35-36測定中,max值差異僅出現(xiàn)在T5、T6和T2、T4之間;Dafi54-55測定中,T2的max值則顯著低于其余處理,T1和T4則顯著低于T3、T5,而T6的max值介于上述兩組處理直接,與之均無顯著差別。2016年,Dafi32-35測定中,T2、T4的max值顯著低于其余處理,而T1、T3、T5的max值差異不顯著,T6的max值介于T1、T2和T5之間,僅顯著低于T5值,與T1、T3差異不顯著。T1~T6的max值4個生長季平均分別為:29.59、26.45、31.44、28.68、32.42和30.90mol/(m2×s)。從max值差異來看,越到生育后期,高水和覆蓋處理(處理5)的優(yōu)勢顯著,而中水覆蓋處理(處理3)的max值與高水覆蓋處理差異均不顯著。中、高水覆蓋處理max在生育期末仍較高,有利于籽粒充實,提高作物產(chǎn)量,這與前人研究結(jié)果一致[25]。

    注:T1~T6分別表示低水覆蓋、低水不覆蓋、中水覆蓋、中水不覆蓋、高水覆蓋和高水不覆蓋處理。圖中不同大寫和小寫字母分別表示處理間在0.1和0.05水平上差異顯著,下同。

    2.2 光合—CO2響應(yīng)曲線參數(shù)差異

    不同滴灌水量和覆蓋處理均不同程度的影響了最大羧化速率cmax和最大電子傳遞效率max(表3)。其中,2013-2015年第二次測定每個處理的CO2響應(yīng)曲線只測定了一條,無法進行均值比較,其余測定日均至少測定兩條曲線,可以進行均值比較及方差分析。

    滴灌水量對cmax值的顯著影響主要表現(xiàn)在抽穗灌漿前、中期的三次測定(Dafi35天以前)中,覆蓋處理對cmax值的顯著影響則主要表現(xiàn)在抽穗灌漿中、后期(Dafi32天以后),研究并未發(fā)現(xiàn)滴灌水量和覆蓋處理的交互作用(表3)。6個處理間的比較結(jié)果顯示,5次測定中有4次處理間差異達到或接近顯著(<0.1,在0.020~0.085之間),僅2015年Dafi35-36測定中處理間cmax值差異不顯著(=0.202)。

    圖2 2013-2016年不同處理冬小麥旗葉最大光合速率Amax比較

    表3 滴灌水量和覆蓋處理及其交互作用對最大羧化速率Vcmax和最大電子傳遞效率Jmax影響的顯著性水平

    圖3顯示了cmax值的處理間多重比較結(jié)果。2013年,Dafi19-20的測定中,低水處理T1、T2的cmax值顯著低于高水處理T5、T6。2014年,Dafi29-30的測定中,T2的cmax值在0.1顯著水平上低于T5、T6,T1的cmax值介于T2和T5、T6之間,且與兩組處理值差異不顯著。2016年,Dafi33-35的測定中,T2的cmax值顯著低于T3、T4、T5、T6,而T1的cmax值介于T2和T3、T4、T5、T6之間,且與兩組處理值差異不顯著,T4的cmax值顯著低于T5;Dafi57-58的測定中,盡管處理間cmax均值差別較大,但處理內(nèi)個體差異(標準差)也較大,處理間差異僅在0.1水平上顯著,其中,只有T2的cmax值顯著低于T5,其余處理的cmax值介于T2和T5之間,且與二者差異不顯著。T1~T6的cmax值4個生長季平均分別為:121.42、107.68、133.42、113.55、141.83和129.95mol/(m2×s)。

    圖3 2013-2016年各處理冬小麥旗葉最大羧化速率Vcmax比較

    滴灌水量對max值的顯著影響在抽穗灌漿前、中期兩次(Dafi19-20和Dafi32-35)觀測到,而覆蓋處理對max值的顯著影響則僅在抽穗灌漿后期末次測定中(Dafi57-58)觀測到,同樣的,研究并未發(fā)現(xiàn)滴灌水量和覆蓋處理對max的交互作用(表3)。6個處理間的比較結(jié)果顯示,5次測定中有3次處理間差異僅為接近顯著(0.05<<0.1)。

    圖4顯示了max值的處理間多重比較結(jié)果。2013年,Dafi19-20的測定中,低水處理T1、T2的max值在0.1顯著水平低于高水處理T5、T6。2016年,Dafi32-35的測定中,T1、T2的max值在0.1顯著水平低于T5、T6,T3、T4的max值介于T1、T2和T5、T6之間,且與兩組處理值差異不顯著;Dafi57-58的測定中,與cmax類似,處理內(nèi)變異較大,可能掩蓋了處理間差異,僅有T2、T4的max值在0.1顯著水平低于T5,其余處理與上述處理均差異不顯著。T1~T6的max值4個生長季平均分別為:246.21、219.24、265.20、248.20、287.04和271.63mol/(m2×s)。從降低灌水量的角度來講,T3的灌溉制度可以作為秸稈覆蓋冬小麥田的滴灌制度,即65%~85%的田間持水量作為土壤水分的上下限控制,該灌溉制度配合田間秸稈覆蓋措施,可以保持冬小麥旗葉較高的光合能力。

    本研究確定的不同處理的cmax值遠高于Wullschleger綜述中對小麥cmax的統(tǒng)計值(83mol/(m2×s))[26],但與近期Müller等對小麥cmax的實測值接近[27],這與近年來小麥品種改良、光合能力提高有關(guān)。cmax的處理間差異一般與旗葉氮含量相關(guān),本研究將在下文詳述。在cmax已知的情況下,光合模型參數(shù)確定時通常采用max與cmax的比值來估算max值,前人研究普遍認為max與cmax的比值較為固定,取值為(2±0.6)[28],本研究各處理的max/cmax比值在1.35~2.68之間變化,均值為2.0,該范圍包含了Evans和Farquhar對小麥光合研究中報道的該比值1.63~2.06[29],且均值與Leuning推薦的C3植物普適值(2.0±0.6)一致[28]?;诒狙芯拷Y(jié)果,筆者推薦無基礎(chǔ)資料研究中冬小麥光合模擬可采用比值2來推算max。

    圖4 2013-2016年不同處理冬小麥最大電子傳遞效率Jmax比較

    2.3 光合參數(shù)與葉片N含量的關(guān)系

    氮素是作物葉綠素及光合相關(guān)酶的重要組成部分,光合速率一般與葉片N含量顯著相關(guān),本研究分析了光合—光強響應(yīng)曲線關(guān)鍵光合參數(shù)max和光合—CO2響應(yīng)曲線關(guān)鍵光合參數(shù)cmax與葉片N含量(mass)的關(guān)系。所有年份、所有測定日期的max和cmax均與mass線性相關(guān)顯著,且不同年份之間相關(guān)關(guān)系的斜率和截距差異不顯著,因此,本研究將所有年份數(shù)據(jù)統(tǒng)一回歸,均達到極顯著水平(<0.0001),mass分別可以解釋max和cmax處理間變異的71.89%和88.76%(圖5)。

    圖5 2013-2016年不同處理冬小麥旗葉光合能力(Amax)及最大羧化速率(Vcmax)與葉片氮質(zhì)量分數(shù)(Nmass)之間的關(guān)系

    本研究中、高水及覆蓋處理的旗葉在生育末期仍能保持較高的max和cmax,與之較高的葉片N含量及葉片衰老的延緩有關(guān)。光合速率一般隨葉片N含量增加而增加,二者的定量關(guān)系在葉片N含量較低時為線性,而葉片N含量較高,且變化范圍較廣時則為指數(shù)增長到最大值的曲線關(guān)系[30]。本研究葉片N含量范圍內(nèi),光合參數(shù)與之關(guān)系仍為線性,未發(fā)現(xiàn)隨著葉片N含量的增加,光合參數(shù)增長變緩的現(xiàn)象。適宜的水分和覆蓋處理增加了冬小麥葉片N含量,高葉片N含量不僅與高光合色素-葉綠素含量相關(guān),而且可提高葉綠體有關(guān)光合碳同化酶類活性,進而提高光合速率[20,31]。

    2.4 產(chǎn)量與光合參數(shù)的關(guān)系

    光合作用是產(chǎn)量形成的基礎(chǔ),本研究分析了產(chǎn)量與光合參數(shù)之間的相關(guān)關(guān)系,產(chǎn)量與max及cmax之間均呈現(xiàn)出顯著的線性相關(guān)關(guān)系(圖6),所有年份產(chǎn)量均與max和cmax線性相關(guān)顯著,且不同年份之間相關(guān)關(guān)系的斜率和截距差異不顯著,因此,本研究將所有年份數(shù)據(jù)統(tǒng)一回歸,均達到極顯著水平(<0.000 1),max和cmax分別能解釋產(chǎn)量變異的83.15%和90.99%。該關(guān)系較直觀的闡明了光合參數(shù)對產(chǎn)量的決定作用,max和cmax較高的處理,產(chǎn)量也相應(yīng)的較高[32]。需要指出的是,本研究2013年度由于返青期初次灌水較晚,冬小麥返青期可能存在干旱脅迫,尤其是低水及不覆蓋處理該年度的產(chǎn)量低于其余年份,這種年際間的產(chǎn)量差異也在光合參數(shù)之間反映出來,并未使2013年度產(chǎn)量與max及cmax之間關(guān)系明顯偏離所有年份的回歸直線。

    圖6 2013-2016年不同處理冬小麥產(chǎn)量(Y)與旗葉最大光合速率(Amax)及最大羧化速率(Vcmax)之間的關(guān)系

    3 結(jié) 論

    本研究通過多年連續(xù)試驗確定了秸稈覆蓋處理和不同滴灌制度下的關(guān)鍵光合參數(shù),并分析了光合參數(shù)處理間差異來源及其與作物產(chǎn)量的關(guān)系。4個生長季的試驗表明,秸稈覆蓋處理和不同滴灌制度均不同程度的影響了旗葉關(guān)鍵光合參數(shù)的值,而二者的交互作用對光合參數(shù)的影響則均未發(fā)現(xiàn)。

    滴灌水量對表觀光量子效率的顯著影響主要表現(xiàn)在抽穗灌漿后期,而覆蓋對影響不顯著;覆蓋和滴灌水量對光合能力max的影響顯著并未表現(xiàn)出明顯季節(jié)差異;滴灌水量對最大羧化速率cmax和最大電子傳遞效率max的顯著影響主要表現(xiàn)在抽穗灌漿前、中期,而覆蓋對cmax和max的顯著影響則主要表現(xiàn)在抽穗灌漿中、后期。

    覆蓋處理和提高滴灌水量有利于提高、max、cmax和max,幾乎所有測定均以覆蓋高水處理(T5)的光合參數(shù)最大,除2016年第一次測定外,高水不覆蓋(T6)和中水覆蓋(T3)的相應(yīng)值略低但與T5差異不顯著,而低水不覆蓋(T2)的相應(yīng)值則均顯著低于T5。從降低灌水量的角度來講,T3的灌溉制度可以作為秸稈覆蓋冬小麥田的滴灌制度,即65%~85%的田間持水量作為土壤水分的上下限控制,該灌溉制度配合田間秸稈覆蓋措施,可以保持冬小麥旗葉較高的光合能力。

    處理間光合參數(shù)的差異可以用葉片N含量mass來解釋,所有年份、所有測定日期的max和cmax均與mass線性相關(guān)顯著。而光合參數(shù)又決定了產(chǎn)量的處理間差異,所有年份的產(chǎn)量均與max和cmax分別線性相關(guān)(2> 0.8)。因此,通過測定不同處理的mass,可以估算其旗葉關(guān)鍵光合參數(shù),進而實現(xiàn)其產(chǎn)量的估算。

    [1] 劉昌明,周長青,張士鋒,等. 小麥水分生產(chǎn)函數(shù)及其效益的研究[J]. 地理研究,2005,24(1):1-10. Liu Changming, Zhou Zhangqing, Zhang Shifeng, et al. Study on water production function and efficiency of wheat[J]. Geographical Research, 2005, 24(1): 1-10. (in Chinese with English abstract)

    [2] 房全孝,王建林,于舜章. 華北平原小麥-玉米兩熟制節(jié)水潛力與灌溉對策[J]. 農(nóng)業(yè)工程學報,2011,27(7):37-44. Fang Quanxiao, Wang Jianlin, Yu Shunzhang. Water-saving potential and irrigation strategies for wheat-maize double cropping system in the North China Plain[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2011, 27(7): 37-44. (in Chinese with English abstract)

    [3] Camp CR. Subsurface drip irrigation: A review[J]. Transactions of the ASABE, 1998, 41(5): 1353-1367.

    [4] 劉立晶,高煥文,李洪文. 玉米-小麥一年兩熟保護性耕作體系試驗研究[J]. 農(nóng)業(yè)工程學報,2004,20(3):70-73. Liu Lijing, Gao Huanwen, Li Hongwen. Conservation tillage for corn-wheat two crops a year region[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2004, 20(3): 70-73. (in Chinese with English abstract)

    [5] 楊永輝,武繼承,張玉亭,等. 耕作與保墑措施對小麥不同生育階段水分利用及產(chǎn)量的影響[J]. 華北農(nóng)學報,2016,31(3):184-190. Yang Yonghui, Wu Jicheng, Zhang Yuting, et al. Effects of tillage, moisture conservation on water use and yield in wheat at different growth stage[J]. Acta agriculturae boreali-sinica 2016, 31(3): 184-190. (in Chinese with English abstract)

    [6] 黃明,吳金芝,李友軍,等. 不同耕作方式對旱作區(qū)冬小麥生產(chǎn)和產(chǎn)量的影響[J]. 農(nóng)業(yè)工程學報,2009,25(1):50-54.Huang Ming, Wu Jinzhi, Li Youjun, et al. Effects of different tillage management on production and yield of winter wheat in dryland[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2009, 25(1): 50-54. (in Chinese with English abstract)

    [7] 劉青林,張恩和,王琦,等. 灌溉與施氮對留茬免耕春小麥耗水規(guī)律,產(chǎn)量和水分利用效率的影響[J]. 草業(yè)學報,2012,21(5):169. Liu Qinglin, Zhang Enhe, Wang Qi, et al. Effect of irrigation and nitrogen supply levels on water consumption, grain yield and water use efficiency of spring wheat on no-tillage with stubble standing farmland[J]. Acta prataculturae sinica, 2012, 21(5): 169. (in Chinese with English abstract)

    [8] Braune Henning, Müller Johannes, Diepenbrock Wulf. Integrating effects of leaf nitrogen, age, rank, and growth temperature into the photosynthesis-stomatal conductance model LEAFC3-N parameterised for barley (L.)[J]. Ecological Modelling, 2009, 220(13/14): 1599-1612.

    [9] Yin X, Struik PC. C3 and C4 photosynthesis models: An overview from the perspective of crop modelling[J]. NJAS Wageningen Journal of Life Sciences, 2009, 57(1): 27-38.

    [10] Müller Johannes, Wernecke Peter, Diepenbrock Wulf. LEAFC3-N: a nitrogen-sensitive extension of the CO2and H2O gas exchange model LEAFC3 parameterised and tested for winter wheat (L.)[J]. Ecological Modelling, 2005, 183: 183-210.

    [11] 江曉東,王蕓,侯連濤,等. 少免耕模式對冬小麥生育后期光合特性的影響[J]. 農(nóng)業(yè)工程學報,2006,22(5):66-69. Jiang Xiaodong, Wang Yun, Hou Liantao, et al. Effects of minimum tillage and no-tillage systems on photosynthetic characteristics at late growth stages of winter wheat[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2006, 22(5): 66-69. (in Chinese with English abstract)

    [12] Pal M, Rao LS, Jain V, et al. Effects of elevated CO2and nitrogen on wheat growth and photosynthesis[J]. Biologia Plantarum, 2005, 49(3): 467-470.

    [13] Foulkes MJ, Hawkesford MJ, Barraclough PB, et al. Identifying traits to improve the nitrogen economy of wheat: Recent advances and future prospects[J]. Field Crops Research, 2009, 114(3): 329-342.

    [14] 張彥群,王建東,龔時宏,等. 滴灌條件下冬小麥施氮增產(chǎn)的光合生理響應(yīng)[J]. 農(nóng)業(yè)工程學報,2015,31(6):170-177. Zhang Yanqun, Wang Jiandong, Gong Shihong, et al. Photosynthetic response of yield enhancement by nitrogen fertilization in winter wheat fields with drip irrigation[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2015, 31(6): 170-177. (in Chinese with English abstract)

    [15] Shangguan Zhouping, Shao Mingan, Dyckmans Jens. Effects of nitrogen nutrition and water deficit on net photosynthetic rate and chlorophyll fluorescence in winter wheat[J]. Journal of Plant Physiology, 2000, 156(1): 46-51.

    [16] 張向前,張賀飛,錢益亮. 不同秸稈覆蓋模式下小麥植株性狀、光合及產(chǎn)量的差異[J]. 麥類作物學報,2016,36(1):120-127. Zhang Xiangqian, Zhang Hefei, Qian Yiliang. Differences in plant traits, photosynthesis and yield of wheat under different straw mulching modes[J]. Journal of triticeae crops. 2016, 36(1): 120-127. (in Chinese with English abstract)

    [17] 馬冬云,郭天財,王晨陽,等. 施氮量對冬小麥灌漿期光合產(chǎn)物積累,轉(zhuǎn)運及分配的影響[J]. 作物學報,2008,34(6):1027-1033. Ma Dongyun, Guo Tiancai, Wang Chenyang, et al. Effects of nitrogen application rates on accumulation, translocation, and partitioning of photosynthate in winter wheat at grain filling stage[J]. Acta agronomica sinica. 2008, 34(6): 1027-1033.(in Chinese with English abstract)

    [18] Evans John R. Photosynthesis and nitrogen relationships in leaves of C3 plants[J]. Oecologia, 1989, 78(1): 9-19.

    [19] Evans John R. Nitrogen and photosynthesis in the flag leaf of wheat (L.)[J]. Plant physiology, 1983, 72(2): 297-302.

    [20] Li Dandan, Tian Mengyu, Cai Jian, et al. Effects of low nitrogen supply on relationships between photosynthesis and nitrogen status at different leaf position in wheat seedlings[J]. Plant Growth Regulation, 2013, 70(3): 257-263.

    [21] 王建東,龔時宏,高占義,等. 滴灌模式對農(nóng)田土壤水氮空間分布及冬小麥產(chǎn)量的影響[J]. 農(nóng)業(yè)工程學報,2009,25(11):68-73. Wang Jiandong, Gong Shihong, Gao Zhanyi, et al. Effects of drip irrigation mode on spatial distribution of soil water and nitrogen and winter wheat yield[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2009, 25(11): 68-73. (in Chinese with English abstract)

    [22] Thornley John HM. Mathematical models in plant physiology[M]. Academic Press (Inc.) London, Ltd., 1976.

    [23] Sharkey Thomas D. Photosynthesis in intact leaves of C3 plants: Physics, Physiology and rate limitations[J]. Botanical Review, 1985, 51(1): 53-105.

    [24] 陸佩玲,于強. 冬小麥光合作用的光響應(yīng)曲線的擬合[J]. 中國農(nóng)業(yè)氣象,2001,22(2):12-14.Lu Peiling, Yu Qiang. Fitting light response curves of photosynthesis of winter wheat[J]. Chinese Journal of Agrometeorology, 2001, 22(2): 12-14. (in Chinese with English abstract)

    [25] 馬東輝,趙長星,王月福,等. 施氮量和花后土壤含水量對小麥旗葉光合特性和產(chǎn)量的影響[J]. 生態(tài)學報,2008,28(10):4896-4901.Ma Donghui, Zhao Zhangxing, Wang Yuefu, et al. Effects of nitrogen fertilizer rate and post-anthesis soil water content on photosynthetic characteristics in flag leaves and yield of wheat[J]. Acta Ecologica Sinica, 2008, 28(10): 4896-4901. (in Chinese with English abstract)

    [26] Wullschleger SD. Biochemical limitations to carbon assimilation in C3 plants-a retrospective analysis of the A/Ci curves from 109 species[J]. Journal of Experimental Botany, 1992, 44(2): 907-920.

    [27] Müller Johannes, Eschenr?der André, Christen Olaf. LEAFC3-N photosynthesis, stomatal conductance, transpiration and energy balance model: Finite mesophyll conductance, drought stress, stomata ratio, optimized solution algorithms, and code[J]. Ecological Modelling, 2014, 290(4): 134-145.

    [28] Leuning R. Temperature dependence of two parameters in a photosynthesis model[J]. Plant Cell & Environment, 2002, 25(9): 1205-1210.

    [29] Evans J, Farquhar GD. Modeling canopy photosynthesis from the biochemistry of C3 chloroplast[M]. In: Boote KJ, Loomis RS, eds. Modelling Crop Photosynthesis-from Biochemistry to Canopy. CSSA special publication no. 19. Madison, WI: American Society of Agronomy, Inc., 1991, 1-15.

    [30] Bindraban PS. Impact of canopy nitrogen profile in wheat on growth[J]. Field Crops Research, 1999, 63(1): 63-77.

    [31] Serrago Román A, Alzueta Ignacio, Savin Roxana, et al. Understanding grain yield responses to source–sink ratios during grain filling in wheat and barley under contrasting environments[J]. Field Crops Research, 2013, 150: 42-51.

    [32] Zhang Yanqun, Wang Jiandong, Gong Shihong, et al. Nitrogen fertigation effect on photosynthesis, grain yield and water use efficiency of winter wheat[J]. Agricultural Water Management, 2017, 179: 277-287.

    Effects of straw mulching and drip irrigation scheduling on photosynthetic characteristic and yield of winter wheat

    Zhang Yanqun1, Wang Jiandong1※, Gong Shihong1, Xu Di1, Sun Jiwen1,2

    (1.100048,;2.255049,)

    Determining the photosynthetic parameters is the basis for understanding the physiological causes of water saving and yield increasing under drip irrigation modes. In this study, intrinsic photosynthetic parameters were determined during four years experiments in winter wheat field treated with six drip irrigation modes (T1 to T6, stands for treatments of lower irrigation amount with straw mulch, lower irrigation amount without straw mulch, moderate irrigation amount with straw mulch, moderate irrigation amount without straw mulch, higher irrigation amount with straw mulch, and higher irrigation amount without straw mulch, respectively). The differences in photosynthetic parameters among treatments were statistically analyzed, and the correlations between photosynthetic parameters and crop yield were analyzed. The results showed that increasing irrigation amount and straw mulch enhanced the apparent photon quantum efficiency (), photosynthetic capacity (max), maximum carboxylation rate (cmax) and maximum electron transfer rate (max), respectively. However, the interaction effects of irrigation amount and mulching on photosynthetic parameters were not significant (> 0.1). The effect of drip irrigation on apparent quantum efficiency () was significant (< 0.05) during the late periods of the heading and grain filling stages, while the effect of straw mulch onwas not significant (> 0.1). The effect of straw mulch and drip irrigation onmaxshow significant (< 0.05) seasonal differences. The effects of drip irrigation oncmaxandmaxwere significant (< 0.05) during the mid-season periods of the heading and grain filling stages, while the effects of straw mulch were significant (< 0.05) during the later periods of the heading and grain filling stages. The results of multiple comparisons showed that the photosynthetic parameters of T5 were the highest (< 0.05). Except for the measurements in 2016, the corresponding values of T6 and T3 were slightly lower but with no significant difference from those of T5 (> 0.1). However, the corresponding values of T2 were significantly lower (< 0.05) than those of T5. These results suggested that, under higher irrigation amount, photosynthetic parameters were not affected significantly (> 0.05) by mulching, but under lower irrigation amount, the parameters significantly (< 0.05) decreased if treated with no mulching. The average values ofmaxduring the four years from T1 to T6 were 29.59, 26.45, 31.44, 28.68, 32.42 and 30.90mol/ (m2·s), respectively. The average values ofcmaxduring the four years from T1 to T6 were 121.42, 107.68, 133.42, 113.55, 141.83 and 129.95mol/ (m2·s), respectively. From the point of view of reducing the amount of irrigation, T3 irrigation scheduling can be used as an optimal strategy for winter wheat under drip irrigation with straw mulch, indicating that 65% to 85% of the field capacity was used as the lower and upper limits of soil moisture control. The management practices of drip irrigation with straw mulch can keep the higher photosynthetic capacity of flag leaf for winter wheat. The differences inmaxandcmaxamong treatments could be explained by the leaf nitrogen content (mass). Grain yield for all years were also linearly related tomaxandcmax, respectively. Therefore, by measuringmass, one can estimatemaxandcmax, thus achieving the yield estimates of winter wheat field in this area. In addition, the photosynthetic characteristic parameters can also be applied to the crop model to improve the accuracy of model prediction under different water and mulch regimes.

    crops; irrigation; photosynthesis; winter wheat; straw mulching; leaf nitrogen content

    10.11975/j.issn.1002-6819.2017.12.021

    S161.4

    A

    1002-6819(2017)-12-0162-08

    2016-11-03

    2017-03-21

    國家“十二五”科技支撐計劃課題(2014BAD12B05);中國水科院基本科研業(yè)務(wù)費項目(ID0145B602017);國家自然科學基金項目(51309250)

    張彥群,女,河北冀州人,博士,高級工程師,主要從事農(nóng)田水循環(huán)與作物高效用水生理調(diào)控研究。Email:zhangyq@iwhr.com

    王建東,男,湖南茶陵人,博士,教授級高級工程師,主要從事節(jié)水灌溉原理與技術(shù)研究。Email:wangjd@iwhr.com

    張彥群,王建東,龔時宏,許 迪,孫繼文. 秸稈覆蓋和滴灌制度對冬小麥光合特性和產(chǎn)量的影響[J]. 農(nóng)業(yè)工程學報,2017,33(12):162-169. doi:10.11975/j.issn.1002-6819.2017.12.021 http://www.tcsae.org

    Zhang Yanqun, Wang Jiandong, Gong Shihong, Xu Di, Sun Jiwen. Effects of straw mulching and drip irrigation scheduling on photosynthetic characteristic and yield of winter wheat[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(12): 162-169. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2017.12.021 http://www.tcsae.org

    猜你喜歡
    羧化旗葉冬小麥
    不同品種小麥灌漿期旗葉光合特性及光合基因表達對臭氧濃度升高的響應(yīng)
    作物學報(2022年9期)2022-07-18 00:58:42
    水分虧缺對小麥芒和旗葉光合特性及蔗糖、淀粉合成的影響
    旗葉衰老產(chǎn)量差異 可作小麥優(yōu)選依據(jù)
    抽穗后不同時期去除旗葉對不同穗型小麥產(chǎn)量的影響
    鈣離子體系下羧化殼聚糖對滑石浮選行為影響
    甘肅冬小麥田
    植物保護(2017年1期)2017-02-13 06:44:34
    冬小麥和春小麥
    中學生(2015年4期)2015-08-31 02:53:50
    機械活化法制備順丁烯二酸酐羧化淀粉的助洗性能
    冬小麥——新冬18號
    冬小麥—新冬41號
    国产一区二区亚洲精品在线观看| 精品午夜福利在线看| 中国国产av一级| 久久国内精品自在自线图片| 日日啪夜夜爽| 成人欧美大片| 91狼人影院| 国产成人freesex在线| 91久久精品国产一区二区三区| 国产精品女同一区二区软件| av在线播放精品| 日本欧美国产在线视频| 麻豆精品久久久久久蜜桃| 中文字幕免费在线视频6| 久久久久久久亚洲中文字幕| 美女被艹到高潮喷水动态| 欧美+日韩+精品| 精品国产露脸久久av麻豆| 人妻夜夜爽99麻豆av| 国产 一区精品| 国产一区二区亚洲精品在线观看| 欧美成人一区二区免费高清观看| 亚洲精品日本国产第一区| 日本黄色片子视频| 国产精品久久久久久久久免| 波多野结衣巨乳人妻| 成人免费观看视频高清| 三级国产精品欧美在线观看| 亚洲精品成人久久久久久| 日日撸夜夜添| 久久精品久久久久久噜噜老黄| 亚洲真实伦在线观看| 久久久午夜欧美精品| 国产极品天堂在线| 亚洲成人精品中文字幕电影| 久久精品久久久久久噜噜老黄| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 人人妻人人澡人人爽人人夜夜| 亚洲色图av天堂| 国产白丝娇喘喷水9色精品| 亚洲成人av在线免费| 三级国产精品欧美在线观看| 亚洲欧美日韩无卡精品| 菩萨蛮人人尽说江南好唐韦庄| 午夜亚洲福利在线播放| 久久99热6这里只有精品| 欧美潮喷喷水| 身体一侧抽搐| 国产亚洲午夜精品一区二区久久 | 久久久久久九九精品二区国产| 一级毛片 在线播放| 欧美zozozo另类| 亚洲成人久久爱视频| 如何舔出高潮| 精品久久久久久久久亚洲| 亚洲精品日韩在线中文字幕| 亚洲精品,欧美精品| 两个人的视频大全免费| 国产 一区精品| av又黄又爽大尺度在线免费看| 日日啪夜夜撸| 欧美bdsm另类| 亚洲欧洲国产日韩| 听说在线观看完整版免费高清| 亚洲综合精品二区| 51国产日韩欧美| 午夜精品国产一区二区电影 | 午夜视频国产福利| 亚洲最大成人中文| 大话2 男鬼变身卡| 亚洲精品国产av蜜桃| 国产 精品1| 美女cb高潮喷水在线观看| 精品人妻熟女av久视频| 国产中年淑女户外野战色| 成年人午夜在线观看视频| 午夜福利视频1000在线观看| 免费黄色在线免费观看| 久久精品国产亚洲av天美| 天堂网av新在线| 欧美日本视频| 亚洲经典国产精华液单| 国产高潮美女av| 久久99热这里只有精品18| 一级毛片电影观看| 卡戴珊不雅视频在线播放| 精品一区二区免费观看| 久久精品熟女亚洲av麻豆精品| 观看免费一级毛片| 亚洲在久久综合| 在线看a的网站| 99视频精品全部免费 在线| 国产熟女欧美一区二区| 欧美日韩综合久久久久久| 婷婷色麻豆天堂久久| 久久久成人免费电影| 国产精品爽爽va在线观看网站| 久久国内精品自在自线图片| 成人二区视频| 男人添女人高潮全过程视频| 亚洲经典国产精华液单| 婷婷色麻豆天堂久久| 久久国产乱子免费精品| 大片电影免费在线观看免费| 六月丁香七月| 秋霞在线观看毛片| 国产综合懂色| 国产成人freesex在线| 国产成人精品久久久久久| 久久国产乱子免费精品| 中文精品一卡2卡3卡4更新| 精品少妇久久久久久888优播| 免费黄色在线免费观看| 高清在线视频一区二区三区| 免费观看在线日韩| 久久久国产一区二区| 国产人妻一区二区三区在| 久热这里只有精品99| 亚洲精品自拍成人| 国产精品久久久久久av不卡| 一级爰片在线观看| 成人毛片a级毛片在线播放| 日本猛色少妇xxxxx猛交久久| 亚洲成人久久爱视频| 91精品一卡2卡3卡4卡| 日本与韩国留学比较| 一级二级三级毛片免费看| 人妻系列 视频| 中文资源天堂在线| 精品久久久久久电影网| 久久久久精品性色| 国产精品不卡视频一区二区| 亚洲天堂av无毛| av卡一久久| 国产大屁股一区二区在线视频| 免费观看av网站的网址| 国产男人的电影天堂91| 人妻少妇偷人精品九色| 国产毛片在线视频| 精品久久久久久电影网| 亚洲国产精品成人久久小说| 欧美激情国产日韩精品一区| 免费av毛片视频| 禁无遮挡网站| 大又大粗又爽又黄少妇毛片口| 午夜福利高清视频| 超碰av人人做人人爽久久| 欧美日韩亚洲高清精品| 夫妻性生交免费视频一级片| 久久久久久久午夜电影| 18禁动态无遮挡网站| 波野结衣二区三区在线| 国产精品人妻久久久影院| 国产片特级美女逼逼视频| 可以在线观看毛片的网站| 蜜臀久久99精品久久宅男| 一级二级三级毛片免费看| 亚洲精品国产av蜜桃| 熟女电影av网| 欧美一区二区亚洲| 在线观看人妻少妇| 神马国产精品三级电影在线观看| 色视频在线一区二区三区| 精品99又大又爽又粗少妇毛片| 精品国产露脸久久av麻豆| 日韩av免费高清视频| 日韩一区二区三区影片| 熟女电影av网| 午夜免费鲁丝| 高清在线视频一区二区三区| 日本与韩国留学比较| 亚洲国产精品成人综合色| 国产 一区 欧美 日韩| 国产永久视频网站| 天天躁日日操中文字幕| 另类亚洲欧美激情| 日日撸夜夜添| 在线免费观看不下载黄p国产| 国产探花极品一区二区| 亚洲国产欧美在线一区| 极品教师在线视频| 国产精品一区二区性色av| 激情 狠狠 欧美| 亚洲久久久久久中文字幕| 亚洲欧美精品自产自拍| 女人十人毛片免费观看3o分钟| 又爽又黄a免费视频| 国内少妇人妻偷人精品xxx网站| 天美传媒精品一区二区| 天堂网av新在线| 汤姆久久久久久久影院中文字幕| 亚洲国产最新在线播放| 国产成人一区二区在线| 国产亚洲91精品色在线| 99热这里只有精品一区| 亚洲激情五月婷婷啪啪| 91精品伊人久久大香线蕉| 中文字幕免费在线视频6| www.av在线官网国产| 只有这里有精品99| 国产视频内射| 中文资源天堂在线| 日产精品乱码卡一卡2卡三| 听说在线观看完整版免费高清| 国产精品国产三级国产专区5o| 欧美老熟妇乱子伦牲交| 国产极品天堂在线| 嫩草影院精品99| 国产精品一二三区在线看| 国产精品99久久久久久久久| av免费观看日本| 一区二区三区四区激情视频| 女的被弄到高潮叫床怎么办| 成人免费观看视频高清| 九九久久精品国产亚洲av麻豆| 欧美xxxx黑人xx丫x性爽| 久久精品国产a三级三级三级| 欧美三级亚洲精品| 97在线视频观看| 久久精品夜色国产| 成人欧美大片| 免费看a级黄色片| 最近中文字幕高清免费大全6| 在线免费观看不下载黄p国产| 人妻 亚洲 视频| 成人毛片60女人毛片免费| 亚洲国产精品专区欧美| 卡戴珊不雅视频在线播放| 麻豆精品久久久久久蜜桃| 国产淫语在线视频| 国产成人精品久久久久久| 久久久国产一区二区| 在线免费十八禁| 成人欧美大片| 人妻夜夜爽99麻豆av| 午夜激情久久久久久久| 男人和女人高潮做爰伦理| 国产色爽女视频免费观看| 大话2 男鬼变身卡| 国产男女内射视频| 下体分泌物呈黄色| 国产探花在线观看一区二区| 日本黄大片高清| 高清av免费在线| 国产精品99久久99久久久不卡 | 国产成人一区二区在线| 麻豆国产97在线/欧美| 国产视频内射| 我的老师免费观看完整版| 国产av国产精品国产| 美女被艹到高潮喷水动态| 少妇人妻 视频| 亚洲精品第二区| 九九在线视频观看精品| 亚洲,一卡二卡三卡| 联通29元200g的流量卡| 免费av毛片视频| 国产精品伦人一区二区| 成人毛片60女人毛片免费| 国产成人免费无遮挡视频| 80岁老熟妇乱子伦牲交| 97超视频在线观看视频| 亚洲人成网站在线播| 亚洲精品自拍成人| 国产片特级美女逼逼视频| 白带黄色成豆腐渣| 我的女老师完整版在线观看| 亚洲国产成人一精品久久久| 日本与韩国留学比较| 亚洲最大成人手机在线| 在线a可以看的网站| 国产精品精品国产色婷婷| 亚洲在线观看片| 免费观看a级毛片全部| 精品一区二区三区视频在线| 最近中文字幕2019免费版| 少妇裸体淫交视频免费看高清| 简卡轻食公司| 国产一区有黄有色的免费视频| 在线a可以看的网站| 七月丁香在线播放| 欧美成人精品欧美一级黄| 久久久久久九九精品二区国产| videossex国产| 久久久久久久久久成人| 日韩强制内射视频| 免费人成在线观看视频色| 精品一区二区免费观看| 亚洲图色成人| 国产一区二区在线观看日韩| 亚洲成人一二三区av| 亚洲精品乱码久久久久久按摩| 久久午夜福利片| 国产伦理片在线播放av一区| 男女边吃奶边做爰视频| 久久这里有精品视频免费| 少妇人妻 视频| 69人妻影院| 亚洲人成网站在线观看播放| 精品久久国产蜜桃| 中文字幕久久专区| 在线观看三级黄色| 夜夜爽夜夜爽视频| 亚洲怡红院男人天堂| 国产精品久久久久久精品古装| av网站免费在线观看视频| 日韩欧美精品免费久久| 特级一级黄色大片| 久久精品综合一区二区三区| 中文字幕久久专区| 2022亚洲国产成人精品| 久久ye,这里只有精品| 尤物成人国产欧美一区二区三区| 在线播放无遮挡| 亚洲国产精品专区欧美| 国产精品熟女久久久久浪| 日日摸夜夜添夜夜添av毛片| av国产免费在线观看| 99久国产av精品国产电影| 欧美变态另类bdsm刘玥| 日本色播在线视频| 久久综合国产亚洲精品| 国产精品一区www在线观看| 亚洲av成人精品一区久久| 日韩av在线免费看完整版不卡| 日本爱情动作片www.在线观看| 国产色爽女视频免费观看| 亚洲精品自拍成人| 69人妻影院| 国产淫语在线视频| 欧美亚洲 丝袜 人妻 在线| 69av精品久久久久久| 蜜桃久久精品国产亚洲av| 日韩,欧美,国产一区二区三区| 国产精品成人在线| 国产免费视频播放在线视频| 好男人在线观看高清免费视频| 亚洲人成网站高清观看| 久久精品人妻少妇| av在线播放精品| 又爽又黄a免费视频| 国内少妇人妻偷人精品xxx网站| 少妇裸体淫交视频免费看高清| 新久久久久国产一级毛片| 边亲边吃奶的免费视频| 成人漫画全彩无遮挡| 免费黄网站久久成人精品| 成人毛片60女人毛片免费| 在线观看美女被高潮喷水网站| 欧美高清性xxxxhd video| 亚洲国产高清在线一区二区三| 色哟哟·www| 成年免费大片在线观看| 久久女婷五月综合色啪小说 | 97超视频在线观看视频| 97精品久久久久久久久久精品| 日韩国内少妇激情av| 亚洲精品色激情综合| 亚洲精品第二区| 大陆偷拍与自拍| 91久久精品国产一区二区三区| 亚洲精品乱码久久久v下载方式| 一个人看视频在线观看www免费| 国产探花极品一区二区| 草草在线视频免费看| 亚洲av一区综合| 国产大屁股一区二区在线视频| 亚洲精品乱码久久久久久按摩| 男女边吃奶边做爰视频| 久热久热在线精品观看| av网站免费在线观看视频| 女人十人毛片免费观看3o分钟| 国产免费又黄又爽又色| av在线观看视频网站免费| 国产高清国产精品国产三级 | 97精品久久久久久久久久精品| 夜夜爽夜夜爽视频| 午夜福利高清视频| 国产精品秋霞免费鲁丝片| 欧美日韩国产mv在线观看视频 | 日韩视频在线欧美| 久久久午夜欧美精品| 日本av手机在线免费观看| 涩涩av久久男人的天堂| 日韩av在线免费看完整版不卡| 欧美激情在线99| 成年免费大片在线观看| 国产有黄有色有爽视频| 一级毛片久久久久久久久女| 久久久a久久爽久久v久久| 26uuu在线亚洲综合色| 亚洲av男天堂| 国产亚洲5aaaaa淫片| 亚洲精品久久午夜乱码| 在线亚洲精品国产二区图片欧美 | av一本久久久久| 一级片'在线观看视频| 在线免费观看不下载黄p国产| 联通29元200g的流量卡| 婷婷色麻豆天堂久久| 综合色av麻豆| 亚洲精品色激情综合| 毛片女人毛片| 亚洲国产精品999| 久久久久久九九精品二区国产| 有码 亚洲区| 久久精品国产亚洲网站| av在线蜜桃| 国产高清国产精品国产三级 | 亚洲精品亚洲一区二区| 亚洲av日韩在线播放| 亚洲在久久综合| 亚洲怡红院男人天堂| 又大又黄又爽视频免费| 国产高清有码在线观看视频| 欧美少妇被猛烈插入视频| 一级毛片久久久久久久久女| 国产爽快片一区二区三区| av在线蜜桃| 一级毛片电影观看| 又爽又黄a免费视频| 22中文网久久字幕| 国产免费一区二区三区四区乱码| 国产精品嫩草影院av在线观看| 白带黄色成豆腐渣| 久久久久精品久久久久真实原创| 亚洲综合精品二区| 亚洲成人久久爱视频| 欧美日韩亚洲高清精品| 一本一本综合久久| 国产日韩欧美在线精品| 岛国毛片在线播放| 亚洲内射少妇av| 噜噜噜噜噜久久久久久91| 久久精品国产a三级三级三级| 波野结衣二区三区在线| 精品久久久久久电影网| 校园人妻丝袜中文字幕| av黄色大香蕉| 久久女婷五月综合色啪小说 | 日本与韩国留学比较| 国产成人免费观看mmmm| 少妇人妻精品综合一区二区| 青春草视频在线免费观看| 在线观看av片永久免费下载| 天美传媒精品一区二区| 一个人观看的视频www高清免费观看| 另类亚洲欧美激情| 国产永久视频网站| av国产精品久久久久影院| 99热这里只有精品一区| 只有这里有精品99| 国产黄频视频在线观看| 黄色怎么调成土黄色| 亚洲精品国产av成人精品| 高清欧美精品videossex| 亚洲经典国产精华液单| 精品久久久久久久人妻蜜臀av| 国产 一区精品| 少妇人妻 视频| 欧美一区二区亚洲| 久久精品夜色国产| 一本一本综合久久| 一级av片app| 中文字幕av成人在线电影| 永久网站在线| 2021天堂中文幕一二区在线观| 国精品久久久久久国模美| 美女cb高潮喷水在线观看| 国产精品国产av在线观看| 欧美另类一区| 天堂网av新在线| 丰满乱子伦码专区| 内射极品少妇av片p| 美女国产视频在线观看| 久久精品熟女亚洲av麻豆精品| 亚洲天堂av无毛| 精品一区二区三区视频在线| tube8黄色片| 女人被狂操c到高潮| 天天躁日日操中文字幕| 亚洲欧美一区二区三区国产| 欧美激情久久久久久爽电影| 夫妻午夜视频| 成人美女网站在线观看视频| av播播在线观看一区| 好男人在线观看高清免费视频| 亚洲av中文字字幕乱码综合| 五月伊人婷婷丁香| 少妇 在线观看| .国产精品久久| 亚洲精品国产色婷婷电影| 国产色爽女视频免费观看| 国产欧美另类精品又又久久亚洲欧美| 下体分泌物呈黄色| 三级国产精品欧美在线观看| 久久久精品欧美日韩精品| 成人欧美大片| 亚洲精品亚洲一区二区| 国产人妻一区二区三区在| 伦理电影大哥的女人| 色吧在线观看| 毛片女人毛片| 三级国产精品欧美在线观看| 大香蕉97超碰在线| 国产精品.久久久| 国产av国产精品国产| 黄色欧美视频在线观看| 亚洲av男天堂| 天堂网av新在线| 精品一区二区免费观看| 亚洲va在线va天堂va国产| 午夜福利视频精品| 直男gayav资源| 久久精品久久久久久噜噜老黄| 亚洲不卡免费看| 午夜福利高清视频| 久久久久性生活片| 婷婷色综合大香蕉| 建设人人有责人人尽责人人享有的 | 亚洲经典国产精华液单| 欧美国产精品一级二级三级 | 在线精品无人区一区二区三 | 大话2 男鬼变身卡| 午夜亚洲福利在线播放| 天天一区二区日本电影三级| 亚洲国产精品专区欧美| 中文字幕亚洲精品专区| 亚洲成人av在线免费| 国产色婷婷99| 亚洲三级黄色毛片| 精品久久久久久电影网| 简卡轻食公司| 两个人的视频大全免费| 日日摸夜夜添夜夜爱| 国产成人精品福利久久| 欧美另类一区| 日韩一区二区视频免费看| 国产精品三级大全| 蜜桃亚洲精品一区二区三区| 日韩一区二区三区影片| 久久久午夜欧美精品| 久久影院123| 秋霞伦理黄片| 18禁在线播放成人免费| 一个人看视频在线观看www免费| 欧美激情久久久久久爽电影| 欧美xxxx黑人xx丫x性爽| 极品少妇高潮喷水抽搐| 国产精品福利在线免费观看| 亚洲美女搞黄在线观看| 精品国产露脸久久av麻豆| 韩国高清视频一区二区三区| 2021天堂中文幕一二区在线观| 国产精品一二三区在线看| 免费不卡的大黄色大毛片视频在线观看| 人人妻人人爽人人添夜夜欢视频 | 大片免费播放器 马上看| 婷婷色av中文字幕| 国产白丝娇喘喷水9色精品| 国产男女超爽视频在线观看| 麻豆国产97在线/欧美| 国产伦精品一区二区三区四那| 亚洲精品自拍成人| 中文欧美无线码| 欧美国产精品一级二级三级 | 亚洲国产日韩一区二区| 精品一区二区免费观看| 青春草国产在线视频| 看免费成人av毛片| 久久久久久国产a免费观看| 特大巨黑吊av在线直播| 亚洲国产高清在线一区二区三| 久久久精品欧美日韩精品| 亚洲人成网站高清观看| 3wmmmm亚洲av在线观看| 国产老妇伦熟女老妇高清| 国产精品99久久久久久久久| 日韩制服骚丝袜av| 久久久久久国产a免费观看| 国产爽快片一区二区三区| 九九在线视频观看精品| 99九九线精品视频在线观看视频| 亚洲人成网站高清观看| 久久精品熟女亚洲av麻豆精品| av天堂中文字幕网| 99视频精品全部免费 在线| 乱系列少妇在线播放| 欧美高清成人免费视频www| 国产成人免费观看mmmm| 最近2019中文字幕mv第一页| 国产永久视频网站| 欧美xxxx性猛交bbbb| 伦精品一区二区三区| 在线观看一区二区三区| 韩国av在线不卡| 免费观看在线日韩| 国产91av在线免费观看| 在线观看国产h片| 国产高清三级在线| 超碰av人人做人人爽久久| 最近中文字幕高清免费大全6| 国产精品女同一区二区软件| 亚洲最大成人中文| 国产精品国产三级国产av玫瑰| 一级黄片播放器| 中文欧美无线码| 精品少妇久久久久久888优播| 久久鲁丝午夜福利片| 国内揄拍国产精品人妻在线| 亚洲精品久久午夜乱码| 伦精品一区二区三区| av免费在线看不卡| 五月天丁香电影| 一个人看视频在线观看www免费|