土壤有效態微量營養元素空間變異特征及對稻米品質的影響

江勝國，詹華明*，劉廣明，盧東琪，鄭愛軍，王志剛，王西玉

（1．天津市地質研究和海洋地質中心，天津 300170；2．天津市農業發展服務中心，天津 300061；3．中國科學院南京土壤研究所，江蘇 南京 210008）

土壤微量營養元素通常指有益于作物生長且需要量微小的元素，主要來源于土壤礦物質，以活動態形式賦存，是鑒定土壤質量的重要因子［1-2］。掌握稻田土壤微量營養元素有效態的空間變異特征和分布規律，研究其與稻谷品質之間的相關性，對指導微肥合理配施和實現稻米提質增產意義重大［3-4］。土壤是非均一性的連續體，屬性隨空間位置變化而發生變異，對其空間變異特性的研究難以定量化，隨著地統計學方法應用到土壤學領域，這一難題逐漸得到解決。20世紀90年代，我國把地統計學應用到土壤領域，取得較多成果［5］。張朝生等［6］認為天津平原區土壤微量元素含量的空間變異具有隨機性和結構性，且具有不同程度的空間自相關性；徐尚平等［7］認為天津土壤中微量元素含量主要受成土母質的影響，空間縱向分異不明顯。很多學者針對不同的土地利用方式也開展了較多土壤微量元素空間分布規律及主要影響因素分析的研究工作［8-10］。但前人針對稻田土壤微量營養元素的空間特征研究較少，且未結合稻米品質進行分析，特別是對天津稻田土壤微量元素空間變異性的研究還未見報道。因此，本次研究運用地統計學方法，結合GIS軟件，把復雜的空間變異特性定量化，系統探討了天津小站稻種植區土壤6種微量營養元素的空間變異特征和分布規律，并結合稻米品質進行了討論。

1 材料與方法

研究區位于天津寶坻區東南部，小站稻產業振興規劃區內（圖1），坐標為：117°46′～117°59′ E、39°41′～39°53′ N（圖1），面積約160 km2，主要包括黃莊鎮及其周邊區域，屬于河流沖積型和濱海型平原地貌，地勢比較平坦，地形趨勢為北高南低，高程2～7 m（大沽高程）。全區屬溫帶濕潤季風型大陸性氣候，全年光、熱、水資源豐富，年日照時數達2566.5 h，太陽總輻射量502.9 kJ，年平均氣溫11.3 ℃，年平均降水量492.2 mm，無霜期平均194 d。全年大于0℃積溫4507.9 ℃，大于10 ℃積溫4116.3 ℃，熱量條件有利于喜溫作物水稻的生長。區內水稻為一年一熟，種植歷史較長，土地利用程度高。土壤類型多為鹽化潮土，質地粘重，為中-重壤質，適合水稻生長。

圖1 研究區位置及采樣點分布圖

1.1 樣品采集與制備

共采集40個樣品（n=40），采集層位為耕作層，采集深度為20 cm。區內每4 km2設1個采樣區，每個采樣區內選定1個代表性點位，然后以代表性點位為中心等距（200 m）選定4個點位，采用“梅花采樣法”多點取樣，最后組合成一個樣品，樣品原始重量大于1000 g。把采集的土樣放置于陰涼通風處，干燥后敲碎成小塊，挑揀出雜質，然后研碎過0.85 mm篩，裝入樣袋貼好標簽，待測。水稻成熟后，選取采樣區代表性地塊采集水稻樣品，品種主要為津原U89，設置重復樣，樣品自然風干，水分降至15%，送至農業農村部稻米及制品質量監督檢驗測試中心進行品質分析。

1.2 分析測定方法

共測定6個微量營養元素指標（Cu、Mn、B、Mo、Ni和S，全量和有效量），5個稻米品質指標（透明度、整精米率、膠稠度、直鏈淀粉和蛋白質）。全Cu、全Mo和全Ni采用電感耦合等離子體質譜法測定；全Mn采用X-射線熒光光譜法測定；全B采用交流電弧-發射光譜法測定；全S采用EDTA間接滴定法測定。有效Cu、Mn和Ni采用DTPA浸提-原子吸收分光光度法測定；有效Mo采用Tammi試劑提取-等離子體質譜法測定；有效B采用沸水提取-等離子體光譜法測定；有效S采用酸性溶液浸提-硫酸鋇比濁法測定。pH值用水浸提樣品，采

用電位法測定。稻谷樣品采用韓國制造的SY-88-TH型礱谷機脫穎殼，用日本制造的CBS300AS型精米機碾成精米，精糙比控制在90%；稻米外觀指標采用日本生產的RJQI20型顆粒評定儀測定；稻米成分采用日本生產的RLTA1082-K型食味分析計測定。

1.3 數據處理與分析

