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.     

Experimental partitioning of trace elements into schreibersite with applications to IIG iron meteorites

Some of the defining characteristics of the IIG iron meteorite group are their high bulk P contents and massive, coarse schreibersite, which have been calculated to make up roughly 11–14 wt% of each specimen. In this study, we produced two data sets to investigate the formation of schreibersites in IIG irons: measurements of trace elements in the IIG iron meteorite Twannberg and experimental determinations of trace element partitioning into schreibersite. The schreibersite‐bearing experiments were conducted with schreibersite in equilibrium with a P‐rich melt and with bulk Ni contents ranging from 0 to 40 wt%. The partitioning behavior for the 20 elements measured in this study did not vary with Ni content. Comparison of the Twannberg measurements with the experimental results required a correction factor to account for the fact that the experiments were conducted in a simplified system that did not contain a solid metal phase. Previously determined solid metal/P‐rich melt partition coefficients were applied to infer schreibersite/solid metal partitioning behavior from the experiments, and once this correction was applied, the two data sets showed broad similarities between the schreibersite/solid metal distribution of elements. However, there were also differences noted, in particular between the Ni and P contents of the solid metal relative to the schreibersite inferred from the experiments compared to that measured in the Twannberg sample. These differences support previous interpretations that subsolidus schreibersite evolution has strongly influenced the Ni and P content now present in the solid metal phase of IIG irons. Quantitative attempts to match the IIG solid metal composition to that of late‐stage IIAB irons through subsolidus schreibersite growth were not successful, but qualitatively, this study corroborates the striking similarities between the IIAB and IIG groups, which are highly suggestive of a possible genetic link between the groups as has been previously proposed.     

The <Emphasis Type="Italic">Wind</Emphasis>/EPACT Proton Event Catalog (1996 – 2016)

We present the finalized catalog of solar energetic proton events detected by the Wind/EPACT instrument over the period 1996?–?2016. Onset times, peak times, peak proton intensity and onset-to-peak proton fluence are evaluated for the two available energy channels, at about 25 and 50 MeV. We describe the procedure utilized to identify the proton events and to relate them to their solar origin (in terms of flares and coronal mass ejections). The statistical relationships between the energetic protons and their origin (linear and partial correlation analysis) are reported and discussed in view of earlier findings. Finally, the different trends found in the first 8 years of Solar Cycles 23 and 24 are discussed.     

Simulating the Propagation of Infrasonic Waves and Estimating the Energy of the Chelyabinsk Meteoroid Explosion Observed on February 15, 2013

Results obtained from simulating the propagation of infrasonic waves from the Chelyabinsk meteoroid explosion observed on February 15, 2013, are given. The pseudodifferential parabolic equation (PDPE) method has been used for calculations. Data on infrasonic waves recorded at the IS31 station (Aktyubinsk, Kazakhstan), located 542.7 km from the likely location of the explosion, have been analyzed. Six infrasonic arrivals (isolated clearly defined pulse signals) were recorded. It is shown that the first “fast” arrival (F) corresponds to the propagation of infrasound in a surface acoustic waveguide. The rest of the arrivals (T1–T5) are thermospheric. The agreement between the results of calculations based on the PDPE method and experimental data is satisfactory. The energy E of the explosion has been estimated using two methods. One of these methods is based on the law of conservation of the acoustic pulse I, which is a product of the wave profile area S/2 of the signal under analysis and the distance to its source E _I [kt] = 1.38 × 10^–10 (I [kg/s])^1.482. The other method is based on the relation between the energy of explosion and the dominant period T of recorded signal E _T [kt] = 1.02 × (T [s]²/σ)^3/2, where σ is the dimensionless distance determining the degree of nonlinear effects during the propagation of sound along ray trajectories. According to the data, the explosion energy E _I,T ranges from 1.87 to 32 kt TNT.     

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.     

Determination of Supermassive Black Hole Spins Based on the Standard Shakura–Sunyaev Accretion Disk Model and Polarimetric Observations

Based on spectropolarimetry for 47 type 1 active galactic nuclei observed with the 6-m BTA telescope, we have estimated the spins of the supermassive black holes at the centers of these galaxies. We have determined the spins based on the standard Shakura–Sunyaev accretion disk model. About 70% of the investigated active galactic nuclei are shown to have Kerr supermassive black holes with a dimensionless spin greater than 0.9.     

Critical Latitude in Tidal Dynamics Using the Kara Sea as an Example

It is well known that, within the linear nonviscous equations of tidal dynamics, the amplitudes of oscillations of the barotropic and baroclinic tidal velocity components unlimitedly increase when approaching the critical latitude. It is also known that the linear equations of tidal dynamics with a constant and specified vertical eddy viscosity indicate the occurrence of significant tidal velocity shears in the near-bottom layer, which are responsible for increasing the baroclinic tidal energy dissipation, the turbulent kinetic energy, and the thickness of the bottom boundary layer. The first circumstance—the growth of the amplitudes of oscillations of the barotropic and baroclinic tidal velocity components—is due to the elimination in the original equations of small terms, which are small everywhere except for the critical latitude zone. The second circumstance—the occurrence of significant tidal velocity shears—is due to the fact that internal tidal waves, which induce the dissipation of the baroclinic tidal energy and the diapycnal diffusion, are either not taken into account or described inadequately. It is suggested that diapycnal diffusion can lead to the degeneration (complete or partial) of tidal velocity shears, with all the ensuing consequences. The aforesaid is confirmed by simulation results obtained using the QUODDY-4 high-resolution three-dimensional finite-element hydrostatic model along the 66.25° E section, which passes in the Kara Sea across the critical latitude.     

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.