Geospatial Modeling of the Tombolo Phenomenon in Sopot using Integrated Geodetic and Hydrographic Measurement Methods

: A tombolo is a narrow belt connecting a mainland with an island lying near to the shore, formed as a result of sand and gravel being deposited by sea currents, most often created as a result of natural phenomena. However, it can also be caused by human activity, as is the case with the Sopot pier—a town located on the southern coast of the Baltic Sea in northern Poland ( ϕ = 54 ◦ 26’N, λ = 018 ◦ 33’E). As a result, the seaﬂoor rises constantly and the shoreline moves towards the sea. Moreover, there is the additional disturbing phenomenon consisting of the rising seaﬂoor sand covering over the waterbody’s vegetation and threatening the city's spa character. Removal of the sand to another place has already been undertaken several times. There is a lack of precise geospatial data about the tombolo’s seaﬂoor course, its size and spatial shape caused by only lowering the seaﬂoor in random places, and the ongoing environmental degradation process. This article presents the results of extensive and integrated geodetic and hydrographic measurements, the purpose of which was to make a 3D model of the phenomena developing in Sopot. The measurements will help determine the size and speed of the geospatial changes. Most of the modern geodetic and hydrographic methods were used in the study of these phenomena. For the construction of the land part of geospatial model, the following were used: photos from the photogrammetric ﬂight pass (unmanned aerial vehicle—UAV), laser scanning of the beach and piers, and satellite orthophotomaps for analysis of the coastline changes. In the sea part, bathymetric measurements were carried out with an unmanned surface vehicle (USV).


Introduction
In the last century, as a result of climate change observed on a global and local scale [1], there has been an increase in sea level. This leads to an increased impact of waves on the sea shore, causing the shore to move inland [2,3]. Therefore, coastal states are increasingly forced to take a number of actions to maintain the current coastline course [4] and protect coastal infrastructure [5]. In parallel, it should be ensured that the steps taken minimize interference with and the impact on the natural environment [6]. Depending on the value of LB*, two basic accumulation forms can be created on the beach: • Salient: When LB* is less than approx. 0.6 to 0.7, a bell-shaped salient in the shoreline will form in the lee of the breakwater. However, parameters other than the breakwater length and distance also influence the accumulation pattern. • Tombolo: When LB* is greater than approx. 0.9 to 1.0, the sand accumulation behind the breakwater will connect the beach to the breakwater in a tombolo formation. But again, parameters other than the breakwater length and distance influence the accumulation pattern. In the case of the marina in Sopot, LB* is 0.78, hence it is likely that a tombolo will eventually develop, which means that in the near future, a strip of land will be created connecting the pier with the shore and protruding above the water.
The purpose of this publication is to present the results of extensive and integrated geodetic and A detached breakwater (beach, coastal, offshore) as a marina in Sopot provides shelter from the waves, whereby the littoral transport behind the breakwater is decreased and the transport pattern adjacent to the breakwater is modified. These characteristics of a breakwater are utilised in different ways for various types of breakwaters by varying relevant parameters. The main parameter determining the future morphological form, which will be shaped as a result of the obstacle in Sopot, Remote Sens. 2020, 12, 737 3 of 18 is the breakwater length relative to breakwater distance to shoreline (LB*) calculated according to the formula [10]: where: LB-breakwater length, d-breakwater distance to shoreline. Depending on the value of LB*, two basic accumulation forms can be created on the beach: • Salient: When LB* is less than approx. 0.6 to 0.7, a bell-shaped salient in the shoreline will form in the lee of the breakwater. However, parameters other than the breakwater length and distance also influence the accumulation pattern. • Tombolo: When LB* is greater than approx. 0.9 to 1.0, the sand accumulation behind the breakwater will connect the beach to the breakwater in a tombolo formation. But again, parameters other than the breakwater length and distance influence the accumulation pattern.
In the case of the marina in Sopot, LB* is 0.78, hence it is likely that a tombolo will eventually develop, which means that in the near future, a strip of land will be created connecting the pier with the shore and protruding above the water.
The purpose of this publication is to present the results of extensive and integrated geodetic and hydrographic measurements that enabled creation of a 3D model of the phenomenon developing in Sopot. This will enable the size and speed of geospatial changes to be determined.

