ALBERT

Description: Rapid assessment of an earthquake’s impact on the affected society is a crucial step in the early phase of disaster management, navigating the need for further emergency response measures. We demonstrate that felt reports collected via the LastQuake service of the European Mediterranean Seismological Center can be utilized to rapidly estimate the probability of a felt earthquake being high impact rather than low impact on a global scale. Our data-driven, transparent, and reproducible method utilizing Bayes’ theorem and kernel density estimation provides results within 10 min for 393 felt events in 2021. Although a separation of high- and low-impact events remains challenging, the cor- rect and unambiguous assessment of a large portion of low-impact events is a key strength of our approach. We consider our method as an inexpensive addition to the pool of earthquake impact assessment tools, one that is fully independent of seismic data and can be utilized in many populated areas on the planet. Although practical deployment of our method remains an open task, we demonstrate the potential to improve disaster management in regions that currently lack expensive seismic instrumentation.

Description: Accelerograms are the primary source for characterizing strong ground motion. It is therefore of paramount interest to have high‐quality recordings free from any nonphysical contamination. Frequently, accelerograms are affected by baseline jumps and drifts, either related to the instrument and/or a major earthquake. In this work, I propose a correction method for these undesired baseline drifts based on segmented linear least squares. The algorithm operates on the integrated waveforms and combines all three instrument components to estimate a model that modifies the baseline to be at zero continuously. The procedure consists of two steps: first a suite of models with variable numbers of discontinuities is derived for all three instrument components. During this process, the number of discontinuities is reduced in a parsimonious way, for example, two very close discontinuities are merged into a single one. In the second step, the optimal model is selected on the basis of the Bayesian information criterion. I exemplify the application on synthetic waveforms with known discontinuities and on observed waveforms from a unified strong‐motion database of the Japan Meteorological Agency (JMA) and the National Research Institute for Earth Science and Disaster Prevention (NIED, Japan) networks for the major events of the 2016 Kumamoto earthquakes. After the baseline jump correction, the waveforms are furthermore corrected for displacement according to Wang et al. (2011). The resulting displacements are comparable to the Interferometric Synthetic Aperture Radar‐derived displacement estimates for the Kumamoto earthquake sequence.

Description: Ground motion with strong‐velocity pulses can cause significant damage to buildings and structures at certain periods; hence, knowing the period and velocity amplitude of such pulses is critical for earthquake structural engineering. However, the physical factors relating the scaling of pulse periods with magnitude are poorly understood. In this study, we investigate moderate but damaging earthquakes (⁠Mw 6–7) and characterize ground‐motion pulses using the method of Shahi and Baker (2014) while considering the potential static‐offset effects. We confirm that the within‐event variability of the pulses is large. The identified pulses in this study are mostly from strike‐slip‐like earthquakes. We further perform simulations using the frequency–wavenumber algorithm to investigate the causes of the variability of the pulse periods within and between events for moderate strike‐slip earthquakes. We test the effect of fault dips, and the impact of the asperity locations and sizes. The simulations reveal that the asperity properties have a high impact on the pulse periods and amplitudes at nearby stations. Our results emphasize the importance of asperity characteristics, in addition to earthquake magnitudes for the occurrence and properties of pulses produced by the forward directivity effect. We finally quantify and discuss within‐ and between‐event variabilities of pulse properties at short distances.

Description: The steady increase of ground-motion data not only allows new possibilities but also comes with new challenges in the development of ground-motion models (GMMs). Data classification techniques (e.g., cluster analysis) do not only produce deterministic classifications but also probabilistic classifications (e.g., probabilities for each datum to belong to a given class or cluster). One challenge is the integration of such continuous classification in regressions for GMM development such as the widely used mixed-effects model. We address this issue by introducing an extension of the mixed-effects model to incorporate data weighting. The parameter estimation of the mixed-effects model, that is, fixed-effects coefficients of the GMMs and the random-effects variances, are based on the weighted likelihood function, which also provides analytic uncertainty estimates. The data weighting permits for earthquake classification beyond the classical, expert-driven, binary classification based, for example, on event depth, distance to trench, style of faulting, and fault dip angle. We apply Angular Classification with Expectation–maximization, an algorithm to identify clusters of nodal planes from focal mechanisms to differentiate between, for example, interface- and intraslab-type events. Classification is continuous, that is, no event belongs completely to one class, which is taken into account in the ground-motionmodeling. The theoretical framework described in this article allows for a fully automatic calibration of ground-motionmodels using large databases with automated classification and processing of earthquake and ground-motion data. As an example, we developed a GMM on the basis of the GMM by Montalva et al. (2017) with data from the strong-motion flat file of Bastías and Montalva (2016) with ∼2400 records from 319 events in the Chilean subduction zone. Our GMMwith the data-driven classification is comparable to the expert-classification-based model. Furthermore, the model shows temporal variations of the between-event residuals before and after large earthquakes in the region.

Description: Surface heat flow is a geophysical variable that is affected by a complex combination of various heat generation and transport processes. The processes act on different lengths scales, from tens of meters to hundreds of kilometers. In general, it is not possible to resolve all processes due to a lack of data or modeling resources, and hence the heat flow data within a region is subject to residual fluctuations. We introduce the REgional HEAT-Flow Uncertainty and aNomaly Quantification (REHEATFUNQ) model, version 2.0.1. At its core, REHEATFUNQ uses a stochastic model for heat flow within a region, considering the aggregate heat flow to be generated by a gamma-distributed random variable. Based on this assumption, REHEATFUNQ uses Bayesian inference to (i) quantify the regional aggregate heat flow distribution (RAHFD) and (ii) estimate the strength of a given heat flow anomaly, for instance as generated by a tectonically active fault. The inference uses a prior distribution conjugate to the gamma distribution for the RAHFDs, and we compute parameters for a uninformed prior distribution from the global heat flow database by Lucazeau (2019). Through the Bayesian inference, our model is the first of its kind to consistently account for the variability in regional heat flow in the inference of spatial signals in heat flow data. Interpretation of these spatial signals and in particular their interpretation in terms of fault characteristics (particularly fault strength) form a long-standing debate within the geophysical community. We describe the components of REHEATFUNQ and perform a series of goodness-of-fit tests and synthetic resilience analyses of the model. While our analysis reveals to some degree a misfit of our idealized empirical model with real-world heat flow, it simultaneously confirms the robustness of REHEATFUNQ to these model simplifications. We conclude with an application of REHEATFUNQ to the San Andreas fault in California. Our analysis finds heat flow data in the Mojave section to be sufficient for an analysis and concludes that stochastic variability can allow for a surprisingly large fault-generated heat flow anomaly to be compatible with the data. This indicates that heat flow alone may not be a suitable quantity to address fault strength of the San Andreas fault.