Esri Binary grid of the bathymetry of Indian River Bay, Delaware, generated from fathometer data acquired in April 2010 during U.S. Geological Survey Field Activity 2010-006-FA (IRB_BATHY, UTM, Zone 18, WGS 84)

Metadata also available as - [Questions & Answers] - [Parseable text] - [XML]

Metadata:

Identification_Information:
Citation:
Citation_Information:
Originator: VeeAnn A. Cross
Publication_Date: 2014
Title:
Esri Binary grid of the bathymetry of Indian River Bay, Delaware, generated from fathometer data acquired in April 2010 during U.S. Geological Survey Field Activity 2010-006-FA (IRB_BATHY, UTM, Zone 18, WGS 84)
Edition: 1
Geospatial_Data_Presentation_Form: raster digital data
Series_Information:
Series_Name: Open-File Report
Issue_Identification: 2011-1039
Publication_Information:
Publication_Place: Woods Hole Coastal and Marine Science Center, Woods Hole, MA
Publisher: U.S. Geological Survey, Coastal and Marine Geology Program
Online_Linkage: https://doi.org/10.3133/ofr20111039
Online_Linkage: Online_Linkage: Larger_Work_Citation:
Citation_Information:
Originator: V.A. Cross
Originator: J.F. Bratton
Originator: H.A. Michael
Originator: K.D. Kroeger
Originator: A. Green
Originator: E. Bergeron
Publication_Date: 2014
Title:
Continuous Resistivity Profiling and Seismic-Reflection Data Collected in April 2010 from Indian River Bay, Delaware
Edition: 1
Series_Information:
Series_Name: Open-File Report
Issue_Identification: 2011-1039
Publication_Information:
Publication_Place: Reston, VA
Publisher: U.S. Geological Survey
Online_Linkage: https://pubs.usgs.gov/of/2011/1039/
Online_Linkage: https://doi.org/10.3133/ofr20111039
Description:
Abstract:
A geophysical survey to delineate the fresh-saline groundwater interface and associated sub-bottom sedimentary structures beneath Indian River Bay, Delaware, was carried out in April 2010. This included surveying at higher spatial resolution in the vicinity of a study site at Holts Landing, where intensive onshore and offshore studies were subsequently completed. The total length of continuous resistivity profiling (CRP) survey lines was 145 kilometers (km), with 36 km of chirp seismic lines surveyed around the perimeter of the bay. Medium-resolution CRP surveying was performed using a 50-meter streamer in a bay-wide grid. Results of the surveying and data inversion showed the presence of many buried paleochannels beneath Indian River Bay that generally extended perpendicular from the shoreline in areas of modern tributaries, tidal creeks, and marshes. An especially wide and deep paleochannel system was imaged in the southeastern part of the bay near White Creek. Many paleochannels also had high-resistivity anomalies corresponding to low-salinity groundwater plumes associated with them, likely due to the presence of fine-grained estuarine mud and peats in the channel fills that act as submarine confining units. Where present, these units allow plumes of low-salinity groundwater that was recharged onshore to move beyond the shoreline, creating a complex fresh-saline groundwater interface in the subsurface. The properties of this interface are important considerations in construction of accurate coastal groundwater flow models. These models are required to help predict how nutrient-rich groundwater, recharged in agricultural watersheds such as this one, makes its way into coastal bays and impacts surface water quality and estuarine ecosystems. For more information on the survey conducted for this project, see https://cmgds.marine.usgs.gov/fan_info.php?fan=2010-006-FA.
Purpose:
The purpose of this dataset is to provide a tide corrected bathymetric grid for the Indian River Bay.
Time_Period_of_Content:
Time_Period_Information:
Range_of_Dates/Times:
Beginning_Date: 20100412
Ending_Date: 20100415
Currentness_Reference: ground condition
Status:
Progress: Complete
Maintenance_and_Update_Frequency: None planned
Spatial_Domain:
Bounding_Coordinates:
West_Bounding_Coordinate: -75.207119
East_Bounding_Coordinate: -75.060608
North_Bounding_Coordinate: 38.621180
South_Bounding_Coordinate: 38.565589
Keywords:
Theme:
Theme_Keyword_Thesaurus: USGS Metadata Identifier
Theme_Keyword: USGS:7991119f-4ffe-4b43-84c2-d30d742ec67b
Theme:
Theme_Keyword_Thesaurus: None
Theme_Keyword: U.S. Geological Survey
Theme_Keyword: USGS
Theme_Keyword: Coastal and Marine Geology Program
Theme_Keyword: CMGP
Theme_Keyword: Woods Hole Coastal and Marine Science Center
Theme_Keyword: WHCMSC
Theme_Keyword: Field Activity Number 2010-006-FA
Theme_Keyword: Info Bank ID K-6-10-DL
Theme_Keyword: bathymetry
Theme_Keyword: Lowrance GPS
Theme_Keyword: Esri binary grid
Theme_Keyword: Open-file Report 2011-1039
Theme:
Theme_Keyword_Thesaurus: ISO 19115 Topic Category
Theme_Keyword: elevation
Theme_Keyword: oceans
Theme:
Theme_Keyword_Thesaurus: USGS Thesaurus
Theme_Keyword: bathymetry
Theme_Keyword: sea-floor characteristics
Theme_Keyword: single-beam echo sounder
Theme_Keyword: geospatial datasets
Place:
Place_Keyword_Thesaurus: None
Place_Keyword: North America
Place_Keyword: North Atlantic Ocean
Place_Keyword: United States
Place_Keyword: Delaware
Place_Keyword: Indian River Bay
Place_Keyword: Indian River Inlet
Place_Keyword: Holts Landing
Place_Keyword: Piney Neck
Place_Keyword: Sussex County
Place_Keyword: Rosedale Beach
Access_Constraints: None.
Use_Constraints:
Not to be used for navigation. The 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.
Point_of_Contact:
Contact_Information:
Contact_Person_Primary:
Contact_Person: VeeAnn A. Cross
Contact_Organization: U.S. Geological Survey
Contact_Position: Marine Geologist
Contact_Address:
Address_Type: mailing and physical address
Address: Woods Hole Coastal and Marine Science Center
Address: 384 Woods Hole Rd.
City: Woods Hole
State_or_Province: MA
Postal_Code: 02543-1598
Contact_Voice_Telephone: (508) 548-8700 x2251
Contact_Facsimile_Telephone: (508) 457-2310
Contact_Electronic_Mail_Address: vatnipp@usgs.gov
Browse_Graphic:
Browse_Graphic_File_Name: Browse_Graphic_File_Description:
Thumbnail GIF image of bathymetry in Indian River Bay. The coastline is included for spatial reference.
Browse_Graphic_File_Type: GIF
Native_Data_Set_Environment:
Microsoft Windows Vista Version 6.1 (Build 7601) Service Pack 1; ESRI ArcCatalog 9.3.1.4095
Data_Quality_Information:
Logical_Consistency_Report: All the data were evaluated with the same criteria.
Completeness_Report:
This grid represents all the valid fathometer readings recorded by the HYPACK navigation software, plus some bathymetric readings recorded by the continuous resistivity profile system when the HYPACK system was not logging data.
Positional_Accuracy:
Horizontal_Positional_Accuracy:
Horizontal_Positional_Accuracy_Report:
The navigation system used was a Lowrance 480M with an LGC-2000 Global Positioning System (GPS) antenna. The antenna was located directly above the fathometer transducer mount point.
Vertical_Positional_Accuracy:
Vertical_Positional_Accuracy_Report:
All collected bathymetry values were collected by the 200 kHz Lowrance fathometer. The fathometer was mounted starboard side aft, directly below the GPS antenna and the resistivity streamer tow point. The Lowrance manufacturer indicates the speed of sound used by the system to convert to depths is 4800 feet/second. All values are assumed to be accurate to within 1 meter. The gridding, grid smoothing, and other processing affects the accuracy of individual cell values as these data were not collected with the intention of generating a grid and have areas of sparse coverage.
Lineage:
Process_Step:
Process_Description:
All the ASCII HYPACK navigation files collected were placed in their own folder on the computer based on day of collection. Each file has the following filename convention: LLL_TTTT.RAW where LLL is the line number and TTTT is the start time (UTC) of data collection in the format HHMM. In order to parse these files to extract navigation and depth information, two scripts were run under a Cygwin operating system. The first script "donav" is an executable shell script that cycles through all the *.RAW files in the folder and calls an AWK script to extract particular lines of information.
donav:
files=`ls *.RAW | cut -d. -f1`
for file in $files
do
	awk -f awkit $file.RAW > $file.navdep
