Search

The Colorado basin in the SW Atlantic is one of the key exploratory frontier basins along the Argentine shelf. The stratigraphy shows a typical rift basin, divided into two major seismic successions, a Jurassic-Lower Cretaceous synrift followed by an Upper Cretaceous to Cenozoic postrift. They represent, respectively, the breakup of Gondwana and passive margin evolution. While the synrift basin formation was associated with the development of structural grabens, the preservation of thick (thousands of meters) postrift sequences are still debated because they cannot be accounted for only by thermal subsidence. In this work, we conducted a multimethodological subsidence analysis across the Colorado basin, including part of the continental slope. We performed a 1D backstripping analysis on 10 hydrocarbon exploration wells and 32 pseudo-wells, covering the whole basin area and a uniform stretching model, which were supported using dynamic topography data from the study region. From the backstripping analyses, considering the six major postrift sedimentary successions described in the literature, we could recognize three subsidence stages. An initial rapid subsidence stage during the post-rift (125–83 Ma), which is restricted to the main depocenters, followed by a stage of generalized and low subsidence focused on the main axis of the basin (83–33 Ma), with an increase in the last 33 my (33–0 Ma). The stretching models exhibit a typical asymptotic curve, depicting significant subsidence at the onset of the postrift stage, gradually diminishing to minimal values within a few tens of millions of years. Upon comparing the backstripping curves with the stretching subsidence models, particularly during the final flat segments, a residual negative value emerges over time (locally exceeding −500 m), indicating considerable residual subsidence during the Cenozoic postrift evolution along the Colorado Basin. This could be interpreted as dynamic subsidence. This residual sinking in the Argentina shelf matches with dynamic topography numerical models in the study region as well as the occurrence of accumulated slabs in the mantle detected by seismic tomography studies. We propose that the accumulation of slabs along this segment of the South American subduction zone might be the drivers of dynamic subsidence, which assisted to the preservation of a large thickness of post-rift strata along the SW Atlantic margin. • Colorado Basin's post-rifting phase revealed three distinct subsidence stages. • Cenozoic residual subsidence linked to dynamic topography detected. • Dynamic topography unveils mantle slab accumulation from South American subduction. • Atlantic post-rift subsidence differs notably from Africa, influenced by dynamic topography from subducted slabs. [ABSTRACT FROM AUTHOR]

The Molucca Sea region is an active arc-arc collision zone whose structure and mantle dynamics are essential for reconstructing the tectonic evolution of Southeast Asia. The key to understanding the tectonic processes of the region is to clarify the morphology of the subducted slabs. However, geometric details of the subducted Molucca Sea slab and other slab segments under this region remain controversial due to the lack of high-resolution mantle tomography. In this study, we jointly invert a large number of local and teleseismic travel-time data recorded at 99 local seismic stations to determine a high-resolution model of P-wave mantle tomography beneath the Molucca Sea region. Our tomography reveals two aseismic slab segments in the deep upper mantle and the mantle transition zone, which may reflect the earlier subducted Molucca Sea slab. The total length of the Molucca Sea slab is estimated to be ∼1800 to 1900 km, being ∼300 km longer than previous estimates. • The Molucca Sea slab has experienced break-offs. • Aseismic slab fragments broken off from the early subducted Molucca Sea slab are revealed in the mantle transition zone. • The maximum length of the Molucca Sea slab is ∼1800 to 1900 km. [ABSTRACT FROM AUTHOR]

We investigate the 3-D P and S wave velocity structure of the East Japan forearc using a large number of high-quality arrival-time data of local earthquakes recorded by a permanent seafloor seismic network (S-net) and the land-based Kiban seismic network. Geometries of the subducting slab boundary and the Conrad and Moho discontinuities under the forearc are refined by referring to results of previous studies including active-source seismic soundings and seismicity located by the S-net. The obtained P and S wave velocity images are quite similar to each other, suggesting that they reflect robust features of the 3-D seismic structure. The subducting Pacific slab exhibits high velocities, whereas significant low-velocity (low-V) anomalies are revealed in the overlying Okhotsk plate beneath the forearc, which may reflect accretionary materials containing abundant fluids. The subducting oceanic crust atop the Pacific slab is revealed clearly as a thin low-V layer extending from the trench axis down to ~100 km depth under the volcanic front. Our results of earthquake relocation using the 3-D velocity model indicate that the double seismic zone in the subducting Pacific slab occurs from the trench axis down to ~180 km depth. The occurrence and focal mechanisms of the double seismic zone are caused by slab deformation during its subduction, such as slab bending or unbending, as well as other processes including hydration, dehydration, compositional variations, and phase changes in the subducting slab. [Display omitted] • Seafloor S-net and Kiban network data are used to study the 3-D forearc structure. • The Conrad, Moho and Pacific slab geometries in the East Japan forearc are refined. • P and S wave velocity images are quite similar, reflecting robust structural features. • The subducting oceanic crust is revealed down to ~100 km depth under the volcanic front. • The intraslab double seismic zone occurs from the trench axis down to ~180 km depth. [ABSTRACT FROM AUTHOR]

