庄河市地理位置：122°29′—123°31′ E，39°25′—40°12′ N（图1）。庄河市地处辽东半岛东侧，大连市东北部，是辽宁沿海经济带的重要节点城市。多年来，庄河市沿海地区城镇化进程加快，地下水开采过量，造成了海水入侵和土壤盐渍化的风险。根据2005年监测的氯离子浓度数据，2005年庄河市的海水入侵呈面状分布，西起明阳镇，东到栗子房乡以及石城岛和王家岛共7个海水入侵地段，合计海水入侵面积达124.5 km²。最严重区域位于黑岛—青堆子地段，入侵范围大，距离长，氯离子浓度最高^[22]；同时，由大连市农业科学研究院在沿海10大乡镇的研究，发现沿海地区各类盐渍土都有不同程度的分布^[23]。研究区域的自然环境概况见表1。

综合考虑庄河市气候、水文地质、地下水资料、海水入侵及土壤盐渍化历史状况，在《海水入侵与土壤盐碱化检测方法（试行）》技术规程指导下，确定地下水和土壤的采样点位置，于2018年4—5月（研究区土壤盐分含量最高的时期）完成采样。利用便携式GPS记录采样点的经纬度。

在海水入侵区、过渡带、未入侵区选取农业灌溉井和居民饮水井为地下水采样点，采样断面垂直于海岸线，每个断面3—5个样品，共采集50个地下水样品，在采集样品前，用皮尺加泡沫浮子物测绳测量水位，并将观测井上沿作为基准点，根据水准测量确定测点高程；随后先用取样水桶采集活水样清洗水样瓶（3次），通过0.45 μm后的滤膜后，在预先经过酸洗的聚乙烯采样瓶封装，贴上标签。运回实验室进行测定，测试指标包括地下水氯离子浓度和潜水矿化度。氯离子浓度采用硝酸银滴定法测定，按照《GB 11896—89 水质 氯化物的测定方法》进行测定；潜水矿化度采用重量法测定，参照《SL79—1994 矿化度的测定方法》进行测定。

土壤样品共采集92个，50个分布在地下水采样点附近，另外42个均匀分布在未采集地下水样品的地区。每个采样点呈S形随机采集3份样品，用土钻分5层（0—5 cm，5—10 cm，10—20 cm，20—40 cm，40—60 cm）采集土样，同层3个点所采的土壤等量均匀混合后作为一个混合代表样本，同时原状土和容重样分为3层（1—20 cm，20—40 cm，40—60 cm）采集土样；采集后的土壤样品在剔除植物根系及大块砂石后，装在聚乙烯样品袋中送回实验室进行测量。测量前将土壤在通风的室内风干，在土样达半干状态时，将大块土壤捏碎，并通过2 mm孔径的筛子过筛。随后配置1∶5的土水悬浮液，充分震荡后测量土壤pH、氯离子浓度（${\mathrm{C}\mathrm{l}}_{\mathrm{s}\mathrm{o}\mathrm{i}\mathrm{l}}^{-}$）、硫酸根离子（${\mathrm{S}\mathrm{O}}_{4}^{2-}$）、全盐含量。其中pH值采用电位法测定，按照《土壤检测土壤pH的测定》（NY/T 1121.2—2006）方法测定；${\mathrm{C}\mathrm{l}}_{\mathrm{s}\mathrm{o}\mathrm{i}\mathrm{l}}^{-}$采用硝酸银滴定法测定，按照《土壤检测 土壤氯离子含量的测定》（NY/T 1121.17—2006）方法测定；${\mathrm{S}\mathrm{O}}_{4}^{2-}$采用EDTA间接滴定法测量，按照《土壤检测 土壤硫酸根离子含量的测定》（NY/T 1121.18—2006）方法测定；全盐含量使用重量法进行测量，按照《土壤检测 土壤水溶性盐总量的测定》（NY/T 1121.16—2006）方法测定；此外，土壤干容重采用称重法进行测量。

应用SPSS、Excel软件，对采样数据进行描述性统计分析，包括平均值、标准差、变异系数等。通过相关指标剔除异常值，分析数据的集中或离散程度。

算术平均法是统计学中最基本、最常用的一种方法，能够直观、简明地评价采样数据的平均水平^[24]。

式中，${X}_{1},{X}_{2},{X}_{n}$为每一个采样点数据，n为采样点个数。

标准差是离均差平方的算术平均数的算术平方根，用$\sigma$表示，在概率统计中作为统计分布程度上的测量依据^[25]。

式中，${x}_{i}$代表第i个采样点处的测量值；${\bar x}$代表样品数据的平均值；n为采样点个数。

