青島中心城區綠地土壤養分特征與評價

李俊翰,高明秀*,李滬波

青島中心城區綠地土壤養分特征與評價

李俊翰1,2,高明秀1,2*,李滬波3

1. 山東農業大學資源與環境學院, 山東 泰安 271018 2. 土肥資源高效利用國家工程實驗室, 山東 泰安 271018 3. 青島市園林環衛技術學校, 山東 青島 266000

城市綠地土壤是城市生態系統的營養基礎，掌握其理化特征及養分狀況對加強城市綠地土壤管理，改善城市綠地質量和環境健康、增進城市生態可持續發展具有重要意義。以青島市東海岸主城區核心區域的城市綠地土壤為研究對象，基于實地調查和室內化驗，分析土壤理化特征，采用主成分分析和模糊隸屬度模型綜合評價土壤養分狀況，基于ArcGIS10.2運用反距離權重法進行空間插值模擬土壤養分空間分布。結果表明：研究區城市綠地砂土、壤土占比分為54.71%和45.29%，部分地段土壤板結。45.28%的綠地土壤pH值偏高；土壤含鹽量較低，呈非鹽漬化（83.02%）和輕度鹽漬化（16.98%），但道路綠地受氯化鈉融雪劑影響含鹽量偏高；土壤養分總體處于中等水平，有機質、堿解氮、有效磷、速效鉀含量的樣點平均值均在三級肥力水平及以上，但部分樣點區域含量較低。綠地類型是影響土壤養分含量的主要因素，海拔高度和海岸距離對研究區綠地土壤理化性質影響不顯著，在人為干擾等強烈作用下質地對養分的影響有悖于一般規律。土壤養分綜合指數介于0.326?~?0.647，總體處中等水平，高、中、低肥力區面積比近似1:1:1。研究區綠地土壤養分指數大致呈現北高南低、西高東低的空間特征，并在中部區域形成一條凸向內陸的近似“M”的弧形高肥力區帶。

城市綠地; 土壤養分; 主成分分析; 養分評價

城市作為人類生產、生活的最主要活動場所和國家重要的組成單元，其生態環境質量的好壞不僅關系到人民生活品質的高低，而且關系到國家的生態文明建設和發展。城市園林綠地是城市生態系統的重要組成部分，在改善環境質量、調節系統平衡、美化城市景觀等方面起到重要作用[1,2]。現階段針對城市園林綠地的研究多集中于景觀設計規劃[3,4]、綠地生態系統結構[5]和功能[6,7]等方面，城市土壤作為城市生態系統的基礎物質承載者，特別是城市園林綠地土壤為城市植被生長提供了必需的附著空間和營養基礎，其質量的高低影響整個生態系統的自我調節能力和系統服務價值的大小[8]。

城市園林綠地土壤主要包括3種類型，即城市綠地土壤、保護地土壤和盆栽土壤。本文選取城市綠地土壤作為研究對象，主要是指公園、苗圃及專用綠地的土壤。這類土壤受降水量、溫濕度等城市自然因素影響較大，同時還受到建筑、堆積、踐踏、挖掘等人為因素干擾，原土層分布已無規律可循，土壤的三相比例發生顯著變化，土壤肥力遭到不同程度的破壞。因此，研究城市綠地土壤的理化特征和養分狀況，可為土壤改良和園林綠化建設提供科學依據，以改善城市生態環境質量，促進城市生態系統的健康、可持續發展。

1 材料與方法

1.1 研究區概況

青島市（35°35′ ~ 37°09′N，119°30′ ~ 121°00′E）地處山東半島東南部，瀕臨黃海，是環渤海經濟圈和山東半島藍色經濟區等多項戰略規劃的重要組成城市，是國家沿海重要的中心城市。本文選取的研究區位于青島東海岸主城區的核心區域，范圍北至青島港鐵路-杭鞍高架路-市北區北界，東至浮山-海青路，南、西至海邊，地理坐標（36°02′ ~ 36°07′N，120°16′ ~ 120°28′E），總面積62.05 ?km2（圖1）。該區是全市的政治、文化、金融中心，青島市人民政府所在地。西臨膠州灣，南臨黃海，地勢北高南低，東高西低，起伏不定，分布有浮山、太平山、青島山等山頭。土壤以棕壤為主并有部分鹽漬土，棕壤多分布于山地丘陵等地勢較高的區域，鹽漬土多分布于沿海低洼區，呈帶狀分布。研究區屬溫帶季風氣候，且有顯著海洋性氣候特點，四季分明，溫度適中，年均降雨量424.6? mm。

