Daoyuan Sun

Prof. Daoyuan Sun

orcid

Research Focus & Metrics

Academic Career Length

11Years (2016-Present)

(Calculated once Scilit publication data is synced: current year − year of first publication + 1)

Representative Publications

View All Publications
7 January 2020 in Nature Communications
Verified
Journal Article
Open access
PDF file
WWWenzhong Wang
Wenzhong Wang
YXYinhan Xu
Yinhan Xu
DSDaoyuan Sun
Daoyuan Sun
SNSidao Ni
Sidao Ni
RWRenata Wentzcovitch
Renata Wentzcovitch
ZWZhongqing Wu
Zhongqing Wu
Seismic heterogeneities detected in the lower mantle were proposed to be related to subducted oceanic crust. However, the velocity and density of subducted oceanic crust at lower-mantle conditions remain unknown. Here, we report ab initio results for the elastic properties of calcium ferrite‐type phases and determine the velocities and density of oceanic crust along different mantle geotherms. We find that the subducted oceanic crust shows a large negative shear velocity anomaly at the phase boundary between stishovite and CaCl2-type silica, which is highly consistent with the feature of mid-mantle scatterers. After this phase transition in silica, subducted oceanic crust will be visible as high-velocity heterogeneities as imaged by seismic tomography. This study suggests that the presence of subducted oceanic crust could provide good explanations for some lower-mantle seismic heterogeneities with different length scales except large low shear velocity provinces (LLSVPs).

19 March 2019 in Geophysical Research Letters
Verified
Journal Article
Open access
PDF file
DSDaoyuan Sun
Daoyuan Sun
DHD. Helmberger
D. Helmberger
VLV. H. Lai
V. H. Lai
MGM. Gurnis
M. Gurnis
JJJennifer Jackson
Jennifer Jackson
HYH.‐Y. Yang
H.‐Y. Yang
At the core‐mantle boundary, most observed ultralow velocity zones (ULVZs) cluster along the edges of the large low shear velocity provinces (LLSVPs) and provide key information on the composition, dynamics, and evolution of the lower mantle. However, their detailed structure near slab‐like structures beneath the mid‐Pacific remains particularly challenging because of the lack of station coverage. While most studies of ULVZs concentrate on SKS‐complexity, here we report on the multipathing of ScS, which expands the sampling for ULVZs. We find the strongest multipathing along a ULVZ patch located just south of Hawaii and the far northeastern edge of the LLSVP, in a zone ~200 km in width and extending 600 km southward. The anomalous ScS travel times and distorted Sdiff waveforms further reveal patches interrupted by observed enhanced D″ indicative of slab‐debris influence on the complexity of the northeastern boundary of the mid‐Pacific LLSVP.

15 February 2022 in Geochemistry, Geophysics, Geosystems
Verified
Journal Article
Open access
PDF file
VLVoon Hui Lai
Voon Hui Lai
DHDon V. Helmberger
Don V. Helmberger
VDVasilije Dobrosavljevic
Vasilije Dobrosavljevic
WWWenbo Wu
Wenbo Wu
DSDaoyuan Sun
Daoyuan Sun
MGMichael Gurnis
Michael Gurnis
Strong waveform complexities, including multipathing of the S diffracted phase and rapid changes in differential ScS‐S times, are observed for multiple deep Fiji earthquakes recorded at the USArray. The complexities occur at the northeastern edge of the Pacific Large Low Shear Velocity Province (LLSVP), about 12 degrees southeast of present‐day Hawaiʻi. Waveform modeling of the multipathing provides good constraints on an ultra‐low velocity zone (ULVZ) with a width of 5 degree located near the inner edge of the LLSVP. Based on the mineralogical‐modeling of the ULVZ as a solid iron‐rich magnesiowüstite‐bearing assemblage with compatible morphology predicted from geodynamical simulations, a ULVZ model with a thickness of 30 km and a shear wave velocity reduction of 18% is preferred. The rapid change in differential ScS‐S travel time is best explained by having both the aforementioned ULVZ and an adjacent high velocity structure near the LLSVP. Furthermore, a low‐velocity plume‐like structure could potentially explain the observed S travel time delay independent of ScS. These seismic features are proposed to be a ULVZ driven toward the edge of the LLSVP while potentially pushed by a subducted slab. This configuration may trigger plume generation due to strong thermal instabilities and is in the same vicinity where mantle flow models place the present‐day Hawaiian plume source. Multiple ScS can potentially be used to verify vertical plume structure in tomographic models but the accuracy of upper mantle structure, which is a key reflection point, needs to be considered.

