Seven Good Reasons for Integrating Terrestrial and Marine Spatial Datasets in Changing Environments

A comprehensive understanding of environmental changes taking place in coastal regions relies on accurate integration of both terrestrial and submerged geo-environmental datasets. However, this practice is hardly implemented because of the high (or even prohibitive) survey costs required for submerged areas and the frequent low accessibility of shallow areas. In addition, geoscientists are used to working on land or at sea independently, making the integration even more challenging. Undoubtedly new methods and techniques of offshore investigation adopted over the last 50 years and the latest advances in computer vision have played a crucial role in allowing a seamless combination of terrestrial and marine data. Although efforts towards an innovative integration of geo-environmental data from above to underwater are still in their infancy, we have identified seven topics for which this integration could be of tremendous benefit for environmental research: (1) geomorphological mapping; (2) Late-Quaternary changes of coastal landscapes; (3) geoarchaeology; (4) geoheritage and geodiversity; (5) geohazards; (6) marine and landscape ecology; and (7) coastal planning and management. Our review indicates that the realization of seamless DTMs appears to be the basic condition to operate a comprehensive integration of marine and terrestrial data sets, so far exhaustively achieved in very few case studies. Technology and interdisciplinarity will be therefore critical for the development of a holistic approach to understand our changing environments and design appropriate management measures accordingly.


Introduction
Ongoing climate changes are producing remarkable impacts worldwide [1]. From the melting of polar ice sheets, and the consequent sea-level rise, to the increasing occurrence of extreme weather events [2], climate changes are notably modifying and/or threatening Earth's environment and ecosystems. Considerable effects have been observed especially in coastal regions [3], which host more than 10% of global population [4]. Coastal erosion, flood risk, increased landslide occurrence and wetland loss are expected to intensify in the coming decades [2], posing serious threats for inhabited areas and environmental assets. A major issue is the uncertainty about the extent and timing of climate-driven impacts, which has often reflected in a "non-immediate" adoption of effective prevention measures and in a non-sustainable coastal management, producing negative socio-economic consequences [2,5].
Coastal zones are the interface between the terrestrial and marine environments and together with nearshore ones are affected by multiple physical processes that originate in both the terrestrial (i.e., fluvial) and marine (i.e., waves and tides) environments. These processes drive geomorphic change, which determines increasingly hazardous conditions (e.g., unstable cliffs, lowlands susceptible to floods) in coastal areas in relation to climate changes and economic development [3,4]. A comprehensive understanding of environmental changes taking place on coastal regions, whose ecosystems are highly productive and very susceptible to changing environmental conditions, relies on an accurate integration of both terrestrial and submerged geo-environmental datasets. This practice is often lacking in coastal management and that still need to be addressed in many regions of the world, where climate change, rising sea levels, tectonic and marine geohazard of different nature are pressing harder year by year and resources need to be more carefully managed (e.g., [6]).
The absence of a seamless spatial data framework prevents in particular a standard procedure for locating and referencing spatial data across the land-marine interface. Different accessibility and investigation costs are the main causes for the huge disparity in size and quality among terrestrial and marine spatial datasets used in environmental research, especially in geomorphological investigation. Landscapes and landforms of terrestrial and marine areas have been traditionally investigated separately, and scientific communities used to working on land-coastal zone included-or at sea independently. On land, geomorphological mapping has been extensively used as a primary method to visualize and analyze Earth surface features ever since early geomorphological research. Taking advantage of the cartographical potentials of Geographic Information Systems (GIS) and the increasing availability of remote sensing tools, data, and products-especially high-resolution aerial and satellite imagery and derived datasets such as Digital Elevation Models (DEMs) and Digital Terrain Models (DTM)-geomorphological maps have been largely produced for many emerged sectors of Earth's surface. Besides being crucial for georisk assessment and land management, detailed geomorphological maps provided reference data for a variety of applied sectors of environmental research, such as landscape ecology, forestry or soil science and spatial planning [7,8]. In the terrestrial domain, the need to share and integrate spatial data for more efficient resource information management has been recognized for over a decade. On the contrary, although the latest developments in submarine acoustic remote sensing [9] have offered increasingly detailed DEMs for the submerged domain of the Earth's surface, only less than 18% of the world's seafloor has been surveyed with a resolution of 30 arc-second, while less than 9% of world's seafloor has been mapped with high-resolution multi-beam sonar [10][11][12][13][14][15]. Undoubtedly, new technologies are revealing more of the seafloor than ever before. The recent development of marine robotics, along with the critical improvement in optical underwater imaging systems (among which the use of hyperspectral cameras in the submarine realm, [16]) and computer vision made it possible to obtain, for the first time, 3D optical reconstruction and a real perception of underwater environments. This allowed for the first time to make more accessible the scientific understanding of submarine processes and environments, not only to scientists but even to the community, with a renewed appreciation of their heterogeneity. Nevertheless, a comprehensive characterization and categorization of submarine landforms and processes is still in its infancy, notwithstanding the important applications for the whole spectrum of submarine geomorphological research [9].
The attention recently paid to the integration of terrestrial and marine spatial datasets for different purposes is notable, although the generation of seamless DEMs and DTMs, based on a reliable integration of onshore, nearshore, and offshore data, is often a challenge for the most part of the coastal regions, especially due to technical issues often associated with the integration of multisource data (i.e., differences in resolution, precision, and accuracy among the datasets available). Addressing this process relies on the development of innovative approaches in using methods and techniques traditionally developed for the investigation of terrestrial environments, for exploring and imaging submarine landforms. A variety of challenging and promising techniques capable of providing a homogeneous and continuum representation of the Earth's surface from the land down to the deep Water 2020, 12, 2221 3 of 35 seafloor has been recently tested. They basically rely on the collection of high-quality data for the integration of multiscale elevation datasets coming from different sources (i.e., bathymetric Light Detecting And Ranging (LiDAR), photogrammetry, and echo-sounders- Figure 1). The application of photogrammetry based on the Structure from Motion technique (using a variety of platforms such as scuba diving, uncrewed surface vehicles-USV-and/or uncrewed aerial vehicles-UAV) has gained in particular new attention as a valuable tool for obtaining bathymetric measurements in very shallow environment, where acoustic devices cannot be safely employed, demanding high costs as in the case of LiDAR surveys [17]. The intertidal and nearshore zones are indeed of critical importance to be surveyed to bridge the gap between the land and the sea and to obtain a reliable integration of marine and terrestrial spatial dataset ( Figure 1).  1). The application of photogrammetry based on the Structure from Motion technique (using a variety of platforms such as scuba diving, uncrewed surface vehicles-USV-and/or uncrewed aerial vehicles-UAV) has gained in particular new attention as a valuable tool for obtaining bathymetric measurements in very shallow environment, where acoustic devices cannot be safely employed, demanding high costs as in the case of LiDAR surveys [17]. The intertidal and nearshore zones are indeed of critical importance to be surveyed to bridge the gap between the land and the sea and to obtain a reliable integration of marine and terrestrial spatial dataset ( Figure 1). This paper aims at showing the reasons why the integration of terrestrial and marine dataset is beneficial for (i) an improved understanding of coastal processes and landforms, (ii) a more effective assessment of coastal risks and impacts, (iii) a sustainable management of coastal environments and associated resources. Progress in the fields of marine acoustic and underwater optical imaging that are contributing to obtain a seamless 3D reconstruction from the land to the sea, are also discussed. Figure 1. Different technological solutions for collecting high-resolution terrestrial and marine data exploiting different platforms: satellite imagery, airborne Light Detecting And Ranging (LiDAR) and bathymetric LiDAR, Uncrewed Aerial Vehicles (UAV)-mounted instruments, and keel-mounted and ROV (Remotely Operated Vehicle)-mounted echo-sounders to obtain high-resolution imagery and derived datasets such as DEMs and DTMs. Note the comparison of different coverages between traditional bathymetric survey vessels (equipped with acoustic systems, such as MBES-MultiBeam EchoSounders) and airborne optical sensors (i.e., bathymetric LiDAR), both in deep and shallow water. The advantages in time and costs of using the latter is obvious, although LiDAR systems provide a lower spatial sampling rate that generates much lower accuracy and precision than that offered by acoustic systems, especially when approaching extinction depths (see Section 3.2 and Table  1 for further details on pros and cons).

