基于灰色關(guān)聯(lián)-結(jié)構(gòu)方程模型的土壤酸化驅(qū)動因子研究*

周碧青 邱龍霞 張黎明 張 秀 陳成榕 邢世和，?

（1 福建農(nóng)林大學(xué)資源與環(huán)境學(xué)院，福州 350002）

（2 土壤生態(tài)系統(tǒng)健康與調(diào)控福建省高校重點(diǎn)實(shí)驗(yàn)室，福州 350002）

（3 School of Environment and Nature Science，Griffith University，Nathan 4111，Australia）

土壤酸堿性是反映耕地質(zhì)量的重要指標(biāo)之一，不同植物生長發(fā)育適宜的土壤酸堿性不同，如棉花、甘藍(lán)、大麥、玉米、小麥等作物適宜在中性至微堿性土壤環(huán)境下生長；水稻、油菜、紫云英、花生等適宜偏酸性土壤環(huán)境生長；而煙草、馬鈴薯等則適宜在酸性土壤環(huán)境下生長［1］。土壤酸化是耕地質(zhì)量退化的一個重要方面，其實(shí)質(zhì)是自然和人為因素共同作用導(dǎo)致土壤中鹽基離子減少，H+和Al3+增加，土壤鹽基飽和度下降，氫飽和度增加的過程［2］。土壤酸化打破了原有的土壤生態(tài)平衡，導(dǎo)致土壤理化和生化性質(zhì)改變，結(jié)構(gòu)性變差，微生物數(shù)量下降，礦質(zhì)養(yǎng)分失衡，重金屬等有毒元素活化，土壤肥力降低，進(jìn)而對作物生長、品質(zhì)和產(chǎn)量產(chǎn)生不良影響［3-4］，揭示區(qū)域耕地土壤酸化的主要驅(qū)動因素，對于科學(xué)阻控區(qū)域耕地土壤酸化、提升區(qū)域耕地質(zhì)量具有十分重要的理論和現(xiàn)實(shí)意義。土壤酸化及其防控已成為國內(nèi)外廣泛關(guān)注的研究熱點(diǎn)，國內(nèi)外至今針對區(qū)域土壤pH變化與空間分異等方面的研究報(bào)道很多且較為深入［5-10］，而有關(guān)土壤酸化影響因素的研究則主要采用常規(guī)統(tǒng)計(jì)分析或線性相關(guān)分析探討單一因子與土壤酸化的關(guān)系［11-15］，然而區(qū)域土壤酸化是諸多外在因素（如降水、施肥、酸雨等）和內(nèi)在因素（如有機(jī)質(zhì)、質(zhì)地、陽離子交換量等）綜合作用的結(jié)果，不同因素對土壤酸化的影響路徑和影響程度均各異，因此現(xiàn)有的相關(guān)研究無法清晰地揭示區(qū)域土壤酸化的多因素綜合影響，也無法闡明區(qū)域土壤酸化影響因素的主次關(guān)系和影響程度。為此，本研究以福建省為研究區(qū)域，在利用1:5萬省域耕地土壤類型空間數(shù)據(jù)庫、1982年36 777個和2016年56 445個耕地表層（0～20 cm）調(diào)查樣點(diǎn)土壤屬性數(shù)據(jù)以及氣象站點(diǎn)相關(guān)氣候要素、酸雨監(jiān)測點(diǎn)降水pH和化肥施用量等數(shù)據(jù)建立省域耕地土壤酸化及其可能影響因素空間數(shù)據(jù)庫的基礎(chǔ)上，以1982—2016年酸化耕地土壤pH變化量為母序列，以土壤酸化的可能影響因素為子序列，借助灰色斜率關(guān)聯(lián)和結(jié)構(gòu)方程分析模型，深入探討1982—2016年間福建省耕地土壤酸化的關(guān)鍵驅(qū)動因素及其影響程度大小，旨在為區(qū)域耕地土壤酸化調(diào)控和土壤培肥保育提供科學(xué)依據(jù)。

