Low-velocity layers within the crust can indicate the presence of melt and lithologic differences with implications for crustal composition and formation. Seismic wave conversions and reverberations across the base of the crust or intracrustal discontinuities, analysed using the receiver function method, can be used to constrain crustal layering. This is commonly accomplished by inverting receiver functions jointly with surface wave dispersion. Recently, the proliferation of model-space search approaches has made this technique a workhorse of crustal seismology. We show that reverberations from shallow layers such as sedimentary basins produce spurious low-velocity zones when inverted for crustal structure with surface wave data of insufficiently high frequency. Therefore, reports of such layers in the literature based on inversions using receiver function data should be re-evaluated. We demonstrate that a simple resonance-removal filter can suppress these effects and yield reliable estimates of crustal structure, and advocate for its use in receiver-function based inversions.

S U M M A R Y Inversion of surface wave data for crustal and upper-mantle structure is a staple of passive seismology, particularly since the advent of techniques enabling surface wave dispersion (SWD) and Rayleigh wave ellipticity measurements from ambient noise. Recent development and application of transdimensional Bayesian (TB) seismic inversion offers an approach to quantify model parameter uncertainties and trade-offs with fewer assumptions than traditional methods. Using synthetic tests, we investigate choices in the implementation of TB for the inversion of SWD and Rayleigh wave ellipticity to constrain the structure of Earth’s continental lithosphere. We focus on three aspects of the inversion: limitation of data sensitivity, assumed scaling among parameters (compressional wave speed, Vp, shear wave speed, Vs, density and radial anisotropy) and parametrization choices. We show that while surface wave data provide relatively strong constraints on the posterior distribution of Vs and, to a lesser extent, Vp, common parametrization choices can potentially bias structure estimates. This is particularly the case for radial anisotropy (ξ ), due to the inability to distinguish variations of Vp and density from those of ξ . Inferred results therefore depend substantially on the parametrization and scaling choices. We illustrate how layered parametrizations can, in the TB framework, recover smoothly varying profiles, and quantify the number of layers recoverable at different levels of measurement uncertainty. Finally, we propose two types of model parametrization for TB inversion involving multiple types of parameters. We demonstrate that by implementing an independent parametrization for different physical quantities, we can avoid imposing identical model geometry across multiple types of model parameters, and obtain better model estimates with reduced trade-offs. We advocate for such a parametrization in TB inversion of radial anisotropy using surface wave data, and when targeting disparate Vp and Vs structures such as those associated with α-β quartz transtion.

Seismic properties and equation-of-state parameters of the liquid iron alloy in the outer core are inferred from normal mode data. Turbulent convection of the liquid iron alloy outer core generates Earth’s magnetic field and supplies heat to the mantle. The exact composition of the iron alloy is fundamentally linked to the processes powering the convection and can be constrained by its seismic properties. Discrepancies between seismic models determined using body waves and normal modes show that these properties are not yet fully agreed upon. In addition, technical challenges in experimentally measuring the equation-of-state (EoS) parameters of liquid iron alloys at high pressures and temperatures further complicate compositional inferences. We directly infer EoS parameters describing Earth’s outer core from normal mode center frequency observations and present the resulting Elastic Parameters of the Outer Core (EPOC) seismic model. Unlike alternative seismic models, ours requires only three parameters and guarantees physically realistic behavior with increasing pressure for a well-mixed homogeneous material along an isentrope, consistent with the outer core’s condition. We show that EPOC predicts available normal mode frequencies better than the Preliminary Reference Earth Model (PREM) while also being more consistent with body wave–derived models, eliminating a long-standing discrepancy. The velocity at the top of the outer core is lower, and increases with depth more steeply, in EPOC than in PREM, while the density in EPOC is higher than that in PREM across the outer core. The steeper profiles and higher density imply that the outer core comprises a lighter but more compressible alloy than that inferred for PREM. Furthermore, EPOC’s steeper velocity gradient explains differential SmKS body wave travel times better than previous one-dimensional global models, without requiring an anomalously slow ~90- to 450-km-thick layer at the top of the outer core.

Abstract. Knowing the location of large-scale turbulent eddies during catastrophic flooding events improves predictions of erosive scour. The erosion damage to the Oroville Dam flood control spillway in early 2017 is an example of the erosive power of turbulent flow. During this event, a defect in the simple concrete channel quickly eroded into a 47 m deep chasm. Erosion by turbulent flow is difficult to evaluate in real time, but near-channel seismic monitoring provides a tool to evaluate flow dynamics from a safe distance. Previous studies have had limited ability to identify source location or the type of surface wave (i.e., Love or Rayleigh wave) excited by different river processes. Here we use a single three-component seismometer method (frequency-dependent polarization analysis) to characterize the dominant seismic source location and seismic surface waves produced by the Oroville Dam flood control spillway, using the abrupt change in spillway geometry as a natural experiment. We find that the scaling exponent between seismic power and release discharge is greater following damage to the spillway, suggesting additional sources of turbulent energy dissipation excite more seismic energy. The mean azimuth in the 5–10 Hz frequency band was used to resolve the location of spillway damage. Observed polarization attributes deviate from those expected for a Rayleigh wave, though numerical modeling indicates these deviations may be explained by propagation up the uneven hillside topography. Our results suggest frequency-dependent polarization analysis is a promising approach for locating areas of increased flow turbulence. This method could be applied to other erosion problems near engineered structures as well as to understanding energy dissipation, erosion, and channel morphology development in natural rivers, particularly at high discharges.