Numerical Modelling of the Anisotropic Mechanical Behaviour of Opalinus Clay at the Laboratory-Scale Using FEM/DEM

The Opalinus Clay (OPA) is an argillaceous rock formation selected to host a deep geologic repository for high-level nuclear waste in Switzerland. It has been shown that the excavation damaged zone (EDZ) in this formation is heavily affected by the anisotropic mechanical response of the material related to the presence of bedding planes. In this context, the purpose of this study is twofold: (i) to illustrate the new developments that have been introduced into the combined finite-discrete element method (FEM/DEM) to model layered materials and (ii) to demonstrate the effectiveness of this new modelling approach in simulating the short-term mechanical response of OPA at the laboratory-scale. A transversely isotropic elastic constitutive law is implemented to account for the anisotropic elastic modulus, while a procedure to incorporate a distribution of preferentially oriented defects is devised to capture the anisotropic strength. Laboratory results of indirect tensile tests and uniaxial compression tests are used to calibrate the numerical model. Emergent strength and deformation properties, together with the simulated damage mechanisms, are shown to be in strong agreement with experimental observations. Subsequently, the calibrated model is validated by investigating the effect of confinement and the influence of the loading angle with respect to the specimen anisotropy. Simulated fracture patterns are discussed in the context of the theory of brittle rock failure and analyzed with reference to the EDZ formation mechanisms observed at the Mont Terri Underground Research Laboratory.     

Micro-seismic monitoring after the shipwreck of the Costa Concordia at Giglio Island (Italy)

A micro-seismic network was used for monitoring the wreck of the Costa Concordia cruise ship, wrecked and run agrounded along the Giglio Island coasts during the night of 13 January 2012, until its removal. The seismic traces were processed by means of real-time and “a posteriori” procedures to detect transients that could be ascribed to wreck movements on the sea bed to integrate this information in an early warning system for assessing the wreck stability. After a first discrimination of the transients using amplitude criteria we proceeded to the localization of the detected signals to focus the attention only on the transients originated in the shipwreck resting area. The results showed that most of the events localized on the wreck were likely related to human work activities or sudden internal brittle failure but not to displacements on the seafloor. Instead, the displacements are associated to the impact on the vessel of great sea storms which approach were well correlated with the increasing seismic noise at low frequency. The carried out procedures based on this unique dataset represent an opportunity to test seismic monitoring techniques also in not usual engineering context to support emergency management activities.     

The distribution of supermassive black holes in the nuclei of nearby galaxies

The growth of supermassive black holes by merging and accretion in hierarchical models of galaxy formation is studied by means of Monte Carlo simulations. A tight linear relation between masses of black holes and masses of bulges arises if the mass accreted by supermassive black holes scales linearly with the mass-forming stars and if the redshift evolution of mass accretion tracks closely that of star formation. Differences in redshift evolution between black hole accretion and star formation introduce a considerable scatter in this relation. A non-linear relation between black hole accretion and star formation results in a non-linear relation between masses of remnant black holes and masses of bulges. The relation of black hole mass to bulge luminosity observed in nearby galaxies and its scatter are reproduced reasonably well by models in which black hole accretion and star formation are linearly related but do not track each other in redshift. This suggests that a common mechanism determines the efficiency for black hole accretion and the efficiency for star formation, especially for bright bulges.     

Preparation and Measurement of Cassiterite for Sn Isotope Analysis

Increased interest in the fractionation of Sn isotopes has led to the development of several techniques for preparing cassiterite (SnO₂, the primary ore of Sn) for isotopic analysis. Two distinct methods have been applied in recent isotopic studies of cassiterite: (a) reduction to tin metal with potassium cyanide (KCN) at high temperature (800 °C), with subsequent dissolution in HCl, and (b) reduction to a Sn solution with hydriodic acid (HI) at low temperature (100 °C). This study compares the effectiveness and accuracy of these two methods and contributes additional methodological details. The KCN method consistently yielded more Sn (> 70% in comparison with < 5%), does not appear to fractionate Sn isotopes at high temperatures over a 2‐hour period and produced consistent Sn isotope values at flux mass ratios of ≥ 4:1 (flux to mineral) with a minimum reduction time of 40 min. By means of a distillation experiment, it was demonstrated that HI could volatilise Sn, explaining the consistently low yields by this method. Furthermore, the distillation generated Sn vapour, which is up to 0.38‰ per mass unit different from the starting material, the largest induced Sn fractionation reported to date. Accordingly, the HI method is not recommended for cassiterite preparation for Sn isotopic analysis.     

