999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

施氮水平對甜高粱主要農藝性狀及其與干物質產量相關關系的影響

2016-09-01 07:16:36渠暉程亮陳俊峰陸曉燕沈益新
草業學報 2016年6期
關鍵詞:產量水平

渠暉,程亮,陳俊峰,陸曉燕,沈益新

(南京農業大學動物科技學院,江蘇 南京 210095)

?

施氮水平對甜高粱主要農藝性狀及其與干物質產量相關關系的影響

渠暉,程亮,陳俊峰,陸曉燕,沈益新*

(南京農業大學動物科技學院,江蘇 南京 210095)

為探求甜高粱干物質產量與主要農藝性狀的關系以及不同施氮水平對兩者相互關系的影響,本研究以大力士甜高粱為材料,采用完全隨機區組設計進行兩年田間試驗,測定不同施氮水平(0,100,200,300,400,500 kg N/hm2)下甜高粱的株高、莖粗、綠葉數、莖葉比、鮮干比、葉面積指數和干物質產量7個參數,并對上述參數進行了相關性分析、多元回歸分析和通徑分析。試驗結果表明,株高,葉面積指數和莖粗均與干物質產量呈極顯著正相關關系(P<0.01),相關系數分別為r=0.921**,r=0.865**,r=0.758**。通徑分析表明,株高是決定干物質產量的主要農藝性狀,莖葉比和葉面積指數與干物質產量的關系對不同施氮水平有不同響應。較低施氮水平(0,100 kg N/hm2)下,莖葉比對干物質產量的直接作用最大(P=0.925),葉面積指數對干物質產量的直接作用最小(P=0.162)。中等施氮水平(200,300 kg N/hm2)下,莖葉比對干物質產量的影響降低(P=0.073),葉面積指數對干物質產量的直接作用增加(P=0.797)。較高施氮水平(400,500 kg N/hm2)下,莖葉比和葉面積指數對干物質產量的直接作用均為負值(P=-0.125,P=-0.040)。施氮量為200~400 kg N/hm2時,干物質產量較高。綜上可知,在較高施氮條件下,合理密植是獲得較高干物質產量的關鍵。當供氮水平較低時,應注意收獲時間的選擇,以便更好地兼顧干物質產量和飼草營養價值。

甜高粱;施氮水平;干物質產量;農藝性狀;通徑分析

氮素是影響飼草產量和營養價值的重要營養元素之一[1-3],是飼草生長發育的主要限制因子。大量研究表明,施用氮肥能顯著提高飼草產量[4]、改善營養價值[5]。對禾本科飼草而言,因其不具備固氮能力,如果土壤中氮含量不足,則只能靠外界供氮來滿足自身生長所需。但不同地區和不同飼草品種對氮肥的響應差異巨大。Marsalis等[6]研究表明,傳統飼用高粱(Sorghumbicolor)在施氮量為218和291 kg N/hm2時,產量最高;而褐色中脈高粱品種的最適施氮量為106~140 kg N/hm2。Almodares等[7]研究認為,高粱產量達到最大值的適宜需氮量為75 kg N/hm2。丁成龍等[8]對美洲狼尾草(Pennisetumamericanum)的研究表明,鮮、干草產量均隨氮量的增加而顯著提高,施氮量在0~450 kg N/hm2范圍內,其鮮、干草產量與施氮量呈線性關系(r=0.9994)。韋家少等[9]研究表明,施氮水平為150和300 kg/hm2時,堅尼草(Panicummaximun)的干草產量顯著高于不施氮處理。曹文俠等[10]連續兩年的施氮試驗表明,高寒地區15齡人工草地的最佳施氮量為150 kg/hm2。因此,根據不同地區氣候條件,土壤特點,結合所選飼草品種,確定最佳的施氮水平,是飼草高產栽培技術的關鍵。

甜高粱具有抗逆性強[11],適應性廣[12-13],生物量大[14]等優良特性,不僅是世界第五大糧食作物,也是優質的飼料作物[15],已成為許多國家優質飼草的重要來源之一。國內外學者對甜高粱選育[16-17],栽培[18-19]和利用[20-21]等做了大量研究。有關施氮對甜高粱產量和品質的影響也有很多研究,但大多集中于干旱半干旱地區[22-23],而在較濕潤氣候條件下的研究鮮見報道,且對產量與其構成因子的相互關系以及不同施氮水平對兩者相互關系影響的研究較少。

我國華東農區氣候潮濕,土壤粘重,地下水位高,許多優質飼草的生長受到限制,探索研究適合當地氣候條件的優質飼草,是推進我國南方草業發展,實現畜牧業現代化的關鍵。本研究以極晚熟甜高粱品種為研究材料,通過測定不同施氮水平下甜高粱株高、莖粗、綠葉數、莖葉比、鮮干比、葉面積指數和干物質產量,就甜高粱干物質產量和主要農藝性狀之間的關系進行相關分析,多元回歸分析和通徑分析,解析不同施氮水平下影響甜高粱干物質產量的主要農藝性狀,以期為氮肥的合理施用和甜高粱的高產、優質栽培提供理論依據。

1 材料與方法

1.1試驗區概況

試驗于2009和2010年在江蘇省常州春暉乳業有限公司農場(31°51′ N,119°49′ E)進行。該試驗地屬于北亞熱帶季風性濕潤氣候,土壤類型為粘壤土。試驗地土壤具體特性及氣象條件見表1。

1.2試驗材料與設計

試驗材料為甜高粱品種大力士,由百綠(北京)國際草業有限公司提供。試驗采用完全隨機區組設計,共設6個氮肥處理(0,100,200,300,400,500 kg N/hm2),3次重復,18個小區。小區面積為12 m2(4 m×3 m),行距30 cm,株距25 cm,每小區種植10行。試驗地播前深翻20 cm,并施入腐熟的牛糞75 t/hm2(干物質含量為25%,總氮含量為0.55% DM)作基肥。播種方式為人工播種,每穴播3~4粒,3葉期按設計株距定苗。氮肥以純氮量換算成尿素分兩次施入,播種時施入1/3尿素,6葉期結合中耕追施余下2/3尿素。耕前噴施2.5 L/hm2阿特拉津并配以播后1周噴施2, 4-D用來防除雜草。兩年的播種時間見表1。

1.3測定內容與方法