为消除测量尺度和量纲的影响，应用变异系数来反映数据离散程度，它是概率分布离散程度的一个归一化量度，其定义为标准差与平均值之比^[26]。

式中，$\mathrm{S}\mathrm{D}$代表样品数据的标准偏差，$\mathrm{M}\mathrm{N}$代表样品数据的算数平均值。

利用基于地统计原理的kriging插值法对地下水和土壤进行探索性空间全局预测分析，总体判断研究区海水入侵与土壤盐渍化状况^[27]。

式中，$Z\left({S}_{i}\right)$代表第i个采样点处的测量值；${W}_{i}$为空间权重，取决于采样点据预测位置的距离以及预测点与其周围的测量值之间的空间关系拟合模型；${S}_{0}$为预测点的位置；N为采样点个数。

在我国，一般用地下水矿化度（M）和氯离子浓度（${\mathrm{C}\mathrm{l}}_{\mathrm{w}\mathrm{a}\mathrm{t}\mathrm{e}\mathrm{r}}^{-}$）两项指标来衡量海水入侵的程度^[28]，根据饮用水卫生标准（GB 5749—2006），M的最大允许限值为1.0 g·L⁻¹，因此将M= 1.0 g·L⁻¹作为咸淡水分界线；但海水入侵陆地含水层造成地下水化学组分发生变化最为显著的是${\mathrm{C}\mathrm{l}}_{\mathrm{w}\mathrm{a}\mathrm{t}\mathrm{e}\mathrm{r}}^{-}$，且${\mathrm{C}\mathrm{l}}_{\mathrm{w}\mathrm{a}\mathrm{t}\mathrm{e}\mathrm{r}}^{-}$相对稳定，监测简单，因此，选择${\mathrm{C}\mathrm{l}}_{\mathrm{w}\mathrm{a}\mathrm{t}\mathrm{e}\mathrm{r}}^{-}$作为衡量海水入侵与否及入侵程度指标。参考《海水入侵监测技术规程（试行）》标准（表2），结合大连市规定的海水入侵标准，将${\mathrm{C}\mathrm{l}}_{\mathrm{w}\mathrm{a}\mathrm{t}\mathrm{e}\mathrm{r}}^{-}$=250 mg·L⁻¹作为海水入侵线，对区域海水入侵程度进行评价。

研究区地下水氯离子浓度和矿化度的范围分别在23.40 — 432.60 mg·L⁻¹和0.22 — 1.98 g·L⁻¹之间（表3），平均值分别为217.67 mg·L⁻¹和0.92 g·L⁻¹，满足饮用水相关标准（GB 5749—2006）。由海水入侵分级标准知（表2），地下水水质属于淡水水质^[29]。但从采样数据散点分布上看（图2），36%的样点的${\mathrm{C}\mathrm{l}}_{\mathrm{w}\mathrm{a}\mathrm{t}\mathrm{e}\mathrm{r}}^{-}$超过大连市规定的海水入侵标准250 mg·L⁻¹，最大超标倍数为73.04%，存在轻微海水入侵。同样的，样品中有40%的样点的M超过1 g·L⁻¹，最大超标倍数达到98%，属于微咸水水质。从地下水氯离子浓度和矿化度的变异系数来看，${\mathrm{C}\mathrm{l}}_{\mathrm{w}\mathrm{a}\mathrm{t}\mathrm{e}\mathrm{r}}^{-}$和M变异程度类似，都属于中等变异性，表明两者的离散程度都相对较大，研究区不同地区不同位置的浓度差异显著。

