Analysis of glacial earthquakes

In 2003, Ekström et al. reported on the detection of a new class of earthquakes that occur in glaciated regions, with the vast majority being in Greenland. The events have a characteristic radiation pattern and lack the high-frequency content typical of tectonic earthquakes. It was proposed that the events correspond to large and sudden sliding motion of glaciers. Here we present an analysis of all 184 such events detected in Greenland between 1993 and 2005. Fitting the teleseismic long-period surface waves to a landslide model of the source, we obtain improved locations, timing, force amplitudes, and force directions. After relocation, the events cluster into seven regions, all of which correspond to regions of very high ice flow and most of which are named outlet glaciers. These regions are Daugaard Jensen Glacier, Kangerdlugssuaq Glacier, Helheim Glacier, the southeast Greenland glaciers, the northwest Greenland glaciers, Rinks Isbrae, and Jakobshavn Isbrae. Event amplitudes range from 0.1 to 2.0 × 10^(14) kg m. Force directions are consistent with sliding in the direction of glacial flow over a period of about 50 s. Each region has a different temporal distribution of events. All glaciers are more productive in the summer, but have their peak activity in different months. Over the study period, Kangerdlugssuaq has had a constant number of events each year, whereas Jakobshavn had most events in 1998–1999, and the number of events in Helheim and the northwest Greenland glaciers has increased substantially between 1993 and 2005. The size distribution of events in Kangerdlugssuaq is peaked above the detection threshold, suggesting that glacial earthquakes have a characteristic size

Earthquakes along Eltanin transform system, SE Pacific Ocean: fault segments characterized by strong and poor seismic coupling and implications for long-term earthquake prediction

Centroid moment tensor solutions are recomputed for 190 earthquakes from 1976 to 2010 along the Heezen, Tharp and Hollister transform faults of the Eltanin system using a 3-D seismic velocity model. The total length of the three en echelon faults is nearly 1000 km; each is characterized by fast long-term rates of displacement of about 80 mm yr^-1. Strike-slip faulting with moment magnitudes Mw up to 6.4 characterizes most of these events. The few involving normal faulting are located up to 40 km on either side of the transforms and involve extension nearly normal to the transforms. This partitioning of slip likely results from changes during the last few million years in the Euler pole for relative motion between the Antarctic and Pacific plates. Some parts of the Heezen and Tharp transforms exhibit strong seismic coupling but others were aseismic at the resolution of our study, Mw > 5.0-5.5. Earthquakes were not found along nearby fast spreading ridges at that resolution. We calculate downdip widths of seismic coupling of about 5 km for four strongly coupled segments from observed moment rates and lengths along strike assuming earthquake activity accounts for the entire plate motion. Major differences in seismic coupling along strike are not in accord with common thermal models of plate cooling but instead are attributed to varying degrees of metamorphism, rock type and effective normal stress and possibly to the presence of short intratransform spreading centres. One 30-42-km-long segment of the Heezen transform that appears to be an isolated well-coupled asperity has ruptured in eight earthquakes of Mw 5.9-6.1 quasi-periodically with a coefficient of variation of 0.26 every 4.0 ± 1.0 yr. Other well-coupled fault segments, which were sites with earthquakes up to Mw 6.39 and fewer events since 1976, have average repeat times of about 7-24 yr. The fast rate of plate motion, maximum size of events and relatively short repeat times make these fault segments a good laboratory for research on quasi-periodic behaviour and earthquake prediction

Simple scaling of catastrophic landslide dynamics

Catastrophic landslides involve the acceleration and deceleration of millions of tons of rock and debris in response to the forces of gravity and dissipation. Their unpredictability and frequent location in remote areas have made observations of their dynamics rare. Through real-time detection and inverse modeling of teleseismic data, we show that landslide dynamics are primarily determined by the length scale of the source mass. When combined with geometric constraints from satellite imagery, the seismically determined landslide force histories yield estimates of landslide duration, momenta, potential energy loss, mass, and runout trajectory. Measurements of these dynamical properties for 29 teleseismogenic landslides are consistent with a simple acceleration model in which height drop and rupture depth scale with the length of the failing slope

Love and Rayleigh phase-velocity maps, 5–40 s, of the western and central USA from USArray data

