Southwest Barents Sea rift basin evolution: comparing results from backstripping and time‐forward modelling

We present tectonic models of progressive basin formation in the south‐west Barents Sea derived as part of the PETROBAR project (Petroleum‐related studies of the Barents Sea region). The basin architecture developed as a multi‐stage rift preceding the creation of the sheared/transtensional margin conjugate to NE Greenland. N‐ to NNE‐striking basins, with sediment thicknesses in places exceeding 15 km, are separated by basement highs. We use two basin analysis approaches, BMT? backstripping and TecMod?time‐forward modelling, to determine stretching factors through time along the profile PETROBAR‐07. This 550 km‐long profile derived from wide‐angle reflection/refraction seismic data acquired in 2007, coincident with deep multichannel seismic reflection data. Detailed stratigraphic analysis of the reflection profile, in concert with a dense grid of 2D profiles tied to wells, provides timing and water depth constraint for the models. Velocity analysis of the wide‐angle data provides constraint on the cumulative crustal stretching. The north‐west trending cross‐section extends from continental craton, at the Varanger Peninsula, to within 16 km of the interpreted continent–ocean boundary. Rifting along the profile was episodic, with four distinct phases of basin formation during the Carboniferous, the Late Permian–Triassic, the Late Jurassic–Early Cretaceous and the Late Cretaceous–Eocene. Collectively, the basins exhibit a general trend of younging, narrowing, and deepening oceanward, suggesting a gradual focusing of rifting prior to final breakup. Cumulative stretching factors derived from BMT and TecMod correlate well with observed crustal thinning, and the two models provide uncertainty bounds for stretching factors for the separate rift phases. In contrast to orthogonally rifted margins, stretching is relatively minor immediately prior to transform breakup, with greater stretching occurring during earlier rift phases.     

Extensional tectonics during the Tyrrhenian back‐arc basin formation and a new morpho‐tectonic map

We present a new tectonic map focused upon the extensional style accompanying the formation of the Tyrrhenian back‐arc basin. Our basin‐wide analysis synthetizes the interpretation of vintage multichannel and single‐channel seismic profiles, integrated with modern seismic images, P‐wave velocity models, and high‐resolution morpho‐bathymetric data. Four distinct evolutionary phases of the Tyrrhenian back‐arc basin opening are further constrained, redefining the initial opening to Langhian/Serravallian time. Listric and planar normal faults and their conjugates bound a series of horst and graben, half‐graben and triangular basins. Distribution of extensional faults, active throughout the basin since Middle Miocene, allows us to define an arrangement of faults in the northern/central Tyrrhenian mainly related to a pure shear which evolved to a simple shear opening. At depth, faults accommodate over a Ductile‐Brittle Transitional zone cut by a low‐angle detachment fault. In the southern Tyrrhenian, normal, inverse and transcurrent faults appear to be related to a large shear zone located along the continental margin of the northern Sicily. Extensional style variation throughout the back‐arc basin combined with wide‐angle seismic velocity models allows to explore the relationships between shallow deformation, faults distribution throughout the basin, and crustal‐scale processes as thinning and exhumation.     

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.     

Spatial distribution of sediment storage types and quantification of valley fill deposits in an alpine basin, Reintal, Bavarian Alps, Germany

Spatial patterns of sediment storage types and associated volumes using a novel approach for quantifying valley fill deposits are presented for a small alpine catchment (17 km²) in the Bavarian Alps. The different sediment storage types were analysed with respect to geomorphic coupling and sediment flux activity. The most landforms in the valley in terms of surface area were found to be talus slopes (sheets and cones) followed by rockfall deposits and alluvial fans and plains. More than two-thirds of the talus slopes are relict landforms, completely decoupled from the geomorphic system. Notable sediment transport is limited to avalanche tracks, debris flows, and along floodplains. Sediment volumes were calculated using a combination of polynomial functions of cross sections, seismic refraction, and GIS modelling. A total of, 66 seismic refraction profiles were carried out throughout the valley for a more precise determination of sediment thicknesses and to check the bedrock data generated from geomorphometric analysis. We calculated the overall sediment volume of the valley fill deposits to be 0.07 km³. This corresponds to a mean sediment thickness of 23.3 m. The seismic refraction data showed that large floodplains and sedimentation areas, which have been developed through damming effects from large rockfalls, are in general characterised by shallow sediment thicknesses (<20 m). By contrast, the thickness of several talus slopes is more than twice as much. For some locations (e.g., narrow sections of valley), the polynomial-generated cross sections resulted in overestimations of up to one order of magnitude; whereas in sections with a moderate valley shape, the modelled cross sections are in good accordance with the obtained seismic data. For the quantification of valley fill deposits, a combined application of bedrock data derived from polynomials and geophysical prospecting is highly recommended.     

