江厚龍，劉國順，楊 超，楊永鋒，胡宏超，劉清華

（1.重慶煙草科學(xué)研究所，重慶 400715；2.河南農(nóng)業(yè)大學(xué)煙草學(xué)院，國家煙草栽培生理生化研究基地，鄭州 450002；3.河南省煙草公司平頂山分公司，河南 平頂山 467000）

土壤作為一個時空連續(xù)的變異體，即便在一個很小范圍內(nèi)也具有高度的空間異質(zhì)性[1-2]。如果忽視土壤養(yǎng)分的這一特性而進行均勻施肥，則會造成高肥區(qū)施肥量過多，低肥區(qū)施肥量過少[3]，從而造成作物產(chǎn)量和品質(zhì)的嚴重下降。因此，根據(jù)土壤養(yǎng)分含量和作物需肥情況來確定施肥量是非常必要的?？緹熓俏覈匾慕?jīng)濟作物，其生長發(fā)育和產(chǎn)質(zhì)量形成受土壤養(yǎng)分和施肥的影響較大[4]。因此，依據(jù)土壤養(yǎng)分情況將煙田劃分為不同肥力水平的區(qū)域進行精準管理是生產(chǎn)優(yōu)質(zhì)煙葉的關(guān)鍵。目前，國內(nèi)外關(guān)于植煙土壤分區(qū)管理方面的研究報道較少，僅見基于土壤養(yǎng)分條件進行煙田養(yǎng)分分區(qū)管理研究的報道[5]?？緹熓且环N特殊的經(jīng)濟作物，就工業(yè)和消費者而言，質(zhì)量比產(chǎn)量更加重要[6]。因此，必要基于烤煙品質(zhì)因素進行分區(qū)精準管理，目前國內(nèi)外尚未報道相關(guān)的研究。另外，在管理分區(qū)方面的報道中，僅對管理分區(qū)進行定義[7-8]，偶有研究建立了施肥模型[9]，但未見將施肥模型應(yīng)用于管理分區(qū)并驗證其效果的研究報道。

筆者運用地理信息系統(tǒng)（GIS）、全球定位系統(tǒng)（GPS）和地統(tǒng)計學(xué)等方法，研究了試驗區(qū)土壤養(yǎng)分的空間分布規(guī)律，定義了基于烤煙品質(zhì)的管理分區(qū)；并通過多年布點試驗建立推薦施肥模型，將模型應(yīng)用于各管理分區(qū)，同時對模型在各分區(qū)的應(yīng)用效果進行驗證。旨在通過推薦施肥減少施肥的盲目性、提高農(nóng)業(yè)生產(chǎn)的科學(xué)水平、促進農(nóng)業(yè)資源高效持續(xù)利用；同時，通過減少農(nóng)業(yè)肥料投入量，降低農(nóng)業(yè)生產(chǎn)對生態(tài)環(huán)境的威脅。本研究結(jié)果將為煙田土壤養(yǎng)分分區(qū)管理和推薦施肥模型提供理論依據(jù)，為精準農(nóng)業(yè)研究開辟新思路。

1 材料與方法

1.1 研究區(qū)概況

本研究在平頂山市郟縣典型煙草種植地塊進行（34°04′36″N，113°05′04″E），面積 4 hm2，以烤煙（中煙100品種）為主要種植作物，土壤類型為褐土。該區(qū)屬暖溫帶大陸性季風(fēng)氣候，年均氣溫14.6℃，≥10 ℃年活動積溫4734.9 ℃，無霜期220 d左右，年均降雨量680 mm。

1.2 土樣采集與分析

本研究于2009年3月上旬整地施肥前進行土樣采集，利用GPS對每樣點進行定位；采用“網(wǎng)格法”取土壤樣品，間隔為20 m。采樣時以網(wǎng)格點為圓心的5 m范圍內(nèi)采集10鉆0～20 cm的土壤組成代表該點的混合樣，共采集樣品111個（圖1）。

