Signature of remnant slabs in the North Pacific from P-wave tomography

A 3-D ray-tracing technique was used in a global tomographic inversion in order to obtain tomographic images of the North Pacific. The data reported by the Geophysical Survey of Russia (1955–1997) were used together with the catalogues of the International Seismological Center (1964–1991) and the US Geological Survey National Earthquake Information Center (1991–1998), and the recompiled catalogue was reprocessed. The final data set, used for following the inversion, contained 523 430 summary ray paths. The whole of the Earth's mantle was parametrized by cells of 2° × 2° and 19 layers. The large and sparse system of observation equations was solved using an iterative LSQR algorithm.
A subhorizontal high-velocity anomaly is revealed just above the 660 km discontinuity beneath the Aleutian subduction zone. This high-velocity feature is observed at latitudes of up to ~70°N and is interpreted as a remnant of the subducted Kula plate, which disappeared through ridge subduction at about 48 Ma. A further positive velocity perturbation feature can be identified beneath the Chukotka peninsula and Okhotsk Sea, extending from ~300 to ~660 km depth and then either extending further down to ~800 km (Chukotka) or deflecting along the 660 km discontinuity (Okhotsk Sea). This high-velocity anomaly is interpreted as a remnant slab of the Okhotsk plate accreted to Siberia at ~55 Ma.     

Sequential integrated inversion of refraction and wide-angle reflection traveltimes and gravity data for two-dimensional velocity structures

A new algorithm is presented for the integrated 2-D inversion of seismic traveltime and gravity data. The algorithm adopts the 'maximum likelihood' regularization scheme. We construct a 'probability density function' which includes three kinds of information: information derived from gravity measurements; information derived from the seismic traveltime inversion procedure applied to the model; and information on the physical correlation among the density and the velocity parameters. We assume a linear relation between density and velocity, which can be node-dependent; that is, we can choose different relationships for different parts of the velocity–density grid. In addition, our procedure allows us to consider a covariance matrix related to the error propagation in linking density to velocity. We use seismic data to estimate starting velocity values and the position of boundary nodes. Subsequently, the sequential integrated inversion (SII) optimizes the layer velocities and densities for our models. The procedure is applicable, as an additional step, to any type of seismic tomographic inversion.
We illustrate the method by comparing the velocity models recovered from a standard seismic traveltime inversion with those retrieved using our algorithm. The inversion of synthetic data calculated for a 2-D isotropic, laterally inhomogeneous model shows the stability and accuracy of this procedure, demonstrates the improvements to the recovery of true velocity anomalies, and proves that this technique can efficiently overcome some of the limitations of both gravity and seismic traveltime inversions, when they are used independently.
An interpretation of field data from the 1994 Vesuvius test experiment is also presented. At depths down to 4.5 km, the model retrieved after a SII shows a more detailed structure than the model obtained from an interpretation of seismic traveltime only, and yields additional information for a further study of the area.     

Joint two-dimensional cross-gradient imaging of magnetotelluric and seismic traveltime data for structural and lithological classification

Collocated magnetotelluric (MT) and seismic profiling is emerging as a necessary combined approach for deep and near-surface imaging but the resulting experimental data are typically interpreted separately since no production programs exist for multidimensional joint inversion of MT and seismic data. We present a joint 2-D inversion approach for imaging collocated MT and seismic refraction data with cross-gradient structural constraints. We describe the main features of the algorithm and first apply it to synthetic data generated for a hypothetical complex geological model. For the synthetic data, we find that the scheme leads to models with remarkable structural resemblance and improved estimates of electrical resistivity and seismic velocity. We apply the scheme to near-surface field data to test the consistency of a previously suggested resistivity–velocity interrelationship and its potential use for subsurface lithofacies discrimination or structural classification. The MT-seismic relationship is found to be in excellent accord with that derived previously for DC resistivity and seismic data set at the test site. Our results suggest that joint MT-seismic cross-gradient imaging leads to improved characterization of heterogeneous geological targets at near-surface to mantle depths.     