由于大力士甜高粱是極晚熟品種,在長江中下游地區不能抽穗,因此按播種后的生長天數進行取樣。播種后46~50 d,大力士進入拔節期后每隔20 d取樣1次,總共取樣4次,具體取樣時間見表1。小區一半用于中期取樣,一半用于收獲測產。取樣時每小區隨機選取長勢一致的3株植株收割地上部分,分離莖葉(包括葉鞘),稱鮮重,用AM-300便攜式葉面積儀測量單株葉面積,計算葉面積指數(leaf area index,LAI)。將樣品帶回實驗室置于105℃鼓風干燥箱下殺青30 min,再于65℃下烘干至恒重,稱重,計算莖葉比,鮮干比。植株各部分粉碎均勻混合后過1 mm 篩,裝入自封袋保存供品質分析。取樣的同時也在每小區內選取長勢一致的5株植株進行掛牌,測量株高(拉直植株,測量從根部到最高處的垂直高度,取平均值,單位:m)、莖粗(以基部以上第三節莖粗的平均值表示,單位:cm)、綠葉數。在測產的一半小區中間收獲長為1.5 m,寬為相鄰3行面積上的植株,測定鮮重,根據鮮干比和收獲面積計算干物質產量(kg/hm2)。

表1 試驗地狀態、土壤特性和取樣時間Table 1 Plot and soil characteristics, and sampling date at each site and year

1.4數據統計

利用Microsoft Excel 2007處理基礎數據及作圖,采用SPSS 19.0軟件對影響干物質產量的各因素進行相關性分析、多元回歸分析和通徑分析。

2 結果與分析

2.1施氮水平對甜高粱農藝性狀的影響

作物的株高隨生育期的推進而增加。從圖1和圖2可以看出,2年試驗中,在300和400 kg N/hm2處理下,甜高粱各生長階段的株高均維持了較高水平,大于其他處理;收獲時,300 kg N/hm2處理的株高最大,分別為3.09 m(2009年)和2.99 m(2010年),500 kg N/hm2處理下的株高在整個生育期均屬于中間水平。在2009年的試驗中,甜高粱株高在第二和第三次取樣間處于拔節中期,其增幅最大,不同施氮水平的平均增幅為49.6%;第三次取樣后增幅趨緩,不同施氮水平下株高的平均增幅為10.8%。2010年試驗第三次取樣后,甜高粱株高依然有較高的增幅,平均增幅為17.0%,這可能是由于2010年平均氣溫高于2009年的緣故。兩年試驗中,不同施氮處理間甜高粱莖粗差異不大。第一次取樣到第二次取樣間,莖粗的增幅較大,而在第二次取樣后,增幅變緩,這一點在2010年的試驗中表現尤為明顯。

各取樣時間段內綠葉數的變化趨勢較為一致,不同處理間的差異大于株高和莖粗,呈現出集中分布的3組(圖1,圖2),這一現象在2010年試驗中表現更為明顯。在300和400 kg N/hm2處理下不同取樣點的平均綠葉數分別為18.0和17.8個,明顯高于其他處理水平;100,200和500 kg N/hm2處理,相互之間差異不大,不同取樣時間的平均綠葉數分別為17.3,17.5和17.1個,明顯高于不施氮處理(16.1個)。第一次取樣時,2009年的綠葉數要低于2010年,可能由于2009年播種時土壤溫度較低,導致苗期的生長較慢,進入拔節期時,表現不如2010年,但在最后一次取樣時,綠葉數與2010年無差異。

隨著生育期的推進,甜高粱植株的莖葉比逐漸增大(圖1,圖2),500 kg N/hm2處理下莖葉比的增速逐漸高于其他施氮處理,這一趨勢在2010年的試驗中表現更為明顯。隨著莖葉比不斷增大,植株鮮干比呈現逐漸減小的趨勢。2009年最后一次采樣時,500 kg N/hm2處理下的鮮干比(3.26);2010年不施氮處理下的鮮干比最大(3.36)。

葉面積指數隨著生育期的推進呈現出階段性增長的趨勢(圖1,圖2),第一次和第二次取樣間,第三次和第四次取樣間,葉面積指數的增幅較大,2009年各施氮水平的平均增幅分別為21.0%和23.7%,2010年的對應增幅分別為23.7%和20.7%。在第二次和第三次取樣間,葉面積指數的增幅較為平緩,兩年試驗各施氮水平的平均增幅分別為5.90%(2009年)和8.85%(2010年),不施氮處理和500 kg N/hm2處理下的葉面積指數,在這一階段幾乎沒有增加。最后一次取樣時,葉面積指數為分別為5.37(2009年)和5.12(2010年)。

圖1 施氮水平對不同生長時期甜高粱農藝性狀的影響(2009年)Fig.1 Effects of nitrogen fertilizer on agronomic traits of sweet sorghum at different growth stage (2009)

2.2甜高粱干物質產量對施氮水平響應的模擬分析

為了分析不同年份間施氮量對甜高粱干物質產量的影響,分別對不同施氮處理下的干物質產量測定值進行了數學模擬,建立了甜高粱干物質產量和施氮水平之間的回歸模型,并利用回歸模型計算了獲得最大干物質產量所需的最優施氮量。圖3表明,施氮量對甜高粱干物質產量的影響呈拋物線變化,隨施氮量的增加呈現先增加后降低的趨勢。通過回歸模型y=Ax2+Bx+C對施氮量的干物質產量進行擬合,決定系數分別為R2=0.77(2009)和R2=0.83(2010),計算出獲得最大干物質產量的施氮量分別為296.9 kg N/hm2(2009年)和323.0 kg N/hm2(2010年)。從擬合曲線可以看出,當施氮量為200~400 kg N/hm2時,干物質產量處于較高水平。2009年試驗結果顯示,當施氮量為200,300和400 kg N/hm2時,對應的干物質產量分別為28640.3,31504.4和29912.6 kg/hm2;2010年對應處理的干物質產量分別為27556.6,31465.1和29893.0 kg/hm2。擬合結果和試驗結果均表明,當施氮量為200~400 kg N/hm2時,大力士甜高粱可獲得較高的干物質產量。

圖2 施氮水平對不同生長時期甜高粱農藝性狀的影響(2010年)Fig.2 Effects of nitrogen fertilizer on agronomic traits of sweet sorghum at different growth stage (2010)

2.3甜高粱干物質產量與主要農藝性狀的相關性、多

圖3 施氮水平對干物質產量的回歸曲線Fig.3 Regression curve of dry matter yield under different nitrogen application rate

