?

保水劑與氮肥協同作用對綠竹抗性生理影響的綜合評價

2019-08-30 07:07張迎輝陳禮光凡莉莉榮俊冬鄭晶晶鄭郁善
西南林業大學學報 2019年5期
關鍵詞:綠竹保水劑沙地

張迎輝 陳禮光 凡莉莉 榮俊冬 鄭晶晶 鄭郁善

( 1. 福建農業職業技術學院,福建 福州 350003;2. 福建農林大學林學院,福建 福州 350012)

我國東南沿海長期以來以木麻黃(Casuarina equisetifolia)、濕地松(Pinus elliottii)等樹種構建防護林體系,存在著林分結構簡單、景觀格局單調、系統穩定性差等缺點[1]。篩選造林樹種,對豐富沿海沙地生物多樣性和提高防護林生態系統的穩定性具有重要的意義[2-3]。綠竹(Bambusa oldhamii)是優良筍用叢生竹種,具有生長快、成材早、發筍期長、產量高等特點,同時表現出良好的抗旱[4]、抗鹽[5]、抗風抗沙[6]能力,具有較高的生態價值和經濟價值。目前,綠竹作為防護林樹種已成功引種到沿海沙地種植,并表現出良好的生長特性[7-8]。研究表明綠竹能夠在低于鹽分質量分數為0.4%的土壤中生存良好[9];沿海沙地引種綠竹能夠改善土壤養分,提高土壤酶活性[10];綠竹與木麻黃混交能夠提高凋落物持水能力,減少地表徑流,具有涵養水源的重要意義[11]。

福建東南沿海土壤多為潮積或風積沙壤土,肥力低,保水性差,夏季出筍旺盛期時易發生短期高溫干旱脅迫,容易造成出筍率低和幼竹死亡,影響竹林生長。土壤干旱影響植物生長發育的各個階段的生理代謝過程[12-13]。保水劑是一種吸水能力特別強的高分子材料,可反復吸水、釋水,還可吸收肥料并緩慢釋放,增加肥效,廣泛用于農業、林業等生產實踐。馬丙堯等[14]認為保水劑—尿素凝膠方式能改善核桃(Juglans regia)苗根際土壤微生態環境,增強抵御干旱脅迫能力,對其生長達到一定的促進作用。宋雙雙等[15]認為在干旱、半干旱區樟子松(Pinus sylvestris)造林時,采用中濃度保水劑和微生物菌肥配施可以改良土壤,提高造林成活率。本研究以沿海沙地綠竹為研究對象,比較不同用量保水劑及氮肥混施后對綠竹葉片生理特性的影響,以期為保水劑在沿海沙地竹林培育中的應用提供理論基礎。

1 試驗地概況

試驗地位于福建省長樂大鶴國有防護林場,地處東經 119°40′43″,北緯 25°57′59″,屬于南亞熱帶海洋性季風氣候,6—9月為夏季高溫干旱,11月到翌年3月為冬春旱季。年降雨量約為1 200 mm,年平均風速為4.2 m/s,年平均氣溫為19.2 ℃,最高氣溫為35.6 ℃,最低氣溫為0 ℃。全年盛行東北風,可達280~300 d。土壤為濱海沙土,有機質含量為2.197 g/kg,水解氮含量7.12~7.88 mg/kg,有效磷含量3.1 mg/kg,速效鉀含量7.64 mg/kg[16],試驗期間試驗地氣溫和土壤含水率調查情況見圖1。天然植被稀少,林分結構簡單,林下常見零星植物有馬纓丹(Lantana camara)、茅莓(Rubus parvifolius)、苞薔薇(Rosa breateata)等。

圖 1 2014年試驗地平均氣溫和土壤含水率Fig. 1 Average temperature and soil moisture content of experimental field in 2014

2 材料與方法

2.1 試驗材料

試驗竹林為2010年4月在沿海沙地木麻黃純林下套種綠竹林,株行距為3 m×3 m,造林面積1.1 hm2。造成林時選擇生長健壯、無病蟲害的2年生母竹,竹稈劈成馬耳型向上切口收集雨水。種植穴規格為80 cm×80 cm×50 cm,每穴放黃心土50 kg,復合肥0.5 kg。

2.2 試驗設計