1 材料與方法

1.1 研究區(qū)概況

福建省地處我國東南沿海，介于23°33′～28°20′N，115°50′～120°40′E；屬亞熱帶海洋性季風(fēng)氣候，年均氣溫14.6～21.3℃，大于等于10℃年活動積溫5 000～7 800℃，年均降水量1 037～2 051 mm；境內(nèi)地形地貌復(fù)雜多樣，以丘陵山地為主，平原、盆地兼之，水熱條件垂直分異明顯［16］。2016年全省耕地總面積13 363 km2，其中水田、旱地和水澆地面積分別占83.00%、13.77%和3.23%，水田主要分布于寧德、福州、南平、三明、漳州和龍巖市，合計(jì)占水田總面積的85.98%；旱地主要分布于寧德、福州、莆田、泉州和漳州市，合計(jì)占旱地總面積的76.03%；水澆地主要分布于福州、莆田、泉州、廈門和漳州市，合計(jì)占水澆地總面積的88.15%。全省耕地土壤共分為9個土類（水稻土、濱海鹽土、潮土、紅壤、赤紅壤、黃壤、風(fēng)沙土、石灰土和紫色土）、15個亞類、31個土屬和81個土種［16］。

1.2 數(shù)據(jù)來源

從省農(nóng)業(yè)廳、國土廳、環(huán)保廳、統(tǒng)計(jì)局和氣象站等部門收集研究所需的福建省耕地、土壤等相關(guān)空間屬性數(shù)據(jù)庫及施肥量資料，包括：（1）2016年耕地利用-土壤類型數(shù)據(jù)庫；（2）數(shù)字高程模型數(shù)據(jù)庫（DEM）；（3）各縣（市、區(qū)）1982年36 777個和2016年56 445個耕地表層土壤調(diào)查樣點(diǎn)數(shù)據(jù)庫（圖1a、圖1b）；（4）氣象站點(diǎn)1982—2016年均溫度和降水量數(shù)據(jù)庫；（5）環(huán)境監(jiān)測點(diǎn)1996—2016降水年均pH數(shù)據(jù)庫；（6）各縣（市、區(qū)）1982—2016年農(nóng)作物年均單位面積施肥量數(shù)據(jù)資料。上述數(shù)據(jù)庫資料比例尺均為1:5萬、高斯-克呂格投影系、西安80坐標(biāo)系和1956年黃海高程系。

圖1 研究區(qū)1982年（a）和2016年（b）耕地土壤調(diào)查樣點(diǎn)分布圖Fig. 1 Distribution of sampling sites in the 1982 (a) and 2016 (b) surveys of cropland soils in Fujian Province

1.3 方法步驟

根據(jù)土壤酸化理論，初選福建省耕地土壤酸化可能影響因素包括年均溫度、年均降水量、海拔、坡度、有機(jī)質(zhì)、CEC、堿解氮、有效磷、砂粒、粉粒、黏粒、降水年均pH和年均單位面積施肥量等。利用上述收集的相關(guān)數(shù)據(jù)庫等資料，借助ArcGIS軟件，采用普通克里格插值法［17］進(jìn)行土壤pH、有機(jī)質(zhì)、CEC、堿解氮、有效磷、砂粒、粉粒、黏粒和年均降水pH空間屬性數(shù)據(jù)的點(diǎn)面推算，采用三維趨勢面分析結(jié)合反距離權(quán)重插值殘差訂正法［18］進(jìn)行年均溫度和年均降水量空間屬性數(shù)據(jù)的點(diǎn)面推算，分別生成福建省1982年和2016年土壤pH以及相關(guān)屬性2 mm×2 mm柵格數(shù)據(jù)庫，采用數(shù)字高程模型（Digital Elevation Model，DEM）數(shù)據(jù)庫生成福建省海拔和坡度2 mm×2 mm柵格數(shù)據(jù)庫。以2016年福建省耕地土種類型矢量圖層分別掩膜上述屬性柵格數(shù)據(jù)圖層，并采用面積加權(quán)平均法分別計(jì)算并建立福建省耕地土壤酸化可能影響因素空間屬性數(shù)據(jù)庫，采用人機(jī)對話賦值法建立福建省耕地單元年均單位面積施肥量空間屬性數(shù)據(jù)庫。研究區(qū)耕地土壤pH及相關(guān)屬性數(shù)據(jù)變化范圍見表1。

