Combined mine tremors source location and error evaluation in the Lubin Copper Mine (Poland)

A modified method of mine tremors location used in Lubin Copper Mine is presented in the paper. In mines where an intensive exploration is carried out a high accuracy source location technique is usually required. The effect of the flatness of the geophones array, complex geological structure of the rock mass and intense exploitation make the location results ambiguous in such mines. In the present paper an effective method of source location and location's error evaluations are presented, combining data from two different arrays of geophones. The first consists of uniaxial geophones spaced in the whole mine area. The second is installed in one of the mining panels and consists of triaxial geophones. The usage of the data obtained from triaxial geophones allows to increase the hypocenter vertical coordinate precision. The presented two-step location procedure combines standard location methods: P-waves directions and P-waves arrival times. Using computer simulations the efficiency of the created algorithm was tested. The designed algorithm is fully non-linear and was tested on the multilayered rock mass model of the Lubin Copper Mine, showing a computational better efficiency than the traditional P-wave arrival times location algorithm. In this paper we present the complete procedure that effectively solves the non-linear location problems, i.e. the mine tremor location and measurement of the error propagation.     

U–Pb SHRIMP zircon geochronology and T–t–d history of the Kampa Dome, southern Tibet

Structural, petrographic and geochronologic studies of the Kampa Dome provide insights into the tectonothermal evolution of orogenic crust exposed in the North Himalayan gneiss domes of southern Tibet. U–Pb ion microprobe dating of zircons from granite gneiss exposed at the deepest levels within the dome yields concordia ²⁰⁶Pb/²³⁸U age populations of 506 ± 3 Ma and 527 ± 6 Ma, with no evidence of new zircon growth during Himalayan orogenesis. However, the granite contains penetrative deformation fabrics that are also preserved in the overlying Paleozoic strata, implying that the Kampa granite is a Cambrian pluton that was strongly deformed and metamorphosed during Himalayan orogenesis. Zircons from deformed leucogranite sills that cross-cut Paleozoic metasedimentary rocks yield concordant Cambrian ages from oscillatory zoned cores and discordant ages ranging from ca. 491–32 Ma in metamict grains. Since these leucogranites clearly post-date the metasedimentary rocks they intrude, the zircons are interpreted as xenocrysts that are probably derived from the Kampa granite. The Kampa Dome formed via a series of progressive orogenic events including regional ~ N–S contraction and related crustal thickening (D₁), predominately top-to-N ductile shearing and crustal extension (D₂), top-to-N brittle–ductile faulting and related folding on the north limb of the dome, localized top-to-S faulting on the southern limb of the dome, and crustal doming (D₃), and continued N–S contraction, E–W extension and doming (D₄). Structural and geochronologic variability amongst adjacent North Himalayan gneiss domes may reflect changes in the magnitude of crustal exhumation along the North Himalayan antiform, possibly relating to differences in the mid-crustal geometry of the exhuming fault systems.     

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.     

Use of sinkhole and specific capacity distributions to assess vertical gradients in a karst aquifer

The carbonate-rock aquifer in the Great Valley, West Virginia, USA, was evaluated using a database of 687 sinkholes and 350 specific capacity tests to assess structural, lithologic, and topographic influences on the groundwater flow system. The enhanced permeability of the aquifer is characterized in part by the many sinkholes, springs, and solutionally enlarged fractures throughout the valley. Yet, vertical components of subsurface flow in this highly heterogeneous aquifer are currently not well understood. To address this problem, this study examines the apparent relation between geologic features of the aquifer and two spatial indices of enhanced permeability attributed to aquifer karstification: (1) the distribution of sinkholes and (2) the occurrence of wells with relatively high specific capacity. Statistical results indicate that sinkholes (funnel and collapse) occur primarily along cleavage and bedding planes parallel to subparallel to strike where lateral or downward vertical gradients are highest. Conversely, high specific capacity values are common along prominent joints perpendicular or oblique to strike. The similarity of the latter distribution to that of springs suggests these fractures are areas of upward-convergent flow. These differences between sinkhole and high specific capacity distributions suggest vertical flow components are primarily controlled by the orientation of geologic structure and associated subsurface fracturing.     

Evaluation of approaches to calculate debris-flow parameters for hazard assessment

