scholarly journals CCREM: New reference Earth model from the global coda-correlation wavefield

2021 ◽  
Author(s):  
Xiaolong Ma ◽  
Hrvoje Tkalčić
Keyword(s):  
Inner Core ◽  
Velocity Structure ◽  
Body Waves ◽  
Earth Model ◽  
Seismic Velocities ◽  
New Paradigm ◽  
Reference Models ◽  
First Order ◽  
Global Seismology ◽  
S Period

The existing Earth reference models have provided an excellent one-dimensional representation of Earth’s properties as a function of its radius and explained many seismic observations in a broad frequency band. However, some discrepancies still exist among these models near the first-order discontinuities (e.g., the core-mantle and the inner-core boundaries) due to different datasets and approaches. As a new paradigm in global seismology, the analysis of coda-correlation wavefield is fundamentally different from interpreting direct observations of seismic phases or free oscillations of the Earth. The correlation features exist in global correlograms due to the similarity of body waves reverberating through the Earth’s interior. As such, there is a great potential to utilize the information stored in the coda-correlation wavefield in constraining the Earth’s internal structure. Here, we deploy the global earthquake-coda correlation wavefield as an independent data source in the 15-50 s period interval to increase the Earth's radial structure constraints. We assemble a dataset of multiple pronounced correlation features and fit both their travel times and waveforms by computing synthetic correlograms through a series of candidate models. Misfit measurements for correlation features are then computed to search for the best-fitting model. The model that provides an optimal representation of the correlation features in the coda-correlation wavefield is CCREM. It displays differences in radial seismic velocities, especially near the first-order discontinuities, relative to previously proposed Earth-reference models. This is the first application of the earthquake-coda correlation wavefield in constraining the whole Earth's radial velocity structure.

scholarly journals Methods for computing internal flattening, with applications to the Earth's structure and geodynamics

1998 ◽  
Vol 132 (3) ◽  
pp. 603-642 ◽  
Author(s):  
C. Denis ◽  
M. Amalvict ◽  
Y. Rogister ◽  
S. Tomecka-Suchoń
Keyword(s):  
Field Theory ◽  
Differential Rotation ◽  
Inner Core ◽  
Static Solution ◽  
Internal Field ◽  
Earth Model ◽  
Earth Structure ◽  
First Order ◽  
The Earth ◽  
Earth Models

SUMMARY After general comments (Section 1) on using variational procedures to compute the oblateness of internal strata in the Earth and slowly rotating planets, we recall briefly some basic concepts about barotropic equilibrium figures (Section 2), and then proceed to discuss several accurate methods to derive the internal flattening. The algorithms given in Section 3 are based on the internal gravity field theory of Clairaut, Laplace and Lyapunov. They make explicit use of the concept of a level surface. The general formulation given here leads to a number of formulae which are of both theoretical and practical use in studying the Earth's structure, dynamics and rotational evolution. We provide exact solutions for the figure functions of three Earth models, and apply the formalism to yield curves for the internal flattening as a function of the spin frequency. Two more methods, which use the general deformation equations, are discussed in Section 4. The latter do not rely explicitly on the existence of level surfaces. They offer an alternative to the classical first-order internal field theory, and can actually be used to compute changes of the flattening on short timescales produced by variations in the LOD. For short durations, the Earth behaves elastically rather than hydrostatically. We discuss in some detail static deformations and Longman's static core paradox (Section 5), and demonstrate that in general no static solution exists for a realistic Earth model. In Section 6 we deal briefly with differential rotation occurring in cylindrical shells, and show why differential rotation of the inner core such as has been advocated recently is incompatible with the concept of level surfaces. In Section 7 we discuss first-order hydrostatic theory in relation to Earth structure, and show how to derive a consistent reference Earth model which is more suitable for geodynamical modelling than are modern Earth models such as 1066-A, PREM or CORE11. An important result is that a consistent application of hydrostatic theory leads to an inertia factor of about 0.332 instead of the value 0.3308 used until now. This change automatically brings ‘hydrostatic’ values of the flattening, the dynamic shape factor and the precessional constant into much better agreement with their observed counterparts than has been assumed hitherto. Of course, we do not imply that non-hydrostatic effects are unimportant in modelling geodynamic processes. Finally, we discuss (Sections 7–8) some implications of our way of looking at things for Earth structure and some current problems of geodynamics. We suggest very significant changes for the structure of the core, in particular a strong reduction of the density jump at the inner core boundary. The theoretical value of the free core nutation period, which may be computed by means of our hydrostatic Earth models CGGM or PREMM, is in somewhat better agreement with the observed value than that based on PREM or 1066-A, although a significant residue remains. We attribute the latter to inadequate modelling of the deformation, and hence of the change in the inertia tensor, because the static deformation equations were used. We argue that non-hydrostatic effects, though present, cannot explain the large observed discrepancy of about 30 days.