表1 研究區(qū)耕地土壤pH及相關(guān)屬性數(shù)據(jù)Table 1 Data of pH and relative properties of cropland soil in Fujian Province

借助ArcGIS軟件，將上述建立的2016年和1982年福建省耕地土壤pH矢量數(shù)據(jù)圖層進(jìn)行減法運(yùn)算，生成各評價(jià)單元pH變化量（△pH）空間數(shù)據(jù)圖層，將2016年土壤pH＜6.5且△pH＜0的耕地土壤歸為酸化土壤，建立1982—2016年福建省耕地土壤酸化空間數(shù)據(jù)庫。在利用灰色斜率關(guān)聯(lián)分析模型（Grey Slope Correlation Model，GSCM）確定耕地土壤酸化主要驅(qū)動因素的基礎(chǔ)上，借助結(jié)構(gòu)方程分析模型（Structural Equation Model，SEM）探討福建省耕地土壤酸化主要驅(qū)動因素的影響路徑和影響程度。

GSCM以關(guān)聯(lián)度系數(shù)為衡量指標(biāo)評判因素之間關(guān)系的密切程度，關(guān)聯(lián)系數(shù)越高表明兩者之間關(guān)系越密切［19］。借助ArcGIS軟件，從福建省1982—2016年耕地土壤酸化及其可能影響因素空間數(shù)據(jù)庫中分別導(dǎo)出各酸化耕地土壤單元的△pH及其相應(yīng)單元的可能影響因素值，分別建立母序列和子序列數(shù)列。將母序列和子序列數(shù)據(jù)導(dǎo)入DPS V15.10軟件中，借助GSCM分析模塊計(jì)算酸化耕地土壤△pH與可能影響因素的關(guān)聯(lián)度系數(shù)，根據(jù)關(guān)聯(lián)系數(shù)高低確定研究區(qū)耕地土壤酸化的主要驅(qū)動因素。

SEM有效地整合了因子分析、回歸分析、路徑分析等多種方法，通過假設(shè)影響路徑，對復(fù)雜數(shù)據(jù)進(jìn)行有效處理，進(jìn)而揭示其因果關(guān)系［20］。在進(jìn)行SEM擬合分析之前，必須驗(yàn)證供試耕地土壤△pH及其主要驅(qū)動因素?cái)?shù)據(jù)的可信度［21］，本研究借助Spss22.0軟件的信度分析模塊，對耕地土壤△pH母序列與主要驅(qū)動因素子序列數(shù)據(jù)進(jìn)行信度分析，結(jié)果表明這些數(shù)據(jù)的信度值為0.783，大于信度理想值標(biāo)準(zhǔn)（0.700）［19］，表明這些供試數(shù)據(jù)具有較好的一致性，可以直接采用SEM進(jìn)行擬合分析。根據(jù)土壤酸化理論，可以假設(shè)年均單位面積施肥量、CEC、黏粒、年均降水量、降水年均pH、有機(jī)質(zhì)等主要驅(qū)動因素對耕地土壤酸化影響包括以下路徑：（1）CEC、降水年均pH、年均單位面積施肥量、黏粒、年均降水量、有機(jī)質(zhì)對土壤△pH有直接影響；（2） 年均降水量通過影響土壤CEC、有機(jī)質(zhì)、黏粒和年均降水pH而間接影響土壤△pH；（3） 黏粒通過影響土壤CEC和有機(jī)質(zhì)而間接影響土壤△pH；（4） 有機(jī)質(zhì)通過影響土壤CEC而間接影響土壤△pH；（5）降水年均pH通過影響土壤有機(jī)質(zhì)、CEC和黏粒而間接影響土壤△pH；（6）年均單位面積施肥量通過影響土壤有機(jī)質(zhì)、CEC、黏粒和降水年均pH而間接影響土壤△pH。根據(jù)上述假設(shè)的6個影響路徑，建立研究區(qū)耕地土壤酸化主要驅(qū)動因素平衡關(guān)系的初始模型，借助Amos21.0的SEM模塊進(jìn)行擬合分析，獲得耕地土壤酸化主要驅(qū)動因素的標(biāo)準(zhǔn)化修正模型以及各顯變量之間的影響路徑及其效應(yīng)，分析探討福建省耕地土壤酸化主要驅(qū)動因素的影響路徑與影響程度。