Many different runout prediction methods can be applied to estimate the mobility of future debris flows during hazard assessment. The present article reviews the empirical, analytical, simple flow routing and numerical techniques. All these techniques were applied to back-calculate a debris flow, which occurred in 1982 at La Guingueta catchment, in the Eastern Pyrenees. A sensitivity analysis of input parameters was carried out, while special attention was paid to the influence of rheological parameters. We used the Voellmy fluid rheology for our analytical and numerical modelling, since this flow resistance law coincided best with field observations. The simulation results indicated that the “basal” friction coefficients rather affect the runout distance, while the “turbulence” terms mainly influence flow velocity. A comparison of the velocity computed on the fan showed that the analytical model calculated values similar to the numerical ones. The values of our rheological parameters calibrated at La Guingueta agree with data back-calculated for other debris flows. Empirical relationships represent another method to estimate total runout distance. The results confirmed that they contain an important uncertainty and they are strictly valid only for the conditions, which were the basis for their development. With regards to the simple flow routing algorithm, this methods could satisfactorily simulate the total area affected by the 1982 debris flow, but it was not able to directly calculate total runout distance and velocity. Finally, a suggestion on how different runout prediction methods can be applied to generate debris-flow hazard maps is presented. Taking into account the definition of hazard and intensity, the best choice would be to divide the resulting hazard maps into two types: “final hazard maps” and “preliminary hazard maps”. Only the use of numerical models provided final hazard maps, because they could incorporate different event magnitudes and they supplied output-values for intensity calculation. In contrast, empirical relationships and flow routing algorithms, or a combination of both, could be applied to create preliminary hazard maps. The present study only focussed on runout prediction methods. Other necessary tasks to complete the hazard assessment can be looked up in the “Guidelines for landslide susceptibility, hazard and risk zoning” included in this Special Issue.     

Assessment of disaster resilience capacity of hillslope communities with high risk for geological hazards

This study presents a novel preparedness assessment method for assessing hazard mitigation and environmental planning of hillslope communities. A professional questionnaire was utilized to weight each indicator. Communities in Hsinchu, Taichung and Nantou counties with debris flow hazards were taken as study samples. Debris flow risk and landslide susceptibility for each community were determined using Geographic Information System (GIS) technology and logistic regression analysis. Thus, a novel risk assessment method for evaluating disaster resilience capacity of hillslope communities was established. This method was then applied to assess casualties caused by Typhoon Herb in 1996 and Typhoon Mindulle in 2004. Additionally, the analytical results generated by this assessment method were discussed with the aim of developing references for implementation of risk analysis, increasing the effectiveness of disaster mitigation, and reducing future loss of life and property.     

松南海坨子—大布苏地区泉四段— 姚家组层序地层研究

根据工区内现有的钻井、测井和地震资料,应用以多级次基准面旋回为参照面的高分辨层序地层学理论与研究方法,对松辽盆地南部海坨子—大布苏地区的泉四段—姚家组地层进行了详细的层序地层研究,共识别出5个长期基准面旋回（相当于三级层序）,9个中期基准面旋回（相当于四级层序）,建立了研究区的层序划分方案,在连井剖面的层序地层对比基础上,建立了该区的地层格架,并对层序构成特征进行了分析,研究了储集层在层序格架内的分布规律,揭示出浊积砂体和三角洲前缘相的分流河道形成于基准面上升期,三角洲前缘相的河口坝、远砂坝和席状砂沉积形成于基准面下降期,指出浊积砂体和三角洲前缘相的远砂坝、席状砂沉积是研究区形成岩性油气藏的有利储集体。     

基于群桩基础沉降估算确定压缩层厚度的探讨

在现行规范中,由于桩基沉降估算模式的不同,压缩层厚度的确定方法也不完全相同。对这些方法进行归纳、总结,并通过具体的工程实例,对由不同估算模式确定的压缩层厚度及沉降估算结果进行探讨与分析,认为对于工程场地地质资料掌握较为翔实的情况下,压缩层厚度可采用变形比法确定,否则,宜根据应力比法确定。     

中国大别山东南缘首次发现大坝陨坑构造

简介首次在大别山东南缘安徽省境内发现的大坝陨坑构造。经初步评价,知陨坑呈椭圆形,长轴呈北北东向,长约19km,短轴长约12km,最大坑深约2km,是一个有中央隆起区的复杂型陨坑。在卫片上陨坑显示环形影像,地貌形态为一洼地。陨坑基座保存尚好,可对它直观和追索陨坑构造边界。形成于230Ma左右的各类撞击变质岩石系列齐全,其中含有柯石英等典型的撞击变质矿物及撞击碎理等超微构造,特别是在陨坑基座内壁普遍发现有鉴别陨坑构造最可靠的标志——干裂自角砾岩,都证明大坝环形影像是一个典型的陨星撞击坑构造。它的发现,具有很大的科学及经济意义,对今后褶皱山区寻找和研究陨坑构造具有示范和指导作用,同时大大丰富了建设大别山世界地质公园的地质依据。无庸置疑此发现将促进大别山旅游业的发展及陨星撞击科普知识的传播。     

Numerical Experiment on Two-Dimensional Line Thermal

The time evolution of a two-dimensional line thermal-a turbulent flow produced by an initial element with signifi-cant buoyancy released in a large water body, is numerically studied with the two-equation k - e model for turbulence closure. The numerical results show that the thermal is characterized by a vortex pair flow and a kidney shaped concentra-tion structure with double peak maxima; the computed flow details and scalar mixing characteristics can be described by self-similar relations beyond a dimensionless time around 10. There are two regions in the flow field of a line thermal: a mixing region where the concentration of tracer fluid is high and the flow is turbulent and rotational with a pair of vortex eyes, and an ambient region where the concentration is zero and the flow is potential and well-described by a model of doublet with strength very close to those given by early experimental and analytical studies. The added virtual mass coeffi-cient of the thermal motion is found to be approximat