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.     

Brine-pool deposition for the Zn–Pb–Cu massive sulphide deposits of the Bathurst mining camp, New Brunswick, Canada. II. Ocean anoxia during mineralisation

Current models of massive sulphide ore genesis in the Bathurst mining camp, New Brunswick, involve settling of sulphide particles from a stagnating, low-salinity hydrothermal plume spreading laterally in an anoxic ocean layer with minimal sulphate content. There is fragmentary evidence of ocean anoxia in the form of local fine lamination in the shales that host some of the deposits but the total organic carbon, S, Fe, U/Th, Ni/Co, V/Ni and V/Cr relationships indicate deposition under oxic or dysoxic conditions. Vanadium and Mn values range from oxic to anoxic and sulphate-reducing to non-sulphate reducing but Mn may be anomalously low due to derivation by erosion of acidic volcanic rocks. The somewhat equivocal physical and chemical data, combined with the likely disturbing effects of penecontemporaneous volcanism, considerably weaken the case for an anoxic bottom layer in a static ocean. The presence of barite with ambient seawater ³⁴S values in Brunswick no. 12 ore, and the abundance of sulphate in modern euxinic basin waters, make a sulphate-free layer unlikely, even if anoxic. Sulphate-bearing, low-salinity fluids mixing with seawater would lead to growth of barite-bearing chimneys and baritic rubble mounds, which are not observed. A model involving brine-pool deposition better explains the major features of the Bathurst ores.     

Delineation of groundwater protection zones by the backward particle tracking method: theoretical background and GIS-based stochastic analysis

The backward particle tracking method, an effective and powerful tool that can be used to delineate groundwater protection zones, is presented. The theoretical background and insights on the applicability of this method are provided. Moreover, the present work enriches the backward particle tracking method with an uncertainty analysis concerning the porosity values, applying a Monte Carlo (MC) approach, coupled with the use of geographical information systems (GIS). As an application example, a wellfield in the Komotini area, Greece, is investigated. The present study may serve as a potential guideline for wellfield delineation, particularly in areas like Greece where lack of data related to the hydrogeological system is often a problem.     

On CO<Subscript>2</Subscript> fluid flow and heat transfer behavior in the subsurface,following leakage from a geologic storage reservoir

Geologic storage of CO₂ is expected to produce plumes of large areal extent, and some leakage may occur along fractures, fault zones, or improperly plugged pre-existing wellbores. A review of physical and chemical processes accompanying leakage suggests a potential for self-enhancement. The numerical simulations presented here confirm this expectation, but reveal self-limiting features as well. It seems unlikely that CO₂ leakage could trigger a high-energy run-away discharge, a so-called “pneumatic eruption,” but present understanding is insufficient to rule out this possibility. The most promising avenue for increasing understanding of CO₂ leakage behavior is the study of natural analogues.     

Structural elements of the Makran region, Oman sea and their potential relevance to tsunamigenisis

The character of convergence along the Arabian–Iranian plate boundary changes radically eastward from the Zagros ranges to the Makran region. This appears to be due to collision of continental crust in the west, in contrast to subduction of oceanic crust in the east. The Makran subduction zone with a length of about 900 km display progressively older and highly deformed sedimentary units northward from the coast, together with an increase in elevation of the ranges. North of the Makran ranges are large subsiding basins, flanked to the north by active volcanoes. Based on 2D seismic reflection data obtained in this study, the main structural provinces and elements in the Gulf of Oman include: (i) the structural elements on the northeastern part of the Arabian Plate and, (ii) the Offshore Makran Accretionary Complex. Based on detailed analysis of these data on the northeastern part of the Arabian Plate five structural provinces and elements—the Musendam High, the Musendam Peneplain, the Musendam Slope, the Dibba Zone, and the Abyssal Plain have been identified. Further, the Offshore Makran Accretionary Complex shown is to consist Accretionary Prism and the For-Arc Basin, while the Accretionary Prism has been subdivided into the Accretionary Wedge and the Accreted/Colored Mélange. Lastly, it is important to note that the Makran subduction zone lacks the trench. The identification of these structural elements should help in better understanding the seismicity of the Makran region in general and the subduction zone in particular. The 1945 magnitude 8.1 tsunamigenic earthquake of the Makran and some other historical events are illustrative of the coastal region’s vulnerability to future tsunami in the area, and such data should be of value to the developing Indian Ocean Tsunami Warning System.     

