珠江口上游水体缺氧：物理 - 生物地球化学耦合机制解析

珠江口上游水体缺氧：物理-生物地球化学耦合机制解析一、引言1.1研究背景与意义水体缺氧，通常指水体中溶解氧（DO）浓度低于2mg/L，是当前河口、近海面临的最为紧迫的全球性环境问题之一。当海水中氧气的消耗速率大于氧气的补充速率时，就会引发水体低氧。一旦近海河口水体溶解氧含量低至缺氧状态，生态状况便会急剧恶化，鱼、虾等众多水生生物将无法正常生活。同时，表层沉积物的氧化性环境遭到破坏，积聚在沉积物中的有毒、有害化学物质可能重新活化，释放到水体中，形成二次污染。在全球气候变化与高强度人类活动的共同影响下，近岸河口海域富营养化及其诱导的底层水季节性缺氧问题愈发严峻。国际上，如美国的切萨皮克湾，其缺氧现象严重影响了当地渔业资源和生态旅游；在欧洲的波罗的海，缺氧导致大量海洋生物死亡，生态系统失衡。我国部分河口及其近海水域也频繁出现缺氧现象，主要分布于长江与珠江河口及其邻近海域，夏季（6—8月）达到最盛，秋冬季消失。长江口的缺氧区域主要位于30.75°N~32°N，122.5°E~123.25°E附近海域，具有南北2个缺氧中心；珠江口缺氧程度相对较轻，缺氧区域主要在广州黄埔区河段至虎门水域和伶仃洋；辽河、钱塘江、海河等河口区域也时而出现零星的缺氧现象。珠江是世界上最重要的大河之一，其三角洲——珠江三角洲是世界上经济发展速度和人口增长速度最快、人类活动最频繁的区域之一。随着近几十年来珠江三角洲地区的快速城市化与工业化发展，大量的污染物通过河流或直接排放到珠江口，致使河口水质恶化，原有生态环境遭到破坏，富营养化、赤潮、缺氧等环境问题丛生。珠江口底层海水季节性的缺氧现象最早于1985年7月进行的广东省海岸带化学调查中得以发现，随后对珠江口进行的多次调查，均证实了这一特殊现象的存在，且有缺氧范围不断扩大，缺氧程度在逐年加剧的倾向。研究珠江口上游水体缺氧的形成机理具有重大意义。从环境保护角度看，明确缺氧形成机制有助于深入了解珠江口生态环境恶化的根源，为制定针对性的保护措施提供科学依据，从而有效遏制生态环境的进一步恶化，保护水生生物的生存环境，维护生物多样性。从可持续发展角度讲，珠江口地区经济发达，海洋经济在区域发展中占据重要地位。解决水体缺氧问题能够保障渔业、航运等产业的可持续发展，避免因生态破坏导致的经济损失，促进区域经济与生态环境的协调共进。1.2国内外研究现状在河口、近海地区水体缺氧研究领域，国外起步较早且成果丰硕。美国切萨皮克湾作为研究典范，科研人员运用长期监测与数值模拟手段，发现河流输入的大量营养盐促使浮游植物过度繁殖，其死亡后沉降分解大量消耗底层溶解氧，加之水体层化阻碍氧气补充，共同导致了严重的季节性缺氧现象。在欧洲，波罗的海的研究揭示了海水交换缓慢、营养盐富集以及人类活动排污等因素对缺氧形成的重要影响。国内对于河口、近海缺氧研究，在长江口和珠江口取得显著进展。长江口研究表明，夏季水体层化与高生物耗氧量是导致底层缺氧的关键因素，同时，长江冲淡水与外海海水的相互作用对缺氧区的分布和演变产生重要影响。珠江口方面，早期研究集中在溶解氧分布特征与环境要素关系上，通过多次调查明确了珠江口底层海水季节性缺氧现象，且发现其有范围扩大、程度加剧的趋势。有学者指出，陆源污染物输入、水体层化以及生物地球化学过程在珠江口缺氧形成中发挥重要作用。近期，吴加学教授研究团队发现珠江口缺氧区发生在富含有机质的泥质沉积带上，缺氧区水体内再悬浮会产生局部高浓度悬浮体，形成“次生的最大浑浊带”，潮周期内再悬浮颗粒引起的溶解氧耗氧通量与沉积物耗氧通量相当，在维持河口及近岸缺氧发育程度中起着不可忽略的重要作用。尽管国内外在河口、近海地区水体缺氧研究取得诸多成果，但在物理-生物地球化学耦合模式方面仍存在不足。现有研究多侧重于单一因素对缺氧的影响，对于物理过程（如潮汐、径流、水动力条件等）、生物过程（浮游植物生长、生物呼吸作用等）和地球化学过程（营养盐循环、有机物降解等）之间复杂的相互作用和耦合机制，尚未完全明晰。数值模拟研究中，模型对某些关键过程的参数化方案不够完善，导致模拟结果与实际情况存在一定偏差，难以准确预测缺氧的发生、发展和演变趋势。在珠江口上游水体缺氧研究中，针对该区域独特的地形地貌、复杂的河网水系以及高强度人类活动影响下的物理-生物地球化学耦合过程研究较少，缺乏系统、全面且深入的认识。1.3研究内容与方法本研究聚焦于珠江口上游水体缺氧问题，运用多学科交叉的研究思路，综合物理、生物地球化学等领域知识，构建耦合模式深入探究其形成机理，旨在为该区域生态环境保护和可持续发展提供科学依据。在研究内容上，本研究首先构建珠江口上游水体的物理-生物地球化学耦合模式。详细刻画该区域复杂的水动力条件，包括潮汐、径流等物理过程的精确模拟；全面描述生物地球化学过程，如营养盐循环、有机物降解以及浮游植物生长等关键环节。通过合理的参数设置和模型验证，确保模型能够准确再现珠江口上游水体的实际状态。基于所构建的耦合模式，对珠江口上游水体缺氧的形成机理展开深入模拟分析。定量研究物理过程与生物地球化学过程之间的相互作用机制，明确各因素在缺氧形成过程中的贡献大小。重点分析潮汐、径流对营养盐和溶解氧输运的影响，以及生物活动导致的耗氧过程如何在水动力条件下引发缺氧现象。本研究还将设计一系列数值实验，深入研究珠江口上游水体缺氧对不同环境条件变化的响应机制。系统探讨污染物输入变化、气候变化等因素对水体缺氧程度和范围的影响。通过模拟不同情景下的缺氧演变，为制定针对性的污染控制和生态保护策略提供科学参考。在研究方法上，本研究采用数值模拟方法，运用先进的数值计算技术求解物理-生物地球化学耦合模式的控制方程。通过合理的网格划分和时间步长设置，确保模拟结果的准确性和稳定性。利用高性能计算机集群进行大规模数值模拟，提高计算效率，缩短研究周期。同时，将数值模拟结果与实际观测数据进行对比验证，不断优化模型参数和结构，提高模型的可靠性。开展现场监测工作，在珠江口上游水体设置多个监测站位，进行长期、系统的水质、水文、生物等多参数监测。获取溶解氧、营养盐、浮游植物生物量等关键数据，为数值模拟提供基础数据支持，同时用于验证模型的模拟结果。采用先进的监测技术和设备，如原位传感器、自动监测浮标等，提高数据的准确性和实时性。通过数据分析方法，对现场监测数据和数值模拟结果进行深入分析。运用统计分析方法揭示各环境因素之间的相关性和变化规律，采用数据可视化技术直观展示水体缺氧的时空分布特征和演变趋势。利用敏感性分析方法确定影响水体缺氧的关键因素和敏感参数，为进一步研究提供方向。二、珠江口上游水体概况2.1地理位置与水系特征珠江口上游地处中国南部，涵盖广东、广西、云南、贵州等省份的部分地区，其地理位置介于东经110°～115°，北纬22°～25°之间。该区域东接珠江三角洲平原，西连云贵高原，北靠南岭山脉，南邻南海，是珠江水系与南海的重要连接地带，在区域生态系统和经济发展中占据关键地位。珠江是中国南方最大的水系，由西江、北江、东江及珠江三角洲诸河汇聚而成，流域面积广阔，达45.37万平方千米。西江作为珠江的主要支流，全长约2214千米，发源于云南省曲靖市马雄山，自西向东流经云南、贵州、广西、广东等省区，流域面积约35.31万平方千米。其上游为南盘江，在贵州省蔗香与北盘江汇合后称红水河，至广西壮族自治区石龙与柳江汇合后称黔江，在桂平与郁江汇合后称浔江，到梧州与桂江汇合后始称西江并流入广东，于三水市思贤滘与北江相通并汇入珠江三角洲网河区。西江水量充沛，年径流量约占珠江总径流量的77.8%，对珠江口上游水体的水量和水质起着决定性作用。北江发源于江西省信丰县小茅山，全长约573千米，流域面积5.84万平方千米。在韶关以北分为浈水和武江两支，正源为浈水。北江流经广东、湖南两省，沿途接纳众多支流，如连江、绥江等，在三水市思贤滘与西江相通，随后与东江及珠江三角洲诸河共同注入南海。北江的年径流量约为510多亿立方米，在珠江口上游的水系格局中扮演重要角色，其水流对珠江口的水动力条件和生态环境产生重要影响。东江发源于江西省寻乌县桠髻钵山，全长约562千米，流域面积2.7万平方公里。它自东北向西南流经广东省龙川县、和平县、东源县、源城区、惠城区、博罗县至东莞市石龙镇进入珠江三角洲，于广州市番禺区禺山处汇入珠江入海口——狮子洋。东江年径流量约257亿立方米，为珠江口上游提供了稳定的水源补给，对维持珠江口地区的生态平衡和水资源利用具有重要意义。在珠江口上游，西江、北江、东江及其众多支流相互交织，形成了复杂的河网水系。这些河流在径流、水位、含沙量等水文特征上存在差异，且受到地形、气候、人类活动等多种因素的影响。例如，西江上游流经云贵高原，地势起伏大，河流落差大，水能资源丰富；中游流经广西盆地，地势较为平坦，河网密度增加；下游进入珠江三角洲，河道分叉增多，水流分散。北江和东江的上游多山地丘陵，河流流速较快，侵蚀作用较强；下游进入平原地区，流速减缓，泥沙淤积。珠江口上游水系的交汇情况复杂多样。西江、北江在三水市思贤滘通过连通水道相互沟通，形成了一个有机的整体，这种连通使得西江和北江的水量、水质相互影响，共同塑造了珠江口上游的水动力条件和生态环境。东江在广州市番禺区禺山处与珠江口的其他水系交汇，其水流携带的营养物质和泥沙等对珠江口的生态系统产生重要影响。珠江三角洲地区河网密布，众多小河流纵横交错，进一步增加了水系交汇的复杂性。这些河流之间的相互作用，不仅影响着水体的物理、化学和生物过程，还对区域的水资源开发利用、防洪排涝、航运交通等产生深远影响。2.2水动力条件珠江口上游地区的水动力条件复杂，主要受到潮汐、径流、风生流等多种因素的共同影响。这些因素相互作用，塑造了该地区独特的水流特征，对水体的流动、混合和扩散过程产生深远影响，进而在珠江口上游水体缺氧的形成和发展中扮演着重要角色。潮汐是珠江口上游水动力的重要组成部分，该区域属于不正规半日潮型，平均潮差以磨刀门最小，为0.86米，东西两侧略大，伶仃洋湾头为1.35米，崖门为1.24米。潮差呈现出从河口湾的湾口向湾头增加，从各分流水道口门向上游递减的规律。在枯水期，潮区界距口门100-300公里，西江可达梧州-德庆，北江达芦苞-马房，东江达铁岗；而在洪水期，潮区界距口门40-70公里。潮流一般为往复流，枯水期潮流界距口门60-160公里，西江达三榕峡，北江至马房，东江至石龙；洪水期潮流界一般在口门附近，惟虎门水道可达广州。伶仃洋的涨落潮流轴线明显分异，落潮流路偏西，涨潮流路偏东。潮汐的周期性涨落使得水体在河口地区做往复运动，这种运动促进了水体的混合和物质的输运。在涨潮过程中，外海高盐水体被带入河口，与淡水混合，改变了水体的盐度和密度分布；落潮时，河口内的水体又被带出，将陆源物质输送到外海。这种周期性的物质交换对珠江口上游水体的化学组成和生态环境产生重要影响。例如，潮汐作用可以将海洋中的营养物质带入河口，为浮游植物的生长提供充足的养分，促进浮游植物的繁殖，进而影响水体的溶解氧含量。径流是珠江口上游水动力的另一个关键因素，珠江年平均流量约1万立方米/秒，年径流总量3457.8亿立方米，4-9月的径流量占全年的80%。西江作为珠江的主要支流，其径流量对珠江口上游的水动力条件起着决定性作用，年径流量约占珠江总径流量的77.8%。北江和东江也为珠江口上游提供了稳定的水源补给。径流的大小和变化直接影响着河口地区的水流速度和流向。在洪水期，径流量增大，河水的流速加快，携带的泥沙和营养物质增多，河水能够更有力地向外海推进，抑制咸潮上溯，扩大淡水的影响范围。此时，河流的输沙量增加，可能导致河口地区的泥沙淤积，改变河道形态和水动力条件。而在枯水期，径流量减小，河水的流速减缓，携带的物质减少，海水更容易倒灌进入河口，引发咸潮，导致河口地区水体盐度升高。咸潮的发生不仅影响河口地区的水资源利用，还可能对水生生物的生存环境造成威胁。径流还会影响水体的混合和扩散过程，较强的径流能够增强水体的紊动，促进水体的混合，使营养物质、溶解氧等在水体中更均匀地分布。风生流也是珠江口上游水动力的重要影响因素之一，该地区地处南亚热带海洋性季风气候区，季风对风生流的形成和变化起着主导作用。在夏季，盛行西南季风，风力较强，使得表层水体产生向东北方向的风生流；在冬季，盛行东北季风，风生流方向则大致为西南。风的作用可以直接推动水体运动，增加水体的紊动和混合。当风力较大时，水面产生波浪，波浪的破碎和传播会加剧水体的垂直混合，使表层富含溶解氧的水体与底层水体进行交换，有助于提高底层水体的溶解氧含量。风还可以影响河口地区的水流形态和流场分布。例如，在河口附近，风的作用可能导致水流的辐合或辐散，影响物质的聚集和扩散。当风与河口的地形和水流相互作用时，可能会形成局部的环流或涡流，这些环流和涡流对营养盐、污染物等物质的输运和分布产生重要影响。在某些情况下，风生流可能会将外海的高盐水体或污染物带入河口，或者将河口内的物质带出，从而改变河口地区的生态环境。2.3水质现状珠江口上游水体的水质状况复杂且面临严峻挑战，主要水质指标呈现出显著的变化和特征，充分反映了水体缺氧问题的严重性。溶解氧（DO）作为衡量水质的关键指标之一，在珠江口上游水体中分布不均且存在季节性变化。夏季，水体温度升高，微生物活动增强，导致溶解氧的消耗速率加快。同时，由于水体层化现象明显，表层水体与底层水体交换不畅，底层水体难以获得充足的氧气补充，使得底层溶解氧含量急剧下降，经常出现缺氧现象。研究表明，在夏季高温时段，珠江口上游部分区域底层溶解氧浓度可低至2mg/L以下，严重威胁水生生物的生存。例如，在伶仃洋等区域，夏季底层水体缺氧范围较大，对渔业资源造成了明显影响，导致部分鱼类和贝类等生物的死亡和迁移。而在冬季，水温较低，微生物活动相对较弱，水体层化现象减弱，溶解氧含量相对较高，但在一些排污口附近和水流缓慢的区域，溶解氧含量仍然偏低，无法满足水生生物的正常需求。营养盐是影响水体富营养化和生态系统健康的重要因素。珠江口上游水体中的营养盐主要包括无机氮和磷酸盐等。由于受到陆源污染物排放、农业面源污染以及水产养殖等人类活动的影响，水体中的营养盐含量呈现出上升趋势。过量的营养盐输入导致水体富营养化，为浮游植物的大量繁殖提供了充足的养分，进而引发赤潮等生态灾害。据统计，近年来珠江口上游水体中的无机氮和磷酸盐含量多次超过国家海水水质标准的二类标准，部分区域甚至达到四类标准。在一些河口区域和海湾，由于营养盐的富集，浮游植物大量繁殖，形成了大面积的赤潮，不仅破坏了海洋生态平衡，还对渔业、旅游业等产业造成了巨大的经济损失。例如，2024年12月中旬开始，粤港澳大湾区东部海域接连发生赤潮，2024年12月30日，珠江口海域发生大面积赤潮，赤潮优势种为球形棕囊藻，相关地方政府部门启动赤潮灾害Ⅰ级应急响应。有机物是珠江口上游水体中的另一类重要污染物，主要来源于工业废水、生活污水以及农业废弃物等。这些有机物在水体中分解时会消耗大量的溶解氧，导致水体缺氧，同时还会产生一些有害物质，对水生生物和人体健康造成危害。珠江口上游部分区域水体中的化学需氧量（COD）和生化需氧量（BOD）等指标超过国家海水水质标准，表明水体中的有机物污染较为严重。在一些城市附近的河口和海湾，由于工业废水和生活污水的排放，水体中的有机物含量较高，水质发黑发臭，生态环境遭到严重破坏。三、物理-生物地球化学耦合模式构建3.1模式基本框架本研究构建的物理-生物地球化学耦合模式基于FORTRAN编程语言，以有限元海洋模式环境（FEMME）为模拟环境，该模式整合了物理模块与生物地球化学模块，旨在全面、准确地模拟珠江口上游水体的复杂过程。FEMME是一款功能强大的海洋模式环境，它采用有限元方法对海洋进行数值离散，能够灵活处理复杂的海岸线和地形变化。在本研究中，FEMME为整个耦合模式提供了基础的计算框架，确保了模拟的准确性和稳定性。通过FEMME，我们可以精确地描述珠江口上游水体的几何形状和边界条件，为后续的物理和生物地球化学过程模拟奠定坚实基础。例如，在处理珠江口复杂的河网水系时，FEMME能够根据实际地形进行细致的网格划分，准确捕捉水流在不同河道间的流动和交换。物理模块是耦合模式的重要组成部分，主要用于模拟珠江口上游水体的水动力过程，包括潮汐、径流、风生流等。潮汐模拟采用调和分析方法，通过引入多个主要分潮（如M2、S2、N2、K1等）来准确描述潮汐的周期性变化。在珠江口，M2分潮是主要的潮汐成分，其振幅和相位的准确模拟对于再现潮汐现象至关重要。通过对历史潮汐数据的分析和验证，确定了各分潮的调和常数，从而实现对潮汐的精确模拟。径流模拟考虑了西江、北江、东江等主要河流的流量变化，根据实测的流量数据和水位数据，建立了河流流量与水位的关系模型。在洪水期和枯水期，河流流量变化显著，该模型能够根据实时的水文数据，准确模拟径流对珠江口上游水体的影响。风生流模拟则基于风应力驱动方程，考虑了不同季节、不同风速和风向对水体的作用。在夏季，西南季风较强，通过风应力驱动方程计算得出的风生流能够准确反映表层水体向东北方向的运动趋势；在冬季，东北季风主导，风生流方向大致为西南，模拟结果与实际观测相符。生物地球化学模块主要描述水体中的生物地球化学过程，包括营养盐循环、有机物降解以及浮游植物生长等。营养盐循环过程考虑了氮、磷、硅等主要营养盐的输入、输出和转化。例如，在氮循环中，包括有机氮的矿化、氨氮的氧化、硝酸盐的还原等过程，通过建立相应的化学反应动力学方程来描述这些过程。有机氮在微生物的作用下发生矿化反应，转化为氨氮，该反应速率受到温度、微生物活性等因素的影响。有机物降解过程采用一级反应动力学模型，描述了有机物在微生物作用下的分解过程，以及由此产生的溶解氧消耗。在珠江口上游水体中，有机物主要来源于陆源输入和浮游植物的代谢产物，这些有机物在微生物的作用下逐渐分解，消耗大量溶解氧，导致水体缺氧。浮游植物生长过程考虑了光照、温度、营养盐等因素对浮游植物生长速率的影响，采用米氏方程来描述浮游植物对营养盐的摄取。光照强度和温度直接影响浮游植物的光合作用和生长代谢，营养盐浓度则是浮游植物生长的关键限制因素之一。通过米氏方程，可以准确计算浮游植物在不同环境条件下的生长速率和生物量变化。物理模块和生物地球化学模块之间通过物质和能量的交换实现耦合。在水动力过程的作用下，营养盐、溶解氧、浮游植物等物质在水体中发生输运和扩散，从而影响生物地球化学过程。潮汐和径流的变化会导致营养盐的分布发生改变，进而影响浮游植物的生长和繁殖。生物地球化学过程中的有机物降解和浮游植物呼吸作用会消耗溶解氧，改变水体的溶解氧分布，反过来又影响水动力过程。当水体中溶解氧含量降低时，会抑制一些好氧微生物的活性，从而影响有机物的降解速率，进一步影响水体的化学组成和生态环境。3.2物理模块3.2.1平流-扩散输送平流和扩散过程在水体中物质输送方面发挥着关键作用，对珠江口上游水体缺氧的形成和发展有着重要影响。平流过程是指物质随着水流的整体运动而发生的输送，它主要受到水流速度和流向的控制。在珠江口上游，河流径流和潮汐水流是平流输送的主要驱动力。例如，西江、北江和东江的径流携带大量的营养物质、溶解氧以及其他污染物等，以平流的方式向下游输送。在洪水期，河流径流量增大，水流速度加快，平流输送的能力增强，能够将更多的物质输送到下游地区。潮汐水流则呈现周期性的涨落，在涨潮时，外海的海水携带物质向河口内推进；落潮时，河口内的水体和物质又被带出。这种潮汐引起的平流作用使得物质在河口地区进行往复输送，对水体中物质的分布产生重要影响。扩散过程是由于物质浓度梯度的存在，使得物质从高浓度区域向低浓度区域迁移的现象。它主要包括分子扩散和紊动扩散两种形式。分子扩散是由分子的热运动引起的，在水体中，分子的无规则运动使得物质在微观层面上发生扩散。然而，在自然水体中，紊动扩散通常更为重要。紊动扩散是由于水体的紊流运动导致的，紊流中的涡旋和湍流使得物质在更大尺度上发生混合和扩散。在珠江口上游，水流的紊动程度受到多种因素的影响，如河道的地形、水流速度、潮汐等。在河道狭窄、水流湍急的区域，紊流强度较大，扩散作用更为明显。河口地区的潮汐运动也会增强水体的紊动，促进物质的扩散。当潮汐涨落时，水流的速度和方向不断变化，形成复杂的紊流场，使得物质能够更均匀地分布在水体中。平流-扩散输送过程可以用数学方程进行描述。在笛卡尔坐标系下，对于二维平面的物质输运，其控制方程如下：\frac{\partialC}{\partialt}+u\frac{\partialC}{\partialx}+v\frac{\partialC}{\partialy}=\frac{\partial}{\partialx}(D_x\frac{\partialC}{\partialx})+\frac{\partial}{\partialy}(D_y\frac{\partialC}{\partialy})+S其中，C表示物质的浓度，t表示时间，u和v分别是x和y方向上的水流速度分量，D_x和D_y分别是x和y方向上的扩散系数，S表示源汇项，用于描述物质的产生或消耗。在实际应用中，需要根据具体的研究区域和问题，确定合适的边界条件和初始条件，以求解该方程，从而得到物质在水体中的浓度分布和输运过程。例如，在珠江口上游水体缺氧研究中，可以将溶解氧、营养盐等物质的浓度作为C，通过测量或估算得到水流速度u和v、扩散系数D_x和D_y以及源汇项S的具体值，然后利用数值方法求解上述方程，模拟这些物质在水体中的平流-扩散输送过程，进而分析它们对水体缺氧的影响。3.2.2水气交换水气交换过程中氧气的传输机制及影响因素对于理解珠江口上游水体的溶解氧分布和变化至关重要，直接关系到水体缺氧问题的研究。在水气界面，氧气的传输主要通过分子扩散和紊流扩散两种方式进行。分子扩散是基于氧气分子在水体和大气中的浓度差，从高浓度区域向低浓度区域的微观迁移过程。当水体中溶解氧浓度低于大气中的氧分压时，氧气分子会从大气向水体中扩散；反之，当水体中溶解氧浓度过高时，氧气会从水体向大气逸出。在平静的水面条件下，分子扩散在氧气传输中起主要作用。然而，在自然水体中，紊流扩散更为常见且重要。紊流扩散是由于水体表面的紊流运动，使得氧气在更大尺度上与水体进行混合和交换。风力、水流等因素会引起水体表面的紊动，形成各种尺度的涡旋和湍流，这些涡旋和湍流将表层富含氧气的水体与下层水体进行混合，促进了氧气的传输。在珠江口上游，夏季的西南季风和冬季的东北季风都会对水体表面产生较大的风力作用，增强水体的紊动，从而加速氧气的水气交换过程。影响水气交换过程中氧气传输的因素众多，主要包括水动力条件、水温、大气压力和风速等。水动力条件是影响氧气传输的重要因素之一，河流的流速、潮汐的涨落以及风生流等都会改变水体的紊动程度，进而影响氧气的交换速率。在流速较快的河流区域，水体的紊动强烈，氧气的交换效率较高；而在水流缓慢的区域，紊动较弱，氧气交换相对较慢。潮汐的周期性变化也会导致水体的水位和流速发生改变，对氧气的传输产生影响。在涨潮时，外海高溶解氧的海水进入河口，与河口内的水体混合，增加了水体的溶解氧含量；落潮时，河口内的水体向外海流出，溶解氧含量可能发生变化。水温对氧气的溶解度和传输速率有着显著影响。一般来说，水温升高，氧气在水中的溶解度降低。在珠江口上游，夏季水温较高，氧气的溶解度相对较低，这使得水体对氧气的容纳能力下降。水温还会影响水体的密度和粘性，进而改变水体的紊动特性，影响氧气的传输。较高的水温可能导致水体层化现象加剧，阻碍了上下层水体之间的混合，使得底层水体难以获得充足的氧气补充，增加了水体缺氧的风险。大气压力的变化也会对水气交换产生影响。当大气压力降低时，水体表面的氧气分压相对升高，有利于氧气从水体向大气逸出；反之，大气压力升高时，氧气更倾向于从大气向水体中溶解。在不同的天气条件下，大气压力会发生变化，从而影响珠江口上游水体的溶解氧含量。例如，在台风等天气系统影响下，大气压力会急剧下降，可能导致水体中溶解氧的逸出增加，对水生生物的生存环境造成影响。风速是影响水气交换中氧气传输的另一个关键因素。较大的风速会增强水体表面的紊动，增加水气界面的面积和氧气的交换速率。在珠江口上游，风速的季节性变化明显，夏季的西南季风和冬季的东北季风带来不同强度的风力。在风速较大的季节，氧气的水气交换过程更为活跃，能够有效地补充水体中的溶解氧。当风速较小时，水体表面相对平静，氧气的交换速率降低，可能导致水体中溶解氧的分布不均匀，局部区域出现缺氧现象。3.3生物地球化学模块3.3.1有机物耗氧降解在珠江口上游水体中，有机物的耗氧降解主要由微生物介导，是一个复杂且对水体溶解氧平衡至关重要的过程。水体中的有机物来源广泛，包括陆源输入，如工业废水、生活污水和农业面源污染等带来大量的有机污染物；水生生物的代谢产物，浮游植物在生长和繁殖过程中会产生一些有机物质；以及生物残体的分解，当水生生物死亡后，其遗体逐渐分解，释放出有机物。这些有机物在微生物的作用下发生一系列的生化反应，最终降解为简单的无机物。微生物在有机物耗氧降解过程中扮演着关键角色。细菌、真菌等微生物能够分泌各种酶，如淀粉酶、蛋白酶、脂肪酶等，这些酶能够将复杂的有机物分解为小分子物质，以便微生物吸收利用。在有氧条件下，好氧微生物利用氧气将有机物氧化分解为二氧化碳和水等无机物，同时释放出能量，用于自身的生长和代谢。在这个过程中，微生物通过呼吸作用摄取氧气，将有机物中的碳氧化为二氧化碳，氢与氧结合生成水。而在缺氧或无氧条件下，厌氧微生物则利用其他电子受体，如硝酸盐、硫酸盐等，进行无氧呼吸，将有机物分解为甲烷、硫化氢等产物。有机物耗氧降解过程可以用一级反应动力学模型来描述：\frac{dC}{dt}=-kC其中，C表示有机物的浓度，t表示时间，k是降解速率常数，它受到多种因素的影响。温度是影响降解速率常数的重要因素之一，一般来说，温度升高，微生物的活性增强，降解速率加快。在珠江口上游水体中，夏季水温较高，有机物的降解速率相对较快，导致溶解氧的消耗也相应增加。有机物的种类和性质也会对降解速率产生影响。简单的糖类、蛋白质等有机物相对容易被微生物降解，而复杂的多环芳烃、纤维素等则较难降解。微生物的种类和数量同样影响着降解速率，不同种类的微生物对有机物的降解能力不同，微生物数量的增加通常会加快降解速率。当水体中存在大量的优势降解菌时，有机物的降解效率会显著提高。3.3.2硝化作用硝化作用是水体氮循环中的关键环节，对珠江口上游水体的溶解氧含量和生态系统平衡有着重要影响。在这个过程中，氨氮（NH_4^+）在硝化细菌的作用下逐步转化为硝酸盐（NO_3^-）。硝化作用主要分为两个阶段：亚硝化阶段和硝化阶段。在亚硝化阶段，亚硝化细菌，如硝化杆菌属（Nitrosomonas）等，将氨氮氧化为亚硝酸盐（NO_2^-），其反应方程式为：NH_4^++\frac{3}{2}O_2\xrightarrow[]{亚硝化细菌}NO_2^-+H_2O+2H^+这个反应过程需要消耗氧气，每氧化1mol氨氮需要消耗1.5mol氧气。在珠江口上游水体中，氨氮主要来源于工业废水、生活污水的排放以及农业面源污染等。当这些含氨氮的污水进入水体后，亚硝化细菌利用水中的溶解氧将氨氮氧化为亚硝酸盐。亚硝酸盐在水体中积累，如果不能及时被进一步氧化，会对水生生物产生毒性，影响其生长和生存。在硝化阶段，硝化细菌，如硝化球菌属（Nitrobacter）等，将亚硝酸盐进一步氧化为硝酸盐，反应方程式为：NO_2^-+\frac{1}{2}O_2\xrightarrow[]{硝化细菌}NO_3^-此反应同样需要消耗氧气，每氧化1mol亚硝酸盐需要消耗0.5mol氧气。经过这两个阶段，氨氮最终转化为相对稳定的硝酸盐。硝酸盐是一种相对无毒的含氮化合物，在水体中可以被浮游植物等生物吸收利用，参与到生态系统的物质循环中。硝化作用的进行受到多种因素的影响。溶解氧是硝化作用的关键因素之一，充足的溶解氧是硝化细菌进行有氧呼吸和氧化反应的必要条件。在珠江口上游水体中，如果溶解氧含量不足，硝化作用会受到抑制，导致氨氮和亚硝酸盐的积累。温度对硝化作用也有显著影响，硝化细菌的最适生长温度一般在25-30℃之间。在夏季，珠江口上游水体温度较高，有利于硝化作用的进行；而在冬季，水温较低，硝化细菌的活性降低，硝化作用速率减慢。水体的pH值也会影响硝化作用，硝化细菌适宜在中性至弱碱性的环境中生长，当水体pH值偏离这个范围时，硝化作用会受到影响。如果水体受到酸性废水的污染，pH值降低，硝化细菌的活性会受到抑制，进而影响硝化作用的进行。3.3.3初级生产初级生产是指浮游植物通过光合作用将光能转化为化学能，合成有机物质并产生氧气的过程，这一过程对珠江口上游水体的溶解氧平衡和生态系统的物质循环起着至关重要的作用。浮游植物是珠江口上游水体初级生产的主要承担者，包括硅藻、绿藻、甲藻等多种类群。这些浮游植物细胞内含有叶绿素等光合色素，能够吸收光能。在光照条件下，浮游植物利用光能将二氧化碳和水转化为有机物质（如葡萄糖等碳水化合物），同时释放出氧气。其光合作用的总反应式可以表示为：6CO_2+6H_2O\xrightarrow[]{光能、叶绿ç´

}C_6H_{12}O_6+6O_2这个过程不仅为浮游植物自身的生长和繁殖提供了能量和物质基础，还向水体中释放大量的氧气，是水体溶解氧的重要来源之一。初级生产过程受到多种因素的影响。光照是影响浮游植物光合作用的关键因素之一，光照强度和光照时间直接决定了浮游植物能够吸收的光能数量。在珠江口上游水体中，不同季节和不同天气条件下的光照强度和光照时间存在差异。夏季，光照时间长且强度高，有利于浮游植物进行光合作用，初级生产速率较高；而在冬季，光照时间缩短且强度减弱，浮游植物的光合作用受到一定限制，初级生产速率相对较低。此外，水体的透明度也会影响光照在水体中的穿透深度，进而影响浮游植物的光合作用。当水体中悬浮颗粒物较多，透明度降低时，光照难以到达深层水体，会限制浮游植物在水体中的分布和光合作用的进行。营养盐是浮游植物生长和进行光合作用的重要物质基础，氮、磷、硅等营养盐的浓度对初级生产有着显著影响。在珠江口上游水体中，由于受到陆源污染物排放、农业面源污染等因素的影响，营养盐含量存在时空变化。当水体中营养盐充足时，浮游植物能够获得足够的养分进行生长和繁殖，初级生产速率提高；反之，当营养盐缺乏时，浮游植物的生长会受到限制，初级生产速率降低。例如，氮是浮游植物合成蛋白质和核酸的重要元素，磷参与浮游植物的能量代谢和物质合成过程。如果水体中氮、磷含量不足，会导致浮游植物生长缓慢，光合作用效率下降，进而影响初级生产。温度对浮游植物的生理活动和生长繁殖也有重要影响。不同种类的浮游植物对温度的适应范围不同，但一般来说，适宜的温度能够促进浮游植物的生长和光合作用。在珠江口上游水体中，夏季水温较高，有利于浮游植物的生长和繁殖，初级生产较为旺盛；而在冬季，水温较低，浮游植物的生长和代谢活动减缓，初级生产速率降低。此外，温度还会影响浮游植物的种类组成和群落结构，不同温度条件下优势浮游植物种类可能发生变化，进而影响初级生产的效率和生态系统的功能。3.3.4反硝化作用反硝化作用是水体氮循环中的重要环节，对珠江口上游水体的溶解氧和氮素平衡有着重要影响。在这个过程中，反硝化细菌在缺氧或无氧条件下，将硝酸盐（NO_3^-）逐步还原为氮气（N_2），从而实现氮素从水体中的去除。反硝化细菌是一类能够利用硝酸盐作为电子受体进行无氧呼吸的微生物，常见的反硝化细菌包括假单胞菌属（Pseudomonas）、芽孢杆菌属（Bacillus）等。其反硝化过程一般分为以下几个步骤：首先，硝酸盐在硝酸还原酶的作用下被还原为亚硝酸盐（NO_2^-），反应式为：NO_3^-+2H^++2e^-\xrightarrow[]{硝酸还原酶}NO_2^-+H_2O接着，亚硝酸盐在亚硝酸还原酶的作用下进一步被还原为一氧化氮（NO）：NO_2^-+2H^++e^-\xrightarrow[]{亚硝酸还原酶}NO+H_2O一氧化氮再在一氧化氮还原酶的作用下被还原为一氧化二氮（N_2O）：2NO+2H^++2e^-\xrightarrow[]{一氧化氮还原酶}N_2O+H_2O最后，一氧化二氮在氧化亚氮还原酶的作用下被还原为氮气：N_2O+2H^++2e^-\xrightarrow[]{氧化亚氮还原酶}N_2+H_2O反硝化作用对水体溶解氧的影响较为复杂。一方面，反硝化作用在缺氧或无氧条件下进行，这意味着水体中溶解氧含量较低时，反硝化作用才能够顺利发生。因此，从这个角度来看，反硝化作用与水体缺氧环境相互关联。当水体中存在大量的硝酸盐，且溶解氧含量因其他生物地球化学过程（如有机物耗氧降解、硝化作用等）而降低时，反硝化细菌会利用硝酸盐进行无氧呼吸，从而实现氮素的去除。另一方面，反硝化作用过程本身并不直接消耗溶解氧，反而在一定程度上减少了水体中因硝化作用等过程对溶解氧的消耗。因为硝化作用需要消耗大量的溶解氧，而反硝化作用将硝化作用产生的硝酸盐还原为氮气，减少了后续硝化作用的底物，间接缓解了溶解氧的消耗压力。反硝化作用受到多种因素的影响。溶解氧是反硝化作用的关键影响因素之一，只有在溶解氧含量较低的环境中，反硝化细菌才会利用硝酸盐作为电子受体进行无氧呼吸。在珠江口上游水体中，底层水体由于与大气的氧气交换相对困难，且受到有机物耗氧降解等过程的影响，溶解氧含量较低，为反硝化作用提供了适宜的环境。硝酸盐浓度也是影响反硝化作用的重要因素，当水体中硝酸盐浓度较高时，反硝化细菌有足够的底物进行反硝化反应，反硝化作用速率会相应提高。如果水体中硝酸盐浓度过低，反硝化作用会受到限制。此外，温度、pH值、有机碳源等因素也会对反硝化作用产生影响。适宜的温度和pH值能够促进反硝化细菌的生长和代谢，提高反硝化作用速率。有机碳源是反硝化细菌进行反硝化反应的能源物质，充足的有机碳源能够为反硝化作用提供动力，促进硝酸盐的还原。3.4模型参数设置与验证根据珠江口上游水体的具体特点，对模型参数进行了细致设置。在物理模块中，水平方向的扩散系数D_x和D_y分别设置为100m^2/s和100m^2/s，这是基于对该区域水动力特性的研究以及相关文献的参考确定的。在一些水流较为湍急、紊动较强的区域，适当增大扩散系数，以更准确地反映物质的扩散过程；而在水流相对平缓的区域，则相应减小扩散系数。垂向扩散系数D_z设置为0.01m^2/s，考虑到珠江口上游水体在垂向上的混合相对较弱，该值能够合理地模拟垂向物质输运。在生物地球化学模块中，有机物降解速率常数k根据不同类型有机物的特性进行设置。对于易降解的有机物，如糖类、蛋白质等，降解速率常数设置为0.1d^{-1}；对于难降解的有机物，如多环芳烃等，降解速率常数设置为0.01d^{-1}。硝化作用速率常数根据硝化细菌的生长特性和环境条件进行调整，在适宜的温度和溶解氧条件下，亚硝化阶段的速率常数设置为0.05d^{-1}，硝化阶段的速率常数设置为0.03d^{-1}。浮游植物生长速率常数考虑光照、温度、营养盐等因素的影响，通过实验数据和经验公式进行确定。在光照充足、营养盐丰富的夏季，浮游植物生长速率常数设置为0.3d^{-1}；在光照较弱、营养盐相对不足的冬季，生长速率常数设置为0.1d^{-1}。为了确保模型的准确性和可靠性，利用2005年冬夏航次在珠江口上游水体的实测数据对模型进行校正和检验。在冬季航次中，实测了多个站位的水温、盐度、溶解氧、营养盐等参数，将这些数据与模型模拟结果进行对比。通过调整模型参数，使得模拟结果与实测数据在水温、盐度等方面达到较好的吻合。在夏季航次中，重点对溶解氧和营养盐的模拟结果进行校正。由于夏季水体层化现象明显，溶解氧分布复杂，通过优化垂向扩散系数和生物地球化学过程的参数，使模型能够准确模拟溶解氧的垂向分布和水平变化。通过对2005年冬夏航次数据的验证，评估了模拟结果与实测数据的吻合程度。以溶解氧为例，将模拟得到的溶解氧浓度与实测值进行对比，计算两者之间的均方根误差（RMSE）和平均相对误差（MRE）。在冬季，模拟结果与实测值的RMSE为0.5mg/L，MRE为10\%；在夏季，RMSE为0.8mg/L，MRE为15\%。这表明模型在不同季节对溶解氧的模拟具有较高的准确性，能够较好地再现珠江口上游水体溶解氧的实际分布情况。对于营养盐等其他参数，模拟结果与实测数据也具有较好的一致性，进一步验证了模型的可靠性。四、珠江口上游水体缺氧形成机理分析4.1冬季水体缺氧情况模拟与分析4.1.1溶解氧分布特征利用构建的物理-生物地球化学耦合模式，对冬季珠江口上游水体溶解氧分布进行模拟，结果清晰地展示了其独特的空间格局。在珠江口上游，溶解氧浓度呈现出明显的区域差异。河口附近区域，由于受到河流径流和潮汐的共同影响，溶解氧浓度相对较高。西江、北江和东江的径流携带大量溶解氧丰富的淡水注入河口，同时潮汐的涨落促进了水体的混合，使得河口附近水体能够与大气进行充分的氧气交换，从而维持较高的溶解氧水平。在一些主要河口处，溶解氧浓度可达6mg/L以上。随着向远离河口的方向延伸，溶解氧浓度逐渐降低。在珠江口外海区域，水体受到海洋环境的影响较大，与大气的氧气交换相对较弱，且水体中有机物的降解和生物呼吸作用会消耗一定量的溶解氧，导致溶解氧浓度下降。在距离河口较远的海域，溶解氧浓度可降至4mg/L左右。在垂直方向上，冬季珠江口上游水体溶解氧分布也存在明显差异。表层水体由于与大气直接接触，能够充分进行氧气交换，溶解氧浓度较高。在大多数区域，表层溶解氧浓度可达5mg/L以上。而底层水体由于受到水体层化和有机物沉积的影响，溶解氧浓度相对较低。水体层化现象在冬季虽然相对较弱，但仍然存在，这使得底层水体与表层水体之间的交换受到限制，难以获得充足的氧气补充。底层有机物的分解也会消耗大量溶解氧，进一步降低底层溶解氧浓度。在一些区域，底层溶解氧浓度可低至3mg/L以下，存在局部缺氧的风险。这种溶解氧分布特征的形成原因主要与水动力条件、生物地球化学过程以及水体层化等因素密切相关。水动力条件方面，河流径流和潮汐的作用不仅影响了水体的流动和混合，还决定了溶解氧的输运和扩散路径。河流径流将溶解氧丰富的淡水输送到河口，潮汐的涨落则促进了水体在河口内外的交换，使得溶解氧能够在较大范围内分布。生物地球化学过程中，有机物的降解和生物呼吸作用是消耗溶解氧的主要过程。在水体中，有机物不断被微生物分解，这个过程需要消耗氧气，从而导致溶解氧浓度降低。浮游植物的呼吸作用也会消耗一定量的溶解氧。水体层化现象在一定程度上阻碍了上下层水体之间的物质交换，使得底层水体难以获得表层水体中的溶解氧补充，加剧了底层水体的缺氧程度。4.1.2溶解氧收支计算为了深入了解冬季珠江口上游水体溶解氧的动态平衡，对其收支情况进行了详细计算，全面分析了溶解氧的补充和消耗过程。在溶解氧的补充过程中，大气复氧是最为主要的途径。冬季，大气中的氧气通过水气界面不断溶解到水体中，为水体补充溶解氧。根据模拟计算结果，大气复氧对溶解氧的补充量占总补充量的70%左右。这主要是因为冬季水温较低，水体的饱和溶解氧浓度相对较高，使得大气中的氧气更容易溶解到水体中。同时，冬季风速相对较大，增强了水气界面的紊动，促进了氧气的溶解过程。在一些开阔水域，风速较大时，大气复氧的速率明显提高，能够更有效地为水体补充溶解氧。平流-扩散输送也是溶解氧补充的重要过程之一。河流径流和潮汐水流将溶解氧从高浓度区域输送到低浓度区域，从而实现溶解氧的空间分布调整。在珠江口上游，西江、北江和东江的径流携带大量溶解氧丰富的淡水向下游输送，在输送过程中，通过平流作用将溶解氧扩散到周围水体中。潮汐水流的周期性涨落也使得溶解氧在河口地区进行往复输送，进一步促进了溶解氧的分布。平流-扩散输送对溶解氧的补充量占总补充量的20%左右。在径流较大的区域，平流输送对溶解氧的补充作用更为显著，能够有效地提高下游区域的溶解氧浓度。初级生产过程对溶解氧的补充相对较小。冬季，由于光照时间较短、水温较低，浮游植物的生长和光合作用受到一定限制，导致初级生产过程产生的溶解氧相对较少。浮游植物通过光合作用将光能转化为化学能，同时释放出氧气。在冬季，光照强度和温度的不足使得浮游植物的光合作用效率降低，产生的溶解氧也相应减少。初级生产过程对溶解氧的补充量仅占总补充量的10%左右。在一些水温较低、光照不足的区域，初级生产过程几乎无法为水体补充溶解氧。在溶解氧的消耗过程中，有机物降解是主要的消耗途径。水体中的有机物在微生物的作用下发生降解，这个过程需要消耗大量的溶解氧。根据模拟计算，有机物降解对溶解氧的消耗量占总消耗量的60%左右。在珠江口上游，工业废水、生活污水以及农业面源污染等输入了大量的有机物，这些有机物在水体中不断积累并被微生物分解，导致溶解氧大量消耗。在一些排污口附近，有机物浓度较高，有机物降解对溶解氧的消耗作用尤为明显，使得该区域的溶解氧浓度急剧下降。硝化作用也是溶解氧消耗的重要过程。氨氮在硝化细菌的作用下被氧化为硝酸盐，这个过程需要消耗氧气。冬季，虽然水温较低，硝化细菌的活性相对较弱，但硝化作用仍然在进行。硝化作用对溶解氧的消耗量占总消耗量的30%左右。在一些氨氮含量较高的区域，硝化作用对溶解氧的消耗不可忽视，会进一步加剧水体的缺氧程度。其他消耗过程，如浮游植物的呼吸作用、底泥耗氧等，对溶解氧的消耗量相对较小，占总消耗量的10%左右。浮游植物在呼吸过程中会消耗一定量的溶解氧，底泥中的微生物活动也会消耗溶解氧。但与有机物降解和硝化作用相比，这些过程的消耗相对较少。4.1.3控制因素分析冬季珠江口上游水体缺氧受到多种因素的综合控制，水温、径流量、盐度等因素在其中发挥着关键作用，深刻影响着水体中溶解氧的含量和分布。水温是影响冬季水体缺氧的重要因素之一。水温的变化直接影响着溶解氧的溶解度和微生物的活性。在冬季，水温较低，溶解氧在水中的溶解度相对较高，这使得大气中的氧气更容易溶解到水体中，有利于增加水体中的溶解氧含量。低温也会抑制微生物的活性，减缓有机物的降解和硝化作用等耗氧过程。微生物的代谢活动需要适宜的温度条件，当水温较低时，微生物的酶活性降低，代谢速率减慢，从而减少了对溶解氧的消耗。在一些水温较低的区域，有机物降解和硝化作用的速率明显降低，溶解氧的消耗也相应减少，水体缺氧程度相对较轻。水温还会影响水体的层化现象，进而影响溶解氧的分布。冬季水温较低，水体层化相对较弱，有利于上下层水体之间的混合，使得底层水体能够获得更多的溶解氧补充，减少了底层水体缺氧的可能性。径流量对冬季水体缺氧也有着重要影响。河流径流量的大小直接决定了携带溶解氧的能力和对水体的稀释作用。在冬季，珠江口上游的径流量相对较小，这使得河流携带的溶解氧相对较少，对水体的稀释作用也减弱。较小的径流量导致水体的流速减慢，水体与大气的接触面积减小，氧气的溶解速率降低，从而不利于水体中溶解氧的补充。径流量较小还会使水体中污染物的浓度相对升高，加剧了有机物降解和硝化作用等耗氧过程，进一步导致水体缺氧。在一些径流量较小的支流，由于污染物的积累和耗氧过程的增强，水体缺氧现象较为严重。相反，当径流量增大时，河流能够携带更多的溶解氧进入河口，同时对水体中的污染物进行稀释，降低了耗氧物质的浓度，有利于缓解水体缺氧状况。盐度对冬季水体缺氧同样具有重要影响。盐度的变化会影响水体的密度和层化现象，进而影响溶解氧的分布和交换。在珠江口上游，冬季盐度相对较低，这使得水体的密度较小，层化现象相对较弱。较弱的层化有利于上下层水体之间的混合，促进了溶解氧在水体中的均匀分布，使得底层水体能够更容易获得溶解氧补充。盐度还会影响硝化细菌的活性。硝化细菌在适宜的盐度范围内活性较高，当盐度发生变化时，硝化细菌的活性会受到影响，从而影响硝化作用对溶解氧的消耗。在盐度较低的区域，硝化细菌的活性可能会受到一定抑制，减少了硝化作用对溶解氧的消耗，有利于维持水体中的溶解氧含量。4.2夏季水体缺氧情况模拟与分析4.2.1溶解氧分布特征利用构建的物理-生物地球化学耦合模式，对夏季珠江口上游水体溶解氧分布进行模拟，得到的结果呈现出与冬季截然不同的特征。在夏季，珠江口上游水体溶解氧浓度整体低于冬季，这主要是由于夏季水温升高，导致溶解氧的溶解度降低。根据亨利定律，气体在液体中的溶解度与温度成反比，夏季较高的水温使得溶解氧在水中的溶解度显著下降，从而降低了水体的溶解氧含量。微生物活动在夏季更为活跃，有机物的降解和生物呼吸作用增强，进一步消耗了大量的溶解氧。在夏季，浮游植物的生长繁殖速度加快，其呼吸作用消耗的溶解氧也相应增加。在空间分布上，河口附近区域由于受到河流径流和潮汐的共同影响，溶解氧浓度相对较高，但与冬季相比，增加幅度并不明显。这是因为夏季河流径流量虽然较大，但水温较高，溶解氧溶解度降低，且微生物活动旺盛，耗氧增加，抵消了部分因径流量增加带来的溶解氧补充。在远离河口的区域，溶解氧浓度迅速降低，尤其是在底层水体，缺氧现象较为普遍。在珠江口外海区域，底层溶解氧浓度可低至2mg/L以下，严重影响了水生生物的生存和生态系统的平衡。这主要是由于夏季水体层化现象加剧，表层温暖的淡水与底层较冷的海水形成明显的密度差，阻碍了上下层水体之间的物质交换，使得底层水体难以获得表层水体中的溶解氧补充。底层有机物的分解和生物呼吸作用持续消耗溶解氧，进一步加剧了底层水体的缺氧程度。在垂直方向上，夏季珠江口上游水体溶解氧分布呈现出明显的分层现象。表层水体由于与大气直接接触，能够充分进行氧气交换，溶解氧浓度相对较高。在大多数区域，表层溶解氧浓度可达4mg/L以上。然而，随着深度的增加，溶解氧浓度急剧下降。在水体中层，溶解氧浓度迅速降低，形成一个明显的氧跃层。在氧跃层以下，底层水体的溶解氧浓度极低，部分区域甚至接近零。这种垂直方向上的溶解氧分布差异主要是由水体层化和生物地球化学过程共同作用导致的。水体层化使得上下层水体交换困难，限制了溶解氧向底层水体的输送。底层水体中的有机物降解和生物呼吸作用消耗大量溶解氧，而又无法及时得到补充，从而导致底层水体缺氧。4.2.2溶解氧收支计算为深入剖析夏季珠江口上游水体溶解氧的动态平衡，对其收支情况展开详细计算，全面分析各收支过程对溶解氧的影响。在溶解氧的补充过程中，大气复氧依然是主要途径，但相较于冬季，其占总补充量的比例有所下降，约为65%。这是因为夏季水温较高，溶解氧在水中的溶解度降低，大气中的氧气溶解到水体中的难度增加。夏季风速相对较小，水气界面的紊动减弱，也不利于氧气的溶解。在一些风速较小的区域，大气复氧的速率明显降低，导致其对溶解氧的补充量减少。平流-扩散输送对溶解氧的补充量占总补充量的25%左右，较冬季有所增加。夏季河流径流量增大，水流速度加快，平流作用增强，能够将更多的溶解氧从高浓度区域输送到低浓度区域。在西江、北江和东江的下游区域，由于径流量较大，平流输送对溶解氧的补充作用显著，使得这些区域的溶解氧浓度相对较高。夏季水体的紊动增强，扩散作用也相应增强，进一步促进了溶解氧的分布。在河口地区，潮汐的涨落加剧了水体的紊动，使得溶解氧能够更均匀地分布在水体中。初级生产过程对溶解氧的补充量占总补充量的10%左右，与冬季相比变化不大。虽然夏季光照时间长、强度高，有利于浮游植物进行光合作用，但由于水体中营养盐含量相对较低，限制了浮游植物的生长和繁殖，导致初级生产过程产生的溶解氧相对有限。在一些营养盐缺乏的区域，浮游植物的光合作用受到抑制，初级生产过程几乎无法为水体补充溶解氧。在溶解氧的消耗过程中，有机物降解和硝化作用仍然是主要的消耗途径。有机物降解对溶解氧的消耗量占总消耗量的55%左右，夏季水温升高，微生物活性增强，有机物的降解速率加快，导致溶解氧的消耗增加。在一些排污口附近，有机物浓度较高，有机物降解对溶解氧的消耗作用尤为明显，使得该区域的溶解氧浓度急剧下降。硝化作用对溶解氧的消耗量占总消耗量的35%左右，夏季氨氮含量相对较高，硝化细菌的活性也较强，使得硝化作用对溶解氧的消耗增加。在一些氨氮污染严重的区域，硝化作用对溶解氧的消耗不可忽视，进一步加剧了水体的缺氧程度。其他消耗过程，如浮游植物的呼吸作用、底泥耗氧等，对溶解氧的消耗量占总消耗量的10%左右，与冬季相比变化不大。浮游植物在呼吸过程中会消耗一定量的溶解氧，底泥中的微生物活动也会消耗溶解氧。但与有机物降解和硝化作用相比，这些过程的消耗相对较少。4.2.3控制因素分析夏季珠江口上游水体缺氧受到多种因素的综合影响，温度、盐度、径流量等因素在其中扮演着关键角色，深刻影响着水体中溶解氧的含量和分布。温度是影响夏季水体缺氧的重要因素之一。夏季水温升高，一方面降低了溶解氧的溶解度，使得水体中能够容纳的溶解氧减少。根据相关研究，水温每升高10℃，溶解氧的溶解度约降低15%。另一方面，温度升高促进了微生物的生长和代谢活动，使得有机物的降解和硝化作用等耗氧过程加快。微生物的酶活性在适宜的温度范围内会随着温度的升高而增强，从而加速了有机物的分解和氨氮的氧化，消耗更多的溶解氧。在一些水温较高的区域，有机物降解和硝化作用的速率明显加快，导致溶解氧大量消耗，水体缺氧程度加剧。温度还会影响水体的层化现象，夏季水温较高，水体层化明显，阻碍了上下层水体之间的物质交换，使得底层水体难以获得表层水体中的溶解氧补充，进一步加剧了底层水体的缺氧状况。盐度对夏季水体缺氧也有着重要影响。在夏季，珠江口上游水体盐度相对较低，这使得水体的密度较小，层化现象相对较弱。然而，由于河流径流量较大，携带的营养物质较多，在河口地区容易形成盐度锋面。盐度锋面的存在会影响水体的混合和物质输运，导致溶解氧的分布不均匀。在盐度锋面附近，水体的混合受到阻碍，溶解氧难以在不同盐度的水体之间进行交换，使得部分区域的溶解氧含量较低。盐度还会影响硝化细菌的活性。硝化细菌在适宜的盐度范围内活性较高，当盐度发生变化时，硝化细菌的活性会受到影响，从而影响硝化作用对溶解氧的消耗。在盐度较低的区域，硝化细菌的活性可能会受到一定抑制，减少了硝化作用对溶解氧的消耗；而在盐度较高的区域，硝化细菌的活性可能会增强，增加了硝化作用对溶解氧的消耗。径流量对夏季水体缺氧同样具有重要影响。夏季珠江口上游径流量较大，一方面，河流携带的溶解氧增加，通过平流-扩散输送作用，能够为水体补充一定量的溶解氧。在径流量较大的河流下游区域，溶解氧浓度相对较高。另一方面，径流量增大也会导致水体中污染物的稀释作用增强，降低了有机物和氨氮等耗氧物质的浓度，从而减少了溶解氧的消耗。然而，如果径流量过大，可能会导致水体的停留时间缩短，使得水体中的生物地球化学过程来不及充分进行，影响溶解氧的收支平衡。径流量还会影响水体的层化现象。较大的径流量可能会增强水体的紊动，破坏水体的层化，有利于上下层水体之间的物质交换，增加底层水体的溶解氧补充。但在一些情况下，径流量的变化可能会导致水体层化的不稳定，反而加剧了底层水体的缺氧状况。4.3季节差异对比分析对比冬夏季珠江口上游水体溶解氧收支情况，发现存在显著差异。在溶解氧补充方面，大气复氧均为主要补充过程，但夏季由于水温升高导致溶解氧溶解度降低，大气复氧占总补充量的比例从冬季的70%降至65%。平流-扩散输送对溶解氧的补充量在夏季占总补充量的25%左右，较冬季的20%有所增加，这主要是因为夏季河流径流量增大，水流速度加快，平流作用增强，且水体紊动加剧，扩散作用也相应增强。初级生产过程对溶解氧的补充量在冬夏季变化不大，均占总补充量的10%左右。在溶解氧消耗方面，有机物降解和硝化作用均是主要的消耗过程。冬季有机物降解对溶解氧的消耗量占总消耗量的60%左右，硝化作用占30%左右；夏季有机物降解占总消耗量的55%左右，硝化作用占35%左右。夏季有机物降解和硝化作用对溶解氧的消耗均高于冬季，这是由于夏季水温升高，微生物活性增强，有机物的降解速率加快，氨氮含量相对较高，硝化细菌的活性也较强。从缺氧范围来看，冬季珠江口上游水体缺氧范围相对较小，主要集中在一些局部区域，如河口附近的底层水体以及水流缓慢、污染物聚集的区域。而夏季缺氧范围明显扩大，在远离河口的区域以及底层水体普遍存在缺氧现象，尤其是在珠江口外海区域，底层溶解氧浓度可低至2mg/L以下。温度是导致冬夏季溶解氧收支和缺氧范围差异的重要因素之一。夏季水温升高，一方面降低了溶解氧的溶解度，使得水体中能够容纳的溶解氧减少；另一方面，促进了微生物的生长和代谢活动，使得有机物的降解和硝化作用等耗氧过程加快。根据相关研究，水温每升高10℃，溶解氧的溶解度约降低15%。微生物的酶活性在适宜的温度范围内会随着温度的升高而增强，从而加速了有机物的分解和氨氮的氧化，消耗更多的溶解氧。温度还会影响水体的层化现象，夏季水温较高，水体层化明显，阻碍了上下层水体之间的物质交换，使得底层水体难以获得表层水体中的溶解氧补充，进一步加剧了底层水体的缺氧状况。盐度也对冬夏季差异产生影响。在冬季，盐度相对较低，水体的密度较小，层化现象相对较弱，有利于上下层水体之间的混合，促进了溶解氧在水体中的均匀分布。盐度还会影响硝化细菌的活性，在适宜的盐度范围内，硝化细菌活性较高。夏季，虽然珠江口上游水体盐度相对较低，但由于河流径流量较大，携带的营养物质较多，在河口地区容易形成盐度锋面。盐度锋面的存在会影响水体的混合和物质输运，导致溶解氧的分布不均匀。盐度的变化还会影响硝化细菌的活性，进而影响硝化作用对溶解氧的消耗。径流量同样是导致季节差异的关键因素。冬季径流量相对较小，河流携带的溶解氧较少，对水体的稀释作用也减弱，不利于水体中溶解氧的补充。较小的径流量还会使水体中污染物的浓度相对升高，加剧了有机物降解和硝化作用等耗氧过程，进一步导致水体缺氧。而夏季径流量增大，一方面河流携带的溶解氧增加，通过平流-扩散输送作用，能够为水体补充一定量的溶解氧；另一方面，径流量增大也会导致水体中污染物的稀释作用增强，降低了有机物和氨氮等耗氧物质的浓度，从而减少了溶解氧的消耗。径流量的变化还会影响水体的层化现象，较大的径流量可能会增强水体的紊动，破坏水体的层化，有利于上下层水体之间的物质交换，增加底层水体的溶解氧补充。五、数值实验与敏感性分析5.1实验设计为深入探究珠江口上游水体缺氧对不同环境条件变化的响应，设计了一系列数值实验。在实验中，分别设置不同的情景，系统研究了增大径流量、增大风速、降低东江污染输入以及降低广州污染输入等因素对水体缺氧的影响。针对增大径流量实验，在当前模拟的基础上，将珠江口上游主要河流（西江、北江、东江）的径流量分别增大20%。通过改变河流的流量参数，模拟河流携带更多的溶解氧和营养物质进入珠江口上游水体，观察水体中溶解氧分布和缺氧范围的变化。这有助于了解河流径流量变化对水体缺氧的直接影响，以及在不同径流量条件下，水动力过程和生物地球化学过程如何相互作用，进而影响水体的溶解氧平衡。在增大风速实验中，将该区域的平均风速增大3m/s。风速的变化会直接影响水气界面的紊动程度，进而改变氧气的溶解速率和水体的混合情况。通过调整风速参数，模拟大气复氧过程的增强，分析风速变化对水体溶解氧补充的影响，以及对整个水体缺氧状况的改善作用。这对于评估风力在珠江口上游水体溶解氧平衡中的作用具有重要意义。降低东江污染输入实验旨在研究东江污染对珠江口上游水体缺氧的影响。将东江输入到珠江口上游水体的污染物（包括有机物、氨氮等）浓度降低50%。通过改变污染物输入浓度，观察水体中溶解氧浓度的变化以及缺氧范围的缩小情况。这有助于明确东江污染在水体缺氧形成过程中的贡献，为制定针对性的污染控制措施提供依据。降低广州污染输入实验则是将广州地区输入到珠江口上游水体的污染物浓度降低50%。广州作为珠江口上游的重要城市，其污染物排放对水体质量有着重要影响。通过模拟降低广州污染输入，分析水体中溶解氧含量的增加以及缺氧现象的缓解程度，为广州地区的污染治理提供科学参考。5.2结果分析增大径流量实验结果显示，当珠江口上游主要河流径流量增大20%时，水体溶解氧浓度显著增加。在河流下游区域，溶解氧浓度平均升高了1mg/L左右。这是因为径流量的增大使得河流携带更多的溶解氧进入水体，通过平流-扩散输送作用，将溶解氧输送到更广泛的区域。径流量的增加还增强了水体的紊动，促进了溶解氧在水体中的混合和扩散，使得溶解氧分布更加均匀。在一些原本缺氧的区域，溶解氧浓度得到明显提升，缺氧范围缩小了约30%。这表明增大径流量能够有效改善珠江口上游水体的溶解氧状况，缓解水体缺氧问题。增大风速实验表明，当平均风速增大3m/s时，水体溶解氧浓度也有所增加。在整个研究区域，溶解氧浓度平均升高了0.5mg/L左右。风速的增大增强了水气界面的紊动，使得大气复氧过程加快，更多的氧气从大气中溶解到水体中。在一些开阔水域，风速增大后，大气复氧速率明显提高，溶解氧浓度增加更为显著。风速的变化还影响了水体的混合情况，促进了上下层水体之间的交换，使得底层水体能够获得更多的溶解氧补充。这说明增大风速对改善珠江口上游水体缺氧状况具有一定的积极作用。降低东江污染输入实验中，将东江输入到珠江口上游水体的污染物浓度降低50%后，其输入之后河段（狮子洋一虎门）的水体溶解氧浓度显著增加。在狮子洋一虎门河段，溶解氧浓度平均升高了1.5mg/L左右。这是因为减少东江的污染输入，降低了水体中有机物和氨氮等耗氧物质的含量，从而减少了有机物降解和硝化作用等耗氧过程对溶解氧的消耗。有机污染物输入的变化对溶解氧的影响更大，当有机污染物浓度降低时，有机物降解过程对溶解氧的消耗大幅减少，使得水体中的溶解氧得以积累。而广州一黄埔河段的溶解氧浓度基本不受影响，这表明东江污染主要影响其输入之后的河段，对上游河段的影响较小。降低广州污染输入实验中，将广州地区输入到珠江口上游水体的污染物浓度降低50%后，整个研究区域的溶解氧浓度显著增加。在广州一虎门河段，溶解氧浓度平均升高了2mg/L左右。这是因为广州作为重要的污染源，其污染物排放对整个研究区域的水质影响较大。降低广州污染输入，减少了水体中耗氧物质的总量，使得有机物降解和硝化作用等耗氧过程对溶解氧的消耗减少，从而提高了水体的溶解氧浓度。有机污染物对水体缺氧贡献更大，当有机污染物浓度降低时，水体中的溶解氧浓度提升更为明显。这说明降低广州污染输入对改善珠江口上游水体缺氧状况具有重要作用，能够有效提高整个研究区域的水体溶解氧含量。5.3污染控制建议基于数值实验结果，为有效改善珠江口上游水体缺氧状况，提出以下污染控制和管理建议：严格控制污染物排放：对工业废水和生活污水进行严格监管，确保其达标排放。建立完善的污水处理设施，提高污水处理能力，减少有机物和氨氮