A new method for the determination of the specific kinetic energy (SKE) released to pyroclastic particles at magmatic fragmentation: theory and first experimental results

Brittle magmatic fragmentation plays a crucial role in explosive eruptions. It represents the starting point of hazardous explosive events that can affect large areas surrounding erupting volcanoes. Knowing the initial energy released during this fragmentation process is fundamental for the understanding of the subsequent dynamics of the eruptive gas-particle mixture and consequently for the forecasting of the erupting column’s behavior. The specific kinetic energy (SKE) of the particles quantifies the initial velocity shortly after the fragmentation and is therefore a necessary variable to model the gas-particle conduit flow and eruptive column regime. In this paper, we present a new method for its determination based on fragmentation experiments and identification of the timings of energy release. The results obtained on compositions representative for basaltic and phonolitic melts show a direct dependence on magma material properties: poorly vesiculated basaltic melts from Stromboli show the highest SKE values ranging from 7.3 to 11.8 kJ/kg, while experiments with highly vesiculated samples from Stromboli and Vesuvius result in lower SKE values (3.1 to 3.8 kJ/kg). The described methodology presents a useful tool for quantitative estimation of the kinetic energy release of magmatic fragmentation processes, which can contribute to the improvement of hazard assessment.     

Seismic wave velocities in the sedimentary cover of Poland: Borehole data compilation

A knowledge of seismic wave velocities in the sedimentary cover is of great importance for interpreting reflection and refraction seismic data, deep seismic soundings and regional and global seismic tomography. This is particularly true for regions characterized by significant thicknesses and a complex sedimentary cover structure. This paper presents the results of an analysis of seismic P-wave velocities in the sedimentary cover of Poland, a complex area of juxtaposition of major tectonic units: the Precambrian East European Craton, the Palaeozoic Platform of Central and Western Europe, and the Alpine orogen represented by the Carpathian Mountains. Based on vertical seismic profiling data from 1188 boreholes, the dependence of velocity versus depth was determined for regional geological units and for successions from the Tertiary and Quaternary to the Cambrian. The data have been approximated by polynomials, and velocity-depth formulas are given down to 6000 m depth. The velocities in the sedimentary cover have been compared with those from other areas in Europe.     

Thermal history and large scale differentiation of the Saturn’s satellite Rhea

Thermal history of Rhea from the beginning of accretion is investigated. We developed a numerical model of convection combined with the parameterized theory. Large scale melting of the satellite’s matter and gravitational differentiation of silicates from ices are included. The results are confronted with observational data from Cassini spacecraft that indicate minor differentiation of the satellite’s interior. We suggest that partial differentiation of the satellite’s interior is accompanied (or followed) by the process of light fraction uprising to the surface. The calculation indicates that the partial differentiation of the matter of the satellite’s interior is possible only for narrow range of parameters. In particular, we found that the time from the formation of CAI (calciumaluminum rich inclusions in chondrites) to the end of accretion of Rhea is in the range of 3–4 My.     

Interaction of geometry and mechanical property of trapezoidal sedimentary basins with incident <Emphasis Type="Italic">SH</Emphasis> waves

This paper investigated the effects of basin geometry and material property on the response of 2D trapezoidal sediment-filled basin to incident plane SH waves. Ten basin configurations with different geometries were developed, and then their seismic responses to both Ricker wavelets and seismic records were simulated by using an explicit finite difference scheme. The definition of deep/shallow basin, the precondition for the observation of prominent surface waves and the influential area of edge effects of the shallow basin were discussed quantitatively in this study. The followings were concluded: in the common velocity contrast range (v _s1/v _s2 < 10), the fundamental frequency a basin with W/H > 3.0 can be estimated approximately by 1D theory. The complexity of peak ground acceleration distribution pattern, the width of the most affected section as well as the amplitude of ground motion in the Edge Region increase with incident frequency. Prominent surface waves can only be observed when the incident wavelength is shorter than the critical wavelength λ _c . The interaction between incident wave and basin dynamic property plays a dominant role on the peak ground acceleration amplitude while the interaction between incident wave and geometry plays a more significant role on the peak ground acceleration distribution. For very shallow basin, different areas along the basin width are affected to different extents. Only a limited area close to the basin edge is influenced significantly. It is more feasible to propose spectral aggravation factor for different surface zones respectively than a uniform constant as a tool to calibrate the 1D-based design spectrum so as to take the basin effects into account.     

