Structure of the crust and upper mantle in Canada

The seismic probing of the crust and upper mantle in Canada started in 1938 and since then has involved many government and university groups using a wide variety of techniques. These have included simple profiling with both wide and narrow station spacing, areal time-term surveys, detailed deep reflection experiments, very long-range refraction studies and the analysis of surface wave dispersion between stations of the Canadian Standard Network.

A review of the published interpretation leads to the general conclusion that:

1. (1) P_n-velocities vary from a value possibly as low as 7.7 km/sec under Vancouver Island to 8.6 km/sec and higher in the extreme eastern part of the shield and some parts of the Atlantic coast.

2. (2) Large areas of Canada have a crustal thickness of 30–40 km, with Vancouver Island, the southwestern Prairies, the Lake Superior basin and parts of the eastern shield of Quebec being thicker. No continental area in Canada is known to have a crust thinner than 29 km.

3. (3) The Riel discontinuity — a deep intra-crustal reflector and sometime refractor, is widely reported in the Prairies and Manitoba. It is not seen to the north in the vicinity of Great Slave Lake, nor in the Hudson Bay, Lake Superior and Maritime regions, nor in the interior of British Columbia. It may be present in some areas of the eastern shield.

4. (4) As experiments have become more detailed, crustal structures of greater complexity have been revealed. The concept that crustal structure becomes simpler with increasing depth is apparently unfounded.

Long-range refraction studies suggest that the Gutenberg P-wave low-velocity channel is poorly developed under the Canadian Shield. The analysis of the dispersion of surface waves, however, suggests that the channel is better developed for S-waves, and is present throughout the country. The lid of the channel is deepest under the central shield and shallowest under the Cordillera.     

The variation of tetrahedral bond lengths in sodic plagioclase feldspars

Multiple linear regression analysis has been applied to the geometric and chemical variables in sodic plagioclases in order to determine their relative effects on individual T-O bond lengths in the Al_1+xSi_3?xO₈ tetrahedral framework. Using data from crystal structure analyses of low and high albite, An₁₆ and An₂₈, and assuming that low albite is completely ordered, 1 $$\begin{gathered} {\text{T}} - {\text{O = 1}}{\text{.568}} + {\text{[(0}}{\text{.122) x (Al content of the T site)]}} \hfill \\ {\text{ }} - {\text{[(0}}{\text{.037) x (}}\Delta {\text{{\rm A}l}}_{{\text{br}}} )] + [0.063){\text{ x }}(\Sigma {\text{[}}q{\text{/(Na,Ca}} - {\text{O)}}^{\text{2}} ])] \hfill \\ {\text{ }} + {\text{[(0}}{\text{.029) x (}} - {\text{1/cosT}} - {\text{O}} - {\text{T)]}} \hfill \\ \end{gathered}$$ where the Al content of a particular tetrahedral (T) site can be estimated from empirically-derived determinative curves, where Δ Al_br is a linkage factor to account for the Al content of adjacent tetrahedral sites, where the formal charge on the (Na_1?xCa_x) atom is q=1+x, and where T-O-T is the inter-tetrahedral angle involving the T-O bond. For sodic plagioclases it is essential to know only the anorthite content and the 2Θ₁₃₁-2Θ_1¯31 spacing (CuK _α radiation) in order to determine the independent variables in this equation and thus to evaluate the individual T-O distances. The 64 individual T-O distances predicted for the four sodic plagioclases by this equation agree well with the observed T-O bond lengths (σ=0.004 Å; r=0.994), and the method has been used by way of example to rationalize the T-O bond lengths in analcime (cf. Ferraris, Jones and Yerkess, 1972).     

An ion-association model for PbO-SiO2 melts: interpretation of thermochemical,conductivity, and density data

