Voluminous Storage and Rapid Magma Ascent Beneath La Palma Revealed by Seismic Tomography

Seismic tomography provides a window into magmatic plumbing systems; however, obtaining sufficient data for ‘real-time’ imaging is challenging. Until now, syn-eruptive tomography has not been successfully demonstrated. For the first time, we obtained high-resolution images of Earth's interior during an ongoing volcanic eruption. We used data from 11,349 earthquakes, most of which during La Palma eruption (19 September-13 December, 2021), to perform travel-time seismic tomography. We present high-precision earthquake relocations and 3D distributions of P and S-wave velocities highlighting the geometry of magma sources. We identified three distinct structures: (1) a shallow localised region (< 3 km) of hydrothermal alteration; (2) spatially extensive, consolidated, oceanic crust extending to ~10 km depth and; (3) a large (> 400 km3) sub-crustal magma-filled rock volume intrusion extending from ~7 to 25 km depth. Our results suggest that this large magma reservoir feeds the La Palma eruption continuously for almost three months. Prior to eruption onset, magma ascended from ~10 km depth to the surface in < 7 days. In the upper 3 km, melt migration is along the western contact between consolidated oceanic crust and altered hydrothermal material. Similar structural weaknesses along the eastern contact could potentially cause new eruptive centres in the future.

Seismic Attenuation Anisotropy in the Southernmost Part of the Taupo Volcanic Zone, North Island, New Zealand

<p>We investigate the mechanisms of seismic anisotropy and attenuation (1/Q) beneath the southernmost part of the Taupo Volcanic Zone (TVZ) by computing variations in S-wave attenuation factors with the direction of wave polarization. We rotate pairs of horizontal components in steps of 22.5◦ from 0◦ to 67.5◦ and into the radial and transverse directions to search for the optimal separation of the attenuation curves and thereby determine an anisotropy symmetry system. The frequency dependence of Q for the rotated S-waves is estimated by means of the non-parametric generalized inversion technique (GIT) of Castro et al. (1990) using shallow earthquakes (< 40 km depth) recorded by GeoNet within 100 km of Mt. Ruapehu. To analyze the effects on computed attenuation properties of source locations, we divide our dataset into two groups: a “TVZ” group containing earthquakes within the TVZ in a distance range of 5–55 km and a “non-TVZ” group containing earthquakes outside the TVZ in a distance range of 5–50 km. To measure Q, we compute the spectral amplitude decay with distance in terms of empirical functions at 20 separate frequencies in the frequency bands 2–10 Hz and 2– 12 Hz for the TVZ and non-TVZ datasets respectively. We construct homogeneous and two-layer Q models for the TVZ dataset based on characteristic features of the attenuation function, while for outside TVZ we only analyse a homogeneous Q model. The homogeneous Q models obtained for the two datasets indicate that S-waves are more attenuated within the TVZ than outside. The homogeneous Q model for the TVZ dataset reveals that the S-wave is anisotropic at high frequencies ( f > 6 Hz) along N–S/E– W directions with the relation QSE ( f ) = (6.15±1.22) f (1.73±0.12) and QSN ( f ) = (4.14± 1.26) f (2.06±0.14), while the non-TVZ dataset shows a weak frequency dependence of attenuation anisotropy at low frequencies in NE–SW/SE–NW directions giving the power law function QSNE ( f ) = (50.93±1.18) f (0.20±0.10) and QSSE ( f ) = (22.60±1.10) f (0.53±0.06). Here, the uncertainty estimates are 95% confidence intervals. To investigate the variation of attenuation anisotropy with depth within the TVZ, we first calculate Q along propagation paths (< 25 km, which corresponds to a maximum turning point depth of 9 km ) and then using paths of 25–55 km length. Small attenuation anisotropy with low attenuation in the N–S direction for the upper crust of TVZ may be related to heterogenous structure as reported by previous studies. Attenuation anisotropy in the northwest direction yielding lower attenuation inferred for the deeper crust suggests the presence of connected melt aligned with the extension direction of TVZ .</p>

Seismic Attenuation Anisotropy in the Southernmost Part of the Taupo Volcanic Zone, North Island, New Zealand