圖1 試驗區(qū)位置及土壤樣點（a）和煙葉烤后樣取樣點（b）分布Fig.1 The location of the study field,spatial distribution pattern of sampling points for soil (a) and flue-cured tobacco (b) (a:111 soil samples; b:39 flue-cured tobacco samples at the corresponding 39 soil samples)

土樣經(jīng)風(fēng)干、去雜、磨碎、過篩后備用。分別測定土壤pH（pH計法）、有機質(zhì)（OM）[10]、活性有機質(zhì)（ASOM）[11]、總氮（TN）[12]、堿解氮（AN）[13]、速效磷（AP）[14]、速效鉀（AK）[15]和陽離子交換量（CEC）[16]等土壤屬性指標。

1.3 試驗方案設(shè)計

2007、2008、2009和2010年分別在平頂山郟縣選取 10個代表性煙田布置大田試驗。每試驗點設(shè)7個處理，3次重復(fù)，隨機區(qū)組排列，小區(qū)面積133.4 m2，各處理施肥量見表1。試驗所施肥料為煙草專用復(fù)合肥m(N):m(P2O5):m(K2O)=10:10:20、硝銨（含N 30%）、重過磷酸鈣（含P2O5為44%）、K2SO4（含K2O為50%）、KNO3（含N為13%，含K2O為 45%）。磷肥作基肥一次性條施；氮肥和鉀肥 70%作基肥，10%作為穴肥，20%在團棵前分 2次追施。供試品種為中煙100。

表1 試驗各處理肥料施用量 kg/hm2Table 1 The fertilizer rates for each treatment

1.4 分區(qū)方法

煙葉樣品于2009年9月獲取?？緹煶墒炱?，從111個土樣點上隨機選取39個（圖1b），采收充分成熟的中部葉進行烘烤，每樣點取 1 kg烤后樣（C2F）。對樣品的外觀質(zhì)量、物理特性和化學(xué)成分及香氣等指標進行測定，并對各樣品進行品吸。將各指標進行量化處理給出分值，并賦予權(quán)重[17]，計算最終的品質(zhì)得分。

1.5 分區(qū)驗證方法

2010年4月在各分區(qū)內(nèi)分別選2個代表性地塊，一塊進行推薦施肥，另一塊常規(guī)施肥（對照）。田間管理按當(dāng)?shù)刈顑?yōu)進行，于9月取各處理充分成熟的中部葉進行烘烤，取不同處理烤后樣（C3F）各2 kg進行分析。

1.6 數(shù)據(jù)處理

本研究用SPSS 10.0軟件進行異常值剔除、描述性統(tǒng)計和逐步回歸分析；用GS+軟件[18]進行地統(tǒng)計分析；模糊聚類分析和最佳分區(qū)數(shù)確定用 MZA軟件[19]進行；利用ArcGIS9.2軟件進行圖層編輯和輸出。

2 結(jié) 果

2.1 土壤養(yǎng)分測試數(shù)據(jù)處理

2.1.1 描述性統(tǒng)計 樣本平均值、標準差、變異系數(shù)、最大值及定義概率等廣泛應(yīng)用于數(shù)據(jù)分析中，是反映數(shù)據(jù)整體分布特征和數(shù)據(jù)分析的第一步，也是對數(shù)據(jù)進一步分析的基礎(chǔ)[20]。對8種土壤屬性進行統(tǒng)計分析，結(jié)果見表2。

由表2可知，該區(qū)域土壤pH變幅為7.39～8.44（=7.93），變幅較小，為普遍認同的烤煙生長pH次適宜區(qū)。OM含量為17.19 g/kg，適宜烤煙生長，可有效提高煙葉香氣質(zhì)和香氣量，減少雜氣和刺激性[21-22]。TN 含量為 0.54～1.40 g/kg，為中等偏低[23]。土壤AN、AP和AK是烤煙生長的主要養(yǎng)分因子[5]，其平均值分別為 74.97、7.39、161.40 mg/kg；AN和AP含量偏低，AK含量為較高水平；AP在該區(qū)內(nèi)分布不均，其變異性較大（CV=43.85%）。許多學(xué)者也報道了相似的結(jié)果[24-26]。表2顯示，大部分土壤養(yǎng)分具有較大的變異性，因此有必要將該區(qū)劃分為不同養(yǎng)分水平的分區(qū)進行精準管理，以減少因肥料因素而造成煙葉產(chǎn)量和品質(zhì)的降低。