We propose that within PbO-SiO₂ melts the component species PbO(1) and SiO₂(1) associate to form a lead orthosilicate complex [Pb₂SiO₄(1)] according to $\text{2PbO(1) + SiO}{}_2\text{(1)}\text{a}\text{?}\text{f Pb}{}_2\text{SiO}{}_4\text{(1)}$ The association constant (K) for the above reaction is near 30 at 900 and 1000°C On the premises that the enthalpy of the above reaction is the only contributor to the enthalpy of mixing and that three melt species [PbO(1), SiO₂(1), Pb₂SiO₄(1)] mix ideally (Raoult's Law obeyed), the partial molar properties and the integral mixing properties of the melt have been calculated and have been found to accurately reproduce the measured thermochemical data.The conductivity data and density data for PbO-SiO₂ melts have been calculated using the association model. The calculated data again are remarkably consistent with the experimental data.Other silicate-melt models have been reinterpreted and that of Toop and Samis is compatible with the proposed model. In fact the Toop and Samis model implicitly assumes the formation of orthosilicate species rather than condensed polymeric species.     

Abundance of Boron in International Geochemical Standards by Prompt-Gamma Neutron Activation Analysis

Eighteen international reference standards, mostly geochemical standards, have been analysed for boron using prompt-gamma neutron activation analysis (PGNAA). There is good agreement between the results for samples with high boron abundances and the few published values using PGNAA and other methods. However, correlation is poor at low boron abundances, especially in those standards with low boron/sodium ratios. Corrections for sodium interference on boron have been critically evaluated and the data corrected using two different methods.     

Mechanism of continental crust subsidence in the Alpine Belt

Formation of deep basins on continental crust in fold belts is often explained by stretching. This mechanism inevitably produces large deformations in the upper crust. No deformations typical of significant stretching were revealed in the predominant part of deep basins on continental crust in the Alpine Belt. This means that these basins were not produced by stretching. Most basins were formed during a short period of time of a few million years. The short duration of the subsidences eliminates thermal relaxation as the mechanism. The space and time relationships between the subsidence and orogeny and the profile of the basin floor exclude thrust loading as a cause of formation for practically all large basins. Gabbro to eclogite transformation is suggested as a mechanism of rapid subsidence. This occurs under the upwelling of hydrous asthenosphere at moderate temperature to the base of the crust. Eclogite sinking into the mantle results in a strong attenuation of the crust and lithosphere, which permits intense subsequent folding. The major part of deep basins in continental crust that formed by rapid subsidence was intensely shortened in the Alpine Belt. Significant crustal shortening did not spread over the cratonic lithosphere.     

Inter‐laboratory Characterisation of Apatite Reference Materials for Chlorine Isotope Analysis

Here we report on a set of six apatite reference materials (chlorapatites MGMH#133648, TUBAF#38 and fluorapatites MGMH#128441A, TUBAF#37, 40, 50) which we have characterised for their chlorine isotope ratios; these RMs span a range of Cl mass fractions within the apatite Ca₁₀(PO₄)₆(F,Cl,OH)₂ solid solution series. Numerous apatite specimens, obtained from mineralogical collections, were initially screened for ³⁷Cl/³⁵Cl homogeneity using SIMS followed by δ³⁷Cl characterisation by gas source mass spectrometry using both dual‐inlet and continuous‐flow modes. We also report major and key trace element compositions as determined by EPMA. The repeatability of our SIMS results was better than ± 0.10‰ (1s) for the five samples with > 0.5 % m/m Cl and ± 0.19‰ (1s) for the low Cl abundance material (0.27% m/m). We also observed a small, but significant crystal orientation effect of 0.38‰ between the mean ³⁷Cl/³⁵Cl ratios measured on three oriented apatite fragments. Furthermore, the results of GS‐IRMS analyses show small but systematic offset of δ³⁷Cl_SMOC values between the three laboratories. Nonetheless, all studied samples have comparable chlorine isotope compositions, with mean 10³δ³⁷Cl_SMOC values between +0.09 and +0.42 and in all cases with 1s ≤ ± 0.25.     

Modeling the natural degradation of earthworks

Stratigraphic, topographic, and ground‐penetrating radar data obtained from a ca. 1800‐year‐old embankment and adjacent ditch at the Hopewell Mound Group (Chillicothe, Ohio) are used to validate the archaeological application of a simple finite‐differences diffusion model employed frequently to assess geomorphic change in natural landscapes. Although diffusion models have been used to describe the topographic degradation of landforms in a variety of geomorphic terrains, the approach has not been applied to ancient earthworks in an archaeological context. The results of this study indicate that a variety of initial earthwork forms can result in the sinusoidal profile apparent on the current landscape. Using the model results to interpret the field data, we suggest that the initial embankment form was steeper and the adjacent ditch was deeper. As a result of natural degradation processes, the earthwork widened and flattened over time. These results have broad implications for any study aimed at: (1) assessing the function of original earthwork forms, (2) determining the formation processes of complicated stratigraphies or artifact assemblages, (3) estimating the time and labor investment required for construction, or (4) identifying the socio‐political structures necessary to build earthworks. © 2005 Wiley Periodicals, Inc.     

Hydrogen isotopic composition of water from fossil micrometeorites in howardites

We have measured the hydrogen isotopic composition (D/H ratios) of the water from 13 carbonaceous chondritic microclasts (CCMs, size <1 mm) trapped in two howardites (Kapoeta and Yamato-793497) early in the evolution of Solar System. The division into tochilinite-rich; magnetite-rich, olivine-poor; magnetite-rich, olivine-rich CCM types is corroborated by the hydrogen isotopic compositions. Both mineralogy and hydrogen isotopic compositions demonstrate that tochilinite-rich CCMs represent CM2 chondritic matter. In contrast, there is no good match between the isotopic and mineralogical properties of the magnetite-rich CCMs and the known groups of carbonaceous chondrites, suggesting that magnetite-rich CCMs represent a new kind of chondritic matter, not yet sampled in meteorite collections. This demonstrates that the view of the asteroid belt revealed by the collection of meteorites is incomplete. The study of (micro)clasts offers a unique opportunity to better decipher the nature and relative abundance of asteroids.The average hydrogen isotopic composition of water belonging to CCMs, D/H = (152.0 ± 4.8) × 10⁻⁶ (1σ_m), is similar to that of Antarctic micrometeorites (AMMs), D/H = (161.2 ± 3.8) × 10⁻⁶ (1σ_m). The similarity, in terms of mineralogy and hydrogen isotopic composition, between CCMs and AMMs demonstrates that the composition of the micrometeorites has not been modified over the whole history of the Solar System. It indicates that the composition of the micrometeorite flux onto Earth has been, and is, dominated by a mixture of CM2-like; magnetite-rich, olivine-poor; magnetite-rich, olivine-rich carbonaceous chondritic matter exemplified by CCMs found in howardites. Because CCMs have not suffered atmospheric entry, they provide an abundant source of pristine micrometeorites.The average D/H ratio of the whole population of CCMs is identical within errors to that of the Earth (149 ± 3 × 10⁻⁶). The match between the CCMs D/H ratio and that of the Earth is especially remarkable because 1) three different populations of CCMs are needed to make the D/H ratio of the Earth; 2) there is no single carbonaceous chondrite group for which a similar match exists. This observation suggests that CCMs population might be representative of the late veneer agent(s) that delivered water to the Earth.     

