环境科学  2023, Vol. 44 Issue (1): 210-218   PDF    
长三角典型城郊流域生物可降解性有机质的分布及影响因素
吴文雨1,2, 马菁晟3, 杨磊4, 李敏4, 唐剑锋1     
1. 中国科学院城市环境研究所城市环境与健康重点实验室, 厦门 361021;
2. 上海应用技术大学生态技术与工程学院, 上海 201418;
3. 江苏省科技镇长团句容团, 句容 212400;
4. 中国科学院生态环境研究中心城市与区域生态国家重点实验室, 北京 100085
摘要: 快速城镇化改变了流域水体中碳的生物地球化学循环过程, 生物可利用性的溶解性有机质组分是其中最为关键的一环, 辨识生物可降解性有机质(BDOM)的时空分布特征及其影响因素对流域水质管控具有重要意义和应用价值.以长三角地区典型城郊流域樟溪为研究对象, 根据流域地形特征、土地利用及人类活动强度布设监测点位, 于2019年分别在雨季和旱季采样, 利用三维荧光平行因子(EEM-PARAFAC)方法结合源汇景观模型研究流域水体BDOM的时空分布特征.结果表明, 流域中生物可降解性有机碳浓度范围在0.57~6.80 mg ·L-1, 且具有较高的时空异质性, 人类活动强度较高的区域水体中BDOM的浓度也相对较高, 且雨季显著高于旱季.EEM-PARAFAC分析结果表明, 流域BDOM主要包括陆源腐殖质(C1)和类蛋白质类(C2)这2种荧光组分.流域BDOM及其陆源腐殖质荧光组分主要受土地利用和人类活动的影响, 其浓度与农业及城镇用地比例和源汇景观负荷比(LWLI)关系密切, 表明城镇化过程中人类活动是影响BDOM分布的重要因素.
关键词: 城郊流域      溶解性有机质(DOM)      生物可降解性有机质(BDOM)      土地利用格局      平行因子分析     
Distribution of Biodegradable Dissolved Organic Matter and Its Affecting Factors in a Typical Peri-urban Watershed in Yangtze River Delta
WU Wen-yu1,2 , MA Jing-sheng3 , YANG Lei4 , LI Min4 , TANG Jian-feng1     
1. Key Laboratory of Urban Environment and Health, Institute of Urban Environment, Chinese Academy of Sciences, Xiamen 361021, China;
2. School of Ecological Engineering, Shanghai Institute of Technology, Shanghai 201418, China;
3. Jurong Group of Jiangsu Science and Technology Mayor Group, Jurong 212400, China;
4. State Key Laboratory of Urban and Regional Ecology, Research Center for Eco-Environmental Sciences, Chinese Academy of Sciences, Beijing 100085, China
Abstract: Dissolved organic matter (DOM) plays key roles in the carbon biogeochemical cycle, and biodegradable dissolved organic matter (BDOM) is one of the key fractions of DOM. Rapid urbanization and intensive human activities substantially influence the distribution of DOM at the watershed scale. Identifying the spatial and temporal variability in BDOM has become an important and urgent issue of water quality control in rapid urbanization areas. However, limited studies have been conducted to explore the role of human activities on the occurrence and distribution of BDOM in peri-urban watersheds. In this study, the spatial and temporal distribution of BDOM and related affecting factors were investigated in a typical peri-urban watershed (Zhangxi watershed) located at Ningbo City in Yangtze River Delta. Water samples were collected in wet and dry seasons in 2019 based on topographic features, land use, and intensity of human activities. The BDOM were characterized by fluorescence excitation-emission matrix and parallel factor analysis (EEM-PARAFAC), and land use patterns were analyzed using the Source-Sink Landscape Model. The results of this study showed that the BDOM concentrations ranged from 0.57 to 6.80 mg ·L-1. Obvious spatial and temporal heterogeneities of BDOM were found at the watershed scale, and significantly higher concentrations of BDOM were observed in the wet season than those in the dry season. Furthermore, relatively high concentrations of BDOM were found in areas with relatively higher intensive human activities. Two fluorescent components (a terrestrial humic-like substance and protein-like substance) were observed using the PARAFAC model. The results of spatial analysis showed that terrestrial humic-like fluorescent components were closely positively correlated with anthropogenic parameters (percentages of agricultural and urban land and ratio of source and sink landscapes). The results showed that the occurrence and distribution of BDOM were strongly influenced by human activities, which could provide scientific guidance for water quality control and related land management in peri-urban aquatic ecosystems.
Key words: peri-urban watershed      dissolved organic matter (DOM)      biodegradable dissolved organic matter (BDOM)      land use pattern      parallel factor analysis     

