Numerical manifold method for dynamic nonlinear analysis of saturated porous media

It is well known that the Babuska–Brezzi stability criterion or the Zienkiewicz–Taylor patch test precludes the use of the finite elements with the same low order of interpolation for displacement and pore pressure in the nearly incompressible and undrained cases, unless some stabilization techniques are introduced for dynamic analysis of saturated porous medium where coupling occurs between the displacement of solid skeleton and pore pressure. The numerical manifold method (NMM), where the interpolation of displacement and pressure can be determined independently in an element for the solution of u–p formulation, is derived based on triangular mesh for the requirement of high accurate calculations from practical applications in the dynamic analysis of saturated porous materials. The matrices of equilibrium equations for the second‐order displacement and the first‐order pressure manifold method are given in detail for program coding. By close comparison with widely used finite element method, the NMM presents good stability for the coupling problems, particularly in the nearly incompressible and undrained cases. Numerical examples are given to illustrate the validity and stability of the manifold element developed. Copyright © 2006 John Wiley & Sons, Ltd.     

CATASTROPHIC ECO-ENVIRONMENTAL CHANGE IN THE SONGNEN PLAIN, NORTHEASTERN CHINA SINCE 1900S

Although the Songnen Plain in the northeastern China was developed relatively late in the temperate zone of the world, its eco-environment has changed greatly. This paper analyzes the changes of land cover and the rates and trends ofdesertification during the past 100 years in the Songnen Plain. According to the macroscopic analysis, we find that the eco-environment in the plain has reached to the threshold of catastrophic change since the 1950s. The Thorn Needle Catastrophic Model was used to determine and validate this conclusion. Human activities, including large-scale construction projects, such as huge dams and dikes, and excessive grazing were the primary factors contributing to regional eco-environmental catastrophe. And irrational reclamation of the wilderness also affected the eco-environmental change. The results reveal the complex human-land interactions.     

Luminescence dating of sand–loess sequences and response of Mu Us and Otindag sand fields (north China) to climatic changes

The sand–loess transition zone in north China is sensitive to climate change, and is an ideal place to investigate past environmental changes. However, past climate change at millennial–centennial timescales in this region has not been well reconstructed because of limited numerical dating. Alternations of sandy loam soils with aeolian sand layers in the Mu Us and Otindag sand fields, which lie along the sand–loess transition zone, indicate multiple intervals of dune activity and stability. This change is probably a response to variations of the East Asian monsoon climate during the late Quaternary. The single aliquot regeneration (SAR) optically stimulated luminescence (OSL) dating protocol, which has been successfully applied to aeolian deposits worldwide, is applied to these two sand fields in this study. The OSL ages provide reliable constraints for reconstruction of past climate changes at suborbital timescale. Sections in both sand fields contain aeolian sand beds recording millennial‐scale episodes of dry climate and widespread dune activation, including episodes at about the same time as Heinrich Event 5 and the Younger Dryas in the North Atlantic region. These results demonstrate the potential of aeolian sediments in semi‐arid north China to record millennial‐scale climatic events, and also suggest that dry–wet climate variation at the desert margin in China may be linked to climatic change elsewhere in the Northern Hemisphere, through atmospheric circulation. This article was published online on 27 November 2008. An error was subsequently identified. This notice is included in the online and print versions to indicate that both have been corrected (16 December 2008). Copyright © 2008 John Wiley & Sons, Ltd.     

A THREE-DIMENSIONAL ELASTIC NESTED-GRID MESO-(β-γ)SCALE ATMOSPHERIC MODEL——PARTI:MODEL DESCRIPTION

A three-dimensional elastic nonhydrostatic mesoscale(β-γ)model with nested-grid is presented.It uses a set of fullequations in terrain-following coordinates as its basic dynamic frame,which is solved with a time-splitting algorithmfor acoustic and gravity waves.The model physical parameterization includes a K-theory subgrid eddy mixing for cloudand free atmosphere,a bulk planetary boundary layer parameterization,and three types of sofisticated cloudmicrophysics schemes with double-parameters for hail-bearing clouds,warm clouds and snowing clouds respectively.The model is designed to be used flexibly for simulations of a variety of meso-and small-scale atmospheric processes,and can be improved as a regional and local operational NWP system in future.     

Development-Conservation Dilemma in the NilgiriMountains of South India

The Nilgiri Mountains of south India is considered unique by anthropologists, geologists, climatologists, botanists as well as tourists. It has remained a subject of constant study and research over the last two centuries. Man-nature balance had continued undisturbed in the Nilgiris for thousands of years until the early 19th century when it became a British colony attracting, in due course, various developmental activities. Subsequently, the Nilgiris and its popular hill stations emerged as favourite places for the British population in India for rest and recuperation, game and for raising commercial plantations. In the process, the traditional indigenous crops were replaced by “English“ vegetables and the natural forests gave way to commercial plantations of coffee, tea and other exotic species of trees.     

Development and Extension of Seawater Desalination by Reverse Osmosis

Seawater desalination has been people‘s fond dream since ancient times,the dream is now becoming a reality.This paper presents a brief development history of reverse osmosis.Much attention was paid to innovative development in mebranes,modules.equipments and applied technology,including asymmetric and composite membranes,spiral-wound element and hollow fiber module,energy recovery equipments and different technological processes.The extension of reverse osmosis,such as desalination,pre-concentration,integrated processes and nanofiltration,is also briefly mentioned.     

中国石油生产的展顾

本文分析了中国石油生产近20年变化的原因,对一些学者认为80年代中国会成为超级产油大国,后又悲观断定80年代中国可能变成石油净进口国作了评述。同时对90年代中国石油生产前景作了初步展望。     

