Long-term changes of CODMn flux into the sea: Retrieval of spectral information from remote sensing images of Oujiang River

  • LIU Yuening , 1, 2, 3 ,
  • GONG Fang , 2, 3, * ,
  • HE Xianqiang 1, 2, 3 ,
  • JIN Xuchen 2, 3
Expand
  • 1. Ocean College, Zhejiang University, Zhoushan 316021, China
  • 2. State Key Laboratory of Satellite Ocean Environment Dynamics, Hangzhou 310012, China
  • 3. Second Institute of Oceanography, MNR, Hangzhou 310012, China

Received date: 2022-08-16

  Revised date: 2022-12-15

  Online published: 2023-04-28

Abstract

As the second largest river in Zhejiang Province, Oujiang River has about 20 billion cubic meters of water flowing into the East China Sea every year. In order to explore the CODMn flux of Oujiang River into the sea and its characteristics of long-term changes based on the correlation between the measured data and the spectral information of Landsat-8 remote sensing image, the Oujiang River CODMn remote sensing inversion model was constructed (the average relative error is 28%) and the inversion results of CODMn in the middle and lower reaches of Oujiang River in 1986-2020 were obtained. In addition, based on rainfall data, the runoff and CODMn flux into the sea was further estimated and then remote sensing monitoring of long-term changes of Oujiang CODMn flux into the sea was realized. The results show that the difference of CODMn between dry and wet seasons in the Oujiang River Estuary is not obvious, but in general, the CODMn in wet season was slightly higher than that in dry season. From 1986 to 2020, the CODMn flux of Oujiang River into the sea fluctuated greatly, with a slight downward trend.

Cite this article

LIU Yuening , GONG Fang , HE Xianqiang , JIN Xuchen . Long-term changes of CODMn flux into the sea: Retrieval of spectral information from remote sensing images of Oujiang River[J]. Journal of Marine Sciences, 2023 , 41(1) : 45 -54 . DOI: 10.3969-j.issn.1001-909X.2023.01.004

0 引言

入海河流是连接陆地和海洋的重要桥梁,人类活动产生的工业废水、生活污水会通过河流输运进入近岸海域。作为水质监测和评估中使用最广泛的参数之一,化学需氧量(Chemical Oxygen Demand, COD)是指水中还原性物质在强氧化剂作用下发生氧化还原反应时消耗氧的量。目前,在水质监测、控制水污染和管理地表水等实际应用中,最常用的指标为高锰酸盐指数(Permanganate Index, CODMn)[1],该指标可表征水体中有机污染程度,常用于自然水体或洁净水体的测定,可有效测量水中可被氧化的还原性物质(有机物)总量。因此,CODMn可间接表征入海河流的有机碳含量,进而有助于估算河流入海有机碳的输运通量[2-4]
相较于传统的基于水化学原理的现场监测方法,遥感技术具有可大范围同步观测、历史可回溯等突出优势。通过卫星遥感影像,结合实测数据,可进一步构建水质反演模型,实现入海河流水质的大范围连续动态监测。自20世纪70年代以来,遥感技术已广泛应用于河流、湖泊水质参数的反演、监测和长时序分析[5]。例如,解启蒙 等[6]通过分析遥感反射率和实测数据之间的相关性,基于Landsat-8卫星数据,建立了具有较高精度的清河水库CODMn遥感反演模型,为水库水质监测提供了可靠数据支撑;王丽艳 等[7]基于实测数据和MODIS遥感数据,估算并分析了呼伦湖COD的时空变化特征,为呼伦湖水环境质量管理提供了科学依据;王歆晖 等[8]利用哨兵2号遥感数据,结合实测数据,构建了基于综合水质指标(包含CODMn等参数)的水质遥感反演模型,为河流水质监管提供了依据;吴铭 等[9]利用1999—2018年的Landsat系列遥感影像数据和实测数据,建立了北极Lena河的水质反演模型,并分析了河流水质长期演变过程。但目前国内外关于污染物遥感反演的研究主要面向内陆湖泊或大型河流等,针对中小河流入海污染物的遥感评估研究相对较少[10]
本研究以瓯江为例,基于Landsat系列卫星遥感数据,结合实测水质数据,构建了瓯江中下游河段CODMn遥感反演模型,分析了瓯江中下游河段CODMn的时空变化特征。在此基础上,进一步结合流域降水数据估算分析了1986—2020年瓯江CODMn入海通量的长期演变特征。

