Deep crustal structure across a young passive margin from wide-angle and reflection seismic data (The Sardinia Experiment) - I. Gulf of Lion's margin

The conjugate margins system of the Gulf of Lion and West Sardinia (GLWS) represents a unique natural laboratory for addressing fundamental questions about rifting due to its landlocked situation, its youth, its thick sedimentary layers, including prominent palaeo-marker such as the MSC event, and the amount of available data and multidisciplinary studies. The main goals of the SARDINIA experiment, were to (i) investigate the deep structure of the entire system within the two conjugate margins: the Gulf of Lion and West Sardinia, (ii) characterize the nature of the crust, and (iii) define the geometry of the basin and provide important constrains on its genesis. This paper presents the results of P-wave velocity modelling on three coincident near-vertical reflection multi-channel seismic (MCS) and wide-angle seismic profiles acquired in the Gulf of Lion, to a depth of 35 km. A companion paper [part II Afilhado et al., 2015] addresses the results of two other SARDINIA profiles located on the oriental conjugate West Sardinian margin. Forward wide-angle modelling of both data sets confirms that the margin is characterised by three distinct domains following the onshore unthinned, 33 km-thick continental crust domain: Domain I is bounded by two necking zones, where the crust thins respectively from 30 to 20 and from 20 to 7 km over a width of about 170 km; the outermost necking is imprinted by the well-known T-reflector at its crustal base; Domain II is characterised by a 7 km-thick crust with anomalous velocities ranging from 6 to 7.5 km/s; it represents the transition between the thinned continental crust (Domain I) and a very thin (only 4-5 km) "atypical" oceanic crust (Domain III). In Domain II, the hypothesis of the presence of exhumed mantle is falsified by our results: this domain may likely consist of a thin exhumed lower continental crust overlying a heterogeneous, intruded lower layer. Moreover, despite the difference in their magnetic signatures, Domains II and III present the very similar seismic velocities profiles, and we discuss the possibility of a connection between these two different domains.

Deep crustal structure across a young passive margin from wide-angle and reflection seismic data (The Sardinia Experiment) - II. Sardinia's margin

Geophysical data acquired on the conjugate margins system of the Gulf of Lion and West Sardinia (GLWS) is unique in its ability to address fundamental questions about rifting (i.e. crustal thinning, the nature of the continent-ocean transition zone, the style of rifting and subsequent evolution, and the connection between deep and surface processes). While the Gulf of Lion (GoL) was the site of several deep seismic experiments, which occurred before the SARDINIA Experiment (ESP and ECORS Experiments in 1981 and 1988 respectively), the crustal structure of the West Sardinia margin remains unknown. This paper describes the first modeling of wide-angle and near-vertical reflection multi-channel seismic (MCS) profiles crossing the West Sardinia margin, in the Mediterranean Sea. The profiles were acquired, together with the exact conjugate of the profiles crossing the GoL, during the SARDINIA experiment in December 2006 with the French R/V L'Atalante. Forward wide-angle modeling of both data sets (wide-angle and multi-channel seismic) confirms that the margin is characterized by three distinct domains following the onshore unthinned, 26 km-thick continental crust : Domain V, where the crust thins from 26 to 6 km in a width of about 75 km; Domain IV where the basement is characterized by high velocity gradients and lower crustal seismic velocities from 6.8 to 7.25 km/s, which are atypical for either crustal or upper mantle material, and Domain III composed of "atypical" oceanic crust.The structure observed on the West Sardinian margin presents a distribution of seismic velocities that is symmetrical with those observed on the Gulf of Lion's side, except for the dimension of each domain and with respect to the initiation of seafloor spreading. This result does not support the hypothesis of simple shear mechanism operating along a lithospheric detachment during the formation of the Liguro-Provencal basin.

Deep crustal structure of the Parana Basin from receiver functions and Rayleigh-wave dispersion: Evidence for a fragmented cratonic root

The deep crustal structure of the Parana Basin of southern Brazil is investigated by analyzing P- and PP-wave receiver functions at 17 Brazilian Lithosphere Seismic Project stations within the basin. The study area can be described as a typical Paleozoic intracratonic basin that hosts one of the largest Large Igneous Province of the world and makes a unique setting for investigating models of basin subsidence and their interaction with mantle plumes. Our study consists of (1) an analysis of the Moho interaction phases in the receiver functions to obtain the thickness and bulk Vp/Vs ratio of the basin`s underlying crust and (2) a joint inversion with Rayleigh-wave dispersion velocities from an independent tomographic study to delineate the detailed S-wave velocity variation with depth. The results of our analysis reveal that Moho depths and bulk Vp/Vs ratios (including sediments) vary between 41 and 48 km and between 1.70 and 1.76, respectively, with the largest values roughly coinciding with the basin`s axis, and that S-wave velocities in the lower crust are generally below 3.8 km/s. Select sites within the basin, however, show lower crustal S-wave velocities slightly above 3.9 km/s suggestive of underplated mafic material. We show that these observations are consistent with a fragmented cratonic root under the Parana basin that defined a zone of weakness for the initial Paleozoic subsidence of the basin and which allowed localized mafic underplating of the crust along the suture zones by Cenozoic magmatism.