done
The AWK script "awkit" simply extracts any line of information in the HYPACK file that contains either the string "GPGGA" or "SDDBT" and writes this information to a new file with the extension navdep.
awkit:
{
if ($0 ~ /GPGGA|SDDBT/) {
	print $0
	}
}
This process step and all subsequent process steps were performed by the same person - VeeAnn A. Cross.
Source_Used_Citation_Abbreviation: *.RAW
Process_Date: 201004
Source_Produced_Citation_Abbreviation: *.navdep
Process_Contact:
Contact_Information:
Contact_Person_Primary:
Contact_Person: VeeAnn A. Cross
Contact_Organization: U.S. Geological Survey
Contact_Position: Marine Geologist
Contact_Address:
Address_Type: mailing and physical address
Address: Woods Hole Coastal and Marine Science Center
Address: 384 Woods Hole Rd.
City: Woods Hole
State_or_Province: MA
Postal_Code: 02543-1598
Contact_Voice_Telephone: (508) 548-8700 x2251
Contact_Facsimile_Telephone: (508) 457-2310
Contact_Electronic_Mail_Address: vatnipp@usgs.gov
Process_Step:
Process_Description:
With the particular lines of interest extracted from the original HYPACK files, additional scripts were run to extract the specific information of interest. The next two scripts run are doholdhypack and awkholdhypack. The doholdhypack is a shell script run under Cygwin which cycles through all the *.navdep files in the folder and calls the awkholdhypack AWK script to process the file, with the results output to *.holdhypack. Of note, if the SDDBT record does not contain a depth value, a value of -9999 is written to the output file - acting as a nodata value.
doholdhypack:
files=`ls *.navdep | cut -d. -f1`
for file in $files
do
	awk -f awkholdhypack $file.navdep > $file.holdhypack
