岩体三维节点偶对分析及其应用

本文讨论了节点偶对分析的三维扩展。节点偶对分析是针对摩擦滑动节理单元的约束边界进行的序列矢量判定。在三维节理单元的分析中,几何约束和相应的力学判据极为复杂,采用节点偶对等效结合力方法才能迎刃而解。文中对三维转向节理模型的分析得到合理的结果,它表明本文所提出的方法可成功地应用于节理岩体的有限元分析。     

Carbonate aeolianites of western Saurashtra,India: experimental decipherment of the depositional mechanisms

Granular carbonate deposits of Late Pleistocene to Early Holocene age, commonly referred to as ‘miliolite limestone’, occur in a linear belt, parallel to the southern coast of Saurashtra, India. In the present study area these carbonate deposits are found in select valleys between ridges and mounds of pyroclastic material present in the Deccan trap plateau. Two different depositional histories have been proposed for these sediments. The presence of marine bioclasts led to the postulation of a marine origin for these deposits. The second school of thought propounded redeposition of the coastal sediments by aeolian processes. Although a few features could not be explained by the proposed aeolian model, critical comparison of these two views favoured the aeolian origin. The mode of occurrence, lithological and structural attributes, and microscopic evidence presented here, also support a possible aeolian origin for these deposits. Experimental observation indicates that these carbonate aeolianites represent backflow deposits, which accumulated because of the flow separation caused by the presence of topographic highs. The conspicuous concave‐up geometry of the deposit conformed to the shape of the separation bulb. In view of the inferred depositional mechanism, the disposition of the deposits and the signature of the palaeoflow direction suggest that the carbonate particles were derived from the north‐western coast of Saurashtra by strong south‐easterly winds. Massive granular carbonates with outsized basement clasts appear to be the product of avalanching of granular material from the higher contours because of oversteepening of the primary deposit.     

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.     

The geochemistry of gem opals as evidence of their origin

Seventy-seven gem opals from ten countries were analyzed by inductively coupled plasma-mass spectrometry (ICP-MS) through a dilution process, in order to establish the nature of the impurities. The results are correlated to the mode of formation and physical properties and are instrumental in establishing the geographical origin of a gem opal. The geochemistry of an opal is shown to be dependant mostly on the host rock, at least for examples from Mexico and Brazil, even if modified by weathering processes. In order of decreasing concentration, the main impurities present are Al, Ca, Fe, K, Na, and Mg (more than 500 ppm). Other noticeable elements in lesser amounts are Ba, followed by Zr, Sr, Rb, U, and Pb. For the first time, geochemistry helps to discriminate some varieties of opals. The Ba content, as well as the chondrite-normalized REE pattern, are the keys to separating sedimentary opals (Ba > 110 ppm, Eu and Ce anomalies) from volcanic opals (Ba < 110 ppm, no Eu or Ce anomaly). The Ca content, and to a lesser extent that of Mg, Al, K and Nb, helps to distinguish gem opals from different volcanic environments. The limited range of concentrations for all elements in precious (play-of-color) compared to common opals, indicates that this variety must have very specific, or more restricted, conditions of formation. We tentatively interpreted the presence of impurities in terms of crystallochemistry, even if opal is a poorly crystallized or amorphous material. The main replacement is the substitution of Si⁴⁺ by Al³⁺ and Fe³⁺. The induced charge imbalance is compensated chiefly by Ca²⁺, Mg²⁺, Mn²⁺, Ba²⁺, K⁺, and Na⁺. In terms of origin of color, greater concentrations of iron induce darker colors (from yellow to “chocolate brown”). This element inhibits luminescence for concentrations above 1000 ppm, whereas already a low content in U (≤ 1 ppm) induces a green luminescence.     

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.     

等离子体质谱法测定玄武岩中微量元素三种样品预处理方法的比较

在测定玄武岩的微量元素时,为了保证得到的微量元素结果反映源区特征,需经预处理去除玄武岩的孔洞中充填的碳酸盐等物质。选取雷州半岛新生代玄武岩三个样品,每一个样品用三种预处理方法进行预处理,等离子体质谱法测定其中的微量元素。分析结果表明,用稀硝酸浸泡样品,会导致玄武岩中大多数微量元素严重丢失,得到的数据失真;而用稀盐酸浸泡去除碳酸盐的方法基本不改变玄武岩中微量元素的含量。     

内蒙古敖包吐萤石矿床的Sr、Nd、Pb同位素地球化学特征

敖包吐萤石矿床是内蒙古北部苏莫查干地区单一萤石矿集区中的一个代表性矿床,产于早二叠世大石寨组火山-沉积岩与早白垩世敖包吐花岗岩的接触带上。文章通过分析该矿床岩、矿石的微量元素和稀土元素,揭示出萤石的成矿作用可分为2个阶段,即交代作用和充填作用。交代作用过程中大石寨组的结晶灰岩可能为萤石的形成提供了部分Ca来源,萤石矿石的稀土元素配分模式与海水基本类似,具有Ce负异常;成矿作用后期主要表现为充填作用,形成颗粒粗大的萤石,表现为重稀土元素富集的特征,并随着萤石的沉淀析出,稀土元素总量逐渐下降,反映出成矿流体经历了较长期的演化过程。各地层单元、花岗岩体和萤石矿石的Sr、Nd、Pb同位素研究表明,萤石的放射性同位素组成具有壳、幔源混合的特点,成矿物质来源具有多源性。早白垩世敖包吐花岗岩可能是萤石中F的主要来源,而大石寨组的结晶灰岩则可能提供了Ca。另外,Pb、Nd同位素的极大不均一性,有可能是成矿流体在运移过程中对艾力格庙群放射性组分的选择性吸收的结果。萤石成矿作用与钾玄岩的时空关系暗示了萤石的成矿过程可能是中国东部岩石圈减薄和下地壳的置换地质事件的结果。在构造转型的过程中,燕山中期富碱的酸性花岗岩浆的活动分异出富含F的成矿流体,与幔源流体混合,沿区域重新活化的深大断裂和大石寨组的层间破碎带上升,交代其间的灰岩透镜体,从而形成敖包吐中型萤石矿床。     

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.     

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.