Geometry and accuracy of reflecting points in bistatic satellite altimetry

The analysis of the time and space distribution of specular (reflecting) points in bistatic altimetry between GPS and CHAMP satellites or SAC-C (taken as examples) is extended from Wagner and Klokočník (2003 J. Geod 77: 128–138). We demonstrate a significantly higher number and density of reflecting points in bistatic altimetry in comparison with traditional monostatic altimetry. After an outline of our older accuracy assessment for the vertical position of the reflecting point, we add a new independent derivation and compare both approaches. We account for orbit errors of both the transmitters (GPS) and receiver (CHAMP) satellites, and the measurement (delay) error. We found that the accuracy of the vertical position of the reflecting point decreases only slowly with increasing off-nadir angle and that the orbit errors must be accounted for if decimeter and better accuracy is required. In this paper, we do not study errors such as state of the ocean, technical parameters of the receiving system, and atmospheric corrections.     

Low-low satellite-to-satellite tracking: a comparison between analytical linear orbit perturbation theory and numerical integration

Low-low satellite-to-satellite tracking (ll-SST) range-rate observations have been predicted by two methods: one based on a linear perturbation theory in combination with the Hill equations, and one based on solving the equations of motion of two low-flying satellites by numerical integration. The two methods produce almost equivalent Fourier spectra of the range-rate observations after properly taking into account a few resonant terms. For a typical GRACE-type configuration, where the two satellites trail each other at a distance of 300 km at an altitude of 460 km, and in the presence of the EGM96 gravity field model, complete to spherical harmonic degree and order 70, the agreement between the Fourier spectra is about 1 mm/s compared to a root-mean-square (RMS) value of more than 220 mm/s for the range-rate signal. The discrepancy of 1 mm/s can be reduced significantly when not taking into account perturbations caused by the J₂ term. Excluding the J₂ term, the agreement between the two methods improves to 0.4 mm/s compared to a RMS value of 6 mm/s for the range-rate signal. These values are 0.01 and 2.3 mm/s when ignoring the spectrum for frequencies below two cycles per orbital revolution, reducing the discrepancy even further to about 0.5% of the signal. The selected linear perturbation theory is thus capable of modeling gravity field induced range-rate observations with very high precision for a large part of the spectrum.     

Outlier detection algorithms and their performance in GOCE gravity field processing

The satellite missions CHAMP, GRACE, and GOCE mark the beginning of a new era in gravity field determination and modeling. They provide unique models of the global stationary gravity field and its variation in time. Due to inevitable measurement errors, sophisticated pre-processing steps have to be applied before further use of the satellite measurements. In the framework of the GOCE mission, this includes outlier detection, absolute calibration and validation of the SGG (satellite gravity gradiometry) measurements, and removal of temporal effects. In general, outliers are defined as observations that appear to be inconsistent with the remainder of the data set. One goal is to evaluate the effect of additive, innovative and bulk outliers on the estimates of the spherical harmonic coefficients. It can be shown that even a small number of undetected outliers (<0.2 of all data points) can have an adverse effect on the coefficient estimates. Consequently, concepts for the identification and removal of outliers have to be developed. Novel outlier detection algorithms are derived and statistical methods are presented that may be used for this purpose. The methods aim at high outlier identification rates as well as small failure rates. A combined algorithm, based on wavelets and a statistical method, shows best performance with an identification rate of about 99%. To further reduce the influence of undetected outliers, an outlier detection algorithm is implemented inside the gravity field solver (the Quick-Look Gravity Field Analysis tool was used). This results in spherical harmonic coefficient estimates that are of similar quality to those obtained without outliers in the input data.     

World-wide synthetic tide parameters for gravity and vertical and horizontal displacements

The response of the Earth’s crust to the direct effect of lunisolar gravitational forcing is known as the body tide. The body tide is superimposed by surface-loading forces due to the pressure of the periodically varying ocean tide acting on the Earth, called ocean tide loading (OTL). Both body tide and OTL can be decomposed into components of the same frequency known as tidal parameters. However, OTL is more complicated than body tides because of the dynamic effects of the ocean. Estimating OTL requires a model of the ocean tides and knowledge of the elastic properties of the solid Earth. Thus, synthetic tide parameters (amplitude factors and phase leads) have been developed here on a world-wide grid for gravity and positional displacements. The body tide contributions were added to the oceanic contribution to provide the Earth tide response. The accuracy and reliability of the synthetic tidal parameters have been estimated by comparing observed gravity and vertical-displacement tide parameters with those interpolated from our synthetic model, which shows good agreement. Tests also indicate that the synthetic tide parameters provide realistic gravimetric and displacements for practical use in tidal prediction.     

