Examining Granger causality between atmospheric parameters and radon

In this paper, we study the relationship between atmospheric parameters (i.e., temperature and humidity) and radon data. We use the linear Granger causality in order to observe possible connections, on short and mid time scale periods, between radon time series and meteorological parameters that strongly influence radon emissions. The analysis suggests radon emission is not affected by these atmospheric parameters on short periods, while there is an evidence of Granger causality on mid periods.     

Geomaterials related to photovoltaics: a nanostructured Fe-bearing kuramite,Cu<Subscript>3</Subscript>SnS<Subscript>4</Subscript>

The successful synthesis of nanoparticles of Fe-bearing kuramite, (Cu,Fe)₃SnS₄, is reported in this study. Nanocrystalline powders were obtained through a mild, environmentally friendly and scalable solvothermal approach, in a single run. The sample was the object of a multidisciplinary investigation, including X-ray diffraction and absorption, scanning electron microscopy and microanalysis, electron paramagnetic resonance, diffuse reflectance and Mössbauer spectroscopy as well as SQUID magnetometry. The nanoparticles consist of pure Fe-bearing kuramite, exhibiting tetragonal structure. The valence state of the metal cations was assessed to be Cu⁺, Sn⁴⁺ and Fe³⁺. The material presents a band gap value of 1.6 eV, which is fully compatible with solar cell applications. The uptake of Fe by nanokuramite opens a compositional field where the physical properties can be tuned. We thus foster the application of Fe-bearing nanokuramite for photovoltaics and energy storage purposes.     

The Maoxian landslide as seen from space: detecting precursors of failure with Sentinel-1 data

Post-event Interferometric Synthetic Aperture Radar (InSAR) analysis on a stack of 45 C-band SAR images acquired by the ESA Sentinel-1 satellites from 9 October 2014 to 19 June 2017 allowed the identification of a clear precursory deformation signal for the Maoxian landslide (Mao County, Sichuan Province, China). The landslide occurred in the early morning of 24 June 2017 and killed more than 100 people in the village of Xinmo. Sentinel-1 images have been processed through an advanced multi-interferogram analysis capable of maximising the density of measurement points, generating ground deformation maps and displacement time series for an area of 460 km² straddling the Minjiang River and the Songping Gully. InSAR data clearly show the precursors of the slope failure in the source area of the Maoxian landslide, with a maximum displacement rate detected of 27 mm/year along the line of sight of the satellite. Deformation time series of measurement points identified within the main scarp of the landslide exhibit an acceleration starting from April 2017. A detailed time series analysis leads to the classification of different deformation behaviours. The Fukuzono method for forecasting the time of failure appear to be applicable to the displacement data exhibiting progressive acceleration. Results suggest that satellite radar data, systematically acquired over large areas with short revisiting time, could be used not only as a tool for mapping unstable areas, but also for landslide monitoring, at least for some typologies of sliding phenomena.     

Location Performance and Detection Magnitude Threshold of the Romanian National Seismic Network

Romania is an earthquake prone area with a few destructive earthquakes per century. The National Institute for Earth Physics carries out the seismic survey of Romania through the Romanian National Seismic Network (RNSN) consisting of 65 real-time seismic stations. Daily reports and monthly bulletins are delivered after routinely analyzing and processing the recorded data. In the present paper we applied the Seismic Network Evaluation through Simulation method for the RNSN configuration as it was in August 2011 to estimate the background noise level, assess the appropriateness of the velocity model adopted in routine location procedure, evaluate the hypocenter location uncertainty and determine the detection magnitude threshold. Areas of greater (southern Romania) and lower (Carpathians and Apuseni Mountains) background noise within the RNSN are identified by mapping the average power of noise in 1–12?Hz frequency range. The statistical study of the P and S phases residual times allow us to assess the appropriateness of the velocity model used in routine location. Both P- and S-wave velocity models can be optimized to improve the quality of the hypocenter location. As shown by our analysis, the RNSN is able to detect and locate earthquakes with M _L magnitude above 2.5 anywhere on the Romanian territory, except the border areas, such as the Crisana–Maramures seismic source zone. Merging data from both sides of the border significantly improves the quality of hypocenter location in these areas.     

Hybrid Inversion Method to Estimate Hydraulic Transmissivity by Combining Multiple-Point Statistics and a Direct Inversion Method

