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.     

Sensitivity analysis of support safety factor for predicting the effects of contributing parameters on roof falls in underground coal mines

Roof falls accounted for 18.18% of all fatal accidents in Indian coal mines, contributing about 35.29% of all fatal accidents in below-ground operations in 2005. The support safety factor, always preferred in support planning and design of underground coal mines, may be an important predictor for roof falls. In this paper, geotechnical data were collected from 14 roof fall incident places in an underground coal mine, located in the Eastern India, which has bord and pillar method of workings. The mean value of probabilistic support safety factor for the case study mine was found to be 1.24. However, the probability, of the estimated support safety factor of less than or equal to one, was found to be 0.246. Sensitivity analysis was conducted to analyze the effects of the contributing parameters on support safety factor and the likelihood of the roof fall. The multi-variate regression analysis was carried out for the data generated by Monte Carlo method to correlate the contributing factors to support safety factor. It ranked gallery width as the first parameter to control the support safety factor.     

The sandstone-hosted stratiform copper mineralization at Mwitapile and its relation to the mineralization at Lufukwe, Lufilian foreland, Democratic Republic of Congo

The Lufilian foreland is a triangular-shaped area located in the SE of the Democratic Republic of Congo and to the NE of the Lufilian arc, which hosts the well-known Central African Copperbelt. The Lufilian foreland recently became an interesting area with several vein-type (e.g., Dikulushi) and stratiform (e.g., Lufukwe and Mwitapile) copper occurrences. The Lufilian foreland stratiform Cu mineralization is, to date, observed in sandstone rock units belonging to the Nguba and Kundelungu Groups (Katanga Supergroup).The Mwitapile sandstone-hosted stratiform Cu prospect is located in the north eastern part of the Lufilian foreland. The host rock for the Cu mineralization is the Sonta Sandstone of the Ngule Subgroup (Kundelungu Group). A combined remote sensing, petrographic and fluid inclusion microthermometric analysis was performed at Mwitapile and compared with similar analysis previously carried out at Lufukwe to present a metallogenic model for the Mwitapile- and Lufukwe-type stratiform copper deposits. Interpretation of ETM+ satellite images for the Mwitapile prospect and the surrounding areas indicate the absence of NE–SW or ENE–WSW faults, similar to those observed controlling the mineralization at Lufukwe. Faults with these orientations are, however, present to the NW, W, SW and E of the Mwitapile prospect. At Mwitapile, the Sonta Sandstone host rock is intensely compacted, arkosic to calcareous with high silica cementation (first generation of authigenic quartz overgrowths). In the Sonta Sandstone, feldspar and calcite are present in disseminated, banded and nodular forms. Intense dissolution of these minerals caused the presence of disseminated rectangular, pipe-like and nodular dissolution cavities. Sulfide mineralization is mainly concentrated in these cavities. The hypogene sulfide minerals consist of two generations of pyrite, chalcopyrite, bornite and chalcocite, separated by a second generation of authigenic quartz overgrowth. The hypogene sulfide minerals are replaced by supergene digenite and covellite. Fluid inclusion microthermometry on the first authigenic quartz phase indicates silica precipitation from an H₂O–NaCl–CaCl₂ fluid with a minimum temperature between 111 and 182 °C and a salinity between 22.0 and 25.5 wt.% CaCl₂ equiv. Microthermometry on the second authigenic quartz overgrowths and in secondary trails related to the mineralization indicate that the mineralizing fluid is characterized by variable temperatures (Th = 120 to 280 °C) and salinities (2.4 to 19.8 wt.% NaCl equiv.) and by a general trend of increasing temperatures with increasing salinities.Comparison between Mwitapile and Lufukwe indicates that the stratiform Cu mineralization in the two deposits is controlled by similar sedimentary, diagenetic and structural factors and likely formed from a similar mineralizing fluid. A post-orogenic timing is proposed for the mineralization in both deposits. The main mineralization controlling factors are grain size, clay and pyrobitumen content, the amount and degree of feldspar and/or calcite dissolution and the presence of NE–SW to ENE–WSW faults. The data support a post-orogenic fluid-mixing model for the Mwitapile- and Lufukwe-type sandstone-hosted stratiform Cu deposits, in which the mineralization is related to the mixing between a Cu-rich hydrothermal fluid, with a temperature up to 280 °C and a maximum salinity of 19.8 wt.% NaCl equiv., with a colder low salinity reducing fluid present in the sandstone host rock. The mineralizing fluid likely migrated upwards to the sandstone source rocks along NE–SW to ENE–WSW orientated faults. At Lufukwe, the highest copper grades at surface outcrops and boreholes were found along and near to these faults. At Mwitapile, where such faults are 2 to 3 km away, the Cu grades are much lower than at Lufukwe. Copper precipitation was possibly promoted by reduction from pre-existing hydrocarbons and non-copper sulfides and by the decrease in fluid salinity and temperature during mixing. Based on this research, new Cu prospects were proposed at Lufukwe and Mwitapile and a set of recommendations for further Cu exploration in the Lufilian foreland is presented.     

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.     