样点${\mathrm{C}\mathrm{l}}_{\mathrm{w}\mathrm{a}\mathrm{t}\mathrm{e}\mathrm{r}}^{-}$和离海距离呈现指数型递减规律，两者的相关系数（R²）达到0.86，具有较高的一致性。随着离海距离的增大，样点${\mathrm{C}\mathrm{l}}_{\mathrm{w}\mathrm{a}\mathrm{t}\mathrm{e}\mathrm{r}}^{-}$逐渐减小的速率越来越小，在离海距离为0—6 km时，${\mathrm{C}\mathrm{l}}_{\mathrm{w}\mathrm{a}\mathrm{t}\mathrm{e}\mathrm{r}}^{-}$最高降幅达240 mg·L⁻¹，最低下降到约250 mg·L⁻¹，而离海距离大于30 km后，${\mathrm{C}\mathrm{l}}_{\mathrm{w}\mathrm{a}\mathrm{t}\mathrm{e}\mathrm{r}}^{-}$基本保持稳定。采样点中，约1/3的采样点超标，属于海水入侵区（${\mathrm{C}\mathrm{l}}_{\mathrm{w}\mathrm{a}\mathrm{t}\mathrm{e}\mathrm{r}}^{-}$=250 mg·L⁻¹），入侵距离大约均在4 km以内。但有些特殊超标样点位于离海距离6 km的位置，入侵距离较远（图3a），这可能是部分沿海地带存在基岩裂隙，导致的海水呈线状入侵引起的。样点M随离海岸线距离的减小呈现逐渐增大的趋势，相关系数（R²）仅为0.25，一致性关系不显著（图3b）。但值得注意的是，在微咸水区，部分样点的离海距离达到10 km以上，最远的甚至达到20 km左右，这些样点的矿化度超标并不是因为海水入侵导致，可能是由于天然本地的矿化度较高导致，且这些点的数量较少，属于微咸水水质的样点几乎均分布在5 km以内，且离海距离为5 km时，样点矿化度位于咸淡水分界线上。海水入侵区和微咸水水质样点均分布在离海5 km以内，可以认为研究区矿化度超标也与海水入侵有关。将地下水氯离子浓度与矿化度进行线性相关性分析（图4）得到相关系数（R²）为0.87，可知地下水咸化与氯离子浓度密切相关，即与海水入侵密切相关。综上，沿海地带地下水氯离子浓度和矿化度升高是由海水入侵引起的，超过饮用标准和灌溉标准。

研究区地下水水位与潜水埋深范围分别为0.80 — 210.90 m和0.90 — 23.10 m，变化范围较大，潜水埋深的变异系数为88%，属于中等变异，离散程度比氯离子浓度和矿化度高，但小于地下水水位的离散水平（表4），表明地下水水位的空间分异性大，随地理位置和地形的不同，水位差异明显。

随着地下水水位的降低，氯离子浓度和矿化度逐渐增大，且地下水水位与${\mathrm{C}\mathrm{l}}_{\mathrm{w}\mathrm{a}\mathrm{t}\mathrm{e}\mathrm{r}}^{-}$的相关性较高（R²=0.70），而与M的相关性一般（R²=0.53）。样点海水入侵区地下水位大多数分布在13 m之内（图5），说明海水入侵集中在地下水水位较低处发生，特别是接近海平面地下水开采量大的地区，这些地区地下水水位相对较低，易导致海水入侵陆地含水层。

地下水水位高低与埋深并无直接关系，两者之间是由地表高程决定的。通过分析得到地表高程越低，潜水埋深就越浅，且两者呈现密切的正相关关系（R²=0.75）。从图6得出，埋深较浅的地方主要在沿海地带。综上，沿海地带地表高程普遍较低，且地下水埋深较浅，较易发生海水入侵，且部分地区满足土壤盐渍化的发生的条件。

研究区土壤颗粒整体呈砂粒（41%）＞粉粒（33%）＞黏粒（26%），从变异系数的角度看，在10% — 100%之间，具有中等强度的变异性。根据土壤质地的USDA标准，研究区土壤质地多为黏壤土（clay loam）及壤土（loam），其次为砂质粘壤土（sandy clay loam），部分地区存在砂质壤土（sand loam）和粉质壤土（silt loam）。研究区土壤容重的均值为1.40 g·cm⁻³（表5），紧实度相对较大，整体上看，土壤容重变异程度较低，属于弱变异水平。除东南沿海以及少部分滨海地区，其余大部分地区土壤的砂粒含量较高，最高可达70%左右，分布在东部行政边界处。粉粒、黏粒含量与砂粒含量的分布几乎相反，在砂粒含量较低的地方，粉粒和粘粒含量都相对较高。

研究区土壤pH值在6.20 — 8.30之间，酸性、中性、碱性土均有分布，土壤pH均值为6.95，属于中性水平；但其标准差与变异系数都较小，分别为0.61和9%（表6），说明土壤的酸碱性较为集中。研究区以中性土为主，占到总样本数的51%，碱性土次之，占比为29%，而酸性土最少，仅占20%（图7），因此，可以说明研究区土壤多非碱化土，少数为碱性和酸性土。