Seismic evidence for a sharp lithospheric base persisting to the lowermost mantle beneath the Caribbean

Broad-band data from South American earthquakes recorded by Californian seismic networks are analysed using a newly developed seismic wave migration method—the slowness backazimuth weighted migration (SBWM). Using the SBWM, out-of-plane seismic P -wave reflections have been observed. The reflection locations extend throughout the Earth's lower mantle, down to the core–mantle boundary (CMB) and coincide with the edges of tomographically mapped high seismic velocities. Modelling using synthetic seismograms suggests that a narrow (10–15 km) low- or high-velocity lamella with about 2 per cent velocity contrast can reproduce the observed reflected waveforms, but other explanations may exist. Considering the reflection locations and synthetic modelling, the observed out-of-plane energy is well explained by underside reflections off a sharp reflector at the base of the subducted lithosphere. We also detect weaker reflections corresponding to the tomographically mapped top of the slab, which may arise from the boundary between the Nazca plate and the overlying former basaltic oceanic crust. The joint interpretation of the waveform modelling and geodynamic considerations indicate mass flux of the former oceanic lithosphere and basaltic crust across the 660 km discontinuity, linking processes and structure at the top and bottom of the Earth's mantle, supporting the idea of whole mantle convection.     

Double-planed deep seismic zone and upper-mantle structure in the Northeastern Japan Arc

Summary. The ScSp wave converted from the ScS wave at the boundary between the descending lithospheric slab and the mantle above it was clearly observed from a nearby deep earthquake with magnitude 7.7 at some stations of the seismic network of Tohoku University which covers the Tohoku District, the northeastern part of Honshu, Japan. By applying the three-dimensional seismic-ray tracing method, the location of this boundary was determined from the difference in arrival time between the ScS and ScSp waves. The result shows that the upper boundary of the descending slab lies exactly on the upper plane of the double-planed deep seismic zone found in the Northeastern Japan Arc.
There is an additional evidence that the boundary is located on the upper plane of the double-planed deep seismic zone. The hypocentre distribution of intermediate-depth earthquakes located by the small-scale seismic-array observation is extremely different from that obtained by the relatively large-scale seismic network. The discrepancy in the distribution of hypocentres of the same earthquake independently located is well explained by the inclined lithospheric slab model derived from the difference in arrival time between the ScS and ScSp waves.
The earthquakes with reverse faulting or with down-dip compressional stresses occur at the upper boundary of the descending slab. Within the descending slab, the earthquakes with down-dip extensional stresses also occur in a very narrow zone from 30 to 40 km below the dipping boundary in the depth range from 50 to about 200 km, and these shocks form the lower plane of the double-planed deep seismic zone.     

Multimode surface waveform tomography of the Pacific Ocean: a closer look at the lithospheric cooling signature

We present a regional surface waveform tomography of the Pacific upper mantle, obtained using an automated multimode surface waveform inversion technique on fundamental and higher mode Rayleigh waves, to constrain the   V_SV   structure down to ∼400 km depth. We have improved on previous implementations of this technique by robustly accounting for the effects of uncertainties in earthquake source parameters in the tomographic inversion. We have furthermore improved path coverage in the South Pacific region by including Rayleigh wave observations from the French Polynesian Pacific Lithosphere and Upper Mantle Experiment deployment. This improvement has led to imaging of vertical low-velocity structures associated with hotspots within the South Pacific Super-Swell region. We have produced an age-dependent average cross-section for the Pacific Ocean lithosphere and found that the increase in   V_SV   with age is broadly compatible with a half-space cooling model of oceanic lithosphere formation. We cannot confirm evidence for a Pacific-wide reheating event. Our synthetic tests show that detailed interpretation of average   V_SV   trends across the Pacific Ocean may be misleading unless lateral resolution and amplitude recovery are uniform across the region, a condition that is difficult to achieve in such a large oceanic basin with current seismic stations.     

Numerical modelling of wavefields in three-dimensional inhomogeneous media