Why Integrating Terrestrial and Marine Datasets?
The generation of a seamless digital terrain model that spans the terrestrial and marine environments is a key process for the mapping, modelling, and forecasting of the impact climatedriven changes to coastal geomorphic processes and environmental responses that may occur. This process has in particular the potential of making achievable a more integrated and holistic approach to the management of the coastal zone and the establishment of a timely disaster response and Figure 1. Different technological solutions for collecting high-resolution terrestrial and marine data exploiting different platforms: satellite imagery, airborne Light Detecting And Ranging (LiDAR) and bathymetric LiDAR, Uncrewed Aerial Vehicles (UAV)-mounted instruments, and keel-mounted and ROV (Remotely Operated Vehicle)-mounted echo-sounders to obtain high-resolution imagery and derived datasets such as DEMs and DTMs. Note the comparison of different coverages between traditional bathymetric survey vessels (equipped with acoustic systems, such as MBES-MultiBeam EchoSounders) and airborne optical sensors (i.e., bathymetric LiDAR), both in deep and shallow water. The advantages in time and costs of using the latter is obvious, although LiDAR systems provide a lower spatial sampling rate that generates much lower accuracy and precision than that offered by acoustic systems, especially when approaching extinction depths (see Section 3.2 and Table 1 for further details on pros and cons). of mitigation measures, especially for those climatically sensitive coastal areas where the changing climate is already severely pressing coastal population and their economy.
Hereafter are seven main distinct reasons why a wider integration of land and sea datasets is of great benefit for the associated implications in the context of global environmental changes. They refer to geomorphological research, and to all those applied sectors of environmental research (e.g., marine spatial planning, archaeology and landscape ecology) that focus on spatial planning and sustainability (Table 1 and Figure 2). 1. Development in the field of geomorphological mapping and coastal morphodynamics thanks to innovative mapping techniques and products for coastal and nearshore environments; 2.
Deeper understanding of Late-Quaternary changes of coastal landscapes and environments, with positive implications for prediction of future risk scenarios; More comprehensive recognition and assessment of geoheritage and geodiversity; 5.
Wider assessment of coastal geohazards and vulnerability in the frame of disaster risk reduction, thanks to the development of models of different processes (e.g., coastal hydrodynamic modelling); 6.
Critical support to other disciplines involved in generating key-data for the sustainable management of marine resources, such as marine and landscape ecology; 7.
Establishment of more sustainable development objectives in planning and management of coastal areas, also in the framework of Integrated Coastal Zone Management.
In detail, integrated research merging terrestrial and submarine data can substantially contribute to: 1. Development in the field of geomorphological mapping and coastal morphodynamics thanks to innovative mapping techniques and products for coastal and nearshore environments; 2. Deeper understanding of Late-Quaternary changes of coastal landscapes and environments, with positive implications for prediction of future risk scenarios; 3. New approaches for geoarchaeological research development in coastal and nearshore environments; 4. More comprehensive recognition and assessment of geoheritage and geodiversity; 5. Wider assessment of coastal geohazards and vulnerability in the frame of disaster risk reduction, thanks to the development of models of different processes (e.g., coastal hydrodynamic modelling); 6. Critical support to other disciplines involved in generating key-data for the sustainable management of marine resources, such as marine and landscape ecology; 7. Establishment of more sustainable development objectives in planning and management of coastal areas, also in the framework of Integrated Coastal Zone Management.  Table 1.  Table 1.

