Transient‐state groundwater flow in various geometries with wells at arbitrary positions: analytical methods

We have developed a method for analytically solving the porous medium flow equation in many different geometries for horizontal (two‐dimensional), homogeneous and isotropic aquifers containing impermeable boundaries and any number of pumping or injection wells located at arbitrary positions within the system. Solutions and results are presented for rectangular and circular aquifers but the method presented here is easily extendible to many geometries. Results are also presented for systems where constant head boundary conditions can be emulated internal to the aquifer boundary. Recommendations for extensions of the present work are briefly discussed. Copyright © 2003 John Wiley & Sons, Ltd.     

A major sulphur isotope event at c. 510 Ma: a possible anoxia–extinction–volcanism connection during the Early–Middle Cambrian transition?

A new approach to constraining seawater δ³⁴S and sulphate concentration using francolite‐bound sulphate reveals an abrupt increase in δ³⁴S to +50‰ around the Early–Middle Cambrian boundary. Such high δ³⁴S values are best explained by increased rates of pyrite burial due to ocean anoxia coupled with an increased sensitivity of the ocean sulphate reservoir to perturbations due to low sulphate concentrations of 500–700 μgL^?1. We argue that the spread of anoxic waters at this time was partly the result of greenhouse warming related to the eruption of the Kalkarindji Large Igneous Province of northern Australia and that it triggered the collapse of early metazoan reef ecosystems during the latest Early Cambrian. Mass extinctions of the last 260 Myr have all coincided with enhanced volcanic activity, while several are also associated with positive shifts in seawater δ³⁴S. Extending this correlation back in time further implicates volcanically induced climate change as a major determining factor in biosphere evolution. Terra Nova, 18, 257–263, 2006     

Particulate organic matter in the coastal and estuarine waters of Goa and its relationship with phytoplankton production

In the coastal and estuarine waters of Goa, particulate organic carbon (POC) varied from 0.52 to 2.51 mg l^?1 and from 0.28 to 5.24 mg l^?1 and particulate phosphorus (PP) varied from 0.71 to 5.18 μg l^?1 and from 0.78 to 20.34 μg l^?1, respectively. The mean values of chlorophyll and primary productivity were 1.94 mg m^?3 and 938.1 mg C m^?2 day^?1 in the coastal waters and 4.3 mg m^?3 and 636.5 mg C m^?1 day^?1 in the estuarine waters, respectively.$\text{POC}\text{chl}$ ratios were low in June and October even when POC values were quite high. The POC in surface waters was linearly correlated with the chlorophyll content. Also PP increased when chlorophyll and primary productivity remained high. The results suggest that the phytoplankton was sharply increasing and contributed to POC and PP content. The percentage of detritus calculated from the intercept values of chlorophyll on POC varied from 46 to 76% depending on season. Results indicate that the major portion of POC and PP during postmonsoon (October–January) is derived from phytoplankton production while the allochthonous matter predominate during monsoon (June–September).     

High-temperature geothermal reservoir possibilities in Hungary

Theoretical studies have shown the possibility of high-temperature ('high enthalpy') geothermal reservoirs in the pre-Tertiary rocks at 4–5 km depth range within the Pannonian Basin. This expectation was proven by the hotwater/steam blowout of Fábiánsebestyén-4 borehole (16.12.85–31.1.86). Exploration efforts carried out during 1987–88 in the broad vicinity of the borehole proved that reservoirs of this type can be found with the combination of seismic reflection, silica-thermometry and magnetotelluric sounding methods. Deliberate prospection should be continued in all suitable areas within the basin, since high enthalpy reservoirs promise profitable operation of geothermal power stations.     

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.     

A micromechanical finite element model for linear and damage‐coupled viscoelastic behaviour of asphalt mixture

