黃三角瀕海區(qū)土壤鹽漬化的驅(qū)動(dòng)力分析及預(yù)測(cè)模擬

趙銘，常春艷，王卓然，趙庚星

（山東農(nóng)業(yè)大學(xué)資源與環(huán)境學(xué)院土肥高效利用國家工程研究中心，山東泰安 271000）

0 引言

土壤鹽漬化是一全球性環(huán)境問題，易導(dǎo)致農(nóng)田退化、土地生產(chǎn)力降低等問題。中國鹽漬土面積約3600萬hm2，占全國耕地面積的6.62%[1]。黃河三角洲自然資源豐富，被譽(yù)為最有開發(fā)潛力的三角洲，但區(qū)域內(nèi)生態(tài)環(huán)境脆弱[2]，土壤鹽漬化現(xiàn)象較為嚴(yán)重，由于長期海水的浸漬作用以及黃河水動(dòng)力的影響，導(dǎo)致地下水埋深較淺、礦化度較高，土壤積鹽問題嚴(yán)重[3]，同時(shí)人類對(duì)土地的不合理開發(fā)利用，使得區(qū)域土壤鹽漬化程度日益加劇，制約著當(dāng)?shù)亟?jīng)濟(jì)社會(huì)的可持續(xù)發(fā)展。

土壤鹽漬化的產(chǎn)生是多種因素綜合作用的結(jié)果，在對(duì)其探究時(shí)，既要理清土壤鹽漬化的主要因素，也要對(duì)不同因素間的作用關(guān)系進(jìn)行分析。蘇春利等[4]探究了地下水因素是河套灌區(qū)土壤鹽漬化產(chǎn)生的主控因素；臧亮等[5]探知了黃驊市土壤含鹽量的主要影響因素為地下水礦化度、地下水埋深和高程。而地理探測(cè)器方法可以很好的用于土壤鹽分影響因素之間的綜合性分析，此方法可以通過探究不同變量之間的空間分異性進(jìn)行背后驅(qū)動(dòng)力的揭示[6-7]，主要包括因子探測(cè)器、風(fēng)險(xiǎn)探測(cè)器、生態(tài)探測(cè)器、交互作用探測(cè)器，目前地理探測(cè)器在不同研究方向中都被廣泛應(yīng)用，如社會(huì)科學(xué)[8]、自然科學(xué)[9]、環(huán)境科學(xué)[10-11]等。

在對(duì)土壤鹽分進(jìn)行預(yù)測(cè)模型的構(gòu)建時(shí)，主要是基于統(tǒng)計(jì)分析方法和機(jī)器學(xué)習(xí)方法來進(jìn)行分析[12-13]，常用的模型主要有多元線性回歸模型、偏最小二乘回歸模型、BP 神經(jīng)網(wǎng)絡(luò)、支持向量機(jī)、隨機(jī)森林[14-16]。如馬國林等[17]基于構(gòu)建的隨機(jī)森林、極限學(xué)習(xí)機(jī)、偏最小二乘回歸模型進(jìn)行了新疆干旱地區(qū)土壤鹽漬化監(jiān)測(cè)，并選出最優(yōu)模型為RF 模型。張智韜等[18]分析了支持向量機(jī)、BP 神經(jīng)模型、隨機(jī)森林模型和多元線性模型的精度效果，得出結(jié)論RF模型可作為土壤鹽分的最佳反演模型。在土壤鹽漬化的模擬研究中，多數(shù)通過水鹽運(yùn)動(dòng)進(jìn)行數(shù)值模擬[19-20]，或者通過模擬鹽分的影響因素變化程度進(jìn)行分析[21]。如黃亞捷等[22]利用Sahysmod 模型模擬未來10 年的土壤水鹽動(dòng)態(tài)變化情況。朱文東[23]通過模擬不同地下水埋深和地下水礦化度的情形來分析土壤鹽漬化的水鹽運(yùn)移規(guī)律。賈春青[24]模擬了地下水位在1.0～2.1 m之間的變化幅度，分析了在不同深度的埋深下土壤鹽分的動(dòng)態(tài)變化情況。由于不同地區(qū)土壤環(huán)境及鹽漬化狀況的差異，土壤鹽分驅(qū)動(dòng)力及其最佳模型并未有很好的兼容性，特別瀕海區(qū)域尚需進(jìn)一步的研究探索。

