新疆土壤剖面盐渍化的时空变化及其对冰川径流变化的响应

《CATENA》：Spatiotemporal changes of soil profile salinization in Xinjiang and their response to glacial runoff variations

　　

研究人员首次在区域尺度上定量分析，将新疆冰川径流与深层土壤盐渍化联系起来。研究人员重建了新疆土壤剖面盐渍化（1960–2022）并分析了其时空演变。研究人员描绘了新疆山地-绿洲-荒漠系统（Mountain-Oasis-Desert System, MODS）中的土壤剖面盐渍化格局。研究人员揭示了剖面盐渍化对冰川径流和多源环境变量的差异化响应。关键发现：盐度热点、垂直格局以及对冰川径流的敏感性——深层（60–100 cm）受影响最大。

**论文解读：新疆土壤剖面盐渍化的时空变化及其对冰川径流变化的响应**

**研究背景**
土壤盐渍化是全球干旱区最普遍的土壤退化形式，严重影响作物生产与生态系统功能。在新疆，山地-绿洲-荒漠系统（MODS）中，盐分随水迁移，冰川融雪水是河流的主要补给源（占总径流量的45%以上）。然而，现有研究多聚焦于表层盐渍化监测，对深层（0–100 cm）剖面的长时间尺度、区域尺度变化及其与冰川径流的关系认识不足。尤其是冰川加速退缩背景下，盐渍化如何响应冰川融水变化尚不明确。因此，揭示土壤剖面盐渍化对冰川径流的宏观响应机制，对干旱区水安全、盐渍化防控及政策制定具有重要科学价值。

**研究内容与结论**
该研究发表在《CATENA》上，以新疆为研究区，利用中国第一、二次冰川编目（CGI-1、CGI-2）数据驱动开放全球冰川模型（OGGM）模拟1960–2022年冰川径流，并结合随机森林模型与多源环境变量（包括气候、地形、土壤属性、植被等）反演同时期0–10、10–20、20–40、40–60、60–100 cm五个深度的土壤电导率（ECe）。随后，采用地理探测器模型分析冰川径流及其他因素对盐渍化的贡献及交互作用。主要结论：新疆土壤盐渍化以表层聚集为主，年均ECe在0–10 cm达74.47 mS/cm，但绿洲区呈现“表层-底部”复合型聚集；过去60年整体呈增加趋势，但绿洲区近30年（1990–2022）有所下降；冰川径流总体增加，但多个流域已过峰值（如JSR、URB等）；冰川径流对盐渍化的影响在60–100 cm深度持续增强（年均q值0.63），且随深度增加，主导因素由温度、地形转为植被和土壤属性，最终变为土壤属性与冰川径流的联合控制；交互作用检测显示，冰川径流与其它变量的共同效应在深层（60–100 cm）最强（q值0.77–0.89）。

**主要技术方法**
1. 随机森林模型：在Google Earth Engine（GEE）平台集成多源环境变量，模拟1960–2022年五个深度（0–10、10–20、20–40、40–60、60–100 cm）的土壤饱和提取液电导率（ECe）。训练数据来自研究团队在新疆多个绿洲（如渭干-库车河三角洲、艾比湖流域、玛纳斯河流域、于田绿洲等）2005–2022年采集的979个（0–10 cm）至451个（60–100 cm）样本，以及世界土壤信息服务（ISRIC）提供的28,666个ECe样本（1980–2014，覆盖山区）。
2. 开放全球冰川模型（OGGM）：基于CGI-1和CGI-2冰川编目数据，模拟21个含冰川流域的冰川径流时空变化，并校准两期数据间的系统偏差。
3. 地理探测器：以流域为单元，采用自然断点法分类环境变量，计算因子检测q值和交互作用，量化冰川径流及多变量对土壤盐渍化的贡献及其动态变化。

