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 km2. 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.

This paper proposes a coastal erosion monitoring system for beach erosion management, which we demonstrate on natural and artificial pocket gravel beaches in Croatia. The approach uses low-cost Structure-from-Motion (SfM) photogrammetric imaging and multi-view stereo (MVS) to produce high-resolution 3D beach models for detecting morphological changes and erosion occurrence.Coastal state indicators, such as the shoreline position and subaerial beach volume, are derived from the 3D models and used to quantify changes between surveys. The method is illustrated through two case studies and, to our knowledge, these are the first repetitive measurements taken on the Croatian eastern Adriatic Coast (CEAC). In case of the natural Brseč beach, beach rotation was found to be a response to natural forcing from waves of various incident directions. For the artificial Dugi Rat beach, which loses sediment every winter and is subsequently re-nourished every spring, monitoring showed that beach nourishment is of limited durability. Both case studies showed that the SfM-MVS technique is suitable for the rapid and frequent acquisition of 3D survey data, from which quantitative coastal indicators can be derived to inform future coastal management interventions. Significantly, this low-cost data acquisition has a great potential for regular beach management survey. The introduction of beach monitoring in Croatia is timely because emerging Integrated CoastalZone Management (ICZM) practices will require data-based approaches. Moreover, rare natural pocket beaches and the ever-increasing number of artificial beaches are extremely vulnerable to natural and man-made changes. Adaptive beach management, based on systematic monitoring data, should be included in the ICZM, and we detail how SfM-MVS-based monitoring can be used at different levels of the ICZM. Implementing robust ICZM monitoring will require broad considerations and consultation with all stakeholders, so we propose that SfM-MVS beach surveys should be initially integrated into the existing monitoring practices for CEAC sea water bathing quality. Extension of the existing database with rapidly-gathered low-cost 3D beach survey data, from a number of targeted beaches, could be used to provide a crucial baseline for the ICZM and strategic coastal monitoring of the CEAC.

Summary On 29 December 2020, a shallow earthquake of magnitude Mw 6.4 struck northern Croatia, near the town of Petrinja, more than 24 hours after a strong foreshock (Ml 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-kilometer en échelon right stepping sections, along a NW-SE striking fault that we call the Petrinja-Pokupsko Fault (PPKF). 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 epicenter. 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.

Recent development of lightweight and small size multi-frequency GNSS receivers allows determination of the precise position of the moving platform and spatial data acquisition without the need for setting up and measuring of ground control points. The main advantage of this approach is a higher operational capacity with reduced time and cost of field measurement. This relates to fieldwork in inaccessible areas with demanding terrain configuration. In this paper development and use of a UAS with direct georeferencing of camera sensor for spatial data acquisition is described, and the possibility of 3D scene reconstruction based on the precise position of the camera with predetermined interior parameters is examined. Modern computer vision-based SfM photogrammetry algorithms are used for determining attitude parameters and reconstruction of the scene. For that purpose, several tests on two different test fields were performed using various system parameters for collecting and analysis of several spatial data sets. The presented results demonstrate a satisfactory accuracy (3.1 cm planar and 6.4 cm spatial) of the system for various applications in geodesy.

Terrestrial laser scanning (TLS) in combination with Unmanned Aircraft System (UAS) and modern computer based photogrammetry is currently the best approach for the acquisition of high-resolution 3D spatial information. Highly realistic 3D spatial data sets are becoming the basis for detailed geological studies, providing a multidisciplinary approach in the study and research of both underground and above ground sites. To emphasize the variety of possible implementations of these state-of-the-art methodologies, four characteristic and yet quite different case studies are presented where such geodetic techniques are successfully employed. The presented case studies demonstrate that TLS and UAS photogrammetry, as non-contact surveying methods, are able to reduce survey time and total project costs. As added value, they provide high-resolution data that can be analyzed in a virtual environment from a sedimentological or structural aspect. Stored digital documentation also allows future multi-temporal spatial data comparison at any timeframe and scale, thus enhancing any target geological data gathering and analyses at the studied sites.