趙明珠 吳婷 唐瑾 馬關(guān)潤 郭鐵英 蕭自位 蘇琳琳 周華 白學慧

摘? 要? 全國第二次土壤普查迄今30年，小?？Х鹊耐寥婪柿Σ粩喔伦兓?，為實現(xiàn)土地資源的有效利用并促進咖啡產(chǎn)業(yè)的可持續(xù)發(fā)展，探尋小粒咖啡的土壤肥力現(xiàn)狀和變化特征尤為重要。2015—2018年間對小?？Х?個主產(chǎn)區(qū)采集的大量土壤數(shù)據(jù)進行分析，研究咖啡地土壤肥力現(xiàn)狀，并與第二次土壤普查數(shù)據(jù)相對比，得出咖啡土壤肥力的變化特征及原因。目前，小?？Х韧寥纏H、有機質(zhì)、堿解氮、有效磷和速效鉀的平均值分別為5.53、32.68 g/kg、128.77 mg/kg、18.01 mg/kg和138.31 mg/kg，主要分布范圍分別為4.5～5.5、>40 g/kg、120～150 mg/kg、10～20 mg/kg和100～150 mg/kg。小粒咖啡土壤肥力在不同海拔高度的分布均存在顯著差異，pH和堿解氮含量隨海拔高度增加而降低，有機質(zhì)和有效磷含量隨海拔高度增加而增加，速效鉀含量隨海拔高度的增加呈先增加后降低。與第二次土壤普查相比，目前小?？Х韧寥莱仕峄厔荩袡C質(zhì)、堿解氮和速效鉀含量明顯降低，而有效磷含量顯著提高。宜耕區(qū)長期頻繁耕作、大量施用化學肥料，小?？Х韧寥婪柿ψ兓@著。建議小?？Х戎鳟a(chǎn)區(qū)應控制磷肥施用量，推進咖啡秸稈、果皮發(fā)酵還田的力度，加快小?？Х韧寥婪柿Φ奶岣卟⑦_到養(yǎng)分供需平衡。

關(guān)鍵詞? 小?？Х?土壤肥力;海拔;全國第二次土壤普查中圖分類號? S571.2? ? ? 文獻標識碼? A

Abstract? The soil fertility of Coffea arabica L. was analyzed by massive data of soil sampled in six main producing regions in Yunnan Province during 2015-2018， and compared with the data obtained during the second national soil survey to explore the laws and causes of the variation of the soil fertility. Presently， the soil fertility of C. arabica L. in Yunnan was pH 4.5-5.5 or 5.53 on average， organic matter 32.68 g/kg， alkalystic N 128.77 mg/kg， available P 18.01 mg/kg， and available K 138.31 mg/kg. The soil fertility had significant difference in different altitudes. pH and the content of alkalystic N decreased， the content of soil organic matter and available P increased， and the content of available K increased first and then decreased with the increase of altitude. Compared with the results from the second national soil survey， the soil pH， the content of organic matter， alkalystic N and available K decreased while available P increased significantly. The soil fertility varied significantly due to long-term frequent tillage and large application of chemical fertilizers in suitable areas. It is suggested that the application amount of phosphate fertilizer should be controlled， and the returning of coffee straw and peel to the field should be promoted to achieve the balance of supply and demand for nutrient.

Keywords? C. arabica L.; soil fertility; altitude; the second national soil survey

DOI? 10.3969/j.issn.1000-2561.2019.04.002

土壤肥力狀況制約著每個農(nóng)業(yè)作物的安全可持續(xù)發(fā)展，通過利用土壤肥力的變化特征可為土壤資源和合理施肥提供有效的科學依據(jù)[1]。作為世界三大飲料作物的咖啡，在中國的新興特色農(nóng)業(yè)中，扮演著舉足輕重的地位。隨著時代的不斷演替與更新，土壤肥力隨著土地的利用方式、施肥水平和管理措施不斷發(fā)生變化，不同的海拔、土壤類型、氣候和生物等因素，致使同一地區(qū)不同的區(qū)域土壤肥力狀況存在較大差異。隨著我國勞動力不斷向城市轉(zhuǎn)移，農(nóng)村的剩余勞動力較為稀缺，咖啡的種植、采摘過程以及后期的加工處理都亟需勞動力。