樣地設置為20 m×20 m,重復3次。每塊樣地10個樣叢,每個樣叢1組處理,共30個樣叢。CK無任何處理,即不施用保水劑和氮肥。保水劑設置3個水平,即A1(30 g/叢)、A2(60 g/叢)、A3(90 g/叢),氮肥設置3個水平,即B1(200 g/叢)、B2(400 g/叢)、B3(600 g/叢),共 9 個混施處理,標記為 A1B1、A1B2、A1B3、A2B1、A2B2、A2B3、A3B1、A3B2、A3B3。距竹蔸 30 cm,溝深 50 cm左右,將保水劑和細沙按1∶2質量比混勻進行環形溝施,再覆蓋沙土并澆水至飽和。于2014年4月施用保水劑,在試驗期間竹林不進行任何灌溉。于2014年7月、8月、9月中旬連續多日晴朗天氣下進行相關指標取樣、測定和分析。保水劑采用河南省農業科學院植物營養與資源環境研究所研制的營養型抗旱保水劑,主要成分為聚丙烯酰胺類物質、腐殖酸有機物質和稀土等。氮肥采用尿素。

2.3 生理指標測定

采樣時間為8:00—9:00,取植株中上部第3~8片完全展開葉為供試材料,用冰袋保鮮,帶回實驗室進行處理。用蒸餾水擦凈葉片表面污物,然后將葉片剪碎、混合均勻,用液氮冷凍后置于-80 ℃冰箱中保存,用于指標測定。葉片葉綠素含量用直接浸提法測定[17];可溶性蛋白含量用考馬斯亮藍法測定[18];可溶性糖含量用蒽酮比色法測定;脯氨酸含量用酸性茚三酮法[18]測定;丙二醛(MDA)含量用硫代巴比妥酸法測定[18];超氧化物歧化酶(SOD)活性用氮藍四唑光化還原法[18]測定;過氧化物酶(POD)活性用愈創木酚法[18]測定;電解質滲透率用電導法測定[18],每個處理重復3次。

2.4 數據處理

采用Excel 2010和SPSS 19.0進行數據統計分析和作圖。用LSD多重比較檢驗綠竹同一時期不同保水劑與氮肥處理下的差異顯著性,顯著性水平為0.05。采用主成分分析法[19]和隸屬函數法[20]對保水劑與氮肥下各項測定指標進行轉換和綜合評價分析。保水劑與氮肥施用效果呈正相關采用隸屬函數公式(公式(1))計算,呈負相關則用反隸屬函數公式(公式(2))計算。

式中:U(Xijk)為第i個處理第j個時期第k項指標的隸屬函數值,且U(Xijk)∈[0,1];Xijk為第i個處理第j個時期第k項指標測定值;Xkmax,Xkmin為所有處理指標中第k項指標的最大和最小測定值。

3 結果與分析

3.1 保水劑與氮肥對綠竹抗性生理的影響

3.1.1 保水劑與氮肥對綠竹葉綠素含量的影響

葉綠素是植物進行光合作用的主要色素,是反映植物在沙地逆境能否順利進行光合作用的重要生理指標。由圖2可知,在保水劑和氮肥處理下,除A2B1外,綠竹葉片中葉綠素含量均高于CK。其中,A3B2、A1B2和A3B1的葉綠素含量最高,分別比CK提高了89.67%(P<0.05)、37.86%和 55.18%(P<0.05)。

圖 2 保水劑與氮肥對葉綠素含量的影響Fig. 2 Effects of water-retaining agent and nitrogen fertilizers mixture on chlorophyll

3.1.2 保水劑與氮肥對綠竹滲透調節物質的影響

由圖3可知,在保水劑和氮肥處理下,綠竹葉片中可溶性糖含量均高于CK,可溶性蛋白和脯氨酸含量基本高于CK,但差異不顯著。其中A2B1、A1B1和A2B2的可溶性糖含量最高,分別比CK提高了58.78%、19.53%和44.06%;A2B2、A2B2和A3B1的可溶性蛋白含量最高,分別比CK提高了147.78%(P<0.05)、152.17%和148.61%;A1B2、A3B2和A3B2的脯氨酸含量最高,分別比CK提高了66.38%,112.95%和33.82%(P<0.05)。在相同保水劑處理下,可溶性糖和可溶性蛋白含量呈現先降低后升高變化趨勢,均在9月表現最高;脯氨酸含量呈現逐漸升高變化趨勢,9月表現最高。

