ALBERT

Description: We study a model of lava flow to determine its thermal and dynamic characteristics from thermal measurements of the lava at its surface. Mathematically this problem is reduced to solving an inverse boundary problem. Namely, using known conditions at one part of the model boundary we determine the missing condition at the remaining part of the boundary. We develop a numerical approach to the mathematical problem in the case of steady-state flow. Assuming that the temperature and the heat flow are prescribed at the upper surface of the model domain, we determine the flow characteristics in the entire model domain using a variational (adjoint) method. We have performed computations of model examples and showed that in the case of smooth input data the lava temperature and the flow velocity can be reconstructed with a high accuracy. As expected, a noise imposed on the smooth input data results in a less accurate solution, but still acceptable below some noise level. Also we analyse the influence of optimization methods on the solution convergence rate. The proposed method for reconstruction of physical parameters of lava flows can also be applied to other problems in geophysical fluid flows.

Description: Studies of discrete volcanic explosions, that usually last less than 2 or 3 min, have suggested that the partitioning of seismic-acoustic energy is likely related to a range of physical mechanisms that depend on magma properties and other physical constraints such as the location of the fragmentation surface. In this paper, we explore the energy partition of a paroxysmal eruptive phase of Tungurahua volcano that lasted for over 4 hr, on 2006 July 14–15, using seismic-acoustic information recorded by stations on its flanks (near field). We find evidence of a linear scaling between seismic and acoustic energies, with time-dependent intensities, during the sustained explosive phase of the eruption. Furthermore, we argue that this scaling can be explained by two different processes: (1) the fragmentation region ultimately acts as the common source of energy producing both direct seismic waves, that travel through the volcanic edifice, and direct acoustic waves coming from a disturbed atmosphere above the summit; (2) the coupling of acoustic waves with the ground to cause seismic waves. Both processes are concurrent, however we have found that the first one is dominant for seismic records below 4 Hz. Here we use the linear scaling of intensities to construct seismic and acoustic indices, which, we argue, could be used to track an ongoing eruption. Thus, especially in strong paroxysms that can produce pyroclastic flows, the index correlation and their levels can be used as quantitative monitoring parameters to assess the volcanic hazard in real time. Additionally, we suggest from the linear scaling that the source type for both cases, seismic and acoustic, is dipolar and dominant in the near field.

Description: Lying below Vatnajökull ice cap in Iceland, Bárðarbunga stratovolcano began experiencing wholesale caldera collapse in 2014 August 16, one of the largest such events recorded in the modern instrumental era. Simultaneous with this collapse is the initiation of a plate boundary rifting episode north of the caldera. Observations using the international constellation of radar satellites indicate rapid 50 cm d –1 subsidence of the glacier surface overlying the collapsing caldera and metre-scale crustal deformation in the active rift zone. Anomalous earthquakes around the rim of the caldera with highly nondouble-couple focal mechanisms provide a mechanical link to the dynamics of the collapsing magma chamber. A model of the collapse consistent with available geodetic and seismic observations suggests that the majority of the observed subsidence occurs aseismically via a deflating sill-like magma chamber.

Description: Deformation experiments on partially molten rocks in simple shear form melt bands at 20° to the shear plane instead of at the expected 45° principal compressive stress direction. These melt bands may play an important role in melt focusing in mid-ocean ridges. Such shallow bands are known to form for two-phase media under shear if strongly non-Newtonian power-law creep is employed for the solid phase, or anisotropy imposed. However laboratory experiments show that shallow bands occur regardless of creep mechanism, even in diffusion creep, which is nominally Newtonian. Here we propose that a couple of forms of two-phase damage allow for shallow melt bands even in diffusion creep.

Description: Progressive deformation of upper mantle rocks via dislocation creep causes their constituent crystals to take on a non-random orientation distribution (crystallographic preferred orientation or CPO) whose observable signatures include shear-wave splitting and azimuthal dependence of surface wave speeds. Comparison of these signatures with mantle flow models thus allows mantle dynamics to be unraveled on global and regional scales. However, existing self-consistent models of CPO evolution are computationally expensive when used with 3-D and/or time-dependent convection models. Here we propose a new method, called ANPAR, which is based on an analytical parametrization of the crystallographic spin predicted by the second-order (SO) self-consistent theory. Our parametrization runs 2–6  x  10 4 times faster than the SO model and fits its predictions for CPO and crystallographic spin with a variance reduction 〉99 per cent. We illustrate the ANPAR model predictions for the deformation of olivine with three dominant slip systems, (010)[100], (001)[100] and (010)[001], for three uniform deformations (uniaxial compression, pure shear and simple shear) and for a corner-flow model of a spreading mid-ocean ridge.

Description: We investigate the capability of Time Reversal Mirror (TRM) algorithm to image local acoustic sources (〈3.5 km) associated with complex, sustained volcanic eruptions. Accurate source localization for volcano infrasound (low-frequency acoustic waves) is often challenging due to pronounced volcanic topography and emergent arrivals of infrasound signals. While the accuracy of the conventional approaches (e.g. triangulation and semblance method) can be severely compromised by the complex volcanic settings, a TRM-based method may have the potential to properly image acoustic sources by the use of full waveform information and numerical modelling of the time-reversed wavefield. We apply the TRM algorithm to a pyroclastic-laden eruption (sustained for ~60 s) at Santiaguito Volcano, Guatemala, and show that an ordinary TRM operation can undergo significant reduction of its focusing power due to strong topographic propagation effects (e.g. reflection and diffraction). We propose a weighted imaging condition to compensate for complicated transmission loss of the time-reversed wavefield and demonstrate that the presented condition significantly improves the focusing quality of TRM in the presence of complex topography. The consequent TRM source images exhibit remarkable agreement with the visual observation of the eruption implying that the TRM method with a proper imaging condition can be used to localize and track acoustic sources associated with complex volcanic eruptions.

