海岛在我国广泛分布，尤其热带岛礁密集成群，具有重要的战略地位。我国热带海岛基本为珊瑚岛礁，表层为形成年代晚、松散未胶结或弱胶结的珊瑚和生物碎屑砂层，其下为珊瑚礁灰岩和珊瑚贝壳灰岩^[1]。这种土壤质地具有粗质多孔、易碎的结构特征，岛上丰富的降水很难产生地表径流^[1]，加之砂土缺乏天然有机物，从而形成了结构性缺水、土地贫瘠、生态脆弱的海岛水土资源环境系统。目前，我国大部分海岛的淡水供应主要依赖远距离运水、海水淡化及通过硬化设施收集并处理的雨水^[2]，对地下水的开发利用不足。拥有稳定的天然淡水资源对海岛的开发建设至关重要，而作为许多海岛地区重要水源之一的地下水^[3]，其对人类的有用性基本上取决于其化学特性。积极开展热带海岛入渗通过其砂质土层的水质特性研究，对于促进海岛资源的开发利用，保障海岛居民用水安全及生态系统维护具有重要的科学意义。

海岛地下水埋深一般较浅^[4]，其水化学成分受诸多因素控制，包括降水组成、地层岩性、水文地质过程以及人类活动等^[5]。表层珊瑚砂土作为联系热带海岛地区水文循环和生物化学循环的纽带，深受海洋环境和人类活动的双重影响，高渗透性能决定了珊瑚砂土层类似于“天然砂漏”，导致其所富含的易溶海洋盐类物质以及人为污染物质可随雨水入渗不断淋溶流失，并迅速到达地下水面，从而改变地下水的化学组成^{[4, 6]}。由于海岛地区普遍缺乏足够的水资源安全管理和健全的污水排放设施，水环境问题呈现出污染源多样化的特征，如工农业活动、人畜粪便和家用化学品流失等造成的地下水污染现象时常发生，使地下水质受到了威胁^[7-9]。有研究表明，我国某些热带海岛地下水因存在色度高、盐度高、有异味等问题而无法使用^[10]。由于海岛珊瑚砂土层受特殊成土母质影响，其入渗水的化学特征是及早发现、治理和防范海岛地下水污染问题的基础，在开发利用海岛淡水资源时需引起重视。

目前，有关入渗通过海岛珊瑚砂土层的水化学特性研究较少，在水文地质条件复杂的海岛地区开展影响因素与过程研究更是鲜见。本研究以我国热带海岛典型珊瑚砂土层为研究对象，通过人工降雨模拟试验，分析渗流通过珊瑚砂土水样的水质变化特征，并选用水化学相图、Visual MINTEQ水文地球化学模拟^[11]、碱基交换指数^[12]计算和Gibbs图解^[13]等方法，探讨通过珊瑚砂土层水样的水化学特征变化及主要控制因素与关键过程，以期为海岛地下水资源的合理开发与利用提供科学保障和参考依据。

1 材料与方法 1.1 砂土特性与机械组成

本研究选择我国典型热带自然珊瑚灰沙岛进行砂土采集。该海岛处于热带海洋性季风气候区，多年平均气温为26.5 ℃，平均降水量为1 500 mm左右，雨季为每年5—10月^[14]。其表层珊瑚砂剖面示意图如图 1所示。采集近海滩地、人工扰动较大的平地、林丛地(近年人工扰动小)3种样地表层砂土(编号分别为D1、D2、D3)，将其带回实验室自然风干、分析并用于开展模拟试验。

表 1为珊瑚砂土的基本理化性质。由此可以看出，珊瑚砂土为高盐碱化土，特别是Ca、Mg含量最高，这主要与其成土母质多为珊瑚、贝壳砂，主要化学成分是CaCO₃^[15]有关。由于珊瑚砂土都还未充分受生物质积累和分解、盐渍和脱盐、磷的富集和淋溶等成土过程^{[14, 16]}，因此其还不能称为严格意义上的土壤。

