A Cost-E ﬀ ective Method to Reproduce the Morphology of the Nearshore and Intertidal Zone in Microtidal Environments

: The photogrammetric method is widely used in coastal areas and in submerged environments. Time-lapse images collected with unmanned aerial vehicles are used to reproduce the emerged areas, while images taken by divers are used to reproduce submerged ones. Conversely, 3D models of natural or human-made objects lying at the water level are severely a ﬀ ected by the di ﬀ erence in refractive index between air and seawater. For this reason, the matching of 3D models of emergent and submerged coasts has been very rarely tested and never used in Earth Sciences. The availability of a large number of time-lapse images, collected at the intertidal zone during previous snorkel surveys, encouraged us to test the merging of 3D models of emerged and submerged environments. Considering the rapid and e ﬀ ective nature of the aforementioned program of swim surveys, photogrammetric targets were not used during image acquisition. This forced us to test the matching of the independent models by recognizing prominent landforms along the waterline. Here we present the approach used to test the method, the instrumentation used for the ﬁeld tests, and the setting of cameras ﬁxed to a specially built aluminum support console and discuss both its advantages and its limits compared to UAVs. 3D models of sea cli ﬀ s were generated by applying structure-from-motion (SfM) photogrammetry. Horizontal time-lapse images, collected with action cameras while swimming parallel to the coastline at nearly constant velocity, were used for the tests. Subsequently, prominent coastal landforms were used to couple the independent models obtained from the emergent and submerged cli ﬀ s. The method was pilot tested in two coastal sites in the north-eastern Adriatic (part of the Mediterranean basin). The ﬁrst site was a 25 m sea wall of sandstone set within a small harbor, while the second site was a 150 m route below plunging limestone cli ﬀ s. The data show that inexpensive action cameras provide a su ﬃ cient resolution to support and integrate geomorphological ﬁeld surveys along rocky coastlines.


Introduction
Rocky coasts are the complex outcome of the work of marine and subaerial processes that result in a variety of landforms [1], in particular in the intertidal zone. This is the area of the coast that is included within the intertidal range or is underwater at high tide and above water level at low tide. Above the intertidal zone is the supratidal zone, spray zone, or splash zone, which is covered by water only during storms [2]. The shape of rocky coasts in the Mediterranean area can be classified as plunging cliffs, sloping coasts, as well as two types of shore platforms according to the shape of the platforms at the toe of the cliff, screes, and composite coasts [3,4]. Rocky coasts represent a key environment to test

Description of the Study Area
Both sites are located on the north-eastern Adriatic coast (Italy and Croatia, northern Mediterranean), as shown in Figure 1A. Site 1, shown in Figure 1B, is a small harbor within a sheltered bay in Italy, south of Trieste and at the border of Slovenia. The survey was performed on the sandstone seawall of a military base. The area is dominated by carbonate and marly-sandstone rocks belonging to

Support Console and Cameras
The support console consisted of a set of 20 × 20 mm aluminum bars, shown in Figure 2, fixed together for their purpose, on which three cameras were attached. The support was hosted on a raft 1.2 m long, used within the Geoswim program, which is well-described by [5]. One GoPro Hero Black 5 was fixed on the support at +0.35 m awl, and up to two GoPRO Hero Black 6s were fixed at − 0.35 m bwl within plastic domes produced by Telesin, shown in Figure 3, in order to preserve the field of view (FOV) and focal length of the camera-lens system [32].

Support Console and Cameras
The support console consisted of a set of 20 × 20 mm aluminum bars, shown in Figure 2, fixed together for their purpose, on which three cameras were attached. The support was hosted on a raft 1.2 m long, used within the Geoswim program, which is well-described by [5]. One GoPro Hero Black 5 was fixed on the support at +0.35 m awl, and up to two GoPRO Hero Black 6s were fixed at −0.35 m bwl within plastic domes produced by Telesin, shown in Figure 3, in order to preserve the field of view (FOV) and focal length of the camera-lens system [32].
We used one camera bwl at site 1, because the sea wall is vertical and homogeneous. Conversely, since the plunging cliffs at site 2 are more jagged, we used two cameras bwl, in order to increase the number of available images. This allowed us to avoid an unpleasant lack of images in the event of technical malfunctions, and to get more images of the same targets of the cliff, in order to reduce shading effects and to choose the best images. The positions of the cameras with respect to the waterline were chosen in order to better cover the intertidal zone and surroundings. The mutual orientation of the cameras bwl and their depth with respect to the waterline were chosen in order to be sure of covering the widest target surface, as shown in Figure 4. The cameras were set in fish-eye mode, or ultra-wide angle, both awl and bwl. The camera inclination can be adjusted, but it was set in a horizontal position and perpendicularly to the swimming direction. The cameras bwl were positioned in the same way, but each was 5 • convergent, as shown in Figures 3 and 4.  We used one camera bwl at site 1, because the sea wall is vertical and homogeneous. Conversely, since the plunging cliffs at site 2 are more jagged, we used two cameras bwl, in order to increase the number of available images. This allowed us to avoid an unpleasant lack of images in the event of technical malfunctions, and to get more images of the same targets of the cliff, in order to reduce shading effects and to choose the best images. The positions of the cameras with respect to the waterline were chosen in order to better cover the intertidal zone and surroundings. The mutual orientation of the cameras bwl and their depth with respect to the waterline were chosen in order to be sure of covering the widest target surface, as shown in Figure 4. The cameras were set in fish-eye mode, or ultra-wide angle, both awl and bwl. The camera inclination can be adjusted, but it was set in a horizontal position and perpendicularly to the swimming direction. The cameras bwl were positioned in the same way, but each was 5° convergent, as shown in Figures 3 and 4.  . Time-lapse images collected in the laboratory corridor. The first line corresponds to the camera above the waterline, while the second and third lines collected images at lower elevations that correspond to the submerged setting. As aforementioned, the latter cameras are slightly convergent.