scholarly journals Thermal Equation of State of Fe3C to 327 GPa and Carbon in the Core

Minerals ◽  
2019 ◽  
Vol 9 (12) ◽  
pp. 744 ◽  
Author(s):  
Suguru Takahashi ◽  
Eiji Ohtani ◽  
Daijo Ikuta ◽  
Seiji Kamada ◽  
Tatsuya Sakamaki ◽  
...  
Keyword(s):  
High Temperature ◽  
Equation Of State ◽  
Inner Core ◽  
Velocity Structure ◽  
Light Element ◽  
Scale Model ◽  
Earth Model ◽  
Thermal Equation ◽  
The Core ◽  
Pressure Scale

The density and sound velocity structure of the Earth’s interior is modeled on seismological observations and is known as the preliminary reference Earth model (PREM). The density of the core is lower than that of pure Fe, which suggests that the Earth’s core contains light elements. Carbon is one plausible light element that may exist in the core. We determined the equation of state (EOS) of Fe3C based on in situ high-pressure and high-temperature X-ray diffraction experiments using a diamond anvil cell. We obtained the P–V data of Fe3C up to 327 GPa at 300 K and 70–180 GPa up to around 2300 K. The EOS of nonmagnetic (NM) Fe3C was expressed by two models using two different pressure scales and the third-order Birch–Murnaghan EOS at 300 K with the Mie–Grüneisen–Debye EOS under high-temperature conditions. The EOS can be expressed with parameters of V0 = 148.8(±1.0) Å3, K0 = 311.1(±17.1) GPa, K0′ = 3.40(±0.1), γ0 = 1.06(±0.42), and q = 1.92(±1.73), with a fixed value of θ0 = 314 K using the KBr pressure scale (Model 1), and V0 = 147.3(±1.0) Å3, K0 = 323.0(±16.6) GPa, K0′ = 3.43(±0.09), γ0 = 1.37(±0.33), and q = 0.98(±1.01), with a fixed value of θ0 = 314 K using the MgO pressure scale (Model 2). The density of Fe3C under inner core conditions (assuming P = 329 GPa and T = 5000 K) calculated from the EOS is compatible with the PREM inner core.

scholarly journals Figures of merit for stellarators near the magnetic axis

2021 ◽  
Vol 87 (1) ◽  
Author(s):  
Matt Landreman
Keyword(s):  
Magnetic Field ◽  
Field Strength ◽  
Magnetic Axis ◽  
Configuration Design ◽  
New Paradigm ◽  
First Order ◽  
Figures Of Merit ◽  
Minor Radius ◽  
Error Field ◽  
Quantities Of Interest

A new paradigm for rapid stellarator configuration design has been recently demonstrated, in which the shapes of quasisymmetric or omnigenous flux surfaces are computed directly using an expansion in small distance from the magnetic axis. To further develop this approach, here we derive several other quantities of interest that can be rapidly computed from this near-axis expansion. First, the $\boldsymbol {\nabla }\boldsymbol {B}$ and $\boldsymbol {\nabla }\boldsymbol {\nabla }\boldsymbol {B}$ tensors are computed, which can be used for direct derivative-based optimization of electromagnetic coil shapes to achieve the desired magnetic configuration. Moreover, if the norm of these tensors is large compared with the field strength for a given magnetic field, the field must have a short length scale, suggesting it may be hard to produce with coils that are suitably far away. Second, we evaluate the minor radius at which the flux surface shapes would become singular, providing a lower bound on the achievable aspect ratio. This bound is also shown to be related to an equilibrium beta limit. Finally, for configurations that are constructed to achieve a desired magnetic field strength to first order in the expansion, we compute the error field that arises due to second-order terms.

Seismic cone penetration test and seismic tomography in permafrost

2004 ◽  
Vol 41 (5) ◽  
pp. 796-813 ◽  
Author(s):  
Anne-Marie LeBlanc ◽  
Richard Fortier ◽  
Michel Allard ◽  
Calin Cosma ◽  
Sylvie Buteau
Keyword(s):  
Seismic Tomography ◽  
Dynamic Properties ◽  
Cone Penetration Test ◽  
The Body ◽  
Body Waves ◽  
Reconstruction Technique ◽  
Seismic Velocities ◽  
Penetration Test ◽  
Cone Penetration ◽  
Seismic Cone Penetration Test