Geomorphological Mapping
Analyzing and representing terrestrial and marine data in an integrated perspective allows us to provide a comprehensive picture of a coastal landscape, taking into account the variety of processes that contribute to landform development through time. The most appropriate tool for the representation of emerged and submerged landscapes and landforms in coastal areas (or, in general, in transition environments) is an integrated geomorphological map that relies on the availability of elevation data from the onshore to the offshore settings. Scientific literature started to report consistent examples toward this effort only during the last decade (e.g., [18][19][20]).
A traditional (terrestrial) geomorphological map generally represents landform genesis, distribution, morphometry, and state of activity, offering a dynamic vision of the investigated coastal landscape (e.g., [23,152,153]). According to the mapping approach used, information on lithology, geological structure, and age of surficial deposits can be included. The value of such maps is unquestionable for both scientific and applied research, aiming at environmental conservation and management, and hazard and risks assessment. The same applies to marine geomorphological maps [9,154].
The main differences between terrestrial and marine geomorphological maps and the hurdle in producing geomorphological maps integrating the representation of land and seafloor features in coastal areas and shallow waters concern: • Representation scale: terrestrial geomorphological maps are more easily drawn at fine scale (e.g., 1:5000) than submarine geomorphological maps, because of the higher ease of data collection on land than underwater; • Standardized terminology and classification schemes: terrestrial landforms are codified at an international level, and their definition and representation are generally shared, either worldwide or country-wide, while for submarine landforms-apart from the main physiographic classification [155]-different terms and classification schemes are used (e.g., see how bedforms  [157]; Wynn and Stow [158]; Stow et al. [159]); • Standardized symbology: standard symbols are codified for terrestrial geomorphological mapping, though they may vary from country to country, while they are lacking for submarine geomorphological mapping and for integrated terrestrial and submarine geomorphological mapping; • Coverage of lithological and chronological information: in the terrestrial environment, the coverage of lithological and chronological data is generally much higher than for the seafloor, thanks to the wider availability of maps, DTMs, scientific literature, and datasets; • Acquiring elevation data in marine regions poses significant challenges and some limitations which make the process more complex, from a technological point of view, than in the emerged system. Customized surveys and dedicated technological solutions are indeed required to apply specific corrections that can address all measurements errors created in particular by hydrodynamics (especially tides and wave motion that must be always severely taken into account in hydrographic survey carried out by mean of echo-sounders) and the physical variability of the water column (which has a strong impact on the sound velocity/refraction of beams, creating at places challenging conditions, such as in the case of fresh water influx at the mouth of a river). Cloud coverage, turbidity, water surface glint and breaking waves can also create challenging environmental and operational condition that may prevail during optical remote sensing surveys in shallow water.
Nevertheless, the scientific community has recently paid increasing attention to the above-mentioned issues and different solutions were introduced to overcome them. First of all, technological improvements in acoustic remote sensing has now made it possible to acquire high-resolution elevation data in shallow environments, ranging from 0.5 m [160] to 0.05 m of resolution [28,161], and also with precision and accuracy close to those achievable for terrestrial landscape investigation. In addition, even the deep environment (deeper than 200 m, [162]) can be now surveyed at high resolution thanks to the progress made by underwater robotics. Some compromises are still needed for shallow-water data acquisition, but recent progress has definitely produced unprecedented results (see Section 3.2 and Table 2 for additional details). On the contrary, the standardization of terminology, classification schemes, and symbology is definitely still an unresolved issue, even though many efforts have been made globally to offer practical and applicable solutions in different contexts. In December 2015, a partnership originated from seabed mapping programs of Norway, Ireland and the UK (MAREANO, INFOMAR, and MAREMAP, respectively-MIM) shares experiences, knowledge, expertise and technology in order to propose advancement in best practices for geological seabed mapping [163]. The development of a standardized geomorphological classification approach, following previous work published by the British Geological Survey [164], has been one of the key activities of the MIM group. At the European level, the EMODnet-Geology Portal (European Marine Observation and Data Network; [165]) is collecting data from European state members on the geology and geomorphology of European seabeds and harmonizing the geomorphological representation through a vocabulary developed by the Federal Institute for Geosciences and Natural Resources [166] and by the British Geological Survey [167]. In addition to research strictly dealing with geological and geomorphological mapping, further inputs have been also provided by organizations extensively involved in providing seafloor maps for several environmental purposes, but where precise representation and definition of submarine landforms were crucial for their applications. The Integrated Ocean and Coastal Mapping (IOCM) of the National Oceanic and Atmospheric Administration (NOAA) developed the Coastal and Marine Ecological Classification Standard (CMECS), a scheme for classifying benthic habitats that includes a geomorphic level for the description of submarine landforms. At the European level (within the frame of different EU funded projects, such as Mapping European Seabed Habitats (MESH), CoralFish [168] and CoCoNet [169] among others), various classification schemes were proposed [170,171] taking inspiration from the CMECS and/or from the analogue European classification scheme EUNIS [172]. Table 2. List of sensors and techniques cited in the text which are relevant for collecting elevation data in the shallow-water environment (i.e., white ribbon). Active sensors are those sensors that control the projection of acoustic or electromagnetic waves (i.e., sound or light) onto the scene, while the term "passive" refers to those systems or techniques that use ambient light to illuminate the scene or techniques such as Structure from Motion (SfM A few attempts have recently been made to overcome the critical issue of defining a standard symbology for submarine and integrated land-sea geomorphological mapping, and most of them addressed the issue at the national level. Within the CARG Project (Italian official Geological Cartography), the Institute for the Protection and Environmental Research (ISPRA) promoted the production of Geological and Geothematic Sheets at the scale 1:50,000 of the Italian territory. In particular, in coastal areas, both terrestrial and marine geology and terrestrial and marine geomorphology were represented [173]. Recently, the Working Group on Coastal Morphodynamics of the Italian Association of Physical Geography and Geomorphology-in agreement with ISPRA-revised and updated the legend for the "Geomorphological Map of Italy" [27] with reference to coastal, marine, lagoon and aeolian landforms, processes and deposits [23]. Within the revision process submerged landforms and deposits down to the continental shelf break were included in the new legend. The first applications of this new legend were provided by Brandolini et al. [25], Furlani et al. [26] and Prampolini et al. [6,24]. Brandolini et al. [25] compared multi-temporal maps of anthropogenic landforms in Genoa city (Liguria region, Italy) to analyze the urban development since the Middle Age to provide a tool for geo-hydrological risk assessment. Furlani et al. [26] surveyed and mapped the coastal stretch of the Conero system (Marche region, Italy), focusing on tidal notches as geomorphological markers for analyzing Late-Holocene tectonic stability of the area. Prampolini et al. [6,24] surveyed and mapped coastal and submarine areas of the Maltese archipelago in order to produce the first integrated geomorphological maps of the islands, constituting the basis for geomorphological reconstruction and coastal hazard assessment in an environment highly sensitive to ongoing climate change ( Figure 3).  -Highly favorable ratio in coverage/unit of time (i.e., wide areas covered in a relatively small period of time) -Source data are public and available (i.e., Sentinel-2) -Source data also offer photorealistic view of the surveyed scene

Late-Quaternary Changes of Coastal Landscapes
Scientific literature referring to integrated analyses of terrestrial and marine geospatial datasets related to coastal areas and continental shelves has proved to be critical in outlining Late-Quaternary coastal changes and paleo-landscapes. With reference to the coverage of lithological and chronological information, it is possible to produce marine geological maps exploiting data coming from marine geological and geophysical surveys. Examples from the Italian coastal and marine areas are given by ISPRA and Gasparo Morticelli et al. [22]. ISPRA, thanks to the CARG Project, has promoted the mapping of the marine geology of the Adriatic Sea, scale 1:250,000 [174], to represent the main geological structures of the Adriatic continental shelf. Gasparo Morticelli et al. [22] published the geological map of Marettimo Island and the submerged surroundings (Egadi Islands, central Mediterranean Sea), scale 1:10,000, integrating data from field surveys and previous marine geological and geophysical surveys, carried out within the CARG Project. Moreover, advances are being made in correlating seabed acoustic backscatter with surface sediments to perform automatic ground discrimination and reduce the number of ground truth stations required to produce robust seabed characterization. Recent availability of well-data from the industry (e.g., ViDEPi [175]) and progress in seismic analysis and seismo-stratigraphic correlation, contributed for a wider geological characterization of the subseafloor too. This way, it is possible to get a wider coverage of lithological and chronological information for submarine regions.

Late-Quaternary Changes of Coastal Landscapes
Scientific literature referring to integrated analyses of terrestrial and marine geospatial datasets related to coastal areas and continental shelves has proved to be critical in outlining Late-Quaternary coastal changes and paleo-landscapes.
Quaternary global sea-level changes have been largely studied worldwide with special attention to the late Pleistocene and Last Glacial Maximum (LGM), when the relative sea level reached its minimum height of ca. 120-130 m below the present sea level ( Figure 4) [176,177].  Consequently, most of the present marine areas down to 130 m of depth, including large portions of the world continental shelves, were emerged during the LGM and modelled by subaerial processes: the paleo-geography of emerged and submerged areas was very different from today, and included land bridges between continental landmasses which do not exist anymore (e.g., the land bridge between Sicily and the Maltese Archipelago- Figure 5-or between Sicily and Egadi Islands of Favignana and Levanzo, in southern Europe [39,42], or the Bering land bridge in the northern Pacific Ocean ( Figure 6) [178]).
The post-glacial sea-level rise covered these landscapes, either drowning [33,37] or re-modelling them through marine processes. At some places, the previously emerged areas hosted prehistoric and historic human settlements-as demonstrated by ruins, pollens, bones, or other remains in the sediments, caves, etc. (cf. Section 2.3) [56]. The post-glacial marine transgression led to the disappearance of vast areas, land bridges, human settlements, and a reconfiguration of geographical boundaries. Indeed, the reconstruction of paleo-landscapes and their geomorphological evolution is the basis for further investigations (e.g., on geohazard assessment, identification, and study of geoarchaeological sites). Consequently, most of the present marine areas down to 130 m of depth, including large portions of the world continental shelves, were emerged during the LGM and modelled by subaerial processes: the paleo-geography of emerged and submerged areas was very different from today, and included land bridges between continental landmasses which do not exist anymore (e.g., the land bridge between Sicily and the Maltese Archipelago- Figure 5-or between Sicily and Egadi Islands of Favignana and Levanzo, in southern Europe [39,42], or the Bering land bridge in the northern Pacific Ocean ( Figure 6) [178]).
The post-glacial sea-level rise covered these landscapes, either drowning [33,37] or re-modelling them through marine processes. At some places, the previously emerged areas hosted prehistoric and historic human settlements-as demonstrated by ruins, pollens, bones, or other remains in the sediments, caves, etc. (cf. Section 2.3) [56]. The post-glacial marine transgression led to the disappearance of vast areas, land bridges, human settlements, and a reconfiguration of geographical boundaries. Indeed, the reconstruction of paleo-landscapes and their geomorphological evolution is the basis for further investigations (e.g., on geohazard assessment, identification, and study of geoarchaeological sites).  Numerous examples of terrestrial and marine data integration for the reconstruction of paleolandscapes are from the Mediterranean Basin [38]. Lo Presti et al. [42] reconstructed the paleogeography of Egadi Islands and relative sea-level variations from the LGM until today. Furlani and Martin [40] reconstructed the paleo-geography of Faraglioni coast (Ustica Island) that was settlement of a Middle Bronze Age village. They combined geomorphological observation made in nearshore and onshore areas. Miccadei et al. [18,19] reconstructed the Late-Quaternary landscape and geomorphological evolution of Tremiti Islands, located north of Gargano promontory, southern Adriatic Sea. Aucelli et al. [36] carried out a multidisciplinary study of submerged ruins of Roman buildings on the Sorrento Peninsula coast (Gulf of Naples). These archaeological remains enabled the reconstruction of the ancient position of both the sea level and the coastline. Rovere et al. [31] analyzed the submarine geomorphology of the offshore between Finale Ligure and Vado Ligure (western Liguria, NW Mediterranean Sea) for the first time, detecting meaningful submarine geomorphological indicators of former sea levels. Micallef et al. [33] and Foglini et al. [37]    Numerous examples of terrestrial and marine data integration for the reconstruction of paleolandscapes are from the Mediterranean Basin [38]. Lo Presti et al. [42] reconstructed the paleogeography of Egadi Islands and relative sea-level variations from the LGM until today. Furlani and Martin [40] reconstructed the paleo-geography of Faraglioni coast (Ustica Island) that was settlement of a Middle Bronze Age village. They combined geomorphological observation made in nearshore and onshore areas. Miccadei et al. [18,19] reconstructed the Late-Quaternary landscape and geomorphological evolution of Tremiti Islands, located north of Gargano promontory, southern Adriatic Sea. Aucelli et al. [36] carried out a multidisciplinary study of submerged ruins of Roman buildings on the Sorrento Peninsula coast (Gulf of Naples). These archaeological remains enabled the reconstruction of the ancient position of both the sea level and the coastline. Rovere et al. [31] analyzed the submarine geomorphology of the offshore between Finale Ligure and Vado Ligure (western Liguria, NW Mediterranean Sea) for the first time, detecting meaningful submarine geomorphological indicators of former sea levels. Micallef et al. [33] and Foglini et al. [37]  To reconstruct ancient landscapes located in coastal and shallow-water environments, geomorphological, archaeological, and geophysical observations of both terrestrial and submarine areas are required, together with absolute dating of terrestrial and marine sediments, remains and landforms [54,61]. Hence, technological advances in acquiring data in nearshore areas and in merging terrestrial and marine spatial datasets (e.g., elevation data) are fundamental for the reconstruction of present emerged and submerged topography and to infer paleo-geography.
Numerous examples of terrestrial and marine data integration for the reconstruction of paleo-landscapes are from the Mediterranean Basin [38]. Lo Presti et al. [42] reconstructed the paleo-geography of Egadi Islands and relative sea-level variations from the LGM until today. Furlani and Martin [40] reconstructed the paleo-geography of Faraglioni coast (Ustica Island) that was settlement of a Middle Bronze Age village. They combined geomorphological observation made in nearshore and onshore areas. Miccadei et al. [18,19] reconstructed the Late-Quaternary landscape and geomorphological evolution of Tremiti Islands, located north of Gargano promontory, southern Adriatic Sea. Aucelli et al. [36] carried out a multidisciplinary study of submerged ruins of Roman buildings on the Sorrento Peninsula coast (Gulf of Naples). These archaeological remains enabled the reconstruction of the ancient position of both the sea level and the coastline. Rovere et al. [31] analyzed the submarine geomorphology of the offshore between Finale Ligure and Vado Ligure (western Liguria, NW Mediterranean Sea) for the first time, detecting meaningful submarine geomorphological indicators of former sea levels. Micallef et al. [33] and Foglini et al. [37] reconstructed Late-Quaternary coastal landscape morphology and evolution of the Maltese archipelago, while Furlani et al. [39] focused their research on marine notches of the Maltese Islands that resulted in confirmation of the slowdown of the Late-Holocene marine transgression.
Examples from northern Europe come, among the others, from southern England and northern France [29], and from Northern Ireland [32]. Bridgland et al. [29] analyzed three fluvial sequences, particularly terrace staircases, from southern England and northern France to reconstruct climate fluctuations and paleo-geography of those areas. Westley et al. [32] mapped the continental shelf of northern coast of Ireland and examined the geomorphology for evidence of past sea-level changes, reconstructed the paleo-geography of the area considering sea-level lowstands of −30, −14 and −6 m. This research allowed the identification of ten areas of high archaeological potential.
The combined mapping of emerged and submerged geomorphological features proved to be functional in analyzing the long-term evolution of coastal landslides. Prampolini et al. [6] showed that coastal block slides along the NW coast of Malta prolong below the sea level, reaching a depth of about 40 m, and Soldati et al. [180] demonstrated by means of cosmogenic nuclide dating that they developed in a subaerial environment-when the coastline was much lower than today-having been submerged only later on, during the post-glacial sea-level rise.

Geoarchaeology
As earlier mentioned, Late-Quaternary sea-level changes have exposed large portion of the present-day continental shelves for long periods of time, resulting in a multitude of archaeological remains lying on the seafloor today [4,45,46,48,49,52]. Therefore, coupling terrestrial and marine datasets can be critical in detecting new archaeological sites in coastal and nearshore areas and for a more comprehensive understanding of already existing ones [60] and references therein. Analyzing coastal archaeological sites can also contribute to the reconstruction of the paleo-geography of ancient landscapes [48], and in particular to infer Late-Holocene relative sea-level oscillations (e.g., [43,51]). Some archaeological remains include functional structures or elements that are unequivocally related to specific elevation of past sea levels, because of their architecture and proximity to the sea. In other cases, an in-depth knowledge of landscape evolution helps inferring about the evolution or the dating of archaeological remains. For example, the proximity of an archaeological site with coastal landforms, whose evolution can be reconstructed, will help in dating and reconstructing the history of the site itself [40] and references therein.
As a matter of fact, early human populations tended to move and expand occupying new territories, in particular during the last glaciation and the early post-glacial period-a period of time characterized by extreme climate fluctuations [56]. In this frame, the most attractive sites for human settlements were coastal lowlands that in some parts of the world, were much more extended than today thanks to sea-level lowstands (e.g., during the LGM, the European land area was 40% wider than presently; [56,181]). During that period, coastal regions were the most densely populated since they profited from more tempered climates that led to enhanced water supplies, and greater ecological diversity. Hence, these areas were sites of prehistoric and historic human settlements, as witnessed by the findings of archaeological remains, pollens, bones, or other ruins in the sediments, caves etc. Then post-glacial sea-level rise led to the flooding of these former territories, redrawing geographical boundaries, and human, plant and animal distributions (cf. Section 2.2) [56,181]. In this context, Harff et al. [56] reported the results achieved within the framework of the SPLASHCOS Project-Submerged Prehistoric Archaeology and Landscapes of the Continental Shelf-(Cooperation in Science and Technology-COST Action TD0902). They succeeded in gathering together experts in geology, archaeology, and climate interested in sea-level changes, paleo-climatology and paleo-geography for reconstructing European submerged landscapes in order to assess their archaeological potential (e.g., [37], cf. Section 2.2).
Geoarchaeological research benefiting from the integration of terrestrial and marine datasets is illustrated with reference to the Mediterranean area by Furlani et al. [50], Aucelli et al. [36,55], Furlani and Martin [40] and Mattei et al. [59]. Furlani et al. [50] studied submerged or partially submerged archaeological structures located along the Maltese coasts providing a first attempt for paleo-environmental reconstruction of the Maltese archipelago from the LGM until today, allowing time and mode of mammal dispersal to the island during the Pleistocene to be inferred. Aucelli et al. [36] analyzed submerged ruins of Roman buildings located along the Sorrento Peninsula coast (Italy) and succeeded in reconstructing sea-level oscillations and coastline changes for the Late-Holocene and tectonic history of the Sorrento Peninsula during the last two millennia. Aucelli et al. [55]  Examples of integration of land-sea datasets for geoarchaeological purposes outside the Mediterranean area are provided by Fisher et al. [47] and Cawthra et al. [57] for South Africa, by Bailey et al. [54] for Saudi Arabia, and by Benjamin et al. [58] and Veth et al. [61] for Western Australia. Fisher et al. [47] developed a conceptual tool that enable correlation of the evolution of human behavior within a dynamic model of changes of paleo-environment. Cawthra et al. [57] analyzed paleo-coastal environments, laying on the present continental shelf, offshore of the Pinnacle Point archaeological locality (Mossel Bay, South Africa). During the Pleistocene, these environments were probably settlement of early-modern humans. Bailey et al. [54] analyzed emerged and submerged landscape of SW Saudi Arabia to study human dispersal in Late Pleistocene and Early Holocene. Finally, Benjamin et al. [58] and Veth et al. [61] analyzed a large amount of data (airborne LiDAR, underwater acoustics, cores and scuba dives observations) acquired in Western Australia for archaeological research.
Coupling land-sea data is now more feasible thanks to modern marine research technologies, integration of large databases and proxy data [60] and references therein, allowing further hidden archaeological sites to be discovered and studied in the near future. Recently, this has proved to be successful is the case of the ancient Roman city of Baia located inside the Bay of Pozzuoli and belonging to the Campi Flegrei Archeological Park (Southern Italy). The site is superbly preserved underwater after having been slowly drowned due to bradyseismic movements which characterize this area near the Vesuvius volcano [182].
Geoarchaeological investigations, particularly in coastal and submerged environments are increasing and are taking advantage of new contributions and new approaches in surveying and collecting data using a combination of acoustic and optical remote sensing sources, to recreate a full picture of the present and old landscapes, validated through field surveys observations and absolute dating evidence (e.g., Uncrewed Surface Vehicle simultaneously acquiring geophysical data and images for photogrammetry and drones equipped with cameras).

Geoheritage and Geodiversity
Integrating terrestrial and marine datasets can be of paramount importance for the assessment of terrestrial and marine sites of geological interest in coastal and shallow-water areas. Both the shore and inner continental shelf show common processes and landforms that should be considered to be a single feature (cf. Sections 2.1 and 2.2) [66].
Geosites-or geodiversity sites (sensu Brilha [183])-are places of a certain value due to human perception or exploitation and include geological elements with high scientific, educational, aesthetic, and cultural importance [71]. Geosites and key geodiversity areas are often protected areas thanks to different directives (e.g., EU Habitat Directive, 1992; OSPAR Convention, 1992; EU Marine Strategy Framework Directive, 2008). Indeed, the importance of preserving geodiversity has been acknowledged mainly thanks to the effect that geodiversity has on biodiversity patterns [184]. Hence, the assessment of geosites and key geodiversity areas enhance the identification of areas that need protection (e.g., Marine Protected Areas-MPAs-Geoparks; cf. Sections 2.6 and 2.7).
Geoheritage and geodiversity have been investigated mainly in terrestrial environments [71] and references therein, while only a few studies on this topic refer to underwater environments (cf. [62][63][64]66]). With reference to the latter, Orrù and Ulzega [62], Orrù et al. [63,64] identified underwater trails for scuba divers in the MPAs of Capo Carbonara and of Capo Caccia (Sardinia, Italy) enhancing the value of the whole underwater environment. Rovere et al. [66] assessed underwater geomorphological heritage in the Bergeggi MPA (Ligurian Sea, Italy) and in the Sigri area (Lesvos Island, Greece); however they considered that a complete approach in studies on geoheritage would take both emerged and submerged landforms of the coastal and nearshore environments into account.
Very few studies deal with the integrated assessment of terrestrial and marine sites of geological interest. This is largely due to technological constraints and, to some extent, to conceptual issues-such as (i) differences in attributes related to geosites in terrestrial and marine environments [185]; (ii) different perception and fruition of abiotic features of the aquatic environment by tourists; and (iii) lack of common schemes and approaches to the identification, assessment, and improvement of submarine geosites. However, attention to these themes is increasing presently.
Coratza et al. [71] identified and assessed the terrestrial and marine geosites of the Portofino Natural Park and MPA (Liguria, Italy), which are internationally known for both terrestrial scenic landforms and quality of the marine ecosystem. They aimed at identifying the most suitable sites for tourist improvement and defining possible connections between terrestrial and marine environments. Finally, they were able to identify a significant number of both terrestrial and marine sites, assessing their scientific value, ecological, cultural, and aesthetic importance, and accessibility, services, and economic potential.
The identification and quantitative assessment of geodiversity in terrestrial areas are already established and several countries have been developing national inventories of geodiversity key areas by means of different methods [183]. Only recently, this discipline has been addressed to the marine environment and its elements of geodiversity. Examples of geodiversity assessment in marine areas come from Scotland [65][66][67][68][69][70], southeast Brazil [67], Hawaiian and Canarian Islands and the New Zealand subduction zone [68]. Brooks et al. [65] and Gordon et al. [69,70] focused on the contribution of geo-conservation within the Marine Protected Area network for Scottish seawaters. Their work is the first systematic assessment of marine geodiversity key areas comparable to the Geological Conservation Review geo-conservation carried out for the terrestrial geology and geomorphology of Great Britain, which was "a world-first project of its type in the systematic assessment of the whole geological heritage of a country, from first principles" [69]. Mansini Maia and Alencar Castro [67] developed a model for characterizing marine geodiversity at a regional scale in the Vitória-Trindade Volcanic Seamount Ridge and its surroundings (SE Brazil). They aimed at supporting Brazilian marine spatial planning regarding geo-conservation of features related to the geological history of Brazil and the most vulnerable habitats. Finally, Veloo [68] developed a geodiversity index for the seafloor and applied it to three study areas (Hawaii, Canary Islands, and the New Zealand subduction zone), considering abiotic factors such as geomorphology, bathymetry, range of slope angle, and light penetration.
Although there is a growing number of studies on terrestrial and marine geodiversity showing that its understanding is essential for several issues, including geo-conservation [186], still there is no specific literature available on the integrated assessment of geodiversity in coastal and shallow-water areas. However, available knowledge and technological tools call for immediate actions in this field, which would be highly beneficial for holistic management and planning in coastal and nearshore areas.