圖 3 保水劑與氮肥對滲透調節物質的影響Fig. 3 Effects of water-retaining agent and nitrogen fertilizers mixture on osmotic adjustment substances

3.1.3 保水劑與氮肥對綠竹MDA含量和電解質滲透率的影響

由圖4可知,在保水劑和氮肥處理下,綠竹葉片中MDA含量和電解質滲透率均低于CK。A2B1、A1B2和A1B2的MDA含量最低,分別比CK降低了 64.15%、95.46%和 93.29%;A3B2處理下電解質滲透率最低,分別比CK降低了63.33%、25.51%和37.50%。在相同保水劑處理下,MDA含量和電解質滲透率呈現先升高后降低變化趨勢,均在8月份表現最高。

圖 4 保水劑與氮肥對MDA含量和電解質滲透率的影響Fig. 4 Effects of water-retaining agent and nitrogen fertilizers mixture on MDA content and electrolyte permeability

3.1.4 保水劑與氮肥對綠竹SOD活性和POD活性的影響

由圖5可知,在保水劑和氮肥處理下,綠竹葉片中SOD活性基本高于CK,POD活性均高于CK。A1B1、A3B3和A2B3的SOD活性最高,分別比 CK提高了 121.24%、728.29%和 63.78%;A3B2、A1B2和A3B1的POD活性最高,分別比CK提高了140.91%、488.76%(P<0.05)和189.40%。在相同保水劑處理下,SOD活性和POD活性呈現先降低后升高變化趨勢,均在9月份表現最高。

圖 5 保水劑與氮肥對SOD和POD活性的影響Fig. 5 Effects of water-retaining agent and nitrogen fertilizers mixture on SOD and POD activity

3.2 綠竹抗性生理指標的相關性分析

由表1可知,保水劑與氮肥處理下,綠竹葉片中8個抗性生理指標間存在一定相關關系。其中,葉綠素含量與可溶性蛋白含量和POD活性成極顯著正相關關系(P<0.01),葉綠素含量與脯氨酸含量呈顯著正相關關系(P<0.05),與電解質滲透率呈顯著負相關關系( P <0.05);可溶性糖含量與脯氨酸含量呈顯著正相關關系(P<0.05);可溶性蛋白含量與POD活性成極顯著正相關關系(P<0.01),與電解質滲透率呈顯著負相關關系(P<0.05);脯氨酸含量與電解質滲透率呈顯著負相關關系(P<0.05);電解質滲透率與POD活性成顯著負相關關系(P<0.05);其他指標之間相關不顯著。

表 1 綠竹抗性生理指標的相關性分析Table 1 Correlation analysis of resistance physiology indexes of B. oldhamii

3.3 綠竹抗性生理指標的主成分分析

由表2可知,8個主成分的累積貢獻率為100%。根據統計學原理,當各主成分累積貢獻率>85%時,即可反映系統的變異信息[19],因此根據涵蓋85%以上信息的主成分確定為主成分個數。在這8個主成分中,前4個主成分的累積貢獻率達92.036%,可以代表原來8個抗性生理指標92.036%的信息,因此,選擇前4個主成分作為評價的綜合指標。第1主成分特征值為4.646,可以反映總變異58.075%的信息。其中可溶性蛋白含量和POD活性指標的載荷矩陣最高??扇苄缘鞍缀康姆e累能提高細胞的保水能力,對細胞的生命物質及生物膜起到保護作用。POD參與植物光呼吸作用,保證植物具有較強的抗氧化能力。保水處理越佳,可溶性蛋白含量越多,POD活性越高,植物的抗性越強。第2主成分特征值為1.442,可以反映總變異18.023%的信息。其中,MDA含量和可溶性蛋白含量有著較高的載荷值。MDA含量積累越少表明細胞膜受損程度越小,細胞內物質外滲程度越低。保水處理越佳,MDA含量越少,可溶性蛋白含量越多。第3主成分特征值為0.684,可以反映總變異8.551%的信息。其中,MDA含量和脯氨酸含量有著較高的載荷值。脯氨酸含量能夠通過滲透調節來降低水勢,維持較高的滲透壓,保證細胞正常生理功能。保水處理越佳,脯氨酸含量越多,MDA含量越少。第4主成分特征值為0.591,可以反映總變異7.386%的信息。其中可溶性糖含量和脯氨酸含量有著較高的載荷值??扇苄蕴菂⑴c滲透調節,并能夠維持蛋白質的穩定。保水處理越佳,可溶性糖含量和脯氨酸含量越多。