元回歸及通徑分析

表2表明,干物質產量(y)與株高(x1)、莖粗(x2)和葉面積指數(x6)呈極顯著正相關關系(P<0.01),相關系數分別為r=0.921,r=0.758和r=0.865;干物質產量(y)與綠葉數(x3)呈顯著正相關關系(P<0.05,r=0.696);株高(x1)與莖粗(x2),綠葉數(x3)和葉面積指數(x6)呈極顯著正相關關系,相關系數分別為r=0.817,r=0.863和r=0.769;莖粗(x2)與綠葉數(x3)呈極顯著正相關關系(r=0.904),與葉面積指數呈顯著正相關關系(r=0.632);干物質產量(y)與各農藝性狀的相關系數大小順序依次為:株高(x1)>葉面積指數(x6)>莖粗(x2)>綠葉數(x3)>鮮干比(x5)>莖葉比(x4),其中鮮干比(x5)和莖葉比(x4)與干物質產量呈負相關,相關系數分別為r=-0.442和r=-0.097。對干物質產量(y)進行正態性檢驗,其顯著水平P=0.811>0.05(表3),服從正態分布,可進行多元回歸分析。以株高(x1)、莖粗(x2)、綠葉數(x3)、莖葉比(x4)、鮮干比(x5)和葉面積指數(x6)為自變量,通過逐步回歸分析,建立上述參數與干物質產量(y)的最優多元回歸方程:y=-43.316+17.744x1+4.086x6(R2=0.908,F=44.572**)。 方程篩選了2個對干物質產量(y)有顯著效應的自變量,分別為株高(x1)和葉面積指數(x6)。這表明,當2個自變量中的一個取值固定在某一水平時,株高(x1)每增加1 m,干物質產量(y)增加17.744 kg;葉面積指數(x6)每增加1%,干物質產量(y)增加40.86 g。決定系數R2=0.908,表明株高(x1)和葉面積指數(x6)對干物質產量(y)的總影響達到90.8%以上。

表2 甜高粱干物質產量與農藝性狀相關性分析Table 2 Correlation analysis of dry matter yield with agronomic trait in sweet sorghum

注:*顯著水平為P<0.05,**顯著水平為P<0.01,x1,株高;x2,莖粗;x3,綠葉數;x4,莖葉比;x5,鮮干比;x6,葉面積指數;y,干物質產量。下同。Note: *significant at the 0.05 level of probability,**significant at the 0.01 level of probability,x1, plant height;x2, stem diameter;x3, green leaf number;x4, stem/leaf ratio;x5, fresh/dry matter ratio;x6, leaf area index;y, dry matter yield. The same below.

表3 正態性檢驗Table 3 Normality test

表4 干物質產量與主要農藝性狀的通徑分析Table 4 Path analysis of dry matter yield and main agronomic traits

在相關分析和多元回歸的基礎上,進一步進行通徑分析,結果表明(表4),各參數對干物質產量(y)的直接通徑系數大小依次為:株高(x1)>葉面積指數(x6)>綠葉數(x3)>莖葉比(x4)>鮮干比(x5)>莖粗(x2)。株高(x1)的直接通徑系數最大(P=0.638)。株高(x1)對干物質產量(y)的間接作用受綠葉數(x3)的負效應影響,但影響較小(P=-0.071),株高通過葉面積指數(x6)對干物質產量(y)產生的間接作用較大(P=0.301),說明株高并不是影響干物質產量(y)的主導因素,而是與葉面積指數(x6)共同對干物質產量(y)產生影響。

葉面積指數(x6)的直接通徑系數為P=0.391。葉面積指數(x6)通過株高(x1)對干物質產量的正向間接作用較高,間接通徑系數為0.491;通過綠葉數(x3)和莖葉比(x4)對干物質產量(y)的負向間接作用相對較小,間接通徑系數分別為-0.036和-0.014。這一結果表明,葉面積指數(x6)對干物質產量(y)所產生的間接正效應可以有效平衡直接負效應。

莖粗(x2)對干物質產量(y)的直接通徑系數較小(P=0.010),通過株高(x1)和葉面積指數(x6)對干物質產量(y)有最大的間接作用,間接通徑系數分別為0.556和0.247。綜合來看,株高(x1)和葉面積指數(x6)是影響甜高粱干物質產量(y)的最重要因素,在栽培中應采取措施增加甜高粱植株株高,并提高其葉面積指數。

2.4不同施氮水平對干物質產量與主要農藝性狀相互關系的影響

不同的施氮水平對各性狀影響不同,進而影響到最終干物質產量的高低,為了進一步了解不同施氮水平下,各農藝性狀之間及其與干物質產量的相互關系,試驗將施氮水平分為低氮(0,100 kg N/hm2)、中氮(200,300 kg N/hm2)和高氮(400,500 kg N/hm2)3組,分別進行多元回歸和通徑分析。

2.4.1低氮水平下干物質產量與主要農藝性狀的多元回歸及通徑分析施氮量為0和100 kg N/hm2水平下,對干物質產量和主要農藝性狀進行相關性分析。結果表明(表5),干物質產量(y)與莖葉比(x4)呈極顯著正相關關系(P<0.01,r=0.740),與鮮干比(x5)呈極顯著負相關關系(r=-0.715);干物質產量(y)與株高(x1)和葉面積指數(x6)呈顯著(P<0.05)正相關關系,相關系數分別為r=0.614和r=0.589;株高(x1)與綠葉數(x3)和葉面積指數(x6)呈極顯著正相關關系,相關系數分別為r=0.737和r=0.802;莖葉比(x4)與鮮干比(x5)呈極顯著負相關關系(r=-0.776)。低氮水平下,干物質產量(y)與各農藝性狀的相關系數大小順序依次為:莖葉比(x4)>鮮干比(x5)>株高(x1)>葉面積指數(x6)>綠葉數(x3)>莖粗(x2),與整體的相關性不同,低氮水平下,莖葉比(x4)和鮮干比(x5)與干物質產量的相關性要大于株高(x1)和葉面積指數(x6)。

表5 低氮水平下甜高粱干物質產量與農藝性狀相關性分析(0,100 kg N/hm2)Table 5 Correlation analysis of dry matter yield with agronomic trait in sweet sorghum under lower nitrogen application rate (0, 100 kg N/ha)

