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.     

Organic carbon deposition and phosphorus accumulation during Oceanic Anoxic Event 2 in Tarfaya, Morocco

With a multi-proxy approach, an attempt was made to constrain productivity and bottom-water redox conditions and their effects on the phosphorus accumulation rate at the Mohammed Plage section on the Tarfaya coast, Morocco, during the Cenomanian-Turonian Anoxic Event (OAE 2). A distinct δ¹³C_org isotope excursion of +2.5‰ occurs close to the top of the section. The unusually abrupt shift of the isotope excursion and disappearance of several planktonic foraminiferal species (e.g. Rotalipora cushmani and Rotalipora greenhornensis) in this level suggests a hiatus of between 40–60 kyrs at the excursion onset. Nevertheless, it was possible to determine both the long-term environmental history as well as the processes that took place immediately prior to and during OAE 2. TOC% values increase gradually from the base of the section to the top (from 2.5% to 10%). This is interpreted as the consequence of a long-term eustatic sea-level rise and subsidence causing the encroachment of less oxic waters into the Tarfaya Basin. Similarly a reduction in the mineralogically constructed ‘detrital index’ can be explained by the decrease in the continental flux of terrigenous material due to a relative sea-level rise. A speciation of phosphorus in the upper part of the section, which spans the start and mid-stages of OAE 2, shows overall higher abundances of P_reactive mass accumulation rates before the isotope excursion onset and lower values during the plateau. Due to the probable short hiatus, the onset of the decrease in phosphorus content relative to the isotope excursion is uncertain, although the excursion plateau already contains lower concentrations. The C_org/P_total and V/Al ratios suggest that this reduction was mostly likely caused by a decrease in the available bottom oxygen content (probably as a result of higher productivity) and a corresponding fall in the phosphorus retention ability of the sediment. Productivity appears to have remained high during the isotope plateau possibly due to a combination of ocean-surface fertilisation via increased aridity (increased K/Al and Ti/Al ratios) and/or higher dissolved inorganic phosphorus content in the water column as a result of the decrease in sediment P retention. The evidence for decreased P-burial has been observed in many other palaeoenvironments during OAE 2. Tarfaya's unique upwelling paleosituation provides strong evidence that the nutrient recycling was a global phenomenon and therefore a critical factor in starting and sustaining OAE 2.     

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.     

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.     

Oxidation of pyrite during extraction of carbonate associated sulfate

The sulfur isotopic composition of carbonate associated sulfate (CAS) has been used to investigate the geochemistry of ancient seawater sulfate. However, few studies have quantified the reliability of δ³⁴S of CAS as a seawater sulfate proxy, especially with respect to later diagenetic overprinting. Pyrite, which typically has depleted δ³⁴S values due to authigenic fractionation associated with bacterial sulfate reduction, is a common constituent of marine sedimentary rocks. The oxidation of pyrite, whether during diagenesis or sample preparation, could thus adversely influence the sulfur isotopic composition of CAS. Here, we report the results of CAS extractions using HCl and acetic acid with samples spiked with varying amounts of pyrite. The results show a very strong linear relationship between the abundance of fine-grained pyrite added to the sample and the resultant abundance and δ³⁴S value of CAS. This data represents the first unequivocal evidence that pyrite is oxidized during the CAS extraction process. Our mixing models indicate that in samples with much less than 1 wt.% pyrite and relatively high δ³⁴S_pyrite values, the isotopic offset imparted by oxidation of pyrite should be much less than ? 4‰. A wealth of literature exists on the oxidation of pyrite by Fe³⁺ and we believe this mechanism drives the oxidation of pyrite during CAS extraction, during which the oxygen used to form sulfate is taken from H₂O, not O₂. Consequently, extracting CAS under anaerobic conditions would only slow, but not halt, the oxidation of pyrite. Future studies of CAS should attempt to quantify pyrite abundance and isotopic composition.     

Application of precise 142Nd/144Nd analysis of small samples to inclusions in diamonds (Finsch,South Africa) and Hadean Zircons (Jack Hills,Western Australia)

¹⁴⁶Sm–¹⁴²Nd and ¹⁴⁷Sm–¹⁴³Nd systematics were investigated in garnet inclusions in diamonds from Finsch (S. Africa) and Hadean zircons from Jack Hills (W. Australia) to assess the potential of these systems as recorders of early Earth evolution. The study of Finsch inclusions was conducted on a composite sample of 50 peridotitic pyropes with a Nd model age of 3.3 Ga. Analysis of the Jack Hills zircons was performed on 790 grains with ion microprobe ²⁰⁷Pb/²⁰⁶Pb spot ages from 3.95 to 4.19 Ga. Finsch pyropes yield 100 × ?¹⁴²Nd = ? 6 ± 12 ppm, ?¹⁴³Nd = ? 32.5, and ¹⁴⁷Sm/¹⁴⁴Nd = 0.1150. These results do not confirm previous claims for a 30 ppm ¹⁴²Nd excess in South African cratonic mantle. The lack of a ¹⁴²Nd anomaly in these inclusions suggests that isotopic heterogeneities created by early mantle differentiation were remixed at a very fine scale prior to isolation of the South African lithosphere. Alternatively, this result may indicate that only a fraction of the mantle experienced depletion during the first 400 Myr of its history. Analysis of the Jack Hills zircon composite yielded 100 × ?¹⁴²Nd = 8 ± 10 ppm, ?¹⁴³Nd = 45 ± 1, and ¹⁴⁷Sm/¹⁴⁴Nd = 0.5891. Back-calculation of this present-day ?¹⁴³Nd yields an unrealistic estimate for the initial ?¹⁴³Nd of ? 160 ?-units, clearly indicating post-crystallization disturbance of the ¹⁴⁷Sm–¹⁴³Nd system. Examination of ^146,147Sm–^142,143Nd data reveals that the Nd budget of the Jack Hills sample is dominated by non-radiogenic Nd, possibly contained in recrystallized zircon rims or secondary subsurface minerals. This secondary material is characterized by highly discordant U–Pb ages. Although the mass fraction of altered zircon is unlikely to exceed 5–10% of total sample, its high LREE content precludes a reliable evaluation of ¹⁴⁶Sm–¹⁴²Nd systematics in Jack Hills zircons.