Application of stress release models to historical earthquakes from North China

Univariate and multivariate stress release models are fitted to historical earthquake data from North China. It is shown that a better fit is obtained by treating separately the Eastern part of the region, including the North China Plain and Bohai Sea, and the Western part of the region, including the Ordos Plateau and its Eastern boundary. Further improvement is obtained by fitting the large events (M7.6) and smaller events in the Western region by different stress release models. The comparisons are made by computing the likelihoods of the fitted models and discounting the number of parameters used by Akaike's AIC criterion. The models are used to develop long-term risk scenarios for the East and West regions.     

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.     

Focal mechanism inversion in the Giudicarie–Lessini seismotectonic region (Southern Alps, Italy): Insights on tectonic stress and strain

A set of 41 focal mechanisms (1989–2006) from P-wave first polarities is computed from relocated seismic events in the Giudicarie–Lessini region (Southern Alps). Estimated hypocentral depths vary from 3.1 to 20.8 km, for duration magnitudes (M_D) in the range 2.7–5.1. Stress and strain inversions are performed for two seismotectonic zones, namely G (Giudicarie) and L (Lessini). This subdivision is supported by geological evidence, seismicity distribution, and focal mechanism types. The available number of data (16 in G, 22 in L) does not make possible any further subdivisions. Seismotectonic zones G and L are undergoing different kinematic regimes: thrust with strike-slip component in G, and strike-slip in L. Principal stress and strain axes in each sub-region show similar orientations. The direction of maximum horizontal compressive stress is roughly perpendicular to the thrust fronts along the Giudicarie Belt in zone G, and compatible with right-lateral strike-slip reactivation of the faults belonging to the Schio-Vicenza system in zone L. On the whole, kinematic regimes and horizontal stress orientations show a good fit with other stress data from focal mechanisms and breakouts and with geodetic strain rate axes.     

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.     

Simultaneous transport of water and solutes under transient unsaturated flow conditions — A case study

The imbalance between incoming and outgoing salt causes salinization of soils and sub-soils that result in increasing the salinity of stream-flows and agriculture land. This salinization is a serious environmental hazard particularly in semi-arid and arid lands. In order to estimate the magnitude of the hazard posed by salinity, it is important to understand and identify the processes that control salt movement from the soil surface through the root zone to the ground water and stream flows. In the present study, Malaprabha sub-basin (up to dam site) has been selected which has two distinct climatic zones, sub-humid (upstream of Khanapur) and semi-arid region (downstream of Khanapur). In the upstream, both surface and ground waters are used for irrigation, whereas in the downstream mostly groundwater is used. Both soils and ground waters are more saline in downstream parts of the study area. In this study we characterized the soil salinity and groundwater quality in both areas. An attempt is also made to model the distribution of potassium concentration in the soil profile in response to varying irrigation conditions using the SWIM (Soil-Water Infiltration and Movement) model. Fair agreement was obtained between predicted and measured results indicating the applicability of the model.     

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.     

华北克拉通破坏与岩石圈减薄

古太古代（约4.0 Ga）时地球上可能只有一个超级大陆, 它的岩石圈厚度高达400 km。在早元古代,这个超级大陆减薄、裂解成十几块,每块中心是太古宙岩石,边缘是元古宙岩石,且各块厚度不等（150~350km）。从元古宙之后这些被称之为稳定克拉通的大陆岩石圈就一直漂游在地幔软流圈之上。中国华北地块就是这些克拉通之一,与众不同的是它在中生代时遭受了第二次破坏,岩石圈厚度从古生代时的180~200 km 减少到现今的80~100 km。本文作者从流变学的视角出发,围绕华北克拉通破坏和岩石圈减薄这一核心问题,从     

江西浒坑钨矿含矿石英脉的地质特征及成矿构造演化

位于赣西武功山地区的江西浒坑钨矿是一个大型石英脉型钨矿床,与赣南典型的“五层楼”模式不同,矿脉主要集中在岩体内接触带。含矿石英脉除了普通的块状构造以外,还普遍存在条带状构造。详细的野外编录和镜下研究表明含矿块状石英脉在形成后遭受不同程度的韧性剪切,从而形成不同成分的分异条带。剪切带以出现新生的面理和线理（拉伸线理）为特征。岩石学研究表明条带状矿石属于糜棱岩,具有丰富的塑性变形显微构造特征。通过对含矿石英脉的构造解析,认为与岩浆岩有关的原生节理系统基本上是同时形成的,沿节理充填的石英脉经历了长期的、多期次的构造运动,从而形成了复杂的石英脉和矿石类型。建立了原生节理形成→石英充填→韧性剪切发生→二期石英充填→脆性碎裂的矿区构造演化序列。韧性剪切作用的动力学背景可能是受后期岩浆侵入的影响,最后的脆性破裂应与区域上存在南北向的压应力有关。     

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.