Summary. Numerical modelling is one of the most efficient methods for an investigation of the relationship between structural features and peculiarities of observed wavefields. It is practically the only method for 2-D and 3-D inhomogeneous media.
An algorithm based on ray theory has been developed for calculations of travel times and amplitudes of seismic waves in 3-D inhomogeneous media with curved interfaces. It was applied for numerical modelling of kinematic and dynamic characteristics of seismic waves propagating in laterally inhomogeneous media.
Travel-time and amplitude patterns were studied in the 2-D and 3-D models of a geosyncline, in which velocity distribution was given by an analytical function of the coordinates. For a more complicated model representing a subducting high-velocity lithospheric plate in a transition zone between oceanic and continental upper mantle, the velocity distribution was given by discrete values on a 2-D non-rectangular grid. It was shown that when a source was placed above the lithospheric plate, a shadow zone appeared along a strike of the structure, i.e. in the direction which is perpendicular to a strong lateral velocity gradient. Travel-time residuals were calculated along the seismological profile for a 3-D velocity distribution in the upper mantle beneath Central Asia, obtained as a result of inversion of travel times by the Backus-Gilbert method. They were found to be in a good agreement with the observed data.     

Active tectonics of the Alpine—Himalayan belt: the Aegean Sea and surrounding regions

Summary. New fault plane solutions, Landsat photographs, and seismic refraction records show that rapid extension is now taking place in the northern and eastern parts of the Aegean sea region. The southern part of the Aegean has also been deformed by normal faulting but is now relatively inactive. In northwestern Greece and Albania there is a band of thrusting near the western coasts adjacent to a band of normal faulting further east. The pre-Miocene geology of the islands in the Aegean closely resembles that of Greece and Turkey, yet seismic refraction shows that the crust is now only about 30 km thick beneath the southern part of the sea, compared with nearly 50 km beneath Greece and western Turkey. These observations suggest that the Aegean has been stretched by a factor of two since the Miocene. This stretching can account for the high heat flow. The sinking slab produced by subduction along the Hellenic Arc may maintain the motions, though the geometry and widespread nature of the normal faulting is not easily explained. The motions in northwestern Greece and Albania cannot be driven in the same way because no slab exists in the area. They may be maintained by blobs of cold mantle detaching from the lower half of the lithosphere, produced by a thermal instability when the lithosphere is thickened by thrusting. Hence generation and destruction of the lower part of the lithosphere may occur beneath deforming continental crust without the production of any oceanic crust.     

Extensional tectonic regimes in the Aegean basins during the Cenozoic

Abstract Kinematics of faults in the Northern Aegean show three extensional tectonic regimes the tensional directions of which trend (1) WNW-ESE, (2) NE-SW and (3) N-S. These were active during the Upper Miocene, Pliocene-Lower Pleistocene and Mid Pleistocene-Present day, respectively. The main characteristics of the stress patterns (1) and (2) on the overall Aegean is tentatively explained by variations of the horizontal lithospheric stress value σ_zz due to the slab push and of the vertical lithospheric stress value σ_zz due to mass heterogeneities. During the Mid Pleistocene-Present, due to the slab push, tectonics were compressional along the arc boundary: σ_zz was σ1. In the Aegean basins, tectonics were extensional, c_2Z was σ1 as a consequence of the thickness of the continental crust and, possibly of an updoming asthenosphere; thus σ_zz became σ2, allowing tension σ3 to be orthogonal to the compression along the arc, i.e. to be roughly parallel to the arc trend. During the Pliocene-Lower Pleistocene, the extensional regime was distinctly different. The tensional directions were roughly radial to the arc. It is suggested that σ_zz was weakly compressional, or eventually tensional, due a seaward migration of the slab so that σ_zz became σ3. In the Northern Aegean, the stress pattern has been also controlled by the westward push of the Anatolian landmass. During the Mid Pleistocene-Present day, this was typically extensional (al was vertical) and the right lateral strike-slip motion on the North Anatolian Fault transformed into a N-S-stretching, E-W-shortening of the Northern Aegean. Dextral strike-slip motions along the North Aegean Trough fault zone were possible on NE-SW-striking faults. During the Pliocene-Lower Pleistocene, normal fault components were higher; however, because the angle between the NE-SW trend of the tensional axis and the strike of the fault zone was acute, dextral strike-slip components were possible on all the faults striking NE-SW to E-W. A clockwise 15^o rotation of Limnos with respect to Samothraki, Thraki and Thassos, suggested by structural data, was probably associated with these dextral motions. The WNW-ESE trending tension during the Upper Miocene indicates that the dextral North Anatolian Fault had not yet merged into the North Aegean Trough fault zone at that time. We propose that the formation of Aegean basins during the Cenozoic was related to the activity of two major Hellenic arcs. The ‘Pelagonian-Pindic Arc’ resulted in the formation of the subsident Aegean basins of Middle Eocene-Lower Miocene age and of the older Northern Aegean orogenic volcanism. The ‘Aegean Arc’ resulted in the formation of the subsident Aegean basins of Middle Miocene to Present day age and of the Southern Aegean orogenic volcanism. Were these arcs associated with a unique subduction zone or with two such zones ? In the first case, the slab is no more than 16 Myr old, in the second it may be as old as 45–50 Myr. The answer depends on the accuracy of the seismic tomography profiles.     

