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.     

Fluid generation and evolution during exhumation of deeply subducted UHP continental crust: Petrogenesis of composite granite–quartz veins in the Sulu belt,China

Composite granite–quartz veins occur in retrogressed ultrahigh pressure (UHP) eclogite enclosed in gneiss at General's Hill in the central Sulu belt, eastern China. The granite in the veins has a high‐pressure (HP) mineral assemblage of dominantly quartz+phengite+allanite/epidote+garnet that yields pressures of 2.5–2.1 GPa (Si‐in‐phengite barometry) and temperatures of 850–780°C (Ti‐in‐zircon thermometry) at 2.5 GPa (~20°C lower at 2.1 GPa). Zircon overgrowths on inherited cores and new grains of zircon from both components of the composite veins crystallized at c. 221 Ma. This age overlaps the timing of HP retrograde recrystallization dated at 225–215 Ma from multiple localities in the Sulu belt, consistent with the HP conditions retrieved from the granite. The ε_Hf(t) values of new zircon from both components of the composite veins and the Sr–Nd isotope compositions of the granite consistently lie between values for gneiss and eclogite, whereas δ¹⁸O values of new zircon are similar in the veins and the crustal rocks. These data are consistent with zircon growth from a blended fluid generated internally within the gneiss and the eclogite, without any ingress of fluid from an external source. However, at the peak metamorphic pressure, which could have reached 7 GPa, the rocks were likely fluid absent. During initial exhumation under UHP conditions, exsolution of H₂O from nominally anhydrous minerals generated a grain boundary supercritical fluid in both gneiss and eclogite. As exhumation progressed, the volume of fluid increased allowing it to migrate by diffusing porous flow from grain boundaries into channels and drain from the dominant gneiss through the subordinate eclogite. This produced a blended fluid intermediate in its isotope composition between the two end‐members, as recorded by the composite veins. During exhumation from UHP (coesite) eclogite to HP (quartz) eclogite facies conditions, the supercritical fluid evolved by dissolution of the silicate mineral matrix, becoming increasingly solute‐rich, more ‘granitic’ and more viscous until it became trapped. As crystallization began by diffusive loss of H₂O to the host eclogite concomitant with ongoing exhumation of the crust, the trapped supercritical fluid intersected the solvus for the granite–H₂O system, allowing phase separation and formation of the composite granite–quartz veins. Subsequently, during the transition from HP eclogite to amphibolite facies conditions, minor phengite breakdown melting is recorded in both the granite and the gneiss by K‐feldspar+plagioclase+biotite aggregates located around phengite and by K‐feldspar veinlets along grain boundaries. Phase equilibria modelling of the granite indicates that this late‐stage melting records P–T conditions towards the end of the exhumation, with the subsolidus assemblage yielding 0.7–1.1 GPa at <670°C. Thus, the composite granite–quartz veins represent a rare example of a natural system recording how the fluid phase evolved during exhumation of continental crust. The successive availability of different fluid phases attending retrograde metamorphism from UHP eclogite to amphibolite facies conditions will affect the transport of trace elements through the continental crust and the role of these fluids as metasomatic agents interacting with the mantle wedge in the subduction channel.     

Impacts of groundwater flow on the evolution of a thermokarst lake in the permafrost–dominated region on the Qinghai–Tibet plateau

As a result of global warming induced permafrost degradation in recent decades, thermokarst lakes in the Qinghai–Tibet plateau (QTP) have been regulating local hydrological and ecological processes. Simulations with coupled moisture–heat numerical models in the Beiluhe basin (located in the hinterland of permafrost regions on the QTP) have provided insights into the interaction between groundwater flow and the freeze–thaw process. A total of 30 modified SUTRA scenarios were established to examine the effects of hydrodynamic forces, permeability, and climate on thermokarst lakes. The results indicate that the hydrodynamic condition variables regulate the permafrost degradation around the lakes. In case groundwater recharges to the lake, a low–temperature groundwater flow stimulates the expansion of the surrounding thawing regions through thermal convection. The thawing rate of the permafrost underlying the lake intensifies when groundwater is discharged from the lake. Under different permeability conditions, spatiotemporal variations in the active layer thickness significantly influence the occurrence of an open talik at the lake bottom. A warmer and wetter climate will inevitably lead to a sharp decrease in the upper limit of the surrounding permafrost, with a continual decrease in the duration of open talik events. Overall, our results underscore that comprehensive consideration of the relevant hydrologic processes is critical for improving the understanding of environmental and ecological changes in cold environments.     