表 2 綠竹葉片生理指標主成分分析Table 2 Principal component results of physiology indexes for B. oldhamii

3.4 綠竹抗性生理指標隸屬函數法的評價結果

根據主成分分析結果,篩選出可溶性蛋白含量、可溶性糖含量、脯氨酸含量和MDA含量作為保水劑與氮肥施用效果綜合評價指標,表明保水劑和氮肥對綠竹滲透調節物質影響較大。其中,可溶性蛋白含量、可溶性糖含量和脯氨酸含量與保水劑與氮肥施用效果呈正相關,用隸屬函數法計算;MDA含量與保水劑與氮肥施用效果呈反相關,用反隸屬函數法計算。用這4個指標的平均隸屬函數值作為保水劑施用效果綜合評判標準,平均值越大,表明保水處理效果越好。保水劑與氮肥處理下綠竹抗性生理指標的綜合評價結果(表 3)為 A2B1>A1B3>A3B2>A2B3>A3B1>A2B2>A3B3>CK>A1B1>A1B2。

表 3 綠竹葉片抗性生理指標的綜合評價Table 3 Comprehensive evaluation results of resistance physiology indexes for B. oldhamii

4 結論與討論

本研究通過對保水劑與氮肥處理下綠竹葉片中8個抗性生理指標的分析表明,不同保水劑與氮肥處理間存在不同程度的差異,但在保水劑與氮肥處理下,各抗性生理指標均有不同程度的提高,與保水劑對仁用杏(Prunus armeniaca)[21]和檳榔(Areca catechu)幼苗[22]生理特性研究結果一致,基本形成了對沿海沙地干旱貧瘠沙地環境的適應性生理特征。保水劑與氮肥的處理提高了綠竹葉片中葉綠素含量,保證光合作用的順利進行,黃慧青等[23]認為在低保水劑用量(10 g/cm2)處理下能夠提高海濱雀稗(Paspalum thunbergii)葉綠素含量,這與本研究結果一致。保水劑與氮肥處理下綠竹在7月、8月、9月中保持較高的葉綠素含量,保證了筍期的光合作用和養分積累,其中除A2B1處理,其他處理均高于CK。

保水劑與氮肥的處理能夠保證綠竹體內主動積累可溶性糖、可溶性蛋白和脯氨酸等滲透調節物質來提高細胞液濃度,降低滲透勢,提高細胞吸水或保水能力。王海珍等[24]研究表明胡楊(Populus euphratica)幼苗以大量積累可溶性糖、維持可溶性蛋白質基本穩定為特征,而灰胡楊(P. pruinosa)幼苗則以積累大量游離脯氨酸來抵抗干旱的環境。關于紅麻(Apocynum venetum)、大麻狀羅布麻(A. cannabinum)[25]、皂莢(Gleditsia sinensis)[26]、油豆角(Phaseolus vulgaris)[27]等研究均表明植物能夠通過提高自身滲透調節物質的積累來抵御一定的干旱脅迫。本研究表明,在7—9月中,綠竹能夠通過積累可溶性糖含量、可溶性蛋白含量和脯氨酸含量,以提高對沙地干旱環境的適應能力??扇苄蕴呛涂扇苄缘鞍淄瑫r作為營養物質供給綠竹筍期生長,因此在出筍盛期(8月),可溶性糖和可溶性蛋白含量積累出現明顯低值。

保水劑與氮肥的處理能夠減輕沙地惡劣環境對植物細胞內活性氧代謝平衡的破環,減少MDA含量、降低電解質滲透率,同時能夠提高細胞內清除活性氧的保護酶SOD和POD活性。周偉偉等[28]對7種景天屬(Sedum)植物研究表明干旱破壞了植物細胞內活性氧產生和清除的動態平衡,MDA含量明顯增加,SOD活性顯著降低,王紅梅等[29]對白刺花(Sophora davidii)研究表明在干旱脅迫下MDA含量逐漸增加,SOD和POD活性隨干旱程度加劇逐漸下降,這與豌豆(Pisum sativum)[30]和辣椒(Capsicum annuum)[31]等研究結果一致。本研究表明,綠竹MDA含量和電解質滲透率在8月積累最多,SOD和POD保護酶活性最低,在保水劑與氮肥施用能夠維持綠竹葉片體內活性氧代謝平衡,降低細胞膜系統損傷。