Magnetotelluric investigations for imaging electrical structure of Garhwal Himalayan corridor,Uttarakhand, India

Magnetotelluric investigations have been carried out in the Garhwal Himalayan corridor to delineate the electrical structure of the crust along a profile extending from Indo-Gangetic Plain to Higher Himalayan region in Uttarakhand, India. The profile passing through major Himalayan thrusts: Himalayan Frontal Thrust (HFF), Main Boundary Thrust (MBT) and Main Central Thrust (MCT), is nearly perpendicular to the regional geological strike. Data processing and impedance analysis indicate that out of 44 stations MT data recorded, only 27 stations data show in general, the validity of 2D assumption. The average geoelectric strike, N70°W, was estimated for the profile using tensor decomposition. 2D smooth geoelectrical model has been presented, which provides the electrical image of the shallow and deeper crustal structure. The major features of the model are (i) a low resistivity (<50Ωm), shallow feature interpreted as sediments of Siwalik and Indo-Gangetic Plain, (ii) highly resistive (> 1000Ωm) zone below the sediments at a depth of 6 km, interpreted as the top surface of the Indian plate, (iii) a low resistivity (< 10Ωm) below the depth of 6 km near MCT zone coincides with the intense micro-seismic activity in the region. The zone is interpreted as the partial melting or fluid phase at mid crustal depth. Sensitivity test indicates that the major features of the geoelectrical model are relevant and desired by the MT data.     

Paleosol at the Archean—Proterozoic contact in NW India revisited: Evidence for oxidizing conditions during paleo-weathering?

A number of fine-grained sericite bearing pelitic, schistose lithologies occur along the Archean (Banded Gneiss Complex)-Proterozoic (Aravalli Supergroup) contact (APC) in the Udaipur valley in NW Indian craton. These Al-rich lithologies (subsequently metamorphosed) have been described as ‘paleosols’, developed over a 3.3 Ga old Archean gneissic basement and are overlain by Paleoproterozoic Aravalli quartzite. The paleosol was developed between 2.5 and 2.1, coincident with the globally recognized Great Oxidation Event (GOE). In previous studies these paleosol sections were interpreted to have developed under reducing environment, however, the finding of a ‘ferricrete’ zone in the upper part of Tulsi Namla section (east of Udaipur) during the present study (in addition to earlier reported lithologies) has led to an alternative suggestion of oxygen-rich conditions during paleosol development. The Tulsi Namla paleosol section shows all the features characteristic of a complete paleosol section described from other Archean cratons. The paleosol includes sericite schist with kyanite as the prevalent Al-silicate in the lower part of profile while chloritoid and Fe-oxides typify the Fe-rich upper part. Alumina has remained immobile during the weathering process while Fe and Mn show a decrease in the lower part of the section and an abrupt rise in the upper part, in the ferricrete zone. The field and geochemical data indicate that the Tulsi Namla section is an in situ weathering profile and at least the upper part shows evidence of oxidizing conditions.     

