Numerical simulation on mining-induced water inrushes related to geologic structures using a damage-based hydromechanical model

A large number of statistics indicate that water inrush has a direct relationship with geological structures such as fault and karst collapse columns. Understanding the mechanism of water inrushes controlled by geologic structures is of vital importance for adopting effective measures to prevent their occurrence. The work begins with formulization of a damage-based hydromechanical model based on elastic damage theory. Next, the model is numerically implemented with finite element method by employing a finite element package called COMSOL Multiphysics, and is also validated against some existing experimental observations. Finally, the model is used to simulate the mining-induced groundwater inrushes when the effect of faults and karst collapse columns is considered in the numerical simulation, and some suggestive conclusions for preventing water inrushes and optimizing underground mining operations are drawn.     

Using monazite and zircon petrochronology to constrain the P–T–t evolution of the middle crust in the Bhutan Himalaya

The growth and dissolution behaviour of accessory phases (and especially those of geochronological interest) in metamorphosed pelites depends on, among others, the bulk composition, the prograde metamorphic evolution and the cooling path. Monazite and zircon are arguably the most commonly used geochronometers for dating felsic metamorphic rocks, yet crystal growth mechanisms as a function of rock composition, pressure and temperature are still incompletely understood. Ages of different growth zones in zircon and monazite in a garnet‐bearing anatectic metapelite from the Greater Himalayan Sequence in NW Bhutan were investigated via a combination of thermodynamic modelling, microtextural data and interpretation of trace‐element chemical ‘fingerprint’ indicators in order to link them to the metamorphic stage at which they crystallized. Differences in the trace‐element composition (HREE, Y, Eu_N/Eu*_N) of different phases were used to track the growth/dissolution of major (e.g. plagioclase, garnet) and accessory phases (e.g. monazite, zircon, xenotime, allanite). Taken together, these data constrain multiple pressure–temperature–time (P–T–t) points from low temperature (<550 °C) to upper amphibolite facies (partial melting, >700 °C) conditions. The results suggest that the metapelite experienced a cryptic early metamorphic stage at c. 38 Ma at <550 °C, ≥0.85 GPa during which plagioclase was probably absent. This was followed by a prolonged high‐T, medium‐pressure (~600 °C, 0.55 GPa) evolution at 35–29 Ma during which the garnet grew, and subsequent partial melting at >690 °C and >18 Ma. Our data confirm that both geochronometers can crystallize independently at different times along the same P–T path and that neither monazite nor zircon necessarily provides timing constraints on ‘peak’ metamorphism. Therefore, collecting monazite and zircon ages as well as major and trace‐element data from major and accessory phases in the same sample is essential for reconstructing the most coherent metamorphic P–T–t evolution and thus for robustly constraining the rates and timescales of metamorphic cycles.     

Constraints on the timing and conditions of high‐grade metamorphism,charnockite formation and fluid–rock interaction in the Trivandrum Block,southern India

Incipient charnockites have been widely used as evidence for the infiltration of CO₂‐rich fluids driving dehydration of the lower crust. Rocks exposed at Kakkod quarry in the Trivandrum Block of southern India allow for a thorough investigation of the metamorphic evolution by preserving not only orthopyroxene‐bearing charnockite patches in a host garnet–biotite felsic gneiss, but also layers of garnet–sillimanite metapelite gneiss. Thermodynamic phase equilibria modelling of all three bulk compositions indicates consistent peak‐metamorphic conditions of 830–925 °C and 6–9 kbar with retrograde evolution involving suprasolidus decompression at high temperature. These models suggest that orthopyroxene was most likely stabilized close to the metamorphic peak as a result of small compositional heterogeneities in the host garnet–biotite gneiss. There is insufficient evidence to determine whether the heterogeneities were inherited from the protolith or introduced during syn‐metamorphic fluid flow. U–Pb geochronology of monazite and zircon from all three rock types constrains the peak of metamorphism and orthopyroxene growth to have occurred between the onset of high‐grade metamorphism at c. 590 Ma and the onset of melt crystallization at c. 540 Ma. The majority of metamorphic zircon growth occurred during protracted melt crystallization between c. 540 and 510 Ma. Melt crystallization was followed by the influx of aqueous, alkali‐rich fluids likely derived from melts crystallizing at depth. This late fluid flow led to retrogression of orthopyroxene, the observed outcrop pattern and to the textural and isotopic modification of monazite grains at c. 525–490 Ma.     