Acquisition Setting and Data Collection
Time-lapse images, shown in Figure 4, were collected at two sites along the north-east Adriatic coast: Site 1-Lazzaretto harbor (Muggia, Trieste, Italy); Site 2-Premantura (Pula, Croatia), as shown in Figure  1. Moreover, a laboratory test at the University of Trieste was carried out in order to calibrate the settings of the cameras.
Below, we describe the procedures developed for collecting time-lapse images in the field and the strategies designed to optimize image quality and overlay.
The GoPROs were set in time-lapse mode, with one shoot every 0.5 s or 1 s bwl, and 1 s awl, as shown in Figure 4. The time-lapse images were collected while swimming parallel to the coast, at an average distance of 1 ± 0.5 m away from shore, at a speed around 1 km/h. In correspondence to topographical obstructions (such as submerged or semi-submerged rocks, boulders, etc.), the route was changed slightly from that planned, but still guaranteed photo coverage by reducing the speed, for example. However, thanks to the characteristics of the software used for reconstructing the 3D models, small differences in swimming speed do not affect the resolution of the models. GoPROs are so-called "actioncams", mainly devoted to creating videos of extreme sport activities. The GoPRO Hero 5 and 6 are fixed-focus still-frame cameras, (or 1080p high-definition video cameras), whose small form-factor and waterproof housing have made them popular with sports enthusiasts (GoPRO technical manuals). Both cameras have a resolution of 12 megapixels (4000 × 3000) in discrete image capture. They have a continuous shooting mode (time-lapse mode) with a nominal maximum rate of two images per second. The nominal focal length of the GoPro camera is 3 mm and the pixel size is 1.73 μm. Assuming no radial distortion, the area covered by a horizontal image taken 1 m from the coastline is 3.0 m × 4.0 m, which . Time-lapse images collected in the laboratory corridor. The first line corresponds to the camera above the waterline, while the second and third lines collected images at lower elevations that correspond to the submerged setting. As aforementioned, the latter cameras are slightly convergent.