由偏度、峰度和Kolmogorov-Smirnov（P＞0.05）檢驗表明，所測土壤指標均為正態(tài)分布或近似符合正態(tài)分布。因此數(shù)據(jù)均能滿足使用地統(tǒng)計分析要求。

表2 土壤養(yǎng)分的描述性統(tǒng)計Table 2 Descriptive statistics of soil nutrients in study area (n=111)

2.1.2 土壤養(yǎng)分和烤煙品質(zhì)空間變異性分析 在忽視各向異性和趨勢效應(yīng)的前提下，擬合了8種養(yǎng)分和烤煙品質(zhì)的半方差函數(shù)模型（表3）。結(jié)果顯示，球狀模型能很好地擬合8種土壤養(yǎng)分，而烤煙品質(zhì)用線性模型能較好擬合，模型能真實地反應(yīng)空間分布情況。CEC、TN和烤煙品質(zhì)的塊金值與基臺值比小于25%，具有較強的空間自相關(guān)性[27]，說明這2種養(yǎng)分及烤煙品質(zhì)的空間變異性主要由結(jié)構(gòu)性變異決定。其他土壤養(yǎng)分的空間變異性均由結(jié)構(gòu)性因素和隨機性因素共同決定。利用半方差函數(shù)模型對土壤養(yǎng)分和烤煙品質(zhì)進行Kriging插值（圖2和圖4A）。由圖可知，烤煙品質(zhì)與土壤養(yǎng)分空間分布極為相似，尤其是與AK和AP。另外，圖2顯示，AK含量屬中等偏上，OM和AP的含量中等偏低，CEC含量表明該研究區(qū)土壤具有中等保肥供肥能力。

表3 土壤屬性半方差函數(shù)理論模型及有關(guān)參數(shù)Table 3 Theoretical models and corresponding parameters for semivariogram of soil properties

圖2 土壤養(yǎng)分空間分布特征圖Fig.2 Spatial distribution maps of soil properties

2.2 土壤養(yǎng)分綜合分區(qū)

2.2.1 分區(qū)數(shù)據(jù)源提取 為了研究影響該區(qū)烤煙品質(zhì)的關(guān)鍵土壤養(yǎng)分因子，利用逐步回歸法進行分析。將土壤 pH、OM、ASOM、TN、AN、AP和CEC作為自變量，烤煙品質(zhì)得分作為因變量，進行逐步回歸分析，得如下回歸方程：

F檢驗顯示，由逐步回歸所建立的方程達到顯著水平（P=0.03），決定系數(shù)R2=0.746，Durbin-Waston 統(tǒng)計量d=1.543，說明該模型擬合度較好。由方程回歸系數(shù)可知，AK、OM、AP和CEC對烤煙品質(zhì)有較大影響，是該區(qū)烤煙品質(zhì)形成的關(guān)鍵土壤養(yǎng)分因子；其他3種土壤屬性對烤煙品質(zhì)無顯著影響。

2.2.2 分區(qū)與驗證 將4種關(guān)鍵養(yǎng)分因子的空間分布圖轉(zhuǎn)化成5 m×5 m柵格圖，通過轉(zhuǎn)點運算提取每個柵格的插值及其坐標，將柵格屬性作為模糊聚類的數(shù)據(jù)源應(yīng)用到 MZA軟件中進行模糊聚類分析。在模糊聚類時軟件計算出每分區(qū)對應(yīng)的 FPI（Fuzziness performance index） 和 NCE（Normalized classification entropy），其值越小分區(qū)效果越好，同時達到最小時對應(yīng)的分區(qū)數(shù)為最佳分區(qū)數(shù)[19]。由圖3可知，該區(qū)域最佳分區(qū)數(shù)為4，分區(qū)結(jié)果見圖4B。圖4顯示，分區(qū)結(jié)果與烤煙品質(zhì)分布具較高吻合度，說明分區(qū)結(jié)果是合理的。