Study Area
In the literature on the subject, the process of modelling three-dimensional sediment transport has been described extensively and in detail in many publications [11,12]. It includes both the impact on the sea shore, as well as the influence on infrastructure and marine structures [13,14]. The diverse research approach to the authors' problem results from different methods of modelling being taken into account: evolution [15], fixed profile [16], coupling area-line [17], area cross-shore-alongshore transport coupling [11,18], diffusion [14], wet-dry [19], hybrid [20], and cut-cell models [19]. It should be noted with reference to [21] that general wave-average area models may also have merits in simulating shoreline evolution for arbitrary geometry. However, the difficulty in the treatment of the shoreward boundary has hindered the use of this type of approach for simulation of shoreline evolution. The major defect is that it is impossible to simulate morphologic change above the wave-average water level, although it is more or less hidden for a macro-tidal environment. A preliminary analysis of the tombolo phenomenon in Sopot was made in study [9]. On this basis, the city of Sopot developed an approximate model for the phenomenon development ( Figure 2) [22].
The developed model, or rather an explanatory figure, is of such a general nature that it prevents referring to coordinates and depth, which makes it impossible to precisely plan dredging works that the city of Sopot has been implementing for 3 years. Hence, direct measurements of the geospatial nature of the phenomenon are necessary from the point of view of assessing its size and planning dredging works for the seafloor (Figure 3) [22,23]. It should be noted that the city's activity in the field of digging sand from the seafloor should be conducted in such a way as to seek restoration sediment transport along the coast and restart natural processes. Therefore, the geospatial distribution of the beach and the seafloor near the pier plays a key role here.
Over the last several years (after the yacht marina was built), the adverse effect of the Sopot tombolo phenomenon on the aquatic environment and humans has been noted ( Figure 3) [22,23]. One of the local internet portals https://sopot.gmina.pl/ explains the cause of the phenomenon: "The emerging tombolo promotes the increasing occurrence of water stagnation and more rapid water heating, which creates more favourable conditions for the blooming of cyanobacteria and other microorganisms" [22]. Ms Magdalena Jachim, the spokeswoman for the Sopot City Hall, speaks in a similar tone: "This is not a thorough deepening, but a straightening of the coastline that will prevent the accumulation of waste and algae brought by the sea to the beach in this place" [24]. research approach to the authors' problem results from different methods of modelling being taken into account: evolution [15], fixed profile [16], coupling area-line [17], area cross-shore-alongshore transport coupling [11,18], diffusion [14], wet-dry [19], hybrid [20], and cut-cell models [19]. It should be noted with reference to [21] that general wave-average area models may also have merits in simulating shoreline evolution for arbitrary geometry. However, the difficulty in the treatment of the shoreward boundary has hindered the use of this type of approach for simulation of shoreline evolution. The major defect is that it is impossible to simulate morphologic change above the waveaverage water level, although it is more or less hidden for a macro-tidal environment. A preliminary analysis of the tombolo phenomenon in Sopot was made in study [9]. On this basis, the city of Sopot developed an approximate model for the phenomenon development ( Figure 2) [22]. The developed model, or rather an explanatory figure, is of such a general nature that it prevents referring to coordinates and depth, which makes it impossible to precisely plan dredging works that the city of Sopot has been implementing for 3 years. Hence, direct measurements of the geospatial nature of the phenomenon are necessary from the point of view of assessing its size and planning dredging works for the seafloor (Figure 3) [22][23]. It should be noted that the city's activity in the field of digging sand from the seafloor should be conducted in such a way as to seek restoration sediment transport along the coast and restart natural processes. Therefore, the geospatial distribution of the beach and the seafloor near the pier plays a key role here.
Over the last several years (after the yacht marina was built), the adverse effect of the Sopot tombolo phenomenon on the aquatic environment and humans has been noted ( Figure 3) [22][23]. One of the local internet portals https://sopot.gmina.pl/ explains the cause of the phenomenon: "The emerging tombolo promotes the increasing occurrence of water stagnation and more rapid water heating, which creates more favourable conditions for the blooming of cyanobacteria and other microorganisms" [22]. Ms Magdalena Jachim, the spokeswoman for the Sopot City Hall, speaks in a similar tone: "This is not a thorough deepening, but a straightening of the coastline that will prevent the accumulation of waste and algae brought by the sea to the beach in this place" [24]. The quoted statements and source materials clearly show that Sopot is likely to lose its status as an international health resort. This is primarily due to the blooming of cyanobacteria and other bacteria, which repeatedly prevents bathing and swimming in this region throughout the year. In the summer time, a swimming prohibition has been repeatedly issued by the Pomeranian Voivodeship Sanitary and Epidemiological Station. In order to determine the geospatial scale of the phenomenon, a photogrammetric flight pass was performed using an unmanned aerial vehicle (UAV) in the vicinity of the Sopot pier.
The measurements were carried out at the end of October and the beginning of November 2018 (a period unfavourable to water blooming). Based on the study results, it can be concluded that the range of cyanobacteria, algae, and other microorganisms occurring in this area is significant. On both orthophotos they are marked with a dark green/black colour. Opposite the Grand Hotel in Sopot, the above-mentioned, 30 m wide blooming area is located along a 200 m stretch of the beach. On the other hand, Figure 4b presents cyanobacteria and algae blooming along the entire length of the yacht marina in Sopot, which stretches for several dozen meters towards the mainland [25]. The quoted statements and source materials clearly show that Sopot is likely to lose its status as an international health resort. This is primarily due to the blooming of cyanobacteria and other bacteria, which repeatedly prevents bathing and swimming in this region throughout the year. In the summer time, a swimming prohibition has been repeatedly issued by the Pomeranian Voivodeship Sanitary and Epidemiological Station. In order to determine the geospatial scale of the phenomenon, a photogrammetric flight pass was performed using an unmanned aerial vehicle (UAV) in the vicinity of the Sopot pier.
The measurements were carried out at the end of October and the beginning of November 2018 (a period unfavourable to water blooming). Based on the study results, it can be concluded that the range of cyanobacteria, algae, and other microorganisms occurring in this area is significant. On both orthophotos they are marked with a dark green/black colour. Opposite the Grand Hotel in Sopot, the above-mentioned, 30 m wide blooming area is located along a 200 m stretch of the beach. On the other hand, Figure 4b presents cyanobacteria and algae blooming along the entire length of the yacht marina in Sopot, which stretches for several dozen meters towards the mainland [25].
Another factor having an adverse effect on the environment is the accumulation of sand in the vicinity of the marina. Despite the positive fact that the Sopot beach is expanding (under reconstruction), in other places the sand reaching the pier is taken away. The University of Gdańsk Professor Leszek Łęczyński, marine geology specialist, stated in an interview for the Gazeta Wyborcza paper that: "For years, we have been observing abrasion of the Orłowo Cliff, in my opinion, the material originating from there is being accumulated exactly near the Sopot pier" [26]. For example, in February 2018 as much as 1200 m 2 of ground originating from this cliff slid down onto a beach in Gdynia Orłowo. The landslide itself was approx. 40 m long [27]. It should also be noted that the Orłowo Cliff is part of the Kępa Redłowska nature reserve. Therefore, this area represents a landscape typical of Gdynia, with natural characteristics primarily associated with the following values: natural (the cliff has a unique multi-layered structure comprising glaciofluvial sands, gravels and silts, glacial till, sand-gravel sediments, etc. [28]), historical (this is the first nature reserve in Pomeranian Voivodeship and one of the oldest in the country), and landscape (from the cliff one can see a panorama of the Hel Peninsula, the harbour in Gdynia, the Sopot pier, etc. [29]).  Another factor having an adverse effect on the environment is the accumulation of sand in the vicinity of the marina. Despite the positive fact that the Sopot beach is expanding (under reconstruction), in other places the sand reaching the pier is taken away. The University of Gdańsk Professor Leszek Łęczyński, marine geology specialist, stated in an interview for the Gazeta Wyborcza paper that: "For years, we have been observing abrasion of the Orłowo Cliff, in my opinion, the material originating from there is being accumulated exactly near the Sopot pier" [26]. For example, in February 2018 as much as 1200 m 2 of ground originating from this cliff slid down onto a beach in Gdynia Orłowo. The landslide itself was approx. 40 m long [27]. It should also be noted that the Orłowo Cliff is part of the Kępa Redłowska nature reserve. Therefore, this area represents a landscape typical of Gdynia, with natural characteristics primarily associated with the following values: natural (the cliff has a unique multi-layered structure comprising glaciofluvial sands, gravels and silts, glacial till, sand-gravel sediments, etc. [28]), historical (this is the first nature reserve in Pomeranian Voivodeship and one of the oldest in the country), and landscape (from the cliff one can see a panorama of the Hel Peninsula, the harbour in Gdynia, the Sopot pier, etc. [29]).
The growing coastal belt also poses a navigational hazard to shipping traffic in this region. The vessels moored to the jetty located in the eastern part of the pier and in the southern part of the yacht marina are at risk of damaging their hull structure. This is due to the fact that for the vast majority of the berthing area, the water depth is actually lower than 1.5 m, and the draft of both sailing and motorised vessels often exceeds 1-1.5 m, thus preventing them from sailing on these waters. It should be added that sailors may not be aware of this fact, because the official electronic navigational chart (ENC) from 2011 ( Figure 5) shows that the isobaths increase linearly (and at equal distances) in relation to the beach coastline, and the water depth for the berthing area ranges from 2 to 5 m. Figure  5 shows a navigational chart (from OpenCPN 5.0 software), which uses official ENC cells used by marine vessels (ships, yachts, etc.). Therefore, it should be stressed that currently no official charts used in navigation include depth changes due to the tombolo phenomenon found near the Sopot pier. Hence, as part of the preliminary study, it was decided to carry out measurements to determine the current topography of the seafloor in this area. The growing coastal belt also poses a navigational hazard to shipping traffic in this region. The vessels moored to the jetty located in the eastern part of the pier and in the southern part of the yacht marina are at risk of damaging their hull structure. This is due to the fact that for the vast majority of the berthing area, the water depth is actually lower than 1.5 m, and the draft of both sailing and motorised vessels often exceeds 1-1.5 m, thus preventing them from sailing on these waters. It should be added that sailors may not be aware of this fact, because the official electronic navigational chart (ENC) from 2011 ( Figure 5) shows that the isobaths increase linearly (and at equal distances) in relation to the beach coastline, and the water depth for the berthing area ranges from 2 to 5 m. Figure 5 shows a navigational chart (from OpenCPN 5.0 software), which uses official ENC cells used by marine vessels (ships, yachts, etc.). Therefore, it should be stressed that currently no official charts used in navigation include depth changes due to the tombolo phenomenon found near the Sopot pier. Hence, as part of the preliminary study, it was decided to carry out measurements to determine the current topography of the seafloor in this area.