Ground motion estimation during the Kashmir earthquake of 8th October 2005

In this article, analytical methods have been used to estimate ground motion during the 8 October 2005, Kashmir earthquake. Peak ground acceleration (PGA) values at several stations in the epicentral region have been estimated by empirical analytical source mechanism models. As an alternate analysis, PGA estimates have also been obtained using the stochastic finite fault seismological model. The estimated PGAs are compared with that obtained from damage values. A PGA contour map in the near-source region is provided. It is found that very near to the epicenter, PGA would have reached more than 1 g. It is demonstrated that empirical analytical models can be effectively used to estimate ground motion due to rupture of active faults.     

陆面过程中稳定水同位素逐日变化的数学模拟——引入稳定同位素效应的CLM实例分析

将稳定同位素效应引入CLM(Community Land Model),并对巴西马瑙斯站在平衡年的稳定水同位素的逐日变化进行模拟和分析.结果表明: 降水、水汽和地表径流中δ18O存在明显的季节变化,并与相应的水量存在显著的负相关关系,但凝结物中δ18O与地面凝结量存在显著的正相关关系,蒸发水汽中δ18O与蒸发量之间无显著的相关关系.受土壤贮水削峰功能的影响,表层土壤和根区水中δ18O的季节变化全无.植被层蒸发水汽中稳定同位素的丰度与大气的干湿程度存在密切联系: 当降水量少时,大气干燥,植被层的蒸发较少,植被蒸发中δ18O较高;当降水量较大时,空气湿润,植被层的蒸发量较大,蒸发中δ18O则较低.植被蒸腾中δ18O的变化与源区水体中δ18O的变化保持一致,尤其是与根区水中的δ18O.由于地下径流直接源自根区水的补充,因此,地下径流中δ18O等于根区水中的δ18O.模拟结果还显示,降水MWL (大气水线)的梯度项和常数项均比全球平均MWL略偏小.尽管主要来自降水的贡献,但地表径流和植被层水体的MWLs与降水MWL存在较大的差异,这一方面与两类水体在蒸发过程中的稳定同位素的富集作用有关,另一方面与CLM模拟的水量有关.大气水汽线与降水的MWL的梯度值相近,说明大气水汽与降水近似处于稳定同位素平衡状态.另外,模拟的地面的凝结线与植被层的凝结线均与全球大气水线相近,且具有非常高相关程度,说明CLM的模拟是合理的.     

南印度洋SST与南亚季风环流年代际变化的研究

利用美国NCEP全球大气再分析资料和JONES全球海表面温度异常(SSTA)资料，分析了南印度洋SSTA和南亚季风环流年代际变化的特征。研究发现，无论是南印度洋副热带海水辐合区的SST还是赤道以北非洲西海岸附近上升运动海区的SST的长期变化趋势，除了准3-5年的变化以外，还存在着明显的年代际的变化。对于全球最显著南亚季风环流的分析表明，南亚季风环流也存在明显的年代际时间尺度的变化。与南太平洋SST的年代际变化相比，南印度洋SST的变化周期要相对短一些。通过分析南半球冷空气年代际活动的特征发现，冷空气与南印度洋SST年代际时间尺度的变化具有密切的联系。     

165°E赤道太平洋溶解氧垂直分布异常分析

根据1986—1990中美热带西太平洋海气相互作用8个航次合作综合调查资料，分析了1987年10月赤道中西太平洋165°E（10°N－6°S）次表层水形成溶解氧最大值的原因。1986－1987年E1Nino衰退时期，该海域赤道附近在E1Nino强盛时期下沉的次表层水开始回升，短时期内形成了类似于中、高纬度海域的理化环境，使浮游生物在混合层内聚集生长，最终导致溶解氧含量在次表层出现最大值和过饱和含量。