1 数据与方法

1.1 研究区域概况

瓯江位于浙江省南部,为浙江省第二大河流,从温州湾入海。其干流全长388 km,流域面积18 165 km2,平均年降水量为1 745.8 mm,年径流量为202.7亿m3。瓯江为典型的山溪性河流,径流年内分配具有明显的季节性,年内径流主要由4—9月的汛期贡献(占全年径流量的78%以上),最大径流量主要出现在6月。按照径流量大小将4—9月划为湿季,其余月份为干季(① 温州市水利电力勘测设计院.温州市第三次水资源调查评价报告.2019.)。瓯江多年平均流量为446 m3/s,径流携带上游各县市乡镇生活污水和工农业废水进入河口海域[11]

1.2 实验数据

1.2.1 实测数据

2021年9月9日—10日(湿季)和2021年12月7日—9日(干季)分别在瓯江中下游河段开展现场调查,站位设置如图1所示。参照《中国海洋生态环境状况公报》[12]对瓯江近海水体的划分,在瓯江入海口划定了一块上承七都岛、下接灵昆岛的“X”型区域,作为瓯江入海水体范围(图1),以其范围内CODMn均值作为瓯江入海CODMn
图1 站位设置

Fig.1 Station setting

实验严格参照《地表水和污水监测技术规范》[13]要求执行。为避免采样点在遥感影像中受到混合像元影响,沿岸采集水样时选择尽可能远离岸边的位置作为实际采样点,同时通过手持GPS仪记录实际采样点坐标。采集的水样于实验当天按照《水和废水监测分析方法》[14]中的“高锰酸盐指数-酸性法”进行分析处理。9月(湿季)和12月(干季)两次调查获得的实测CODMn样本数分别为27个和15个。其中,靠近口门处的OJ01、OJ02和OJ27号点位受海水影响,盐度偏高,“高锰酸盐指数-酸性法”不再适用于上述点位水样[15],建模时筛除了相应数据。
为分析同一点位CODMn在干、湿季节的差异,选取两次采样均覆盖到的站点进行对比(表1),结果表明,研究区域同一采样点的CODMn在干、湿季差异并不显著,但总体上来说,湿季的CODMn略大于干季。
表1 瓯江中下游河段现场调查CODMn结果

Tab.1 CODMn results of field survey in the middle and lower reaches of Oujiang River

站点名称 CODMn/(mg·L-1)
9月湿季航次 12月干季航次
OJ10 1.26 1.10
OJ09 1.37 0.80
OJ13 0.98 1.10
OJ14 1.63 1.10
OJ08 1.89 1.00
OJ15 1.69 1.00
OJ16 2.01 1.40
OJ17 2.28 1.10
OJ19 2.20 1.30
OJ20 3.12 3.10
OJ03 5.47 5.20
最大值 6.67 6.20
最小值 0.98 0.80
平均值 2.84 2.21
标准差 1.60 1.87

1.2.2 遥感数据

