A smoothed particle hydrodynamics modelling of soil–water mixing and resulting changes in average strength

Soil–water interaction is a pivotal process in many underwater geohazards such as underwater landslides where soil sediments gradually evolve into turbidity currents after interactions with ambient water. Due to the large deformations, multiphase interactions and phase changes this involves, investigations from numerical modelling of the transition process have been limited so far. This study explores a simple numerical replication of such soil–water mixing with respect to changes in average strength using smoothed particle hydrodynamics (SPH). A uniform viscoplastic model is used for both the solid-like and fluid-like SPH particles. The proposed numerical solution scheme is verified by single-phase dam break tests and multiphase simple shear tests. SPH combinations of solid-like and fluid-like particles can replicate the clay–water mixture as long as the liquidity index of the solid-like particles is larger than unity. The proposed numerical scheme is shown to capture key features of an underwater landslide such as hydroplaning, water entrainment and wave generation and thus shows promise as a tool to simulate the whole process of subaquatic geohazards involving solid–fluid transition during mass transport.     

Analytical prediction of aggregates' effects on the ITZ volume fraction and Young's modulus of concrete

In order to evaluate analytically the ITZ volume fraction (f_ITZ) in concrete, a three phase model is proposed for the random concrete microstructure using the Voronoï tessellation. Within this model, the ITZ local thickness is a statistical variable depending on the local paste thickness available between each couple of neighbouring aggregates. The f_ITZ is found to not exceed 7% for typical concretes. Then, the concrete Young's modulus is predicted analytically using a four‐phase generalized self consistent model but in which the proposed f_ITZ is considered. It is found that the concrete Young's modulus increases when increasing aggregates volume fraction, aggregates maximum size and the proportion of coarse aggregates and when decreasing the ITZ thickness and Young's modulus. Finally, the validity of the proposed model is discussed based on a comparison between its predictions and three sets of experimental results related to normal and high strength concretes taken from literature. Copyright © 2016 John Wiley & Sons, Ltd.     

Experimental determination of REE partition coefficients between zircon,garnet and melt: a key to understanding high‐T crustal processes

The partitioning of rare earth elements (REE) between zircon, garnet and silicate melt was determined using synthetic compositions designed to represent partial melts formed in the lower crust during anatexis. The experiments, performed using internally heated gas pressure vessels at 7 kbar and 900–1000 °C, represent equilibrium partitioning of the middle to heavy REE between zircon and garnet during high‐grade metamorphism in the mid to lower crust. The D_REE (zircon/garnet) values show a clear partitioning signature close to unity from Gd to Lu. Because the light REE have low concentrations in both minerals, values are calculated from strain modelling of the middle to heavy REE experimental data; these results show that zircon is favoured over garnet by up to two orders of magnitude. The resulting general concave‐up shape to the partitioning pattern across the REE reflects the preferential incorporation of middle REE into garnet, with D_Gd (zircon/garnet) ranging from 0.7 to 1.1, D_Ho (zircon/garnet) from 0.4 to 0.7 and D_Lu (zircon/garnet) from 0.6 to 1.3. There is no significant temperature dependence in the zircon–garnet REE partitioning at 7 kbar and 900–1000 °C, suggesting that these values can be applied to the interpretation of zircon–garnet equilibrium and timing relationships in the ultrahigh‐T metamorphism of low‐Ca pelitic and aluminous granulites.     

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.     

Snow model sensitivity analysis to understand spatial and temporal snow dynamics in a high‐elevation catchment

In this paper, we addressed a sensitivity analysis of the snow module of the GEOtop2.0 model at point and catchment scale in a small high‐elevation catchment in the Eastern Italian Alps (catchment size: 61 km²). Simulated snow depth and snow water equivalent at the point scale were compared with measured data at four locations from 2009 to 2013. At the catchment scale, simulated snow‐covered area (SCA) was compared with binary snow cover maps derived from moderate‐resolution imaging spectroradiometer (MODIS) and Landsat satellite imagery. Sensitivity analyses were used to assess the effect of different model parameterizations on model performance at both scales and the effect of different thresholds of simulated snow depth on the agreement with MODIS data. Our results at point scale indicated that modifying only the “snow correction factor” resulted in substantial improvements of the snow model and effectively compensated inaccurate winter precipitation by enhancing snow accumulation. SCA inaccuracies at catchment scale during accumulation and melt period were affected little by different snow depth thresholds when using calibrated winter precipitation from point scale. However, inaccuracies were strongly controlled by topographic characteristics and model parameterizations driving snow albedo (“snow ageing coefficient” and “extinction of snow albedo”) during accumulation and melt period. Although highest accuracies (overall accuracy = 1 in 86% of the catchment area) were observed during winter, lower accuracies (overall accuracy < 0.7) occurred during the early accumulation and melt period (in 29% and 23%, respectively), mostly present in areas with grassland and forest, slopes of 20–40°, areas exposed NW or areas with a topographic roughness index of ?0.25 to 0 m. These findings may give recommendations for defining more effective model parameterization strategies and guide future work, in which simulated and MODIS SCA may be combined to generate improved products for SCA monitoring in Alpine catchments.     

