Coastal Multibeam Bathymetry and Backscatter Data Collected in May 2023 from Rockaway Peninsula, New York

Acquisition: An ERS consisting of multibeam bathymetry and acoustic backscatter was conducted using two Teledyne Reson SeaBat T50-P multibeam echosounders, in dual-head configuration. The pair of Mills Cross (configuration) transmit and receive arrays, were placed side-by-side within a bracket that oriented them at opposing 30-degree angles (relative to horizontal). The bracket was pole-mounted on the starboard gunnel of the research vessel (R/V) Sallenger, a 26-foot fiberglass vessel. The sonar arrays were oriented athwartships (primary and secondary arrays facing outward and down to port and starboard, respectively) and located approximately 1-m below the waterline when deployed. An Applanix Pos MV Wavemaster II V5 inertially aided navigation system was used for timing, positioning, and attitude. The Inertial Measurement Unit (IMU) was mounted atop the sonar bracket just aft of the pole. Two Trimble GA830 Global Navigation Satellite System (GNSS) antennas were affixed atop the cabin approximately amidships of the vessel, two meters apart, port-starboard. An AML Oceanographic Micro-X sound velocity sensor was used for speed of sound readings at the sonar head. A Sontek CastAway CTD was used for measuring sound velocity profiles. Hydrographic survey acquisition and processing software, HYPACK/HYSWEEP (version 2022) managed data inputs from sensors and produced two files (.RAW and .HSX). Data were collected in WGS84 [G1150] Universal Transverse Mercator (UTM) Zone 18 North (18N) coordinates.

Processing: Positioning and attitude data were post-processed in Applanix POSPac MMS (POSPac) software version 8.8 Service Pack 1. Raw position files were imported into POSPac and processed with the IN-Fusion PP-RTX processing mode utilizing Trimble CenterPoint RTX technology, creating smoothed best estimate of trajectory (SBET) files and dynamic error files (SMRMSG) for each day of the survey. Bathymetric data were post-processed in CARIS Hydrographic processing software HIPS/SIPS version 11.4.28. HYPACK acquisition files (.HSX) were first imported into CARIS and then post-processed navigation files (SBET and SMRMSG) were imported and applied. Additionally, Global Positioning System (GPS) tide (GPS tidal corrections utilize GPS height information from real time survey observations to obtain an ellipsoid height) was computed and sound velocity profiles were applied. Roll calibration was also performed, as necessary. TPU was calculated from real time values recorded in the HSX and the SMRMSG files. Data were filtered to exclude any points with RESON Quality flags of 0, 1, and 2 (these flags are assigned to individual soundings at the time of acquisition and are a metric of pass/fail in regard to both brightness and collinearity: 0 = fail both, 1 = pass brightness only, 2 = pass collinearity only), as well as soundings that had a TPU greater than stated IHO Special Order standards (IHO, 2008). Refraction artifacts were addressed using Refraction Editor. Further editing was performed in Subset Editor. A 1-m bathymetric CUBE surface was computed and a visual and statistical inspection of the data effecting the surface was performed and edited as needed to serve as quality control (Calder and Wells, 2007). Small gaps in the continuous surface were caused by schools of baitfish obscuring the signal during acquisition as well as limitations in swath width and planned line navigation. The final surface was exported as an x,y,z ASCII file using the Base Surface to ASCII function.

Backscatter Mosaic Creation: Acoustic Intensity Backscatter data (.7k) were loaded in CARIS Hydrographic processing software HIPS/SIPS version 11.4.28 using the SIPS Backscatter tool. The final CUBE bathymetry surface serves as input for geometric corrections and average temperature and salinity data (derived from the Sontek CastAway CTD measurements taken during acquisition) are utilized by the program to calculate local absorption and transmission loss. These inputs, along with the following settings were entered into the tool, and then processing was initiated [Weighted Average = blending, image type = beam average, weighting = swath, and an Angle-Varying gain = 200]. During processing, the tool automatically creates a pattern file (.bbp) which is used to correct for variations in the beam pattern of the system. Upon completion, a 1-m backscatter mosaic with intensity values ranging from -42.2189 to 20.5824 was created. After visual inspection, the mosaic was exported as a 32-bit float GeoTIFF (.tiff) using the 'surface to GeoTIFF' function.

Datum Transformation: Final x,y,z data were converted using NOAA's VDatum v.4.3 software conversion tool from the WGS84 [G1150] realization to NAD83 [2011] NAVD88, using the GEOID12B model (orthometric height).