2 結(jié) 果

2.1 耕地土壤酸化與空間差異性

研究結(jié)果表明（表2），1982—2016年的34年間，福建省耕地土壤pH平均降幅達(dá)0.34個單位，其中pH降低0～0.3、0.3～1和大于等于1個單位的耕地土壤分別占耕地總面積的33.45%、31.79%和5.43%；全省強(qiáng)酸性、酸性耕地土壤面積分別增加4.26%和2.00%，微酸性和中性耕地土壤面積分別減少5.46%和1.48%；全省有75.30%的中性耕地土壤轉(zhuǎn)為微酸性、酸性和強(qiáng)酸性，55.72%的微酸性耕地土壤轉(zhuǎn)為酸性和強(qiáng)酸性；耕地土壤pH降幅均值高于全省平均水平的有福州、龍巖、莆田、泉州和廈門市，酸化耕地面積較大的地級市包括福州、龍巖、南平、寧德、泉州和漳州市，合計(jì)占全省酸化耕地總面積的84.70%?？梢姡?982—2016年福建省有70.67%的耕地土壤發(fā)生不同程度酸化，且土壤酸化的區(qū)域差異較為明顯。

表2 1982—2016年福建省耕地土壤△pH均值和酸化面積Table 2 Mean of △pH of cropland soil and area of soil acidification in Fujian Province during 1982 and 2016

2.2 耕地土壤酸化的主要驅(qū)動因素

以土壤pH變化量（△pH）為母序列，將可能的影響因子作為子序列進(jìn)行GSCM分析，結(jié)果表明（表3），年均單位面積施肥量、CEC、黏粒與△p H的關(guān)聯(lián)程度最高，關(guān)聯(lián)系數(shù)介于0.884～0.954，其次為年均降水量、降水年均pH、有機(jī)質(zhì)、堿解氮及有效磷，與△pH的關(guān)聯(lián)系數(shù)絕對值|R|介于0.609～0.669，而砂粒、年均溫度、坡度及粉粒與△pH的關(guān)聯(lián)程度相對較低，關(guān)聯(lián)系數(shù)絕對值|R|介于0.561～0.596。其中年均單位面積施肥量、年均降水量、堿解氮、有效磷、砂粒、年均溫度和坡度等對△pH呈負(fù)向影響，即上述因素?cái)?shù)值越大，耕地土壤酸化越嚴(yán)重；而CEC、黏粒、降水年均pH、有機(jī)質(zhì)和粉粒等對△pH呈正向影響，即上述因素?cái)?shù)值越大，耕地土壤酸化越弱。由于耕地土壤堿解氮和有效磷含量與年均單位面積施肥量存在密切關(guān)系，故其對耕地土壤酸化的影響可以通過年均單位面積施肥量來反映。因此，可以確定關(guān)聯(lián)系數(shù)絕對值|R|＞0.620的年均單位面積施肥量、CEC、黏粒、年均降水量、降水年均pH、有機(jī)質(zhì)等6個因子為福建省耕地土壤酸化的主要驅(qū)動因子。