圖3 不同分區(qū)數(shù)所對應(yīng)的FPI和NCE值Fig.3 FPI and NCE shown as a function of the number of management zone.

圖4 基于煙葉質(zhì)量的土壤養(yǎng)分分區(qū)管理圖Fig.4 Soil nutrient management zone map base on the tobacco-quality

為了驗證基于煙葉品質(zhì)的管理分區(qū)能否有效地量化作物生長的關(guān)鍵土壤養(yǎng)分，進而有效地進行精準管理，對各分區(qū)的土壤養(yǎng)分和煙葉品質(zhì)進行分析，結(jié)果見表4、5。表4顯示，各土壤養(yǎng)分和煙葉品質(zhì)在不同分區(qū)間具顯著差異，說明分區(qū)致使土壤養(yǎng)分和煙葉品質(zhì)在分區(qū)間的差異性趨于增大。由表5可知，分區(qū)后土壤養(yǎng)分含量和煙葉品質(zhì)在分區(qū)內(nèi)趨于一致。因此，認為該分區(qū)結(jié)果是合理的，可作為精準管理依據(jù)。

表4 各管理分區(qū)內(nèi)土壤養(yǎng)分和煙葉品質(zhì)的差異性檢驗Table 4 One-way variance analysis of soil nutrients and tobacco quality in the four management zones.

表5 管理分區(qū)內(nèi)土壤養(yǎng)分和煙葉品質(zhì)得分的變異系數(shù) %Table 5 CV of soil nutrient and tobacco quality in the four management zones %

2.3 施肥模型建立

根據(jù)4年10點大田試驗數(shù)據(jù)，計算土壤養(yǎng)分效率、肥料利用率、需肥量，綜合考慮烤煙的品質(zhì)條件及專家意見來確定目標產(chǎn)量。然后利用土壤基礎(chǔ)肥力，結(jié)合土壤有效養(yǎng)分校正系數(shù)和肥料利用率，通過公式便可計算出不同目標產(chǎn)量下氮、磷、鉀的推薦施肥量，得到如下推薦施肥模型：

式中，Y為目標產(chǎn)量，SN、SP、SK分別為堿解氮、速效磷和速效鉀的土壤測定值。

2.4 施肥模型的分區(qū)驗證

2.4.1 推薦施肥量的計算 在4個分區(qū)內(nèi)分別取土樣測定基礎(chǔ)肥力。利用上述推薦施肥模型和基礎(chǔ)肥力計算各分區(qū)的推薦施肥量（表6）。

2.4.2 常規(guī)化學(xué)成分分析 朱尊權(quán)認為用“可用性”來評價煙葉的工業(yè)使用價值更為合適[28]，而化學(xué)成分又反映了煙葉可用性的優(yōu)劣[29]。由表7可知，推薦施肥區(qū)煙葉的總氮（除分區(qū)2外）[30]、煙堿（除分區(qū)3外）[31]、鉀（除分區(qū)3外）[4]、氯[31]和還原糖[32]及水溶性總糖[32]含量均較對照區(qū)更滿足優(yōu)質(zhì)煙葉的要求。各區(qū)域煙葉的化學(xué)成分均較協(xié)調(diào)，但推薦施肥區(qū)的糖類與煙堿、總氮與煙堿及鉀和氯之間的較對照區(qū)更加協(xié)調(diào)[30,33-34]。因此，推薦施肥區(qū)煙葉更能滿足優(yōu)質(zhì)煙葉要求，說明本試驗所得的推薦施肥模型和推薦施肥量是比較合理的。

2.4.3 致香物質(zhì)含量分析 致香物質(zhì)含量是評價煙葉內(nèi)在質(zhì)量的核心內(nèi)容和重要指標之一[35]。本研究將致香物質(zhì)簡單劃分為苯丙氨酸類、棕色化產(chǎn)物、類西柏烷降解產(chǎn)物、類胡蘿卜素降解產(chǎn)物、新植二烯5類[36]進行分析，結(jié)果見表8。

