Improving the assessment of ICESat water altimetry accuracy accounting for autocorrelation

Given that water resources are scarce and are strained by competing demands, it has become crucial to develop and improve techniques to observe the temporal and spatial variations in the inland water volume. Due to the lack of data and the heterogeneity of water level stations, remote sensing, and especially altimetry from space, appear as complementary techniques for water level monitoring. In addition to spatial resolution and sampling rates in space or time, one of the most relevant criteria for satellite altimetry on inland water is the accuracy of the elevation data. Here, the accuracy of ICESat LIDAR altimetry product is assessed over the Great Lakes in North America. The accuracy assessment method used in this paper emphasizes on autocorrelation in high temporal frequency ICESat measurements. It also considers uncertainties resulting from both in situ lake level reference data. A probabilistic upscaling process was developed. This process is based on several successive ICESat shots averaged in a spatial transect accounting for autocorrelation between successive shots. The method also applies pre-processing of the ICESat data with saturation correction of ICESat waveforms, spatial filtering to avoid measurement disturbance from the land–water transition effects on waveform saturation and data selection to avoid trends in water elevations across space. Initially this paper analyzes 237 collected ICESat transects, consistent with the available hydrometric ground stations for four of the Great Lakes. By adapting a geostatistical framework, a high frequency autocorrelation between successive shot elevation values was observed and then modeled for 45% of the 237 transects. The modeled autocorrelation was therefore used to estimate water elevations at the transect scale and the resulting uncertainty for the 117 transects without trend. This uncertainty was 8 times greater than the usual computed uncertainty, when no temporal correlation is taken into account. This temporal correlation, corresponding to approximately 11 consecutive ICESat shots, could be linked to low transmitted ICESat GLAS energy and to poor weather conditions. Assuming Gaussian uncertainties for both reference data and ICESat data upscaled at the transect scale, we derived GLAS deviations statistics by averaging the results at station and lake scales. An overall bias of −4.6 cm (underestimation) and an overall standard deviation of 11.6 cm were computed for all lakes. Results demonstrated the relevance of taking autocorrelation into account in satellite data uncertainty assesment.     

Observation on height changes of chinese tide gauges by GPS

IntroductionBecause of global temperature rising and glacierthawing,the sea level has risen about 10-25 cmin the past century[1]. The sea level change isone of the main ai ms in global change monito-ring. Presently the basic tool used in sea levelmonitori…     

星载加速度计的动力法校准

通过模拟轨道,利用动力法对星载加速度计校准进行了研究。用不同参考重力场模型获得的校准参数分别重新积分卫星轨道,与模拟轨道的差异最大可达m级。结果表明,参考重力场模型误差对加速度计校准的影响不可忽略,在校准加速度计的同时应解算重力场模型,以削弱模型误差的影响。     

导航用海洋重力异常图的孔斯曲面重构方法

本文基于地球重力物理场连续的特性,将计算机构图中的孔斯（Coons）曲面建模引入到导航用海洋重力异常图的加密重构中,建立了不同边界曲线的双一次Coons曲面和双三次C1Coons曲面重力异常模型。通过对实验区数据计算分析,结果表明基于Coons曲面建立的重力异常模型精度均高于双线性内插模型,且边界条件数越多,模型精度越高;同时利用现有全球SRTM30plus海深数据,借鉴于间接内插的思想,以布格异常基于Coons曲面进行重力异常图的重构,与空间重力异常直接建模相比,精度整体提高了25%。     

基于卫星加速度恢复地球重力场的去相关滤波法

基于加速度法恢复地球重力场时,卫星加速度是由卫星轨道数值微分得到,而数值微分会放大高频误差,进而降低了重力场解算结果的精度.针对数值微分导出的加速度误差具有有色噪声的特性,提出利用去相关算法构造白化滤波器对加速度有色噪声进行滤波处理,并根据去相关的基本原理分别构造了基于三点差分和ARMA模型的白化滤波器.采用不同噪声背景的CHAMP卫星模拟轨道数据进行解算,结果表明:基于去相关滤波解算的重力场模型精度均要比等权解算的重力场模型精度高,初步验证了去相关滤波方法的有效性.     

利用引力梯度数据计算径向梯度的优化方法

根据重力梯度观测各分量的方差及协方差信息,提出了利用GOCE梯度数据计算径向重力梯度的优化方法。首先给出了径向重力梯度的计算方法,并深入分析了误差传播规律,通过建立相应的条件极值问题,给出了计算径向重力梯度最优组合因子的方法;通过模拟数据验证了本文所提出的优化因子的优越性。实际数据计算表明:相对于传统方法,采用优化组合因子可使反演所得引力位模型的累积大地水准面精度在250阶时提高约2 cm。由于径向重力梯度不仅可以用于地球引力场模型的求解,也可直接应用于地球物理问题的讨论,因此本文所提出的优化方法也可对部分地球动力学问题的讨论提供方便。     

