Seismic moment tensors from synthetic rotational and translational ground motion: Green's functions in 1-D versus 3-D

Seismic moment tensors are an important tool and input variable for many studies in the geosciences. The theory behind the determination of moment tensors is well established. They are routinely and (semi-) automatically calculated on a global scale. However, on regional and local scales, there are still several difficulties hampering the reliable retrieval of the full seismic moment tensor. In an earlier study, we showed that the waveform inversion for seismic moment tensors can benefit significantly when incorporating rotational ground motion in addition to the commonly used translational ground motion. In this study, we test, what is the best processing strategy with respect to the resolvability of the seismic moment tensor components: inverting three-component data with Green’s functions (GFs) based on a 3-D structural model, six-component data with GFs based on a 1-D model, or unleashing the full force of six-component data and GFs based on a 3-D model? As a reference case, we use the inversion based on three-component data and 1-D structure, which has been the most common practice in waveform inversion for moment tensors so far. Building on the same Bayesian approach as in our previous study, we invert synthetic waveforms for two test cases from the Korean Peninsula: one is the 2013 nuclear test of the Democratic People’s Republic of Korea and the other is an Mw  5.4 tectonic event of 2016 in the Republic of Korea using waveform data recorded on stations in Korea, China and Japan. For the Korean Peninsula, a very detailed 3-D velocity model is available. We show that for the tectonic event both, the 3-D structural model and the rotational ground motion, contribute strongly to the improved resolution of the seismic moment tensor. The higher the frequencies used for inversion, the higher is the influence of rotational ground motions. This is an important effect to consider when inverting waveforms from smaller magnitude events. The explosive source benefits more from the 3-D structural model than from the rotational ground motion. Nevertheless, the rotational ground motion can help to better constraint the isotropic part of the source in the higher frequency range.SD thanks Bayrische Forschungsallianz (BayFOR) for funding her internship at the Australian National University (ANU, grant: BayIntAn-LMU-2017-66). Further funding comes from the European Research Council (advanced grant to HI: ROMY, number: 339991)

Global quieting of high-frequency seismic noise due to COVID-19 pandemic lockdown measures

Human activity causes vibrations that propagate into the ground as high-frequency seismic waves. Measures to mitigate the COVID-19 pandemic caused widespread changes in human activity, leading to a months-long reduction in seismic noise of up to 50%. The 2020 seismic noise quiet period is the longest and most prominent global anthropogenic seismic noise reduction on record. While the reduction is strongest at surface seismometers in populated areas, this seismic quiescence extends for many kilometers radially and hundreds of meters in depth. This provides an opportunity to detect subtle signals from subsurface seismic sources that would have been concealed in noisier times and to benchmark sources of anthropogenic noise. A strong correlation between seismic noise and independent measurements of human mobility suggests that seismology provides an absolute, real-time estimate of population dynamics

Large expert-curated database for benchmarking document similarity detection in biomedical literature search

Document recommendation systems for locating relevant literature have mostly relied on methods developed a decade ago. This is largely due to the lack of a large offline gold-standard benchmark of relevant documents that cover a variety of research fields such that newly developed literature search techniques can be compared, improved and translated into practice. To overcome this bottleneck, we have established the RElevant LIterature SearcH consortium consisting of more than 1500 scientists from 84 countries, who have collectively annotated the relevance of over 180 000 PubMed-listed articles with regard to their respective seed (input) article/s. The majority of annotations were contributed by highly experienced, original authors of the seed articles. The collected data cover 76% of all unique PubMed Medical Subject Headings descriptors. No systematic biases were observed across different experience levels, research fields or time spent on annotations. More importantly, annotations of the same document pairs contributed by different scientists were highly concordant. We further show that the three representative baseline methods used to generate recommended articles for evaluation (Okapi Best Matching 25, Term Frequency-Inverse Document Frequency and PubMed Related Articles) had similar overall performances. Additionally, we found that these methods each tend to produce distinct collections of recommended articles, suggesting that a hybrid method may be required to completely capture all relevant articles. The established database server located at https://relishdb.ict.griffith.edu.au is freely available for the downloading of annotation data and the blind testing of new methods. We expect that this benchmark will be useful for stimulating the development of new powerful techniques for title and title/abstract-based search engines for relevant articles in biomedical research.Peer reviewe