A high-resolution geophysical investigation of sediment distribution controlled by catchment size and tides in a multi-basin turbid outwash fjord: Simpson Bay,Prince William Sound,Alaska

Surficial sediment distribution within Simpson Bay is a function of antecedent bedrock and recently deposited glacial geology, as well as active physical processes both within Simpson Bay and Prince William Sound (PWS). Simpson Bay is a turbid, outwash fjord located in northeastern PWS, Alaska. Freshwater from heavy precipitation, and the melting of high alpine glaciers enter the bay through bay head rivers and small shoreline creeks. The catchment has a high watershed/basin surface area ratio (∼8:1), and easily erodible bedrock that contribute to high sediment loads. The system can be divided into three discrete basins, each with specific morphologic and circulatory characters. Side scan sonar, swath bathymetry, and seismic profiles reveal that bathymetric highs are areas of outcropping glacial surfaces. High backscatter coupled with surface grab samples reveal these surfaces to be composed of coarse sediment and bedrock outcrops. Bathymetric lows are areas of low backscatter, and grab samples reveal these areas to be ponded deposits of organic-rich estuarine muds. The data provide evidence of terminal morainal bank systems, and glacial grounding line deposits at the mouth of the bay and rocky outcrops were identified as subsurface extensions of aerial rocky promontories. Radioisotope analyses of short cores reveal that the bay has an average accumulation rate of approx. 0.5 cm year⁻¹, but that this varies in function of the watershed/basin surface area ratios of the different basins. The interaction of tidal currents and sediment source drives sediment distribution in Simpson Bay. Hydrographic data reveal high spatial variability in surface and bottom currents throughout the bay. Subsurface currents are tide dominated, but generally weak (5–20 cm s⁻¹), while faster currents are found along shorelines, outcrops, and bathymetric highs. Bathymetric data reveal steep slopes with little to no modern sediment throughout the bay, suggesting lack of deposition due to tidal currents.     

Spectroscopic analysis (FTIR, Raman) of water in mafic and intermediate glasses and glass inclusions

Micro-Raman spectroscopy, even though a very promising technique, is not still routinely applied to analyse H₂O in silicate glasses. The accuracy of Raman water determinations critically depends on the capability to predict and take into account both the matrix effects (bulk glass composition) and the analytical conditions on band intensities. On the other hand, micro-Fourier transform infrared spectroscopy is commonly used to measure the hydrous absorbing species (e.g., hydroxyl OH⁻ and molecular H₂O) in natural glasses, but requires critical assumptions for the study of crystal-hosted glasses. Here, we quantify for the first time the matrix effect of Raman external calibration procedures for the quantification of the total H₂O content (H₂O_T = OH⁻ + H₂O_m) in natural silicate glasses. The procedures are based on the calibration of either the absolute (external calibration) or scaled (parameterisation) intensity of the 3550 cm⁻¹ band. A total of 67 mafic (basanite, basalt) and intermediate (andesite) glasses hosted in olivines, having between 0.2 and 4.8 wt% of H₂O, was analysed. Our new dataset demonstrates, for given water content, the height (intensity) of Raman H₂O_T band depends on glass density, reflectance and water environment. Hence this matrix effect must be considered in the quantification of H₂O by Raman spectroscopy irrespective of the procedure, whereas the parameterisation mainly helps to predict and verify the self-consistency of the Raman results. In addition, to validate the capability of the micro-Raman to accurately determine the H₂O content of multicomponent aluminosilicate glasses, a subset of 23 glasses was analysed by both micro-Raman and micro-FTIR spectroscopy using the band at 3550 cm⁻¹. We provide new FTIR absorptivity coefficients (ε₃₅₅₀) for basalt (62.80 ± 0.8 L mol⁻¹ cm⁻¹) and basanite (43.96 ± 0.6 L mol⁻¹ cm⁻¹). These values, together with an exhaustive review of literature data, confirm the non-linear decline of the FTIR absorptivity coefficient (ε₃₅₅₀) as the glass depolymerisation increases. We demonstrate the good agreement between micro-FTIR and micro-Raman determination of H₂O in silicate glasses when the matrix effects are properly considered.     