**研究结果**
**4.1 精度验证**
通过十折交叉验证，表层（0–10 cm）ECe模拟精度最高（RMSE=4.29 mS/cm，R2=0.71），深层（60–100 cm）R2为0.52，所有深度均达到可接受范围。

**4.2 土壤盐度时空变化**
新疆剖面ECe年均值排序：0–10 cm（74.47 mS/cm）> 40–60 cm（27.93）> 10–20 cm（22.85）> 60–100 cm（22.62）> 20–40 cm（21.11）。过去60年各深度均呈波动上升趋势，增速分别为0.48、0.19、0.10、0.19、0.02 mS/cm·a?1。绿洲区呈现先升（1960–1990）后降（1990–2022）趋势，而荒漠区持续上升。空间上，三个盐碱热点区：艾比湖-天山北坡-玛纳斯湖周边、吐鲁番-哈密盆地、罗布泊周围。

**4.3 冰川径流时空变化**
1960–2022年新疆冰川径流总量介于18.32–31.92 Gt，平均23.86 Gt；多数流域在1960–1990年呈下降趋势，之后上升。叶尔羌河年均径流最大（4.48 Gt），其次为渭干河（3.66 Gt）、伊犁河（2.67 Gt）等。冰川物质平衡持续负值，平均减少4.46 mm w.e.·a?1。

**4.4 土壤剖面盐渍化对冰川径流变化的响应**
在非山区（绿洲+荒漠），大多数流域的剖面盐度随冰川径流增加而增加，但绿洲区盐度与径流呈反向关系（1960–1990径流减少时盐度增加，1990–2022径流增加时盐度下降）。因子检测显示，冰川径流对0–60 cm盐渍化的影响在1960–1990年增强，之后减弱；但对60–100 cm深度的影响持续增强（年均q值0.63，2013年达0.79）。交互作用检测表明，所有因子对均表现为双因子或非线性增强，且冰川径流与其它变量的交互作用在60–100 cm深度最强（q值0.77–0.89）。

**讨论与结论**
讨论部分指出，研究结果与全球及中亚地区盐渍化格局一致，绿洲区灌溉深层的盐分淋洗导致60–100 cm出现“底部聚集”现象。冰川径流模拟与已有研究（如Huss & Hock, 2018；Rounce et al., 2020）吻合，部分流域（如JSR、URB等）已过峰值，而伊犁河、阿克苏河等将在2034–2060年达到峰值。深层盐渍化受土壤属性与冰川径流联合控制，且交互效应随深度增强。局限性包括：缺乏实测径流验证、CGI-1定位误差、山区土壤样本缺失等。

**研究结论**
（1）新疆土壤盐渍化主要表现为表层聚集，年均ECe为74.47（0–10 cm）、27.93（40–60 cm）、22.85（10–20 cm）、22.62（60–100 cm）、21.11 mS/cm（20–40 cm）。绿洲区呈现“表层-底部”复合型聚集，荒漠区自上而下递减。过去60年各深度盐度总体上升，年增量分别为0.48、0.19、0.10、0.19、0.02 mS/cm·a?1，但绿洲区近30年（1990–2022）呈下降趋势。（2）过去60年新疆冰川径流总体增加，但多数流域在1960–1990年下降，若干流域（如JSR、URB、ETNR、TPB）已过峰值（2010–2014年）。叶尔羌河年均径流最大（4.48 Gt），其次为渭干河（3.66 Gt）、伊犁河（2.67 Gt）等。（3）冰川径流对盐渍化的影响在1960–1990年持续增强，之后减弱，但在60–100 cm深度持续增强（年均q值0.63）。随着深度增加，主导因素从温度、地形转变为植被、土壤属性，最终变为土壤属性与冰川径流联合控制。交互作用显示，所有因子对均增强，冰川径流与其它变量的交互在深层最强（q值0.77–0.89）。该研究揭示了新疆土壤剖面盐渍化的时空格局及其对冰川径流的响应机制，为水资源管理、盐渍化防控和生态保护提供了科学依据。