基于TERRA和AQUA双星MODIS影像协同反演陆地上空气溶胶光学厚度的研究

基于TERRA和AQUA双星MODIS影像协同反演陆地上空气溶胶光学厚度的研究一、引言1.1研究背景与意义1.1.1研究背景气溶胶作为分散在大气中的小颗粒物质，涵盖了沙尘、烟尘、雾霾等多种类型，广泛存在于地球大气环境中。其对环境与气候有着多方面的重要影响，在环境领域，气溶胶是导致近年来城市及郊区频繁出现“灰霾”天气的重要原因之一。当气溶胶颗粒大量增加时，会显著降低大气能见度，干扰人们的正常出行和生活，如2013年我国中东部地区多次出现的严重雾霾天气，使得多个城市的能见度降至极低水平，对交通、旅游业等造成了巨大冲击。同时，气溶胶中的可吸入颗粒物（PM10）和可入肺颗粒物（PM2.5）对人体健康危害巨大。这些细小颗粒能够随着呼吸进入人体呼吸道和肺部，引发呼吸道疾病、心血管问题等，长期暴露在高浓度气溶胶环境中，还会增加患肺癌等严重疾病的风险。从气候角度来看，气溶胶对太阳辐射有着吸收和散射作用，进而影响地球的能量收支平衡。一方面，部分气溶胶能够散射太阳辐射，将部分太阳辐射反射回宇宙空间，使得到达地球表面的太阳辐射减少，从而产生冷却效应；另一方面，一些具有吸光性的气溶胶，如黑碳等，能够吸收太阳辐射，增加大气的能量，起到加热作用。气溶胶还在云的形成和性质方面扮演关键角色。它可以作为云凝结核，影响云的微物理过程，改变云的光学特性、云量和降水效率。例如，在某些地区，气溶胶浓度的增加会导致云滴数量增多、尺寸减小，从而使云层的反照率增加，进一步影响地球的辐射平衡和气候。光学厚度是衡量气溶胶浓度的关键参数，通常用AOD（AerosolOpticalDepth）来表示，它反映了气溶胶对太阳辐射的吸收和散射作用所产生的光学效应。AOD的空间分布和变化蕴含着丰富的大气环境信息，对环境监测和气候研究具有重要意义。准确获取陆地上空气溶胶光学厚度，能够帮助我们更精确地了解大气污染的程度和范围，为制定有效的污染治理措施提供科学依据。在全球气候变化研究中，气溶胶光学厚度的变化是评估气候变化的重要指标之一，有助于深入理解气候变化的机制和趋势。MODIS（ModerateResolutionImagingSpectroradiometer）是美国NASA卫星上的一种先进的多光谱遥感探测仪器，它能够获取地表和大气的高分辨率遥感影像，覆盖了陆地、海洋和云层等多种目标。TERRA和AQUA是两颗搭载着MODIS仪器的卫星，它们顺序飞行于地球轨道上，并且可以实现同步或交替观测同一地区。这种观测模式使得它们能够提供更加全面和准确的地球观测数据，为陆地上空气溶胶光学厚度的反演研究提供了丰富的数据来源。通过对TERRA和AQUA双星MODIS影像的协同分析，可以充分利用两颗卫星的观测优势，提高气溶胶光学厚度反演的精度和可靠性，为深入研究气溶胶的时空分布特征和变化规律奠定坚实基础。1.1.2研究意义本研究利用TERRA和AQUA双星MODIS影像协同反演陆地上空气溶胶光学厚度，具有多方面的重要意义。在大气环境研究领域，准确的气溶胶光学厚度数据能够帮助我们更全面、深入地了解大气污染的时空分布特征。通过分析不同地区、不同时间的气溶胶光学厚度，我们可以清晰地掌握大气污染的严重程度和变化趋势，识别出污染高发区域和时段。这为进一步研究大气污染的形成机制提供了关键数据支持，有助于揭示污染源的排放特征、传输路径以及气象条件等因素对污染扩散的影响，从而为制定针对性的大气污染治理策略提供科学依据。在气候预测方面，气溶胶作为影响地球气候的重要因子，其光学厚度的精确测量对于提高气候模型的准确性至关重要。气候模型中对气溶胶的描述和参数化是影响模型模拟和预测能力的关键环节，而准确的气溶胶光学厚度数据能够为模型提供更真实的输入参数。通过本研究获取的高精度气溶胶光学厚度数据，可以改善气候模型对气溶胶辐射强迫的模拟，进而提高气候预测的精度，使我们能够更准确地预测未来气候变化趋势，为应对气候变化的决策制定提供有力支持。从环境监测角度来看，卫星遥感技术具有时效性高、覆盖面广、分辨率高等优势，能够实现对陆地大面积区域的气溶胶光学厚度的快速反演。利用TERRA和AQUA双星MODIS影像协同反演气溶胶光学厚度，可以进一步提高气溶胶数据的时空分辨率，并降低误差。这使得我们能够实时、动态地监测大气环境质量的变化，及时发现大气污染事件，为环境保护部门的环境监管和应急响应提供及时、准确的信息。研究气溶胶的季节变化和区域分布特征，还可以为治理雾霾、改善空气质量提供科学依据和决策支持，有助于制定合理的环境政策和措施，促进区域可持续发展。1.2国内外研究现状1.2.1国外研究进展国外在利用双星MODIS影像反演气溶胶光学厚度方面开展了大量的研究工作，取得了一系列重要成果。早在21世纪初，随着MODIS传感器的广泛应用，研究人员就开始探索利用其数据进行气溶胶光学厚度的反演。Kaufman等科学家通过大量的飞机试验，发现对于植被密集且具有较低反照率的地表，2.13μm近红外通道反照率与0.47μm、0.66μm可见光通道反射率之间存在良好的相关性，并成功将这一结论应用于MODIS的气溶胶反演算法中。这一开创性的工作为后续基于MODIS数据的气溶胶反演研究奠定了基础。随着研究的深入，许多学者致力于改进反演算法以提高气溶胶光学厚度的反演精度。例如，Levy等对MODIS气溶胶反演算法进行了优化，提出了改进的查找表（LUT）方法，通过更准确地考虑地表反射率和气溶胶模型等因素，有效降低了反演误差。此外，一些研究还关注不同地区气溶胶特性的差异对反演结果的影响。在对沙漠地区的研究中，发现该地区气溶胶以沙尘为主，其光学特性与其他地区有所不同，因此需要针对沙尘气溶胶的特点对反演算法进行调整和优化，以获得更准确的气溶胶光学厚度数据。近年来，随着多源数据融合技术的发展，国外研究人员开始尝试将MODIS数据与其他卫星数据或地面观测数据相结合，以进一步提高气溶胶反演的精度和可靠性。如将MODIS数据与MISR（Multi-angleImagingSpectroRadiometer）数据进行融合，利用MISR的多角度观测优势，获取更丰富的气溶胶信息，从而改善气溶胶光学厚度的反演效果。还有研究将MODIS数据与地面太阳光度计观测数据相结合，通过对比和验证，对反演算法进行校准和改进，使得反演结果与地面观测数据具有更好的一致性。在应用研究方面，国外学者利用双星MODIS影像反演得到的气溶胶光学厚度数据，开展了广泛的研究。在气候变化研究领域，通过分析长时间序列的气溶胶光学厚度数据，研究气溶胶对全球辐射平衡和气候变暖的影响机制；在空气质量监测方面，利用高分辨率的气溶胶光学厚度数据，实时监测城市和区域的空气质量状况，为空气污染预警和治理提供科学依据。1.2.2国内研究进展国内在双星MODIS影像反演气溶胶光学厚度领域也取得了显著的研究成果。在早期阶段，国内学者主要是对国外已有的反演算法进行引进和应用，并结合我国的实际情况进行适应性改进。毛节泰等利用北京大学地面多波段太阳光度计的观测数据，对MODIS卫星遥感气溶胶的方法进行了对比研究，结果表明二者具有较好的相关性，为国内开展相关研究提供了重要的参考依据。随着研究的不断深入，国内学者在反演算法的改进和创新方面取得了一系列突破。张维刚、李龙剑等针对MODIS数据特点，提出了一种基于遗传算法优化神经网络的气溶胶光学厚度反演方法。该方法通过遗传算法对神经网络的权值和阈值进行优化，提高了神经网络的收敛速度和反演精度，有效改善了传统反演算法在复杂地表条件下的反演效果。此外，一些研究还关注到气溶胶的时空分布特征与区域环境因素之间的关系。通过对我国不同地区气溶胶光学厚度的长期监测和分析，发现气溶胶的分布受到地形、气象条件以及人类活动等多种因素的综合影响。在京津冀地区，由于工业排放和机动车尾气排放量大，加上特殊的地形条件不利于污染物的扩散，导致该地区气溶胶光学厚度较高，大气污染问题较为严重。在应用方面，国内研究人员将双星MODIS影像反演的气溶胶光学厚度数据广泛应用于环境监测、气候研究和城市规划等领域。在环境监测中，利用气溶胶光学厚度数据对大气污染进行实时监测和预警，为环保部门制定污染治理措施提供科学依据；在气候研究中，分析气溶胶对区域气候的影响，为气候变化的预测和应对提供支持；在城市规划中，结合气溶胶分布特征，合理规划城市布局和产业发展，减少大气污染对城市居民的影响。1.3研究目标与内容1.3.1研究目标本研究旨在充分利用TERRA和AQUA双星MODIS影像的优势，通过协同反演的方法，实现对陆地上空气溶胶光学厚度的高精度反演。具体目标包括：构建适用于双星MODIS影像的气溶胶光学厚度协同反演算法，提高反演精度，降低反演误差；获取长时间序列、高时空分辨率的陆地上空气溶胶光学厚度数据，为研究气溶胶的时空分布特征和变化规律提供数据支持；分析气溶胶光学厚度与气象因素、地表特征等环境要素之间的相互关系，揭示气溶胶的形成、传输和演化机制。1.3.2研究内容本研究的主要内容涵盖数据获取与预处理、气溶胶光学厚度反演算法构建以及反演结果分析与验证三个关键方面。在数据获取与预处理阶段，从NASA官网下载TERRA和AQUA卫星的MODISL1B原始影像数据，这些数据涵盖了多个波段，包含丰富的地表和大气信息。下载时需筛选出研究区域内的影像，并确保数据的完整性和准确性。利用ENVI等专业遥感图像处理软件，对下载的MODIS影像进行辐射校正，将原始的数字量化值（DN）转换为具有物理意义的辐射亮度值，消除传感器自身特性和大气辐射传输过程对辐射量测量的影响，使不同时间、不同条件下获取的影像数据具有可比性。进行大气校正，去除大气分子、气溶胶、水汽等对太阳辐射的散射和吸收作用，获取真实的地表反射率信息。同时，对影像进行几何校正，消除因卫星轨道、姿态、地球曲率等因素导致的影像几何变形，确保影像中地物的位置精度。此外，还需进行云检测和云掩膜处理，识别并去除影像中的云覆盖区域，避免云对气溶胶光学厚度反演的干扰。气溶胶光学厚度反演算法构建是研究的核心内容。基于MODIS气溶胶反演产品MOD04（TERRA卫星）和MYD04（AQUA卫星）的算法原理，深入研究暗像元法、深蓝算法等经典反演算法在双星MODIS影像中的适用性。暗像元法利用植被、水体等低反照率地表的特性，通过查找表法建立卫星观测反射率与气溶胶光学厚度之间的关系；深蓝算法则针对高反照率地表，如沙漠、干旱地区等，通过特定的波段组合和假设条件来反演气溶胶光学厚度。针对研究区域的特点，如地形复杂、地表类型多样等，对现有算法进行改进和优化。考虑不同地表类型的反射率特性差异，采用更准确的地表反射率模型；结合研究区域的气象数据，如相对湿度、风速等，对气溶胶模型进行调整，以提高反演算法对复杂环境的适应性。利用地面太阳光度计观测数据，如AERONET（AErosolROboticNETwork）站点数据，对改进后的反演算法进行验证和校准。通过对比卫星反演结果与地面观测数据，评估反演算法的精度和可靠性，进一步优化算法参数，提高反演精度。在反演结果分析与验证方面，利用反演得到的陆地上空气溶胶光学厚度数据，采用空间分析方法，如克里金插值、泰森多边形等，分析气溶胶光学厚度的空间分布特征，绘制气溶胶光学厚度空间分布图，揭示其在不同地形、地貌和土地利用类型下的分布规律。运用时间序列分析方法，分析气溶胶光学厚度的时间变化趋势，包括日变化、季节变化和年际变化等，探讨其变化的驱动因素，如气象条件的季节性变化、人类活动的周期性等。结合气象数据，如温度、湿度、气压、风速等，以及地表特征数据，如土地利用类型、植被覆盖度等，通过相关性分析、多元线性回归等统计方法，研究气溶胶光学厚度与这些环境要素之间的相互关系，揭示气溶胶的形成、传输和演化机制。将反演结果与其他卫星反演产品，如VIIRS（VisibleInfraredImagingRadiometerSuite）气溶胶产品，以及地面监测站点的实测数据进行对比分析，评估本研究反演结果的准确性和可靠性，进一步验证反演算法的有效性。二、TERRA和AQUA双星MODIS影像概述2.1TERRA和AQUA卫星介绍2.1.1TERRA卫星特性TERRA卫星是美国国家航空航天局（NASA）地球观测系统（EOS）计划中的第一颗卫星，于1999年12月18日成功发射。其主要任务是实现从单系列极轨空间平台上对太阳辐射、大气、海洋和陆地进行综合观测，获取有关海洋、陆地、冰雪圈和太阳动力系统等信息，进而实现对大气和地球环境变化的长期观测和研究的总体战略目标。TERRA卫星运行于太阳同步近极地轨道，轨道高度约为705千米，轨道倾角为98.2°，重复周期为16天。这样的轨道设计使得TERRA卫星能够对地球表面进行全面且系统的观测，每天可获得4条过境轨道资料，能有效覆盖地球大部分区域。卫星搭载了多种先进的探测仪器，其中中分辨率成像光谱仪（MODIS）是其核心仪器之一。MODIS具有36个中等分辨率水平（0.25μm-1μm）的光谱波段，覆盖了从可见光到热红外的广泛光谱范围，光谱分辨率高，能够提供丰富的地表和大气信息。其空间分辨率分别为250米（波段1-2）、500米（波段3-7）和1000米（波段8-36）。不同分辨率的波段适用于不同的研究和应用需求，例如250米分辨率的波段可用于对城市区域、农田等较小尺度地物的精细观测；500米分辨率波段有助于对较大范围的植被、水体等进行监测；1000米分辨率波段则可用于宏观的全球尺度的研究，如气候变化、生态系统监测等。MODIS每1-2天可对地球表面观测一次，这使得它能够及时捕捉到地球表面的动态变化，为实时监测和研究提供了有力支持。除MODIS外，TERRA卫星还搭载了云与地球辐射能量系统（CERES）、多角度成像光谱仪（MISR）、先进星载热辐射与反射辐射计（ASTER）和对流层污染测量仪（MOPITT）等仪器。这些仪器各自具有独特的观测能力和科学目标，与MODIS相互配合，实现了对地球系统多方面的综合观测。CERES主要用于测量地球的辐射平衡，研究云层和辐射的相互作用，为理解地球的能量收支和气候变化提供关键数据；MISR通过多角度观测，获取地表和云层的反射特性，有助于提高对气溶胶、云特性等的反演精度；ASTER用于高分辨率的地形和地表温度测量，其14个光谱波段可提供丰富的地表信息，在地质勘探、生态环境研究等领域发挥重要作用；MOPITT则专注于对流层中一氧化碳等污染物的测量，对于研究大气污染传输和化学过程具有重要意义。2.1.2AQUA卫星特性AQUA卫星同样是NASA地球观测系统计划中的重要成员，于2002年5月4日成功发射。它的主要使命是研究地球的水循环和气候系统，通过对大气、海洋和陆地表面的观测，获取与水相关的各种信息，如海洋表面温度、海冰、降水和土壤湿度等。AQUA卫星与TERRA卫星在轨道和观测模式上既有相似之处，也存在一些差异。AQUA卫星也运行在太阳同步近极地轨道，轨道高度同样约为705千米，这使得它与TERRA卫星能够在相近的高度对地球进行观测。但AQUA卫星的过境时间与TERRA不同，它是下午星，从南向北于地方时13：30左右通过赤道。这种时间上的差异使得两颗卫星能够在不同时段对同一地区进行观测，从而获取更全面的地球表面信息。例如，对于一些在一天中变化明显的现象，如大气边界层的演变、地表能量平衡的日变化等，通过TERRA和AQUA卫星的不同时段观测，可以更完整地了解其变化过程。在仪器搭载方面，AQUA卫星也装载了MODIS，其MODIS仪器与TERRA卫星上的MODIS具有相同的技术参数和观测能力，同样拥有36个光谱波段，空间分辨率也分为250米、500米和1000米。这使得两颗卫星的MODIS数据具有可比性和互补性。AQUA卫星还搭载了先进微波扫描辐射计-EOS（AMSR-E）等仪器。AMSR-E使用微波波段进行观测，主要用于测量海洋表面温度、海冰、降水和土壤湿度等参数。微波具有较强的穿透能力，能够穿透云层和一定深度的地表，获取其他波段难以探测到的信息。在测量土壤湿度时，微波可以探测到地表以下一定深度的水分含量，为农业、水资源管理等领域提供重要的数据支持。2.2MODIS传感器及其影像特点2.2.1MODIS传感器参数MODIS传感器具有独特且丰富的参数设置，为地球观测提供了全面而精细的数据。在波段设置方面，它拥有36个离散的光谱波段，光谱范围从0.4微米的可见光波段一直延伸到14.4微米的热红外波段，实现了全光谱覆盖。这些波段被精心设计用于探测不同的地球物理和生物地球化学过程。其中，波段1（0.62-0.67μm）和波段2（0.841-0.876μm）属于可见光和近红外波段，对植被的反射特性非常敏感，常用于计算归一化植被指数（NDVI），以评估植被的生长状况和健康程度。在植被生长茂盛的区域，近红外波段的反射率较高，而可见光波段的反射率相对较低，通过NDVI的计算可以清晰地反映出植被的分布和生长态势。波段3（0.459-0.479μm）、波段4（0.545-0.565μm）和波段5（1.23-1.25μm）等波段在监测大气气溶胶、云特性以及水体颜色等方面发挥着重要作用。波段3和波段4对大气中的气溶胶粒子的散射和吸收特性敏感，可用于反演气溶胶光学厚度，评估大气污染程度；波段5则有助于区分不同类型的水体，如识别海洋中的浮游植物分布和水体的浑浊程度。在分辨率方面，MODIS具有三种不同的空间分辨率，以满足不同尺度的研究需求。其中，波段1和波段2的空间分辨率为250米，这使得它能够对较小尺度的地物进行精细观测，如城市中的建筑物、农田的边界以及小型水体等。在城市研究中，可以利用这两个波段的高分辨率影像，准确地识别城市的建筑类型、道路网络和绿地分布，为城市规划和管理提供详细的数据支持。波段3-7的分辨率为500米，适用于对中等尺度区域的监测，如大面积的森林、草原和湖泊等。通过这些波段的影像，可以对森林的覆盖范围、草原的植被变化以及湖泊的面积和水质变化进行有效的监测。波段8-36的分辨率为1000米，适合用于宏观尺度的全球观测，如研究全球气候变化、生态系统的分布和演化等。在全球气候变化研究中，利用1000米分辨率的波段数据，可以分析全球植被覆盖的变化趋势、冰川的退缩和扩张以及海洋表面温度的变化等。MODIS的扫描宽度达到2330公里，这意味着它在一次观测中能够覆盖非常广阔的区域，大大提高了数据获取的效率。结合其每天1-2次对地球表面的观测频率，使得MODIS能够及时捕捉到地球表面的动态变化，为实时监测和研究提供了有力支持。2.2.2MODIS影像优势MODIS影像在气溶胶监测领域展现出多方面的显著优势。首先，其具有较高的时间分辨率，每天可对地球表面观测1-2次。这使得它能够实时追踪气溶胶的动态变化，及时捕捉到气溶胶浓度的突然升高或降低等异常情况。在工业排放集中的地区，可能会出现突发的污染物排放事件，导致气溶胶浓度瞬间增加。MODIS影像能够快速监测到这种变化，为环境监管部门提供及时的预警信息，以便采取相应的措施进行应对。相比之下，一些其他卫星传感器的观测频率较低，可能无法及时发现这些短期的气溶胶变化，从而影响对大气污染事件的及时处理。MODIS影像拥有丰富的光谱信息，其36个波段覆盖了从可见光到热红外的广泛光谱范围。不同类型的气溶胶在不同波段上具有独特的光谱特征，通过对这些光谱信息的分析，可以实现对气溶胶类型的准确识别。沙尘气溶胶在近红外波段具有较高的反射率，而有机气溶胶和黑碳气溶胶在可见光和近红外波段的吸收特性则有所不同。利用MODIS影像的多波段数据，通过建立合适的反演模型，可以有效地识别出不同类型的气溶胶，进而深入研究它们的来源、传输和演化规律。这对于理解大气污染的形成机制和制定针对性的污染治理策略具有重要意义。在空间覆盖范围上，MODIS影像也具有明显优势。其扫描宽度达2330公里，能够在一次观测中覆盖大面积的区域。这使得它适用于对大尺度区域的气溶胶分布进行监测，如整个大洲或大洋的气溶胶分布情况。在研究跨区域的大气污染传输时，MODIS影像可以提供全面的气溶胶分布信息，帮助研究人员了解气溶胶在不同地区之间的传输路径和扩散范围。相比高分辨率但小范围的影像，MODIS影像能从宏观角度展示气溶胶的整体分布格局，为研究区域乃至全球尺度的气溶胶问题提供了宏观视角。MODIS数据还具有全球免费获取的优势。这使得全球范围内的科研人员、政府机构和环保组织等都能够方便地获取和使用这些数据。无论是发达国家还是发展中国家，都可以利用MODIS数据开展气溶胶监测和研究工作，促进了全球范围内对大气污染问题的共同关注和研究。免费的数据获取政策降低了研究成本，提高了数据的可及性和共享性，有利于推动气溶胶监测技术的发展和应用，促进全球大气环境的保护和治理。2.3在气溶胶反演中的应用基础2.3.1气溶胶光学特性探测原理MODIS影像探测气溶胶光学特性主要基于其对太阳辐射的散射和吸收原理。当太阳辐射进入地球大气层后，会与气溶胶粒子发生相互作用。气溶胶粒子的大小、形状和化学成分等因素决定了其对不同波长太阳辐射的散射和吸收能力。MODIS传感器通过探测不同波段的太阳辐射在经过气溶胶层后的反射、散射和透射情况，获取气溶胶的光学信息。MODIS拥有36个光谱波段，这些波段覆盖了从可见光到热红外的广泛光谱范围，为探测气溶胶光学特性提供了丰富的信息。在可见光波段（如0.47μm、0.66μm等波段），气溶胶对太阳辐射的散射作用较为明显，且不同类型的气溶胶在这些波段的散射特性存在差异。沙尘气溶胶由于其粒子较大，在可见光波段的散射能力较强，且散射光的强度和方向与粒子的大小、形状密切相关；而黑碳气溶胶等吸光性气溶胶则在这些波段表现出较强的吸收特性，导致反射光的强度降低。通过分析MODIS在这些波段接收到的反射辐射，可以反演气溶胶的光学厚度和粒子大小等参数。例如，利用MODIS在0.47μm和0.66μm波段的反射率数据，结合辐射传输模型，可以建立气溶胶光学厚度与反射率之间的关系，从而反演出气溶胶光学厚度。在近红外波段（如0.86μm、1.24μm等波段），气溶胶的散射和吸收特性也与可见光波段有所不同。一些气溶胶在近红外波段的散射能力相对较弱，但对某些成分的气溶胶，如有机气溶胶，其在近红外波段的吸收特性可能更为突出。MODIS在近红外波段的观测数据可以提供额外的信息，用于进一步区分气溶胶的类型和特性。结合可见光和近红外波段的数据，可以更准确地识别气溶胶的成分和来源。通过对比不同波段的反射率和吸收特性，可以判断气溶胶中是否含有沙尘、黑碳、有机物质等成分，进而了解气溶胶的来源是自然源（如沙尘排放）还是人为源（如工业排放、生物质燃烧）。在热红外波段（如10.7μm、11.0μm等波段），气溶胶对地球表面发射的热辐射也有一定的吸收和散射作用。虽然热红外波段主要用于探测地表温度和云特性等，但气溶胶在该波段的影响也不容忽视。在一些污染严重的地区，气溶胶的存在会改变热红外辐射的传输路径和强度，从而影响对地表温度的准确测量。通过分析MODIS在热红外波段的观测数据，可以间接获取气溶胶对热辐射的影响信息，为研究气溶胶与地表能量平衡的相互关系提供依据。2.3.2双星协同观测的优势TERRA和AQUA双星协同观测在气溶胶光学厚度反演中具有多方面的显著优势，能够有效提高反演的精度和可靠性。双星协同观测大大提高了时间分辨率。TERRA卫星是上午星，从北向南于地方时10：30左右通过赤道；AQUA卫星是下午星，从南向北于地方时13：30左右通过赤道。这种不同的过境时间使得两颗卫星能够在一天内的不同时段对同一地区进行观测。在一天中，气溶胶的浓度和分布可能会发生变化，例如在城市地区，早晨由于人类活动的增加，如机动车出行、工业生产启动等，气溶胶浓度可能会逐渐升高；而在下午，随着大气扩散条件的变化和部分污染物的沉降，气溶胶浓度可能会有所降低。通过TERRA和AQUA卫星的不同时段观测，可以更全面地捕捉到气溶胶在一天内的变化情况，从而获取更准确的气溶胶光学厚度时间序列数据。相比单颗卫星观测，双星协同观测能够提供更丰富的时间维度信息，有助于研究气溶胶的短期变化规律和动态过程。双星观测还可以增加数据的空间覆盖范围。尽管TERRA和AQUA卫星的轨道相似，但它们在观测过程中会存在一定的观测角度差异。这种角度差异使得两颗卫星能够从不同视角获取同一地区的影像，从而覆盖更广泛的区域。在山区等地形复杂的地区，由于地形起伏和阴影的影响，单颗卫星可能无法完全观测到所有区域。但双星协同观测时，不同的观测角度可以弥补这一不足，使得更多的区域能够被观测到，减少观测盲区。通过对双星影像的融合处理，可以获得更完整的气溶胶光学厚度空间分布信息，提高空间分辨率和精度。双星协同观测还可以提高数据的可靠性和准确性。由于卫星观测过程中可能会受到各种因素的干扰，如云层覆盖、大气噪声等，导致观测数据存在误差。TERRA和AQUA卫星对同一地区的观测数据可以相互验证和补充。当一颗卫星的观测数据受到云层影响时，另一颗卫星在不同时间的观测数据可能是晴空条件下的，从而可以利用这颗卫星的数据进行反演或对受影响的数据进行校正。通过对比两颗卫星的观测数据，还可以评估数据的质量和可靠性，识别出异常数据并进行剔除或修正。利用双星观测数据进行交叉验证，可以有效降低反演误差，提高气溶胶光学厚度反演的准确性。三、陆地上空气溶胶光学厚度反演原理与方法3.1反演基本原理3.1.1辐射传输理论基础辐射传输理论是研究电磁辐射在介质中传播过程的理论，它在气溶胶反演中起着至关重要的基础作用。在地球大气环境中，太阳辐射从太阳表面发出，穿过大气层到达地球表面，在这个过程中，辐射与大气中的各种成分，包括气溶胶、气体分子、云等发生相互作用。辐射传输理论的核心是辐射传输方程（RTE，RadiativeTransferEquation），它描述了辐射强度在空间和方向上的变化规律。对于平面平行大气模型，辐射传输方程可以表示为：\frac{dI(\tau,\mu)}{d\tau}=-I(\tau,\mu)+\frac{\omega_0}{4\pi}\int_{-1}^{1}\int_{0}^{2\pi}p(\mu,\mu',\varphi-\varphi')I(\tau,\mu')d\varphi'd\mu'+\frac{S(\tau,\mu)}{4\pi}其中，I(\tau,\mu)是在光学厚度为\tau处，沿天顶角余弦为\mu方向的辐射强度；\omega_0是单次散射反照率，表示气溶胶粒子散射辐射与吸收辐射的相对比例，\omega_0的值介于0（完全吸收）和1（完全散射）之间，不同类型的气溶胶具有不同的单次散射反照率，例如，黑碳气溶胶的单次散射反照率较低，表明其吸收能力较强，而硫酸盐气溶胶的单次散射反照率相对较高，散射作用更为显著；p(\mu,\mu',\varphi-\varphi')是散射相函数，用于描述气溶胶粒子散射辐射的方向分布特性，它与气溶胶粒子的大小、形状和折射率等因素密切相关，不同形状的气溶胶粒子，如球形、非球形等，其散射相函数会有明显差异，通过研究散射相函数，可以了解气溶胶对辐射的散射方向偏好，从而更好地理解辐射在大气中的传输过程；S(\tau,\mu)是源函数，包括直接太阳辐射、大气散射辐射和地表反射辐射等，它综合考虑了辐射的各种来源，在不同的大气条件和地表特征下，源函数的组成和大小会发生变化，例如，在晴朗天气下，直接太阳辐射在源函数中占主导地位，而在多云天气或高反射率地表区域，大气散射辐射和地表反射辐射的贡献可能会增加。在气溶胶反演中，辐射传输方程用于建立卫星观测到的辐射信号与气溶胶光学厚度之间的定量关系。通过测量卫星接收到的不同波段的辐射强度，并结合已知的大气和地表参数，如大气成分、地表反射率等，求解辐射传输方程，可以反演出气溶胶的光学厚度等参数。由于辐射传输方程的复杂性，通常需要采用数值方法进行求解，如离散坐标法、累加法等。离散坐标法将辐射传输方程在空间和方向上进行离散化，通过求解离散化后的方程组得到辐射强度的数值解；累加法是一种逐步累加辐射传输过程中各种相互作用的方法，通过多次迭代计算，逐步逼近辐射传输方程的精确解。3.1.2气溶胶与辐射的相互作用气溶胶对太阳辐射的吸收和散射作用是其影响地球大气能量平衡和气候的关键机制，也是陆地上空气溶胶光学厚度反演的重要物理基础。气溶胶对太阳辐射的吸收作用主要取决于其化学成分和粒径大小。一些气溶胶成分，如黑碳、有机碳等，具有较强的吸光性。黑碳是一种由不完全燃烧产生的含碳颗粒物，它在可见光和近红外波段具有较高的吸收系数。当太阳辐射照射到含有黑碳的气溶胶上时，部分辐射能量被黑碳吸收，转化为热能，从而使气溶胶粒子温度升高，进而影响大气的热力学状态。气溶胶的粒径大小也会影响其吸收能力，一般来说，粒径较小的气溶胶粒子对短波辐射的吸收效率较高，而粒径较大的气溶胶粒子对长波辐射的吸收相对较强。在城市污染区域，由于存在大量的细颗粒物，其中包含一定量的黑碳和有机碳，这些细颗粒物对太阳辐射的吸收作用较为显著，会导致到达地面的太阳辐射减少，同时使大气温度升高，加剧城市热岛效应。气溶胶对太阳辐射的散射作用同样与气溶胶的粒径、形状和折射率等因素密切相关。当太阳辐射与气溶胶粒子相遇时，会发生散射现象，使辐射方向发生改变。根据气溶胶粒子粒径与辐射波长的相对大小，散射可分为瑞利散射、米氏散射和非球形粒子散射等类型。瑞利散射发生在气溶胶粒子粒径远小于辐射波长的情况下，其散射强度与波长的四次方成反比，因此对短波辐射的散射作用更为明显。在晴朗的天空中，由于大气中主要是气体分子和少量的细颗粒物，其粒径远小于可见光波长，瑞利散射使得天空呈现出蓝色，因为蓝光的波长较短，更容易被散射。米氏散射适用于气溶胶粒子粒径与辐射波长相当的情况，其散射特性较为复杂，散射强度与粒子的大小、形状、折射率以及辐射波长等因素都有关系。沙尘气溶胶粒子的粒径通常在微米量级，与可见光波长相当，当太阳辐射遇到沙尘气溶胶时，主要发生米氏散射，使得天空变得浑浊，颜色也会发生变化。对于非球形气溶胶粒子，如不规则形状的矿物颗粒、烟尘粒子等，其散射特性不能简单地用瑞利散射或米氏散射来描述，需要考虑粒子的形状因子等因素，采用更复杂的散射模型进行分析。非球形粒子的散射相函数具有独特的分布特征，会导致辐射在不同方向上的散射强度分布不均匀。气溶胶的吸收和散射作用对太阳辐射的传输和分布产生了重要影响。一方面，气溶胶的散射作用使得太阳辐射在大气中发生多次散射，部分辐射被散射回宇宙空间，部分辐射则被散射到其他方向，从而改变了辐射的传输路径和强度分布。这种散射作用会降低到达地面的直接太阳辐射强度，但增加了散射辐射的成分，使地面接收到的辐射更加均匀。在阴天或雾霾天气中，由于气溶胶浓度较高，散射作用强烈，直接太阳辐射大幅减弱，天空显得较为昏暗，但散射辐射使得周围环境仍有一定的亮度。另一方面，气溶胶的吸收作用则直接消耗了太阳辐射能量，减少了到达地面的总辐射量。吸收的辐射能量转化为热能，会加热气溶胶粒子和周围大气，影响大气的温度分布和垂直运动。在一些污染严重的地区，气溶胶的大量吸收会导致地面气温降低，而大气温度升高，形成逆温层，进一步抑制了污染物的扩散，加剧了大气污染。3.2反演算法3.2.1暗目标法暗目标法是一种在陆地上空气溶胶光学厚度反演中广泛应用的经典方法，其核心在于巧妙地利用特定地表类型的低反射率特性，从而有效扣除地表反射对卫星观测信号的影响，进而实现对气溶胶光学厚度的准确反演。该方法的关键假设是在植被茂密的区域以及水体等特定地表类型处，地表反射率相对较低，并且在近红外波段与可见光波段之间存在较为稳定的相关性。在植被茂密的地区，植被的叶绿素对近红外光具有强烈的反射作用，而对可见光的吸收相对较强，使得该区域在近红外波段的反射率较高，而在可见光波段的反射率较低。这种特性使得植被成为暗目标法中理想的地表类型选择。水体表面对太阳辐射的反射也具有一定的规律性，在某些波段的反射率相对较低，同样符合暗目标法的应用条件。在实际应用中，首先需要从MODIS影像中准确识别出这些暗目标区域。这通常借助于归一化植被指数（NDVI）等指标来实现。NDVI通过计算近红外波段与红光波段反射率的差值与和值之比，能够有效地突出植被的分布范围和生长状况。当NDVI值大于某一阈值时，可将该区域初步判定为植被茂密的暗目标区域。对于水体的识别，则可以利用水体在近红外波段的低反射率特性，通过设定合适的反射率阈值来筛选出水体区域。在识别出暗目标区域后，根据辐射传输理论，建立卫星观测反射率与气溶胶光学厚度之间的关系。由于暗目标区域的地表反射率相对稳定且较低，因此可以将其作为已知量，通过辐射传输方程来求解气溶胶光学厚度。辐射传输方程描述了太阳辐射在大气中传输过程中与气溶胶、大气分子等相互作用后到达卫星传感器的辐射强度变化。在考虑了大气的散射、吸收以及地表反射等因素后，通过对辐射传输方程的求解，可以得到气溶胶光学厚度与卫星观测反射率之间的定量关系。然而，暗目标法在实际应用中也存在一定的局限性。该方法对地表类型的要求较为苛刻，仅适用于植被茂密和水体等低反射率地表区域，对于沙漠、干旱地区等高反射率地表，暗目标法的假设不再成立，无法准确反演气溶胶光学厚度。大气条件的变化，如相对湿度、风速等，会影响气溶胶的光学特性和分布，从而对暗目标法的反演精度产生影响。在高湿度环境下，气溶胶粒子可能会吸湿增长，改变其粒径分布和光学特性，使得基于固定假设的暗目标法难以准确描述气溶胶与辐射的相互作用。3.2.2查找表法查找表法是陆地上空气溶胶光学厚度反演中的一种重要方法，其原理是通过预先建立的查找表，快速获取与卫星观测数据相对应的气溶胶光学厚度值，从而实现反演过程。查找表的构建是该方法的关键环节。通常利用辐射传输模型，结合大量的模拟计算来生成查找表。辐射传输模型能够精确地描述太阳辐射在大气中的传输过程，包括与气溶胶、大气分子、云层等的相互作用。在构建查找表时，需要考虑多种因素，如不同的气溶胶模型、地表反射率、太阳天顶角、卫星观测角等。对于气溶胶模型，常见的有大陆型、海洋型、城市型等，不同的气溶胶模型具有不同的光学特性和粒径分布，会对辐射传输产生不同的影响。地表反射率也是一个重要因素，不同的地表类型，如植被、水体、沙漠等，其反射率在不同波段存在显著差异，需要分别进行考虑。太阳天顶角和卫星观测角决定了辐射传输的路径和角度，对卫星接收到的辐射强度有重要影响。以6S（SecondSimulationoftheSatelliteSignalintheSolarSpectrum）辐射传输模型为例，在构建查找表时，首先设定一系列不同的气溶胶光学厚度值，从0到较大的污染浓度对应的光学厚度范围。针对每个气溶胶光学厚度值，结合不同的地表反射率、太阳天顶角和卫星观测角组合，利用6S模型计算出对应的卫星观测反射率。将这些计算结果整理成查找表，其中查找表的行和列分别对应不同的输入参数（如太阳天顶角、卫星观测角、地表反射率等），表格中的值则是相应条件下的卫星观测反射率和对应的气溶胶光学厚度。在实际反演过程中，当获取到MODIS影像的观测数据后，首先提取影像中的太阳天顶角、卫星观测角以及地表反射率等信息。然后，根据这些信息在预先构建的查找表中进行查找，找到与之最匹配的记录。该记录中对应的气溶胶光学厚度值即为反演得到的结果。在查找过程中，通常采用线性插值等方法来提高查找的精度，以获取更准确的气溶胶光学厚度值。如果查找表中没有与观测数据完全匹配的记录，可以通过对相邻记录进行线性插值，根据观测数据与相邻记录的参数差异，计算出相应的气溶胶光学厚度值。查找表法的优点在于计算速度快，能够快速实现气溶胶光学厚度的反演。一旦查找表构建完成，在反演过程中只需进行简单的查找和插值操作，无需进行复杂的辐射传输模型计算，大大提高了反演效率。查找表法也存在一定的局限性。查找表的构建依赖于辐射传输模型的准确性和输入参数的合理性，如果模型存在误差或输入参数与实际情况不符，会导致查找表的精度下降，进而影响反演结果的准确性。对于复杂的大气环境和多变的地表条件，查找表可能无法涵盖所有的情况，使得在某些特殊情况下反演结果的可靠性降低。三、陆地上空气溶胶光学厚度反演原理与方法3.3双星MODIS影像协同反演步骤3.3.1数据获取与预处理本研究所需的TERRA和AQUA卫星的MODIS影像数据均从NASA官网的LAADSDAAC（LandAtmosphereNearreal-timeCapabilityforEOS）数据平台获取。该平台提供了丰富的MODIS数据产品，其中L1B级数据是经过初步处理的原始影像数据，包含了传感器获取的原始辐射亮度信息，是进行气溶胶光学厚度反演的基础数据。在下载数据时，需根据研究区域的地理位置和时间范围进行筛选。研究区域若为中国华北地区，时间范围设定为2020年1月至2020年12月，在LAADSDAAC平台上设置相应的经纬度范围和时间筛选条件，确保下载的数据能够准确覆盖研究区域和时间段，且数据质量良好，无明显的数据缺失或异常值。获取数据后，利用ENVI（TheEnvironmentforVisualizingImages）软件对MODIS影像进行预处理，这是确保后续反演结果准确性的关键步骤。辐射校正旨在将MODIS影像的原始数字量化值（DN）转换为具有物理意义的辐射亮度值。在ENVI软件中，通过调用MODIS数据自带的定标参数文件，结合传感器的辐射响应函数，进行辐射定标操作。对于MODIS影像，其定标参数文件包含了增益系数、偏移量等关键信息，利用这些参数，通过公式L=G\timesDN+B（其中L为辐射亮度值，G为增益系数，DN为原始数字量化值，B为偏移量），可以将每个像元的DN值转换为辐射亮度值，消除传感器自身特性对辐射量测量的影响，使不同时间、不同条件下获取的影像数据具有可比性。大气校正则是为了去除大气分子、气溶胶、水汽等对太阳辐射的散射和吸收作用，获取真实的地表反射率信息。在ENVI软件中，选用6S（SecondSimulationoftheSatelliteSignalintheSolarSpectrum）辐射传输模型进行大气校正。该模型能够精确地模拟太阳辐射在大气中的传输过程，考虑了大气成分、气溶胶类型、太阳天顶角、卫星观测角等多种因素对辐射的影响。在进行大气校正时，需要输入大气模式、气溶胶模式、地表类型等参数。根据研究区域的地理位置和季节特点，选择合适的大气模式，如中纬度夏季模式（MidlatitudeSummer）用于夏季数据的处理，中纬度冬季模式（MidlatitudeWinter）用于冬季数据的处理；根据该区域常见的气溶胶类型，选择相应的气溶胶模式，如大陆型气溶胶模式（Continental）适用于污染较为严重的大陆地区；根据影像中地物的主要类型，确定地表类型参数。通过运行6S模型，对辐射亮度值进行校正，得到地表反射率影像。几何校正也是预处理的重要环节，其目的是消除因卫星轨道、姿态、地球曲率等因素导致的影像几何变形，确保影像中地物的位置精度。在ENVI软件中，利用地面控制点（GCP，GroundControlPoint）进行几何校正。通过在影像和参考地图或高分辨率影像上选取同名地物点作为GCP，如明显的道路交叉点、建筑物拐角等，一般选取不少于20个均匀分布的GCP。利用这些控制点，采用多项式变换等方法建立影像坐标与地理坐标之间的转换关系，对影像进行重采样和几何纠正，使影像的地理坐标与实际地理位置相符。云检测和云掩膜处理同样不可或缺，云的存在会严重干扰气溶胶光学厚度的反演，必须识别并去除影像中的云覆盖区域。在ENVI软件中，采用基于阈值法的云检测算法。该算法利用云在多个波段的反射率和亮度温度特征，设定相应的阈值来识别云像元。在可见光波段，云的反射率通常较高；在热红外波段，云的亮度温度较低。通过设定可见光波段反射率阈值大于0.4，热红外波段亮度温度阈值小于270K等条件，将满足条件的像元识别为云像元。对识别出的云像元进行掩膜处理，将其从影像中剔除，得到无云的影像数据，为后续的气溶胶光学厚度反演提供纯净的数据基础。3.3.2影像配准与融合双星影像的配准是协同反演的关键步骤，其目的是使TERRA和AQUA卫星获取的MODIS影像在空间位置上精确对齐，确保后续数据融合和分析的准确性。在进行影像配准时，首先在ENVI软件中打开TERRA和AQUA卫星的MODIS影像。选择影像中明显且稳定的地物特征作为同名点，如大型水体的边界、山脉的轮廓、城市的主要街道等。通过人工交互的方式，在两幅影像上精确标记同名点，一般选取30-50个同名点，以保证配准的精度。利用这些同名点，采用多项式变换模型进行影像配准。多项式变换模型能够通过数学函数描述两幅影像之间的几何变换关系，如一次多项式可描述平移和旋转，二次多项式可进一步描述缩放和扭曲等变形。在配准过程中，根据同名点的分布和数量，选择合适的多项式阶数，通过最小二乘法求解多项式的系数，实现影像的几何变换和配准。配准完成后，通过计算同名点的均方根误差（RMSE，RootMeanSquareError）来评估配准精度，一般要求RMSE小于半个像元，以确保配准后的影像能够准确对齐。影像融合是将配准后的TERRA和AQUA卫星MODIS影像进行合并，以充分利用两颗卫星的观测信息，提高数据的时空分辨率和信息量。本研究采用加权平均融合法进行影像融合。该方法基于TERRA和AQUA卫星的过境时间和观测条件，为两幅影像的每个像元分配不同的权重。由于TERRA卫星是上午星，AQUA卫星是下午星，它们在不同时段对同一地区进行观测，反映了不同时刻的大气和地表状况。在白天，太阳辐射较强，气溶胶的分布和特性可能随时间发生变化。根据研究区域的气候特点和观测经验，对于上午时段的TERRA卫星影像，赋予权重0.6；对于下午时段的AQUA卫星影像，赋予权重0.4。通过加权平均计算，得到融合后的影像，其像元值为I_{fused}=0.6\timesI_{TERRA}+0.4\timesI_{AQUA}（其中I_{fused}为融合后影像的像元值，I_{TERRA}为TERRA卫星影像的像元值，I_{AQUA}为AQUA卫星影像的像元值）。这样，融合后的影像既包含了上午时段气溶胶的信息，又融合了下午时段的观测数据，能够更全面地反映气溶胶的时空变化特征。加权平均融合法还能有效降低因单颗卫星观测误差或数据缺失带来的影响，提高数据的可靠性和稳定性。3.3.3气溶胶光学厚度计算在完成数据获取与预处理以及影像配准与融合后，利用改进的暗像元法结合查找表法来计算气溶胶光学厚度。在融合后的MODIS影像中，根据归一化植被指数（NDVI）和地表反射率特征，准确识别出暗像元区域。对于植被覆盖区域，当NDVI大于0.5且近红外波段（如2.1μm波段）的地表反射率小于0.1时，将该像元判定为暗像元。在水体区域，根据水体在可见光和近红外波段的低反射率特性，通过设定合适的反射率阈值，如可见光波段反射率小于0.05，近红外波段反射率小于0.1，筛选出水体暗像元。针对研究区域的地形和气候特点，对暗像元法中的查找表进行优化。在构建查找表时，充分考虑地形起伏对太阳辐射传输的影响，利用数字高程模型（DEM，DigitalElevationModel）数据，如SRTM（ShuttleRadarTopographyMission）数据，获取研究区域的地形高度信息。对于不同海拔高度的区域，调整查找表中的大气参数，如大气压强、水汽含量等，因为随着海拔升高，大气压强降低，水汽含量减少，这些因素会影响气溶胶与辐射的相互作用。根据研究区域的气候特征，如相对湿度、风速等，对气溶胶模型进行调整。在高湿度地区，气溶胶粒子可能会吸湿增长，改变其粒径分布和光学特性，因此在查找表中，针对不同的相对湿度范围，设置不同的气溶胶粒径分布和光学参数。在计算气溶胶光学厚度时，对于每个暗像元，提取其在MODIS影像中的多个波段反射率，如0.47μm、0.66μm、0.86μm等波段。根据暗像元所在的地理位置和时间，结合研究区域的气象数据，确定太阳天顶角、卫星观测角、大气模式等参数。利用这些参数，在优化后的查找表中进行查找和插值计算，得到该暗像元对应的气溶胶光学厚度值。采用双线性插值法，当查找表中没有与当前暗像元参数完全匹配的记录时，通过对相邻四个记录的参数和对应的气溶胶光学厚度值进行双线性插值，计算出该暗像元的气溶胶光学厚度。通过对影像中所有暗像元的计算，得到研究区域的气溶胶光学厚度分布影像。四、案例研究与数据分析4.1研究区域选择4.1.1区域概况本研究选取京津冀地区作为案例研究区域。京津冀地区位于华北平原北部，地理坐标大致在北纬36°05′-42°40′，东经113°27′-119°50′之间，总面积约21.6万平方千米。该地区是中国的政治、文化和国际交往中心，也是重要的经济区域，涵盖北京市、天津市以及河北省的石家庄、保定、秦皇岛等11个地级市。从地形地貌来看，京津冀地区地形复杂多样，西北部为山地和高原，属于燕山山脉和太行山脉的一部分，地势较高，海拔在1000米以上，这些山地在防风固沙、水源涵养方面发挥着重要作用。东南部为广阔的华北平原，地势平坦，自然条件优越，生产建设成本较低，是该地区经济社会活动的主要区域，集中了大量的城市、工业和农业用地。该地区还拥有丰富的海岸资源，秦皇岛市、天津市、唐山市、沧州市东临渤海，拥有大陆岸线长641千米，海岛岸线长199千米，天津港是京津冀城市群的海上门户，秦皇岛港是世界最大能源输出港之一，京唐港、黄骅港条件优越，曹妃甸是优良的深水港址，为对外开放和物流航运提供了重要保障。在气候方面，京津冀地区属于温带大陆性季风气候，四季分明。春季干旱多风，气温回升快，蒸发旺盛，降水稀少，平均降水量仅占全年的10%-15%，常出现沙尘天气，对大气环境质量产生较大影响。夏季高温多雨，降水集中在7-8月，平均降水量占全年的60%-70%，但降水分布不均，局部地区易出现暴雨洪涝灾害。高温高湿的天气条件有利于气溶胶的吸湿增长和化学反应，导致气溶胶浓度升高。秋季天高气爽，气温适中，降水逐渐减少，大气扩散条件相对较好，但在某些年份，由于冷空气活动较弱，也会出现静稳天气，不利于污染物的扩散。冬季寒冷干燥，受大陆冷气团控制，盛行西北风，平均降水量仅占全年的5%左右。冬季供暖期的到来，使得燃煤等能源消耗大幅增加，大量污染物排放到大气中，加上低温、高湿度和弱风等不利气象条件，容易形成雾霾天气，导致气溶胶浓度急剧上升。京津冀地区是中国经济发展的重要区域之一，也是人口密集区。近年来，随着城市化进程的加速和工业的快速发展，该地区的大气污染问题日益突出。大量的工业排放、机动车尾气排放、燃煤排放以及建筑扬尘等，使得该地区的气溶胶浓度一直处于较高水平。在工业排放方面，钢铁、化工、建材等行业是主要的污染源，这些行业排放的污染物中含有大量的二氧化硫、氮氧化物、颗粒物等，是气溶胶的重要组成部分。机动车保有量的快速增长也导致尾气排放成为大气污染的重要来源之一，汽车尾气中含有碳氢化合物、一氧化碳、氮氧化物和颗粒物等污染物，在阳光照射下，这些污染物会发生一系列复杂的光化学反应，形成二次气溶胶，进一步加重了大气污染。4.1.2选择理由选择京津冀地区作为研究区域，主要基于以下几个方面的考虑。京津冀地区是中国大气污染较为严重的区域之一，长期受到雾霾等大气污染问题的困扰，对该地区气溶胶光学厚度的研究具有重要的现实意义。该地区的大气污染问题不仅影响当地居民的身体健康，还对区域经济发展和生态环境造成了严重影响。准确获取该地区的气溶胶光学厚度，深入了解气溶胶的时空分布特征和变化规律，对于制定有效的大气污染治理措施，改善区域空气质量具有重要的指导作用。京津冀地区的经济活动活跃，城市化进程快速推进，人类活动对大气环境的影响显著。大量的工业企业、机动车和居民生活排放的污染物，使得该地区的气溶胶来源复杂多样。研究该地区气溶胶光学厚度与人类活动之间的关系，有助于揭示人类活动对大气环境的影响机制，为合理规划城市发展、调整产业结构提供科学依据。在城市化进程中，城市扩张、工业布局和交通规划等因素都会对大气污染产生影响，通过对气溶胶光学厚度的研究，可以评估这些因素对大气环境的影响程度，为城市可持续发展提供决策支持。该地区的地形和气候条件复杂，对气溶胶的传输和扩散产生重要影响。西北部的山地阻挡了冷空气的南下，使得东南部平原地区在冬季容易出现静稳天气，不利于污染物的扩散；夏季的东南季风又会将海洋上的水汽和污染物输送到该地区，增加了大气污染的复杂性。研究地形和气候因素对气溶胶光学厚度的影响，有助于深入理解气溶胶的形成、传输和演化机制，为提高大气污染预测和预警能力提供理论支持。在不同的地形和气候条件下，气溶胶的传输路径和扩散范围会有所不同，通过对这些因素的研究，可以更准确地预测大气污染的发展趋势，及时采取有效的防控措施。京津冀地区拥有较为完善的地面监测站点和丰富的气象数据资源，为研究提供了充足的数据支持。这些地面监测站点能够实时监测大气污染物的浓度、气象要素等信息，与卫星遥感数据相结合，可以对气溶胶光学厚度的反演结果进行验证和校准，提高研究的准确性和可靠性。地面监测站点的实测数据可以作为参考标准，用于评估卫星反演结果的精度，通过对比分析，可以发现反演算法中存在的问题，进一步优化算法，提高反演精度。气象数据则可以用于分析气象条件对气溶胶光学厚度的影响，为研究气溶胶的形成和演化机制提供必要的信息。四、案例研究与数据分析4.2数据获取与处理4.2.1数据来源本研究的数据主要来源于TERRA和AQUA双星的MODIS影像，这些影像数据均从美国国家航空航天局（NASA）的官网数据平台获取。该平台提供了丰富且全面的MODIS数据产品，涵盖了不同分辨率和处理级别的数据，为气溶胶光学厚度的反演研究提供了坚实的数据基础。具体来说，选用的是MODISL1B级数据，这是经过初步处理的原始影像数据，包含了传感器获取的原始辐射亮度信息，其保留了最原始的观测数据特征，对于后续进行高精度的气溶胶光学厚度反演至关重要。在获取TERRA和AQUA双星MODIS影像数据时，根据研究区域京津冀地区的地理位置（北纬36°05′-42°40′，东经113°27′-119°50′），在数据下载平台上精确设置经纬度范围，确保下载的数据能够完整覆盖京津冀地区。同时，为了研究气溶胶光学厚度的时间变化特征，选择了2020年1月至2020年12月这一时间段内的影像数据。在下载过程中，仔细检查数据的完整性和质量，剔除数据缺失或存在明显异常值的影像，保证数据的可靠性和可用性。除了MODIS影像数据，还收集了辅助数据以提高反演精度和分析气溶胶光学厚度与其他环境因素的关系。从美国国家海洋和大气管理局（NOAA）的相关数据平台获取了研究区域的气象数据，包括温度、湿度、气压、风速和风向等信息。这些气象数据对于理解大气环境状况以及其对气溶胶分布和光学特性的影响至关重要。温度的变化会影响气溶胶粒子的挥发性和化学反应速率，湿度的增加可能导致气溶胶粒子吸湿增长，改变其粒径分布和光学性质。通过获取这些气象数据，可以在反演过程中更准确地考虑大气条件的影响，提高气溶胶光学厚度反演的精度。为了准确识别地表类型和地形信息，收集了研究区域的土地利用类型数据和数字高程模型（DEM）数据。土地利用类型数据来源于中国科学院资源环境科学数据中心，该数据将研究区域的土地利用类型划分为耕地、林地、草地、建设用地、水域和未利用地等多种类型。通过这些数据，可以在反演过程中针对不同的地表类型选择合适的反演算法和参数，提高反演结果的准确性。DEM数据则采用了美国地质调查局（USGS）提供的SRTM（ShuttleRadarTopographyMission）数据，其空间分辨率为90米，能够精确反映研究区域的地形起伏情况。地形因素对气溶胶的传输和扩散有着重要影响，在山区，地形的阻挡和抬升作用会改变气溶胶的运动轨迹和浓度分布。利用DEM数据，可以在反演算法中考虑地形对太阳辐射传输和大气运动的影响，进一步优化反演结果。4.2.2数据处理流程在获取TERRA和AQUA双星MODIS影像数据及相关辅助数据后，利用ENVI软件进行一系列的数据处理操作，以确保数据质量和满足后续反演分析的要求。首先进行辐射校正，目的是将MODIS影像的原始数字量化值（DN）转换为具有物理意义的辐射亮度值。在ENVI软件中，调用MODIS数据自带的定标参数文件，结合传感器的辐射响应函数进行辐射定标。MODIS数据的定标参数文件包含了每个波段的增益系数和偏移量等关键信息，通过公式L=G\timesDN+B（其中L为辐射亮度值，G为增益系数，DN为原始数字量化值，B为偏移量），可以将每个像元的DN值准确转换为辐射亮度值。经过辐射校正后，消除了传感器自身特性对辐射量测量的影响，使得不同时间、不同条件下获取的影像数据具有可比性。大气校正是数据处理的重要环节，其主要目的是去除大气分子、气溶胶、水汽等对太阳辐射的散射和吸收作用，从而获取真实的地表反射率信息。在ENVI软件中，选用6S（SecondSimulationoftheSatelliteSignalintheSolarSpectrum）辐射传输模型进行大气校正。该模型能够精确模拟太阳辐射在大气中的传输过程，充分考虑了大气成分、气溶胶类型、太阳天顶角、卫星观测角等多种因素对辐射的影响。在进行大气校正时，需要输入一系列参数，根据研究区域的地理位置和季节特点，选择合适的大气模式，如中纬度夏季模式用于夏季数据的处理，中纬度冬季模式用于冬季数据的处理。根据该区域常见的气溶胶类型，选择相应的气溶胶模式，如大陆型气溶胶模式适用于污染较为严重的大陆地区。还需根据影像中地物的主要类型，确定地表类型参数。通过运行6S模型，对辐射亮度值进行校正，最终得到地表反射率影像。几何校正也是必不可少的步骤，其作用是消除因卫星轨道、姿态、地球曲率等因素导致的影像几何变形，确保影像中地物的位置精度。在ENVI软件中，利用地面控制点（GCP，GroundControlPoint）进行几何校正。通过在影像和参考地图或高分辨率影像上选取同名地物点作为GCP，如明显的道路交叉点、建筑物拐角、河流交汇点等。为了保证校正精度，一般选取不少于20个均匀分布的GCP。利用这些控制点，采用多项式变换等方法建立影像坐标与地理坐标之间的转换关系，对影像进行重采样和几何纠正，使影像的地理坐标与实际地理位置相符。在选择控制点时，尽量选择在不同地形和地物类型上，以提高几何校正的精度。在山区选择山峰、山谷等特征点，在平原地区选择农田边界、道路交叉口等。通过这种方式，可以更好地纠正因地形起伏和卫星观测角度不同导致的几何变形。云检测和云掩膜处理同样至关重要，云的存在会严重干扰气溶胶光学厚度的反演，必须识别并去除影像中的云覆盖区域。在ENVI软件中，采用基于阈值法的云检测算法。该算法利用云在多个波段的反射率和亮度温度特征，设定相应的阈值来识别云像元。在可见光波段，云的反射率通常较高；在热红外波段，云的亮度温度较低。通过设定可见光波段反射率阈值大于0.4，热红外波段亮度温度阈值小于270K等条件，将满足条件的像元识别为云像元。对识别出的云像元进行掩膜处理，将其从影像中剔除，得到无云的影像数据，为后续的气溶胶光学厚度反演提供纯净的数据基础。在云检测过程中，还可以结合纹理分析等方法，进一步提高云像元识别的准确性。云的纹理通常较为均匀，与周围地物的纹理特征有明显差异，通过分析影像的纹理信息，可以更准确地识别出云像元，减少误判。4.3反演结果与分析4.3.1气溶胶光学厚度空间分布通过对TERRA和AQUA双星MODIS影像协同反演，得到了京津冀地区2020年气溶胶光学厚度的空间分布结果（图1）。从空间分布来看，气溶胶光学厚度呈现出明显的区域差异。在京津冀地区的东南部平原区域，包括河北的保定、石家庄、衡水以及天津的部分地区，气溶胶光学厚度普遍较高，大部分区域的AOD值在0.8-1.2之间。这主要是由于这些地区是人口密集和经济活动高度活跃的区域，工业生产、机动车尾气排放以及燃煤等人类活动产生了大量的气溶胶污染物。保定和石家庄是重要的工业城市，拥有众多的钢铁、化工、建材等企业，这些企业在生产过程中会向大气中排放大量的颗粒物和挥发性有机物，这些污染物在大气中经过复杂的物理和化学过程，形成二次气溶胶，导致气溶胶浓度升高。该地区的城市化进程快速推进，城市建设和交通发展也加剧了气溶胶的排放。大量的建筑工地产生的扬尘以及机动车保有量的持续增加，使得交通尾气排放成为气溶胶的重要来源之一。相比之下，京津冀地区的西北部山区，如张家口和承德的部分地区，气溶胶光学厚度较低，AOD值大多在0.4以下。这些地区地形以山地为主，植被覆盖度较高，生态环境相对较好，人类活动强度较小。山区的植被具有吸附和净化空气中污染物的作用，能够有效降低气溶胶的浓度。茂密的森林可以通过叶片表面的吸附作用和植物的生理代谢过程，吸收空气中的颗粒物和有害气体，减少气溶胶的含量。该地区的工业活动相对较少，主要以农业和旅游业为主，污染物排放源相对较少，使得气溶胶光学厚度维持在较低水平。在城市区域，气溶胶光学厚度也存在明显的空间差异。以北京市为例，中心城区的气溶胶光学厚度普遍高于郊区。在北京市的核心区域，如东城区、西城区和朝阳区的部分地区，由于人口密集、交通拥堵以及高楼大厦的阻挡，不利于气溶胶的扩散，导致气溶胶浓度较高，AOD值在1.0-1.5之间。而在北京市的郊区，如延庆区、怀柔区等，植被覆盖度相对较高，人口密度和经济活动强度较低，气溶胶光学厚度相对较低，AOD值大多在0.6-0.8之间。城市的热岛效应也会对气溶胶的分布产生影响，中心城区的高温环境会导致大气边界层不稳定，有利于气溶胶的垂直扩散，但同时也会吸引周边地区的污染物向中心城区汇聚，进一步增加气溶胶浓度。4.3.2时间变化特征对2020年京津冀地区气溶胶光学厚度的时间变化特征进行分析，发现其存在明显的季节变化和日变化规律。从季节变化来看（图2），春季（3-5月）和冬季（12-2月）的气溶胶光学厚度相对较高，而夏季（6-8月）和秋季（9-11月）相对较低。春季气溶胶光学厚度较高主要是由于该季节多沙尘天气。京津冀地区受北方沙尘源地的影响，当冷空气南下时，会将沙尘颗粒携带至该地区，导致气溶胶浓度急剧增加。2020年4月的一次沙尘天气过程中，京津冀地区的气溶胶光学厚度明显升高，部分地区的AOD值超过了1.5。春季气温回升快，蒸发旺盛，土壤干燥，地表沙尘容易被风吹起，也增加了气溶胶的来源。冬季气溶胶光学厚度高的原因主要是供暖期的到来。在冬季，京津冀地区大量使用煤炭等化石燃料进行供暖，燃煤排放的大量污染物，如二氧化硫、氮氧化物和颗粒物等，是气溶胶的重要组成部分。冬季的气象条件不利于污染物的扩散，低温、高湿度和弱风等条件容易形成逆温层，抑制了大气的垂直运动，使得污染物在近地面层积聚，导致气溶胶浓度升高。在2020年12月的一些时段，由于连续的静稳天气，京津冀地区的气溶胶光学厚度持续上升，部分城市出现了严重的雾霾天气。夏季气溶胶光学厚度较低主要得益于有利的气象条件。夏季降水较多，雨水对气溶胶粒子具有冲刷和清除作用，能够有效降低气溶胶浓度。在一次强降雨过程后，京津冀地区的气溶胶光学厚度会明显下降。夏季太阳辐射强，大气对流活动旺盛，有利于气溶胶的垂直扩散和稀释。高温使得大气边界层高度增加，污染物能够在更大的空间范围内扩散，从而降低了近地面的气溶胶浓度。秋季气溶胶光学厚度相对较低，一方面是因为秋季天气晴朗，大气扩散条件较好；另一方面，秋季植被仍然保持一定的覆盖度，对气溶胶有一定的吸附和净化作用。随着秋季的推进，气温逐渐降低，大气稳定性增加，气溶胶浓度也会逐渐升高。在日变化方面，气溶胶光学厚度呈现出明显的早晚高、中午低的特征。早晨（6-9时），由于夜间大气边界层稳定，污染物容易积聚，加上机动车出行等人类活动的增加，导致气溶胶浓度升高。随着太阳辐射增强，大气边界层逐渐抬升，对流活动加剧，气溶胶在垂直方向上得到扩散，浓度逐渐降低，在中午（12-14时）达到最低值。下午（15-18时），随着太阳辐射减弱，大气边界层逐渐稳定，气溶胶浓度又开始回升。晚上（19-2时），大气边界层进一步稳定，人类活动虽然减少，但污染物的扩散条件变差，气溶胶浓度维持在较高水平。在2020年5月10日的观测中，北京城区早晨的气溶胶光学厚度为0.9，中午降至0.6，下午又回升至0.8。4.3.3与其他数据对比验证为了验证TERRA和AQUA双星MODIS影像协同反演得到的气溶胶光学厚度的准确性，将反演结果与AERONET地基站点数据进行对比分析。选取了京津冀地区的3个AERONET站点，分别为北京站点、香河站点和天津站点，这些站点能够代表不同的地理位置和环境条件。在对比过程中，首先对AERONET站点数据进行筛选，选择与MODIS影像观测时间相近（时间差在1小时以内）的数据进行对比。利用相关系数（R）和均方根误差（RMSE）来评估反演结果与地基站点数据之间的相关性和误差。相关系数R越接近1，表示两者之间的相关性越好；均方根误差RMSE越小，表示反演结果与地基站点数据之间的差异越小。对比结果表明（表1），在3个站点上，反演结果与AERONET地基站点数据之间具有较高的相关性，相关系数R均在0.8以上。北京站点的相关系数R达到了0.85，均方根误差RMSE为0.12；香河站点的相关系数R为0.83，均方根误差RMSE为0.15；天津站点的相关系数R为0.81，均方根误差RMSE为0.18。这说明本研究采用的双星MODIS影像协同反演方法能够较为准确地反演陆地上空气溶胶光学厚度，反演结果与地基站点实测数据具有较好的一致性。在某些情况下，反演结果与地基站点数据之间仍存在一定的差异。在一些复杂的气象条件下，如强风、暴雨等，大气中的气溶胶分布和光学特性会发生快速变化，可能导致反演结果与地基站点数据之间的偏差增大。地基站点