采集的砂土样中>2 mm的颗粒含量相对较高，为保持其原始状态，仅剔除极个别巨砾。参照SL237—1999《土工试验规程》^[17]，按照以下粒组^[18]：>2、2 ~ 0.5、0.5 ~ 0.25、0.25 ~ 0.125、0.125 ~ 0.063 mm和＜0.063 mm，对样品分次取样筛分，精度0.01 g，每组样品筛分前后样品总量误差必须小于1%。通过颗粒分析，发现3种样地砂土均以粒径>0.25 mm的粗粒为主，粒径＜0.25 mm的细粒含量较少。根据GBJ145—90《土的分类标准》^[19]和GB50007—2011《建筑地基基础设计规范》^[20]，将D1、D2和D3分别归类为粗砂、砾砂和中砂。

1.2 试验设计与样品分析

根据热带海岛表层为珊瑚砂土层特点，本试验装置(图 2)采用高40 cm，内部直径5 cm的有机玻璃土柱，设计珊瑚砂柱长为30 cm，采用分层装填法，每装一层(约5 cm)捣实一次，并计算容重，以确保砂柱装填均匀。D1、D2和D3砂柱装填容重分别约为1.40、1.53和1.35 g/cm³，介于测定的珊瑚砂堆积密度1.35 ~ 1.55 g/cm³。砂柱上下各铺一层玻纤棉，防止水流冲散砂样。利用校正好的蠕动泵(BL100-pump)，开展连续人工降雨(雨强为60 mm/h)模拟试验，并分别于5、10、20、30、40、50、60、80、100、120、160、200、240、300、420、600、840、1 320、1 800、2 280、2880 min采集水样，以分析水质的时间变化特征。试验设置3组重复，每批试验按照相同步骤做两个空白样，以减少样品误差。

由于该区域降水主要靠对流雨和台风，人类活动对降水的污染可忽略不计，雨水水质纯净，因此，试验用水选用去离子水，其主要理化指标见表 2。

根据GB/T 5750—2006《生活饮用水标准检验方法》^[21]收集和储存水样。采用台式水质分析仪(innoLab 20P和innoLab 20C)测定水样的pH、电导率(EC)、盐度(SAL)、总溶解固体(TDS)和水温(T)；水的碱度(HCO₃^-和CO₃^2-)采用碱度计(德国Aquanmerck公司)滴定测定，精度0.1 mmol/L；阴离子(F^-、Cl^-、NO₂^-、SO₄^2-、Br^-、NO₃^-和H₂PO₄^-)采用离子色谱法(ics-2100)进行测定，检出限为0.001 mg/L；阳离子(Na⁺、K⁺、Mg²⁺、Ca²⁺)及其他微量元素(Al、As、B、Ba、Cd、Co、Cr、Cu、Fe、Li、Mn、Mo、Ni、P、Pb、Se、Si、S、Sr、V和Zn)采用电感耦合等离子体光谱仪(ICP-OES: Optima 5300DV，美国PerkinElmer公司)完成检测，检出限为0.001 mg/L，仪器1 h内相对标准偏差≤0.5%。检测时先对每批试验的两个空白样进行测定，结果取均值。测定样品时，每20个样重复1次，以验证结果的稳定性。最后将测定结果扣除空白作为样品值。为了从整体上反映人工降雨条件下入渗通过不同砂土层的水质状况，将每一砂柱在不同时间采集的所有水样的实测数据求平均值。

1.3 研究方法

本研究选用水化学相图、Visual MINTEQ水文地球化学模拟、碱基交换指数计算以及Gibbs图解等方法揭示入渗通过砂土表层水的水化学特征及主要控制过程。

1) Visual MINTEQ水文地球化学模型。该模型基于化学平衡常数和较全面的热力学反应数据库，定量计算水溶液中矿物质的形态、溶解度以及固相和溶解相的平衡^[22]，并可简单快速地预测所有输入离子的固相饱和指数(saturation index, SI)，从而指示地下水中哪种溶质发生了溶解和沉淀反应。为评价入渗通过珊瑚砂层水样的潜在化学反应，本研究采用地球化学模型Visual MINTEQ 3.00计算了水中矿物相的饱和指数。以25℃被测水样中离子的摩尔浓度作为输入成分，输入分析水样各自对应的pH，运行并获得输出文件。与固体相比，若SI>0，溶液过饱和；若SI=0，溶液处于平衡状态；若SI＜0，则溶液不饱和。

2) 碱基交换指数。碱基交换指数(indices of base exchange，IBE)^[23]主要用于理解发生离子交换的化学反应，计算方法如下：