因此，本研究選取黃河三角洲墾利區(qū)為研究區(qū)，篩選出土壤鹽分主要驅(qū)動(dòng)因子，構(gòu)建鹽分的最佳預(yù)測(cè)模型，進(jìn)而對(duì)土壤鹽分進(jìn)行情景模擬分析，為黃河三角洲瀕海區(qū)土壤鹽漬化的治理和開發(fā)利用提供科學(xué)依據(jù)。

1 材料與方法

1.1 研究區(qū)概況

墾利區(qū)(118°15′—119°19'E，37°24′—38°10′N)隸屬于山東省東營市（圖1），位于黃河入?？冢瑬|瀕渤海，區(qū)域?qū)儆谂瘻貛Т箨懶约撅L(fēng)氣候，四季分明，春季干旱多風(fēng)，夏季炎熱多雨，秋季涼爽溫和，冬季寒冷干燥。地勢(shì)自西南到東北呈扇狀分布，地貌類型主要有微斜平地、海灘地、河灘高低，植被有天然植被和人工植被，天然植被多為耐鹽性植被。土地利用類型主要有耕地、林地、草地等，主要農(nóng)作物有小麥、玉米和棉花，土壤類型多為濱海鹽化潮土，土壤質(zhì)地以輕壤和砂壤為主，毛細(xì)管作用強(qiáng)烈，加上區(qū)域地下水埋深較淺且礦化度較高，同時(shí)受到海水浸漬的影響，土壤鹽漬化普遍。

圖1 研究區(qū)地理位置及采樣點(diǎn)分布

1.2 數(shù)據(jù)采集與處理

1.2.1 數(shù)據(jù)采集綜合考慮研究區(qū)地理位置、植被類型、土地利用方式等因素，在墾利區(qū)布設(shè)觀測(cè)點(diǎn)位47 個(gè)，如圖1所示，每個(gè)點(diǎn)位有多個(gè)觀測(cè)樣點(diǎn)，使用EC110便攜式鹽分計(jì)，對(duì)土壤表層(0～15 cm)電導(dǎo)率進(jìn)行測(cè)定，并采集土壤樣本，同時(shí)手持GPS 定位儀記錄實(shí)地坐標(biāo)，一共獲得108組樣點(diǎn)數(shù)據(jù)。

1.2.2 數(shù)據(jù)處理對(duì)土壤樣本進(jìn)行室內(nèi)自然風(fēng)干、研磨過篩，采用烘干法測(cè)定土壤鹽分含量，建立土壤電導(dǎo)率Eco(dS/m)和土壤含鹽量St(g/kg)之間關(guān)系模型(1)[25]。利用模型進(jìn)行電導(dǎo)率轉(zhuǎn)換，最終得到所有樣點(diǎn)土壤含鹽量數(shù)據(jù)。

1.3 研究方法

1.3.1 土壤鹽漬化程度分級(jí)根據(jù)土壤含鹽量(g/kg)的高低，將研究區(qū)土壤鹽漬化程度分為非鹽漬化土(0～1 g/kg)；輕度鹽漬化土(1～2 g/kg)；中度鹽漬化土(2～4 g/kg)；重度鹽漬化土(4～10 g/kg)；鹽土(＞10 g/kg)。

1.3.2 土壤鹽分驅(qū)動(dòng)因子的篩選基于前期研究與實(shí)地調(diào)查，綜合選取了蒸發(fā)量、降水量、地下水埋深、地下水礦化度、土壤粘粒含量、相對(duì)高程、植被覆蓋度、距海遠(yuǎn)近8 個(gè)與土壤鹽分有關(guān)的影響因子，通過地理探測(cè)器方法分析土壤鹽分與其相關(guān)因子之間的關(guān)系，從而進(jìn)行驅(qū)動(dòng)力的篩選。

地理探測(cè)器，通常用于探測(cè)兩個(gè)變量之間可能會(huì)發(fā)生的因果關(guān)系，主要分為因子探測(cè)器、風(fēng)險(xiǎn)探測(cè)器、生態(tài)探測(cè)器和交互作用探測(cè)器四部分。本研究從因子探測(cè)和交互作用探測(cè)兩方面進(jìn)行分析，其中因子探測(cè)器用于分析不同的影響因子對(duì)土壤鹽分的解釋力大小，用q值來表示，值越大說明解釋能力越強(qiáng)；交互作用探測(cè)器用于分析不同影響因子之間的相互作用力，用q值來衡量，因子間的關(guān)系有非線性減弱、單因子非線性減弱、雙因子增強(qiáng)、獨(dú)立和非線性增強(qiáng)5種。

1.3.3 土壤鹽分預(yù)測(cè)模擬

（1）土壤鹽分預(yù)測(cè)模型構(gòu)建。在進(jìn)行模型的構(gòu)建時(shí)，將研究區(qū)的樣本數(shù)據(jù)分為兩組，隨機(jī)選取2/3樣本數(shù)據(jù)用于建模集，1/3 用于驗(yàn)證集，建模集與驗(yàn)證集數(shù)據(jù)的基本特征也保持一致，使二者都能反映整體的鹽分狀況。

以篩選出的土壤鹽分驅(qū)動(dòng)因子為輸入變量，以土壤鹽分為輸出變量，分別采用統(tǒng)計(jì)分析中的多元線性回歸(Multiple Linear Regression，MLR)、偏最小二乘回歸(Partial Least Squares Regression，PLSR)以及機(jī)器學(xué)習(xí)中的BP 神經(jīng)網(wǎng)絡(luò)(Back Propagation Neural Network，BPNN)、支持向量機(jī)(Support Vector Machine，SVM)4 種方法進(jìn)行模型構(gòu)建。通過決定系數(shù)R2(Coefficient of determination) 和RMSE(Root mean square error)驗(yàn)證模型的預(yù)測(cè)精度，選用精度最高的模型作為土壤鹽分最佳預(yù)測(cè)模型，其中R2表示模型的構(gòu)建與驗(yàn)證的穩(wěn)定性，其數(shù)值與1越相近，說明模型穩(wěn)定性越強(qiáng)；RMSE用來檢驗(yàn)?zāi)Ｐ偷念A(yù)估效果，其數(shù)值越小，說明模型精度越高。同時(shí)統(tǒng)計(jì)分析模型通過SPSS 19 軟件實(shí)現(xiàn)，機(jī)器學(xué)習(xí)模型通過MATLAB R2016a 軟件完成。

