An investigation of upper mantle heterogeneity beneath the Archaean and Proterozoic crust of western Canada from Lithoprobe controlled-source seismic experiments

Observations of upper mantle reflectivity at numerous locations around the world have been linked to the presence of a heterogeneous distribution of rock types within a broad layer of the upper mantle. This phenomenon is observed in wide-angle reflection data from Lithoprobe's Alberta Basement Transect [the SAREX and Deep Probe experiments of 1995] and Trans-Hudson Orogen Transect [the THoRE experiment of 1993]. SAREX and Deep Probe image the Archaean lithosphere of the Hearne and Wyoming Provinces, whereas THoRE images the Archaean and Proterozoic lithosphere of the Trans-Hudson Orogen and neighbouring areas.Finite-difference synthetic seismograms are used to constrain the position and physical properties of the reflective layer. SAREX/Deep Probe modelling uses a 2-D visco-elastic finite-difference routine; THoRE modelling uses a pseudospectral algorithm. In both cases, the upper mantle is parameterized in terms of two media. One medium is the background matrix; the other is statistically distributed within the first as a series of elliptical bodies. Such a scheme is suitable for modelling: (1) variations in lithology (e.g., a peridotite matrix with eclogite lenses) or (2) variations in rheology (e.g., lenses of increased strain within a less strained background).The synthetic seismograms show that the properties of heterogeneities in the upper mantle do not change significantly between the two Lithoprobe transects. Beneath the Trans-Hudson Orogen in Saskatchewan, the layer is best modelled to lie at depths between 80 and 150 km. Based on observations from perpendicular profiles, anisotropy of the heterogeneities is inferred. Beneath the Precambrian domains of Alberta, 400 km to the west, upper mantle heterogeneities are modelled to occur between depths of 90 and 140 km. In both cases the heterogeneous bodies within the model have cross-sectional lengths of tens of kilometers, vertical thicknesses less than 1 km, and velocity contrasts from the background of − 0.3 to − 0.4 km/s. Based on consistency with complementary data and other results, the heterogeneous layer is inferred to be part of the continental lithosphere and may have formed through lateral flow or deformation within the upper mantle.     

Contemporary vertical crustal deformation in Cascadia

Data from 29 tide gauges and 113 pairs of first and second order leveling lines are analyzed to determine the pattern of vertical deformation in the Pacific Northwest of the United States. The data span nearly 100 years and represent the interseismic elastic deformation related to the great earthquake cycle. Uplift rates calculated from leveling surveys are adjusted to a model surface in the tidal reference frame using a robust, weighted, linear, least square technique. Rapid uplift occurs in two distinct broad regions along the coast separated by a narrow zone of slow subsidence. Vertical deformation rates range from > 4 mm/year of uplift on southern Vancouver Island to > 2 mm/year of subsidence in northern coastal Oregon. The deformation pattern is consistent with the results of previous studies and subduction models.     

Equation of state of MgGeO3 perovskite to 65 GPa: comparison with the post-perovskite phase

The equation of state of MgGeO₃ perovskite was determined between 25 and 66 GPa using synchrotron X-ray diffraction with the laser-heated diamond anvil cell. The data were fit to a third-order Birch–Murnaghan equation of state and yielded a zero-pressure volume (V ₀) of 182.2 ± 0.3 Å³ and bulk modulus (K ₀) of 229 ± 3 GPa, with the pressure derivative (K′₀ = (?K ₀/?P) _T ) fixed at 3.7. Differential stresses were evaluated using lattice strain theory and found to be typically less than about 1.5 GPa. Theoretical calculations were also carried out using density functional theory from 0 to 205 GPa. The equation of state parameters from theory (V ₀ = 180.2 Å³, K ₀ = 221.3 GPa, and K′₀ = 3.90) are in agreement with experiment, although theoretically calculated volumes are systematically lower than experiment. The properties of the perovskite phase were compared to MgGeO₃ post-perovskite phase near the observed phase transition pressure (～65 GPa). Across the transition, the density increased by 2.0(0.7)%. This is in excellent agreement with the theoretically determined density change of 1.9%; however both values are larger than those for the (Mg,Fe)SiO₃ phase transition. The bulk sound velocity change across the transition is small and is likely to be negative [?0.5(1.6)% from experiment and ?1.2% from theory]. These results are similar to previous findings for the (Mg,Fe)SiO₃ system. A linearized Birch–Murnaghan equation of state fit to each axis yielded zero-pressure compressibilities of 0.0022, 0.0009, and 0.0016 GPa^?1 for the a, b, and c axis, respectively. Magnesium germanate appears to be a good analog system for studying the properties of the perovskite and post-perovskite phases in silicates.