An official website of the United States government
SUMMARY Improving the resolution of seismic anelastic models is critical for a better understanding of the Earth’s subsurface structure and dynamics. Seismic attenuation plays a crucial role in estimating water content, partial melting and temperature variations in the Earth’s crust and mantle. However, compared to seismic wave-speed models, seismic attenuation tomography models tend to be less resolved. This is due to the complexity of amplitude measurements and the challenge of isolating the effect of attenuation in the data from other parameters. Physical dispersion caused by attenuation also affects seismic wave speeds, and neglecting scattering/defocusing effects in classical anelastic models can lead to biased results. To overcome these challenges, it is essential to account for the full 3-D complexity of seismic wave propagation. Although various synthetic tests have been conducted to validate anelastic full-waveform inversion (FWI), there is still a lack of understanding regarding the trade-off between elastic and anelastic parameters, as well as the variable influence of different parameter classes on the data. In this context, we present a synthetic study to explore different strategies for global anelastic inversions. To assess the resolution and sensitivity for different misfit functions, we first perform mono-parameter inversions by inverting only for attenuation. Then, to study trade-offs between parameters and resolution, we test two different inversion strategies (simultaneous and sequential) to jointly constrain the elastic and anelastic parameters. We found that a sequential inversion strategy performs better for imaging attenuation than a simultaneous inversion. We also demonstrate the dominance of seismic wave speeds over attenuation, underscoring the importance of determining a good approximation of the Hessian matrix and suitable damping factors for each parameter class.
more »
« less
Teleseismic P ‐Wave Attenuation Beneath the Southeastern United States
Abstract Seismic attenuation is an important parameter for characterizing subsurface morphology and thermal structure. In this study, we useP‐wave amplitude spectra from 588 teleseismic events recorded by 477 broadband seismic stations in the southeastern United States to examine the spatial variations of seismic attenuation in the crust and upper mantle. The resulting seismic attenuation parameter (∆t*) measurements obtained using the spectral ratio technique reveal a zone of relatively low attenuation in the Gulf of Mexico Coastal Plain and the southwestern terminus of the Piedmont province. Spatial coherency analysis of the ∆t* observations suggests that the center of the low attenuation layer is located within the uppermost mantle at about 70 km depth. This low attenuation anomaly lies along the suture zone between Laurentia and Gondwana and approximately coincides with the east‐west trending Brunswick magnetic anomaly. The origin of this low attenuation anomaly can be attributed to low attenuation bodies in the form of remnant lithospheric fragments in the deep crust and the uppermost mantle. The contribution of scattering to the observed ∆t* is estimated by calculating the ratio of amplitude on the transverse and vertical components in theP‐wave window. Relative to the rest of the study area, the Gulf of Mexico Coastal Plain demonstrates weaker scattering which is suggestive of a more homogenous crustal and uppermost mantle structure.
more »
« less
- Award ID(s):
- 1919789
- PAR ID:
- 10449192
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Geochemistry, Geophysics, Geosystems
- Volume:
- 22
- Issue:
- 6
- ISSN:
- 1525-2027
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
SUMMARY We present models of crustal and uppermost mantle structure beneath the Hawaiian Swell and surrounding region. The models were derived from ambient-noise intermediate-period Rayleigh-wave phase velocities and from seafloor compliance that were estimated from continuous seismic and pressure recordings collected during the Hawaiian Plume-Lithosphere Undersea Mantle Experiment (PLUME). We jointly inverted these data at the locations of over 50 ocean-bottom instruments, after accounting for variations in local bathymetry and sediment properties. Our results suggest that the crystalline crust is up to 15 km thick beneath the swell and up to 23 km thick closer to the islands. Anomalously thick crust extends towards the older seamounts, downstream of Hawaii. In a second region, anomalies immediately to the south of Hawaii may be associated with the leading edge of the shallow Hawaiian magma conduit. In a third region, thickened crust to the immediate west of Hawaii may be related to Cretaceous seamounts. Low seismic velocities identified in the uppermost mantle to the northeast of Hawaii may be linked to the Molokai fracture zone and may be manifest of complex non-vertical pathways of melt through the upper lithosphere. Velocity anomalies decrease in amplitude towards the surface, suggesting that melt becomes focused into conduits at depths between 20 and 40 km that escape the resolution capabilities of our data set.more » « less
-
Abstract We image seismic attenuation near the Hikurangi trench offshore New Zealand, using ocean bottom and land‐based seismometers, revealing high attenuation above a recurring shallow slow‐slip event and within the subducting Hikurangi Plateau. The Hikurangi subduction margin east of the North Island, New Zealand is the site of frequent shallow slow slip events. Overpressured fluids are hypothesized to lead to slow slip at shallow depths close to the oceanic trench. Seismic attenuation, energy loss of seismic waves, can be used to detect high temperatures, melt, the presence of fluids, and fractures. We use local earthquake P‐ and S‐waves from 180 earthquakes to invert fort*, and subsequently invert for Qp and Qs, offshore the North Island directly above the area of slow slip. We image Qp and Qs to ∼25 km depth, increasing resolution of previously identified coastal lowQ(100–300), and finding a new region of even higher attenuation (Qp and Qs < 50–100) directly above the shallow slow slip event of 2014–2015, beneath the offshore seismic array. This highest attenuation is downdip of a subducting seamount, and is spatially correlated with a high seismic reflectivity zone and Vp/Vs > 1.85, all of which provide evidence for the presence of fluids. The Qp and Qs is low at the trench (<50–100) and in the subducting plate (100–200), suggesting that seismic wave scattering due to faults, fractures, and the inherent heterogeneous composition of the Hikurangi Plateau, a large igneous province, plays a role in seismic attenuation.more » « less
-
Abstract The eastern North American passive margin was modified by Mesozoic rifting. Seismic data from recent deployment of onshore and offshore stations offer a unique opportunity for studying the signature of syn‐rifting and postrifting in lithospheric structures. Using full‐wave ambient noise tomography, we construct a new seismic velocity model for the lithosphere of the southeastern United States. Our model confirms an oceanic‐continental transitional crust over a ∼70 km wide zone across the coastline. Our model reveals (a) a patch of lower‐than‐average mantle lithospheric velocities underlying this transitional crust and (b) a low‐velocity column in the mantle lithosphere beneath the Virginia volcanoes. We propose that anomaly 1 represents cooled enriched mantle that underplated the thinning crust during the initial stages of rifting around 230 Ma. Anomaly 2 likely has a more recent origin in the Eocene and may result from an asthenospheric upwelling induced by a localized lithospheric delamination.more » « less
-
null (Ed.)SUMMARY Interfaces are important part of Earth’s layering structure. Here, we developed a new model parametrization and iterative linearized inversion method that determines 1-D crustal velocity structure using surface wave dispersion, teleseismic P-wave receiver functions and Ps and PmP traveltimes. Unlike previous joint inversion methods, the new model parametrization includes interface depths and layer Vp/Vs ratios so that smoothness constraint can be conveniently applied to velocities of individual layers without affecting the velocity discontinuity across the interfaces. It also allows adding interface-related observation such as traveltimes of Ps and PmP in the joint inversion to eliminate the trade-off between interface depth and Vp/Vs ratio and therefore to reduce the uncertainties of results. Numerical tests show that the method is computationally efficient and the inversion results are robust and independent of the initial model. Application of the method to a dense linear array across the Wabash Valley Seismic Zone (WVSZ) produced a high-resolution crustal image in this seismically active region. The results show a 51–55-km-thick crust with a mid-crustal interface at 14–17 km. The crustal Vp/Vs ratio varies from 1.69 to 1.90. There are three pillow-like, ∼100 km apart high-velocity bodies sitting at the base of the crust and directly above each of them are a low-velocity anomaly in the middle crust and a high-velocity anomaly in the upper crust. They are interpreted to be produced by mantle magmatic intrusions and remelting during rifting events in the end of the Precambrian. The current diffuse seismicity in the WVSZ might be rooted in this ancient distributed rifting structure.more » « less
