Hydrological trends and the evolution of catchment research in the Alptal valley,central Switzerland

When the observation of small headwater catchments in the pre-Alpine Alptal valley (central Switzerland) started in the late 1960s, the researchers were mainly interested in questions related to floods and forest management. Investigations of geomorphological processes in the steep torrent channels followed in the 1980s, along with detailed observations of biogeochemical and ecohydrological processes in individual forest stands. More recently, research in the Alptal has addressed the impacts of climate change on water supply and runoff generation. In this article, we describe, for the first time, the evolution of catchment research at Alptal, and present new analyses of long-term trends and short-term hydrologic behaviour. Hydrometeorological time series from the past 50 years show substantial interannual variability, but only minimal long-term trends, except for the ~2°C increase in mean annual air temperature over the 50-year period, and a corresponding shift towards earlier snowmelt. Similar to previous studies in larger Alpine catchments, the decadal variations in mean annual runoff in Alptal's small research catchments reflect the long-term variability in annual precipitation. In the Alptal valley, the most evident hydrological trends were observed in late spring and are related to the substantial change in the duration of the snow cover. Streamflow and water quality are highly variable within and between hydrological events, suggesting rapid shifts in flow pathways and mixing, as well as changing connectivity of runoff-generating areas. This overview illustrates how catchment research in the Alptal has evolved in response to changing societal concerns and emerging scientific questions.     

Holocene permafrost history and cryostratigraphy in the High‐Arctic Adventdalen Valley,central Svalbard

This paper presents the history and cryostratigraphy of the upper permafrost in the High‐Arctic Adventdalen Valley, central Svalbard. Nineteen frozen sediment cores, up to 10.7 m long, obtained at five periglacial landforms, were analysed for cryostructures, ice, carbon and solute contents, and grain‐size distribution, and were ¹⁴C‐ and OSL‐dated. Spatial variability in ice and carbon contents is closely related to the sedimentary history and mode of permafrost aggradation. In the valley bottom, saline epigenetic permafrost with pore ice down to depths of 10.7 m depth formed in deltaic sediments since the mid‐Holocene; cryopegs were encountered below 6 m. In the top 1 to 5 m, syngenetic and quasi‐syngenetic permafrost with microlenticular, lenticular, suspended and organic‐matrix cryostructures developed due to loess and alluvial sedimentation since the colder late Holocene, which resulted in the burial of organic material. At the transition between deltaic sediments and loess, massive ice bodies occurred. A pingo developed where the deltaic sediments reached the surface. On hillslopes, suspended cryostructure on solifluction sheets indicates quasi‐syngenetic permafrost aggradation; lobes, in contrast, were ice‐poor. Suspended cryostructure in eluvial deposits reflects epigenetic or quasi‐syngenetic permafrost formation on a weathered bedrock plateau. Landform‐scale spatial variations in ground ice and carbon relate to variations in slope, sedimentation rate, moisture conditions and stratigraphy. Although the study reveals close links between Holocene landscape evolution and permafrost history, our results emphasize a large uncertainty in using terrain surface indicators to infer ground‐ice contents and upscale from core to landform scale in mountainous permafrost landscapes.     

Three-dimensional finite element simulation of fracture propagation in rock specimens with pre-existing fissure(s) under compression and their strength analysis

A FORTRAN program, consistent with the commercially available finite element (FE) code ABAQUS, is developed based on a three-dimensional (3D) linear elastic brittle damage constitutive model with two damage criteria. To consider the heterogeneity of rock, the developed FORTRAN program is used to set the stiffness and strength properties of each element of the FE model following a Weibull distribution function. The reliability of the program is assessed against available experimental results for granite cylindrical specimens with a throughgoing, flat and inclined fissure. The calibration procedure of the material parameters is explained in detail, and it is shown that the compressive to tensile strength ratio can have a substantial influence on the failure response of the specimens. Numerical simulations are conducted for models with different levels of heterogeneity. The results show a smaller load bearing capacity for models with less homogeneity, representing gradual coalescence of fully damaged elements forming throughout the models during loading. The maximum load bearing capacity is studied for various combinations of inclination angles of two centrally aligned, throughgoing and flat fissures of equal length embedded in cylindrical models under uniaxial and multiaxial loading conditions. The key role of the compressive to tensile strength ratio is highlighted by repeating certain simulations with a lower compressive to tensile strength ratio. It is proven that the peak loads of the rock models with sufficiently small compressive to tensile strength ratios containing two throughgoing fissures of equal length are similar, provided that the minimum inclination angles of the models are the same. The results are presented and discussed with respect to the existing experimental findings in the literature, suggesting that the numerical model applied in this study can provide useful insight into the failure behaviour of rock-like materials.     

Time–space fractional governing equations of transient groundwater flow in confined aquifers: Numerical investigation

