A High-Resolution Digital Bathymetric Elevation Model Derived from ICESat-2 for Adam’s Bridge

This data descriptor elaborates the details of a high-resolution digital bathymetric elevation model generated for the region, namely, Adam’s Bridge, which encompasses a chain of shoals between Rameswaram Island, off the southeastern coast of Tamil Nadu, India, and Mannar Island, off the northwestern coast of Sri Lanka. The proposed dataset has taken advantage of the photon penetrability in the shallow waters by the green laser of ICESat-2 LiDAR to derive the seabed topography. Seafloor depths from ~0.2 million geolocated photons of ICESat-2 for the study area were accrued and interpolated to generate a 10 m digital bathymetric elevation model. Adam’s Bridge, an isthmus and submerged reefal assemblage in shallow and super-shallow waters, is a feature of scientific curiosity. Our dataset has the potential to enhance the understanding of Adam’s Bridge structure by providing substantial information to reconstruct its evolution.


Background & Summary
Adam's Bridge is an isthmus of length ~30 km connecting Dhanushkodi, the southeastern point of Rameswaram Island, off the southeastern coast of Tamil Nadu, India, and Talai Mannar, the western end of Mannar Island, off the northwestern coast of Sri Lanka [1][2][3][4][5][6][7][8][9][10][11] ; the extent along with its surroundings is shown in Fig. 1 comprising high-resolution satellite imagery.The southern part of the Adam's Bridge has the Gulf of Mannar, an arm of the Indian Ocean, and Palk Strait, an inlet of the Bay of Bengal, in its northern direction.
On this heavily submerged isthmus, loose and fine sandbanks are exposed at irregular intervals above the water level.Generally, these sandbanks exhibit seasonal shifts; moreover, these sandbanks have no presence of rocks 1,2,10 .Most of this isthmus crest (the central position of the submerged ridge) is within the depth of 1 m of the water currents with great rapidity, but in some parts, sudden channels of depth 2-3 m exist; however, this depth is more when the channels are narrower 10 .Studies concerning the material strata through borings done by various surveyors yielded sand up to a depth of 7-9 m, and further depths yielded Holocene conglomerate, beds of limestone, calcareous sandstone, and occasional corals 1,7,8,10,12,13 .Next to this second strata, a layer of sand followed by hard rocks or continuous compact formation was detected 12,13 .Such a type of reefal assemblage is a matter of scientific curiosity, and further understanding of its morphological structure can provide substantial information to reconstruct its evolution.
Besides the interest in the unique geological structure of Adam's Bridge, right from the 19 th century, it has remained at the center of controversies and debates surrounding the reasons for attributing its name and origin (man-made or geological formation)  . From he ages, with reference to the Hindu mythology called Ramayana, this geographic extent, in most South Asian countries, is referred to as Ram Setu [23][24][25][26][27]46,47 .Major James

