融合高分一号、风云四号及葵花-8的光合有效辐射遥感方法及产品研究

光合有效辐射是植被生产力模型的重要输入参数，也是陆地生态系统模型及生物地化模型等的重要特征参量。目前的全球或区域尺度光合有效辐射产品在起伏地表情形的精度及可用性尚存在不足，高分系列卫星的发射，为高空间分辨率光合有效辐射遥感产品反演提供了可能。本文结合高分一号卫星(GF-1)反演的地表反照率产品，气溶胶光学厚度产品，日本新一代静止卫星Himawari-8及国产风云四号卫星(FY-4)的云光学厚度产品及美国的SNPP卫星反演的气溶胶光学厚度产品(用于填补GF气溶胶产品的空值区)，发展了一种基于参数化模型的光合有效辐射遥感反演方法。晴空条件下，主要考虑气溶胶和瑞利散射对光合有效辐射的衰减；云天情形，则主要考虑云光学厚度参数对入射辐射的影响，以球形粒子的Mie散射理论和平面平行大气辐射传输原理计算总透过率；针对起伏地表，分别计算地形对坡面入射角度的影响、对直接辐射的遮蔽以及对散射辐射的增强/衰减效应，从而得到高分卫星光合有效辐射遥感反演模型及产品。利用河北怀来试验站、西南大学柑研所及黑河流域地表过程综合观测网的连续观测数据对光合有效辐射产品进行了对比验证，其相关系数为0.87，平均偏差为1.5...     

扩展多波长偏振测量的太阳—天空辐射计观测网

太阳辐射和大气气溶胶的地基遥感观测在气候变化、大气环境和对地观测等领域有重要意义。结合多波长偏振测量技术在气溶胶观测中的应用,本文介绍了在中国典型区域建立的太阳—天空辐射计观测网SONET(Sunsky radiometer Oberservation NETwork)。主要内容包括:(1)多波长偏振太阳—天空辐射计的参数、基本原理和功能,以及据此建立的中国区域太阳—天空辐射计观测网及其站点的基本情况;(2)太阳直射、太阳主平面和平纬圈几何条件下,太阳—天空辐射和偏振的观测模式及其测量参数;(3)针对气溶胶、水汽等大气成分观测发展的太阳直射辐射、天空光辐亮度、线偏振度等观测量的定标和数据处理方法,并给出在晴朗和灰霾等天气下的SONET观测结果示例;(4)基于在SONET典型站点连续的较长期遥感观测,初步建立了城市、乡村和灰霾情况下的整层大气气溶胶参数模型。在扩展多波长偏振测量的基础上,SONET还发展了自动观测、快速数据处理、定期标定和维护等观测网运行关键技术,可为典型大气气溶胶模型构建、卫星气溶胶产品反演和验证、大气环境遥感监测等应用提供支撑。     

暗目标法的Himawari-8静止卫星数据气溶胶反演

Himawari-8(H8)是由日本气象厅发射的新一代静止气象卫星,可实现10 min/次的高频次对地观测,搭载的AHI(Advanced Himawari Imager)传感器设置有与MODIS暗目标气溶胶反演算法所需的类似波段。本文参考暗目标算法构建了针对该卫星传感器的陆地气溶胶反演算法:首先,通过基于地基站点观测数据的精确大气校正,统计得到短波红外与可见光波段的地表反射率比值关系,将此作为先验知识用于地—气解耦时的反射率估计;然后,初步假设大陆型气溶胶类型,利用辐射传输模型建立查找表;最后,通过模拟与卫星观测的表观反射率误差最小实现气溶胶光学厚度反演解算。选取2016年5月覆盖京津冀地区的观测数据进行测试,将反演结果与对应时间的MODIS气溶胶光学厚度产品进行对比验证,空间分布趋势一致、相关性较高,相关系数R达到0.852;通过与地基观测网AERONET站点实测数据对比验证,所有站点的相关系数R~2均大于0.88,精度较高。利用反演的高时间分辨率产品,分析了京津冀地区的大气空间分布和日变化情况,结果表明:采用暗目标法对H8静止卫星陆地气溶胶光学厚度反演具有一定的潜力和可行性,能反映气溶胶的高时间变化信息,有望成为大气环境污染变化监测新的重要手段。     

气溶胶大气辐射校正区域性分析