本研究选择美国地质调查局公布的陆地卫星(Landsat)高空间分辨率(30 m)多光谱遥感数据进行分析,遥感数据从美国地质调查局网站(https://earthexplorer.usgs.gov/)下载获取。本研究拟选择1986—2020年瓯江河道上无云覆盖或云量小于30%的Landsat影像数据,有243幅Landsat Level-1级遥感数据符合上述条件(表2)。其中,数据选择时确保每年干、湿季至少有1幅遥感影像可用,影像基本覆盖了各个月份、季节、年份,具有一定的代表性,以满足分析不同水文条件下CODMn时空动态分布状况的需要。
表2 研究使用的遥感影像数据

Tab.2 Remote sensing data used in the study

卫星名称 遥感器 时间跨度 图幅数量/幅
Landsat-5 TM 1986-01-07—2011-10-27 134
Landsat-7 ETM+ 1999-10-02—2003-01-14 21
Landsat-7 ETM+(修复条带后) 2003-07-25—2021-01-15 42
Landsat-8 OLI 2013-04-23—2021-12-09 46
本研究中,建模时使用2021年9月4日和2021年12月9日的Landsat-8卫星影像,用于匹配现场实测的CODMn数据。
本研究使用ENVI 5.3软件对Landsat系列遥感数据进行预处理,主要包括辐射定标和大气校正。其中,辐射定标通过Radiometric Calibration工具实现,将辐亮度值输出为BIL格式;大气校正通过FLAASH大气校正算法实现,并根据干、湿季节选择不同的气溶胶模型。基于获得的遥感反射率数据,进一步提取采样点的遥感反射率信息,提取时以采样点所在像元为中心,选择该中心像元周围3×3窗口范围内反射率值最低的像元作为该采样点对应的遥感反射率信息进行取样,以降低混合像元影响。

1.2.3 降水数据

本研究使用的水文数据来源于PENG et al[16]提供的基于CRU(Climatic Research Unit)气候数据集和World Clim气候数据集生成的中国地区月均降水数据。经过区域筛选和计算处理,获得瓯江流域月尺度降水量信息,并在此基础上计算得到瓯江流域年均降水量。

1.3 遥感模型的构建

1.3.1 相关性分析

对采样点临近像元各波段的反射率与实测CODMn数据进行敏感性分析。使用SPSS 25软件分析拟合单波段、多波段组合与CODMn之间的相关性,寻找皮尔逊相关系数较高的波段组合,以便后续得到拟合度更高的回归曲线和反演模型。

1.3.2 反演模型

以相关性较高的波段组合作为自变量,CODMn实测值作为因变量,建立回归方程,得到瓯江中下游河段CODMn的遥感反演模型。
经过相关性分析并对比决定系数(R2),最终筛选得到如表3所示的3种R2较高的波段组合形式和遥感反演模型(上述模型均已通过显著性检验)。
表3 部分波段组合、模型形式及模型决定系数

Tab.3 Partial band combination results and R2 of the models

波段组合 模型形式 R2 p n
ln(Blue+NIR)/ln(Red) Y=816.827-1 484.774X+675.798X2 0.523 <0.001 42
ln(Green×NIR)/ln(Red+NIR) Y=3 368.372-2 742.496X+269.527X3 0.577 <0.001 42
ln(Blue×NIR)/ln(Red+NIR) Y=1 706.314-1 400.998X+140.141X3 0.611 <0.001 42
比较了各模型在瓯江中下游河段的CODMn反演误差,最终确定使用ln(Blue+NIR)/ln(Red)的波段组合形式构建遥感反演模型:
Y=816.827-1 484.774X+675.798X2
基于上述遥感反演模型,结合遥感数据和水体光谱信息,可进一步反演瓯江中下游河段CODMn,从而获取CODMn时空变化特征。

1.3.3 反演模型验证

通过温州环境监测站(七都岛,2017年)和中国环境监测总站(龙湾,2018—2021年)提供的CODMn实测数据对模型进行精度验证。
基于前文构建的定量遥感反演模型,获得各个站点的CODMn反演值,并与实测值进行对比分析。验证结果表明,反演模型的平均绝对误差为0.63 mg/L,平均相对误差为28%,在可接受范围内。

1.4 CODMn入海通量计算方法

利用下式进行CODMn入海通量估算:
T=C×Q
式中:T为入海通量,C为河流入海处CODMn(即图1中“X”型区域的反演值均值),Q为河流入海径流量。
相关文献表明,瓯江径流量与流域降水量之间存在显著线性关系[17]。因此,本研究使用LIU et al[18]的研究结果估算瓯江径流量(R2=0.98,p<0.005):
Qi=0.061 6×P-22.87
Q=Qi×S
式中:Qi为地表年径流量,P为年降水量,Q为河流入海径流量,S为流域面积。
利用1998—2020年《浙江省水资源公报》[19]中钱塘江入海径流量监测数据对径流量估算公式进行验证。验证结果(图2)显示,径流量估算值与监测值之间的绝对误差均值为69.09亿m3,相对误差均值为13%。误差较大的年份有2015年、2019年和2020年,对应的相对误差分别为22%、25%和23%。整体误差小于15%,在可接受范围内。
图2 钱塘江入海径流量计算值与实际监测值

Fig.2 Comparison of estimated and in situ runoff data of Qiantang River

受降水量年际差异的影响,若仅从CODMn入海通量变化角度分析瓯江入海水质演变会存在一定误差。因此,本研究引入“比值”的概念:将瓯江1986—2020年各年的CODMn入海通量与入海径流量之比作为当年瓯江入海水体水质状况优劣的参考依据,以探讨瓯江流域水质状况的年际变化。该概念可以理解为“单位入海水体中所含CODMn”,即消除径流量年际差异后的瓯江入海水体CODMn变化情况。

2 结果与分析

2.1 瓯江CODMn空间分布特征

基于遥感反演模型,获得1986—2020年瓯江中下游河段干、湿季CODMn长时序反演结果(图3图4)。在此基础上绘制了1986—2020年瓯江中下游河段干、湿季CODMn反演值的多年平均空间分布状况(图5)。
图3 瓯江中下游河段CODMn空间分布(1986—2020年干季)

Fig.3 Spatial distribution of CODMn in the middle and lower reaches of Oujiang River (dry season 1986-2020)

图4 瓯江中下游河段CODMn空间分布(1986—2020年湿季)

Fig.4 Spatial distribution of CODMn in the middle and lower reaches of Oujiang River (wet season 1986-2020)

图5 瓯江中下游河段CODMn多年平均空间分布(1986—2020年)

Fig.5 Multi-year average spatial distribution of CODMn in the middle and lower reaches of Oujiang River (1986-2020)

从遥感反演的CODMn空间分布状况(图3图4)看,瓯江中下游河段CODMn在各年份的空间分布没有明显差异,结合图1可以得出,瓯江中下游河段CODMn空间分布整体呈现温溪—梅岙段低、梅岙—龙湾/黄华段高的特征。
从瓯江中下游河段CODMn多年平均空间分布(图5)可见,瓯江中下游河段CODMn的多年平均值在干、湿季的差异并不显著,并呈现出高度的一致性。温溪—梅岙段森林植被覆盖程度较高,人类活动干扰程度低,因此温溪—梅岙段整体水质状况良好,水体清澈,CODMn总体约为1 mg/L;随着河道途经的乡村、城市数量增加,人类活动产生的生活污水、废水排入瓯江,水体中的CODMn逐渐升高,至灵昆岛附近区域,CODMn已增至约3 mg/L;但由于海水的稀释作用,CODMn在入海口附近略有下降。

2.2 瓯江CODMn入海通量估算

根据公式(2)、(4)估算瓯江CODMn年入海通量,结果见图6
图6 1986—2020年瓯江CODMn年入海通量及其与入海径流量的比值

Fig.6 The annual CODMn flux of Oujiang River into the sea from 1986 to 2020 and the ratio of CODMn flux into the sea with the runoff

结果显示,1986—2020年瓯江CODMn入海通量年际变化较大,整体呈现波动变化趋势。CODMn入海通量与径流量比值的变化趋势总体与CODMn入海通量变化趋势一致。1986—1995年,瓯江入海水体CODMn变化波动较大;1996—2011年,CODMn开始趋于稳定,CODMn的年际差异趋于平缓;2012年之后,CODMn又开始趋于波动变化。但总体来说,瓯江水质在1986—2020年之间的变化并不显著,并有稳中向好的变化趋势。
结合公报数据[19-21],对照CODMn入海通量与径流量比值发现:1987年、2005年、2013年、2014年等年份由于受台风影响,瓯江沿岸发生诸多洪涝灾害,使得近海水质显著下降;2001年,瓯江流域废水、污水排放量同比增长0.69亿 t,水质下降,该年份对应的比值出现极值;2004年降水量明显低于常年平均降水量,较多年平均降水量减少了12.8%,属于枯水年,径流量较小,水质相较2003年有所下降,比值也有所上升。
在没有人为干预的理想情况下,河流径流量和CODMn入海通量通常有着很好的相关性[20],但随着人类生产、生活、建设活动的加剧,对瓯江中下游河段CODMn的影响逐渐加强。结合图6,可进一步发现:1986—1994年瓯江CODMn入海通量与径流量比值的平均值约为0.34,处于相对较高的水平,其中,最小值为0.16(1986年),最大值为0.54(1991年),标准差为0.12,在该时段内瓯江入海水质指标CODMn呈现出较大的波动变化。1994年后,随着人们对环境问题的重视程度逐步提升以及各项环保措施的颁布和实施,瓯江水质有所提升,CODMn入海通量和径流量的比值整体呈现下降趋势。1995—2020年,瓯江CODMn入海通量与径流量比值的平均值降为0.29,标准差为0.07,其中,最小值为0.14(1995年),最大值为0.44(2014年),该时期瓯江入海水质状况相对稳定,与1986—1994年相比,水质得到改善。2018—2020年,瓯江CODMn入海通量与径流量比值逐年下降,到2020年已降至0.24左右。总地来说,瓯江CODMn入海总量呈现出波动降低的趋势。

3 结论

本文以瓯江作为研究对象,基于Landsat系列卫星遥感数据,结合实测数据,构建了瓯江CODMn遥感反演算法,并进一步结合降水量数据估算径流量和CODMn入海通量,进而实现了瓯江CODMn入海通量长时序变化遥感监测,研究结果表明。
1)瓯江中下游河段CODMn在干、湿季的差异并不显著,但在大部分年份,湿季的CODMn略大于干季。
2)1986—2020年瓯江中下游河段CODMn多年平均空间分布显示:温溪—梅岙河段CODMn普遍较低,保持在1 mg/L左右;梅岙—灵都岛河段,CODMn逐渐累积升高,到灵昆岛附近,CODMn升至3 mg/L左右。
3)1986—2020年瓯江CODMn入海通量长时序反演结果表明,瓯江年均CODMn入海通量整体波动较大,呈现降低—升高—降低的演变趋势。
4)1986—1994年瓯江CODMn入海通量与径流量比值的平均值约为0.34,处于相对较高的水平;1994年后,瓯江CODMn入海通量和径流量的比值开始下降,1995—2020年比值均值降低至0.29左右,说明瓯江CODMn入海总量呈现波动降低的趋势。
[1]
YANG B, LIU Y P, OU F P, et al. Temporal and spatial analysis of COD concentration in East Dongting Lake by using of remotely sensed data[J]. Procedia Environmental Sciences, 2011, 10: 2703-2708.