Measurement Concept
In terms of measurement, the area where the tombolo phenomenon occurs is specific, placing high demands on surveyors and hydrographers regarding the measurement precision, as well as a need to use various advanced measurement techniques, such as:

•
Land measurements (geodesy) enable the acquisition of geospatial data of a beach section adjacent to the sea. Measurements of this type can be carried out in two ways. In order to obtain a precise 3D model of the beach adjacent to the shore, the best solution is the use of the terrestrial laser scanning (TLS) [30][31][32]. However, this is a time-consuming method, it requires expensive measurement equipment and advanced methods of point cloud processing. The predicted positioning accuracy (3D) is 2-3 cm, for distances up to 100 m ( Figure 6a). An alternative to TLS is the use of photogrammetric methods, where the beach model is developed based on a sequence Remote Sens. 2020, 12, 737 6 of 18 of photos from the UAV [33][34][35]. The predicted model accuracy (3D) depends on the positioning accuracy of the global navigation satellite systems (GNSSs) used and the type of photogrammetric camera used, ranging from 5 cm to 1 m ( Figure 6b). • Sea surveys (hydrography) for depths below 1 m require the use of an atypical hydrographic vessel that has a very small draft and is equipped with a specialized echo sounder dedicated to ultra-shallow waters, the range of which should be at least 30 cm. The predicted accuracy of this solution allows a depth accuracy of 2-5 cm root mean square (RMS) (Figure 7a). An alternative solution to the use of a manned hydrographic vessel is the use of an unmanned surface vehicle (USV), which is increasingly used in bathymetric measurements [36][37][38]. The predicted accuracy of bathymetric measurements is identical to that of manned hydrographic measurements ( Figure 7b). • Sea surveys (hydrography) for depths: 0-0.6 m are specific because this part of the waterbody could not, for obvious reasons, be measured by a manned vessel. Therefore, no actual bathymetric data that could determine the spread of the tombolo phenomenon were available in this area. Currently, bathymetric measurements in ultra-shallow waters and on the coastline are carried out by either tachymetric [39] or geodetic method with a person stepping into the sea to a predefined depth using a GNSS receiver operating in real time [40][41][42].