A crustal seismic velocity model for the UK, Ireland and surrounding seas

A regional model of the 3-D variation in seismic P -wave velocity structure in the crust of NW Europe has been compiled from wide-angle reflection/refraction profiles. Along each 2-D profile a velocity–depth function has been digitised at 5 km intervals. These 1-D velocity functions were mapped into three dimensions using ordinary kriging with weights determined to minimise the difference between digitised and interpolated values. An analysis of variograms of the digitised data suggested a radial isotropic weighting scheme was most appropriate. Horizontal dimensions of the model cells are optimised at 40 × 40 km and the vertical dimension at 1 km. The resulting model provides a higher resolution image of the 3-D variation in seismic velocity structure of the UK, Ireland and surrounding areas than existing models. The construction of the model through kriging allows the uncertainty in the velocity structure to be assessed. This uncertainty indicates the high density of data required to confidently interpolate the crustal velocity structure, and shows that for this region the velocity is poorly constrained for large areas away from the input data.     

Upper mantle structure beneath the Mid-Atlantic Ridge north of the Azores based on observations of compressional wave

Summary. Four seismic refraction profiles have been interpreted which serve to indicate the structure of the lithosphere near the Mid-Atlantic Ridge close to the Azores. An east–west profile which crosses the ridge axis yields a crustal structure. Although energy is propagated across the ridge axis within the crust the axial region marks a clear barrier to propagation within the mantle. A profile parallel to the axis (4 my isochron) shows, below a 7.6 km/s layer, a low-velocity zone underlain by an 8.3 km/s refractor 9 km below the sea bed. On profies normal to the ridge axis higher velocities, which are observed on lines shot towards the ridge, can be attributed to this refractor if it has a dip of several degrees away from the ridge. On another profile parallel to the axis (9 my isochron) a velocity of about 8.3 km/s is only found to exist much deeper at about 30 km depth. These observations are interpreted in the light of seismic refraction results recently obtained by Lewis & Snydsman and of quantitative petrological models, such as that of Bottinga & Allègre. A velocity model based on Bottinga & Allègre's model allows us to understand our results qualitatively. In particular the two 8.3 km/s refractors at 9 and 30 km depth correspond to two different residual peridotite layers. The upper layer contains 1.5–2 wt per cent water and as the lithosphere moves away from the ridge axis the temperature in this layer becomes low enough to start hydration reactions. These cause the low-velocity zone observed at 4 my and the total disappearance of the shallow level refractor before 9 my.     

New seismic data support Cenozoic rifting in George VI Sound, Antarctic Peninsula

New multifold seismic-reflection and wide-angle reflection/refraction data across George VI Sound, Antarctic Peninsula, show the presence of graben and horst structures indicating an extensional origin. The data suggest that rocks of an accretionary complex and fore-arc basin underlie the Sound and are in faulted contact along its eastern boundary with volcanic and plutonic rocks of the associated Mesozoic arc of western Palmer Land. A cover of possible syn- and postglacial Cenozoic deposits drapes the structures. The combination of new seismic, synthetic-aperture radar and previously acquired data suggests subduction-related rifting in the Sound was segmented, with opening in the south predominately by normal extension whilst in the north, dextral transtension predominated.     

Migration - why doesn't it work for deep continental data?

Summary. Conventional migration of deep seismic reflection data often produces disappointingly poor results even when the original unmigrated data are of high quality and are relatively noise free. Surprisingly, deep data often subjectively appear to be best migrated at velocities which are up to 50% less than appropriate interval velocities derived from crustal refraction experiments or directly from stacking velocities. The explanation for this behaviour is that near surface features distort and attenuate the seismic wave field and produce apparent discontinuities in deep reflections. Since these discontinuities are spurious they are not associated with the appropriate large diffractions which real discontinuities at depth would produce. During the process of migration reflections are invented in order to cancel out the missing diffractions thereby producing a "smiley" section which appears overmigrated. Since the lateral extent of individual smiles increases with increasing migration velocity, increasing two-way-time, and increasing seismic wavelength, the effect" is almost unnoticeable on conventional shallow seismic data but is overwhelming for deep crustal data.     

