Numerical modeling for analyzing thermal surface anomalies induced by underground coal fires

Coal seams burning underneath the surface are recognized all over the world and have drawn increasing public attention in the past years. Frequently, such fires are analyzed by detecting anomalies like increased exhaust gas concentrations and soil temperatures at the surface. A proper analysis presumes the understanding of involved processes, which determine the spatial distribution and dynamic behavior of the anomalies.In this paper, we explain the relevance of mechanical and energy transport processes with respect to the occurrence of temperature anomalies at the surface. Two approaches are presented, aiming to obtain insight into the underground coal fire situation: In-situ temperature mapping and numerical simulation. In 2000 to 2005, annual temperature mapping in the Wuda (Inner Mongolia, PR China) coal fire area showed that most thermal anomalies on the surface are closely related to fractures, where hot exhaust gases from the coal fire are released. Those fractures develop due to rock mechanical failure after volume reduction in the seams. The measured signals at the surface are therefore strongly affected by mechanical processes.More insight into causes and effects of involved energy transport processes is obtained by numerical simulation of the dynamic behavior of coal fires. Simulations show the inter-relation between release and transport of thermal energy in and around underground coal fires. Our simulation results show a time delay between the coal fire propagation and the observed appearance of the surface temperature signal. Additionally, the overall energy flux away from the burning coal seam into the surrounding bedrock is about 30-times higher than the flux through the surface. This is of particular importance for an estimation of the energy released based on surface temperature measurements. Finally, the simulation results also prove that a fire propagation rate estimated from the interpretation of surface anomalies can differ from the actual rate in the seam.     

The palaeoenvironment associated with a partial Iguanodon skeleton from the Upper Weald Clay (Barremian, Early Cretaceous) at Smokejacks Brickworks (Ockley, Surrey, UK), based on palynomorphs and ostracods

In 2001 a partial skeleton of an Iguanodon was discovered in the Upper Weald Clay (Barremian, Early Cretaceous) at Smokejacks Brickworks near Ockley, Surrey, UK. When the dinosaur was excavated, a detailed stratigraphic section was logged and 25 samples taken for palynological and micropalaeontological (ostracod and megaspore) analysis, including a detailed sample set of the dinosaur bed itself. Qualitative and quantitative analysis of the palynoflora revealed rich and well-preserved non-marine assemblages of pollen and spores, including early angiosperms, and freshwater green algae. Four types of angiosperm pollen are described and assigned to the genus Retimonocolpites Pierce, 1961, but left in open nomenclature. Some marine elements such as dinoflagellate cysts are identified as the result of reworking of Middle and Upper Jurassic sediments. The pollen/spore assemblages depict a vegetational change from principally gymnosperm-dominated assemblages at the base to principally pteridophyte-dominated assemblages at the top of the section. The dinosaur bed shows a pteridophyte-dominated assemblage, with a significantly high amount of the freshwater green alga Scenedesmus novilunaris He Cheng-quan et al., 1992. Samples close to the dinosaur bed yielded the first useful ostracod finds from Smokejacks Brickworks: well-preserved assemblages containing Cypridea clavata (Anderson, 1939), Damonella cf. pygmaea (Anderson, 1941), Stenestroemia cf. cressida Anderson, 1971 and Stenestroemia sp. A, and fragments and damaged valves of a thin-shelled ostracod, possibly belonging to Mantelliana Anderson, 1966. Those identified as Cypridea clavata show a wide range of morphological variety and in opposition to Anderson's (1967, 1985) taxonomic scheme, which would assign them to up to five different taxa, they are considered to be intraspecific variants of a single species. The possibilities and limitations of age determination of the Wealden sediments using palynomorphs and ostracods are discussed; distinct forms of early angiosperm pollen, together with the ostracod fauna, are consistent with an early Barremian age. Pollen and spores are discussed in terms of their parent plants and the reconstruction of vegetation and palaeoclimate. Palynology and ostracods give evidence for temporary freshwater conditions at the time when the Iguanodon died and the carcase was buried.     

Ultrahigh pressure macro diamonds from Copeton (New South Wales, Australia), based on Raman spectroscopy of inclusions

