Seizmičnost Hrvatske u razdoblju 2006–2016

During the ten-year period from 2006 to 2015 a total of 36 733 earthquakes were located in Croatia and its surrounding areas, with 37 main events registering magnitudes from 4.0 to 4.9. Seismically the most active was the coastal part of Croatia confined to two seismically distinguished areas. The NW domain was seismically less active, with almost 10000 located events (seven were of magnitude ML ≥ 4.0), among which were the three strongest events that occurred in Croatia during the observed period. Two of them occurred in the Senj picentral area, the first one on 5 February 2007 at 8:30 UTC (ML = 4.9, Imax = VII °MSK) and the second one on 30 July 2013 at 12:58 UTC, (ML = 4.8, Imax = VI °MSK). The third event occurred near Kornati Islands on 18 July 2007 at 10:54 UTC (ML = 4.8). The SE domain experienced the highest number of earthquakes (over 19 000 located events, with 24 events of magnitude ML ≥ 4.0, among which the strongest one was of magnitude ML = 4.9 with the epicentre in Bosnia and Herzegovina near the Croatian border). The seismicity in the continental part of Croatia was weak-to-moderate, with earthquakes of magnitudes ML ≤ 4.1. Focal mechanisms were obtained for 31 earthquakes with magnitudes ML ≥ 4.0, and individual earthquakes have also been macroseismically analysed. Low current moment release rates for both regions (continental and coastal) as compared to long-term averages, indicate the regions are currently in the strain accumulation phase.Tijekom desetogodišnjeg razdoblja od 2006. do 2015. godine u Hrvatskoj i okolnim područjima locirano je 36 733 potresa, od toga 37 glavnih potresa s magnitudama u rasponu od 4,0 do 4,9. Seizmički najaktivniji je bio priobalni dio Hrvatske karakteriziran s dva različita područja seizmičke aktivnosti. U sjeverozapadnom priobalnom odručju, koje je u odnosu na ostatak obalnog područja bilo slabije seizmički aktivno, locirano je gotovo 10000 potresa, od čega je sedam potresa bilo magnitude ML ≥ 4,0, a među njima su i tri najjača potresa koja su se dogodila u Hrvatskoj tijekom promatranog razdoblja. Dva potresa imala su epicentar u Senjskom epicentralnom području. Prvi potres se dogodio 5. veljače 2007. u 8:30 UTC (ML = 4,9, Imax = VII °MSK), a drugi 30. srpnja 2013. u 12:58 UTC, (ML = 4,8, Imax = VI °MSK). Treći potres se dogodio u Kornatskom arhipelagu 18. srpnja 2007. u 10:54 UTC (ML = 4,8). U jugoistočnom obalnom području dogodio se značajno najveći broj zabilježenih potresa (više od 19000 lociranih potresa, od čega 24 potresa magnitude ML ≥ 4.0, među kojima je najjači bio magnitude ML= 4,9 s epicentrom u Bosni i Hercegovini u blizini hrvatske granice). Seizmičnost u kontinentalnom dijelu Hrvatske bila je slaba do umjerena s potresima magnitude ML ≤ 4,1. Žarišni mehanizmi potresa izračunati su za 31 potres s magnitudama ML ≥ 4,0. Potresi su makroseizmički analizirani. U usporedbi s dugoročnim prosjekom, trenutna relativno niska seizmička aktivnost mjerena oslobođenim seizmičkim momentom u jedinici vremena ukazuje da se oba dijela Hrvatske nalaze se u fazi akumuliranja tektonskih deformacija

Magnus expansion for a chirped quantum two-level system

We derive a Magnus expansion for a frequency chirped quantum two-level system. We obtain a time-independent effective Hamiltonian which generates a stroboscopic time evolution. At lowest order the according dynamics is identical to results from using a rotating wave approximation. We determine, furthermore, also the next higher order corrections within our expansion scheme in correspondence to the Bloch-Siegert shifts for harmonically driven systems. Importantly, our scheme can be extended to more complicated systems, i.e. even many-body systems.Comment: 4 pages, 1 figur