在低氮水平下(0和100 kg N/hm2),對干物質產量(y)進行正態性檢驗,其顯著水平P=0.561>0.05(表6),服從正態分布,可進行多元回歸分析。以株高(x1)、莖粗(x2)、綠葉數(x3)、莖葉比(x4)、鮮干比(x5)和葉面積指數(x6)為自變量,通過逐步回歸分析法,建立各參數與干物質產量(y)的最優多元回歸方程:y=-12.039+5.966x1+3.295x2+4.253x4(R2=0.912,F=27.475**)。方程篩選了3個對干物質產量(y)有顯著效應的自變量,分別為株高(x1)、莖粗(x2)和莖葉比(x4)。這表明,當3個自變量中的一個取值固定在某一水平時,株高(x1)每增加1 m,干物質產量(y)增加5.966 kg;莖粗(x2)每增加1 cm,干物質產量(y)增加3.295 kg;莖葉比(x4)每增加1%,干物質產量(y)增加42.53 g,決定系數R2=0.912,所以株高(x1)、莖粗(x2)和莖葉比(x4)對干物質產量(y)的總影響達到91.2%以上。

表6 低氮水平下正態性檢驗(0,100 kg N/hm2)Table 6 Normality test under lower nitrogen application rate (0, 100 kg N/ha)

低氮水平下對干物質產量和各農藝性狀進行通徑分析,討論各農藝性狀對干物質產量的直接作用和間接作用。表7表明,各參數對干物質產量(y)的直接通徑系數大小依次為:莖葉比(x4)>莖粗(x2)>株高(x1)>綠葉數(x3)>鮮干比(x5)>葉面積指數(x6)。

莖葉比(x4)的直接通徑系數最大(P=0.925)。莖葉比(x4)對干物質產量(y)的間接作用受莖粗(x2)、綠葉數(x3)和鮮干比(x5)的負效應影響,間接作用分別為-0.028,-0.063,-0.148,三者總的間接作用達-0.239,對莖葉比的影響較大,但因莖葉比(x4)的直接作用更大,所以在各參數中,是影響干物質產量(y)的最重要因素。

表7 低氮水平下干物質產量與主要農藝性狀的通徑分析(0,100 kg N/hm2)Table 7 Path analysis of dry matter yield and main agronomic traits under lower nitrogen application rate (0, 100 kg N/ha)

莖粗(x2)的直接通徑系數為0.293。莖粗(x2)通過株高(x1)、綠葉數(x3)和葉面積指數(x6)對干物質產量產生正向間接作用,間接通徑系數分別為0.032,0.043和0.01;通過莖葉比(x4)和鮮干比(x5)對干物質產量(y)產生負向間接作用,間接通徑系數分別為-0.090和-0.036,說明莖粗(x2)對干物質產量(y)所產生的間接正效應可以有效平衡直接負效應。

株高(x1)的直接通徑系數為0.246,對干物質產量(y)產生正向影響,但作用小于莖葉比(x4)和莖粗(x2)。株高(x1)僅通過鮮干比(x5)對干物質產量(y)產生負效應,且影響較小(P=-0.062),除此以外,通過其他參數均對干物質產量產生正向影響。綜合來看,低氮水平下,除株高(x1)外,莖葉比(x4)和莖粗(x2)是影響干物質產量(y)的重要因素。

2.4.2中氮水平下干物質產量與主要農藝性狀的多元回歸及通徑分析施氮為200和300 kg N/hm2水平下,對干物質產量和主要農藝性狀進行相關性分析。結果表明(表8),干物質產量(y)與株高(x1)和葉面積指數(x6)呈顯著正相關關系(P<0.05),相關系數分別為r=0.615和r=0.687。株高(x1)和莖粗(x2)呈顯著正相關關系(r=0.696);與鮮干比(x5)呈顯著負相關關系(r=-0.697)。莖粗(x2)與綠葉數(x3)呈極顯著正相關關系(P<0.01,r=0.778)。綠葉數(x3)與鮮干比(x5)呈顯著負相關關系(r=-0.643)。中氮水平下,干物質產量(y)與各農藝性狀的相關系數大小順序依次為:葉面積指數(x6)>株高(x1)>鮮干比(x5)>莖粗(x2)>綠葉數(x3)>莖葉比(x4),其中,鮮干比(x5)和莖葉比(x4)與干物質產量(y)呈負相關關系,相關系數分別為r=-0.477和r=-0.003。與低氮水平的相關性不同,中氮水平下,葉面積指數(x6)和株高(x1)與干物質產量的相關性要大于鮮干比(x5)和莖葉比(x4)。

中氮水平下(200和300 kg N/hm2),對干物質產量(y)進行正態性檢驗,其顯著水平P=0.239>0.05(表9),服從正態分布,可進行多元回歸分析,按低氮中的方法通過逐步回歸分析法,建立各參數與干物質產量(y)的最優多元回歸方程:y=-38.443+27.889x1-24.469x2+7.308x6(R2=0.839,F=13.847**)。方程篩選了3個對干物質產量(y)有顯著效應的自變量,分別為株高(x1)、莖粗(x2)和葉面積指數(x6)。這表明,當3個自變量中的一個取值固定在某一水平時,株高(x1)每增加1 m,干物質產量(y)增加27.889 kg;莖粗(x2)每增加1 cm,干物質產量(y)減少24.469 kg;葉面積指數(x6)每增加1%,干物質產量(y)增加73.08 g,決定系數R2=0.839,所以株高(x1)、莖粗(x2)和葉面積指數(x6)對干物質產量(y)的總影響達到83.9%以上。

中氮水平下對干物質產量和各農藝性狀進行通徑分析,討論各農藝性狀對干物質產量的直接作用和間接作用。表10表明,各參數對干物質產量(y)的直接通徑系數大小依次為:株高(x1)>葉面積指數(x6)>莖粗(x2)>莖葉比(x4)>鮮干比(x5)>綠葉數(x3)。

株高(x1)的直接通徑系數最大(P=0.861)。株高(x1)對干物質產量(y)的間接作用受莖粗(x2)、綠葉數(x3)和莖葉比(x4)的負效應影響,間接作用分別為-0.470,-0.013,-0.028,三者總的間接作用達-0.511,株高(x1)通過葉面積指數(x6)和鮮干比(x5)對干物質產量(y)產生正向間接作用,間接系數分別為0.048和0.217,兩者總的間接作用小于負向的間接作用,株高(x1)主要通過直接作用對干物質產量(y)產生影響,再次表明,中氮水平下,株高(x1)是影響干物質產量(y)的重要性狀之一。