Deep seismic soundings along Hirapur-Mandla profile, central India

Summary. The crustal depth section along Hirapur-Mandla profile has been computed in two steps from Deep Seismic Sounding (DSS) data. The shallow section up to the crystalline basement is derived by inverting first arrival refraction travel times. The upper Vindhyan sediments (velocity 4.5 km s⁻¹) have a maximum thickness of about 1.5 km at Bakshaho. The lower Vindhyan sediments (velocity 5.4 km s⁻¹) were deposited north of Narmada-Son lineament between Katangi and Narsinghgarh in a graben developed in crystalline basement. The thickness of the lower Vindhyans increases from north to south towards Katangi and the depth to the basement reaches 5.5 km near Jabera. The depth to the Moho boundary varies from 39.5 km near Tikaria to 45 km at Narsinghgarh. The narrow block between Katangi and Jabalpur forms a horst feature which represents the Narmada-Son lineament forming the southern boundary of the Vindhyan basin. Two-dimensional ray tracing was performed generating travel time curves from various shot points which were matched with observed travel time data.     

Pre-drilling reflection survey of the Black Forest, SW Germany

Summary. The unified seismic exploration program, consisting of 345 km of deep reflection profiling, a 200 km refraction profile, an expanding spread profile and near-surface high resolution reflection meaasurements, revealed a strongly differentiated crust beneath the Black Forest. The highly reflective lower crust contains numerous horizontal and dipping reflectors at depths of 13-14 km down to the crust-mantle boundary (Moho). The Moho appears as a flat horizontal first order discontinuity at a relatively shallow level of 25–27 km above a transparent upper mantle. From modelling of synthetic near-vertical and wide-angle seismograms using the reflectivity method the lower crust is supposed to be composed of laminae with an average thickness of about 100 m and velocity differences of greater than 10% increasing from top to bottom. The upper crust is characterised by mostly dipping reflectors, associated with bivergent underthrusting and accretion tectonics of Variscan age and with extensional faults of Mesozoic age. A bright spot at 9.5 km depth is characterised by low velocity material suggesting a fluid trap. It appears on all of the three profiles in the centre of the intersection region. The upper crust seems to be decoupled from the lowest crust by a relatively transparent zone which is' also identified as a low-velocity zone. This low velocity channel is situated directly above the laminated lower crust. The laminae in the Rhinegraben area are displaced vertically to greater depths indicating an origin before Tertiary rift formation and a subsidence of the whole graben wedge.     

Geological and geophysical investigation of the Mid-Cayman Spreading Centre: seismic velocity measurements and implications for the constitution of layer 3

Summary . In this paper we present laboratory measurements of compressional and shear wave velocities of a diverse suite of gabbroic rocks collected from the walls of the Mid-Cayman Spreading Centre with DSRV Alvin. The degree of deformation and alteration affecting these gabbros is quite variable, and we believe that they are typical of plutonic rocks emplaced at shallow levels (upper portion of seismic layer 3 and shallower) of the oceanic crust. The compositional and textural variations are reflected in the wide range of laboratory velocities which span most of the range of seismic velocities reported for oceanic and ophiolite rock samples including basalts, gabbros, ultramafics, and their altered derivatives. Based upon the laboratory velocities and the geological setting of the Mid-Cayman gabbros, it is argued that no unique lithology, except anhydrous peridotite, can be unequivocally identified in the oceanic lithosphere from seismic velocity data alone. Furthermore, these data allow for the possibility of considerable lithologic heterogeneity within portions of the oceanic crust at the scale of a few centimetres to a few hundred metres. Such heterogeneities would go unrecognized because seismic refraction studies mask these variations resulting in a picture of apparent uniformity.     

Inversion of seismic refraction data in planar dipping structure