Two high-resolution multi-offset vertical seismic profile (VSP) surveys were carried out in a permafrost mound near Umiujaq in northern Quebec, Canada, while performing seismic cone penetration tests (SCPT) to study the cryostratigraphy and assess the body waves velocities and the dynamic properties of warm permafrost. Penetrometer-mounted triaxial accelerometers were used as the VSP receivers, and a swept impact seismic technique (SIST) source generating both compressional and shear waves was moved near the surface following a cross configuration of 40 seismic shot-point locations surrounding each of the two SCPTs. The inversion of travel times based on a simultaneous iterative reconstruction technique (SIRT) provided tomographic images of the distribution of seismic velocities in permafrost. The Young's and shear moduli at low strains were then calculated from the seismic velocities and the permafrost density measured on core samples. The combination of multi-offset VSP survey, SCPT, SIST, and SIRT for tomographic imaging led to new insights in the dynamic properties of permafrost at temperatures close to 0 °C. The P- and S-wave velocities in permafrost vary from 2400 to 3200 m/s and from 900 to 1750 m/s, respectively, for a temperature range between –0.2 and –2.0 °C. The Young's modulus varies from 2.15 to 13.65 GPa, and the shear modulus varies from 1.00 to 4.75 GPa over the same range of temperature.Key words: permafrost, seismic cone penetration test, vertical seismic profiling, seismic tomography, dynamic properties.

Quaternary sodic and potassic intraplate volcanism in northeast China controlled by the underlying heterogeneous lithospheric structures

Geology ◽  
2021 ◽  
Author(s):  
Xingli Fan ◽  
Qi-Fu Chen ◽  
Yinshuang Ai ◽  
Ling Chen ◽  
Mingming Jiang ◽  
...  
Keyword(s):  
Upper Mantle ◽  
Lithospheric Mantle ◽  
Northeast China ◽  
Velocity Structure ◽  
Seismic Velocity ◽  
Subcontinental Lithospheric Mantle ◽  
Seismic Velocities ◽  
S Wave ◽  
Crust And Upper Mantle ◽  
Potassic Volcanism

The origin and mantle dynamics of the Quaternary intraplate sodic and potassic volcanism in northeast China have long been intensely debated. We present a high-resolution, three-dimensional (3-D) crust and upper-mantle S-wave velocity (Vs) model of northeast China by combining ambient noise and earthquake two-plane wave tomography based on unprecedented regional dense seismic arrays. Our seismic images highlight a strong correlation between the basalt geochemistry and upper-mantle seismic velocity structure: Sodic volcanoes are all characterized by prominent low seismic velocities in the uppermost mantle, while potassic volcanoes still possess a normal but thin upper-mantle “lid” depicted by high seismic velocities. Combined with previous petrological and geochemical research findings, we propose that the rarely erupted Quaternary potassic volcanism in northeast China results from the interaction between asthenospheric low-degree melts and the overlying subcontinental lithospheric mantle. In contrast, the more widespread Quaternary sodic volcanism in this region is predominantly sourced from the upwelling asthenosphere without significant overprinting from the subcontinental lithospheric mantle.

scholarly journals New measurements of long-period radial modes using large earthquakes

2020 ◽  
Vol 224 (2) ◽  
pp. 1211-1224
Author(s):  
S Talavera-Soza ◽  
A Deuss
Keyword(s):  
Inner Core ◽  
Cross Coupling ◽  
Earth Model ◽  
Cylindrical Anisotropy ◽  
Long Period ◽  
Centre Frequency ◽  
The Earth ◽  
Radial Anisotropy ◽  
Crucial Information ◽  
Quality Factor Q

SUMMARY Radial modes, nS0, are long-period oscillations that describe the radial expansion and contraction of the whole Earth. They are characterized only by their centre frequency and quality factor Q, and provide crucial information about the 1-D structure of the Earth. Radial modes were last measured more than a decade ago using only one or two earthquakes. Here, we measure radial modes using 16 of the strongest and deepest earthquakes of the last two decades. By introducing more earthquake data into our measurements, we improve our knowledge of 1-D attenuation, as we remove potential earthquake bias from our results. For mode 0S0, which is dominated by compressional energy, we measure a Q value of 5982, much higher than previously measured, and requiring less bulk attenuation in the Earth than previously thought. We also show that radial modes cross-couple (resonate) strongly to their nearest spheroidal mode due to ellipticity and inner core cylindrical anisotropy. Cross-coupling improves the fit between data and synthetics, and gives better estimates of the centre frequency and attenuation value of the radial modes. Including cross-coupling in our measurements results in a systematic shift of the centre frequencies of radial modes towards the Preliminary Reference Earth Model. This shift in centre frequencies, has implications for the strength of the radial anisotropy present in the uppermost inner core, with our cross-coupling results agreeing with lower values of anisotropy than the ones inferred from just measuring the modes in self-coupling (isolation). Furthermore, cross-coupling between radial modes and angular-order two modes provides constraints on cylindrical inner core anisotropy, that will help us improve our knowledge of the 3-D structure of the inner core.