（2）土壤鹽分情景模擬。以篩選出的土壤鹽分驅(qū)動(dòng)因子為輸入變量，以土壤鹽分為輸出變量，通過選出的最佳預(yù)測(cè)模型進(jìn)行土壤鹽分的模擬預(yù)測(cè)，主要分為以下步驟：a.將模型應(yīng)用于研究區(qū)的點(diǎn)位數(shù)據(jù)，從而得到各點(diǎn)位的鹽分預(yù)測(cè)值；b.通過ARCGIS 軟件將土壤鹽分值進(jìn)行空間插值分析；c.將插值結(jié)果圖根據(jù)土壤鹽漬化的分級(jí)標(biāo)準(zhǔn)進(jìn)行重分類；d.統(tǒng)計(jì)不同鹽漬化的占比情況。

2 結(jié)果與分析

2.1 土壤鹽分空間分布特征

通過ArcGIS軟件對(duì)土壤鹽分進(jìn)行插值，根據(jù)鹽漬化程度等級(jí)進(jìn)行劃分，所得墾利區(qū)土壤鹽分空間分布圖（圖2），可以看出，整體上研究區(qū)的鹽漬化程度呈現(xiàn)出北高南低、東高西低的分布，東部沿海地區(qū)較易受到海水的入侵，導(dǎo)致鹽土在這些區(qū)域分布廣泛；重度鹽漬土主要分布在中東部地區(qū)；中度鹽漬土主要分布在西南和中部地區(qū)，這些區(qū)域種植農(nóng)作物多，較多的植被覆蓋使得鹽漬化程度較輕；研究區(qū)的輕度鹽漬土分布較少。