黑龙江省耕地非农化的空间格局与重心曲线分析

为了揭示耕地非农化的空间格局与重心曲线,为耕地资源保护和利用政策制定提供决策支撑,该文以黑龙江省1980、1990、2000、2010、2015年土地利用数据为基础,运用GIS空间分析方法,引入重心模型分析黑龙江省耕地非农化的空间格局及重心曲线。研究结果表明:1980—2015年黑龙江省耕地非农化发展不平衡,4个时期耕地非农化面积呈波动增长的趋势;耕地非农化等级空间分布特征明显,等级较高的城市集中分布在西部地区;耕地非农化速度等级的重心曲线在空间分布上,具有从中心向边缘延伸的特点;耕地非农化围绕中心地区进行,其空间扩散路径在图形上呈现出一个回旋状。     

基于DMSP-OLS和珞珈一号夜光遥感影像的粤港澳大湾区城市空间格局演变

夜光遥感影像数据可有效反映城市空间格局变化。本文基于1992—2012年的DSMP-OLS夜光遥感影像和2018年的珞珈一号遥感影像,利用分层阈值法提取粤港澳大湾区内各城市建成区;通过计算平均灯光强度、平均灯光增长速率、城市建成区面积、城市建成区增长速率、城市重心、城市重心偏移距离等一系列指数,揭示区内各城市的空间格局演变过程。研究结果表明:①1992—2018年,粤港澳大湾区的城市规模大幅增长,沿珠江口两侧形成了以澳门、广州、深圳和香港为核心的倒“U”形城市群,并呈辐射状向周边扩张。②以珠江口为界,粤港澳大湾区东部各个城市的发展水平整体高于西部各个城市,广州、深圳、香港等核心城市发展水平明显高于江门、肇庆、惠州等外围城市。③1992—2018年,粤港澳大湾区建成区的增长速率由小变大,最后逐渐趋于稳定,2002—2007年是城市扩张最迅猛时期。④1992—2018年,粤港澳大湾区的各城市重心迁移方式表现为3种类型:持续向区域中心迁移;持续向相邻城市邻接区迁移;持续向海洋方向迁移。大部分城市的重心迁移方向呈“震荡”特征。     

Spaceborne gravitational field determination by means of locally supported wavelets

In space-borne gravitational field determination, two challenges are inherent. First, the continuation of the data down to the surface of the Earth is an ill-posed problem, requiring therefore regularization techniques. Second huge data sets result requiring efficient numerical methods. In this paper, we show how locally supported wavelets on the sphere can be developed by means of a spherical version of the so-called up function. By construction, the corresponding scaling functions and wavelets are infinitely smooth, so that they can be used for regularization purposes. In particular, we show how the ill-posed pseudo-differential equations coming from satellite missions can be regularized by efficient numerical schemes using locally supported wavelets. These methods seem in particular to be interesting for regional gravity field modelling.     

Ellipsoidal geoid computation

Modern geoid computation uses a global gravity model, such as EGM96, as a third component in a remove–restore process. The classical approach uses only two: the reference ellipsoid and a geometrical model representing the topography. The rationale for all three components is reviewed, drawing attention to the much smaller precision now needed when transforming residual gravity anomalies. It is shown that all ellipsoidal effects needed for geoid computation with millimetric accuracy are automatically included provided that the free air anomaly and geoid are calculated correctly from the global model. Both must be consistent with an ellipsoidal Earth and with the treatment of observed gravity data. Further ellipsoidal corrections are then negligible. Precise formulae are developed for the geoid height and the free air anomaly using a global gravity model, given as spherical harmonic coefficients. Although only linear in the anomalous potential, these formulae are otherwise exact for an ellipsoidal reference Earth—they involve closed analytical functions of the eccentricity (and the Earths spin rate), rather than a truncated power series in e². They are evaluated using EGM96 and give ellipsoidal corrections to the conventional free air anomaly ranging from –0.84 to +1.14 mGal, both extremes occurring in Tibet. The geoid error corresponding to these differences is dominated by longer wavelengths, so extrema occur elsewhere, rising to +766 mm south of India and falling to –594 mm over New Guinea. At short wavelengths, the difference between ellipsoidal corrections based only on EGM96 and those derived from detailed local gravity data for the North Sea geoid GEONZ97 has a standard deviation of only 3.3 mm. However, the long-wavelength components missed by the local computation reach 300 mm and have a significant slope. In Australia, for example, such a slope would amount to a 600-mm rise from Perth to Cairns.