Measurement Concept
In terms of measurement, the area where the tombolo phenomenon occurs is specific, placing high demands on surveyors and hydrographers regarding the measurement precision, as well as a need to use various advanced measurement techniques, such as: • Land measurements (geodesy) enable the acquisition of geospatial data of a beach section adjacent to the sea. Measurements of this type can be carried out in two ways. In order to obtain a precise 3D model of the beach adjacent to the shore, the best solution is the use of the terrestrial laser scanning (TLS) [30][31][32]. However, this is a timeconsuming method, it requires expensive measurement equipment and advanced methods of point cloud processing. The predicted positioning accuracy (3D) is 2-3 cm, for distances up to 100 m (Figure 6a). An alternative to TLS is the use of photogrammetric methods, where the beach model is developed based on a sequence of photos from the UAV [33][34][35]. The predicted model accuracy (3D) depends on the positioning accuracy of the global navigation satellite systems (GNSSs) used and the type of photogrammetric camera used, ranging from 5 cm to 1 m (Figure 6b).

Measurement Concept
In terms of measurement, the area where the tombolo phenomenon occurs is specific, placing high demands on surveyors and hydrographers regarding the measurement precision, as well as a need to use various advanced measurement techniques, such as: • Land measurements (geodesy) enable the acquisition of geospatial data of a beach section adjacent to the sea. Measurements of this type can be carried out in two ways. In order to obtain a precise 3D model of the beach adjacent to the shore, the best solution is the use of the terrestrial laser scanning (TLS) [30][31][32]. However, this is a timeconsuming method, it requires expensive measurement equipment and advanced methods of point cloud processing. The predicted positioning accuracy (3D) is 2-3 cm, for distances up to 100 m ( Figure 6a). An alternative to TLS is the use of photogrammetric methods, where the beach model is developed based on a sequence of photos from the UAV [33][34][35]. The predicted model accuracy (3D) depends on the positioning accuracy of the global navigation satellite systems (GNSSs) used and the type of photogrammetric camera used, ranging from 5 cm to 1 m (Figure 6b). • Sea surveys (hydrography) for depths below 1 m require the use of an atypical hydrographic vessel that has a very small draft and is equipped with a specialized echo sounder dedicated to ultra-shallow waters, the range of which should be at least 30 cm. The predicted accuracy of this solution allows a depth accuracy of 2-5 cm root mean square (RMS) (Figure 7a). An alternative solution to the use of a manned hydrographic vessel is the use of an unmanned surface vehicle (USV), which is increasingly used in  bathymetric measurements [36][37][38]. The predicted accuracy of bathymetric measurements is identical to that of manned hydrographic measurements (Figure 7b).
(a) (b) • Sea surveys (hydrography) for depths: 0-0.6 m are specific because this part of the waterbody could not, for obvious reasons, be measured by a manned vessel. Therefore, no actual bathymetric data that could determine the spread of the tombolo phenomenon were available in this area. Currently, bathymetric measurements in ultra-shallow waters and on the coastline are carried out by either tachymetric [39] or geodetic method with a person stepping into the sea to a predefined depth using a GNSS receiver operating in real time [40][41][42].

Geodetic Measurements
Creating a 3D beach model can be done on the basis of two methods: using UAV photos or taking measurements based on geodetic laser scanning [31]. The technique of measuring and creating a 3D beach model based on TLS compared to the method using a UAV is much more time consuming, but allows theoretically higher measurement accuracy.
In view of the elongated surface of the beach under measurement, it was necessary to plan and arrange an appropriate number of sites. The measurements were taken with a Trimble TX8 laser scanner without the photo-taking option. Hence, the obtained point clouds had only colours resulting from the calculated laser beam reflection intensity. To cover the assumed study area, it was necessary to establish 27 sites located at a distance of approx. 60 m from each other ( Figure 8). In view of the small number of characteristic objects to be used for the recording of a point cloud in the field at a later time, spherical tags (located in the sand in a manner ensuring the stability of their position) were used. The tags had to be located a relatively short distance from the neighboring measurement sites. Therefore, they were located halfway between the measurement sites or in their immediate vicinity. Such an approach ensured that a relatively large set of points on each spherical tag's surface were obtained during the measurements. This enabled the precise fitting of spheres into the set of points and the determination of their midpoints which, in the recording process, were the points of adjustment for particular local point cloud systems.