表8 中氮水平下甜高粱干物質產量與農藝性狀相關性分析(200,300 kg N/hm2)Table 8 Correlation analysis of dry matter yield with agronomic trait in sweet sorghum under moderate nitrogen application rate (200, 300 kg N/ha)

葉面積指數(x6)的直接通徑系數為0.797。葉面積指數(x6)通過株高(x1)、莖葉比(x4)和鮮干比(x5)對干物質產量產生正向間接作用,間接通徑系數分別為0.234,0.023和0.015;通過莖粗(x2)和綠葉數(x3)對干物質產量(y)產生負向間接作用,間接通徑系數分別為-0.371和-0.010,說明葉面積指數(x6)對干物質產量(y)所產生的間接正效應可以有效平衡直接負效應。

表9 中氮水平下正態性檢驗(200,300 kg N/hm2)Table 9 Normality test under moderate nitrogen application rate (200, 300 kg N/ha)

表10 中氮水平下干物質產量與主要農藝性狀的通徑分析(200,300 kg N/hm2)Table 10 Path analysis of dry matter yield and main agronomic traits under moderate nitrogen application rate (200, 300 kg N/ha)

莖粗(x2)的直接通徑系數為-0.675,對干物質產量(y)產生負向影響,莖粗(x2)通過株高(x1)、鮮干比(x5)和葉面積指數(x6)對干物質產量(y)產生較大的正向影響,間接系數分別為0.599,0.035和0.438,這說明,中氮水平下莖粗(x2)主要通過株高(x1)和葉面積指數(x6)的間接作用影響干物質產量(y)。綜合來看,中氮水平下,株高(x1)和葉面積指數(x6)是影響干物質產量(y)的主要性狀,相比于低氮水平,莖葉比(x4)的影響下降,莖粗(x2)則對干物質產量(y)有負向影響。說明隨著施氮水平的提高,葉面積指數(x6)對干物質產量的重要性提高。

2.4.3高氮水平下干物質產量與主要農藝性狀的多元回歸及通徑分析施氮為400和500 kg N/hm2水平下,對干物質產量和主要農藝性狀進行相關性分析。結果表明(表11),干物質產量(y)與株高(x1)呈極顯著正相關關系(P<0.01,r=0.932),與鮮干比(x5)呈極顯著負相關關系(r=-0.884),與葉面積指數(x6)呈顯著正相關關系(P<0.05,r=0.650)。株高(x1)與鮮干比(x5)呈極顯著負相關關系(r=-0.792)。鮮干比(x5)與葉面積指數(x6)呈極顯著負相關關系(r=-0.752)。高氮水平下,干物質產量(y)與各農藝性狀的相關系數大小順序依次為:株高(x1)>鮮干比(x5)>葉面積指數(x6)>綠葉數(x3)>莖葉比(x4)>莖粗(x2),其中,莖粗(x2)與干物質產量(y)呈負相關關系(r=-0.166)。與中氮水平的相關性不同,高氮水平下,鮮干比(x5)與干物質產量的相關性要大于葉面積指數(x6)。

高氮水平下(400和500 kg N/hm2),對干物質產量(y)進行正態性檢驗,其顯著水平P=0.275>0.05(表12),服從正態分布,可進行多元回歸分析,按低氮處理的方法,通過逐步回歸分析,建立各參數與干物質產量(y)的最優多元回歸方程:y=8.653+12.986x1-5.592x5(R2=0.925,F=55.884**)。方程篩選了2個對干物質產量(y)有顯著效應的自變量,分別為株高(x1)和鮮干比(x5)。這表明,當2個自變量中的一個取值固定在某一水平時,株高(x1)每增加1 m,干物質產量(y)增加12.986 kg;鮮干比(x5)每增加1%,干物質產量(y)減少55.92 g;決定系數R2=0.925,所以株高(x1)和鮮干比(x5)對干物質產量(y)的總影響達到92.5%以上。

表11 高氮水平下甜高粱干物質產量與農藝性狀相關性分析(400,500 kg N/hm2)Table 11 Correlation analysis of dry matter yield with agronomic trait in sweet sorghum under higher nitrogen application rate (400, 500 kg N/ha)

高氮水平下對干物質產量和各農藝性狀進行通徑分析,討論各農藝性狀對干物質產量的直接作用和間接作用。表13表明,各參數對干物質產量(y)的直接通徑系數大小依次為:株高(x1)>鮮干比(x5)>葉面積指數(x6)>綠葉數(x3)>莖葉比(x4)>莖粗(x2)。

株高(x1)的直接通徑系數最大(P=0.562)。株高(x1)對干物質產量(y)的間接作用受綠葉數(x3)和葉面積指數(x6)的負效應影響,間接作用分別為-0.094和-0.022,兩者的影響相對較小,株高(x1)通過其他參數對干物質產量(y)產生正向間接作用,影響最大的是鮮干比(x5),間接系數為0.408。說明在高氮情況下株高(x1)和鮮干比(x5)共同對干物質產量(y)產生影響,株高(x1)不再具有決定性的作用。

鮮干比(x5)的直接通徑系數為-0.515。鮮干比(x5)通過株高(x1)和莖葉比(x4)對干物質產量產生負向間接作用,間接通徑系數分別為-0.445和-0.032;通過莖粗(x2)和綠葉數(x3)和葉面積指數(x6)對干物質產量(y)產生正向間接作用,間接通徑系數分別為0.015,0.063和0.030,說明鮮干比(x5)株高(x1)的間接作用和自身的直接作用共同影響對干物質產量(y)所產生負效應。

葉面積指數(x6)的直接通徑系數為-0.040,對干物質產量(y)產生負向影響,葉面積指數(x6)通過株高(x1)、莖粗(x2)和鮮干比(x5)對干物質產量(y)產生較大的正向影響,間接系數分別為0.311,0.019和0.387;通過綠葉數(x3)和莖葉比(x4)對干物質產量(y)產生的間接負向作用較小,間接系數分別為-0.020和-0.007。這說明,高氮水平下葉面積指數(x6)對干物質產量(y)的直接作用小于通過株高(x1)、莖粗(x2)和鮮干比(x5)對干物質產量(y)產生的間接作用。綜合來看,在高氮水平下,株高依然是影響干物質產量的重要因素,而葉面積指數對干物質產量的影響低于中氮水平,并開始有負向的直接影響。

