We present a new gravity map between 45°-70° W and 5°-40° S integrating open source terrestrial gravity data of Argentina with the South American Gravity Model 2004 (SAGM04), a 5 min-arc resolution gravity model. The Bouguer anomaly map reveals a 2,000 km long linear gravity feature from 15° S to 30° S at longitude 55° W, with a steep horizontal gradient separating two gravity domains. The eastern domain is the Paraná basin, with NE-SW trending Bouguer anomalies of-80 mGal in average. The western domain comprises the Chaco-Paraná, Chaco-Tarija and Pantanal basins, with circular positive anomalies of up to 20 mGal in amplitude. Previous seismic studies mapped a thinner crust of less than 35 km in the western domain and the present gravity models indicate a 10 to 20 kg/m³ denser crust. On the other hand, the eastern domain has a thicker crust of more than 40 km. Seismic tomography models also show P-and S-wave velocity reduction in the western domain whereas high-velocity characterises the Paraná basin. These geophysical data indicate that the gravity gradient marks a transition between two distinct lithospheres. The gravity gradient is associated with a tectonic feature referred to as the Western Paraná suture/shear zone. Granites of 530 to 570 Ma ages, located parallel or over the gravity gradient, suggest a Neoproterozoic to Early Cambrian age suture/shear zone, thus approximately synchronous and parallel to the Pampean belt. Sediment corrected residual gravity map and its vertical derivative allow us to define the limits of the Rio Apa, Rio de la Plata and Rio Tebicuary cratons. Their eastern and western limits are the Western Paraná suture and the Pampean belt, respectively. This study unravels Precambrian tectonic elements concealed by the Phanerozoic sedimentary basins adding new constraints for the amalgamation history of SW Gondwana. Research Highlights Gravity map reveals the Neoproterozoic 2,000 km long Western Paraná suture/shear zone Geophysical delimitation of the Rio Apa, Rio Tebicuary and Rio de la Plata cratons Gravity anomalies of the Amazonian and Rio Apa cratons are distinct New tectonic features of SW Gondwana final amalgamation revealed Key Words SW Gondwana; intracontinental basins; cratons; sutures; gravity Research Highlights Gravity map reveals the Neoproterozoic 2,000 km long Western Paraná suture/shear zone Geophysical delimitation of the Rio Apa, Rio Tebicuary and Rio de la Plata cratons Gravity anomalies of the Amazonian and Rio Apa cratons are distinct New tectonic features of SW Gondwana final amalgamation revealed

We estimate density and P‐wave velocity perturbations in the mantle beneath the southeastern South America plate from geoid anomalies and P‐wave traveltime residuals to constrain the structure of the lithosphere underneath the Paraná Magmatic Province (PMP) and conterminous geological provinces. Our analysis shows a consistent correlation between density and velocity anomalies. The P‐wave speed and density are 1% and 15 kg/m3 lower, respectively, in the upper mantle under the Late Cretaceous to Cenozoic alkaline provinces, except beneath the Goiás Alkaline Province (GAP), where density (+20 kg/m3) and velocity (+0.5%) are relatively high. Underneath the PMP, the density is higher by about 50 kg/m3 in the north and 25 kg/m3 in the south, to a depth of 250 − 300 km. These values correlate with high‐velocity perturbations of +0.5% and +0.3%, respectively. Profiles of density perturbation versus depth in the upper mantle are different for the PMP and the adjacent Archean São Francisco (SFC) and Amazonian (AC) cratons. The Paleoproterozoic PMP basement has a high‐density root. The density is relatively low in the SFC and AC lithospheres. A reduction of density is a typical characteristic of chemically depleted Archean cratons. A more fertile Proterozoic and Phanerozoic subcontinental lithospheric mantle has a higher density, as deduced from density estimates of mantle xenoliths of different ages and composition. In conjunction with Re‐Os isotopic studies of the PMP basalts, chemical and isotopic analyses of peridodite xenoliths from the GAP in the northern PMP, and electromagnetic induction experiments of the PMP lithosphere, our density and P‐wave speed models suggest that the densification of the PMP lithosphere and flood basalt generation are related to mantle refertilization. Metasomatic refertilization resulted from the introduction of asthenospheric components from the mantle wedge above Proterozoic subduction zones, which surrounded the Paraná lithosphere. The high‐density PMP lithosphere is presently gravitationally unstable and prone to delamination.

