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.

On 29 December 2020, a shallow earthquake of magnitude M_w 6.4 struck northern Croatia, near the town of Petrinja, more than 24 hr after a strong foreshock (M_L 5). We formed a reconnaissance team of European geologists and engineers, from Croatia, Slovenia, France, Italy and Greece, rapidly deployed in the field to map the evidence of coseismic environmental effects. In the epicentral area, we recognized surface deformation, such as tectonic breaks along the earthquake source at the surface, liquefaction features (scattered in the fluvial plains of Kupa, Glina and Sava rivers), and slope failures, both caused by strong motion. Thanks to this concerted, collective and meticulous work, we were able to document and map a clear and unambiguous coseismic surface rupture associated with the main shock. The surface rupture appears discontinuous, consisting of multi-kilometre en échelon right stepping sections, along a NW–SE striking fault that we call the Petrinja-Pokupsko Fault. The observed deformation features, in terms of kinematics and trace alignments, are consistent with slip on a right lateral fault, in agreement with the focal solution of the main shock. We found mole tracks, displacement on faults affecting natural features (e.g. drainage channels), scarplets and more frequently breaks of anthropogenic markers (roads, fences). The surface rupture is observed over a length of ∼13 km from end-to-end, with a maximum displacement of 38 cm, and an average displacement of ∼10 cm. Moreover, the liquefaction extends over an area of nearly 600 km² around the epicentre. Typology of liquefaction features include sand blows, lateral spreading phenomenon along the road and river embankments, as well as sand ejecta of different grain size and matrix. Development of large and long fissures along the fluvial landforms, current or ancient, with massive ejections of sediments is pervasive. These features are sometimes accompanied by small horizontal displacements. Finally, the environmental effects of the earthquake appear to be reasonably consistent with the usual scaling relationships, in particular the surface faulting. This rupture of the ground occurred on or near traces of a fault that shows clear evidence of Quaternary activity. Further and detailed studies will be carried out to characterize this source and related faults in terms of future large earthquakes potential, for their integration into seismic hazard models.

There is growing evidence that many large coastal boulder deposits found on the exposed rocky ocean shores were deposited by extreme storm waves rather than by catastrophic tsunamis, as previously thought. In addition, before the first discovery in the northern Adriatic a few years ago, such deposits were not expected in relatively shallow semi-enclosed inland basins. Here we report on a large coastal coarse-clast deposit on the central Adriatic island of Mana, which also contains numerous large storm boulders that weigh up to several tons. Large solitary boulders are also located outside of the deposit, closer to the partly submerged sea cliff and the wave impact. The erosion of the cliff top and displacement of the carbonate bedrock fragments began when the extreme waves inundated the lowermost part of the cliff edge, probably during the late Holocene sea-level rise. The UAS photogrammetry-based fragmentation analysis of the storm deposit and the calculated fractal dimension value indicate that the material was fragmented by multiple high-energy events. A comparison of the available photographs indicates that displacements of the most exposed solitary boulders probably occurred during Vaia, the last extreme storm that hit the Adriatic on the 29th of October 2018. However, the modeled maximum wave height south of Mana during the peak of the storm would be insufficient to move these boulders. Yet local geomorphology probably further influenced the increase in wave height that, in combination with specific geological features, caused displacements of the boulders. There is a shorter fetch affecting Mana Island with respect to the northern Adriatic boulder field in southern Istria. Thus, such an active local erosion of the generally stable eastern Adriatic karstic coast depends on the extreme storms that have a weaker impact in the central than in the northern Adriatic.

Plitvice Lakes National Park is the largest national park in Croatia and also the oldest from 1949. It was added to the UNESCO World Natural Heritage List in 1979, due to the unique physicochemical and biological conditions that have led to the creation of 16 named and several smaller unnamed lakes, which are cascading one into the next. Previous scientific research proved that the increased amount of dissolved organic matter (pollution) stops the travertine processes on Plitvice Lakes. Therefore, this complex, dynamic but also fragile geological, biological and hydrological system required a comprehensive limnological survey. Thirteen of the sixteen lakes mentioned above were initially surveyed from the air by an unmanned aircraft equipped with a survey grade GNSS and a full frame high-resolution full-screen camera. From these recordings, a georeferenced, high-resolution orthophoto was generated, on which the following surveys by a multibeam sonar depended. It is important to mention that this was the first time that these lakes had ever been surveyed both with the multibeam sonar technique and with such a high-resolution camera. Due to the fact that these thirteen lakes are difficult to reach and often too shallow for a boat-mounted sonar, a special autonomous surface vehicle was developed. The lakes were surveyed by the autonomous surface vehicle mounted with a multibeam sonar to create detailed bathymetric models of the lakes. The missions were planned for the surface vehicle based on the orthophoto from the preliminary studies. A detailed description of the methodology used to survey the different lakes is given here. In addition, the resulting high-resolution bathymetric maps are presented and analysed together with an overview of average, maximum depths and number of data points. Numerous interesting depressions, which are phenomena consistent with previous studies of Plitvice Lakes, are noted at the lake beds and their causes are discussed. This study shows the huge potential of remote sensing technologies integrated into autonomous vehicles in terms of much faster surveys, several orders of magnitude more data points (compared to manual surveys of a few decades ago), as well as data accuracy, precision and georeferencing.

Earthquakes and related coseismic effects at the surface, both primary and secondary, such as liquefaction and lateral spreading, can impact humans due to induced economic or social disruptions (e.g. slope, bridge and building foundation failures, flotation of buried structures). In this respect, it results of primary interest to map liquefaction induced evidences soon after an earthquake. On the 29th December 2020, a major earthquake (Mw 6.4) occurred in Croatia, close to Petrinja, 45 km south of Zagreb, generating widespread liquefaction and lateral spreading phenomena in a radius of approximately 20 km from the epicentre. A European team of researchers (geologists and engineers), in strict collaboration with the Croatian Geological Survey, performed field reconnaissance campaigns with the aim to provide a detailed identification and characterization of the primary and secondary geological and geotechnical coseismic effects induced by the Croatian earthquakes. Specifically with reference to the liquefaction phenomena, the Working Group integrated the data collected directly in the field with those from remote survey by drone aerial photos acquired in the post-event immediate. The adopted process allowed the collection of the liquefaction record with the highest possible completeness both in terms of pattern and distribution of the phenomena. The database includes several detailed case studies typified by the following characteristics: (1) liquefaction occurring on alluvial plain sites (Kupa river, Sava river and Glina river); (2) blows made by sand and/or gravel with local presence of shells and armored mud balls; (3) lateral spreading phenomena along road and river embankments; (4) sand ejecta of different grain size and matrix, even at the same site; (5) sand and/or gravel ejecta along fault traces. The characteristics of these features are discussed with reference to the alluvial setting and tectonic context. All together, the detailed survey of these recent liquefaction features will assist to build new empirical relations, to update the existing ones and to mitigate the effects of future earthquakes recognizing liquefaction prone areas for a correct land use planning, as for seismic microzonation studies.