Mining of Cenozoic alluvial deposits at Copeton and Bingara (Eastern Australia) has produced two million macrodiamonds (0.25 ct median size). Raman spectroscopy is used to identify included minerals within uncut Copeton diamonds, with sealed chamber remnant pressures of 31.7 to 35.6 kbar for coesite, 13.6 and 22.7 kbar for clinopyroxene, and 7.6 kbar for grossular garnet. Assuming elastic behaviour, these values generate inclusion entrapment PT loci which intersect, restricting diamond formation conditions: from 250 °C, 43 kbar to 800 °C, 52 kbar. Larger than error (± 100 °C and ± 4 kbar), this range shows a systematic variation in inclusion composition with diamond zoning and N properties. Published research shows 1) Copeton and Bingara diamonds are unique, and 2) modern alluvium in the Bingara district carries mantle-formed garnet, captured by post-tectonic alkali basalt from an extensive diamondiferous ultrahigh pressure (UHP) terrane that stalled at depth because it is dominated by mafic eclogite. The combined Raman and geological results indicate two sets of subduction UHP diamond formation conditions/protolith are required, firstly cooler oceanic slab and secondly including higher temperature continental crust. The Copeton and Bingara stones are UHP macrodiamonds, and Carboniferous ⁴⁰Ar/³⁹Ar age dates on clinopyroxene inclusions should be interpreted as ages of crystallisation, representing the termination of subduction. The characteristic features of ruptured inclusions and etched percussion marks on Copeton and Bingara diamond indicate volcanic delivery to the earth's surface. Alluvial deposits elsewhere in Eastern Australia may carry similar diamond along with diamond of different origin.     

The effect of extracellular polymers (EPS) on the proton adsorption characteristics of the thermophile Bacillus licheniformis S-86

This study investigated proton adsorption to an extracellular polymeric substance (EPS) producing bacterial strain, Bacillus licheniformis S-86, in order to characterise and quantify the contribution made by EPS to cell surface reactivity. Potentiometric titrations were conducted using both untreated cells and cells from which the EPS layer had been extracted. Surface-complexation modelling indicated the presence of four different functional groups in both untreated and EPS-free cells. These sites are assigned to phosphodiester, (pK_a 3.3–3.4), carboxylic (pK_a 5.3–5.4), phosphoryl/ (pK_a 7.4–7.5) and hydroxyl/amine (pK_a 9.9–10.1) type groups. The pK_a values for the four groups were very similar for untreated and EPS-free cells, indicating no qualitative difference in composition, but site concentrations in the untreated cells were statistically found to be significantly higher than those in the EPS-free cells for the pK_a 3.3–3.4 and pK_a 9.9–10.1 sites. Infrared analysis provided supporting evidence that site 2 is carboxylic in nature but did not reveal any difference in IR absorption between the native and EPS-free cells. Dissolved organic carbon (DOC) analysis conducted during this study indicated that DOC release by cells is significant, and that the EPS layer is the major contributor.     

Wind wave-forced fine sediment erosion during the slack water periods in Hangzhou Bay,China

The shallow water wave simulation model-SWAN incorporated with a simple fine sediment erosion model is applied to Hangzhou Bay, China, to model the horizontal distribution of the maximum bottom orbital velocity and corresponding fine sediment erosion rates induced by: (1) southeasterly steady winds (5, 20 and 30 m/s), (2) southwesterly steady winds (5 and 20 m/s); (3) northwesterly steady winds (5 and 20 m/s); (4) east-southeasterly steady winds (5 and 20 m/s); (5) easterly steady winds (5 and 20 m/s) under closed and unclosed boundaries; and (6) unsteady winds during the slack water periods. Results suggest: (1) the steady wind wave-induced maximum bottom orbital velocities and corresponding fine sediment erosion rates generally increased with the increasing steady winds; (2) closed and unclosed boundary conditions had more significant influences on modeled fine sediment erosion rates under 5 m/s easterly steady winds than 20 m/s; and (3) steady and unsteady wind wave-induced maximum bottom currents could be significant in eroding fine sediment bed in Hangzhou Bay. The results show implications for geomorphology, sedimentology, coastal erosion, and environmental pollution mitigation in Hangzhou Bay.     

A design methodology for rock slopes susceptible to wedge failure using fracture system modelling

This paper demonstrates how the use of fracture system modelling can be linked to limit equilibrium analysis of rock slopes susceptible to wedge failure. The use of fracture systems highlights some of the limitations inherent in traditional structural data analysis and representation. Consequently it allows for more comprehensive input data that can be used for stability analysis of rock slopes. In particular the developed methodology addresses important issues such as spatial variability and wedge size distributions. The paper introduces a series of guidelines for interpretation of the results of rock slopes. The proposed techniques arguably result in an improved level of confidence in the design of rock slopes susceptible to wedge failure.     