圖2 土壤鹽分空間分布圖

2.2 土壤鹽分驅(qū)動(dòng)因子分析

（1）因子探測(cè)器。表1是因子探測(cè)器的分析結(jié)果，通過顯著性檢驗(yàn)P值可以看出，蒸發(fā)量、降水量、地下水埋深、地下水礦化度、地表高程、植被覆蓋度、距海遠(yuǎn)近和粘粒含量8 個(gè)影響因子均通過了顯著性檢驗(yàn)，說明各影響因子均對(duì)土壤鹽分的分布和土壤鹽漬化的產(chǎn)生有著一定的控制作用。各因子對(duì)土壤鹽分的影響力大小為地下水礦化度(0.762)＞植被覆蓋度(0.737)＞地下水埋深(0.730)＞距海遠(yuǎn)近(0.627)＞粘粒含量(0.613)＞地表高程(0.609)＞降水量(0.419)＞蒸發(fā)量(0.364)，說明不同因子對(duì)土壤鹽漬化的影響存在差異性。因子探測(cè)器的結(jié)果表明，土壤鹽漬化是在蒸發(fā)量、降水量、地下水埋深、地下水礦化度、地表高程、植被覆蓋度、距海遠(yuǎn)近和粘粒含量這8 個(gè)因素的共同影響下發(fā)生的，其中地下水礦化度、植被覆蓋度、地下水埋深是土壤鹽漬化的主導(dǎo)因素；距海遠(yuǎn)近、粘粒含量和地表高程是土壤鹽漬化的重要因素；蒸發(fā)量和降水量對(duì)土壤鹽漬化的影響力較小。

表1 因子探測(cè)器結(jié)果

（2）交互作用探測(cè)器。表2 是交互作用探測(cè)器的分析結(jié)果，可以看出不同因子間的交互作用影響力均大于單一影響因子對(duì)土壤鹽漬化的影響力，均表現(xiàn)為雙因子增強(qiáng)，說明因子間的交互作用對(duì)土壤鹽漬化的影響力顯著增加。交互作用探測(cè)器的結(jié)果表明，地下水礦化度與植被覆蓋度交互作用的影響力最高，達(dá)到了94.2%，地下水礦化度與地下水埋深、距海遠(yuǎn)近交互作用的影響力也分別高達(dá)93.2%、93.0%，地下水礦化度與蒸發(fā)量、降水量、地表高程、距海遠(yuǎn)近和粘粒含量之間的交互作用力也均在70%以上，同時(shí)地下水礦化度、地下水埋深、植被覆蓋度和距海遠(yuǎn)近這4個(gè)因子之間的兩兩交互作用力也較強(qiáng)。由此可以看出地下水礦化度是土壤鹽漬化的最主要影響因子，植被覆蓋度、地下水埋深、距海遠(yuǎn)近是土壤鹽漬化的重要影響因子。

表2 交互作用探測(cè)器結(jié)果

結(jié)合因子探測(cè)器和交互作用探測(cè)器的分析結(jié)果可以看出，土壤鹽漬化主要受到地下水、植被和海水的影響，由此，將地下水礦化度、地下水埋深、植被覆蓋度和距海遠(yuǎn)近作為土壤鹽分的主要驅(qū)動(dòng)因子。