<p>We investigate the mechanisms of seismic anisotropy and attenuation (1/Q) beneath the southernmost part of the Taupo Volcanic Zone (TVZ) by computing variations in S-wave attenuation factors with the direction of wave polarization. We rotate pairs of horizontal components in steps of 22.5◦ from 0◦ to 67.5◦ and into the radial and transverse directions to search for the optimal separation of the attenuation curves and thereby determine an anisotropy symmetry system. The frequency dependence of Q for the rotated S-waves is estimated by means of the non-parametric generalized inversion technique (GIT) of Castro et al. (1990) using shallow earthquakes (< 40 km depth) recorded by GeoNet within 100 km of Mt. Ruapehu. To analyze the effects on computed attenuation properties of source locations, we divide our dataset into two groups: a “TVZ” group containing earthquakes within the TVZ in a distance range of 5–55 km and a “non-TVZ” group containing earthquakes outside the TVZ in a distance range of 5–50 km. To measure Q, we compute the spectral amplitude decay with distance in terms of empirical functions at 20 separate frequencies in the frequency bands 2–10 Hz and 2– 12 Hz for the TVZ and non-TVZ datasets respectively. We construct homogeneous and two-layer Q models for the TVZ dataset based on characteristic features of the attenuation function, while for outside TVZ we only analyse a homogeneous Q model. The homogeneous Q models obtained for the two datasets indicate that S-waves are more attenuated within the TVZ than outside. The homogeneous Q model for the TVZ dataset reveals that the S-wave is anisotropic at high frequencies ( f > 6 Hz) along N–S/E– W directions with the relation QSE ( f ) = (6.15±1.22) f (1.73±0.12) and QSN ( f ) = (4.14± 1.26) f (2.06±0.14), while the non-TVZ dataset shows a weak frequency dependence of attenuation anisotropy at low frequencies in NE–SW/SE–NW directions giving the power law function QSNE ( f ) = (50.93±1.18) f (0.20±0.10) and QSSE ( f ) = (22.60±1.10) f (0.53±0.06). Here, the uncertainty estimates are 95% confidence intervals. To investigate the variation of attenuation anisotropy with depth within the TVZ, we first calculate Q along propagation paths (< 25 km, which corresponds to a maximum turning point depth of 9 km ) and then using paths of 25–55 km length. Small attenuation anisotropy with low attenuation in the N–S direction for the upper crust of TVZ may be related to heterogenous structure as reported by previous studies. Attenuation anisotropy in the northwest direction yielding lower attenuation inferred for the deeper crust suggests the presence of connected melt aligned with the extension direction of TVZ .</p>

Changes in soil conditions before and after earthquakes at a repetitive soil liquefaction site in Taiwan

Beishih Village of Hsinhua Township in southern Taiwan is a unique location for studying soil liquefaction. Soil liquefaction was observed at the same site after earthquakes in 1946, 2010, and 2016, each of which had a Richter magnitude greater than six. This recurrence provides an opportunity for analyzing soil condition variations resulting from soil liquefaction. Seismic data sets were collected in 2011, 2014, 2016, and 2017. We used seismic refraction tomography and the multichannel analysis of surface waves to estimate P- and S-wave velocities. In S-wave velocity profiles, low shear velocity zones were located beneath sand volcanoes shortly after two earthquakes and disappeared 4 years after a 2010 earthquake. However, the P-wave velocity is less sensitive to soil condition changes, possibly because groundwater obscures the effect of soil liquefaction on velocity profiles. In addition, we used seismic wave velocities to determine the importance of soil properties such as Poisson’s ratio, shear modulus, and porosity to identify the cause of the low shear velocity zone. Notably, although porosity decreased after soil grain rearrangement, sand and clay mixing increased the Poisson’s ratio, reducing the shear modulus of the soil. In addition, a soil layer between 2 and 7 m and a deeper layer below 10 m that resulted in sand volcanoes were both liquefied. We also considered how the evaluation of soil liquefaction potential could be affected by long-term variations in soil conditions and changes resulting from liquefaction. The factor of safety was used to evaluate the liquefaction potential of the site. The results revealed that the assessment conducted long after the earthquake underestimated risk because the soil developed shear strength after the earthquake.

Ultrasonic and seismic constraints on crystallographic preferred orientations of the Priestley Glacier shear margin, Antarctica

Crystallographic preferred orientations (CPOs) are particularly important in controlling the mechanical properties of glacial shear margins. Logistical and safety considerations often make direct sampling of shear margins difficult and geophysical measurements are commonly used to constrain the CPOs. We present here the first direct comparison of seismic and ultrasonic data with measured CPOs in a polar shear margin. The measured CPO from ice samples from a 58 m deep borehole in the left lateral shear margin of the Priestley Glacier, Antarctica, is dominated by horizontal c-axes aligned sub-perpendicular to flow. A vertical seismic profile experiment with hammer shots up to 50 m away from the borehole, in four different azimuthal directions, shows velocity anisotropy of both P-waves and S-waves. Matching P-wave data to the anisotropy corresponding to CPO models defined by horizontally aligned c-axes gives two possible solutions for c-axis azimuth, one of which matches the c-axis measurements. If both P-wave and S-wave data are used, there is one best fit for azimuth and intensity of c-axis alignment that matches well the measurements. Azimuthal P-wave and S-wave ultrasonic data recorded in the laboratory on the ice core show clear anisotropy that matches that predicted from the CPO of the samples. With good quality data, azimuthal increments of 30° or less will constrain well the orientation and intensity of c-axis alignment. Our experiments provide a good framework for planning seismic surveys aimed at constraining the anisotropy of shear margins.