Relationship between footwall composition,crustal contamination,and fluid–rock interaction in the Platreef,Bushveld Complex,South Africa

In the northern limb of the 2.06-Ga Bushveld Complex, the Platreef is a platinum group elements (PGE)-, Cu-, and Ni-mineralized zone of pyroxenite that developed at the intrusion margin. From north to south, the footwall rocks of the Platreef change from Archaean granite to dolomite, hornfels, and quartzite. Where the footwall is granite, the Sr-isotope system is more strongly perturbed than where the footwall is Sr-poor dolomite, in which samples show an approximate isochron relationship. The Nd-isotope system for samples of pyroxenite and hanging wall norite shows an approximate isochron relationship with an implied age of 2.17 ± 0.2 Ga and initial Nd-isotope ratio of 0.5095. Assuming an age of 2.06 Ga, the ɛNd values range from −6.2 to −9.6 (ave. −7.8, n = 17) and on average are slightly more negative than the Main Zone of the Bushveld. These data are consistent with local contamination of an already contaminated magma of Main Zone composition. The similarity in isotope composition between the Platreef pyroxenites and the hanging wall norites suggests a common origin. Where the country rock is dolomite, the Platreef has generally higher plagioclase and pyroxene δ ¹⁸O values, and this indicates assimilation of the immediate footwall. Throughout the Platreef, there is considerable petrographic evidence for sub-solidus interaction with fluids, and the Δ _{plagioclase–pyroxene} values range from −2 to +6, which indicates interaction at both high and low temperatures. Whole-rock and mineral δD values suggest that the Platreef interacted with both magmatic and meteoric water, and the lack of disturbance to the Sr-isotope system suggests that fluid–rock interaction took place soon after emplacement. Where the footwall is granite, less negative δD values suggest a greater involvement of meteoric water. Consistently higher values of Δ _{plagioclase–pyroxene} in the Platreef pyroxenites and hanging wall norites in contact with dolomite suggest prolonged interaction with CO₂-rich fluid derived from decarbonation of the footwall rocks. The overprint of post crystallization fluid–rock interaction is the probable cause of the previously documented lack of correlation between PGE and sulfide content on the small scale. The Platreef in contact with dolomite is the focus of the highest PGE grades, and this suggests that dolomite contamination played a role in PGE concentration and deposition, but the exact link remains obscure. It is a possibility that the CO₂ produced by decarbonation of assimilated dolomite enhanced the process of PGE scavenging by sulfide precipitation.     

Magnetotelluric investigations for imaging electrical structure of Garhwal Himalayan corridor,Uttarakhand, India

Magnetotelluric investigations have been carried out in the Garhwal Himalayan corridor to delineate the electrical structure of the crust along a profile extending from Indo-Gangetic Plain to Higher Himalayan region in Uttarakhand, India. The profile passing through major Himalayan thrusts: Himalayan Frontal Thrust (HFF), Main Boundary Thrust (MBT) and Main Central Thrust (MCT), is nearly perpendicular to the regional geological strike. Data processing and impedance analysis indicate that out of 44 stations MT data recorded, only 27 stations data show in general, the validity of 2D assumption. The average geoelectric strike, N70°W, was estimated for the profile using tensor decomposition. 2D smooth geoelectrical model has been presented, which provides the electrical image of the shallow and deeper crustal structure. The major features of the model are (i) a low resistivity (<50Ωm), shallow feature interpreted as sediments of Siwalik and Indo-Gangetic Plain, (ii) highly resistive (> 1000Ωm) zone below the sediments at a depth of 6 km, interpreted as the top surface of the Indian plate, (iii) a low resistivity (< 10Ωm) below the depth of 6 km near MCT zone coincides with the intense micro-seismic activity in the region. The zone is interpreted as the partial melting or fluid phase at mid crustal depth. Sensitivity test indicates that the major features of the geoelectrical model are relevant and desired by the MT data.     

Characteristics of spectral aerosol optical depths over India during ICARB

Spectral aerosol optical depth (AOD) measurements, carried out regularly from a network of observatories spread over the Indian mainland and adjoining islands in the Bay of Bengal and Arabian Sea, are used to examine the spatio-temporal and spectral variations during the period of ICARB (March to May 2006). The AODs and the derived Ångström parameters showed considerable variations across India during the above period. While at the southern peninsular stations the AODs decreased towards May after a peak in April, in the north Indian regions they increased continuously from March to May. The Ångström coefficients suggested enhanced coarse mode loading in the north Indian regions, compared to southern India. Nevertheless, as months progressed from March to May, the dominance of coarse mode aerosols increased in the columnar aerosol size spectrum over the entire Indian mainland, maintaining the regional distinctiveness. Compared to the above, the island stations showed considerably low AODs, so too the northeastern station Dibrugarh, indicating the prevalence of cleaner environment. Long-range transport of aerosols from tshe adjoining regions leads to remarkable changes in the magnitude of the AODs and their wavelength dependencies during March to May. HYSPLIT back-trajectory analysis shows that enhanced long-range transport of aerosols, particularly from the west Asia and northwest coastal India, contributed significantly to the enhancement of AOD and in the flattening of the spectra over entire regions; if it is the peninsular regions and the island Minicoy are more impacted in April, the north Indian regions including the Indo Gangetic Plain get affected the most during May, with the AODs soaring as high as 1.0 at 500 nm. Over the islands, the Ångström exponent (α) remained significantly lower (～1) over the Arabian Sea compared to Bay of Bengal (BoB) (～1.4) as revealed by the data respectively from Minicoy and Port Blair. Occurrences of higher values of α, showing dominance of accumulation mode aerosols, over BoB are associated well with the advection, above the boundary layer, of fine particles from the east Asian region during March and April. The change in the airmass to marine in May results in a rapid decrease in α over the BoB.