done
The AWK script awkholdhypack:
BEGIN {
FS = ","
}
{
FS = ","
depth = -9999
if ($1 ~ /GPGGA/)
	{
	utctime = $2
	latdeg = substr($3,1,2)
	latmin = substr($3,3,6)
	declat = latdeg + (latmin/60)
	londeg = substr($5,1,3)
	lonmin = substr($5,4,6)
	declon = -1 * (londeg + (lonmin/60))
	if (NR==1) {
		holddepth = -9999
		}
	else {
		printf("%s, %9.6f, %9.6f, %5.1f\n", holdutctime, holddeclon, holddeclat, holddepth)
	}
	holdutctime = utctime
	holdutcdate = utcdate
	holddeclon = declon
	holddeclat = declat
	holddepth = -9999
	}
if ($1 ~ /SDDBT/)
	{
	if ($4 != "")
		{
		depthreal = $4
		holddepth = depthreal
		}
	else
		{
		depthreal = -9999
		holddepth = -9999
		}
	}
}
END {
printf("%s, %9.6f, %9.6f, %5.1f\n", holdutctime, holddeclon, holddeclat, holddepth)
}
Source_Used_Citation_Abbreviation: *.navdep
Process_Date: 201004
Source_Produced_Citation_Abbreviation: *.holdhypack
Process_Step:
Process_Description:
With each individual HYPACK file processed, all the individual files on a given day were concatenated into a single file. Because the file naming convention starts with the line number, and surveying wasn't done in numerical order of lines, the appending of the files needs to be in the order the files were acquired. The time is in the filename, but another additional file written during acquisition was a LOG file containing the order of the lines. By editing this file and converting it to a shell script, the individual files can be concatenated in chronological order on each day. Running the shell script for each day under Cygwin successfully concatenated the files. What follows is the example from Julian day 103:
cat 201_1213.holdhypack \
202_1414.holdhypack \
203_1441.holdhypack \
204_1600.holdhypack > jd103hypack.csv
Source_Used_Citation_Abbreviation: *.holdhypack
Process_Date: 201004
Source_Produced_Citation_Abbreviation: jd103hypack.csv
Source_Produced_Citation_Abbreviation: jd104hypack.csv
Source_Produced_Citation_Abbreviation: jd105hypack.csv
Process_Step:
Process_Description:
Using VI editor, under Cygwin, the resulting CSV file for each day was edited to add the header line "gpstime, longitude, latitude, depth_m". With this header line added, the CSV file can be imported to ArcMap 9.2 as an event theme using Tools - Add XY data and defining the projection as Geographic, WGS84. This event theme is then converted to a point shapefile by right mouse click - data - export and saving a shapefile for each HYPACK event theme.
Source_Used_Citation_Abbreviation: jd103hypack.csv
Source_Used_Citation_Abbreviation: jd104hypack.csv
Source_Used_Citation_Abbreviation: jd105hypack.csv
Process_Date: 201004
Source_Produced_Citation_Abbreviation: jd103hypack.shp
Source_Produced_Citation_Abbreviation: jd104hypack.shp
Source_Produced_Citation_Abbreviation: jd105hypack.shp
Process_Step:
Process_Description:
The table of this shapefile was edited in ArcMap 9.2 to add two additional attributes: gpsdate and jday. These attributes were populated with the appropriate information for each day of data collection.
Source_Used_Citation_Abbreviation: jd103hypack.shp
Source_Used_Citation_Abbreviation: jd104hypack.shp
Source_Used_Citation_Abbreviation: jd105hypack.shp
Process_Date: 201004
Source_Produced_Citation_Abbreviation: jd103hypack.shp
Source_Produced_Citation_Abbreviation: jd104hypack.shp
Source_Produced_Citation_Abbreviation: jd105hypack.shp
Process_Step:
Process_Description:
Only the valid depth values are needed to create the bathymetry grid, so within ArcMap 9.2 under the shapefile properties, apply a definition query on the shapefile: "depth_m" <> -9999. Then use XTools Pro 5.2 to export the shapefile to a text file. Fields for export: gpstime, longitude, latitude, depth, jday. Save as type text, ANSI encoding.
Source_Used_Citation_Abbreviation: jd103hypack.shp
Source_Used_Citation_Abbreviation: jd104hypack.shp
Source_Used_Citation_Abbreviation: jd105hypack.shp
Process_Date: 201004
Source_Produced_Citation_Abbreviation: jd103_exp.txt
Source_Produced_Citation_Abbreviation: jd104_exp.txt
Source_Produced_Citation_Abbreviation: jd105_exp.txt
Process_Step:
Process_Description:
Use an AWK script to reformat each text file into the format needed for the MATLAB scripts that will be used for tide corrections. The AWK script awkbathy_fmt:
BEGIN {
FS=","
}
{
FS=","
hr=substr($1,1,2)
min=substr($1,3,2)
sec=substr($1,5,2)
jday = $5
longitude = $2
latitude = $3
depth = $4
printf("%s\t%s\t%s\t%s\t%9.6f\t%9.6f\t%3.1f\n", jday, hr, min, sec, longitude, latitude, depth)
}
Source_Used_Citation_Abbreviation: jd103_exp.txt
Source_Used_Citation_Abbreviation: jd104_exp.txt
Source_Used_Citation_Abbreviation: jd105_exp.txt
Process_Date: 201004
Source_Produced_Citation_Abbreviation: jd103_exp.dat
Source_Produced_Citation_Abbreviation: jd104_exp.dat
Source_Produced_Citation_Abbreviation: jd105_exp.dat
Process_Step:
Process_Description:
One day of equipment testing (JD102 - April 12, 2010) acquired resistivity data, but no HYPACK data. Additionally, in the western portion of the bay on JD105 (April 15) the HYPACK system was not recording, but navigation with bathymetry was being recorded by the CRP system during line L21F1 and L22F1. Although the test line on JD102 was not fully processed, the GPS file from that day was processed in the same manner as the other GPS files from the CRP data collection. An example of how the GPS files are handled can be read in the metadata for jd103gps_bestdepth.shp available from https://pubs.usgs.gov/of/2011/1039/html/ofr2011-1039-catalog.html. The processed shapefile from JD102, as well as selected records from JD105 (between GPS times 133421 and 140157) were exported using XTools Pro v 5.2 with the fields being exported - gpstime; longitude, latitude, depth_m and jday.
Source_Used_Citation_Abbreviation: jd102.shp
Source_Used_Citation_Abbreviation: jd105gps_spatjoin.shp
Process_Date: 201004
Source_Produced_Citation_Abbreviation: tstbathy.txt (JD102 resistivity portion)
Source_Produced_Citation_Abbreviation: jd105gps_spatjoin_extradepths_exp.txt
Process_Step:
Process_Description:
Each of the resistivity acquired bathymetry values were processed with the same AWK script (awkbathy_fmt) to generate files in the appropriate format for the MATLAB tide processing scripts.
awkbathy_fmt
BEGIN {
FS=","
}
{
FS=","
hr=substr($1,1,2)
min=substr($1,3,2)
sec=substr($1,5,2)
jday = $5
longitude = $2
latitude = $3
depth = $4
printf("%s\t%s\t%s\t%s\t%9.6f\t%9.6f\t%3.1f\n", jday, hr, min, sec, longitude, latitude, depth)
}
Source_Used_Citation_Abbreviation: tstbathy.txt (JD102 resistivity portion)
Source_Used_Citation_Abbreviation: jd105gps_spatjoin_extradepths_exp.txt
Process_Date: 201004
Source_Produced_Citation_Abbreviation: tstbathy.dat
Source_Produced_Citation_Abbreviation: bathy.dat (JD105 formatted bathymetry)
Process_Step:
Process_Description:
All of the HYPACK formatted bathymetry data and the resistivity data from JD102 were concatenated together in Cygwin using the command: cat tstbathy.dat jd103_hyp.txt jd104_hyp.txt jd105_hyp.txt > bathy.dat. The data from JD105 resistivity collection was handled separately since that gap wasn't noticed until later.
Source_Used_Citation_Abbreviation: tstbathy.dat
Source_Used_Citation_Abbreviation: jd103_hyp.txt
Source_Used_Citation_Abbreviation: jd104_hyp.txt
Source_Used_Citation_Abbreviation: jd105_hyp.txt
Process_Date: 201004
Source_Produced_Citation_Abbreviation: bathy.dat
Process_Step:
Process_Description:
In order to tide correct the bathymetry, data from tide stations in the area was required. There are two tide stations in the Indian River Bay area. One is from Rosedale Beach, DE, and the other is from the Indian River Inlet. The Rosedale Beach tides are available from: https://waterdata.usgs.gov/usa/nwis/uv?site_no=01484540 and the Indian River Inlet tides are available from: https://waterdata.usgs.gov/de/nwis/uv/?site_no=01484683&. In both cases the tab-delimited format was downloaded for all the days covering the CRP data acquisition. For both tide stations it's important to note that the time recorded in the tide files is in local time - Eastern Standard Time (EST), even though at the time of interest (April 12-15, 2010) local time is actually EDT (Eastern Daylight Time). The offset between EST and UTC is -5 hours (to go from EST to UTC). Other pertinent information from the website for the Rosedale Beach tide station: LOCATION.--Lat 38deg35min29.5sec, long. 75deg12min44.7sec, Sussex County, DE, Hydrologic Unit 02060010, on left bank attached to a privately owned fishing pier, at Seals Point, 1.9 miles west of Oak Orchard. DRAINAGE AREA.--Not determined. PERIOD OF RECORD.--April 1991 to current year. GAGE.--Water-stage recorder. Datum of gage is 0.0 ft above National Geodetic Vertical Datum of 1929. REMARKS.--This is a tidal station. Discharge are not determined at this location. U.S. Geological Survey satellite collection platform at station.
Other pertinent information from the website for the Indian River Inlet tide station: LOCATION.--Lat 38deg36min35.4sec, long. 75deg04min04.8sec, Sussex County, DE, Hydrologic Unit 02060010, 0.3 mi northwest of the Indian River Inlet, 0.2 mi west of State Highway 1, 4.9 mi north of Bethany Beach and at the Indian River Coast Guard station. DRAINAGE AREA.--Not determined. PERIOD OF RECORD.--June 1988 to June 1989, April 1991 to current year. GAGE.--Water-stage recorder. Datum of gage is 0.0 ft above National Geodetic Vertical Datum of 1929. Mean Low-Low Water: Gage Datum + 1.51 ft. REMARKS.--This is a tidal station. Discharge are not determined at this location. U.S. Geological Survey satellite collection platform at station.
Process_Date: 201004
Source_Produced_Citation_Abbreviation: rosedale_beach.txt
Source_Produced_Citation_Abbreviation: indianriver_inlet.txt
Process_Step:
Process_Description:
Using VI under Cygwin removed the header in the text file of the tide station data.
Source_Used_Citation_Abbreviation: rosedale_beach.txt
Source_Used_Citation_Abbreviation: indianriver_inlet.txt
Process_Date: 201004
Source_Produced_Citation_Abbreviation: rosedale_beach_noheader.txt
Source_Produced_Citation_Abbreviation: indianriver_inlet_noheader.txt
Process_Step:
Process_Description:
The tide data needs to be reformatted to work with the MATLAB processing scripts. In order to do this, process the files through an AWK script awktide_fmt:
{
year = substr($3,1,4)
month = substr($3,6,2)
day = substr($3,9,2)
hr = substr($4,1,2)
min = substr($4,4,2)
utchr = hr + 5
tideval = $6
tideval_meters = tideval * 0.3048
printf("1\t%s\t%s\t%02d\t%02d\t%02d\t%02d\t%02.4f\t%2.3f\n", $2, year, >month, day, utchr, min, tideval_meters, tideval)
}
The first element written to the output file is the number 1. The value of this number isn't important, but the MATLAB script does expect the first element to be a number. Another important piece of information about this AWK script is that 5 is added to the hours in the tide file to convert from EST to UTC to be compatible with the bathymetry data. Also, the original tide data was recorded in feet, but the work I'm doing is in meters. So that calculation is also taken care of in the AWK script.
Source_Used_Citation_Abbreviation: rosedale_beach_noheader.txt
Source_Used_Citation_Abbreviation: indianriver_inlet_noheader.txt
Process_Date: 201004
Source_Produced_Citation_Abbreviation: rd_tides.dat
Source_Produced_Citation_Abbreviation: ir_tides.dat
Process_Step:
Process_Description:
MATLAB version 7.5.0.342 (R2007b) was used to process the tide stations and bathymetry. The commands issued in MATLAB rely on MATLAB m files written by the USGS in Woods Hole, MA. The series of scripts used are: julian.m, linterp.m, gregorian.m, s2hms.m, and hms2h.m. There are actually two separate bathymetry files (both called bathy.dat) in separate folders that are being processed. One is the HYPACK data plus the CRP bathymetry data from Julian day 102. The other bathy.dat file represents the bathymetry extracted from Julian day 105 when the HYPACK system was not recording. So this series of commands is actually run twice - in the two separate folders. The following MATLAB commands were issued to process the data including comments explaining some of the processes. The comment lines are prefaced with the "%".
%Step 1. Outside of Matlab, format the ASCII tide gauge data.
%Load the tide data.
load ir_tides.dat
load rd_tides.dat
greg_ir=ir_tides(:,3:7);  %year month day hour minute needs to be in columns 3,4,5,6,7
greg_ir(:,6)=zeros(length(greg_ir(:,1)),1);  %adding secs bc the tide file doesn't have any
greg_rd=rd_tides(:,3:7);
greg_rd(:,6)=zeros(length(greg_rd(:,1)),1);
%then establish the julian timebase for the observed tide data
jd_ir=julian(greg_ir);
jd_rd=julian(greg_rd);
h_ir=ir_tides(:,8);  %observed tide in meters (in this case above NGVD29)
h_rd=rd_tides(:,8);
%Step 2: outside of Matlab I reformatted my bathy data.  The original data was exported
%from a shapefile using XTools and then needed to be reformatted for my needs here
load bathy.dat %there are 2 bathy.dat file in 2 separate folders
jd0=julian([2010 1 1 0 0 0]);  %starting point for yearday
%setting the julian time base for the ship stuff
%this assumes column 1=jd, 2=hour, 3=min, 4=sec
jdtrack=jd0+bathy(:,1)-1+(bathy(:,2)+bathy(:,3)/60+bathy(:,4)/3600)/24;
x=bathy(:,5);  %longitude in decimal degrees in column 5
y=bathy(:,6);  %latitude in decimal degrees in column 6
htrack=bathy(:,7);
%Step 3: interpolate observed tide data onto the ship track timebase using a linear
%interpolation. Splining the tide data could result in overshoots
tidetrack_ir=linterp(jd_ir,h_ir,jdtrack);
tidetrack_rd=linterp(jd_rd,h_rd,jdtrack);
h_corrected_ir=htrack-tidetrack_ir;
h_corrected_rd=htrack-tidetrack_rd;
%Step 4: output the information I want into a text file.  I decided I basically want my
%original bathymetry data, plus the tide value and the corrected value
%I have to use the ' to get the columns flipped so the append thing works
alltidestuff=[bathy'; tidetrack_ir'; h_corrected_ir';tidetrack_rd';h_corrected_rd'];
outfile=['veetest'];
fid=fopen(outfile,'w');
fprintf(fid,'%03d:%02d:%02d:%02d, %9.6f, %9.6f, %3.1f, %3.1f, %3.1f, %3.1f, %3.1f \n',alltidestuff);
fclose(fid);
The result is a comma-delimited text file with 8 columns of information: jdtime, longitude, latitude, depth_m, ir_tide, ir_cor, rd_tide, rd_cor The definition of this header information is as follows:
jdtime: julian day and time from the GPS
longitude: longitude from the GPS
latitude: latitude from the GPS
depth_m: depth in meters of the fathometer (as recorded by the GPS)
ir_tide: interpolated tide values at the GPS time from the Indian River Inlet tide station
ir_cor: the corrected depth in meters of the fathometer based on the Indian River tide value (depth_m - ir_tide)
rd_tide: interpolated tide values at the GPS time from the Rosedale Beach tide station
rd_cor: the corrected depth in meters of the fathometer based on the Rosedale tide value (depth_m - rd_tide)
The header line was added to the veetest file and saved under a new file name: hypbathytides.csv and jd105_moredepths.csv
Source_Used_Citation_Abbreviation: bathy.dat
Process_Date: 201004
Source_Produced_Citation_Abbreviation: hypbathytides.csv
Source_Produced_Citation_Abbreviation: jd105_moredepths.csv
Process_Step:
Process_Description:
Each of the resulting CSV files were added to ArcMap 9.2 as event themes: Tools - Add XY Data. The projection was defined as Geographic, WGS84. These event themes were then converted to shapefiles - right mouse click on the event theme - Data Export.
Source_Used_Citation_Abbreviation: hypbathytides.csv
Source_Used_Citation_Abbreviation: jd105_moredepths.csv
Process_Date: 201004
Source_Produced_Citation_Abbreviation: hypbathytides.shp
Source_Produced_Citation_Abbreviation: jd105_moredepths.shp
Process_Step:
Process_Description:
The resulting shapefiles are geographic, in order to do a weighted tide influence on each point, the point shapefiles are projected to UTM, Zone 18 so that distance values are in meters. Using ArcMap 9.2 - ArcToolbox - Data Management Tools - Projections and Transformations - Feature. Input each geographic shapefile, and output a new shapefile in a UTM, Zone 18, WGS84 projection - no transformation necessary.
Source_Used_Citation_Abbreviation: hypbathytides.shp
Source_Used_Citation_Abbreviation: jd105_moredepths.shp
Process_Date: 201004
Source_Produced_Citation_Abbreviation: hypbathytides_utm18.shp
Source_Produced_Citation_Abbreviation: jd105_moredepths_utm18.shp
Process_Step:
Process_Description:
The tide station locations are known, so that latitude/longitude information was added to a comma-delimited text file and added to ArcMap 9.2 as an event theme and converted to a shapefile. Within ArcMap 9.2, the Geographic, WGS84 shapefile was projected to a UTM, Zone 18, WGS84 shapefile.
Source_Used_Citation_Abbreviation: tidesta_locs.csv
Process_Date: 201004
Source_Produced_Citation_Abbreviation: tidestations.shp
Source_Produced_Citation_Abbreviation: tidestations_utm18.shp
Process_Step:
Process_Description:
The UTM tide station shapefile and the two shapefiles with the bathymetry and tide information were added to ArcMap 9.2. One tide station (Indian River Inlet) is at the eastern end of the study area, while the other tide station (Rosedale Beach) is at the western end of the study area. Instead of just using the adjustments from one tide station, a weighted tide adjustment was applied based on the distance from the tide stations. In order to get the distance to each tide station in the bathymetry shapefiles, the Near tool will be used (ArcToolbox - Analysis Tools - Proximity - Near). Prior to using the tool, add the following attributes to the shapefile (as type double): dist_inlet; dist_rose; wt_tide; wttide_dep. With these attributes added, apply a definition query to the tide station shapefile so only one station is present. The Near tool will only report the distance to the closest tide station if two are available. To get the distance to each tide station so that a weighted tide value can be calculated based on distance to each tide station, the Near tool will need to be run twice - with only one station available each time. So with one tide station visible, use the near tool with the following parameters: input-hypbathytides_utm18; near features - tidestations_utm18. leave everything else at default (don't use a search radius). Then copy the "near_dist" attribute (automatically created by the Near tool) to the appropriate distance attribute created for that station. Repeat the same process for the other tide station. Repeat the same process with the jd105_moredepths_utm18.shp (add the attributes and run the Near tool with the tide stations).
Source_Used_Citation_Abbreviation: hypbathytides_utm18.shp
Source_Used_Citation_Abbreviation: jd105_moredepths_utm18.shp
Source_Used_Citation_Abbreviation: tidestations_utm18.shp
Process_Date: 201004
Source_Produced_Citation_Abbreviation: hypbathytides_utm18.shp
Source_Produced_Citation_Abbreviation: jd105_moredepths_utm18.shp
Process_Step:
Process_Description:
In ArcMap 9.2, after using the Near tool and copying the distances to the tide stations into the appropriate fields, use the field calculator on the "wt_tide" attribute.

wt_tide = ((1-([dist_rose]/( [dist_rose] + [dist_inlet])))* [rd_tide]) + ((1-( [dist_inlet]/( [dist_rose] + [dist_inlet])))* [ir_tide])

This should give a weighted tide value for each point based on the distance to each tide station. Now apply the tide correction: tideadded = [depth_m] - [wt_tide]
Source_Used_Citation_Abbreviation: hypbathytides_utm18.shp
Source_Used_Citation_Abbreviation: jd105_moredepths_utm18.shp
Process_Date: 201004
Source_Produced_Citation_Abbreviation: hypbathytides_utm18.shp
Source_Produced_Citation_Abbreviation: jd105_moredepths_utm18.shp
Process_Step:
Process_Description:
Use topo2raster to generate a new grid. ArcMap 9.2 - ArcToolbox - Spatial Analyst Tools - Interpolation - Topo to Raster.
Inputs:
bounding_poly; ; boundary
jd105_moredepths_utm18; wttide_dep; PointElevation
hypbathytides_utm18; wttide_dep; PointElevation
output: bestbathy
cell size: 25
drainage enforcement: no_enforce
primary type of input data: spot
Everything else left to defaults.

The bounding polygon was a quick polygon created around the survey area to limit the extent of extrapolation in areas where no data were available.
Source_Used_Citation_Abbreviation: jd105_moredepths_utm18
Source_Used_Citation_Abbreviation: hypbathytides_utm18
Process_Date: 201004
Source_Produced_Citation_Abbreviation: bestbathy
Process_Step:
Process_Description:
Then filter the resulting grid: ArcToolbox - Spatial Analyst Tools - Neighborhood - Filter
Input: bestbathy
Output: bestbat_fil
Filter type: low
Check the box next to ignore NODATA in calculations.
Source_Used_Citation_Abbreviation: bestbathy
Process_Date: 201004
Source_Produced_Citation_Abbreviation: bestbat_fil
Process_Step:
Process_Description:
Used Spatial Analyst toolbar in ArcMap 9.2 - Raster Calculator and multiplied the grid bestbat_fil by -1. Then right mouse clicked on the calculation grid - Data - Make permanent with the output: bbat_filneg
Source_Used_Citation_Abbreviation: bestbat_fil
Process_Date: 201009
Source_Produced_Citation_Abbreviation: bbat_filneg
Process_Step:
Process_Description: The bathymetry grid was renamed to irb_bathy in ArcCatalog 9.2.
Source_Used_Citation_Abbreviation: bbat_filneg
Process_Date: 201012
Source_Produced_Citation_Abbreviation: irb_bathy
Process_Step:
Process_Description:
Edits to the metadata were made to fix any errors that MP v 2.9.36 flagged. This is necessary to enable the metadata to be successfully harvested for various data catalogs. In some cases, this meant adding text "Information unavailable" or "Information unavailable from original metadata" for those required fields that were left blank. Other minor edits were probably performed (title, publisher, publication place, etc.). Attempted to modify http to https where appropriate. Moved the minimal source information provided to make it the first process step. The metadata date (but not the metadata creator) was edited to reflect the date of these changes. The metadata available from a harvester may supersede metadata bundled within a download file. Compare the metadata dates to determine which metadata file is most recent.
Process_Date: 20171002
Process_Contact:
Contact_Information:
Contact_Organization_Primary:
Contact_Organization: U.S. Geological Survey
Contact_Person: VeeAnn A. Cross
Contact_Position: Marine Geologist
Contact_Address:
Address_Type: Mailing and Physical
Address: 384 Woods Hole Road
City: Woods Hole
State_or_Province: MA
Postal_Code: 02543-1598
Contact_Voice_Telephone: 508-548-8700 x2251
Contact_Facsimile_Telephone: 508-457-2310
Contact_Electronic_Mail_Address: vatnipp@usgs.gov
Process_Step:
Process_Description: USGS Thesaurus keywords added to the keyword section.
Process_Date: 20180720
Process_Contact:
Contact_Information:
Contact_Organization_Primary:
Contact_Organization: U.S. Geological Survey
Contact_Person: VeeAnn A. Cross
Contact_Position: Marine Geologist
Contact_Address:
Address_Type: Mailing and Physical
Address: 384 Woods Hole Road
City: Woods Hole
State_or_Province: MA
Postal_Code: 02543-1598
Contact_Voice_Telephone: 508-548-8700 x2251
Contact_Facsimile_Telephone: 508-457-2310
Contact_Electronic_Mail_Address: vatnipp@usgs.gov
Process_Step:
Process_Description:
Added keywords section with USGS persistent identifier as theme keyword.
Process_Date: 20200908
Process_Contact:
Contact_Information:
Contact_Organization_Primary:
Contact_Organization: U.S. Geological Survey
Contact_Person: VeeAnn A. Cross
Contact_Position: Marine Geologist
Contact_Address:
Address_Type: Mailing and Physical
Address: 384 Woods Hole Road
City: Woods Hole
State_or_Province: MA
Postal_Code: 02543-1598
Contact_Voice_Telephone: 508-548-8700 x2251
Contact_Facsimile_Telephone: 508-457-2310
Contact_Electronic_Mail_Address: vatnipp@usgs.gov
Spatial_Data_Organization_Information:
Direct_Spatial_Reference_Method: Raster
Raster_Object_Information:
Raster_Object_Type: Grid Cell
Row_Count: 246
Column_Count: 510
Vertical_Count: 1
Spatial_Reference_Information:
Horizontal_Coordinate_System_Definition:
Planar:
Grid_Coordinate_System:
Grid_Coordinate_System_Name: Universal Transverse Mercator
Universal_Transverse_Mercator:
UTM_Zone_Number: 18
Transverse_Mercator:
Scale_Factor_at_Central_Meridian: 0.999600
Longitude_of_Central_Meridian: -75.000000
Latitude_of_Projection_Origin: 0.000000
False_Easting: 500000.000000
False_Northing: 0.000000
Planar_Coordinate_Information:
Planar_Coordinate_Encoding_Method: row and column
Coordinate_Representation:
Abscissa_Resolution: 25.000000
Ordinate_Resolution: 25.000000
Planar_Distance_Units: meters
Geodetic_Model:
Horizontal_Datum_Name: D_WGS_1984
Ellipsoid_Name: WGS_1984
Semi-major_Axis: 6378137.000000
Denominator_of_Flattening_Ratio: 298.257224
Vertical_Coordinate_System_Definition:
Altitude_System_Definition:
Altitude_Datum_Name: National Geodetic Vertical Datum of 1929
Altitude_Resolution: 0.1
Altitude_Distance_Units: meters
Altitude_Encoding_Method: Implicit coordinate
Entity_and_Attribute_Information:
Overview_Description:
Entity_and_Attribute_Overview:
This grid is in meters in the vertical datum of NGVD29. Based on the NOAA website (http://www.ngs.noaa.gov/cgi-bin/VERTCON/vert_con.prl accessed October 2010) the conversion from NGVD29 to NAVD88 is -0.238 meters at both the Indian River Inlet tide station and the Rosedale Beach tide station. This indicates a direct offset of the whole grid will convert to NAVD88. So taking this grid and subtracting 0.238 meters will yield a grid with the vertical datum NAVD88.
Entity_and_Attribute_Detail_Citation: U.S. Geological Survey
Distribution_Information:
Distributor:
Contact_Information:
Contact_Person_Primary:
Contact_Person: VeeAnn A. Cross
Contact_Organization: U.S. Geological Survey
Contact_Position: Marine Geologist
Contact_Address:
Address_Type: mailing and physical address
Address: Woods Hole Coastal and Marine Science Center
Address: 384 Woods Hole Rd.
City: Woods Hole
State_or_Province: MA
Postal_Code: 02543-1598
Contact_Voice_Telephone: (508) 548-8700 x2251
Contact_Facsimile_Telephone: (508) 457-2310
Contact_Electronic_Mail_Address: vatnipp@usgs.gov
Resource_Description: Downloadable Data
Distribution_Liability:
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.
Standard_Order_Process:
Digital_Form:
Digital_Transfer_Information:
Format_Name: AIG
Format_Version_Number: 9.2
Format_Information_Content:
The WinZip (v. 14.0) file contains the Esri binary grid - data contained in the grid folder and the info folder, as well as the associated metadata files.
Transfer_Size: 0.25
Digital_Transfer_Option:
Online_Option:
Computer_Contact_Information:
Fees: None.
Technical_Prerequisites:
This WinZip file contains data available in Esri binary grid format. The user must have software capable of uncompressing the WinZip file and reading/displaying the grid.
Metadata_Reference_Information:
Metadata_Date: 20200908
Metadata_Contact:
Contact_Information:
Contact_Person_Primary:
Contact_Person: VeeAnn A. Cross
Contact_Organization: U.S. Geological Survey
Contact_Position: Marine Geologist
Contact_Address:
Address_Type: mailing and physical address
Address: Woods Hole Coastal and Marine Science Center
Address: 384 Woods Hole Rd.
City: Woods Hole
State_or_Province: MA
Postal_Code: 02543-1598
Contact_Voice_Telephone: (508) 548-8700 x2251
Contact_Facsimile_Telephone: (508) 457-2310
Contact_Electronic_Mail_Address: vatnipp@usgs.gov
Metadata_Standard_Name: FGDC Content Standards for Digital Geospatial Metadata
Metadata_Standard_Version: FGDC-STD-001-1998
Metadata_Time_Convention: local time

This page is <https://cmgds.marine.usgs.gov/catalog/whcmsc/open_file_report/ofr2011-1039/irb_bathymeta.html>
Generated by mp version 2.9.50 on Tue Sep 21 18:20:35 2021