《Geochemistry, Geophysics, Geosystems》：Subseafloor Temperature Variability During Ridge-Flank Hydrothermal Circulation: Implications for a Moderate-Temperature Earth System and Other Ocean Worlds
编辑推荐：
海底热液循环驱动着固体地球与其海洋之间的生物地球化学交换，并可能在其他海洋世界中发挥类似作用。本研究评估了胡安·德富卡洋脊东翼约3.5 Ma海底下热液循环过程中的海底温度变化幅度。这一经过充分研究的、中等温度系统（上部洋壳中约65°C–70°C）代表了地球上的
海底热液循环驱动着固体地球与其海洋之间的生物地球化学交换，并可能在其他海洋世界中发挥类似作用。本研究评估了胡安·德富卡洋脊东翼约3.5 Ma海底下热液循环过程中的海底温度变化幅度。这一经过充分研究的、中等温度系统（上部洋壳中约65°C–70°C）代表了地球上的暖端成员，其缓慢流动跨越数十公里。补给区与排泄区之间的适度侧向流动伴随着小尺度对流，导致局部温度和流动条件的变化，即使在简单的层状模型中也是如此。流动和温度场随时间变化，海底变化的幅度围绕准稳态条件振荡。研究人员计算了与不同微生物栖息地相关的温度范围内地壳含水层的比例。在排泄流体约为70°C的模拟中，含水层中65%–90%的区域温度处于40°C–80°C。相同的模拟在适用于木卫二的低重力条件下运行时，含水层中50%–70%的温度为40°C–80°C，30%–50%的温度大于80°C。流体追踪表明，由于次级对流，流体可以以0.2°C–10°C/千年的速率加热和冷却。基于这些模拟的地球化学建模表明，对于这种端元情况，流体中某些营养物质和有机分子的浓度可能因这种温度变化而出现高达一个数量级的变化。这些发现突显了水岩反应变化的含义，进而揭示了可居住的海底生态位可能具有瞬时性质。
研究背景与科学问题
海底热液循环是驱动固体地球与海洋之间物质与能量交换的关键地质过程。在地球上，高温黑烟囱喷口虽视觉壮观且对局部元素收支影响显著，但从全球尺度看，洋脊侧翼的低至中等温度热液系统贡献了更大的热量和质量通量，并为微生物活动创造了有利条件。然而，以往对洋脊侧翼系统的数值模拟多聚焦于流体排放温度、流速以及与海底热流观测相匹配的整体特征，而对补给区和排泄区之间地下走廊内流体流动和热状态的空间与时间变异性关注不足。这种变异性直接影响水岩反应、矿物稳定性、化学动力学以及沿热液路径的微生物过程，因而亟需系统评估。
本研究以胡安·德富卡洋脊东翼约3.5 Ma洋壳上的热液虹吸系统为典型案例。该区域沉积速率极高，数百米厚的低渗透沉积层覆盖于火山基底之上，迫使流体通过火山露头进行补给和排泄。该系统是地球上表征最完善的洋脊侧翼热液系统之一，其流体排放温度约65°C–70°C，被视为中等温度系统的暖端成员。已有三维数值模拟成功再现了其观测到的排放速率和热流分布，但这些模拟结果尚未被用于深入剖析地下含水层内部的精细温度变化格局及其对生物地球化学过程的潜在影响。
研究内容与主要结论
研究人员通过对已有三维耦合流体-热量输运数值模拟结果的再分析，系统量化了洋脊侧翼热液虹吸系统中地下含水层温度的时空变异幅度，并探讨了这种变异性对微生物适居性和有机地球化学示踪物的意义。研究发现，即使补给区和排泄区相距约50 km，在整体侧向流动的背景下，含水层内部仍发育有小尺度对流细胞，导致局部温度在空间和时间上持续波动。在排放温度约72°C的模拟中，含水层中65%–90%的区域温度处于嗜热菌适宜的范围（40°C–80°C）。在适用于木卫二低重力的模拟中，该比例降至50%–70%，同时有30%–50%的区域温度超过80°C，进入超嗜热菌范围。流体粒子追踪显示，由于次级对流，流体可经历高达0.2°C–10°C/千年的升温或降温速率。结合地球化学建模，研究人员进一步揭示，±2°C–12°C的温度波动即可导致某些营养元素（如硫）浓度产生近一个数量级的变化，而钙、镁、铁等元素的浓度变化也可达数倍。此外，研究人员还探索了利用苯甲胺-苯甲醇-氨体系等有机分子作为地下温度示踪物的可行性，发现在较高温路径（如深部含水层案例）中，这类化合物可快速响应并锁定热液信号，即使在排放后迅速冷却至接近冰点的海水温度，其比值也能长期保存。
关键技术方法
研究人员采用有限元热质输运代码FEHM进行三维瞬态热液流动与温度场模拟。模型域设计基于先前成功复现胡安·德富卡洋脊东翼露头间热液流动特征的模拟方案，包含450 m沉积层、600 m渗透性地壳含水层及低渗透基底。模拟运行至准稳态后继续演化0.5–1.0 Myr，期间按特定时间步长输出流动和温度场数据。研究人员选取补给-排泄之间圆柱形区域评估含水层温度分布，并划分四个微生物相关温度区间。通过后处理粒子追踪技术，基于动态稳态速度场采用Runge-Kutta积分法计算流体粒子路径，获取沿程温度与流速演化。在地球化学方面，研究人员采用EQ3/6软件进行平衡态水岩反应模拟，并结合Arrhenius关系推导的温度依赖反应速率常数与热力学软件计算的平衡常数，建立可逆化学反应动力学模型，预测有机示踪物浓度演变。
研究结果
3.1 含水层与流体温度的物理预测
对于浅部含水层的地球模型，含水层中65%–90%的区域温度处于40°C–80°C，3%–9%低于40°C，其余3%–31%处于80°C–120°C。高渗透率案例因对流更强，更多热量被带出，故高于80°C的区域比例较低。深部含水层案例排放温度达114°C，温度分布更均匀，约2%超过120°C。特定位置的温度振荡幅度在高渗透率案例中约为±15°C，在低渗透率或深部案例中可达±20°C–30°C。温度变化速率大多为0.2°C–0.8°C/千年，唯低重力低渗透率案例仅0.05°C–0.07°C/千年。低热输入案例同样出现次级对流，温度横向变幅约±7°C，表明即使在地球最古老洋壳对应的低热流条件下，温度变异性依然存在。流体粒子的温度变化速率可远高于固定位置，常见短期速率达5°C–10°C/千年。流速在高渗透率案例中为0.5–6 m/yr，低渗透率案例相应降低。
3.2 地球化学预测的示例规模
以碳质球粒陨石成分为岩石、水岩质量比0.1为例，模拟Ea-1中±2°C–12°C的温度波动导致钙浓度变化约2倍，硫浓度变化近一个数量级，磷和碳几乎不变，镁和铁变化约4倍。对于有机示踪物，在深部含水层案例Ea-4的一条流体路径上，初始低温段反应缓慢，但当温度超过125°C后反应迅速，示踪物很快锁定当地条件。动力学模型显示，0.5 Myr于95°C或约1000年于145°C即可使该套示踪物响应不同反应条件。若流体排放后快速冷却至接近冰点，热液信号将被永久保存。
讨论与结论
讨论要点
讨论部分指出，温度是控制海水与洋壳化学反应速率和程度的主要因素之一，水岩比和流体/岩石成分也起重要作用。中等温度波动（如±15°C）即可显著影响营养元素浓度，进而影响微生物可利用的化学能和细胞丰度。氢气的生成随温度在10°C–120°C范围内可相差1–2个数量级。对于其他海洋世界，尽管缺乏板块构造和现代火山活动，但潮汐应力、过去岩浆活动以及可能的沉积物覆盖等因素，使得其海底可能存在类似地球洋脊侧翼的异质性结构，从而聚焦流体补给和排泄。在木卫二高渗透率模拟中，温度变异性与地球案例相似，但低渗透率下对流细胞无法形成。有机示踪物在木卫二案例中因平均温度偏低（75°C–80°C）而难以在短于2 Myr的时间内记录信号，但更高温系统或更低温度响应的示踪物（如乙醇-乙烯体系）可能适用。研究也承认其局限性，包括假设含水层均质各向同性、未考虑断裂和裂缝、未纳入反应输运过程等。未来工作应扩展至非均质介质、离散裂隙网络模型，并实现全瞬态粒子追踪以捕捉系统失稳事件。
研究结论翻译
本研究初步评估了由热驱动的水热系统中浅层含水层内可能出现的海底温度变异性，其中排放流体表现出中等温度（70°C–110°C）。结果为理解补给区和排泄区相隔数十公里时可能发生的时空变异性尺度提供了新的认识。本研究所探索的模型基于先前的数值模拟，这些模拟成功再现了胡安·德富卡洋脊东翼观测到的排放特征。研究表明，即使在相对简单的层状模型中，侧向流动也与小尺度对流共存，导致含水层温度在空间和时间上持续变化。对于所测试的大多数地球模型，含水层的大部分区域（65%–90%）处于40°C–80°C，这一温度范围支持嗜热微生物群落。在适用于木卫二低重力的模拟中，虽然整体流速减慢，但温度变异性仍然显著，且在足够高的渗透率下，相当一部分含水层温度超过80°C。通过将流体粒子温度演化与地球化学模型相链接，研究人员证明，即使适度的温度变化（±2°C–12°C）也可能导致某些营养元素和有机分子浓度发生数倍至近一个数量级的变化。此外，研究人员展示了苯甲胺-苯甲醇-氨体系等有机分子可作为地下温度的有效示踪物，特别是在较高温度的热液系统中，这些信号可能在流体排放后得以长期保存。这些发现强调了水岩反应的潜在变异性，以及因此产生的可居住海底生态位的瞬时性质，无论是对地球还是对其他海洋世界都具有重要意义。