Acquisition Setting and Data Collection
Time-lapse images, shown in Figure 4, were collected at two sites along the north-east Adriatic coast: Site 1-Lazzaretto harbor (Muggia, Trieste, Italy); Site 2-Premantura (Pula, Croatia), as shown in Figure 1. Moreover, a laboratory test at the University of Trieste was carried out in order to calibrate the settings of the cameras.
Below, we describe the procedures developed for collecting time-lapse images in the field and the strategies designed to optimize image quality and overlay.
The GoPROs were set in time-lapse mode, with one shoot every 0.5 s or 1 s bwl, and 1 s awl, as shown in Figure 4. The time-lapse images were collected while swimming parallel to the coast, at an average distance of 1 ± 0.5 m away from shore, at a speed around 1 km/h. In correspondence to topographical obstructions (such as submerged or semi-submerged rocks, boulders, etc.), the route was changed slightly from that planned, but still guaranteed photo coverage by reducing the speed, for example. However, thanks to the characteristics of the software used for reconstructing the 3D models, small differences in swimming speed do not affect the resolution of the models. GoPROs are Remote Sens. 2020, 12, 1880 6 of 21 so-called "actioncams", mainly devoted to creating videos of extreme sport activities. The GoPRO Hero 5 and 6 are fixed-focus still-frame cameras, (or 1080p high-definition video cameras), whose small form-factor and waterproof housing have made them popular with sports enthusiasts (GoPRO technical manuals). Both cameras have a resolution of 12 megapixels (4000 × 3000) in discrete image capture. They have a continuous shooting mode (time-lapse mode) with a nominal maximum rate of two images per second. The nominal focal length of the GoPro camera is 3 mm and the pixel size is 1.73 µm. Assuming no radial distortion, the area covered by a horizontal image taken 1 m from the coastline is 3.0 m × 4.0 m, which corresponds to a ground sampling distance (GSD) of 1.7 mm [25]. The GoPROs were set to collect time-lapse images every second, or 3600 images/hour. Swimming at 1 km/h (about 0.3 m/s), the camera shoots every 0.3 m, as shown in Figure 5. Overlap exceeds 90% at the aforementioned selected swimming speeds. The system is mainly intended for corridor mapping along the waterline with a single strip. If needed, a camera inside the dome in the bow, shown in Figure 6A,C,D, at the water level can be used to film the survey operations.
The distortion index in underwater photography was calibrated using a print target made using a regular W/B patch grid, shown in Figure 6B.  Overlap exceeds 90% at the aforementioned selected swimming speeds. The system is mainly intended for corridor mapping along the waterline with a single strip. If needed, a camera inside the dome in the bow, shown in Figure 6A,C,D, at the water level can be used to film the survey operations. Overlap exceeds 90% at the aforementioned selected swimming speeds. The system is mainly intended for corridor mapping along the waterline with a single strip. If needed, a camera inside the dome in the bow, shown in Figure 6A,C,D, at the water level can be used to film the survey operations.
The distortion index in underwater photography was calibrated using a print target made using a regular W/B patch grid, shown in Figure 6B.  The distortion index in underwater photography was calibrated using a print target made using a regular W/B patch grid, shown in Figure 6B.
The grid, or patch size, represents the relative size of each box compared to the neighbor. Images of the target were collected in every part of the printed patch, summarizing the entire surface of interest, in order to obtain a perfect correction model. Images were georeferenced and assigned on the map using the integrated GPS tags of the camera awl. The GoPRO GPS precision is not provided by the company, but usually it is 1 m under open sky, and worsens close to obstacles, such as sea caves, etc. GPS data were used only to accelerate the process of tie point generation, since the software compares only nearby images [Technical Manual of Metashape, 2019].

Post-Processing of Time-Lapse Images
Time-lapse images were included in the photogrammetric process to obtain 3D models. They have become a key concept in geosciences for modelling topographic surfaces, and for the monitoring of spatial and temporal changes in natural environments [33]. A wide range of commercial and open source software tools enable the generation of 3D models of complex geomorphological features, including by geoscientists and other non-experts [34]. In particular, SfM is a photogrammetric range imaging technique for estimating three-dimensional structures from two-dimensional image sequences. SfM offers Earth scientists a tool for building high-resolution models at a range of scales and for detection purposes. SfM allows the use of consumer-grade digital cameras and highly automated data processing, and it is based upon algorithms for the efficient and automatic orientation of large image datasets without further acquisition information [34].
3D structures can be resolved from a series of overlapping, offset images. The automatic camera coordinates estimation algorithms simultaneously allow the operators to identify 3D camera location and scene geometry [20]. Automatic internal and external orientation are important innovations in the SfM method. Digital photogrammetry works using recursive analysis to determine the external orientation from the inner one, using the images' metadata. Normally, photographs are distorted in some way with respect to a perfect monoscopic projection. The distortion of the images allows analysts to determine where and which points are not in the right position on the photo (internal orientation of the camera). These parameters change with the focal length and, therefore, a pixel situated in two different positions in two different images can be represented with a different projection. Metadata are processed by software, so it is also possible to obtain a correction profile starting from them and computing the average error from every single image of a dataset.
Color corrections were performed using Photoshop CC 2019, while images were processed using Agisoft Metashape 1.5.0 workflow [35], in order to create the mesh models.
Tie points (TPs) were created using Agisoft Metashape. They represent the number of recognized features and define the total population of common points for each image. The software also creates the connections between them. Then, the dense cloud process defines the model at a more detailed level, creating up to 93 times more points than the tie point cloud. The final 3D model is a set of multiple triangles linking the points. The model can be displayed using a color film, called texture, over each mesh triangle.
The emergent and underwater models were obtained independently from images collected awl and bwl, thus each model has its own arbitrary datum. In order to join up the emergent and submerged part of the coastline, at the end of the procedure the independent models were oriented and jointed using a common coordinates system.
Menna et al. [11] suggested joining up the separated models using calibrated Orientation Devices (ODs), such as rigid bars. We did not use this technique as much as the image dataset was collected during very rapid expeditions. For this reason, the models were linked up using prominent features of the coastline instead of ODs, as at site 2, shown in Figure 4. A 2 m-long aluminum bar was used as a vertical reference to define the model orientation, to determine the z axis, and to scale the models.

