程曉梅 王晉峰 孫楠 蒯雁 范志勇 張久權(quán) 徐明崗 李建華

摘? 要：探明煙田土壤速效養(yǎng)分時(shí)空變異特征，為烤煙高產(chǎn)優(yōu)質(zhì)的施肥管理提供科學(xué)依據(jù)。應(yīng)用地統(tǒng)計(jì)學(xué)等方法，分析了大理州祥云和彌渡兩個(gè)典型植煙縣1982、2012和2022年煙田土壤堿解氮、有效磷和速效鉀的時(shí)空演變特征。1982—2022年祥云土壤堿解氮含量呈先增后降趨勢(shì)，彌渡土壤堿解氮含量3個(gè)年度間無(wú)顯著差異且以III級(jí)為主（100～150 mg/kg），空間分布均勻。煙區(qū)土壤有效磷和速效鉀含量均呈不同程度上升趨勢(shì)，彌渡北部和祥云中西部高，其他區(qū)域低。近40年土壤速效養(yǎng)分塊金效應(yīng)由<25%上升至50%，表明養(yǎng)分空間變異的影響因素中地形等自然因素影響減弱，施肥等人為因素影響增強(qiáng)。土壤堿解氮含量一直處于中等水平（100～150 mg/kg），磷鉀肥施用導(dǎo)致有效磷和速效鉀含量持續(xù)升高。根據(jù)土壤速效養(yǎng)分豐缺狀況，祥云中北部、西部及彌渡北部區(qū)域煙田應(yīng)“穩(wěn)氮減磷鉀”，其余區(qū)域需注意“穩(wěn)氮控磷鉀”。

關(guān)鍵詞：大理州；煙田；速效養(yǎng)分；時(shí)空變異

中圖分類(lèi)號(hào)：S572.01? ? ? ? ? ? ? ? ? ? 文獻(xiàn)標(biāo)識(shí)碼：A? ? ? ? ? ? ? ? ? ? 文章編號(hào)：1007-5119（2023）02-0015-07

Abstract： The spatial-temporal variation characteristics of soil available nutrients were studied to provide scientific basis for fertilization management of high yield and high quality in tobacco fields. The spatio-temporal evolution characteristics of soil available nitrogen， available phosphorus and available potassium in two typical counties of Xiangyun and Midu in Dali prefecture in 1982， 2012 and 2022 were analyzed by using geostatistics. From 1982 to 2022， the average content of soil available nitrogen in Xiangyun county increased first and then decreased. There was no significant difference in soil available nitrogen between three years in Midu， which was mainly grade III （100-150 mg/kg） and evenly distributed spatially. The contents of soil available phosphorus and available potassium in tobacco-growing areas showed an upward trend to various degrees， which were high in the north of Midu and the central and western parts of Xiangyun， and low in other regions. In recent 40 years， the nugget effect of soil available nutrients increased from <25% to 50%， indicating that the influence of natural factors such as topography on the spatial variation of nutrients was weakened， and human factors such as fertilization was enhanced. The content of soil available nitrogen was at a medium level （100-150 mg/kg）， and the contents of available phosphorus and available potassium increased continuously due to the application of phosphorus and potassium fertilizers. According to the abundance and deficiency of soil available nutrients， for the tobacco fields in north central and west of Xiangyun and north of Midu we should 'stabilize nitrogen and reduce phosphorus and potassium'， and for the rest areas we need to pay attention to 'stabilize nitrogen and control phosphorus and potassium'.

Keywords： Dali Prefecture; tobacco field; soil available nutrients; spatial-temporal variability

烤煙是我國(guó)重要的經(jīng)濟(jì)作物，云南是我國(guó)最大煙葉生產(chǎn)省域[1]。該區(qū)域?yàn)閭鹘y(tǒng)農(nóng)業(yè)發(fā)展區(qū)，由于其豐富的光熱條件和充沛的水資源，生產(chǎn)的煙葉品質(zhì)獨(dú)特、氣味純凈，烤煙產(chǎn)業(yè)已然成為該區(qū)域農(nóng)業(yè)經(jīng)濟(jì)的重要支柱[2]。土壤速效養(yǎng)分是烤煙營(yíng)養(yǎng)元素的直接來(lái)源，很大程度上決定了煙葉產(chǎn)量與品質(zhì)[3]。隨著烤煙種植及煙田管理的不斷發(fā)展，煙田土壤養(yǎng)分發(fā)生了很大的變化[4-5]。因此系統(tǒng)了解典型區(qū)域煙田土壤速效養(yǎng)分時(shí)空變異特征及分布狀況有利于科學(xué)掌握土壤肥力，實(shí)現(xiàn)平衡施肥。