An Eocene/Oligocene blueschist‐/greenschist facies P–T loop from the Cycladic Blueschist Unit on Naxos Island,Greece: Deformation‐related re‐equilibration vs. thermal relaxation

Geothermobarometric and geochronological work indicates a complete Eocene/early Oligocene blueschist/greenschist facies metamorphic cycle of the Cycladic Blueschist Unit on Naxos Island in the Aegean Sea region. Using the average pressure–temperature (P–T) method of thermocalc coupled with detailed textural work, we separate an early blueschist facies event at 576 ± 16 to 619 ± 32°C and 15.5 ± 0.5 to 16.3 ± 0.9 kbar from a subsequent greenschist facies overprint at 384 ± 30°C and 3.8 ± 1.1 kbar. Multi‐mineral Rb–Sr isochron dating yields crystallization ages for near peak‐pressure blueschist facies assemblages between 40.5 ± 1.0 and 38.3 ± 0.5 Ma. The greenschist facies overprint commonly did not result in complete resetting of age signatures. Maximum ages for the end of greenschist facies reworking, obtained from disequilibrium patterns, cluster near c. 32 Ma, with one sample showing rejuvenation at c. 27 Ma. We conclude that the high‐P rocks from south Naxos were exhumed to upper mid‐crustal levels in the late Eocene and early Oligocene at rates of 7.4 ± 4.6 km/Ma, completing a full blueschist‐/greenschist facies metamorphic cycle soon after subduction within c. 8 Ma. The greenschist facies overprint of the blueschist facies rocks from south Naxos resulted from rapid exhumation and associated deformation/fluid‐controlled metamorphic re‐equilibration, and is unrelated to the strong high‐T metamorphism associated with the Miocene formation of the Naxos migmatite dome. It follows that the Miocene thermal overprint had no impact on rock textures or Sr isotopic signatures, and that the rocks of south Naxos underwent three metamorphic events, one more than hitherto envisaged.     

Transient seepage analysis in zoned anisotropic soils based on the scaled boundary finite‐element method

The scaled boundary finite‐element method, a semi‐analytical computational scheme primarily developed for dynamic stiffness of unbounded domains, is applied to the analysis of unsteady seepage flow problems. This method is based on the finite‐element technology and gains the advantages of the boundary element method as well. Only boundary of the domain is discretized, no fundamental solution is required and singularity problems can be modeled rigorously. Anisotropic and non‐homogeneous materials satisfying similarity are modeled with no additional efforts. In this study, firstly, formulation of the method for the transient seepage flow problems is derived followed by its solution procedures. The accuracy, simplicity and applicability of the method are demonstrated via four numerical examples of transient seepage flow – three of them are available in the literature. Homogenous, non‐homogenous, isotropic and anisotropic material properties are considered to show the versatility of the technique. Excellent agreement with the finite‐element method is observed. The method out‐performs the finite‐element method in modeling singularity points. Copyright © 2014 John Wiley & Sons, Ltd.     

Assessment of flow paths and confluences for saltwater intrusion in a deltaic river network