In order to model non‐Fickian transport behaviour in groundwater aquifers, various forms of the time–space fractional advection–dispersion equation have been developed and used by several researchers in the last decade. The solute transport in groundwater aquifers in fractional time–space takes place by means of an underlying groundwater flow field. However, the governing equations for such groundwater flow in fractional time–space are yet to be developed in a comprehensive framework. In this study, a finite difference numerical scheme based on Caputo fractional derivative is proposed to investigate the properties of a newly developed time–space fractional governing equations of transient groundwater flow in confined aquifers in terms of the time–space fractional mass conservation equation and the time–space fractional water flux equation. Here, we apply these time–space fractional governing equations numerically to transient groundwater flow in a confined aquifer for different boundary conditions to explore their behaviour in modelling groundwater flow in fractional time–space. The numerical results demonstrate that the proposed time–space fractional governing equation for groundwater flow in confined aquifers may provide a new perspective on modelling groundwater flow and on interpreting the dynamics of groundwater level fluctuations. Additionally, the numerical results may imply that the newly derived fractional groundwater governing equation may help explain the observed heavy‐tailed solute transport behaviour in groundwater flow by incorporating nonlocal or long‐range dependence of the underlying groundwater flow field.     

Predicting quantiles of water quality from catchment characteristics

Water quality is often highly variable both in space and time, which poses challenges for modelling the more extreme concentrations. This study developed an alternative approach to predicting water quality quantiles at individual locations. We focused on river water quality data that were collected over 25 years, at 102 catchments across the State of Victoria, Australia. We analysed and modelled spatial patterns of the 10th, 25th, 50th, 75th and 90th percentiles of the concentrations of sediments, nutrients and salt, with six common constituents: total suspended solids (TSS), total phosphorus (TP), filterable reactive phosphorus (FRP), total Kjeldahl nitrogen (TKN), nitrate-nitrite (NO_x), and electrical conductivity (EC). To predict the spatial variation of each quantile for each constituent, we developed statistical regression models and exhaustively searched through 50 catchment characteristics to identify the best set of predictors for that quantile. The models predict the spatial variation in individual quantiles of TSS, TKN and EC well (66%–96% spatial variation explained), while those for TP, FRP and NO_x have lower performance (37%–73% spatial variation explained). The most common factors that influence the spatial variations of the different constituents and quantiles are: annual temperature, percentage of cropping land area in catchment and channel slope. The statistical models developed can be used to predict how low- and high-concentration quantiles change with landscape characteristics, and thus provide a useful tool for catchment managers to inform planning and policy making with changing climate and land use conditions.     

Dynamic modelling of retrogressive landslides with emphasis on the role of clay sensitivity

This paper presents a detailed numerical study of the retrogressive failure of landslides in sensitive clays. The dynamic modelling of the landslides is carried out using a novel continuum approach, the particle finite element method, complemented with an elastoviscoplastic constitutive model. The multiwedge failure mode in the collapse is captured successfully, and the multiple retrogressive failures that have been widely observed in landslides in sensitive clays are reproduced with the failure mechanism, the kinematics, and the deposition being discussed in detail. Special attention has been paid to the role of the clay sensitivity on each retrogressive failure as well as on the final retrogression distance and the final run‐out distance via parametric studies. Moreover, the effects of the viscosity of sensitive clays on the failure are also investigated for different clay sensitivities.     

Petrochronology of oxidized granulites from southern Peru

Sapphirine–quartz granulites from the Cocachacra region of the Arequipa Massif in southern Peru record early Neoproterozoic ultrahigh‐temperature metamorphism. Phase equilibrium modelling and zircon petrochronology are used to quantify timing and pressure–temperature (P–T) conditions of metamorphism. Modelling of three magnetite‐bearing sapphirine–quartz samples indicates peak temperatures of >950°C at ~0.7 GPa and a clockwise P–T evolution. Elevated concentrations of Al in orthopyroxene are also consistent with ultrahigh‐temperature conditions. Neoblastic zircon records ages of c. 1.0–0.9 Ga that are interpreted to record protracted ultrahigh‐temperature metamorphism. Th/U ratios of zircon of up to 100 reflect U‐depleted whole‐rock compositions. Concentrations of heavy rare earth elements in zircon do not show systematic trends with U–Pb age but do correlate with variable whole‐rock compositions. Very large positive Ce anomalies in zircon from two samples probably relate to strongly oxidizing conditions during neoblastic zircon crystallization. Low concentrations of Ti‐in‐zircon (<10 ppm) are interpreted to result from reduced titania activities due to the strongly oxidized nature of the granulites and the sequestration of titanium‐rich minerals away from the reaction volume. Whole‐rock compositions and oxidation state have a strong influence on the trace element composition of metamorphic zircon, which has implications for interpreting the geological significance of ages retrieved from zircon in oxidized metamorphic rocks.     

Active capture and anomalous adsorption: New mechanisms for the incorporation of heavy noble gases

Abstract— Active capture is a new process for the incorporation of large quantities of heavy noble gases into growing surfaces. Adsorption in the conventional sense involves surface bonding by polarization (Van der Waals forces). What is referred to as “anomalous adsorption” of heavy noble gases involves chemical bonds and can occur when other (more chemically active) species are not available to preempt sites with unfilled bonds. Anomalous adsorption has been observed under conditions of fracture, vacuum deposition and ionizing radiation. Active capture depends upon anomalous adsorption to retain noble gases on a surface long enough to be captured in a growing surface film as it is deposited. The fundamental principle may be the impingement onto the growing film with sufficient energy to liberate surface electrons (work function energy of a few electronvolts) so that they are retained by anomalous adsorption long enough to be entrapped in the growing surface. Trapping efficiencies of ?1% have been observed for Kr and Xe in laboratory experiments, implying a fundamentally new mechanism for the incorporation of heavy noble gases onto surfaces. It may play a role in explaining the large concentrations of planetary noble gases contained in phase‐Q.