Question 1 of 5: Verifying a Bathymetric LiDAR System

Question 1: Verifying a Bathymetric LiDAR System (25 marks)

Question text not reproduced: the examination questions are © Engineers and Geoscientists BC. Open the official past paper (linked at the top of this page) to read the question, then follow the worked solution below.

An airborne bathymetric LiDAR (ALB) fires a pulsed green laser (532 nm) that penetrates the water column together with a coincident near-infrared (1064 nm) pulse that reflects off the surface; depth is recovered from the time separation of the surface and bottom returns, $d=\tfrac{1}{2}\,v_w\,\Delta t$ (with $v_w\approx0.75c$ in water and a geometric correction for the off-nadir scan angle and refraction at the interface). Verification is an independent-check exercise: you compare the system against ground truth that is more accurate than the system itself, over the full range of conditions the specification claims.

(a) Verifying accuracy and depth penetration. Accuracy verification is done against independent, higher-order reference depths.

Given. A worked check point at depth $d=15\ \text{m}$, an IHO S-44 Order 1a specification ($a=0.5\ \text{m}$, $b=0.013$), and a survey-time Secchi depth $Z_S=6\ \text{m}$. Find. the depth-accuracy tolerance the residuals must satisfy and the expected depth penetration.

1. Fly a calibration / verification site of known bathymetry. Choose an area already surveyed by acoustic multibeam or a leveled reference, spanning shallow to deep water and a range of bottom types and slopes. Include a hard, flat, well-defined surface (e.g. a concrete ramp or a leveled bar) whose height is known from geodetic leveling.
2. Reference every sounding to a common vertical datum. Reduce both the LiDAR and the reference depths to the same chart datum (via GNSS-derived aircraft trajectory + a geoid/SEP model, or a tide model), so the comparison is not corrupted by tide or trajectory error.
3. Compute depth differences and compare against the IHO tolerance. Form the LiDAR-minus-reference residuals at each check point; their mean tests for a bias (a sound-speed/refraction or timing offset) and their spread tests precision. The allowable total vertical uncertainty follows the IHO S-44 rule $$\text{TVU}=\sqrt{a^{2}+(b\,d)^{2}}$$ For Order 1a ($a=0.5\ \text{m},\ b=0.013$) at a check depth $d=15\ \text{m}$: $$\text{TVU}=\sqrt{0.5^{2}+(0.013\times 15)^{2}}=\boxed{0.54\ \text{m}}$$ so 95% of the residuals at 15 m must fall within $\pm0.54$ m for the system to meet Order 1a.
4. Fly reciprocal and cross lines. Overlapping and opposing-direction lines expose direction-dependent errors (scanner-angle refraction, heading/attitude bias); the crossover differences at line intersections are an internal accuracy check independent of any reference survey.

Depth penetration is verified by flying progressively deeper water of known depth until the bottom return can no longer be reliably detected, and recording the greatest depth at which a valid bottom is still resolved. Penetration is governed by water clarity, so it must be reported relative to clarity — the practical rule is that maximum penetration is roughly two-to-three Secchi depths:

5. Tie penetration to a clarity measurement. Measure the Secchi depth $Z_S$ at survey time. With $Z_S=6\ \text{m}$ the expected penetration is $$d_{\max}\approx (2\text{ to }3)\,Z_S = \boxed{12\text{ to }18\ \text{m}}$$ Confirm that the deepest reliably-detected known depth falls in this band, and report penetration together with the diffuse attenuation coefficient $K_d$ (or Secchi depth) so the number is meaningful.

(b) Verifying hazard-detection ability. Hazard detection is the ability to resolve a small least-depth over an obstruction (a shoal, wreck, or boulder) — a spatial-resolution and detection question, not an accuracy question. Verify it empirically: deploy known targets of known size and least depth (calibrated spheres, a sunken frame, or a mapped small shoal) and confirm the system detects them and reports their least depth correctly. Vary target size and depth to establish the smallest cube/least-depth the system can capture at the operational point spacing and altitude (this is the "capability to detect features of size X" specification). Because a bathymetric LiDAR footprint is metres wide and the sounding density is coarser than multibeam, small isolated hazards can fall between shots — so the verification must be at the operational spot spacing, and any feature detection claim is a statistical one (probability of detection versus feature size).

(c) Seafloor classification with LiDAR. Only in a limited, coarse sense. The green-laser waveform does carry some information about the bottom: the amplitude and shape of the bottom return depend on bottom reflectance and roughness, so a high-reflectance sandy bottom returns a stronger, sharper pulse than dark mud or vegetation, and some systems classify substrate crudely from bottom-return intensity or the presence of a volume (vegetation) echo. But this is unreliable for true classification: the return amplitude is confounded by water-column attenuation, depth, scan angle and slope, all of which change the signal independently of the bottom type. For dependable seafloor classification you use acoustic backscatter (multibeam or sidescan), whose intensity is a direct function of the sediment's acoustic impedance and roughness. So: LiDAR can give a coarse, qualitative substrate hint, but it is not a substitute for acoustic backscatter classification.