Continuous data recorded on more than 1600 USArray seismic stations operating in the western and central US between 2006 and 2012 are used to map phase velocities of Love and Rayleigh waves at short periods (5–40 s) using a noise-correlation technique. Vertical and transverse records for all station pairs separated by less than 600 km are cross correlated in 4-h-long segments, and the resulting spectra are stacked for the time period of station operation. Dispersion curves are determined from the locations of zeros in the real component of the correlation spectra using a method based on that of Aki (1957). Phase-velocity maps expanded on a 0.25°-by-0.25° pixel grid are estimated by inversion of the phase-velocity measurements. Comparison with predicted phase-velocity maps based on the crustal model CRUST2.0 combined with the mantle model ND08 shows good agreement at the longer periods. Strong slow anomalies (greater than 25%) in the short-period maps are geographically correlated with basins and regions of thick sedimentary cover. The strength of these anomalies is not well predicted by existing crustal-velocity models

Quantify and account for field reference errors in forest remote sensing studies

Field inventoried data are often used as references (ground truth) in forest remote sensing studies. However, the reference values are affected by various kinds of errors, which tend to make the reported accuracies of the remote sensing-based predictions worse than they are. The more accurate the remote sensing techniques are becoming, the more pronounced this problem will be. This paper addresses the impact of uncertainties in field reference data due to measurement errors, model errors, and position errors when evaluating the accuracy of biomass predictions from airborne laser scanning at plot level. We present novel theoretical analysis methods that take the interactions of the error sources into account. Further, an error characterization model (ECM) is used to describe the error structure of the remote sensing-based predictions, and we show how the parameters of the ECM can be adjusted when field references contain errors. We also show how root mean square error (RMSE) estimates can be adjusted. Based on data from Scandinavian forests, we conclude that the field reference errors have an impact on the remote sensing-based predictions. By accounting for these errors the RMSE of the remote sensing-based predictions was reduced by 6-18%. The most influential sources of error in the field references were found to be the residual errors of the allometric biomass model and the field plot position errors. Together, these two sources accounted for 97% of the variance while measurement errors and biomass model parameter uncertainties were negligible in our study

A comparison of approaches to the prediction of surface wave amplitude

A controlled experiment is performed to investigate how assumptions and simplifications in the measurement and analysis of surface wave amplitudes affect inferred attenuation variations in the mantle. Synthetic seismograms are generated using a spectral-element method for 42 earthquakes, 134 receiver locations and two earth models, both of which contain 3-D elastic properties and 1-D anelastic properties. Fundamental-mode Rayleigh-wave amplitudes are measured at periods of 50, 75 and 125 s for 4749 paths. The amplitudes are measured with respect to a reference waveform based on 1-D Earth structure, and thus amplitude observations that are not equal to unity can be attributed to differences in the computation of the spectral-element and reference waveforms or to uncertainties in the amplitude measurements themselves. Calculation of earthquake source excitation in the 3-D earth model versus the 1-D earth model has a significant effect on the amplitudes, especially at shorter periods, and variations in the average amplitude for each event are well explained by the effect of Earth structure at the event location on the source excitation. The effect of local Earth structure at the receiver location on the amplitude is, for most paths, much smaller than for the source amplitude. After correcting for source and receiver effects on amplitude, the remaining signal is compared to predictions of elastic focusing effects using the great-circle ray approximation, exact ray theory (ERT) and finite-frequency theory (FFT). We find that, for the earth models we have tested, ERT provides the best fit at 50 s, and FFT is most successful at 75 and 125 s, indicating that the broad zone of surface wave sensitivity cannot be neglected for the longer periods in our experiment. The bias introduced into attenuation models by focusing effects, which is assessed by inverting the measured amplitudes for 2-D attenuation maps, is most important at high spherical-harmonic degrees. Unaccounted-for scattering of seismic energy may slightly (<5 per cent) raise average global attenuation values at short periods but has no detectable effect at longer periods. The findings of this study also provide a set of guidelines for handling source, receiver and focusing effects that can be applied to surface wave amplitudes measured for the real Earth

The European Upper Mantle as Seen by Surface Waves