Saltwater intrusion is a serious issue in estuarine deltas all over the world due to rapid urban sprawl and water shortage. Therefore, detecting the major flow paths or locations at risk of saltwater intrusion in estuarine ecosystems is important for mitigating saltwater intrusion. In this paper, we introduce a centrality index, the betweenness centrality (BC), to address this problem. Using the BC as the weighted attribute of the river network, we identify the critical confluences for saltwater intrusion and detect the preferential flow paths for saltwater intrusion through the least‐cost‐path algorithm from a graph theory approach. Moreover, we analyse the responses of the BC values of confluences calculated in the river network to salinity. Our results show that the major flow paths and critical confluences for saltwater intrusion in a deltaic river network can be represented by the least cost paths and the BC values of confluences, respectively. In addition, a significant positive correlation between the BC values of confluences and salinity is determined in the Pearl River Delta. Changes in the salinity can produce significant variation in the BC values of confluences. Therefore, freshwater can be diverted into these major flow paths and critical confluences to improve river network management under saltwater intrusion. Copyright © 2015 John Wiley & Sons, Ltd.     

Numerical simulation of percussive drilling

This paper presents a novel dynamical model to analyze the long‐term response of a percussive drilling system. This departs from existing approaches that usually consider a single activation and bit/rock interaction cycle for the analysis of the process performance. The proposed model integrates the axial dynamics of an elastic piston and an elastic drill bit, a motion‐dependent pressure law to drive the piston, and a generalized bit/rock interaction law representative of the dynamic indentation taking place at the bit/rock interface. It applies to down‐the‐hole percussive drilling as well as top‐hole, with minor modifications. The model does not account for the angular motion or the hole cleaning, however. The model is first formulated mathematically; then, a finite‐dimensional approximation is proposed for computations. Numerical analyses of the model response, for a low‐size down‐the‐hole percussive system, follow. The period‐1 stationary response for the reference configuration is studied in detail, and parametric analyses assessing the influence on the rate of penetration of the bit/rock interaction parameters, the feed force, and the percussive activation parameters are conducted. These analyses reveal that the multiscale nature of the process is well captured by the model and recover expected trends for the influence of the parameters. They also suggest that a significant increase of the penetration rate can be achieved by increasing the percussive frequency. Copyright © 2015 John Wiley & Sons, Ltd.     

Nine thousand years of upper montane soil/vegetation dynamics from the summit of Caratuva Peak,Southern Brazil

Biodiversity loss, climate change, and increased freshwater consumption are some of the main environmental problems on Earth. Mountain ecosystems can reduce these threats by providing several positive influences, such as the maintenance of biodiversity, water regulation, and carbon storage, amongst others. The knowledge of the history of these environments and their response to climate change is very important for management, conservation, and environmental monitoring programs. The genesis of the soil organic matter of the current upper montane vegetation remains unclear and seems to be quite variable depending on location. Some upper montane sites in the very extensive coastal Sea Mountain Range present considerable organic matter from the late Pleistocene and other from only the Holocene. Our study was carried out on three soil profiles (two cores in grassland and one in forest) on the Caratuva Peak of the Serra do Ibitiraquire (a sub-range of Sea Mountain Range – Serra do Mar) in Southern Brazil. The δ¹³C isotopic analyses of organic matter in soil horizons were conducted to detect whether C₃ or C₄ plants dominated the past communities. Complementarily, we performed a pollen analysis and ¹⁴C dating of the humin fraction to obtain the age of the studied horizons. Except for a short and probably drier period (between 6000 and 4500 cal yr BP), C₃ plants, including ombrophilous grasses and trees, have dominated the highlands of the Caratuva Peak (Pico Caratuva), as well as the other uppermost summits of the Serra do Ibitiraquire, since around 9000 cal yr BP. The Caratuva region represents a landscape of high altitude grasslands (campos de altitude altomontanos or campos altomontanos) and upper montane rain/cloud forests with soils that most likely contain some organic matter from the late Pleistocene, as has been reported in Southern and Southeastern Brazil for other sites. However, our results indicate that the studied deposits (near the summit) are from the early to late Holocene, when somewhat wetter and warmer conditions (since around 9000 cal yr BP) enabled a stronger colonization of the ridge of Pico Caratuva by mainly C₃ plants, especially grassland species. However, at the same time, even near the summit, the soil core from the forest site already presented the current physiognomy (or a shrubby/elfin or successional forest), indicating that the colonization of the neighboring uppermost saddles and valleys were probably populated mainly by upper montane forest species.