The influence of coupled physical swelling and chemical reactions on deformable geomaterials

Coupled thermo‐hydro‐mechanical‐chemical modelling has attracted attention in past decades due to many contemporary geotechnical engineering applications (e.g., waste disposal, carbon capture and storage). However, molecular‐scale interactions within geomaterials (e.g., swelling and dissolution/precipitation) have a significant influence on the mechanical behaviour, yet are rarely incorporated into existing Thermal‐Hydro‐Mechanical‐Chemical (THMC) frameworks. This paper presents a new coupled hydro‐mechanical‐chemical constitutive model to bridge molecular‐scale interactions with macro‐physical deformation by combining the swelling and dissolution/precipitation through an extension of the new mixture‐coupling theory. Entropy analysis of the geomaterial system provides dissipation energy, and Helmholtz free energy gives the relationship between solids and fluids. Numerical simulation is used to compare with the selected recognized models, which demonstrates that the swelling and dissolution/precipitation processes may have a significant influence on the mechanical deformation of the geomaterials.     

Including spatial distribution in a data‐driven rainfall‐runoff model to improve reservoir inflow forecasting in Taiwan

Multi‐step ahead inflow forecasting has a critical role to play in reservoir operation and management in Taiwan during typhoons as statutory legislation requires a minimum of 3‐h warning to be issued before any reservoir releases are made. However, the complex spatial and temporal heterogeneity of typhoon rainfall, coupled with a remote and mountainous physiographic context, makes the development of real‐time rainfall‐runoff models that can accurately predict reservoir inflow several hours ahead of time challenging. Consequently, there is an urgent, operational requirement for models that can enhance reservoir inflow prediction at forecast horizons of more than 3 h. In this paper, we develop a novel semi‐distributed, data‐driven, rainfall‐runoff model for the Shihmen catchment, north Taiwan. A suite of Adaptive Network‐based Fuzzy Inference System solutions is created using various combinations of autoregressive, spatially lumped radar and point‐based rain gauge predictors. Different levels of spatially aggregated radar‐derived rainfall data are used to generate 4, 8 and 12 sub‐catchment input drivers. In general, the semi‐distributed radar rainfall models outperform their less complex counterparts in predictions of reservoir inflow at lead times greater than 3 h. Performance is found to be optimal when spatial aggregation is restricted to four sub‐catchments, with up to 30% improvements in the performance over lumped and point‐based models being evident at 5‐h lead times. The potential benefits of applying semi‐distributed, data‐driven models in reservoir inflow modelling specifically, and hydrological modelling more generally, are thus demonstrated. Copyright © 2012 John Wiley & Sons, Ltd.     

What is the source of baseflow in agriculturally fragmented catchments? Complex groundwater/surface‐water interactions in three tributary catchments of the Wabash River,Indiana, USA

Some conceptual models suggest that baseflow in agriculturally fragmented watersheds may contain little, if any, groundwater. This has critical implications for stream quality and ecosystem functioning. Here, we (a) identify the sources and flowpaths contributing to baseflow using ²²²Rn and ⁸⁷Sr/⁸⁶Sr and (b) quantify mean apparent ages of groundwater and baseflow using multiple isotopic tracers (CFC, SF₆, ³⁶Cl, and ³H) in 4 small (0.08 to 0.64 km²) tributary catchments to the Wabash River in Indiana, USA. ²²²Rn activities and ⁸⁷Sr/⁸⁶Sr ratios indicate that baseflow in 3 catchments is sourced primarily from groundwater; baseflow in the fourth is dominated by a source similar to agricultural run‐off. CFC‐12 data indicate that springs in 1 catchment are discharging significant proportions of water that recharged between 1974 (42 ± 2 years) and 1961 (55 ± 2 years). Those same springs have ³⁶Cl/Cl ratios between 1,381.08 ± 29.37 (×10^?15) and 1,530.64 ± 27.65 (×10^?15) indicating that a substantial proportion of the discharge likely recharged between 1975 (41 years) and 1950 (66 years). Groundwater samples collected from streambed mini‐piezometers in a separate catchment have CFC‐12 concentrations indicating that a large proportion of the recharge occurred between 1948 (68 ± 2 years) and 1950 (66 ± 2 years). Repeat sampling conducted in September 2015 after above‐average summer rainfall did not show significant decreases in mean apparent age. The relatively old ages observed in 3 of the catchments can be explained by geological complexities that are likely present in all 4 catchments, but overwhelmed by flow from the shallow phreatic aquifer in the fourth catchment.     

Pressure and temperature conditions for crystallization of metamorphic allanite and monazite in metapelites: a case study from the Miyar Valley (high Himalayan Crystalline of Zanskar,NW India)

