On 28 December 2020, seismic activity in the wider Petrinja area strongly intensified after a period of relative seismological quiescence that had lasted more than 100 years (since the well-known M5.8 Kupa Valley earthquake of 1909, which is known based on the discovery of the Mohorovičić discontinuity). The day after the M5 foreshock, a destructive M6.2 mainshock occurred. Outcomes of preliminary seismological, geological and SAR image analyses indicate that the foreshocks, mainshock and aftershocks were generated due to the (re)activation of a complex fault system—the intersection of longitudinal NW–SE right-lateral and transverse NE–SW left-lateral faults along the transitional contact zone of the Dinarides and the Pannonian Basin. According to a survey of damage to buildings, approximately 15% of buildings were very heavily damaged or collapsed. Buildings of special or outstanding historical or cultural heritage significance mostly collapsed or became unserviceable. A preliminary analysis of the earthquake ground motion showed that in the epicentral area, the estimated peak ground acceleration PGA values for the bedrock ranged from 0.29 to 0.44 g. In the close Petrinja epicentral area that is characterized by the superficial deposits, significant ground failures were reported within local site effects. Based on that finding and building damage, we assume that the resulting peak ground acceleration (PGA_site) values were likely between 0.4 and 0.6 g depending on the local site characteristics and the distance from the epicentre.

On 22 March 2020, Zagreb was struck by an M5.5 earthquake that had been expected for more than 100 years and revealed all the failures in the construction of residential buildings in the Croatian capital, especially those built in the first half of the 20th century. Because of that, extensive seismological, geological, geodetic and structural engineering surveys were conducted immediately after the main shock. This study provides descriptions of damage, specifying the building performances and their correlation with the local soil characteristics, i.e., seismic motion amplification. Co-seismic vertical ground displacement was estimated, and the most affected area is identified according to Sentinel-1 interferometric wide-swath data. Finally, preliminary 3D structural modeling of the earthquake sequence was performed, and two major faults were modeled using inverse distance weight (IDW) interpolation of the grouped hypocenters. The first-order assessment of seismic amplification (due to site conditions) in the Zagreb area for the M5.5 earthquake shows that ground motions of approximately 0.16–0.19 g were amplified at least twice. The observed co-seismic deformation (based on Sentinel-1A IW SLC images) implies an approximately 3 cm uplift of the epicentral area that covers approximately 20 km². Based on the preliminary spatial and temporal analyses of the Zagreb 2020 earthquake sequence, the main shock and the first aftershocks evidently occurred in the subsurface of the Medvednica Mountains along a deep-seated southeast-dipping thrust fault, recognized as the primary (master) fault. The co-seismic rupture propagated along the thrust towards northwest during the first half-hour of the earthquake sequence, which can be clearly seen from the time-lapse visualization. The preliminary results strongly support one of the debated models of the active tectonic setting of the Medvednica Mountains and will contribute to a better assessment of the seismic hazard for the wider Zagreb area.

Boulder detachment from the seafloor and subsequent transport and accumulation along rocky coasts is a complex geomorphological process that requires a deep understanding of submarine and onshore environments. This process is especially interesting in semi-enclosed shallow basins characterized by extreme storms, but without a significant tsunami record. Moreover, the response of boulder deposits located close to the coast to severe storms remains, in terms of accurate displacement measurement, limited due to the need to acquire long-term data such as ongoing monitoring datasets and repeated field surveys. We present a multidisciplinary study that includes inland and submarine surveys carried out to monitor and accurately quantify the recent displacement of coastal boulders accumulated on the southernmost coast of the Premantura (Kamenjak) Promontory (Croatia, northern Adriatic Sea). We identified recent boulder movements using unmanned aerial vehicle digital photogrammetry (UAV-DP). Fourteen boulders were moved by the waves generated by a severe storm, named Vaia, which occurred on 29 October 2018. This storm struck Northeast Italy and the Istrian coasts with its full force. We have reproduced the storm-generated waves using unstructured wave model Simulating WAves Nearshore (SWAN), with a significant wave height of 6.2 m in front of the boulder deposit area. These simulated waves are considered to have a return period of 20 to 30 years. In addition to the aerial survey, an underwater photogrammetric survey was carried out in order to create a three-dimensional (3D) model of the seabed and identify the submarine landforms associated with boulder detachment. The survey highlighted that most of the holes can be considered potholes, while only one detachment shape was identified. The latter is not related to storm Vaia, but to a previous storm. Two boulders are lying on the seabed and the underwater surveys highlighted that these boulders may be beached during future storms. Thus, this is an interesting example of active erosion of the rocky coast in a geologically, geomorphologically, and oceanologically predisposed locality.

Seismically induced ground motion at a site is generally influenced by the seismic source, the propagation path and the local site conditions. Over the last several decades, researchers have consistently asserted that for near-site attenuation, the spectral parameter κ is subject primarily to the site conditions. In this research, we estimated the parameter κ based on the acceleration amplitude spectrum of shear waves from local earthquakes recorded by seismological stations situated in the western part of Croatia from the slope of the high-frequency part. The spatial distribution of κ values is comparable with seismological, geophysical and geological features, with the published coda-Q values for each station as well as with the isoseismal maps for selected stronger earthquakes in the study area. The complex pattern of longitudinal and transversal major late-orogenic fault zones dissecting early-orogenic thin-skinned tectonic cover in the Kvarner area and the shallow depth to the Moho (Mohorovičić discontinuity) in the Adriatic foreland (southern Istria) are probably responsible for a significant part of wave attenuation and for the anisotropy of attenuation. Regional near-surface attenuation distribution and modelled macroseismic fields point to the conclusion that attenuation properties of rocks in the northwestern External Dinarides are far from isotropic, and the most likely anisotropy sources are the preferential orientations of cracks and fractures under the local tectonic stress field, trapping of waves along major faults (waveguides), and/or attenuation within the fault zones. These results are important for gaining further insight into the attenuation of near-surface crust layers in the northwestern External Dinarides and the associated Adriatic foreland as well as in similar geotectonic settings.