Influence of random heterogeneity of shear wave velocity on sliding mass response and seismic deformations of earth slopes

Soil shear wave velocity has been recognized as a governing parameter in the assessment of the seismic response of slopes. The spatial variability of soil shear wave velocity can influence the seismic response of sliding mass and seismic displacements. However, most analyses of sliding mass response have been carried out by deterministic models. This paper stochastically investigates the effect of random heterogeneity of shear wave velocity of soil on the dynamic response of sliding mass using the correlation matrix decomposition method and Monte Carlo simulation(MCS). The software FLAC 7.0 along with a Matlab code has been utilized for this purpose. The influence of statistical parameters on the seismic response of sliding mass and seismic displacements in earth slopes with different inclinations and stiffnesses subject to various earthquake shakings was investigated. The results indicated that, in general, the random heterogeneity of soil shear modulus can have a notable impact on the sliding mass response and that neglecting this phenomenon could lead to underestimation of sliding deformations.     

A Nonlinear Model of Surf Beat

A nonlinear short-wave-averaged (surf beat) model is presented. The model is based on that of Roelvink (1993), but the numerical techniques used in the solution are based on the so-called weighted-averaged flux (WAF) method (eg Watson et al., 1992), with time-operator splitting used for the treatment of some of the source terms. This method allows a small number of computational points to be used, and is particularly efficient in modelling breaking long waves. The short-wave (or primary-wave) energy equation is solved using a more traditional Lax-Wendroff technique. Results of validation indicate that the model performs satisfactorily in most respects.     

The problem of underpowered rivers

This study has hypothesized that for many rivers the trade-off between flow accumulation and the decrease in slope along channel length means that stream power increases downstream and, moreover, that given the low slope angles in headwater and low-order streams, they would have insufficient stream power to erode let alone transport sediment. The study considered the stream power profile, the particle travel distances and the application of the Hjulström curve based on the velocity profile of nine, large UK catchments. The study showed that:

1. Some rivers never showed a maximum in their longitudinal stream power profile, implying that some rivers never develop a deposition zone before they discharge at the tidal limit.
2. Particle travel distances during a bankfull discharge event showed that for some rivers 91% of the upper main channel would not be cleared of sediment. Furthermore, while some rivers could transport a 2 mm particle their entire length in one bankfull event, for another river it would take 89 such events.
3. The Hjulström curve shows that for three of the study rivers the upper 20 km of the river was not capable of eroding a 2 μm particle.
4. The study has shown that for all rivers studied, erosion is focused downstream and deposition upstream. Many UK rivers have a dead zone where, on time scales of the order of centuries, no erosion or transport occurs and erosion only occurs in the lower courses of the channel where discharge rather than slope dominates – we propose these as underpowered rivers.

© 2020 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd     

Upstairs-downstairs: supercontinents and large igneous provinces,are they related?

There is a correlation of global large igneous province (LIP) events with zircon age peaks at 2700, 2500, 2100, 1900, 1750, 1100, and 600 and also probably at 3450, 3000, 2000, and 300 Ma. Power spectral analyses of LIP event distributions suggest important periodicities at 250, 150, 100, 50, and 25 million years with weaker periodicities at 70–80, 45, and 18–20 Ma. The 25 million year periodicity is important only in the last 300 million years. Some LIP events are associated with granite-forming (zircon-producing) events and others are not, and LIP events at 1900 and 600 Ma correlate with peaks in craton collision frequency. LIP age peaks are associated with supercontinent rifting or breakup, but not dispersal, at 2450–2400, 2200, 1380, 1280, 800–750, and ≤200 Ma, and with supercontinent assembly at 1750 and 600 Ma. LIP peaks at 2700 and 2500 Ma and the valley between these peaks span the time of Neoarchaean supercraton assemblies. These observations are consistent with plume generation in the deep mantle operating independently of the supercontinent cycle and being controlled by lower-mantle and core-mantle boundary thermochemical dynamics. Two processes whereby plumes can impact continental assembly and breakup are (1) plumes may rise beneath supercontinents and initiate supercontinent breakup, and (2) plume ascent may increase the frequency of craton collisions and the rate of crustal growth by accelerating subduction.     

A non-hydrostatic model for the numerical study of landslide-generated waves

This paper proposes and demonstrates a two-layer depth-averaged model with non-hydrostatic pressure correction to simulate landslide-generated waves. Landslide (lower layer) and water (upper layer) motions are governed by the general shallow water equations derived from mass and momentum conservation laws. The landslide motion and wave generation/propagation are separately formulated, but they form a coupled system. Our model combines some features of the landslide analysis model DAN3D and the tsunami analysis model COMCOT and adds a non-hydrostatic pressure correction. We use the new model to simulate a 2007 rock avalanche-generated wave event at Chehalis Lake, British Columbia, Canada. The model results match both the observed distribution of the rock avalanche deposit in the lake and the wave run-up trimline along the shoreline. Sensitivity analyses demonstrate the importance of accounting for the non-hydrostatic dynamic pressure at the landslide-water interface, as well as the influence of the internal strength of the landslide on the size of the generated waves. Finally, we compare the numerical results of landslide-generated waves simulated with frictional and Voellmy rheologies. Similar maximum wave run-ups can be obtained using the two different rheologies, but the frictional model better reproduces the known limit of the rock avalanche deposit and is thus considered to yield the best overall results in this particular case.     