Dynamics of Saxothuringian subduction channel/wedge constrained by phase‐equilibria modelling and micro‐fabric analysis

Subduction and exhumation dynamics can be investigated through analysis of metamorphic and deformational evolution of associated high‐grade rocks. The Erzgebirge anticline, which forms at the boundary between the Saxothuringian and Teplá‐Barrandian domains of the Bohemian Massif, provides a useful study area for these processes owing to the occurrence of numerous meta‐basites preserving eclogite facies assemblages, and coesite and diamond bearing quartzofeldspathic lithologies indicating subduction to deep mantle depths. The prograde and retrograde evolution of meta‐basite from the Czech portion of the Erzgebirge anticline has been constrained through a combination of thermodynamic modelling and conventional thermobarometry. Garnet growth zoning indicates that the rocks underwent burial and heating to peak conditions of 2.6 GPa and at least 615 °C. Initial exhumation occurred with concurrent cooling and decompression resulting in the growth of amphibole and zoisite poikiloblasts overgrowing and including the eclogite facies assemblage. The development of clinopyroxene–plagioclase–amphibole symplectites after omphacite and Al‐rich rims on matrix amphibole indicate later heating at the base of the lower crust. Omphacite microstructures, in particular grain size analysis and lattice‐preferred orientations, indicate that the prograde evolution was characterized by a constrictional strain geometry transitioning into plane strain and oblate fabrics during exhumation. The initial constrictional strain pattern is interpreted as being controlled by competing slab pull and crustal buoyancy forces leading to necking of the subducting slab. The transition to plane strain and flattening geometries represents transfer of material from the subducting lithosphere into a subduction channel, break‐off of the dense slab and rebound of the buoyant crustal material.     

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.     

Modelling and experimental study of coupled exchange of p,p′‐DDE and colloids in a flume simulated stream–streambed system

Stream–subsurface exchange plays a significant role in the fate and transport of contaminants in streams. It has been modelled explicitly by considering fundamental processes such as hydraulic exchange, colloid filtration, and contaminant interactions with streambed sediments and colloids. The models have been successfully applied to simulate the transport of inorganic metals and nutrients. In this study, laboratory experiments were conducted in a recirculating flume to investigate the exchange of a hydrophobic organic contaminant, p,p′‐dichloro‐diphenyl‐dichloroethane (DDE), between a stream and a quartz sand bed. A previously developed process‐based multiphase exchange model was modified by accounting for the p,p′‐DDE kinetic adsorption to and desorption from the bed sediments/colloids and was applied to interpret the experimental results. Model input parameters were obtained by conducting independent small‐scale batch experiments. Results indicate that the immobilization of p,p′‐DDE in the quartz sand bed can occur under representative natural stream conditions. The observed p,p′‐DDE exchange was successfully simulated by the process‐based model. The model sensitivity analysis results show that the exchange of p,p′‐DDE can be sensitive to either the sediment sorption/desorption parameters or colloidal parameters depending on the experimental conditions tested. For the experimental conditions employed here, the effect of colloids on contaminant transport is expected to be minimal, and the stream–subsurface exchange of p,p′‐DDE is dominated by the interaction of p,p′‐DDE with bed sediment. The work presented here contributes to a better mechanistic understanding of the complex transport process that hydrophobic organic contaminants undergo in natural streams and to the development of reliable, predictive models for the assessment of impacted streams. Copyright © 2015 John Wiley & Sons, Ltd.     

Hammond Hill Research Catchment: Supporting hydrologic investigations of rooting zone and vegetation water dynamics under climate change

The Hammond Hill Research Catchment (HH) is a small (120 ha), temperate, second order tributary to Six Mile Creek, Cayuga Lake, and the Great Lakes (42.42°, −76.32°). The HH has been monitored since January 2017 for the purpose of understanding how recent infiltration mixes with antecedent soil water on hillslope forest floors and the spatial and temporal patterns of Root Water Uptake (RWU) by temperate northeastern US tree species (eastern hemlock [Tsuga canadensis], American beech [Fagus grandifolia], and sugar maple [Acer saccharum]). These data are informing us about the hydrologic consequences of anticipated tree species composition change and supporting the development of more refined ecohydrological models. The glaciated catchment is underlain by a shallow confining siltstone layer (1–1.5 m depth) and densely covered with an approximately 60 year old regrowth mixed species forest of hemlock, beech, and other deciduous tree species common to the northeastern US. Current datasets from the HH include precipitation snow water equivalent, discharge, and associated isotopic water compositions, δ²H & δ¹⁸O. Measurements of (top 10 cm) soil water content, as well as bulk soil water and hemlock and beech xylem isotopic compositions are made at several locations across a topographic wetness gradient. The near-term role of the HH is to support an understanding of the environmental and ecological drivers of plant RWU competition. All data from the HH are publicly available.