Modelling Freshwater Resources at the Global Scale: Challenges and Prospects

Quantification of spatially and temporally resolved water flows and water storage variations for all land areas of the globe is required to assess water resources, water scarcity and flood hazards, and to understand the Earth system. This quantification is done with the help of global hydrological models (GHMs). What are the challenges and prospects in the development and application of GHMs? Seven important challenges are presented. (1) Data scarcity makes quantification of human water use difficult even though significant progress has been achieved in the last decade. (2) Uncertainty of meteorological input data strongly affects model outputs. (3) The reaction of vegetation to changing climate and CO₂ concentrations is uncertain and not taken into account in most GHMs that serve to estimate climate change impacts. (4) Reasons for discrepant responses of GHMs to changing climate have yet to be identified. (5) More accurate estimates of monthly time series of water availability and use are needed to provide good indicators of water scarcity. (6) Integration of gradient-based groundwater modelling into GHMs is necessary for a better simulation of groundwater–surface water interactions and capillary rise. (7) Detection and attribution of human interference with freshwater systems by using GHMs are constrained by data of insufficient quality but also GHM uncertainty itself. Regarding prospects for progress, we propose to decrease the uncertainty of GHM output by making better use of in situ and remotely sensed observations of output variables such as river discharge or total water storage variations by multi-criteria validation, calibration or data assimilation. Finally, we present an initiative that works towards the vision of hyperresolution global hydrological modelling where GHM outputs would be provided at a 1-km resolution with reasonable accuracy.     

Earth’s Electromagnetic Environment

The natural spectrum of electromagnetic variations surrounding Earth extends across an enormous frequency range and is controlled by diverse physical processes. Electromagnetic (EM) induction studies make use of external field variations with frequencies ranging from the solar cycle which has been used for geomagnetic depth sounding through the 10\(^{-4}\)–10\(^4\) Hz frequency band widely used for magnetotelluric and audio-magnetotelluric studies. Above 10\(^4\) Hz, the EM spectrum is dominated by man-made signals. This review emphasizes electromagnetic sources at \(\sim\)1 Hz and higher, describing major differences in physical origin and structure of short- and long-period signals. The essential role of Earth’s internal magnetic field in defining the magnetosphere through its interactions with the solar wind and interplanetary magnetic field is briefly outlined. At its lower boundary, the magnetosphere is engaged in two-way interactions with the underlying ionosphere and neutral atmosphere. Extremely low-frequency (3 Hz–3 kHz) electromagnetic signals are generated in the form of sferics, lightning, and whistlers which can extend to frequencies as high as the VLF range (3–30 kHz).The roughly spherical dielectric cavity bounded by the ground and the ionosphere produces the Schumann resonance at around 8 Hz and its harmonics. A transverse resonance also occurs at 1.7–2.0 kHz arising from reflection off the variable height lower boundary of the ionosphere and exhibiting line splitting due to three-dimensional structure. Ground and satellite observations are discussed in the light of their contributions to understanding the global electric circuit and for EM induction studies.     

Present-Day Crustal Vertical Motion Around the Ordos Block Constrained by Precise Leveling and GPS Data