應用SPSS 20.0進行數據的描述性統計分析、相關性研究以及正態分布檢驗；應用GS+9.0進行微量營養元素有效態的半方差函數及擬合參數的分析；應用ArcGIS 10.2地統計分析模塊中的克里格插值法進行空間分布規律的研究。

1.3.1 地統計分析

本次選取6種微量營養元素，運用半方差變異函數方法進行空間變異特性的研究。該方法首先要對數據進行半方差函數分析，選取適合的半方差函數模型，再運用克里格法進行空間插值。半方差函數主要的屬性模型包括：球狀模型、指數模型、高斯模型和線狀模型，表1為模型中的主要參數及反應的信息情況［11］?？死锔癫逯捣ㄊ窃趯ν寥涝紨祿治龅幕A上，結合半方差函數分析對數據優化處理的方法之一。該方法可以對研究區內未調查取樣部位的土壤信息進行線性無偏最優估計［10］，并可實現點信息向面信息的轉化，更直觀地呈現出元素的空間分布規律。

表1 半方差函數模型主要參數及反映的信息情況

1.3.2 稻米品質分析

本次選取透明度、整精米率、膠稠度、直鏈淀粉和蛋白質5個指標（表2），應用SPSS 22.0進行微量營養元素Cu、Mn、B、Mo、Ni和S有效態與稻米品質的相關性研究。分析認為，區內稻米透明度和整精米率達到了我國農業行業食用稻品種品質標準（NY/T 593-2013）所規定的一等級，而膠稠度和直鏈淀粉（干基）為二等級。蛋白質含量較高，均值為7.56%，有相關研究認為蛋白質含量在6%～7%質量最優［12］。

表2 主要稻米指標值的描述性統計特征

2 結果與分析

2.1 土壤微量營養元素有效態數據統計分析

通過文獻檢索，整理了土壤中S、Mn、B、Cu、Mo和Ni有效態含量評價標準［13-17］（表3）。對比數據統計的結果（表4），可知研究區有效S、Mn和B含量極為豐富，有效Cu和Ni含量較高，有效Mo含量中等水平。有效S含量均值210.10 mg/kg，處于極高水平，遠超出全國最高均值（25 mg/kg），這與研究區特殊的硫酸鹽-氯化物鹽漬土類型有關；有效Mn含量均值41.73 mg/kg，處于極高水平；有效B含量均值2.29 mg/kg，處于極高水平；有效Cu含量均值5.68 mg/kg，處于高水平；有效Mo含量均值0.20 mg/kg，處于中等水平；有效Ni含量均值0.55 mg/kg，處于高水平。從變異系數角度分析，有效S、Mn和Mo的變異程度較高，有效B、Cu和Ni的變異程度較低，說明有效S、Mn和Mo的含量分布不均，具有相對較高的空間異質性。

表3 土壤有效態微量元素含量評價標準 （mg/kg）

表4 研究區6種土壤有效態微量元素含量統計特征

2.2 空間變異性研究

元素數據的描述性統計只能概略說明含量的整體變化趨勢，不能具體表明空間變異的結構特征、分布規律及相關性、獨立性、隨機性和結構性等［9］，因此需要應用地統計學方法進一步研究。

2.2.1 空間結構特征

根據半方差函數分析的結果（表5，圖2），選取了擬合效果最好的模型參數（R2接近于1，RSS最?。?。結果顯示，有效S和Cu擬合效果最好的模型為指數模型，有效Mn、B和Ni為高斯模型，有效Mo為球狀模型。從塊金值分析，有效Mo和Ni主要受結構性因素影響，而有效B和Cu的變異在一定程度上還受隨機性因素的影響。有效S和Mn的基臺值較高，說明兩種元素的總變異程度高。6種元素的塊基比都＜25%，說明樣本間的變異總體主要受結構性因素影響，空間自相關性強烈［6］。有效Cu和S變程值最大，說明在較大的范圍內具有空間相關性，而有效Ni變程值最小，說明在較小的范圍內具有空間相關性。有效S、Mn、Mo和Ni的決定系數值最接近于1，擬合效果最理想，說明相關參數的選擇能較好地反映這些元素的空間結構特征。

表5 土壤有效態微量元素半方差函數及擬合參數

圖2 土壤有效態微量元素半方差函數

2.2.2 空間分布規律

依據半方差函數分析的結果，借助ArcGIS 10.2地統計分析模塊，選取普通克里格插值法繪制了6種土壤微量營養元素有效態的空間分布規律圖（圖3）。

圖3 土壤有效態微量元素含量空間分布規律