UAV Comparison
To validate the models, we performed a UAV survey using a DJI Spark drone (DJI, Nanshan District, Shenzhen, China). This UAV is categorized as a micro drone, according to the weight classification proposed by [36]. Micro drones are widely used in various fields of geoscience and offer major benefits for their maneuverability during flight operations and landing phases. In addition, these drones can provide high-resolution images from reduced flight times [37,38]. Unmanned aerial vehicles-digital photography (UAV-DP)-enabled the reconstruction of high-resolution point clouds, orthophotographs and 3D models of the sea cliffs studied. These models were compared with that produced using horizontal images collected awl from the raft.
Shots were taken every 2 s and the drone was piloted manually 25 m above sea level and approximately 15 m from the plunging cliffs with the camera orientation set at approximately 45 • towards the terrain in order to avoid a lack of coverage on the vertical structure [10].
The model was produced through the processing of 247 images acquired by the drone in February 2020, with the same procedure followed with horizontal time-lapse images.

Weather Data during the Surveying Days
Weather data for the surveying days were obtained, respectively from the meteorological station at Muggia (ARPA FVG), 5 km from site 1, and from the station of the Croatian Meteorological and Hydrological Service of Pula, 13 km from site 2. The data are reported in Table 1. No rainfall events were observed, while the sea and air temperatures conformed to seasonal data and ranged between 25.8 • C and 27.5 • C. The highest wind speed blew from the second quadrant, shown in Table 1. The tide elevations during the surveys were measured and corrected using the tide gauge in Trieste. No correction was applied for site 1, because it is located close to the gauge. A 25 min advance was applied for site 2, considering the 70% tide amplitude in Trieste [39,40]. At site 1, the tide was −5 cm below mean sea level, and at site 2 during the June survey, the tide was 0 cm above mean sea level, while at site 2 during the UAV survey, the tide was +5 cm above mean sea level.

Case Studies
Data from site 1 are described in Section 4.1, site 2 in Section 4.2, and drone acquisition in site 2 is described in Section 4.3.

Site 1
The acquisition of time-lapse images was carried out on the 16 August 2019 along a test route of 25 m, shown in Table 1, at site 1, on the southernmost sector of a seawall, shown in Figure 1B. The acquisition area is located in a very sheltered bay. The seawall is composed of regular sandstone blocks ranging in length from 0.10 m to 0.5 m.
Remote Sens. 2020, 12, 1880 9 of 21 The sea and weather conditions were perfect at the time the survey was undertaken, shown in Table 2, with a very calm sea, this being a sheltered site. The raft was very stable during the surveying activities, without significant rolling and pitching. Moreover, the good sea conditions allowed the raft to be pushed parallel to the coastline. A total of 511 images were collected, 298 above the waterline and 213 below it at a rate of one per second, shown in Table 2. Table 2 reports all the acquisition parameters set for this test. The emerged and underwater acquisitions were conducted using the same parameters and following the same route, but a GoPro Hero 5 was used awl and a GoPro Hero 6 was used bwl, shown in Figure 7.

Site 1
The acquisition of time-lapse images was carried out on the 16 th August 2019 along a test route of 25 m, shown in Table 1, at site 1, on the southernmost sector of a seawall, shown in Figure 1B. The acquisition area is located in a very sheltered bay. The seawall is composed of regular sandstone blocks ranging in length from 0.10 m to 0.5 m.
The sea and weather conditions were perfect at the time the survey was undertaken, shown in Table  2, with a very calm sea, this being a sheltered site. The raft was very stable during the surveying activities, without significant rolling and pitching. Moreover, the good sea conditions allowed the raft to be pushed parallel to the coastline. A total of 511 images were collected, 298 above the waterline and 213 below it at a rate of one per second, shown in Table 2. Table 2 reports all the acquisition parameters set for this test. The emerged and underwater acquisitions were conducted using the same parameters and following the same route, but a GoPro Hero 5 was used awl and a GoPro Hero 6 was used bwl, shown in Figure 7.

Site 2
The acquisition of time-lapse images was carried out on the 21 June 2019 on a route of 150 m, shown in Figure 4, along plunging cliffs in the southeasternmost sector of the Premantura peninsula (Croatia), east of the entrance to a small sea cave, shown in Figure 1C. The acquisition area is exposed to southerly winds. The plunging cliffs are composed of horizontal decametric limestone beds.
The sea conditions were calm during the surveying activities, being a sheltered site, shown in Figure 1, with a weak SE wind, shown in Table 1, so the raft was stable, without significant rolling and pitching during the surveying activities. The good sea conditions allowed the raft to be pushed parallel to the coastline. However, some submerged blocks and outcrops prevented the perfect following of a route parallel to the coastline in the western sector.
Field acquisition was carried out using three cameras: one GoPro Hero 5 was set awl while two convergent GoPro Hero 6s were set bwl. In total, 715 images were collected above the waterline and 1453 from each camera below it, shown in Table 3. Images were georeferenced and localized on the map, using the integrated GPS tags of the camera awl, shown in Figure 8.