Geohazards
Coastal and shallow-water environments are threatened by different kinds of geohazards that can produce significant impacts on (i) the economy, due to possible reduction in tourism and disruption of urbanized areas, (ii) on landscapes, due to possible severe morphological changes both onshore and underwater, and (iii) on ecosystems, due to possible loss of sensitive habitats [187][188][189]. In this context, the availability of terrestrial and marine spatial datasets is fundamental to get a full picture of coastal geohazards. A combined analysis of terrestrial and marine processes should be considered to be a necessary step in geohazard assessment in coastal environments. In the past few years, there has been an increase of published papers in this field of research. Here we briefly present a review of literature with special reference to tsunami and storm waves, volcanic eruptions, coastal landslides, coastal inundation, and erosion due to sea-level rise.
Recent progress in hydrodynamic modelling and simulation produced considerable results in topics such as tsunami and storm wave hazards in coastal environments, mostly because bathymetric data is a crucial parameter in nearshore wave and hydrodynamic modelling [28,190]. De Jongh and van Opstal [95] proposed an interesting combined analysis of topography and bathymetry in Mozambique to model tsunami and storm surges impacts on land. One of the most common effects to storm and tsunami, after the flooding of low-lying areas, is the detachment of boulders of a variety of sizes (from decametric to metric) from the seafloor. The knowledge of both nearshore and coastal geomorphology is fundamental to develop models reconstructing the height of the wave necessary to produce such a detachment and boulders' possible path on land. Examples of tsunami or storm deposits and models on the waves that caused them come from the Mediterranean Sea, in particular from the Istrian coast (e.g., [116]), the Apulian coasts (e.g., [99] and references therein), the Maltese archipelago (e.g., [104,105]) and the Greek coasts (e.g., [81,113]).
De Gange et al. [88] illustrated the effects of volcanic eruptions on coastal and marine environments, such as spreading of volcanic ashes and pyroclasts, which can affect also terrestrial and marine habitats, earthquakes and landslides (Figure 7). which would be highly beneficial for holistic management and planning in coastal and nearshore areas.