通過對推薦施肥區(qū)和對照區(qū)烤后煙葉中這幾類致香物質(zhì)含量分析發(fā)現(xiàn)，除分區(qū)3和4的推薦施肥區(qū)棕色化產(chǎn)物類含量低于對照外，其他各推薦施肥區(qū)煙葉各類香氣物質(zhì)含量及總致香物質(zhì)含量均較對照高。說明推薦施肥有利于提高烤后煙葉致香物質(zhì)含量，改善煙葉的品質(zhì)。因此，可認為本研究結(jié)果是合理的。

表6 各分區(qū)推薦施肥量Table 6 The recommended fertilizer rates for management zones

表7 烤后煙葉化學(xué)成分Table 7 The chemical component contents of cured tobacco leaves

3 討 論

研究表明，該區(qū)域土壤養(yǎng)分具有較高的異質(zhì)性，尤其是土壤AP的變異性最大，因此有必要對該區(qū)域進行分區(qū)，實施精準管理。地統(tǒng)計分析顯示，8種土壤養(yǎng)分均能用球狀模型進行較好擬合，各養(yǎng)分的空間變異性均受人為因素和結(jié)構(gòu)性因素共同作用（除 TN和 CEC外），本結(jié)果與 Wang[37]和Young[38]等報道相似?；貧w分析表明，土壤 OM、AP、AK和 CEC是該區(qū)域烤煙品質(zhì)的關(guān)鍵養(yǎng)分因子，這與劉國順[4]提出的土壤氮、磷、鉀等養(yǎng)分對煙葉品質(zhì)有著較大的貢獻一致。利用這4種關(guān)鍵養(yǎng)分因子進行模糊聚類分析可將該區(qū)劃分為4個不同肥力水平；檢驗顯示，烤煙品質(zhì)和土壤養(yǎng)分在各分區(qū)內(nèi)趨于一致，分區(qū)間具較大差異，說明用此方法來實現(xiàn)煙田養(yǎng)分的差異化管理是可行的；陳彥[39]、李翔[40]和李艷[41]等也曾利用類似方法進行了農(nóng)田養(yǎng)分分區(qū)管理研究。

表8 不同處理對煙葉香氣物質(zhì)成分含量的影響Table 8 The effect of different treatments on aroma component contents in tobacco leaves

本研究僅以中煙100品種為研究材料，獲得的養(yǎng)分需求量對于其他烤煙品種有一定局限性[42]（尤其是氮肥）。再者，由于受到土壤、降水、灌溉、施肥和氣候等因素影響，造成養(yǎng)分的揮發(fā)、淋溶、徑流和禁錮，影響了肥料利用率和土壤養(yǎng)分貢獻率。所以，本研究所得推薦施肥模型參數(shù)存在一定局限性，而且各參數(shù)年際間還有一定程度的變化。因此，本研究所獲得的推薦施肥模型需經(jīng)多年多點試驗進行不斷地修正和完善。孫克剛等[43]和潘劍君等[44]也曾運用這種方法來研究施肥模型，并得到較為理想的效果。

將推薦施肥模型運用到管理分區(qū)中顯示，運用推薦施肥模型來確定施肥量比常規(guī)施肥量更有利于形成化學(xué)成分諧調(diào)、香氣物質(zhì)含量高的優(yōu)質(zhì)烤煙。本研究的推薦施肥有利于優(yōu)質(zhì)烤煙生產(chǎn)、減少農(nóng)田肥料投入、降低農(nóng)業(yè)成本、緩解農(nóng)業(yè)對環(huán)境壓力。因此，可認為本研究的推薦施肥模型是合理的。本文僅從化學(xué)成分和致香物質(zhì)來評價推薦施肥模型是充分的，今后可用烤煙外觀質(zhì)量、物理特性、化學(xué)成分和安全性等對施肥模型進行綜合評價。

[1]Goovaerts P.Geostatistical tools for characterizing the spatial variability of microbiological and physicochemical soil properties.Biol[J].Fertil.Soils.1998,27:315-334.