表13 高氮水平下干物質產量與主要農藝性狀的通徑分析(400,500 kg N/hm2)Table 13 Path analysis of dry matter yield and main agronomic traits under higher nitrogen application rate (400, 500 kg N/ha)

3 討論

兩年試驗中,甜高粱各農藝性狀在施氮量為300和400 kg N/hm2綜合表現較好,各性狀隨施氮量的變化趨勢基本一致。通過對施氮量和干物質產量進行二次回歸擬合得出,施氮量在200~400 kg N/hm2時,可獲得較高的干物質產量。500 kg N/hm2處理下的干物質產量分別為27132.6 kg/hm2(2009年)和26928.8 kg/hm2(2010年),與200 kg N/hm2處理下的干物質產量相近,但投入產出比低。

飼用作物的生物產量受多個農藝性狀的影響。不同品種和栽培措施使得農藝性狀對生物產量構成的重要性不盡相同。譚禾平等[25]通過綜合分析和最優線性回歸選擇后發現,株高、穗粗、穗長、禿尖長、百粒重和出籽率是影響玉米產量的最有效因子。王亞玲等[26]通過分析苜蓿(Medicagosativa)產量構成因素發現,對產量貢獻率最大的是生長高度。楊偉光等[27]研究表明,分蘗數對甜高粱生物產量的影響最大,其次才是株高。向清華等[28]的研究表明,莖重和主莖直徑對紫花苜蓿干物質產量有較大作用。馮西博等[29]研究了青稞(Hordeumvulgare)農藝性狀在不同施氮水平下對產量的貢獻后發現,不論施氮量是否合理,每穗結實小穗數和千粒重均是影響青稞產量的關鍵因子。因此,明確農藝性狀與生物產量的關系,可以更好地為高產品種的選育提供參考。

本試驗通過對低氮、中氮和高氮3種施氮水平下的多元回歸和通徑分析發現,不論施氮水平是否合理,株高始終是影響甜高粱干物質產量的關鍵因素。低氮水平,干物質產量的主要影響因素為株高、莖粗和莖葉比;中氮水平,干物質產量的主要影響因素為株高、莖粗和葉面積指數;高氮水平,干物質產量的主要影響因素為株高、鮮干比。莖粗在低氮水平下對干物質產量為正向直接作用,而在中氮水平下對干物質產量的影響主要以負向直接作用(P=-0.675)和通過株高(P=0.599)和葉面積指數(P=0.438)的正向間接作用為主,且莖粗與干物質產量呈較弱的正相關關系(r=0.378),說明中氮條件下,莖粗并不能直接提高干物質產量,而是莖粗的增加會削弱其他性狀的表現,進而影響干物質產量的增加。隨著施氮量的提高,葉面積指數對干物質產量的影響增大,但隨著施氮量的進一步提高,高氮水平下,葉面積指數對干物質產量的直接影響為負效應,且與干物質產量的相關性低于鮮干比,說明高氮水平下,葉面積指數成了干物質產量的限制因子。這可能是由于高水平氮肥促進了植株葉面積的生長,使葉面積指數過大,葉片之間相互遮陰,造成群體間通風不暢,中下部葉片光照不足,導致整個植株光合效率的下降,最終影響了光合產物的合成。因此在甜高粱栽培時,施肥量不宜過高;土壤貧瘠需要較多肥料時,應仔細考慮種植密度,合理安排群體結構;施氮量較低時,莖葉比是影響干物質產量的主要因素,延遲收獲可以提高莖葉比,有利于干物質產量的提高,但過高的莖葉比會降低飼草營養價值[30]。因此施氮量較低時,甜高粱收獲時間的選擇是平衡飼草產量和品質的關鍵。

4 結論

兩年試驗結果表明,200~400 kg N/hm2是較適宜的施氮量,可獲得較大的干物質產量。多元回歸和通徑分析表明,株高和葉面積指數是影響干物質產量的最重要因素,但在高氮(≥400 kg N/hm2)水平下,葉面積指數對干物質產量有負向直接作用。低氮(≤100 kg N/hm2)水平下,莖葉比對干物質產量的正向影響增大,在此施氮水平下如果將甜高粱作為飼草利用,栽培時應特別注意選擇適宜的收獲時間,以便更好地調整干物質產量和營養價值間的平衡關系。

References:

[1]Elia A, Santamaria P, Serio F. Nitrogen nutrition, yield and quality of spinach. Journal of the Science of Food and Agriculture, 1998, 76(3): 341-346.

[2]Ferguson R, Hergert G, Schepers J,etal. Site-specific nitrogen management of irrigated maize. Soil Science Society of America Journal, 2002, 66(2): 544-553.

[3]Samonte S, Wilson L, Medley J,etal. Nitrogen utilization efficiency: relationships with grain yield, grain protein and yield-related traits in rice. Agronomy Journal, 2006, 98(1): 168-176.

[4]Stewart W, Dibb D, Johnston A,etal. The contribution of commercial fertilizer nutrients to food production. Agronomy Journal, 2005, 97(1): 1-6.

[5]Reddy B, Sanjana Reddy P, Bidinger F,etal. Crop management factors influencing yield and quality of crop residues. Field Crops Research, 2003, 84(1): 57-77.

[6]Marsalis M, Angadi S, Contreras-Govea F. Dry matter yield and nutritive value of corn, forage sorghum, and BMR forage sorghum at different plant populations and nitrogen rates. Field Crops Research, 2010, 116(1): 52-57.

[7]Almodares A, Taheri R, Hadi M,etal. The effect of nitrogen and potassium fertilizers on the growth parameters. Pakistan Journal of Biological Sciences, 2006, 9(12): 2350-2353.

[8]Ding C L, Gu H R, Bai S J,etal. Influence of different rate of fertilizer application and plant densities on the yield of pennisetum glaucum. Grassland of China, 1999, (5): 12-14.

[9]Wei J S, Liu G D, Cai B Y. A study on effect of cutting period and fertilizing to yield and quality of panicum maximum. Acta Agrestia Sinica, 2002, 10(2):139-143.

[10]Cao W X, Liu H D, Li W,etal. Nitrogen fertilizer effects on a 15 year old mixed species alpine pasture after two years. Acta Prataculturae Sinica, 2015, 24(9): 130-137.

[11]Vasilakoglou I, Dhima K, Karagiannidis N,etal. Sweet sorghum productivity for biofuels under increased soil salinity and reduced irrigation. Field Crops Research, 2011, 120(1): 38-46.