Late Palaeozoic hydrocarbon migration through the Clair field, West of Shetland, UK Atlantic margin

Geochemical analysis of bitumen- and hydrocarbon-bearing fluid inclusions from the Devonian-Carboniferous Clair field indicates that the reservoirs contain a mixture of oils from different marine and lacustrine sources. Reconstruction of the Clair field oil-charge history using fluid inclusion petrography show that oil-charging occurred at times of K-feldspar, quartz and calcite cementation. Temperature–composition–time data yielded from the integration of fluid inclusion microthermometry with high-resolution Ar–Ar dating, date hydrocarbon-bearing K-feldspar overgrowths at 247 ± 3.3 Ma. These data show that in order for oil to be trapped within primary fluid inclusions in K-feldspar overgrowths, hydrocarbon migration throughout the UK Atlantic margin must have been taking place during the Late Palaeozoic and as such, current industry oil-play models based solely on oil charging from Jurassic-Cretaceous marine sources are clearly incomplete and need revision. Apatite fission track analysis and vitrinite reflectance data were used to reconstruct thermal burial histories and assess potential oil generation from Middle Devonian lacustrine source rocks. Thermal history data from wells along The Rona Ridge adjacent to the Clair field show that the Palaeozoic section was heated to greater than 100 °C at some time between 270 and 230 Ma, confirming that Devonian source rocks were mature and expelling oil during the Late Palaeozoic at the time that authigenic K-feldspar overgrowths were growing in the Clair field.     

Quantitative textural analysis of packings of elongate crystals

The spatial distribution of grains in a solidifying igneous rock controls the physical properties of the crystal mush, and in turn is controlled by the rate of crystal growth and accumulation. A predominant non-spherical habit for igneous minerals brings into question the use of spherical particles in reference packings used for quantification of spatial distribution. Furthermore, variations of crystal clustering/ordering with length scale require spatial statistics which take into account the distribution of particles beyond nearest neighbours. Using random close packings of spherocylinders, we demonstrate the importance of aspect ratio for the aggregation index (usually known as R) and show that packings of spherical particles have more structure than packings of rods. The spatial distribution functions demonstrate that the plagioclase grains in the colonnade from the Holyoke basalt are clustered on a length scale of 0.5 mm. Understanding the controls on grain spatial distribution in igneous rocks will depend on the application of these techniques to well-understood environments.     

The effects of metamorphism on O and Fe isotope compositions in the Biwabik Iron Formation, northern Minnesota

The Biwabik Iron Formation of Minnesota (1.9 Ga) underwent contact metamorphism by intrusion of the Duluth Complex (1.1 Ga). Apparent quartz–magnetite oxygen isotope temperatures decrease from ∼700°C at the contact to ∼375°C at 2.6 km distance (normal to the contact in 3D). Metamorphic pigeonite at the contact, however, indicates that peak temperatures were greater than 825°C. The apparent O isotope temperatures, therefore, reflect cooling, and not peak metamorphic conditions. Magnetite was reset in δ¹⁸O as a function of grain size, indicating that isotopic exchange was controlled by diffusion of oxygen in magnetite for samples from above the grunerite isograd. Apparent quartz–magnetite O isotope temperatures are similar to calculated closure temperatures for oxygen diffusion in magnetite at a cooling rate of ∼5.6°C/kyr, which suggests that the Biwabik Iron Formation cooled from ∼825 to 400°C in ∼75 kyr at the contact with the Duluth Complex. Isotopic exchange during metamorphism also occurred for Fe, where magnetite–Fe silicate fractionations decrease with increasing metamorphic grade. Correlations between quartz–magnetite O isotope fractionations and magnetite–iron silicate Fe isotope fractionations suggest that both reflect cooling, where the closure temperature for Fe was higher than for O. The net effect of metamorphism on δ¹⁸O–δ⁵⁶Fe variations in magnetite is a strong increase in δ¹⁸O_Mt and a mild decrease in δ⁵⁶Fe with increasing metamorphic grade, relative to the isotopic compositions that are expected at the low temperatures of initial magnetite formation. If metamorphism of Iron Formations occurs in a closed system, bulk O and Fe isotope compositions may be preserved, although re-equilibration among the minerals may occur for both O and Fe isotopes. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.