2.3 耕地土壤酸化主要驅(qū)動因素的影響路徑

利用Amos 21.0軟件對建立的研究區(qū)耕地土壤酸化主要驅(qū)動因素平衡關(guān)系的初始模型進(jìn)行擬合，結(jié)果表明年均單位面積施肥量和降水年均pH對黏粒影響的P分別為0.06和0.08，表明這兩條假設(shè)影響路徑系數(shù)為零的概率達(dá)到了6%和8%，大于5%的顯著性水平，即其顯著性檢驗(yàn)未通過，說明初始模型擬合指數(shù)較差，需要對模型進(jìn)行修正。通過對這兩條顯著性較差的假設(shè)影響路徑進(jìn)行釋放，然后對修正后的模型進(jìn)行第二次擬合，獲得平衡關(guān)系的標(biāo)準(zhǔn)化修正模型及其路徑系數(shù)見圖2。第二次擬合結(jié)果各驅(qū)動因素影響路徑系數(shù)為零的概率均小于5%，表明修正模型的擬合效果達(dá)到顯著性水平，其擬合指數(shù)見表4。從表4可見，修正模型的x2/df為0.569，相對應(yīng)的P為0.451（＞0.05），其他擬合指數(shù)也均符合評判標(biāo)準(zhǔn)，故認(rèn)為該修正模型是合理的。

從修正模型的路徑系數(shù)可見（圖2），年均單位面積施肥量、年均降水量、降水年均pH、黏粒、有機(jī)質(zhì)和CEC對土壤△pH有直接影響，路徑系數(shù)分別為-0.70、-0.40、0.38、0.30、0.22和0.23；年均降水量通過影響土壤CEC、有機(jī)質(zhì)和黏粒而間接影響土壤△pH，路徑系數(shù)分別為-0.29、-0.46和-0.29；黏粒通過影響有機(jī)質(zhì)和CEC而間接影響土壤△pH，路徑系數(shù)分別為0.34和0.23；有機(jī)質(zhì)通過影響CEC而間接影響土壤△pH，路徑系數(shù)為0.50；降水年均pH通過影響土壤有機(jī)質(zhì)和CEC而間接影響土壤△pH，路徑系數(shù)分別為-0.40和-0.23；土壤黏粒含量通過影響土壤有機(jī)質(zhì)和CEC而間接影響土壤△pH，路徑系數(shù)分別為-0.40和-0.23；年均單位面積施肥量通過影響土壤有機(jī)質(zhì)、CEC和降水年均pH而間接影響土壤的△pH，路徑系數(shù)分別為0.20、0.33和-0.27。根據(jù)結(jié)構(gòu)方程模型路徑系數(shù)分析的原則，在有顯著性關(guān)系的前提下，路徑系數(shù)絕對值大的數(shù)列代表外生顯變量對內(nèi)生顯變量的影響程度較大，反之則小［21］?？梢?，年均單位面積施肥量、年均降水量、降水年均pH和黏粒對福建省耕地土壤酸化具有較顯著的直接影響；而年均降水量和黏粒通過影響有機(jī)質(zhì)、有機(jī)質(zhì)通過影響CEC、降水年均pH通過影響有機(jī)質(zhì)以及年均單位面積施肥量通過影響CEC而對福建省耕地土壤酸化具有較顯著的間接影響。

表3 福建省耕地土壤pH變化量（△ pH）與可能影響因子的GSCM分析Table 3 Grey slope correlation analysis between △ pH and potential impact factor of the cropland soil in Fujian

表4 平衡關(guān)系標(biāo)準(zhǔn)化修正模型擬合指數(shù)Table 4 Fitting indices of the standardized modification model for balanced relation