[12]Ibrahim Y M. Comparative analysis of forage sorghum for best cutting time. Annals of Arid Zone, 1999, 38(1): 75-78.

[13]Mccaughey W P, Therrien M C, Mabon R. Forage sorghum in southern Manitoba. Canadian Journal of Plant Science, 1996, 76(1): 123-125.

[14]Rooney W L, Blumenthal J, Bean B,etal. Designing sorghum as a dedicated bioenergy feedstock. Biofuels, Bioproducts and Biorefining, 2007, 1(2): 147-157.

[15]Chang J H, Xia X Y, Zhang L,etal. Analysis of the resistance gene to the sorghum aphid,Melanaphissacchari, with SSR marker inSorghumbicolor. Acta Prataculturae Sinica, 2006, 15(2): 113-118.

[16]Sun S J, Wang Y, Zheng G C,etal. The feasibility research of breeding and application of Sorghum-Sudan grass hybrids. Acta Prataculturae Sinica, 1999, 8(3): 39-39.

[17]Burke J, Franks C, Burow G,etal. Selection system for the stay-green drought tolerance trait in sorghum germplasm. Agronomy Journal, 2010, 102(4): 1118-1122.

[18]Li C X, Feng H S. A study on the adaptability of sweet sorghum planted in different altitudinal areas of the Qinghai plateau. Acta Prataculturae Sinica, 2013, 22(3): 51-59.

[19]Marsalis M A, Angadi S, Contreras-Govea F E,etal. Harvest timing and byproduct addition effects on corn and forage sorghum silage grown under water stress. Research Bulletin BL-799[C]. Las Cruces, NM, USA: Agricultural Experiment Station, New Mexico State University, 2009.

[20]Zhang S J, Amerjan O, Xue X Z,etal. Quality analysis on different sweet sorghum silages in Southern Xinjiang compared with a corn silage. Acta Prataculturae Sinica, 2014, 23(3): 232-240.

[21]Howell T, Evett S, Tolk J,etal. Evapotranspiration of corn and forage sorghum for silage. Proceedings of World Environmental and Water Resources Congress[C]. US: Hawaii, 2008: 1-14.

[22]Afzal M, Ahmad A U H, Zamir S I. Performance of multicut forage sorghum under various sowing methods and nitrogen application rates. The Journal of Animal and Plant Sciences, 2013, 23(1): 232-239.

[23]Zegada-Lizarazu W, Monti A. Are we ready to cultivate sweet sorghum as a bioenergy feedstock? A review on field management practices. Biomass and Bioenergy, 2012, (40): 1-12.

[24]Du J J, Chen Z W. Method of path analysis by using SPSS linear regression program. Bulletin of Biology, 2010, (2): 4-6.

[25]Tan H P, Wang G Y, Hu X N,etal. Multiple regression and path analysis of effective factors affecting maize yield. Acta Agriculturae Zhejiangensis, 2006, 18(4): 238-240.

[26]Wang Y L, Li X F, Shi S L,etal. Analysis and evaluation on the production performance components of different alfalfa. Chinese Journal of Grassland, 2007, 29(5): 8-15.

[27]Yang W G, Yang F, Gao C F,etal. Studies on the genetic correlation of important traits of sweet sorghum. Journal of Jilin Agricultural University, 1995, 17(1): 29-31.

[28]Xiang Q H, Deng R, Zhang D H,etal. Application of path analysis on economical character of alfalfa. Guizhou Journal of Animal Husbandry and Veterinary Medicine, 2010, 34(5): 44-45.

[29]Fu X B, Wang J J, Du Jie-ciren,etal. Analysis of the contribution of the agronomic characters of Tibet spring naked barley under different N fertilizers and seeding dates to yield. Barley and Cereal Sciences, 2009, (2): 1-5.

[30]Zhang D Q, Zhou G Y, Wen D Z,etal. Effect of cutting frequency on yield and quality of forage grasses. Journal of Tropical and Subtropical Botany, 2000, (S1): 43-51.

[8]丁成龍, 顧洪如, 白淑娟, 等. 不同施肥量, 密度對美洲狼尾草產量的影響. 中國草地, 1999, (5): 12-14.

[9]韋家少, 劉國道, 蔡碧云. 刈割周期, 施肥量對堅尼草產草量和質量的影響. 草地學報, 2002, 10(2): 139-143.

[10]曹文俠, 劉皓棟, 李文, 等. 連續兩年施氮對15齡混作禾草草地的改良效果研究. 草業學報, 2015, 24(9): 130-137.

[15]常金華, 夏雪巖, 張麗, 等. 高梁抗蚜基因的遺傳分析和SSR標記定位. 草業學報, 2006, 15(2): 113-118.

[16]孫守鈞, 王云, 鄭根昌, 等. 高粱—蘇丹草雜交種的應用可行性研究及其選育. 草業學報, 1999, 8(3): 39-39.

[18]李春喜, 馮海生. 甜高粱在青海高原不同海拔生態區的適應性研究. 草業學報, 2013, 22(3): 51-59.

[20]張蘇江, 艾買爾江·吾斯曼, 薛興中, 等. 南疆玉米和不同糖分甜高梁的青貯品質分析. 草業學報, 2014, 23(3): 232-240.

[24]杜家菊, 陳志偉. 使用SPSS線性回歸實現通徑分析的方法. 生物學通報, 2010, (2): 4-6.

[25]譚禾平, 王桂躍, 胡賢女, 等. 影響玉米產量效應因子的多元回歸與通徑分析. 浙江農業學報, 2006, 18(4): 238-240.

[26]王亞玲, 李曉芳, 師尚禮, 等. 紫花苜蓿生產性能構成因子分析與評價. 中國草地學報, 2007, 29(5): 8-15.

[27]楊偉光, 楊福, 高春福, 等. 甜高粱主要農藝性狀相關性研究. 吉林農業大學學報, 1995, 17(1): 29-31.

[28]向清華, 鄧蓉, 張定紅, 等. 通徑分析在紫花苜蓿農藝性狀分析中的應用. 貴州畜牧獸醫, 2010, 34(5): 44-45.

[29]馮西博, 王金菊, 趙慧芬, 等. 西藏春青稞農藝性狀在不同施氮量和播期水平下對產量貢獻的分析. 大麥與谷類科學, 2009, (2): 1-5.