UAV Imagery at Site 2
The acquisition of time-lapse images with UAV was carried out on 20 February 2020 at site 2, shown in Figure 1C, using the parameters of acquisition described in Table 4. The acquisition points with the drone at site 2 are shown in Figure 9, together with the dense cloud model.

UAV Imagery at Site 2
The acquisition of time-lapse images with UAV was carried out on 20 th February 2020 at site 2, shown in Figure 1C, using the parameters of acquisition described in Table 4. The acquisition points with the drone at site 2 are shown in Figure 9, together with the dense cloud model.
The weather was sunny and the sea very calm, as shown in Table 1. Table 4. Acquisition parameters with UAV at Premantura (Croatia).

Data Post-Processing and Production of 3D Models
Time-lapse images and UAV images were used to produce 3D models of the intertidal and nearshore zones. The following sections show the models obtained and describe the analysis parameters of each case study at sites 1 and 2.

D Model at Site 1
A total of 413,025 TPs were obtained from 298 images, using high-quality analysis, shown in Table  5. The alignment of the emergent and submerged parts of the seawall was done using the groups of markers previously described.  The weather was sunny and the sea very calm, as shown in Table 1.

Data Post-Processing and Production of 3D Models
Time-lapse images and UAV images were used to produce 3D models of the intertidal and nearshore zones. The following sections show the models obtained and describe the analysis parameters of each case study at sites 1 and 2.

3D Model at Site 1
A total of 413,025 TPs were obtained from 298 images, using high-quality analysis, shown in Table 5. The alignment of the emergent and submerged parts of the seawall was done using the groups of markers previously described.  Figure 10C,D show that the sandstone blocks are easily distinguishable awl. The dark brown color refers to the high intertidal band, while above it is the splash zone. Intertidal organisms can also be distinguished, in particular limpets, such as Patella Cerulea, and oysters. The resolution of the bwl model, shown in Figure 10E,F, also permits us to distinguish the organisms. In Figure 10F, at the bottom right, brown algae are visible. The sandstone blocks are less distinguishable bwl because they are partially covered by organisms. The awl/bwl transition is influenced by the intersection of the respective models and by wave movements.
Remote Sens. 2020, 12, 1880 13 of 23 Figure 10C,D show that the sandstone blocks are easily distinguishable awl. The dark brown color refers to the high intertidal band, while above it is the splash zone. Intertidal organisms can also be distinguished, in particular limpets, such as Patella Cerulea, and oysters. The resolution of the bwl model, shown in Figure 10E,F, also permits us to distinguish the organisms. In Figure 10F, at the bottom right, brown algae are visible. The sandstone blocks are less distinguishable bwl because they are partially covered by organisms. The awl/bwl transition is influenced by the intersection of the respective models and by wave movements.

3D Model at Site 2 with Horizontal and Drone Images
At Kamenjak, Premantura, the survey was carried out with good sea and weather conditions, as testified by the weather data, shown in Table 1, and therefore the boundary conditions of the test were excellent. Table 6 summarizes the total number of TPs and points in the Dense Cloud. Due to the large number of images, the reference preselection function was used during the alignment process. This forced the software to match only nearby images, therefore reducing the time required for alignment. The distortion decreased using the top and bottom of the iron bar used to define the z axis. The awl model was positioned using the GPS tags of the GoPRO images, while the bwl model was scaled and aligned to the awl model using prominent landforms along the coastline.
Part of the jointed model, obtained from the alignment of independent emergent and submerged models, is shown in Figure 11. The matching process was difficult, mainly because images that lay at sea level did not perfectly coincide with each other.

