高分五号大气痕量气体差分吸收光谱仪对流层NO2柱浓度遥感反演研究

搭载于“高分五号”（GF-5）卫星上的痕量气体差分吸收光谱仪（EMI）是一台星下观测的高光谱载荷,测量紫外和可见光光谱范围的地球后向散射辐射,设计用于大气痕量气体的探测。本研究基于EMI在VIS1通道的实测光谱,利用差分光学吸收光谱（DOAS）方法进行了对流层NO₂柱浓度反演,展示了基于EMI载荷的对流层NO₂柱浓度反演结果,与同类载荷产品进行交叉验证,并利用地基观测结果进行了地基验证。研究表明,EMI反演结果与OMI、TROPOMI具有较好的空间分布一致性和较低的相对偏差,与TROPOMI具有较好的时间变化一致性。地面验证结果表明EMI NO₂反演结果具有较高的精度。本研究证明了EMI在全球NO₂监测方面的能力,可以为中国后续污染气体探测载荷的设计和反演算法的开发提供参考。     

基于OMI的亚洲地区对流层NO2高分辨率反演产品POMINO v2.1及其与其他产品的定量对比

对流层二氧化氮（NO₂）是一种重要的痕量污染气体。现有基于OMI卫星探测器、覆盖亚洲地区的NO₂公开产品QA4ECV、OMNO2和POMINO受到广泛使用,然而对于这3个产品的差异的定量认识仍然不足。我们将前期开发的POMINO产品进行了改进和优化,更新至v2.1并将反演区域扩大至覆盖东亚、东南亚和南亚大部分地区,随后定量分析了QA4ECV、OMNO2 v4和POMINO v2.1对流层NO₂垂直柱浓度在2015年—2020年在不同采样条件下的异同。结果显示,POMINO版本的更新对自身NO₂柱浓度反演结果的整体影响较小（<10%）。当3个产品均基于POMINO v2.1进行一致采样时,产品之间在整个亚洲区域的平均差异约为10%,在京津冀等污染地区的差异最高可达40%。3个产品均显示,京津冀地区NO₂柱浓度在五年间下降了约30%,而长三角地区的变化趋势较小。当3个产品进行分别采样时,POMINO v2.1的有效数据量比其他两个产品增加了11%—44%,特别是更好地保留了重污染情形下的NO₂数据,从而降低了采样引起的对NO₂污染水平的系统性低估。本研究对于NO₂卫星产品差异的定量分析有助于认识氮氧化物污染状况以及排放和影响评估。     

2005年-2014年京津冀对流层NO2柱浓度时空变化及影响因素

基于臭氧监测仪OMI对流层NO₂柱浓度产品研究了京津冀城市群2005年-2014年NO₂时空变化及影响因素:(1)10年柱浓度年均增长率为3.35%,且年度波动较大。忽略2008年国家奥运会的环境治理所引起的变化,2005年-2011年NO₂呈逐渐升高趋势;2012年-2014年呈逐渐降低趋势,以2014年下降最为显著。(2)呈西北低东南高的趋势。燕山-太行山山系以北的承德和张家口市浓度较低,山系以南主要有北京-天津-唐山与石家庄-邢台-邯郸两个污染中心。(3)京津冀北部三面环山不利于NO₂的扩散,夏季丰富的降水对NO₂具有显著湿沉降作用。(4)通过相关性分析、文献及国家政策印证等方法,发现地区产业及能源结构很大程度上决定了地区的污染来源。北京市10年来第三产业一直处于主导且稳步提高,煤炭消耗量低,但汽车保有量增加了1.5倍,主要来源为机动车尾气排放;天津市第二产业比第三产业比重略高,煤炭消耗量是北京的两倍之余,但汽车保有量仅是北京市的一半,由此可知工业排放和机动车是共同来源;河北省第二产业比重很高,燃煤量占京津冀地区的80.6%,河北省工业排放是NO₂的主要来源,但近几年随着机动车保有量的剧增,其尾气排放分担率不可小觑。     

奥运大气环境遥感