地統(tǒng)計(jì)學(xué)與地理信息系統(tǒng)（GIS）的結(jié)合是當(dāng)前研究土壤養(yǎng)分時(shí)空分布的主要方法[6-8]，可在不同尺度范圍內(nèi)將土壤屬性與地理數(shù)據(jù)有機(jī)結(jié)合，更加高效直觀地反映出研究區(qū)土壤養(yǎng)分變異特征和分布狀況[9]。GIS技術(shù)具有管理空間不均勻分布資源的功能，可以探究土壤養(yǎng)分時(shí)空變異，能夠直接反映區(qū)域內(nèi)煙田施肥策略及耕作方式等農(nóng)業(yè)制度的合理性，為區(qū)域煙葉生產(chǎn)中養(yǎng)分的科學(xué)管理提供參考[10]。目前，眾多研究已利用地統(tǒng)計(jì)學(xué)分析了不同區(qū)域煙田土壤養(yǎng)分的時(shí)空變異特征。鄒靜等[11]以安順市西秀區(qū)煙田土壤為研究對(duì)象，結(jié)果顯示土壤堿解氮含量略高，空間分布上無(wú)明顯規(guī)律，有效磷含量適宜，速效鉀含量較低，均呈明顯西低東高規(guī)律。譚智勇等[12]對(duì)文山州基本煙田進(jìn)行研究，該區(qū)土壤堿解氮含量整體呈由北向南、由東向西增加趨勢(shì)，土壤有效磷空間變化不顯著，土壤有效鉀含量整體上呈由中部向北部及南部增加趨勢(shì)。郭迎新等[2]研究了2011—2020年洱海流域煙田土壤養(yǎng)分時(shí)空變異特征，表明有效磷含量上升，高值區(qū)以斑塊狀分布在洱海北部、東部和西部，速效鉀含量豐富，高值區(qū)呈片狀分布在洱海北部和東部。綜合來(lái)看，我國(guó)煙田土壤速效養(yǎng)分存在分布不均衡、空間變異較大等特點(diǎn)，不同植煙區(qū)土壤養(yǎng)分含量時(shí)間變異特征及空間分布格局存在較大差異?？h域尺度的研究不僅可以具體地分析當(dāng)?shù)責(zé)熖锿寥浪傩юB(yǎng)分的時(shí)空變異特征，而且可以精準(zhǔn)結(jié)合當(dāng)?shù)貧夂驐l件、地形地貌、施肥習(xí)慣等因素為煙草種植提供科學(xué)施肥指導(dǎo)，對(duì)當(dāng)?shù)剞r(nóng)業(yè)生產(chǎn)具有重要指導(dǎo)意義[13-14]。有學(xué)者對(duì)重慶石柱縣、云南羅平縣、湘西花垣縣等地?zé)熖锿寥浪傩юB(yǎng)分空間變異進(jìn)行研究并提出區(qū)劃施肥，為縣域煙田因土施肥提供現(xiàn)實(shí)參考[15-17]。然而，目前已有的關(guān)于縣域尺度煙田土壤速效養(yǎng)分的研究大多集中于短時(shí)期內(nèi)的時(shí)空變化，關(guān)于時(shí)空分布長(zhǎng)期動(dòng)態(tài)變化的定量研究較少，缺乏對(duì)云南省典型植煙區(qū)合理化施肥策略的精確指導(dǎo)。

祥云和彌渡縣地處大理州紅花大金元特色煙葉“黃金走廊”上的核心區(qū)域，生產(chǎn)的烤煙以其獨(dú)特的“清、香、甜、潤(rùn)”而聞名全國(guó)，是許多知名品牌卷煙的首選原料，該區(qū)域植煙面積較大，具有典型性。本文探討1982—2022年大理州祥云和彌渡縣2個(gè)典型縣域煙田土壤堿解氮、有效磷和速效鉀時(shí)空變異特征及分布格局，為大理州乃至云南和全國(guó)煙田養(yǎng)分管理及平衡施肥提供參考。

1? 材料與方法

1.1? 研究區(qū)概況

大理白族自治州地處云南省中部偏西（東經(jīng)98°52′～101°03′，北緯24°41′～26°42′），平均海拔2090 m，年均溫15.8 ℃，年均日照時(shí)數(shù)2332 h，年均降水量835.7 mm。煙田土壤主要以水稻土、紅壤、紫色土、潮土和黃壤為主。全州常年烤煙種植面積約3.23萬(wàn)hm2，其中祥云縣和彌渡縣煙田面積大，分別為5600 hm2和3100 hm2。