We present a high-resolution 3-D model of P wave tomography of the crust and mantle down to 1600 km depth under Southeast Asia (SE Asia), which is determined by inverting more than 1.3 million P wave arrival time data of local, regional and teleseismic events recorded at 1241 seismic stations deployed in SE Asia and mainland China. Strong and weak smoothing inversions are conducted to reveal different characteristic scales of velocity anomalies in the mantle. Both seismic and aseismic portions of the subducted plates in different parts of the circular subduction system in SE Asia are clearly revealed as high-velocity bodies in our tomographic images. The intraplate volcanoes in Borneo are underlain by significant low-velocity anomalies in the upper mantle, whereas a clear high-velocity anomaly is visible in the mantle transition zone and the lower mantle. We deem that hot and wet upwelling flows in the big mantle wedge above the stagnant slab in the mantle transition zone may have caused the intraplate volcanism in Borneo. A large low-velocity anomaly consisting of a flat head and a twisted tail extending to the bottom of our model is imaged beneath the southeastern Asian basalt province. The low-velocity anomaly possibly represents a hot zone whose shape may have been modified by mantle wind associated with a plate rotation. The hot zone has a lower-mantle origin and probably contains some small mantle upwellings inside, with the Hainan plume being a major trunk. • Aseismic portions of the subducted slabs beneath SE Asia are revealed. • Hot upwellings in the big mantle wedge caused intraplate volcanism in Borneo. • The Hainan plume has a complex geometry due to mantle wind. [ABSTRACT FROM AUTHOR]

The motivation for this work was to model the crustal thickness, and the geometry of the Mohorovičić (Moho) discontinuity, the crust-upper mantle boundary around the Japanese archipelago. The western Pacific Ocean is one of the most eminent subduction zones in the world. The subduction zone complexity is due to alteration within the geometry of the stagnant slab. We compiled a three-dimensional density distribution model for subduction plates around the Japan Islands, and estimated the depth to Moho using the constrained non-linear inversion of gravity data derived from the latest Gravity field and Ocean Circulation Explorer, or GOCE, model known as GOCO06s. The study had three steps, starting with the data corrections followed by the inversion, and ending with the modeling and comparing of the estimated results with previous models and studies. In this study, we used available receiver function results, seismic tomography models, global crustal thickness, and geological and tectonic information as constraints. The model estimated the crustal thickness, surface geometry of the Moho, and geometry of the stagnant slab in the Kuril Trench, Japan Trench, Nankai Trough, and Izu-Bonin Trough. In gravity inversion, we used 34 km and 600 kg m−3 as the Moho depth reference and density contrast, respectively. Results showed that the crustal thickness ranged from 14 to 43 km with a − 0.78 mean misfit for Moho depths derived from seismic P-wave analysis. The calculated gravity data had a 0.16 mGal misfit with the observed data. The normalized analytical signal helped to sharply delineate the different plates and their edges. Three cross-sections were conducted to model the stagnant slab along the subduction zone of the Eastern Eurasian plate, western Pacific plate, Okhotsk plate, and Philippine Sea plate. The eastern area of the entire region (Pacific plate) had a low crust thickness and high density compared to the western area (Eurasian plate), which had a high crust thickness and low density. • 3D gravity inversion of the Northwestern Pacific Plate at the Japan Islands. • Estimating the Moho depths around the Japan Islands. • 3D density modeling of the Stagnant Slabs and Subduction plates. [ABSTRACT FROM AUTHOR]

Although it is well understood that mantle flow in a subduction setting can develop seismic anisotropy, many tomographic studies still assume an isotropic mantle. In this study, we explore how unaccounted-for seismic anisotropy contributes to P-wave anomalies, focusing on apparent sub-slab low velocities. We use data from the westernmost Mediterranean, where the Alboran slab has subducted beneath the Gibraltar Arc. Our isotropic P-wave velocity model shows a low-velocity anomaly at depths of 75–225 km to the west of the high-velocity-anomaly that represents the Alboran slab. We explore whether this low-velocity anomaly could be an anisotropy-induced artifact as a significant shear-wave-splitting (SWS) along the arc is observed. We consider five different hypothetical anisotropy models representing different mantle configurations. The anisotropy models are based on 1) the SWS observations, 2) the depth distribution of the low-velocity anomaly and 3) geodynamic modeling of the development of entrained and toroidal mantle flow in generic subduction zones. Given the underdetermined nature of an anisotropic inversion, we include the anisotropy models as a priori constraints on the tomography. The tomographic results including anisotropy models have smaller travel time residuals and isotropic anomalies in the region of interest while fitting the SWS observation better than the isotropic model. We propose that unaccounted-for seismic anisotropy can produce artifacts in isotropic tomography of real datasets, and this is at least partially the case for the low-velocity anomaly below the Alboran slab. Further, we show that not taking this into account can lead to erroneous conclusions about mantle temperature and the presence of melt. From the insights that our anisotropy models provide, we also suggest that toroidal mantle flow is currently dominant below and around the Alboran slab. • Seismic tomography considering seismic anisotropy in the westernmost Mediterranean. • Anisotropic models improve fit to travel time and SWS observations. • Best-fit model includes only anisotropy from toroidal flow. • Results show isotropic models overestimate the magnitude of sub-slab low-velocities. • Inclusion of anisotropy affects interpretations of temperature and melt content. [ABSTRACT FROM AUTHOR]