结合耦合传输-反应模型解析DGT通量：土壤/沉积物中溶解配体与有限吸附容量对化学有效性的影响

推荐：土壤和沉积物中污染物的生物有效性评估是环境科学的关键问题。薄层梯度扩散技术 (DGT) 作为一种原位被动采样工具，在此类评估中应用广泛。然而，在非搅拌介质（如土壤/沉积物）中，DGT累积数据的解读缺乏精准的数学工具。为此，研究人员开发了开放源代码程序 SSOCO，并推导了分析解，以更准确地预测DGT在不同参数下的累积行为，并揭示其物理化学意义。该研究为利用DGT评估污染物（如重金属、有机污染物、营养盐等）在复杂环境介质中的动态形态和有效浓度提供了更强大的理论框架和实用工具。

想象一下，我们如何能精准评估土壤或河流沉积物中，有多少重金属、抗生素或营养盐是真正能被生物（比如植物、微生物）利用，从而可能被吸收、富集，甚至进入食物链的？这并非一个简单的“浓度是多少”的问题，而是一个关于“有多少是真正可用的”的动态难题。在环境科学与地球化学领域，薄层梯度扩散技术 (Diffusive Gradients in-Thin-films, DGT) 就像一位沉默的侦探，能通过原位被动采样，捕捉目标分析物在环境介质中的动态通量，从而反映其“生物有效性”。它已被成功应用于评估从重金属（如铜、镉、铅）、有机污染物（如抗生素、多环芳烃）、放射性核素到营养盐（如磷酸盐、硝酸盐）等多种物质的有效性，在监测沉积物修复效果、研究生物地球化学氧化还原循环或评估施肥效果等方面至关重要。

然而，当DGT这位侦探深入到静止的土壤或沉积物中执行任务时，解读其收集到的“证据”（即分析物在结合胶中的累积量）变得异常复杂。与在流动水体中可以快速达到稳定状态不同，在非搅拌介质中，由于缺乏对流补充，分析物在DGT装置周围的扩散区域会随时间不断扩展，导致其累积量与时间并非简单的线性关系，也无法用简单的解析表达式将累积量与扩散系数、孔隙度等关键物理化学参数直接关联。这使得从DGT测量结果中提取精确的、定量的环境信息变得困难重重。尽管现有的一维DIFS (DGT Induced Fluxes in Soils and Sediments) 模型等工具被广泛使用，但它们存在局限性，例如无法准确模拟短期部署，且通常未考虑溶解配体形成部分可动络合物、土壤吸附位点可能饱和等复杂但现实的情景。因此，开发一个能更全面、更真实地模拟土壤/沉积物中复杂相互作用，并能够处理有限吸附容量和部分可动络合物的新模型，成为提升DGT技术解释能力和应用价值的关键。

为解决上述挑战，由 Josep Galceran、Guanlei Li、Jordi Sans-Duñó、Carlos Rey-Castro、Jaume Puy、Yue Gao 和 Joan Cecília 组成的研究团队在《Geoderma》上发表了一项重要研究。他们开发了一套强大的数学工具和开源代码，旨在更准确地解释和预测当DGT设备部署在土壤或沉积物中时，目标分析物的累积行为。这项工作的核心是建立并验证了一个名为SSOCO (Sediments and SOils with COmplexation) 的一维数值模型，该模型明确包含了土壤和凝胶的孔隙度差异，并首次将土壤吸附位点描述为遵循朗缪尔等温线的、浓度有限的固定配体，从而能够模拟土壤吸附位点的潜在饱和效应。同时，模型还考虑了目标分析物与溶解配体、树脂结合位点之间的动力学结合/解离过程。此外，研究人员还为分析物单纯在两种介质中扩散这一简单情况，推导出了一个精确的解析解（可通过提供的电子表格计算），该解可进一步扩展到分析物被吸附和/或与溶液中其他物种（如形成可动络合物）相互作用的某些情况。通过比较，他们指出了广泛使用的DIFS模拟程序在短部署时间下精度有限的问题。这些工具共同构成一个强大的框架，使得从DGT测量中提取复杂土壤/沉积物系统中更真实、更有意义的信息成为可能，即使在非理想“完美吸附”或“配体过量”条件下亦然。

本研究主要采用了理论建模、数值模拟与解析推导相结合的技术方法。首先，研究人员基于质量传输和反应动力学原理，构建了一个一维数学模型（SSOCO模型），该模型耦合了扩散过程以及分析物与溶解配体、土壤吸附位点、DGT树脂之间的动力学反应。模型通过有限元方法在Fortran中实现为开源代码。其次，针对仅有分析物扩散的简化情况，通过求解偏微分方程组，推导出了累积量随时间变化的精确解析解表达式。最后，利用开发的数值程序和解析解，系统地研究了不同特征参数（如扩散系数、孔隙度、络合/吸附动力学常数、土壤吸附容量等）对DGT累积量和再供给参数R的影响，并与现有模型（DIFS）的结果进行了对比验证。

3.1. 无反应时M扩散的分析解

为建立基础，研究首先考虑了最简单的“扩散情形”，即目标金属M不与土壤颗粒或任何配体发生吸附或络合反应。在此情况下，M仅在材料扩散层（凝胶）和土壤孔隙水两个介质中扩散，分别具有扩散系数 D_M^g和 D_M^s。通过求解相应的扩散方程，结合完美的吸附边界条件（在树脂/凝胶界面浓度为零）和初始条件，研究人员推导出了累积量 n_M(t) 和再供给参数 R_diff(t) 关于时间的精确解析表达式。这个解析解避免了数值模拟中可能出现的数值误差，并且作为更复杂情况下的基准和扩展基础。