我國從20世紀30年代就開始在全國進行第二次土壤普查，但與上次土壤普查已有30年之隔，目前海拔與土壤肥力的相關(guān)性研究已有諸多報道，取得豐富成果，地形、植被結(jié)構(gòu)、溫度、水分和光照等環(huán)境因子隨著海拔的上升表現(xiàn)出一定的規(guī)律性，但由于研究的作物、區(qū)域差異，所得出的結(jié)論也不盡相同。焦?jié)櫚驳萚2]對隴南白龍江流域油橄欖園土壤肥力研究指出土壤肥力總體呈現(xiàn)了隨海拔升高而降低的趨勢，華北山區(qū)海拔越低土壤有效磷含量越高[3]，武夷山、廬山茶園和納板河自然保護區(qū)土壤有機質(zhì)含量隨海拔高度升高而升高[4-5]，薛沛沛等[6]對大崗山不同海拔毛竹林土壤肥力研究也指出海拔與土壤肥力呈正相關(guān)關(guān)系。云南省的普洱、臨滄、德宏、保山、文山和怒江作為小?？Х鹊闹鳟a(chǎn)區(qū)，每個產(chǎn)區(qū)的海拔不一，隨之的土壤肥力也各不相同。隨著近年來我國對咖啡行業(yè)的不斷重視，小?？Х鹊姆N植面積和產(chǎn)量不斷上升，化肥施用量也隨之升高，有機肥攝入量卻不斷萎縮，導致土壤理化性質(zhì)不斷發(fā)生改變，土壤質(zhì)量降低。

本文通過對云南咖啡地土壤樣品的測定、數(shù)據(jù)分析統(tǒng)計，明確小?？Х韧寥婪柿ΜF(xiàn)狀，并與全國第二次土壤普查的結(jié)果比較，得出土壤肥力的變化特征，為咖啡地合理施肥、促進土壤養(yǎng)分平衡提供一定的理論依據(jù)。

1? 材料與方法

1.1? 研究區(qū)概況

云南省西南部的普洱、臨滄、怒江，西部的德宏、保山以及東南部的文山都屬于云南的熱區(qū)，同時是生產(chǎn)小?？Х鹊狞S金地帶，擁有著得天獨厚的自然條件優(yōu)勢：位于北緯21～25°的亞熱帶氣候，年均氣溫15～25 ℃，海拔500～2000 m（高海拔促進咖啡優(yōu)良品質(zhì)的形成），年平均降雨量豐沛700～1800 m，雨熱同期，日照充足（促進咖啡生長），晝夜溫差大（15～20 ℃，利于咖啡養(yǎng)分積累），干濕分明（利于咖啡采收），基本全年無霜（咖啡為喜熱作物，不耐寒），成為孕育優(yōu)質(zhì)咖啡的溫床。普洱、臨滄、德宏和保山并稱云南小粒咖啡的四大核心主產(chǎn)區(qū)。截至目前，云南全省的小?？Х确N植面積達11.8萬hm2，年產(chǎn)量13.9萬t。根據(jù)小?？Х戎鳟a(chǎn)區(qū)每個采樣點的海拔高度，劃分為3個海拔梯度，分別為低海拔（500～1000 m）地區(qū)、中海拔（1000～1500 m）地區(qū)和高海拔（>1500 m）地區(qū)。

1.2? 土樣采集與分析

2015—2018年間在云南省小粒咖啡主產(chǎn)區(qū)：普洱、臨滄、德宏、保山、文山和怒江地區(qū)，采用GPS定位。每個點采集3份土樣（0～20 cm），采用“梅花形”采樣法采集[7]，共采集684份土壤樣品，土樣分布詳細情況如表1所示。

采用2.5∶1水土比，pH計測定土壤pH;重鉻酸鉀容量法-外加熱法測定土壤有機質(zhì)（OM）;1.0 mol/L NaOH堿解擴散法測定土壤堿解氮（AN）;0.03 mol/L NH4F-0.025 mol/L HCl浸提-比色法測定土壤有效磷（AP）;1.0 mol/L NH4OAc浸提-火焰光度法測定土壤速效鉀（AK）[8]。

1.3? 數(shù)據(jù)分析