Geohazards
Coastal and shallow-water environments are threatened by different kinds of geohazards that can produce significant impacts on (i) the economy, due to possible reduction in tourism and disruption of urbanized areas, (ii) on landscapes, due to possible severe morphological changes both onshore and underwater, and (iii) on ecosystems, due to possible loss of sensitive habitats [187][188][189]. In this context, the availability of terrestrial and marine spatial datasets is fundamental to get a full picture of coastal geohazards. A combined analysis of terrestrial and marine processes should be considered to be a necessary step in geohazard assessment in coastal environments. In the past few years, there has been an increase of published papers in this field of research. Here we briefly present a review of literature with special reference to tsunami and storm waves, volcanic eruptions, coastal landslides, coastal inundation, and erosion due to sea-level rise.
Recent progress in hydrodynamic modelling and simulation produced considerable results in topics such as tsunami and storm wave hazards in coastal environments, mostly because bathymetric data is a crucial parameter in nearshore wave and hydrodynamic modelling [28,190]. De Jongh and van Opstal [95] proposed an interesting combined analysis of topography and bathymetry in Mozambique to model tsunami and storm surges impacts on land. One of the most common effects to storm and tsunami, after the flooding of low-lying areas, is the detachment of boulders of a variety of sizes (from decametric to metric) from the seafloor. The knowledge of both nearshore and coastal geomorphology is fundamental to develop models reconstructing the height of the wave necessary to produce such a detachment and boulders' possible path on land. Examples of tsunami or storm deposits and models on the waves that caused them come from the Mediterranean Sea, in particular from the Istrian coast (e.g., [116]), the Apulian coasts (e.g., [99] and references therein), the Maltese archipelago (e.g., [104,105]) and the Greek coasts (e.g., [81,113]).
De Gange et al. [88] illustrated the effects of volcanic eruptions on coastal and marine environments, such as spreading of volcanic ashes and pyroclasts, which can affect also terrestrial and marine habitats, earthquakes and landslides (Figure 7).  Examples of landslides triggered by volcanic eruptions or, in general volcanic setting, are from Mount Etna (Sicily, Italy) and oceanic and insular volcanoes. Urlaub et al. [115] analyzed the deformations of the southeastern flank of Etna volcano that is sliding into the Ionian Sea and carried out the first long-term seafloor displacement monitoring campaign. Oceanic and insular volcanoes commonly experience giant landslides with relevant run-out (i.e., debris avalanches), able to create huge depositional areas in the offshore and even deep domain. Examples come from Canary Islands [76,77,101], Hawaii Islands [72,75,78], Stromboli Island ( Figure 8) [72,75,78], Stromboli Island ( Figure 8) [83,85,91], Lipari Island [110] etc. As for more common coastal landslides, there are case studies showing how the integration of land and sea datasets can be beneficial for landslide hazard assessment. In the Calabria region, seismic-induced landslides originated on land and reached the seafloor [94,100]. De Blasio and Mazzanti [87] produced a few-centimeters resolution DTLM (Digital Terrain and Lacustrine Model) and a DTMM (Digital Terrain and Marine Model) of two Italian sites in Latium and Calabria affected by coastal rock falls in order to model the falling of material into the water. Casalbore et al. [90] analyzed terrestrial and marine DTMs both pre-and post-hyperpycnal flows at Fiumara (Western Messina Strait, Italy) aiming at detecting the morphological evidence of the event on the seabed and to assess flash flood occurrence a posteriori. Another example of integration of terrestrial and marine datasets for landslide analysis is from the Nice landslide (Ligurian Sea, NW Mediterranean) subsequent to the 1979 catastrophe of the Nice International Airport (NE France) that caused a 2-3 m high tsunami, generated by a landslide that progressively turned into a debris flow and, then, in a turbidity current [73,80]. Several studies analyzed morphology, stratigraphy, geotechnics of the landslide and surrounding terrestrial and marine areas, providing numerical models of the phenomenon [74] and evaluating the possibility of future collapses and related impacts on the environment and human activities [89]. Coastal and marine mass wasting can also be related to past and present sea-level rise due to climate change. This is the case of the Vasto landslide (Abruzzo, Italy), a rotational slide continuing under the sea level, whose geomorphological evolution, and past and historical reactivations have been reconstructed by Della Seta et al. [96]. Another example comes from the Maltese coastal block slides that developed during the last glaciation and were then influenced by the successive sea-level oscillations [180].
Among hazards induced by the ongoing climate change, coastal inundation (especially along stretches of coast affected by subsidence; [108]) and coastal erosion triggered by extreme meteorological events, and sea-level rise are the most reported in the literature. In particular, several papers concern (i) the quantification of sea-level rise (e.g., [2,93,107,111]; Figure 9), (ii) the general impacts of sea-level rise on coastal environments [97,102] and (iii) the assessment of coastal exposure and coastal erosion (e.g., [98,114,117,118]). Terrestrial and marine datasets are differently analyzed and integrated, although bathymetry is often taken into consideration, especially in those works presenting predictive models on sea-level rise and coastal inundation. As for more common coastal landslides, there are case studies showing how the integration of land and sea datasets can be beneficial for landslide hazard assessment. In the Calabria region, seismic-induced landslides originated on land and reached the seafloor [94,100]. De Blasio and Mazzanti [87] produced a few-centimeters resolution DTLM (Digital Terrain and Lacustrine Model) and a DTMM (Digital Terrain and Marine Model) of two Italian sites in Latium and Calabria affected by coastal rock falls in order to model the falling of material into the water. Casalbore et al. [90] analyzed terrestrial and marine DTMs both pre-and post-hyperpycnal flows at Fiumara (Western Messina Strait, Italy) aiming at detecting the morphological evidence of the event on the seabed and to assess flash flood occurrence a posteriori. Another example of integration of terrestrial and marine datasets for landslide analysis is from the Nice landslide (Ligurian Sea, NW Mediterranean) subsequent to the 1979 catastrophe of the Nice International Airport (NE France) that caused a 2-3 m high tsunami, generated by a landslide that progressively turned into a debris flow and, then, in a turbidity current [73,80]. Several studies analyzed morphology, stratigraphy, geotechnics of the landslide and surrounding terrestrial and marine areas, providing numerical models of the phenomenon [74] and evaluating the possibility of future collapses and related impacts on the environment and human activities [89]. Coastal and marine mass wasting can also be related to past and present sea-level rise due to climate change. This is the case of the Vasto landslide (Abruzzo, Italy), a rotational slide continuing under the sea level, whose geomorphological evolution, and past and historical reactivations have been reconstructed by Della Seta et al. [96]. Another example comes from the Maltese coastal block slides that developed during the last glaciation and were then influenced by the successive sea-level oscillations [180].
Among hazards induced by the ongoing climate change, coastal inundation (especially along stretches of coast affected by subsidence; [108]) and coastal erosion triggered by extreme meteorological events, and sea-level rise are the most reported in the literature. In particular, several papers concern (i) the quantification of sea-level rise (e.g., [2,93,107,111]; Figure 9), (ii) the general impacts of sea-level rise on coastal environments [97,102] and (iii) the assessment of coastal exposure and coastal erosion (e.g., [98,114,117,118]). Terrestrial and marine datasets are differently analyzed and integrated, although bathymetry is often taken into consideration, especially in those works presenting predictive models on sea-level rise and coastal inundation.