Inversion methods that rely on measurements of the hydraulic head h cannot capture the fine-scale variability of the hydraulic properties of an aquifer. This is particularly true for direct inversion methods, which have the further limitation of providing only deterministic results. On the other hand, stochastic simulation methods can reproduce the fine-scale heterogeneity but cannot directly incorporate information about the hydraulic gradient. In this work, a hybrid approach is proposed to join a direct inversion method (the comparison model method, CMM) and multiple-point statistics (MPS), for determination of a hydraulic transmissivity field T from a map of a reference hydraulic head \(h^\mathrm {(ref)}\) and a prior model of the heterogeneity (a training image). The hybrid approach was tested and compared with pure MPS and pure CMM approaches in a synthetic case study. Also, sensitivity analysis was performed to test the importance of the acceptance threshold \(\delta \), a simulation parameter that allows one to tune the influence of \(h^\mathrm {(ref)}\) on the final results. The transmissivity fields T obtained using the hybrid approach take into account information coming from the hydraulic gradient while simultaneously reproducing some of the fine-scale features provided by the training image. Furthermore, many realizations of T can be obtained thanks to the stochasticity of MPS. Nevertheless, it is not straightforward to exploit the correlation between the T maps provided by the CMM and the prior model introduced by the training image, because the former depends on the boundary conditions and flow settings. Another drawback is the growing number of simulation parameters introduced when combining two diverse methods. At the same time, this growing complexity opens new possibilities that deserve further investigation.     

Elastic properties of six silicate garnet end members from accurate ab initio simulations

The elastic properties of six silicate garnet end members, among the most important rock-forming minerals, are investigated here for the first time via accurate ab initio theoretical simulations. The Crystal program is used, which works within periodic boundary conditions and allows for all-electron basis sets to be adopted. From the computed elastic tensor, Christoffel’s equation is solved along a set of crystallographic directions in order to fully characterize the seismic wave velocity anisotropy in such materials. Polycrystalline isotropic aggregate elastic properties are derived from the computed single-crystal data via the Voigt-Reuss-Hill averaging procedure. Transferability of the elastic properties from end members to their solid solutions with different chemical compositions is also addressed.     

Earthquake‐induced turbidites triggered by sea level oscillations in the Upper Cretaceous and Paleocene of Italy

Earthquakes are widely recognized as triggers for turbidites, submarine debris flows and slumps. In tectonically active areas, surprisingly small changes in stress can trigger seismic events, implying that past sea level changes may be important in controlling the timing of seismicity and the occurrence of turbidites. We apply this idea in an analysis of turbidites from the upper Cretaceous–Paleocene Scaglia Rossa Formation of the Umbria‐Marches region of Italy. These turbidites are composed of resedimented foraminiferal tests derived from fluidizing deep‐water (∼1500 m), pelagic sediments; seismic triggering is the most likely triggering mechanism given this setting and composition. The timing of these turbidites (and associated synsedimentary slumps), constrained by biostratigraphy and magnetostratigraphy, reveals an unusual, non‐random temporal pattern that appears to correlate well with proposed eustatic fluctuations. This correlation between turbidites and eustatic fluctuations leads us to suggest that stress and pore fluid pressure changes associated with changing sea level may trigger periods of increased seismicity in the geological past.     

Time evolution of complex networks: commuting systems in insular Italy

The aim of this paper is to study the dynamics of the commuting system of two insular regions of Italy, Sardinia and Sicily, inspected as complex networks. The authors refer to a 20-year time period and take into account three census data sets about the work and study-driven inter-municipal origin-destination movements of residential inhabitants in 1981, 1991 and 2001. Since it is likely that the number of municipalities (in this case, the vertices of the system) does not display sharp variations, the authors direct the study to the variation of the properties emerging through both a topological and a weighted network representation of commuting in the time periods indicated.     

The use of iron chemical analysis of podzols to date the Late Pleistocene–Holocene deglaciation history of the Central Italian Alps

Podzols that have developed on glacial and periglacial features provide the opportunity to reconstruct glacial evolution after the Last Glacial Maximum (LGM) using different soil indices. Analysing 17 soils classified as podzol, we used the crystallinity ratio of free iron oxides (CRF) on both the A and Bs horizons, and absolute ages for the same landforms containing the soil profile, to create dating curves. Two equations were generated: age = 4566.9 × ln (CRF) + 1760 (1), and age = 3907 × ln (CRF) + 3508.2 (2). The reliability of the curves was evaluated with the Fe_o/Fe_d ratio, and with the difference of ages calculated using both equations. Equation (2) is considered more reliable because the A horizon may be influenced by new pedogenesis on the pre-existing podzol, leading to the development of a new type of soil. By dating the soils, we reconstructed the glacial history of the three main upper branches of the LGM Adda Glacier in the Central Italian Alps, specifically the Stelvio Pass area (ST), Gavia Pass area (GV), and the Val Viola valley (VV). Seven glacial advances were identified at 16.7–14.7 ka (phase I), 12.3 ka (phase II), 11 ka (phase III), 10–9.7 ka (phase IV), 9 ka (phase V), 7.5 ka (phase VI) and 5.3 ka (phase VII). The first five phases are chronologically similar to the main Late Pleistocene–Early Holocene phases recorded in the Central European Alps. The last two Holocene phases, which are both longer in duration than the Little Ice Age, are recorded in ST and GV. Interestingly, these phases generally are not recorded in the rest of the Central European Alps, where the late Holocene glaciers were smaller than their present size.