城郊是连接城镇与乡村的关键过渡地带, 城郊地区土地利用多样、景观结构复杂、生态系统结构和功能独特[1].城镇化梯度的三元结构(城镇核心区、城镇郊区和乡村/自然地区)的边界及其承载的环境过程也随城镇扩张不断改变[2, 3].一方面, 城郊为城镇提供关键的原材料、能源和农产品等, 使得城郊成为生态系统服务功能冲突最为强烈的.另一方面, 大量的城市工业污染和废弃物进入城郊, 使得城郊地区生态系统服务的供需冲突最为强烈[4, 5].溶解性有机质(dissolved organic matter, DOM)是一类结构复杂的有机混合物且极其不稳定, 是无机碳和生命形态碳的关键纽带[6].生物可降解性有机质(biodegradable dissolved organic matter, BDOM)是指水体中可被细菌降解成CO2或用于维持自身生长繁殖的有机物, 为微生物提供碳源和能源, 是DOM重要的组成部分[7].环境中的细菌群落结构与功能也与DOM的生物降解性密切相关[8, 9].具有不同生物降解性的DOM会引起细菌多样性, 群落结构和代谢的不同反应, 反映在碳消耗、生长和呼吸速率和酶的降解方面[10~12].不稳定的蛋白质和脂肪族化合物能够被细菌快速利用, 而稳定的材料只能被特定的基团降解[10~12].以生物可降解溶解性有机碳(biodegradable dissolved organic carbon, BDOC)可作为饮用水生物稳定性指标来衡量BDOM, 其代表异养微生物降解的溶解性有机碳[13].BDOC已广泛应用于地表水有机物[14]、天然有机物[15]和海水[16]的生物降解性和转化研究, 但目前对于实际水体中BDOM的特点和组成分析报道较少, 探索生物可降解性的过程, 可以更好地理解DOM动力学以及生态系统中能量和物质循环[17].

人类活动对城郊生态系统DOM的含量和性质造成了剧烈影响, 打破了原有的碳循环过程[18].有研究表明, 植被类型和土地利用类型可对流域内DOM和BDOM造成重要影响[7].Zhou等[19]的研究发现, 人为的干扰、经济发展水平的提高以及城市和农业土地使用加剧, 会增加流域中有机碳的负荷或改变其组成和生物可利用性.近年来, 随着我国快速的城镇化过程, 城郊地区已成为生态系统服务功能冲突最激烈的区域[5].不同的人类活动会造成独特的DOM组成, 城郊剧烈的人类活动直接改变了DOM的含量和性质, Zhang等[20]对嘉陵江南充段的研究发现, 其有色可溶性有机物(colored dissolved organic matter, CDOM)的来源和组成主要受城市化发展的影响; 对于东北吉林省水域的研究表明, 城镇化发展过程中的人为活动促进了DOM数量的增加, 导致DOM组分向类蛋白质成分转化[21].源汇景观格局分析方法将景观格局与生态过程有机结合, 可定量表征不同景观类型及其空间配置对特定生态过程的影响[22], 已广泛应用于生态环境领域.例如, 周添惠等[23]利用源汇景观格局分析方法研究了黄土高原渭河源流域径流及泥沙特征, 结果表明在长时间尺度上源汇景观格局演变改变了流域径流及泥沙特征.在海河流域的研究中也表明, 源汇景观格局指数与总氮的空间变异特征密切关联[24].李敏等[25]在樟溪流域中的研究也显示, 源汇景观负荷比指数(LWLI)与水体抗生素浓度呈正相关关系, 流域LWLI越大其水体抗生素浓度越高.

然而现有的研究多致力于流域DOM的分布及影响因素, 对城郊流域BDOM的分布特性及其与源汇景观格局的关系的研究还很薄弱.因此, 本文以宁波城郊关键带樟溪流域为研究区域, 通过对流域水体的采样分析, 分析不同城镇化梯度下流域水体BDOM的分布特性, 揭示城郊流域BDOM浓度变化的驱动因子, 以期为快速城镇化背景下城郊流域的水环境质量和安全保障提供科学支撑.

1 材料与方法 1.1 研究区概况

樟溪流域(29°44′~29°51′N, 121°13′~121°21′E) 位于浙江省宁波市(图 1), 总面积89.7 km2.该流域属亚热带季风气候, 有明显的雨季和旱季, 雨季一般在4~9月, 累计降雨量为1 035 mm, 旱季一般在10~12月和1~3月, 累计降水量为782 mm[4].研究区内樟溪为奉化江支流, 全长14.5 km, 上游为饮用水源地皎口水库, 流经章水、龙观和鄞江这3个乡镇.研究区位于四明山区向宁绍平原的过渡地带, 主要土地利用类型有城镇建设用地、农田、园地、林地和水库等, 其土地利用格局受城镇化影响较大.农田主要以贝母、花生和蔬菜种植为主, 林地主要为亚热带常绿阔叶林, 园地主要为经济林木苗圃[26].

图 1 研究区土地利用和采样点分布示意 Fig. 1 Land use and water sampling sites of study area

基于流域的土地利用特征、水文路径和地形因素, 根据人为干扰程度, 将研究区划分为4个城市化梯度明显的子流域:子流域4>子流域1>子流域2>子流域3, 设置了16个监测断面用于地表水的长期定位监测, 其中采样点W1~W9属于子流域1, W11~W12属于子流域2, W14属于子流域3, 流域点位靠近水库, 为森林覆盖率较高的自然生态系统(96.3%), W10、W13、W15和W16属于子流域4, 相比其它子流域其农业用地和建设用地比例较高, 为耕地(35.8%)和城镇(10.3%)[5].

1.2 样品采集