1.2? 樣品采集與分析

根據(jù)祥云和彌渡縣烤煙種植情況，于2012年和2022年分別采集289個(gè)（173和116個(gè)）和1126個(gè)（700和426個(gè)）耕層（0～20 cm）土壤樣品（圖1）。隨著時(shí)間變化和技術(shù)發(fā)展，土壤樣點(diǎn)采集更加密集。土樣采集時(shí)間選在烤煙前茬作物收獲后未進(jìn)行翻耕和施肥時(shí)，以反映取樣點(diǎn)的真實(shí)理化性狀，同時(shí)使用GPS定位，記錄經(jīng)緯度等信息。取樣時(shí)，每取樣點(diǎn)以對(duì)角線法采20個(gè)點(diǎn)，土壤均勻混合后四分法留取1 kg，作為1個(gè)樣品備用。土樣經(jīng)自然風(fēng)干、研磨、過(guò)篩后測(cè)定土壤速效養(yǎng)分，堿解氮采用堿解擴(kuò)散法測(cè)定；有效磷采用鉬銻抗比色法測(cè)定；速效鉀采用原子吸收分光光度法測(cè)定。1982年的數(shù)據(jù)來(lái)源于祥云和彌渡縣第二次土壤普查資料，共51個(gè)（24和27個(gè)）煙田耕層土壤數(shù)據(jù)。

1.3? 土壤速效養(yǎng)分分級(jí)

依據(jù)《大理州植煙土壤分析與評(píng)價(jià)》土壤養(yǎng)分分級(jí)標(biāo)準(zhǔn)[18]，對(duì)煙田土壤速效養(yǎng)分含量進(jìn)行分級(jí)，共分為5個(gè)等級(jí)，分別為V（低）、IV（較低）、III（中等）、II（較高）、I（高）（表1）。

1.4? 數(shù)據(jù)處理

利用SPSS 17.0進(jìn)行統(tǒng)計(jì)描述和正態(tài)檢驗(yàn)，Minitab 20進(jìn)行Box-Cox變換，使之服從正態(tài)分布。

采用ArcGIS 10.6地統(tǒng)計(jì)學(xué)軟件中的Geostatistical Analyst模塊的普通Kriging插值法繪制土壤速效養(yǎng)分空間分布圖，運(yùn)用重分類(lèi)（Reclassify）計(jì)算每個(gè)等級(jí)面積比例。半方差函數(shù)模型計(jì)算和最優(yōu)理論模型擬合使用GS+9.0完成，計(jì)算塊金值、基臺(tái)值、塊金效應(yīng)、變程、殘差、決定系數(shù)等參數(shù)，模型選取以決定系數(shù)越大、殘差越小為最佳。塊金效應(yīng)即塊金值與基臺(tái)值之比，表示空間異質(zhì)性程度或土壤屬性的空間依賴(lài)性。當(dāng)比值<25%時(shí)，表示變量具有強(qiáng)烈的空間自相關(guān)性，空間變異主要由結(jié)構(gòu)因素引起；當(dāng)比值在25%～75%范圍內(nèi)，表示變量具有中等水平的空間自相關(guān)性，空間變異由結(jié)構(gòu)因素和隨機(jī)因素共同作用；當(dāng)比值>75%時(shí)，表示變量具有較弱的空間自相關(guān)性，變異主要由隨機(jī)因素引起。結(jié)構(gòu)因素包括地形、氣候、土壤母質(zhì)等自然因素，隨機(jī)因素主要包括施肥、耕作、種植制度和灌溉等人為活動(dòng)。

2? 結(jié)? 果

2.1? 土壤速效養(yǎng)分含量統(tǒng)計(jì)特征

表2示出，40年間祥云植煙土壤堿解氮含量呈先增后降趨勢(shì)，1982—2012年年均增速1.1 mg/kg，2012—2022年年均下降速率為1.2 mg/kg，至2022年降為119.1 mg/kg；空間變異系數(shù)逐漸降低。有效磷在1982—2012年間由V級(jí)上升為III級(jí)，年均增速為0.5 mg/kg，2012—2022年年均增速有所提高（1.6 mg/kg）；3次取樣有效磷空間變異系數(shù)均屬中等，其中1982年變異性最強(qiáng)。速效鉀含量在1982—2022年間呈上升趨勢(shì)，增加了1.3倍；均屬中等變異。