针对使用6S模型中标准气溶胶模型进行大气辐射校正时结果较差的问题,该文提出使用实测气溶胶模型代替其进行计算。利用AERONET数据产品分析了香河、榆中等部分站点的气溶胶成分及其区域性分布特点,证明了中国气溶胶的成分具有明显的区域特性;并对6S模型中标准气溶胶模型与香河地基观测气溶胶光学参数进行比较,结果表明,两者具有较大的差异;最后使用不同的气溶胶模型校正Landsat5TM影像,结果表明,采用实测气溶胶数据进行大气辐射校正结果较优。     

灰霾污染状况下气溶胶组分及辐射效应的遥感估算

本文利用地基遥感观测数据和气溶胶模型,分析了北京地区2013年1月份灰霾期间气溶胶中不同成分的光学与物理化学特性,并使用大气辐射传输模式估算了气溶胶中不同成分对直接辐射强迫的贡献。结果显示,在严重灰霾时期,大气中的硫酸盐等水溶性气溶胶光学厚度（440nm）显著增加,对气溶胶光学厚度的贡献达到70%以上,同时黑碳气溶胶对总光学厚度（440nm）的贡献也呈现增长趋势,是晴朗天气下的8倍以上。严重灰霾（28日）时气溶胶的直接辐射强迫显著增加,但气溶胶中不同成分对辐射效应的贡献不同,黑碳成分在气溶胶对大气层的辐射强迫中占主要贡献（57%）,而水溶性成分则在气溶胶对地表的辐射强迫中占主要贡献（60%）。     

Landsat 8地表温度反演及验证—以黑河流域为例

地表温度是区域和全球尺度地表物理过程的一个重要参数,目前已有的地表温度产品空间分辨率较低,缺乏高空间分辨率的地表温度产品。Landsat系列卫星提供了大量免费的高空间分辨率遥感数据,然而对应的高空间分辨率地表温度产品还未见到,为了获取长时间序列的高空间分辨率地表温度数据,针对Landsat 8 TIRS数据提出了一个物理单通道地表温度反演算法。该算法首先利用ASTER全球地表发射率产品(ASTER GED)结合Landsat 8地表反射率产品计算Landsat 8影像的地表发射率,然后利用快速辐射传输模型RTTOV结合MERRA大气廓线数据对热红外影像进行大气校正,最后利用物理单通道地表温度反演算法得到地表温度。利用黑河流域HiWATER试验2013年—2015年15个站点的实测地表温度数据对本文方法和普适性单通道算法进行了验证,同时对验证站点的空间异质性进行了分析。结果表明,本文方法和普适性单通道算法估算的地表温度整体精度均较高,能够获取高精度、高空间分辨率的地表温度数据,可以服务于城市热岛效应、地表蒸散发估算等相关研究。     

内蒙古辐射校正场特性评价与应用潜力分析

利用2010年6月在内蒙古贡格尔辐射校正场(简称内蒙古辐射校正场)实测的地表反射率和大气气溶胶数据,分析了该辐射校正场的地表反射率和大气气溶胶特性,并与敦煌辐射校正场相关数据进行了对比分析.结果表明:内蒙古辐射校正场地表反射率的相对方差小于5％,表现出较好的均一性;在天气晴朗的情况下,各通道的大气气溶胶光学厚度可达0....     

利用云下阴影实现陆地上空气溶胶和地表反射率的同时反演—理论方法和模拟

卫星探测信号包含大气中分子和粒子的散射贡献以及地表反射的贡献,在陆地上空二者的贡献相当,并且陆地地表反射率在时间和空间上极度不均一,因此,很难区分二者的各自贡献从而定量提取大气气溶胶和地表反射率,陆地上空气溶胶的反演也一直是一个极具挑战性的课题.而高分辨率卫星资料如TM5的可见光通道能够很好地区分云和云下阴影,如果云是不透光的,在阴影上空,卫星信号仅包含大气散射贡献和地表漫反射贡献,而在邻近的非阴影区上空,卫星探测信号还包含地表直接反射的贡献,根据这个原理,利用辐射传输模式分析了阴影区和非阴影区上空卫星探测的辐射量差别与地表反射率和大气气溶胶的关系,提出一种利用云下阴影来同时提取阴影上空大气气溶胶和地表反射率的单波长反演方案,并对气溶胶单次散射反照率,散射相函数,测量精度以及地表反射率的不均一性进行了敏感性分析.     

MODIS产品估算复杂地形下的光合有效辐射