样品采集于2019年6月(雨季)和12月(旱季).使用pH < 2的硫酸处理过的5 L棕色玻璃瓶, 在每个河流采样点中间水面30 cm深处收集3个平行水样(样瓶用水样预冲洗3次), 在密封无气泡状态下4℃低温黑暗保存运回, 并在1 d之内完成实验培养[7].

1.3 实验室培养测定BDOC

参照文献[7, 27]进行DOM的培养.首先加入3 mL过1.2 μm GF/F滤膜作为接种微生物到血清瓶中, 然后移取30 mL过0.45 μm GF/F滤膜的水样到血清瓶中进行培养, 再分别加入100 μL的NH4NO3和K2HPO4溶液, 以降低营养盐对微生物活性的限制, 用等体积去离子水加营养盐设置空白对照实验.所有样品放入恒温培养箱25℃条件黑暗培养, 每24 h打开瓶盖并摇晃一次, 以确保有氧条件; 分别在第0、2、7、14、21和28 d取样, 用0.45 μm GF/F滤膜(预先在马弗炉中450℃灼烧5 h除去有机质)过滤, 并立即测定DOC的浓度, 最终结果取3个平行水样的平均值.由于0.45 μm滤膜重新过滤几乎可以过滤掉全部的微生物量, 所以最后测得DOC降解程度包括了微生物矿化作用和同化作用两部分的量, 以上所有操作均在无菌超净台完成.DOC浓度采用总有机碳(TOC V-CPH, Elementar, Germany)进行.

1.4 DOM的光谱特性