The textural and chemical evolution of allanite and monazite along a well‐constrained prograde metamorphic suite in the High Himalayan Crystalline of Zanskar was investigated to determine the P–T conditions for the crystallization of these two REE accessory phases. The results of this study reveals that: (i) allanite is the stable REE accessory phase in the biotite and garnet zone and (ii) allanite disappears at the staurolite‐in isograd, simultaneously with the occurrence of the first metamorphic monazite. Both monazite and allanite occur as inclusions in staurolite, indicating that the breakdown of allanite and the formation of monazite proceeded during staurolite crystallization. Staurolite growth modelling indicates that staurolite crystallized between 580 and 610 °C, thus setting the lower temperature limit for the monazite‐forming reaction at ~600 °C. Preservation of allanite and monazite inclusions in garnet (core and rim) constrains the garnet molar composition when the first monazite was overgrown and subsequently encompassed by the garnet crystallization front. Garnet growth modelling and the intersection of isopleths reveal that the monazite closest to the garnet core was overgrown by the garnet advancing crystallization front at 590 °C, which establishes an upper temperature limit for monazite crystallization. Significantly, the substitution of allanite by monazite occurs in close spatial proximity, i.e. at similar P–T conditions, in all rock types investigated, from Al‐rich metapelites to more psammitic metasedimentary rocks. This indicates that major silicate phases, such as staurolite and garnet, do not play a significant role in the monazite‐forming reaction. Our data show that the occurrence of the first metamorphic monazite in these rocks was mainly determined by the P–T conditions, not by bulk chemical composition. In Barrovian terranes, dating prograde monazite in metapelites thus means constraining the time when these rocks reached the 600 °C isotherm.     

Severe soil erosion during a 3‐day exceptional rainfall event: combining modelling and field data for a fallow cereal field

Exceptional rainfall events cause significant losses of soil, although few studies have addressed the validation of model predictions at field scale during severe erosive episodes. In this study, we evaluate the predictive ability of the enhanced Soil Erosion and Redistribution Tool (SERT‐2014) model for mapping and quantifying soil erosion during the exceptional rainfall event (~235 mm) that affected the Central Spanish Pyrenees in October 2012. The capacity of the simulation model is evaluated in a fallow cereal field (1.9 ha) at a high spatial scale (1 × 1 m). Validation was performed with field‐quantified rates of soil loss in the rills and ephemeral gullies and also with a detailed map of soil redistribution. The SERT‐2014 model was run for the six rainfall sub‐events that made up the exceptional event, simulating the different hydrological responses of soils with maximum runoff depths ranging between 40 and 1017 mm. Predicted average and maximum soil erosion was 11 and 117 Mg ha^?1 event^?1, respectively. Total soil loss and sediment yield to the La Reina gully amounted to 16.3 and 9.0 Mg event^?1. These rates are in agreement with field estimations of soil loss of 20.0 Mg event^?1. Most soil loss (86%) occurred during the first sub‐event. Although soil accumulation was overestimated in the first sub‐event because of the large amount of detached soil, the enhanced SERT‐2014 model successfully predicted the different spatial patterns and values of soil redistribution for each sub‐event. Further research should focus on stream transport capacity. Copyright © 2014 John Wiley & Sons, Ltd.     

The stiff upper LIP: investigating the High Arctic Large Igneous Province

The Canadian Arctic Islands expose a complex network of dykes and sills that belong to the High Arctic Large Igneous Province (HALIP), which intruded volatile‐rich sedimentary rocks of the Sverdrup Basin (shale, limestone, sandstone and evaporite) some 130 to 120 million years ago. There is thus great potential in studying the HALIP to learn how volatile‐rich sedimentary rocks respond to magmatic heating events during LIP emplacement. The HALIP remains, however, one of the least well known LIPs on the planet due to its remote location, short field season, and harsh climate. A Canadian–Swedish team of geologists set out in summer 2015 to further explore HALIP sills and their sedimentary host rocks, including the sampling of igneous and meta‐sedimentary rocks for subsequent geochemical analysis, and high pressure‐temperature petrological experiments to help define the actual processes and time‐scales of magma–sediment interaction. The research results will advance our understanding of how climate‐active volatiles such as CO_2, SO₂ and CH₄ are mobilised during the magma–sediment interaction related to LIP events, a process which is hypothesised to have drastically affected Earth's carbon and sulphur cycles. In addition, assimilation of sulphate evaporites, for example, is anticipated to trigger sulphide immiscibility in the magma bodies and in so doing could promote the formation of Ni‐PGE ore bodies. Here we document the joys and challenges of ‘frontier arctic fieldwork’ and discuss some of our initial observations from the High Arctic Large Igneous Province.     

Middle Paleolithic occupation on a Marine Isotope Stage 5 lakeshore in the Nefud Desert,Saudi Arabia

Major hydrological variations associated with glacial and interglacial climates in North Africa and the Levant have been related to Middle Paleolithic occupations and dispersals, but suitable archaeological sites to explore such relationships are rare on the Arabian Peninsula. Here we report the discovery of Middle Paleolithic assemblages in the Nefud Desert of northern Arabia associated with stratified deposits dated to 75,000 years ago. The site is located in close proximity to a substantial relict lake and indicates that Middle Paleolithic hominins penetrated deeply into the Arabian Peninsula to inhabit landscapes vegetated by grasses and some trees. Our discovery supports the hypothesis of range expansion by Middle Paleolithic populations into Arabia during the final humid phase of Marine Isotope Stage 5, when environmental conditions were still favorable.