Summary. A new method is presented for the direct inversion of seismic refraction data in dipping planar structure. Three recording geometries, each consisting of two common-shot profiles, are considered: reversed, split, and roll-along profiles. Inversion is achieved via slant stacking the common-shot wavefield to obtain a delay time–slowness (tau– p ) wavefield. The tau– p curves from two shotpoints describing the critical raypath of refracted and post-critically reflected arrivals are automatically picked using coherency measurements and the two curves are jointly used to calculate velocity and dip of isovelocity lines iteratively, thereby obtaining the final two-dimensional velocity model.
This procedure has been successfully applied to synthetic seismograms calculated for a dipping structure and to field data from central California. The results indicate that direct inversion of closely-spaced refraction/wide-aperture reflection data can practically be achieved in laterally inhomogeneous structures.     

Basement depth and sedimentary velocity structure in the northern Arabian platform, eastern Syria

Basement depth in the Arabian plate beneath eastern Syria is found to be much deeper than previously supposed. Deep-seated faulting in the Euphrates fault system is also documented. Data from a detailed 300 km long reversed refraction profile, with offsets up to 54 km, are analysed and interpreted, yielding a velocity model for the upper 9 km of continental crust. The interpretation integrates the refraction data with seismic-reflection profiles, well logs and potential field data, such that the results are consistent with all available information. A model of sedimentary thicknesses and seismic velocities throughout the region is established. Basement depth on the north side of the Euphrates is interpreted to be around 6 km, whilst south of the Euphrates basement depth is at least 8.5 km. Consequently, the potentially hydrocarbon-rich pre-Mesozoic section is shown, in places, to be at least 7 km thick. The dramatic difference in basement depth on adjacent sides of the Euphrates graben system may suggest that the Euphrates system is a suture/shear zone, possibly inherited from Late Proterozoic accretion of the Arabian plate. Gravity modelling across the southeast Euphrates system tends to support this hypothesis. Incorporation of previous results allows us to establish the first-order trends in basement depth throughout Syria     

Seismic anisotropy within the uppermost mantle of southern Germany

This paper presents an updated interpretation of seismic anisotropy within the uppermost mantle of southern Germany. The dense network of reversed and crossing refraction profiles in this area made it possible to observe almost 900 traveltimes of the P_n phase that could be effectively used in a time-term analysis to determine horizontal velocity distribution immediately below the Moho. For 12 crossing profiles, amplitude ratios of the P_n phase compared to the dominant crustal phase were utilized to resolve azimuthally dependent velocity gradients with depth. A P -wave anisotropy of 3–4 per cent in a horizontal plane immediately below the Moho at a depth of 30 km, increasing to 11 per cent at a depth of 40 km, was determined. For the axis of the highest velocity of about 8.03 km s⁻¹ at a depth of 30 km a direction of N31°F was obtained. The azimuthal dependence of the observed P_n amplitude is explained by an azimuth-dependent sub-Moho velocity gradient decreasing from 0.06 s⁻¹ in the fast direction to 0 s⁻¹ in the slow direction of horizontal P -wave velocity. From the seismic results in this study a petrological model suggesting a change of modal composition and percentage of oriented olivine with depth was derived.     

Receiver functions at the stations of the German Regional Seismic Network (GRSN)

Several years of broad-band teleseismic data from the GRSN stations have been analysed for crustal structure using P -to- S converted waves at the crustal discontinuities. An inversion technique was developed which applies the Thomson-Haskell formalism for plane waves without slowness integration. The main phases observed are Moho conversions, their multiples in the crust, and conversions at the base of the sediments. The crustal thickness derived from these data is in good agreement with results from other studies. For the Gräfenberg stations, we have made a more detailed comparison of our model with a previously published model obtained from refraction seismic experiments. The refraction seismic model contains boundaries with strong velocity contrasts and a significant low-velocity zone, resulting in teleseismic waveforms that are too complicated as compared to the observed simple waveforms. The comparison suggests that a significant low-velocity zone is not required and that internal crustal boundaries are rather smooth.     

Combined seismic reflection and refraction profiling in southwest Germany - detailed velocity mapping by the refraction survey

