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.     

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     

40Ar/39Ar and oxygen isotope studies of polymetamorphism from Tinos Island, Cycladic blueschist belt, Greece

Abstract Petrological, oxygen isotope and ⁴⁰Ar/³⁹Ar studies were used to constrain the Tertiary metamorphic evolution of the lower tectonic unit of the Cyclades on Tinos. Polyphase high-pressure metamorphism reached pressures in excess of 15 kbar, based on measurements of the Si content in potassic white mica. Temperatures of 450–500° C at the thermal peak of high-pressure metamorphism were estimated from critical metamorphic assemblages, the validity of which is confirmed by a quartz–magnetite oxygen isotope temperature of 470° C. Some ⁴⁰Ar/³⁹Ar spectra of white mica give plateau ages of 44–40 Ma that are considered to represent dynamic recrystallization under peak or slightly post-peak high-pressure metamorphic conditions. Early stages in the prograde high-pressure evolution may be documented by older apparent ages in the high-temperature steps of some spectra. Eclogite to epidote blueschist facies mineralogies were partially or totally replaced by retrograde greenschist facies assemblages during exhumation. Oxygen isotope thermometry of four quartz–magnetite pairs from greenschist samples gives temperatures of 440–470° C which cannot be distinguished from those deduced for the high-pressure event. The exhumation and overprint is documented by decreasing ages of 32–28 Ma in some greenschists and late-stage blueschist rocks, and ages of 30–20 Ma in the lower temperature steps of the Ar release patterns of blueschist micas. Almost flat parts of Ar–Ar release spectra of some greenschist micas gave ages of 23–21 Ma which are assumed to represent incomplete resetting caused by a renewed prograde phase of greenschist metamorphism. Oxygen isotope compositions of blueschist and greenschist facies minerals show no evidence for the infiltration of a δ¹⁸O-enriched fluid. Rather, the compositions indicate that fluid to rock ratios were very low, the isotopic compositions being primarily controlled by those of the protolith rocks. We assume that the fundamental control catalysing the transformation of blueschists into greenschists and the associated resetting of their isotopic systems was the selective infiltration of metamorphic fluid. A quartz–magnetite sample from a contact metamorphic skarn, taken near the Miocene monzogranite of Tinos, gave an oxygen isotope temperature of 555° C and calculated water composition of 9.1%. The value of δ¹⁸O obtained from this water is consistent with a primary magmatic fluid, but is lower than that of fluids associated with the greenschist overprint, which indicates that the latter event cannot be directly related to the monozogranite intrusion.     

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.     

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.     

A dated cave sediment record of Pleistocene transgressions on Berry Head,Southwest England

Berry Head, a limestone headland in Torbay, southwest England, exhibits a series of subaerial marine erosion platforms and raised beaches spanning an altitude range of 97 m. Solution caves on the headland show preferred horizontal development at elevations that are correlated with the marine erosion platforms, and developed in a marine/freshwater mixing zone whose position was controlled by high sea-level still-stands. Corbridge Cave in Berry Head Quarry lies below the raised beaches in Torbay, and contains evidence of three marine transgressions in the form of fine-grained marine ponding deposits with a marine microfauna. Uranium-series dating of intercalated speleothems indicates that a transgression during oxygen isotope stage 5e reached an elevation of 5.8 m OD, while an earlier transgression (probably during stage 7) reached at least 7.2 m OD. These findings are used to constrain possible interpretations of the aminostratigraphy of raised beaches in Southwest Britain, and a correlation of the Unnamed Stage of Bowen, Sykes, Reeves, Miller, Andrews, Brew and Hare with oxygen isotope stage 5e is proposed.     

Seismotectonic significance of the 29 January 1989 Etne earthquake, SW Norway

The integrated analysis of geological, seismological and field observations with lineament data derived from satellite images allows the identification of a possible seismogenic fault zone for an earthquake which occurred near Etne in southwestern Norway, on 29 February 1989. The hypocentre of the earthquake was located at the mid-crust at a depth of 13.8±0.9 km which is typical of small intraplate earthquakes. The Etne earthquake occurred as a result of normal faulting with a dextral strike-slip component on a NW–SE trending fault. Available geological and lineament data indicate correlation of the inferred seismogenic fault with the NW–SE trending Etne fault zone. An aeromagnetic anomaly related to the Etne fault zone forms a regional feature intersecting both Precambrian basement and allochthonous Caledonian rocks. Based on these associations the occurrence of the Etne event is ascribed to the reactivation of a zone of weakness along the Etne fault zone. Slope-instabilities developed in the superficial deposits during the Etne event demonstrate the existence of potentially hazardous secondary-effects of such earthquakes even in low seismicity areas such as southwestern Norway.     

Optimization of grid-drilling using computer simulation

A computer simulation method has been developed to find efficient drilling grids for mineral deposits. A well-known ore deposit is used as a model to develop an efficient pattern for undiscovered ore bodies in the same area or in other prospects where similar geometry is suspected. The model for this study is the Austinville, Virginia deposit, a Mississippi Valley-type deposit composed of 17 ore bodies totaling 34 million short tons (30 million metric tons). The method employs a computer program that simulates drilling the model deposit with different patterns, including various levels of follow-up drilling. Follow-up holes are drilled in fences at one half the original spacing around holes in the grid that show ore-grade mineralization. Each pattern is drilled 100 times from random starting locations to provide a range of outcomes of drilling, including the best, worst, and most likely. For this study, patterns of 100 drill holes were composed of 10 fences spaced 1000–5000 feet (305–1524 m) apart, each with 10 holes spaced 200–1000 feet (61–305 m) apart. In all, 25 grids were used with zero to three levels of follow-up drilling. The 600/2000 grid, with drill holes spaced 600 feet (183 m) apart in fences spaced 2000 feet (610 m) apart, was compared with the 200/5000 grid because they represented contrasting outcomes. The 600/2000 grid penetrated many ore bodies consistently but with few multiple hits to individual ore bodies; whereas the 200/5000 grid inconsistently penetrated few ore bodies with many multiple hits. The 600/2000 grid was more efficient than the 200/5000 grid at hitting large ore bodies of 1,000,000 short tons or greater (900,000 metric tons or greater) and was made more effective by adding one cycle of follow-up drilling. The 600/2000 grid had a 97% chance of hitting one or more large ore bodies with at least one drill hole per ore body, and the 200/5000 grid had a 64% chance. Once hit, there was an 82% chance that the largest ore body would be penetrated by three or more holes when using the 600/2000 grid and an 88% chance using the 200/5000 grid.