U.S. Geological Survey 2016 1.0 raster digital data Open-File Report 2016-1168 Woods Hole Coastal and Marine Science Center, Woods Hole, Massachusetts U.S. Geological Survey, Coastal and Marine Geology Program https://doi.org/10.3133/ofr20161168 http://pubs.usgs.gov/of/2016/1168/GIS_catalog/raster/bathy/2013-003-FA_Bath5m.zip http://pubs.usgs.gov/of/2016/1168/ofr2016-1168-data_catalog.html Seth D. Ackerman Laura L. Brothers David S. Foster Brian D. Andrews Wayne E. Baldwin William C. Schwab 2016 1.0 Open-File Report 2016-1168 Reston, VA U.S. Geological Survey http://dx.doi.org/10.3133/ofr20161168 http://pubs.usgs.gov/of/2016/1168/ These data were collected under a cooperative agreement between the Massachusetts Office of Coastal Zone Management (CZM) and the U.S. Geological Survey (USGS), Coastal and Marine Geology Program, Woods Hole Coastal and Marine Science Center (WHCMSC). Initiated in 2003, the primary objective of this program is to develop regional geologic framework information for the management of coastal and marine resources. Accurate data and maps of sea floor geology are important first steps toward protecting fish habitat, delineating marine resources, and assessing environmental changes due to natural or human impacts. The project is focused on the inshore waters of coastal Massachusetts, primarily in water depths of 5 to 30 meters (m) deep. Data collected for the mapping cooperative have been released in a series of USGS Open-File Reports (http://woodshole.er.usgs.gov/project-pages/coastal_mass/). The geophysical data were collected during a survey in 2013 during USGS Field Activity 2013-003-FA (http://cmgds.marine.usgs.gov/fan_info.php?fa=2013-003-FA) and cover approximately 185 square kilometers of the inner continental shelf. These data are used to define sea-floor morphology as part of the Massachusetts Seafloor Mapping Project. This gridded bathymetry layer and depth-colored hillshaded image represents approximately 185 square kilometers of interferometric sonar (Systems Engineering and Assessment Ltd. SWATHPlus-M) bathymetry data collected by the U.S. Geological Survey during USGS survey 2013-003-FA in 2013 offshore of Massachusetts, south of Martha's Vineyard and north of Nantucket. An additional approximately 50 square kilometers of interferometric bathymetry data was collected north of Nantucket during the same survey. 20130522 20130611 ground condition of survey dates: 20130522-20130611; see Completeness_Report for more information Complete None planned -70.810609 -70.029975 41.390816 41.192901 USGS Metadata Identifier USGS:a2c7d2f2-3ba2-4dad-967f-1b34f8a184ab None U.S. Geological Survey USGS Coastal and Marine Geology Program CMGP Woods Hole Coastal and Marine Science Center WHCMSC seafloor sea floor Marine Geology field activity 2013-003-FA field activity number 13003 USGS CMG InfoBank ID U-03-13-MA raster GeoTIFF bathymetry hillshaded bathymetry shaded relief bathymetry swath bathymetry interferometry SEA SWATHPlus SEA Submetrix M/V Scarlett Isabella ISO 19115 Topic Category oceans elevation imageryBaseMapsEarthCover geoscientificInformation USGS Thesaurus bathymetry sea-floor characteristics interferometric sonar marine geophysics geospatial datasets None North America United States Atlantic Ocean Massachusetts Marthas Vineyard Martha's Vineyard Nantucket Nantucket Sound Cape Cod None Seafloor Sea floor Surface None 2013 None These data are NOT to be used for navigation. Mariners should refer to the appropriate nautical chart. Public domain data from the U.S. Government are freely redistributable with proper metadata and source attribution. Please recognize the U.S. Geological Survey as the originator of the dataset. Seth Ackerman U.S. Geological Survey Geologist mailing and physical address

384 Woods Hole Rd.

