An official website of the United States government
Dynamic topography refers to vertical deflections of Earth’s surface from viscous flow within the mantle. Here we investigate how past subduction history affects present dynamic topography. We assimilate two plate reconstructions into TERRA forward mantle convection models to calculate past mantle states and predict Earth’s present dynamic topography; a comparison is made with a database of observed oceanic residual topography. The two assimilated plate reconstructions ‘Earthbyte’ and ‘Tomopac’ show divergent subduction histories across an extensive deep-time interval within Pacific-Panthalassa. We find that introducing an alternative subduction history perturbs our modelled present-day dynamic topography on the same order as the choice of radial viscosity. Additional circum-Pacific intra-oceanic subduction in Tomopac consistently produces higher correlations to the geoid (more than 20% improvement). At spherical harmonic degrees 1–40, dynamic topography models with intra-oceanic subduction produce universally higher correlations with observations and improve fit by up to 37%. In northeast Asia, Tomopac models show higher correlations (0.46 versus 0.18) to observed residual topography and more accurately predict approximately 1 km of dynamic subsidence within the Philippine Sea plate. We demonstrate that regional deep-time changes in subduction history have widespread impacts on the spatial distribution and magnitude of present-day dynamic topography. Specifically, we find that local changes to plate motion histories can induce dynamic topography changes in faraway regions located thousands of kilometres away. Our results affirm that present-day residual topography observations provide a powerful, additional constraint for reconstructing ancient subduction histories.
more »
« less
NW Pacific‐Panthalassa Intra‐Oceanic Subduction During Mesozoic Times From Mantle Convection and Geoid Models
Abstract Pacific‐Panthalassa plate tectonics are the most challenging on Earth to reconstruct during the Mesozoic and Cenozoic eras due to extensive subduction, which has resulted in large (>9,000 km length) unconstrained gaps between the Pacific and Laurasia (now NE Asia) back to the Early Jurassic. We build four contrasted NW Pacific‐Panthalassa global plate reconstructions and assimilate their velocity fields into global geodynamic models. We compare our predicted present mantle structure, synthetic geoid and dynamic topography to Earth observations. P‐wave tomographic filtering of predicted mantle structures allows for more explicit comparisons to global tomography. Plate reconstructions that include intra‐oceanic subduction in NW Pacific‐Panthalassa fit better to the observed geoid and residual topography, challenging popular models of Andean‐style subduction along East Asia. Our geodynamic models predict significant SE‐ward lateral slab advections within the NW Pacific basin lower mantle (∼2,500 km from Mesozoic times to present) that would confound “vertical slab sinking”‐style restorations of imaged slabs and past subduction zone locations.
more »
« less
- Award ID(s):
- 1848327
- PAR ID:
- 10382369
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Geochemistry, Geophysics, Geosystems
- Volume:
- 23
- Issue:
- 11
- ISSN:
- 1525-2027
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract The plate tectonic history of the hypothesized “proto‐South China Sea” (PSCS) ocean basin and surrounding SE Asia since Cenozoic times is controversial. We implement four diverse proto‐South China Sea plate reconstructions into global geodynamic models to constrain PSCS plate tectonics and possible slab locations. Our plate reconstructions consider the following: southward versus double‐sided PSCS subduction models; earlier (Eocene) or later (late Oligocene) initiation of Borneo counterclockwise rotations; and larger or smaller reconstructed Philippine Sea plate sizes. We compare our modeling results against tomographic images by accounting for mineralogical effects and the finite resolution of seismic tomography. All geodynamic models reproduce the tomographically imaged Sunda slabs beneath Peninsular Malaysia, Sumatra, and Java. Southward PSCS subduction produces slabs beneath present Palawan, northern Borneo, and offshore Palawan. Double‐sided PSCS subduction combined with earlier Borneo rotations uniquely reproduces subhorizontal slabs under the southern South China Sea (SCS) at ~400 to 700 km depths; these models best fit seismic tomography. A smaller Philippine Sea (PS) plate with a ~1,000‐km‐long restored Ryukyu slab was superior to a very large PS plate. Considered together, our four end‐member plate reconstructions predict that the PSCS slabs are now at <900 km depths under present‐day Borneo, the SCS, the Sulu and Celebes seas, and the southern Philippines. Regardless of plate reconstruction, we predict (1) mid‐Cenozoic passive return‐flow upwellings under Indochina; and (2) late Cenozoic downwellings under the SCS that do not support a deep‐origin “Hainan plume.” Modeled Sundaland dynamic topography strongly depends on the imposed plate reconstructions, varying by almost 1 km.more » « less
-
Abstract Simulating present‐day solid Earth deformation and volatile cycling requires integrating diverse geophysical data sets and advanced numerical techniques to model complex multiphysics processes at high resolutions. Subduction zone modeling is particularly challenging due to the large geographic extent, localized deformation zones, and the strong feedbacks between reactive fluid transport and solid deformation. Here, we develop new workflows for simulating 3‐dimensional thermal‐mechanical subduction and patterns of volatile dehydration at convergent margins, adaptable to include reactive fluid transport. We apply these workflows to the Hikurangi margin, where recent geophysical investigations have offered unprecedented insight into the structure and processes coupling fluid transport and solid deformation across broad spatiotemporal scales. Geophysical data sets constraining the downgoing and overriding plate structure are collated with the Geodynamic World Builder, which provides the initial conditions for forward simulations using the open‐source geodynamic modeling software code ASPECT. We systematically examine how plate interface rheology and hydration of the downgoing plate and upper mantle influence Pacific–Australian convergence and seismic anisotropy. Models prescribing a plate boundary viscosity of Pa s best reproduce observed plate velocities, and changing the configuration of the Pacific–Australia plate boundary directly influences the modeled plate motions. Models considering hydrated olivine fabrics best reproduce observations of seismic anisotropy. Predicted patterns of slab dehydration and mantle melting correlate well with observations of seismic attenuation and arc volcanism. These results suggest that hydration‐related rheological heterogeneity and related fluid weakening may strongly influence slab dynamics. Future investigations integrating coupled fluid transport and global mantle flow will provide insight into the feedbacks between subduction dynamics, fluid pathways, and arc volcanism.more » « less
-
null (Ed.)The configuration of mid-ocean ridges subducted below North America prior to Oligocene time is unconstrained by seafloor isochrons and has been primarily inferred from upper-plate geology, including near-trench magmatism. However, many tectonic models are permitted from these constraints. We present a fully kinematic, plate tectonic reconstruction of the NW Cordillera since 60 Ma built by structurally unfolding subducted slabs, imaged by mantle tomography, back to Earth’s surface. We map in three-dimensions the attached Alaska and Cascadia slabs, and a detached slab below western Yukon (Canada) at 400−600 km depth that we call the “Yukon Slab.” Our restoration of these lower plates within a global plate model indicates the Alaska slab accounts for Pacific-Kula subduction since ca. 60 Ma below the Aleutian Islands whereas the Cascadia slab accounts for Farallon subduction since at least ca. 75 Ma below southern California, USA. However, intermediate areas show two reconstruction gaps that persist until 40 Ma. We show that these reconstruction gaps correlate spatiotemporally to published NW Cordillera near-trench magmatism, even considering possible terrane translation. We attribute these gaps to thermal erosion related to ridge subduction and model mid-ocean ridges within these reconstruction gap mid-points. Our reconstructions show two coeval ridge-trench intersections that bound an additional “Resurrection”-like plate along the NW Cordillera prior to 40 Ma. In this model, the Yukon slab represents a thermally eroded remnant of the Resurrection plate. Our reconstructions support a “northern option” Farallon ridge geometry and allow up to ∼1200 km Chugach terrane translation since Paleocene time, providing a new “tomographic piercing point” for the Baja-British Columbia debate.more » « less
-
Abstract The Earth's long‐ and intermediate‐wavelength geoid anomalies are surface expressions of mantle convection and are sensitive to mantle viscosity. While previous studies of the geoid provide important constraints on the mantle radial viscosity variations, the mantle buoyancy in these studies, as derived from either seismic tomography or slab density models, may suffer significant uncertainties. In this study, we formulate 3‐D spherical mantle convection models with plate motion history since the Cretaceous that generate dynamically self‐consistent mantle thermal and buoyancy structures, and for the first time, use the dynamically generated slab structures and the observed geoid to place important constraints on the mantle viscosity. We found that non‐uniform weak plate margins and strong plate interiors are critical in reproducing the observed geoid and surface plate motion, especially the net lithosphere rotation (i.e., degree‐1 toroidal plate motion). In the best‐fit model, which leads to correlation of 0.61 between the modeled and observed geoid at degrees 4–12, the lower mantle viscosity is ∼1.3–2.5 × 1022 Pa⋅s and is ∼30 and ∼600–1,000 times higher than that in the transition zone and asthenosphere, respectively. Slab structures and the geoid are also strongly affected by slab strength, and the observations prefer moderately strong slabs that are ∼10–100 times stronger than the ambient mantle. Finally, a thin weak layer below the 670‐km phase change on a regional scale only in subduction zones produces stagnant slabs in the mantle transition zone as effectively as a weak layer on a global scale.more » « less