Middle East - North Africa and the millennium development goals : implications for German development cooperation

          Closed-loop controlled combustion is a promising technique to improve the overall performance of internal combustion engines and Diesel engines in particular. In order for this technique to be implemented some form of feedback from the combustion process is required. The feedback signal is processed and from it combustionrelated parameters are computed. These parameters are then fed to a control process which drives a series of outputs (e.g. injection timing in Diesel engines) to control their values. This paper’s focus lies on the processing and computation that is needed on the feedback signal before this is ready to be fed to the control process as well as on the electronics necessary to support it. A number of feedback alternatives are briefly discussed and for one of them, the in-cylinder pressure sensor, the CA50 (crank angle in which the integrated heat release curve reaches its 50% value) and the IMEP (Indicated Mean Effective Pressure) are identified as two potential control variables. The hardware architecture of a system capable of calculating both of them on-line is proposed and necessary feasibility size and speed considerations are made by implementing critical blocks in VHDL targeting a flash-based Actel ProASIC3 automotive-grade FPGA

Improved finite-source inversion through joint measurements of rotational and translational ground motions: a numerical study

With the prospects of seismic equipment being able to measure rotational ground motions in a wide frequency and amplitude range in the near future, we engage in the question of how this type of ground motion observation can be used to solve the seismic source inverse problem. In this paper, we focus on the question of whether finite-source inversion can benefit from additional observations of rotational motion. Keeping the overall number of traces constant, we compare observations from a surface seismic network with 44 three-component translational sensors (classic seismometers) with those obtained with 22 six-component sensors (with additional three-component rotational motions). Synthetic seismograms are calculated for known finite-source properties. The corresponding inverse problem is posed in a probabilistic way using the Shannon information content to measure how the observations constrain the seismic source properties. We minimize the influence of the source receiver geometry around the fault by statistically analyzing six-component inversions with a random distribution of receivers. Since our previous results are achieved with a regular spacing of the receivers, we try to answer the question of whether the results are dependent on the spatial distribution of the receivers. The results show that with the six-component subnetworks, kinematic source inversions for source properties (such as rupture velocity, rise time, and slip amplitudes) are not only equally successful (even that would be beneficial because of the substantially reduced logistics installing half the sensors) but also statistically inversions for some source properties are almost always improved. This can be attributed to the fact that the (in particular vertical) gradient information is contained in the additional motion components. We compare these effects for strike-slip and normal-faulting type sources and confirm that the increase in inversion quality for kinematic source parameters is even higher for the normal fault. This indicates that the inversion benefits from the additional information provided by the horizontal rotation rates, i.e., information about the vertical displacement gradient

The Ahar-Varzeghan earthquake doublet (M<inf>w</inf> 6.4 and 6.2) of 11 August 2012: Regional seismic moment tensors and a seismotectonic interpretation

On 11 August 2012 an earthquake doublet (Mw 6.4 and 6.2) occurred near the city of Ahar, northwest Iran. Both events were only 6 km and 11 minutes apart, producing a surface rupture of about 12 km in length. Historical and modern seismicity has so far been sparse in this area. Spatially, the region represents a transitional zone between different tectonic domains, including compression in Iran, westward extrusion of the Anatolian plate, and thrusting beneath the Caucasus. In this study, we inverted the surface waveforms of the two mainshocks and 11 aftershocks (Mw ≥4:3) to obtain regional seismic moment tensors. The earthquakes analyzed can be grouped into pure strike slip (including the first mainshock) and oblique reverse mechanisms (including the second mainshock). The sequence provides information about faulting mechanisms at the spatial scale of the entire rock volume affected by the earthquake doublet, including coinciding deformation on minor faults (sub)parallel to the main fault and Riedel shears. It occurred on a so far unknown fault structure, which we call the Ahar fault. Alongside the seismological data, we used geological maps, satellite images, and digital elevation data to analyze the geomorphology of the region. Our analysis suggests that the adjacent North Tabriz fault, which accomodates up to 7 mm=yr of rightlateral strike-slip faulting, does not compensate the entire lateral shear strain, and that part of it is compensated farther north. Combined, our results suggest a temporally and spatially complex style of deformation (reverse and strike slip) overprinting older reverse deformation