Crustal Thinning From Orogen to Back-Arc Basin: The Structure of the Pannonian Basin Region Revealed by P-to-S Converted Seismic Waves

We present the results of P-to-S receiver function analysis to improve the 3D image of the sedimentary layer, the upper crust, and lower crust in the Pannonian Basin area. The Pannonian Basin hosts deep sedimentary depocentres superimposed on a complex basement structure and it is surrounded by mountain belts. We processed waveforms from 221 three-component broadband seismological stations. As a result of the dense station coverage, we were able to achieve so far unprecedented spatial resolution in determining the velocity structure of the crust. We applied a three-fold quality control process; the first two being applied to the observed waveforms and the third to the calculated radial receiver functions. This work is the first comprehensive receiver function study of the entire region. To prepare the inversions, we performed station-wise H-Vp/Vs grid search, as well as Common Conversion Point migration. Our main focus was then the S-wave velocity structure of the area, which we determined by the Neighborhood Algorithm inversion method at each station, where data were sub-divided into back-azimuthal bundles based on similar Ps delay times. The 1D, nonlinear inversions provided the depth of the discontinuities, shear-wave velocities and Vp/Vs ratios of each layer per bundle, and we calculated uncertainty values for each of these parameters. We then developed a 3D interpolation method based on natural neighbor interpolation to obtain the 3D crustal structure from the local inversion results. We present the sedimentary thickness map, the first Conrad depth map and an improved, detailed Moho map, as well as the first upper and lower crustal thickness maps obtained from receiver function analysis. The velocity jump across the Conrad discontinuity is estimated at less than 0.2 km/s over most of the investigated area. We also compare the new Moho map from our approach to simple grid search results and prior knowledge from other techniques. Our Moho depth map presents local variations in the investigated area: the crust-mantle boundary is at 20–26 km beneath the sedimentary basins, while it is situated deeper below the Apuseni Mountains, Transdanubian and North Hungarian Ranges (28–33 km), and it is the deepest beneath the Eastern Alps and the Southern Carpathians (40–45 km). These values reflect well the Neogene evolution of the region, such as crustal thinning of the Pannonian Basin and orogenic thickening in the neighboring mountain belts

Shear-wave velocity structure beneath the Dinarides from the inversion of Rayleigh-wave dispersion

Highlights • Rayleigh-wave phase velocity in the wider Dinarides region using the two-station method. • Uppermost mantle shear-wave velocity model of the Dinarides-Adriatic Sea region. • Velocity model reveals a robust high-velocity anomaly present under the whole Dinarides. • High-velocity anomaly reaches depth of 160 km in the northern Dinarides to more than 200 km under southern Dinarides. • New structural model incorporating delamination as one of the processes controlling the continental collision in the Dinarides. The interaction between the Adriatic microplate (Adria) and Eurasia is the main driving factor in the central Mediterranean tectonics. Their interplay has shaped the geodynamics of the whole region and formed several mountain belts including Alps, Dinarides and Apennines. Among these, Dinarides are the least investigated and little is known about the underlying geodynamic processes. There are numerous open questions about the current state of interaction between Adria and Eurasia under the Dinaric domain. One of the most interesting is the nature of lithospheric underthrusting of Adriatic plate, e.g. length of the slab or varying slab disposition along the orogen. Previous investigations have found a low-velocity zone in the uppermost mantle under the northern-central Dinarides which was interpreted as a slab gap. Conversely, several newer studies have indicated the presence of the continuous slab under the Dinarides with no trace of the low velocity zone. Thus, to investigate the Dinaric mantle structure further, we use regional-to-teleseismic surface-wave records from 98 seismic stations in the wider Dinarides region to create a 3D shear-wave velocity model. More precisely, a two-station method is used to extract Rayleigh-wave phase velocity while tomography and 1D inversion of the phase velocity are employed to map the depth dependent shear-wave velocity. Resulting velocity model reveals a robust high-velocity anomaly present under the whole Dinarides, reaching the depths of 160 km in the north to more than 200 km under southern Dinarides. These results do not agree with most of the previous investigations and show continuous underthrusting of the Adriatic lithosphere under Europe along the whole Dinaric region. The geometry of the down-going slab varies from the deeper slab in the north and south to the shallower underthrusting in the center. On-top of both north and south slabs there is a low-velocity wedge indicating lithospheric delamination which could explain the 200 km deep high-velocity body existing under the southern Dinarides

