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.     

Mineralogy and composition of Kuroko deposits from northeastern Honshu and their possible modern analogues from the Izu-Ogasawara (Bonin) Arc south of Japan: Implications for mode of formation

The Kuroko deposits of NE Honshu are a key type deposit for the study of volcanogenic massive sulfide deposits. However, these deposits have not been studied in detail since the early 1980's and knowledge of their mode of formation is now dated. In this study, we present the analysis of 12 samples of the Kuroko deposits, 12 samples of submarine hydrothermal minerals from the Sunrise deposit and 6 samples from Suiyo Seamount, both of which are located on the Izu-Ogasawara (Bonin) Arc, for 27 elements. For the Kuroko deposit, Cd>Sb>Ag>Pb>Hg>As>Zn>Cu are highly enriched, Au>Te>Bi>Ba>Mo are moderately enriched, In>Tl are somewhat enriched and Fe is not significantly enriched relative to the average continental crust. Within each of these deposits, a similar pattern of element associations is apparent: Zn–Pb with As, Sb, Cd, Ag, Hg, Tl and Au; Fe–Cu–Ba with As, Sb, Ag, Tl, Mo, Te and Au; Si–Ba with Ag and Au; CaSO₄. The enrichment of the chalcophilic elements in these deposits is consistent with hydrothermal leaching of these elements from the host rocks which are dominantly rhyolite–dacite in the case of the Kuroko deposits, rhyolite in the case of the Sunrise deposit and dacite–rhyolite in the case of the Suiyo Seamount deposit. However, this pattern of element enrichment is also similar to that observed in fumarolic gas condensates from andesitic volcanoes. This suggests that there may be a significant magmatic contribution to the composition of the hydrothermal fluids responsible for the formation of the Kuroko deposits, although it is not yet possible to quantify the relative contributions of these two sources of elements.The compositional data show that Sunrise and Suiyo Seamount deposits are much closer compositionally to the Kuroko deposits from NE Honshu than are the submarine hydrothermal deposits from the JADE site in the Okinawa Trough which contain, on average, significantly higher concentrations of Pb, Zn, Sb, As and Ag than each of these deposits. In spite of the greater similarity in tectonic setting of the Hokuroku Basin in which the Kuroko deposits formed to the Okinawa Trough (intracontinental rifted back-arc basin) compared to Myojin Knoll and Suiyo Seamount (active arc volcanoes), it appears that submarine hydrothermal deposits from Myojin Knoll and Suiyo Seamount are closer analogues of the Kuroko deposit than are those from the Okinawa Trough. The present data are consistent with the magmatic hydrothermal model for the formation of Kuroko-type deposits as formulated by Urabe and Marumo [Urabe, T., Marumo, K., 1991. A new model for Kuroko-type deposits of Japan. Episodes 14, 246–251].     

Mapping Proterozoic unconformity-related uranium deposits in the Rockhole area, Northern Territory, Australia using landsat ETM+

The Rockhole area, Northern Territory, Australia, hosts a number of Proterozoic unconformity-related uranium deposits. The geology of the area features within Paleoproterozoic rocks of the Pine Creek Orogen, near the unconformity with overlying platform cover sandstone of the Paleo- to Mesoproterozoic McArthur Basin. Landsat Enhanced Thematic Mapper plus (ETM+) data was used in the Rockhole area to assist in mapping geological structures and lithology, and to identify anomalous concentrations of ferrous minerals, the product of alteration, which can be indicators of buried uranium mineralization. Several image-processing procedures were applied to the ETM+ data to identify, isolate and enhance mineralogical information as simple and complex false color composites. ETM+ 754 shown as red green and blue respectively was the best simple image. Overall, complex images based on Principal Component Analysis proved to be the most useful products. Sandstone, shale and siltstone, the target lithologies, Koolpin Formation, the target stratigraphic unit, and bleaching pattern due to the removal of iron(II) compounds, the target alteration pattern, were confidently mapped to provide information required by the mineral emplacement model, which ultimately identified areas of likely uranium mineralization. Thus the contrasting behavior of the two principle oxidation states of uranium and iron can be utilized to map/delineate bleached alteration zones associated with economic concentrations of uranium using multispectral sensors like Landsat or better hyperspectral sensors.     

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.     

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.     

Evaluation of coal bed methane potential of coal seams of Sawang Colliery,Jharkhand, India