提出基于MODIS大气和地表产品反演高精度复杂地形下的光合有效辐射模型。该方法是对普通的大气传输模型进行简化, 在模型中分别考虑了削减光合有效辐射的5个主要因子: 瑞利散射、臭氧和水汽吸收、气溶胶散射以及地表和大气间的多次反射。从3个方面考虑了复杂地形对于光合有效辐射的影响: (1) 坡度和坡向; (2) 天空可视角; (3) 邻近地形的附加辐射。提取2006年晴空条件下的估算结果与中国生态网络山东禹城站和长白山站的实测值的比较, 得到相关系数分别为0.924和0.9, 表明模型反演结果能够较好地解译实际光合有效辐射的变化。通过比较考虑地形和不考虑地形2种情况得到的光合有效辐射, 定量地分析了地形对光合有效辐射的影响。     

地表反射率产品支持的GF-1 PMS气溶胶光学厚度反演及大气校正

GF-1卫星PMS(GF-1 PMS)数据具有高空间分辨率、短重访周期的特点,可以在地表类型识别、参数提取中发挥重要作用。但由于缺少2.1μm附近的短波红外波段,使得气溶胶反演时地表反射率的精确确定非常困难,从而导致其高精度大气校正难以开展,限制了该数据的应用。本文提出了一种地表反射率数据支持的气溶胶反演方法,用于GF-1PMS数据的大气校正。其基本思想是:使用现有的地表反射率数据集为GF-1PMS数据提供地表反射率,用于确定GF-1PMS图像中浓密植被像元(DDV)的分布,基于确定的浓密植被像元反演气溶胶光学厚度(AOD),并用于大气校正。这里使用的地表反射率数据集为合成的无云MODIS地表反射率产品,对GF-1PMS数据做了空间尺度的转换。为降低两类数据配准误差对地表反射率确定的影响,提出了使用区域NDVI分布百分比匹配的方法,回避了像元的直接匹配,为GF-1PMS数据提供DDV的空间分布。为验证该方法的有效性,利用北京、太湖两个AERONET站点观测的气溶胶光学厚度对气溶胶反演结果进行精度验证,结果表明,气溶胶反演算法精度较高,稳定性较强。AOD反演结果应用于北京和敦煌地区的GF-1PMS数据大气校正,获得的地表反射率与地面实测的地表反射率的误差低于0.015,且大气校正后影像对比度明显提高。     

Preliminary validation of GLASS-DSSR products using surface measurements collected in arid and semi-arid regions of China

Global Land Surface Satellite-downward surface shortwave radiation (GLASS-DSSR) products have been routinely produced from 2008–2010 based on an improved look-up table algorithm, which explicitly accounts for the variations of cloud optical depth, water vapor content, and elevation. In this study, we validated and assessed the accuracy of these products in arid and semiarid regions of China. Toward this goal, observation data-sets provided by the Arid and Semiarid Region Collaborative Observation Project as well as four other metrological sites were collected, chosen, and preprocessed for the final validation. Due to the possible effect of spatial collocation and the strong adjacency pixel effect in instantaneous products, we used a more sophisticated validating scheme in order to reduce the impacts from these effects as much as possible. Evidences indicate that the GLASS-DSSR products are considerably accurate over most parts of arid and semiarid regions in China, but in complex terrain areas the products might need further refinements. The R² at all sites (except Naqu) was larger than 0.8 with a root mean square error (RMSE) range of about in 90–130 W/m². Linear regression analyses suggest that GLASS-DSSR products tend to overestimate DSSR in the interval of low surface-measured values and symmetrically underestimate DSSR in the interval of high values. This systematic error may result from inappropriate assumptions about clouds and aerosol loadings over the regions in the operational algorithm.     

Influence of the atmospheric masses on the gravitational field of the earth

Seasonal and latitude dependent corrections to the gravity and height anomalies are developed in order to account for the neglect of the atmospheric masses outside the geold, when using Stokes’ equation. It is shown that the atmospheric correction to gravity at sea level is almost constant, equal to0.871 mgals with a variation of2 μ gals whereas the height anomaly correction varies between −0.1 cm and −1.3 cm. Further, when the combined latitudinal/seasonal dependence is neglected in the atmospheric corrections, the maximum error introduced is of the order of40 μ gals for the gravity corrections and0.7 cm for the height anomaly corrections.     

Interpretation of the tidal residuals during the 11 July 1991 total solar eclipse