[30]張德強, 周國逸, 溫達志, 等. 刈割時間和次數對牧草產量和品質的影響. 熱帶亞熱帶植物學報, 2000, (S1): 43-51.

*Effects of nitrogen on agronomic traits and dry matter yield and their relationship in sweet sorghum

QU Hui, CHENG Liang, CHEN Jun-Feng, LU Xiao-Yan, SHEN Yi-Xin*

CollegeofAnimalScienceandTechnology,NanjingAgriculturalUniversity,Nanjing210095,China

This study was conducted to determine the relationship between dry matter (DM) yield and agronomic traits of sweet sorghum (Sorghumbicolorcv. ‘Hunnigreen’) using correlation and path coefficient analysis under different rates of nitrogen (0, 100, 200, 300, 400, 500 kg N/ha). The experiment was set up in eastern China during 2009-2010 and 7 agronomic traits of sorghum were evaluated. Results showed that DM yield was positively correlated (P<0.01) with plant height (r=0.921**), leaf area index (LAI,r=0.865**) and stem diameter (SD,r=0.758**). Plant height and LAI produced the strongest path coefficient (P=0.638,P=0.391) for DM yield. Plant height was the best indicator of DM yield using stepwise regression analysis. The contribution of stem∶leaf ratio (S∶L ratio) and LAI to DM yield differed under different rates of N fertilizer. S∶L ratio and LAI had positive path coefficient (P=0.925,P=0.162) for DM yield under lower N fertilizer rate (0, 100 kg N/ha). LAI had a positive path coefficient (P=0.797) for DM yield under intermediate N fertilizer rates (200, 300 kg N/ha), but had a negative path coefficient for DM yield under higher N fertilizer rates (400, 500 kg N/ha). The Effect of S∶L ratio on DM yield decreased with increasing N rate. DM yield was highest at N rates of 200-400 kg N/ha. Plant population was important for higher DM yield at higher rates of N fertilizer. Under lower N fertilizer rates, harvest date was a key issue in terms of the trade-off between DM yield and nutritive value of sweet sorghum.

sweet sorghum; nitrogen fertilizer rate; dry matter yield; agronomic traits; path analysis

10.11686/cyxb2015400

http://cyxb.lzu.edu.cn

2015-09-01;改回日期:2015-12-10基金項目:國家科技支撐計劃南方優質飼草高效生產加工利用關鍵技術研究與集成示范項目(2011BAD17B03)資助。

渠暉(1984-),男,內蒙古烏蘭察布人,在讀博士。E-mail: qhyulin@hotmail.com

Corresponding author. E-mail: yxshen@njau.edu.cn

渠暉,程亮,陳俊峰,陸曉燕,沈益新. 施氮水平對甜高粱主要農藝性狀及其與干物質產量相關關系的影響. 草業學報, 2016, 25(6): 13-25.

QU Hui, CHENG Liang, CHEN Jun-Feng, LU Xiao-Yan, SHEN Yi-Xin. Effects of nitrogen on agronomic traits and dry matter yield and their relationship in sweet sorghum. Acta Prataculturae Sinica, 2016, 25(6): 13-25.

猜你喜歡
產量水平
張水平作品
2022年11月份我國鋅產量同比增長2.9% 鉛產量同比增長5.6%
提高玉米產量 膜下滴灌有效
今日農業(2021年14期)2021-11-25 23:57:29
世界致密油產量發展趨勢
海水稻產量測評平均產量逐年遞增
今日農業(2020年20期)2020-11-26 06:09:10
作家葛水平
火花(2019年12期)2019-12-26 01:00:28
2018年我國主要水果產量按省(區、市)分布
加強上下聯動 提升人大履職水平
人大建設(2019年12期)2019-05-21 02:55:32
2018上半年我國PVC產量數據
聚氯乙烯(2018年9期)2018-02-18 01:11:34
老虎獻臀
主站蜘蛛池模板: 国产精品福利导航| 久久久精品无码一区二区三区| 在线亚洲小视频| 波多野结衣一二三| 欧美国产日产一区二区| 人妻少妇乱子伦精品无码专区毛片| 欧美日韩国产精品综合 | 一级成人欧美一区在线观看 | 欧美精品在线视频观看| 天堂成人在线| 久久精品66| 一级看片免费视频| 99爱在线| 91小视频在线观看免费版高清| 免费国产在线精品一区 | 亚洲无码免费黄色网址| 国产尤物视频在线| 免费观看欧美性一级| 欧美综合区自拍亚洲综合天堂 | 不卡视频国产| 亚洲人成网址| 国产美女无遮挡免费视频| 国产激情影院| 在线观看av永久| 亚洲国产中文精品va在线播放| 国产欧美在线观看精品一区污| 中文成人无码国产亚洲| 狠狠做深爱婷婷久久一区| 欧美精品二区| 亚洲综合精品香蕉久久网| 国产欧美视频在线观看| 精品综合久久久久久97超人| 国产精选自拍| 91福利免费| 依依成人精品无v国产| 亚洲色图欧美在线| 91久久夜色精品国产网站| 2020久久国产综合精品swag| 久热中文字幕在线观看| 成人免费黄色小视频| 国产91视频免费| 亚洲日韩在线满18点击进入| 五月天福利视频| 老司机aⅴ在线精品导航| 国产91丝袜在线播放动漫| 亚洲日本www| 99久久免费精品特色大片| 国产成人在线无码免费视频| 精品无码国产自产野外拍在线| 欧美日韩午夜| 福利在线免费视频| 欧美精品H在线播放| 婷婷色中文| 亚洲丝袜第一页| 四虎永久在线视频| 国产日韩欧美中文| 国产丝袜一区二区三区视频免下载| 天天视频在线91频| 91精选国产大片| 一级一毛片a级毛片| 国产日韩丝袜一二三区| 国产精鲁鲁网在线视频| 国产精品性| 亚洲综合狠狠| 最新日韩AV网址在线观看| 国产一区二区影院| 国产91av在线| 欧美日韩在线成人| 国产自在线播放| 欧美日韩久久综合| 极品国产在线| 91综合色区亚洲熟妇p| 婷婷色丁香综合激情| 就去吻亚洲精品国产欧美| 国产精品手机在线观看你懂的| 在线a视频免费观看| 久久中文字幕av不卡一区二区| 欧美在线精品怡红院| 日韩欧美在线观看| 91精品福利自产拍在线观看| 中文字幕在线视频免费| 欧美一区二区精品久久久|