The transition from subduction to continental collision: crustal structure in the North Canterbury region, New Zealand

The North Canterbury region marks the transition from Pacific plate subduction to continental collision in the South Island of New Zealand. Details of the seismicity, structure and tectonics of this region have been revealed by an 11-week microearthquake survey using 24 portable digital seismographs. Arrival time data from a well-recorded subset of microearthquakes have been combined with those from three explosions at the corners of the microearthquake network in a simultaneous inversion for both hypocentres and velocity structure. The velocity structure is consistent with the crust in North Canterbury being an extension of the converging Chatham Rise. The crust is about 27 km thick, and consists of an 11 km thick seismic upper crust and 7 km thick seismic lower crust, with the middle part of the crust being relatively aseismic. Seismic velocities are consistent with the upper and middle crust being composed of greywacke and schist respectively, while several lines of evidence suggest that the lower crust is the lower part of the old oceanic crust on which the overlying rocks were originally deposited.
The distribution of relocated earthquakes deeper than 15 km indicates that the seismic lower crust changes dip markedly near 43S. To the south-west it is subhorizontal, while to the north-east it dips north-west at about 10. Fault-plane solutions for these earthquakes also change near 43S. For events to the south, P -axes trend approximately normal to the plate boundary (reflecting continental collision), while for events to the north, T -axes are aligned down the dip of the subducted plate (reflecting slab pull). While lithospheric subduction is continuous across the transition, it is not clear whether the lower crust near 43S is flexed or torn.     

Three-dimensional lithospheric modelling at NORSAR: linearity of the method and amplitude variations from the anomalies

Summary. A modification of the Aid et al . technique for three-dimensional lithospheric modelling is used to find smoothly varying models for the P -wave velocity structure beneath NORSAR. The method includes ray tracing and calculation of geometrical spreading in the anomalies. The results of linear inversion of the travel-time data compare well with those of previous investigators. The assumption of linearity, which removes the need to ray trace through the anomalies, is tested with iterative solutions for both synthetic and real data. A model with an rms velocity perturbation of 3 per cent, extending to 120 km depth, is found to be reasonably linear. In fact the procedure leads to two models which satisfy the same amount of the real data but which differ by far more than the standard errors. However, these differences are not significant once the imperfect resolution is accounted for by using the total estimation error of the stochastic inverse.
The depth of major anomalies appears to be greater than the array diameter and is therefore not well constrained. Comparing the geometrical spreading produced by these models with the amplitude variations observed at the array indicates that structure deeper than 120 km but shallower than 200 km makes an important contribution to the observations. None of the models used can produce variations as large as those in the amplitude data. For deep, essentially two-dimensional, anomalies the fit to these data is much better for sources to the NE of the array than for sources in other quadrants.     