DOI

[2]
夏圣斌, 张慧霖, 朱世凯. 昆山吴淞江流域地表水中TOC与高锰酸盐指数的相关性分析[J]. 智慧农业导刊, 2021, 1(3):40-42.

XIA S B, ZHANG H L, ZHU S K. Correlation analysis between TOC and permanganate index in surface water of Wusong River Basin in Kunshan[J]. Journal of Smart Agriculture, 2021, 1(3): 40-42.

[3]
陈坤, 李捷, 徐庆彤, 等. 渤海西部陆域河流CODMn与TOC浓度分布及相关性研究[J]. 湖北农业科学, 2016, 55(8):1926-1930.

CHEN K, LI J, XU Q T, et al. CODMn and TOC concentration distribution and correlation studies of different rivers in the west of Bohai Sea[J]. Hubei Agricultural Sciences, 2016, 55(8): 1926-1930.

[4]
陈光, 刘廷良, 孙宗光. 水体中TOC与COD相关性研究[J]. 中国环境监测, 2005, 21(5):9-12.

CHEN G, LIU T L, SUN Z G. Research of the correlation between TOC and COD in water[J]. Environmental Monitoring in China, 2005, 21(5): 9-12.

[5]
杨一鹏, 王桥, 王文杰, 等. 水质遥感监测技术研究进展[J]. 地理与地理信息科学, 2004, 20(6):6-12.