3.2. 快速络合和/或吸附极限情况下的扩展

研究表明，上述纯扩散情形的解析解可以扩展到两种重要的极限情况。当溶解络合物ML的形成/解离动力学非常快（即完全可动）时，以及当分析物M与土壤吸附位点S的吸附/解吸动力学非常快时，系统可以分别用等效的、具有“表观”扩散系数的纯扩散问题来近似描述。对于快速络合，有效扩散系数是游离M和络合物ML的扩散系数按其平衡浓度加权的平均值。对于快速吸附（假设土壤位点过量），有效扩散系数需除以一个因子 (1 + K_ads')，其中K_ads' 是条件吸附常数，这反映了固相吸附对总体传输的缓冲作用。这些扩展使得解析解能够应用于更广泛的实际场景。

3.3. 再供给参数R的近似表达式及其峰值

再供给参数 R = c_DGT/ c_T,M* 是DGT数据解读中的关键指标，其中c_T,M* 是孔隙水中总溶解分析物浓度。研究发现，R随时间并非单调变化，而是呈现先上升后下降的峰值行为。通过对解析解进行简化，研究人员推导出了R峰值出现的时间 t_peak和峰值大小 R_peak的近似表达式。例如，t_peak大致与 (δ^g)²/ D_M^g成正比，而 R_peak总是小于1。这些近似公式为了解和控制DGT测量中R的行为提供了直观的指导。

3.4. 长期部署行为与有效孔隙水浓度

对于足够长的部署时间，分析物的累积量 n_M(t) 表现出与时间的平方根成正比的渐近行为。这一特性使得可以从长期DGT累积数据中，估算出一个“有效”的土壤溶液浓度 c_E。研究表明，c_E不仅取决于总溶解浓度 c_T,M*，还受到土壤和凝胶的孔隙度比值 (φ_s/φ_g)、扩散系数比值以及络合/吸附动力学的影响。在快速平衡的极限下，c_E可以简化为一个与上述参数相关的表达式，这为利用DGT长期测量来反推介质中的有效浓度提供了理论依据。

3.5. SSOCO模型的应用：参数影响研究

利用新开发的SSOCO数值程序，研究团队系统探讨了多个关键参数对DGT响应的影响。例如，他们发现土壤孔隙度 φ_s的增大会降低R值，因为增大的孔隙空间稀释了扩散通量。土壤吸附强度 (通过K_D或K_ads' 体现) 的增加会显著降低R，并延缓R峰值的出现，因为更多的分析物被固相滞留。溶解配体的存在及其与M形成的络合物的可动性会显著影响R：如果络合物是完全可动且扩散较快，它可能提供额外的传输途径从而提高R；如果动力学较慢（部分可动），则可能限制供应。特别重要的是，SSOCO模型能够模拟土壤吸附位点有限（非过量）的情况。研究表明，当DGT的萃取导致局部土壤吸附位点趋于饱和时，R值会比位点无限过剩的假设下预测的要高，这突显了考虑有限吸附容量在污染严重或吸附位点有限的土壤中的重要性。

3.6. 与DIFS模型的比较

研究将SSOCO模型和新的解析解的结果与广泛使用的DIFS模型进行了比较。发现在短时间部署（例如小于24小时）下，DIFS预测的累积量和R值与精确解存在显著偏差，且DIFS预测的R峰值出现时间可能不准确，峰值高度也可能被高估。这主要是由于DIFS在数值离散化上的一些近似处理所致。因此，对于需要高精度解读短期DGT实验，或研究早期动力学过程的工作，使用SSOCO或本文的解析解更为可靠。

结论与意义

本研究通过发展SSOCO数值模型和一系列解析解，显著增强了对DGT在土壤和沉积物中测量结果的解释和预测能力。主要结论包括：

1. 1.
   提供了精确的数学工具：为分析物在两层介质中纯扩散的情形提供了严格的解析解，该解可作为基准并可扩展到快速络合/吸附的极限情况。
2. 2.
   揭示了关键参数的影响规律：通过参数研究，定量揭示了土壤孔隙度、吸附强度、络合可动性、吸附位点有限性等对DGT响应（累积量、R值）的影响，提供了更深入的物理化学见解。
3. 3.
   提出了新的概念和近似：推导了长期部署下累积量与时间的平方根关系，并定义了有效浓度c_E的概念；给出了R峰值时间和大小的近似公式，便于实验设计和数据快速解读。
4. 4.
   开发了更先进的模拟工具：SSOCO作为开源代码，其创新在于同时考虑了孔隙度差异、溶解配体的部分可动性、树脂结合动力学，并首次引入了基于朗缪尔等温线的有限土壤吸附容量，从而能更真实地模拟复杂环境介质中的情况。
5. 5.
   指出了现有工具的不足：通过与DIFS的比较，明确了其在短时模拟中的局限性，为DGT用户在选择模型时提供了重要参考。

这项研究的重要意义在于，它将DGT数据的解读从传统的、往往依赖于简化假设和数值近似的框架，推向了一个更严谨、更灵活、更能适应复杂真实环境的新阶段。所开发的开源工具和解析解不仅有助于环境科学家和地球化学家更准确地评估污染物和营养盐的生物有效性、动态形态和迁移通量，也为优化DGT实验设计（如部署时间选择）、发展新的DGT构型（如考虑有限吸附容量的影响）以及将DGT与其他模型耦合应用提供了坚实的理论基础和实用手段。最终，这将提升我们在土壤修复、环境风险评价、农业养分管理和全球生物地球化学循环研究等方面的能力和精度。