A palaeogeographical reconstruction of the South American and African continents back to anomaly C34 (84 Ma) brings together the Rio Grande Rise (RGR) and the central portion of the Walvis Ridge (WR), thus the RGR–WR aseismic ridges may have a common origin. If the construction of the RGR–WR basaltic plateau took place mainly between 89 and 78 Ma, as indicated by the ages of the basalts sampled by DSDP wells, then the basaltic magmas are the result of an ‘on-ridge’ volcanism. Once separated, the normal sea-floor spreading and thermal subsidence of the RGR and WR ridges continued until approximately 47 Ma when an Eocene magmatism took place in the RGR. In the WR, a younger volcanism is observed in the Guyot Province. The available geochemical and isotope data of the WR–RGR basalts do not indicate the participation of the continental crust melting component. Incompatible trace element ratios and isotope signatures of the basalts from the RGR–WR ridges are distinct from the present-day Tristan da Cunha alkaline rocks, and are nearly identical to the high-Ti Paraná Magmatic Province (PMP) tholeiites (133–132 Ma). Both the high-Ti PMP and the WR–RGR basalts are characterized by moderate initial 87Sr/86Sr and low 206Pb/204Pb isotope ratios [Enriched Mantle I (EMI) mantle component], suggesting melting from a common source, with significant participation of sub-continental lithospheric mantle (SCLM). A three-dimensional (3D) flexural modelling of the RGR and WR was conducted using ETOPO1 digital topography/bathymetry and EGM2008-derived free-air anomalies as a constraint. The best fit between the observed and calculated free-air anomalies was obtained for an elastic plate with elastic plate thickness (Te) of less than 5 km, consistent with an ‘on-ridge’ initial construction of the RGR–WR. The modelling of the crust–mantle interface depths indicates a total crustal thickness of up to 30 km in the RGR–WR. Flexural analysis reinforces the geological evidence that RGR was constructed during two main magmatic episodes, the tholeiitic basalts in the Santonian–Conician times and the alkaline magmatism in the Eocene. Geochemical and geophysical evidence, which rules out the classical deep-mantle plume model in explaining the generation of basalts of these volcanic provinces, is presented. Finally, three models to explain the geochemical and isotope signatures of RGR–WR basalts are reviewed: (1) thermal erosion of SCLM owing to edge-driven convection; (2) melting of fragmented or detached SCLM and lower crust; and (3) thermal erosion at the base of the SCLM with lateral transport of enriched components by mantle flow.

[1] We developed a three-dimensional scheme to invert geoid anomalies aiming to map density variations in the mantle. Using an ellipsoidal-Earth approximation, the model space is represented by tesseroids. To assess the quality of the density models, the resolution and covariance matrices were computed. From a synthetic geoid anomaly caused by a plume tail with Gaussian noise added, the inversion code was able to recover a plausible solution about the density contrast and geometry when it is compared to the synthetic model. To test the inversion algorithm in a natural case study, geoid anomalies from the Yellowstone Province (YP) were inverted. From the Earth Gravitational Model 2008 expanded up to degree 2160, lower crust-and mantle-related negative geoid anomalies with amplitude of approximately 70 m were obtained after removing long-wavelength components (>5400 km) and crustal effects. We estimated three density models for the YP. The first model, the EDM-1 (estimated density model), uses a starting model with density contrast equal to 0. The other two models, the EDM-2 and EDM-3, use an initial density derived from two S-velocity models for the western United States, the Dynamic North America Models of S Waves by Obrebsky et al. (2011) and the Northwestern United States Teleseismic Tomography of S Waves (NWUS11-S) by James et al. (2011). In these three models, a lower and an upper bound for the density solution was also imposed as a priori information. Regardless of the initial constraints, the inversion of the residual geoid indicates that the lower crust and the upper mantle of the YP have a predominantly negative density contrast (~À50 kg/m 3 ) relative to the surrounding mantle. This solution reveals that the density contrast extends at least to 660 km depth. Regional correlation analysis between the EDM-1 and NWUS11-S indicates an anticorrelation (coefficient of~À0.7) at 400 km depth. Our study suggests that the mantle density derived from the inversion of geoid could be integrated with seismic velocity models to image mantle anomalous features beyond the depth limit of investigation achieved combining gravity and seismic tomography.

No abstract

