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.     

Wavelet Analysis of Space Solar Telescope Images

The scientific satellite SST (Space Solar Telescope) is an important research project strongly supported by the Chinese Academy of Sciences. Every day, SST acquires 50 GB of data (after processing) but only 10GB can be transmitted to the ground because of limited time of satellite passage and limited channel volume. Therefore, the data must be compressed before transmission. Wavelets analysis is a new technique developed over the last 10 years, with great potential of application. We start with a brief introduction to the essential principles of wavelet analysis, and then describe the main idea of embedded zerotree wavelet coding, used for compressing the SST images. The results show that this coding is adequate for the job.     

SAMPLING REMOTELY SENSED IMAGERY FOR STORAGE,RETRIEVAL, AND RECONSTRUCTION*

A major difficulty in remote sensing is handling the many data from sensors aboard aircraft and satellites. In this paper we identify an optimal procedure for sampling remotely sensed data before their storage or on their retrieval. The procedure depends on spatial correlation in the scene and uses kriging to estimate values that have been lost. An example in which data from an airborne multispectral scanner could be diminished to only about one tenth without serious loss of precision illustrates the method.     

Surface Photometry of Comet Hale-Bopp Carbon Images

A very long series of photographic observations of the comet Hale-Bopp has been made during January–April 1997 at the double astrograph (400/2000) of the Main Astronomical Observatory (Kyiv, Ukraine). Some of the cometary photos were obtained with two wide-band filter combinations. One of these combinations isolates C₂ emission, another — the nearby dust continuum. The images were digitized by means of AMDPH-XY machine and then calibrated following the standard procedure. After subtraction of the dust continuum the distribution of surface brightness in the C₂ emission coma of comet Hale-Bopp was studied. We found an asymmetric brightness distribution both pre- and post-perihelion. On 21.77 April 1997 a secondary brightness peak is found at the distance of 1.03 × 10⁵ km from the nucleus. It is possible that this peak is related to the extended source of the C₂ molecules. This revised version was published online in July 2006 with corrections to the Cover Date.     

Post-solidification cooling and the age of Io's lava flows

The modeling of thermal emission from active lava flows must account for the cooling of the lava after solidification. Models of lava cooling applied to data collected by the Galileo spacecraft have, until now, not taken this into consideration. This is a flaw as lava flows on Io are thought to be relatively thin with a range in thickness from ∼1 to 13 m. Once a flow is completely solidified (a rapid process on a geological time scale), the surface cools faster than the surface of a partially molten flow. Cooling via the base of the lava flow is also important and accelerates the solidification of the flow compared to the rate for the ‘semi-infinite’ approximation (which is only valid for very deep lava bodies). We introduce a new model which incorporates the solidification and basal cooling features. This model gives a superior reproduction of the cooling of the 1997 Pillan lava flows on Io observed by the Galileo spacecraft. We also use the new model to determine what observations are necessary to constrain lava emplacement style at Loki Patera. Flows exhibit different cooling profiles from that expected from a lava lake. We model cooling with a finite-element code and make quantitative predictions for the behavior of lava flows and other lava bodies that can be tested against observations both on Io and Earth. For example, a 10-m-thick ultramafic flow, like those emplaced at Pillan Patera in 1997, solidifies in ∼450 days (at which point the surface temperature has cooled to ∼280 K) and takes another 390 days to cool to 249 K. Observations over a sufficient period of time reveal divergent cooling trends for different lava bodies [examples: lava flows and lava lakes have different cooling trends after the flow has solidified (flows cool faster)]. Thin flows solidify and cool faster than flows of greater thickness. The model can therefore be used as a diagnostic tool for constraining possible emplacement mechanisms and compositions of bodies of lava in remote-sensing data.     

Direct mid-IR images of Gl 2591,W51-IRS2 and S140

We present the results of direct mid-IR images of the luminous young stellar objects GL 2591, S140-IRS1, and W51-IRS2. The sources show an extended mid-IR emission at the limit of our spatial resolution indicating the presence of dense circumstellar dust disks.     