2.3 土壤鹽分預(yù)測(cè)模型

以地下水礦化度、地下水埋深、植被覆蓋度和距海遠(yuǎn)近為輸入數(shù)據(jù)，以土壤含鹽量為輸出數(shù)據(jù)，構(gòu)建了BPNN、SVM、MLR 和PLSR 模型，不同模型的預(yù)測(cè)精度結(jié)果如表3 所示，以驗(yàn)證集樣本數(shù)據(jù)對(duì)模型進(jìn)行驗(yàn)證的結(jié)果如圖3所示。由表3可以看出，4個(gè)模型的建模集R2均在0.6 以上，其中BPNN 模型的建模精度最高，R2為0.8847，RMSE為1.1350；SVM 模型次之，R2為0.8397，RMSE為1.1112，二者均表現(xiàn)出了較強(qiáng)的預(yù)測(cè)能力；PLSR模型的建模精度最差。BPNN模型的驗(yàn)證集R2為0.7999，RMSE為1.1204，模型表現(xiàn)出較強(qiáng)的穩(wěn)定性，依次是SVM模型、MLR模型、PLSR模型。從模型驗(yàn)證結(jié)果來看（圖3），BPNN 模型的驗(yàn)證效果與SVM 模型、MLR 模型和PLSR 模型相比結(jié)果最佳，其實(shí)際值與預(yù)測(cè)值的R2最高，為0.8263，其余3個(gè)模型的R2均在0.8以下。

表3 土壤鹽分預(yù)測(cè)模型評(píng)價(jià)指標(biāo)表

圖3 土壤鹽分預(yù)測(cè)模型驗(yàn)證結(jié)果

經(jīng)對(duì)比分析，BPNN 模型在4 個(gè)模型中的擬合程度最好，其穩(wěn)定性最強(qiáng)，模型精度最高，由此可選用BPNN模型作為土壤鹽分的最佳預(yù)測(cè)模型。

2.4 土壤鹽分預(yù)測(cè)結(jié)果

通過BP 神經(jīng)網(wǎng)絡(luò)最佳預(yù)測(cè)模型對(duì)墾利區(qū)表層土壤鹽分進(jìn)行預(yù)測(cè)，得到的預(yù)測(cè)結(jié)果與實(shí)測(cè)結(jié)果面積統(tǒng)計(jì)如表4 所示。輕度鹽漬土占比最少，實(shí)測(cè)面積比例為0.38%，預(yù)測(cè)結(jié)果面積比例為0.27%；中度鹽漬土實(shí)測(cè)面積比例為27.91%，預(yù)測(cè)結(jié)果面積比例為25.72%；重度鹽漬土實(shí)測(cè)面積比例為23.90%，預(yù)測(cè)結(jié)果面積比例為25.73%；鹽土占比最多，實(shí)測(cè)面積比例為47.81%，預(yù)測(cè)結(jié)果面積比例為48.29%。

表4 不同鹽漬化程度面積統(tǒng)計(jì) %

實(shí)測(cè)與預(yù)測(cè)結(jié)果面積統(tǒng)計(jì)的結(jié)果顯示，不同程度的鹽漬土面積比例差別較小，因此使用BP神經(jīng)網(wǎng)絡(luò)模型對(duì)土壤鹽分進(jìn)行預(yù)測(cè)分析是可行的，同時(shí)預(yù)測(cè)效果良好。

2.5 土壤鹽分情景模擬

通過上述土壤鹽分驅(qū)動(dòng)因子分析結(jié)果可以看出，地下水因素對(duì)土壤鹽漬化的產(chǎn)生起到關(guān)鍵性影響作用，因此本研究模擬地下水的變化情景，同時(shí)借助上述最佳預(yù)測(cè)模型-BP神經(jīng)網(wǎng)絡(luò)模型對(duì)土壤鹽漬化狀況進(jìn)行分析討論。

在地下水的變化情景下，將水位變化幅度限定為0.5 m，分別設(shè)定對(duì)照組、地下水位下降0.5 m、地下水位上升0.5 m 3種情景。