Precise leveling data observed in the period of 1970–2014 around the Ordos block were collected and processed to estimate present-day crustal vertical movement. Vertical rates of 6 GPS sites were employed as a priori constraints to define the reference frame. The velocity field shows that the interior of the Ordos block moves upward at a rate of 3 mm/a as a stable block. With respect to the central Ordos, the grabens and rifts around the Ordos block are undergoing subsidence, while the northeastern and southwestern Ordos uplift at the average rates of 2 and 1 mm/a, respectively. To the southeastern margin of the Ordos block, the Weihe and southern Shanxi grabens are subsiding at the rates of 4–6 mm/a. The subsidence of the Shanxi graben indicates that the graben is experiencing extensional movement on a long timescale. To the northwestern margin of the Ordos block, the Hetao and Yinchuan rifts are subsiding at the rates of 2–3 mm/a. A 2-D buried faulting model is used to infer the normal or reverse dip-slip rates. Our solution shows that most of the normal slip rates along the faults in the grabens and rifts are ~2 mm/a.     

Source parameters for the 2013–2015 earthquake sequence in Nógrád county,Hungary

Between 2013 June and 2015 January, 35 earthquakes with local magnitude M _L ranging from 1.1 to 4.2 occurred in Nógrád county, Hungary. This earthquake sequence represents above average seismic activity in the region and is the first one that was recorded by a significant number of three-component digital seismographs in the county. Using a Bayesian multiple-event location algorithm, we have estimated the hypocenters of 30 earthquakes with M _L ≥1.5. The events occurred in two small regions of a few squared kilometers: one to the east of Érsekvadkert and the other at Iliny. The uncertainty of the epicenters is about 1.5–1.7 km in the E-W direction and 1.8–2.1 km in the N-S direction at the 95 % confidence level. The estimated event depths are confined to the upper 3 km of the crust. We have successfully estimated the full moment tensors of 4 M _w ≥3.6 earthquakes using a probabilistic waveform inversion procedure. The non-double-couple components of the retrieved moment tensor solutions are statistically insignificant. The negligible amount of the isotropic component implies the tectonic nature of the investigated events. All of the analyzed earthquakes have strike-slip mechanism with either right-lateral slip on an approximately N-S striking or left-lateral movement on a roughly E-W striking nodal plane. The orientations of the obtained focal mechanisms are in good agreement with the main stress pattern published for the epicentral region. Both the P and T principal axes are horizontal, and the P axis is oriented along a NE-SW direction.     

Earthquake dynamics on circular faults: a review 1970–2015

We review the development of earthquake dynamics taken from the point of view of the origin of seismic radiation instead of the detailed study of rupture propagation on complex surfaces. Many features of seismic radiation can be explained by simple models that serve as elementary canonical problems. Some of these properties are very well known like the fact that at low frequencies, the seismic spectrum is proportional to the seismic moment. At high frequencies, on the other hand, radiation is generated by the motion of the rupture front, in particular stopping phases and geometrical obstacles (barriers). A rupture front moving at constant speed does not produce far-field radiation. For many practical applications, for determining source size and stress drop, for example, it is not necessary to determine geometrical details of the rupture. For such cases, a simple circular crack model is quite sufficient. An improvement on this method is to do dynamic inversion on simple, elliptical-shaped sources and letting the rupture start arbitrarily from a point on the fault. This problem can be solved nowadays with finite differences and a variety of inversion techniques.     

Three-dimensional admittance analysis of lithospheric elastic thickness over the Louisville Ridge

Using bathymetry and altimetric gravity anomalies, a 1° 9 1° lithospheric effective elastic thickness(Te) model over the Louisville Ridge and its adjacent regions is calculated using the moving window admittance technique. For comparison, three bathymetry models are used: general bathymetric charts of the oceans, SIO V15.1,and BAT_VGG. The results show that BAT_VGG is more suitable for calculating T e than the other two models. T e along the Louisville Ridge was re-evaluated. The southeast of the ridge has a medium Te of 10–20 km, while Te increases dramatically seaward of the Tonga-Kermadec trench as a result of the collision of the Pacific and IndoAustralian plates.