Woods Hole MA 02543-1598 USA 508-548-8700 x2315 508-457-2310 sackerman@usgs.gov http://pubs.usgs.gov/of/2016/1168/GIS_catalog/browse_img/big/2013-003-FA_Bath5m.jpg Depth-colored hillshade image of bathymetry JPEG Microsoft Windows Vista Version 6.1 (Build 7601) Service Pack 1; ESRI ArcCatalog 9.3.1.4095 This bathymetric surface represents interpolated data at 5-meter per pixel resolution. Bathymetric data were processed to account for gaps that occurred along-track and between adjacent lines. Quality control and data processing were conducted to remove spurious points and reduce speed of sound artifacts (refraction). Soundings from the surveys were processed and edited using Computer Aided Resource Information System (CARIS) Hydrographic Information Processing System (HIPS; vers. 8.0, 8.1 and 9.0). Although the soundings were edited, small data spikes may still exist. While much effort was devoted to cleaning the data and minimizing survey artifacts, some artifacts may still be seen in the final gridded bathymetric surface including an along track artifact at nadir (the ship track) and at the far edges of a survey line where adjacent swaths overlap. These artifacts are especially noticeable in areas of little local relief. Nearly all quality main survey line bathymetric data (generally shore parallel) that were collected within the survey areas south of Martha's Vineyard and north of Nantucket were incorporated in this bathymetric surface. Tielines and transit lines were generally excluded from this bathymetric surface, especially if equal or higher quality main survey line data were available. The tieline and transit line data were used to verify tide corrections during post-processing. Line L1f1 was heavily edited due to navigation errors and issues with staying on the survey line. The line was restarted as L1f2, and much of the end of line L1f1 was discarded in swath editing. Time gaps occurred during this field activity where no surveying was done due to bad weather, rough seas and stops in port for crew changes. Planned surveying was paused and no main survey line data were collected during these times: >JD143 15:20 - JD146 20:06 (no data were collected on JD144 or JD145) >JD149 15:40 - JD151 18:59 (no data were collected on JD150) >JD153 00:01 - JD153 14:32 >JD155 18:16 - JD156 02:06 >JD158 06:12 - JD160 19:15 (no data were collected on JD159) USGS Field Activity 2013-003-FA: Navigation was acquired with a Coda Octopus F185R Real-time Kinematic Global Positioning System (RTK-GPS); which is accurate to less than 1 meter (typically on the order of centimeters) horizontally. The SWATHPlus transducers were mounted on a rigid pole, approximately 2.17 m below the water line, along the port side of the M/V Scarlett Isabella. The RTK-GPS antenna was located on the same pole approximately 4.6 m above the sea surface. These offsets were corrected within the SWATHPlus acquisition software. USGS Field Activity 2013-003-FA: Vertical accuracy of the raw data based on system specifications may approximate 1% of water depth, 0.1 to 0.4 meters within the survey area. However, overall vertical accuracies on the order of 50-cm or better are assumed based on the following considerations: The Coda Octopus F185R Attitude and Positioning system, used to correct for vessel roll, pitch, heave, and yaw, has a theoretical vertical accuracy of a few mm. Tidal offsets were corrected to MLLW using Real-Time Kinematic GPS (RTK-GPS) heights in post-processing (see process steps). Corrections were sent via cellular modem and VHF radio to the ship from the base station on land. The primary reference station is located at the USGS Marine Operations Facility (MOF) in Falmouth, MA. USGS field tests (unpublished) using submerged targets suggest that the vertical accuracy of the RTK-GPS tidal correction is less than 30cm. Refraction artifacts were minimized by acquiring a range of sound velocity profiles with an ODIM MVP30 moving vessel profiler (until ~0550UTC on JD 156) and a hand-casted AML Oceanographic MinosX sound velocimeter (for the duration of survey 2013-003-FA). Sound velocity data are entered into the SWATHPlus acquisition software to model the sound velocity structure of the water column. Gridding algorithms and cell sizes for these data could introduce errors as great as 3m along the edges of the data, but gridding-induced errors are likely much smaller (less than 20-cm) for most of the survey area. Changes in ship draft due to water and fuel usage were not considered. U.S. Geological Survey Unpublished Material disc 20130522 20130611 ground condition during the survey dates: 20130522-20130611; see Completeness_Report for more information raw data USGS Field Activity 2013-003-FA: A 234 kHz Systems Engineering and Assessment Ltd. (SEA) SWATHplus interferometric sonar was mounted on the port side of the M/V Scarlett Isabella. Data were acquired during May and June 2013. USGS Configuration: Survey: Survey lines were run at an average speed of 5 knots and were generally spaced 100 m apart to obtain overlapping swaths of data and full sea-floor coverage. Sonar: The SEA SWATHplus-M operates at a frequency of 234 kHz and a variable range (increased or decreased manually depending on water depth). The system was operated with a transmit power: 80 percent, a transmit length: 43 cycles, 4096 samples per channel, ping rate of 15 per second. Speed of Sound: Sound-velocity profiles were collected approximately every 2 hours using an ODIM MVP30 moving vessel profiler or a hand-casted AML Oceanographic MinosX sound velocimeter. Tides: Tides were corrected during post-processing using two-minute interval averages of the RTK tide data recorded in the "TID" message string of the HYPACK navigation file for each survey line. VDatum transformation grid “MENHMAgome01_8301” was used within the HYPACK software for the kinematic tide correction process during 2013-003-FA to offset water level heights from the WGS84 ellipsoid to MLLW. See the HYPACK navigation files metadata for additional information "2013-003-FA_HYPACK_meta.xml" available from the larger work citation. Remaining artifacts: While much effort was devoted to cleaning the data and minimizing survey artifacts, some artifacts may still be seen in the final gridded bathymetric surface including an along track artifact at nadir (the ship track) and at the far edges of a survey line where adjacent swaths overlap. These artifacts are especially noticeable in areas of little local relief. 1. Raw to Processed conversion for USGS field activity 2013-003-FA: Each raw SWATHPlus bathymetric sonar file (SXR) was converted to SWATHPlus processed file (SXP) using SEA SWATHPlus Swath Processor (ver. 3.07.17). During the conversion process, sound velocity profiles were used to minimize potential refraction artifacts from fluctuations in the speed of sound within the water column. Several bathymetric filters were applied to remove erroneous soundings and reduce the density of the data. Predicted tidal information was merged into the processed file (SXP) during this conversion, however final tidal corrections were applied later in the processing flow. Bathymetric filtering typically included low amplitude (100%), range (0-4m), box (3-50m depth, 1.5-75m horizontal), median (window size 5), alongtrack 1 (depth difference of 5-m, window size 5-m, and learn rate of 0.7), alongtrack 2 (depth difference of 1.5-m, window size 1m, and learn rate of 0.9), and mean filters (0.25m). These represent the typical parameters used for the majority of the survey but slight adjustments were required for some areas. This process step and all subsequent process steps were performed by the same person - Seth Ackerman. 2013 Seth Ackerman U.S. Geological Survey Geologist mailing and physical address