采用Excel軟件和SPSS軟件進行數(shù)據(jù)整理與分析。土壤養(yǎng)分分級按照全國第二次土壤普查[9]時指定的標準，利于二者間的分析比較。

2? 結(jié)果與分析

2.1? 土壤pH分析

土壤pH是評價土壤質(zhì)量高低的重要指標之一，其大小制約土壤養(yǎng)分的含量及其有效性，從而影響土壤肥力的形成[10]。目前云南省小粒咖啡主產(chǎn)區(qū)的土壤pH范圍4.21～8.11（表2），呈正態(tài)分布（圖1A），平均值為5.53，變異系數(shù)為13.75%，

較低。土壤強酸、酸性、微酸、中性、堿性和強堿水平的分布頻率分別為2.19%、55.12%、31.14%、9.65%、1.90%和0.00%，小?？Х戎鳟a(chǎn)區(qū)的3個海拔高度的土壤pH：低海拔>中海拔>高海拔，三者間存在顯著差異，表現(xiàn)為土壤pH隨著海拔的升高逐漸降低的趨勢。

小粒咖啡最適宜種植的土壤pH為5.5～6.5，主產(chǎn)區(qū)仍有57.31%的土壤酸性過強，11.55%的土壤堿性過強。低海拔地區(qū)21.43%的土壤過酸現(xiàn)象，比第二次土壤普查時增加15.09%;41.84%的土壤過堿，比第二次土壤普查時減少18.73%（圖1B）。中海拔地區(qū)60.91%的土壤過酸現(xiàn)象，比第二次土壤普查時增加9.72%;6.91%的土壤過堿，比第二次土壤普查時增加6.12%（圖1C）。高海拔地區(qū)100%的土壤過酸，比第二次土壤普查時增加49.86%（圖1D）。

2.2? 土壤的有機質(zhì)含量分析

有機質(zhì)作為衡量土壤肥力高低的關(guān)鍵因素之一，其含量高低影響土壤的質(zhì)量，同時制約著作物的生長發(fā)育。目前云南咖啡主產(chǎn)區(qū)的土壤有機質(zhì)含量4.32～79.63 g/kg（表2），呈正態(tài)分布（圖2A），平均含量為32.68 g/kg，變異系數(shù)為51.30%，屬中等變異。云南咖啡主產(chǎn)區(qū)的土壤有機質(zhì)含量極缺、很缺、缺乏、中等、豐富和很豐富水平的分布頻率分別為0.58%、3.95%、25.15%、21.93%、14.91%和33.48%，3個海拔范圍內(nèi)的有機質(zhì)含量高海拔>中海拔>低海拔，三者間存在顯著差異，隨著海拔的升高有機質(zhì)含量逐漸升高。

小?？Х戎鳟a(chǎn)區(qū)的土壤有機質(zhì)由第二次土壤普查時的50.51 g/kg，減少到目前的32.68 g/kg，降幅35.30%。低海拔地區(qū)32.65%的土壤存在有機質(zhì)缺乏現(xiàn)象，比第二次土壤普查時增加32.65%（圖2B）。中海拔地區(qū)30.18%的土壤存在有機質(zhì)缺乏現(xiàn)象，比第二次土壤普查時增加29.39%（圖2C）。高海拔地區(qū)13.89%的土壤存在有機質(zhì)缺乏，比第二次土壤普查時增加13.89%（圖2D）。

2.3? 土壤的堿解氮含量分析

土壤堿解氮含量作為農(nóng)業(yè)活動中衡量氮肥施用量的主要指標，其高低反映出土壤的短期供氮能力及作物的養(yǎng)分吸收狀況。云南咖啡主產(chǎn)區(qū)的土壤堿解氮含量為9.90～286.62 mg/kg，呈正態(tài)分布（圖3A），平均為128.77屬中等變異（表2）。云南咖啡主產(chǎn)區(qū)的土壤堿解氮，極缺、很缺、缺乏、中等、豐富和很豐富水平的分布頻率分別為2.19%、2.19%、13.45%、23.39%、29.68%和29.09%，3個海拔范圍內(nèi)的堿解氮含量依次為：低海拔≈中海拔>高海拔，高海拔與中低海拔間差異顯著，整體呈現(xiàn)出隨著海拔的升高堿解氮量逐漸降低的趨勢。中海拔地區(qū)的堿解氮變幅最大，為9.90～286.62 g/kg，與整個產(chǎn)區(qū)的變幅一致。