NIOM Analysis of Earthquake Observation Records in a Rockfill Dam

ABSTRACT We applied the normalized input–output minimization method (a method developed for the analysis of propagation times in vertical array records) to long-term earthquake observation records from Aratozawa Dam (in Kurihara, Miyagi prefecture, Japan), spanning the period from July 1992 to December 2019 to determine the propagation velocity of seismic waves in the embankment, and investigated changes in soil properties. As a result, we showed that (1) the velocities of S and P waves in the upper section were 449 and 993 m/s, respectively, prior to the strong earthquake motions derived from earthquake records from January 1997 through October 2001, whereas 608 and 1538, respectively, in the lower section, (2) in the Iwate–Miyagi Nairiku earthquake, the S-wave velocity in the upper section decreased to 158 m/s in the principal shock, and (3) in subsequent minor earthquakes the propagation velocity increased more or less in proportion with the logarithm of the number of elapsed days, requiring three years or longer to return to the initial value, (4) although similar changes were observed in the Great East Japan earthquake of 2011, the reduction in propagation velocity that remained after the principal shock was smaller than in the case of the Iwate–Miyagi Nairiku earthquake, and it was judged that there were no large effects on the dam body such as those that occurred in the Iwate–Miyagi Nairiku earthquake, and furthermore (5) in the principal shock of the Iwate–Miyagi Nairiku earthquake, the shear modulus in the upper part of the dam body decreased from 400 to 50 MPa (with a maximum shear strain of 10−3), resulting in more pronounced changes than in the lower section, whereas the damping ratio increased by at least 10% in the lower section during the principal shock of the Iwate–Miyagi Nairiku earthquake, resulting in much greater changes than in the upper section.

Pressure-dependent topological superconductivity on the surface of FeTe0.5Se0.5

FeTe1-xSex is a family of iron-based superconductors with its critical temperature (Tc) dependent on the composition of Se. A well-known Tc is 14.5 K for x = 0.45, which exhibits an s-wave superconducting gap between the topological superconducting surfaces states. Exchange interaction between the electrons has been proposed as the mechanism behind the formation of Cooper pairs for the sample of FeTe0.5Se0.5. In this article we provide further proof that exchange interaction, and hence the associated Tc, depends on the applied pressure on FeTe0.5Se0.5. Using density functional calculations for electrons and phonons and the Bardeen-Cooper-Schrieffer (BCS) theory for superconductivity, we found that Tc and superconducting gap for FeTe0.5Se0.5 soars under increasing compression, consistent with the results of experiment.

Elastic least-squares reverse time migration of steeply dipping structures using prismatic reflections

The migration of prismatic reflections can be used to delineate steeply dipping structures, which is crucial for oil and gas exploration and production. Elastic least-squares reverse time migration (ELSRTM), which considers the effects of elastic wave propagation, can be used to obtain reasonable subsurface reflectivity estimations and interpret multicomponent seismic data. In most cases, we can only obtain a smooth migration model. Thus, conventional ELSRTM, which is based on the first-order Born approximation, considers only primary reflections and cannot resolve steeply dipping structures. To address this issue, we develop an ELSRTM framework, called Pris-ELSRTM, which can jointly image primary and prismatic reflections in multicomponent seismic data. When Pris-ELSRTM is directly applied to multicomponent records, near-vertical structures can be resolved. However, the application of imaging conditions established for prismatic reflections to primary reflections destabilizes the process and leads to severe contamination of the results. Therefore, we further improve the Pris-ELSRTM framework by separating prismatic reflections from recorded multicomponent data. By removing artificial imaging conditions from the normal equation, primary and prismatic reflections can be imaged based on unique imaging conditions. The results of synthetic tests and field data applications demonstrate that the improved Pris-ELSRTM framework produces high-quality images of steeply dipping P- and S-wave velocity structures. However, it is difficult to delineate steep density structures because of the insensitivity of the density to prismatic reflections.

Innovation in the computer industry

A learning-curve approach has been employed to study inventiveness in the computer industry. The appearance of new models as well as the appearance of new computer manufacturers have been following typical logistic S-curves over the last 26 years of computer history. Innovation seems to come in S-waves, and even though many of today's major manufacturers are close to exhausting their maximum innovation potential in their present S-wave, the overall computer market is rather “young”. In contrast, the personal computer market is already beyond its maturity phase. Remarkably invariant over the 26 years considered, and through the complete range of computer sizes, remains the fact that for every five new computer models appearing on the market, there is one new computer company emerging also.