Characteristics of spectral aerosol optical depths over India during ICARB

Spectral aerosol optical depth (AOD) measurements, carried out regularly from a network of observatories spread over the Indian mainland and adjoining islands in the Bay of Bengal and Arabian Sea, are used to examine the spatio-temporal and spectral variations during the period of ICARB (March to May 2006). The AODs and the derived Ångström parameters showed considerable variations across India during the above period. While at the southern peninsular stations the AODs decreased towards May after a peak in April, in the north Indian regions they increased continuously from March to May. The Ångström coefficients suggested enhanced coarse mode loading in the north Indian regions, compared to southern India. Nevertheless, as months progressed from March to May, the dominance of coarse mode aerosols increased in the columnar aerosol size spectrum over the entire Indian mainland, maintaining the regional distinctiveness. Compared to the above, the island stations showed considerably low AODs, so too the northeastern station Dibrugarh, indicating the prevalence of cleaner environment. Long-range transport of aerosols from tshe adjoining regions leads to remarkable changes in the magnitude of the AODs and their wavelength dependencies during March to May. HYSPLIT back-trajectory analysis shows that enhanced long-range transport of aerosols, particularly from the west Asia and northwest coastal India, contributed significantly to the enhancement of AOD and in the flattening of the spectra over entire regions; if it is the peninsular regions and the island Minicoy are more impacted in April, the north Indian regions including the Indo Gangetic Plain get affected the most during May, with the AODs soaring as high as 1.0 at 500 nm. Over the islands, the Ångström exponent (α) remained significantly lower (～1) over the Arabian Sea compared to Bay of Bengal (BoB) (～1.4) as revealed by the data respectively from Minicoy and Port Blair. Occurrences of higher values of α, showing dominance of accumulation mode aerosols, over BoB are associated well with the advection, above the boundary layer, of fine particles from the east Asian region during March and April. The change in the airmass to marine in May results in a rapid decrease in α over the BoB.     

Aerosol mass loading over the marine environment of Arabian Sea during ICARB: Sea-salt and non-sea-salt components

Mass loading and chemical composition of atmospheric aerosols over the Arabian Sea during the pre-monsoon months of April and May have been studied as a part of the Integrated Campaign for Aerosols, gases and Radiation Budget (ICARB). These investigations show large spatial variabilities in total aerosol mass loading as well as that of individual chemical species. The mass loading is found to vary between 3.5 and 69.2 μg m^?3, with higher loadings near the eastern and northern parts of Arabian Sea, which decreases steadily to reach its minimum value in the mid Arabian Sea. The decrease in mass loading from the coast of India towards west is estimated to have a linear gradient of 1.53 μg m^?3/° longitude and an e^?1 scale distance of ～2300 km. SO ₄ ^2? , Cl^? and Na⁺ are found to be the major ionic species present. Apart from these, other dominating watersoluble components of aerosols are NO ₃ ^? (17%) and Ca²⁺ (6%). Over the marine environment of Arabian Sea, the non-sea-salt component dominates accounting to ～76% of the total aerosol mass. The spatial variations of the various ions are examined in the light of prevailing meteorological conditions and airmass back trajectories.     

A technical note on seismic microzonation in the central United States

Microzonation is an effort to evaluate and map potential hazards found in an area, urban area in particular, that could be induced by strong ground shaking during an earthquake. These hazards include: ground motion amplification, liquefaction, and slope failure. The microzonation maps, depicting ground-motion amplification, liquefaction, and landslide potentials, can be produced if the ground motion on bedrock (input) and the site conditions are known. These maps, in combination with ground-motion hazard maps (on bedrock), can be used to develop a variety of hazard mitigation strategies such as seismic risk assessment, emergency response and preparedness, and land-use planning. However, these maps have certain limitations that result from the nature of regional mapping, data limitations, generalization, and computer modeling. These microzonations show that when strong ground shaking occurs, damage is more likely to occur, or be more severe, in the higher hazard areas. The zones shown on the hazard maps should not serve as a substitute for site-specific evaluations.