圖 1 研究區位置和樣點分布圖

1.2 土壤樣品采集與化驗

1.2.1 土壤樣品采集樣品于2016年12月采集，在研究區內的公園、廣場、街頭綠地、道路沿線參照綠地分布、街區分隔情況進行布點，使樣點遍布研究區，采樣點如圖1所示。采樣時，根據綠地分布情況，選擇有代表性的樣點，在100~200 m2左右范圍內，采用“S”形方法采樣；采樣深度為0~30? cm，將所采土樣進行多點（5~10點）混合，然后用四分法，對角線取樣品1 ?kg左右裝袋密封。共采集土壤樣品53個。根據中華人民共和國行業標準《城市綠地分類標準（CJJ/T 85—2002）》，結合研究區綠地實際情況，將綠地分為公園綠地、街旁綠地、道路綠地、單位綠地4種類型。單位綠地指位于某單位內部，由單位自行管理的綠地，取4個樣本；道路綠地指道路廣場用地內的綠地，包括行道樹綠帶、分車綠帶、交通島綠地、交通廣場和停車場綠地等，并將道路防護綠地（團島路防護林）歸入道路綠地，共取12個樣本；公園綠地指向公眾開放，以游憩為主要功能，兼具生態、美化、防災等作用的綠地，共取19個樣本；街旁綠地指位于城市道路用地之外，相對獨立成片的綠地，包括街道廣場綠地、小型沿街綠化用地等，其中“街道廣場綠地”與“道路綠地”中的“廣場綠地”不同，“街道廣場綠地”位于道路紅線之外，而“廣場綠地”在城市規劃的道路廣場用地（道路紅線范圍）內，共取18個樣本。公園綠地以原生土為主，而單位綠地、道路綠地、街旁綠地以客土為主。

1.2.2 土壤樣品化驗土壤樣品測試分析按《土壤分析技術規范》[9]的規定，參考《土壤農化分析》[10]的方法進行。土壤質地測定采用過1 mm孔篩風干土樣，加分散劑（酸性土壤加0.5 ?mol·L-1NaOH，堿性土壤加0.5 ?mol·L-1六偏磷酸鈉，中性土壤加0.5 ?mol·L-1草酸鈉），用吸管法測定；土壤pH采用水土比2.5:1，酸度計法；土壤含鹽量采用水土比5:1，烘干稱重法；有機質采用油浴加熱重鉻酸鉀—硫酸溶液氧化，硫酸亞鐵滴定測定；堿解氮采用堿解擴散法；有效磷采用0.5? mol·L-1NaHCO3浸提，鉬銻抗比色法；速效鉀采用1 ?mol·L-1NH4AC浸提，火焰光度計法測定。

1.3 數據來源與處理

DEM數字高程數據來自地理空間數據云平臺（http://www.gscloud.cn），數據日期2009年，采用通用橫軸墨卡托（UTM）投影，1984年世界大地（WGS84）坐標系統，空間分辨率30 ?m。借助ArcGIS10.2軟件平臺進行采樣點海拔高度和海岸距離數據提取，運用反距離權重（Inverse Distance Weighting，IDW）空間插值功能進行指標的空間分布模擬[11]。數據統計與相關分析采用SPSS軟件。

1.4 城市綠地土壤養分綜合評價

選取pH值、有機質、堿解氮、有效磷和速效鉀5項指標進行青島城市綠地土壤養分評價，采用主成分分析和模糊隸屬度模型進行綜合評價[12]。

1.4.1 評價指標權重運用主成分分析法確定評價指標權重[13,14]，排除人為因素干擾，體現指標權重的客觀性。具體步驟如下：（1）運用SPSS軟件進行主成分分析，提取主成分方差貢獻率和累計貢獻率，并計算各主成分權重；（2）提取成份得分系數矩陣，結合主成分權重，計算各評價指標的因子得分；（3）對因子得分進行標準化處理（式1），確定為評價指標權重。

其中，EIW是評價指標權重，a是各評價指標的因子得分，為指標個數。