Description: The expansion of offshore renewable energy infrastructure and the need for trans-continental shelf power transmission require the use of submarine high-voltage (HV) cables. These cables have maximum operating surface temperatures of up to 70 °C and are typically buried 1–2 m beneath the seabed, within the wide range of substrates found on the continental shelf. However, the heat flow pattern and potential effects on the sedimentary environments around such anomalously high heat sources in the near-surface sediments are poorly understood. We present temperature measurements from a 2-D laboratory experiment representing a buried submarine HV cable, and identify the thermal regimes generated within typical unconsolidated shelf sediments—coarse silt, fine sand and very coarse sand. We used a large (2 x 2.5 m 2 ) tank filled with water-saturated spherical glass beads (ballotini) and instrumented with a buried heat source and 120 thermocouples to measure the time-dependent 2-D temperature distributions. The observed and corresponding Finite Element Method simulations of the steady state heat flow regimes and normalized radial temperature distributions were assessed. Our results show that the heat transfer and thus temperature fields generated from submarine HV cables buried within a range of sediments are highly variable. Coarse silts are shown to be purely conductive, producing temperature increases of 〉10 °C up to 40 cm from the source of 60 °C above ambient; fine sands demonstrate a transition from conductive to convective heat transfer between cf. 20 and 36 °C above ambient, with 〉10 °C heat increases occurring over a metre from the source of 55 °C above ambient; and very coarse sands exhibit dominantly convective heat transfer even at very low ( cf. 7 °C) operating temperatures and reaching temperatures of up to 18 °C above ambient at a metre from the source at surface temperatures of only 18 °C. These findings are important for the surrounding near-surface environments experiencing such high temperatures and may have significant implications for chemical and physical processes operating at the grain and subgrain scale; biological activity at both microfaunal and macrofaunal levels; and indeed the operational performance of the cables themselves, as convective heat transport would increase cable current ratings, something neglected in existing standards.

Description: We present a strategy to thoroughly investigate the effects of prominent topography on the surface tilt due to a spherical pressure source. We use Etna's topography as a case of study and, for different source positions, we compare the tilt fields calculated through (i) a 3-D boundary element method and (ii) analytical half-space solutions. We systematically determine (i) the source positions leading to the strongest tilt misfits when numerical and analytical results are compared and (ii) the surface areas where the strongest distortions in the tilt field are most likely to be observed. We also demonstrate that, under critical circumstances, in terms of respective positions of pressure source and observation points, results of inversion procedures aimed at retrieving the source parameters can be misleading, if tilt data are analysed using models that do not account for topography.

Description: The agreement between shear wave velocities for the Earth's inner core observed from seismology with those derived from mineral physics is considerably worse than for any other region of the Earth. Furthermore, there is still debate as to the phase of iron present in the inner core, particularly when alloying with nickel and light elements is taken into account. To investigate the extent to which the mismatch between seismology and mineral physics is a function of either crystal structure and/or the amount of nickel present, we have used ab initio molecular dynamics simulations to calculate the elastic constants and seismic velocities ( V p and V s ) of face centred cubic (fcc) iron at Earth's inner core pressures (360 GPa) and at temperatures up to ~7000 K. We find that V p for fcc iron (fcc-Fe) is very similar to that for hexagonal close packed (hcp) iron at all temperatures. In contrast, V s for fcc-Fe is significantly higher than in hcp-Fe, with the difference increasing with increasing temperature; the difference between V s for the core (from seismology) and V s for fcc-Fe exceeds 40 per cent. These results are consistent with previous work at lower temperatures. We have also investigated the effect of 6.5 and 13 atm% Ni in fcc-Fe. We find that Ni only slightly reduces V p and V s (e.g. by 2 per cent in V s for 13 atm% Ni at 5500 K), and cannot account for the difference between the velocities observed in the core and those of pure fcc-Fe. We also tried to examine pre-melting behaviour in fcc-Fe, as reported in hcp-Fe by extending the study to very high temperatures (at which superheating may occur). However, we find that fcc-Fe spontaneously transforms to other hcp-like structures before melting; two hcp-like structures were found, both of hexagonal symmetry, which may most easily be regarded as being derived from an hcp crystal with stacking faults. That the structure did not transform to a true hcp phase is likely as a consequence of the limited size of the simulation box (108 atoms). At 360 GPa, in pure fcc-Fe, we find that the transition from fcc to the hcp-like structure occurs at 7000 K, whereas in the Ni bearing system, the transition occurs at higher temperature (7250 K). This reinforces previous work showing that fcc-Fe might transform to hcp-Fe just before melting, and that Ni tends to stabilize the fcc structure with respect to hcp.

Description: Viscoelastic behaviour of materials in nature is observed in post-event deformations due to seismic or volcanic activities. In this paper, by adopting the correspondence principle, we propose an inelastic model to predict first the Laplace-domain response of a transversely isotropic viscoelastic half-space due to a shear or tensile fault of polygonal shape. The displacement and stress fields in the time domain are then obtained using an efficient and accurate algorithm for the inverse Laplace transform. Numerical examples are presented to validate the proposed solution and to show the viscoelastic displacement and stress fields due to a strike-slip, dip-slip and tensile fault of rectangular shape. The obtained results indicate that both viscoelasticity and transverse isotropy play significant roles in the viscoelastic response of the half-space due to faults, which could be used as benchmarks for the future numerical analysis of realistic post-seismic or volcanic event.