Geodetic Measurements
Creating a 3D beach model can be done on the basis of two methods: using UAV photos or taking measurements based on geodetic laser scanning [31]. The technique of measuring and creating a 3D beach model based on TLS compared to the method using a UAV is much more time consuming, but allows theoretically higher measurement accuracy.
In view of the elongated surface of the beach under measurement, it was necessary to plan and arrange an appropriate number of sites. The measurements were taken with a Trimble TX8 laser scanner without the photo-taking option. Hence, the obtained point clouds had only colours resulting from the calculated laser beam reflection intensity. To cover the assumed study area, it was necessary to establish 27 sites located at a distance of approx. 60 m from each other (Figure 8). In view of the small number of characteristic objects to be used for the recording of a point cloud in the field at a later time, spherical tags (located in the sand in a manner ensuring the stability of their position) were used. The tags had to be located a relatively short distance from the neighboring measurement sites. Therefore, they were located halfway between the measurement sites or in their immediate vicinity. Such an approach ensured that a relatively large set of points on each spherical tag's surface were obtained during the measurements. This enabled the precise fitting of spheres into the set of points and the determination of their midpoints which, in the recording process, were the points of adjustment for particular local point cloud systems. Combined point clouds coming from all stations formed a numerical representation of the environment. The connected point clouds contained measurement noise in the form of wrongly registered measurement points. Developing a point cloud often required manual removal of these measurement errors (Figure 9). Combined point clouds coming from all stations formed a numerical representation of the environment. The connected point clouds contained measurement noise in the form of wrongly registered measurement points. Developing a point cloud often required manual removal of these measurement errors (Figure 9). Combined point clouds coming from all stations formed a numerical representation of the environment. The connected point clouds contained measurement noise in the form of wrongly registered measurement points. Developing a point cloud often required manual removal of these measurement errors (Figure 9). The photogrammetric flight pass was performed using a DJI Mavic Pro drone independently for each of the adopted parts. The data were recorded using a Pix4D Capture mobile application. During the mission, 621 photographs were taken at an altitude of 60 m above the drone take-off level (i.e., the beach). The obtained average ground sample distance (GSD) coefficient value amounted to 2.25 cm/pixel (px) at a camera resolution of 4000 x 3000 effective px. The overlap parameter defining the degree of photograph overlapping in transverse and longitudinal directions was determined to be 80%. The main parameters of the camera are: •

Hydrographic Measurements
The marine vessel that was used in the research was the hydrographic survey motorboat Navigator One AMG (Figure 7a), a manned vessel intended for carrying out hydrographic surveys, with a length of 5.5 m, width of 2.55, draft of 0.3 m, a crew of 6, coastal buoyancy class, and a weight of 780 kg. It is equipped with an outboard motor with a power of 115 hp (Suzuki) and the Lowrance HDS Carbon navigation system with the StructureScan 3D module. The photogrammetric flight pass was performed using a DJI Mavic Pro drone independently for each of the adopted parts. The data were recorded using a Pix4D Capture mobile application. During the mission, 621 photographs were taken at an altitude of 60 m above the drone take-off level (i.e., the beach). The obtained average ground sample distance (GSD) coefficient value amounted to 2.25 cm/pixel (px) at a camera resolution of 4000 x 3000 effective px. The overlap parameter defining the degree of photograph overlapping in transverse and longitudinal directions was determined to be 80%. The main parameters of the camera are:

Hydrographic Measurements
The marine vessel that was used in the research was the hydrographic survey motorboat Navigator One AMG (Figure 7a), a manned vessel intended for carrying out hydrographic surveys, with a length of 5.5 m, width of 2.55, draft of 0.3 m, a crew of 6, coastal buoyancy class, and a weight of 780 kg. It is equipped with an outboard motor with a power of 115 hp (Suzuki) and the Lowrance HDS Carbon navigation system with the StructureScan 3D module.
The single beam echo sounder (SBES) was calibrated using the tare method, which involves adjusting the average speed of sound in water at known depths to echo sounder indications. The basic method for measuring the speed of sound is the use of a CTD (conductivity, temperature, depth) sensor or a sound velocity profiler (SVP). An alternative is the measurement of water temperature and the determination of the speed of sound based on an empirical equation when the salinity is known. Certain echo sounders determine the speed of sound in water based on the temperature by means of a sensor located in the transducer.
Oceanographic measurements which are part of hydrographic surveys include water level observations and measurements of the speed of sound in water [43]. The water level can be observed on a gauging station which is situated in a marina or port and is referenced to the geodetic height system used in a particular country. In Poland, according to the current regulations [44], water levels are referenced to the height system of Kronstadt (PL-KRON86-NH) (until 2019) and Amsterdam (PL-EVRF2007-NH). For the area of Sopot marina surveys, the closest gauging stations are Gdańsk Port Północny and Gdynia. Interpolation between observations at distant points carries the risk that in the location of bathymetric measurements, the actual water level deviates from the interpolated level due to e.g., wind causing water lifting. Another solution is to periodically measure the water level above the sea level using a GNSS receiver and a geoid model.
During the surveys, the water levels recorded at the Gdańsk Port Północny (at a distance of 9.7 km from Sopot) and Gdynia (at a distance of 7.8 km from Sopot) gauging stations were identical and were adopted for the reduction of depth measurement to the chart datum. The methodology for considering water level measurements during hydrographic surveys is described in detail in [45,46].
Sounding profiles have been planned for the marina and pier areas, taking into account the following assumptions: • Area includes a rectangle 800 m long, calculated along the coast and symmetrical to the pier and 500 m wide corresponding to the pier length, excluding the marina ( Figure 10). Oceanographic measurements which are part of hydrographic surveys include water level observations and measurements of the speed of sound in water [43]. The water level can be observed on a gauging station which is situated in a marina or port and is referenced to the geodetic height system used in a particular country. In Poland, according to the current regulations [44], water levels are referenced to the height system of Kronstadt (PL-KRON86-NH) (until 2019) and Amsterdam (PL-EVRF2007-NH). For the area of Sopot marina surveys, the closest gauging stations are Gdańsk Port Północny and Gdynia. Interpolation between observations at distant points carries the risk that in the location of bathymetric measurements, the actual water level deviates from the interpolated level due to e.g. wind causing water lifting. Another solution is to periodically measure the water level above the sea level using a GNSS receiver and a geoid model.
During the surveys, the water levels recorded at the Gdańsk Port Północny (at a distance of 9.7 km from Sopot) and Gdynia (at a distance of 7.8 km from Sopot) gauging stations were identical and were adopted for the reduction of depth measurement to the chart datum. The methodology for considering water level measurements during hydrographic surveys is described in detail in [45][46].
Sounding profiles have been planned for the marina and pier areas, taking into account the following assumptions: • Area includes a rectangle 800 m long, calculated along the coast and symmetrical to the pier and 500 m wide corresponding to the pier length, excluding the marina (Figure 10).