Investigation of the central Adriatic lithosphere structure with the AlpArray-CASE seismic experiment

The tectonics of the Adriatic microplate is not well constrained and remains controversial, especially at its contact with the Dinarides, where it acts as the lower plate. While the northern part of the Adriatic microplate will be accu- rately imaged within the AlpArray project, its central and southern parts de- serve detailed studies to obtain a complete picture of its structure and evolution. We set up the Central Adriatic Seismic Experiment (CASE) as a AlpArray Complementary Experiment with a temporary seismic network to provide high- quality seismological data as a foundation for research with state-of-the-art methods and high-precision seismic images of the controversial area. The inter- national AlpArray-CASE project involves four institutions: the Department of Earth Sciences and the Swiss Seismological Service of ETH Zürich (CH), the Department of Geophysics of the Faculty of Science at the University of Zagreb (HR), the Republic Hydrometeorological Service of the Republic of Srpska (BIH) and Istituto Nazionale di Geofisica e Vulcanologia (I). The established temporary seismic network will be operational for at least 18 months. It combines existing permanent and temporary seismic stations operated by the involved institutions together with newly deployed temporary seismic stations, installed in November and December 2016, managed by ETH Zürich and INGV: five in Croatia, four in Bosnia and Herzegovina and one in Italy. We present our scientific aims and network geometry as well as the newly deployed stations sites and settings. In particular, the new stations show favourable noise level (power spectral density estimates). The new network improves considerably the theoretical ray coverage for ambient noise tomography and the magnitude threshold shown in the Bayesian magnitude of completeness threshold map

Seismicity of Croatia in the period 2006-2015

During the ten-year period from 2006 to 2015 a total of 36 733 earthquakes were located in Croatia and its surrounding areas, with 37 main events registering magnitudes from 4.0 to 4.9. Seismically the most active was the coastal part of Croatia confined to two seismically distinguished areas. The NW domain was seismically less active, with almost 10000 located events (seven were of magnitude ML ≥ 4.0), among which were the three strongest events that occurred in Croatia during the observed period. Two of them occurred in the Senj picentral area, the first one on 5 February 2007 at 8:30 UTC (ML = 4.9, Imax = VII °MSK) and the second one on 30 July 2013 at 12:58 UTC, (ML = 4.8, Imax = VI °MSK). The third event occurred near Kornati Islands on 18 July 2007 at 10:54 UTC (ML = 4.8). The SE domain experienced the highest number of earthquakes (over 19 000 located events, with 24 events of magnitude ML ≥ 4.0, among which the strongest one was of magnitude ML = 4.9 with the epicentre in Bosnia and Herzegovina near the Croatian border). The seismicity in the continental part of Croatia was weak-to-moderate, with earthquakes of magnitudes ML ≤ 4.1. Focal mechanisms were obtained for 31 earthquakes with magnitudes ML ≥ 4.0, and individual earthquakes have also been macroseismically analysed. Low current moment release rates for both regions (continental and coastal) as compared to long-term averages, indicate the regions are currently in the strain accumulation phase.Tijekom desetogodišnjeg razdoblja od 2006. do 2015. godine u Hrvatskoj i okolnim područjima locirano je 36 733 potresa, od toga 37 glavnih potresa s magnitudama u rasponu od 4,0 do 4,9. Seizmički najaktivniji je bio priobalni dio Hrvatske karakteriziran s dva različita područja seizmičke aktivnosti. U sjeverozapadnom priobalnom odručju, koje je u odnosu na ostatak obalnog područja bilo slabije seizmički aktivno, locirano je gotovo 10000 potresa, od čega je sedam potresa bilo magnitude ML ≥ 4,0, a među njima su i tri najjača potresa koja su se dogodila u Hrvatskoj tijekom promatranog razdoblja. Dva potresa imala su epicentar u Senjskom epicentralnom području. Prvi potres se dogodio 5. veljače 2007. u 8:30 UTC (ML = 4,9, Imax = VII °MSK), a drugi 30. srpnja 2013. u 12:58 UTC, (ML = 4,8, Imax = VI °MSK). Treći potres se dogodio u Kornatskom arhipelagu 18. srpnja 2007. u 10:54 UTC (ML = 4,8). U jugoistočnom obalnom području dogodio se značajno najveći broj zabilježenih potresa (više od 19000 lociranih potresa, od čega 24 potresa magnitude ML ≥ 4.0, među kojima je najjači bio magnitude ML= 4,9 s epicentrom u Bosni i Hercegovini u blizini hrvatske granice). Seizmičnost u kontinentalnom dijelu Hrvatske bila je slaba do umjerena s potresima magnitude ML ≤ 4,1. Žarišni mehanizmi potresa izračunati su za 31 potres s magnitudama ML ≥ 4,0. Potresi su makroseizmički analizirani. U usporedbi s dugoročnim prosjekom, trenutna relativno niska seizmička aktivnost mjerena oslobođenim seizmičkim momentom u jedinici vremena ukazuje da se oba dijela Hrvatske nalaze se u fazi akumuliranja tektonskih deformacija

