The stability of a remediated bed in Hamilton Harbour, Lake Ontario, Canada

Abstract

In situ measurements of lakebed sediment erodibility were made on three sites in Hamilton Harbour, Lake Ontario, using the benthic flume Sea Carousel. Three methods of estimating the surface erosion threshold (τ c(0)) from a Carousel time series were evaluated: the first method fits measures of bed strength to eroded depth (the failure envelope) and evaluates threshold as the surface intercept; the second method regresses mean erosion rate (E m) with bed shear stress and solves for the floc erosion rate (E f) to derive the threshold for E m = E f = 1 × 10 -5 kg m -2 s -1; the third method extrapolates a regression of suspended sediment concentration (S) and fluid transmitted bed shear stress (τ o) to ambient concentrations. The first field site was undisturbed (C) and acted as a control; the second (W) was disturbed through ploughing and water injection as part of lakebed treatment, whereas the third site (OIP) was disturbed and injected with an oxidant used for remediation of contaminated sediment. The main objectives of this study were: (1) to evaluate the three different methods of deriving erosion threshold; (2) to compare the physical behaviour of lacustrine sediments with their marine estuarine counterparts; and (3) to examine the effects of ploughing and chemical treatment of contaminated sediment on bed stability. Five deployments of Sea Carousel were carried out at the control site. Mean erosion thresholds for the three methods were: τ c(0) = 0·5 (±0·06), 0·27 (±0·01) and 0·34 (±0·03) Pa respectively. Method 1 overpredicted bed strength as it was insensitive to effects in the surface 1-2 mm, and the fit of the failure envelope was also highly subjective. Method 2 exhibited a wide scatter in the data (low correlation coefficients), and definition of the baseline erosion rate (E f) is largely arbitrary in the literature. Method 3 yielded stable (high correlation coefficients), reproducible and objective results and is thus recommended for evaluation of the erosion threshold. The results of this method correlated well with sediment bulk density and followed the same trend as marine counterparts from widely varying sites. Mass settling rates, expressed as a decay constant, k, of S(t), were strongly related to the maximum turbidity at the onset of settling (S max) and were also in continuity with marine counterparts. Thus, it appears that differences in salinity had little effect on mass settling rates in the examples presented, and that biological activity dominated any effects normally attributable to changes in salinity. Bedload transport of eroded aggregates (2-4 mm in diameter) took place by rolling below a mean tangential flow velocity (U y) of 0·32 ms -1 and by saltation at higher velocities. Mass transport as bedload was a maximum at U y = 0·4 ms -1, although bedload never exceeded 1% of the suspended load. The proportion of material moving as bedload was greatest at the onset of erosion but decreased as flow competence increased. Given the low bulk density and strength of the lakebed sediment, the presence of a bedload component is notable. Bedload transport over eroding cohesive substrates should be greater in estuaries, where both sediment density and strength are usually higher. Significant differences between the ploughed and control sites were apparent in both the erosion rate and the friction coefficient (∅), and suggest that bed recovery after disruption is rapid (< 24 h). τ c(0) increased linearly with time after ploughing and recovered to the control mean value within 3 days. The friction coefficient was reduced to zero by ploughing (diagnostic of fluidization), but increased linearly with time, regaining control values within 6 days. No long-term reduction in bed strength due to remediation was apparent.