$ IBE = \frac{{C{l^ - } - (N{a^ + } + {K^ + })}}{{S{O_4}^{2 - } + HC{O_3}^ - + C{O_3}^{2 - }}} $

(1)

式中：Cl^-指所有卤族元素(F、Cl、Br、I、At、Ts)之和；Na⁺和K⁺代表所有碱金属元素(Li、Na、K、Rb、Cs、Fr)之和，以质量浓度计。当水中IBE＜0时，表示地下水中的Ca²⁺或Mg²⁺与含水层介质中的Na⁺和K⁺发生交换；如果IBE>0，则表示地下水与其围岩之间存在反向离子交换。

3) Gibbs比。通过Gibbs图解可以直观反映水样中TDS与占优势的阴阳离子比值的关系，进而描述水化学组成以及地层岩性与水化学之间的关系^[24]。Gibbs比的计算方法如下：

$ Gibbs比（阳离子）=\frac{N{\rm{a}}^{+}+{K}^{+}}{N{\rm{a}}^{+}+{K}^{+}+C{a}^{2+}} $

(2)

$ \rm{Gibbs}比（阴离子）=\frac{{\rm{Cl}}^{\rm{-}}}{{\rm{Cl}}^{\rm{-}}+{\rm{HCO}}_{3}^{\rm{-}}} $

(3)

当(Na⁺+K⁺)/(Na⁺+K⁺+Ca²⁺)或Cl^-/(Cl^-+HCO₃^-)比值大于0.5，且TDS含量较低时，水化学组分主要受大气降水影响；当(Na⁺+K⁺)/(Na⁺+K⁺+Ca²⁺)或Cl^-/(Cl^-+HCO₃^-)比值小于0.5，且TDS含量中等时，水化学组分主要受水岩相互作用影响；当(Na⁺+K⁺)/(Na⁺+K⁺+Ca²⁺)或Cl^-/(Cl^-+HCO₃^-)比值大于0.5，且TDS含量较大时，水化学组分主要受蒸发结晶作用影响^[25]。

2 结果与讨论 2.1 入渗通过珊瑚砂土层水样pH和EC动态变化

本研究结果表明，入渗通过不同样地砂土层水样的pH随降雨历时呈先快速上升、后缓慢变化至平稳的趋势(图 3)。入渗通过D1、D2和D3的水样pH分别介于7.92 ~ 8.81、7.86 ~ 9.30和7.82 ~ 9.35，与试验用水(pH 6.78)和大气降水(pH 6.0)pH均值^[26]相比碱性增加，甚至部分水样超过了世卫组织(WHO)^[27]规定的饮用水pH安全限值(6.5 ~ 8.5)。以上结果说明珊瑚砂土为高盐碱化土。

同时研究发现，降雨初期，入渗通过D3(林丛地)的水样pH最小。这与近年来该样地受人为干扰小且林下植被修复后土壤微生物的改善有关，如频繁的土壤微生物活动有助于改善土壤盐碱化状况，并可降低其pH^{[14, 28]}。

通过重复试验，不同时间入渗通过D1、D2和D3的水样电导率变异系数分别介于8.55% ~ 44.68%、8.46% ~ 23.53% 和2.41% ~ 28.67%，根据变异性分级标准^[29]，均属弱或中等变异性，说明试验重现性较好。人工降雨条件下，入渗通过不同样地砂土层的水样电导率呈先剧烈后缓慢、并在短时间内达到平稳状态的下降趋势(图 4)。降雨初期，入渗通过D1、D2、D3的水样电导率分别为48 496.7、89 900、77 600 μS/cm，并各自在300、200、240 min时降到了1 682.3、1 428.7、1 952.3 μS/cm，分别在初始水平基础上降低了96.53%、98.41% 和99.87%；累计降雨48 h后，入渗通过D1、D2、D3水样的电导率分别为131.4、111.6、98.4 μS/cm。由筛分结果可知，组成砂土的颗粒中黏粒含量极少，加之珊瑚砂土有机质含量低，导致它对各种盐分离子的吸附性能较差。因此，累积降雨条件下，珊瑚砂土层中的盐分不断淋溶流失，淋溶水质均可由盐水变为淡水。然而，影响珊瑚砂土层含盐量大小的因素有很多，比如气候、植被、水文、人类活动等，是各种环境因子影响的综合反映。我国热带海岛地区常年蒸发旺盛，尤其在持续干季，埋藏较浅(一般小于2 m)^[4]的海岛地下水会不断向地表运移，严重时会引发含水层疏干，导致咸水反渗盐化珊瑚砂，进而对岛上人类和动植物的生存环境构成威胁。