通過地下水位變化情景下的土壤鹽漬化程度以及面積統(tǒng)計(jì)可知（圖4，表5），當(dāng)?shù)叵滤幌陆禃r(shí)，輕度鹽漬土在研究區(qū)分布最少，其變化最小，面積比例較對(duì)照組增加了0.22%，由0.39%變化到0.61%；中度鹽漬土在中東部地區(qū)變化顯著，同時(shí)呈現(xiàn)減少的趨勢(shì)，面積比例由27.76%減少到22.30%；中部地區(qū)黃河兩岸多為重度鹽漬化，西南地區(qū)的部分鹽土也轉(zhuǎn)變?yōu)橹囟塞}漬土，重度鹽漬土的面積比例由23.70%變化到了39.98%；鹽土仍集中分布在沿海地區(qū)，同時(shí)部分鹽土轉(zhuǎn)換到重度鹽漬土水平，面積比例由48.15%減少到38.11%。

表5 地下水位變化情景下的土壤鹽漬化面積統(tǒng)計(jì) %

圖4 地下水位變化情景下的土壤鹽漬化程度

當(dāng)?shù)叵滤簧仙龝r(shí)，輕度鹽漬土呈現(xiàn)微弱變化，面積比例從0.39%變化到了0.37%；中度鹽漬土變化顯著，面積比例從27.76%減少到了12.99%，重度鹽漬土面積比例由23.70%增加到46.21%，二者主要體現(xiàn)在西南地區(qū)和中部地區(qū)，中度鹽漬土逐漸變化為重度鹽漬土；南部沿海區(qū)域的鹽土面積略有減少，從鹽土水平變化到了重度鹽漬土水平，面積比例由48.15%變化到40.43%，同時(shí)鹽土仍集中在沿海分布。

整體來看，當(dāng)?shù)叵滤幌陆禃r(shí)，鹽漬化程度有所改善，主要表現(xiàn)為鹽土面積大幅減少，改變了最初以鹽土所主導(dǎo)的鹽漬化現(xiàn)象。地下水位的上升，加重了土壤鹽漬化程度，主要體現(xiàn)在中、重度鹽漬化水平上，使得中度鹽漬化面積下降，重度鹽漬化面積大幅升高。

3 結(jié)論

（1）地理探測(cè)器中因子探測(cè)結(jié)果表明，各因子對(duì)土壤鹽分的影響力大小為地下水礦化度＞植被覆蓋度＞地下水埋深＞距海遠(yuǎn)近＞粘粒含量＞地表高程＞降水量＞蒸發(fā)量；交互作用探測(cè)結(jié)果表明地下水因素、植被覆蓋度、距海遠(yuǎn)近之間的交互作用影響力較強(qiáng)，由此選定土壤鹽分的主要驅(qū)動(dòng)因子為地下水礦化度、地下水埋深、植被覆蓋度、距海遠(yuǎn)近。

（2）4種預(yù)測(cè)模型的精度評(píng)價(jià)表明，BP神經(jīng)網(wǎng)絡(luò)模型的建模集R2為0.8847，RMSE為1.1350，驗(yàn)證集R2為0.7999，RMSE為1.1204，因此選用BP神經(jīng)網(wǎng)絡(luò)模型作為土壤鹽分的最佳預(yù)測(cè)模型。土壤鹽分預(yù)測(cè)結(jié)果分析顯示，鹽漬化預(yù)測(cè)面積結(jié)果與實(shí)際面積結(jié)果較為一致：輕度鹽漬土預(yù)測(cè)結(jié)果為0.27%，實(shí)測(cè)面積比例為0.38%；中度鹽漬土預(yù)測(cè)結(jié)果為25.72%，實(shí)測(cè)面積比例為27.91%；重度鹽漬土預(yù)測(cè)結(jié)果為25.73%，實(shí)測(cè)面積比例為23.90%；鹽土預(yù)測(cè)結(jié)果為48.29%，實(shí)測(cè)面積比例為47.81%。