1.4.2 評價指標隸屬度運用隸屬度函數對各評級指標進行量綱歸一化處理，消除因為指標不同所帶來的數據差異性。其中，pH值適用于拋物線型函數（式2），此類指標存在一個最佳適宜范圍，指標值偏離程度越大對土壤肥力的負向影響越大。有機質、堿解氮、有效磷和速效鉀適用于“S”型函數（式3），該指標在一定范圍內與土壤肥力呈現正相關關系，偏離此范圍對土壤肥力狀況影響很小。

其中，為隸屬函數，為評價指標的實際測定值，1、2、3和4分別為評價指標的臨界值。根據全國第二次土壤普查推薦的土壤肥力分級標準，結合研究區實際情況，確定指標臨界值。pH值以微酸性下限值（5.5）、中性下限值（6.5）、中性上限值（7.5）和弱堿性上限值（8.5）為臨界值；有機質、堿解氮、有效磷和速效鉀等養分指標，分別以二級上限值和五級下限值作為臨界值。

1.4.3 土壤養分綜合評價指數依據土壤肥力評級指標權重和隸屬度計算研究區土壤養分綜合評價指數。其中，是土壤養分綜合評價指數，EIW和DMF分別為第項指標的權重值和隸屬度值。

2 結果與分析

2.1 土壤物理性質特征

土壤物理性質主要通過野外采樣時觀察土壤表觀狀況和室內土壤質地測定進行描述分析。土壤質地對土壤理化性質和養分供應具有顯著的影響作用，主要表現在對土壤水分、養分、空氣、吸附性、耕性及植物生長的影響上[15,16]。青島城市綠地土壤質地總體上偏砂，顆粒較粗。在全部53個樣品中，砂土25個、占47.17%，壤砂土4個、占7.55%；砂壤土23個、占43.40%，壤土只有1個、占1.89%。由圖2可知，在水平方向上，青島市園林土壤質地分布較為復雜，各種質地交錯分布，但總體表現為沿海地區多砂質土壤，內陸地區多壤質土壤，土壤顆粒由沿海向內陸逐漸變小的變化規律。在垂直方向上，砂質土壤多分布在海拔較低的沿海區域，壤質土壤多分布在海拔較高的山丘公園、坡地等區域。

在調查中發現，砂質土土層較薄、粗骨物質多，通氣透水性好保水力弱，有利于好氣性微生物的活動。因含水量少，熱容量較小，所以晝夜溫差變化大，土溫變化快，這對于某些植物生長不利。壤土和部分砂壤土人為壓實嚴重。這導致表層以下好氣性微生物活性降低，土壤熟化緩慢、密實程度高，破壞通透性良好的團粒結構，形成理化性能較差的密實、板結的片狀或塊狀結構，出現板結堅硬、透氣性能差、循環不暢現象，阻礙園林植物根系的正常伸展[17]。

圖 2 土壤質地分布圖

2.2 土壤化學性質特征

土壤化學性質是影響土壤肥力的重要因素，青島城市綠地土壤化學性質描述性統計見表1。由表1可見，根據標準差判斷，各指標組內個體間離散程度大小依次是：速效鉀>堿解氮>有效磷>有機質> pH>含鹽量。根據變異系數判斷標準（<0.1弱變異性，0.1<<1中等變異性，>1強變異性），所有項目的變異系數均大于0.1小于1，屬于中等程度變異。各指標變異程度大小依次為：含鹽量>有效磷>堿解氮>速效鉀>有機質> pH。

表 1 青島城市綠地土壤化學性質描述性統計表

2.2.1 土壤酸堿性土壤pH平均值7.23，區間范圍4.79~8.98。53個樣本中有24個（45.28%）超過7.50，最大值高達8.98，分析原因，其一可能是因為園林綠地受到擾動大，特別是建筑垃圾含量高（部分建筑垃圾往往含有較多的石灰等堿性物質）引起內陸棕壤的酸堿性不正常；其二，可能是沿海的部分鹽漬土在較長時間的施肥影響和淡水淋洗后鹽分減少而堿性離子增加導致[18]。pH高于8.0，堿性較高，不利于園林植物生長，應注意通過施用酸性肥料或調理劑進行調控，使土壤pH值接近中性（6.5?~?7.5）促進植物生長。