Geodetic Measurement Results
A photogrammetric model was developed using the Pix4d Mapper Pro software which enables photographic data processing and generation based on three-dimensional models and orthophotomaps. A median of 21,128 nodal points on a single photograph was obtained, which indicates a relatively high number of characteristic points and areas, considering the frequently little diversified coverage of the area surface. Figure 11 shows a 3D model built based on the UAV photos [25].
A comparative picture of point clouds is illustrated in Figure 12, obtained from measurements using TLS and UAV methods.
To assess the accuracy of the 3D beach model made, it is necessary to compare it to the state geodetic coordinate system used in Poland (Figure 13) [47].
A photogrammetric model was developed using the Pix4d Mapper Pro software which enables photographic data processing and generation based on three-dimensional models and orthophotomaps. A median of 21,128 nodal points on a single photograph was obtained, which indicates a relatively high number of characteristic points and areas, considering the frequently little diversified coverage of the area surface. Figure 11 shows a 3D model built based on the UAV photos [25]. A comparative picture of point clouds is illustrated in Figure 12, obtained from measurements using TLS and UAV methods.  indicates a relatively high number of characteristic points and areas, considering the frequently little diversified coverage of the area surface. Figure 11 shows a 3D model built based on the UAV photos [25]. A comparative picture of point clouds is illustrated in Figure 12, obtained from measurements using TLS and UAV methods.  To assess the accuracy of the 3D beach model made, it is necessary to compare it to the state geodetic coordinate system used in Poland (Figure 13) [47]. In the discussed TLS case study, the detailed inaccuracy components of the applied laser scanner measurement are equal to [47]: • Inaccuracy of the Trimble Tx8 laser scanner-1 cm horizontally and 1.5 cm vertically.
• Inaccuracy of the Trimble R10 GNSS receiver used to determine the coordinates of reference points-1 cm horizontally and 1.5 cm vertically. In the discussed TLS case study, the detailed inaccuracy components of the applied laser scanner measurement are equal to [47]: • Inaccuracy of the Trimble Tx8 laser scanner-1 cm horizontally and 1.5 cm vertically.

•
Inaccuracy of the Trimble R10 GNSS receiver used to determine the coordinates of reference points-1 cm horizontally and 1.5 cm vertically.

•
Inaccuracy of the PL-2000 system-2 cm horizontally and vertically.
A cloud of TLS points after transformation to the PL-2000 system, including the data collected at the determined area (terrain), can be a reference set for the cloud of points obtained by the UAV method in the form of terrain digital photos. A series of photos, processed by the appropriate program, enables creation of a space projection of the investigated area (terrain), in which besides the X-Y coordinates, the Z coordinate is determined in relation to each point of the cloud. On the basis of the reference points coordinates, the obtained local coordinates are transformed in the PL-2000 system. In practice, the sequence of operations realized on the basis of photo data does not allow achievement of a satisfactory level of accuracy of the data on the investigated area [47].
In the UAV method of terrain imaging, the picture processing operation resulted in the generation of a cloud containing 20,666,253 points, which translates into a density of approx. 117 points per m 2 . Numerical data processing was carried out using a high-parameter workstation (16 GB RAM, i7-6600U, GTX 1070) which lasted for 3 hours 42 minutes. The following values of errors for particular coordinates were obtained: 2.83 m (X), 4.08 m (Y) and 9.81 m (Z) [47].
In the discussed UAV case study, the detailed inaccuracy components of the applied unmanned aerial vehicle are equal to [47]: • Camera px resolution-2.25 cm.