Methods
Study area.The extent of the study area for which the proposed dataset was generated is shown in Fig. 1b and marked with a white box.The presence of the Adam's Bridge, as seen in optical satellite data, is shown in Fig. 1c,where mostly, in its entire form, is a submerged ridge with occasional exposure to sand cays at irregular intervals.From the tip of Dhanushkodi, Adam's Bridge's general trending long axis is oriented from west-northwest to east-southeast in secondary intercardinal directions, entirely perpendicular to the predominant wave approach directions of both the Gulf of Mannar and the Palk Strait.The general direction of Adam's Bridge for the first 17 km of its length is about east-south-east; later, the direction is gradually curved towards the north of the east and finally touches Talaimannar Island at its east.The study area includes neighbouring stretches of Adam's Bridge like Rameswaram Island (on the Indian side), part of Mannar Island (on the Sri Lankan side), and the surrounding waters of the Gulf of Mannar and Palk Strait; the rationale for extending the study area beyond the core area of interest, i.e., Adam's Bridge, is to accommodate the relationship between the location of interest and the distance to surrounding sample points during the digital seafloor surface generation.
Data resources and methodology.The primary data resource used to generate DBEM for Adam's Bridge is the depth information of the seafloor retrieved using water-penetrated ICESat-2 photons.The ICESat-2 mission operates in a non-sun-synchronous orbit (thus, both the day-time and night-time acquisitions are possible) from an average altitude of 496 km with a temporal resolution of ~91 days, during which 1387 unique ground www.nature.com/scientificdata/track acquisitions will be possible for the global extent 73 .The along-track sampling interval of ICESat-2 photons is 0.7 m (spatial resolution).For the extent of Adam's Bridge, there exist seven ICESat-2 ground tracks (shown in Fig. 2), and thus, data acquisitions from these available ground tracks were considered in this research.From these seven reference ground tracks, one hundred thirty-three tracks of ICESat-2 data acquisitions were available between October 2018 and October 2023 for the study area considered in this research.
ICESat-2, equipped with a solo sensor called Advanced Topographic Laser Altimeter System (ATLAS), uses a 532 nm wavelength (green) laser operating at a pulse repetition frequency of 10 kHz.A single incident of laser from ATLAS encounters a diffractive optical element before splitting into six beams (organized as three pairsleft, near-nadir, and right); of which, within each pair, one beam (termed as a strong beam) has four times the energy of the other (termed as a weak beam) with a separation of 90 m between them 73 .Figure 3a illustrates the concept of ICESat-2's multi-beam data acquisition over a part of Adam's Bridge.The left and right paired beams are 3.3 km apart from the nadir most paired beam.
Once the laser beams from ICESat-2 hit the earth's surface, reflected photons recorded by the photon-counting telescope mounted on the ATLAS sensor yield the range measurements.Data related to the ICESat-2 platform (like position, orientation, attitude, and orbital velocity), laser pointing vectors, pulse emission timing, and the range measurements will be assimilated to generate the geodetic position (latitude and longitude) and ellipsoidal heights for each geolocated photon 73 .At ground stations, science teams of ICESat-2 will process and distribute various levels of data products through a web portal at https://nsidc.org/data/icesat-2 maintained by the National Snow and Ice Data Centre (NSIDC).Global geolocated photon data available as a Level-2A product under the nomenclature ATL03 contain heights above the WGS84 ellipsoid (ITRF2014 reference frame), geodetic latitude and longitude, and other relevant attributes for all the photons downlinked by the ATLAS instrument onboard the ICESat-2 81 .
For our research, ATL03 data was downloaded using the OpenAltimetry application 82 , which is available at https://openaltimetry.earthdatacloud.nasa.gov/data/icesat2/.The OpenAltimetry application, a NASA-funded collaboration between NSIDC, Scripps Institution of Oceanography, the EarthScope Consortium, and the University of California San Diego, is a web-based cyberinfrastructure platform that allows users to locate, visualize, and download ICESat-2 surface elevation data and photon clouds for any location on Earth 82 .Photon data acquisitions from the ICESat-2 are available in Comma Separated Values (CSV) format; thus, the data can be analyzed in any electronic Spreadsheet application like Microsoft Excel or Google Sheets.Even though ATL03 is a non-imaging product (tabular data), the presence of attributes like geodetic latitude and longitude in the data will enable us to infer and analyze in the spatial domain using Geographic Information System (GIS) software like ESRI ArcGIS (https://www.esri.com/)or QGIS (https://www.qgis.org/).
Along-track 2D profiles generated from a sequence of ICESat-2 photon data enable visualization of Earth's surface variations.Figure 4a shows a subset of ICESat-2 geolocated photons overlaid on high-resolution satellite imagery for the extent of the study area.Figure 4b shows the accumulated geolocated photons in the 2D profile chart, with the X-axis having the along-track latitude acquired by the ICESat-2 ground track and the Y-axis having the elevation (ellipsoidal height) in meters.www.nature.com/scientificdata/During ICESat-2's data acquisition over the water bodies, most photons will get reflected from the water surface.However, depending on the optical properties of the water, some of the photons will return from the water column and seafloor [76][77][78][83][84][85] (shown in Fig. 4b). The difuse attenuation coefficient for downwelling irradiance, K d (λ) (in m −1 ), where λ is the wavelength of light, is a measure of how light dissipates with depth in water; it indicates how strongly light intensity at a specified wavelength is attenuated within the water column and generally computed from remote sensing techniques 86,87 .Open access and daily temporal resolution-based data like K d ( 490) is available from the Sentinel-3 A/B OLCI level-2 series of data product services and helps characterize the transparency of water [88][89][90] .For successful seafloor detection up to a depth of ~40 m from the ICESat-2 photons, the recommended value 76,91 of K d (490) should be less than 0.12 m −1 .In our research, acquisitions during high turbid load in the water were avoided by referring to the K d (490) of Sentinel-3 A/B; during this crosschecking procedure, the overlap period between ICESat-2 and Sentinel-3 A/B acquisitions is within +/− 24 hours.Similarly, for water depth-related studies using ICESat-2, it is recommended to prefer night-time acquisitions because the background noise caused by solar spectral radiation significantly impacts the seafloor detection performance of LiDAR 92,93 .Thus, only those acquisitions of ICESat-2 obtained during night-time were considered to generate DBEM for the study area.
Abiding by prerequisite conditions like preferring night-time acquisitions and omitting the data acquired during turbid load periods, 66 tracks are found qualified from the available 133 tracks, and these tracks comprise 396 strong and weak beams of along-track data.All these 396 beams (shown in Fig. 3b) of data were processed to classify the returned photons from the water surface, water column, land, and seafloor using Density-Based Spatial Clustering of Applications with Noise (DBSCAN) [94][95][96][97] followed by manual correction using localized statistical algorithms 98,99 to eliminate the outliers.Typically, the DBSCAN algorithm, by taking parameters like radius and threshold of minimum points, will make clusters of points as individual units/classes within the said radius when the density of points exceeds the pre-set threshold of minimum points.Also, the DBSCAN algorithm identifies the points as outliers in low-density regions (refer to Fig. 4b,where outliers exist above the water surface and below the seabed).Figure 4c shows the classified photons based on their interaction with the surface feature for a subset of the ICESat-2 beam.and need applying refraction correction to retrieve the actual depth 76 ; this is because there will be a change in the speed of light that occurs at the air-water interface due to the difference of refractive indices of air and sea-water.Avoiding the refraction correction for the depths from the photons returned from the seafloor includes vertical www.nature.com/scientificdata/errors in the bathymetric estimation due to apparent depth values.Figure 4c shows the difference in depths due to the apparent seafloor and the refraction-corrected seafloor.For the DBEM generated for the extent of Adam's Bridge, we implemented the refraction correction method 76,85 based on Eq. 1. www.nature.com/scientificdata/ In Eq. 1, Depth C is the refraction-corrected depth of the seafloor, and Depth apparent is the apparent depth (default depth) of the photons returned from the seafloor.n 1 ≈ 1.00029 and n2 ≈ 1.34116 are the refractive indices of air and water bodies, respectively.
Only those photons that returned from the seafloor and land were considered for further processing.Photons returned from the water surface and water column were discarded during the computation of bathymetry as they do not have any role in estimating the seafloor depths.Figure 5 shows profiles generated for two different acquisitions of ICESat-2 beams at the extent of Adam's Bridge.It is evident that from Fig. 5c,d, both the exposed sandbanks (above the mean sea level) and the seafloor depths are seen in the profiles; thus, to understand the intricate details of Adam's Bridge, elevation values (of those features that are above the mean sea level) and the sea floor depth prove to be vital; else, a simple digital bathymetric model that holds only depth values (below the mean sea level) may not give information about the exposed features of the study area.
Figure 6 shows the profiles for two acquisitions of ICESat-2 beams over Pamban Island and Dhanushkodi, which include the land elevation and the seafloor depth.In general, depending on the surface features over the land, ICESat-2 photons may return from the bare earth as well as from the canopy/built-up features 85,103,104 .Generally, photons returning from the ground are denser in the horizontal direction than in the vertical direction.In contrast, the photons returning from the canopy have a high density in the vertical direction (refer to Fig. 6c).Using this as a hint, algorithms like DBSCAN can easily distinguish the photons returning from the ground and canopy 105,106 .In this research, we have classified the photons returned from the ground and surface features using DBSCAN, and finally, towards generating the DBEM for the study area, only those photons that returned from the bare earth were considered at the land parts, and the photons falling on the canopy/built-up were discarded (refer to Fig. 6c,d).
Approximately 0.2 million ICESat-2 photons representing the values of the seafloor depths and land elevations were collected as a part of the data aggregation process; from this, a few well-distributed points were reserved as checkpoints for technical validation of the output (discussed in the subsequent sections).To generate a continuous surface, the point database containing the accumulated depths and the land elevation values from the ICESat-2 photons must undergo interpolation 107 .Interpolation is a mathematical process to predict unknown values using the surrounding measured values 108 .However, the pattern of data collection through ICESat-2 is along-track, and solely using the along-track points in the interpolation methods may not yield a continuous surface as gaps between individual tracks impact the point distribution.The maximum gap in the study area is ~400 m in our collection of along-track points.Thus, additional points are needed to compensate for the fill-gap issue.For this, ~4.9 thousand points of depth values were extracted from available Electronic Navigation Charts (ENCs) 109,110 mentioned in Table 1.The available ENCs for the study area accord with the International Hydrographic Organization (IHO) standard called S-57, which is a data transfer standard used for the exchange of digital hydrographic data between national hydrographic offices and its distribution to manufacturers, mariners, and other data users 111 .These ENCs are processed using the tools available in ESRI's ArcGIS software 112 .Earlier, the procedure for using ENCs in GIS platforms was described by Hui et al. 113 and Taylor et al. 114 .The procedure to retrieve the sounding depth values and the associated geodetic coordinates from the ENCs includes reading the values from a data layer titled Geographic Reference Record.During the retrieval of depth values from the ENCs, it was ensured, especially for the gap areas, that a minimum of 2 points were retrieved successfully.A few points from the pool of sounding depths acquired from the ENCs were reserved as checkpoints for technical validation.
Similarly, for the extent of having the land part (like Rameshwaram and Talai Mannar Islands), orthometric elevation values from an open-access bare-earth model called Forest And Buildings removed Copernicus DEM (FABDEM) 115 were considered.FABDEM, available at 30 m spatial resolution, is the first kind of global DEM to represent the elevations of near bare-earth 116 .Figure 7 shows the distribution of total points considered from various sources used to generate DBEM for the study area.The details about the data sources are mentioned in Table 1. Figure 8 shows the schematic representation of the methodology implemented towards generating DBEM for Adam's Bridge; the methodology is partly considered from the earlier researchers' works 100 that have attempted to derive bathymetry from the ICESat-2 geolocated photons.
Interpolation is essential in generating bathymetric surfaces from well-distributed points; the process enables estimating the values in areas lacking direct measurements.Most GIS software has implemented various interpolation methods 117 .However, there is no single choice in selecting the optimal interpolation technique, as most of them can significantly impact the accuracy based on the surface type, data density, and point distribution 107 .Interpolation techniques are generally classified into two categories: deterministic and geostatistical 117 .Deterministic interpolation techniques create surfaces based on measured points or mathematical formulas; methods such as Inverse Distance Weight (IDW), natural neighbour, and spline/trend fitting fall in the deterministic category.In contrast, geostatistical interpolation techniques, such as various forms of Kriging, are based on statistics and are used for more advanced prediction surface modeling that also includes some measure of the certainty or accuracy of predictions 117 .In our experiment, we generated surfaces using various methods available in the GIS software, namely, ESRI's ArcGIS Pro 10.3 version (http://www.esri.com); the methods are listed in Table 2. Technical validation for all the surfaces generated using the available interpolation methods was performed using the reserved checkpoints described earlier.The study considered the DBEM with less Root Mean Square Error (RMSE) obtained with the IDW method.Earlier, researchers mentioned that IDW is considered a highly adaptable estimation method as it is the best to reconstruct natural surfaces given dense and well-distributed points in www.nature.com/scientificdata/ the study area 118,119 .As in our case, high-density and well-distributed points for the study area are available, and due to this, the IDW method generated a superior DBEM compared to other methods (refer to Table 2).
Determining the resolution (i.e., the cell/grid size) of the bathymetric elevation model depends on the density of the points used to generate the model.Attempting to generate a high-resolution surface model using less-density points will accumulate artifacts in the output; thus, the higher the density of points, the fewer the artifacts in the high-resolution surface model.Hu 120 , through his experiments, suggested a standard method for determining the cell size of the raster-based bathymetry/elevation model through Eq. 2.