The Cloud Ice Mountain Experiment (CIME) 1998: experiment overview and modelling of the microphysical processes during the seeding by isentropic gas expansion

The second field campaign of the Cloud Ice Mountain Experiment (CIME) project took place in February 1998 on the mountain Puy de Dôme in the centre of France. The content of residual aerosol particles, of H₂O₂ and NH₃ in cloud droplets was evaluated by evaporating the drops larger than 5 μm in a Counterflow Virtual Impactor (CVI) and by measuring the residual particle concentration and the released gas content. The same trace species were studied behind a round jet impactor for the complementary interstitial aerosol particles smaller than 5 μm diameter. In a second step of experiments, the ambient supercooled cloud was converted to a mixed phase cloud by seeding the cloud with ice particles by the gas release from pressurised gas bottles. A comparison between the physical and chemical characteristics of liquid drops and ice particles allows a study of the fate of the trace constituents during the presence of ice crystals in the cloud.In the present paper, an overview is given of the CIME 98 experiment and the instrumentation deployed. The meteorological situation during the experiment was analysed with the help of a cloud scale model. The microphysics processes and the behaviour of the scavenged aerosol particles before and during seeding are analysed with the detailed microphysical model ExMix. The simulation results agreed well with the observations and confirmed the assumption that the Bergeron–Findeisen process was dominating during seeding and was influencing the partitioning of aerosol particles between drops and ice crystals. The results of the CIME 98 experiment give an insight on microphysical changes, redistribution of aerosol particles and cloud chemistry during the Bergeron–Findeisen process when acting also in natural clouds.     

Seismic displacement demand prediction in non-linear domain: Optimization of the N2 method

In Europe, computation of displacement demand for seismic assessment of existing buildings is essentially based on a simplified formulation of the N2 method as prescribed by Eurocode 8(EC8). However, a lack of accuracy of the N2 method in certain conditions has been pointed out by several studies. This paper addresses the assessment of effectiveness of the N2 method in seismic displacement demand determination in non-linear domain. The objective of this work is to investigate the accuracy of the N2 method through comparison with displacement demands computed using non-linear timehistory analysis(NLTHA). Results show that the original N2 method may lead to overestimation or underestimation of displacement demand predictions. This may affect results of mechanical model-based assessment of seismic vulnerability at an urban scale. Hence, the second part of this paper addresses an improvement of the N2 method formula by empirical evaluation of NLTHA results based on EC8 ground-classes. This task is formulated as a mathematical programming problem in which coefficients are obtained by minimizing the overall discrepancy between NLTHA and modified formula results. Various settings of the mathematical programming problem have been solved using a global optimization metaheuristic. An extensive comparison between the original N2 method formulation and optimized formulae highlights benefits of the strategy.     

Remote sensing of Mercury's magnetospheric plasma environment via energetic neutral atoms imaging

The potentiality of exploring the plasma dynamical behavior of Mercury via energetic neutral atom (ENA) analysis is briefly discussed. The analysis of the Hermean ENA would give clues to the global magnetospheric configuration, clarifying crucial issues like the existence of radiation belts. In order to demonstrate the feasibility of ENA imaging in the Hermean environment, an empirical model of the equatorial proton distributions, originally developed for the Earth's magnetosphere, is here modified in order to be applicable to the expected conditions, according to the present knowledge of the planet. Simulations of the ENA radiation are then presented and discussed. It is shown that the simulated radiation is significantly high, especially when pointing the field of view towards the Hermean surface. The way of determining the proton flux distributions and the cross-tail electric field characteristics by using the ENA signal is discussed. Finally, a preliminary design of an ENA detector suitable for the ESA mission to Mercury ‘BepiColombo’ is proposed and discussed.