Crustal thickness estimation beneath the southern central Andes at 30°S and 36°S from S wave receiver function analysis

We have used the S wave receiver function (SRF) technique to investigate the crustal thickness beneath two seismic profiles from the CHARGE project in the southern central Andes. A previous study employing the P wave receiver function method has observed the Moho interface beneath much of the profiles. They found, however, that the amplitude of the P to S conversion was diminished in the western part of the profiles and have attributed it to a reduction of the impedance contrast at the Moho due to lower crustal ecologitization. With SRF, we have successfully detected S to P converted waves from the Moho as well as possible conversions from other lithospheric boundaries. The continental South American crust reaches its maximum thickness of ∼70 km (along 30°S between 70°W and 68.5°W) beneath the Principal Cordillera and the Famatina system and becomes thinner towards the Sierras Pampeanas with a thickness of ∼40 km. Negative phases, possibly related to the base of the continental and oceanic lithosphere, can be recognized in the summation traces at different depths. By comparing our results with data obtained from previous investigations, we are able to further constrain the thickness of the crust and lithosphere beneath the central Andes.     

Crustal tomographic imaging of a transitional continental rift: the Ethiopian rift

In this study we image crustal structure beneath a magmatic continental rift to understand the interplay between crustal stretching and magmatism during the late stages of continental rifting: the Main Ethiopian Rift (MER). The northern sector of this region marks the transition from continental rifting in the East African Rift to incipient seafloor spreading in the southern Red Sea and western Gulf of Aden. Our local tomographic inversion exploits 172 broad-band instruments covering an area of 250 × 350 km of the rift and adjacent plateaux. The instruments recorded a total of 2139 local earthquakes over a 16-month period. Several synthetic tests show that resolution is good between 12 and 25 km depth (below sea level), but some horizontal velocity smearing is evident along the axis of the Main Ethiopian Rift below 16 km. We present a 3-D P -wave velocity model of the mid-crust and present the first 3-D Vp / Vs model of the region. Our models show high P -wave velocities (6.5 km s⁻¹) beneath the axis of the rift at a depth of 12–25 km. The presence of high Vp / Vs ratios (1.81–1.84) at the same depth range suggest that they are cooled mafic intrusions. The high Vp / Vs values, along with other geophysical evidence, suggest that dyking is pervasive beneath the axis of the rift from the mid-crustal depths to the surface and that some portion of partial melt may exist at lower crustal depths. Although the crustal stretching factor across the Main Ethiopian Rift is ∼1.7, our results indicate that magma intrusion in narrow zones accommodates a large proportion of extensional strain, with similarities to slow-spreading mid-ocean ridge processes.     

Seismotectonic implications of relocated aftershock sequences in Iran and Turkey

Summary. Six aftershock sequences in Iran and Turkey are relocated using existing teleseismic data. Two of these are in the Zagros mountains where local fieldwork has failed to detect subcrustal seismicity but published teleseismic locations show depths greater than 100 km. All apparently deep events are shown to be small and badly recorded with poor depth resolution. There is thus no evidence for active lithospheric subduction in the Zagros.
Relocations of other sequences in Iran and Turkey are used with fault plane solutions, satellite photographs and surface faulting to provide new insight on the geometry of faulting and crustal deformation of those regions. Linear seismic trends from these sequences are shown to cut older geological structures and do not always bear a simple relation to surface faulting. In such cases aftershock activity may be on primary buried faults whose behaviour is not simply revealed in surface structure and deformation.
A linearized inversion scheme is used to investigate the trade-off between resolution and uncertainty in the hypocentral parameters. The ultimate resolution of teleseismic locations is shown to be limited by the quality of arrival time data.     

P-wave traveltime and polarization tomography of VSP data