Observations of gravity and atmospheric pressure variations during the total solar eclipse of 11 July 1991 in Mexico City are presented. An LCR-G402 gravimeter equipped with a feedback system and a digital data acquisition system scanned gravity and pressure every second around the totality. On the pressure record an oscillation, starting at the totality, with a peak to peak amplitude of 0.5 hPa and a periodicity of 40 to 50 min, can clearly be seen. This oscillation results from the thermal shock wave produced by the Moon shadow travelling at supersonic speed. At the 0.1 μGal (1 nm · s⁻²) level all gravity perturbations are explained by the atmospheric pressure effect. Received: 10 February 1995 / Accepted: 7 June 1996     

Determination of Geopotential of Local Vertical Datum Surface

1　ＩｎｔｒｏｄｕｃｔｉｏｎＥａｃｈｃｏｕｎｔｒｙｏｒｅａｃｈｇｒｏｕｐｏｆｃｏｕｎｔｒｉｅｓｓｅｌｅｃｔｓｍｅａｎｓｅａｌｅｖ ｅｌａｔａｄｅｆｉｎｅｄｔｉｄｅｇａｕｇｅｏｒａｔａｇｒｏｕｐｏｆｇａｕｇｅｓｆｏｒｉｔｓｖｅｒｔｉ ｃａｌｄａｔｕｍｓｕｒｆａｃｅ .Ｉｔｉｓｒｅａｌｉｚｅｄ ,ｈｏｗｅｖｅｒ,ｔｈａｔｔｈｅｌｏｃａｌｍｅａｎｓｅａｌｅｖｅｌｉｓｕｓｕａｌｌｙｄｅｐａｒｔｅｄｆｒｏｍｔｈｅｇｅｏｉｄ ,ｗｈｉｃｈｓｈｏｕｌｄｂｅｔｈｅｉｄｅａｌｄａｔｕｍｓｕｒｆａｃｅｆｏｒｈｅｉｇｈｔ,ｏｗ…     

Calibration of an optical equation to analyse the atmospheric turbidity and water quality of an estuarine environment

Mapping of oceanic productive zones through an optical sensor on board a satellite requires a detailed knowledge about the interaction of solar flux with the atmosphere and water column. Hence the radiative transfer on February 14, 2002 over estuaries of Goa has been examined through a couple radiative transfer model (CRTM). This has been carried out in the context of an optical sensor when it looks down at nadir. The forcing of extraterrestrial solar irradiance with atmosphere results in its absorption and scattering due to aerosol particles and gaseous constituents and rayleigh particles. The optical spectrum has been subjected to atmospheric forcing and radiation has been calculated for every 1 nm. Special reference is given to two wavelengths (443 and 510 nm), which are important for ocean optics. The atmospheric components pertaining to these wavelengths are 0.04193 w/m²/nm/sr and 0.0225 w/m²/nm/sr (rayleigh) and 0.0121 w/m²/nm/sr and 0.015 w/m²/nm/sr (aerosol). Within the optical frame work (400 –700 nm), the effect of aerosol overrides that of rayleigh from 545 nm onwards. The specular reflectance at the sea surface contributes 1.725 x l0^-12 w/m²/nm/sr and 1.947x 10^-12 w/m²/nm/sr while the respective signature from water column is 0.418x 10^-3w/m²/nm/sr and 0.536x 10^-3 w/m²/nm/sr. The simulated water signature has been compared with radiances measured by a radiometer (Satlantic radiometer). The correlation (R²) is found to be around 0.97. The sensitivity analysis further reveals the respective wavelengths of the optical spectrum affected by atmospheric and oceanic constituents. Moreover, the effect of change in solar zenith angle on water signature has also been carried out. This could be applied to invert satellite data so as to retrieve coastal water constituents and productivity.     

Spatial and Temporal Variability of Net Primary Productivity (NPP) over Terrestrial Biosphere of India Using NOAA-AVHRR Based GloPEM Model