23 April 2016 in Journal of Geophysical Research: Solid Earth
Verified
Journal Article
Open access
PDF file
DSDaoyuan Sun
Daoyuan Sun
DHDon Helmberger
Don Helmberger
MMMeghan S. Miller
Meghan S. Miller
JJJennifer M. Jackson
Jennifer M. Jackson
D″ represents one of the most dramatic thermal and compositional layers within our planet. In particular, global tomographic models display relatively fast patches at the base of the mantle along the circum‐Pacific which are generally attributed to slab debris. Such distinct patches interact with the bridgmanite (Br) to post‐bridgmanite (PBr) phase boundary to generate particularly strong heterogeneity at their edges. Most seismic observations for the D″ come from the lower mantle S wave triplication (Scd). Here we exploit the USArray waveform data to examine one of these sharp transitions in structure beneath Alaska. From west to east beneath Alaska, we observed three different characteristics in D″: (1) the western region with a strong Scd, requiring a sharp δVs = 2.5% increase; (2) the middle region with no clear Scd phases, indicating a lack of D″ (or thin Br‐PBr layer); and (3) the eastern region with strong Scd phase, requiring a gradient increase in δVs. To explain such strong lateral variation in the velocity structure, chemical variations must be involved. We suggest that the western region represents relatively normal mantle. In contrast, the eastern region is influenced by a relic slab that has subducted down to the lowermost mantle. In the middle region, we infer an upwelling structure that disrupts the Br‐PBr phase boundary. Such an interpretation is based upon a distinct pattern of travel time delays, waveform distortions, and amplitude patterns that reveal a circular‐shaped anomaly about 5° across which can be modeled synthetically as a plume‐like structure rising about 400 km high with a shear velocity reduction of ~5%, similar to geodynamic modeling predictions of upwellings.

3 September 2025 in Nature
Verified
Journal Article
Open access
PDF file
HBHuixing Bi
Huixing Bi
DSDaoyuan Sun
Daoyuan Sun
NSNingyu Sun
Ningyu Sun
ZMZhu Mao
Zhu Mao
MDMingwei Dai
Mingwei Dai
DHDouglas Hemingway
Douglas Hemingway
For rocky planets, the presence of a solid inner core has notable implications on the composition and thermal evolution of the core and on the magnetic history of the planet 1–3 . On Mars, geophysical observations have confirmed that the core is at least partially liquid 4–7 , but it is unknown whether any part of the core is solid. Here we present an analysis of seismic data acquired by the InSight mission, demonstrating that Mars has a solid inner core. We identify two seismic phases, the deep core-transiting phase, PKKP, and the inner core boundary reflecting phase, PKiKP, indicative of the inner core. Our inversions constrain the radius of the Martian inner core to about 613 ± 67 km, with a compressional velocity jump of around 30% across the inner core boundary, supported by additional inner-core-related seismic phases. These properties imply a concentration of distinct light elements in the inner core, segregated from the outer core through core crystallization. This finding provides an anchor point for understanding the thermal and chemical state of Mars. Moreover, the relationship between inner core formation and the Martian magnetic field evolution could provide insights into dynamo generation across planetary bodies.

24 February 2022 in Nature Communications
Verified
Journal Article
Open access
PDF file
JLJiewen Li
Jiewen Li
DSDaoyuan Sun
Daoyuan Sun
DBDan J. Bower
Dan J. Bower
Ultra-low velocity zones (ULVZs) are localized small-scale patches with extreme physical properties at the core-mantle boundary that often gather at the margins of Large Low Velocity Provinces (LLVPs). Recent studies have discovered several mega-sized ULVZs with a lateral dimension of ~900 km. However, the detailed structures and physical properties of these ULVZs and their relationship to LLVP edges are not well constrained and their formation mechanisms are poorly understood. Here, we break the degeneracy between the size and velocity perturbation of a ULVZ using two orthogonal seismic ray paths, and thereby discover a mega-sized ULVZ at the northern edge of the Pacific LLVP. The ULVZ is almost double the size of a previously imaged ULVZ in this region, but with half of the shear velocity reduction. This mega-sized ULVZ has accumulated due to stable mantle flow converging at the LLVP edge driven by slab-debris in the lower mantle. Such flow also develops the subvertical north-tilting edge of the Pacific LLVP.