Stochastic modelling of groundwater flow and solute transport in aquifers

This paper presents an introductory overview of recently developed stochastic theories for tackling spatial variability problems in predicting groundwater flow and solute transport. Advantages and limitations of the theories are discussed. Lastly, strategies based on the stochastic approaches to predict solute transport in aquifers are recommended.     

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.     

Metamorphic phase relations in orthopyroxene-bearing granitoids: implication for high-pressure metamorphism and prograde melting in the continental crust

In this work, the factors controlling the formation and preservation of high-pressure mineral assemblages in the metamorphosed orthopyroxene-bearing metagranitoids of the Sandmata Complex, Aravalli-Delhi Mobile Belt (ADMB), northwestern India have been modelled. The rocks range in composition from farsundite through quartz mangerite to opdalite, and with varying K₂O, Ca/(Ca + Na)_rock and FeO_tot + MgO contents. A two stage metamorphic evolution has been recorded in these rocks.
An early hydration event stabilized biotite with or without epidote at the expense of magmatic orthopyroxene and plagioclase. Subsequent high-pressure granulite facies metamorphism (∼15 kbar, ∼800 °C) of these hydrated rocks produced two rock types with contrasting mineralogy and textures. In the non-migmatitic metagranitoids, spectacular garnet ± K-feldspar ± quartz corona was formed around reacting biotite, plagioclase, quartz and/or pyroxene. In contrast, biotite ± epidote melting produced migmatites, containing porphyroblastic garnet incongruent solids and leucosomes.
Applying NCKFMASHTO T–M (H₂O) and P–T pseudosection modelling techniques, it is demonstrated that the differential response of these magmatic rocks to high-pressure metamorphism is primarily controlled by the scale of initial hydration. Rocks, which were pervasively hydrated, produced garnetiferous migmatites, while for limited hydration, the same metamorphism formed sub-solidus garnet-bearing coronae. Based on the sequence of mineral assemblage evolution and the mineral compositional zoning features in the two metagranitoids, a clockwise metamorphic P–T path is constrained for the high-pressure metamorphic event. The finding has major implications in formulating geodynamic model of crustal amalgamation in the ADMB.     

Biological and granulometric controls on sedimentary organic matter of an intertidal mudflat

A variety of measures of organic matter concentration and quality were made on samples collected from the top few mm of intertidal mudflat sediment over the course of a year, in order to assess the relative importance of biological and sedimentological influences on sedimentary organic matter. Winter and summer were times of relatively fine-grained sediment accumulation, caused by biological deposition or stabilization processes and resulting in higher organic matter concentrations. Stable carbon isotope and Br:C ratios indicated a planktonic source of bulk organic matter. Ratios of organic carbon to specific surface area of the sediments were consistent with an organic monolayer coverage of sediment grains. Correction for changing grain size during the year showed no change in the organic concentration per unit surface area, in spite of organic matter inputs by in situ primary production, buildup of heterotroph biomass and mucus coatings, and biodeposition of organic-rich seston. There were also no indications of changes in bulk organic quality, measured as hydrolyzable carbohydrates and amino acids, in response to these biological processes. It is concluded that biological processes on a seasonal time scale affect the bulk organic matter of these sediments via a modulation of grain size rather than creation or decay of organic matter.     

Elastic solutions for stresses in a transversely isotropic half‐space subjected to three‐dimensional buried parabolic rectangular loads

In many areas of engineering practice, applied loads are not uniformly distributed but often concentrated towards the centre of a foundation. Thus, loads are more realistically depicted as distributed as linearly varying or as parabola of revolution. Solutions for stresses in a transversely isotropic half‐space caused by concave and convex parabolic loads that act on a rectangle have not been derived. This work proposes analytical solutions for stresses in a transversely isotropic half‐space, induced by three‐dimensional, buried, linearly varying/uniform/parabolic rectangular loads. Load types include an upwardly and a downwardly linearly varying load, a uniform load, a concave and a convex parabolic load, all distributed over a rectangular area. These solutions are obtained by integrating the point load solutions in a Cartesian co‐ordinate system for a transversely isotropic half‐space. The buried depth, the dimensions of the loaded area, the type and degree of material anisotropy and the loading type for transversely isotropic half‐spaces influence the proposed solutions. An illustrative example is presented to elucidate the effect of the dimensions of the loaded area, the type and degree of rock anisotropy, and the type of loading on the vertical stress in the isotropic/transversely isotropic rocks subjected to a linearly varying/uniform/parabolic rectangular load. Copyright © 2002 John Wiley & Sons, Ltd.     

The influence of a non‐associated flow rule on the calculation of the factor of safety of soil slopes

This paper presents a method that incorporates a non‐associated flow rule into the limit analysis to investigate the influence of the dilatancy angle on the factor of safety for the slope stability analysis. The proposed method retain's the advantage of the upper bound method, which is simple and has no stress involvement in the calculation of the energy dissipation and the factor of safety. Copyright © 2001 John Wiley & Sons, Ltd.     

Progress in Earthquake Science and Technology in China: Review and Prospects (Ⅰ)

