Estimation of Sea Ice Thickness from SARAL/AltiKa in Drifting Orbit Phase

Abstract

Intra and inter-annual variations in the sea ice thickness are highly sensitive indicators of climatic variations undergoing in the earth’s atmosphere and oceans. This paper describes the method of estimating sea ice thickness using radar waveforms data acquired by SARAL/Altika mission during its drifting orbit phase from July 2016 onwards yielding spatially dense data coverage. Based on statistical analysis of return echoes, classification of the surface has been carried out in three different types, viz. floe, lead and mixed. Time delay correction methods were suitably selected and implemented to make corrections in altimetric range measurements and thereby freeboard. By assuming hydrostatic equilibrium, freeboard data were converted into sea ice thickness. Results show that sea ice thickness varies from 4 to 5?m near ice shelves and 1 to 2.5?m in the marginal sea ice regions. Freeboard and sea ice thickness estimates were also validated using NASA’s Operation Ice Bridge (OIB) datasets. Freeboard measurements show very high correlation (0.97) having RMSE of 0.13. Overestimation of approximately 1–2?m observed in the sea ice thickness, which could be attributed to distance between AltiKa footprint and OIB locations. Moreover, sensitivity analysis shows that snow depth and snow density over sea ice play crucial role in the estimation of sea ice thickness.     

Separation of Waves by Interpolation Method

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

Optimizing interpolation of shoot density data from a Posidonia oceanica seagrass bed

A case study on the optimization of Posidonia oceanica density interpolation, using a data set from a large meadow at Porto Conte Bay (NW Sardinia, Italy), is presented. Ordinary point kriging, cokriging and a weighted average based on inverse square distance were used to interpolate density data measured in 36 sampling stations. The results obtained from different methods were then compared by means of a leave‐one‐out cross‐validation procedure. The scale at which interpolation was carried out was defined on the basis of the Hausdorff dimension of the variogram. Optimizing spatial scale and data points search strategy allowed obtaining more accurate density estimates independently of the interpolation method.     

重力异常垂向一阶导数的一种简便算法

利用拉格朗日插值导出了重力异常垂向一阶导数的计算公式，给出了地面以及地面之上不同高度的求导系数。该公式可以计算地面上的重力垂向一阶导数，还可以直接计算地面之上任意高度上的垂向一阶导数。鉴于该公式除系数不同外，与上延公式完全相同，因此，程序设计尤为简单。使用本公式对模型数据和实际资料进行了处理，证明了本算法的实用性。     

单中继星单低轨卫星的联合定轨精度分析

针对由单中继星和单低轨卫星组成的联合定轨系统,给出了系统内不同轨道卫星摄动项的选取方案和卫星间的可见性判别模型。在模拟出含有白噪声的四程测距观测数据文件的基础上,研究了测距精度和采样弧段对联合定轨中高、低轨卫星定轨精度的影响。得出如下结论:联合定轨更有利于对低轨卫星的轨道改进;同样的采样时间条件下,测距噪声越小定轨精度越高,并且采样时间越短它的影响越明显;同样的测距噪声条件下,所用资料的采样时间越长精度越高,但当测距精度很高时,TDRS达到最好定轨精度所需的采样时间相应有所缩短。并在定出轨道后进行了轨道预报,分析了轨道预报的趋势及精度。     

一种基于启发式算法的等高线局部内插方法

在基于启发式内插等高线算法的基础上提出了一种局部内插方法。首先利用Douglas-Peucker算法提取等高线的特征点,根据特征点判断等高线之间的相似性程度,找出导致等高线出现异常的特征点;然后将相似性程度很低的两条等高线自动分解为若干简单等高线再进行内插。有效地解决了局部弯曲很大、马鞍型地貌等复杂等高线的内插问题。算法已经在以Microstation为平台的数字制图系统中实现并逐渐实用化。     

基于精密星历的计算卫星位置的方法

分析了拉格朗日多项式插值算法,提出了将Neville算法实现于精密星历的插值,数值试验表明,其克服了函数插值的龙格现象,有很好的逼近效果.     

基于钻孔数据的三维地质建模与可视化研究

三维地质模型能够完整准确地表达复杂地质现象的边界条件及地质体内含的各种地质构造,直观地再现地质单元的空间展布及其相互关系,最大限度地提高地质分析的直观性和准确性,因此三维地质建模已引起地质、采矿、岩土工程等诸多领域的日益重视。笔者在综合考察多种建模方法的基础上,结合工程实践,提出了一种实用的工程地质三维建模实施方法。该方法以钻孔资料和地质剖面数据为三维地质建模的源数据,允许用户手工修正模型剖面地层分界线,从而实现对建模过程的干预和控制,克服了单纯依靠钻孔数据建模结果不精确且难以修正的问题。最后,通过一个研究实例展示了该方法的实际建模效果。     

一种新型非线性滤波理论及其应用

研究一种新型的非线性滤波理论,即Unscented卡尔曼滤波（UKF）,同时为了获得更高的计算效率和确保协方差阵的非负定性,研究了平方根UKF。将UKF和平方根UKF应用到星载GPS卫星定轨中,实际算例表明UKF和平方根UKF的性能要优于常用的推广卡尔曼滤波的性能。     

Interpolation by regularized spline with tension: I. Theory and implementation

Bivariate and trivariate functions for interpolation from scattered data are derived. They are constructed by explicit minimization of a general smoothness functional, and they include a tension parameter that controls the character of the interpolation function (e.g., for bivariate case the surface can be tuned from a membrane to a thin steel plate), Tension can be applied also in a chosen direction, for modeling of phenomena with a simple type of anisotropy. The functions have regular derivatives of all orders everywhere. This makes them suitable for analysis of surface geometry and for direct application in models where derivatives are necessary. For processing of large datasets (thousands of data points), which are now common in geosciences, a segmentation algorithm with a flexible size of overlapping neighborhood is presented. Simple examples demonstrating flexibility and accuracy of the functions are presented.On leave from the Department of Physical Geography and Cartography, Comenius University, Mlynská dolina, Bratislava, Czechoslovakia.On leave from the Institute of Physics, Dúbravská cesta 9, Bratislava, Czechoslovakia.