[2]Iqbal J,Thomasson J A,Jenkins J N,et al.Spatial variability analysis of soil physical properties of alluvial soils[J].Soil Sci.Am.J.,2005,69:1338-1350.

[3]Ferguson R B,Hergert G W,Schepers J S,et al.Site-specific nitrogen management of irrigated maize:yield and soil residual nitrate effects[J].Soil Sci.Am.J.2002,66:544-553.

[4]劉國順.煙草栽培學(xué)[M].北京：中國農(nóng)業(yè)出版社，2003.

[5]Wang X Z,Liu G S,Hu H C,et al.Determination of management zones for a tobacco field based on soil fertility[J].Comput.Electron.Agric.2009,65:168-175.

[6]韓錦峰.煙草栽培生理[M].北京：中國農(nóng)業(yè)出版社，1986.

[7]Li Y,Shi Z,Li F.Delineation of Site-Specific Management zones based on temporal and spatial variability of soil electrical conductivity[J].Pedosphere.2007,17:156-164.

[8]Moral F J,Terro'n J M,Marques da Silva J R.Delineation of management zones using mobile measurements of soil apparent electrical conductivity and multivariate geostatistical techniques[J].Soil Till.Res.,2010,106:335-343.

[9]郭熙，賴錦春，趙小敏，等.基于 GIS 丘陵土壤分區(qū)的水稻施肥配方研究[J].中國農(nóng)業(yè)科學(xué)，2011，44（2）：307-315.

[10]Nelson D W,Sommers L E.Total carbon,organic carbon and organic matter.In:Methods of soil analysis[J].Part 2.Agron.Monogr.,1982 :539-579.

[11]Lefroy RDB,Blair G,Strong W M.Changes in soil organic matter with cropping as measured by organic carbon fractions and 13C natural isotope abundance[J].Plant Soil J.,1993:155-156,399-402.

[12]Bremner J M,Mulvaney C S.Total nitrogen.In:Methods of Soil Analysis [M].Madison,WI,Am.Soc.Argon.,1984:595-624.

[13]鮑士旦.土壤農(nóng)化分析[M].北京：農(nóng)業(yè)出版社，2005.

[14]Olsen S R,Cole C V,Watanabe F S,et al.Estimation of available phosphorus in soils by extraction with sodium bicarbonate[M].Washington,DC.USDA Circ.939,1954.

[15]Richards L A.Diagnosis and improvement of saline and alkaline soil[M].Washington,DC.USDA Hand Book.No.60.U.S.Govt.Print.Office.,1954:60.

[16]Chapman H D.Cation exchange capacity.In:methods of soil analysis[M].Madison,WI,American institute agronomy,1965:891-901.

[17]Wu D X,Yuan Z Y,Yan K Y,et al.Establishment on quality indexes of different grades of flue-cured tobacco[J].Tobacco Science & Technology,2001,12:9-15.

[18]Robertson G P.GS+geostatistics for environmental science user manual[M].Gamma Design Software.Version 3.1.Plainwell,MI,1998.

[19]Fridgen J J,Kitchen N R,Sudduth K A,et al.Management zone analyst (MZA):software for subfield management zone delineation[J].A gron J.,2004,96:100-108.

[20]Facchinelli A,Sacchi E,Mallen L.Multivariate statistical and GIS-based approach to identify heavy metal sources in soil[J].Environmental Pollution,2001,114(3):313-324.

[21]任永浩，韓錦峰.不同根際pH值下烤煙生長反應(yīng)的比較研究[J].中國煙草，1995（3）：1-5.

[22]強繼業(yè)，朱海平，周振春，等.云南省部分地區(qū)烤煙適宜 pH值范圍的緩沖研究[J].中國生態(tài)農(nóng)業(yè)學(xué)報，2005，13（2）：149-l51.

[23]全國土壤普查辦公室.中國土壤普查技術(shù)[M].北京：農(nóng)業(yè)出版社，1990.

[24]高玉蓉，許紅衛(wèi)，周斌.稻田土壤養(yǎng)分的空間變異性研究[J].土壤通報，2005，36（6）：822-825.