384 Woods Hole Rd.

Woods Hole MA 02543-1598 USA 508-548-8700 x2315 508-457-2310 sackerman@usgs.gov 2. CARIS processing for USGS field activity 2013-003-FA: A new CARIS HIPS project (ver. 8.0) was created for this field activity with projection information set to Universal Transverse Mercator (UTM) Zone 19, WGS84. Each SWATHPlus processed file (sxp) was imported to the new CARIS project using the Import/Conversion Wizard. A 5 meter resolution Bathymetric and Statistical Error (BASE) Surface was created from the files for each Julian day. The BASE surface for each day was reviewed for any inconsistencies or data anomalies. Navigation was reviewed and edited as needed using the navigation editor tool. Beam-to-beam slopes and across track angle filters were applied to the soundings line by line. The refraction editor was used to adjust sound speed values in areas where velocimeter data did not adequately correct depth profiles obviously influenced by local anomalies in speed of sound through the water column. Preliminary processing was done during the survey and additional processing was done post-survey with CARIS HIPS and SIPS versions 8.0, 8.1 and 9.0. This process step was done between 2013 and 2015. 2015 3. Prepare and apply RTK Tide data for USGS field activity 2013-003-FA: A tidal correction file was created during post-processing by applying a two-minute moving median filter to the RTK tide data recorded in the "TID" message string of the raw HYPACK navigation files. The TID messages from device 1 (DEV 1; which was the NovAtel DL-V3 GPS receiver) were extracted using a shell script called web_novTid_movmed and a Python script called tid_movmed.py. Only HYPACK extracted tide data were used for this survey. RTK tide corrections were loaded and applied to the soundings line by line in CARIS, resulting in bathymetric values referenced to MLLW. Each line was re-merged and the preliminary BASE surfaces were recomputed to incorporate the new tide information. Additional swath editing was done in CARIS using the swath and area based editors to remove any additional inconsistencies and minimize survey artifacts. This process step was done between 2013 and 2015. 2015 2013_003MasterTide.tid 4. Final BASE surface creation and export: Three separate 5-meter BASE surfaces were created for the (1) east-west trending lines south of Martha's Vineyard; (2) the north-south trending lines south of Martha's Vineyard and east of Noman's Land; and (3) the area north of Nantucket utilizing all the quality bathymetric soundings data from the survey. An interpolation algorithm was run on each of these BASE surfaces in CARIS HIPS/SIPS to fill small gaps between lines. The BASE surfaces were then combined in CARIS HIPS using the Tools -> Surfaces -> Combine tool (overlapping areas results in the shoalest depth in the resulting surface). The final surface was then exported from CARIS as 32-bit GeoTIFF and as hillshaded GeoTIFF image. Except for setting the sun angle to 40 degrees above the horizon and the vertical exaggeration to 5x, the default hillshade settings were used to create the hillshaded GeoTIFF image. This process step was done between 2015 and 2016. 2016 Bathymetric images (Gridded bathymetric GeoTIFF image - 2013-003-FA_Bath5m_MLLW.tif and depth-colored bathymetric shaded-relief image - 2013-003-FA_BathHlshd5m.tif) USGS Thesaurus keywords added to the keyword section. 20180720 U.S. Geological Survey VeeAnn A. Cross Marine Geologist Mailing and Physical