Damage Evaluation and Seismic Assessment of a Typical Historical Unreinforced Masonry Building in the Zagreb 2020 Earthquake: A Case Study—Part I

The city of Zagreb, the national capital and economic hub of Croatia, is situated in a seismically active region and hosts a significant array of historical buildings, from the medieval to Austro-Hungarian periods. These buildings possess varying but generally high degrees of vulnerability to seismic loading. This was highlighted in the Zagreb earthquake of 22 March 2020, emphasizing the need for seismic retrofitting in order to preserve this architectural heritage. In this paper, the seismic capacity of one such unreinforced masonry building is considered through a number of analysis methods, including response spectrum, pushover, and out-of-plane wall failure analyses. Given the advantages and disadvantages of the individual methods, their applicability and value in a seismic analysis is considered. Ambient vibration measurements before and after the Zagreb 2020 earthquake, used for model calibration, are also presented. Conclusions are drawn from each individual analysis and later compared. In conclusion, no single analysis method considers all relevant failure modes, and a combination of nonlinear static or dynamic analysis and out-of-plane analysis is recommended. Due to the large volume of the material, it is published in two parts, with ground motion record selection, dynamic analysis, and a comparison of the results published in part two

Transversely isotropic lower crust of Variscan central Europe imaged by ambient noise tomography of the Bohemian Massif

The recent development of ambient noise tomography, in combination with the increasing number of permanent seismic stations and dense networks of temporary stations operated during passive seismic experiments, provides a unique opportunity to build the first high-resolution 3-D shear wave velocity (vS) model of the entire crust of the Bohemian Massif (BM). This paper provides a regional-scale model of velocity distribution in the BM crust. The velocity model with a cell size of 22 km is built using a conventional two-step inversion approach from Rayleigh wave group velocity dispersion curves measured at more than 400 stations. The shear velocities within the upper crust of the BM are ∼0.2 km s−1 higher than those in its surroundings. The highest crustal velocities appear in its southern part, the Moldanubian unit. The Cadomian part of the region has a thinner crust, whereas the crust assembled, or tectonically transformed in the Variscan period, is thicker. The sharp Moho discontinuity preserves traces of its dynamic development expressed in remnants of Variscan subductions imprinted in bands of crustal thickening. A significant feature of the presented model is the velocity-drop interface (VDI) modelled in the lower part of the crust. We explain this feature by the anisotropic fabric of the lower crust, which is characterised as vertical transverse isotropy with the low velocity being the symmetry axis. The VDI is often interrupted around the boundaries of the crustal units, usually above locally increased velocities in the lowermost crust. Due to the north-west–south-east shortening of the crust and the late-Variscan strike-slip movements along the north-east–south-west oriented sutures preserved in the BM lithosphere, the anisotropic fabric of the lower crust was partly or fully erased along the boundaries of original microplates. These weakened zones accompanied by a velocity increase above the Moho (which indicate an emplacement of mantle rocks into the lower crust) can represent channels through which portions of subducted and later molten rocks have percolated upwards providing magma to subsequently form granitoid plutons