彌渡土壤堿解氮含量在3次取樣間（115.0～129.1 mg/kg）變化不顯著，均處于III級(jí)。有效磷含量在1982—2022年呈不斷增加趨勢(shì)，1982—2012年年均增加1.1 mg/kg，2012—2022年年均增加1.7 mg/kg。速效鉀含量1982—2022年呈增加趨勢(shì)，1982—2012年提升較緩，年均增加3.0 mg/kg，2012—2022年提升較快，年均增加15 mg/kg；其中1982年空間分布呈強(qiáng)變異，2012和2022年均呈中等變異。

2.2? 土壤速效養(yǎng)分半方差函數(shù)分析

表3為3次取樣煙田土壤速效養(yǎng)分最優(yōu)半方差函數(shù)模型及參數(shù)。1982年土壤堿解氮塊金效應(yīng)為15.0%，該時(shí)期堿解氮主要受結(jié)構(gòu)因素影響，2012和2022年塊金效應(yīng)為40.9%和50.0%，受結(jié)構(gòu)因素和隨機(jī)因素共同影響。1982和2012年有效磷與速效鉀塊金效應(yīng)均小于25%，具有較強(qiáng)的空間自相關(guān)性；2022年土壤有效磷與速效鉀塊金效應(yīng)均為49.9%，屬中等程度空間自相關(guān)?？梢?jiàn)近40年土壤有效磷和速效鉀結(jié)構(gòu)性減弱，隨機(jī)性增強(qiáng)，烤煙種植過(guò)程中施肥等農(nóng)藝措施影響增強(qiáng)。

2.3? 土壤速效養(yǎng)分空間分布格局

由圖2可見(jiàn)，1982—2022年祥云縣土壤堿解氮含量以III級(jí)為主，呈大片狀分布于全區(qū)，下莊鎮(zhèn)、普棚鎮(zhèn)及鹿鳴鄉(xiāng)部分區(qū)域堿解氮含量隨時(shí)間推移上升趨勢(shì)明顯。彌渡縣土壤堿解氮含量1982年以III級(jí)為主（58.6%）；2012年和2022年III級(jí)面積進(jìn)一步擴(kuò)大，占比為75%以上。

由圖3可見(jiàn)，祥云縣土壤有效磷含量1982年以IV（45.8%）和V級(jí)（45.6%）為主；2012年III級(jí)以上等級(jí)占比增加至85.1%，空間分布相對(duì)均勻；2022年III級(jí)及以上等級(jí)占比增加至94.7%，其中祥城鎮(zhèn)、下莊鎮(zhèn)、禾甸鎮(zhèn)和普棚鎮(zhèn)大部分區(qū)域I級(jí)和II級(jí)面積明顯增加。彌渡有效磷含量1982年以V級(jí)為主，空間分布均勻；2012年以III級(jí)（69.7%）為主，II級(jí)次之（19.4%）；2022年III級(jí)及以上等級(jí)面積占比增加至92.1%，其中彌渡北部部分區(qū)域有效磷含量超過(guò)80 mg/kg。

由圖4可見(jiàn)，祥云速效鉀含量1982年以IV級(jí)（73.6%）為主；2012年III級(jí)及以上等級(jí)占比為57.8%；2022年IV等級(jí)面積明顯下降，III級(jí)及以上等級(jí)占比增加至94.1%，其中I和II等級(jí)占比62.8%，主要分布在祥城鎮(zhèn)、禾甸鎮(zhèn)、下莊鎮(zhèn)和云南驛鎮(zhèn)等大部分區(qū)域。彌渡速效鉀含量1982年以IV級(jí)為主（37.1%），V級(jí)次之（33.8%），III級(jí)及以上等級(jí)主要分布于米甸鎮(zhèn)和牛街彝族鄉(xiāng)等區(qū)域；2012年速效鉀含量大幅提升，III級(jí)面積占比42.7%；2022年I級(jí)和II級(jí)面積占比高達(dá)85.0%，呈片狀分布于全區(qū)，III級(jí)及以下等級(jí)面積占比僅為15.0%，零星分布于新街鎮(zhèn)及牛街鄉(xiāng)等部分區(qū)域。

總體來(lái)看，大理兩個(gè)縣域煙田堿解氮、有效磷和速效鉀含量的時(shí)空變異明顯。40年間彌渡北部及祥云中部堿解氮含量上升，空間分布更加均勻。有效磷和速效鉀的空間分布趨勢(shì)相似，呈彌渡北部和祥云中西部高其他區(qū)域低的分布格局。

3? 討? 論