Site preservation along an active meandering and avulsing river: The Red River,Arkansas

The low-gradient Red River is a rapidly migrating, sinuous stream with easily erodible banks. Avulsion is common at many scales, from individual meander bends that are cut off to major sections of the river that form multiple, complex meander belts. The present meander belt can be subdivided into mappable landforms—termed phases—that are associated with river courses of different ages and thus associated with archeological sites of different ages. Within the study area two phases are present. The younger Modern meander belt phase has formed within the past 0.2–0.3 ky, precluding preservation of prehistoric archaeological sites. Any protohistoric artifacts that may have been preserved in this meander belt phase would be deeply buried because as much as 2 m of the vertical accretion sediment has accumulated between artificial levees in <0.1 ky and 1–2 m of sediment has accumulated beyond the artificial levees in <0.2 ky. Archeological site preservation in this highly mobile fluvial end member can be used as a predictor for other, similar streams. A large prehistoric site is preserved on an older (0.5–1 kya) Late Prehistoric meander belt phase associated with an abandoned river course. In the study area a Fourche Maline 7 period (A.D. 800–900) through Caddo IV period (ca. A.D. 1500–1700) archeological site (3MI3/30) is preserved on this slightly higher altitude portion of the flood plain. At locations proximal to the river, the site may be buried by overbank sediment 0.4 m thick, but at more distant locations the site is at the surface or only buried by thin overbank sediment because of low sedimentation rates (0.04 cm yr⁻¹) over the span of a millennium. Sites, such as 3MI3/30, that are occupied contemporaneous with overbank sedimentation may be stratified; however, localized erosion and removal of some archeological material may occur where channelized flow crosses the natural levee. © 1998 John Wiley & Sons, Inc.     

Geochemically induced shifts in catabolic energy yields explain past ecological changes of diffuse vents in the East Pacific Rise 9°50'N area

The East Pacific Rise (EPR) at 9°50'N hosts a hydrothermal vent field (Bio9) where the change in fluid chemistry is believed to have caused the demise of a tubeworm colony. We test this hypothesis and expand on it by providing a thermodynamic perspective in calculating free energies for a range of catabolic reactions from published compositional data. The energy calculations show that there was excess H₂S in the fluids and that oxygen was the limiting reactant from 1991 to 1997. Energy levels are generally high, although they declined in that time span. In 1997, sulfide availability decreased substantially and H₂S was the limiting reactant. Energy availability dropped by a factor of 10 to 20 from what it had been between 1991 and 1995. The perishing of the tubeworm colonies began in 1995 and coincided with the timing of energy decrease for sulfide oxidizers. In the same time interval, energy availability for iron oxidizers increased by a factor of 6 to 8, and, in 1997, there was 25 times more energy per transferred electron in iron oxidation than in sulfide oxidation. This change coincides with a massive spread of red staining (putative colonization by Fe-oxidizing bacteria) between 1995 and 1997. For a different cluster of vents from the EPR 9°50'N area (Tube Worm Pillar), thermodynamic modeling is used to examine changes in subseafloor catabolic metabolism between 1992 and 2000. These reactions are deduced from deviations in diffuse fluid compositions from conservative behavior of redox-sensitive species. We show that hydrogen is significantly reduced relative to values expected from conservative mixing. While H₂ concentrations of the hydrothermal endmember fluids were constant between 1992 and 1995, the affinities for hydrogenotrophic reactions in the diffuse fluids decreased by a factor of 15 and then remained constant between 1995 and 2000. Previously, these fluids have been shown to support subseafloor methanogenesis. Our calculation results corroborate these findings and indicate that the 1992-1995 period was one of active growth of hydrogenotrophic communities, while the system was more or less at steady state between 1995 and 2000.