Deep Crustal Structure of the North-West African Margin from Combined Wide-angle and Reflection Seismic Data (MIRROR Seismic Survey)

The structure of the Moroccan and Nova Scotia conjugate rifted margins is of key importance for understanding the Mesozoic break-up and evolution of the northern central Atlantic Ocean basin. Seven combined multichannel reflection (MCS) and wide-angle seismic (OBS) data profiles were acquired along the Atlantic Moroccan margin between the latitudes of 31.5° and 33° N during the MIRROR seismic survey in 2011, in order to image the transition from continental to oceanic crust, to study the variation in crustal structure and to characterize the crust under the West African Coast Magnetic Anomaly (WACMA). The data were modeled using a forward modeling approach. The final models image crustal thinning from 36 km thickness below the continent to approximately 8 km in the oceanic domain. A 100 km wide zone characterized by rough basement topography and high seismic velocities up to 7.4 km/s in the lower crust is observed westward of the West African Coast Magnetic Anomaly. No basin underlain by continental crust has been imaged in this region, as has been identified north of our study area. Comparison to the conjugate Nova Scotian margin shows a similar continental crustal thickness and layer geometry, and the existence of exhumed and serpentinized upper mantle material on the Canadian side only. The oceanic crustal thickness is lower on the Canadian margin.

Deep crustal electromagnetic structure of Bhuj earthquake region (India) and its implications

The existence of fluids and partial melt in the lower crust of the seismically active Kutch rift basin (on the western continental margin of India) owing to underplating has been proposed in previous geological and geophysical studies. This hypothesis is examined using magnetotelluric (MT) data acquired at 23 stations along two profiles across Kutch Mainland Uplift and Wagad Uplift. A detailed upper crustal structure is also presented using twodimensional inversion of MT data in the Bhuj earthquake (2001) area. The prominent boundaries of reflection in the upper crust at 5, 10 and 20 km obtained in previous seismic reflection profiles correlate with conductive structures in our models. The MT study reveals 1-2 km thick Mesozoic sediments under the Deccan trap cover. The Deccan trap thickness in this region varies from a few meters to 1.5 km. The basement is shallow on the northern side compared to the south and is in good agreement with geological models as well as drilling information. The models for these profiles indicate that the thickness of sediments would further increase southwards into the Gulf of Kutch. Significant findings of the present study indicate 1) the hypocentre region of the earthquake is devoid of fluids, 2) absence of melt (that is emplaced during rifting as suggested from the passive seismological studies) in the lower crust and 3) a low resistive zone in the depth range of 5-20 km. The present MT study rules out fluidsand melt (magma) as the causative factors that triggered the Bhuj earthquake. The estimated porosity value of 0.02% will explain 100-500 ohm·m resistivity values observed in the lower crust. Based on the seismic velocities and geochemical studies, presence of garnet is inferred. The lower crust consists of basalts - probably generated by partial melting of metasomatised garnet peridotite at deeper depths in the lithosphere - and their composition might be modified by reaction with the spinel peridotites.

Crustal structure and reflectivity of the Swiss Alps from 3-Dimensional seismic modeling. 2. Penninic Nappes

Surface geological mapping, laboratory measurements of rock properties, and seismic reflection data are integrated through three-dimensional seismic modeling to determine the likely cause of upper crustal reflections and to elucidate the deep structure of the Penninic Alps in eastern Switzerland. Results indicate that the principal upper crustal reflections recorded on the south end of Swiss seismic line NFP20-EAST can be explained by the subsurface geometry of stacked basement nappes. In addition, modeling results provide improvements to structural maps based solely on surface trends and suggest the presence of previously unrecognized rock units in the subsurface. Construction of the initial model is based upon extrapolation of plunging surface. structures; velocities and densities are established by laboratory measurements of corresponding rock units. Iterative modification produces a best fit model that refines the definition of the subsurface geometry of major structures. We conclude that most reflections from the upper 20 km can be ascribed to the presence of sedimentary cover rocks (especially carbonates) and ophiolites juxtaposed against crystalline basement nappes. Thus, in this area, reflections appear to be principally due to first-order lithologic contrasts. This study also demonstrates not only the importance of three-dimensional effects (sideswipe) in interpreting seismic data, but also that these effects can be considered quantitatively through three-dimensional modeling.

Upper crustal structure of Madeira Island revealed from ambient noise tomography