本研究表明，1982—2022年間，祥云和彌渡縣煙田土壤堿解氮含量整體處于III級(jí)，屬中等范圍。祥云縣堿解氮呈先升后降趨勢(shì)，這是由于20世紀(jì)90年代化學(xué)氮肥用量逐年增多[19]所致。2012年之后，隨著煙草的規(guī)模化種植，為保證煙草產(chǎn)量和質(zhì)量，煙草公司制定了“控氮、穩(wěn)磷、增鉀”原則，因此堿解氮含量又有所下降，但氮素含量整體相對(duì)適宜。而彌渡縣3次取樣堿解氮含量變化不顯著（115.0～129.1 mg/kg），這可能與該縣植煙區(qū)煙農(nóng)的施肥習(xí)慣及有機(jī)肥的合理配施等有關(guān)。據(jù)調(diào)查，祥云縣煙區(qū)農(nóng)戶(hù)習(xí)慣施氮量（90 kg/hm2）高于彌渡（60 kg/hm2）。彌渡北部及祥云中部等低氮區(qū)域土壤堿解氮含量逐漸增加至中等水平，中等肥力區(qū)域保持相對(duì)穩(wěn)定，這與大理重視煙草產(chǎn)業(yè)發(fā)展，合理的施肥指導(dǎo)政策緊密相關(guān)。1982年堿解氮以強(qiáng)自相關(guān)為主，2012年以來(lái)自相關(guān)減弱，說(shuō)明烤煙種植過(guò)程中受人為管理影響增強(qiáng)。近年來(lái)煙田推行使用含氮復(fù)合肥和有機(jī)肥以及煙田管理水平的提高，土壤肥力變得更加均勻。針對(duì)祥云和彌渡煙田現(xiàn)狀，大部分區(qū)域應(yīng)以“穩(wěn)氮”為主，祥云米甸鎮(zhèn)及彌渡牛街鄉(xiāng)部分低氮區(qū)應(yīng)適當(dāng)增施氮肥，確保煙葉的正常生長(zhǎng)發(fā)育。

1982—2022年大理州2個(gè)典型縣域煙田土壤有效磷含量均呈上升趨勢(shì)。由于第二次土壤普查發(fā)現(xiàn)大面積農(nóng)田土壤磷含量較低，增施磷肥對(duì)作物增產(chǎn)明顯，農(nóng)民加大了磷肥施用量[20]，而烤煙對(duì)磷肥的當(dāng)季利用率僅為15%～18%，導(dǎo)致土壤中磷素大量累積，煙田磷含量不斷提升。祥云中西部和彌渡北部有效磷含量高于其他區(qū)域，因祥云中西部和彌渡北部地勢(shì)較平緩，多為壩地，有相對(duì)優(yōu)越的灌溉設(shè)施及管理?xiàng)l件，植煙面積大，促使煙田肥料的投入量高于其他區(qū)域，而且地勢(shì)平緩，磷不易隨水分流失。土壤有效磷水平過(guò)高，會(huì)導(dǎo)致煙葉葉片厚而粗糙，油分少?gòu)椥孕?，粗筋暴梗，并引發(fā)環(huán)境污染[21-22]。針對(duì)大理州煙田土壤現(xiàn)狀，應(yīng)以“減磷”為主，尤其對(duì)祥云中北部、西部及彌渡北部高磷區(qū)應(yīng)實(shí)施減磷措施，其他中磷區(qū)應(yīng)控制磷肥施用，而彌渡牛街鄉(xiāng)南部區(qū)域則可適當(dāng)增施磷肥。

大理2個(gè)典型縣域煙田土壤速效鉀含量均呈顯著上升趨勢(shì)。煙葉鉀含量是衡量煙葉品質(zhì)的一個(gè)重要指標(biāo)，能提高烤煙的燃燒性，減少煙葉中有害物質(zhì)如焦油等含量，提升煙葉香氣、吃味、色澤等[23]。土壤鉀素的供應(yīng)狀況直接影響著煙葉含鉀量[24]。20世紀(jì)80年代，隨著社會(huì)進(jìn)步和化肥種類(lèi)的完善，烤煙生產(chǎn)注重多養(yǎng)分協(xié)調(diào)，硫酸鉀、硝酸鉀等鉀肥施用量逐年增加[19]。為追求高品質(zhì)高效益，大理煙草公司制定“增鉀”原則，部分煙區(qū)存在超標(biāo)準(zhǔn)施鉀現(xiàn)象[2，25]，導(dǎo)致煙區(qū)土壤速效鉀含量不斷增加。鉀肥長(zhǎng)期大量施用導(dǎo)致當(dāng)季利用率下降[26]，土壤鉀素積累增加。研究區(qū)土壤速效鉀的空間分布表現(xiàn)為彌渡北部和祥云中西部高其他區(qū)域低的分布格局。祥云彌渡為云南老牌煙區(qū)，植煙年限較長(zhǎng)，各區(qū)域農(nóng)戶(hù)在農(nóng)業(yè)生產(chǎn)過(guò)程中采用不同的施肥習(xí)慣，且祥云中西部和彌渡北部多為壩地，有機(jī)肥推廣施用較易，投入較多，且該區(qū)域地勢(shì)平緩，鉀流失少，土壤殘留高，因此鉀素含量高。針對(duì)鉀含量較高的祥云中西部、北部和彌渡北部、中南部等區(qū)域宜實(shí)施“減鉀”措施，確?？緹熒a(chǎn)投入減量不減收，其他區(qū)域宜實(shí)施“穩(wěn)鉀”措施。