This study presents a finite element (FE) micromechanical modelling approach for the simulation of linear and damage‐coupled viscoelastic behaviour of asphalt mixture. Asphalt mixture is a composite material of graded aggregates bound with mastic (asphalt and fine aggregates). The microstructural model of asphalt mixture incorporates an equivalent lattice network structure whereby intergranular load transfer is simulated through an effective asphalt mastic zone. The finite element model integrates the ABAQUS user material subroutine with continuum elements for the effective asphalt mastic and rigid body elements for each aggregate. A unified approach is proposed using Schapery non‐linear viscoelastic model for the rate‐independent and rate‐dependent damage behaviour. A finite element incremental algorithm with a recursive relationship for three‐dimensional (3D) linear and damage‐coupled viscoelastic behaviour is developed. This algorithm is used in a 3D user‐defined material model for the asphalt mastic to predict global linear and damage‐coupled viscoelastic behaviour of asphalt mixture. For linear viscoelastic study, the creep stiffnesses of mastic and asphalt mixture at different temperatures are measured in laboratory. A regression‐fitting method is employed to calibrate generalized Maxwell models with Prony series and generate master stiffness curves for mastic and asphalt mixture. A computational model is developed with image analysis of sectioned surface of a test specimen. The viscoelastic prediction of mixture creep stiffness with the calibrated mastic material parameters is compared with mixture master stiffness curve over a reduced time period. In regard to damage‐coupled viscoelastic behaviour, cyclic loading responses of linear and rate‐independent damage‐coupled viscoelastic materials are compared. Effects of particular microstructure parameters on the rate‐independent damage‐coupled viscoelastic behaviour are also investigated with finite element simulations of asphalt numerical samples. Further study describes loading rate effects on the asphalt viscoelastic properties and rate‐dependent damage behaviour. Copyright © 2006 John Wiley & Sons, Ltd.     

Radon variations in a Hungarian village

 A steady radon exhalation is assumed in most publications. In a village of North-East Hungary, however, high radon concentrations have been measured, differing strongly in neighbouring houses and varying in time, due to the interplay of geochemical phenomena. Received: 20 November 1995 · Accepted: 18 June 1996     

Recognition of Significant Surface Soil Geochemical Anomalies Via Weighted 3D Shortest-Distance Field of Subsurface Orebodies: A Case Study in the Hongtoushan Copper Mine,NE China

Natural Resources Research - Quantitative prediction of concealed mineralization is always confronted with difficulties in comprehensive analysis between 2D and 3D data and between qualitative and...     

Predicting saturation of gas hydrates using pre-stack seismic data, Gulf of Mexico

A promising method for gas hydrates exploration incorporates pre-stack seismic inversion data, elastic properties modeling, and seismic interpretation to predict saturation of gas hydrates (Sgh). The technology can be modified slightly and used for predicting hydrate concentrations in shallow arctic locations as well. Examples from Gulf of Mexico Walker Ridge (WR) and Green Canyon (GC) protraction areas illustrate how Sgh was derived and used to support the selection of well locations to be drilled for gas hydrates in sand reservoirs by the Chevron-led Joint Industry Project (JIP) Leg II cruise in 2009. Concentrations of hydrates were estimated through the integration of seismic inversion of carefully conditioned pre-stack data, seismic stratigraphic interpretation, and shallow rock property modeling. Rock property trends were established by applying principles of rock physics and shallow sediment compaction, constrained by regional geological knowledge. No nearby sonic or density logs were available to define the elastic property trends in the zone of interest. Sgh volumes were generated by inverting pre-stack data to acoustic and shear impedance (PI and SI) volumes, and then analyzing deviations from modeled impedance trends. In order to enhance the quality of the inversion, we stress the importance of maximizing the signal to noise ratio of the offset data by conditioning seismic angle gathers prior to inversion. Seismic interpretation further plays an important role by identifying false anomalies such as hard, compact strata, which can produce apparent high Sgh values, and by identifying the more promising strata and structures for containing the hydrates. This integrated workflow presents a highly promising methodology, appropriate for the exploration of gas hydrates.