The monitoring of terrestrial carbon dynamics is important in studies related with global climate change. This paper presents results of the inter-annual variability of Net Primary Productivity (NPP) from 1981 to 2000 derived using observations from NOAA-AVHRR data using Global Production Efficiency Model (GloPEM). The GloPEM model is based on physiological principles and uses the production efficiency concept, in which the canopy absorption of photosynthetically active radiation (APAR) is used with a conversion “efficiency” to estimate Gross Primary Production (GPP). NPP derived from GloPEM model over India showed maximum NPP about 3,000 gCm⁻²year⁻¹ in west Bengal and lowest up to 500 gCm⁻²year⁻¹ in Rajasthan. The India averaged NPP varied from 1,084.7 gCm⁻²year⁻¹ to 1,390.8 gCm⁻²year⁻¹ in the corresponding years of 1983 and 1998 respectively. The regression analysis of the 20 year NPP variability showed significant increase in NPP over India (r = 0.7, F = 17.53, p < 0.001). The mean rate of increase was observed as 10.43 gCm⁻²year⁻¹. Carbon fixation ability of terrestrial ecosystem of India is increasing with rate of 34.3 TgC annually (t = 4.18, p < 0.001). The estimated net carbon fixation over Indian landmass ranged from 3.56 PgC (in 1983) to 4.57 PgC (in 1998). Grid level temporal correlation analysis showed that agricultural regions are the source of increase in terrestrial NPP of India. Parts of forest regions (Himalayan in Nepal, north east India) are relatively less influenced over the study period and showed lower or negative correlation (trend). Finding of the study would provide valuable input in understanding the global change associated with vegetation activities as a sink for atmospheric carbon dioxide.     

Incorporating dynamic factors for improving a GIS‐based solar radiation model

Solar radiation has been a major input to agricultural, hydrological, and ecological modeling. However, solar radiation is usually influenced by three groups of dynamic factors: sun–earth position, terrain, and atmospheric effects. Therefore, an integrated approach to accurately consider the impacts of those dynamic factors on solar radiation is essential to estimate solar radiation over rugged terrain. In this study, a spatial and temporal gap‐filling algorithm was proposed to obtain a seamless daily MODIS albedo dataset. A 1 km‐resolution digital elevation model was used to model the impact of local topography and shading by surrounding terrain on solar radiation. A sunshine‐based model was adopted to simulate radiation under the influence of clouds. A GIS‐based solar radiation model that incorporates albedo, shading by surrounding terrain, and variations in cloudiness was used to address the spatial variability of these factors in mountainous terrain. Compared with other independent solar radiation products, our model generated a more reliable solar radiation product over rugged terrain, with an R² of 0.88 and an RMSE of 2.55 MJ m^?2 day^?1. The improved solar radiation products and open source app can be used further in practice or scientific research.     

A New GIS-based Solar Radiation Model and Its Application to Photovoltaic Assessments

The solar radiation model r.sun is a flexible and efficient tool for the estimation of solar radiation for clear‐sky and overcast atmospheric conditions. In contrast to other models, r.sun considers all relevant input parameters as spatially distributed entities to enable computations for large areas with complex terrain. Conceptually the model is based on equations published in the European Solar Radiation Atlas (ESRA). The r.sun model was applied to estimate the solar potential for photovoltaic systems in Central and Eastern Europe. The overcast radiation was computed from clear‐sky values and a clear‐sky index. The raster map of the clear‐sky index was computed using a multivariate interpolation method to account for terrain effects, with interpolation parameters optimized using a cross‐validation technique. The incorporation of terrain data improved the radiation estimates in terms of the model's predictive error and the spatial pattern of the model outputs. Comparing the results of r.sun with the ESRA database demonstrates that integration of the solar radiation model and the spatial interpolation tools in a GIS can be especially helpful for data at higher resolutions and in regions with a lack of ground measurements.     

Modeling the Global Solar Radiation Under Cloudy Sky Using Meteosat Second Generation High Resolution Visible Raw Data

This paper presents an improved model for mapping both hourly and daily values of the global radiation reaching the earth’s surface from satellite data over Algeria. An empirical model with a large spectral band using Linke atmospheric turbidity was chosen to modeling the clear-sky global radiation. The extinction of the solar radiation (absorption plus scattering) is described by an exponential decay function of the cloud albedo and cloud index. To compute these two parameter, we resolve the sun-pixel-satellite radiative model equation using the following hypothesizes. The first one supposes the total transmittance of the atmospheric in the pixel-satellite direction as the ratio of the direct radiation in real sky to the same quantity in case of a very clean sky. The second one estimates that the noise energy reflected back towards the radiometer by molecules of gases in the upper levels of the atmosphere represents about 5 % of the energy measured by the radiometer. The model performance was validated for five radiometric stations across the country for three selected months in 2013 and satisfactory results were obtained and discussed.