D Model at Site 2 with Horizontal and Drone Images
At Kamenjak, Premantura, the survey was carried out with good sea and weather conditions, as testified by the weather data, shown in Table 1, and therefore the boundary conditions of the test were excellent. Table 6 summarizes the total number of TPs and points in the Dense Cloud. Due to the large number of images, the reference preselection function was used during the alignment process. This forced the software to match only nearby images, therefore reducing the time required for alignment. The distortion decreased using the top and bottom of the iron bar used to define the z axis. The awl model was positioned using the GPS tags of the GoPRO images, while the bwl model was scaled and aligned to the awl model using prominent landforms along the coastline.
Part of the jointed model, obtained from the alignment of independent emergent and submerged models, is shown in Figure 11. The matching process was difficult, mainly because images that lay at sea level did not perfectly coincide with each other.  The final drawing up of the 3D model of the seawall is shown in Figure 11A,B. A zoom of the textured model, shown in Figure 11C-F, shows the high-resolution and the color range quality Figure 11. shows a sector of the plunging cliffs studied at site 2. The morphology of the cliff is tabular and indented following the horizontal limestone beds. Color belts of biological origin, well described by Schneider [41], show the brownish belt, the white upper intertidal belt, and the grayish supratidal belt. Color belts fit well with the trend of the geological beds, shown in Figure 11A, that range in width from 2 cm to 15 cm. Awl, the surface of the cliff is bare, with very few organisms (such as gastropods, etc.), while bwl, limestone beds are almost completely covered by organisms. Bwl, the cliff is less indented than awl because organisms cover the limestone beds. Joints, fractures, and rock-mass subdivision are easily recognizable, and together with geological beds, characterize the surface morphology of the cliff. The latter is vertically plunging into the sea, while the tidal notch is lacking, shown in Figure 11B, neither above nor below the sea level. A roughly 10 cm-indentation affects the mean sea level, probably due to the selective abrasion of limestone beds. This confirms the surveys of Furlani et al. [6], as the authors have not found tidal notches in southern Istria. In Figure 11C, a white fresh surface is clearly visible at the top center of the photo and this is due to recent wave abrasion. These white surfaces are common at site 2 and the surroundings [25]. The same considerations go for the remaining part of the studied cliffs.
However, the drone model does not allow the modeling of the submerged part of the coast and the resolution is lower than that obtained with horizontal images, shown in Table 7. In Figure 12, we report the distribution of distances between points of UAV-DP images and horizontal time-lapse images. Errors are commonly lower than 5 cm, and most of them are lower than 2.5 cm. The textured model with drone images is very scenic, but the resolution is lower than the model obtained with horizontal images, in particular close to the water level, shown in Figures 13 and 14. The textured model with drone images is very scenic, but the resolution is lower than the model obtained with horizontal images, in particular close to the water level, shown in Figures 13 and 14. The textured model with drone images is very scenic, but the resolution is lower than the model obtained with horizontal images, in particular close to the water level, shown in Figures 13 and 14.

Discussion
Rocky coasts, and in particular plunging cliffs, represent a reliable environment to test photogrammetric procedures to build 3D models of intertidal and nearshore zones. Plunging cliffs are particularly suitable to test the alignment of coastal models because cliff walls are often vertical and almost homogeneous from a geomorphological viewpoint both awl and bwl. Moreover, plunging cliffs are particularly interesting for the recurrent presence of sea level markers developing in the intertidal zone, such as marine notches and other specific landforms [42,43] that would be useful to model. Studies on rocky coasts would greatly benefit from 3D reconstructions, in which the submerged and emergent parts of the intertidal zone can be matched in a single model. Moreover, this work illustrates the possibility of using images from the Geoswim archive [5], collected using low-cost equipment with rapid acquisition, to reconstruct the intertidal morphology of long stretches of rocky coastlines.
The SfM technique has been usually restricted to the independent reconstruction of emergent or underwater objects, respectively through drone/terrestrial photogrammetry or submerged photogrammetric surveys. The resulting 3D representations are very useful for geologists and geomorphologists, as well as archaeologists, as they provide a model to interpret the morphometric parameters of coastal features, such as coastal landforms, geological bedding, etc. [19]. Vertical sea cliffs in submerged environments have also been successfully studied using photogrammetry by combining sideways multibeam echosounder (MBES) and remotely operated vehicle (ROV) videos [44], but

