劉華君,陳友強,林 明,白曉山,鄧超宏,潘竟海,李承業(yè)
(新疆農(nóng)業(yè)科學院經(jīng)濟作物研究所,烏魯木齊 830091)
【研究意義】目前我國籽用西瓜主要產(chǎn)區(qū)集中在新疆地區(qū),其中塔額盆地縣市加上兵團第九師為我國最大的商品籽用西瓜生產(chǎn)基地,新疆生態(tài)區(qū)屬于大陸性干旱氣候,晝夜溫差大,年均降水量少,氣候冷涼,為籽用西瓜提供了良好的生長環(huán)境,新疆各地州均有大面積種植,據(jù)調(diào)查發(fā)現(xiàn)2008年阿勒泰地區(qū)種植面積4×104hm2(60萬畝)[1]。針對新疆不同生態(tài)環(huán)境下,挖掘、選育生育期適當、經(jīng)濟產(chǎn)量優(yōu)良的優(yōu)質(zhì)籽用西瓜種質(zhì)資源,對實現(xiàn)籽瓜大面積、優(yōu)質(zhì)高效生產(chǎn)具有重要意義?!厩叭搜芯窟M展】研究發(fā)現(xiàn)作物生長發(fā)育是其產(chǎn)量的重要衡量指標[2-5]。俄有浩等[2]研究發(fā)現(xiàn)小麥生長發(fā)育后期,溫度升高不能促進其提前成熟,有效積溫對生育期的影響較大;紀海波等[3]通過對西瓜種質(zhì)資源主要性狀的表型多樣性進行分析,發(fā)現(xiàn)單果質(zhì)量的變異系數(shù)最大,種子百粒重、單瓜種子數(shù)、果實橫徑的遺傳多樣性指數(shù)均較大,分別為2.02、2.01、2.01,而莖斷面形狀的遺傳多樣性指數(shù)最小。趙勝杰等[4]利用SSR標記對無籽西瓜進行指紋圖譜構(gòu)建及遺傳多樣性分析,發(fā)現(xiàn)核心SSR標記適用于構(gòu)建無籽西瓜品種的DNA指紋圖譜,我國主栽無籽西瓜品種遺傳相似度較高,亟需加強種質(zhì)創(chuàng)新。呂尊富等[5]利用WheatGrow、CERES模型將小麥的生育期模型擴大,進一步量化中國小麥主產(chǎn)區(qū)的區(qū)域生產(chǎn)力。有關(guān)從生育期、經(jīng)濟產(chǎn)量等方面分析作物生態(tài)適應性已有大量報道,龐丹波[6],張曉杰[7],吳舒致等[8]基于主成分分析、隸屬函數(shù)法對作物種質(zhì)資源進行綜合評價;徐守俊等[10]利用主成分分析將參試品種分成2大類,第一類大豆品種有籽粒具有較高的鉬富集能力;第二類豆莢獲取鉬能力較強?!颈狙芯壳腥朦c】新疆不同生態(tài)區(qū)差異較大,對不同籽瓜的生育期、產(chǎn)量性狀要求各不相同,但目前有關(guān)新疆生態(tài)區(qū)籽瓜優(yōu)良品系篩選方法報道較少。研究選育出適宜新疆不同生態(tài)區(qū)的優(yōu)良籽瓜種質(zhì)資源?!緮M解決的關(guān)鍵問題】基于新疆生態(tài)區(qū)籽瓜生長發(fā)育時期、經(jīng)濟產(chǎn)量的不同,利用主成分分析對參試品種進行綜合評價,為加快優(yōu)質(zhì)籽瓜選育進展提供一定理論依據(jù)。
試驗分別在新疆農(nóng)業(yè)科學院瑪納斯農(nóng)業(yè)試驗站和新疆昌吉市農(nóng)業(yè)技術(shù)推廣中心試驗基地進行,這兩個地點代表新疆北疆天山北坡沖積平原生態(tài)類型區(qū),屬于較為典型的溫帶大陸性干旱半干旱氣候,年平均氣溫6.8 ℃,年平均降雨量230 mm,年日照時數(shù)2 833 h,有效積溫3 450℃,年均蒸發(fā)量1 787 mm,平均無霜期180 d。列出參試品種(系)11個,均來源于新疆農(nóng)科院經(jīng)作所,其列出特性。表1
表1 參試品種(系)來源及特性
Table 1 Source and characteristics of tested varieties (lines)
品種(系)Varieties(lines)皮色Skincolor籽色Seedcolor大小Size1601青網(wǎng)紅中1602青皮黑小1603白皮黑中1604白碎黑大1607白皮黑中1608青皮黑小1609青網(wǎng)紅中1612白碎黑大1613白皮黑中1614白碎黑大1617青網(wǎng)紅中1618青皮黑小
1.2.1 試驗設計
于2016~2017年各試點分別按照當?shù)剡m宜的播種時間播種,采用隨機區(qū)組排列, 3 次重復,行距1 m,株距25 cm,播種時結(jié)合整地施45 000 kg/hm2農(nóng)家肥和450 kg/hm2磷酸二銨作底肥,三葉期間苗,定苗,留苗4×104株/hm2,6月中旬澆頭水,隨水適當追施尿素,全生育期灌水6次,8月下旬停水。6月下旬和7月上旬分別噴施氟硅唑防治白粉病,9月初收獲。
1.2.2 生育期統(tǒng)計及測定指標
生育期統(tǒng)計指標有出苗期(X1)、伸蔓期(X2)、始花期(X3)、盛花期(X4)、坐果期(X5)、定果期(X6)、成熟期(X7)、全生育期(X8);產(chǎn)量測定指標有粒數(shù)(X9)、粒重(X10)、百粒重(X11)。
所有數(shù)據(jù)用Excel 2010 和SPSS 19.0 統(tǒng)計分析軟件進行相關(guān)性分析及主成分分析。
研究表明,參試品種(系)不同生育時期差異較大,其中出苗期時間最長的是1603,平均11.3 d,其次是1604,平均10.67 d;時間最短的是1607、1613,分別為7.57 d;伸蔓期時間最長的是1602,平均21.67 d;其次是1603,平均為21.00 d;時間最短的是1608,平均17.67 d;盛花期時間最長的是1604,平均9 d,其次是1608,平均8.67 d;時間最短的是1601,平均5.67 d;成熟期較其他生育時期,時間最長,均在50 d左右,其中1608成熟期平均52.33 d,1618成熟期平均50.33 d;參試品種中全生育期時間最長的為1612,平均119.33 d,其次是1614,平均119.00 d;生育期最短的是1607、1608,均為112.00 d。表2
研究表明,其中1601、1617粒數(shù)最多,分別360.0、360.67個;1601、1617粒重最大,分別為62.01、61.67 g;其次是1609粒數(shù)為359.67個,粒重為60.67 g;1607的粒數(shù)最少81.33 個,粒重16.98 g;品系1614、1613、1604的百粒重最高,分別平均為24.53、24.52和24.46 g;百粒重最低的是1608,平均為14.33 g。表3
表3 不同籽瓜品種(系)經(jīng)濟產(chǎn)量
Table 3 The economic yield of different seed-watermelon lines (varieties)
品種(系)Lines(Varieties)粒數(shù)Theseednumbe(個)粒重Grainweight(g)百粒重100grainweight(g)1601360.00±0.58a62.01±0.57a16.88±0.05c1602355.67±0.33b51.18±0.43b14.4±0.15d160381.33±0.67d17.64±0.32d21.28±0.3b1604138.33±0.33c33.83±0.6c24.46±0.32a160781.33±0.67d16.98±0.56d20.69±0.59b1608356.33±0.88b51.51±0.49b14.33±0.13d1609359.67±0.33a60.67±0.34a16.92±0.06c1612137.67±0.88c33.17±0.6c24.52±0.28a161382.00±0.58d17.31±0.66d21.23±0.33b1614138.33±0.33c34.17±0.44c24.53±0.24a1617360.67±0.33a61.67±0.33a16.63±0.32c1618354.67±1.33b51.18±0.43b14.32±0.25d
注:小寫字母表示不同處理間差異5%水平顯著
Note: lower case letters indicate that the difference between different treatments is 5% level significant
對不同籽瓜品種生育期與經(jīng)濟產(chǎn)量的相關(guān)性進行分析,研究表明,伸蔓期(X2)與盛花期(X4)呈極顯著負相關(guān),相關(guān)系數(shù)-0.75;出苗期(X1)與坐果期(X5)呈顯著正相關(guān),相關(guān)系數(shù)0.70;全生育期(X8)與坐果期(X5)呈極顯著正相關(guān),相關(guān)系數(shù)0.76;與定果期(X6)呈顯著正相關(guān),相關(guān)系數(shù)0.61;粒數(shù)(X9)與始花期(X3)呈顯著負相關(guān),相關(guān)系數(shù)-0.64;粒重(X10)與始花期(X3)呈顯著負相關(guān),相關(guān)系數(shù)-0.62;粒重(X10)與粒數(shù)(X9)呈顯著正相關(guān),相關(guān)系數(shù)0.96;百粒重(X11)與全生育期(X8)呈顯著正相關(guān),相關(guān)系數(shù)為0.62,與粒數(shù)(X9)、粒重(X10)呈極顯著負相關(guān),相關(guān)系數(shù)分別為-0.089、-0.74。說明籽瓜不同時期的長短,對其全生育期、經(jīng)濟產(chǎn)量影響較大,高產(chǎn)籽瓜的始花期較短,且后期成熟期時間較長。表4
對參試品種的11個指標進行主成分分析,得到各提取成分的初始化特征值,研究表明,參試品種的11個指標的累積貢獻率達到100%,說明11個指標對籽瓜的綜合性狀均有貢獻。根據(jù)主成分對應特征值大于1的原則提取主成分,提取到4 個主成分,其中第1主成分的貢獻率為32.742%,第2主成分貢獻率高達53.418%,前4個主成分累積貢獻率為81.620%,說明前4個主成分代表了籽瓜主要性狀的85%左右的信息,因此,可以用前4個主成分來反映11個原始指標的基本信息。表5
表4 不同籽瓜品種生育期、經(jīng)濟產(chǎn)量相關(guān)系數(shù)
Table 4 The correlation coefficient of growth period and economic yield of different seed melon varieties
指標IndexX1X2X3X4X5X6X7X8X9X10X11X11X2-0.131X3-0.320.221X40.09-0.75**0.201X50.70*-0.25-0.120.231X60.070.150.53-0.080.341X70.13-0.52-0.300.290.43-0.061X80.560.040.400.260.76**0.61*0.301X9-0.210.13-0.64*-0.47-0.10-0.450.45-0.401X10-0.080.07-0.62*-0.480.11-0.270.57-0.210.96**1**X110.40-0.220.560.400.440.62*-0.190.62*-0.89**-0.74**1
注:*在 0.05水平(雙側(cè))上顯著相關(guān),**在0.01 水平(雙側(cè))上顯著相關(guān)
Note:*significant correlation at 0.05 level (bilateral),**significant correlation at 0.01 level (bilateral)
表5 不同參試品種生育期、經(jīng)濟產(chǎn)量主成分分析
Table 5 Principal component analysis of growth period and economic yield of different
成分Component特征值Eigenvalue方差Variance(%)累積貢獻率Cumulativecontributionrate(%)13.60232.74232.74222.27420.67653.41831.85216.83470.25241.25011.36881.62050.7526.83588.45460.5104.63793.09270.4524.11197.20380.2111.91999.12290.0760.69199.813100.0180.16599.979110.0020.021100.000
籽瓜生育期、經(jīng)濟產(chǎn)量成分矩陣反映看主要指標在主成分上的載荷,研究表明,決定第一主成分的主要指標有粒數(shù)(X9)、百粒重(X11)、粒重(X10)、始花期(X3),說明主成分1是由上述指標組成的一個綜合指標,其中粒數(shù)的載荷值最大,為0.913,其次是百粒重、粒重、始花期,說明籽瓜進入生殖生長期越早,經(jīng)濟產(chǎn)量越高;第2主成分貢獻率20.676%,其中成熟期(0.744)、坐果期(0.720)、全生育期(0.633)、出苗期(0.534)為相對較高的正向載荷,說明第2主成分為生育期因子其決定性作用,為生育期因子;第3主成分中貢獻率為16.834%,其中載荷較高且符號為正的指標有伸蔓期、定果期、全生育期、坐果期,其中伸蔓期載荷較高(0.813),說明第3主成分為營養(yǎng)生長決定因子;第4主成分貢獻率為11.368%,其中載荷較高且符號為正的指標有成熟期(0.387)、盛花期(0.273),說明在籽瓜花期、成熟期對其綜合性狀影響較大。表6
表6 不同參試品種生育期和經(jīng)濟產(chǎn)量成分矩陣
Table 6 The component matrix of growth period and economic yield of different varieties
指標Index成份 Component1234X90.9130.3390.0430.054X110.9130.0760.101-0.030X100.8130.4960.1310.078X30.644-0.2140.2130.126X7-0.0660.744-0.3870.387X50.1680.7200.331-0.198X80.5970.6330.3680.017X2-0.197-0.2230.813-0.135X40.5220.144-0.6250.273X10.3170.534-0.103-0.700X60.2390.1790.5650.673
研究表明,參試品種(系)中,第1主成分經(jīng)濟產(chǎn)量得分最高的是1614、1612,均為1.20,其次是1613(0.87);第2主成分生育期因子得分最高的是1614(1.17),其次是1612(0.97);第3主成分營養(yǎng)生長決定因子得分最高的是1618(1.11),其次是1601(0.81);第4主成分得分最高的是1613(1.66),其次是1617(1.32)。而綜合得分最高的是1612(0.65),其次是1614(0.63),說明1612、1614品系經(jīng)濟產(chǎn)量、生育期性狀優(yōu)良,具有較高的產(chǎn)量潛力,1618的經(jīng)濟產(chǎn)量較低,但其營養(yǎng)生長特性優(yōu)良,具有提高產(chǎn)量的潛力。表7
表7 不參試品種(系)生育期和經(jīng)濟產(chǎn)量各主成分得分及綜合得分
Table 7 The principal component scores and comprehensive scores of growth period and economic yield of non tested varieties (lines)
品種(系)Varieties(lines)主成分因子Componentfactor1234綜合得分Generalscore等級Grade1601-1.140.130.81-0.40-0.2691602-1.06-0.360.110.30-0.371116031.22-0.440.21-1.580.17516040.780.450.28-0.540.33316070.29-2.150.24-0.09-0.32101608-1.070.61-2.58-0.76-0.74121609-1.070.290.54-0.07-0.21816121.200.970.37-0.060.65116130.87-1.29-1.001.660.04616141.201.17-0.240.280.6321617-0.360.720.141.320.2141618-0.88-0.101.11-0.06-0.137
生育期適當、經(jīng)濟產(chǎn)量高是選擇優(yōu)良品種(系)的重要因素,發(fā)掘出生育期適當?shù)钠贩N(系)是品種(系)實現(xiàn)高產(chǎn)的基礎[9、10],國內(nèi)外研究者對小麥[11]、花生[12、13]、玉米[14-16]、棉花等[17-18]作物的生育期、經(jīng)濟產(chǎn)量進行了大量研究,但是對籽瓜的生育期、經(jīng)濟產(chǎn)量的綜合分析則鮮見報道,研究通過對籽瓜生育期、經(jīng)濟產(chǎn)量等11個指標進行相關(guān)性分析,發(fā)現(xiàn)籽瓜百粒重與全生育期呈顯著正相關(guān),粒數(shù)與始花期呈顯著負相關(guān);粒重與始花期呈顯著負相關(guān),說明籽瓜較高的經(jīng)濟產(chǎn)量與其各個生長發(fā)育時期緊密相關(guān)。研究發(fā)現(xiàn)作物的生育期、經(jīng)濟產(chǎn)量性狀等11個指標品種(系)之間存在明顯的遺傳差異,但可以用4個主成分來反映11個原始指標的基本信息。
根據(jù)主成分綜合得分情況,對參試品種進行評價,綜合得分越高說明該品種(系)生育期適當、經(jīng)濟產(chǎn)量優(yōu)良,具有較高的產(chǎn)量潛力,其中1614、1612,第1、2主成分得分最高,第3、4主成分得分較低,但綜合得分分別為第1、第2,選育優(yōu)良籽瓜品種(系)時,不能單獨從經(jīng)濟產(chǎn)量、營養(yǎng)生長狀況、生育期進行評價,而應該對其進行綜合全面評價。
籽瓜生長發(fā)育過程中經(jīng)濟產(chǎn)量與其各個生長發(fā)育時期緊密相關(guān),百粒重與全生育期呈顯著正相關(guān),粒數(shù)與始花期呈顯著負相關(guān),粒重與始花期呈顯著負相關(guān)。參試品種(系)1614、1612綜合評價分別排名第1名、第2名,其中,品種(系)1614生育期119 d、單瓜粒數(shù)138.33個、籽粒重34.17 g、百粒重24.53 g,品種(系)1612生育期119.33 d、單瓜粒數(shù)137.67個、籽粒重33.17 g、百粒重24.52 g,兩品種(系)綜合性狀優(yōu)良,具有豐產(chǎn)潛力,可作為重點材料加以利用。
參考文獻(References)
[1] 潘竟海,陳友強,林明,等. 黑籽瓜產(chǎn)量性狀的相關(guān)與通徑分析[J]. 新疆農(nóng)業(yè)科學,2013,50(8):1 406-1 410.
PAN Jing-hai, CHEN You-qiang, LIN Ming, et al. (2013). White Xiaoshan, correlation and path analysis of yield traits of black seed melon [J].XinjiangAgriculturalSciences, 50(8): 1,406-1,410. (in Chinese)
[2] 俄有浩,霍治國,馬玉平, 等. 中國北方春小麥生育期變化的區(qū)域差異性與氣候適應性[J]. 生態(tài)學報,2013,33(19):6 295-6 302.
E You-hao, HUO Zhi-guo, Ma Yu-ping, et al. (2013). Regional variation and climatic adaptability of spring wheat growth stages in northern China [J].JournalofEcology, 33(19):6,295-6,302. (in Chinese)
[3] 紀海波,張玉鑫,李玉明, 等.西瓜種質(zhì)資源主要性狀的表型多樣性[J].西北農(nóng)林科技大學學報(自然科學版),2013,41(8):155-160.
JI Hai-bo, ZHANG Yu-xin, LI Yu-ming, et al. (2013). The phenotypic diversity of main characters of Watermelon Germplasm Resources [J].JournalofNorthwestAgricultureandForestryUniversity(NaturalScienceEdition) , 41(8): 155-160. (in Chinese)
[4] 趙勝杰,朱紅菊,路緒強,等.無籽西瓜品種SSR指紋圖譜構(gòu)建及遺傳多樣性分析[J].植物遺傳資源學報,2013,14(6):1 142-1 146.
ZHAO Sheng-jie, ZHU Hong-ju, LU Xu-qiang, et al. (2013). Seedless watermelon varieties SSR fingerprinting construction and genetic diversity analysis [J].JournalofPlantGeneticResources, 14(6): 1,142-1,146. (in Chinese)
[5] 呂尊富,劉小軍,湯亮, 等. 基于WheatGrow和CERES模型的區(qū)域小麥生育期預測與評價[J]. 中國農(nóng)業(yè)科學,2013,46(6):1 136-1 148.
Lü Zun-fu, LIU Xiao-jun, TANG Liang, et al. (2013). Regional Prediction and Evaluation of Wheat Phenology Based on the WheatGrow and CERES Models [J].ChineseAgriculturalScience, 46(6):1,136-1,148. (in Chinese)
[6] 龐丹波,李生寶,潘占兵, 等. 基于主成分分析和隸屬函數(shù)的紫花苜蓿引種初步評價[J]. 西南農(nóng)業(yè)學報,2015,28(6):2 815-2 819.
PANG Dan-bo, LI Sheng-bao, PAN Zhan-bing, et al. (2015). Preliminary evaluation of Alfalfa Introduction Based on principal component analysis and membership function [J].SouthwestAgriculturalJournal, 28(6): 2,815-2,819. (in Chinese)
[7] 張曉杰,姜慧芳,任小平, 等. 中國花生核心種質(zhì)的主成分分析及相關(guān)分析[J]. 中國油料作物學報,2009,31(3):298-304.
ZHANG Xiao-jie, JIANG Hui-fang, REN Xiao-ping, et al. (2009). Principal component analysis and correlation analysis of peanut core collection in China [J].JournalofChineseJournalofOilCrops, 31(3): 298-304. (in Chinese)
[8] 吳舒致,黎裕. 谷子種質(zhì)資源的主成分分析和圖論主成分分類[J]. 西北農(nóng)業(yè)學報,1997,(2):49-51, 53.
WU Shu-zhi, LI Yu. (1997). Principal component analysis and principal component classification of graph theory for Foxtail Millet Germplasm [J].NorthwestAgriculturalJournal, (2): 49-51, 53. (in Chinese)
[9] Tan, K. C., Lim, H. S., & Jafri, M. Z. M. (2016). Prediction of column ozone concentrations using multiple regression analysis and principal component analysis techniques: a case study in peninsular malaysia.AtmosphericPollutionResearch, 7(3): 533-546.
[10] 徐守俊,孫學成,胡承孝, 等. 不同大豆品種鉬富集因素的主成分分析和聚類分析[J]. 華北農(nóng)學報,2016,31(S1):399-404.
XU Shou-jun, SUN Xue-cheng, HU Cheng-xiao, et al. (2016). Principal component analysis and cluster analysis of molybdenum enrichment factors in different soybean varieties [J].ActaSinicaSinica, 31(S1): 399-404. (in Chinese)
[11] 李桂萍,張根生,巴青松,等. 雜種小麥品質(zhì)性狀的性狀相關(guān)和主成分分析[J]. 浙江農(nóng)業(yè)學報,2016,28(9):1 447-1 453.
LI Gui-ping, ZHANG Gen-sheng, BA Qin-song, et al. (2016). Quality traits of wheat hybrids trait correlation and principal component analysis [J].ZhejiangJournalofAgriculturalSciences, 28(9): 1,447-1,453. (in Chinese)
[12] 李玉發(fā),竇忠玉,梁軍,等. 花生主要農(nóng)藝性狀的遺傳變異及相關(guān)性和主成分分析[J]. 遼寧農(nóng)業(yè)科學,2013,(3):11-14.
LI Yu-fa, DOU Zhong-yu, LIANG Jun, et al. (2013). Ho Chinese. Main Agronomic Traits of peanut and principal component analysis of genetic variation and correlation of [J].LiaoningAgriculturalSciences, (3): 11-14. (in Chinese)
[13] 殷冬梅,張幸果,王允,等. 花生主要品質(zhì)性狀的主成分分析與綜合評價[J]. 植物遺傳資源學報,2011,12(4):507-512, 518.
YIN Dong-mei, ZHANG Xing-guo, WANG Yun, et al. (2011). plant genetic resources and comprehensive evaluation of principal component analysis of main quality traits in peanut, [J].JournalofPlantGeneticResources, 12(4): 507-512, 518. (in Chinese)
[14] 徐云姬,顧道健,秦昊,等. 玉米灌漿期果穗不同部位籽粒碳水化合物積累與淀粉合成相關(guān)酶活性變化[J]. 作物學報,2015,41(2):297-307.
XU Yun-ji, GU Dao-jian, QIN Hao, et al. (2015). Changes of carbohydrate accumulation and enzyme activities related to starch synthesis in different parts of ear during maize filling stage [J].ActaSinica, 41(2): 297-307. (in Chinese)
[15] 陳傳永,王榮煥,趙久然,等. 不同生育時期遮光對玉米籽粒灌漿特性及產(chǎn)量的影響[J]. 作物學報,2014,40(9):1 650-1 657.
CHEN Chuan-yong, WANG Rong-huan, ZHAO Jiu-ran, et al. (2014). Effects of shading at different growth stages on Grain Filling Characteristics and yield of maize [J].ActaAgronomicaSinica, 40(9): 1,650-1,657. (in Chinese)
[16] 陳群,耿婷,侯雯嘉,等. 近20年東北氣候變暖對春玉米生長發(fā)育及產(chǎn)量的影響[J]. 中國農(nóng)業(yè)科學,2014,47(10):1 904-1 916.
CHEN Qun, GENG Ting, HOU Wen-jia, et al. (2014). Effects of climate warming in Northeast China on growth and yield of Spring Maize in recent 20 years [J].ChineseAcademyofAgriculturalSciences, 47(10): 1,904-1,916. (in Chinese)
[17] 鄧忠,翟國亮,呂謀超,等. 施肥策略對新疆棉花產(chǎn)量、品質(zhì)與水氮利用的影響[J/OL]. 排灌機械工程學報,2017,(10):1-5.
DENG Zhong, ZHAI Guo-liang, Lü Mou-chao, et al. (2017).Effects of fertilization strategies on cotton yield, quality and water and nitrogen use in Xinjiang [J/OL].JournalofIrrigationandDrainageMachineryEngineering, (10): 1-5. (in Chinese)
[18] 崔日鮮,陳明利. 山東棉花生育期降水量和作物需水量的時空變化特征研究[J]. 青島農(nóng)業(yè)大學學報(自然科學版),2014,31(4):255-265.
CUI Ri-xian, CHEN Ming-li. (2014). Shandong cotton growing period precipitation and evapotranspiration temporal variation characteristics [J].JournalofQingdaoAgriculturalUniversity(NaturalScienceEdition) , 31(4): 255-265. (in Chinese)