The coal seams of Sawang Colliery, East Bokaro Coalfields are bituminous to sub-bituminous in nature and categorized as high gaseous seams (degree II to degree III level). These seams have the potential for coal bed methane (CBM) and their maturity increases with increasing depth, as a result of enhanced pressure-temperature conditions in the underground. The vitrinite maceral group composition of the investigated coal seams ranges from 62.50–83.15%, whereas the inertinite content varies from 14.93–36.81%. The liptinite content varies from 0.66% to 3.09%. The maximum micro-pores are confined within the vitrinite group of macerals. The coal seams exhibit vitrinite reflectance values (Ro% calculated) from 0.94% (sample CG-97) to 1.21% (sample CG-119). Proximate analyses of the investigated coal samples reveal that the moisture content (M%) ranges from 1.28% to 2.98%, whereas, volatile matter (VM%) content is placed in the range of 27.01% to 33.86%. The ash content (A%) ranges from 10.92% to 30.01%. Fixed carbon (FC%) content varies from 41.53% to 55.93%. Fuel ratio variation shows a restricted range from 1.53 to 1.97. All the coal samples were found to be strongly caking and forming coke buttons. The present study is based on the adsorption isotherm experiments carried out under controlled P-T conditions for determination of actual gas adsorption capacity of the coal seams. This analysis shows that the maximum methane gas adsorbed in the coal sample CG-81 is 17 m³/t (Std. daf), at maximum pressure of 5.92 MPa and experimental temperature of 30°C. The calculated Langmuir regression parameters P_L and V_L range from 2.49 to 3.75 MPa and 22.94 to 26.88 m³/t (Std. daf), respectively.     

Proposed seismic hazard maps of Sumatra and Java islands and microzonation study of Jakarta city,Indonesia

This paper presents the development of spectral hazard maps for Sumatra and Java islands, Indonesia and microzonation study for Jakarta city. The purpose of this study is to propose a revision of the seismic hazard map in Indonesian Seismic Code SNI 03-1726-2002. Some improvements in seismic hazard analysis were implemented in the analysis by considering the recent seismic activities around Java and Sumatra. The seismic hazard analysis was carried out using 3-dimension (3-D) seismic source models (fault source model) using the latest research works regarding the tectonic setting of Sumatra and Java. Two hazard levels were analysed for representing 10% and 2% probability of exceedance (PE) in 50 years ground motions for Sumatra and Java. Peak ground acceleration contour maps for those two hazard levels and two additional macrozonation maps for 10% PE in 50 years were produced during this research. These two additional maps represent short period (0.2 s) and long-period (1.0 s) spectra values at the bedrock. Microzonation study is performed in order to obtain ground motion parameters such as acceleration, amplification factor and response spectra at the surface of Jakarta. The analyses were carried out using nonlinear approach. The results were used to develop contour of acceleration at the surface of Jakarta. Finally, the design response spectra for structural design purposes are proposed in this study.     

The WRF model performance for the simulation of heavy precipitating events over Ahmedabad during August 2006

The summer monsoon season of the year 2006 was highlighted by an unprecedented number of monsoon lows over the central and the western parts of India, particularly giving widespread rainfall over Gujarat and Rajasthan. Ahmedabad had received 540.2mm of rainfall in the month of August 2006 against the climatological mean of 219.8mm. The two spells of very heavy rainfall of 108.4mm and 97.7mm were recorded on 8 and 12 August 2006 respectively. Due to meteorological complexities involved in replicating the rainfall occurrences over a region, the Weather Research and Forecast (WRF-ARW version) modeling system with two different cumulus schemes in a nested configuration is chosen for simulating these events. The spatial distributions of large-scale circulation and moisture fields have been simulated reasonably well in this model, though there are some spatial biases in the simulated rainfall pattern. The rainfall amount over Ahmedabad has been underestimated by both the cumulus parameterization schemes. The quantitative validation of the simulated rainfall is done by calculating the categorical skill scores like frequency bias, threat scores (TS) and equitable threat scores (ETS). In this case the KF scheme has outperformed the GD scheme for the low precipitation threshold.     

Implications of Kali-Hindon inter-stream aquifer water balance for groundwater management in western Uttar Pradesh

The Kali-Hindon inter-stream region extends over an area of 395 km² within the Ganga-Yamuna interfluve. It is a fertile tract for sugarcane cultivation. Groundwater is a primary resource for irrigation and industrial purposes. In recent years, over-exploitation has resulted in an adverse impact on the groundwater regime. In this study, an attempt has been made to calculate a water balance for the Kali-Hindon inter-stream region. Various inflows and outflows to and from the aquifer have been calculated. The recharge due to rainfall and other recharge parameters such as horizontal inflow, irrigation return flow and canal seepage were also evaluated. Groundwater withdrawals, evaporation from the water table, discharge from the aquifer to rivers and horizontal subsurface outflows were also estimated. The results show that total recharge into the system is 148.72 million cubic metres (Mcum), whereas the total discharge is 161.06 Mcum, leaving a deficit balance of −12.34 Mcum. Similarly, the groundwater balance was evaluated for the successive four years. The result shows that the groundwater balance is highly sensitive to variation in rainfall followed by draft through pumpage. The depths to water level are shallow in the canal-irrigated northern part of the basin and deeper in the southern part. The pre-monsoon and post-monsoon water levels range from 4.6 to 17.7 m below ground level (bgl) and from 3.5 to 16.5 m bgl respectively. It is concluded that the groundwater may be pumped in the canal-irrigated northern part, while withdrawals may be restricted to the southern portion of the basin, where intense abstraction has led to rapidly falling water table levels.