2.2 入渗通过珊瑚砂土层水样微量元素浓度变化

微量元素是水质综合评价的重要内容之一。虽然微量元素在水体中的含量甚微，但若在水质检测和评估时忽略它们的存在，就有可能给人类健康带来潜在风险^[30]。本研究从检测出的20种微量元素中选取在WHO饮用水标准^[27]中有规定限值的14种微量元素(Al、As、B、Ba、Cd、Cr、Cu、Fe、Mn、Mo、Ni、Pb、Se和Zn)以及P元素浓度进行了对比分析。由表 3可知，本研究入渗通过不同珊瑚砂土层的水样中除Ba、Cu、Cd(D3)外，其余微量元素浓度最大值均超出了饮用水标准值，但只有D3水样中的Al元素浓度最小值始终未降到饮用水标准值；从微量元素浓度均值来看，入渗通过3种类型样地水样中的Ba、Cd、Cu、Mn、Mo、Ni、Se、Zn、B(D3)元素浓度都达到了饮用水标准，相比之下，微量元素Al、As、B(D1和D2)、Cr、Fe、Pb较难溶于水，应该引起重视。此外，入渗通过D1、D2和D3的水样中P元素浓度最大值分别为2.062、0.877、1.778 mg/L，明显高于大气降水值(0.02 mg/L)^[31]。龚子同等^[31]研究发现，海岛土壤浅层地下水中P元素含量为0.2 ~ 4.9 mg/L，本研究结果在其范围内，说明热带海岛入渗通过珊瑚砂土层的水中P元素含量增加。这与热带海岛频繁的鸟类活动密切相关，由于海鸟粪中富含P等元素，海岛珊瑚砂经过脱盐阶段，在生物积累、淋溶、磷素富集、积盐脱盐过程等作用下，形成了以磷质石灰土为主的珊瑚砂土，因此鸟粪中富含的P元素是砂土中富含磷矿的主要来源^{[31, 15]}。

2.3 入渗通过珊瑚砂土层水样主要离子化学特征

我国热带海岛远离大陆，受热带海洋性气候和海水等因素影响，珊瑚砂土层盐碱化程度高。由表 4可知，入渗通过D1、D2、D3的水样TDS变化范围分别为24 246.67 ~ 65.70、44 933.30 ~ 55.80、38 766.67 ~ 49.20 mg/L，平均值分别为5 175.98、6 004.45、4 903.61 mg/L。当TDS＜600 mg/L时，水质最好，当TDS>1 000 mg/L时，饮用水水质较差^[27]。因此，入渗通过珊瑚砂土层的水样TDS，不适合饮用。由各离子浓度均值可推断出，整体上不同珊瑚砂土中的主要阳离子和阴离子含量大小分别为：Na⁺>Mg²⁺>Ca²⁺>K⁺和Cl^->SO₄^2->HCO₃^->NO₃^->CO₃^2-> Br^->F^-。根据WHO标准^[27]，入渗通过D1、D2和D3水中的Na⁺、Ca²⁺、Mg²⁺、F^-、Cl^-、SO₄^2-、Br^-和NO₃^-最小浓度值都在饮用水标准范围之内，但最大值远超过饮用水标准。虽然K⁺、CO₃^2-和HCO₃^-没有规定允许限值，但浓度太高也会影响水质，严重时会危及人体健康。因此，入渗通过珊瑚砂土层的水样TDS主要受Na⁺、Ca²⁺、Mg²⁺、Cl^-、SO₄^2-、Br^-、HCO₃^-和F^-等控制，而这些离子恰好是海水的主要组分，这说明海水混合作用深刻影响着海岛表层珊瑚砂的化学组成。同时，珊瑚砂土中碳酸盐类矿物(如文石、方解石)的溶解也是Ca²⁺、Mg²⁺和HCO₃^-等的主要自然来源；而淋溶于水中的NO₃^-和SO₄^2-等营养盐类物质主要是人为原因所致。