The presence of anisotropy requires that tomographic methods be generalized to account for anisotropy. This generalization allows geological structure to be correctly imaged and allows the anisotropic parameters to be estimated. Use of isotropic inversion for imaging anisotropic structures gives systematic trends in the traveltime and polarization residuals. However, due to the limited directional coverage, the traveltimes along may not be sufficient to study the anisotropic properties of the structure. Polarizations can provide independent information on the structure. Traveltime and polarization inversion are applied to synthetic examples simulating VSP experiments. Transverse isotropy and 1-D structure are assumed. Plots of traveltime and polarization residuals are an important tool to detect the anomalies due to the presence of anisotropy. For receivers located in anisotropic layers, polarization residuals display consistent anomalies of several degrees. The synthetic examples show that even the simple 1-D problem is difficult, when using direct arrivals only. Large a posteriori errors in anisotropic parameters are obtained by traveltime inversion in layers where available incidence angles are less than 45°. Resolution of the tomographic image of VSP data is greatly improved by a combination of traveltime and polarization information. In order to obtain accurate inversion results, the measurement error of polarization data should be kept to within a few degrees.     

New insight into the crust and upper mantle structure under Alaska

To better understand the seismic structure of the subducting Pacific plate under Alaska, we determined the three-dimensional P-wave velocity structure to a depth of approximately 200 km beneath Alaska using 438,146 P-wave arrival times from 10,900 earthquakes. In this study an irregular grid parameterization was adopted to express the velocity structure under Alaska. The number of grid nodes increases from north to south in the study area so that the spacing between grid nodes is approximately the same in the longitude direction. Our results suggest that the subducting Pacific slab under Alaska can be divided into three different parts based on its geometry and velocity structure. The western part has features similar to those in other subduction zones. In the central part a thick low-velocity zone is imaged at the top of the subducting Pacific slab beneath north of the Kenai Peninsula, which is believed to be most likely the oceanic crust plus an overlying serpentinized zone and the coupled Yakutat terrane subducted with the Pacific slab. In the eastern part, significant high-velocity anomalies are visible to 60–90 km depth, suggesting that the Pacific slab has only subducted down to that depth.     

Tectonic evolution of the Alboran Sea basin

The Alboran Sea is an extensional basin of Neogene age that is surrounded by highly arcuate thrust belts. Multichannel seismic (MCS) reflection profile data suggest the basin has a complex tectonic fabric that includes extensional, compressional and strike-slip structures. The early Miocene history appears to be dominated by graben formation with border faults that are in large part contemporaneous with thrust movements in the external zones of the Betic and Rif mountains. Extension appears to have continued into the late Miocene although the main movements were probably completed by the time of the Messinian ‘salinity crisis’. The Pliocene and younger history of the basin is dominated by infilling of the Messinian topography, gentle subsidence, and extensional, compressional and strike-slip movements. There is evidence from the sea-floor morphology and seismicity patterns that the basin is actively deforming in response to present-day plate motions. Backstripping of well data in the basin margin suggests that the initial extensional event was accompanied by crustal and lithospheric thinning. The depth to Moho inferred from backstripping is greater than the depth expected based on seismic and gravity modelling, suggesting that backstripping underestimates the true amount of thinning. One explanation is that some of the thinning occurred while the crust was above sea level, perhaps as a result of either crustal thickening, or a period of lithospheric heating and thinning, prior to rifting. We found that a model with a ‘normal’ crustal thickness of 31.2 km, a lithospheric thickness of 50 km, and β= 1.4 predicts 0.8 km of initial uplift. These parameters fit the well subsidence data and bring the backstripped Moho into better agreement with the seismic and gravity Moho. The origin of such a thin lithosphere is not constrained by the data, but we believe that it may be a result of the detachment of a cold lithospheric ‘root’ that formed during pre-Neogene collisional orogeny in the region.     

Magnetotelluric image of the crust and upper mantle in the backarc of the northwestern Argentinean Andes