Summary. Within the framework of site survey studies of the Deep Drilling Program of the Fedral Republic of Germany (KTB), coincident deep-seismic reflection and refraction experiments in the Black Forest, southwest Germany, were carried out. The simultaneous interpretation of the reflection and the refraction data reveals in particular both a strong velocity reduction in the upper crust and a laterally varying laminated structure of the lower curst. Additional refraction lines result in a three-dimensional crustal model which shows two distinct crustal types of different seismic properties. These crustal types seem to correlate with the major geologic units of Southwest Germany. Variations of Poisson's ratio derived from clearly recorded shear wave data show a similar trend.     

Physical properties and structure of the lower crust revealed by one- and two-dimensional modelling

Summary. Data from a refraction and a reflection seismic survey in the Black Forest, southwest Germany, are used for extensive one- and two-dimensional modelling. The data are available along approximately the same line, and therefore the same piece of crust is probed by two seismic methods. We utilize this favorable circumstance for detailed model calculations concerning both data sets. Lower crustal properties vary on the scale of a wavelength and thus full solutions of the elastic equations are required: the Reflectivity Method for the evaluation of refraction seismograms and numerical solutions of the acoustic wave equation for the reflection response. Details of the geometry and physical properties of the lamination are derived. Vertical layering on a scale of 100 m is found; horizontal extent of reflecting elements is in the range of a few hundreds of meters; rocks with velocities between 5.6 and 7.2 km/s constitute the lower crust.     

Operations and main results of the ECORS project in France

Summary. The ECORS project was launched in 1983 with the aim of studying fundamental mechanisms of geodynamics in France. The ECORS deep seismic profiles have concentrated on a few structures of major significance: outer zone of the Variscan orogen in northern France, the basin formed during the Bay of Biscay's opening and transects of recent orogenic ranges (Pyrenees and Alps). The seismic profiles have been carried out with all the available modern techniques of industry and completed wherever possible by additional geophysical surveys (magnetism, gravity, MT, wide-angle and refraction seismics) and geological surveys. The first results already shed new light on major geodynamic phenomena such as variations in the frontal Variscan detachment, lower crust formation, crustal behaviour during orogenesis and variations in the formation of cratonic basins.     

Detailed S-wave structure in the Dublin Basin and its northern margin

Summary. The seismic structure has been measured to a depth of about 3 km along a 30 km seismic profile in east central Ireland. This profile is unusual in that it is the S -wave velocity—depth structure that has been measured to a degree of precision more normally associated with P -wave results. One reason for this is that the sources used were quarry blasts which generated strong S -waves and short-period surface waves but rather weak P -waves.
The results show a layer of Carboniferous limestone with shear velocity 2.65 km⁻¹ s overlying a layer with a velocity of 3.06 km s⁻¹. This second layer was interpreted as Lower Palaeozoic strata (Silurian/Ordovician) since this velocity was evident in an inlier seen at the surface at the northern end of the line. A third refraction horizon, shear velocity 3.45 km s⁻¹ and displaying a basinal structure, was also recognized. This may be Cambrian or Precambrian basement.     

Crustal structure from an OBS survey of the Nootka fault zone off western Canada

Summary. The Nootka fault zone is the boundary between the small Explorer and Juan de Fuca plates which are situated between the America and Pacific plates off western Canada. To investigate the crustal structure in the region, three explosive/large airgun refraction lines were shot into three ocean bottom seismometers (OBSs) with three-component geophone assemblies. In this phase of the study, P -wave velocity—depth models are interpreted by comparison of the travel time and amplitude characteristics of the observed data with theoretical seismograms computed using a WKBJ algorithm. The interpretation gives relatively consistent results for the upper crust. However, the structure of the lower crust is significantly different among the various profiles. Upper mantle velocities range from 7.5 to 8.3 kms⁻¹ and the sub-bottom crustal thickness vanes from 6.4 to 11 km. Nevertheless, these seismic models are consistent in general terms with oceanic crustal models represented by ophiolite complexes. Some aspects of the differences among profiles can be explained by consideration of a recent tectonic model for the development of the fault zone. This requires, within a 1 Myr time interval, variations in the process of crustal formation at the ridge, crustal 'maturing', or both. The abnormally thick crust near a spreading centre may result in part from the complex interaction of the Juan de Fuca and Explorer plates with the larger and older America and Pacific plates. Upper mantle velocity variations are consistent with the concept of velocity anisotropy. The different record sections show that seismic energy is attenuated for ray paths traversing the Nootka fault zone.