High-Resolution Crustal S-wave Velocity Model and Moho Geometry Beneath the Southeastern Alps: New Insights From the SWATH-D Experiment

We compiled a dataset of continuous recordings from the temporary and permanent seismic networks to compute the high-resolution 3D S-wave velocity model of the Southeastern Alps, the western part of the external Dinarides, and the Friuli and Venetian plains through ambient noise tomography. Part of the dataset is recorded by the SWATH-D temporary network and permanent networks in Italy, Austria, Slovenia and Croatia between October 2017 and July 2018. We computed 4050 vertical component cross-correlations to obtain the empirical Rayleigh wave Green’s functions. The dataset is complemented by adopting 1804 high-quality correlograms from other studies. The fast-marching method for 2D surface wave tomography is applied to the phase velocity dispersion curves in the 2–30 s period band. The resulting local dispersion curves are inverted for 1D S-wave velocity profiles using the non-perturbational and perturbational inversion methods. We assembled the 1D S-wave velocity profiles into a pseudo-3D S-wave velocity model from the surface down to 60 km depth. A range of iso-velocities, representing the crystalline basement depth and the crustal thickness, are determined. We found the average depth over the 2.8–3.0 and 4.1–4.3 km/s iso-velocity ranges to be reasonable representations of the crystalline basement and Moho depths, respectively. The basement depth map shows that the shallower crystalline basement beneath the Schio-Vicenza fault highlights the boundary between the deeper Venetian and Friuli plains to the east and the Po-plain to the west. The estimated Moho depth map displays a thickened crust along the boundary between the Friuli plain and the external Dinarides. It also reveals a N-S narrow corridor of crustal thinning to the east of the junction of Giudicarie and Periadriatic lines, which was not reported by other seismic imaging studies. This corridor of shallower Moho is located beneath the surface outcrop of the Permian magmatic rocks and seems to be connected to the continuation of the Permian magmatism to the deep-seated crust. We compared the shallow crustal velocities and the hypocentral location of the earthquakes in the Southern foothills of the Alps. It revealed that the seismicity mainly occurs in the S-wave velocity range between ∼3.1 and ∼3.6 km/s

Evidence for radial anisotropy in the lower crust of the Apennines from Bayesian ambient noise tomography in Europe

Probing seismic anisotropy of the lithosphere provides valuable clues on the fabric of rocks. We present a 3-D probabilistic model of shear wave velocity and radial anisotropy of the crust and uppermost mantle of Europe, focusing on the mountain belts of the Alps and Apennines. The model is built from Love and Rayleigh dispersion curves in the period range 5–149 s. Data are extracted from seismic ambient noise recorded at 1521 broad-band stations, including the AlpArray network. The dispersion curves are first combined in a linearized least squares inversion to obtain 2-D maps of group velocity at each period. Love and Rayleigh maps are then jointly inverted at depth for shear wave velocity and radial anisotropy using a Bayesian Monte Carlo scheme that accounts for the trade-off between radial anisotropy and horizontal layering. The isotropic part of our model is consistent with previous studies. However, our anisotropy maps differ from previous large scale studies that suggested the presence of significant radial anisotropy everywhere in the European crust and shallow upper mantle. We observe instead that radial anisotropy is mostly localized beneath the Apennines while most of the remaining European crust and shallow upper mantle is isotropic. We attribute this difference to trade-offs between radial anisotropy and thin (hectometric) layering in previous studies based on least-squares inversions and long period data (>30 s). In contrast, our approach involves a massive data set of short period measurements and a Bayesian inversion that accounts for thin layering. The positive radial anisotropy (VSH > VSV) observed in the lower crust of the Apennines cannot result from thin layering. We rather attribute it to ductile horizontal flow in response to the recent and present-day extension in the region