DB31T+1691-2026城市绿地碳汇遥感监测技术要求

ICS65.020.403131DB31/T1691—2026前言 2规范性引用文件 3术语和定义 4监测方法与内容 监测方法 监测内容 5数据选择与处理 遥感影像数据 气象数据 绿地植被类型数据 6碳汇计算与制图 植被净初级生产力计算 植被碳汇计算 结果统计 空间分布制图 附录A（规范性）净初级生产力计算 6A.1总则 A.2净初级生产力计算 附录B（资料性）不同植被类型植被指数参考值 9附录C（资料性）不同植被类型最大光能利用率参考值 10参考文献 DB31/T1691—2026本文件按照GB/T1.1—2020《标准化工作导则第1部分：标准化文件的结构和起草规则》的规定起草。请注意本文件的某些内容可能涉及专利。本文件的发布机构不承担识别专利的责任。本文件由上海市绿化和市容管理局提出并组织实施。本文件由上海市园林绿化标准化技术委员会归口。本文件起草单位：上海市园林科学规划研究院、上海市绿化管理指导站、上海市测绘院。本文件主要起草人：张浪、张桂莲、邢璐琪、林勇、郑谐维、朱春玲、张雯、易扬、张希金、张艳华、仲启铖、林奕成、王歆晖、张琪。11范围GB/T15968遥感影像平面图制作规范城市绿地植被urbangreenspacevegetation城市绿地内各种植被类型的总称，包括常绿阔叶林、落叶阔叶林、常绿针叶林、落叶针叶林、针阔混交林、竹林、灌丛和草地等。城市绿地植被碳汇carbonsinkofurbangreenspacevegetation像元pixel2空间分辨率spatialresolution波谱spectrum波谱分辨率spectralresolution绿色植物在单位时间、单位面积内固定的有机物总量减去植物自养呼吸消耗后的净积累量，也称净注：一般把400nm～700nm的太阳辐射称为光合有效辐射。3植被对波长在400nm~700nm间太阳辐射能量的吸收比率，描述太阳光在冠层辐射传输过程中植被吸收比例。一系列利用遥感影像不同谱段数据的线性或非线性组合而形成的能反映绿色植物的生长状况和分布的特征指数。利用气象数据、遥感影像数据、绿地植被类型数据，基于光能利用率模型，核算城市绿地植被碳汇总量，制作城市绿地植被碳汇空间分布图。5.1.1基本要求基于遥感影像数据开展城市绿地碳汇监测，主要数据来源包括但优于10m;蓝波段)、近红外(760nm~1100nm)、短波红外(1550nm～1750nm)波段；5.1.2.1影像预处理：包括辐射定标、大气校正、几何校正、图像融合、图像增强(去噪声)、裁剪、镶嵌等。4NPP(x,t)城市绿地植被碳汇总量计算见公式(2)。DB31/T1691—2026总——监测期内，城市绿地碳汇总量，单位为克二氧化碳（gCO2——监测周期内的监测间隔总数，如以一年为监测周期，监测间隔为季度，则监测间隔总——监测周期内的第个监测间隔，如以一年为监测周期，监测间隔为季度，则监测间隔为第季度；——城市绿地空间数据的像元总数；——城市绿地空间数据的第个像元；,——像元在时期（月度/季度/年度）的碳汇量，单位为克二氧化碳（gCO2）。空间分布制图利用空间数据处理工具对城市绿地碳汇空间分布进行专题图制作，图面信息包括但不限于：——行政界线、监测区；——图面注记（图名、图例、监测时段、数字和线划比例尺、指北针、公里网格等——图廓整饰（内图廓、外图廓、坐标注记），按照GB/T15968、CH/T3019执行。DB31/T1691—2026（规范性）净初级生产力计算A.1总则基于光能利用率模型，利用遥感数据、气象数据等，根据植被吸收的光合有效辐射和光能利用率，计算净初级生产力。A.2净初级生产力计算A.2.1净初级生产力城市绿地监测范围内，净初级生产力计算见公式（A.1）。,=,×(,)·······································(A.1),——在监测间隔期，像元的净初级生产力，单位为克碳（gC(,)——在监测间隔期，像元吸收的光合有效辐射，单位为兆焦（MJ,——在监测间隔期，像元的实际光能利用率，单位为克碳每兆焦（gC/MJ）。A.2.2吸收光和有效辐射A.2.2.1吸收光合有效辐射计算见公式（A.2）。(,)=(,)×0.5×(,)································(A.2)(,)——在监测间隔期，像元吸收的光合有效辐射，单位为兆焦（MJ,——在监测间隔期，像元的太阳总辐射量，单位为兆焦（MJ0.5——常数，植被所能利用的太阳有效辐射占太阳总辐射的比例；(,)——在监测间隔期，像元的植被层对入射光合有效辐射的吸收比。A.2.2.2光合有效辐射的吸收比计算见公式（A.3～A.7）。,=,,+,,/2···························(A.3),=+··················(A.4),=+······················(A.5),=,,((,,))···········································(A.6),=···············································(A.7),——在监测间隔期，像元的植被层对入射光合有效辐射的吸收比；DB31/T1691—2026,——在监测间隔期，像元的归一化植被指数；,——第植被类型最小NDVI，取值参见表B.1；或对应各植被类型NDVI值5%处的百分位数；——常数，取值为0.95；——常数，取值为0.001；,——第植被类型最大NDVI，取值参见表B.1；或对应各植被类型NDVI值95%处的百分位数；,——在监测间隔期，像元的比值植被指数；,——第植被类型最小SR，取值参见表B.1；或基于各植被类型NDVI值5%处的百分位数和SR公式计算得到；,——第植被类型最大SR，取值参见表B.1；或基于各植被类型NDVI值95%处的百分位数和SR公式计算得到；,——在监测间隔期，像元的近红外波段反射率数据；(,)——在监测间隔期，像元的红波段反射率数据。A.2.3光能利用率A.2.3.1光能利用率计算见公式（A.8）。,=1,×2,×,×,·······························(A.8),——在监测间隔期，像元的实际光能利用率，单位为克碳每兆焦（gC/MJ1,、2,——在监测间隔期，像元的温度胁迫指数；,——在监测间隔期，像元的水分胁迫指数；,——第植被类型的最大光能利用率，单位为克碳每兆焦（gC/MJ参考值参见表C.1。A.2.3.2温度胁迫指数计算见公式（A.9～A.10）。1,=0.8+0.02×−0.0005×[2]·························(A.9)2,=1.184+1x[.×.3((−10−0(+,)(],))]······························(A.10)1,、2,——在监测间隔期，像元的温度胁迫指数；——研究区域内一年之中NDVI值达到最高时当月的平均气温，单位为摄氏度(,)——在监测间隔期，像元的平均气温，单位为摄氏度(℃)。A.2.3.3水分胁迫指数计算见公式（A.11～A.12）。,=··············································(A.11)=,,((,,))···········································(A.12),——在监测间隔期，像元的水分胁迫指数；——地表水分指数，描述植被叶片水分含量；DB31/T1691—2026——单个像元生长期内的LSWI最大值；,——在监测间隔期，像元的近红外波段反射率数据；,——在监测间隔期，像元的短波红外波段反射率数据。DB31/T1691—2026（资料性）不同植被类型植被指数参考值表B.1给出了不同植被类型植被指数参考值。表B.1不同植被类型植被指数参考值DB31/T1691—2026（资料性）不同植被类型最大光能利用率参考值表C.1给出了不同植被类型的最大光能利用率参考值。表C.1不同植被而类型最大光能利用率参考值DB31/T1691—2026[1]GB/T14950—2009摄影测量与遥感术语[2]GB/T34815—2017植被生态质量气象评价指数[3]GB/T41280—2022卫星遥感影像植被覆盖度产品规范[4]GB/T41281—2022光合有效辐射遥感产品真实性检验[5]GB/T51346—2019城市绿地规划标准[6]CJJ/T85—2017城市绿地分类标准[7]LY/T3253—2021林业碳汇计量监测术语[8]朱文泉,潘耀忠,何浩,于德永,扈海波.中国典型植被最大光利用率模拟[J].科学通