•
Inaccuracy of UAV position: X, Y-a few meters, Z-several meters.
To estimate the quality of the UAV method measurement, it is proposed to use the TLS data cloud as reference data. These data are obtained in the same area as the data in the UAV data cloud. Then, it is possible to estimate the systematic X-Y-Z coordinates errors components in the UAV cloud, as well as their dispersion (noise) in relation to reference values ( Figure 14) [47]. Dedicated software is currently under construction to allow for the determination of the uncertainty characteristics of X-Y-Z coordinates collected in large data files obtained using the UAV method [47].
Results of the executed measurements of the terrain geometry are treated as random variables, that is, they are subject to statistical rules and the probability calculus is used to assess them. The term "metrological aspects" in this paper concerns the estimation of uncertainty in relation to the measurements of the sea shore and the seafloor under investigation [48]. Measurement results obtained with use of the aforementioned methods defined as TLS and UAV, respectively, create the series (clouds) of measurement points. Each of them has different references in the X-Y-Z coordinates system. Consequently, it is possible to estimate an uncertainty of determination of space coordinates of each measurement point in cloud by applying the methodology of type B uncertainty determining [49]. Uncertainty of type B concerns mainly the estimation of inaccuracy of measurement Dedicated software is currently under construction to allow for the determination of the uncertainty characteristics of X-Y-Z coordinates collected in large data files obtained using the UAV method [47].
Results of the executed measurements of the terrain geometry are treated as random variables, that is, they are subject to statistical rules and the probability calculus is used to assess them. The term "metrological aspects" in this paper concerns the estimation of uncertainty in relation to the measurements of the sea shore and the seafloor under investigation [48]. Measurement results obtained with use of the aforementioned methods defined as TLS and UAV, respectively, create the series (clouds) of measurement points. Each of them has different references in the X-Y-Z coordinates system. Consequently, it is possible to estimate an uncertainty of determination of space coordinates of each measurement point in cloud by applying the methodology of type B uncertainty determining [49]. Uncertainty of type B concerns mainly the estimation of inaccuracy of measurement instrumentation. Its calculation requires determination of the absolute error threshold value ∆X t and assumption of the given function of probability distribution for this error. Then, for a given component, the type B standardized uncertainty u s (x) is expressed as for normal distribution [47]: or for uniform distribution: where: ∆X t -maximum permissible error threshold value. Taking into account that estimation of the type B measurement procedure may concern a few inaccuracy components characterized by the standard uncertainties u s1 (x), u s2 (x), ..., and u sn (x), we can calculate the complex uncertainty u c (x) as [47]: where: u s1 (x), u s2 (x), ..., and u sn (x) correspond to all sources of errors in the considered measurement procedure.
Then, for a given component, the type B extended uncertainty U B (x) is expressed as [47]: where: k-coefficient of extension (coverage factor). First of all, the detailed uncertainty components can be analysed and calculated for the method TLS and related to its measurement instrumentation. Taking into consideration known instrumental and method data, the following error components are considered [47]:
The uncertainty should be estimated for all three coordinates (X, Y, and Z). The uncertainty budget is shown in Table 1 for X or Y coordinates and in Table 2 for Z coordinate [47]. Assuming the normal distribution of complex uncertainties for all coordinates, the extended uncertainty (5) u TLS is 2.82 cm (X or Y coordinate) and 3.36 cm (Z coordinate) (p = 0.95) [47].

Hydrographic Measurement Results
Hydrographic measurements enabled creation of the bathymetric chart of the tombolo phenomenon in Sopot. This map was developed based on measurements realised in October 2018. In addition, bathymetric data (isobaths) from the 2011 ENC were processed and supplemented with a point cloud from the beach and pier (Figure 15).  Assuming the normal distribution of complex uncertainties for all coordinates, the extended uncertainty (5) uTLS is 2.82 cm (X or Y coordinate) and 3.36 cm (Z coordinate) (p = 0.95) [47].

Hydrographic Measurement Results
Hydrographic measurements enabled creation of the bathymetric chart of the tombolo phenomenon in Sopot. This map was developed based on measurements realised in October 2018. In addition, bathymetric data (isobaths) from the 2011 ENC were processed and supplemented with a point cloud from the beach and pier (Figure 15).   Figure 15a, it should be noted that the berthing area located in the eastern part of the pier and in the southern part of the yacht marina is not suitable for navigation due to the very shallow depths (1-1.5 m) occurring there. It should also be stressed that the authors carried out bathymetric   Figure 15a, it should be noted that the berthing area located in the eastern part of the pier and in the southern part of the yacht marina is not suitable for navigation due to the very shallow depths (1-1.5 m) occurring there. It should also be stressed that the authors carried out bathymetric measurements of the areas in the vicinity of the Sopot pier using a manned vessel Navigator One AMG, which has one of the shallowest drafts possible (30 cm). Nevertheless, the vessel's bottom repeatedly scraped the sand on the seafloor, which posed a risk of damage to both the vessel and the measuring equipment. These measurements clearly proved that sailing on such vessels (even with the shallowest draft possible) in ultra-shallow waters is risky and poses a hazard of damaging the measurement systems [50,51].
The shape of the isobaths from measurements of the tombolo phenomenon (Figure 15a) indicates that as a result of reduced sand transport along the coast, sand began to accumulate on the right side of the pier in Sopot. The result of this state is that the seafloor has significantly increased over the entire length of the waterbody connecting the marina with the shore. At the inner marina, the depths decreased by 1.5-2 m. Halfway between the marina and the shore, the seafloor rose by about 1 m. The width of the peninsula, which will be built in the future, can be estimated at about 90-100 m, at the water surface.
The degree of advancement of the progressive tombolo phenomenon can also be estimated by comparing two digital terrain models (DTMs) [52] from 2011 and 2018 (measured). As a result of their comparison, the volume sand deposited over 8 years on the pier area in Sopot was determined. Figure 16 presents a map in which depth differences between the tombolo DTM models developed on the basis of measurements and the ENC data are shown in colour. The volume sand deposited is 31,963,450 m 3 . This means that to remove all the sand accumulated with the tombolo phenomenon, 2663 large dump trucks with a capacity of 12 m 3 would have to be used. decreased by 1.5-2 m. Halfway between the marina and the shore, the seafloor rose by about 1 m. The width of the peninsula, which will be built in the future, can be estimated at about 90-100 m, at the water surface.
The degree of advancement of the progressive tombolo phenomenon can also be estimated by comparing two digital terrain models (DTMs) [52] from 2011 and 2018 (measured). As a result of their comparison, the volume sand deposited over 8 years on the pier area in Sopot was determined. Figure  16 presents a map in which depth differences between the tombolo DTM models developed on the basis of measurements and the ENC data are shown in colour. The volume sand deposited is 31,963,450 m 3 . This means that to remove all the sand accumulated with the tombolo phenomenon, 2663 large dump trucks with a capacity of 12 m 3 would have to be used. The measurements and bathymetric charts made it possible to determine the areas of dredging works that are planned by the city of Sopot. It is noteworthy that the width of the depth changes due to sand accumulation is about 300 m.