[25]梁中龍，甘海華，戴軍.廣州岑村農(nóng)用地土壤養(yǎng)分空間變異性研究[J].華南農(nóng)業(yè)大學(xué)學(xué)報，2004，25（3）：22-25.

[26]趙軍，孟凱，隋躍宇，等.海倫黑土有機碳和速效養(yǎng)分空間異質(zhì)性分析[J].土壤通報，2005，36（4）：488-492.

[27]Gao X B,Pei J H,Xie W X.A study of weighting exponent m in a fuzzy c-means algorithm[J].Chinese Journal of Electronics,2000,28(4):80-83.

[28]朱尊權(quán).煙葉的可用性與卷煙的安全性[J].煙草科技，2000（8）：3-6.

[29]趙銘欽，蘇長濤，姬小明，等.不同成熟度對烤后煙葉物理性狀、化學(xué)成分和中性香氣成分的影響[J].華北農(nóng)學(xué)報，2008，23（3）：146-150.

[30]胡國松，趙元寬，曹志洪，等.我國主要產(chǎn)煙省烤煙元素組成和化學(xué)品質(zhì)評價[J].中國煙草學(xué)報，1997，3（3）：36-43.

[31]王瑞新.煙草化學(xué)[M].北京：中國農(nóng)業(yè)出版社，2003.

[32]金聞博，戴亞，橫田平，等.煙草化學(xué)[M].北京：清華大學(xué)出版社，1993.

[33]章新軍，黎妍妍，許自成，等.河南烤煙外觀與內(nèi)在質(zhì)量的綜合評價[J].安徽農(nóng)業(yè)科學(xué)，2007，35（7）：1953-1954，1959.

[34]周冀衡，朱小平，王彥亭，等.煙草生理與生物化學(xué)[M].合肥：中國科學(xué)技術(shù)大學(xué)出版社，1996.

[35]Weeks W W.Chemistry of tobacco constituents influencing flavor and aroma [J].Rec.Adv.Tob.Sci.,1985(11):175-200.

[36]劉國順，江厚龍，劉淑端，等.基于烤煙品質(zhì)確定煙田的養(yǎng)分管理分區(qū)[J].植物營養(yǎng)與肥料學(xué)報，2011，17（4）：1028-1035.

[37]Wang Y D,Feng N N,Li T X,et al.Spatial variability of soil cation exchange capacity in hilly tea plantation soils under different sampling scales[J].Agricultural Sciences in China,2008,7(1):96-103.

[38]Young M H,Caldwell T G,Meadows D G,et al.Variability of soil physical and hydraulic properties at the Mojave global change facility,Nevada:implications for water budget and evapotranspiration [J].J.Arid Environ,2009,73:733-744.

[39]陳彥，呂新.基于模糊 C-均值聚類法的綠洲農(nóng)田精確管理分區(qū)研究[J].生態(tài)學(xué)報，2008，28（7）：3067-3074.

[40]李翔，潘瑜春，馬景宇，等.基于多種土壤養(yǎng)分的精準管理分區(qū)方法研究[J].土壤學(xué)報，2007，44（1）：14-20.

[41]李艷，史舟，吳次芳，等.基于多源數(shù)據(jù)的鹽堿地精確農(nóng)作管理分區(qū)研究[J].農(nóng)業(yè)工程學(xué)報，2007，23（8）：84-89.

[42]嚴小龍，張福鎖.植物營養(yǎng)遺傳[M].北京：中國農(nóng)業(yè)出版社，1997：78-106.

[43]孫克剛，楊稚娟，王守剛，等.油菜目標產(chǎn)量施肥數(shù)學(xué)模型及氮磷鉀肥推薦措施研究[J].土壤肥料，2006（1）：21-24.

[44]潘劍君，房世波，孫維俠，等.利用土壤信息系統(tǒng)進行水稻施氮量決策初探—以江蘇省大豐市為例[J].南京農(nóng)業(yè)大學(xué)學(xué)報，2000，23（3）：53-56.