A revised evolutionary history of armadillos (Dasypus) in North America based on ancient mitochondrial DNA

The large, beautiful armadillo, Dasypus bellus, first appeared in North America about 2.5 million years ago, and was extinct across its southeastern US range by 11 thousand years ago (ka). Within the last 150 years, the much smaller nine‐banded armadillo, D. novemcinctus, has expanded rapidly out of Mexico and colonized much of the former range of the beautiful armadillo. The high degree of morphological similarity between these two species has led to speculation that they might be a single, highly adaptable species with phenotypical responses and geographical range fluctuations resulting from environmental changes. If this is correct, then the biology and tolerance limits for D. novemcinctus could be directly applied to the Pleistocene species, D. bellus. To investigate this, we isolated ancient mitochondrial DNA from late Pleistocene‐age specimens of Dasypus from Missouri and Florida. We identified two genetically distinct mitochondrial lineages, which most likely correspond to D. bellus (Missouri) and D. novemcinctus (Florida). Surprisingly, both lineages were isolated from large specimens that were identified previously as D. bellus. Our results suggest that D. novemcinctus, which is currently classified as an invasive species, was already present in central Florida around 10 ka, significantly earlier than previously believed.     

A long‐term study of stable isotopes as tracers of processes governing water flow and quality in a lowland river basin: the upper Thames,UK

A long‐term study of O, H and C stable isotopes has been undertaken on river waters across the 7000‐km² upper Thames lowland river basin in the southern UK. During the period, flow conditions ranged from drought to flood. A 10‐year monthly record (2003–2012) of the main River Thames showed a maximum variation of 3‰ (δ¹⁸O) and 20‰ (δ²H), although interannual average values varied little around a mean of –6.5‰ (δ¹⁸O) and –44‰ (δ²H). A δ²H/δ¹⁸O slope of 5.3 suggested a degree of evaporative enrichment, consistent with derivation from local rainfall with a weighted mean of –7.2‰ (δ¹⁸O) and –48‰ (δ²H) for the period. A tendency towards isotopic depletion of the river with increasing flow rate was noted, but at very high flows (>100 m³/s), a reversion to the mean was interpreted as the displacement of bank storage by rising groundwater levels (corroborated by measurements of specific electrical conductivity). A shorter quarterly study (October 2011–April 2013) of isotope variations in 15 tributaries with varying geology revealed different responses to evaporation, with a well‐correlated inverse relationship between Δ¹⁸O and baseflow index for most of the rivers. A comparison with aquifer waters in the basin showed that even at low flow, rivers rarely consist solely of isotopically unmodified groundwater. Long‐term monitoring (2003–2007) of carbon stable isotopes in dissolved inorganic carbon (DIC) in the Thames revealed a complex interplay between respiration, photosynthesis and evasion, but with a mean interannual δ¹³C‐DIC value of –14.8 ± 0.5‰, exchange with atmospheric carbon could be ruled out. Quarterly monitoring of the tributaries (October 2011–April 2013) indicated that in addition to the aforementioned factors, river flow variations and catchment characteristics were likely to affect δ¹³C‐DIC. Comparison with basin groundwaters of different alkalinity and δ¹³C‐DIC values showed that the origin of river baseflow is usually obscured. The findings show that long‐term monitoring of environmental tracers can help to improve the understanding of how lowland river catchments function. Copyright © NERC 2015. Hydrological Processes © 2015 John Wiley & Sons, Ltd.     

Pre‐ and post‐impoundment nitrogen in the lower Missouri River

