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 (τ₀) 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.     

The electrical structure of the Slave craton

The Slave craton in northwestern Canada, a relatively small Archean craton (600×400 km), is ideal as a natural laboratory for investigating the formation and evolution of Mesoarchean and Neoarchean sub-continental lithospheric mantle (SCLM). Excellent outcrop and the discovery of economic diamondiferous kimberlite pipes in the centre of the craton during the early 1990s have led to an unparalleled amount of geoscientific information becoming available.

Over the last 5 years deep-probing electromagnetic surveys were conducted on the Slave, using the natural-source magnetotelluric (MT) technique, as part of a variety of programs to study the craton and determine its regional-scale electrical structure. Two of the four types of surveys involved novel MT data acquisition; one through frozen lakes along ice roads during winter, and the second using ocean-bottom MT instrumentation deployed from float planes.

The primary initial objective of the MT surveys was to determine the geometry of the topography of the lithosphere–asthenosphere boundary (LAB) across the Slave craton. However, the MT responses revealed, completely serendipitously, a remarkable anomaly in electrical conductivity in the SCLM of the central Slave craton. This Central Slave Mantle Conductor (CSMC) anomaly is modelled as a localized region of low resistivity (10–15 Ω m) beginning at depths of 80–120 km and striking NE–SW. Where precisely located, it is spatially coincident with the Eocene-aged kimberlite field in the central part of the craton (the so-called “Corridor of Hope”), and also with a geochemically defined ultra-depleted harzburgitic layer interpreted as oceanic or arc-related lithosphere emplaced during early tectonism. The CSMC lies wholly within the NE–SW striking central zone defined by Grütter et al. [Grütter, H.S., Apter, D.B., Kong, J., 1999. Crust–mantle coupling; evidence from mantle-derived xenocrystic garnets. Contributed paper at: The 7th International Kimberlite Conference Proceeding, J.B. Dawson Volume, 1, 307–313] on the basis of garnet geochemistry (G10 vs. G9) populations.

Deep-probing MT data from the lake bottom instruments infer that the conductor has a total depth-integrated conductivity (conductance) of the order of 2000 Siemens, which, given an internal resistivity of 10–15 Ω m, implies a thickness of 20–30 km. Below the CSMC the electrical resistivity of the lithosphere increases by a factor of 3–5 to values of around 50 Ω m. This change occurs at depths consistent with the graphite–diamond transition, which is taken as consistent with a carbon interpretation for the CSMC.

Preliminary three-dimensional MT modelling supports the NE–SW striking geometry for the conductor, and also suggests a NW dip. This geometry is taken as implying that the tectonic processes that emplaced this geophysical–geochemical body are likely related to the subduction of a craton of unknown provenance from the SE (present-day coordinates) during 2630–2620 Ma. It suggests that the lithospheric stacking model of Helmstaedt and Schulze [Helmstaedt, H.H., Schulze, D.J., 1989. Southern African kimberlites and their mantle sample: implications for Archean tectonics and lithosphere evolution. In Ross, J. (Ed.), Kimberlites and Related Rocks, Vol. 1: Their Composition, Occurrence, Origin, and Emplacement. Geological Society of Australia Special Publication, vol. 14, 358–368] is likely correct for the formation of the Slave's current SCLM.     

Petrological and geochronological constraints on high pressure, high temperature metamorphism in the Snowbird tectonic zone, Canada