介绍了中国科学院在遥感应用研究所楼顶建立的空气质量超级监测站, 并主要介绍了2007 年和2008 年遥感监测的结果。监测结果表明, 2007 年好运测试期间北京机动车实行单双号限制, 使北京市气溶胶光学厚度和NO₂ 的柱浓度下降30%左右；2008 年从6 月份开始, 由于北京市的控制措施使NO₂ 柱浓度有了比较大的下降。卫星遥感监测表明, 通过有效的环境控制可以降低北京市及其周边的污染浓度。     

融合星地多源数据资料的长三角地区高分辨率时空无缝PM2.5浓度数据

大气污染物浓度全方位动态监测是进行区域大气污染精细化防控的重要前提。为开展长三角地区小时分辨率PM_2.5浓度无缝制图,本研究通过耦合AOD缺失信息重建与多模数据融合技术,建立了一套能够有效集成卫星遥感、地面观测、数值模拟等多源异构数据资料的近地面PM_2.5浓度无缝制图方案,并据此生产了2015年—2020年长三角地区小时分辨率无缝PM_2.5浓度格点数据产品。结果表明：本研究生产的PM_2.5浓度无缝格点产品与国控站点观测数据的交叉验证相关系数达0.9,平均偏差不超过10 μg·m^-3。较于空间分布不均且相对稀疏的站点观测PM_2.5浓度资料,面域无缝PM_2.5浓度格点数据更能有效揭示长三角地区PM_2.5污染的时空变化特征;在2015年—2020年研究期内,其平均下降速率超过3 μg·m^-3·a^-1。本研究发展的PM_2.5浓度无缝制图方法和生产的相关数据产品有望为区域灰霾污染防控和PM_2.5暴露健康风险评估研究提供方法参考和基础数据支撑。     

基于卫星遥感的近地面PM2.5浓度反演进展

细颗粒物PM_2.5（Fine Particulate Matter）是影响空气质量和公共健康的关键因素之一。高时空分辨率的PM_2.5数据是公共健康风险评估和流行病学研究的基本需求。相较于地面站点,卫星遥感技术具有连续观测、宽覆盖和低成本的优势,基于卫星气溶胶光学厚度AOD（Aerosol Optical Depth）反演PM_2.5质量浓度的方法已成为热点。本研究概述了卫星AOD产品反演PM_2.5浓度的原理,介绍了用于PM_2.5反演的主要卫星AOD产品及其反演精度;总结了现有的PM_2.5估算方法及其优缺点,指出目前PM_2.5反演研究存在的问题;提出未来PM_2.5反演方向主要集中在高时空分辨率的PM_2.5浓度重建、基于激光雷达数据的三维PM_2.5浓度反演及PM_2.5化学组分反演等方向。比例因子法、物理机理模型和统计模型这3种方法都能在不同时期不同程度地准确估算PM_2.5浓度,代表了那个时期较为前沿的研究热点,但比例因子法和物理机理模型因其自身的局限性而应用较少,而统计模型因其独特的时间或时空异质性的可描述性和强大的非线性描述能力的优势而被广泛应用并不断改进。目前PM_2.5反演研究存在的问题主要有3种：（1）卫星AOD的非随机缺失问题造成估算的PM_2.5数据缺失;（2）反演模型的精度问题;（3）PM_2.5的化学成分估算问题。基于此,本文为了准确揭示近地面PM_2.5的时空变化趋势,提高基于卫星AOD产品的近地面PM_2.5反演研究的准确性,提出了几点未来的研究方向：首先,新型的高空间（如风云四号、高分五号）、高时间分辨率（Himawari-8/9）卫星AOD产品在PM_2.5的精细化估算研究上具有很大优势,这对于高时空分辨率的PM_2.5浓度重建具有重要意义;其次,随着大气探测技术的发展,星载、机载及地基激光雷达都能够获取垂直分布信息,搭载在无人机上的颗粒物传感器可实现PM_2.5垂直方向上的监测,将其与光学遥感卫星数据及地面监测数据结合,可实现三维的PM_2.5浓度反演;最后,PM_2.5化学组分信息对于分析污染成因、暴露特征等尤其重要,其时空变化趋势研究是一个重要的发展方向,然而,地面PM_2.5组分观测站网仍不完善,如何克服卫星遥感估算中对地面站网的依赖,实现PM_2.5化学成分的高精度反演需要进一步研究。通过本研究,有助于进一步了解不同PM_2.5估算方法的原理机制及其优缺点,为基于卫星AOD产品反演近地面PM_2.5浓度的新的发展方向提供启示,提升近地面PM_2.5浓度反演的精度及时空分辨率。     