Additional Research Results
In order to determine the scale of the tombolo phenomenon, it was decided to analyse the coastline course in Sopot over the last 10 years. To this end, aerial photographs from the website https://mapa.trojmiasto.pl/ provided by the MGGP Aero (years 2008, 2011, 2014, and 2016) and satellite photographs from Google Earth Pro, a popular geographical photo and information sharing application, were used. Vectorisation and cartographic adjustment of the photographs in terms of the inclination angle, rotation, and scale helped to identify considerable morphological changes in the coast in the vicinity of Sopot ( Figure 17) [25]. The measurements and bathymetric charts made it possible to determine the areas of dredging works that are planned by the city of Sopot. It is noteworthy that the width of the depth changes due to sand accumulation is about 300 m.

Additional Research Results
In order to determine the scale of the tombolo phenomenon, it was decided to analyse the coastline course in Sopot over the last 10 years. To this end, aerial photographs from the website https://mapa.trojmiasto.pl/ provided by the MGGP Aero (years 2008, 2011, 2014, and 2016) and satellite photographs from Google Earth Pro, a popular geographical photo and information sharing application, were used. Vectorisation and cartographic adjustment of the photographs in terms of the inclination angle, rotation, and scale helped to identify considerable morphological changes in the coast in the vicinity of Sopot ( Figure 17) [25].
Moreover, the conducted analyses enabled both the calculation of the Sopot beach area in selected years and the determination of its area variability degree over the years (Table 3) [25].  Moreover, the conducted analyses enabled both the calculation of the Sopot beach area in selected years and the determination of its area variability degree over the years (Table 3) [25].

Discussion
As a result of the construction of the marina in Sopot, there was a change in the morphological forms of the coast around the pier on a sandy beach of 800 m length. Accumulating sediment transport along the beach has caused the seafloor to rise, which has lead to a tombolo in the proceeding years (a narrow belt connecting the mainland with an island lying near the shore formed as a result of sand and gravel being deposited by sea currents). In addition to coastline changes, this phenomenon has led to increased vegetation of this area, threatening the city's spa character. This publication presents research results aimed at integrating measurement methods and systems of modern geodesy and hydrography for geospatial modeling of the tombolo phenomenon.
As part of the study, TLS and UAV beach measurements were carried out with an analysis of their accuracy, coastline changes were determined based on the analysis of satellite images, bathymetric measurement of the waterbody adjacent to the sea was carried out, and bathymetric charts of the tombolo phenomenon were made and they were compared with data from the ENC. Finally, the volume and area of accumulating sand in the developing tombolo was determined.

Discussion
As a result of the construction of the marina in Sopot, there was a change in the morphological forms of the coast around the pier on a sandy beach of 800 m length. Accumulating sediment transport along the beach has caused the seafloor to rise, which has lead to a tombolo in the proceeding years (a narrow belt connecting the mainland with an island lying near the shore formed as a result of sand and gravel being deposited by sea currents). In addition to coastline changes, this phenomenon has led to increased vegetation of this area, threatening the city's spa character. This publication presents research results aimed at integrating measurement methods and systems of modern geodesy and hydrography for geospatial modeling of the tombolo phenomenon.
As part of the study, TLS and UAV beach measurements were carried out with an analysis of their accuracy, coastline changes were determined based on the analysis of satellite images, bathymetric measurement of the waterbody adjacent to the sea was carried out, and bathymetric charts of the tombolo phenomenon were made and they were compared with data from the ENC. Finally, the volume and area of accumulating sand in the developing tombolo was determined.
The analysis results can be used by the local government in Sopot to solve this important problem for the town, as well as develop a systemic approach to sand removal that is very expensive. It is recommended that dredging works be carried out only in selected places and in such a way that restores sediment transport along the beach, because its complete cessation would accelerate the process of tombolo development, as well as have a significant impact on the waterbody's vegetation.
Based on Figure 17, it can be concluded that the municipal beach area has been steadily increasing, and the coastline to the east of the pier is increasingly approaching the yacht marina. It can also be seen that in the years 2008-2018, the coastal belt area in this place increased by approx. 13,062 m 2 , and the coastline shifted towards the sea by 53 m [25]. This has happened mainly as a result of the construction of an additional breakwater protecting the pier against storms. Therefore, it is necessary to conduct continuous cartographic monitoring in this area.
It should also be emphasized that hydrographic measurements with a manned vessel in such shallow water are difficult and can often be unachievable due to the small depths (below 1 m). Hence, future work has focused on the use of the hydrographic USV, which was successfully carried out in December 2019.
The proposed research methodology involving integration of geodetic and hydrographic methods aimed at creating a 3D tombolo model can be successfully used to assess other phenomena of this type.