The upper deck of the East Athabasca mylonite triangle (EAmt), northern Saskatchewan, Canada, contains mafic granulites that have undergone high P–T metamorphism at conditions ranging from 1.3 to 1.9 GPa, 890–960 °C. Coronitic textures in these mafic granulites indicate a near‐isothermal decompression path to 0.9 GPa, 800 °C. The Godfrey granite occurs to the north adjacent to the upper deck high P–T domain. Well‐preserved corona textures in the Godfrey granite constrain igneous crystallization and early metamorphism in the intermediate‐pressure granulite field (Opx + Pl) at 1.0 GPa, 775 °C followed by metamorphism in the high pressure granulite field (Grt + Cpx + Pl) at 1.2 GPa, 860 °C. U–Pb geochronology of zircon in upper deck mafic granulite yields evidence for events at both c. 2.5 Ga and c. 1.9 Ga. The oldest zircon dates are interpreted to constrain a minimum age for crystallization or early metamorphism of the protolith. A population of 1.9 Ga zircon in one mafic granulite is interpreted to constrain the timing of high P–T metamorphism. Titanite from the mafic granulites yields dates ranging from 1900 to 1894 Ma, and is interpreted to have grown along the decompression path, but still above its closure temperature, indicating cooling following the high P–T metamorphism from c. 960–650 °C in 4–10 Myr. Zircon dates from the Godfrey granite indicate a minimum crystallization age of 2.61 Ga, without any evidence for 1.9 Ga overgrowths. The data indicate that an early granulite facies event occurred at c. 2.55–2.52 Ga in the lower crust (c. 1.0 GPa), but at 1.9 Ga the upper deck underwent high P–T metamorphism, then decompressed to 0.9–1.0 GPa. Juxtaposition of the upper deck and Godfrey granite would have occurred after or been related to this decompression. In this model, the high P–T rocks are exhumed quickly following the high pressure metamorphism. This type of metamorphism is typically associated with collisional orogenesis, which has important implications for the Snowbird tectonic zone as a fundamental boundary in the Canadian Shield.     

Eastern Laurentia in Rodinia: constraints from whole-rock Pb and U/Pb geochronology

Whole-rock Pb isotopic signatures and U/Pb geochronology refute a Rodinian correlation of northeastern Laurentia and proto-Andean Amazonia. According to this previously proposed model, the Labrador–Scotland–Greenland Promontory (LSGP) of northeastern Laurentia collided with the proto-Andean margin of Amazonia, at the Arica Embayment, during the Grenville/Sunsás Orogeny (ca. 1.0 Ga). Links between the two margins were based upon the correlation of the LSGP with Arequipa-Antofalla Basement (AAB), a Proterozoic block along the proto-Andean margin of Amazonia adjacent to the Arica Embayment. Specifically, similarities in 1.8–1.0 Ga basement rocks in both regions suggested that the AAB was originally a piece of the LSGP. Furthermore, similarities in unique, post-collisional, but pre-rift, glacial sedimentary sequences also supported a link between the AAB and LSGP.Tests of these apparent similarities fail to support correlation of the AAB and the LSGP and, thus, eliminate a direct link between northeastern Laurentia and southwestern Amazonia in Rodinia. However, Pb isotopic compositions and U/Pb geochronology provide the basis for two new correlations, namely, (1) the ca. 1.3–1.0 Ga basement in the central and southern Appalachians may be an allochthonous block that was transferred to Laurentia from Amazonia at ca. 1.0 Ga, and (2) an allochthonous AAB may be a piece of the Kalahari Craton that was transferred to Amazonia at ca. 1.0 Ga. Based on these new correlations and a previously proposed Grenvillian connection between southern Laurentia (Llano) and Kalahari, we propose that Amazonia may have collided with a contiguous southeastern Laurentia/Kalahari margin at ca. 1.0 Ga.     

A statistical approach to evaluating distance metrics and analog assignments for pollen records

The modern analog technique typically uses a distance metric to determine the dissimilarity between fossil and modern biological assemblages. Despite this quantitative approach, interpretation of distance metrics is usually qualitative and rules for selection of analogs tend to be ad hoc. We present a statistical tool, the receiver operating characteristic (ROC) curve, which provides a framework for identifying analogs from distance metrics. If modern assemblages are placed into groups (e.g., biomes), this method can (1) evaluate the ability of different distance metrics to distinguish among groups, (2) objectively identify thresholds of the distance metric for determining analogs, and (3) compute a likelihood ratio and a Bayesian probability that a modern group is an analog for an unknown (fossil) assemblage. Applied to a set of 1689 modern pollen assemblages from eastern North America classified into eight biomes, ROC analysis confirmed that the squared-chord distance (SCD) outperforms most other distance metrics. The optimal threshold increased when more dissimilar biomes were compared. The probability of an analog vs no-analog result (a likelihood ratio) increased sharply when SCD decreased below the optimal threshold, indicating a nonlinear relationship between SCD and the probability of analog. Probabilities of analog computed for a postglacial pollen record at Tannersville Bog (Pennsylvania, USA) identified transitions between biomes and periods of no analog.     