资源三号02星激光测高数据质量分析

近几年来,中国的对地观测卫星激光测高技术发展较快,资源三号02星上搭载了国内首台对地观测的卫星激光测高试验载荷,在后续的高分七号、陆地生态系统碳监测卫星上均装备业务化应用的激光测高仪,开展国产对地观测卫星激光测高数据质量分析研究非常必要,能有效填补国内空白并推动相关技术发展。本文对资源三号02星载荷的激光测高数据进行了严密几何处理,重点从夜间观测、卫星小角度侧摆、海面区域等不同条件下的数据质量进行了分析,对数据的平均可用率进行了统计分析。试验表明,资源三号02星激光测高仪能获得约30%有效测高数据,在夜间数据有效利用率略有提高,在大型水面也能获得有效数据,经后处理后在平坦地区验证的绝对高程精度优于1.0 m,部分点高程精度优于0.5 m,硬件本身的测距精度和卫星姿态测量误差是主要的误差源。相关结论对于后续国产卫星激光测高载荷的研制以及数据后处理和应用具有参考价值。     

重大公共卫生事件影响下广州气溶胶分布变化

通过暗像元法对广州市2020-01-03、2020-01-31、2020-02-22、2020-02-24、2020-03-15、2020-04-09、2020-04-27的MODIS影像进行气溶胶光学厚度反演,分析新型冠状病毒肺炎疫情期间广州市气溶胶光学厚度在时空上的变化情况。结果表明广州市全域气溶胶光学厚度在1月31日达到最低,期间正值广东省启动重大突发公共卫生事件一级响应,城市人为活动大幅减少。2月22日及以后,全市气溶胶光学厚度总体上升,此时城市进入复工复产阶段。研究表明利用卫星遥感手段快速反演气溶胶光学厚度,可在重大突发公共卫生事件发生时监测城市人为活动情况,辅助政府决策和宏观调控。     

北京城市下垫面大气CO2反演算法

大气温室气体监测仪GMI（Greenhouse gases Monitor Instrument）是高分五号（GF-5）卫星载荷之一，主要用于全球温室气体含量监测和碳循环研究。高精度反演是卫星大气CO₂遥感的基本需求。地表反射率影响卫星遥感辐射量及辐射传输过程中的地气耦合过程，严重制约着CO₂的反演精度，针对GMI开发高精度的大气CO₂反演算法，地表反射是一个需要重点考虑的因素。城市是CO₂重要的发射源，且城市下垫面存在明显的二向反射特性，加上城市大气条件不良，复杂的地气耦合效应存在这都考验反演算法的准确性和鲁棒性。本文针对北京城市地区，利用2011年—2016年共5年的MODIS（MODerate-resolution Imaging Spectroradiometer）地表二向反射分布函数BRDF（Bidirectional Reflectance Distribution Function）数据，构建了适合利用单次观测数据反演的BRDF模型，并提出一种同时反演地表BRDF参数和大气CO₂含量的算法。结果表明在550 nm波长处气溶胶光学厚度AOD（Aerosol Optical Depth）小于0.4时，大部分GMI模拟数据的反演误差控制在0.5%（~2 ppm）内。利用GOSAT （Greenhouse gases Observing SATellite）实测数据的反演结果与修正后的日本国立环境研究所NIES（National Institute for Environmental Studies）反演结果进行对比，其平均误差为1.25 ppm，相关性达到0.85。本算法满足GMI数据在北京城市区域高精度CO₂反演的需求，并使得反演高值气溶胶区域数据成为可能，增加了GMI观测数据的利用率。     

