基于近岸系列分层竖井的海底地下水排泄及其营养盐输入研究
2021-04-01吴自军王富康崔振昂AliyaZidan朱虹霓王一清
吴自军,王富康,崔振昂,高 龙,Aliya Zidan,朱虹霓,王一清
(1.同济大学海洋地质国家重点实验室,上海200092;2.广州海洋地质调查局,广东广州,510075;3.广西海洋地质调查研究院,广西北海,536000)
海岸带是海陆地之间相互连接、相互作用的过渡地带。由于海岸带滨海含水层与海洋自由相通,地下水在由陆向海输送过程中与近岸海水发生交换混合作用,最终经海底释放进入近岸水体,该过程被认为是海底地下水排泄(submarine groundwater discharge,SGD)。由于受陆地和海洋驱动力的双重作用,SGD由2个分量组成:一个是来自陆地的淡水排泄,为大气降水入渗补给地下形成的地下水,包括潜水和承压水,其向海底渗流的驱动力主要为水力梯度;另一个为再循环咸水排泄,即海水在潮汐、波浪、咸淡水密度差等作用下通过海陆界面进入近海含水层后又回归到海洋[1]。
海底地下水排泄作为海岸带陆海相互作用的典型过程,不仅是全球水循环的重要组成部分,更是多种溶解化学物质从陆地向海洋输送的一个重要的通道,在全球海洋物质循环及其收支平衡中发挥着不容忽视的作用[2]。诸多研究显示,海底地下水排泄为近海提供了重要的营养盐、微量元素、碳等物质,沿岸海域海水富营养化、赤潮暴发、重金属和有机物污染,甚至近岸海洋生态系统的失衡,皆与海底地下水排泄密切相关[3]。
广西北海市近岸海域水体面临着富营养化的威胁。据统计,在1995年至2005年间,北部湾海域共报告了赤潮灾害18次。诸多研究认为赤潮发生的原因主要是由陆源径流、工农业废水、城市生活污水等携带高含量营养盐排入近海所致[4],而陆源地下水排泄的影响与贡献却长期被忽视。为调查陆地地下水与临近海水的相互作用与影响,广州海洋地质调查局在广西北海市南部冯家江西岸钻取了系列分层地下水竖井,其分别位于海水入侵区、咸淡水分界线和未入侵区(淡水区),这为研究陆源地下水营养盐分布特征、地球化学过程及其对近海的环境影响提供了绝好条件。本研究通过野外现场监测和采样,分析系列分层竖井不同层位地下水的氯离子、硫酸根离子、营养盐、溶解无机碳(DIC)及其碳同位素(δ13C-DIC)组成等,查明近岸地下水水化学组成及咸化特征,重点探讨不同竖井及不同层位的地下水营养盐含量的空间变化及其制约因素,并讨论地下水携带的营养盐输入对近海生态系统的潜在影响。
1 研究区地质及环境概况
北海市地处广西壮族自治区南部,紧邻北部湾北岸,在构造上属南康盆地的西南段,其地形呈向西突出的半岛状,即北、西、南三面临海(图1a,图1b)。北海市除西面为基岩残丘外,总体上主要为河海混合松散层堆积的滨海平原。受地质构造的制约,南康盆地构成了一个相对独立的水文地质单元,其主要含水层为地表覆盖数米至十多米厚的中更新统北海组潜水含水层,下伏第四系下更新统湛江组承压含水层(I承压水)、上第三系尚村组承压含水层(II承压水)及黄牛岭组承压含水层(Ⅲ承压水)。含水层以粗砂、砾砂为主,呈多层产出,层间为粘土、粉质粘土相隔,粘性土往往尖灭缺失或相变为砂性土,构成岩性天窗,使各含水层间的水力联系较为密切。北海组潜水含水层普遍高于海平面,承压含水层则向南、北方向延伸入海。补给来源主要是大气降水,还有运河季节性渗漏补给[5]。
改革开放后,随着北海市社会经济迅速发展,水资源需求激增,地下水超采现象十分普遍,致使地下水水位急剧下降,最终在采水点形成降落漏斗,引发海水入侵等环境问题[6]。另外,上世纪90年代中后期,由于海水高位养殖等工农业活动和污水入渗地下,加剧了沿海地带地下水环境的恶化。为监测北海市禾塘村水源地地下水受海水的影响,广州海洋地质调查局基于水源地开采现状以及物探测量成果,于2015年在北海市南部冯家江西岸钻取了系列分层地下水监测竖井,竖井断面编号为GXHSRQ-Ⅱ(以下简写为GX断面)(图1c,图1d)。系列分层地下水竖井由3组监测井构成,其分别设置在海水入侵区(GX1井组)、咸淡水过渡带(GX2井组)和未入侵区(GX3井组),三个井组的位置距离海岸线由近及远,每组井之间相距约500 m。同时,每组井又分别钻取3口不同含水层位的监测井,分别监测潜水层(深度5.1~9.1 m)、Ⅰ承压水层(深度48.6~60.6 m)和Ⅱ承压水层(深度60.6~69.6 m)。系列分层地下竖井采用井组编号与-1、-2、-3组合的形式表示,如GX-1井组的潜水井代号为GX-1-1)。监测井位置及钻孔深度见表1。
图1 广西北海近岸系列分层竖井位置与剖面示意Fig.1 Location and section of stratified vertical wells in coastal areas of Beihai,Guangxi
2 采样与分析
2.1 样品采集
分别于2017年10月,2018年11月在广西北海近岸GX断面的系列分层竖井进行地下水的监测和采样。在连续监测和采样前使用发电机带动潜水泵连续抽取井水20 min,以获取新鲜的地下水。随后使用有机玻璃采水器降至井底以上约1 m深处采集未与大气接触的地下水。采集的水样现场用孔径0.22μm的滤膜过滤,过滤后的水样分别进行5项营养盐(包括NH4-N、NO2-N、NO3-N、SiO3-Si和PO4-P)、氯离子、硫酸根离子、溶解无机碳(DIC)等含量的测试。用于营养盐测试的样品盛装于150 ml聚乙烯采样瓶,并加入10μl饱和氯化汞;用于Cl-、SO42-测试的样品置于150 ml聚乙烯采样瓶保存;用于DIC测试的样品从采样器底部的硅胶管中引入到10 ml的玻璃顶空瓶中,并溢出约采样瓶体积的2/3以去除气泡,加入2滴饱和氯化汞抑制微生物活动,所有采集的样品立即送至广西海洋地质调查研究院于4℃低温保存。样品采集结束后,将YSI水质多参数监测仪(CTD)放置在各水井离底约1 m深处,放置时间为24~25 h,测定的参数主要有深度、盐度和溶解氧等。
表1 近岸系列分层竖井位置与钻孔深度Tab.1 Locations and depths of stratified vertical wells
2.2 样品分析
营养盐分析采用AA3连续流动分析仪(Bran-Luebee,德国),其中NH4-N的测定采用靛酚蓝法,PO4-P的测定采用磷钼蓝法,SiO3-Si的测定采用硅钼蓝法,NO2-N和NO3-N的测定采用磺胺/N(1-萘乙酸)-乙二胺法,分析精度分别为99.69%、98.02%、99.05%、99.05%、99.41%。
地下水DIC及碳稳定同位素测定利用MAT 253同位素比质谱仪。前处理装置为Gas Bench II(美国Thermo Fisher公司),包括:GC-PAL自动进样器,PoraPlot Q(2 5m×0.32 mm),恒温样品盘(控制温度±0.1℃)。测定水中溶解无机碳同位素的反应流程为注入磷酸、充氦气、注入样品、反应达到平衡后进行样品测试。DIC的测试精度≤±0.15%;δ13C值的测试精度≤±0.15‰。
地下水Cl-、Na+、K+含量测定使用离子色谱法(离子色谱仪ICS-1500,美国Dionex公司);Ca2+、Mg2+离子的浓度采用乙二胺四乙酸二钠滴定法测定;HCO3-、等离子的浓度采用双指示剂中和法测定。
3 结果与讨论
3.1 近岸系列分层地下水的水化学类型及咸化特征
3.1.1水化学类型
地下水的水化学类型和阴阳离子当量可在一定程度上指示地下水来源及影响因素。图2(表示集中离子在溶液中的质量浓度百分比)显示GX1-1属于Cl-HCO3-Mg-Na型水,HCO3-毫克当量占阴离子毫克当量的27.16%;GX2-1属于HCO3-Ca型,HCO3-毫克当量占阴离子毫克当量的82.27%;GX3-1属于Cl-Na-Ca型,HCO3-毫克当量占阴离子毫克当量的22.88%,说明GX2-1和GX3-1具有地表水特征,而GX1-1明显兼具地表水和海水特征,表明其受到来自咸水的影响。GX1-2与GX2-2均为Cl-Na型水,且Cl-毫克当量和Na+毫克当量均分别占阴、阳离子毫克当量总量的90%和75%以上,HCO3-毫克当量仅为阴离子毫克当量总量的1%左右,指示GX1-2和GX2-2化学组成受氯钠离子控制,受海水影响明显。GX3-2为HCO3-Na-Ca型,说明其未受氯钠离子控制,井水未发生咸化。GX1-3与GX1-2,GX2-2化学组成特征基本一致,其水化学类型为Cl-Na型,且Cl-毫克当量与Na+毫克当量分别为阴阳离子毫克当量总量的90.72%与76.16%,具有与海水相似的离子组成;GX2-3与GX3-3分属HCO3-Na-Ca型和HCO3-Cl-Na型,与当地的II承压水特征一致[7]。
3.1.2 咸化特征
图3显示系列分层地下水竖井的I承压水与II承压水的水深、盐度、溶解氧含量的变化。GX1是最靠海方向的地下水井组,该井组GX1-2与GX1-3的水深日间变化均表现出明显的周期性波动,最大波动幅度近1 m,且其波动规律与近海潮汐涨落对应,表明其与近岸海水联系紧密且可能遭遇海水入侵。GX2和GX3井组的I承压和II承压水由于距离海洋较远,其水深日间变化微弱且缺乏周期性规律,水位波动变化与潮汐涨落无明显对应。I承压水中,GX1-2盐度值平均值为4.02;GX2-2盐度值在监测期间随时间波动明显,最大值3.83,最小值2.29。可见,GX1和GX2井组的I承压水盐度值均明显高于一般天然淡水,表明其受海水入侵发生咸化;GX3-2的盐度值较低,平均值0.08,说明其未受海水影响。GX1-3盐度值在监测期间盐度平均值则为6.70,为全部监测井最高,GX2-3和GX3-3盐度值则几乎无波动,其平均值分别为0.10和0.03,表明II承压水仅GX1井组被咸化。
Cl-化学行为保守,其含量高低可指示地下水的咸化程度。一般认为当地下水中的Cl-高于250 mg·L-1时,即认为地下水发生咸化[8]。在潜水层位中,GX1-1的Cl-含量最高,平均为386.48 mg·L-1,远大于250 mg·L-1,说明其被严重咸化;GX2-1和GX3-1的Cl-含量平均分布为46.60和28.93 mg·L-1,均远低于GX1-1,显示GX2和GX3井组的潜水未受咸水入侵影响。在I承压水层位中,GX1井组Cl-含量均值为643.87 mg·L-1,较GX1-1更高,表明其咸化更为严重;GX2-2平均Cl-含量为355.21 mg·L-1,亦高于250 mg·L-1警戒线;而GX3-2平均Cl-含量仅为11.21 mg·L-1,指示其未受咸水影响。在II承压水层位,GX1-3平均Cl-含量为1 271.64 mg·L-1,为全部监测层位最高,咸化最为严重,GX2-3平均Cl-含量为49.87 mg·L-1,该值与GX2-1井水相当,但要略高于淡水,说明可能受到较为微弱的咸水影响。GX3-3平均Cl-仅6.24 mg·L-1,基本代表淡水地下水组分,因此该井水可认为未受到咸化的影响。
图2 广西北海近岸系列分层竖井地下水化学类型piper图(单位:%)Fig.2 Piper trilinear nomograph of hydrochemical type of water in stratified vertical wells in coastal areas of Beihai,Guangxi(unit:%)
图3 广西北海近岸系列分层地下水井的水深、盐度、溶解氧质量浓度的变化Fig.3 Daily variation of depths,salinity and dissolved oxygen of stratified vertical wells in coastal areas of Beihai,Guangxi
3.1.3 井水咸化的影响因素
总体而言,GX系列监测井呈现咸化程度海边高,内陆低,潜水与II承压水咸化规模小,而I承压水咸化程度大的特征,同时各分层地下水又呈现各自独特的咸化特点。
(1)潜水:GX1井组的潜水咸化程度远远大于GX2和GX3井组潜水,主要原因受附近的高位养殖池塘的影响。野外调查发现,GX1井组附近存在一废弃的养殖鱼塘。20世纪90年代,北海地区流行高位咸水养殖,即通过人工在海平面之上的近海岸处开挖池塘,并灌注海水用于水产养殖。近年来,高位养殖方式虽被取缔,但大量咸水池塘仍存留于地表。由于北海地区地表盖层土质松散,且高位海水养殖池塘底部通常缺乏有效的防渗漏设施,故GX1井组的潜水极易因附近池塘水下渗而发生咸化[9]。同时,GX1井组潜水水化学特征指示其与地表水联系较密切,且此处潜水埋藏较浅,标程较高,因此难以发生海水入侵。
(2)I承压水:Cl-含量与盐度监测数据表明I承压水仅有GX3井组为淡水,GX1和GX2井组I承压水均被咸化,表明半咸水-淡水界线位于GX2和GX3井组之间,并且较潜水、II承压水更往内陆延伸。水化学证据显示GX1和GX2井组I承压水受氯钠离子控制,具有与海水相似的离子组成,说明I承压水咸水来源为单一海水,水深变化也反映了GX1井组I承压水与海水联系密切,同样指示其咸化是由海水入侵造成的。地下水是北海市的主要水源,自90年代以来,北海市地下水开采量急剧升高,造成地下水水位大幅下降,形成了数个明显的地下水降落漏斗,被认为是海水侵入陆地含水层的主要原因。
(3)II承压水:仅有位于靠海的GX1-3井水遭遇海水入侵,咸化明显,而GX2和GX3井组的II承压水均为淡水。GX1-3的水化学特征及水深变化规律与GX1、GX2井组的I承压水一致,说明其咸化的因素同样是由海水入侵引起的。与I承压水相比,II承压水的半咸水-淡水界限退至GX1和GX2井组之间,咸化范围相较I承压水有所缩小,这可能是由于II承压水埋藏深,开采难度大,开采量远低于I承压水,故地下水的水力梯度得以维持。
3.2 近岸系列分层竖井的地下水营养盐含量间分布与影响因素
3.2.1 SiO3-Si
在潜水层位,GX1-1的SiO3-Si含量最高,平均为241.08μmol·L-1;GX2-1和GX3-1的SiO3-Si含量平均分别为127.97和42.61μmol·L-1,这显示了在北海近岸系列分层地下水竖井的潜水层位中,SiO3-Si含量具有由内陆向海升高的趋势。在I承压水层位中,GX1-2、GX2-2和GX3-2井水的SiO3-Si含量分别为169.43,121.97和114.90μmol·L-1。可见,I承压水SiO3-Si含量亦显示出与潜水层位相似的变化趋势。在II承压水层位中,GX1-3平均SiO3-Si含量为272.08μmol·L-1,为全部井组的最大值。随着离海距离的增大,井水的SiO3-Si含量减小,GX2-3平均SiO3-Si含量降低至188.83μmol·L-1,GX3-3平均含量则减为121.94μmol·L-1。整体而言,广西北海近岸系列分层竖井的井水SiO3-Si含量表现出明显的“内陆低,海边高”的分布特征(图4)。
3.2.2 PO4-P
在潜水层位,GX1-1的PO4-P平均含量仅为0.12 μmol·L-1,部分水样中PO4-P浓度低于检出限;GX2-1的PO4-P平均含量为0.53μmol·L-1,较GX1-1略高,而GX3-1的PO4-P平均含量则更低(0.03μmol.L-1)。在I承压水层位,GX1-2和GX2-2的PO4-P平均含量均为0.04μmol·L-1,GX3-2的PO4-P平均含量为略高于前者,为0.07μmol·L-1。在II承压水层位,GX1,GX2,GX3井组的II承压水的PO4-P平均含量分别为0.07,0.05和0.04μmol·L-1,亦处于低位。总体而言,GX系列监测井的PO4-P含量整体偏低(图4),暗示地下水中存在强烈的PO4-P消除途径。
3.2.3 NH4-N
在潜水层位,GX1-1的NH4-N平均浓度为271.48 μmol·L-1,明显高于GX2-1(82.11μmol·L-1)和GX3-1(0.90μmol·L-1)的NH4-N平均含量,指示NH4-N含量在潜水层位呈现由海向陆递减趋势。在I承压水层位,GX1-2的NH4-N平均含量约为19.63μmol·L-1,较同井组潜水显著减少;GX2-2的NH4-N平均含量约为44.56μmol·L-1,GX3-2的NH4-N平均含量约为0.59 μmol·L-1。II承压水中NH4-N含量整体偏低,且各井组差别不大,其中GX1-3、GX2-3和GX3-3的NH4-N平均浓度分别为0.49、0.68和0.70μmol·L-1。总体而言,GX系列监测井中NH4-N含量整体表现出“表层多,下层少”的趋势(图4),暗示NH4-N可能主要来自地表。在含量最高且变化最明显的潜水层,NH4-N则表现出“海边高,内陆少”的分布规律,指示NH4-N可能更多向海边一侧累积。
3.2.4 NOx-N
NO2-N作为N循环的中间产物,其在地下水中的含量远低于NO3-N,又因为二者地球化学行为较一致,故可将NO3-N与NO2-N合并以NOx-N来表示。NOx-N含量在潜水层位中呈现出由陆向海递减的趋势,其中GX1-1的NOx-N平均含量为3.74μmol·L-1,GX2-1的NOx-N平均含量升至127.46μmol·L-1,GX3-1的NOx-N平均含量则高达604.22μmol·L-1。在I承压水层位,GX1-2的NOx-N平均含量为68.46μmol·L-1,GX2-2的NOx-N平均含量为70.97μmol·L-1,GX3-2的NOx-N平均含量为110.10μmol·L-1,I承压水亦存在NOx-N含量自内陆向海边的降低现象,但变化程度不如潜水层显著。在II承压水层位,GX1和GX2井组II承压水NOx-N含量差别不大,其平均含量分别为39.40和39.73μmol·L-1,略高于GX3-3的NOx-N的平均含
量(29.03μmol·L-1)。总体而言,除GX1-1外,GX系列监测井的NOx-N含量分布呈“表层高,下层低”趋势,在同层位,潜水层中NOx-N含量变化最为明显,呈现出“海边低,内陆高”的趋势,且I承压水与II承压水各井组间的NOx-N含量分布较为均衡,变化趋势不明显。
图4 广西北海近岸系列分层竖井的地下水营养盐随Cl-含量的变化关系Fig.4 Nutrients(PO4-P,SiO3-Si,NO x-N and NH4+-N)versus Cl-concentrations in stratified vertical wells in coastal areas of Beihai,Guangxi
3.3 营养盐含量变化的影响因素
3.3.1 SiO3-Si
广西北海近岸系列分层竖井的地下水中SiO3-Si含量表现出明显的“内陆低,海边高”的分布特征,此现象可能是咸水混入导致的“盐效应”造成。“盐效应”通常发生在海岸带地下河口环境中,因咸淡水混合而导致含硅矿物溶解加剧的现象[10]。由于地下水环境中的SiO3-Si主要来自于含硅矿物的溶解,因此“盐效应”会使得地下水环境中的SiO3-Si出现富集。“盐效应”产生的必要条件有二,首先为地下水环境中应混入大量的NaCl或碱金属、碱土金属离子。NaCl的混入会促进含水层中非定形硅溶解度的增加[11],而大量碱金属或碱土金属电解质离子混入,可显著提升石英的溶解速率[12]。在广西北海近岸系列分层竖井中,GX1井组由于靠近海岸线及高位养殖池塘,海水侵入及咸池塘水的下渗将大量氯钠离子,碱金属或碱土金属离子带入地下水环境,具体表现为GX1井组地下水溶解性总固体(total dissolved solids,TDS)的显著升高。在潜水层位,YSI水质多参数分析仪监测到GX1-1的TDS值高达3 063.4 mg·L-1,GX2-1和GX3-1则分别仅有300.9及69.6 mg·L-1;在I承压水层位,GX1-2和GX2-2均受到海水入侵,其TDS值分别为906.1及1 294.1 mg·L-1,GX3-2的TDS值则仅为60.2 mg·L-1;在II承压水层位,GX1-3的TDS值为4 322.1 mg·L-1,为全部监测井最高,GX2-3和GX3-3的TDS值则分别仅为89.4及26.1 mg·L-1。“盐效应”产生的第二个必要条件为当地含水层中需含有大量的硅质矿物,以作为地下水中SiO3-Si稳定输入源。据相关水文地质资料显示,GX系列监测井所主要钻探的下更新统湛江组(Q1z)地层与上新统尚村组(N2sh)地层主要岩性均为砂砾与粘土互层[5],富含石英、长石等高硅矿物,具备向地下水中供给大量SiO3-Si的潜力。
3.3.2 PO4-P
近岸系列分层竖井地下水中PO4-P分布特征是PO4-P含量普遍偏低,但变化趋势不明显。一般认为,此现象的原因主要是铁锰氧化物或其他颗粒有效吸附地下水中溶解态的PO4-P,导致PO4-P从水体中被吹扫清除。在可以吸附PO4-P的颗粒中,又以铁氧化物最常见且吸附能力最强,因此常以地层中出现的铁质氧化物集合作为PO4-P被强烈吸附作用的直接指示证据[13]。相关水文地质资料显示,广西北海近岸系列分层竖井主要钻探的地层为下更新统湛江组(Q1z),地层中赋存有丰富氧化铁质夹层[5],因此可作为PO4-P被铁氧化物吸附继而从水体中被清除的直接证据。同时,除GX1井组II承压水外,GX系列监测井溶解氧含量普遍较高(图3),有利于水中的二价铁离子的氧化,从而加剧对水中PO4-P的吸附清除作用。
3.3.3 NH4-N
NH4-N在GX系列监测井中的分布特征为“表层多,下层少”,指示NH4-N可能主要来自地表。在潜水层位,GX1-1中的NH4-N明显高于GX2-1和GX3-1,表明GX1井组附近的地表存在NH4-N输入源。3.1节提及,GX1-1咸化的主要原因是附近的高位养殖池塘水的下渗。野外调查发现,高位养殖池塘因富含有机质而发生了严重的富营养化,诱发藻类过盛生长而使水体变绿,池塘底泥则发黑变臭。因此,池塘水及底泥中的有机质发生分解,可将有机氮转化NH4-N并随池塘水下渗至潜水层,导致附近潜水井中的NH4-N含量增大。
3.3.4 NOx-N
近岸系列分层竖井地下水中的NOx-N的分布特征与NH4-N相似,均为“表层多,下层少”,指示GX系列分层竖井地下水中N可能来自于地表,但与NH4-N相反的是,NOx-N在潜水层中呈现出自内陆至海边削减趋势。海岸带地区微生物活动强烈,NOx-N的去除作用多与微生物参与的反硝化作用有关。反硝化作用是一种由反硝化菌参与的将硝酸盐中的N通过一系列中间产物还原为N2或N2O的生物地球化学过程,此过程通常以有机碳作电子供体,并生成DIC[14]。为了追溯地下水DIC的来源,可将DIC浓度与其对应的DIC浓度和碳同位素值的乘积做线性回归,所得的斜率即为反映“添加”的碳的δ13C值[15]。图5a显示潜水层“添加”的DIC同位素值为-16.84%。一般认为,地下水DIC的主要来源为大气、海水以及有机质的矿化分解,其中现代大气的δ13C值约为-7.5‰,表层海水的δ13C值约为0~1‰,陆源有机质δ13C约为-25‰~-30‰[16]。根据DIC端元同位素特征,可推断地下水DIC主要来自陆源有机质的分解作用。同时,潜水组井中NOx-N含量与DIC含量呈反向变化趋势(图5b),暗示DIC很大可能是由有机质氧化与NOx-N还原过程产生的,即反硝化过程有关。
图5 北海近岸系列分层竖井中潜水层DIC与DIC×δ13 C-DIC及NO x-N与DIC的关系Fig.5 DIC versus DIC×δ13 C-DIC,and NO x-N versus DIC in unconfined aquifer in coastal areas of Beihai,Guangxi
3.4 地下水排泄携带的营养盐对近海生态影响
只要含水层与海洋有水力连接,且水位在海平面以上,就会存在地下淡水向海洋排泄,同时将其中的营养盐物质输入海洋。一般认为,这些营养盐物质进入海洋后,将会对近海海洋环境产生两方面的影响,即促进海洋初级生产力的提升与影响海水中的营养物质平衡。
海洋富营养化与赤潮的暴发是海洋初级生产力过度升高引发的环境问题。北海附近海域曾多次遭受赤潮侵袭,Kaiser等[17]认为以南流江为代表的地表径流是北海海域营养盐的主要来源,也是造成该海域赤潮的主要原因。但调查发现,南流江河水中的部分营养盐浓度远低于近岸系列分层竖井地下水中的浓度,例如,南流江河水中的SiO3-Si含量约47μmol·L-1,而GX系列监测井GX1井组潜水、I承压水、II承压水中SiO3-Si平均含量分别为241.08,169.43和272.08μmol·L-1,远远高于河水中的含量;南流江河水中的NH4-N含量约47μmol·L-1,而GX系列监测井GX1井组潜水,I承压水中NH4-N平均含量分别为271.48和82.11μmol·L-1,亦比河水高出很多。尽管没测定采样区附近的海底地下水排泄通量,但Chen等[18]利用222Rn示踪剂测得采样区附近北部湾茅尾海SGD排泄通量为2.6×107~4.9×107m3·d-1,远大于当地河流得输入通量。故地下水高浓度营养盐及排泄通量是广西北海近海营养盐输入的重要途径,是诱发附近海域水体富营养化的重要推手。
地下水排泄不仅增加近海水体的营养盐浓度,同时也会改变水体N/P比值。海洋中浮游生物对N、P元素的利用比例约为16(Redfield比值),常规海水中N/P比略低于16。因此在正常情况下,浮游生物将处于弱N限制状态。但是由于P祛除的效率往往比N高,因此海岸带地下水的N/P比往往高于16[19]。GX系列监测井N/P比值的范围为141.01~4 326.04,平均值为595.85,远高于Redfield比值(图6)。长期大量高N/P比地下水排泄进入海洋,将可能迫使海洋中N/P比偏离Redfield比值,破坏海水中的营养物质平衡,促进浮游生物营养盐结构由N限制转向P限制。
图6 GX系列监测井氮磷比Fig.6 N/P ratio in groundwater of GX wells
4 结论
以广西北海市南部近海分层地下水竖井(分别位于海水入侵区、咸淡水分界线和未入侵区)为研究对象,通过现场水质监测和采样分析,主要结论如下:
(1)近岸系列竖井含水层整体上表现为靠近海边咸化明显,其中I承压水层咸化程度较大,而潜水与II承压水层咸化不明显;系列分层竖井含水层咸化程度差异主要受海水入侵及海水高位养殖的影响。
(2)不同含水层营养盐空间分布特征表现为SiO3-Si含量自内陆至海岸方向上升,各含水层的PO4-P含量皆偏低;潜水层中的NH4-N和NOx-N含量整体高于承压水,其中NH4-N在近海岸的竖井潜水层(GX1-1)含量最高,而NOx-N自内陆至海岸明显减小。地下含水层营养盐空间分布受不同过程影响,其中SiO3-Si含量变化可能的原因是咸水带来的碱金属离子造成含水层溶解加剧;地下水PO4-P含量低是受含水层铁氧化物吸附清除作用所致;NH4-N和NOx-N主要来源地表,潜水层NOx-N可能受反硝化作用影响,致使其自内陆至海岸呈减小趋势。
(3)研究区地下水含水层的N/P比平均值为595.85,高N/P比值地下水排泄进入海洋,将促进浮游生物营养盐结构由N限制转向P限制,势必会对海洋正常营养盐结构和浮游生物造成影响。
中国有18 000 km的绵长海岸带,随着沿海地区人口的不断增长和经济的快速发展,沿海海域赤潮发生的频率呈现出不断上升的态势。在探寻富营养化及赤潮的原因时,需要关注陆地地下水排泄对近海营养盐的输送及其对生态环境的影响。
作者贡献声明:
吴自军:设计了研究思路,制定了研究方案,指导了研究过程,并对论文进行了修改;王富康:撰写了论文并参与了样品采集和实验工作;崔振昂、高龙:为野外工作提供了协助;Aliya Zidan、朱虹霓、王一清:参与了样品采集工作。