Reaction mechanism for the replacement of calcite by dolomite and siderite: implications for geochemistry, microstructure and porosity evolution during hydrothermal mineralisation

Carbonate reactions are common in mineral deposits due to CO₂-rich mineralising fluids. This study presents the first in-depth, integrated analysis of microstructure and microchemistry of fluid-mediated carbonate reaction textures at hydrothermal conditions. In doing so, we describe the mechanisms by which carbonate phases replace one another, and the implications for the evolution of geochemistry, rock microstructures and porosity. The sample from the 1.95 Moz Junction gold deposit, Western Australia, contains calcite derived from carbonation of a metamorphic amphibole—plagioclase assemblage that has further altered to siderite and dolomite. The calcite is porous and contains iron-rich calcite blebs interpreted to have resulted from fluid-mediated replacement of compositionally heterogeneous amphiboles. The siderite is polycrystalline but nucleates topotactically on the calcite. As a result, the boundaries between adjacent grains are low-angle boundaries (<10°), which are geometrically similar to those formed by crystal–plastic deformation and recovery. Growth zoning within individual siderite grains shows that the low-angle boundaries are growth features and not due to deformation. Low-angle boundaries develop due to the propagation of defects at grain faces and zone boundaries and by impingement of grains that nucleated with small misorientations relative to each other during grain growth. The cores of siderite grains are aligned with the twin planes in the parent calcite crystal showing that the reactant Fe entered the crystal along the twin boundaries. Dolomite grains, many of which appear to in-fill space generated by the siderite replacement, also show alignment of cores along the calcite twin planes, suggesting that they did not grow into space but replaced the calcite. Where dolomite is seen directly replacing calcite, it nucleates on the Fe-rich calcite due to the increased compatibility of the Fe-bearing calcite lattice relative to the pure calcite. Both reactions are interpreted as fluid-mediated replacement reactions which use the crystallography and elemental chemistry of the calcite. Experiments of fluid-mediated replacement reactions show that they proceed much faster than diffusion-based reactions. This is important when considering the rates of reactions relative to fluid flow in mineralising systems.     

Zircon M127 – A Homogeneous Reference Material for SIMS U–Pb Geochronology Combined with Hafnium,Oxygen and,Potentially, Lithium Isotope Analysis

In this article, we document a detailed analytical characterisation of zircon M127, a homogeneous 12.7 carat gemstone from Ratnapura, Sri Lanka. Zircon M127 has TIMS‐determined mean U–Pb radiogenic isotopic ratios of 0.084743 ± 0.000027 for ²⁰⁶Pb/²³⁸U and 0.67676 ± 0.00023 for ²⁰⁷Pb/²³⁵U (weighted means, 2s uncertainties). Its ²⁰⁶Pb/²³⁸U age of 524.36 ± 0.16 Ma (95% confidence uncertainty) is concordant within the uncertainties of decay constants. The δ¹⁸O value (determined by laser fluorination) is 8.26 ± 0.06‰ VSMOW (2s), and the mean ¹⁷⁶Hf/¹⁷⁷Hf ratio (determined by solution ICP‐MS) is 0.282396 ± 0.000004 (2s). The SIMS‐determined δ⁷Li value is ?0.6 ± 0.9‰ (2s), with a mean mass fraction of 1.0 ± 0.1 μg g^?1 Li (2s). Zircon M127 contains ~ 923 μg g^?1 U. The moderate degree of radiation damage corresponds well with the time‐integrated self‐irradiation dose of 1.82 × 10¹⁸ alpha events per gram. This observation, and the (U–Th)/He age of 426 ± 7 Ma (2s), which is typical of unheated Sri Lankan zircon, enable us to exclude any thermal treatment. Zircon M127 is proposed as a reference material for the determination of zircon U–Pb ages by means of SIMS in combination with hafnium and stable isotope (oxygen and potentially also lithium) determination.     

Timing of advance and basal condition of the Laurentide Ice Sheet during the last glacial maximum in the Richardson Mountains,NWT

This study presents new ages for the northwest section of the Laurentide Ice Sheet (LIS) glacial chronology from material recovered from two retrogressive thaw slumps exposed in the Richardson Mountains, Northwest Territories, Canada. One study site, located at the maximum glacial limit of the LIS in the Richardson Mountains, had calcite concretions recovered from aufeis buried by glacial till that were dated by U/Th disequilibrium to 18,500 cal yr BP. The second site, located on the Peel Plateau to the east yielded a fossil horse (Equus) mandible that was radiocarbon dated to ca. 19,700 cal yr BP. These ages indicate that the Peel Plateau on the eastern flanks of the Richardson Mountains was glaciated only after 18,500 cal yr BP, which is later than previous models for the global last glacial maximum (LGM). As the LIS retreated the Peel Plateau around 15,000 cal yr BP, following the age of the Tutsieta phase, we conclude that the presence of the northwestern margin of the LIS at its maximum limit was a very short event in the western Canadian Arctic.