2.4 耕地土壤酸化主要驅(qū)動因素的影響效應(yīng)

圖2 福建省耕地土壤酸化主要驅(qū)動因子平衡關(guān)系的標(biāo)準(zhǔn)化修正模型Fig. 2 Standardized modification model for balanced relations between key impact factors of the cropland soil acidification in Fujian

結(jié)構(gòu)方程模型分析獲得的耕地土壤酸化主要驅(qū)動因素的影響效應(yīng)結(jié)果表明，年均單位面積施肥量對△pH的直接效應(yīng)最為顯著，路徑系數(shù)絕對值高達(dá)0.70，同時(shí)，大量施用化肥還會通過影響土壤有機(jī)質(zhì)、CEC及降水pH等因素間接影響土壤△pH，其間接效應(yīng)系數(shù)分別為0.20、0.33和-0.27，年均單位面積施肥量的總效應(yīng)絕對值高達(dá)0.86，故長期持續(xù)大量施用尿素、氯化銨、硫酸銨和過磷酸鈣等酸性或生理酸性肥料是造成福建省耕地土壤酸化的最主要人為因素。年均降水量對△pH的直接效應(yīng)為-0.40，間接效應(yīng)為-0.35，總效應(yīng)絕對值也高達(dá)0.75，故多雨的氣候條件是導(dǎo)致1982—2016年福建省耕地土壤酸化的最主要自然因素。降水年均pH對△pH的直接效應(yīng)為0.38，間接效應(yīng)為0.11，總效應(yīng)達(dá)到0.49，故酸雨是導(dǎo)致福建省1982—2016年耕地土壤酸化的另一最主要人為因素。表7結(jié)果還表明，福建省耕地土壤△pH與黏粒、有機(jī)質(zhì)和CEC的直接效應(yīng)分別為0.30、0.22和0.23，間接效應(yīng)分別為0.13、0.18和0，總效應(yīng)分別為0.43、0.30和0.23，故黏粒、有機(jī)質(zhì)和CEC是減緩2008—2016年福建省耕地土壤酸化的主要內(nèi)在因子。

3 討 論

3.1 關(guān)鍵外在驅(qū)動因素對耕地土壤酸化的影響

Guo等［9］研究指出大量施用化肥加速中國農(nóng)田土壤酸化進(jìn)程，使過去20年來我國主要農(nóng)田土壤pH均值下降約0.5個單位，相當(dāng)于土壤酸量（H+）在原有基礎(chǔ)上增加2.2倍。長期以來，福建省單位面積化肥施用量一直位居全國前列，近31年來全省年化肥施用量持續(xù)增長，2016年化肥施用量為1985年的2.52倍，并以氮、磷肥占優(yōu)勢，分別占年化肥施用量的61%和17%，且以尿素、氯化銨、硫酸銨和過磷酸鈣等生理酸性或酸性化肥為主。耕地土壤酸化面積較大的南平、泉州、龍巖和漳州等市年均單位面積施肥量介于508～1 980 kg·hm-2，是全省年均單位面積施肥量（434 kg·hm-2）的1.17倍～4.56倍［22］，長期不合理施用酸性或生理酸性肥料，或直接向土壤輸入磷酸而導(dǎo)致土壤酸化，或因植物喜好吸收而使SO4-2在土壤中殘留并與作物代換吸收釋放出的H+結(jié)合形成硫酸而導(dǎo)致土壤酸化［23］，或因在土壤中發(fā)生硝化作用釋放H+而加速土壤酸化［24］，故年均單位面積施肥量必然對福建省耕地土壤酸化產(chǎn)生顯著的直接影響，其直接影響效應(yīng)高達(dá)-0.70。此外，有研究表明偏施化肥而忽視有機(jī)肥施用，造成土壤有機(jī)質(zhì)含量下降以及鹽基離子補(bǔ)充不足和淋失數(shù)量增加［25］，并顯著提高土壤鋁、鐵的活性及其含量［26］，進(jìn)而加劇土壤酸化；不合理施用氮肥所致的氮素反硝化作用，致使土壤向大氣中排放的含氮化合物增加，大氣氮沉降也成比例提高，也會加速土壤酸化［12］，故大量施用化肥還會通過影響土壤有機(jī)質(zhì)、CEC及降水pH等因素而間接影響土壤pH，其間接影響效應(yīng)分別達(dá)0.20、0.33和-0.27，致使年均單位面積施肥量在耕地土壤酸化主要驅(qū)動因素中位居首位，總影響效應(yīng)高達(dá)0.86。因此，長期持續(xù)大量施用尿素、氯化銨、硫酸銨和過磷酸鈣等酸性或生理酸性肥料是福建省耕地土壤酸化首要的外在驅(qū)動因素。