Discussion
Rocky coasts, and in particular plunging cliffs, represent a reliable environment to test photogrammetric procedures to build 3D models of intertidal and nearshore zones. Plunging cliffs are particularly suitable to test the alignment of coastal models because cliff walls are often vertical and almost homogeneous from a geomorphological viewpoint both awl and bwl. Moreover, plunging cliffs are particularly interesting for the recurrent presence of sea level markers developing in the intertidal zone, such as marine notches and other specific landforms [42,43] that would be useful to model. Studies on rocky coasts would greatly benefit from 3D reconstructions, in which the submerged and emergent parts of the intertidal zone can be matched in a single model. Moreover, this work illustrates the possibility of using images from the Geoswim archive [5], collected using low-cost equipment with rapid acquisition, to reconstruct the intertidal morphology of long stretches of rocky coastlines.
The SfM technique has been usually restricted to the independent reconstruction of emergent or underwater objects, respectively through drone/terrestrial photogrammetry or submerged photogrammetric surveys. The resulting 3D representations are very useful for geologists and geomorphologists, as well as archaeologists, as they provide a model to interpret the morphometric parameters of coastal features, such as coastal landforms, geological bedding, etc. [19]. Vertical sea cliffs in submerged environments have also been successfully studied using photogrammetry by combining sideways multibeam echosounder (MBES) and remotely operated vehicle (ROV) videos [44], but intertidal and nearshore zones remain a neglected set. In fact, as aforementioned, the transition zone has rarely been investigated because of inherent differences in the optical settings of the cameras that preclude the common processing of images collected awl and bwl. Recently, [11,14], with very good results, modelled several natural and human-made structures lying or floating at the waterline. Their method can be very useful and reliable for small objects, or short sectors of coast. Moreover, such a task can be accomplished in a relatively short time and with a high degree of precision. The availability of a large number of time-lapse images collected awl and bwl, collected within the Geoswim program [5], directed our efforts at testing the matching of independent models obtained from images collected awl and bwl avoiding the use of metal bars or other targets to perform the matching. This is particularly important in the case of surveys conducted within campaigns along large stretches of coast. In these cases, the time available for the collection of images, field calibration, and target spreading across the study area is very limited. On the contrary, the Geoswim program provides a large number of time-lapse images, which are collected by simply swimming parallel to the coastline. Thus, these images can be used to reproduce the intertidal and nearshore zone and can help provide information for several coastal frameworks.
Photogrammetric models are affected by several types of errors [34], such as the distance of acquisition, and consequently, the variability of the scale of acquisition, camera calibration, image network geometry, surface texture and lighting conditions, ground control point (GCP) characteristics, height of the camera awl and bwl, contrast, light reflection, and photos discarded due to water drops or waves over the dome. Taking all these errors into consideration, we tested the possibility of using these models for rocky coast studies. Figures 10 and 11 show that the alignment fits well, at least when models involve limited sectors of coastline. In this case, the resolution can be up to centimetric. The models obtained in our tests proved to be complete, fast, and quite easy to obtain, and aligned above and below the water interface.
The test performed at site 1, shows that the resolution of the resulting model is sufficient to measure the size of the sandstone blocks on the wall, to count marine shells, and recognize the species of organisms ranging in size from 1 to 2 cm, shown in Figure 10A,B. The alignment was done through rotating the models by matching prominent landforms recognized on the coast, such as two pipes that are recognizable on the left, shown in Figure 11B. Moreover, the colors relative to different heights of the intertidal zone [43] can easily be distinguished and can be used for ecological studies.
The test performed at site 2 was carried out using both horizontal images collected from the raft and 45 • tilted images collected from the drone, along a 150 m route, as shown in the jointed model in Figure 13. The model obtained from horizontal time-lapse images, shown in Figure 11 defines a good resolution for geomorphological studies from the centimetric scale up.
The comparison between the resolution of the models, shown in Table 7, shows that the resolution for bwl models is higher (0.846 mm/pixel) than awl models (1.990 mm/pixel), but the resolution of the drone model is significantly lower (7.000 mm/pixel), because images were collected further away from the cliff. Bwl models can be significantly affected by turbidity, which reduces the resolution and significantly changes the colors of the texturized model, as shown in Figure 10. Moreover, as the morphological complexity of the cliff increases, the models can produce data gaps, shown in Figure 11A, which can reduce the reproduced rock surface. However, this problem affects both drone and horizontal time-lapse acquisitions. For this reason, it is fundamental to collect images with very good sea and weather conditions. The visual comparison of the models, shown in Figure 14, has shown that, using the same camera resolution, awl models are more detailed than drone models, because of the higher distance of the drone camera from the intertidal zone.
The swimming survey was carried out a few meters from the cliff, while the UAV surveys were significantly further away. In particular, images collected from the swim surveys are more detailed on vertical surfaces, because of an inherent difficulty in collecting images with drones from vertical surfaces. Conversely, UAV-DP acquisition is very reliable in modelling horizontal, or slightly inclined rock surfaces, shown in Figure 15A. For example, swim survey images cannot model the top of the stack, shown in Figure 15B, but the model must be integrated with UAV-DP images. What is more, underwater, common commercial drone modelling is completely unreliable, while horizontal images permit the collection of images, even if with a reduced resolution due to the lower visibility bwl. The differences in light conditions of the two models are due to the software analysis of images collected during two surveying days. However, this fact does not affect the resolution of the 3D models. To summarize the pros and cons of the discussed approach, we present a SWOT matrix, shown in Table 8. Swim surveys are particularly effective when modelling plunging cliffs, in particular the intertidal and nearshore zone, when awl and bwl the general morphology is homogeneous, shown in Table 8.
Conversely, the acquisition of horizontal images on low-lying rocky coasts or shore platforms is completely impossible, in particular in the submerged sections [5]. Unfortunately, the collection of images in the latter segment of sea cliffs can be strongly influenced by sea conditions. High turbidity can significantly reduce the number of reliable images, because it reduces the resolution of underwater photos creating shaded images with low visibility. At its worst, with very rough seas, turbidity can completely prevent the building of 3D models. Low light levels also affect the quality of images, because part of the coastline may be lacking image coverage, therefore, high values of local ruggedness decrease the feasibility of building effective models. Furthermore, the sectors of the coast without direct light are all but impossible to model using SfM, because the software cannot align the models. The same goes for backlit images. Local topographic obstacles along the swimming direction may force the operator to To summarize the pros and cons of the discussed approach, we present a SWOT matrix, shown in Table 8. Swim surveys are particularly effective when modelling plunging cliffs, in particular the intertidal and nearshore zone, when awl and bwl the general morphology is homogeneous, shown in Table 8. Table 8. The SWOT matrix. This analysis of the strengths allows for the definition of the characteristics of the method that give an advantage over other approaches of coastal modelling; the weaknesses allow for the definition of the characteristics of the approach that need improvements; the opportunities define the goals, or the elements that the project could exploit to its advantage, while the threats analyze the overall obstacles to using the method.