Marine and Landscape Ecology
In the last few years, attention has been focused on the so-called "mapping from ridge to reef" approach [12] in order to investigate the connectivity between upland watersheds, intertidal zones and shallow coastal areas, including the influence that coastal (or riverine) and submarine morphological features can have on habitats' distribution.
As a matter of fact, habitats can be defined as physically distinct areas characterized by specific physical, chemical and biological properties (and oceanographic properties as well for benthic habitats) and hosting distinct species or communities of species. Among the physical components, seabed morphology and its geomorphological significance can have a remarkable control on ecological processed and associated biota [191,192]. Indeed, it is known that a wide variety of terrain attributes (e.g., orientation, slope angle, roughness), substrate type and chemical and oceanographic variables deeply affect species distribution and in turn biodiversity, providing surrogates used to identify places that deserve protection [68]. In submarine environments, different landforms are usually associated with specific benthic habitats as discussed in Harris and Baker [129] and some species are defined as ecosystem engineers providing themselves typical submarine landforms or geomorphic proxies for habitat detection, even in the deep submarine environments (cf. [193][194][195][196][197][198]). In tropical coastal environments detailed and accurate representations of topography and bathymetry are essential for habitat mapping [199,200], since they are required for modelling nearshore hydrodynamics, sediment transport and reef evolutionary processes [21] and references therein, [201]. Finally, we must consider that terrestrial and marine ecosystems are linked by freshwater inputs (e.g., rivers discharge) that supply sediments, nutrient exchange and larval transport, and pollutants [121,124].
In this framework, the acquisition of reliable base maps, in terms of elevation data, for both onland and marine environments constitutes the basis for any studies aiming at analyzing seabed landscapes. Benthic habitat mapping means "plotting the distribution and extent of habitats to create a map with complete coverage of the seabed showing distinct boundaries separating adjacent habitats" as stated within the MESH project [202]. Benthic habitat and, more generally, habitat mapping practices constitute a basic tool for habitats conservation as part of an ecosystem approach [168] (cf. Section 2.7).
The growing interest in seafloor mapping, habitat mapping and development of an integrated management of coastal and marine environment fostered large use of abiotic surrogates to represent biodiversity [126], as witnessed also by the international GeoHab (Marine Geological and Biological Habitat Mapping) community [127,203] and at the European level by the MAREANO program in Norway [204][205][206], and the MAREMAP in the UK [207,208].
Important contributions including habitat maps generated by coupling terrestrial and underwater geospatial datasets are from Hogrefe et al. [121], McKean et al. [122,123], Vierling et al., [125], Leon et al. [21], Marchese et al. [127] and Prampolini et al. [6,128]. The latter exploited latest technological advancement, among which the LiDAR-derived elevation data of the Maltese coast and of its seafloor down to a depth of ca. 50 m for mapping both geomorphological features and habitats.

Marine and Landscape Ecology
In the last few years, attention has been focused on the so-called "mapping from ridge to reef" approach [12] in order to investigate the connectivity between upland watersheds, intertidal zones and shallow coastal areas, including the influence that coastal (or riverine) and submarine morphological features can have on habitats' distribution.
As a matter of fact, habitats can be defined as physically distinct areas characterized by specific physical, chemical and biological properties (and oceanographic properties as well for benthic habitats) and hosting distinct species or communities of species. Among the physical components, seabed morphology and its geomorphological significance can have a remarkable control on ecological processed and associated biota [191,192]. Indeed, it is known that a wide variety of terrain attributes (e.g., orientation, slope angle, roughness), substrate type and chemical and oceanographic variables deeply affect species distribution and in turn biodiversity, providing surrogates used to identify places that deserve protection [68]. In submarine environments, different landforms are usually associated with specific benthic habitats as discussed in Harris and Baker [129] and some species are defined as ecosystem engineers providing themselves typical submarine landforms or geomorphic proxies for habitat detection, even in the deep submarine environments (cf. [193][194][195][196][197][198]). In tropical coastal environments detailed and accurate representations of topography and bathymetry are essential for habitat mapping [199,200], since they are required for modelling nearshore hydrodynamics, sediment transport and reef evolutionary processes [21] and references therein, [201]. Finally, we must consider that terrestrial and marine ecosystems are linked by freshwater inputs (e.g., rivers discharge) that supply sediments, nutrient exchange and larval transport, and pollutants [121,124].
In this framework, the acquisition of reliable base maps, in terms of elevation data, for both on-land and marine environments constitutes the basis for any studies aiming at analyzing seabed landscapes. Benthic habitat mapping means "plotting the distribution and extent of habitats to create a map with complete coverage of the seabed showing distinct boundaries separating adjacent habitats" as stated within the MESH project [202]. Benthic habitat and, more generally, habitat mapping practices constitute a basic tool for habitats conservation as part of an ecosystem approach [168] (cf. Section 2.7).
The growing interest in seafloor mapping, habitat mapping and development of an integrated management of coastal and marine environment fostered large use of abiotic surrogates to represent biodiversity [126], as witnessed also by the international GeoHab (Marine Geological and Biological Habitat Mapping) community [127,203] and at the European level by the MAREANO program in Norway [204][205][206], and the MAREMAP in the UK [207,208].
Important contributions including habitat maps generated by coupling terrestrial and underwater geospatial datasets are from Hogrefe et al. [121], McKean et al. [122,123], Vierling et al., [125], Leon et al. [21], Marchese et al. [127] and Prampolini et al. [6,128]. The latter exploited latest technological advancement, among which the LiDAR-derived elevation data of the Maltese coast and of its seafloor down to a depth of ca. 50 m for mapping both geomorphological features and habitats. Leon et al. [21] produced a seamless and high-resolution DEM of the fringing reef system of Lizard Island in northern Great Barrier Reef (Australia), merging multisource 3D models (topographic and bathymetric LiDAR data, passive optical remote sensing data, nautical charts, and single-beam and multi-beam echo-sounder data reaching 30 m b.s.l.). McKean et al. [122,123] tested the high-resolution Experimental Advanced Airborne Research LiDAR (EAARL), a new technology for cross-environment surveys of channels and floodplains, to acquire elevation-depth data of a channel in the Bear Valley Creek (Idaho, USA), and map its landforms and habitats.
Combining terrestrial and bathymetric LiDAR allows reconstructing a 3D view of terrestrial habitats (e.g., St. Joe Woodlands and sagebrush-steppe ecosystem in Idaho (USA) as showed by Vierling et al. [125]). Hogrefe et al. [121] combined depths derived from IKONOS satellite imagery and sonar data to produce a seamless DEM of Tutuila Island (American Samoa) that can be used for evaluating the assessment of human population and land use practices on coral reefs. Finally, worthy of note are the recent studies reporting first applications of photogrammetric technique to UAV imagery to map coral habitats in tropical coastal environments [127,209].