由图 5可知，在连续48 h的降雨过程中，入渗通过不同珊瑚砂土层的水样中主要离子浓度都呈现出先剧烈下降、后缓慢变化并逐渐趋于平稳状态的趋势。本研究发现，入渗通过D1的水样离子浓度出现波动上升现象，是因为珊瑚砂虽然属于粗颗粒砂土，但其具有颗粒大小不一，且内、外孔隙丰富的微观结构特征，这就导致砂柱内必然会有优先流产生，从而造成水中溶解物的浓度会出现峰值。根据WHO标准^[27]，入渗通过D1(图 5A、D)、D2(图 5B、E)和D3(图 5C、图 5F)的水样中阴阳离子浓度都达到饮用水标准所需降雨历时分别为2 880、2 880、600 min，其中阳离子浓度达到饮用水标准的先后顺序为Ca²⁺、Mg²⁺、Na⁺，而阴离子没有统一规律。不同的是，入渗通过D1(图 5A)、D2(图 5B)和D3(图 5C)的水样中占优势的阳离子分别在累积降雨1 800、2 280、1 800 min后由Na⁺变为Ca²⁺，且Ca²⁺浓度基本保持不变，而其余阳离子浓度继续保持略微下降趋势。入渗通过D1(图 5D)、D2(图 5E)和D3(图 5F)的水样中占优势的阴离子分别早在累计降雨100、120、80 min及以后由Cl^-变为HCO₃^-或CO₃^2-。其中，入渗通过D1的水中始终未检测到CO₃^2-，入渗通过D2和D3的水样分别在累计降雨840和240 min后检测到了CO₃^2-。此外，D1(图 5D)、D2(图 5E)和D3(图 5F)分别在累计降雨淋洗2 280、2 280、240 min之后，在水样中也检测不到Br^-。

2.4 入渗通过珊瑚砂土层水样水化学类型

利用水化学相图可以判断水样主要离子组成变化，体现不同水样的化学组成特征，从而辨别其控制端元。如图 6所示，降雨初期，入渗通过不同珊瑚砂土层的水样化学类型主要为Cl-Na水。随着降雨的不断冲刷，珊瑚砂土层所含盐分不断被淋失，水化学类型也在不断发生变化。其中，入渗通过D1的水样水化学类型变化情况为：Cl-Na水→HCO₃·Cl-Na水→HCO₃·Cl-Na·Mg水→HCO₃-Na·Mg水→HCO₃-Na·Ca·Mg水；入渗通过D2和D3水样水化学类型变化情况为：Cl-Na水→HCO₃·Cl-Na水→HCO₃-Na水→HCO₃-Na·Mg水→HCO₃-Na·Ca·Mg水。总体来说，人工降雨条件下，入渗通过不同珊瑚砂土层水样水化学类型都由Cl-Na水转变为HCO₃-Na·Ca·Mg水。入渗通过砂土层水样水化学类型发生这一变化，主要与珊瑚砂土的成岩环境、成土母质以及所含矿物溶解度密切相关。珊瑚砂土为海相沉积物，盐碱化程度高，主要矿物组分为富含钙镁的碳酸盐岩，特别是海水作用导致土层中的Na⁺和Cl^-含量最高。在降雨淋溶作用下，珊瑚砂土层中的各种盐类依次溶出，最先是氯化物溶出，形成高矿化度的Cl-Na水；当渗入水继续冲淡岩层时，氯化物不断被消耗，则此时可形成矿化度较低的HCO₃·Cl-Na水或HCO₃-Na水；当氯化物几乎全部被渗入水从珊瑚砂土层中驱替后，继续渗入则会发生碳酸盐的溶解和淋滤作用，从而形成低矿化度的HCO₃-Na·Mg水或HCO₃-Na·Ca·Mg水。

2.5 入渗通过珊瑚砂土层水样水化学控制过程