目前学者多以易测的饱和土壤浸出液电导率（EC）为指标^[30－33]，还有学者利用钠吸附比（SAR）、交换钠离子百分比（ESP）等指标^[34]来计算土壤溶液的离子之间的占比，进而分析盐渍化的水平。研究区土壤全盐含量最小值仅为0.1%，而最大值可达到0.8%，土壤结构极差；盐分含量的平均水平为0.4%，超过了部分植物正常生长的限值（表7）。从${\mathrm{C}\mathrm{l}}^{-}$/${\mathrm{S}\mathrm{O}}_{4}^{2-}$化学计量比上看，其范围在0.30 — 3.65之间，硫酸盐土、氯化物-硫酸盐土以及硫酸盐-氯化物型土壤均有分布。从变异强度上看，两项指标均属于中等变异性，而${\mathrm{C}\mathrm{l}}^{-}$/${\mathrm{S}\mathrm{O}}_{4}^{2-}$变异性相对更高。但从采样数据散点分布上看，研究区非盐渍化（38%）＞中度盐渍化（23%）＞轻度盐渍化（21%）＞重度盐渍化较少（18%）（图8），土壤盐渍化类型为硫酸盐型（46%）＞硫酸盐-氯化物型（34%）＞硫酸盐型（21%）。对比土壤样品盐渍化类型与程度，发现盐渍化土壤主要分布在氯化物-硫酸盐土以及硫酸盐-氯化物型土壤上，且在硫酸盐-氯化物型土壤上，盐渍化程度更为严重，多为重度盐渍土和中度盐渍土。土壤全盐含量与离海距离呈现对数型递减规律，两者的相关系数（R²）达到0.58。随着离海距离的增大，土壤全盐含量趋向于稳定。离海距离在10 km之内，土壤全盐含量下降趋势显著，随后则趋于稳定（图8）。综上，土壤中${\mathrm{C}\mathrm{l}}_{\mathrm{w}\mathrm{a}\mathrm{t}\mathrm{e}\mathrm{r}}^{-}$相对较高的土壤，其发生盐渍化的可能性较大，而土壤中氯离子浓度可能来由含盐地下水的毛细上升作用^[35]，沿海地带地下水氯离子浓度较高主要由海水入侵导致。因此研究区的滨海盐渍土也极可能由于海水入侵导致。

从庄河区土壤盐渍化空间分布看，盐渍化土壤多分布在沿海地带（图9），面积达1150 km²，约占陆地总面积的28%，内陆地区几乎全为非盐渍化土壤，特别是东北部的山区，这可能与该处地势高，雨量充分，蒸发强度小，没有造成盐渍化的外界条件有关。海水入侵集中在3块独立的区域，并没有连成片（图9），最大入侵距离达到5.21 km，入侵最大面积达到17.2 km²，总面积达29.81 km²，占陆地总面积的0.73%，3处区域均出现了不同程度的土壤盐渍化问题。

从空间分布来看，研究区沿海地带的地下水埋深很浅，几乎都在3 m以下（图10a），地下水浓度、矿化度含量普遍较高，存在不同程度的微咸水水质。参考《海水入侵监测技术规程（试行）》的标准并依据大连市规定的判断海水入侵的标准，以${\mathrm{C}\mathrm{l}}_{\mathrm{w}\mathrm{a}\mathrm{t}\mathrm{e}\mathrm{r}}^{-}$=250 mg·L⁻¹作为海水入侵的临界线，可以看出沿海地带均有不同程度的海水入侵，东南沿海地带入侵程度相对于其他地方更为严重。海水入侵线与咸淡水分界线几乎是重叠的，离海距离均在5 km左右，再次证实了沿海地区矿化度超标与海水入侵的密切相关性。对照土壤全盐含量空间分布（图10b），沿海土壤盐分远远高于内陆，内陆地区主是非盐渍化土，以中性氯化物—硫酸盐型土壤为主，特别是研究区东北部的山区，这可能与该处地势高，雨量充分，蒸发强度小等外界条件相关。而东南沿海地带的全盐含量高于其他沿海地带，以硫酸盐-氯化物型土壤为主，进一步表明了沿海地带土壤中氯离子浓度相对较高。从浓度变化来看，研究区的地下水盐离子浓度呈沿海向内陆逐渐递减的规律，土壤盐分也由沿海向内陆呈层状递减，沿海地区地下水盐分浓度明显高于内陆地区，且全盐含量高的区域地下水盐分浓度均在500 mg·L⁻¹以上。可见，研究区地下水浓度和土壤全盐含量在空间分布上呈现一致性的规律，而地下水浓度过高可能是土壤盐渍化的风险源之一，这与其他沿海地带研究呈现的规律是一致的^[36]。

