Simulating methane emission from a Chinese rice field as influenced by fertilizer and water level

Natural and agricultural wetlands are considered to be the major sources of global atmospheric methane (CH₄). A one‐dimensional model was developed to simulate methane emission and used to examine the influence of various physical processes on the rate of methane emission. Three processes involved in the methane emission are implemented in the model: production, reoxidation and transport. Three transport pathways were considered: diffusion across water–air or soil–air interfaces, ebullition and diffusion through plants. These pathways are influenced by soil properties, plant growth, water‐table conditions, temperature and external inputs (e.g. fertilizer). The model was used to examine the seasonal variation of the methane emission at a rice field in Hunan, China, which was observed during a field experiment for consecutive (early and late) rice seasons in 1992. The observed seasonal variations of methane emission, and role of plants in transporting methane to the atmosphere, are captured by the model simulation. Further model applications were conducted to simulate effects of fertilizer and water‐level condition on the methane emission. The results indicate that unfermented organic fertilizer produces a higher methane emission rate than mineral fertilizer. The simulations with treatments of a deep‐water covering and constant moisture reduced the methane emission. The rice field study provides a framework for further development of the model towards simulations based on spatially distributed variables (e.g. water table, soil temperature and vegetation) at a regional scale. Copyright © 2003 John Wiley & Sons, Ltd.     

Tensile stresses around boreholes due to transient fluid flow

We present solutions for the effective stress induced by gas flow through a porous solid into a borehole resulting from sudden pressure reduction. Tensile effective stress that exceeds the strength of the solid will lead to borehole failure. This has applications to the intentional creation of cavities, relevant to the efficient recovery of coalbed methane, and the avoidance of borehole stability problems in conventional gas production.     

Estimation of the sub-iron-iron flux ratio in cosmic rays over alice springs using a balloon flight experiment with a CR-39 detector

The sub-iron-iron flux ratio in cosmic rays at an atmospheric depth of 9.8 g cm^–2 has been estimated using a balloon-borne CR-39 (HCB) stack launched from Alice Springs for 32 hours. The recovered and chemically etched plates were analysed optically and the measured etch pit diameters yielded the sub-iron-iron flux ratio at the flight altitude. The sub-iron-iron flux ratio has been corrected for the top of the atmosphere by considering the loss of heavy ions due to nuclear interaction and fragmentation. The present result has been compared with the result expected from the source composition derived by Protheroeet al. as well as other authors.     

The estimation of the shape of an array using a hidden Markov model

A hidden Markov model (HMM) technique for the estimation of the shape of a towed array is presented. It is assumed that there is a far-field source radiating sound containing possibly weak spectral lines. The technique uses either the Fourier coefficients at a given frequency computed from a single time block or the maximal eigenvector of a sample spectral covariance matrix. The technique is illustrated using several simulations. The results of these simulations indicate that the HMM technique yields shape and bearing estimates more accurate than those provided by a maximum-likelihood array shape estimation technique     

Late Oligocene-Miocene syn-/-late-orogenic successions in Western and Central Mediterranean Chains from the Betic Cordillera to the Southern Apennines*

The identification of syn- and late-orogenic flysch deposits, extending from the Betic Cordillera to the Southern Apennines, assists in the reconstruction of the tectonic-sedimentary evolution of the perimediterranean chains. A microplate was located between the European and African Plates during the Late Jurassic–Early Cretaceous, bordered northwards by the Piemontese Ocean and southwards by another (North Africa ‘Flysch’ Basin or Maghrebian) Ocean. The Piemontese Ocean and the northern margin of the microplate were structured from the Late Cretaceous to the Eocene to create an Eo-alpine Chain. The southern margin of the microplate was deformed in the Aquitanian, when the internal areas of the Maghrebian Ocean were characterized by syn-orogenic flysch deposits. This episode culminated with metamorphism (25–22 Ma) and nappe emplacement, which destroyed the former palaeogeography and created an orogenic belt (AlKaPeCa). Afterwards, a lower Burdigalian late-orogenic cycle started in the deformed area, which as a result of the opening of the Algero-Provençal Basin, caused the fragmentation of the AlKaPeCa, its thrusting on the ‘Flysch’ Basin and the collision with the North Africa and South Iberia Margins. These latter were folded and thrusted, the ‘Flysch’ Units pushed over the External Domain and also back-thrusted. Langhian late-orogenic deposits suture the new tectonic features. Finally, the whole orogen was thrust onto the foredeep during the Middle–Late Miocene.     

Transient response of a circular cavity in a poroelastic medium

This paper considers the transient response of a pressurized long cylindrical cavity in an infinite poroelastic medium. To obtain transient solutions, Biot's equations for poroelastodynamics are specialized for this problem. A set of exact general solutions for radial displacement, stresses, pore pressure and discharge are derived in the Laplace transform space by using analytical techniques. Solutions are presented for three different types of prescribed transient radial pressures acting on the surface of a permeable as well as an impermeable cavity surface. Time domain solutions are obtained by inverting Laplace domain solutions using a reliable numerical scheme. A detailed parametric study is presented to illustrate the influence of poroelastic material parameters and hydraulic boundary conditions on the response of the medium. Comparisons are also presented with the corresponding ideal elastic solutions to portray the poroelastic effects. It is noted that the maximum radial displacement and hoop stress at the cavity surface are substantially higher than the classical static solutions and differ considerably from the transient elastic solutions. Time histories and radial variations of displacement, hoop stress, pore pressure and fluid discharge corresponding to a cavity in two representative poroelastic materials are also presented.     

Shakedown in layered pavements under moving surface loads

Elastic–plastic deformations in pavements consisting of layers of different frictional materials are investigated. The upper bound, kinematic shakedown theorem is used to obtain estimates of the critical shakedown loads. Fully general three-dimensional deformations are considered. The influence of the loading distribution, interactions between loads, and the effect of varying the thickness, stiffness and strength of the layers are explored. Consequences of this investigation for particular existing designs of flexible pavements are investigated. In particular, it is found that the strength of the subgrade has no effect on the magnitude of the critical shakedown load.     

Upscaling of elastic properties for large scale geomechanical simulations

Large scale geomechanical simulations are being increasingly used to model the compaction of stress dependent reservoirs, predict the long term integrity of under‐ground radioactive waste disposals, and analyse the viability of hot‐dry rock geothermal sites. These large scale simulations require the definition of homogenous mechanical properties for each geomechanical cell whereas the rock properties are expected to vary at a smaller scale. Therefore, this paper proposes a new methodology that makes possible to define the equivalent mechanical properties of the geomechanical cells using the fine scale information given in the geological model. This methodology is implemented on a synthetic reservoir case and two upscaling procedures providing the effective elastic properties of the Hooke's law are tested. The first upscaling procedure is an analytical method for perfectly stratified rock mass, whereas the second procedure computes lower and upper bounds of the equivalent properties with no assumption on the small scale heterogeneity distribution. Both procedures are applied to one geomechanical cell extracted from the reservoir structure. The results show that the analytical and numerical upscaling procedures provide accurate estimations of the effective parameters. Furthermore, a large scale simulation using the homogenized properties of each geomechanical cell calculated with the analytical method demonstrates that the overall behaviour of the reservoir structure is well reproduced for two different loading cases. Copyright © 2004 John Wiley & Sons, Ltd.