Magnetotelluric data from the backarc of the Central Andes in NW Argentinawere re-examined by employing impedance tensor decomposition and 2-D inversion and modelling techniques. The data in the period range of 50–15 000 s were collected on a profile of 220 km length reaching from the Eastern Cordillera across the Santa Barbara System to the Andean foreland of the Argentinean Chaco.
After a dimensionality analysis, data from most sites were treated as regional 2-D. The exception was the eastern section of the profile, where the magnetotelluric transfer functions for periods ≤ 1000 s reflect a 3-D earth. Application of two tensor decomposition schemes yielded a regional strike direction of N–S, which is the azimuth of the Central Andean mountain chains. Several 2-D models were obtained by pseudo- and full 2-D Occam inversion schemes. Special emphasis was placed on the inversion of phase data to reduce the influence of static shifts in the apparent resistivity data. The smooth inversion models all show a good conductor at depth. A final model was then calculated using a finite element forward algorithm.
The most prominent feature of the resulting model is a conductor which rises from depths of 180 km below the Chaco region to 80 km beneath the Santa Barbara System and the Eastern Cordillera. Its interpretation as a rise of the electrical asthenosphere is supported by seismic attenuation studies. Magnetotelluric results, surface heat-flow distribution in the area, and the electrical properties of crustal and mantle rocks suggest that the upper mantle is predominantly ductile beneath the Eastern Cordillera and the western Santa Barbara System. This generally agrees with anelastic seismic attenuation models of the area and is useful in discriminating between models of Q quality factor distribution.     

The Acre Basin basement (NW Brazil) and the transition from the intracratonic to retroarc foreland basin system

The understanding of the crustal transition between orogenic zones and cratonic portions in distal regions of foreland basins has received increasing attention, but the analysis is often hampered by the sedimentary cover. Despite the peculiar location of the Acre Basin, specifically between the Amazonian Craton and the sub-Andean zone, local basement studies are still scarce due to lacking seismic data and exploratory wells. Therefore, this work aims to map basement depths, estimate crustal compositions and identify the main depocenters, structures and limits of Acre Basin using an integrated analysis to understand better the region lithospheric evolution, its relationship with the Amazonian Craton and its positioning within the Andean orogeny. For that, we used well, 2D seismic reflection, airborne and ground gravity and magnetic data as well as the EMG2008. Tilt Depth estimates indicate basement depths between 500 and 7800 m and larger sedimentary thicknesses in the northern portion. Additionally, we modelled groups of potential sources between 0.1 and 22 km and Moho depths between 26 and 37 km. Compositionally, the upper crust consists dominantly of meta-sedimentary and low-grade metamorphic rocks and granites, indicating that the sub-Andean and Acre Basins share a similar basement. Thus, there are indications that the basement of the Acre Basin is essentially formed by the Sunsás province in the Amazonian Craton. However, local differences in basement depth, magnetic susceptibility and exploratory potential led to the subdivision into Divisor and Xapuri sub-basins, north and south of the Fitzcarrald Arch, respectively. Finally, it was possible to establish the limits of the Andean orogeny influence in the Acre Basin and delimit the area of the Western Amazon Foredeep installed during the Neogene.     

Crustal structure of the Newfoundland rifted continental margin from constrained 3-D gravity inversion

The rifting history of the Atlantic continental margin of Newfoundland is very complex and so far has been investigated at the crustal scale primarily with the use of 2-D seismic surveys. While informative, the results generated from these surveys cannot easily be interpreted in a regional sense due to their sparse sampling of the margin. A 3-D gravity inversion of the free air data over the Newfoundland margin allows us to generate a 3-D density anomaly model that can be compared with the seismic results and used to gain insight into regions lacking seismic coverage. Results of the gravity inversion show good correspondence with Moho depths from seismic results. A shallowing of the Moho to 12 km depth is resolved on the shelf at the northern edge of the Grand Banks, in a region poorly sampled by other methods. Comparisons between sediment thickness and crustal thickness show deviations from local isostatic compensation in locations which correlate with faults and rifting trends. Such insights must act as constraints for future palaeoreconstructions of North Atlantic rifting.