We present the first image of the Madeira upper crustal structure, using ambient seismic noise tomography. 16 months of ambient noise, recorded in a dense network of 26 seismometers deployed across Madeira, allowed reconstructing Rayleigh wave Green's functions between receivers. Dispersion analysis was performed in the short period band from 1.0 to 4.0 s. Group velocity measurements were regionalized to obtain 20 tomographic images, with a lateral resolution of 2.0 km in central Madeira. Afterwards, the dispersion curves, extracted from each cell of the 2D group velocity maps, were inverted as a function of depth to obtain a 3D shear wave velocity model of the upper crust, from the surface to a depth of 2.0 km. The obtained 3D velocity model reveals features throughout the island that correlates well with surface geology and island evolution. (C) 2015 Elsevier B.V. All rights reserved.

Upper crustal structure of Madeira Island revealed from ambient noise tomography

We present the first image of the Madeira upper crustal structure, using ambient seismic noise tomography. 16 months of ambient noise, recorded in a dense network of 26 seismometers deployed across Madeira, allowed reconstructing Rayleigh wave Green's functions between receivers. Dispersion analysis was performed in the short period band from 1.0 to 4.0 s. Group velocity measurements were regionalized to obtain 20 tomographic images, with a lateral resolution of 2.0 km in central Madeira. Afterwards, the dispersion curves, extracted from each cell of the 2D group velocity maps, were inverted as a function of depth to obtain a 3D shear wave velocity model of the upper crust, from the surface to a depth of 2.0 km. The obtained 3D velocity model reveals features throughout the island that correlates well with surface geology and island evolution. (C) 2015 Elsevier B.V. All rights reserved.

Gravity Modelling of the Crustal Structure Below Bavali Shear Zone And Adjacent Major Plutons Of Northern Kerala

The primary aim of the present study is to acquire a large amount of gravity data, to prepare gravity maps and interpret the data in terms of crustal structure below the Bavali shear zone and adjacent regions of northern Kerala. The gravity modeling is basically a tool to obtain knowledge of the subsurface extension of the exposed geological units and their structural relationship with the surroundings. The study is expected to throw light on the nature of the shear zone, crustal configuration below the high-grade granulite terrain and the tectonics operating during geological times in the region. The Bavali shear is manifested in the gravity profiles by a steep gravity gradient. The gravity models indicate that the Bavali shear coincides with steep plane that separates two contrasting crustal densities extending beyond a depth of 30 km possibly down to Moho, justifying it to be a Mantle fault. It is difficult to construct a generalized model of crustal evolution in terms of its varied manifestations using only the gravity data. However, the data constrains several aspects of crustal evolution and provides insights into some of the major events.

Genetic algorithm inversion of the average 1D crustal structure using local and regional earthquakes

Knowing the best 1D model of the crustal and upper mantle structure is useful not only for routine hypocenter determination, but also for linearized joint inversions of hypocenters and 3D crustal structure, where a good choice of the initial model can be very important. Here, we tested the combination of a simple GA inversion with the widely used HYPO71 program to find the best three-layer model (upper crust, lower crust, and upper mantle) by minimizing the overall P- and S-arrival residuals, using local and regional earthquakes in two areas of the Brazilian shield. Results from the Tocantins Province (Central Brazil) and the southern border of the Sao Francisco craton (SE Brazil) indicated an average crustal thickness of 38 and 43 km, respectively, consistent with previous estimates from receiver functions and seismic refraction lines. The GA + HYPO71 inversion produced correct Vp/Vs ratios (1.73 and 1.71, respectively), as expected from Wadati diagrams. Tests with synthetic data showed that the method is robust for the crustal thickness, Pn velocity, and Vp/Vs ratio when using events with distance up to about 400 km, despite the small number of events available (7 and 22, respectively). The velocities of the upper and lower crusts, however, are less well constrained. Interestingly, in the Tocantins Province, the GA + HYPO71 inversion showed a secondary solution (local minimum) for the average crustal thickness, besides the global minimum solution, which was caused by the existence of two distinct domains in the Central Brazil with very different crustal thicknesses. (C) 2010 Elsevier Ltd. All rights reserved.

Gravity Data and Crustal Structure of North-Central Brazil

The Bouguer gravity anomaly of the northwest Ceará state in north-central Brazil was separated into its regional and residual components which were interpreted separately. By assuming that the sources of the regional anomalies are the depth variations of the crust-mantle interface, the mapping of these variations permited identifying crustal thickening zones which may be related to regional structures. The gravity residual sources coincide with occurrences of high-grade rocks (granulites) associated to medium-grade gneisses. Besides, the major strike-slip zones present significant signatures in the gravity data. This geophysical interpretation is compatible with the interpretation that the tectonic framework of the area is related to two crustal blocks conjoined by an A-type suture. The blocks are displaced along an oblique ramp with dextral movement, which played an important role in uplifting high-grade rocks from the lower crust to upper crustal levels. The suture zone corresponds to an imbricated compressive system dipping to the east and complicated by late dextral strike-slip shear zones.