表 5列出了基于Visual MINTEQ模型和碱基交换指数方程的计算结果。就钙、镁矿物而言，入渗通过不同砂土层的水样中都含有大量文石[CaCO₃]、方解石[CaCO₃]、白云石(有序)[CaMg(CO₃)₂]、碳酸钙镁石[CaMg₃(CO₃)₄]、菱镁矿[MgCO₃]和球霰石[CaCO₃，多晶]，其中，入渗通过D1的水样中矿物SI平均值分别为1.86、2.01、4.73、5.84、1.58和1.44；入渗通过D2的水样中矿物SI平均值分别为1.86、2.00、4.68、5.69、1.53和1.44；入渗通过D3的水样中矿物SI平均值分别为1.82、1.96、4.66、5.72、1.55和1.40，这充分说明组成砂土的主要矿物组分是文石、白云石和方解石等碳酸盐矿物，主要化学组分为CaCO₃^[32]。根据SI最小值(表 5)可知，以上6种矿物在入渗通过不同砂土层的所有水样中都为饱和状态，甚至严重过饱和(根据SI最大值)，这主要与珊瑚砂土本身的理化性质有关。一方面，珊瑚砂土主要由弱(无)胶结的粗粒松散沉积物组成，渗透性强，吸附性能较差；另一方面，其化学成分深受不同生物及矿物组分制约^[33]，例如文石质的造礁岩生物(如珊瑚、仙掌藻、多孔螅、多毛类等)和方解石质生物(如牡蛎、浮游有孔虫等)的Ca含量很高，一般>37%；镁方解石质生物(如珊瑚藻、软珊瑚古针、棘皮动物、苔藓虫、底栖有孔虫等)含Mg高，Ca含量稍低，约27% ~ 36%。因此，在充分淋溶作用下，这些矿物会在短时间内被大量析出，从而导致其SI值出现严重过饱和现象。

入渗通过D1(图 7A)、D2(图 7B)和D3(图 7C)的水样中，以上几种矿物的SI值随降雨历时整体上呈减小趋势。不同的是，入渗通过D1和D2的水样中矿物SI值在降雨初期出现再次增大的现象，这主要与降水在多孔珊瑚砂层中的渗流机制和饱和状态密切相关。筛分试验结果表明，D1、D2和D3分别为粗砂、砾砂和中砂。相比较而言，D2颗粒最粗，渗入水通透性最好，而D1和D3则对水流具有较大的阻滞作用。珊瑚砂土颗粒大小不一会导致砂柱内部孔隙分布不均，流经较大孔隙的水流(也称优先流)会最先溶出部分以上矿物，待整个砂柱被入渗水完全浸润饱和时，同类矿物会再次被析出，从而导致这些矿物的SI值在整体减小的过程中出现再次增大现象。

除碳酸盐以外，珊瑚砂土层还含有少量硅酸盐。基于硅酸盐和碳酸盐风化循环反应，文石、方解石、白云石、碳酸钙镁石、菱镁矿和球霰石可以在水中沉淀出来，其基本化学反应如下：①二氧化碳在水中溶解形成碳酸：CO₂ + H₂O → H₂CO₃；②碳酸在风化过程中分解硅酸盐岩石并向地下水释放阳离子：H₂CO₃ + H₂O + 硅酸盐矿物→ HCO₃^- + 阳离子(K⁺、Na⁺、Ca²⁺、Mg²⁺ 等)+ 黏土；③溶解物之间的反应产生新的固体，在风化环境中比原始基岩矿物更稳定。例如，在珊瑚砂地层地下水中，Ca²⁺和HCO₃^–反应生成方解石：Ca²⁺ + 2HCO₃^- → CaCO₃ + CO₂ + H₂O。

基于Visual MINTEQ模型计算结果，硅酸盐矿物的溶解和文石、方解石、白云石等碳酸盐矿物的沉淀可能是影响入渗通过珊瑚砂土层水样化学组成的重要机制之一。而其他钙、镁矿物都不饱和，水纤菱镁矿[Mg₂CO₃(OH)₂·3H₂O]、水镁石[Mg(OH)₂]、泻利盐[MgSO₄·7H₂O]、萤石[CaF₂]、石膏[CaSO₄]、水菱镁矿[Mg₅(CO₃)₄(OH)₂·4H₂O]、石灰[CaO]、三水菱镁矿[Mg₅(HCO₃)(OH)·2H₂O]、方镁石[MgO]和氢氧钙石[Ca(OH)₂]等矿物的SI值均小于0(表 5)。另外，岩盐[NaCl]、芒硝[Na₂SO₄·10H₂O]、泡碱[Na₂CO₃·10H₂O]、无水芒硝[Na₂SO₄]和水碱[Na₂CO₃·H₂O]等钠矿物以及岩盐[KCl]也不饱和。因此，碳酸盐矿物(不包括碳酸钠和碳酸钾)的沉淀以及氯化物和硫酸盐矿物的溶解度深刻影响着入渗通过珊瑚砂土层水样的化学组成。