Atmospheric deposition of Be and Be in New Zealand rain (1996-98)

Beryllium isotope concentrations were determined in monthly rainfall collections at three sites across New Zealand (36 to 45° S), from October 1996 to November 1998. At the northern sites of Leigh (near Auckland) and Gracefield (near Wellington), ⁷Be and ¹⁰Be concentrations are relatively constant at 1.2 to 1.4 × 10⁷ atoms kg⁻¹ rain and 2.1 to 2.6 × 10⁷ atoms kg⁻¹ rain, respectively. These concentrations correspond to annual flux rates of ∼15 and ∼27 × 10⁹ atoms m⁻² y⁻¹, respectively. At the southern site of Dunedin, concentrations are similar to those at the northern sites, but flux rates are significantly lower at ∼ 9 and ∼19 × 10⁹ atoms m⁻² y⁻¹, respectively, because of lower average rainfall east of the Southern Alps mountain range. These results are broadly similar to those reported by Brown et al. (1989) and Knies et al. (1994) for rain from midlatitude sites in the USA sampled from 1986 to 1994. Mean ⁷Be/¹⁰Be ratios for New Zealand (0.47 to 0.61) are, however, significantly lower than for the USA (0.69 to 0.78), due in part to the addition of ¹⁰Be from re-suspended dust. Subtraction of the dust component increases the New Zealand ⁷Be/¹⁰Be ratios to 0.70 (Leigh), 0.65 (Gracefield) and 0.50 (Dunedin). The adjusted results provide evidence for transfer of older stratospheric air to the troposphere in late-spring-summer, an effect which is strongest in the north. The overall reduction of ⁷Be/¹⁰Be from north to south implies an increase in residence time from ∼ 80 to ∼100 d for Be isotopes in the atmosphere above New Zealand.     

UV Laser Ablation ICP-MS: Some Applications in the Earth Sciences

This study reports a series of applications of UV laser ablation ICP-MS in the geological sciences. The advantages and disadvantages of the PQ "S" option and the use of nitrogen in the carrier gas are discussed. A general problem common to all ablation techniques is the calibration technique and experiments involving synthetic calibration samples are covered. Zircon geochemistry and geochronology by LA-ICP-MS are discussed and data are presented for REE, Hf and U for a standard zircon (91500) as well as a series of zircons from Zimbabwe. The potential of using Ce and Eu anomalies in petrologic studies is illustrated by zircons from a fractionated gabbroic-granite in the Urals. The potential of the LA-ICP-MS method to utilise standard X-ray fluorescence glass discs is demonstrated as a useful semi-quantitative tool in determining REE patterns. LA-ICP-MS is a powerful tool in the analysis of the platinum group elements (and Re) and some examples are given in the successful application of the technique to partitioning in iron meteorites.     

Shear-wave anisotropy: spatial and temporal variations in time delays at Parkfield, Central California

Shear-wave splitting is analysed on data recorded by the High Resolution Seismic Network (HRSN) at Parkfield on the San Andreas fault, Central California, during the three-year period 1988-1990. Shear-wave polarizations either side of the fault are generally aligned in directions consistent with the regional horizontal maximum compressive stress, at some 70° to the fault strike, whereas at station MM in the immediate fault zone, shear-wave polarizations are aligned approximately parallel to the fault. Normalized time delays at this station are found to be about twice as large as those in the rock mass either side. This suggests that fluid-filled cracks and fractures within the fault zone are elastically or seismically different from those in the surrounding rocks, and that the alignment of fault-parallel shear-wave polarizations are associated with some fault-specific phenomenon.
Temporal variations in time delays between the two split shear-waves before and after a M_L = 4 earthquake can be identified at two stations with sufficient data: MM within the fault zone and VC outside the immediate fault zone. Time delays between faster and slower split shear waves increase before the M_L = 4 earthquake and decrease near the time of the event. The temporal variations are statistically significant at 68 per cent confidence levels. Earthquake doublets and multiplets also show similar temporal variations, consistent with those predicted by anisotropic poroelasticity theory for stress modifications to the microcrack geometry pervading the rock mass. This study is broadly consistent with the behaviour observed before three other earthquakes, suggesting that the build-up of stress before earthquakes may be monitored and interpreted by the analysis of shear-wave splitting.