主成分分析表明可溶性蛋白含量、可溶性糖含量、脯氨酸含量和MDA含量可作為衡量保水劑和氮肥處理下綠竹抗性生理能力的指標,但是從單一指標無法判斷保水劑與氮肥施用對植物抗性生理影響的強弱。目前,隸屬函數法已應用于保水劑施用效果的綜合評價[32],本研究同樣利用隸屬函數法對不同保水劑和氮肥處理進行綜合評價,用可溶性蛋白含量、可溶性糖含量、脯氨酸含量和MDA含量4項指標隸屬度的平均值作為保水劑和氮肥處理優劣綜合評價標準,綜合評 價 結 果 為 A2B1>A1B3>A3B2>A2B3>A3B1>A2B2>A3B3>CK>A1B1>A1B2。在沿海沙地特殊生境中,水肥條件是制約沿海防護林樹種生長發育重要因素,提高土壤肥力并維持土壤含水量是有效保證沿海沙地造林保存率的關鍵,適量的保水劑與氮肥施用能夠提高植物對雨水的儲藏和利用能力,緩解高溫干旱對植物的傷害??赡苡捎诒K畡┡c氮肥的互作效應受氮肥水溶液離子強度制約[33-34],再加上沿海沙地夏季高溫影響較大,使保水劑功能受到抑制,造成A1B1和A1B2的效果低于CK,有關保水劑和氮肥相互作用原理還有待于做進一步研究和探討。

[ 參 考 文 獻 ]

[1]林武星. 沿海木麻黃基干林帶更新方式和更新效果研究 [J]. 江西農業大學學報, 2016, 22(4): 576-580.

[2]林武星. 閩南沿海沙地引種樹種的評價 [J]. 中南林學院學報, 2005, 25(1): 42-45.

[3]譚芳林, 徐俊森, 林武星, 等. 福建濱海沙地造林樹種的適應性與選擇研究 [J]. 林業科學, 2003, 39(S1):100-105.

[4]杜瀾, 陳亮, 謝錦忠, 等. 干旱脅迫及復水對綠竹容器苗生理特征的影響 [J]. 林業科學研究, 2018, 31(4):164-171.

[5]鄭容妹, 鄭郁善, 張梅, 等. 鹽分脅迫對沿海綠竹光合作用及葉綠素的影響 [J]. 竹子學報, 2002, 21(4):76-80.

[6]林愛玉, 涂志華, 上官保國, 等. 沿海沙地引種竹子和木麻黃固沙功能比較研究 [J]. 福建林學院學報,2013, 33(1): 28-33.

[7]張梅, 鄭郁善, 陳禮光. 濱海沙地竹子引種試驗初報 [J]. 西南林學院學報, 2007, 27(1): 48-50.

[8]陳禮光, 鄭郁善, 姚慶端, 等. 沿海沙地新造綠竹林生物量結構 [J]. 福建林學院學報, 2002, 22(3): 249-252.

[9]鄭郁善, 鄭容妹, 姚慶端. 鹽分脅迫對沿海沙地綠竹生理生化的效應 [J]. 福建農林大學學報(自然科學版), 2003, 32(4): 478-481.

[10]涂志華, 潘瑞, 盧仁鋒, 等. 沿海沙地8個竹種根際土壤的酶活性 [J]. 福建農林大學學報(自然科學版),2012, 41(3): 243-247.

[11]鄭晶晶, 蔡錳柯, 林宇, 等. 6種不同沿海防護混交林凋落葉持水性能比較 [J]. 水土保持通報, 2015, 35(5):111-116.

[12]李大培, 李生萍, 高向倩, 等. 安康‘紫仁’核桃響應干旱脅迫的生理評價研究 [J]. 西南林業大學學報(自然科學), 2018, 38(1): 202-206.

[13]馬曉華, 鄭堅, 張旭樂, 等. 干旱脅迫對大花紫薇與毛萼紫薇生理生化特征的影響 [J]. 西南林業大學學報(自然科學), 2018, 38(5): 205-209.

[14]馬丙堯, 劉方春, 馬海林, 等. 保水劑-尿素凝膠對核桃苗根際土壤微生態環境影響 [J]. 中國農學通報,2018, 34(27): 55-59.