地下水与其围岩环境之间的离子交换对于了解地下水中的溶质来源也很重要^[33]。本研究结果表明，珊瑚砂中富含Na、K、Ca、Mg等组分。入渗通过D1、D2和D3水样的IBE值变化范围分别为2.1 ~ 0.01、2.6 ~ 0.01、3.12 ~ -0.05，平均值分别为0.73、0.7、0.93(表 5)。入渗通过不同砂土层水样的IBE值基本为正(图 7D)，表明在降雨入渗过程中，已经淋溶于水体中的Na⁺和K⁺与珊瑚砂土层中还未彻底溶解的Ca²⁺和Mg²⁺的交换很普遍。同时发现，入渗通过不同珊瑚砂土层的水样，其IBE值都随降雨历时的延长在逐渐减小，并趋于0，表明阳离子交换作用在减弱，这是因为珊瑚砂土中溶解能力相对较强的海洋盐类物质(Na⁺和K⁺)随降雨冲刷最先被淋失，降雨后期其淋溶水化学类型几乎完全受制于矿物沉淀和溶解过程，水化学组分继而被碳酸盐矿物中溶解出来的Ca²⁺和Mg²⁺等替代。由此可见，若是在珊瑚砂土含水层的一个完整的地下水循环系统中，阳离子交换作用也是影响其地下水化学类型的主要控制机制之一。

由图 8可知，D1、D2和D3分别在降雨300、240、120 min之前，其淋溶水样主要受蒸发结晶作用影响，(Na⁺+K⁺)/(Na⁺+K⁺+Ca²⁺)和Cl^-/(Cl^-+HCO₃^-)比值都接近于1，TDS>1 000 mg/L，其中最大的TDS值接近海水特性；之后水岩相互作用增强，具有中等偏高的(Na⁺+K⁺)/(Na⁺+K⁺+Ca²⁺)比值和中等偏低的Cl^-/(Cl^-+HCO₃^-)比值，且TDS＜1 000 mg/L。由此可以确定，热带海岛入渗通过珊瑚砂土层的水化学成分深受蒸发结晶、水岩相互作用和海水混合作用的影响，这主要与海岛所处的海洋环境和干湿季分明的气候密切相关。研究结果可以进一步揭示，当海岛地区正值雨季时，海岛珊瑚砂土层中的各种海洋易溶盐类物质随降雨入渗不断被淋失，水–岩相互作用也越来越明显，这与人工降雨中后期珊瑚砂土层淋溶水中的矿物沉淀占优势结果相一致。

3 结论

1) 本研究入渗通过D1、D2和D3水样的TDS最大值分别为24 246.67、44 933.3、38 766.67 mg/L，表明热带海岛表面珊瑚砂土层盐碱化程度普遍较高。但在人工降雨条件下，珊瑚砂土层中的盐分不断淋溶流失，其淋溶水化学类型均可由高矿化度的Cl-Na水转变为低矿化度的HCO₃-Na·Ca·Mg水。

2) 相比其他微量元素，入渗通过珊瑚砂土层水样中超标微量元素Al、As、B(D1和D2)、Cr、Fe、Pb的浓度难以下降，需引起重视。

3) 控制入渗通过珊瑚砂土层水样水质的水化学机制主要包括蒸发结晶作用和岩石沉淀作用，以及降雨入渗过程中Na⁺与K⁺、Ca²⁺与Mg²⁺离子交换作用的影响。其中，本研究入渗通过不同珊瑚砂土层的水中都含有过饱和的文石[CaCO₃]、方解石[CaCO₃]、白云石[CaMg(CO₃)₂]、碳酸钙镁石[CaMg₃(CO₃)₄]、菱镁矿[MgCO₃]和球霰石[CaCO₃，多晶]等碳酸钙镁矿物，这充分说明热带海岛珊瑚砂土本身富含这些矿物。