Softening the Lower Crust: Modes of Syn-Transport Transposition Around and Adjacent to a Deep Crustal Granulite Nappe, Parry Sound Domain, Grenville Province, Ontario, Canada

The Parry Sound domain is a granulite nappe-stack transported cratonward during reactivation of the ductile lower and middle crust in the late convergence of the Mesoproterozoic Grenville orogeny. Field observations suggest the following with respect to the ductile sheath: (1) Formation of a carapace of transposed amphibolite facies gneiss derived from and enveloping the western extremity of the Parry Sound domain and separating it from high-strain gneiss of adjacent allochthons. This ductile sheath formed dynamically around the moving granulite nappe through the development of systems of progressively linked shear zones. (2) Transposition initiated by hydration (amphibolization) of granulite facies gneiss by introduction of fluid along cracks accompanying pegmatite emplacement. Shear zones nucleated along pegmatite margins and subsequently linked and rotated. The source of the pegmatites was most likely subjacent migmatitic and pegmatite-rich units or units over which Parry Sound domain was transported. Comparison of gneisses of the ductile sheath with high-strain layered gneiss of adjacent allochthons show the mode of transposition of penetratively layered gneiss depended on whether or not the gneiss protoliths were amphibolite or granulite facies tectonites before initiation of transposition, resulting in, e.g., folding before shearing, no folding before shearing, respectively. Meter-scale truncation along high-strain gradients at the margins of both types of transposition-related shear zones observed within and marginal to Parry Sound domain mimic features at kilometer scales, implying that apparent truncation by transposition originating in a manner similar to the ductile sheath may be a common feature of deep crustal ductile reworking. Citation: Culshaw, N., C. Gerbi, and J. Marsh (2010), Softening the lower crust: Modes of syn-transport transposition around and adjacent to a deep crustal granulite nappe, Parry Sound domain, Grenville Province, Ontario, Canada, Tectonics, 29, TC5013, doi:10.1029/2009TC002537.

Upper crustal structure of Madeira Island revealed from ambient noise tomography

We present the first image of the Madeira upper crustal structure, using ambient seismic noise tomography. 16 months of ambient noise, recorded in a dense network of 26 seismometers deployed across Madeira, allowed reconstructing Rayleigh wave Green's functions between receivers. Dispersion analysis was performed in the short period band from 1.0 to 4.0 s. Group velocity measurements were regionalized to obtain 2D tomographic images, with a lateral resolution of 2.0 km in central Madeira. Afterwards, the dispersion curves, extracted from each cell of the 2D group velocity maps, were inverted as a function of depth to obtain a 3D shear wave velocity model of the upper crust, from the surface to a depth of 2.0 km. The obtained 3D velocity model reveals features throughout the island that correlates well with surface geology and island evolution.

Response of a multi-domain continental margin to compression: Study from seismic reflection-refraction and numerical modelling in the Tagus Abyssal Plain

The effects of the Miocene through Present compression in the Tagus Abyssal Plain are mapped using the most up to date available to scientific community multi-channel seismic reflection and refraction data. Correlation of the rift basin fault pattern with the deep crustal structure is presented along seismic line IAM-5. Four structural domains were recognized. In the oceanic realm mild deformation concentrates in Domain I adjacent to the Tore-Madeira Rise. Domain 2 is characterized by the absence of shortening structures, except near the ocean-continent transition (OCT), implying that Miocene deformation did not propagate into the Abyssal Plain, In Domain 3 we distinguish three sub-domains: Sub-domain 3A which coincides with the OCT, Sub-domain 3B which is a highly deformed adjacent continental segment, and Sub-domain 3C. The Miocene tectonic inversion is mainly accommodated in Domain 3 by oceanwards directed thrusting at the ocean-continent transition and continentwards on the continental slope. Domain 4 corresponds to the non-rifted continental margin where only minor extensional and shortening deformation structures are observed. Finite element numerical models address the response of the various domains to the Miocene compression, emphasizing the long-wavelength differential vertical movements and the role of possible rheologic contrasts. The concentration of the Miocene deformation in the transitional zone (TC), which is the addition of Sub-domain 3A and part of 3B, is a result of two main factors: (1) focusing of compression in an already stressed region due to plate curvature and sediment loading; and (2) theological weakening. We estimate that the frictional strength in the TC is reduced in 30% relative to the surrounding regions. A model of compressive deformation propagation by means of horizontal impingement of the middle continental crust rift wedge and horizontal shearing on serpentinized mantle in the oceanic realm is presented. This model is consistent with both the geological interpretation of seismic data and the results of numerical modelling.