In the article the author looks back the hard development course and great progress in earth quake science and technology in China during the last near a half of century and expounds the following 3 aspects: (1) The strong desire of the whole society to mitigate seismic disasters and reduce the effect of earthquakes on social-economic live is a great driving force to push forward the development of earthquake science and technology in China; (2) To better ensure people‘ s life and property, sustainable economic development, and social stability is an essential purpose to drive the development of earthquake science and technology in China; and (3) To insist on the dialectical connection of setup of technical system for seismic monitoring with the scientific research of earthquakes and to better handle the relation between crucial task, current scientif ic level, and the feasibility are the important principles to advance the earthquake science and technology in China. Some success and many setbacks in earthquake disaster mitigation consistently enrich our knowledge regarding the complexity of the conditions for earthquake occurrence and the process of earthquake preparation, promote the reconstruction and modernization of technical system for earthquake monitoring, and deepen the scientific research of earthquakes. During the last 5 years, the improvement and modernization of technical system for earthquake monitoring have clearly provided the technical support to study and practice of earthquake prediction and pre caution, give prominence to key problems and broaden the field of scientific research of earth quakes. These have enabled us to get some new recognition of the conditions for earthquake oc currence and process of earthquake preparation, characteristics of seismic disaster, and mecha nism for earthquake generation in China‘s continent. The progress we have made not only en courages us to enhance the effectiveness of earthquake disaster mitigation, but also provides a basis for accelerating further development of earthquake science and technology in China in the new century, especially in the 10th five-year plan. Based on the history reviewed, the author sets forth a general requirement for develop ment of earthquake science and technology in China and brings out 10 aspects to be stressed and strengthened at present and in the future. These are: upgrade and setup of the network of digitized seismic observation; upgrade and setup of the network for observation of seismic pre cursors; setup of the network for observation of strong motion; setup of the laboratories for ex periment on seismic regime; establishment of technical system for seismic information, emer gency command and urgent rescue; research on short-term and imminent earthquake predic tion; research on intermediate- and long-term earthquake prediction; research on attenuation of seismic ground motion, mechanism for seismic disaster, and control on seismic disaster; ba sic research fields related to seismology and geoscience. We expect that these efforts will signifi cantly elevate the level of earthquake science and technology in China to the advanced interna tional level, improve theories, techniques, and methods for earthquake precaution and predic tion, and enhance the effectiveness of earthquake disaster mitigation.     

Late Quaternary submarine bedforms and ice‐sheet flow in Gerlache Strait and on the adjacent continental shelf,Antarctic Peninsula

Geophysical data from Gerlache Strait, Croker Passage, Bismarck Strait and the adjacent continental shelf reveal streamlined subglacial bedforms that were produced at the bed of the Antarctic Peninsula Ice Sheet (APIS) during the last glaciation. The spatial arrangement and orientation of these bedforms record the former drainage pattern and flow dynamics of an APIS outlet up‐flow, and feeding into, a palaeo‐ice stream in the Western Bransfield Basin. Evidence suggests that together, they represent a single ice‐flow system that drained the APIS during the last glaciation. The ice‐sheet outlet flowed north/northeastwards through Gerlache Strait and Croker Passage and converged with a second, more easterly ice‐flow tributary on the middle shelf to form the main palaeo‐ice stream. The dominance of drumlins with low elongation ratios suggests that ice‐sheet outlet draining through Gerlache Strait was comparatively slower than the main palaeo‐ice stream in the Western Bransfield Basin, although the low elongation ratios may also partly reflect the lack of sediment. Progressive elongation of drumlins further down‐flow indicates that the ice sheet accelerated through Croker Passage and the western tributary trough, and fed into the main zone of streaming flow in the Western Bransfield Basin. Topography would have exerted a strong control on the development of the palaeo‐ice stream system but subglacial geology may also have been significant given the transition from crystalline bedrock to sedimentary strata on the inner–mid‐shelf. In the broader context, the APIS was drained by a number of major fast‐flowing outlets through cross‐shelf troughs to the outer continental shelf during the last glaciation. Copyright © 2004 John Wiley & Sons, Ltd.     

Low-cost volcano surveillance from space: case studies from Etna, Krafla, Cerro Negro, Fogo, Lascar and Erebus

 Satellite data offer a means of supplementing ground-based monitoring during volcanic eruptions, especially at times or locations where ground-based monitoring is difficult. Being directly and freely available several times a day, data from the advanced very high resolution radiometer (AVHRR) offers great potential for near real-time monitoring of all volcanoes across large (3000×3000 km) areas. Herein we describe techniques to detect and locate activity; estimate lava area, thermal flux, effusion rates and cumulative volume; and distinguish types of activity. Application is demonstrated using data for active lavas at Krafla, Etna, Fogo, Cerro Negro and Erebus; a pyroclastic flow at Lascar; and open vent systems at Etna and Stromboli. Automated near real-time analysis of AVHRR data could be achieved at existing, or cheap to install, receiving stations, offering a supplement to conventional monitoring methods. Received: 21 January 1997 / Accepted: 3 April 1997