DOM的光谱表征包括紫外-可见光光谱和荧光光谱, 采用紫外分光光度计(UV3600, Shimadzu, Japan)扫描200~600 nm吸光度, 采用荧光分光光度计(F-4600, Hitachi, Japan)进行水体DOM的三维荧光光谱测定, 配以1 cm石英比色皿.以超纯水为实验空白, 对水样进行三维荧光扫描.仪器光源为150 W氙灯, 光电倍增电压为800 V, 激发和发射狭缝宽度均设置为5 nm, 相应时间0.1 s, 激发波长(Ex)为220~450 nm, 扫描间隔5 nm; 发射波长(Em)为260~600 nm, 扫描间隔1 nm; 扫描速度为2 400 nm ·min-1.所有水样的三维荧光均需减去空白光谱, 以消除拉曼散射的影响[28, 29]; 并将发射波长等于激发波长和发射波长等于2倍激发波长包围外的区域光谱值赋值为零, 以修正瑞利散射的影响[30].使用MATLAB 2014a(Math Works, Natick, USA)和DOMFluor工具箱执行平行因子分析(parallel factor analysis, PARAFAC)[31, 32], 使用OpenFluor网站(http://www.openfluor.org)的在线光谱库来验证PARAFAC的结果[4].

1.5 源汇景观模型

源汇景观模型(source-sink landscape model, SSLM)是中国科学院生态环境研究中心城市与区域生态国家重点实验室开发的空间分布式评价模型[26], 是基于ArcGIS开发的toolbox模型插件.SSLM根据研究区的土地利用特征和DOM的排放性质, 将建设用地、农田和林地划分为“源”景观, 将河流和水库划分为“汇”景观.通过SSLM可以快捷地进行流域范围划分, 提取汇水区范围内的DEM和土地利用, 输出距流域出水口距离、流域坡度并通过海拔、坡度和距离等参数绘制不同源和汇景观的洛伦兹曲线, 并计算出源汇景观负荷比指数(LWLI).对于海拔高程和距离要素, LWLI越大, 景观在高程或距离要素上的分布越有助于推动子流域出口点的环境效应, 而对于坡度要素, LWLI越大, 景观在坡度要素上的分布越有助于阻碍子流域出口点的环境效应[25].

1.6 数据处理

本文对实验室培养测定BDOM用到的两种关系见式(1)和式(2)[7]

(1)
(2)

式中, BDOCtt时刻的微生物可降解溶解性有机碳浓度(mg ·L-1), DOCt0为样品初始时刻溶解性有机碳浓度(mg ·L-1), DOCtt时刻的溶解性有机碳浓度(mg ·L-1), BDOCA为BDOC的降解率(%).本文研究利用SPSS 25.0和EXCEL 2017对数据进行处理和描述性统计, 采用Origin 2017处理土地利用与DOC、BDOC的相关关系; 使用ArcGIS 10.5对研究区土地利用数据进行处理, 源汇景观负荷比指数由SSLM进行计算; 使用R语言Corrplot对DOC和BDOC和源汇景观负荷比指数以及土地利用进行矩阵相关分析.

2 结果与分析 2.1 BDOM的生物利用时间动力学

DOM浓度经常用溶解性有机碳(dissolved organic carbon, DOC) 浓度来进行量化, 所以利用DOC浓度变化情况来解释BDOM的降解动力学过程[33].樟溪流域水样的DOM降解动力学过程如图 2所示, 在培养过程中, 水样的DOC浓度逐渐下降直至达到平衡.整个降解动力学过程可用一级降解动力学方程拟合.

图 2 典型样本DOC浓度的变化和生物降解动力学模型 Fig. 2 Variation in DOC concentration and biodegradation kinetic model for a typical sample

基于一级降解动力学理论, 拟合典型土地利用样本数据随时间的变化, 采用生物降解动力学模型, 表述如式(3)[34]:

(3)

式中, t为降解时间(d), CDOC, tt时刻DOC生物降解率(%), CDOC, eq为达到平衡时DOC生物降解所占的质量分数(%), k为降解系数(d-1).对不同典型土地利用样本DOC生物降解所占的质量分数随时间变化进行曲线拟合, 得到降解系数k和达到平衡时CDOC, eq的浓度, 动力学模型结果显示所有样品在28 d之内均已达到平衡[19], 因此本研究以最终培养28 d降解后损失的DOM为BDOM.

2.2 樟溪流域BDOC的区域分布特征

结果表明, 培养完成后樟溪流域的ρ(DOC)在0.87~3.29 mg ·L-1之间, 其中雨季ρ(DOC)为0.98~3.29 mg ·L-1, 平均值为(1.85±0.66) mg ·L-1, 旱季ρ(DOC)为0.87~2.13 mg ·L-1, 平均值为(1.51±0.36) mg ·L-1, 雨季浓度显著高于旱季(P < 0.001). ρ(BDOC)在0.57~6.80 mg ·L-1之间, 雨季ρ(BDOC)范围为1.12~6.80 mg ·L-1, 平均值为(3.69±1.70) mg ·L-1, 在旱季ρ(BDOC)范围为0.57~2.17mg ·L-1, 平均值为(1.41±0.49) mg ·L-1, BDOC浓度在旱、雨季存在极显著性差异(P < 0.001).在各子流域之间BDOC浓度也存在极显著的差异(P=0.001), 旱雨季中BDOC浓度在各子流域中的分布趋势保持一致, 中位值由高到低依次为:子流域4>子流域1>子流域2>子流域3[图 3(a)].

图 3 旱雨季BDOC浓度和BDOC降解率的空间分布 Fig. 3 Spatial distribution of BDOC concentrations and BDOC degradation rate in dry and rainy seasons

樟溪水样的BDOCA年增长率平均值为(56.2±12.4)%, 其中雨季样品的BDOCA增长率平均值为(65.1±6.6)%, 旱季增长相对较低[(47.3±10.3)%], 雨季和旱季存在极显著差异(P < 0.001), 这是因为水温对水样的影响, 旱季多低温(< 5℃)会降低水样中的微生物活性, 减缓DOC的降解, 而雨季温度相对较高(>20℃)促进了水样中的细菌生长与活动, 更易造成DOC的损失与降解[35]; 各子流域之间BDOCA值也存在极显著差异[图 3(b)], 其中子流域4的BDOCA值最高, 雨季和旱季分别为(65.0±9.3)%和(47.3±10.3)%; 子流域3则最低, 雨季和旱季分别为61.5%和35.2%.同时, 各流域BDOCA值与BDOC浓度在子流域中的分布趋势一致, 随着人类活动程度增加而升高[36, 37], 表明溶解有机质的可利用性越强.

2.3 DOM的荧光特性

利用三维荧光-平行因子模型分析得到2个组分:陆源腐殖质组分C1(Ex=230, Em=425 nm)[38, 39], 蛋白质类组分C2(Ex=225(275), Em=332 nm)[40].本研究结果与EEM-PARAFAC模型数据库OpenFluor有很好的匹配(相似性>95%).

培养前后樟溪水样的陆源腐殖质C1无显著差异(P>0.05), 且C1组分的Fmax值(0.45±0.18)略低于蛋白质类C2组分(0.49±0.15); C2组分反而平均降低了0.17±0.2, 各样点之间具有显著差异性(P < 0.05).4个子流域中组分C1的荧光强度与BDOCA分布一致:子流域4>子流域1>子流域2>子流域3; 组分C2的荧光强度分布为:子流域3>子流域4>子流域1>子流域2 [图 4(a)]; 这可能是因为夏季垂钓事件频繁发生, 水库人为影响因素较大, 所以蛋白质类C2组分相对较高[41]; 而以森林为主的子流域2荧光强度最低, Tang等[4]的研究结果表明, 森林等自然流域所含的荧光组分相较于城区等强度较低.培养后各子流域陆源腐殖质C1和蛋白质类C2分布与培养前C1组分一致[图 4(b)].而组分C2的降解幅度显著高于C1, 相较于其它子流域, 子流域4的组分C2降解量最高.该结果与樟溪流域整体变化趋势一致, 表明在培养过程中, 腐殖质类物质含量相对较低且相对稳定[42].而微生物活动也会优先与类蛋白质发生反应[19, 43], 同时C1和C2组分在不同子流域内的分布情况表示人为活动和土地利用比例会影响河流DOM的化学组成和浓度变化[44].

P < 0.05和P < 0.01表示显著性 图 4 PARAFAC降解前与降解后荧光组分箱形图比较 Fig. 4 Comparison of box diagrams of fluorescent fractions before and after PARAFAC degradation

3 讨论

源汇景观模型通过输入流域数字高程模型(digital elevation model, DEM)数据、土地利用数据和流域出口位置计算得出基于距离的景观负荷比指数(LWLI.D), 基于坡度的景观负荷比指数(LWLI.S), 源汇景观负荷比指数(LWLI), 基于海拔的景观负荷比指数(LWLI.E)等指标. 然后进一步通过源汇模型指数, 将源汇土地利用比例和景观负荷比指数与其流域水文过程结合起来, 可描述BDOM在流域内的时空分布规律.利用源汇景观计算樟溪流域的土地利用格局可间接表征流域内人类活动强度.结果表明樟溪流域培养28 d后的BDOC浓度与源汇模型中“农业+建设用地”占比有显著正相关关系[R2=0.60, P < 0.05, 图 5(a)], 与林地用地占比显著负相关[R2=0.58, P < 0.05, 图 5(b)].在城市化梯度的研究中也发现, 城镇流域中BDOC浓度显著高于森林溪流(P < 0.001)[45]; 这表明随着流域内人类活动强度(农业和城镇用地)的增加, 流域水体中BDOC浓度有明显的上升趋势[46].土地利用方式与人类活动关系密切, 而人为干扰较少的林地中BDOC浓度则较低[47].樟溪流域DOM浓度分布的结果也表明, 流域内未进行充分处理的生活污水等进入水体也可能增加BDOC浓度[5].这些发现说明, 流域城市化导致河流DOM浓度变化影响河流生态系统的新陈代谢, 并最终影响河流中有机碳的命运[45].

(a)和(b)分别表示BDOC浓度与“农业+建设用地”和森林用地占比的相关性分析; (c)和(d)分别表示BDOM荧光强度(BDOMC)与“农业+建设用地”和森林用地占比的相关性分析 图 5 BDOC浓度和BDOM与不同土地利用类型比例的相关性分析 Fig. 5 Correlation analysis of BDOC concentrations and BDOM to different land use percentages

BDOM与土地利用的相关性和BDOC浓度一致, 与源汇模型中“农业+建设用地”占比有显著正相关关系[R2=0.47, P < 0.01, 图 5(c)], 与林地用地占比显著负相关[R2=0.51, P < 0.01, 图 5(d)].这是因为在林地中腐殖质的来源主要为陆源植物的凋落; 而在农业和建设用地中, 人为活动会影响水体中的DOM组成, 这在一定程度上会促成微生物的生长[25], 在内源微生物的降解和部分陆源输入的情况下[48], 导致BDOM荧光强度随着“农业+建设用地”占比的增加而上升.

图 5(c)图 5(d)可知, BDOM与“农业+建设用地”占比和BDOC浓度相一致, 也存在着显著正相关(R2=0.47, P < 0.01), 且BDOM与森林用地占比存在显著负相关(R2=0.51, P < 0.01).由此也可表明人类活动强的土地利用类型是水体BDOM和BDOC的重要来源.这是因为与城镇相邻的河流, 其沿岸分布了一定面积的农田, 由人类活动产生的DOM会随地表径流进入水体, 这在一定程度上会造成微生物的生长, 容易发生藻类暴发等, 严重污染水体环境[30].

BDOC和BDOM与源汇景观格局特征的相关矩阵分析结果显示(图 6), 林地和水库用地占比与BDOCA和BDOM均呈负相关, 农业用地占比与BDOC28(雨季)和BDOMC1有较高的正相关(R2=0.68和R2=0.65), 建设用地占比与BDOMC1来源密切相关(R2=0.77).此外, BDOCA、BDOC28和BDOM均与源汇景观负荷比指数呈正相关关系, 表明源汇景观负荷比指数越大的流域, 水体BDOM的浓度越高, 但景观负荷比指数和土地利用比例对BDOM、BDOC雨季和旱季变化的影响程度也各有不同.BDOMC1相比BDOMC2与源汇景观负荷比指数(LWLI)关系更密切, BDOMC1受距离景观负荷比指数(LWLI.D)和坡度景观负荷比指数(LWLI.S)影响较大, 而BDOMC2受源汇景观负荷比指数影响较小, 这可能与蛋白质类组分C2排放主要来自于下游人口密集活动相关, 从上游到下游“源”景观(城镇和农田)面积逐渐增大, 而“汇”(林地和水库)景观面积逐渐减小, 在地理空间上, 也影响了BDOM在上下游区域的可降解性浓度.除此之外, BDOCA (雨季)和BDOC28(雨季)与景观负荷比指数(LWLI)的正相关性较高, 表明雨季相比旱季于景观格局分布上更容易影响流域BDOC和BDOM的迁移转化.这可能因为雨季流域地表水和土壤各物质发生更复杂的相互交换过程要较旱季更为强烈[7].

蓝色方块表示正相关性, 红色方块表示负相关性, 大小和颜色深浅表示相关性的强弱; 1.林地用地占比, 2.水库用地占比, 3.农业用地占比, 4.建设用地占比, 5.LWLI, 6.LWLI.E, 7.LWLI.D, 8. LWLI.S, 9.BDOCA(旱季), 10.BDOC28(旱季), 11. BDOCA(雨季), 12. BDOC28(雨季), 13. BDOMC1, 14. BDOMC2 图 6 不同季节下BDOC、BDOM与源汇土地利用比例和景观负荷比指数相关矩阵分析 Fig. 6 Correlation matrix analysis of BDOC, BDOM, and source-sink land use ratio and landscape load ratio index in different seasons

综上, 可利用源汇景观模型探究景观格局对BDOM的影响, 通过合理的景观分布和格局优化来预防水体生态环境污染和控制水体污染物, 有助于城郊流域土地利用格局对BDOM生物可利用性影响研究, 对降低城乡复合生态系统水体水质污染、提升水质净化等多种生态系统服务功能具有重要意义.

4 结论

(1) 流域ρ(BDOC)为0.57~6.80 mg ·L-1, 在空间上BDOC浓度分布为:子流域4>子流域1>子流域2>子流域3, 在时间上BDOC浓度在旱雨季存在极显著差异(P < 0.001).

(2) 流域中BDOM的荧光组分通过EEM-PARAFAC得到2个组分, 蛋白质类C2组分略高于陆源腐殖质类C1的Fmax值, C1组分与BDOC空间分布一致.

(3) 流域内土地利用类型与BDOC和组分C1的BDOM的分布密切相关, 人类活动强度越大的区域其浓度越高.

(4) 流域源汇景观格局特征与BDOC和BDOM显著相关, “源”景观(城镇和农田)影响了BDOM在流域内的空间分布, 源汇景观格局影响了流域碳的环境行为和过程.

参考文献
[1] Zhao F K, Yang L, Chen L D, et al. Soil contamination with antibiotics in a typical peri-urban area in eastern China: Seasonal variation, risk assessment, and microbial responses[J]. Journal of Environmental Sciences, 2019, 79: 200-212. DOI:10.1016/j.jes.2018.11.024
[2] Cilliers S S, Bredenkamp G J. Vegetation of road verges on an urbanisation gradient in Potchefstroom, South Africa[J]. Landscape and Urban Planning, 2000, 46(4): 217-239. DOI:10.1016/S0169-2046(99)00057-2
[3] Wang Y Q, Moskovits D K. Tracking fragmentation of natural communities and changes in land cover: applications of Landsat data for conservation in an urban landscape (Chicago Wilderness)[J]. Conservation Biology, 2001, 15(4): 835-843. DOI:10.1046/j.1523-1739.2001.015004835.x
[4] Tang J F, Wang W D, Yang L, et al. Variation in quantity and chemical composition of soil dissolved organic matter in a peri-urban critical zone observatory watershed in Eastern China[J]. Science of the Total Environment, 2019, 688: 622-631. DOI:10.1016/j.scitotenv.2019.06.270
[5] Tang J F, Wang W D, Yang L, et al. Seasonal variation and ecological risk assessment of dissolved organic matter in a peri-urban critical zone observatory watershed[J]. Science of the Total Environment, 2020, 707. DOI:10.1016/j.scitotenv.2019.136093
[6] Massicotte P, Asmala E, Stedmon C, et al. Global distribution of dissolved organic matter along the aquatic continuum: across rivers, lakes and oceans[J]. Science of the Total Environment, 2017, 609: 180-191. DOI:10.1016/j.scitotenv.2017.07.076
[7] 马小亮, 刘桂民, 吴晓东, 等. 青藏高原多年冻土区典型植被下河流溶解性有机碳的生物可利用性[J]. 环境科学, 2018, 39(5): 2086-2094.
Ma X L, Liu G M, Wu X D, et al. Bioavailability of dissolved organic carbon in rivers for typical vegetation types in the Permafrost Regions on the Qinghai-Tibet Plateau[J]. Environmental Science, 2018, 39(5): 2086-2094. DOI:10.13227/j.hjkx.201709280
[8] Kellerman A M, Kothawala D N, Dittmar T, et al. Persistence of dissolved organic matter in lakes related to its molecular characteristics[J]. Nature Geoscience, 2015, 8(6): 454-457. DOI:10.1038/ngeo2440
[9] Logue J B, Stedmon C A, Kellerman A M, et al. Experimental insights into the importance of aquatic bacterial community composition to the degradation of dissolved organic matter[J]. The ISME Journal, 2016, 10(3): 533-545. DOI:10.1038/ismej.2015.131
[10] Baña Z, Ayo B, Marrasé C, et al. Changes in bacterial metabolism as a response to dissolved organic matter modification during protozoan grazing in coastal Cantabrian and Mediterranean waters[J]. Environmental Microbiology, 2014, 16(2): 498-511. DOI:10.1111/1462-2920.12274
[11] Pérez M T, Sommaruga R. Differential effect of algal-and soil-derived dissolved organic matter on alpine lake bacterial community composition and activity[J]. Limnology and Oceanography, 2006, 51(6): 2527-2537. DOI:10.4319/lo.2006.51.6.2527
[12] Roiha T, Peura S, Cusson M, et al. Allochthonous carbon is a major regulator to bacterial growth and community composition in subarctic freshwaters[J]. Scientific Reports, 2016, 6(1). DOI:10.1038/srep34456
[13] 杜嘉丹, 董秉直. 太湖微污染水中生物可降解有机碳的特点分析[J]. 中国给水排水, 2018, 34(13): 43-47.
Du J D, Dong B Z. Characteristics of biodegradable dissolved organic carbon in Taihu Lake[J]. China Water & Wastewater, 2018, 34(13): 43-47. DOI:10.19853/j.zgjsps.1000-4602.2018.13.009
[14] Labanowski J, Feuillade G. Combination of biodegradable organic matter quantification and XAD-fractionation as effective working parameter for the study of biodegradability in environmental and anthropic samples[J]. Chemosphere, 2009, 74(4): 605-611. DOI:10.1016/j.chemosphere.2008.09.063
[15] Stutter M I, Richards S, Dawson J J C. Biodegradability of natural dissolved organic matter collected from a UK moorland stream[J]. Water Research, 2013, 47(3): 1169-1180. DOI:10.1016/j.watres.2012.11.035
[16] Lønborg C, Davidson K, Álvarez-Salgado X A, et al. Bioavailability and bacterial degradation rates of dissolved organic matter in a temperate coastal area during an annual cycle[J]. Marine Chemistry, 2009, 113(3-4): 219-226. DOI:10.1016/j.marchem.2009.02.003
[17] Gu N T, Song Q B, Yang X L, et al. Fluorescence characteristics and biodegradability of dissolved organic matter (DOM) leached from non-point sources in southeastern China[J]. Environmental Pollution, 2020, 258. DOI:10.1016/j.envpol.2019.113807
[18] Chantigny M H. Dissolved and water-extractable organic matter in soils: a review on the influence of land use and management practices[J]. Geoderma, 2003, 113(3-4): 357-380. DOI:10.1016/S0016-7061(02)00370-1
[19] Zhou L, Zhou Y Q, Tang X M, et al. Biodegradable dissolved organic carbon shapes bacterial community structures and co-occurrence patterns in large eutrophic Lake Taihu[J]. Journal of Environmental Sciences, 2021, 107: 205-217. DOI:10.1016/j.jes.2021.02.011
[20] Zhang L Q, Zhu X H, Huang X, et al. Effects of land use and nutrients on the characteristics of dissolved organic matter in the Nanchong Section of Jialing River, China in December 2019[J]. Water Supply, 2022, 22(2): 1863-1875. DOI:10.2166/ws.2021.307
[21] Lyu L, Liu G, Shang Y X, et al. Characterization of dissolved organic matter (DOM) in an urbanized watershed using spectroscopic analysis[J]. Chemosphere, 2021, 277. DOI:10.1016/j.chemosphere.2021.130210
[22] 孙然好, 孙龙, 苏旭坤, 等. 景观格局与生态过程的耦合研究: 传承与创新[J]. 生态学报, 2021, 41(1): 415-421.
Sun R H, Sun L, Su X K, et al. A review on the coupling research of landscape patterns and ecological processes: consistency and innovation[J]. Acta Ecologica Sinica, 2021, 41(1): 415-421.
[23] 周添惠, 杨磊, 赵方凯, 等. 渭河源流域源汇景观演变对径流泥沙的影响[J]. 生态学报, 2022, 42(1): 58-66.
Zhou T H, Yang L, Zhao F K, et al. Impacts of source and sink landscape dynamics on runoff and sediment in the source watershed of the Weihe River[J]. Acta Ecologica Sinica, 2022, 42(1): 58-66.
[24] 孙然好, 陈利顶, 王伟, 等. 基于"源""汇"景观格局指数的海河流域总氮流失评价[J]. 环境科学, 2012, 33(6): 1784-1788.
Sun R H, Chen L D, Wang W, et al. Correlating landscape pattern with total nitrogen concentration using a location-weighted sink-source landscape index in the Haihe River Basin, China[J]. Environmental Science, 2012, 33(6): 1784-1788.
[25] 李敏, 唐剑锋, 陈利顶, 等. 城郊流域源汇景观格局与水体抗生素的关系[J]. 环境科学, 2020, 41(5): 2264-2271.
Li M, Tang J F, Chen L D, et al. Relationship between source-sink landscape pattern and antibiotics in surface water in peri-urban watershed[J]. Environmental Science, 2020, 41(5): 2264-2271. DOI:10.13227/j.hjkx.201911114
[26] 陈利顶, 傅伯杰, 赵文武. "源""汇"景观理论及其生态学意义[J]. 生态学报, 2006, 26(5): 1444-1449.
Chen L D, Fu B J, Zhao W W. Source-sink landscape theory and its ecological significance[J]. Acta Ecologica Sinica, 2006, 26(5): 1444-1449. DOI:10.3321/j.issn:1000-0933.2006.05.020
[27] Vonk J E, Tank S E, Mann P J, et al. Biodegradability of dissolved organic carbon in permafrost soils and aquatic systems: a meta-analysis[J]. Biogeosciences, 2015, 12(23): 6915-6930. DOI:10.5194/bg-12-6915-2015
[28] 马琦琦, 李刚, 魏永. 城郊关键带土壤中溶解性有机质的光谱特性及其时空变异[J]. 环境化学, 2020, 39(2): 455-466.
Ma Q Q, Li G, Wei Y. Spectral characteristics and spatiotemporal variation of DOM in Peri-urban Critical Zone[J]. Environmental Chemistry, 2020, 39(2): 455-466.
[29] Murphy K R, Hambly A, Singh S, et al. Organic matter fluorescence in municipal water recycling schemes: toward a unified PARAFAC model[J]. Environmental Science & Technology, 2011, 45(7): 2909-2916.
[30] Chen W, Westerhoff P, Leenheer J A, et al. Fluorescence excitation? emission matrix regional integration to quantify spectra for dissolved organic matter[J]. Environmental Science & Technology, 2003, 37(24): 5701-5710.
[31] Stedmon C A, Bro R. Characterizing dissolved organic matter fluorescence with parallel factor analysis: a tutorial[J]. Limnology and Oceanography: Methods, 2008, 6(11): 572-579. DOI:10.4319/lom.2008.6.572
[32] 何杰, 李学艳, 林欣, 等. 光谱特征法辨识不同污染景观河道中溶解性有机物的组分与来源[J]. 环境科学学报, 2021, 41(3): 1000-1010.
He J, Li X Y, Lin X, et al. Spectral feature method was used to identify the components and sources of dissolved organic matter in different polluted landscape channels[J]. Acta Scientiae Circumstantiae, 2021, 41(3): 1000-1010.
[33] Baker A, Curry M. Fluorescence of leachates from three contrasting landfills[J]. Water Research, 2004, 38(10): 2605-2613. DOI:10.1016/j.watres.2004.02.027
[34] Tang J F, Li X H, Yu X W, et al. Application of triolein-embedded cellulose acetate membrane (TECAM) passive sampler to study phase distribution of hexabromocyclododecanes (HBCDs) in sediment[J]. Environmental Science and Pollution Research, 2018, 25(21): 20823-20832. DOI:10.1007/s11356-018-2254-1
[35] So S H, Choi I H, Kim H C, et al. Seasonally related effects on natural organic matter characteristics from source to tap in Korea[J]. Science of the Total Environment, 2017, 592: 584-592. DOI:10.1016/j.scitotenv.2017.03.063
[36] Cioce D M, Aitkenhead-Peterson J A. Biodegradable dissolved organic carbon in urban and remnant soils in south-central Texas, USA[J]. Geoderma, 2015, 245-246: 52-55. DOI:10.1016/j.geoderma.2015.01.013
[37] Zhou Y Q, Yao X L, Zhou L, et al. How hydrology and anthropogenic activity influence the molecular composition and export of dissolved organic matter: observations along a large river continuum[J]. Limnology and Oceanography, 2021, 66(5): 1730-1742. DOI:10.1002/lno.11716
[38] Painter S C, Lapworth D J, Woodward E M S, et al. Terrestrial dissolved organic matter distribution in the North Sea[J]. Science of the Total Environment, 2018, 630: 630-647. DOI:10.1016/j.scitotenv.2018.02.237
[39] Amaral V, Graeber D, Calliari D, et al. Strong linkages between DOM optical properties and main clades of aquatic bacteria[J]. Limnology and Oceanography, 2016, 61(3): 906-918. DOI:10.1002/lno.10258
[40] Murphy K R, Ruiz G M, Dunsmuir W T M, et al. Optimized parameters for fluorescence-based verification of ballast water exchange by ships[J]. Environmental Science & Technology, 2006, 40(7): 2357-2362.
[41] Goździejewska A M, Skrzypczak A R, Koszałka J, et al. Effects of recreational fishing on zooplankton communities of drainage system reservoirs at an open-pit mine[J]. Fisheries Management and Ecology, 2020, 27(3): 279-291. DOI:10.1111/fme.12411
[42] 曹昌丽, 梁梦琦, 何桂英, 等. 城镇化河流溶解性有机质的荧光特性与水质相关性: 以宁波市北仑区芦江为例[J]. 环境科学, 2018, 39(4): 1560-1567.
Cao C L, Liang M Q, He G Y, et al. Fluorescent dissolved organic matter and its correlation with water quality in a urban river: A case study of the Lujiang River in Beilun, Ningbo[J]. Environmental Science, 2018, 39(4): 1560-1567.
[43] Bai L L, Cao C C, Wang C H, et al. Toward quantitative understanding of the bioavailability of dissolved organic matter in freshwater lake during cyanobacteria blooming[J]. Environmental Science & Technology, 2017, 51(11): 6018-6026.
[44] Mattsson T, Kortelainen P, Räike A. Export of DOM from boreal catchments: impacts of land use cover and climate[J]. Biogeochemistry, 2005, 76(2): 373-394.
[45] Hosen J D, McDonough O T, Febria C M, et al. Dissolved organic matter quality and bioavailability changes across an urbanization gradient in headwater streams[J]. Environmental Science & Technology, 2014, 48(14): 7817-7824.
[46] Mei Y, Hornberger G M, Kaplan L A, et al. The delivery of dissolved organic carbon from a forested hillslope to a headwater stream in southeastern Pennsylvania, USA[J]. Water Resources Research, 2014, 50(7): 5774-5796.
[47] Nguyen T M H, Le T P Q, Van Hoang V, et al. Seasonal and spatial variation in biodegradability of organic carbon along the Red River, Vietnam[J]. Carbon Management, 2021, 12(5): 549-557.
[48] 高凤, 邵美玲, 唐剑锋, 等. 城镇流域水体-沉积物中溶解性有机质的荧光特性及影响因素: 以宁波市小浃江为例[J]. 环境科学, 2019, 40(9): 4009-4017.
Gao F, Shao M L, Tang J F, et al. Fluorescence characteristics and influencing factors of dissolved organic matter (DOM) in water and sediment of urban Watershed: a case study of Xiaojia River in Ningbo City[J]. Environmental Science, 2019, 40(9): 4009-4017.