總體來(lái)看，隨機(jī)性因素如施肥等成為該區(qū)當(dāng)前土壤速效氮磷鉀的主要影響因素。從當(dāng)前煙田土壤速效養(yǎng)分來(lái)看，氮含量相對(duì)適宜，部分區(qū)域磷鉀豐富，持續(xù)采取“控氮提鉀”措施會(huì)造成土壤養(yǎng)分失衡[26]，應(yīng)根據(jù)土壤氮、磷、鉀養(yǎng)分狀況，及時(shí)調(diào)整施肥策略。

4? 結(jié)? 論

1982—2022年40年間，祥云土壤堿解氮含量呈現(xiàn)先增后降趨勢(shì)，彌渡3次采樣土壤堿解氮含量無(wú)顯著差異且均以III級(jí)為主（100～150 mg/kg）；兩縣土壤有效磷和速效鉀含量均表現(xiàn)為不同程度上升趨勢(shì)，有效磷含量高于20 mg/kg 和速效鉀含量高于150 mg/kg的面積增加至90%以上，空間上呈彌渡北部和祥云中西部高其他區(qū)域低的分布特點(diǎn)。當(dāng)前煙田土壤速效養(yǎng)分氮含量相對(duì)適宜，部分區(qū)域磷鉀豐富，祥云中北部、西部及彌渡北部區(qū)域應(yīng)“穩(wěn)氮減磷鉀”，其他區(qū)域需注意“穩(wěn)氮控磷鉀”。

參考文獻(xiàn)

[1]孫計(jì)平，吳照輝，李雪君，等. 21世紀(jì)中國(guó)烤煙種植區(qū)域及主栽品種變化分析[J]. 中國(guó)煙草科學(xué)，2016，37（3）：86-92.

SUN J P， WU Z H， LI X J， et al. Analysis of regional variation and major varieties of flue-cured tobacco planted in China in the Twenty-first Century[J]. Chinese Tobacco Science， 2016， 37（3）： 86-92.

[2]郭迎新，陳永亮，苗琪，等. 洱海流域植煙土壤養(yǎng)分時(shí)空變異特征及肥力評(píng)價(jià)[J]. 中國(guó)農(nóng)業(yè)科學(xué)，2022，55（10）：1987-1999.

GUO Y X， CHEN Y L， MIAO Q， et al. Spatial-temporal variability of soil nutrients and assessment of soil fertility in Erhai Lake Basin[J]. Scientia Agricultura Sinica， 2022， 55（10）： 1987-1999.

[3]翟?xún)?yōu)雅，阮志，張立新，等. 施用秸稈對(duì)烤煙土壤酶活性及烤煙氮和鉀含量的影響[J]. 中國(guó)煙草科學(xué)，2014，35（5）：45-49.

ZHAI Y Y， RUAN Z， ZHANG L X， et al. Effects of straw application on soil enzyme activities and nitrogen and potassium contents in flue-cured tobacco[J]. Chinese Tobacco Science， 2014， 35（5）： 45-49.

[4]李倩，索炎炎，景延秋，等. 南陽(yáng)植煙土壤養(yǎng)分時(shí)空演變特征及肥力評(píng)價(jià)[J]. 中國(guó)煙草科學(xué)，2022，43（5）：22-30.

LI Q， SUO Y Y， JING Y Q， et al. Spatial-temporal variability of soil nutrients and assessment of soil fertility in Nanyang Tobacco planting areas[J]. Chinese Tobacco Science， 2022， 43（5）： 22-30.

[5]張恒，黃鶯，劉明宏，等. 基于空間插值法的遵義煙區(qū)植煙土壤養(yǎng)分時(shí)空變化[J]. 中國(guó)煙草科學(xué)，2020，41（3）：36-43.