小?？Х戎鳟a(chǎn)區(qū)的土壤堿解氮由30年前的159.92 mg/kg，減少到目前的128.77 mg/kg，降幅19.48%。低海拔地區(qū)14.29%土壤存在堿解氮缺乏現(xiàn)象，比第二次土壤普查時減少1.56%（圖3B）。中海拔地區(qū)18.36%的土壤存在堿解氮缺乏現(xiàn)象，比第二次土壤普查時增加17.57%（圖3C）。高海拔地區(qū)19.44%的土壤存在堿解氮缺乏現(xiàn)象，比第二次土壤普查時增加19.44%（圖3D）。

2.4? 土壤的有效磷含量分析

土壤有效磷能定量反映土壤對作物的當季供磷能力高低，對指導科學施肥具有積極意義。云南咖啡主產(chǎn)區(qū)土壤有效磷含量為0.49～92.06 mg/kg，呈正態(tài)分布（圖4A），平均為18.01 mg/kg，變異系數(shù)為97.20%，是咖啡土壤肥力中變異系數(shù)最高的指標（表2）。云南咖啡主產(chǎn)區(qū)的土壤有效磷含量極缺、很缺、缺乏、中等、豐富和很豐富的分布頻率分別為14.33%、10.53%、15.94%、29.82%、18.71%和10.67%，3個海拔范圍內(nèi)的有效磷含量依次為：高海拔>中海拔>低海拔，三者之間存在顯著差異，隨著海拔的升高有效磷含量逐漸升高。中海拔地區(qū)的有效磷變幅最大為0.49～92.06 mg/kg，與整個產(chǎn)區(qū)的變幅一致。

小?？Х戎鳟a(chǎn)區(qū)的土壤有效磷含量由30年前的9.93 mg/kg，增加到目前的18.01 mg/kg，增幅81.37%。低海拔地區(qū)38.78%土壤有效磷缺乏現(xiàn)象，比第二次土壤普查時減少45.73%（圖4B）。中海拔地區(qū)40.18%土壤有效磷缺乏現(xiàn)象，比第二次土壤普查時減少15.77%（圖4C）。高海拔地區(qū)55.56%土壤有效磷缺乏現(xiàn)象，比第二次土壤普查時減少27.82%（圖4D）。

2.5? 土壤的速效鉀含量分析

土壤的速效鉀含量作為評價土壤對作物的當季供鉀能力高低的重要指標，對指導農(nóng)業(yè)生產(chǎn)活動中合理施用鉀肥具有重要意義。云南咖啡主產(chǎn)區(qū)的土壤速效鉀含量為43.34～317.49 mg/kg，呈正態(tài)分布（圖5A），平均為138.31 mg/kg，變異系

數(shù)為28.94%，屬較低變異（表2）。云南咖啡主產(chǎn)區(qū)的土壤速效鉀含量極缺、很缺、缺乏、中等、豐富和很豐富水平的分布頻率分別為0.00%、0.29%、17.54%、46.49%、30.12%和5.56%，3個海拔范圍內(nèi)的速效鉀含量依次為：中海拔>高海拔≈低海拔，中海拔與高低海拔間差異顯著，隨著海拔的升高速效鉀含量先升高后降低。中海拔地區(qū)的速效鉀變幅最大為43.34～317.49 mg/kg，與整個產(chǎn)區(qū)的變幅一致。

小?？Х戎鳟a(chǎn)區(qū)的土壤速效鉀含量由30年前的156.59 mg/kg，減少到目前的138.31 mg/kg，降幅為11.67%。低海拔地區(qū)25.51%土壤速效鉀缺乏現(xiàn)象，比第二次土壤普查時增加25.51%（圖5B）。中海拔地區(qū)16.55%的土壤速效鉀缺乏，比第二次土壤普查時增加8.61%（圖5C）。高海拔地區(qū)16.67%的土壤速效鉀缺乏，比第二次土壤普查時增加16.67%（圖5D）。

3? 討論