Large water‐sample sets collected from 1899 through 1902, 1907, and in the early 1950s allow comparisons of pre‐impoundment and post‐impoundment (1969 through 2008) nitrogen concentrations in the lower Missouri River. Although urban wastes were not large enough to detectably increase annual loads of total nitrogen at the beginning of the 20^th century, carcass waste, stock‐yard manure, and untreated human wastes measurably increased ammonia and organic‐nitrogen concentrations during low flows. Average total‐nitrogen concentrations in both periods were about 2.5 mg/l, but much of the particulate‐organic nitrogen, which was the dominant form of nitrogen around 1900, has been replaced by nitrate. This change in speciation was caused by the nearly 80% decrease in suspended‐sediment concentrations that occurred after impoundment, modern agriculture, drainage of riparian wetlands, and sewage treatment. Nevertheless, bioavailable nitrogen has not been low enough to limit primary production in the Missouri River since the beginning of the 20^th century. Nitrate concentrations have increased more rapidly from 2000 through 2008 (5 to 12% per year), thus increasing bioavailable nitrogen delivered to the Mississippi River and affecting Gulf Coast hypoxia. The increase in nitrate concentrations with distance downstream is much greater during the post‐impoundment period. If strategies to decrease total‐nitrogen loads focus on particulate N, substantial decreases will be difficult because particulate nitrogen is now only 23% of total nitrogen in the Missouri River. A strategy aimed at decreasing particulates also could further exacerbate land loss along the Gulf of Mexico, which has been sediment starved since Missouri River impoundment. In contrast, strategies or benchmarks aimed at decreasing nitrate loads could substantially decrease nitrogen loadings because nitrates now constitute over half of the Missouri's nitrogen input to the Mississippi. Ongoing restoration and creation of wetlands along the Missouri River could be part of such a nitrate‐reduction strategy. Published 2013. This article is a U.S. Government work and is in the public domain in the USA.     

Bound and mobile soil water isotope ratios are affected by soil texture and mineralogy,whereas extraction method influences their measurement

Questions persist about interpreting isotope ratios of bound and mobile soil water pools, particularly relative to clay content and extraction conditions. Interactions between pools and resulting extracted water isotope composition are presumably related to soil texture, yet few studies have manipulated the bound pool to understand its influence on soil water processes. Using a series of drying and spiking experiments, we effectively labelled bound and mobile water pools in soils with varying clay content. Soils were first vacuum dried to remove residual water, which was then replaced with heavy isotope-enriched water prior to oven drying and spiking with heavy isotope-depleted water. Water was extracted via centrifugation or cryogenic vacuum distillation (at four temperatures) and analysed for oxygen and hydrogen isotope ratios via isotope ratio mass spectrometry. Water from centrifuged samples fell along a mixing line between the two added waters but was more enriched in heavy isotopes than the depleted label, demonstrating that despite oven drying, a residual pool remains and mixes with the mobile water. Soils with higher clay + silt content appeared to have a larger bound pool. Water from vacuum distillation samples have a significant temperature effect, with high temperature extractions yielding progressively more heavy isotope-enriched values, suggesting that Rayleigh fractionation occurred at low temperatures in the vacuum line. By distinctly labelling bound and mobile soil water pools, we detected interactions between the two that were dependent on soil texture. Although neither extraction method appeared to completely extract the combined bound and mobile (total water) pool, centrifugation and high temperature cryogenic vacuum distillations were comparable for both δ²H and δ¹⁸O of soil water isotope ratios.     

The potential for agricultural land use change to reduce flood risk in a large watershed

Effects of agricultural land management practices on surface runoff are evident at local scales, but evidence for watershed‐scale impacts is limited. In this study, we used the Soil and Water Assessment Tool model to assess changes in downstream flood risks under different land uses for the large, intensely agricultural, Raccoon River watershed in Iowa. We first developed a baseline model for flood risk based on current land use and typical weather patterns and then simulated the effects of varying levels of increased perennials on the landscape under the same weather patterns. Results suggest that land use changes in the Raccoon River could reduce the likelihood of flood events, decreasing both the number of flood events and the frequency of severe floods. The duration of flood events were not substantially affected by land use change in our assessment. The greatest flood risk reduction was associated with converting all cropland to perennial vegetation, but we found that converting half of the land to perennial vegetation or extended rotations (and leaving the remaining area in cropland) could also have major effects on reducing downstream flooding potential. We discuss the potential costs of adopting the land use change in the watershed to illustrate the scale of subsidies required to induce large‐scale conversion to perennially based systems needed for flood risk reduction. Copyright © 2013 John Wiley & Sons, Ltd.