HELPFUL
Strengths: Opportunities: • The method represents a useful way to model rocky coast topography in microtidal environments, using low-cost instruments; • The acquisition and elaboration times are relatively short; • The 3D models can be improved with subsequent surveys in different time periods, or by using a drone; • The method permits the study of wide sections of coastline; • Modelization processes need to build independent underwater and submerged photogrammetry.
• Data collected can implement a database of rocky coast images; • The models are created in scale and can be manipulated later; • The models can be compared to each other over time; • A wide range of applications in the field of marine environments, such as coastal erosion processes, biological and ecological conservation, and the survey of human-made intertidal structures, such as artificial quays and port structures, are suited to this approach.

Weakness:
Threats: • Images collected in shaded environments, can be strongly affected by a lack of light; • Images are strongly affected by the sun's lighting; • The alignment process could require several steps of pre-process elaboration, such as color correction, etc.; • Many kilometers of rocky coast are surveyed through rapid campaigns, so image acquisition must be optimized in relation to other data; • Low-lying rock coasts are difficult to survey using this approach, because cameras are set horizontally. In this case, UAV-DP outputs might be necessary to add complementary information.
• Data collected during bad sea conditions may prove unreliable for the SfM method, in particular, in the intertidal rather than in the nearshore zone; • Raft roll due to waves can produce blurry and incomplete images; • Images collected bwl can be affected by turbidity, which can reduce or prevent modelling processes; • The pattern of coastal landforms, and details of such structures, might be recognizable only at a certain distance from the coast, which should be maintained during the entire survey.
Conversely, the acquisition of horizontal images on low-lying rocky coasts or shore platforms is completely impossible, in particular in the submerged sections [5]. Unfortunately, the collection of images in the latter segment of sea cliffs can be strongly influenced by sea conditions. High turbidity can significantly reduce the number of reliable images, because it reduces the resolution of underwater photos creating shaded images with low visibility. At its worst, with very rough seas, turbidity can completely prevent the building of 3D models. Low light levels also affect the quality of images, because part of the coastline may be lacking image coverage, therefore, high values of local ruggedness decrease the feasibility of building effective models. Furthermore, the sectors of the coast without direct light are all but impossible to model using SfM, because the software cannot align the models. The same goes for backlit images. Local topographic obstacles along the swimming direction may force the operator to instantly change the image direction of acquisition, in order to achieve the best coverage of the coastline under study, but the final resolution of the model will be reduced. Alongside problems concerning image quality, even a lack of GPS signal may reduce the accuracy of the model, in particular within sea caves.
3D models of the intertidal and nearshore zone are expected to have a wide range of applications in the fields of Earth Sciences, such as geomorphology, geology, biology, ecology, as well as surveys of harbors and other human-made structures.

Conclusions
The approach discussed here suggests that the intertidal and nearshore zones can be successfully described by 3D models obtained using time-lapse images, following precise procedures and under good weather and sea conditions. The merged models allow for the analysis of morphological, but also morphometric features of the intertidal and nearshore zone, where several coastal landforms, including intertidal notches, biogenic or coastal karst holes, furrows and channels, joints, faults, potholes, blocks or boulders of marine or subaerial origin, stacks and sea caves, etc., occur.
Independent models of the emergent and submerged part can be generated using common SfM software and then merged using prominent landforms located on sea cliffs along the surveyed route.
The tests highlighted that the generation of 3D models of the intertidal zone are strongly conditioned by several parameters, such as light transmission in and between air and water, weather and sea conditions, the rapid changes in the water level due to waves, as well as by water turbidity, which may reduce or even prevent 3D modelling. To improve the accuracy of the models, further photographic campaigns in different tidal heights, such as high and low water, can be performed in order to supplement data at the waterline. The comparison with UAV-DP outputs showed that drones permit the mapping of large stretches of coastline and are very reliable for horizontal and sub-horizontal rock surfaces, but the resolution and detail of 3D models collected from swim surveys are higher in the intertidal and nearshore zone. In situations of an abundance of flat surfaces near the water level, it is necessary to integrate the models with UAV-DP outputs. Funding: Funding for this study was provided by the University of Trieste for field activities and instruments. The Falck family partially funded its publication.