Search

The determination of earthquake source parameters is an important task in seismology. For many applications, it is also valuable to understand the uncertainties associated with these determinations, and this is particularly true in the context of earthquake early warning (EEW) and hazard mitigation. In this paper, we develop a framework for probabilistic moment tensor point source inversions in near real time. Our methodology allows us to find an approximation to p(m|d), the conditional probability of source models (m) given observations (d). This is obtained by smoothly interpolating a set of random prior samples, using Mixture Density Networks (MDNs)—a class of neural networks which output the parameters of a Gaussian mixture model. By combining multiple networks as ‘committees’, we are able to obtain a significant improvement in performance over that of a single MDN. Once a committee has been constructed, new observations can be inverted within milliseconds on a standard desktop computer. The method is therefore well suited for use in situations such as EEW, where inversions must be performed routinely and rapidly for a fixed station geometry. To demonstrate the method, we invert regional static GPS displacement data for the 2010 MW 7.2 El Mayor Cucapah earthquake in Baja California to obtain estimates of magnitude, centroid location and depth and focal mechanism. We investigate the extent to which we can constrain moment tensor point sources with static displacement observations under realistic conditions. Our inversion results agree well with published point source solutions for this event, once the uncertainty bounds of each are taken into account. [ABSTRACT FROM AUTHOR]

How do body-wave traveltimes constrain the Earth's radial (1-D) seismic structure? Existing 1-D seismological models underpin 3-D seismic tomography and earthquake location algorithms. It is therefore crucial to assess the quality of such 1-D models, yet quantifying uncertainties in seismological models is challenging and thus often ignored. Ideally, quality assessment should be an integral part of the inverse method. Our aim in this study is twofold: (i) we show how to solve a general Bayesian non-linear inverse problem and quantify model uncertainties, and (ii) we investigate the constraint on spherically symmetric P-wave velocity (VP) structure provided by body-wave traveltimes from the EHB bulletin (phases Pn, P, PP and PKP). Our approach is based on artificial neural networks, which are very common in pattern recognition problems and can be used to approximate an arbitrary function. We use a Mixture Density Network to obtain 1-D marginal posterior probability density functions (pdfs), which provide a quantitative description of our knowledge on the individual Earth parameters. No linearization or model damping is required, which allows us to infer a model which is constrained purely by the data.We present 1-D marginal posterior pdfs for the 22 VP parameters and seven discontinuity depths in our model. P-wave velocities in the inner core, outer core and lower mantle are resolved well, with standard deviations of ∼0.2 to 1 per cent with respect to the mean of the posterior pdfs. The maximum likelihoods of VP are in general similar to the corresponding ak135 values, which lie within one or two standard deviations from the posterior means, thus providing an independent validation of ak135 in this part of the radial model. Conversely, the data contain little or no information on P-wave velocity in the D′′ layer, the upper mantle and the homogeneous crustal layers. Further, the data do not constrain the depth of the discontinuities in our model. Using additional phases available in the ISC bulletin, such as PcP, PKKP and the converted phases SP and ScP, may enhance the resolvability of these parameters. Finally, we show how the method can be extended to obtain a posterior pdf for a multidimensional model space. This enables us to investigate correlations between model parameters. [ABSTRACT FROM AUTHOR]

SUMMARY What makes a seismogram look like a seismogram? Seismic data sets generally contain waveforms sharing some set of visual characteristics and features-indeed, seismologists routinely exploit this when performing quality control 'by hand'. Understanding and harnessing these characteristics offers the prospect of a deeper understanding of seismic waveforms, and opens up many potential new techniques for processing and working with data. In addition, the fact that certain features are shared between waveforms suggests that it may be possible to transform the data away from the time domain, and represent the same information using fewer parameters. If so, this would be a significant step towards making fully non-linear tomographic inversions computationally tractable. Hinton & Salakhutdinov showed that a particular class of neural network, termed 'autoencoder networks', may be used to find lower-dimensional encodings of complex binary data sets. Here, we adapt their work to the continuous case to allow the use of autoencoders for seismic waveforms, and offer a demonstration in which we compress 512-point waveforms to 32-element encodings. We also demonstrate that the mapping from data to encoding space, and its inverse, are well behaved, as required for many applications. Finally, we sketch a number of potential applications of the technique, which we hope will be of practical interest across all seismological disciplines, and beyond. [ABSTRACT FROM AUTHOR]

The ability to handle large amounts of data automatically is essential for any major tomographic inversion. As part of this process, it is necessary to differentiate between high-quality seismograms, and those that are unusable due to noise or other errors. This quality assessment is traditionally made visually; however, the sheer quantity of data in a modern tomographic data set makes this approach unfeasible. It is therefore necessary to develop techniques for automating this quality assessment process. We demonstrate that a simple neural network, trained to recognize the frequency-domain characteristics of high- and low-quality data, can successfully distinguish the two classes in unseen data. We demonstrate that the resulting clean data sets are of sufficient quality to allow full-waveform determination of event focal mechanisms and hypocentral parameters. The process we outline allows the rapid creation of a high-quality data set for seismic tomography. Depending on application, this may be suitable for use without further refinement. In some circumstances, a further visual inspection may remain desirable to ensure the data set is noise-free; however, a significant benefit will still derive from the reduction in number of traces to be examined. This will enable full-waveform inversion using significantly larger data sets than has hitherto been possible. The selection strategy relies only on measurements made from the seismogram, and on rough estimates of hypocentral location—the final data set does not depend on any a priori assumptions regarding earth structure or wave propagation. Our focus has been on data selection for seismic tomography, but the approach is general and may find application across a wide range of seismic investigations. An automated system is of interest wherever large data sets must be handled, or where time is of the essence—such as in earthquake hazard assessment. [ABSTRACT FROM AUTHOR]