ZHANG H， HUANG Y， LIU M H， et al. Temporal and spatial variations of soil nutrients in Zunyi tobacco growing area based on spatial interpolation method[J]. Chinese Tobacco Science， 2020， 41（3）： 36-43.

[6]YAN Y B， YANG Y. Uncertainty assessment of spatiotemporal distribution and variation in regional soil heavy metals based on spatiotemporal sequential Gaussian simulation[J]. Environmental pollution， 2023， 322： 121243.

[7]SUN T， TONG W J， CHANG N J， et al. Estimation of soil organic carbon stock and its controlling factors in cropland of Yunnan Province， China[J]. Journal of Integrative Agriculture， 2022， 21（5）： 1475-1487.

[8]SUN X L， ZHAO Y G， WU Y J， et al. Spatio-temporal change of soil organic matter content of Jiangsu Province， China， based on digital soil maps[J]. Soil Use and Management， 2012， 28（3）： 318-328.

[9]ZHANG Q， YANG Z P， LI Y， et al. Spatial variability of soil nutrients and GIS-based nutrient management in Yongji County， China[J]. International Journal of Geographical Information Science， 2010， 24（7）： 965-981.

[10]張隆偉，伍仁軍，王昌全，等. 四川涼攀煙區(qū)植煙土壤有效銅和有效鋅空間變異特征[J]. 中國(guó)煙草科學(xué)，2014，35（3）：1-6.

ZHANG L W， WU R J， WANG C Q， et al. Spatial variability characteristics of available Cu and Zn in tobacco planting soils in Liangshan and Panzhihua， Sichuan Province[J]. Chinese Tobacco Science， 2014， 35（3）： 1-6.

[11]鄒靜，李熙全，代飛，等. 喀斯特山區(qū)植煙土壤肥力綜合評(píng)價(jià)及空間分布特征——以安順西秀區(qū)為例[J]. 山地農(nóng)業(yè)生物學(xué)報(bào)，2021，40（2）：29-35.

ZOU J， LI X Q， DAI F， et al. Comprehensive evaluation of tobacco-growing soil fertility in Karst Mountain Are-take Xixiu district of Anshun as an example[J]. Journal of Mountain Agriculture and Biology， 2021， 40（2）： 29-35.

[12]譚智勇，周冀衡，李強(qiáng)，等. 文山州土壤養(yǎng)分空間分布與變化規(guī)律[J]. 礦物學(xué)報(bào)，2016，36（4）：600-604.

TAN Z Y， ZHOU J H， LI Q， et al. Spatial distribution and variation characteristics of soil nutrients in Wenshan cultivated soil[J]. Acta Mineralogica Sinica， 2016， 36（4）： 600-604.

[13]張隆偉，伍仁軍，王昌全，等. 四川涼攀煙區(qū)植煙土壤有效銅和有效鋅空間變異特征[J]. 中國(guó)煙草科學(xué)，2014，35（3）：1-6.

ZHANG L W， WU R J， WANG C Q， et al. Spatial variability characteristics of available Cu and Zn in tobacco planting soils in Liangshan and Panzhihua， Sichuan Province[J]. Chinese Tobacco Science， 2014， 35（3）： 1-6.

[14]郭婷，李宏光，李文，等. 湘南稻作煙區(qū)土壤有效硼的時(shí)空變異及影響因素分析[J]. 煙草科技，2019，52（7）：27-34.

GUO T， LI H G， LI W， et al. Spatial-temporal variability and factors to influence boron availability in tobacco-planting soils of rotating paddy-upland fields in southern Hunan[J]. Tobacco Science & Technology， 2019， 52（7）： 27-34.

[15]陳潔菲，查宇璇，楊超，等. 重慶市石柱縣煙田土壤肥力演變與施肥區(qū)劃[J]. 土壤，2021，53（6）：1207-1214.

CHEN J F， ZHA Y X， YANG C， et al. Evolution and fertilization zoning of tobacco-growing soil fertility of Shizhu County， Chongqing City[J]. Soils， 2021， 53（6）： 1207-1214.

[16]王育軍，周冀衡，孫書(shū)斌，等. 云南省羅平縣煙區(qū)土壤肥力適宜性評(píng)價(jià)及養(yǎng)分時(shí)空變異特征[J]. 土壤，2015，47（3）：515-523.

WANG Y J， ZHOU J H， SUN S B， et al. Evaluation of soil fertility suitability and spatial/temporal variability of nutrient contents of tobacco-planting soils in Luoping County of Yunnan， China[J]. Soils， 2015， 47（3）： 515-523.

