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.     

Fracture energy release and size effect in borehole breakout

The paper presents a simple approximate analytical solution of the remote stresses that cause the collapse of a borehole or other circular cylindrical cavity in an infinite elastic space. Regions of parallel equidistant splitting cracks are assumed to form on the sides of the cavity. Their boundary is assumed to be an ellipse of a growing horizontal axis, the other axis remaining equal to the borehole diameter. The slabs of rock between the splitting cracks are assumed to buckle as slender columns, and their post-critical stress is considered as the residual stress in the cracked rock. The buckling of these slab columns is assumed to be resisted not only by their elastic bending stiffness but also shear stresses produced on rough crack faces by relative shear displacements. The energy release from the infinite medium caused by the growth of the elliptical cracking region is evaluated according to Eschelby's theorem. This release is set equal to the energy dissipated by the formation of all the splitting cracks, which is calculated under the assumption of constant fracture energy. This yields the collapse stress as a function of the elastic moduli, fracture energy, ratio of the remote principal stresses, crack shear resistance characteristic and borehole diameter. The collapse stress as a function of crack spacing is found to have a minimum, and the correct crack spacing is determined from this minimum. For small enough diameters, the crack spacing increases as the (4/5)-power of the borehole diameter, while for large enough diameters a constant spacing is approached. In contrast to plastic solutions, the breakout stress exhibits a size effect, such that for small enough diameters the breakout stress decreases as the (? 2/5)-power of the borehole diameter, while for large enough diameters a constant limiting value is approached. Finally, some numerical estimates are given and the validity of various simplifying assumptions made is discussed.     

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     

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.     

Geometry and kinematics of recent deformation in the Mondy–Tunka area (south-westernmost Baikal rift zone, Mongolia–Siberia)

We used satellite imagery and field data to investigate the south‐westernmost Baikal rift zone. We focus our study in the Mondy and Ikhe Ukhgun valleys, site of an Mw = 6.9 seismic event in 1950. Surface deformations are observed along the E–W‐trending Mondy strike‐slip fault and along the Ikhe Ukhgun thrust. The Mondy fault system is 80 km long and is composed of four segments 10–15 km long. These segments are characterized by subvertical planes with left‐lateral movements. The Ikhe Ukhgun thrust is 20 km long, dips 40° to the south and shows reverse movement with a left‐lateral component. These observations are consistent with the present‐day regional NNE–SSW compression and with the focal mechanism of the 1950 Mondy earthquake that was recently re‐evaluated. These features, like those observed in the Tunka basin, demonstrate a recent change of regional strain regime from transtension to transpression that we place before the Late Pleistocene.     

Preliminary organic petrography of Palaeozoic rocks from the Grand Banks,Newfoundland

Samples of lower Palaeozoic bedrock from the Grand Banks of Newfoundland were examined using reflected light microscopy and Rock-Eval pyrolysis. These samples contained organic material which included bitumen (structureless organic material), acritarchs, chitinozoa, scolecodonts and graptolites. The reflectance of the organic material increased from acritarchs towards graptolites with bitumen showing a wide range of reflectance. Multiple phases of oil migration through one sample were inferred by examination of the bitumens: the lowest reflecting being more recent than the highest reflecting bitumen. Combined reflectance and Rock-Eval analysis indicate that these samples are mature to overmature. The two oldest samples (Arenig-Llanvirn, Ordovician) are assessed as mature and as having potential for generation of liquid hydrocarbons, but must be considered as mainly gas prone.     

Hydrodynamic control of phytoplankton in low salinity waters of the James River estuary,Virginia, U.S.A.

Autotrophic biomass and productivity as well as nutrient distributions and phytoplankton cell populations in the James River estuary, Virginia, were quantified both spatially and temporally over a 17-month period. Emphasis was placed on the very low salinity region of the estuary in order to gain information on the fate of freshwater phytoplankters. Differing amounts of freshwater plant biomass are advected into the estuary as living material, DOC or POC and the demonstrated variability of this input must play an important role in marine biogeochemical cycling.Late summer and fall maxima in both chlorophyll a and the photosynthetic production of particulate organic carbon in very low salinity regions were inversely correlated with river discharge.During periods of low river discharge greater than 50% of the chlorophyll a biomass measured at 0‰ disappeared within a narrow range of salinity (0–2‰). Cell enumeration data suggest that species introduced from the freshwater end-member tend to comprise the bulk of the biomass removed. Confounding factors, which may contribute to the regulation of both the abundance and species of phytoplankters mid-river, include the flocculation of colloidal material with phytoplankton cells, the presence of the turbidity maximum and the growth of endemic phytoplankton populations.An inverse relationship exists between the phytoplankton abundance in very low salinity waters and the abundance of biomass measured in the lower portion of the river (estuary). Thus, autotrophic production in the fresh and very low salinity areas may indirectly regulate the onset on the spring bloom in the estuary by controlling the amount of nutrients available.