This portion of the USGS data release presents bathymetry data collected during surveys performed in the Columbia River littoral cell, Washington and Oregon in 2016 (USGS Field Activity Number 2016-663-FA). Bathymetry data were collected using four personal watercraft (PWCs) equipped with single-beam sonar systems and global navigation satellite system (GNSS) receivers. The sonar systems consisted of an Odom Echotrac CV-100 single-beam echosounder and 200 kHz transducer with a 9 degree beam angle. Raw acoustic backscatter returns were digitized by the echosounder with a vertical resolution of 1.25 cm. Depths from the echosounders were computed using sound velocity profiles measured using a YSI CastAway CTD during the survey. Positioning of the survey vessels was determined at 5 to 10 Hz using Trimble R7 GNSS receivers. Output from the GNSS receivers and sonar systems were combined in real time on the PWC by a computer running HYPACK hydrographic survey software. Navigation information was displayed on a video monitor, allowing PWC operators to navigate along survey lines at speeds of 2 to 3 m/s. Survey-grade positions of the PWCs were achieved with a single-base station and differential post-processing. Positioning data from the GNSS receivers were post-processed using Waypoint Grafnav to apply differential corrections from a GNSS base station with known horizontal and vertical coordinates relative to the North American Datum of 1983. Orthometric elevations relative to the NAVD88 vertical datum were computed using National Geodetic Survey Geoid12a offsets. Bathymetric data were merged with post-processed positioning data and spurious soundings were removed using a custom Graphical User Interface (GUI) programmed with the computer program MATLAB. The average estimated vertical uncertainty of the bathymetric measurements is 10 cm. The final point data from the PWCs are provided in a comma-separated text file and are projected in cartesian coordinates using the Washington State Plane South, meters coordinate system.