Detailed multi-beam bathymetry, sub-bottom acoustic, and side-scan sonar observations of submerged canyons with tufa barriers were used to characterize the Zrmanja River karst estuary on the eastern Adriatic coast, Croatia. This unique karst environment consists of two submerged karst basins (Novigrad Sea and Karin Sea) that are connected with river canyons named Novsko Ždrilo and Karinsko Ždrilo. The combined use of high-resolution geophysical data with legacy topography and bathymetry data in a GIS environment allowed for the description and interpretation of this geomorphological setting in relation to the Holocene sea-level rise. The tufa barriers had a predominant influence on the Holocene flooding dynamics of the canyons and karst basins. Here, we describe the possible river pathways from the basins during the lowstand and the formation of a lengthening estuary during the Holocene sea-level rise. Based on the analyzed morphologies and the relative sea-level curve for the Adriatic Sea, the flooding of the Novsko Ždrilo occurred 9200 years before present (BP) and Karinsko Ždrilo was flooded after 8400 years BP. The combination of high-resolution geophysical methods gave an accurate representation of the karst estuarine seafloor and the flooding of semi-isolated basins due to sea-level rise.

The upper part of the Krka River estuary and Prokljan Lake are a specific example of a well-stratified estuarine environment in a submerged river canyon. Here, we reconstructed the geomorphological evolution of the area and classified the data gathered in the study, integrating multibeam echosounder data, backscatter echosounder data, side-scan sonar morpho-bathymetric surveys, and acoustic sub-bottom profiling, with the addition of ground-truthing and sediment analyses. This led to the successful classification of the bottom sediments using the object-based image analysis method. Additional inputs to the multibeam echosounder data improved the segmentation of the seafloor classification, geology, and morphology of the surveyed area. This study uncovered and precisely defined distinct geomorphological features, specifically submerged tufa barriers and carbonate mounds active during the Holocene warm periods, analogous to recent tufa barriers that still exist and grow in the upstream part of the Krka River. Fine-grained sediments, classified as estuarine sediments, hold more organic carbon than coarse-grained sediments sampled on barriers. A good correlation of organic carbon with silt sediments allowed the construction of a prediction map for marine sedimentary carbon in this estuarine/lake environment using multibeam echosounder data. Our findings highlight the importance of additional inputs to multibeam echosounder data to achieve the most accurate results.

Koločep bay is a 30 km elongated karst basin located between the Elafiti Islands and the mainland with a NW–SE general direction. The bay lies within the seismically active southern Dalmatia zone. Irregular grid sub-bottom profiles and two legacy reflection seismic profiles have been used to determine the overall morphology of the bay and to establish the seismic stratigraphy of its sedimentary infill. Three major seismic–stratigraphic units have been identified in the upper ~50 m of the ~120-meter-thick sedimentary sequence that lies above the karstified limestone bedrock. The karst polje basin was flooded due to sea-level rise before 12.1 kyr BP. The morphology of the bay implies complex influences of Late Glacial meltwater discharges, aeolian sand deposition, the existence of paleo–ponor/karst spring zones and buried Pleistocene river channels. The Pleistocene seismotectonic units are deformed in the NW and SE parts of the basin. The central part of the basin has no signatures of intensive tectonic activity during the Holocene. A major erosion event was identified that led to the formation of a basin within the older sedimentary infill. In the southern part of the basin, we have evidence of Holocene tectonic activity with the formation of erosional scarps on the seafloor of the bay.

This study focuses on the analysis of sediment core retrieved from the deepest part (25 m) of Pirovac Bay. A long sedimentary sequence (7.45 m) supplemented by a shorter sediment core (1.45 m) from a shallower part of the bay was analyzed for sedimentological, mineralogical, geochemical, and micropaleontological (ostracod) parameters. The sediment thickness above the underlying karst paleorelief (karstic bedrock) is up to 12 m. Sediments recorded a transition from a freshwater to a marine environment starting from post-Neapolitan Yellow Tuff tephra sedimentation. First, the floodplain developed in Pirovac Bay, with intermittent pools and ponds, followed by wetland environment. The formation of a shallow freshwater paleolake during the Middle Holocene at 10 cal kyr BP was enabled by the rising sea level and high freshwater input from the karstified underground from the adjacent Lake Vrana (Biograd na Moru). The onset of marine intrusions through the karstified underground is evident with formation of a brackish lake in the Pirovac Bay basin. Marine transgression and flooding of the bay occurred at 7.3 cal kyr BP, evidenced by the geochemical and ostracod parameters, providing crucial insights into the dynamics of coastal inundation under past climate change. Intriguingly, freshwater ostracod species were still present in the marine sediments, brought into the bay from Lake Vrana through surficial canal Prosika and groundwater discharge (numerous estavelles) along the northeastern shores of the bay, proving their mutual influence. This submerged Holocene freshwater paleolake, reported here for the first time, underlines the sensitivity of coastal karst systems to the rise in sea level and serves to stress how important understanding of these processes is for effective management in coastal zone and climate change adaptation strategies. The findings provided evidence supporting the existence of coastal marine basins as freshwater lakes prior to being flooded by seawater as a consequence of the Holocene post-glacial sea level rise.

A 7.4 m long sediment core has been retrieved from the central part of Lake Vrana on the island of Cres to reconstruct the paleoenvironmental conditions. Lake Vrana is the deepest freshwater lake in Croatia, located in the karst region of the eastern Adriatic coast. A dated sediment sequence in Lake Vrana of 4.4 m has spanned the past 16.4 kyr, and it featured a dynamic sediment deposition until the beginning of the Holocene, including strong sediment input and supply to the lake by runoff sediments of dolomitic origin from the catchment in the period 16.4–14.4 cal kyr BP. High organic carbon content, which originates from mixed terrestrial and aquatic origins in the periods 14.4–13.3 cal kyr BP and 12.7–11.7 cal kyr BP, indicates fluctuating lake levels in shallow water environments during the Late Glacial to Holocene transition. The Holocene sequence indicates the development of more stable conditions and continuous sediment deposition, characterized by an increasing trend of siliciclastic sediments delivered into the lake during the early Holocene (11.7–10 cal kyr BP) and dominantly from 8 to 4.4 cal kyr BP, indicating enhanced input and erosion, which coincides with the humid and pluvial period recorded in the central Mediterranean region. It is followed by sediments with high organic carbon content between 4.4 and 1.6 cal kyr BP, which points to higher lake productivity. Calcite sedimentation prevailed between 1.6 to 0.4 cal kyr BP, indicating stable deeper-lake conditions. Predominantly, siliciclastic sediments from 0.4 to 0.1 cal kyr BP pointed to erosion during the Little Ice Age (LIA), with enhanced precipitation and sediment discharge from the catchment. The re-establishment of calcite sedimentation has been observed over the last 100 years.