福建省地處亞熱帶氣候區(qū)，年均氣溫15.8～21.7 ℃，最熱月均溫28 ℃，年均降水量高達(dá)1 015～1 923 mm，大于等于10 ℃積溫高達(dá)5 000～7 800 ℃［16］。高溫多雨的氣候條件，致使全省土壤形成和發(fā)育過程脫硅富鋁化作用和有機(jī)質(zhì)礦化作用強(qiáng)烈，鹽基物質(zhì)大量淋失，有機(jī)質(zhì)含量總體不高，陽離子代換量和鹽基飽和度降低，氫飽和度上升；此外，多雨的氣候條件易引發(fā)水土流失而使黏粒大量淋失，從而降低土壤對酸的緩沖性，故多雨的氣候條件必然對耕地土壤酸化產(chǎn)生較顯著的直接影響，并通過影響土壤有機(jī)質(zhì)、CEC和黏粒而對耕地土壤酸化產(chǎn)生間接影響。本研究結(jié)果表明，福建省年均降水量對耕地土壤△pH的直接影響效應(yīng)為-0.40，間接影響效應(yīng)為-0.35，總影響效應(yīng)絕對值高達(dá)0.75，在耕地土壤酸化主要驅(qū)動因素中位居第二，致使降水量較大的閩北地區(qū)（南平市）發(fā)生酸化的耕地土壤面積較大，占全省酸化耕地土壤總面積的20.41%。因此，多雨的氣候條件是福建省耕地土壤發(fā)生不同程度酸化主要的外在驅(qū)動因素。

人類生活、生產(chǎn)活動以及土壤氮素反硝化作用產(chǎn)生的SO2、NOx等酸性氣體，在大氣中經(jīng)過一系列化學(xué)反應(yīng)形成H2SO4和HNO3，這些強(qiáng)酸隨著降水（即酸雨）落到地表或滲入土壤，直接向土壤輸入H+而加速土壤酸化［27］。福建省屬于全國的酸雨敏感區(qū)，酸雨較嚴(yán)重的區(qū)域主要位于閩北、閩西和閩南，其中龍巖、泉州等城市的酸雨頻率均＞50%［28］。根據(jù)酸雨監(jiān)測點(diǎn)資料，南平、泉州、漳州和福州等市的降水年均pH分別為5.2、5.1、4.6和4.7，較全省降水年均pH（6.0）低0.9～1.4個單位，故較強(qiáng)的酸雨必然對這些區(qū)域耕地土壤酸化產(chǎn)生較顯著的直接影響，表現(xiàn)為福建省降水年均pH對耕地土壤△ pH的直接影響效應(yīng)達(dá)0.38，致使南平、龍巖、泉州、漳州和福州等市的耕地土壤發(fā)生酸化的耕地面積較大，合計(jì)占全省酸化耕地總面積的74.07%。有研究表明，模擬酸雨對有機(jī)碳礦化有激發(fā)作用而影響土壤有機(jī)質(zhì)含量［29］，且明顯促進(jìn)土壤礦物風(fēng)化作用而影響土壤CEC［30］，故酸雨通過影響土壤有機(jī)質(zhì)和CEC可間接影響土壤酸化。本研究結(jié)果表明，福建省降水年均pH通過影響有機(jī)質(zhì)和CEC而間接影響土壤△pH的間接影響效應(yīng)分別為-0.40和-0.23，致使降水年均pH對全省耕地土壤酸化的影響程度次于年均降水量而位居第三，總影響效應(yīng)達(dá)0.49。因此，酸雨也必然成為福建省耕地土壤酸化的另一主要外在驅(qū)動因素。