（3）地下水變化模擬情景結(jié)果表明，地下水位的適當(dāng)下降改善了土壤鹽漬化狀況，在對(duì)照組、地下水位下降0.5 m、地下水上升0.5 m 3種情景下，輕度鹽漬土面積比例分別為0.39%、0.61%、0.37%，中度鹽漬土面積比例分別為27.76%、22.30%、12.99%，重度鹽漬土面積比例分別為23.70%、38.98%、46.21%，鹽土面積比例分別為48.15%、38.11%、40.43%。

本研究對(duì)黃三角瀕海區(qū)進(jìn)行分析探究，篩選出土壤鹽分的主要驅(qū)動(dòng)因子，構(gòu)建了最佳預(yù)測(cè)模型，并根據(jù)設(shè)定的情景對(duì)土壤鹽分進(jìn)行模擬分析，為黃河三角洲土壤鹽漬化的調(diào)控及防治提供了依據(jù)。

4 討論

土壤鹽分含量的高低受到多種因素的共同作用，除自然因素外，人類活動(dòng)也會(huì)在一定程度上對(duì)鹽漬化產(chǎn)生影響，通過篩選土壤鹽分的主要驅(qū)動(dòng)因子發(fā)現(xiàn)，地下水因素對(duì)土壤鹽漬化的影響較大，影響因子的篩選結(jié)果與馬玉蕾[26]、周在明[27]、阮本清[28]等結(jié)果相似。氣候、地形、土壤質(zhì)地等這些自然因素不會(huì)輕易改變，但人為因素對(duì)土壤鹽分的影響綜合體現(xiàn)在地下水水位和水質(zhì)方面[29]，因此要加強(qiáng)地下水的有效調(diào)控與動(dòng)態(tài)檢測(cè)，對(duì)水資源進(jìn)行合理的利用，實(shí)現(xiàn)對(duì)土壤鹽漬化的高效治理。

在黃河三角洲土壤鹽分預(yù)測(cè)模型的探究中，機(jī)器學(xué)習(xí)模型的精度顯著高于統(tǒng)計(jì)分析模型[30-32]，其中BP神經(jīng)網(wǎng)絡(luò)模型表現(xiàn)出較高的模擬精度，這與前人的探究結(jié)果相一致[33-35]。通過此方法對(duì)研究區(qū)的土壤鹽漬化現(xiàn)狀進(jìn)行分析，得到該區(qū)域的土壤鹽漬化預(yù)測(cè)分布圖，該圖符合研究區(qū)土壤鹽漬化分布的實(shí)際狀況，因此BP 神經(jīng)網(wǎng)絡(luò)模型為快速準(zhǔn)確的對(duì)黃河三角洲瀕海區(qū)土壤鹽漬化預(yù)測(cè)分析提供可能。但是不同地區(qū)的土壤鹽漬化狀況不同，本研究所建立的最佳預(yù)測(cè)模型是否在黃河三角洲其他地區(qū)所適用，還需進(jìn)一步的驗(yàn)證。

對(duì)于地下水位變化的情景模擬結(jié)果可以看出，當(dāng)?shù)叵滤幌陆岛?，地下水補(bǔ)給的蒸發(fā)量減少，鹽分在土壤中的積累作用減弱，能夠在一定程度上減輕土壤鹽漬化程度，因此，情景化模擬可以對(duì)鹽漬化發(fā)展趨勢(shì)進(jìn)行有效的探究預(yù)測(cè)。但當(dāng)前對(duì)于鹽漬土的影響因子進(jìn)行情景化模擬的研究較少，鑒于鹽漬化的影響因子有多種，本研究只對(duì)地下水這一個(gè)因素進(jìn)行了模擬，在之后的研究中，對(duì)其他因子以及因子間的綜合作用也要進(jìn)行相應(yīng)的情景模擬分析。