Data source
Remarks/Usage Details of data availability
• Preferred acquisitions: January to May (low turbid period) and mostly night-time acquisitions.
• Refraction-corrected photons returned from the seafloor were used to compute the bathymetry.
• Elevations from the bare earth returned photons were used to estimate the digital elevation model for the land part.www.nature.com/scientificdata/ Where s is the estimated raster cell size, and n is the number of points in the minimum area of density (A) within the extent of the point distribution.From the point distribution map obtained in our experiment, it is observed that for every 200 sq.m, a minimum of 2 points representing seafloor depths/land elevations exist.Thus, the output pixel's cell size was kept to 10 m during the interpolation stage based on Eq. 2.

Data Records
The DBEM for Adam's Bridge generated from this study is accessible at figshare 121

technical Validation
Technical validation of the DBEM generated in this research was performed using qualitative and quantitative methods.Qualitative assessment of any modeled surface (like a Digital Elevation Model -DEM, Digital Surface Model -DSM, Digital Terrain Model -DTM, and DBEM) can be done using visual methods in a three-dimensional (3D) viewer, which in general will be available in satellite-based image processing or GIS software having advanced capabilities.However, the use of visual methods for quality checking depends on the expertise and experience of the human resource 122 .3D perspective-based visual analysis of digital surface features is highly useful for perceiving intricate details of seabed topography 123 .Additionally, 3D perspective views with directional lighting effects and coloring schemes enhance the impressions of relief data 124 .
During the visual analysis to assess the quality of the DBEM, it was compared with the free and open accessible global bathymetric data sources like GEBCO (latest version titled GEBCO_2023 grid) and GMRT bathymetric datasets.Figure 9 illustrates 3D perspective views of the Adam's Bridge and its surroundings that were generated using two free and open accessible global bathymetric datasets and the DBEM generated in this research.GEBCO_2023 grid, which is of 450 m spatial resolution, exhibited a relatively flat surface, especially at the extent of Adam's Bridge (refer to Fig. 9a); the reason may be attributed to its coarser resolution, which resulted in giving less details with respect to the seabed topography.Similarly, the digital bathymetry from GMRT, which is available at 100 m spatial resolution, not only failed to give a relief impression of the Adam's Bridge but also accounted for artifacts in the form of spikes and large sinks (refer to Fig. 9b); the reasons for the errors can be because GMRT lacks any sounding data near the vicinity of Adam's Bridge 72 .Visual impressions from the DBEM generated from our research with 10 m spatial resolution exhibited a higher degree of detail for the entire study area than the other two bathymetric models.Also, it has resulted in a more realistic representation of the Adam's Bridge as a submarine continuation of Dhanushkodi and Talaimannar Island (refer to Fig. 9c).Moreover, at regular intervals of Adam's Bridge, sudden narrow channels with depths varying between 2 to 3 m are seen which are not evident in the other bathymetric datasets.These narrow channels permit the exchange of water waves between the Gulf of Mannar and the Palk Strait.Importantly, from the crest line of Adam's Bridge, the narrow channels are accompanied by perpendicular ridges (refer to Fig. 9c), which are nullified in the GEBCO and GMRT bathymetric datasets.
By applying the processing scheme mentioned in Fig. 8, we generated the high-resolution DBEMs for the Adam's Bridge using various interpolation methods, in which ~0.2 million points are from the ICESat-2 photons contributing to the surface generation.During this process, 400 points of depth and elevation values from ICESat-2 photons were reserved as checkpoints towards quality checking of the output DBEMs (refer to Fig. 7).From these reserved points, 300 depth values were used to check the accuracy of the DBEM over the vicinity of Adam's Bridge (extent having water), and 100 points were used for quality checking over the extent containing the islands.Additionally, 200 checkpoints from ENCs were used to assess the accuracy for the extent of the Gulf of Mannar and Palk Strait, where the depth is more than 40 m (refer to Fig. 7).RMSE, a statistical formula, was used to quantify the vertical accuracies for the surfaces generated using various interpolation methods and is based on Eq.  number of observations.Table 2 summarizes the results with the model accuracies for the seafloor and terrains for all the surfaces generated using various interpolation methods.From Table 2, it is evident that the RMSE of the surface generated using IDW interpolation method performed better than all other interpolation methods with the error being less than 0.79 m over the extent of the study area.