scholarly journals Constraining 1-D inner core attenuation through measurements of strongly coupled normal mode pairs

2020 ◽  
Vol 223 (1) ◽  
pp. 612-621 ◽  
Author(s):  
S Talavera-Soza ◽  
A Deuss
Keyword(s):  
Normal Mode ◽  
Inner Core ◽  
Strong Dependence ◽  
Normal Modes ◽  
Cross Coupling ◽  
Body Waves ◽  
Spectral Amplitude ◽  
Cylindrical Anisotropy ◽  
Strongly Coupled ◽  
Mode Pair

SUMMARY We measured inner core normal mode pair 10S2–11S2, which cross-couples strongly for 1-D structure and is sensitive to shear wave velocity, and find that our measurements agree with a strongly attenuating inner core. In the past, this mode pair has been used to try to resolve the debate on whether the inner core is strongly or weakly attenuating. Its large spectral amplitude in observed data, possible through the apparent low attenuation of 10S2, has been explained as evidence of a weakly attenuating inner core. However, this contradicted body waves and other normal modes studies, which resulted in this pair of modes being excluded from inner core modelling. Modes 10S2 and 11S2 are difficult to measure and interpret because they depend strongly on the underlying 1-D model used. This strong dependence makes these modes change both their oscillation characteristics and attenuation values under a small 1-D perturbation to the inner core model. Here, we include this effect by allowing the pair of modes to cross-couple or resonate through 1-D structure and treat them as one hybrid mode. We find that, unlike previously thought, the source of 10S2 visibility is its strong cross-coupling to 11S2 for both 1-D elastic and anelastic structure. We also observe that the required 1-D perturbation is much smaller than the 2 per cent vs perturbation previously suggested, because we simultaneously measure 3-D structure in addition to 1-D structure. Only a 0.5 per cent increase in inner core vs or a 0.5 per cent decrease in inner core radius is required to explain 10S2–11S2 observations and a weakly attenuating inner core is not needed. In addition, the 3-D structure measurements of mode 10S2 and its cross-coupling to 11S2 show the typical strong zonal splitting pattern attributed to inner core cylindrical anisotropy, allowing us to add further constrains to deeper regions of the inner core.

Crustal velocity structure of the southern Nechako basin, British Columbia, from wide-angle seismic traveltime inversion1This article is one of a series of papers published in this Special Issue on the theme of New insights in Cordilleran Intermontane geoscience: reducing exploration risk in the mountain pine beetle-affected area, British Columbia.

2011 ◽  
Vol 48 (6) ◽  
pp. 1050-1063 ◽  
Author(s):  
A.L. Stephenson ◽  
G.D. Spence ◽  
K. Wang ◽  
J.A. Hole ◽  
K.C. Miller ◽  
...  
Keyword(s):  
British Columbia ◽  
Mountain Pine Beetle ◽  
Velocity Structure ◽  
Volcanic Rocks ◽  
Sedimentary Rocks ◽  
Velocity Model ◽  
Seismic Profile ◽  
Wide Angle ◽  
Seismic Velocities ◽  
Coast Mountains

In the BATHOLITHSonland seismic project, a refraction – wide-angle reflection survey was shot in 2009 across the Coast Mountains and Interior Plateau of central British Columbia. Part of the seismic profile crossed the Nechako Basin, a Jurassic–Cretaceous basin with potential for hydrocarbons within sedimentary strata that underlies widespread volcanic rocks. Along this 205 km-long line segment, eight large explosive shots were fired into 980 seismometers. Forward and inverse modelling of the traveltime data were conducted with two independent methods: ray-tracing based modelling of first and secondary arrivals, and a higher resolution wavefront-based first-arrival seismic tomography. Material with velocities less than 5.0 km/s is interpreted as sedimentary rocks of the Nechako Basin, while velocities from 5.0–6.0 km/s may correspond to interlayered sedimentary and volcanic rocks. The greatest thickness of sedimentary rocks in the basin is found in the central 110 km of the profile. Two sub-basins were identified in this region, with widths of 20–50 km and maximum sedimentary depths of 2.5 and 3.3 km. Such features are well-defined in the velocity model, since resolution tests indicate that features with widths greater than ∼13 km are reliable. Beneath the sedimentary rocks, seismic velocities increase more slowly with depth — from 6.0 km/s just below the basin to 6.3 km/s at ∼17 km in depth, and then to 6.8–7.0 km/s at the base of the crust. The Moho is found at a depth of 33.5–35 km beneath the profile, and mantle velocities are high at 8.05–8.10 km/s.

Sign in / Sign up

Export Citation Format

Share Document