土壤盐渍化问题是由自然、地理环境和人类活动综合作用的结果，受气候、地形、土体结构和地下水埋深、地下水开采等多种因素的影响。通过对庄河市土壤盐渍化影响因子相关性进行分析，可知土壤全盐含量与海拔高程、地貌类型、土壤质地、潜水埋深都呈负相关性，仅与潜水矿化度呈现正相关性。其中，土壤全盐量与潜水矿化度具有较大的正相关关系，相关系数为0.568，反映了盐分的积累和地下水的含盐量及地下水盐上行关系密切（表8）。由前面分析可知，海水入侵导致了庄河市沿海地带地下水氯离子浓度和矿化度超标，进一步表明了庄河市沿海地带的土壤盐渍化的形成主要根源在于海水入侵，且海水入侵贡献率可能达到50%以上。

在此基础上，利用SPSS对影响因子进行聚类分析（图11），将影响因子分为3类：全盐含量、潜水矿化度、潜水埋深为一类，反映土壤盐分主要受到潜水矿化度和潜水埋深的影响，体现影响水盐运移的主要过程；地貌类型和高程为一类，反映土壤的水文地质过程；土壤质地为一类，反映土壤本身性质对盐分运移的影响。聚类分析结果同相关性分析结果一致，全盐含量与潜水矿化度具有显著的相关性，表明海水入侵对土壤盐分有着直接且明显的影响，是土壤盐分积聚的主要驱动因素。同时沿海部分地区潜水埋深相对较浅也是导致盐渍化的主要驱因。

地下水矿化度、埋深与土壤全盐含量较高的相关性，进一步间接证实了庄河市海水入侵与土壤盐渍化的密切关联性。首先，沿海地带受到一定程度的海水入侵，导致了地下水的咸化，同时沿海地带地势较平，海水的倒灌导致了陆域的排水不畅，又进一步导致了地下水埋深浅化，加之地下水资源的开采，加剧了海水入侵的发生。因此，呈现出地下水埋深越浅，盐份含量越高的现象。同时，庄河市的气候属于典型的蒸发量大于降雨量型，旱季气温相对较高时，含盐地下水会通过毛管作用上升至土壤表层，水分蒸发后盐分便积聚起来，随着时间的推移，盐分积累就造成了土壤盐渍化，可见气候因素亦成为盐渍化问题形成的助力。因此，庄河市沿海地带海水入侵，带来了地下水盐分的变化，在气候因素的助推下，导致水生态和土壤生态环境失衡，从而导致了土壤盐渍化。

（1） 庄河市沿海地带地下水氯离子与矿化度的超标由海水入侵引起，区域内氯离子浓度和矿化度分别在23.40 — 432.60 mg·L⁻¹和0.22 — 1.98 g·L⁻¹之间，空间分异性大；且离海距离越近，地下水水位越低，海水入侵越严重；沿海地带地下水埋深较内陆浅，海水入侵的风险增加，导致部分地带达到了土壤盐渍化的发生条件。

（3）从空间插值分析结果看，庄河市沿海地带的地下水浓度与土壤盐渍化程度在空间分布上呈现较好一致性。氯离子浓度与离海距离呈现出指数型递减规律，相关系数为0.86，且递减的速率随着离海距离的增大不断减小，最后趋于稳定。以250 mg·L⁻¹为海水入侵临界值，最大入侵距离为6 km，土壤全盐含量与离海距离呈现对数型层状递减规律，相关系数为0.58，递减速率在离海距离约10 km时趋于平缓，离海距离约20 km后土壤全盐含量趋于稳定。由于地质原因存在空间差异，微咸水水质、重度盐渍化土均分布在沿海地带，以东南沿海地带最为明显。

（4）庄河市海水入侵是造成沿海地带土壤盐渍化的主要控制因素，贡献率达到50%以上；沿海地带较浅的地下水埋深和地下水位为海水入侵的发生提供有利条件，海水入侵通过影响地下水氯离子浓度和矿化度，使地下水咸化，在旱季气温相对较高时，含盐地下水可以通过毛管作用上升至土壤表层，水分蒸发后盐分积聚，造成了土壤盐渍化。

（5）土壤盐渍化阻碍了庄河市的水资源开发利用，使其沿海地带的土壤生态失衡，从而对经济发展产生影响。因而引起庄河市土壤盐渍化的直接诱因是海水入侵。因此，治理土壤盐渍化，要结合庄河市的实际情况同步采取防治海水入侵的措施。应控制地下水开采时间和间隔、开采井布局及密度，调控地下水水位和地下水埋深，加强水资源的管理配置，全面建设节水型社会，发展节水农业，提高节水技术，建设降低万元产值耗水量的相关工程。同时，可结合海绵城市建设，充分利用城市雨水系统选择恰当的“渗、蓄、排”措施改善沿海地带的灌—排形式，避免沿海地带挖沙、盐场建设等活动高强度发展，降低海水入侵的几率，防止咸水下渗，重视沿海湿地建设，为海水入侵竖起屏障。