小?？Х戎鳟a(chǎn)區(qū)的土壤肥力狀況具有明顯的海拔差異性。成土母質(zhì)、地形地貌和生產(chǎn)活動等不同的生態(tài)環(huán)境顯著影響著土壤的肥力特性[11]。低海拔地區(qū)，土壤多為黃土、赤土和紅灰泥土，該區(qū)域土體中含有石英顆粒，有機質(zhì)含量低，土壤中的碳酸鈣含量增多，降低有效磷含量[12]。中海拔地區(qū)，土壤以赤土、黃棕土為主，該區(qū)域土壤偏酸性，保水保肥能力較強，土體熟化程度較高，有機質(zhì)含量高[12]，是小粒咖啡的主要生產(chǎn)地段。高海拔地區(qū)，土壤大多數(shù)為黃紅壤，質(zhì)地偏黏，該區(qū)域多為新開墾的林地，有機質(zhì)含量豐富，酸性強。高海拔地區(qū)咖啡生長期長[13]，利于生物堿和酚性物質(zhì)合成[14]，咖啡杯品厚度佳，成為生產(chǎn)精品咖啡的黃金地帶。

土壤酸堿度不僅影響土壤中的微生物活性，而且制約著土壤養(yǎng)分的形成、轉(zhuǎn)化、吸收與利用，成為判別土壤肥力高低的主要指標之一[15]。30年內(nèi)，土壤pH下降0.5個單位，主要原因是長期過量施用氮肥[16]。土壤酸化現(xiàn)象在高海拔地區(qū)比較顯著，土壤過酸會影響咖啡對土壤養(yǎng)分的吸收與利用，阻礙根系生長，大量咖啡收獲加速土壤礦物風化并過度消耗堿性礦物養(yǎng)分庫，最終導致土壤pH不斷降低[17]。可結(jié)合當?shù)匦夂蛱攸c進行相應的養(yǎng)分優(yōu)化管理和種植調(diào)整[3]，采用咖啡與橡膠、龍眼、牛油果、西番蓮的復合種植模式，充分利用作物的空間分布特征。施用石灰石和生物炭混合改良劑、生物調(diào)理劑，土壤pH改善至5.0～5.5，土壤有機質(zhì)、速效磷、堿解氮分別提高了25%～39%、38%～74%、7%～13%[18]。

土壤有機質(zhì)是組成土壤肥力的重要指標，是肥力形成的實質(zhì)，同時影響并制約著土壤的理化性質(zhì)[2]。自全國第二次土壤普查以來，云南熱區(qū)的耕地面積不斷減少，優(yōu)質(zhì)耕地鹽堿化、沙化或占用，土壤肥力差的土地被開墾為耕地，致使整體咖啡產(chǎn)區(qū)的耕地質(zhì)量不斷下降，重化肥輕有機肥也是關(guān)鍵因素之一。海拔升高改變了微生物群落結(jié)構(gòu)[3]，土壤有機質(zhì)含量隨海拔高度升高而升高[4-5]，高海拔地區(qū)的土壤活性有機碳含量和分配比例較高[19]，與本研究結(jié)果一致。

土壤堿解氮是氮肥轉(zhuǎn)化成作物可吸收利用的主要形式，其含量高低表征氮肥的豐缺和利用率狀況?？Х鹊氐牡适┯昧繀T乏，氮肥在分解過程中釋放CO2和有機酸，CO2部分被植物吸收利用，部分溶解于土壤產(chǎn)生各種酸，加快土壤酸化進程，但促進難溶物質(zhì)溶解，增加土壤的有效養(yǎng)分含量[15]。土壤有效磷缺乏導致植物的咖啡受到限制，過多則影響質(zhì)量和品質(zhì)的形成[20]。低海拔地區(qū)多為石灰性母質(zhì)，所以磷養(yǎng)分貧瘠[15， 21]。小?？Х茸钸m宜的氮磷鉀配比為25∶9∶19，咖啡對肥料的需求量是氮肥最多，其次是鉀肥，最少是磷肥。從生態(tài)友好型農(nóng)業(yè)發(fā)展的角度看，土壤磷含量需要保持在既能保證咖啡的高產(chǎn)、穩(wěn)產(chǎn)和優(yōu)質(zhì)，又能最大限度降低環(huán)境的污染破壞力[22]。土壤中的鉀不會發(fā)生形態(tài)的轉(zhuǎn)化，與易轉(zhuǎn)化的氮和隨雨水淋失易被土壤固定的磷不一致，因此能夠真實反映咖啡地的施肥現(xiàn)狀[15]，可采用磷素衡量監(jiān)控技術(shù)[23]，提高磷肥的有效利用率。中、高海拔地區(qū)土壤質(zhì)地較為黏重，土壤速效鉀含量高[24]?？Х葘︹浀男枨罅績H次于氮，但農(nóng)民長期重施氮磷肥輕鉀肥的習慣，導致消耗土壤中大量的鉀素，同時得不到及時的補充。中國普遍存在鉀肥缺乏現(xiàn)象，80%以上的鉀素存于作物秸稈中，因此咖啡果皮、咖啡渣、修剪的枝干等發(fā)酵腐熟還田作為一種重要的速效鉀資源，可有效增加土壤的有機質(zhì)和速效養(yǎng)分含量[25-26]。因地制宜，堅持咖啡帶走多少施用多少的原則，保證土壤養(yǎng)分的動態(tài)平衡。

