Three dimensional simulation of fluid flow in X‐ray CT images of porous media

A numerical scheme is developed in order to simulate fluid flow in three dimensional (3‐D) microstructures. The governing equations for steady incompressible flow are solved using the semi‐implicit method for pressure‐linked equations (SIMPLE) finite difference scheme within a non‐staggered grid system that represents the 3‐D microstructure. This system allows solving the governing equations using only one computational cell. The numerical scheme is verified through simulating fluid flow in idealized 3‐D microstructures with known closed form solutions for permeability. The numerical factors affecting the solution in terms of convergence and accuracy are also discussed. These factors include the resolution of the analysed microstructure and the truncation criterion. Fluid flow in 2‐D X‐ray computed tomography (CT) images of real porous media microstructure is also simulated using this numerical model. These real microstructures include field cores of asphalt mixes, laboratory linear kneading compactor (LKC) specimens, and laboratory Superpave gyratory compactor (SGC) specimens. The numerical results for the permeability of the real microstructures are compared with the results from closed form solutions. Copyright © 2004 John Wiley & Sons, Ltd.     

NORTHERN AND SOUTHERN HEMISPHERIC POLAR LOWS——A COMPARATIVE STUDY

Northern Hemispheric“polar lows”,i.e.maritime,small,but fairly intense cyclonic systems have been studied formore than two decades.Recently researchers have turned their attention towards the occurrence of polar lows and othersmall-scale cyclonic disturbances on the Southern Hemisphere.In this paper a short account of our present knowledgeof Northern Hemispheric polar lows will be given,followed by some preliminary results of an investigation of the natureand occurrence of polar lows on the Southern Hemisphere.     

江苏岸外辐射沙脊群东沙稳定性研究

东沙是江苏岸外辐射沙脊群中的第二大沙洲,具有独特的地形地貌和水动力条件,对它进行稳定性研究为揭示整个辐射沙洲及其邻近岸滩的动态演变都非常有益。通过利用多年遥感卫片资料、1998年取得的现场水文泥沙观测资料和东沙滩面表层沉积物资料等,对东沙的地形地貌特征、沉积特征和东沙两侧潮汐通道的水流泥沙特征等进行了详细分析。研究结果表明,东沙的沙脊偏于西侧,西侧滩面较窄、高程较高且岸线较为顺直,东侧滩面较宽、高程较低且岸线较为破碎;西洋和陈家坞槽均处于冲刷状态,净输沙的主要方向为输向槽外或输向条子泥;东沙近三十年来面积有所缩小且有外围向中央收缩的趋势,尤其以向东、向南方向的迁移最为明显。     

西太平洋暖池区热带波动的天气效应

为了探讨西太平洋暧地区热带波动的天气效应，利用1980年2－9月140°E日本静止卫星纬度时间剖面卫星云图，分析了5d和14d左右周期热带波动云的演变特征，井综合分析了14d周期的云系演变型式与流场的关系，为预报热带中期天气变化规律提供了依据；热带波动中30－60d大气低频振荡的云量变化最显著，北半球夏半年热带波动的天气压主要在10°N－0°，各半年在10°N－10°S，超过这个区域热带云量的港分布型式就有明显的变化。     

GNSS three carrier ambiguity resolution using ionosphere-reduced virtual signals

This paper presents a general modeling strategy for ambiguity resolution (AR) and position estimation (PE) using three or more phase-based ranging signals from a global navigation satellite system (GNSS). The proposed strategy will identify three best “virtual” signals to allow for more reliable AR under certain observational conditions characterized by ionospheric and tropospheric delay variability, level of phase noise and orbit accuracy. The selected virtual signals suffer from minimal or relatively low ionospheric effects, and thus are known as ionosphere-reduced virtual signals. As a result, the ionospheric parameters in the geometry-based observational models can be eliminated for long baselines, typically those of length tens to hundreds of kilometres. The proposed modeling comprises three major steps. Step 1 is the geometry-free determination of the extra-widelane (EWL) formed between the two closest L-band carrier measurements, directly from the two corresponding code measurements. Step 2 forms the second EWL signal and resolves the integer ambiguity with a geometry-based estimator alone or together with the first EWL. This is followed by a procedure to correct for the first-order ionospheric delay using the two ambiguity-fixed widelane (WL) signals derived from the integer-fixed EWL signals. Step 3 finds an independent narrow-lane (NL) signal, which is used together with a refined WL to resolve NL ambiguity with geometry-based integer estimation and search algorithms. As a result, the above two AR processes performed with WL/NL and EWL/WL signals respectively, either in sequence or in parallel, can support real time kinematic (RTK) positioning over baselines of tens to hundreds of kilometres, thus enabling centimetre-to-decimentre positioning at the local, regional and even global scales in the future.     

Permanent Scatterers for landslide investigations: outcomes from the ESA-SLAM project

Within the SLAM project (Service for Landslide Monitoring), launched in 2003 by the European Space Agency (ESA) the Permanent Scatterers (PS) technique, a multi-image interferometric approach, coupled with the interpretation of aerial-photos and optical satellite images, was carried out for landslide investigations. The PS analysis was applied at a regional scale as support for landslide inventory mapping and at local scale for the monitoring of single well-known slope movements. For the integration of the PS measurements within a landslide inventory the Arno river basin (Italy) was chosen as test site for the presence of a high number of mass movements (to date about 300 areas at high landslide risk and more than 27,000 individual landslides mapped by the institutional authorities). About 350 SAR images have been interferometrically processed by means of the PS technique, with the detection of about 600,000 PS. The use of optical images contributed spatial meaning to the point-wise information provided by the PS, making it easier to identify terrain features related to slope instability and the landslide boundaries. Here we describe the employed methodology and its impact in the updating of a preexisting landslide inventory. 6.8% of the total number of landslides were characterized by ground displacement measurements from the PS: 6.1% of already mapped landslides and 0.8% of new unstable areas detected through the PS analysis. Moreover, most of the PS are located in urban areas, showing that the proposed methodology is suitable for landslide mapping in areas with a quite high density of urbanization, but that over vegetated areas it still suffers from the limitations induced by the current space-borne SAR missions (e.g. temporal de-correlation). On the other hand, the use of InSAR for the monitoring of single slow landslides threatening built-up areas has provided satisfactory results, allowing the measurement of superficial deformations with high accuracy on the landslide sectors characterized by a good radar reflectivity and coherence.     

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.