3.2 關(guān)鍵內(nèi)在驅(qū)動因素對耕地土壤酸化的影響

土壤抵抗酸堿變化的能力與土壤膠體數(shù)量、組成及CEC密切相關(guān)［1］。有機(jī)質(zhì)、黏粒和CEC高的耕地土壤對酸緩沖能力強(qiáng)，土壤越不易酸化，反之，則越易發(fā)生酸化。此外，土壤黏粒高低制約著土壤通氣性和微生物活動，進(jìn)而影響土壤有機(jī)質(zhì)積累，而土壤CEC主要取決于黏粒和有機(jī)質(zhì)含量與組成［1］，故土壤黏粒、有機(jī)質(zhì)和CEC必然對耕地土壤酸化產(chǎn)生較顯著的直接影響，且黏粒通過影響有機(jī)質(zhì)、有機(jī)質(zhì)通過影響CEC也必然對耕地土壤酸化產(chǎn)生間接影響。本研究結(jié)果表明，福建省黏粒、有機(jī)質(zhì)和CEC對耕地土壤△ pH的直接影響效應(yīng)分別為0.30、0.22和0.23，黏粒和有機(jī)質(zhì)對耕地土壤△ pH的間接影響效應(yīng)分別為0.13和0.18，黏粒、有機(jī)質(zhì)和CEC對耕地土壤酸化總影響效應(yīng)分別為0.43、0.30和0.23，致使福建省耕地土壤發(fā)生酸化程度的高低因土壤黏粒、有機(jī)質(zhì)和CEC不同而差異明顯，如位于龍巖市西北部（長汀縣）的耕地土壤CEC、黏粒和有機(jī)質(zhì)含量均較低，分別為7.88 cmol·kg-1、14.09%和28.02 g·kg-1，致使該區(qū)域耕地土壤酸化較為嚴(yán)重，而位于三明市北部（沙縣）的耕地土壤CEC、黏粒和有機(jī)質(zhì)含量均較高，分別為9.23 cmol·kg-1、15.64%和30.81g·kg-1，致使該區(qū)域耕地土壤酸化程度則相對較弱。因此，黏粒、有機(jī)質(zhì)和CEC成為減緩福建省耕地土壤酸化最主要的內(nèi)在驅(qū)動因素。

4 結(jié) 論

采用灰色斜率關(guān)聯(lián)分析和結(jié)構(gòu)方程分析模型相結(jié)合的方法，可以清晰地揭示區(qū)域耕地土壤酸化的多因素綜合影響，且可以深入闡明區(qū)域耕地土壤酸化主要驅(qū)動因素及其影響路徑和影響效應(yīng)。年均單位面積施肥量、年均降水量、降水年均pH、有機(jī)質(zhì)、黏粒和CEC對福建省耕地土壤酸化程度影響的總效應(yīng)分別為-0.86、-0.75、0.49、0.43、0.30和0.23?？梢?，福建省耕地土壤酸化是自然和人為多種因素交互作用所致，其中長期大量施用酸性或生理酸性化肥、亞熱帶豐沛的年降水量和較嚴(yán)重的酸雨是加速福建省耕地土壤酸化的三大主要外因，而黏粒、有機(jī)質(zhì)和CEC則是減緩福建省耕地土壤酸化的三大主要內(nèi)因。