咖啡主產(chǎn)區(qū)由低海拔向高海拔土壤酸堿度逐漸降低，交換性酸增加，酸性增強;土壤有機質(zhì)、有效磷含量隨海拔的升高而增加[27]，土壤速效鉀含量隨著海拔的升高呈現(xiàn)先升后降的趨勢，與焦?jié)櫚驳萚2]的研究結(jié)果不一致，可能與咖啡主產(chǎn)區(qū)所處的干熱河谷小氣候特征在不同海拔下的地形、植被結(jié)構(gòu)、溫度、水分和光照等環(huán)境異質(zhì)性息息相關(guān)。

4? 結(jié)論

（1）目前小粒咖啡主產(chǎn)區(qū)的土壤pH平均為5.53，屬酸性土壤;有機質(zhì)平均含量為32.68 g/kg，屬豐富水平;堿解氮平均含量為128.77 mg/kg，屬豐富水平;有效磷平均含量為18.01 mg/kg，屬中等水平;速效鉀平均含量為138.31 mg/kg，屬中等水平。

（2）土壤pH、有機質(zhì)、堿解氮、有效磷和速效鉀含量的主要分布范圍為4.5～5.5、>40 g/kg、120～150 mg/kg、10～20 mg/kg和100～150 mg/kg，比例分別為55.12%、33.48%、29.68%、29.82%和46.49%;全國第二次土壤普查時的主要分布范圍為5.5～6.5、>40 g/kg、>150 mg/kg、10～20 mg/kg和150～200 mg/kg，所占比例分別為44.91%、59.62%、58.14%、29.42%和62.23%。

（3）不同海拔高度比較，土壤pH和堿解氮含量表現(xiàn)為低海拔>中海拔>高海拔，有機質(zhì)、有效磷含量表現(xiàn)為高海拔>中海拔>低海拔，速效鉀含量表現(xiàn)為中海拔>高海拔≈低海拔地區(qū)。

參考文獻

[1] 全國農(nóng)業(yè)技術(shù)服務中心. 耕地質(zhì)量演變趨勢研究[M]. 北京： 中國農(nóng)業(yè)科技出版社， 2008： 1-15.

[2] 焦?jié)櫚玻?李朝周， 趙? 陽， 等. 海拔對隴南白龍江流域油橄欖園土壤肥力的影響[J]. 生態(tài)學雜志， 2018， 37（2）： 360-365.

[3] 郭永龍， 劉友兆， 王利環(huán). 華北山區(qū)不同海拔臺地不同土地利用方式下土壤肥力及脫鹽趨勢[J]. 水土保持學報， 2012， 26（6）： 131-134.

[4] 程? 浩， 張厚喜， 黃智軍， 等. 武夷山不同海拔高度土壤有機碳含量變化特征[J]. 森林與環(huán)境學報， 2018， 38（2）： 135-141.

[5] 申佳艷， 李小英， 袁? 勇， 等. 納板河自然保護區(qū)土壤酶對不同海拔、坡向的響應[J]. 水土保持研究， 2018， 25（1）： 111-119， 125.

[6] 薛沛沛， 王? 兵， 牛? 香. 大崗山不同海拔毛竹林土壤肥力的灰色關(guān)聯(lián)度分析[J]. 浙江農(nóng)業(yè)學報， 2013， 25（6）： 1354-1359.