[17]胡瑞文，黎娟，向德明，等. 湘西州花垣煙區(qū)植煙土壤肥力適宜性評(píng)價(jià)及影響因素[J]. 云南農(nóng)業(yè)大學(xué)學(xué)報(bào)（自然科學(xué)），2018，33（3）：520-528.

HU R W， LI J， XIANG D M， et al. Suitability evaluation and influencing factors of soil fertility in tobacco growing areas of Huayuan County， Xiangxi Prefecture[J]. Journal of Yunnan Agricultural University（Natural Science）， 2018， 33（3）： 520-528.

[18]李再光. 大理州植煙土壤分析與評(píng)價(jià)[M]. 昆明：云南人民出版社，2014.

LI Z G. Analysis and evaluation of tobacco planting soil in Dali Prefecture[M]. Kunming： Yunnan People's Press， 2014.

[19]云南省大理州煙草專(zhuān)賣(mài)局，云南省煙草大理州公司. 大理州煙草志[M]. 昆明：云南民族出版社，2005.

Dali Prefecture Tobacco Monopoly Bureau of Yunnan Province， Yunnan Tobacco Dali Prefecture Company. Tobacco of Dali Prefecture[M]. Kunming： Yunnan Ethnic Publishing House， 2005.

[20]陳朝陽(yáng)，吳平，陳星峰，等. 南平市植煙土壤氮、磷、鉀養(yǎng)分狀況與演變趨勢(shì)[J]. 中國(guó)農(nóng)學(xué)通報(bào)，2011，27（25）：68-76.

CHEN C Y， WU P， CHEN X F， et al. Dynamic variations of nitrogen， phosphorus and potassium in Nanping tobacco-growing soil[J]. Chinese Agricultural Science Bulletin， 2011， 27（25）： 68-76.

[21]楊智棋. 攀西地區(qū)植煙土壤有效磷空間變異特征及其對(duì)煙葉品質(zhì)的影響研究[D]. 雅安：四川農(nóng)業(yè)大學(xué)，2014.

YANG Z Q. Panxi tobacco planting area soil spatial variability of available P and its effects on the quality of tobacco[D]. Yaan Sichuan Agricultural University， 2014.

[22]劉建國(guó)，羅建新，楊磊，等. 湖南植煙土壤磷素狀況分析[J]. 河南農(nóng)業(yè)大學(xué)學(xué)報(bào)，2020，54（1）：125-132.

LIU J G， LUO J X， YANG L， et al. Analysis of phosphorus status in tobacco planting soils in Hunan[J]. Journal of Henan Agricultural University， 2020， 54（1）： 125-132.

[23]邢云霞，王蒙，劉世亮，等. 氮鉀肥施用方法對(duì)烤煙品質(zhì)及經(jīng)濟(jì)效益的影響[J]. 中國(guó)水土保持科學(xué)，2014，12（6）：90-96.

XING Y X， WANG M， LIU S L， et al. Effect of N and K fertilizer application methods on quality and economic benefits of Nicotiana tabacum L[J]. Science of Soil and Water Conservation， 2014， 12（6）： 90-96.

[24]季璇，陳熙卓，宋文靜，等. 施鉀量對(duì)不同生態(tài)區(qū)烤煙生長(zhǎng)和鉀肥利用率的影響[J]. 中國(guó)煙草科學(xué)，2020，41（5）：36-42.

JI X， CHEN X Z， SONG W J， et al. Effects of potassium application rate on growth and potassium use efficiency of flue-cured tobacco in different ecological regions[J]. Chinese Tobacco Science， 2020， 41（5）： 36-42.

[25]賀彩平，劉峰，石孝均，等. 重慶市巫山縣煙草施肥現(xiàn)狀調(diào)查分析[J]. 中國(guó)農(nóng)學(xué)通報(bào)，2013，29（7）：179-184.

HE C P， LIU F， SHI X J， et al. Investigation and analysis of fertilizer application on tobacco in Wushan County of Chongqing[J]. Chinese Agricultural Science Bulletin， 2013， 29（7）： 179-184.

[26]肖怡，馮浪，王鵬，等. 鹽源縣植煙土壤速效養(yǎng)分空間變異特征研究[J]. 浙江農(nóng)業(yè)學(xué)報(bào)，2017，29（8）：1329-1335.

XIAO Y， FENG L， WANG P， et al. Study on spatial variability of soil available nutrients in tobacco planting soil in Yanyuan County[J]. Acta Agriculturae Zhejiangensis， 2017， 29（8）： 1329-1335.