Groundwater flow regime under natural conditions as inferred from past evidence and contemporary field observations in a semi-arid basin: Cuenca de la Independencia, Guanajuato, México

This study was carried out in the Cuenca de la Independencia, a semi-arid basin in Central Mexico. The objective is to describe the main features of a groundwater flow regime under natural conditions, based on groundwater discharge manifestations. Information obtained from paleoecological, paleontological, archaeological and historical data suggests that, prior to heavy development (starting in the 1950s), the hydrogeologic regime was characterized by a larger groundwater availability in a more humid and colder climate. Manifestations associated to groundwater discharges are springs, lagoons, wetlands, saline soils, chalcedony deposits, phreatophytes, thermalism, and artesianism. The different types of manifestations and their position in the basin indicate the influence of groundwater flow systems hierarchically nested, forming concentric zones at ground level. The groundwater flow regime corresponds to a classical gravity-induced flow system with generation of local, intermediate and regional patterns. Integrating several types of data to establish the flow geometry and its dynamics has proven a useful tool to increase understanding of the original groundwater regimes. This approach can also be applied in other over-exploited semi-arid basins.     

Zebra dolomitization as a result of focused fluid flow in the Rocky Mountains Fold and Thrust Belt, Canada

Discordant zebra dolomite bodies occur locally in the Middle Cambrian Cathedral and Eldon Formations of the Main Ranges of the Canadian Rocky Mountains Fold and Thrust Belt. They are characterized by alternating dark grey (a) and white (b) bands, forming an ‘abba’ diagenetic cyclicity. These bands developed parallel to both bedding and cleavage. Dark grey (a) bands consist of fine (< 300 μm) non-planar crystalline impure dolomite. The white (b) bands are composed of coarse (up to several millimetres) milky-white pure saddle dolomites (b1) which are often covered by pore-lining zoned dolomite (b2). The b phases often possess a saddle-shaped morphology. In contrast to the replacement origin of the a dolomite, the zoned b2 dolomite rims are interpreted as a cement formed in open cavities. The b1 dolomite is interpreted as the result of recrystallization with diagenetic leaching of non-carbonate components. All the zebra dolomites studied are (nearly) stoichiometric and are characterized by enriched Na and depleted Sr concentrations. Fe and Mn concentrations in these dolomites differ depending on the sample locality. Fluid inclusion data indicate that the dolomites formed from relatively hot (T_H = 130–200 °C), saline (20–23 wt% CaCl₂ eq.) fluids. A diagenetic high temperature origin is also supported by depleted δ¹⁸O values (−20 to −14‰ VPDB). A contribution of ⁸⁷Sr-enriched fluids is reflected in the ⁸⁷Sr/⁸⁶Sr values (0·7091–0·7123). Zebra dolomite development is explained by focused fluid flow, which exploited areas of structural weaknesses (e.g. basin-platform, rim areas, faults, etc.). Expulsion of hot basinal brines in a tectonically active regime generated overpressures, which explains the development of secondary porosity during zebra dolomitization as well as the intra-zebra fracturing at decimetre to micrometre scale.     

Mathematical model for hydraulically aided electrokinetic remediation of aquifer and removal of nonanionic copper

One of the most cost-effective in situ technologies for soil and groundwater (i.e., aquifer) remediation is electrokinetic remediation. In electrokinetic remediation, electromigration due to electric field is combined with hydromigration due to hydraulic flow by purge water to remove pollutants from aquifers through the pore water. This study aims at investigating theoretically the role of electromigration (as active movement) of pollutants and the role of hydromigration (as passive movement) of pollutants in electrokinetic remediation, and making it clear that the control variables for electrokinetic remediation are the applied voltage and the hydraulic flow rate. These aims are pursued by construction of a mathematical model based on physico-chemical considerations and by model simulations of the electrokinetic remediation applied to the virtual aquifer polluted by heavy metals of copper sulfate. According to numerical simulations with the model: (1) heavy metal (nonanionic copper) is removed from the upstream anode region and accumulated in the downstream cathode region; (2) to carry away the heavy metal outside the aquifer (global removal), hydromigration by purge water flow is essential; and (3) electromigration contributes mainly to the redistribution of heavy metals within the aquifer (local removal and local accumulation).     

A size classification for debris flows

A 10-fold classification for debris flow size is proposed based on total volume, peak discharge and area inundated by debris. Size classes can be used for regional overview studies where detailed site investigations are either unnecessary, too costly or where the highest hazard and risk creeks need to be identified for further study. They are also useful to compare the regional impact between affected areas and the effects of rainstorms, and they allow lay-people to obtain an understanding of debris flow magnitude and consequences. Finally, different size classes allow the estimation of travel times to points of interest based on empirically derived equations. It is proposed that agencies concerned with debris flows should establish a documentation of debris flow size according to this classification, which serves as a data base for hazard and risk planning.     

泥石流运动衰减动力学初步研究

我国公路泥石流病害严重，泥石流淤埋公路构建筑物是一类常见的公路泥石流病害类型。泥石流衰减动力学是防治泥石流淤埋病害的重要关键技术，也是泥石流运动学、动力学研究的核心问题之一。本文作者运用泥沙运动力学及流体力学原理，初步建立了泥石流固相颗粒和液相浆体的能量衰减条件，把泥石流衰减模式概化为两类，即能量抑制衰减和能量自由衰减；通过泥石流沉积模型试验，得到了不同粘度泥石流体的沉积扇变化形态，随着泥石流体粘度的增大，沉积扇边缘变陡、扩展范围变小、纵轴线长度减小等结论与实际情况吻合；初步建立了泥石流能量衰减速率计算方法。研究成果为防治公路泥石流病害奠定了基础。     