2.2.2 土壤含鹽量53個樣點中非鹽漬化土壤樣點有44個，占83.02%，說明青島雖然是海濱城市，但海岸為砂石海岸，距離海邊一定距離后土壤無自然鹽漬化的條件。輕度鹽漬化樣點9個、占16.98%，多數園林植物生長不受影響。無中、重度鹽漬化土壤和鹽土。值得注意的是，有5個樣點土壤含鹽量超過1.50 ?g/kg，占9.43%。這5個樣點有的分布于海邊，受到海水影響含鹽量偏高。但靠內陸的個別樣點鹽分含量異常，達1.90 ?g/kg以上，接近中度鹽漬化水平，分析原因，這幾個樣點均位于橋梁或斜坡附近，可能是因為冬季降雪，城市管理中使用工業鹽作為融雪劑，而后在這些地方積累所致。此外，園林土壤過量施用化肥，也可能使得園林土壤中無機鹽離子增多，pH值升高，鹽漬化加重，這會導致植物形成生理干旱，產生鹽離子毒害，抑制園林植物的生長和發育。

2.2.3 土壤養分指標有機質、堿解氮、有效磷和速效鉀是土壤養分的重要組成成分，是衡量土壤肥力的重要指標。根據全國第二次土壤普查推薦的土壤肥力分級標準，青島城市綠地土壤肥力分級結果見表2。有機質含量介于4.62~43.19 g/kg之間，平均值26.86 g/kg，屬于三級肥力標準。低于6 g/kg的有2個樣點，均為公園綠地，但僅有少量喬木和灌木，無地被草本植物；高于30 g/kg的有24個樣點，占45.28%。堿解氮含量介于16.17~371.88 mg/kg之間，平均值為105.71 mg/kg，屬于三級肥力標準。低于30 mg/kg的有2個樣點，均為公園綠地，但僅有少量喬木和灌木，無地被草本植物；120~150 mg/kg的有9個樣點，占16.98%；高于150 mg/kg的有6個樣點，占11.32%。有效磷含量介于1.86~65.06 mg/kg之間，平均值為20.87 mg/kg，屬于二級肥力標準。低于3 mg/kg的有1個樣點，3~5 mg/kg的有3個樣點，分別為街頭或道路綠地、山地；20~40 mg/kg的有15個樣點，占28.30%；高于40 mg kg-1的有8個樣點，占15.09%。速效鉀含量介于32.61~397.15 mg/kg之間，平均含量141.51 mg/kg，屬于三級肥力標準。低于50 mg/kg的有3個樣點，分別為街頭或道路綠地、公園綠地；50 ~ 100 mg kg-1的有15個樣點，占比33.96%；150~200 mg/kg的有6個樣點，占11.32%；高于200 mg/kg的有11個樣點，占20.75%。總體而言，土壤養分含量居于中等水平，部分樣點區域含量較低，需要注意施肥補充。

表 2 青島城市綠地土壤養分分級統計

2.3 土壤化學性質的影響因素分析

2.3.1 不同綠地類型對土壤化學性質的影響對不同綠地類型下的土壤化學性質進行對比分析（表3），可見，不同綠地類型由于其位置、管理方式等的不同，對土壤化學性質有明顯的影響。道路綠地含鹽量最高，平均0.83 g/kg，雖然整體上處于非鹽漬化水平，但個別樣點鹽漬化已達中度邊緣。北方城市冬季常以工業鹽（氯化鈉）作為融雪劑，氯化鈉在道路兩側、特別是橋梁或斜坡下方積累，使道路綠地土壤含鹽量升高。單位綠地含鹽量最低，但pH最高，可能是單位綠地往往建筑垃圾較多，或者一般單位對綠地的管理僅重視澆水，在降低土壤鹽分的同時，pH有所升高。公園綠地含鹽量平均為0.58 g/kg，個別樣點處于海邊，可能受到海潮影響，浪花濺入帶進鹽分所致。從土壤養分含量來看，各類型綠地養分含量平均值均在三級以上，與青島城市綠地有較好的施肥管理有關。有機質含量平均值均三級以上肥力水平，其中道路綠地最高，可能與道路綠地一般處于圍欄之中，較少人為踐踏擾動、枯枝落葉積累較多有關；公園綠地最低，主要原因是公園綠地人為干擾較多，地被植物人為踐踏較多、枯枝落葉常有清掃，有機質積累相對較少。公園綠地堿解氮含量最高，道路綠地與之相近。公園綠地一般管理較其他類型綠地多，施肥以氮肥為主，但多為丘陵、海岸公園，地形起伏較大，存在土壤和養分流失問題，故有效磷、速效鉀含量并不是最高的。街旁綠地有效磷、速效鉀含量最高，因為街旁綠地一般地勢起伏較小，水土流失較輕，再就是街旁綠地中有幾處為近年建成的街道廣場綠地，磷鉀肥施用時間較近，土壤中保存較多。