Occurrence and evolution of the Xiaotangshan hot spring in Beijing, China

Thermal groundwater occurs in bedrock aquifers consisting of the dolomite of the Wumishan Group of the Jixianin System and the Cambrian carbonate in the Xiaotangshan geothermal field near the northern margin of the North China Plain, China. The hot water in the geothermal field of basin-type discharges partly in the form of the Xiaotangshan hot spring under natural conditions. The hot water has TDS of less than 600 mg/L and is of Na·Ca-HCO₃ type. The geothermal water receives recharge from precipitation in the mountain area with elevation of about 500 m above sea level to the north of the spring. Thermal groundwater flows slowly south and southeast through a deep circulation with a residence time of 224 years estimated with the Ra–Rn method. The Xiaotangshan hot spring dried up in the middle of the 1980s owing to the increasing withdrawal of the hot water in the geothermal field in the past decades. The water level of the geothermal system still falls continually at an annual average rate of about 2 m, although water temperature changes very little, indicating that the recharge of such a geothermal system of basin-type is limited. Over-exploitation has a dramatic impact on the geothermal system, and reduction in exploitation and reinjection are required for the sustainable usage of the hot water.     

延边东部五道沟岩群的单颗粒锆石SHRIMP U-Pb 年代学及其地质意义

通过对延边地区东部五道沟岩群黑云阳起石片岩的单颗粒锆石SHRIMP U-Pb年龄测定.获得21个单点年龄,其中包括7组谐和年龄和1组不谐和年龄.谐和年龄分别为1347.8 Ma、844.8 Ma、340.2～313.7 Ma(平均值为323±23 Ma,N=4,MSWD=0.23,置信度=0.88)、292.9～288.3 Ma(平均值=291±25 Ma,N=3,MSWD=0.031.置信度=0.86)、279.2～266.2 Ma(平均值为279±28 Ma,N=4,MSWD=0.031,置信度=0.86)、127.4～124.2 Ma(206Pb/208U年龄平均值为126.5±3.7 Ma,N=5,MSWD=0.12;置信度=O.97)、116.1～106.3 Ma(平均值为115±39 Ma,N=2,MSWD=1.2置信度=0.27),不谐和年龄的下交点年龄为(451±120)Ma、上交点年龄为(1811±400)Ma(MSWD=7.2);这一结果表明:阳起石片岩的原岩主要是来自中元古代、新元古代和早古生代的碎屑物,指示五道沟群的沉积成岩作用发生在石炭世(323±23 Ma),变质作用发生在晚二叠世(291±25 Ma),之后在279.2～266.2 Ma、126.5～106.3 Ma先后受两次岩浆用和蚀变作用的改造.     

Images of the seabed in the Faroe-Shetland Channel from commercial 3D seismic data

A large number of 3D deep seismic surveys in the Faroe-Shetland Channel gives continuous coverage over most of the region. These surveys were designed primarily to image depths in excess of 4 km, use low frequency sources and are recorded at low temporal sample rates. However, commercial 3D data can generate highly detailed images of the seabed due to the high spatial sample rate, typically 12.5 m. This is particularly true in waters below 200 m. Despite geophysical artefacts, the images reveal that there are a number of sedimentary processes at work adjacent to and within this channel. On the West Shetland Shelf, iceberg scouring and moraines reflect the impact of glaciation. On the West Shetland slope there is clear evidence for down-slope processes, such as debris flows, linear erosion channels, basal fans and (one case) slope failure. Along-slope processes are also active as indicated by the presence of sediment waves and contourite mounds. On the floor of the basin, polygonal cracking can be observed. The most spectacular feature appears to be the Judd Deeps, a system of cliffs approximately 200 m high and 40 km across. Traditionally, seabed investigation has been performed using high-resolution surveys. This study shows that deep exploration data can also provide useful images of the seafloor.     

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

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