YANG Y P, WANG Q, WANG W J, et al. Application and advances of remote sensing techniques in determining water quality[J]. Geography and Geo-Information Science, 2004, 20(6): 6-12.

[6]
解启蒙, 林茂森, 杨国范, 等. 清河水库水体高锰酸盐指数遥感反演模型研究[J]. 中国农村水利水电, 2017(10):57-61.

XIE Q M, LIN M S, YANG G F, et al. Study on remote sensing inversion model of permanganate index in Qinghe Reservoir[J]. China Rural Water and Hydropower, 2017(10): 57-61.

[7]
王丽艳, 史小红, 孙标, 等. 基于MODIS数据遥感反演呼伦湖水体COD浓度的研究[J]. 环境工程, 2014, 32(12):103-108.

WANG L Y, SHI X H, SUN B, et al. Determination of COD concentration of water in Hulun Lake based on MODIS data[J]. Environmental Engineering, 2014, 32(12): 103-108.

[8]
王歆晖, 田华, 季铁梅, 等. 哨兵2卫星综合水质指标的河流水质遥感监测方法[J]. 上海航天(中英文), 2020, 37(5):92-97,104.

WANG X H, TIAN H, JI T M, et al. Remote sensing monitoring method for comprehensive water quality index in rivers based on Sentinel-2 satellite[J]. Aerospace Shanghai (Chinese & English), 2020, 37(5): 92-97, 104.