表 3 不同綠地類型土壤化學性質對比

2.3.2 不同質地類型對土壤化學性質的影響對不同質地類型下的土壤化學性質進行對比分析（表4），pH值由小到大依次為：砂壤土<砂土<壤砂土<壤土。含鹽量由小到大依次為砂壤土<壤砂土<砂土<壤土，與pH值有較高的一致性。有機質、堿解氮、有效磷和速效鉀都表現為砂土含量最高，壤土含量最低。其中有機質和有效磷表現為壤砂土含量低于砂壤土含量，堿解氮和速效鉀表現為壤砂土含量高于砂壤土含量。綜合來看，堿解氮和速效鉀在一定程度上較為明顯的表現出質地越輕含量越高的特征，有機質和有效磷含量總體上也表現出砂土比壤土高的現象。

表 4 不同質地類型土壤化學性質對比

2.3.3 樣點指標自相關分析通過對青島城市綠地土壤指標的相關性分析，各指標之間存在著一定程度的相關關系，相關系數結果如表5所示。pH值與各指標之間都存在著負相關關系，特別是與堿解氮存在極顯著的負相關關系，與其它指標之間的負相關關系并不明顯；含鹽量與pH值、有機質、海拔高度表現出不顯著的負相關關系，與速效鉀、海岸距離表現出不顯著的正相關關系，而與堿解氮呈顯著正相關，與有效磷呈極顯著正相關；有機質與pH值、含鹽量呈不顯著負相關，與其他指標呈不顯著正相關；堿解氮與pH值呈不顯著負相關，與其他指標呈現正相關關系，特別的與含鹽量、有效磷呈顯著正相關，與速效鉀呈極顯著正相關；有效磷與pH值、海拔高度、海岸距離呈現不顯著的負相關關系，與有機質呈不顯著正相關，與堿解氮呈顯著正相關，與含鹽量速效鉀呈極顯著正相關；速效鉀與pH值、海拔高度呈不顯著負相關，與含鹽量、有機質、海岸距離呈不顯著正相關，與堿解氮、有效磷呈極顯著正相關；海拔高度與海岸距離之間存在極顯著正相關。海拔高度和海岸距離與其他指標之間的相關性不顯著，說明二者并不是影響青島城市綠地土壤養分的重要影響因素。

表 5 青島城市綠地土壤指標相關系數

備注：**顯著性水平<0.01，即極顯著水平；*顯著性水平<0.05，即顯著水平。

Note: **Very significance level<0.01; *significance level<0.05.

2.4 土壤養分綜合評價

選取pH值、有機質、堿解氮、有效磷、速效鉀5項指標進行青島城市綠地土壤養分的主成分分析。通過Bartlett球度檢驗，各指標的相關系數的相伴概率為0.00，適合進行主成分分析，共提取5個成分。參考相關文獻[19,20]，一般選取特征值大于1，或累計貢獻率大于80%的主成分進行綜合分析。本研究中第3主成分特征值0.882，累計貢獻率達到83.117%，故選取前3個主成分描述各項指標的數據信息特征。由表6可知，第1主成分反映堿解氮、有效磷和速效鉀的綜合狀況；第2主成分反映pH值狀況；第3主成分主要反映有機質的狀況，3個主成分能夠反映所有指標的數據信息。

表 6 青島城市綠地土壤指標的因子載荷矩陣

