Open File Report
RV INLAND SURVEYER CRUISE IS-98
Lake Tahoe, California and Nevada
This report is preliminary and has not been reviewed for conformity with U.S. Geological Survey editorial standards or with the North American Stratigraphic Code. Use of trade, product, or firm names in this report is for descriptive purposes only and does not imply endorsement by the U.S. Government.
1 U.S. Geological Survey, MS 999 345 Middlefield Road, Menlo
Park, CA 94025
THE BATHYMETRY OF LAKE TAHOE, CALIFORNIA-NEVADA
The major objective of cruise IS-98 was to map the bathymetry of Lake Tahoe, California-Nevada (Fig. 1) to fulfill a commitment made during the Lake Tahoe Presidential Forum in 1997. The only existing bathymetry of Lake Tahoe, collected in 1923, was recently compiled by Rowe and Stone (1997), but the data density is inadequate for the level of scientific studies ongoing and anticipated in the near future for Lake Tahoe. Recent advances in marine multibeam-sonar capabilities now permit a cost-effective way, to precisely map the bathymetry of large areas of the ocean floor with 100% coverage. Cruise IS-98 applied this state-of-the-art ocean technology to Lake Tahoe. The newest of these high-resolution multibeam mapping systems also simultaneously collects backscatter (similar to sidescan sonar) imagery that results in a complimentary and co-registered data set that is related to the distribution of lake-floor materials and textures. The two types of maps that resulted from this cruise provide the multidisiplinary Lake Tahoe research community an unprecedented set of base maps upon which to build their studies. This report describes the high-resolution multibeam mapping system used at Lake Tahoe, outlines the data-processing steps used to produce the maps, and includes the daily log of the cruise.
There are several different brands of high-resolution multibeam mapping systems that are appropriate for shallow-water surveys. After a review of the currently available systems, we chose to use the Kongsberg Simrad EM1000 system because of its proven track record, its demonstrated portability, its ability to map large areas at high speed without compromising data quality and, most importantly, its ability to simultaneously produce high-resolution side-scan sonar imagery. For the Lake Tahoe survey, we used an EM1000 system owned and operated by C&C Technologies, Inc. (Lafayette, LA) installed aboard the 26-ft Inland Surveyor (Fig. 2).
An overview of high-resolution multibeam mapping systems can be found in Hughes-Clarke, et al. (1996). The Simrad EM1000 system operates at a frequency of 95 kHz from a semi-circular transducer that is mounted on a rigidly attached ram on the bow of Inland Surveyor. The system was designed to operate in several modes through a range of depths from 5 to approximately 800-m water depth. The shallow (ultrawide) mode, used to maximum depths of about 200 m, forms 60 receive beams with a spacing of 2.5° distributed across track and 3° wide along track. This geometry generates a 150° swath that covers 7.4 times the water depth (Fig. 3). The wide mode, used to maximum depths of about 600 m, uses 48 formed beams with a spacing of 2.5 ° that results in a swath of 120° or 3.4 times the water depth. The deep mode, for water depths greater than 500 m, uses 48 formed beams with a spacing of 1.25° that generates a swath of 60° or 1.1 times the water depth. There are options within these modes for beam distribution (equiangular or equidistant) and pulse lengths. The specific options used for the Lake Tahoe survey are discussed in the data processing section below.
Most conventional vertical-incidence echo sounders determine the time of arrival of the returned pulse (and thus the depth) by detecting the position of the sharp leading edge of the returned echo, a technique called amplitude detection. In multibeam sonars, where the angle of incidence increases for each consecutive receive beam to either side of the vertical (nadir), a returned echo loses its sharp leading edge and the accurate determination of depth becomes inaccurate. To address this problem, the EM1000 multibeam system uses an interferometric principle in which each beam is split, through electronic beamforming, into two "halfbeams" and the phase difference between these "halfbeams" is calculated to provide a measure of the angle of arrival of the echo. The point at which the phase is zero (i.e., where the wavefront of the returned echo is normal to the receive-beam bore) is determined for each beam and provides an accurate measure of the range to the seafloor. Both amplitude and phase detection are recorded for each beam and then the system software picks the "best" detection method for each beam, based on a number of quality measurements, and uses this method to calculate depth.
The EM1000 also provides quantitative seafloor-backscatter data that can be displayed in a sidescan-sonar-like image. The backscatter images can be used to gain insight into the spatial distribution of seafloor properties. A time series of echo amplitudes from each beam is recorded at 0.2- to 2.0-ms sampling rate, depending on the water depth. The echo amplitudes are sampled a much faster rate than the beam spacing and can be processed from beam-to-beam to produce a backscatter image with the theoretical resolution of the sampling interval (15 cm at 0.2 ms). The amplitude information can be placed in its geometrically correct position relative to the across-track profile because the angular direction of each range sample is known. The EM1000 software corrects the amplitude time series for gain changes, propagation losses, predicted beam patterns and for the insonified area (with the simplifying assumptions of a flat seafloor and Lambertian scattering). Subsequent processing (see below) uses real seafloor slopes and applies empirically derived beam-pattern corrections to produce a quantitative estimate of seafloor backscatter across the swath.
In addition to the multibeam sonar array, a multibeam survey requires a careful integration of a number of ancillary systems. These include: (1) a differentially corrected Global Positioning System (DGPS); (2) an accurate measure the heave, pitch, roll, and heading of the vessel, all to better than 0.01° and the transformation these measurements to estimates of the motion of the transducer at the times of transmission and reception (motion sensor); (3) a method to precisely determine the sound-speed structure of the water column, using measurements of temperature, salinity and depth with a CTD and the calculation of sound velocity profiles [SVP].
The Lake Tahoe survey was navigated with inertial positioning provided by an Applied Analytic POS/MV motion sensor as well as a Trimble DGPS. The POS/MV records vehicle motion at 100 Hz with an accuracy of 0.01°. USGS benchmarks were located at two sites on the north end of Lake Tahoe and their positions were determined with a DGPS. These provided a local calibration for a position accuracy for the survey of ±1 m.
Sound-velocity profiles were calculated several times each day so that ray-tracing techniques could be used to determine the effect of acoustic refraction in the water, principally caused by variations in water temperature. A SeaBird CTD was deployed the first day of operations to get a good reference SVP. Sippican T-7 expendable bathythermographs (0 to 600-m water depth) were routinely collected anytime it was suspected that the water structure had changed during the course of the survey. An additional sound-velocity profiler is installed at the transducer array to determine the speed of sound in water directly at the transducer. All the SVP data are fed directly into the Simrad EM1000 processor for instantaneous beamforming and raytracing of the individual beams.
Raw EM1000 data telegrams were acquired over an ethernet on a shipboard workstation. The data stream include is shown in Table 1. A number of ancillary data sources were also acquired by C&C Technologies (Table 2).
The accurate reduction of swath bathymetric data critically depends on a proper knowledge of the geometry and relative positions of the sonar transducer with respect to the motion sensor, the ship, and the positioning-system antennae (Fig. 4). C&C Technologies, using standard surveying techniques, measured these values before the survey began (Table 3). The transducer is at 1.3-m depth (vessel stationary). All values are relative to the transducer.
Three Vertical Reference Units were available onboard the Inland Surveyor to measure the pitch, roll, and yaw. The POS/MV, a full DGPS-aided inertial navigational system, was used as the primary attitude sensor and directly interfaced to the Simrad EM1000 and provided 100-Hz measurements of boat attitude to 0.01°. The POS/MV was chosen because it has high accuracy specifications, it is insensitive to long period horizontal accelerations, and it provides an inertial position solution that is reliable through DGPS outages of periods of less than about 1 or two minutes.
There are several operational modes available for the EM1000. The differences in the modes are a function of pulse length, beam spacing, and angular sector. The pulse length controls the amount of energy transmitted into the water column. The system can be operated in an "equiangular" (EA) mode in which the beams are spaced at equal angles apart, resulting in a non-linear (increasing spacing away from nadir) spacing of sonar footprints on the seafloor. The system can also be operated in an "equidistant" (EDBS) mode in which the beams are spaced such that the sonar footprints are equally spaced in the across-track profile. The EDBS geometry is achieved by generating variable beam-angular spacings. Although EDBS has advantages in data handling (i.e., provides even sounding density), there are two limitations. The beams in the 140° and 150° modes are spaced wider than their beam widths and results in incomplete coverage that produces a striping close to nadir. This problem disappears as the swath width closes to ~120°. However, the second limitation occurs because of attitude uncertainties and imperfect refraction models that can result in sounding errors that grow with angle from the vertical. Consequently, these limitations render the outermost beams less reliable than for the EA mode.
For the reasons stated above, the EA mode was used for the Lake Tahoe survey. In the EA mode, the EM1000 was operated with a 0.2 ms pulse length in waters less than 150 m deep, that provided a 150° swath. In waters deeper than 150 m, the EM1000 was operated with a 0.7 ms pulse length that generated a 120° swath.
All bathymetric data were adjusted through Kongsberg Simrad software for (1) transducer draft, (2) static roll, pitch and gyro misalignments, (3) pitch on transmit, (4) heave at transmit and reception, (5) roll at reception, (6) refracted ray path, and (7) beam steering at transducer interface. Post-logging transformations included (1) transformation of navigation from antenna to transducer, (2) correction for positioning to sonar time shifts, (3) lake surface elevation, and (4) any unaccounted for static attitude misalignments.
The Kongsberg Simrad EM1000 provides a backscatter-intensity time series for the bottom insonification period for each of the 60 individual beams. The corrections applied by the shipboard recording system are listed in the table below in Table 4.
A set of required transformations are performed by specialized software written by the Ocean Mapping Group at the University of New Brunswick. These transformations include conversion of each beam time series to a horizontal range equivalent, splicing the 60 beam traces together to produce one full slant-range corrected trace, and removal of residual beam-pattern effects. Although the system software corrects for average beam pattern, there are ± 2 dB ripples in the average beam pattern that vary from transducer to transducer. The beam pattern of the C & C EM1000 may be slightly different each time the stave connectors are reconnected in the transceiver.
Our processing approach to backscatter was to stack several thousand pings to view the angular variation of received backscatter intensity as a function of beam angle. Inherent in this function is both the transmit and receive sensitivities, as well as the mean angular response of the seafloor. We then invert this function to minimize the beam pattern and angular variations.
Kongsberg Simrad uses a variable gain within 15° of vertical to reduce logged dynamic range. The sidescan data at this stage have had a Lambertian response backed out and the beam pattern has been corrected with respect to the vertical and all receive beams have been roll stabilized. Consequently, corrections have been made for variations in the beam-forming amplifiers but not variations in the stave sensitivities of the physical array.
Additional transformations were required to produce calibrated backscatter measurements. These include (1) removal of Lambertian model, (2) true lake-floor slope correction, (3) refracted ray path correction, (4) residual beam pattern correction, and (5) aspherical-spreading corrections.
Despite the careful measurements of transducer alignments and offsets, the true geometry of the installed system can only be determined through the determination of the self-consistency of seafloor measurements. To facilitate such a determination, we conducted a series of "patch tests" in Lake Tahoe whereby the system was run back and forth across various features on the lake floor to determine if there were residual roll, pitch, heading, or timing offsets that required corection factors.
A full patch test procedure was started prior to data collection on August 2 to calibrate any time delay, and gyro misalignment. The static adjustments were estimated from the patch test are listed in Table 5.
The 1-Hz DGPS data were logged with the Kongsberg Simrad EM1000 software. The Simrad Bottom Detection Unit (BDU time stamps the depth and sidescan telegrams and was slaved to a shipboard SUN Sparc 20 that itself was synched to the GPS 1 PPS. The navigation telegrams were externally stamped by the Trimble 4000 GPS receiver. The receiver antenna positions were shifted to the transducer according to the X and Y offsets using the POS/MV output.
Every 1-Hz navigation fix was checked for gross time and/or distance jumps by graphical examination during data processing. Outliers were interactively interrogated for time, flagged and rejected (or re-accepted). All navigation jumps greater than 20 s were automatically flagged as uninterpolatable.
Land-based data processing (Fig. 5) consisted of (1) the editing the 1-Hz navigation fixes to flag bad fixes; (2) examining each ping of each beam to flagg outliers, bad data, etc.; (3) merging the depth and backscatter data with the cleaned navigation; (4) reducing all depth values to mean lake level (1899 m); (5) performing additional refraction corrections, if necessary, for correct beam raytracing; (6) separating out the amplitude measurements for conversion to backscatter; (7) gridding depth and backscatter into a geographic projection at the highest resolution possible with water depth; (8) regridding individual subareas of bathymetry and backscatter into final georeferenced map sheets; (9) gridding and contouring the bathymetry; and (10) generation of the final maps. Nearly finalized maps were completed in the field one day after the completion of the cruise and the final maps that accompany this report were completed one week after the end of the cruise.
All swath bathymetric data were passed through a number of simple filters that flagged all beams that form a spike when compared to their neighbors for which the across-track slope anomalies exceed 25° and flagged all beams that form a spike when compared to their neighbors for which the along-track slope anomalies exceed 55°. These two criteria removed the worst of the spikes. All beams without immediate neighbors fore and aft were rejected and all beams without immediate neighbors port and starboard were rejected. The two criteria removed isolated beams but the remaining outliers had to be examined and subjectively examined using the UNB/OMG SwathEditor tool in a line by line mode.
The single biggest limitation on the quality of sounding data is water-column refraction. Refraction-related anomalies grow non-linearly with beam angle and the resulting artifacts can create short-wavelength topographic features that may be misinterpreted as seabed geology. There was some fear voiced that the suspected strong lake stratification would present a problem for the beam steering and ray tracing of individual beams. This did not prove to be the case and the thermocline and density structure, as repeatedly measured by CTD and XBT, remaimed remarkable constant over the 16 days of the survey. A representative plot of a water-velocity profile is shown in Figure 6. In fact, no additional emperical refraction corrections were necessary during processing. If all of the alignments were correctly determined, Kongsberg Simrad states that depth "...accuracies better than 0.5% of water depth are achieved across the entire swath."
A lake-elevation station, located on the Coast Guard pier in Tahoe City, is maintained by the USGS. Readings are taken every four hours so that a detailed time series was available for the elevation of the lake surface relative to mean sea level. The lake elevation on August 2nd was 1899.08 m and on August 17th it was 1899.02 m. A lake elevation of 1899 m was used as the reference datum for all depth determinations. However, it appears that there may be a 0.43-m offset between the datum used for the elevation and the NGVD datum (Gerald Rockwell, personal commun., 1998). The lake level dropped only 6 cm over the duration of the mapping survey, but because the drop and the potential datum offset are smaller than the resolution of the system, so we did not adjust the depth determinations for these offsets. Once each lake depth determination was calculated, the value was converted to meters above mean sea level. This allows the entire Tahoe bathymetry to be merged with the regional DEM and to be viewed as a topographic feature.
The large overview maps of backscatter and shaded relief that accompany this report were generated from larger-scale subarea maps (Fig. 7 and Table 6). The higher-resolution subarea maps were regridded at 10 m to produce the maps of the entire lake (Figs. 7 and 8). This regridding loses resolution in the shallower areas but allows the entire area to be mapped. The detailed subarea maps were produced at the maximum resolution allowable for the data.
Contour maps represent the more traditional method of displaying bathymetry. The contours were derived from the gridded elevations. The resultant contours were smoothed by a 3-point running average in the overview maps. Even at the original contour grid, more than 90% of the data have to be discarded so as to only show some chosen contour interval. A much better representation of bathymetry, using 100% of the data is a shaded-relief map.
A shaded-relief map (Fig. 8) is a pseudo-sun-illumination of a topographic surface using the Lambertian scattering law (equation 1),
SI = K * cos phi (Eq. 1)
where SI is the pseudo-sun intensity, K is a constant that allows for even background, and phi is the angle between the pseudo sun and the bathymetric surface.
The backscatter map (Fig. 9) is a representation of the amount of acoustic energy, at 95 kHz, that is scattered back to the receiver from the lake floor. The Kongsberg Simrad EM1000 system has been calibrated at the factory and all gains, power levels, etc. that are applied during signal generation and detection are recorded for each beam and are used to adjust the amplitude value prior to recording. Consequently, the backscatter is calibrated to an absolute reflectance of the seabed. However, the amount of energy, measured in decibels (dB), is some complex function of constructional and destructional interference caused by the interaction of an acoustic wave with a volume of sediment or, in the case of hard rock, the rock material. The backscatter from sediment is volume reverberation to at least 10 cm caused by seabed and subsurface interface roughnesses above the Rayleigh criteria (a function of acoustic wave length), the composition of the sediment, and its bulk properties (water content, bulk density, etc.). Although, it is not possible to determine a unique geological facies from the backscatter value, reasonable predictions can be made from the backscatter based on the local geology.
It can not be stressed too strongly that one of the great advantages of this survey is that the bathymetry is completely georeferenced with the backscatter. Consequently, each pixel on the map has a latitude, longitude, depth, and backscatter value assigned to it.
August 3, 1998 (980803)
A second surprise was that the 1900-m elevation of Lake Tahoe rendered the engines of Inland Surveyor less efficient than was predicted. The result of this inefficiency was that a maximum speed of only ~8 kts could be maintained. A new set of propellers was ordered to increase speed.
C&C located two USGS benchmarks around the lake and calibrated their navigation systems with ease. The Coast Guard measures the lake elevation (to the nearest millimeter) once a day and agreed to provide me with those measurements. These data provided our "tides" for producing corrected elevations for the lake floor.
August 4, 1998 (980804)
August 5, 1998 (980805)
The sediments in the chutes between headlands and in small channels have such low backscatter (~-55 dB) that the EM1000 could not acquire bottom detect in the areas where the wall of the channel/chute meets the floor. The result was data dropouts and we could find no solution to this problem. It is a minor problem, but an interesting one. The new propellers arrrived and were installed at the end of the day.
August 6, 1998 (980806)
August 7, 1998 (980807)
August 8, 1998 (980808)
August 9, 1998 (980809)
August 10, 1998 (980810)
August 11, 1998 (980811)
August 12, 1998 (980812)
The afternoon was spent at a "stake-holders conference" where we presented the images on a large-screen TV connected to the Silicon Graphics workstation and answered questions. Toward the conclusion of this conference, we asked the stake holders what they wanted us to do with our remaining survey time. A consensus was reached that they wanted us to complete the mapping to the 20-m isobath.
While the conference was in session the Inland Surveyor returned to the lake and spent the afternoon continuing the survey.
August 13, 1998 (980813)
August 14, 1998 (980814)
August 15, 1998 (980815)
August 16, 1998 (980816)
August 17, 1998 (980817)
August 18, 1998 (980818)
The data from August 17th were processed and the final field maps were completed. The computers were packed and the field program ended.
Hughes Clarke, J.E., Mayer, L.A., and Wells, D.E., 1996, Shallow-water imaging multibeam somars: A new tool for investigating seafloor processes in the coastal zone and on the continental shelf. Marine Geophysical Researches, 18: 607-629.
Rowe, T.G. and Stone, J.C., 1997, Selected hydrographic features of the Lake Tahoe Basin, California and Nevada. U.S. Geol. Survey Open-File Rept. 97-384.
U.S. Department of the Interior
|| U.S. Geological Survey
|| National Mapping Information
|| Western Regional Geography|
Page Contact: M.Gould, firstname.lastname@example.org
Last Modified: Wednesday, 08-Mar-2017 16:01:50 EST (tm)
Privacy Statement || Disclaimer || FOIA || Accessibility || Web Team || Feedback || Search