大气氮沉降是大气中各类含氮组分通过干、湿过程向下垫面传输沉积的现象，是全球氮循环的关键环节，其核心理论梳理可为后续研究明确科学逻辑与作用边界。
大气中活性氮分为气态与颗粒态两类，不同组分理化性质差异直接决定其迁移转化路径：排放进入对流层的含氮污染物，会经历光化学反应、气粒转化、云化学转化等过程，改变活性氮的存在形态分布，直接影响后续沉降的发生概率与速率。
干沉降是含氮组分在湍流扩散、布朗运动、重力作用下直接沉积的过程，强度与湍流条件、下垫面特征、组分粒径高度相关；湿沉降是含氮组分经云内冲刷、云下清除随降水沉降的过程，通量由降水强度、组分浓度共同决定。
传统地面观测仅能获取单点离散的大气氮沉降数据，无法支撑大区域连续空间的氮沉降特征研究，而卫星遥感具备大范围同步观测、空间连续覆盖、可追溯历史变化的天然优势，是区域尺度估算的核心支撑。
当前可用遥感数据源分为两类：极轨卫星可提供千米级对流层活性氮柱浓度逐日产品，匹配干沉降估算的浓度输入需求；气象卫星可提供高分辨率降水产品，满足湿沉降对降水输入的要求。现有卫星产品仍存在精度缺陷，本研究将以极轨卫星产品为核心，结合气象卫星降水产品，引入地面观测校准系统误差，充分发挥遥感技术优势。
本研究选取我国东部典型城市群为核心研究区，全域面积约18.3万平方公里，覆盖平原农业区、低山丘陵生态区与城镇密集区三类功能单元，兼具高强度排放与多样化下垫面特征，可充分检验模型适配性、校准模型参数。统一采用WGS-84坐标系与UTM分度带投影，消除空间匹配误差。 研究时段为2018-2022年共5个完整自然年，以月为基础统计单元，可获取连续稳定的Sentinel-5P TROPOMI观测数据，满足氮沉降时序演变分析需求。 本研究关注干、湿沉降中的气态二氧化氮、气态硝酸、气态氨气、颗粒态铵盐与硝酸盐五类活性氮组分，剔除低贡献惰性氮气，兼顾估算完整性与计算效率。
卫星遥感反演大气活性氮浓度，依赖不同含氮组分对特定波段辐射的特征吸收原理，原始观测信号受传感器偏差、大气效应与仪器噪声干扰，必须经规范化校正才能保障反演精度。 先开展辐射定标，本研究结合卫星配套星上定标参数与交叉定标结果，将传感器原始DN值转换为高精度大气顶部辐射亮度，消除传感器响应衰减与漂移误差。再基于MODTRAN6模型结合实地观测开展大气校正，剔除大气干扰与条带噪声，校正后数据信噪比提升18%，为后续反演提供可靠基础。
下垫面属性是大气氮沉降过程的核心调控因子，本研究整合三类下垫面基础数据：采用研究区同期30米分辨率Landsat地表覆盖分类产品，结合Google Earth高分辨率影像修正模糊分类区域；采用30米分辨率SRTM数字高程模型，提取海拔、坡度、坡向核心参数并计算局地粗糙度。所有数据统一投影至WGS-84坐标系，分类数据采用最邻近法、连续参数采用三次卷积，重采样至1千米分辨率匹配Sentinel-5P氮浓度栅格，生成匹配性良好的下垫面属性数据集。
大气氮干湿沉降估算对气象要素的时空精度要求较高，干沉降速率计算依赖近地面风速、大气边界层稳定度等参数，湿沉降估算则与降水特征直接相关，多源气象数据存在空间不连续、系统偏差不一致等问题，需通过处理生成连续可靠的输入数据集。
本研究整合研究区国家级气象站逐时原位观测、欧洲中期天气预报中心再分析格点数据、FY卫星反演气象产品三类数据，以集合卡尔曼滤波算法按1日时间窗口滚动开展时序同化，用高精度站点观测动态更新格点与卫星产品，补充缺测数据。
同化后利用独立站点观测构建分要素偏差校正方程，最终重采样得到研究区2018-2022年空间连续的1km分辨率逐日气象栅格数据集，为后续估算提供可靠输入。
从卫星遥感光谱信号中提取与近地面大气氮污染物高度关联的特征参量，是构建干沉降遥感估算模型的核心前提，直接决定氮浓度反演与干沉降通量估算的精度。本研究以Sentinel-5P TROPOMI观测数据为核心，该传感器在紫外-可见光波段设置多个大气活性氮观测通道，可捕捉不同氮组分的特征吸收信号，为特征参量提取提供波段基础。
针对二氧化氮提取对流层柱浓度差分斜柱量、大气质量校正因子，针对氨气提取210-225nm波段连续吸收斜率、差分吸收光学厚度，针对颗粒态氮提取气溶胶光学厚度波段差值、气溶胶体积中值半径。提取后通过最大互信息系数剔除低关联参量，最终保留12个核心特征参量，建立随机森林回归映射模型，对二氧化氮、氨气浓度的拟合决定系数分别达0.78、0.72，满足估算精度要求。
干沉降速率是表征大气含氮组分向下垫面沉积效率的核心参数，受大气边界层特征、下垫面属性与组分性质共同调控，传统参数化方案采用固定阻抗默认值，忽略动态变化与精细差异，易引发系统性测算偏差。
本优化对干沉降速率公式
中的三类阻力分别修正：引入逐日大气稳定度参数实现逐栅格动态计算，结合组分扩散特性修正，依托高分辨率下垫面分类与动态遥感参数分季节修正，优化后测算偏差降低约16%。
本研究采用未参与模型训练与特征筛选的地面独立观测数据验证干沉降估算结果，避免验证结果虚高。选取研究区12个覆盖2018-2022年、涵盖三类典型下垫面的观测站点作为验证样本。空间维度分三类排放区域验证，所有分层相对偏差均控制在25%以内，校准后高海拔区域偏差降低7.2个百分点；时间维度分季节、月度验证，整体决定系数达0.71，校准后冬季氨气估算偏差降低8.1个百分点，精度满足研究需求。
降水是大气氮湿沉降的核心驱动因子，其反演精度直接决定湿沉降估算可靠性，传统地面站点仅能获取离散单点数据，难以覆盖地形复杂、站点稀疏区域的降水空间分异。本研究整合GPM微波观测与FY-3D多光谱红外观测数据，引入地形参数修正复杂地形反演误差，填补薄云缺测信息空白，采用贝叶斯模型分配动态权重融合两类反演结果，最终生成研究区2018-2022年1km分辨率逐日降水栅格，整体平均相对误差控制在12.8%以内，满足估算需求。
大范围开展雨相氮浓度原位观测受限于成本与维护条件，现有观测多为单点离散分布，无法支撑全域连续湿沉降估算，雨相氮浓度的空间分布估算一直是湿沉降研究的核心瓶颈。依托卫星获取的大气柱氮要素遥感观测结果，结合活性氮迁移转化规律，可建立以柱浓度为核心输入的雨相氮浓度间接反演方法。
大气柱氮浓度与雨相氮浓度存在稳定统计关联，本研究以预先提取的各活性氮特征参量为输入，结合近地面氮浓度映射得到逐栅格氮组分含量，分别构建气态、颗粒态活性氮与雨相氮浓度的关联模型，引入降水强度、相对湿度、云光学厚度构建参数化方案，计算两类氮组分对雨相氮浓度的贡献。
累加得到总雨相氮浓度后，利用站点实测数据优化模型参数，最终模型拟合决定系数达0.69，可满足区域估算精度要求，填补了大范围雨相氮浓度观测空白。
本研究基于优化后的降水时空分布栅格与反演得到的雨相氮浓度栅格，完成时空匹配后搭建湿沉降估算模型：以月为统计单元，逐日降水累加得到月降水量，逐日雨相氮浓度按降水量加权得到月均值，采用1公里分辨率WGS-84坐标系UTM投影的逐像元对齐匹配，得到多时间尺度湿沉降通量结果。采用研究区17个独立站点2018-2022年实测数据校验，整体拟合决定系数达0.72，修正山地偏差后精度满足研究需求。
基于干沉降估算模型输出的2018-2022年逐月氮干沉降通量栅格数据，结合空间趋势分析识别研究区氮干沉降的空间分异特征，结果显示研究区多年平均氮干沉降通量为，整体呈现显著的梯度分异规律：高值区集中在中南部城镇密集区与西部河流沿岸农业平原，核心城区通量可达，北部低山丘陵生态保护区为低值区，通量普遍低于，核心排放源空间分布是分异的核心驱动。
时间维度上，受产业调整与减排政策影响，研究区氮干沉降通量五年累计降幅达13.8%，整体呈逐年下降趋势；季节波动呈现“夏高冬低”特征，夏季占全年通量的32.1%，冬季仅占18.7%。
基于2018-2022年逐月氮湿沉降通量栅格数据，结合降水格局分异特征分析可得：研究区氮湿沉降多年平均通量为，占总沉降约40%，空间呈现“南部高、北部西部低”的分异特征，高值集中在南部山前迎风坡与东南部沿湖水网区，降水空间差异是该分异格局的核心驱动因子。时间上整体呈小幅波动下降趋势，季节分异呈“夏高冬低”特征，丰水年通量较枯水年高出48.7%，降水主导湿沉降时序变化。
2018-2022年研究区大气氮总沉降多年平均值为，干沉降占比约60.0%，是区域氮沉降的主导贡献类型，与研究区降水天数少、近地面氮浓度高的特征匹配。
活性氮总沉降中还原态氮占比约52.6%，略高于氧化态氮，和区域近半数面积为农业区、氨排放占比超五成直接相关。不同功能区组分结构差异显著：城镇区氧化态氮占比达56.8%，农业区还原态氮占比升至58.7%，生态涵养区两类占比接近。氧化态氮占比五年间从45.1%升至49.3%，需强化农业氨排放管控。
传统区域氮沉降生态效应评估，大多依赖地面离散监测站点的人工原位观测获取数据，这类观测方式受站点布设成本、覆盖范围、监测人力等多重条件限制，仅能获取零散点位尺度的氮输入水平，无法覆盖山地、偏远林区、连片农田等人类活动难以长期抵达的区域，更难以从全域层面清晰刻画氮沉降通量的空间梯度差异，也无法精准捕捉不同生态系统氮饱和临界阈值的空间分异特征。当研究人员需要开展大尺度区域氮沉降生态效应评估时，只能依靠点位数据插值推算全域氮输入水平，插值过程中不可避免会引入系统误差，最终得到的评估结果往往与实际情况存在明显偏差，难以满足精细化生态环境管理的需求。
本研究依托构建的卫星遥感估算模型，输出了覆盖研究区域的1千米分辨率连续网格化氮沉降通量产品，不仅同步包含干沉降、湿沉降两种沉降途径的总氮通量数据，还可以区分不同形态活性氮的沉降贡献，能够为森林生态系统氮饱和风险识别、农田生态系统氮素盈余核算、流域水生态系统富营养化来源解析等多类研究，提供全覆盖的核心输入数据支撑。这类连续网格化的沉降数据，可以精准匹配行政区划、自然流域、生态功能区等不同层级管理单元的实际氮输入水平，有效解决传统观测数据覆盖不足带来的评估偏差问题，从基础数据层面提升氮污染风险区划、区域氮素平衡核算、水体富营养化综合治理的准确性，为不同层级生态环境管理部门的氮污染防控政策落地，提供扎实可靠的数据支撑。
本研究输出的高分辨率氮沉降空间分布产品，可清晰识别不同层级管控单元的氮沉降总负荷、干湿组分结构与不同污染源贡献的定量差异，为区域大气氮污染减排的分区分类精准管控提供扎实的数据支撑，破解传统管控模式依赖区域平均沉降量制定统一减排目标，导致管控措施针对性不足、资源错配的突出问题。当前我国大气氮污染管控已从传统的总量控制导向，全面转向精细化分区管控的新阶段，精准识别高氮沉降负荷风险区域与优先管控单元，是提升整体管控效率、降低全社会减排成本的核心前提。依托本研究产出的1千米分辨率连续覆盖的氮沉降通量空间数据，可结合现有行政区划框架与各类生态管控单元边界，精准计算每个管控单元的平均氮沉降负荷、超过生态临界负荷的面积占比，量化不同单元的氮污染风险等级，为管控优先级的科学划分提供明确的量化依据。研究区城镇密集区氧化态氮沉降占比超过五成，且集中了全域约70%的氮氧化物人为排放，农业平原区还原态氮沉降占比接近六成，农业源氨排放对沉降的贡献超过60%，基于沉降组分结构的空间分异特征，可针对不同功能区域制定差异化的减排方向：城镇开发单元优先推进工业深度脱硝、移动源氮氧化物减排，农业生产单元优先推进化肥减量增效、规模化畜禽养殖氨排放管控，避免统一减排要求带来的管控错位与成本浪费。对于跨界传输影响显著的生态涵养保护区域，可基于氮沉降空间分布数据明确上游排放对下游保护区的贡献占比，为建立跨区域氮污染联防联控机制提供量化支撑，清晰划分不同行政区的减排责任。针对氮沉降负荷已超过生态临界阈值的水源涵养区、自然保护区，可依托空间分布数据精准划定优先管控范围，匹配针对性减排措施，保障管控措施落地在高风险区域，提升生态环境质量改善的整体效率。模型成果可直接与区域“三线一单”生态环境管控体系对接，将氮沉降负荷、组分特征纳入生态管控单元的核心污染属性，为差异化管控措施制定提供精准数据支撑，推动区域大气氮污染管控从粗放式总量调控向精细化分区分类管控转型。
本研究构建的大气氮干湿沉降遥感估算模型以我国东部典型城市群为研究区完成参数校准，不同气候带、下垫面类型区域的氮排放结构、气象背景、下垫面属性存在显著差异，需针对目标区域特征完成通用化适配调整，才能保障模型输出精度，形成可复制推广的标准化业务化应用方案。
模型适配调整的核心环节围绕核心输入参数的本地化校准展开，针对目标区域的下垫面属性差异，调整干沉降速率计算的冠层阻力参数区间：若目标区域为高纬度寒温带针叶林区，需结合针叶林常年恒定的叶面积指数修正冠层阻力默认值，避免落叶林季节波动参数适配错误；若目标区域为西北干旱半干旱荒漠草原区，需结合低植被覆盖度的裸露地表特征，重新标定城镇未利用地混合下垫面的表面阻抗系数，匹配稀疏植被覆盖下干沉降过程的实际特征；针对云贵高原等喀斯特复杂地形区域，需结合区域数字高程模型重新计算地形校正系数，优化降水反演的地形抬升与雨影效应修正参数，保障湿沉降估算的地形适配性。
针对目标区域的气象特征差异，调整气象要素的同化校准方案：若目标区域云量常年偏高，卫星遥感活性氮观测存在大量缺测，需引入再分析氮浓度产品补充缺测数据，同时优化特征参量映射模型的训练样本，提升云覆盖条件下浓度反演的稳定性；若目标区域为季风气候区，降水季节波动极强，需重新标定雨相氮浓度反演的降水强度依赖系数，匹配不同降水强度下氮组分冲刷比例的变化特征。
完成参数适配后，依托目标区域已有的地面观测站点完成模型二次校准，输出符合业务化应用要求的标准化沉降产品，形成包含数据预处理、模型参数设置、结果校准全流程的可复制推广方案，可直接供不同区域的氮沉降业务化监测工作使用，整体适配调整流程可在1个月内完成，具备较高的推广效率。
