SUMMARY: The Eastern-Europe region (EER), is a complex geotectonic area that captures part of the Alpine-Himalayan Orogen, the subduction of multiple NeoTethys Branches and part of the East European Craton. It is one of the most exciting geological areas in Europe due to a diversity of tectonic processes acting within it: extensional basin evolution, oceanic subduction, post-collisional volcanism, as well as active crustal deformation associated with the push of the Adria plate or the pull of the actively detaching Vrancea slab. This makes EER an excellent natural laboratory to study the behaviour of the lithosphere–asthenosphere system in a heterogeneous tectonic setting. To investigate the lateral heterogeneity and physical properties of the crust in the EER, we use ambient seismic noise data recorded by the vertical components of broad-band stations that have been operational between 1999 and 2020 in Eastern Europe and surrounding regions. We used this significant amount of data and the latest processing techniques of the ambient seismic noise field based on the continuous wavelet transform to compute cross-correlations between various station pairs, turning every available seismic station into a virtual source. The coda of the interstation cross-correlograms were used to determine coda quality factors (Qc) of Rayleigh waves in four different period ranges (3.0–5.0, 5.0–10.0, 10.0–20.0 and 20.0–30.0 s) and to invert them in the 2-D space, constructing the highest resolution attenuation tomography of the region. Our results reveal high attenuation features throughout the northeast Pannonian region, the Bohemian Massif, the East Carpathians and the Moesian Platform. Nevertheless, our findings do not emphasize a close correlation between the depth of sedimentary basins and attenuation features identified at longer periods. In addition, Qc variations are larger at short periods, indicating higher heterogeneity in the uppermost crust of Eastern Europe. Our findings demonstrate the higher efficiency of noise correlation approaches relative to earthquake data analyses investigating Qc at low frequencies.

Romanian seismicity is mainly confined to the Eastern Carpathians Arc bend (ECAB), where strong subcrustal earthquakes (magnitude up to 7.9) are generated in a narrow lithospheric body descending into the mantle. The seismic activity in the overlying crust is spread over a larger area, located mostly toward the outer side of the ECAB. It is significantly smaller than subcrustal seismicity, raising controversies about possible upper mantle-crust coupling. A significant earthquake sequence took place in the foreland of the ECAB triggered on 22 November 2014 by a mainshock of magnitude 5.7 (the greatest instrumentally recorded earthquake in this region) located in the lower crust. The mainshock triggered a significant increase in the number of small-magnitude events spread over an unusually large area in the ECAB. The paper’s goal is to compute the source parameters of the earthquakes that occurred during the aforementioned sequence, by empirical application of Green’s function and spectral ratio techniques. Fault plane solutions are determined using multiple methods and seismicity evolution at regional scale is investigated. Our results highlight a still active deformation regime at the edge of the EE Craton, while the source parameters reveal a complex fracture of the mainshock and a very high-stress drop.

The Vrancea seismic zone (VSZ), located in Romania, at the sharp bend of the southeast Carpathians, is an anomalous intraplate seismic nest releasing the largest strain in continental Europe. The last powerful earthquake in Vrancea occurred in 1977 causing significant damage to the densely populated cities in southeast Romania. The seismic infrastructure was underdeveloped at that time and the earthquake was not well recorded locally. Using ambient seismic noise recorded at modern broadband seismic stations around Vrancea in the year 2020, and the approximate moment tensor solutions of significant earthquakes Mw>6 since 1977, we reconstruct ground-motion waveforms recorded by modern seismometers decades after their occurrence. The virtual earthquake approach uses empirical Green’s functions from ambient noise cross correlations between pairs of seismic stations and adds the signatures of a model earthquake: double couple mechanism, buried source, and a realistic earth model in the epicentral area. Ambient noise data capture the signatures of complex subsurface structures and the 3D path effects, providing a more physically accurate representation of the seismic wavefield than purely synthetic data. Our new results demonstrate the viability of this innovative method and provide a unique opportunity for more accurate seismic hazard analysis in intraplate seismic zones with insufficient instrumental data.

Seismic velocity is the geophysical property that has a key role in characterizing dynamic processes and the state of the stress around the faults, providing valuable information regarding the change in the tectonic regime. The stress in the crust is an important indicator of the possible occurrence of a major earthquake, and the variation of seismic velocities, in time, can provide a clearer picture on the tectonic processes taking place in the region. In the crust, velocities change before, during, and after earthquakes through several mechanisms related to fault deformations, pore pressure, stress changes, and recovery processes. In this study, we investigate the possible correlation between the changes of seismic velocities (Vp/Vs) in time and the occurrence of moderate size crustal and intermediate depth earthquakes from the Vrancea region. Our findings show that there are no significant variations in Vp/Vs for the intermediate depth earthquakes, while crustal events have decreased seismic activity prior to the main earthquake and no high Vp/Vs anomalies. Our results indicate key aspects, and such analyses should be carried out in real-time to continuously explore any unusual pattern pointed out by the seismic velocity changes. Vp/Vs and their standard errors can also be used to describe seismic activity patterns that shape the tectonic evolution of the area.