资源一号02D高光谱影像内陆水体叶绿素a浓度反演

2019-09-12成功发射的资源一号02D卫星(ZY-102D)搭载了新一代可见短波红外高光谱相机AHSI(Advanced Hyperspectral Imager),其丰富的细分波段和较高的空间分辨率在内陆湖库水质监测方面具有较大潜力,但数据可用性有待分析和验证.本研究以中国华东和华北平原的典型富营养湖库(太湖、...     

Increased accuracy of geostatistical prediction of nitrogen dioxide in the United Kingdom with secondary data

Five techniques were used to map nitrogen dioxide (NO₂) concentrations in the United Kingdom. The methods used to predict from point data, collected as part of the UK NO₂ diffusion tube network, were local linear regression (LR), inverse distance weighting (IDW), ordinary kriging (OK), simple kriging with a locally varying mean (SKlm) and kriging with an external drift (KED). These techniques may be divided into two groups: (i) those that use only a single variable in the prediction process (IDW, OK) and (ii) those that make use of additional variables as a part of prediction (LR, SKlm and KED). Nitrous oxides emission data were used as secondary data with LR, SKlm and KED. It was concluded that SKlm provided the most accurate predictions based on the summary statistics of prediction errors from cross-validation.     

Evaluating Hyperspectral Imager for the Coastal Ocean (HICO) data for seagrass mapping in Indian River Lagoon,FL

Differentiation between benthic habitats, particularly seagrass and macroalgae, using satellite data is complicated because of water column effects plus the presence of chlorophyll-a in both seagrass and algae that result in similar spectral patterns. Hyperspectral imager for the coastal ocean data over the Indian River Lagoon, Florida, USA, was used to develop two benthic classification models, Slope_RED and Slope_NIR. Their performance was compared with iterative self-organizing data analysis technique and spectral angle mapping classification methods. The slope models provided greater overall accuracies (63–64%) and were able to distinguish between seagrass and macroalgae substrates more accurately compared to the results obtained using the other classifications methods.     

CO2 concentration and its spatiotemporal variation in the troposphere using multi-sensor satellite data,carbon tracker,and aircraft observations

Satellite-based atmospheric CO₂ observations have provided a great opportunity to improve our understanding of the global carbon cycle. However, thermal infrared (TIR)-based satellite observations, which are useful for the investigation of vertical distribution and the transport of CO₂, have not yet been studied as much as the column amount products derived from shortwave infrared data. In this study, TIR-based satellite CO₂ products – from Atmospheric Infrared Sounder, Tropospheric Emission Spectrometer (TES), and Thermal And Near infrared Sensor for carbon Observation – and carbon tracker mole fraction data were compared with in situ Comprehensive Observation Network for Trace gases by AIrLiner (CONTRAIL) data for different locations. The TES CO₂ product showed the best agreement with CONTRAIL CO₂ data resulting in R² ~ 0.87 and root-mean-square error ~0.9. The vertical distribution of CO₂ derived by TES strongly depends on the geophysical characteristics of an area. Two different climate regions (i.e., southeastern Japan and southeastern Australia) were examined in terms of the vertical distribution and transport of CO₂. Results show that while vertical distribution of CO₂ around southeastern Japan was mainly controlled by horizontal and vertical winds, horizontal wind might be a major factor to control the CO₂ transport around southeastern Australia. In addition, the vertical transport of CO₂ also varies by region, which is mainly controlled by anthropogenic CO₂, and horizontal and omega winds. This study improves our understanding of vertical distribution and the transport of CO₂, both of which vary by region, using TIR-based satellite CO₂ observations and meteorological variables.     

Spatiotemporal mixed effects modeling for the estimation of PM2.5 from MODIS AOD over the Indian subcontinent

ABSTRACT

The physical processes associated with the constituents of the troposphere, such as aerosols have an immediate impact on human health. This study employs a novel method to calibrate Aerosol Optical Depth (AOD) obtained from the MODerate resolution Imaging Spectrometer (MODIS – Terra satellite) for estimating surface PM_2.5 concentration. The Combined Deep Blue Deep Target daily product from the MODIS AOD data acquired across the Indian Subcontinent was used as input, and the daily averaged PM_2.5pollution level data obtained from 33 monitoring stations spread across the country was used for calibration. Mixed Effect Models (MEM) is a linear model to deal with non-independent data from multiple levels or hierarchy using fixed and random effects of dependent parameters. MEM was applied to the dataset obtained for the period from January to August 2017. The MEM considers a fixed and random component, where the random components model the daily variations of the AOD – PM_2.5 relationships, site-specific adjustment parameters, temporal (meteorological) variables such as temperature, and spatial variables such as the percentage of agricultural area, forest cover, barren land and road density with the resolution of 10 km × 10 km. Estimation accuracy was improved from an R² value of 0.66 from our earlier study (when PM_2.5 was modeled against only AOD and site-specific parameters) toR² value of 0.75 upon the inclusion of spatiotemporal (meteorological) variables with increased % within Expected Error from 18% to 35%, reduced Mean Bias Error from 3.22 to 0.11 and reduced RMSE from 29.11 to 20.09. We also found that spline interpolation performed better than IDW and Kriging inefficiently estimating the PM_2.5 concentrations wherever there were missing AOD data. The estimated minimum PM_2.5 is 93 ± 25μg/m³ which itself is in the upper limit of the hazardous level while the maximum is estimated as 170 ± 70μg/m³. The study has thus made it possible to determine the daily spatial variations of PM_2.5 concentrations across the Indian subcontinent utilizing satellite-based AOD data.     

Application of geographical information system for mapping/monitoring seismological hazards in Pakistan

This study attempts to use the geographic information system (GIS) technique to map and understand the tectonics and crustal structures of Pakistan. Maps of surficial tectonic features and seismological parameters including Moho depth, P_n velocity and P_g velocity are complied. Based on the seismological data-set of the country the earthquake hazard map of Pakistan is also presented by applying regression technique on seismological, geological and topographical parameters. A case study of 8 October 2005 earthquake is used to validate the hazard map. It is envisaged that the developed GIS database would help policy-makers and scientists in natural hazard evaluation, seismic risk assessment and understanding of earthquake occurrences in Pakistan.     

GF-5 AIUS切高校正算法研究与反演验证

利用红外掩星获取温压及其大气成分廓线,需要获得精确的指向信息,针对高分五号大气环境甚高分辨率探测仪（GF-5 AIUS）一级光谱信息,本文首先分析了不同切高的光谱变化情况,在此基础上,利用查找表方法进行切高校正。在低切高段（10—20 km）,选取N₂吸收波段,模拟6—40 km高度上的透过率光谱,通过查找表校正方法对一级数据中的低层切高进行校正。在高切高段（20—100 km）,模拟6—90 km切高上的透过率光谱,利用O₃吸收波段（20—90 km）进行校正。校正算法采取全局最小均方根误差的方法进行统计实验,确定最小校正半径为15 km。最后选取O₃、HCl、N₂O共3种大气成分的反演波段,进行校正后的透过率数据对比验证与反演产品精度的验证。结果表明：在O₃、HCl、N₂O共3种成分的反演通道上,校正后的观测透过率与模拟的透过率数据一致性较好,两者最大偏差小于0.1,切高校正算法得到的结果与模拟透过率吻合度较好;同时基于Aura MLS（Microwave Limb Sounder）、ACE-FTS（Atmospheric Chemistry Experiment Fourier-Transform Spectrometer）Level 2产品进行了O₃、HCl、N₂O单廓线交叉验证。结果表明：GF-5 AIUS反演的单廓线O₃与MLS、ACE-FTS相对偏差在30 km以下小于50%,在30—60 km小于20%;单廓线HCl与MLS、ACE-FTS相对偏差整体小于20%;单廓线N₂O与MLS、ACE-FTS相对偏差30 km以下小于20%,30—60 km相对偏差较大。     

利用遥感数据评价燃煤电厂空气质量

卫星观测数据可以评价燃煤电厂的空气质量等级。NO₂、SO₂ 和烟尘是燃煤电厂排放的主要污染物,本文利用卫星遥感观测的NO₂、SO₂和气溶胶光学厚度AOD(Aerosol Optical Depth)开展燃煤电厂空气质量评价。以中国华北地区为实验区,分析对比了3种污染物不同时间分辨率和空间分辨率的污染状况,确定了单因子的5级分级标准,根据燃煤电厂排放污染物的权重不同,提出了评价近地表空气质量状况的模型。本文综合考虑3种污染因子来反映电厂空气质量,有利于提高评价的准确性以及反应信息的全面性。结果表明,该模型能正确反映不同地区电厂的空气质量特点。     

Application of hyperspectral remote sensing reflectance data to photochemical rate calculations in the Duplin River,a tidal river on the coast of Georgia,USA

Estuaries are photochemically dynamic environments with high carbon loads and relatively small areas. The small area poses problems for large-scale satellite-based remote sensing calculations, where the resolution is too coarse to distinguish land from water. Airborne remote sensing instruments have the potential to reveal the dynamics of these areas with fine-scale resolution. In June 2006, hyperspectral remote sensing imagery, using an AISA Eagle instrument, was collected over the tidal Duplin River, Georgia, USA. A dark-water updated version of the SeaUV algorithm was applied to the AISA remote sensing image to determine diffuse attenuation constants in the ultraviolet and calculate surface photochemical production rates of two inorganic products – carbon monoxide (CO) and carbon dioxide (CO₂). For an average day in June at the study site, the modeled photoproduction rates for CO₂ and CO averaged ~7 × 10^?1 nmol C/day/cm³ and ~3.5 × 10^?2 nmol C/day/cm³, respectively.     

Implications of ocean color in the upper water thermal structure at NINO3.4 region: a sensitivity study for optical algorithms and ocean color variabilities

Chlorophyll a (Chl-a) has been the most commonly used biomass metric in biological oceanographic processes. Although limited to two-dimensional surfaces, remote-sensing tools have been successfully providing the most recent state of marine phytoplankton biomass to better understand bottom-up processes initiating daily marine material cycles. In this exercise, ocean color products with various time-scales, derived from Sea-Viewing Wide Field-of-View Sensor (SeaWiFS), were used to investigate how their bio-optical properties affect the upper-ocean thermal structure in a global ocean modeling framework. This study used a ¼-degree Hybrid Coordinate Ocean Model forced by hourly atmospheric fluxes from the Climate Forecast System Reanalysis at National Oceanic Atmospheric Administration. Three numerical experiments were prepared by combining two ocean color products – downwelling diffuse attenuation coefficients (K_dPAR) and chlorophyll a (Chl-a) – and two shortwave radiant flux algorithms. These three runs are: (1) KparCLM, based on a 13-year long-term climatological K_dPAR derived from SeaWiFS; (2) ChlaCLM, based on a 13-year long-term Chl-a derived from SeaWiFS; and (3) ChlaID, which uses the inter-annual time-series of monthly-mean SeaWiFS Chl-a product. The KparCLM experiment uses a Jerlov-like two-band scheme; whereas, both ChlaCLM and ChlaID use a two-band scheme that considers inherent (absorption (a) and backscattering (b_b) coefficients) and apparent optical properties (downwelling attenuation coefficient (K_d) and solar zenith angle (θ, varying 0–60°)). It is found that algorithmic differences in optical parameterizations have a bigger impact on the simulated temperatures in the upper-100 m of the eastern equatorial Pacific, NINO3.4 region, than other parts of the ocean. Overall, the K_dPAR-based approach estimated relatively low surface temperatures compared to those estimated from the chlorophyll-based method. In specific, this cold bias, pronounced in the upper 20–30 m, is speculated to be due to optical characteristics of the algorithm and K_dPAR products, or due to nonlinear hydrodynamical processes involving displacement of mixed-layer depth. Comparisons between each experiment against Global Ocean Data Assimilation System (GODAS; Behringer and Xue 2004) analyses find that KparCLM-based simulations have lower mean differences and variabilities with higher cross-correlation coefficients compared to ChlaCLM- and ChlaID-based experiments.