[15]宋雙雙, 孫保平, 張建鋒. 保水劑和微生物菌肥對半干旱區造林和土壤改良的影響 [J]. 水土保持學報,2018, 32(3): 334-339.

[16]章桂芝. 沿海沙地竹林施肥效應與經營管理措施研究[D]. 福州: 福建農林大學, 2014.

[17]高俊鳳. 植物生理學[M]. 北京: 高等教育出版社,2006.

[18]李合生. 植物生理生化實驗原理和技術[M]. 北京: 高等教育出版社, 2000.

[19]溫延臣, 李燕青, 袁亮, 等. 長期不同施肥制度土壤肥力特征綜合評價方法 [J]. 農業工程學報, 2015, 31(7):91-99.

[20]吳曉鳳, 倪沛, 楊濤, 等. 10種菊科植物的抗旱性與抗鹽性評價 [J]. 生態學雜志, 2018, 37(7): 1959-1968.

[21]龐海穎, 牛東偉, 李彥慧. 施用保水劑對仁用杏抗旱生理特性的影響 [J]. 河北農業大學學報, 2017, 40(4):43-49.

[22]李佳, 劉立云, 李艷, 等. 保水劑對干旱脅迫檳榔幼苗生理特征的影響 [J]. 南方農業學報, 2018, 49(1):104-108.

[23]黃慧青, 周林濤, 安勐潁, 等. 保水劑對海濱雀稗抗旱性的影響 [J]. 西南農業學報, 2016, 29(8): 1828-1833.

[24]王海珍, 徐雅麗, 張翠麗, 等. 干旱脅迫對胡楊和灰胡楊幼苗滲透調節物質及抗氧化酶活性的影響 [J]. 干旱區資源與環境, 2015, 29(12): 125-130.

[25]王東清, 李國旗, 蘇德喜. 干旱脅迫對兩種羅布麻滲透調節物質積累和保護酶活性的影響 [J]. 干旱區資源與環境, 2012, 26(12): 177-181.

[26]陳麗培, 劉瑞霞, 楊玉珍, 等. 干旱脅迫對刺槐、皂莢葉片滲透調節物質含量及保護酶活性的影響 [J]. 河南農業科學, 2017, 46(10): 122-127.

[27]Farhangi A S, Torabianb S. Antioxidant enzyme and osmotic adjustment changes in bean seedlings as affected by biochar under salt stress [J]. Ecotoxicology and Environmental Safety, 2017, 137: 64-70.

[28]周偉偉, 王雁, 杜靜. 干旱脅迫對景天屬植物生理生化特性的影響 [J]. 林業科學研究, 2009, 22(6):829-834.

[29]王紅梅, 包維楷, 李芳蘭. 不同干旱脅迫強度下白刺花幼苗葉片的生理生化反應 [J]. 應用與環境生物學報, 2008, 14(6): 757-762.

[30]顧建勤, 閆志利, ??×x, 等. 干旱脅迫及復水對豌豆苗期保護酶活性及膜脂過氧化的影響 [J]. 干旱地區農業研究, 2012, 30(3): 159-164.

[31]劉佳, 郁繼華, 徐秉良, 等. 干旱氣候條件下水分脅迫對辣椒葉片生理特性的影響 [J]. 核農學報, 2012,26(8): 1197-1203.

[32]原慧芳, 肖榮才, 黃菁, 等. 東試早柚對不同保水處理的生理響應及綜合評價 [J]. 浙江農業學報, 2016,28(4): 586-594.

[33]趙和瓊, 谷麗萍, 鄭科, 等. 不同集水保水措施對金沙江干熱河谷區林地土壤儲水量的影響 [J]. 西部林業科學, 2018, 47(6): 85-89.

[34]楊磊, 蘇文強. 化肥對保水劑吸水性能的影響 [J]. 東北林業大學學報, 2004, 32(5): 37-38.

猜你喜歡
綠竹保水劑沙地
能鉆過柔軟沙地的蛇形機器人
保水劑吸水保水特性的試驗分析
呼倫貝爾沙地實現良性逆轉
沙地迷宮
綠竹俠
保水劑在山地核桃園中的應用
林業保水劑及其在園林綠化中應用
風滾草
國風·衛風·淇奧
千年指尖過,綠竹仍猗猗
91香蕉高清国产线观看免费-97夜夜澡人人爽人人喊a-99久久久无码国产精品9-国产亚洲日韩欧美综合