蓖麻生长农业气象条件分析及栽培技术初探

根据对蓖麻发育期的观测，结合当地气象条件，分析了阿克苏种植蓖麻的适应性及栽培中应采取的趋利避害的技术措施。     

淮河流域土壤湿度异常的时空分布特征及其与气候异常关系的初步研究

利用中国东部1990~2000年旬平均土壤湿度、降水和气温观测资料，通过对0~50 cm层次土壤湿度进行旋转主分量分析 (REOF)，重点分析了淮河流域土壤湿度的时空分布特征, 并初步研究了土壤湿度与前期、同期和后期不同时段降水与气温的关系。发现春季以30 cm为界，30 cm以上各层土壤湿度异常的第一旋转空间模态十分相似, 其大值中心主要位于淮河流域，而30 cm以下 (30~50 cm) 各层的第二旋转空间模态与之亦十分类似, 因此称该模态为“淮河型”，而夏季和秋季虽然该模态也很显著, 但特征不如春季突出。该模态在各层次土壤中具有明显的持续性特征，均存在40旬左右的显著周期；并与前期和同期降水（气温）呈显著正 (负) 相关关系，与约半年后的降水 (气温) 呈负 (正) 相关关系。     

上海“8.5”特大暴雨的成因和特点

利用WSR－88D多普勒雷达资料，通过分析2001年8月5日到8月6日产生于上海地区的特大暴雨过程的雷达回波特征，研究探讨这次过程的降水成因和降水特点。发现副热带高压边缘偏南气流的风速辐合、螺旋雨带中强回波短带的合并加强、风的垂直切变产生的辐合以及热带系统的高效率降水作用是这次暴雨过程的降水成因和降水特点。并提出了预报这类暴雨系统需要注意的要点：强降水回波的发生和发展加强往往和风场中的局地强辐合区相联系，利用多普勒雷达连续演变的风场资料，可及时了解辐合区的形成和演变过程，提前作出强降水可能性的预测。     

长江三角洲地区水和热通量的时空变化特征及影响因子

文中利用改进的K B模式和牛顿扩散方法及 196 1年以来的长江三角洲 (2 8～ 33°N ,118～ 12 3°E)地区的 4 8个测站的常规气象资料 ,估计了该地区近 4 0a来的蒸散量和感热通量。结合该地区的气温、太阳辐射等气候资料和 196 0年以来该区域土地资源利用变化等有关信息对该地区的潜热通量和感热通量的时 空间变化特征及其可能成因进行了综合分析。结果认为该地区自 2 0世纪 70年代开始平均蒸散量有逐渐减小的趋势 ,与 1980年相比 ,1998年区域年平均蒸散量减小了 2 4mm。区域平均感热通量与蒸散量相比在此期间变化并不明显。通过对该地区的云量、太阳辐射及土地利用变化资料分析认为 ,造成该地区平均蒸散量减少趋势的的原因之一是用于蒸发的能量即太阳辐射的减少 ,而造成太阳辐射减少的可能原因为云量及大气透明度的变化所至 ;原因之二是该地区地表覆盖条件的改变。近 2 0a来 ,该地区的水田、旱地及水域面积占总面积的比率分别减少 1.35 3% ,4 .4 4 2 %和2 .5 97% ,而城镇建设、工矿及其它建设用地面积则增加 3.345 %。耕地及水面的减小和城镇及建设用地面积的增加从整体上使区域平均蒸发量减少。     

波包传播诊断的理论基础和计算方法

文中提出了一种新的利用实际观测资料研究波包传播的诊断方法 (WPD) ,并通过理想资料和实际资料的计算证明了该方法的可用性和广泛性     

冷暖事件对大气能量循环和纬向平均环流影响的模拟研究

利用中国科学院大气物理研究所大气科学与地球流体力学数值模拟国家重点实验室新发展的GOALS 5全球海 陆 气耦合模式研究了暖事件 (ElNi no)和冷事件 (LaNina)对大气能量循环和纬向平均环流的影响 ,并用观测资料进行了对比分析。结果表明 :对于纬向平均资料来说 ,冷、暖事件在热带和副热带地区的大气环流相关量的反相变化特征非常清晰 ,中高纬度地区并不明显。此外 ,还发现 ,暖事件时定常涡动的经向热通量的变化是北半球对流层热带外地区温度异常的主要原因 ,而瞬变波的影响则起抵消作用。冷事件时定常波和瞬变波相互抵消的局地特征也依然存在 ,但瞬变波的影响有所增强。