[9]
吴铭, 黄珏, 宫丽娇, 等. 长时序Landsat的北极Lena河DOC浓度变化及驱动力分析[J]. 遥感学报, 2021, 25(3):830-845.

WU M, HUANG J, GONG L J, et al. Analysis of DOC concentration variation and driving forces in the Arctic River Lena based on long-term Landsat time series[J]. National Remote Sensing Bulletin, 2021, 25(3): 830-845.

[10]
TOPP S N, PAVELSKY T M, JENSEN D, et al. Research trends in the use of remote sensing for inland water quality science: Moving towards multidisciplinary applications[J]. Water, 2020, 12(1): 169.

DOI

[11]
李彦伟, 刘菁, 侯庆志, 等. 瓯江江心屿河段分流比演变及其影响因素[J]. 人民长江, 2020, 51(11):22-27,133.

LI Y W, LIU J, HOU Q Z, et al. Evolution and influence factors of diversion ratio of Jiangxinyu reach in Oujiang River[J]. Yangtze River, 2020, 51(11): 22-27, 133.

[12]
中华人民共和国生态环境部. 2021年中国海洋生态环境状况公报[R/OL].[2022-06-30]. https://www.mee.gov.cn/hjzl/sthjzk/jagb/.

Ministry of Ecology and Environment of the People’s Republic of China. Bulletin on China’s marine ecological environment in 2021[R/OL].[2022-06-30]. https://www.mee.gov.cn/hjzl/sthjzk/jagb/.

[13]
国家环境保护总局. 地表水和污水监测技术规范:HJ/T 91—2002[S]. 北京: 中国环境出版社, 2002.

State Environmental Protection Administration. Technical specifications requirements for monitoring of surface water and waste water: HJ/T 91—2002[S]. Beijing: China Environment Publishing House, 2002.

[14]
国家环境保护总局《水和废水监测分析方法》编委会. 水和废水监测分析方法[M]. 第四版. 北京: 中国环境科学出版社, 2002.

The Editorial Committee of water and wastewater monitoring and analysis method of State Environmental Protection Administration. Water and wastewater monitoring and analysis method[M]. Fourth Edition. Beijing: China Envi-ronmental Science Press, 2002.

[15]
LI Z R, SHENG Y Q, SHI W J, et al. Influence of salinity on COD measurements in coastal water management[J]. Desalination and Water Treatment, 2016, 57(39): 18338-18345.

DOI

[16]
PENG S Z, DING Y X, LIU W Z, et al. 1 km monthly temperature and precipitation dataset for China from 1901 to 2017[J]. Earth System Science Data, 2019, 11(4): 1931-1946.

DOI

[17]
王东. 广通河流域降水径流关系判断与变异特性研究[J]. 甘肃水利水电技术, 2021, 57(11):4-7.

WANG D. Study on the relationship between precipitation and runoff and its variation characteristics in Guangtong River Basin[J]. Gansu Water Resources and Hydropower Technology, 2021, 57(11): 4-7.

[18]
LIU D, BAI Y, HE X Q, et al. Changes in riverine organic carbon input to the ocean from mainland China over the past 60 years[J]. Environment International, 2020, 134:105258.

DOI

[19]
浙江省水文管理中心. 浙江省水资源公报[R/OL].1998—2020[2022-06-30]. https://www.zjsw.cn/pages/nav.jsp?catId=1029.

Zhejiang Provincial Hydrological Management Center. Zhejiang water resources bulletin[R/OL]. 1998-2020[2022-06-30]. https://www.zjsw.cn/pages/nav.jsp?catId=1029.

[20]
浙江省生态环境厅. 浙江省环境状况公报[R/OL].1997—2020[2022-06-30]. http://sthjt.zj.gov.cn/col/col1201912/index.html.

Zhejiang Provincial Department of Ecological Environment. Zhejiang environmental status bulletin[R/OL]. 1997-2020[2022-06-30]. http://sthjt.zj.gov.cn/col/col1201912/index.html.

[21]
温州市生态环境局. 温州市环境状况公报[R/OL].2005—2020[2022-06-30]. http://sthjj.wenzhou.gov.cn/col/col1229561926/index.html.

Wenzhou Ecological Environment Bureau. Wenzhou environ-mental status bulletin[R/OL]. 2005-2020 [2022-06-30]. http://sthjj.wenzhou.gov.cn/col/col1229561926/index.html.

Outlines

/