作为我国固有蓝色国土,南海及其岛屿是我国海上丝绸之路的重点建设区域,也是维护区域安全的战略要地。近年来,我国创造性地对珊瑚礁砂进行绞吸、耙吸等方式通过管道运输到礁盘,以实现对南海现有岛礁扩容,短时间内完成珊瑚岛基础设施建设。同时,通过珊瑚岛稳态工程、植被构建与深度绿化、生态淡水涵养等措施,推进建设生态活岛和宜居岛取得显著成效。这有效地巩固了海防,维护了我国海洋权益,满足建设海上安全通道的基本需求。
作为珊瑚岛淡水重要贮存空间,珊瑚砂中的水分是联系降雨、入渗、蒸散发、补给地下淡水资源等地表和地下介质间水文过程的关键环节,是建设生态岛的制约因子之一。珊瑚岛淡水来源以大气降水为主,降水入渗进入珊瑚砂体后形成稳定的地下淡水体
[1]。珊瑚礁砂作为一种主要成分为碳酸盐的海相砂土类介质,具有多孔隙、易破碎和胶结弱等特性
[2]。珊瑚砂颗粒在填筑过程中受自重和水力作用下发生分选。粗颗粒砂运移距离短,先沉积,而粉粒和黏粒随水流迁移沉积,再经机械振冲作用促进粗细分离、紧实,最后形成黏粒致密层。经过反复填压施工,地层剖面形成了厚度不等、质地差异显著的层状结构。剖面珊瑚砂物理特性,如质地、体积质量、颗粒机械组成和导水率,决定了珊瑚砂水文物理性质,影响着珊瑚砂持水性和导水性,进而改变珊瑚砂剖面中的水分动态、运移速率与方向,增强珊瑚砂水文过程异质性。可见,珊瑚岛地下淡水形成与演化,不仅取决于礁盘地质结构、大气降水和蒸散发等因子,还与珊瑚砂质地、剖面结构及其水力学性质等密切相关。因此,开展珊瑚礁砂剖面水力学性质、水文与砂输移特征的研究,有助于明确珊瑚岛致密层对地下淡水形成与演化等海岛水文过程作用机理。
国内外学者针对土壤剖面夹黏层性质特征及其效应展开系列研究。结果表明土壤在发育演化进程中因受到水文、地质、气象和人类活动等作用,剖面中分布有导水率差异显著的细质与粗质土层。例如,河流冲积平原区域因黏砂沉积差异形成的剖面
[3];农地因长期耕作形成的压实致密层或犁底层
[4];荒漠绿洲区的夹黏土层
[5]。有关土壤剖面黏土夹层的研究目前多集中在其基本物理性质、水分再分配、水文和盐分运移过程等方面
[6-8]。土壤剖面由于夹黏层孔隙分布变化,一方面明显减弱水分向下入渗,补给深层地下水,但可显著提高夹黏层以上的土壤水分体积分数,供给植被水分生长需求,如在荒漠绿洲区中出现夹黏土层区域的植被覆盖度高于相同水热状况下不含夹黏土层的区域
[9];另一方面,夹黏土层能减少因强烈蒸散发作用引起的地下水位上升作用,如在农田黏土层具有抑制深层水分提升能力,减少表土盐分积聚
[10-11]。可见,土壤剖面中夹黏层对水分、养分和盐分物质运移具有明显的双重影响。相关研究为珊瑚岛剖面致密层生态水文与工程效应研究提供了参考借鉴。
因建成过程特殊性,地下礁砂剖面分布有厚度不等的珊瑚砂致密层,这将使得珊瑚岛淡水循环过程在水平、垂直方向存在异质性,从而影响珊瑚岛礁地下淡水形成、演化速率和淡水储量等过程。目前,关于珊瑚岛的研究主要关注珊瑚砂工程力学特性、理化性质及海岛植被恢复等方面
[12-15],但针对珊瑚岛剖面致密层水力学性质及其水文与珊瑚砂输移过程的研究较为薄弱,仍需进一步开展探究。基于此,本研究以热带珊瑚砂剖面为研究对象,通过原位开挖珊瑚砂剖面调查与采样,测定珊瑚砂剖面粗颗粒层与致密层水文物理性质和水分特征曲线,分析剖面致密层导水和持水能力。同时利用模拟降雨试验,分析不同水文条件下致密层在降雨过程中径流、珊瑚砂输移特征,明确珊瑚砂致密层对坡面水文和输砂作用,探讨致密层对珊瑚岛地表水文过程、地下淡水形成与演化的影响,以期为深入了解其淡水循环过程与机理、植被恢复提供理论支撑。
1 材料与方法
1.1 研究区概况
研究区位于海南省三沙市南沙群岛中某一由钙质沉积物形成的珊瑚岛,剖面结构分层明显。土壤类型参照《第三次全国土壤普查暂行土壤分类系统(试行)》
[16]命名为石灰扰动人为新成土,归属于砂质碳酸盐型-石灰扰动人为新成土土属。该地属于典型热带海洋湿润气候,年平均气温为28.4 ℃,最高月平均气温30.3 ℃,最低月平均气温27 ℃。年平均降水量为1 808 mm,降水年内分布不均;最大月降水可达724 mm,最小月降水甚至为0 mm,降水量月际分布不均。岛上主要植被有木麻黄(
Casuarina equisetifolia Forst.)、榄仁(
Terminalia catapppa L.)、海滨木巴戟(
Morinda citrifolia L.)、草海桐(
Scaevola sericea Vahl)和狗牙根(
Cynodon dactylon (L.) Pers)等灌木和草本。
1.2 野外取样与室内分析
本研究于2023年9月在南沙新成珊瑚岛开展剖面调查。因此,本研究选取有致密层分布的代表性样点进行珊瑚砂剖面开挖,分别在0~10、10~20、20~40、40~60、60~80 cm深度处取样。每层取3个环刀(100 cm3)样品,分别用于测定珊瑚砂容重、珊瑚砂饱和导水率、水分特征曲线。同时,在相同深度处收集约50 g珊瑚砂样品,装入塑封袋用于测试珊瑚砂颗粒的机械组成。
珊瑚砂体积质量分数采用烘干称重法
[5]测定;饱和导水率采用定水头法
[17]测定;珊瑚砂颗粒机械组成采用马尔文激光粒度仪测定样品的黏粒、粉粒和砂粒百分比质量分数
[5],按照土壤颗粒分级系统对珊瑚砂进行质地分类。珊瑚砂水分特征曲线采用压力膜仪测定,每个样品测定9个吸力值(10、20、40、50、80、150、300、580、780 kPa)对应的珊瑚砂水分体积分数
[18]。采用van Genuchten模型对吸力值-水分体积分数关系进行特征曲线拟合,获得该模型相关参数
[19]。依照土壤水分特征曲线在33、1 500 kPa吸力值时土壤水分体积分数对应为土壤田间持水量和凋萎系数。根据各珊瑚砂水分特征曲线计算珊瑚砂有效水分体积分数
[17]。
1.3 室内模拟降雨试验
室内模拟降雨于2024年5月9日至6月1日在中国科学院华南植物园科研区(广东省广州市)进行。试验槽采用可调节坡度不锈钢材质,其规格为1.6 m(长)×0.5 m(宽)×0.5 m(深),表层和深20 cm处分别设有地表径流、壤中流收集口,试验钢槽见
图1。为确保水分能自由入渗,钢槽底部均匀分布有7排直径为5 mm的排水孔,每排3个排水孔。试验材料所用的珊瑚砂为珊瑚条经粉碎过筛后制成。根据致密层颗粒级配分析和前人
[20]研究结果,珊瑚岛礁致密层颗粒粒径<0.075 mm。同时,根据珊瑚岛剖面实地调查结果显示,剖面中含有多层厚度为10~50 cm粒径<0.075 mm的致密层。本研究为模拟致密层分布的坡面,致密层厚度设置为20 cm。本试验将粉碎珊瑚砂分别通过0.85、0.15、0.075 mm孔径。待珊瑚砂分筛完成后,先将少于0.075 mm的珊瑚砂装填至底层20 cm,浇水至饱和,让其自然沉降紧实,然后在上层20 cm装填0.85~0.15 mm粗质珊瑚砂。均质坡面则装填40 cm厚0.85~0.15 mm粗质珊瑚砂(
图1)。因新成珊瑚岛较平整,本研究试验槽坡度设定为3°。同时,考虑到热带珊瑚岛在枯水期水分能完全下渗;丰水期连续降雨导致下层珊瑚砂水分饱和及海水涨潮涌至珊瑚砂体形成壤中流等不同水文条件。因此,本试验对应设置3种不同水文条件,分别为自由入渗(free drainage)、饱和(saturation)和渗流(seepage)条件。珊瑚砂坡面饱和条件采用马氏瓶控制,马氏瓶进水管跟试验槽底部导水管连通,让试验槽20~40 cm间珊瑚砂充分饱和,在降雨过程中始终保持连通,直至降雨结束。渗流条件则是试验槽底部导水管连接1个自由水面高于20 cm的导水管,降雨试验进行时导水管连通水泵进行持续供水,试验槽深20 cm处形成1个有地下渗流界面
[21-22]。
试验模拟降雨装置为下喷式降雨机,喷头为Fulljet GW型号,喷头高度为3.0 m,有效降雨面积为3.0 m×3.0 m(
图2)。降雨均匀度为80%,降雨强度为60~140 mm/h。暴雨和大暴雨事件是影响珊瑚砂体水分与输移的重要因素之一,根据珊瑚岛气象站监测数据显示,最大日降雨量可达116 mm
[23]。因此,本试验设计降雨强度为120 mm/h,降雨历时为1 h,场降水总量为120 mm,每个处理重复1场降水。
降雨试验开始后,当收集口有连续径流时,记录起始时间及接取径流样。径流-珊瑚砂流失过程用1.0 L的塑料瓶收集,取样周期为3 min。当降雨结束后,进行样品称量、静置24 h,倒掉上清液后将珊瑚砂冲洗到铝盒中,放入105 ℃烘箱,待烘干后称重。由此可计算场降水引起的径流总量、珊瑚砂输移总量等特征
[24-25]。
1.4 统计分析
本研究采用IBM SPSS Statistics 20软件进行单因素方差分析和珊瑚砂水分特征曲线Van Genuchten模型拟合,采用SigmaPlot 14.0软件制图。
2 结果与分析
2.1 珊瑚砂剖面质地与水文物理性质差异
与自然形成的热带珊瑚岛相比,珊瑚岛剖面因构筑材料特殊性、成岛时间短,其剖面结构分布有间断沉积的黏粒层(
图1)。该黏粒致密层厚度约10~50 cm,黏粒质量分数约74.0%,粉粒和砂粒质量分数分别为17.2%和8.8%,质地类型为重黏土(
表1)。
剖面其他珊瑚砂层以粉粒和砂粒为主,其平均质量分数分别为47.5%和47.0%,黏粒仅占5.5%,质地类型为壤土、粉壤土和砂壤土(
表1)。此类珊瑚砂层体积质量为1.01~1.55 g/cm
3,平均体积质量为1.18 g/cm
3;而致密层体积质量约为1.68 g/cm
3,显著高于其他质地类型珊瑚砂层。
珊瑚砂质地和体积质量等性质差异必然影响其导水和持水性能。珊瑚砂饱和导水率为0.23~1.43 cm/min,平均值为0.69 cm/min;而珊瑚砂致密层饱和导水率为0.04~0.06 cm/min,仅为其他质地珊瑚砂饱和导水率的8%(
表1)。粉质或壤质珊瑚砂由于粉粒或砂粒质量分数大,颗粒间孔隙度大,水分吸力作用弱,其水分体积分数在低水吸力阶段随水吸力增加而陡降,尤其是砂质类型珊瑚砂表现最明显(
图2)。随着水吸力持续增加,珊瑚砂水分体积分数会缓慢降低,直至相对稳定。致密层黏粒质量分数高,具较大比表面积,毛管孔隙作用强,在低水吸力阶段水分体积分数减少相对缓慢,而在高水吸力阶段水分体积分数持续降低后趋于平缓。
因此,在珊瑚岛剖面的致密层含砂量少,颗粒细腻,体积质量大,导水能力弱,保水性能好;而其他质地类型珊瑚砂以粗颗粒为主,呈较弱持水和较强导水能力。
2.2 珊瑚砂坡面产流过程对剖面结构与水文条件的响应
均质与有致密层分布的2类剖面特征的坡面在模拟降雨条件下径流-珊瑚砂输移过程见
图3。结果表明,均质珊瑚砂(粒径为0.85~0.15 mm)与深层有致密层(粒径<0.075 mm)分布的珊瑚砂坡面对降雨产流响应过程差异明显。
由于珊瑚砂颗粒较大,全部入渗转化为珊瑚砂水分体积分数和向深层淋溶,2个坡面均无地表径流产生。在均质坡面,降水均向深层入渗,其淋溶速率总体呈现为快速增长至稳定产流(2.91~3.10 mm/min)变化特征。而在致密层分布的珊瑚砂坡面,降水在向深层运动过程中因致密层导水能力锐减引起水分滞留,导致密层界面形成横向径流与珊瑚砂流失。随降雨持续进行,珊瑚砂水分体积分数逐渐达到饱和,初始壤中径流速率由0.55~0.92 mm/min逐渐增大至稳定阶段(3.10~3.13 mm/min),其壤中径流速率与均质珊瑚砂坡面淋溶速率相近。整体上珊瑚砂坡面的降水分配主要转为壤中径流,平均占降雨总量的87.2%,而淋溶至深层的水量仅占12.8%。
近地表水文条件是影响珊瑚砂体水文过程的重要因子之一。由
图3可知,均质和有致密层分布的珊瑚砂坡面在自由入渗、饱和与渗流等3种水文条件下壤中流过程、径流总量发生明显分异。结果表明,均质坡面在自由入渗条件下入渗径流总量(131.88 mm)与饱和水文条件下壤中径流总量(135.07 mm)相近。均质坡面在饱和条件下其实质形成1个滞水界面,但入渗水流运动由垂直方向转变为横向壤中径流。
均质坡面在渗流条件下(速率为0.9 mm/min),壤中径流平均速率和总量显著增加,分别为饱和水文条件下的1.7、2.0倍(图54)。因致密层珊瑚砂具有较低导水能力,一方面可减少水分向下淋溶漏失,促进形成横向壤中径流;另一方面,也可降低由深层至浅层的渗流补水。
致密层分布的坡面在自由入渗、饱和和渗流水文条件下,其壤中径流量分别为116.99、132.73、143.31 mm。在饱和水文条件下,2类坡面壤中径流总量相近。而在自由入渗条件下,致密层降低约87%淋溶径流量;在渗流水文条件下,致密层降低约92%壤中径流量。因此,珊瑚砂剖面中致密层对坡面降水再分配具有调控作用,在垂直方向上降低降水向深层淋溶,削减渗流作用。
2.3 珊瑚砂流失过程对剖面结构与水文条件的响应
由
图3和
图4可知,坡面结构与近地表水文条件对珊瑚砂颗粒输移流失过程和总量均具有重要影响。在均质坡面,因全部降水向深层入渗而未监测到珊瑚砂颗粒随入渗水流在垂直方向上的迁移流失。由于在本试验钢槽底部淋溶液出口处垫有细孔塑料膜,阻滞细颗粒珊瑚砂垂向流失。对于有致密层分布的坡面,全部降水通过高入渗能力的粗颗粒珊瑚砂层并渗流至致密层,雨水沿松散层与致密层界面向下坡汇集形成壤中径流,引起珊瑚砂颗粒流失。由
图4可知,在初始阶段,壤中径流含沙量为0.66~5.76 g/L,均处于较低水平。随着降雨进行,珊瑚砂逐渐达到饱和,水分体积分数和壤中径流连通性增强,壤中径流速率迅速增加,其携带的珊瑚砂粒能力增强,含砂量增至168.53~191.41 g/L,造成珊瑚砂体内颗粒位移与流失。自由入渗、饱和和渗流3种水文条件对均质和有致密层分布的坡面壤中径流含砂率变化过程与输移总量影响见
图3和
图4。
相比于自由入渗条件,2类坡面在饱和与渗流水文条件下,降雨过程中含砂率波动变化更加剧烈。含砂率峰值均出现在降雨后期,同时出现多次峰值。一方面由于珊瑚砂体积含水量持续增加和水文连通性提升,横向壤中径流速率增强使侵蚀搬运能力增加;另一方面,由于珊瑚砂从剖面内部不断流失,造成砂体失稳崩塌,增加可搬运的砂体来源。整体上,2类坡面在饱和条件下珊瑚砂流失量相近,但在渗流条件下,均质坡面珊瑚砂流失量高于致密层坡面珊瑚砂流失量的8.2%。均质珊瑚砂坡面从自由入渗转换至饱和条件,导致形成壤中径流与珊瑚砂流失。当由饱和条件转换为渗流条件时,珊瑚砂流失量增幅为47.5%。致密层分布坡面在饱和条件时,珊瑚砂流失量为自由入渗条件下的1.5倍;由水分饱和条件演变为渗流条件时,珊瑚砂流失量增加约31.1%。
3 讨 论
3.1 剖面致密层与水文条件对珊瑚砂体水文-砂输移作用
珊瑚岛剖面致密层的形成,导致体积质量、导水率等水文物理性质差异显著,与土壤剖面中夹黏层的持水特性相一致
[5]。因此,珊瑚砂致密层显著改变剖面水力结构,质地极细的珊瑚砂颗粒形成的滞水层,为侧向壤中径流形成提供一个界面,导致截然不同的壤中径流和珊瑚砂输移过程。显著的坡面水文响应表现在降水入渗的横向和垂直再分配特征,珊瑚砂体内出现壤中径流。鞠远江等
[26]对不同级配珊瑚砂在水力梯度条件下渗透过程及其颗粒分析结果显示,珊瑚砂因胶结弱、结构松散,导致细颗粒易随渗流迁移分布,可改变剖面渗透特性;且随渗流强度、珊瑚砂体级配及其紧实程度相关。本研究结果表明短时间堆积且固结作用弱的珊瑚砂颗粒因横向壤中径流作用发生迁移,进而引起珊瑚砂体孔隙形成,增强珊瑚砂输移风险,与方乾等
[27]、李焱秋等
[28]在喀斯特地区水土流失过程相似。由于喀斯特地区土壤层厚度较薄,降水一部分能入渗到岩土界面或基岩裂缝引起土壤侵蚀;另一部分形成坡面地表径流和侵蚀,即喀斯特地区同时存在地表与地下水土流失。不同的是,本研究中坡面表层20 cm珊瑚砂颗粒较粗,入渗能力高于试验降雨强度;同时导水能力较弱的致密层起到类似喀斯特区基岩的作用,导致降水易形成壤中径流,出现地下珊瑚砂颗粒流失现象。因此,地层剖面中质地细腻且导水性弱的致密层,作为滞水层有助于促进横向壤中径流形成,进而影响珊瑚岛砂体水分横向、垂向再分配等淡水水文过程
[29-30]。
总之,由于近地表水文条件变化,不仅改变均质珊瑚砂体内渗流方向,增加横向水分再分配,而且显著加剧珊瑚砂流失总量。致密层分布的坡面在自由入渗、饱和水文情形下加剧珊瑚砂流失,但珊瑚砂致密层能削减渗流条件下砂体流失量。因此,致密层剖面横向壤中径流可引起珊瑚砂颗粒迁移,同时地下水文条件可加剧珊瑚砂输移强度。此特殊地表与地下双重水文过程与珊瑚砂流失输移机制增加了珊瑚岛结构生态效应的复杂性。基于野外海岛调查和模拟试验结果,本研究提出了海岛珊瑚砂颗粒流失与迁移模式。珊瑚砂体内水文-砂颗粒流失速率主要受珊瑚砂粒径、剖面结构、包气带水文条件、地下水位和降雨量等因子影响;而珊瑚砂颗粒流失总量则主要受珊瑚砂粒径、珊瑚砂体紧实固结程度、工程措施和植物措施等因子影响。然而,珊瑚砂体壤中径流及其引起的珊瑚砂颗粒输移如何影响海岛地下淡水水文过程、裂缝发育演变等方面仍是不清楚。同时,本研究仅考虑了一种级配组成的珊瑚砂与剖面致密层在强降雨强度下水文过程和珊瑚砂流失,仍有复杂的珊瑚砂剖面结构及其在多个影响因子驱动下珊瑚砂流失特征未被研究。因此,需要进一步开展试验与监测珊瑚岛剖面结构对珊瑚岛水文过程、地下淡水形成与蓄积等生态效应的影响,以便更好地明确珊瑚岛淡水循环过程机理。
3.2 剖面结构异质性的水文-珊瑚砂输移效应展望
珊瑚岛地下淡水透镜体的形成是解决淡水资源短缺、维持生态系统平衡、支撑人与自然生态系统的关键生源要素之一。天然形成的热带珊瑚岛,地层结构和地表生态系统相对稳定,形成稳定的地下淡水透镜体,其形成演变过程主要受降水、土壤结构、地形、水文地质和植被等多要素、多过程耦合作用
[31-33]。而新成珊瑚岛地层剖面结构与天然珊瑚岛存在显著差异,天然珊瑚岛的地下淡水透镜体形成和演变理论、模型难以解析新成珊瑚岛地下淡水循环机理。因此,天然降水能否形成有效、稳定的地下淡水资源是构建绿色、宜居和生态岛的根本保障。面对海岛淡水形成与演化机理以及生产、生活和生态需水等难题,需要在3个方面开展深入研究。
1)珊瑚砂颗粒间孔隙大,致密层分布导致珊瑚砂剖面结构空间异质性强。珊瑚砂体内水分运动可致颗粒流失、再沉积等物理搬运侵蚀和岩溶化学侵蚀
[34-36]。致密层能阻滞水分下渗、延长水分与珊瑚砂溶蚀作用;层间横向流促进珊瑚砂迁移流失,导致形成大孔隙
[20]。在长期岩溶作用下可能形成孔隙溶洞
[37]。因此,面对珊瑚岛降雨量大、气温高、良好的水热条件,需明确珊瑚岛地层结构对补给地下淡水和供给植被生长能力的影响,进一步开展有关热带珊瑚岛地下淡水水文过程及其溶蚀作用强度等研究。
2)珊瑚岛粗质颗粒砂持水能力有限,致密层会造成局部区域地下淡水补给减少,可提高其剖面上层珊瑚砂的水分,引起水文过程浅层化。而植被为适应珊瑚砂水分分布状况,其所需水分大多来源于表层,导致植被根系分布在近地表层,将不利于抗阻海上台风。因此,面对热带珊瑚岛淡水资源缺乏、土壤贫瘠,且高温、高风速和高盐雾等自然条件,需进一步筛选具有较强适应性和生态功能的植被种类,如深根系和主根发达、茎干粗壮、抗风能力强、耐盐雾等生理生态特征;同时,研发改良珊瑚砂性质的技术,优化植被恢复措施,从而促进热带珊瑚岛植被恢复。
3)植被恢复和生态重建是生态岛礁建设的核心内容之一,而地下淡水资源的形成和稳定是实现生态岛礁建设的基本前提。因此,面对珊瑚岛生产、生活用水、生态需水等多重压力,需进一步关注岛礁植被生长动态和植被水分利用特征,制定合理的植被恢复措施;亟待开展海岛地下淡水蓄积过程与机理研究,系统评估和预测海岛淡水资源生态承载力演变,研发促进地下淡水透镜体形成的技术方法和更高效、环保的节水措施,为海岛淡水资源利用与管理提供参考。
4 结 论
1)以黏粒为主要颗粒组成的致密层体积质量高,具有较低导水能力,可阻滞水分向深处入渗,促进珊瑚砂体内形成横向壤中径流。
2)因珊瑚岛剖面结构分异明显,形成时间短,胶结性弱,细颗粒珊瑚砂极易在径流作用下发生横向位移与流失。
3)珊瑚砂近地表的水文条件,如饱和或渗流条件,可加剧珊瑚砂迁移流失。珊瑚砂颗粒迁移将进一步重塑剖面珊瑚砂的孔隙度和导水率等水力性质。