研究區微量營養元素有效態含量分布規律性較強，為面狀展布，除了有效S外（雙向異性），都在一定方向呈遞增或遞減的趨勢，這與前期數據描述性統計分析結果基本一致。有效S、Cu和Mo連續分布，基本呈北高南低的分布規律；有效Mn和Ni含量基本呈南高北低的分布規律；有效B含量呈東高西低的分布規律。在八里莊和小傾店區域有效S含量有雙向遞增的趨勢，這與農田的水稻種植年限有關，該區2018年實施的旱改水工程，水稻種植年限短，土壤的鹽漬化程度相對較高。有效Cu在局部區域呈異常相對高的斑塊狀，這可能由人為活動對土壤造成的污染所致［11］。有效Mn分布相對復雜，局部呈島狀，最高值分布在南部里自沽、東白莊和小傾店3個區域。有效B面狀分布規律性最為明顯，最大值分布在東部的八里莊和小傾店區域，最小值分布在西部的北里自沽和苑洪橋區域。

2.3 微量營養元素與稻米品質的相關性

稻米品質與肥料的施用關系較為密切［18］，研究認為微量營養元素肥的調控可產生稻米品質的差異（表6，圖4）。透明度是稻米外觀品質中最直觀的指標，與有效S呈顯著正相關（P＜0.05）。整精米率是稻米加工品質最重要的指標之一，與有效Mo呈極顯著負相關（P＜0.01）。膠稠度、直鏈淀粉和蛋白質與稻米蒸煮食味品質關系較為密切，其中，膠稠度與有效Cu呈顯著正相關（P＜0.05），與有效Mo呈極顯著正相關（P＜0.01），與有效Mn呈顯著負相關（P＜0.05）；直鏈淀粉與有效Mo呈顯著負相關（P＜0.05）；蛋白質與有效B呈顯著正相關（P＜0.05）。

表6 土壤有效態微量元素含量與稻米品質指標的相關性

圖4 有效態微量元素含量與稻米品質指標的相關性圖

3 討論

研究區耕地土壤為不同程度的鹽化潮土，成土母質主要來源于河流沉積物和海湖沉積物。土壤中各種微量營養元素初始含量豐富，有效態含量相對較高（表4）。與全國第二次土壤普查［19］的數據對比，區內有效S、Mn、B、Cu、Ni和Mo含量較豐富，說明了多年的作物種植和精耕細作顯著提升了土壤的肥力水平。這6種元素的塊基比值（表5）都小于25%，說明空間自相關性強烈［6］，空間變異主要受結構性因素的影響，受人為耕作和管理等因素的影響較小。土壤各元素分布規律性較強，多呈面狀展布，而有效S卻在兩個區域出現高值，有雙向異性的特征。分析認為，有效S受天津濱海區硫酸鹽-氯化物類型鹽漬化土壤內在屬性的影響較大。有機質和pH值對區內稻田土壤微量營養元素空間變異特征的影響起主導作用，土壤中有機質與微量營養元素有效態多呈正相關，而與pH值呈負相關［20］。

研究認為，研究區有效S與透明度呈顯著正相關，有效Mo與整精米率呈顯著負相關，有效Cu和Mo與膠稠度呈顯著正相關，有效Mo與直鏈淀粉呈顯著負相關，有效B與蛋白質呈顯著正相關。由數據分析可知，研究區土壤S、Mo、Cu、B元素的含量較為豐富，水稻的有效吸收量較好。從而，使得小站稻米的品質指標表現優異。前人對天津小站稻米的品質指標進行了深入研究，且與稻花香米和中日泰優質稻米在品質上進行了對比［12］，也證明了小站稻米的透明度、整精米率、膠稠度和蛋白質指標都相對較高。然而，較高的蛋白質含量對小站稻米的食味品質影響較大，前人研究認為，稻米蛋白質含量在6%～7%之間最好，超過7%食味就會降低［21］。研究區稻米蛋白質含量平均為7.56%，最高可達8.4%，嚴重影響了食味品質。通過實地調查，研究區內稻田在前期底肥和水稻生長期3次追肥的施加中，都有較多的氮肥，且土壤B 含量豐富，能提高作物的固氮能力［22-23］，使得水稻的有效吸收量也較高，直接導致了蛋白質含量的增加。因此，未來還需要進一步研究合理的施肥方法，在提升小站稻米食味品質的同時達到節肥增效的目的。

4 結論

本次研究對指導天津小站稻種植區土壤的施肥工作和進一步提升稻米品質具有較高的參考價值。研究區土壤有效S、Mn和B含量極為豐富，有效Cu和Ni含量較高，有效Mo含量為中等水平；有效S、Mn和Mo的變異程度較高，具有相對較強的空間異質性；有效S和Cu擬合效果最好的模型為指數模型，有效Mn、B和Ni為高斯模型，有效Mo為球狀模型；6種元素的塊基比值都＜25%，有效Cu和S變程值最大，有效Ni變程值最小，有效S、Mn、Mo和Ni的決定系數值最接近于1，擬合效果最理想；元素空間分布規律性較強，有效S、Cu和Mo呈北高南低的分布規律，有效Mn和Ni呈南高北低的分布規律，有效B呈東高西低的分布規律；土壤中S、Cu、B、Mo的有效量與稻米透明度、精米率、蛋白質、膠稠度和直鏈淀粉等品質指標具有顯著相關性。