利用本文構建的土壤養分綜合評價指數模型計算研究區城市綠地土壤養分指數大小，基于arcgis10.2軟件平臺，運用反距離權重（Inverse Distance Weighting，IDW）空間插值法進行空間插值，結果如圖3。其中，53個樣點綜合評價指數最大值0.916，最小值0.192，平均值0.579，方差0.187。如圖3所示，土壤養分綜合評價指數大致呈現北高南低、西高東低的特征，并在中部區域形成一條凸向內陸的近似“M”的弧形高肥力區帶。采用自然裂點分級法將青島城市綠地土壤養分水平分為五級：極低[0,0.455）、低[0.455,0.544）、中[0.544,0.612）、高[0.612,0.699）、極高[0.699,∞）。其中，極高土壤養分區面積4.31 ?km2，占比6.96%；高土壤養分區面積17.84 ?km2，占比28.79%；中等土壤養分區面積18.74 ?km2，占比30.24%；低土壤養分區面積17.19 ?km2，占比27.74%；極低土壤養分區面積3.89? km2，占比6.27%。高、中、低肥力區面積比近似1:1:1。

圖 3 青島城市綠地土壤養分綜合評價得分

對不同質地和綠地類型的土壤養分評價，結果見表7。總體來說，土壤質地越輕土壤養分綜合評價指數越高，但并不嚴格按照此趨勢變化，土壤養分綜合評價指數從小到大依次是：壤土<壤砂土<砂壤土<砂土。不同的園林綠地類型，土壤養分也有所不同，土壤養分綜合評價指數從小到大依次是：公園<道路<單位<街旁。單位綠地和街旁綠地地形平坦，大多為新建綠地，施肥時間短，土壤保存多，且人工干預，特別是對綠地的培肥養護較多，肥力相對較高；公園綠地多山丘、坡地，地勢起伏不定，存在一定程度的水肥流失現象，土壤養分相對較低。

表 7 不同質地和綠地類型的土壤養分評價

3 討論

土壤養分是土壤肥力的重要物質基礎，是物理、化學和生物性質的綜合反映與作用結果，同時受自然因素和人類活動的影響，表現出明顯的空間分異特征[21,22]。研究區土壤質地從沿海低海拔到內陸高海拔表現為由砂質土向壤質土逐漸變化的規律。一般而言，土壤質地越輕，顆粒越粗，保水保肥能力越弱[23,24]。而青島城市綠地土壤總體表現為砂質土各養分含量高于壤質土各養分含量，并沒遵循土壤含量與質地之間的一般性變化規律，究其原因與人為擾動影響較大有關。自然環境中，土壤中氮、磷、鉀和有機質等各養分主要來源于枯枝落葉、植物枯萎根系和動物遺體等在微生物參與下的分解、利用。城市環境中，人為的清掃活動極大地削弱了這種土壤養分的自然來源途徑，人工施加肥料和土壤調理劑成為土壤養分的主要來源，并且以人類意志為轉移，打破了自然界的土壤養分分布規律。施工、廢棄物堆棄和垃圾掩埋等人類活動也會對土壤物理、化學和生物性質產生影響，促使土壤養分含量向有悖于自然規律的方向發生改變[25-27]。

相關研究表明，土壤養分是隨海拔變化的，氮、磷、鉀和有機質含量與海拔高度具有顯著的相關關系[28-30]，但本研究中二者相關性并不顯著。其一，研究區以低海拔丘陵和沿海平原為主，地勢較為平坦，海拔變化范圍僅為0 ~ 318m，變幅較小，不能很好地體現出土壤理化性質的海拔梯度變化特征；其二，城市土壤受人為因素干擾強烈，以不同城市綠地類型為尺度施加不同程度、不同作用方向人工干預，破壞了其原有的自然空間變化規律。

本文研究了青島東海岸主城區的核心區域的城市綠地土壤理化性質養分特征，為綠地土壤管理提供了依據，但僅是一個時間斷面，今后尚需連續監測，分析其連續性演變特征。并且，運用反距離空間插值法進行空間模擬還存在受樣點方位影響大、預測精度低和局部細節不明顯等局限性[31]，今后還需深入研究具有更高精度的空間預測方法[32]。

4 結論

（1）研究區綠地砂土、壤土占比分為54.71%和45.29%，部分地段存在土壤板結現象。研究區綠地45.28%的土壤pH偏高，不利于園林植物生長。土壤含鹽量較低，以非鹽漬化土（83.02%）和輕度鹽漬化土（16.98%）為主，無中、重度鹽漬化土壤和鹽土。土壤養分總體處于中等水平，有機質、堿解氮、有效磷、速效鉀含量的樣點平均值均在三級肥力水平及以上，但部分樣點區域含量較低，需要注意施肥補充，并控制水土流失；

（2）離海岸距離遠近和海拔高度不同的樣點土壤養分含量沒有明顯的變化規律，二者對研究區綠地土壤養分的影響不顯著。不同綠地類型因所處微地形的差異和所受人為因素的干擾不同，土壤養分存在較明顯差異。不同質地土壤養分含量存在一定差異，但受人為干擾等影響強烈而有悖于一般壤土比砂土肥力高的規律；

（3）土壤養分綜合評價結果顯示，研究區土壤養分綜合指數介于0.326?~?0.647，總體處于中等水平，高、中、低肥力區面積比近似1:1:1。從空間上看，研究區綠地土壤養分指數大致呈現北高南低、西高東低的特征，并在中部區域形成一條凸向內陸的近似“M”的弧形高肥力區帶。

[1] 毛齊正,羅上華,馬克明,等.城市綠地生態評價研究進展[J].生態學報,2012,32(17):5589-5600

[2] 蘇泳嫻,黃光慶,陳修治,等.城市綠地的生態環境效應研究進展[J].生態學報,2011,31(23):7287-7300

[3] 趙警衛,徐紋艷,林保英.基于公眾審美的城市綠地植物景觀設計依據[J].城市問題,2018,37(1):41-45,80

[4] 鄧位,李翔.英國城市綠地標準及其編制步驟[J].國際城市規劃,2017,33(6):20-26

[5] 黃玉梅,黃勝嵐,張健,等.成都市溫江區城市綠地不同植物配置下中小型土壤動物群落特征[J].中國科學院大學學 報,2018,35(1):33-41

[6] 張緒良,徐宗軍, 張朝暉,等. 青島市城市綠地生態系統的環境凈化服務價值[J].生態學報, 2011, 31(9):2576-2584

[7] Lin BB, Egerer M, Liere H,Local and landscape-scale land cover affects microclimate and water use in urban gardens[J]. Science of the Total Environment, 2018,47(1):610-611:570-575

[8] 張甘霖,趙玉國,楊金玲,等.城市土壤環境問題及其研究進展[J].土壤學報,2007,44(5):925-933

[9] 全國農業技術推廣服務中心.土壤分析技術規范[M].北京:中國農業出版社,2006

[10] 鮑士旦.土壤農化分析[M].北京:中國農業科技出版社,2000

[11] 劉光孟,汪云甲,王允.反距離權重插值因子對插值誤差影響分析[J].中國科技論文,2010,5(11):879-884

[12] 王舒甜,張金池,鄭丹揚,等.鐘山風景區土壤環境對人為踩踏擾動的響應[J].林業科學,2017,53(8):9-16

[13] 黃安,楊聯安,杜挺,等.基于主成分分析的土壤養分綜合評價[J].干旱區研究,2014,31(5):819-825

[14] 林海明,杜子芳.主成分分析綜合評價應該注意的問題[J].統計研究,2013,30(8):25-31

[15] 李志鵬,趙業婷,常慶瑞.渭河平原縣域農田土壤速效養分空間特征[J].干旱地區農業研究,2014(2):163-170

[16] 葛楠楠,石蕓,楊憲龍,等.黃土高原不同土壤質地農田土壤碳、氮、磷及團聚體分布特征[J].應用生態學 報,2017,28(5):1626-1632

[17] 閔祥宇.采煤塌陷地復墾土壤緊實特征及其土壤環境效應[D].泰安:山東農業大學,2018

[18] 王志剛,趙永存,廖啟林,等.近20年來江蘇省土壤pH值時空變化及其驅動力[J].生態學報,2008,28(2):720-727

[19] 劉毓,趙鳳蓮,劉紅權,等.濟南市城市綠地土壤肥力綜合評價及管理[J].山東農業大學學報(自然科學 版),2016,47(3):366-371

[20] 蘭忠明,王飛.多元統計分析在三明煙區土壤肥力評價中的運用[J].中國農學通報,2006,22(3):263-266

[21] 張晗,趙小敏,歐陽真程,等.江西省不同農田利用方式對土壤養分狀況的影響[J].水土保持研究,2018,25(6):53-60

[22] 梁斌,齊實.北京山區土壤養分空間變化特征研究[J].土壤,2018,50(4):769-777

[23] 葛楠楠,石蕓,楊憲龍,等.黃土高原不同土壤質地農田土壤碳、氮、磷及團聚體分布特征[J].應用生態學 報,2017,28(5):1626-1632

[24] 譚長強,彭玉華,申文輝,等.廣西都安地區5種森林類型土壤機械組成及其肥力比較[J].生態科學,2017,36(2):119-125

[25] 郭彥海,高國龍,王慶,等.典型平原地區生活垃圾焚燒廠周邊土壤中重金屬賦存形態分布及生物有效性評價 [J/OL].環境科學研究.[2018-10-22].http://kns.cnki.net/kcms/detail/11.1827.X.20181022.0946.002.html

[26] 葉紅云,熊黑鋼,張芳,等.基于CWT的人類不同干擾程度下干旱區土壤有機質含量估算研究[J].激光與光電子學進 展,2019,56(5):1-10

[27] 葉飛,吳勝軍,姜毅,等.人類活動對三峽消落帶土壤亞硝酸鹽型甲烷厭氧氧化菌群落的影響[J].環境科學學 報,2018,38(8):3266-3277

[28] 俞月鳳,何鐵光,宋同清,等.桂西北喀斯特地區石灰土養分空間變異特征[J].生態學報,2018,38(8):1-9

[29] 王云全,吳學群,彭占偉.區域耕地土壤養分空間變異研究及影響因素分析——以昌寧縣為例[J].熱帶作物學 報,2015,36(4):753-758

[30] 謝紅花,李超,錢曄,等.云南烏蒙山區土壤養分空間變異及海拔梯度分布規律[J].中國農學通報,2019,35(8):52-59

[32] 李正泉,吳堯祥.顧及方向遮蔽性的反距離權重插值法[J].測繪學報,2015,44(1):91-98

[33] 陸訪儀,趙永存,黃標,等.海倫市耕層土壤有機質含量空間預測方法研究[J].土壤通報,2012,43(3):662-667

Soil Nutrient Characteristics and Evaluation of Green Space in Qingdao Central Urban Area

LI Jun-han1,2, GAO Ming-xiu1,2*, LI Hu-bo3

1.271018,2.271018,3.266000,

The soil of urban green space is the nutritional basis of the urban ecosystem. It is of great significance to master the physical and chemical characteristics and the nutrient condition of the urban green space, to improve the urban green space quality，the environmental health and the sustainable development of the urban ecosystem. Taking the urban greenbelt soil in the core area of the main urban area of the Qingdao east coast as the research object, the physical and chemical characteristics of the soil were analyzed with the field investigation and laboratory test, and the soil nutrient status was comprehensively evaluated by principal component analysis and fuzzy membership model. Built on ArcGIS10.2, the spatial distribution of soil nutrients was simulated by spatial interpolation with inverse distance weight method. The results showed that the proportion of sandy soil and loam in the urban greenbelt in the study area was 54.71% and 45.29%, and the soil in some sections was consolidated, 45.28% of the green space soil has a high pH value; the soil salt content was low, showing non-salinization (83.02%) and mild salinization (16.98%). Road greenbelt salt content is high under the influence of sodium chloride snow melt agent. Soil nutrients were generally at the medium level, the average values of organic matter, alkali-hydrolyzable nitrogen, available phosphorus and available potassium were all at three levels of fertility and above, but the contents of some samples were relatively low. The type of green space is the primary factor affecting the soil nutrient content, while the effects of altitude and coastal distance on the physical and chemical properties of green soil in the study area were not significant. The effect of texture on nutrients is contrary to the general rule under the persuasive effect of human disturbance. The soil nutrient comprehensive index was between 0.326 and 0.647, and the overall level was medium. The area ratio of high, medium and low fertilizer area is about 1:1:1. The soil nutrient index of the greenbelt in the study area showed the spatial characteristics of soaring in the north and the west, low in the south and the east. There was a high arc fertility zone with convex inland approximately "M" in the central region.

Urban green space; soil nutrient; principal component analysis; nutrient evaluation

S151.9; S158.2

A

1000-2324(2019)04-0550-09

2018-10-18

2018-11-05

國家科技支撐計劃項目(2013BAD05B06-5);山東省重大科技創新工程項目(2017CXGC0301);山東省高校“雙一流”獎補資金項目(SYL2017XTTD02)

李俊翰(1993-),男,碩士研究生,主要研究方向為土地利用變化與生態安全. E-mail:junhanl@163.com

Author for correspondence. E-mail:mxgao@sdau.edu.cn