Coastal Planning and Management
Land and sea interaction is increasingly perceived as relevant in the context of planning and management of terrestrial and sea areas, since most of the activities occurring in the marine environment are also connected with the terrestrial vicinities. The interdependence of land and offshore systems drives the need for integration between terrestrial and marine planning systems, considering driver issues that cross the land/sea boundary [138]. Among others, changes in both landscapes and seascapes due to urbanization and anthropogenic activities represent a key element to consider within any planning processes.
The land-sea interactions and related processes constitute a key element of the Mid-Term Strategy 2016-2021 of UN Environment /MAP adopted with Decision IG. 22/1 [147], and correspond to the first objective of both the Mediterranean Strategy for Sustainable Development (MSSD) 2016-2025, adopted with Decision IG 22/2 [148], and the Sustainable Development Goals 14 and 15 [149,150]. Indeed, the goals of "Life below water" (SDG 14) and "Life on land" (SDG 15) are strictly interconnected [145].
In this context, few diverse approaches facilitate land-sea planning system integration. Among them there are the Integrated Coastal Zone Management (ICZM) and the Maritime Spatial Planning (MSP) tools, coupled with an Ecosystem-Based Approach (EBA) [137].
The ICZM initiatives provide a support to integrated and holistic planning and management of the coastal areas, including both the land (inland limit decided by the countries) and marine (territorial seas) components (Art. 1 of the Protocol on ICZM in the Mediterranean [151]). The importance of considering land and sea space as a whole within the ICZM process is re-affirmed by some of the Protocol's objectives and principles as for example the following: "Ensure preservation of the integrity of the coastal ecosystems, landscape and geomorphology" (Art. 5; objective d). Given the definition of the coastal zone provided by the Protocol, this integrity can be preserved only if the land and marine parts of the landscape are considered together.
The "Conceptual Framework for MSP in the Mediterranean" (Barcelona Convention, December 2017 [146]) foster this integration also facilitating the introduction of MSP into ICZM in the framework of the Barcelona Convention Protocols. A step by step methodology for the implementation of the MSP following common principles in the Mediterranean has been designed thanks to the existing guiding documents ( Figure 10) [132,134,[141][142][143].
Recent examples of integrated approaches to terrestrial and marine spatial planning occur in the design of a network of MPAs using models such as Marxan ("marine reserve design using spatially explicit annealing"), the most widely used software at global level for conservation planning and designed for solving conservation planning issues in landscapes and seascapes [135].
Following an ecosystem-based approach within the MSP and ICZM leads to an evolution within the different planning and management actions, taking into account the need to embrace multidisciplinary approaches and to advocate cross-realm connectivity [144]. Recent examples of integrated approaches to terrestrial and marine spatial planning occur in the design of a network of MPAs using models such as Marxan ("marine reserve design using spatially explicit annealing"), the most widely used software at global level for conservation planning and designed for solving conservation planning issues in landscapes and seascapes [135].
Following an ecosystem-based approach within the MSP and ICZM leads to an evolution within the different planning and management actions, taking into account the need to embrace multidisciplinary approaches and to advocate cross-realm connectivity [144].
These imply the integration of spatial data stored in a GIS, including relevant Earth Observation services and the characterization of marine habitats and seabed landscapes, especially as premise for representing coastal and maritime space [136].
MSP practice highlights the necessity of having a strong data and knowledge base among its principles [210]. Hence, the integration of terrestrial and marine spatial dataset constitutes a fundamental element to be able to undertake a new interdisciplinary approach for an integrated analysis of marine ecosystems and common maritime space [137].

Advances in Data Collection Technology and Data Processing Methodology
The accomplishment of the outcomes reported in the literature and examples listed in previous chapters is mostly due to the widespread recent availability of new technologies and software that enable scientists to acquire geo-environmental spatial data that were unrecoverable before the 1970s in the submarine environments. Their integration with terrestrial data has become feasible especially with the generation of data format and products suitable for implementation into GIS platforms that in turn made possible to handle and analyze complex and heterogeneous datasets from the onshore to the offshore zone, as shown in most of the previously cited works. These imply the integration of spatial data stored in a GIS, including relevant Earth Observation services and the characterization of marine habitats and seabed landscapes, especially as premise for representing coastal and maritime space [136].
MSP practice highlights the necessity of having a strong data and knowledge base among its principles [210]. Hence, the integration of terrestrial and marine spatial dataset constitutes a fundamental element to be able to undertake a new interdisciplinary approach for an integrated analysis of marine ecosystems and common maritime space [137].

Advances in Data Collection Technology and Data Processing Methodology
The accomplishment of the outcomes reported in the literature and examples listed in previous chapters is mostly due to the widespread recent availability of new technologies and software that enable scientists to acquire geo-environmental spatial data that were unrecoverable before the 1970s in the submarine environments. Their integration with terrestrial data has become feasible especially with the generation of data format and products suitable for implementation into GIS platforms that in turn made possible to handle and analyze complex and heterogeneous datasets from the onshore to the offshore zone, as shown in most of the previously cited works.
The expansion of GIScience can be dated back to the 1990s. As soon as marine datasets became accurately "geo-referenced", thanks to our ability in obtaining geographical positions at sea through the development of GPS, their structure, format and way of representation moved immediately toward the form of geospatially enabling the data to create maps and 3D scenes of the marine environment. Since "marine" GIS has evolved adapting a technology originally designed for land-based applications, the integration of marine and terrestrial datasets has been quite immediate. This especially happened for those marine studies dealing with coastal environments or applications that benefit from the investigation of spatial relationships within and between marine and terrestrial dataset-such as measuring and monitoring the seascape or modelling/predict future scenarios [211]. An important trigger in applying this new approach for coastal and marine data visualization and analysis came undoubtedly from industry, especially in the fields of hydrographic surveying and production of nautical charts and publications.
As soon as marine and terrestrial elevation datasets started to be implemented and harmonized in a continuum dataset, the British Geological Survey introduced the new term "white ribbon" in the hydrographic sector, to designate the information gap of elevation data in the shallow area formed by the intertidal and nearshore zones, meaning the interface between land and sea. Covering the white ribbon with high-resolution bathymetric data became soon a challenge in all first attempts devoted to integrate marine and terrestrial spatial datasets (e.g., [121,124,[212][213][214][215]), and the scientific community soon realized both the relevance and the issues to be addressed in carrying out topo-bathymetric surveys [15,28,190,[216][217][218][219][220][221]. Most of the difficulties in getting elevation data in the white ribbon are caused by the water depth: it is generally too shallow for traditional bathymetric surveys (because of the draft and the need to submerge the echo-sounders keeping them at a certain distance from the seafloor to obtain an efficient coverage, and because of the unsafe conditions generated by the common occurrence of rocky outcrops and/or waves) and too deep for traditional optical land-based survey methods. Shallow water is in addition more expensive to be surveyed than deeper ones since MBES seafloor coverage is narrowed as water becomes shallower, requiring the vessel to spend excessive time in shallower areas due to the need to run very close sur vey lines to achieve adequate coverage ( Figure 1). Nonetheless, the intertidal and nearshore zones are of extreme importance to investigate for all the reasons listed in the previous sections. This challenge favored the development of new advanced methods and techniques to improve the capability of obtaining shallow-water bathymetric data, and especially: • The production of new advanced acoustic systems designed for obtaining depth measurements in shallow water; • The application of cutting-edge visualization technology to images and data collected with optical sensors to obtain elevation data from shallow areas (i.e., underwater photogrammetry, image derived bathymetry, LiDAR, laser scanning).
Finally, given all the technological and economic difficulties mentioned above in mapping the seafloor, both in deep and shallow water, sharing data is increasingly appreciated and encouraged by several research funding programs. The goal "map once, use many times" supports the creation of national (underpinned by governments), regional and international repositories of bathymetric data. Examples of regional repositories are the European EMODnet (www.emodnet-bathymetry.eu) and Baltic Sea Bathymetry Database (http://data.bshc.pro), while at international/global level, GEBCO (https://www.gebco.net/) is a repository of world bathymetry, which is also updated thanks to local portals (e.g., from EMODnet), and currently under update thanks to GEBCO Seabed 2030 Project, aiming at completing mapping of the world ocean by 2030.

Shallow-Water Acoustic Systems
Most of the interest in surveying shallow water comes from industry (oil and gas companies, port and harbor authorities and maritime engineering among others) and academics and has grown rapidly in recent years. This has pushed manufacturers to both produce MBES systems (i.e., "beamforming system"- [222]) adapted for fast mobilization on smaller vessels (easy-to-use and quick-to-deploy) and explore new innovations in swath bathymetry systems, developing novel swath sonar technology to reach greater seafloor coverage (up to 15 times the depth), such as (i) the interferometric echo-sounders, also known as Phase Differencing Bathymetric Systems (PDBS) [222,223] or (ii) the multi-phase echo-sounders (MPES- [224]).
Technological developments in beamforming system have especially affected the geometry and the performance of the transmit array and sounding frequencies, refining the capability of the systems in offering a wider coverage (up to 7 times the water depth) and narrower acoustic beam (reaching accuracies that have been shown to exceed the IHO-International Hydrographic Organization-Standards for Hydrographic Surveys). The broader coverage (i.e., swath) is obtained using multi-transducer multi-beam products. Two sonar heads (i.e., transducers), for instance, can easily achieve double the coverage, by simply adjusting the angle of the heads.
Developments in shallow-water swath bathymetry systems involved all aspects of the "seafloor mapping system", including all ancillary sensors and software involved in the survey to provide the so-called "integrated survey system", namely the GPS/GNSS positioning systems, the motion sensors, and sound velocity recording sensors. The goal is to simplify installation and calibration procedure and make the shallow-water MBES systems perfect for use on vessels of opportunity, small survey launches, and even Autonomous Surface Vessels (ASVs) or USVs. ASV/USV are vehicles that can navigate and collect data from the surface of the water without a crew. ASV/USV are currently produced to remotely control data acquisition especially in shallow marine water, rivers and channels.

Optical Sensor for Underwater Imaging and Mapping
Latest technological developments in underwater 3D reconstruction, based on airborne active optical sensor, has given rise to a wide range of new systems and techniques such as the LiDAR systems, Structured Light (SL), Laser Stripe (LS), Laser Line Scanning (LLS), Stereo Vision (SV) and SfM [225][226][227][228].
LiDAR systems have substantially improved shallow seafloor mapping in coastal environments. Using infrared laser pulses, topographic LiDAR systems achieve a very high-resolution mapping performance (i.e., meter to sub-meter point spacing) with sub-meter vertical accuracy. Bathymetric LiDAR systems, can even use both infrared and blue-green laser pulses, to simultaneously acquire depth measurements down to~70 m below Mean Sea Level (MSL), according to water turbidity, typically <40 m is achieved in most applications. Bathymetric LiDAR systems have been the first active optical sensors that provided elevation data from the nearshore areas, allowing surveying shallow seafloor with much more efficiency in terms of coverage and required time (Figure 1). The LiDAR technology can presently be integrated also to terrestrial or surface or underwater platforms, carrying out ROV-based LiDAR inspection surveys, benefiting from increased spatial and temporal resolution, and greater accuracy [229]. LLS can now be used just like a multi-beam although the technology is slightly different. Both subsea LiDAR and Laser scanners generate a relative point cloud (referenced by flow and bearing measurements) with a resolution of even millimeters, i.e., much higher than any acoustics-based system. The high resolution offered by subsea LiDAR or laser and the need to operate under lighting conditions determine, however, a limited range, strongly regulated by the environmental conditions. The resolution and accuracy typical of LiDAR/submarine laser systems require indeed clear water with good visibility [230]. Thus, they cannot be employed to scan those nearshore areas characterized by high turbidity such as the ones close to river discharge or with sediment/pollutant moved by water movements [231,232].
Hence, starting from the 1990s, active sensors based on underwater acoustic (e.g., multi-beam or interferometric echo-sounders) and light signals LiDAR/LLS systems improved substantially the capacity of obtaining elevation data in underwater environments, despite expensive techniques, especially for small scale surveys. These instruments directly provide a point cloud of bathymetric measurements that can generate a DTM with sub-meter resolution. LiDAR systems can even combine onshore topographic and nearshore bathymetric mapping obtaining detailed emerged and submerged surfaces in a single acquisition ( Figure 11) [233].
With the availability of high-resolution images collected by satellite remote sensors (or even by drones), 3D underwater models can now be generated using also passive optical systems, at least for those areas in which visible light can penetrate down to the seafloor. Satellite-Derived Bathymetry (SDB) is indeed the most recently developed method of surveying shallow waters. Different companies developed ad-hoc algorithms since the 1990s to convert the information collected by satellite sensor into bathymetric data. SDB is based on the connection between the seabed reflected energy and the depth of the water [234]. The method, using dedicated computational algorithms, basically exploits, for each pixel of the satellite image where the seafloor is visible, the statistical relationship between the depth of water and the type and intensity of energy detected by the sensor. Since SDB can estimates the water depth of the seafloor up to the extent of light penetration into the water medium (i.e., around 20-30 m under optimal conditions), water transparency is the main limiting factor. Atmospheric absorption, sun glint, high substrate heterogeneity, algal blooms, suspended sediment, or waves, can also all limit SDB performance. measurements that can generate a DTM with sub-meter resolution. LiDAR systems can even combine onshore topographic and nearshore bathymetric mapping obtaining detailed emerged and submerged surfaces in a single acquisition ( Figure 11) [233]. Figure 11. Integration of terrestrial and marine elevation data derived from bathymetric LiDAR: example from the NW coast of Malta, Central Mediterranean Sea (modified after Prampolini et al. [6]).
With the availability of high-resolution images collected by satellite remote sensors (or even by drones), 3D underwater models can now be generated using also passive optical systems, at least for those areas in which visible light can penetrate down to the seafloor. Satellite-Derived Bathymetry (SDB) is indeed the most recently developed method of surveying shallow waters. Different companies developed ad-hoc algorithms since the 1990s to convert the information collected by satellite sensor into bathymetric data. SDB is based on the connection between the seabed reflected energy and the depth of the water [234]. The method, using dedicated computational algorithms, basically exploits, for each pixel of the satellite image where the seafloor is visible, the statistical relationship between the depth of water and the type and intensity of energy detected by the sensor. Since SDB can estimates the water depth of the seafloor up to the extent of light penetration into the water medium (i.e., around 20-30 m under optimal conditions), water transparency is the main limiting factor. Atmospheric absorption, sun glint, high substrate heterogeneity, algal blooms, suspended sediment, or waves, can also all limit SDB performance.
Finally, with the advent of underwater camera systems, progress in deep-sea robotics, and the increased number of videos and images being captured underwater, researchers began to obtain optical 3D reconstruction of recorded scenes with (sub)centimetric resolution, employing numerous techniques, among which the use of stereo cameras and the principle of photogrammetry even underwater (i.e., SfM technique). Optical underwater imaging is emerging as a key technology for a variety of oceanography applications [227] and [235] among others.
However, very few techniques employing photorealistic seafloor imagery take critically into account the extent to which scattering affects the scenes captured under daylight in shallow water or using active illumination in deep water. In most cases reported in the literature, it is implicitly assumed that light is neither absorbed nor scattered by the medium in which the source and scene are immerged (as it happens in pure air [236]). However, the major challenge facing optical imaging in these applications is the severe degradation of image quality caused by scattering generated by impurities and organisms. SfM has been also applied to UAV-based RGB imagery, on coastal waters [237,238] among others. UAV imagery processed with SfM techniques offers a low-cost alternative to established shallow seafloor mapping techniques providing also important visual information with Figure 11. Integration of terrestrial and marine elevation data derived from bathymetric LiDAR: example from the NW coast of Malta, Central Mediterranean Sea (modified after Prampolini et al. [6]).
Finally, with the advent of underwater camera systems, progress in deep-sea robotics, and the increased number of videos and images being captured underwater, researchers began to obtain optical 3D reconstruction of recorded scenes with (sub)centimetric resolution, employing numerous techniques, among which the use of stereo cameras and the principle of photogrammetry even underwater (i.e., SfM technique). Optical underwater imaging is emerging as a key technology for a variety of oceanography applications [227] and [235] among others.
However, very few techniques employing photorealistic seafloor imagery take critically into account the extent to which scattering affects the scenes captured under daylight in shallow water or using active illumination in deep water. In most cases reported in the literature, it is implicitly assumed that light is neither absorbed nor scattered by the medium in which the source and scene are immerged (as it happens in pure air [236]). However, the major challenge facing optical imaging in these applications is the severe degradation of image quality caused by scattering generated by impurities and organisms. SfM has been also applied to UAV-based RGB imagery, on coastal waters [237,238] among others. UAV imagery processed with SfM techniques offers a low-cost alternative to established shallow seafloor mapping techniques providing also important visual information with the generation of an orthomosaic for the surveyed scene [17]. Nevertheless, water refraction introduces severe errors when UAVs imagery is used for bathymetric applications. Although the application of photogrammetric procedures on images captured directly in the water medium (in-water) needs only a thorough calibration to correct the effects of refraction, in instance where the image acquisition occurs through-water (two-media), the sea surface undulations caused by waves [239,240] and the magnitude of refraction that can change at each point of every image, lead to uncertain results [241,242].
Overall, it is clear that no single optical imaging system can meet all the needs of underwater 3D reconstruction. The different systems cover very different spatial scales, resolutions and accuracy, being suitable for different applications. Furthermore, it is important to underline the lack of systematic studies to precisely compare the performance of different sensors in relation to the same scenario and under identical conditions. However, technology and computer vision are definitely on the way of addressing all pitfalls of the mentioned applications, to obtain 3D optical reconstructions more reliable over multiple spatial scales, through innovative sensors and data processing. Finally, it should be emphasized that the possibility of obtaining accurate photorealistic 3D reconstructions, also allows the use of interactive tools for visualization and exploration in 3 dimensions, designed to support the interpretation and analysis of the obtained spatial data. These large, complex and multicomponent spatial datasets can indeed be used to develop innovative learning tools for environmental sciences, presenting new worlds of interactive exploration to a multitude of users [243].

Conclusions and Future Perspectives
Most terrestrial landscapes and landforms have always been investigated for several purposes and benefited from high accessibility to surveyed areas and high-resolution data. On the contrary, the submarine environment has struggled in being represented with the same resolution and coverage as its terrestrial counterpart because of its remoteness and/or its limited accessibility and the high costs imposed by the need to use expensive infrastructures and sophisticated technologies (especially research vessels, underwater robotics and a multitude of sensors). Nevertheless, remote sensing tools operating from satellite, aerial platforms and vessels or autonomous vessels and drones have contributed to obtain elevation data for both land and sea areas with a comparable resolution, establishing submarine geomorphology as a field of research that is also remarkably contributing to marine environmental management, with an increase in many associated applicable research.
This paper has pointed out seven good reasons to pursuit such a comprehensive and homogeneous integration of terrestrial and submarine datasets, showing the outputs of relevant research in this field. The interest in producing integrated land-sea geomorphological maps is now at its beginning, even though it would be the basis for further applied research. The integrated assessment of geoheritage and geodiversity in coastal and marine environments has been the subject of a very limited number of papers so far. A much larger number of papers refers instead on the coupling of terrestrial and marine spatial datasets aiming at reconstructing paleo-landscapes in coastal and marine areas and outlining their geomorphological evolution, supporting also the identification and study of archaeological sites. The field of application that has mostly benefited so far from the integration of terrestrial and marine datasets is the integrated assessment of geohazards in coastal and marine areas, with special reference to tsunami and storms, coastal and marine landslides and sea-level-rise-related hazards. Attention to the interrelations between land and sea and their effects on marine habitats is being paid, but still only few works show combined investigation in terrestrial and nearshore environments. Finally, in the last ten years, stakeholders have pointed out the need for an integrated planning and management of marine and terrestrial areas. During this span of time, discussion has extended from MSP alone to ICZM of coastal areas thanks to several directives and plans developed and applied at international level (i.e., European Union, United Nations).
In the near future, it is likely that technological improvements will allow an increasingly easier accessibility to the "white ribbon" and a better integration of terrestrial and marine spatial datasets that strongly relies on the realization of seamless DTMs of land and sea areas. This will enhance further and much wider research in transition environments based on interdisciplinary approaches. It is also desirable that a standardization and/or harmonization of data will be soon achieved, to adopt common terminology and classification schemes for both terrestrial and submarine geomorphological features.
The resulting increased awareness of the interconnection between landscapes of terrestrial and marine areas calls for a holistic approach to better understand environmental changes taking place on Earth and to consequently design appropriate management measures. Scientists and stakeholders typically working on terrestrial and marine areas separately will hopefully understand the benefits of coupling terrestrial and marine investigation and activities, which is of paramount importance for environmental protection and enhancement.