We derive a global, three-dimensional tomographic model of horizontally and vertically polarized shear velocities in the upper mantle. The model is based on a recently updated global database of Love- and Rayleigh-wave fundamental-mode phase-anomaly observations, with a good global coverage and a particularly dense coverage over Europe and the Mediterranean basin (broadband stations from the Swiss and German seismic networks). The model parameterization is accordingly finer within this region than over the rest of the globe. The large-scale, global structure of our model is very well correlated with that of earlier shear-velocity tomography models, based both on body- and surface-wave observations. At the regional scale, within the region of interest, correlation is complicated by the different resolution limits associated to different databases (surface waves, compressional waves, shear waves), and, accordingly, to different models; while a certain agreement appears to exist for what concerns the grand tectonic features in the area, heterogeneities of smaller scale are less robustly determined. Our new model is only one step towards the identification of a consensus model of European/Mediterranean upper-mantle structure: on the basis of the findings discussed here, we expect that important improvements will soon result from the combination, in new tomographic inversions, of fundamental-mode phase-anomaly data like ours with observations of surface-wave overtones, of body-wave travel times, of ambient "noise”, and by accounting for an a-priori model of crustal structure more highly resolved than the one employed her

Cavity-free vacuum-Rabi splitting in circuit quantum acoustodynamics

Artificial atoms coupled to surface acoustic waves (SAWs) have played a crucial role in the recent development of circuit quantum acoustodynamics (cQAD). In this paper, we have investigated the interaction of an artificial atom and SAWs beyond the weak coupling regime, focusing on the role of the interdigital transducer (IDT) that enables the coupling. We find a parameter regime in which the IDT acts as a cavity for the atom, rather than an antenna. In other words, the atom forms its own cavity. Similar to an atom coupled to an explicit cavity, this regime is characterized by vacuum-Rabi splitting, as the atom hybridizes with the phononic vacuum inside the IDT. This hybridization is possible because of the interdigitated coupling, which has a large spatial extension, and the slow propagation speed of SAWs. We work out a criterion for entering this regime from a model based on standard circuit-quantization techniques, taking only material parameters as inputs. Most notably, we find this regime hard to avoid for an atom on top of a strong piezoelectric material, such as LiNbO$_3$. The SAW-coupled atom on top of LiNbO$_3$ can thus be regarded as an atom-cavity-bath system. On weaker piezoelectric materials, the number of IDT electrodes need to be large in order to reach this regime.Comment: 11 pages, 5 figure

Propagating phonons coupled to an artificial atom

Quantum information can be stored in micromechanical resonators, encoded as quanta of vibration known as phonons. The vibrational motion is then restricted to the stationary eigenmodes of the resonator, which thus serves as local storage for phonons. In contrast, we couple propagating phonons to an artificial atom in the quantum regime, and reproduce findings from quantum optics with sound taking over the role of light. Our results highlight the similarities between phonons and photons, but also point to new opportunities arising from the unique features of quantum mechanical sound. The low propagation speed of phonons should enable new dynamic schemes for processing quantum information, and the short wavelength allows regimes of atomic physics to be explored which cannot be reached in photonic systems.Comment: 30 pages, 6 figures, 1 tabl

Global observation of vertical-CLVD earthquakes at active volcanoes

Some of the largest and most anomalous volcanic earthquakes have non-double-couple focal mechanisms. Here, we investigate the link between volcanic unrest and the occurrence of non-double-couple earthquakes with dominant vertical tension or pressure axes, known as vertical compensated-linear-vector-dipole (vertical-CLVD) earthquakes. We determine focal mechanisms for 313 target earthquakes from the standard and surface wave catalogs of the Global Centroid Moment Tensor Project and identify 86 shallow 4.3 ≤ MW ≤ 5.8 vertical-CLVD earthquakes located near volcanoes that have erupted in the last ~100 years. The majority of vertical-CLVD earthquakes occur in subduction zones in association with basaltic-to-andesitic stratovolcanoes or submarine volcanoes, although vertical-CLVD earthquakes are also located in continental rifts and in regions of hot spot volcanism. Vertical-CLVD earthquakes are associated with many types of confirmed or suspected eruptive activity at nearby volcanoes, including volcanic earthquake swarms as well as effusive and explosive eruptions and caldera collapse. Approximately 70% of all vertical-CLVD earthquakes studied occur during episodes of documented volcanic unrest at a nearby volcano. Given that volcanic unrest is underreported, most shallow vertical-CLVD earthquakes near active volcanoes are likely related to magma migration or eruption processes. Vertical-CLVD earthquakes with dominant vertical pressure axes generally occur after volcanic eruptions, whereas vertical-CLVD earthquakes with dominant vertical tension axes generally occur before the start of volcanic unrest. The occurrence of these events may be useful for identifying volcanoes that have recently erupted and those that are likely to erupt in the future