Experimental evidences and numerical modelling of debris flow initiated by channel runoff

Debris flow initiation by channel bed mobilization is a common process in high mountainous areas. Initiation is more likely to occur at the outlet of small, steeply sloping basins where concentrated overland flow feeds an ephemeral channel incised in slope deposits. Such geological conditions are typical of the Dolomite region (Italian Alps), which is characterized by widespread debris flow activity triggered by severe summer thunderstorms. Real-time data and field observations for one of these catchments (Acquabona catchment, Belluno, Italian Alps) were used to characterize the hydrological response of the initiation area to rainfalls of varying intensity and duration. The observed behaviour was then reproduced by means of a simple hydrological model, based on the kinematic wave assumption, to simulate the generation of channel runoff. The model is capable of predicting the observed hydrological response for a wide range of rainfall impulses, thus providing a physical basis for the understanding of the debris flow triggering threshold.     

Temporal and spatial variations in water flow and sediment load in Narmada River Basin, India: natural and man-made factors

The Narmada River flows through the Deccan volcanics and transports water and sediments to the adjacent Arabian Sea. In a first-ever attempt, spatial and temporal (annual, seasonal, monthly and daily) variations in water discharge and sediment loads of Narmada River and its tributaries and the probable causes for these variations are discussed. The study has been carried out with data from twenty-two years of daily water discharge at nineteen locations and sediment concentrations data at fourteen locations in the entire Narmada River Basin. Water flow in the river is a major factor influencing sediment loads in the river. The monsoon season, which accounts for 85 to 95% of total annual rainfall in the basin, is the main source of water flow in the river. Almost 85 to 98% of annual sediment loads in the river are transported during the monsoon season (June to November). The average annual sediment flux to the Arabian Sea at Garudeshwar (farthest downstream location) is 34.29×10⁶ t year⁻¹ with a water discharge of 23.57 km³ year⁻¹. These numbers are the latest and revised estimates for Narmada River. Water flow in the river is influenced by rainfall, catchment area and groundwater inputs, whereas rainfall intensity, geology/soil characteristics of the catchment area and presence of reservoirs/dams play a major role in sediment discharge. The largest dam in the basin, namely Sardar Sarovar Dam, traps almost 60–80% of sediments carried by the river before it reaches the Arabian Sea.     

Mixing of event and pre‐event water in a shallow Entisol in sloping farmland based on isotopic and hydrometric measurements,SW China

Water percolation and flow processes in subsurface geologic media play an important role in determining the water source for plants and the transport of contaminants or nutrients, which is essential for water resource management and the development of measures for pollution mitigation. During June 2013, the dynamics of the rainwater, soil water, subsurface flows and groundwater in a shallow Entisol on sloping farmland were monitored using a hydrometric and isotopic approach. The results showed that effective mixing of rainwater and soil water occurred in hours. The rebound phenomenon of δD profiles in soils showed that most isotope‐depleted rainwater largely bypassed the soil matrix when the water saturation in the soil was high. Preferential‐flow, which was the dominant water movement pattern in the vadose zone, occurred through the whole soil profile, and infrequent piston‐flow was mainly found at 20–40 cm in depth. The interflow in the soil layer, composed of 75.2% rainwater, was only generated when the soil profile had been saturated. Underflow in the fractured mudrock was the dominant flow type in this hillslope, and outflow was dominated by base flow (groundwater flow) with a mean contribution of 76.7%. The generation mechanism of underflow was groundwater ridging, which was superimposed upon preferential‐flow composed mainly of rainwater. The quick mixing process of rainwater and soil water and the rapid movement of the mixture through preferential channels in the study soil, which shows a typical bimodal pore size distribution, can explain the prompt release of pre‐event water in subsurface flow. Water sources of subsurface flows at peak discharge could be affected by the antecedent soil water content, rain characteristics and antecedent groundwater levels. Copyright © 2016 John Wiley & Sons, Ltd.     

Investigation into the effect of fault properties on wave transmission

The estimation of wave transmission across the fractured rock masses is of great importance for rock engineers to assess the stability of rock slopes in open pit mines. Presence of fault, as a major discontinuity, in the jointed rock mass can significantly impact on the peak particle velocity and transmission of blast waves, particularly where a fault contains a thick infilling with weak mechanical properties. This paper aims to study the effect of fault properties on transmission of blasting waves using the distinct element method. First, a validation study was carried out on the wave transmission across a single joint and different rock mediums through undertaking a comparative study against analytical models. Then, the transmission of blast wave across a fault with thick infilling in the Golgohar iron mine, Iran, was numerically studied, and the results were compared with the field measurements. The blast wave was numerically simulated using a hybrid finite element and finite difference code which then the outcome was used as the input for the distinct element method analysis. The measured uplift of hanging wall, as a result of wave transmission across the fault, in the numerical model agrees well with the recorded field measurement. Finally, the validated numerical model was used to study the effect of fault properties on wave transmission. It was found that the fault inclination angle is the most effective parameter on the peak particle velocity and uplift. Copyright © 2017 John Wiley & Sons, Ltd.