384 Woods Hole Road

Woods Hole MA 02543-1598 508-548-8700 x2251 508-457-2310 vatnipp@usgs.gov Crossref DOI link was added as the first link in the metadata. 20191118 U.S. Geological Survey VeeAnn A. Cross Marine Geologist Mailing and Physical

384 Woods Hole Road

Woods Hole MA 02543-1598 508-548-8700 x2251 508-457-2310 vatnipp@usgs.gov Added keywords section with USGS persistent identifier as theme keyword. 20200908 U.S. Geological Survey VeeAnn A. Cross Marine Geologist Mailing and Physical

384 Woods Hole Road

Woods Hole MA 02543-1598 508-548-8700 x2251 508-457-2310 vatnipp@usgs.gov Raster Pixel 4183 13008 Universal Transverse Mercator 19 0.999600 -69.000000 0.000000 500000.000000 0.000000 row and column 5.000000 5.000000 meters D_WGS_1984 WGS_1984 6378137.000000 298.257224 Mean lower low water 0.1 meters Explicit depth coordinate included with horizontal coordinates There are no attributes associated with a GeoTIFF however the pixel values of the 32-bit image are depth values. The "no data" value is set to 3.402823e+038. Data values represent depth in meters referenced to mean lower low water (MLLW). The color-ramp of the hillshade bathymetry image is full visual spectrum from red to purple representing shallow (less than 10 meters) to deep (approximately 30 meters), respectively. The "no data" value in the hillshaded image is RGB value: 0,0,0. U.S. Geological Survey Seth Ackerman U.S. Geological Survey Geologist mailing and physical address

384 Woods Hole Rd.

Woods Hole MA 02543-1598 USA 508-548-8700 x2315 508-457-2310 sackerman@usgs.gov USGS data release 2013-003-FA bathymetry from south of Martha's Vineyard and north of Nantucket (2013-003-FA_Bathy.zip). The zip file contains a folder with the the 32-bit Floating Point GeoTIFF image (2013-003-FA_Bath5m_MLLW.tif) with a world file (2013-003-FA_Bath5m_MLLW.tfw). A depth-colored hillshade bathymetry image (2013-003-FA_BathHlshd5m.tif) with a world file (2013-003-FA_BathHlshd5m.tfw) and metadata files in three standard formats are also included in the zip file. The pixel values in the depth-colored hillshade bathymetry image DO NOT correspond to depths. The pixel values of the 32-bit Floating Point GeoTIFF image are depth values. Neither the U.S. Government, the Department of the Interior, nor the USGS, nor any of their employees, contractors, or subcontractors, make any warranty, express or implied, nor assume any legal liability or responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, nor represent that its use would not infringe on privately owned rights. The act of distribution shall not constitute any such warranty, and no responsibility is assumed by the USGS in the use of these data or related materials. Any use of trade, product, or firm names is for descriptive purposes only and does not imply endorsement by the U.S. Government. GeoTIFF CARIS HIPS/SIPS 9.0 GeoTIFF 32-bit floating point GeoTIFF image and depth-colored hillshaded GeoTIFF image The zip file contains a folder with a 32-bit GeoTIFF floating point GeoTIFF image at 5-meter resolution, a depth-colored hillshaded GeoTIFF image and the associated metadata and world files. Use 7zip, WinZip, Peazip, pkUnzip, Keka or other uncompressing software 41 http://pubs.usgs.gov/of/2016/1168/GIS_catalog/raster/bathy/2013-003-FA_Bath5m.zip http://pubs.usgs.gov/of/2016/1168/ofr2016-1168-data_catalog.html http://dx.doi.org/10.3133/ofr20161168 Data can be downloaded via the World Wide Web (WWW) None The ZIP file contains a 32-bit Floating Point GeoTIFF image compressed using the DEFLATE compression format with a world file. To utilize these data an image processing or GIS software package capable of viewing a 32-bit Floating Point GeoTIFF image with DEFLATE compression. Standard image viewing software cannot translate a 32-bit floating point GeoTIFF image. The depth-colored hillshade GeoTIFF image is 32-bit GeoTIFF with LZW compression which can be viewed in any image processing or GIS software package capable of viewing a 32-bit GeoTIFF image with LZW compression. 20240318 Seth Ackerman U.S. Geological Survey Geologist mailing and physical address

384 Woods Hole Rd.

Woods Hole MA 02543-1598 USA 508-548-8700 x2315 508-457-2310 whsc_data_contact@usgs.gov The metadata contact email address is a generic address in the event the person is no longer with USGS. (updated on 20240318) FGDC Content Standards for Digital Geospatial Metadata FGDC-STD-001-1998 local time