Summary Long-period (T > 10 s) shear-wave reflections between the surface and reflecting boundaries below seismic stations are useful for studying phase transitions in the mantle transition zone (MTZ) but shear-velocity heterogeneity and finite-frequency effects complicate the interpretation of waveform stacks. We follow up on a recent study by Shearer & Buehler (2019) (SB19) of the top-side shear-wave reflection Ssds as a probe for mapping the depths of the 410-km and 660-km discontinuities beneath the USArray. Like SB19, we observe that the recorded Ss410s-S and Ss660s-S traveltime differences are longer at stations in the western US than in the central-eastern US. The 410-km and 660-km discontinuities are about 40–50 km deeper beneath the western US than the central-eastern US if Ss410s-S and Ss660s-S traveltime differences are transformed to depth using a common-reflection point (CRP) mapping approach based on a 1-D seismic model (PREM in our case). However, the east-to-west deepening of the MTZ disappears in the CRP image if we account for 3-D shear-wave velocity variations in the mantle according to global tomography. In addition, from spectral-element method synthetics, we find that ray theory overpredicts the traveltime delays of the reverberations. Undulations of the 410-km and 660-km discontinuities are underestimated when their wavelengths are smaller than the Fresnel zones of the wave reverberations in the MTZ. Therefore, modeling of layering in the upper mantle must be based on 3-D reference structures and accurate calculations of reverberation traveltimes.

Regional waveforms of deep‐focus Tonga‐Fiji earthquakes indicate anomalous traveltime differences (ScS2‐ScS) and amplitude ratios (ScS2/ScS) of the phases ScS and ScS2. The correlation between the ScS2‐ScS delay time and the ScS2/ScS amplitude ratio suggests that shear wave apparent Q in the mantle below the Tonga‐Fiji region is highest when shear wave velocities are lowest. This observation is unexpected if temperature variations were responsible for the seismic anomalies. Using spectral element method waveform simulations for four tomographic models, we demonstrate that focusing and scattering of shear waves by long‐wavelength 3‐D heterogeneity in the mantle may overwhelm the signal from intrinsic attenuation in long‐period ScS2/ScS amplitude ratios. The tomographic models reproduce the trends in recorded ScS2‐ScS difference times and ScS2/ScS amplitude ratios. Although they cannot be ruled out, variations in shear wave attenuation (i.e., the quality factor Q) are not necessary to explain the data.

Summary Based on new data from permanent and temporary networks, we present fundamental mode Rayleigh wave group velocity maps at periods of 10-150 s related to the lithosphere beneath South America. We analyse waveform data from 1043 earthquakes, from 2002 to 2019, which were recorded by 282 stations. To isolate fundamental mode Rayleigh waves, a phase-matched filter is applied, and the measurements of group velocity are obtained from multiple filter analysis techniques. Thus, we obtain 17838 paths, covering most of the South American continent, which reach their maximum at the period of 30 s and decrease for both shorter and longer periods. We calculate average dispersion curves and probability density distribution of all measured curves to check the consistency of our dataset. Then, regionalised group velocity maps are obtained by iteratively combining the fast marching method and the subspace inversion method. The resolution of our models is assessed by checkerboard tests, which show that the synthetic group velocities are well recovered, despite some amplitude and smearing effects in some portions of the model, probably owing to regularisation and uneven raypath coverage. Compared to previous group velocity studies for South America, our models present better resolution, mainly for shorter periods. Our maps of 10 and 20 s, for example, show an excellent correlation with the sedimentary thickness (CRUST1.0) and topography density (UNB$\_$TopoDens). Regions of exposed basement and high-density are related to fast group velocities, while sedimentary basins and low-densities are observed as areas of slow group velocities. We identify small-scale fast group velocity heterogeneities that may be linked to the Rio Apa and Rio Tebicuary cratons as well as to the geochronological provinces of the Amazonian Craton. The most striking feature of our map at 40 s is a fast group velocity structure with the same NE trend of the Transbrasiliano lineament, a Neoproterozoic megashear fault that crosses a large part of the South American continent. Our long-period maps sample lithospheric depths, revealing that cratonic areas of South America, such as the Amazonian and São Francisco cratons, correlate well with fast group velocities. Another interesting feature is the presence of a strong group velocity gradient between the Paraná and Chaco-Paraná basins, which nearly coincides with the location of the Western Paraná Suture, a continental-scale gravity discontinuity. From our group velocity maps, we estimate 1D S-velocity depth profiles at ten locations in South America: Chaco-Tarija Basin, Borborema Province (BP), Amazonian Craton, Paraná Basin, Tocantins Province, Acre Basin (AcB), Altiplano-Puna Volcanic Complex, Mantiqueira Province (MP), Parnaába Basin, and São Francisco Craton. Most of our inverted S-velocity profiles show good agreement with the SL2013sv model at lithospheric depths, except the BP, AcB, and MP profiles. Particularly for the BP, a low shear-wave velocity, from about 75 to 150 km depth, is a feature that is not present in the SL2013sv model and was probably resolved in our model because of our denser raypath coverage. This decreased S-velocity may be due to a lithospheric thinning beneath the BP, as already pointed out by previous studies.