高祥照， 馬常寶， 杜? 森. 測土配方施肥技術(shù)[M]. 北京： 中國農(nóng)業(yè)出版社， 2005： 1-7， 14-20.

鮑士旦. 土壤農(nóng)化分析（第3版） [M]. 北京： 中國農(nóng)業(yè)出版社， 2000： 30-107.

全國土壤普查辦公室. 中國土壤[M]. 北京： 中國農(nóng)業(yè)出版社， 1998： 843-984.

Aciego P J C， Brookes P C. Relationships between soil pH and microbial properties in a UK arable soil[J]. Soil Biology and Biochemistry， 2008， 40： 1856-1861.

焉? 莉， 王? 寅， 馮國忠， 等. 吉林省農(nóng)田土壤肥力現(xiàn)狀及變化特征[J]. 中國農(nóng)業(yè)科學， 2015， 48（23）： 4800-4810.

陳彥清， 楊建宇， 鄖文聚， 等. 國家尺度上基于地形因子的光溫及氣候生產(chǎn)潛力修正算法[J]. 中國農(nóng)業(yè)科學， 2016， 49（11）： 2082-2092.

李? 強， 張? 翀， 任志遠. 近15年黃土高原植被物候時空變化特征分析[J]. 中國農(nóng)業(yè)科學， 2016， 49（22）： 4352-4365.

方洪生， 周迎春， 蘇有健. 海拔高度對茶園環(huán)境及茶葉品質(zhì)的影響[J]. 安徽農(nóng)業(yè)科學， 2014， 42（20）： 6573-6575.

張紅桔， 馬閃閃， 趙科理， 等. 山核桃林地土壤肥力狀況及其空間分布特征[J]. 浙江農(nóng)林大學學報， 2018， 35（4）： 664-673.

Guo J H， Liu X J， Zhang Y， et al. Significant acidification in major Chinese croplands[J]. Science， 2010， 327（5968）： 1008-1010.

Lesturgez G， Poss R， Noble A， et al. Soil acidification without pH drop under intensive cropping systems in Northeast Thailand[J]. Agriculture Ecosystem and Environment， 2006， 114（2/4）： 239-248.

張寶林. 廬山云霧茶茶園土壤酸度和肥力改良措施研究[D]. 南昌： 南昌航空大學， 2017.

向成華， 欒軍偉， 駱宗詩， 等. 川西沿海拔梯度典型植被類型土壤活性有機碳分布[J]. 生態(tài)學報， 2010， 30（4）： 1025-1034.

Hinsinger P. Bioavailability of soil inorganic P in the rhizosphere as affected by root-induced chemical changes： a review[J]. Plant Soil， 2001， 237（2）： 173-195.

Braschi I， Ciavatta C， Giovannini C， et al. Combined effect of water and organic matter on phosphorus availability in calcareous soils[J]. Nutr Cycl Agroecosys， 2003， 67（1）： 67-74.

胡? 斌， 和樹莊， 陳春瑜， 等. 滇池流域土壤氮磷分布特征及關(guān)鍵影響因素研究[J]. 土壤學報， 2012， 49（6）： 1178-1184.

賈良良， 孫彥銘， 劉克桐， 等. 河北省不同生態(tài)區(qū)農(nóng)田土壤肥力現(xiàn)狀及變化特征[J]. 土壤通報， 2018， 49（2）： 367-376.

Horst W J， Kamh M， Jibrin J M， et al. Agronomic measures for increasing P availability to crops[J]. Plant and Soil， 2001， 237： 211-223.

Liu C， Lu M， Cui J， et al. Effects of straw carbon input on carbon dynamics in agricultural soils： a meta- analysis[J]. Global Change Biology， 2014， 20（5）： 1366-1381.

Yu C J， Qin J G， Xu J， et al. Straw combustion in circulating fluidized bed at low-temperature： Transformation and distribution of potassium[J]. Canadian Journal of Chemical Engineering， 2010， 88（5）： 874-880.

劉? 穎， 宮淵波， 李? 瑤， 等. 川西高寒灌叢草地不同海拔梯度土壤化學計量特征[J]. 四川農(nóng)業(yè)大學學報， 2018， 36（2）： 167-174.