Usage Notes
This dataset, available in GeoTIFF format, can be opened, visualized, and further used to derive additional terrain/surface characteristics (ex., slope, aspect, and contours.)with the help of satellite-image processing or GIS software.However, for 3D visualization, the software should be equipped to support 3D viewer.Free and open-source GIS software, namely QGIS, can be used to view this DBEM in 3D viewer.Alternatively, commercial-off-the-shelf (COTS) GIS software like ESRI ArcGIS (http://www.esri.com)can be used with advanced interactive features.Figure 10 is a typical view for the extent of Adam's Bridge generated using ArcScene module of ESRI ArcGIS software.Figure 11 shows 2D and 3D perspective view of the Adam's Bridge.Figure 12 shows a wire-mesh mode of this dataset for a zoomed extent of the Adam's Bridge generated in the ArcScene module of ESRI ArcGIS software (http://www.esri.com);here, the z-exaggeration was set to 200 times during the visualization to amplify the visual intricacies of Adam's Bridge structure.The proposed DBEM can

Fig. 1
Fig. 1 Adam's Bridge and its surroundings as seen in satellite imagery.(a) Map showing southern India and Sri Lanka.(b) The extent marked with the white box comprises part of India, including Pamban/Rameswaram Island, with its eastern point named Dhanushkodi, Adam's Bridge, and Talai Mannar, the western end of Mannar Island, Sri Lanka.(c) A zoomed extent showing the exposed sandbanks of Adam's Bridge.The satellite imagery used in the maps is from the web mapping services of the Sentinel-2 cloudless layer for 2021 by EOX (https://s2maps.eu/and https://esa.maps.eox.at/).

Fig. 2
Fig.2Available reference ground tracks of ICESat-2 mission over the extent of Adams' Bridge.The satellite imagery used in the map is from web mapping services of the Sentinel-2 cloudless layer for 2021 by EOX (https://s2maps.eu/and https://esa.maps.eox.at/), and the source of ICESat-2 ground tracks is https://icesat-2.gsfc.nasa.gov/science/specs.

Fig. 3
Fig. 3 Beams of ICESat-2 acquisitions.(a) Six beam configuration of ICESat-2 organized as three pairs -left, near-nadir, and right.Within each pair, one strong beam and a weak beam have a separation of 90 m between them.The left and right paired beams are 3.3 km apart from the nadir most paired beam.(b) Qualified beams (both strong and weak) over the extent of Adam's Bridge used to generate a digital bathymetric elevation model.

Fig. 4
Fig. 4 Illustration of ICESat-2's data acquisition over water bodies.(a) A subset of ICESat-2 single beam acquired over the Adam's Bridge.(b) The two-dimensional profile generated from the ICESat-2 photons for an extent of Adam's Bridge -primarily, the photons were returned from the water surface, water column, and seafloor.A small patch of exposed sandbank is also evident from the profile.(c) 2-dimensional profile showing classified photons.Differences in depths computed from the photons of apparent seafloor and refractioncorrected seafloors can be perceived.

Fig. 5
Fig. 5 Illustration of seafloor variation at Adam's Bridge using the profiles from ICESat-2 photons.(a) ICESat-2 beam acquired over the Adam's Bridge (at the tail-end of Dhanushkodi) and overlaid on the satellite data.(b) Subset of the ICESat-2 beam acquired over the Adam's Bridge (at the head of Talai Mannar) and overlaid on the satellite data.(c) 2-dimensional profile showing the refraction-corrected seafloor for the beam corresponding to (a).(d) 2-dimensional profile showing the seafloor variation for the beam corresponding to (b).Both the profiles shown in (c) and (d) signify the presence of Adam's Bridge's structure, which gradually rises from a depth of nearly 8 m to the surface of the water level.

Fig. 6
Fig. 6 Illustration of ICESat-2 photons acquired over the extent of land and seawater surrounding Adam's Bridge.(a) ICESat-2 beam acquired over Rameshwaram Island and overlaid on the satellite data.(b) ICESat-2 beam acquired over Dhanushkodi and overlaid on the satellite data.(c) 2-dimensional profile showing the seafloor and the land/canopy for the beam shown in (a).(d) 2-dimensional profile showing the seafloor and the exposed land part at the Dhanushkodi with respect to the beam shown in (b).

Fig. 7
Fig. 7 Map showing the distribution of points representing the seafloor depths and land elevations over the Adam's Bridge and its surroundings.Nearly 0.2 million points were accrued towards creating a database primarily acquired from the ICESat-2 geolocated photons returned from the seafloor and land, sounding depths from ENCs, and land elevations from a digital bare-earth model, FABDEM.Few points that were reserved as checkpoints for technical validation of the output are also shown in the map.

Fig. 8
Fig. 8 Schematic representation of methodology implemented to generate a digital bathymetric elevation model for Adam's Bridge and its surroundings using ICESat-2 photon data.

Fig. 9
Fig. 9 Comparison of free and open-access bathymetric data sources with the high-resolution bathymetric data generated using ICESat-2 seafloor returned photons.(a,b) GEBCO_2023 grid with a spatial resolution of 450 m and GMRT bathymetry with a spatial resolution of 100 m, respectively.(c) Digital Bathymetric Elevation data of 10 m resolution generated using ICESat-2 seafloor returned photons.These 3D perspectives views were generated using ESRI's ArcScene Ver.10.8.1 software (https://www.esri.com/).

3 .H
(Depth or Elevation) (Depth or Elevation) In Eq. 3, (Depth or Elevation) DBEM is the set of depth or elevation values obtained from the modeled DBEMs and (Depth or Elevation) checkpoints are the set of depth or elevation values of the reserved checkpoints.n is the

Fig. 10 A
Fig. 10 A 10 m Digital Bathymetric Elevation Model generated using ICESat-2 photons for the extent of Adam's Bridge and its surroundings.This 3D perspective view was generated using the ArcScene module of ESRI ArcGIS Ver.10.8.1 software (https://www.esri.com/).

Fig. 11
Fig. 11 3D perspective generated with a 10 m Digital Bathymetric Elevation Model (DBEM) generated using ICESat-2 photons for the extent of Adam's Bridge.(a) Extent showing the Adam's Bridge in high-resolution satellite imagery.The satellite imagery used in the map is from the web mapping services of the Sentinel-2 cloudless layer for 2021 by EOX (https://s2maps.eu/and https://esa.maps.eox.at/).(b) Perspective view for the extent of Adam's Bridge and its surroundings.This 3D perspective view was generated using the ArcScene module of ESRI ArcGIS Ver.10.8.1 software (https://www.esri.com/).

Fig. 12
Fig.12 Wire-mesh mode of the 10 m Digital Bathymetric Elevation Model generated using ICESat-2 photons for the extent of Adam's Bridge and its surroundings.This 3D perspective view was generated using the ArcScene module of ESRI ArcGIS Ver.10.8.1 software (https://www.esri.com/).

Table 1 .
• ICESat-2 Approximately 0.2 million ICESat-2 photons representing the seafloor's depths and terrain elevations were collected as a part of data collection.Details These charts are issued by hydrographic offices, as per the International Hydrographic Organisation's (IHO) standards, specifications, and symbol sets.Towards densification of elevation points for the extent of the area having the land part (like Rameshwaram and Talai Mannar Islands), elevation values from FABDEM were used.https://data.bris.ac.uk of the data sources used to generate a digital bathymetric elevation model for the Adam's Bridge.
. The DBEM is in GeoTIFF file format with a cell size of 10 m (spatial resolution) projected in the Universal Transverse Mercator (UTM) coordinate system with the Zone 44 N [EPSG:32644].The pixel values of the DBEM represent orthometric depth values for the extent with water and elevation values for the extent with peninsula/land/islands in meters.

Table 2 .
Details of accuracies of the surfaces generated using various interpolation method.