Seismic Reflection Experiment Research Articles

Integration of in‐mine seismic and GPR surveys to gain advanced knowledge of Bushveld Complex orebodies

AbstractImproving the exploration of deep‐seated mineral deposits and assessing the stability of the mine pillars require that geophysical techniques are deployed in a fast and cost‐effective manner with minimal environmental impact. This research presents results from in‐mine reflection seismic experiments and a ground penetrating radar (GPR) survey conducted at the Maseve platinum mine, South Africa. The research aims to develop and implement methods to image platinum group metal (PGM) deposits and geological structures near mine tunnels and assess the stability of pillars. The seismic experiments were conducted using a sledgehammer source (10 lb), conventional cabled geophones (14 Hz), and a landstreamer with 4.5 Hz vertical component geophones. The GPR survey was conducted using a Noggin 500 GPR system with 500 MHz centre frequency. An image of the underlying orebody and geological structures down to 100 m from the mine tunnel floor (∼500 m below ground surface) was produced. We correlated the coherent reflections beneath the tunnel floor with a known Upper Group (UG2) PGM orebody. The final seismic section shows that the UG2 mineralisation is dissected by near‐vertical dykes, faults and fractures. These structures, faults in particular, are interpreted to have been active post‐mineralisation, implying that they may have contributed to the current complex geometry of the deposit. Four GPR profiles were collected around a stability pillar adjacent to the seismic lines. The radargram sections were processed to improve the signal‐to‐noise ratio (S/N). The results show different patterns of fracturing and stress‐induced structures. These fractures were shown to be sub‐vertical and, possibly, constitute complex micro‐structures within the pillar, which could compromise the pillar stability and integrity. The study demonstrates that in‐mine seismic and GPR surveys can be cost‐effective and valuable for mineral exploration.

Cretaceous magmatic underplating and delamination beneath continental SE Brazil and their tectonic implications: Evidence from the PABBRISE wide-angle reflection and refraction seismic profile

Controlled-source seismology provides robust methods to obtain information on the velocity structure of the crust and uppermost mantle. We present the first systematic wide-angle reflection and refraction (WARR) seismic experiment for continental southeastern Brazil: the NW-SE-oriented and ca. 700 km-long PABBRISE (PAraná Basin, Brasília and RIbeira orogens Seismic Experiment) profile. Using ray tracing, we obtained a two-dimensional P-wave seismic velocity model of the crust and uppermost mantle. Our model shows that Moho depths vary between ca. 31 and 45 km, with crustal thinning towards the continental margin. The modeled velocity structure sheds new light on the present-day lithospheric morphology of the region, with the most striking features noted within the lower crust and uppermost mantle. We show that beneath the Paraná Basin, parts of the lower crust and uppermost mantle were pervasively intruded due to magmatic underplating. We modeled a ca. 200 km-long and 11 km-thick (maximum thickness), high-velocity (7.20–7.80 km/s) lower crustal (HVLC) body in this region. Outside the Paraná Basin, near the continental margin, comparatively thinner crust (∼31 km) with low mean P-wave velocity (6.20 km/s) and a low velocity (<7.90 km/s) uppermost mantle probably are the result of delamination. Magmatic underplating and delamination can be linked with the breakup of Pangea and the opening of the South Atlantic Ocean during the Cretaceous.

Toward a Practical Appraisal for Waveform Tomography of Band- and Offset-Limited Marine Seismic Data

We present a generalized workflow to retrieve high-resolution P-wave velocity ( <formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"> <tex>$V_{ P}$</tex> </formula> ) models of complex Earth's subsurface structures from traditional marine near-vertical seismic reflection experiments. These records have typically offsets too short to map refraction phases and lack low-frequency information. The workflow is composed of three steps: 1) downward continuation (DC) of seismic records to the seafloor to recover diving wave information; 2) travel-time tomography (TTT) of first arrivals obtained from DC data, to retrieve a kinematically correct model; and 3) full-waveform inversion (FWI) of the original streamer dataset, starting with the model obtained with TTT and sequentially introducing higher wavenumber details into the model. We show that the TTT allows overcoming the issues associated with the nonlinearity intrinsic to FWI. We also disentangle envelope and phase from the waveform to choose the objective function most suitable for FWI. We assess the accuracy of initial models and predict the quality of the FWI results by quantifying the early arrival cycle skipping between original and simulated data. The efficiency of the workflow is tested with a challenging synthetic target model, containing vertical boundaries with strong velocity contrasts and velocity inversions embedded in a checkerboard-like pattern. We show that workflow steps 1) and 2) provide a TTT model that is not cycle skipped at the frequencies available in most marine seismic experiments and thus allow step 3) FWI to obtain high-resolution <formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"> <tex>$V_{ P}$</tex> </formula> models of the subsurface using band- and offset-limited field datasets, traditionally collected in marine airgun and streamer acquisitions.

Near-surface structure of the Sodankylä area in Finland, obtained by passive seismic interferometry

Abstract. Controlled-source seismic exploration surveys are not always possible in nature-protected areas. As an alternative, the application of passive seismic techniques in such areas can be proposed. In our study, we show results of passive seismic interferometry application for mapping the uppermost crust in the area of active mineral exploration in northern Finland. We utilize continuous seismic data acquired by the Sercel Unite wireless multichannel recording system along several profiles during XSoDEx (eXperiment of SOdankylä Deep Exploration) multidisciplinary geophysical project. The objective of XSoDEx was to obtain a structural image of the upper crust in the Sodankylä area of northern Finland in order to achieve a better understanding of the mineral system at depth. The key experiment of the project was a high-resolution seismic reflection experiment. In addition, continuous passive seismic data were acquired in parallel with reflection seismic data acquisition. Due to this, the length of passive data suitable for noise cross-correlation was limited from several hours to a couple of days. Analysis of the passive data demonstrated that dominating sources of ambient noise are non-stationary and have different origins across the XSoDEx study area. As the long data registration period and isotropic azimuthal distribution of noise sources are two major conditions for empirical Green function (EGF) extraction under the diffuse field approximation assumption, it was not possible to apply the conventional techniques of passive seismic interferometry. To find the way to obtain EGFs, we used numerical modelling in order to investigate properties of seismic noise originating from sources with different characteristics and propagating inside synthetic heterogeneous Earth models representing real geological conditions in the XSoDEx study area. The modelling demonstrated that scattering of ballistic waves on irregular shape heterogeneities, such as massive sulfides or mafic intrusions, could produce a diffused wavefield composed mainly of scattered surface waves. In our study, we show that this scattered wavefield can be used to retrieve reliable EGFs from short-term and non-stationary data using special techniques. One of the possible solutions is application of “signal-to-noise ratio stacking” (SNRS). The EGFs calculated for the XSoDEx profiles were inverted, in order to obtain S-wave velocity models down to the depth of 300 m. The obtained velocity models agree well with geological data and complement the results of reflection seismic data interpretation.

Sparse 3D reflection seismic survey for deep-targeting iron oxide deposits and their host rocks, Ludvika Mines, Sweden

Abstract. Many metallic mineral deposits have sufficient physical property contrasts, particularly density, to be detectable using seismic methods. These deposits are sometimes significant for our society and economic growth and can help to accelerate the energy transition towards decarbonization. However, their exploration at depth requires high-resolution and sensitive methods. Following a series of 2D seismic trials, a sparse, narrow source–receiver azimuth, 3D seismic survey was conducted in the Blötberget mine, in central Sweden, covering an area of approximately 6 km2 for deep-targeting iron oxide deposits and their host rock structures. The survey benefited from a collaborative work by putting together 1266 seismic recorders and a 32 t vibrator, generating 1056 shot points in a fixed geometry setup. Shots were fired at every 10 m where possible, and receivers were placed at every 10–20 m. Notable quality data were acquired despite the area being dominated by swampy places as well as by built-up roads and historical tailings. The data processing had to overcome these challenges for the static corrections and strong surface waves in particular. A tailored for hardrock setting and processing workflow was developed for handling such a dataset, where the use of mixed 2D and 3D refraction static corrections was relevant. The resulting seismic volume is rich in terms of reflectivity, with clear southeast-dipping reflections originating from the iron oxide deposits extending vertically and laterally at least 300 m beyond what was known from available boreholes. As a result, we estimate potential additional resources from the 3D reflection seismic experiment on the order of 10 Mt to be worth drilling for detailed assessments. The mineralization is crosscut by at least two major sets of northwest-dipping reflections interpreted to dominantly be normal faults and to be responsible for much of the lowland in the Blötberget area. Moreover, these post-mineralization faults likely control the current 3D geometry of the deposits. Curved and submerged reflections interpreted from folds or later intrusions are also observed, showing the geological complexity of the study area. The seismic survey also delineates the near-surface expression of a historical tailing as a by-product of refraction static corrections, demonstrating why 3D seismic data are so valuable for both mineral exploration and mine planning applications.

Effective utilization of seismic reflection technique with moderate cost in uranium exploration

ABSTRACTIn participation with numerous industrial partners, the Seismic Laboratory of the University of Saskatchewan has conducted a variety of active seismic reflection experiments; both on the west and east sides of the Athabasca Basin. Results of the investigations at Shea Creek, McArthur River and Keefe Lake illustrate that the seismic investigations deliver effective, highly relevant primary structural images of the subsurface, with resolution that no other geophysical technique can match. Correlation of similar seismic signatures, on several distant but inter‐related seismic sections, allowed spatial extension of promising exploration target zones previously unrecognized. Within the three‐dimensional seismic volume, comparable reflectivity patterns defined the complex areal distribution of mineralization‐related fault systems. Beyond these novel contributions, extended analysis of seismic signal attributes (amplitude and frequency), optical televiewer, and full‐wave sonic data offer detailed lithological characterization, including alteration zones, clay content, as well as porosity and fracture density information. Although these structural and geologically relevant anomalies are primary indicators of mineralization, presenting novel exploration advantages, the seismic method is still not a standard component of the Athabasca Basin exploration approach, due to the negative perception that ‘it is very expensive’. Comparing the costs of all geophysical techniques to the cost of a single logged drill hole illustrates that the results of a properly designed seismic data acquisition program not only leads to more effective planning of a drilling program, but also would lead to a much quicker recognition of the major mineralized zone(s), and a reduction in the number of required exploration boreholes. This integrated approach to exploration would then translate into a significant reduction of the total exploration expenditures. Unquestionably, the drilling of boreholes provides the most explicit, reliable information to a certain depth, but only within a very small area. Directly connecting the borehole information to seismic results extends the local reliable data; permitting reduction of the number of boreholes to create accurate two‐dimensional or three‐dimensional subsurface images and reduction of the expenditures of the total exploration program.

Passive processing of active nodal seismic data: estimation of &amp;lt;i&amp;gt;V&amp;lt;/i&amp;gt;&amp;lt;sub&amp;gt;P&amp;lt;/sub&amp;gt;∕&amp;lt;i&amp;gt;V&amp;lt;/i&amp;gt;&amp;lt;sub&amp;gt;S&amp;lt;/sub&amp;gt; ratios to characterize structure and hydrology of an alpine valley infill

Abstract. The advent of cable-free nodal arrays for conventional seismic reflection and refraction experiments is changing the acquisition style for active-source surveys. Instead of triggering short recording windows for each shot, the nodes are continuously recording over the entire acquisition period from the first to the last shot. The main benefit is a significant increase in geometrical and logistical flexibility. As a by-product, a significant amount of continuous data might also be collected. These data can be analyzed with passive seismic methods and therefore offer the possibility to complement subsurface characterization at marginal additional cost. We present data and results from a 2.4 km long active-source profile, which have recently been acquired in western Colorado (US) to characterize the structure and sedimentary infill of an over-deepened alpine valley. We show how the “leftover” passive data from the active-source acquisition can be processed towards a shear wave velocity model with seismic interferometry. The shear wave velocity model supports the structural interpretation of the active P-wave data, and the P-to-S-wave velocity ratio provides new insights into the nature and hydrological properties of the sedimentary infill. We discuss the benefits and limitations of our workflow and conclude with recommendations for the acquisition and processing of similar datasets.

Evaluation of the 2D reflection seismic method toward the exploration of thrust-controlled mineral deposits in southwestern Fujian Province, China

The southwestern region of the Fujian Province is one of the major ore districts in China. The current model states that mineral deposition is highly controlled by thrust structure, which means that there may be concealed deposits located deep within overlapping thrust areas. Reflection seismology, which has great depth penetration and higher resolution than other geophysical methods, has great potential to delineate complex structures and give some clues to mineralization. In 2015, an experimental 2D reflection seismic survey called “Fujian 2D” was conducted in this region. Data were acquired along a 13.8 km length, with a source interval of 60 m, and 691 identical receivers with an equal spacing of 20 m were used to record data for each source. Due to topographical restrictions caused by the source environment, the mass or position of some shots was changed. Despite the restrictions, the average fold number reached 64 for a 10 km distance along the middle of the survey line. During the data processing procedure, conventional technologies involving static correction, noise elimination, deconvolution, and iterative velocity analyses were applied. After the prestack time migration failed to obtain a high-quality imaging result, rugged prestack depth migration (PSDM) was introduced that resulted in a better quality image of the subsurface structure and which included near-surface parts of the thrusts. In addition, P- and S-wave velocities and density data were determined from two borehole cores. Forward modeling and imaging found that the Permian marble hosting the mineral deposits has lower velocity than the surrounding rocks, where contacts give rise to strong reflections. The final rugged PSDM also clearly delineated the thrust bodies and magma intrusion zones. Combining this forward modeling with the known geology of the investigated site, the Fujian 2D reflection seismic experiment demonstrates great potential for unveiling the main elements controlling mineral deposition, such as tectonic structure, stratigraphic contacts, and lithology. Our experimental results demonstrate that reflection seismology has a wide range of applications for future mineral exploration at greater depths.

Crustal structure of the Mendeleev Rise and the Chukchi Plateau (Arctic Ocean) along the Russian wide-angle and multichannel seismic reflection experiment “Arctic-2012”

We present a seismic and density model for the crust and the uppermost mantle of the Arctic Ocean off-shore Chukotka down to a 40 km depth along a 740-km long latitudinal (at ca. 77°N) “Arctic-2012” wide-angle/MCS profile. Joint seismic and gravity modeling indicates significant differences in the crustal velocity and density structure of the northeastern Vilkitsky Trough, the Mendeleev Rise, the Chukchi Basin, and the Chukchi Plateau.The Vilkitsky Trough and the Chukchi Basin have a thin crust (23 km and 18 km, correspondingly), 6–8 km thick sedimentary cover, 3–6 km thick upper/middle crust (with the smallest thickness of 3–4 km beneath the Chukchi Basin), and 9–10 km thick lower crust. The uppermost mantle of the Chukchi Basin has a high density (3.27–3.31 g/cm3) and a low velocity (Vp ∼ 7.8 km/s), which we explain by 5–10% serpentinization of mantle peridotite at a 22–35 km depth as a result of crustal hyperextension and seawater penetration.The Chukchi Plateau and the Mendeleev Rise have a thick crust (28–29 km and 33–34 km, correspondingly), underlain by a normal mantle (Vp ∼ 8.0 km/s). The Chukchi Plateau has a 2‐4 km thick sedimentary cover, a thick (15–18 km) upper/middle crust with low-Vp, low-density lenses interpreted as magmatic intrusions, and a 9–12 km thick lower crust. The Mendeleev Rise has a 3–7 km thick sedimentary cover (most of which is formed by metasediments with a possible presence of volcanic rocks), a 7–8 km thick upper/middle crust, and a thick (20 km) lower crust which includes a 3–4 km thick high-velocity (Vp ∼ 7.3 km/s) underplated magmatic material. The high density anomaly (at depths >35 km) below the Mendeleev Rise is interpreted as an eclogitic body in the upper mantle lithosphere.Seismic Vp and Vp/Vs structure of the crust along the “Arctic-2012” profile indicates its continental nature: a 3–18 km thick upper/middle crustal layer with Vp ∼ 6.0–6.8 km/s and Vp/Vs ∼ 1.70–1.73 typical of felsic-intermediate continental upper crust is present along the entire profile. Strong variability of the crustal structure along the profile reflects its significant modification by metamorphism and magmatism, possibly related to the High-Arctic Large Igneous Province and localized lithosphere extension beneath the Chukchi Basin.

The derivation of an anisotropic velocity model from a combined surface and borehole seismic survey in crystalline environment at the COSC-1 borehole, central Sweden

The Scandinavian Caledonides provide a well-preserved example of a Palaeozoic continent–continent collision, where surface geology in combination with geophysical data provides information about the geometry of parts of the Caledonian structure. The project COSC (Collisional Orogeny in the Scandinavian Caledonides) investigates the structure and physical conditions of the orogen units and the underlying basement with two approximately 2.5 km deep cored boreholes in western Jämtland, central Sweden. In 2014, the COSC-1 borehole was successfully drilled through a thick section of the Seve Nappe Complex. This tectonostratigraphic unit, mainly consisting of gneisses, belongs to the so-called Middle Allochthons and has been ductilely deformed and transported during the collisional orogeny. After the drilling, a major seismic survey was conducted in and around the COSC-1 borehole with the aim to recover findings on the structure around the borehole from core analysis and downhole logging. The survey comprised both seismic reflection and transmission experiments, and included zero-offset and multiazimuthal walkaway Vertical Seismic Profile (VSP) measurements, three long offset surface lines centred on the borehole, and a limited 3-D seismic survey. In this study, the data from the multiazimuthal walkaway VSP and the surface lines were used to derive detailed velocity models around the COSC-1 borehole by inverting the first-arrival traveltimes. The comparison of velocities from these tomography results with a velocity function calculated directly from the zero-offset VSP revealed clear differences in velocities for horizontally and vertically travelling waves. Therefore, an anisotropic VTI (transversely isotropic with vertical axis of symmetry) model was found that explains first-arrival traveltimes from both the surface and borehole seismic data. The model is described by a vertical P-wave velocity function derived from zero-offset VSP and the Thomsen parameters ε = 0.03 and δ = 0.3, estimated by laboratory studies and the analysis of the surface seismic and walkaway VSP data. This resulting anisotropic model provides the basis for further detailed geological and geophysical investigations in the direct vicinity of the borehole.

Using large dynamite shots to image the structure of the Moho from deep seismic reflection experiment between the Sichuan basin and Qinling orogen

The Qinling orogen was formed as a result of the collision between the North and South China blocks. The Qinling orogen represents the location at which the southern and northern parts of the Chinese mainland collided, and it’s also the intersection of the Central China orogen and the north-south tectonic belt. There is evidence of strong deformation in this orogen, and it has had a long and complex geological history. We investigated the structure of the Moho in the southern Qinling orogen using large dynamite shot imaging techniques. By integrating the analysis of the single-shot and the move-out corrections profile, we determined the structure of the Moho beneath the northern Dabashan thrust belt and the southern Qinling orogen, including the mantle suture beneath Fenghuang mountain. The Moho is divided into two parts by the mantle suture zone beneath Fenghuang mountain: (1) from Ziyang to Hanyin, the north-dipping Moho is at about 45–55 km depth and the depth increases rapidly; and (2) from Hanyin to Ningshan, the south-dipping Moho is at about 40–45 km depth and shallows slowly. The mantle suture is located beneath Fenghuang mountain, and the Moho overlaps at this location: the shallower Moho is connected to the northern part of China, and the deeper Moho is connected to the southern part. This may indicate that the lithosphere in the Sichuan basin subducts to the Qinling block and that the subduction frontier reaches at least as far as Fenghuang mountain.

PREDICTED SEDIMENTARY SECTION OF SUBGLACIAL LAKE VOSTOK

In early February 2012, the drill hole at the Vostok Station encountered theLakeVostokwater. This step is important to study the lake composition including possible microbial life and to model subglacial environments however, the next ambitious target of the Vostok Drilling Project is sampling of bottom sediments, which contain the unique record of ice sheet evolution and environmental changes in centralAntarcticafor millions of years. In this connection, the forecast of sedimentary succession based on existing geophysical data, study of mineral inclusions in the accretion ice cores and tectonic models is important task. Interpretation of Airborne geophysical data suggests thatLakeVostokis the part of spacious rift system, which exists at least from Cretaceous. Reflection and refraction seismic experiments conducted in the southern part ofLakeVostokshow very thin (200–300 m) stratified sedimentary cover overlying crystalline basement with velocity of 6.0–6.2 km/s. At present, deposition in southernLakeVostokis absent and similar conditions occurred likely at least last3 m.y. when ice sheet aboveLakeVostokchanged insignificantly. It can be also inferred that from the Late Miocene the rate of deposition inLakeVostokwas extremely low and so the most of sedimentary section is older being possibly of Oligocene to early to middle Miocene age when ice sheet oscillated and deposition was more vigorous. If so, the sampling of upper few meters of this condensed section is very informative in terms of history of Antarctic glaciation. Small thickness of sedimentary cover raises a question about existence of lake (rift) depression during preglacial and early glacial times.

Large explosive shot gathers along the SinoProbe deep seismic reflection profile and Moho depth beneath the Qiangtang terrane in central Tibet

Owing to extreme topography, rapid velocity and thickness variation of near-surface layer and strong seismic attenuation through the thickest crust of the Earth, it is not easy to acquire reflective Moho in central Tibet. However, with help from the SinoProbe deep seismic reflection experiment, large explosive of 1000 kg seismic sources have been tentatively detonated in the Qiangtang terrane and good quality data were acquired. Compared with multifold data, the single-fold record from single shot gather can also show clear Moho image in the Qiangtang terrane. Moho reflection appears at ∼ 24 s TWT (∼ 75.1 km) in the northernmost Lhasa terrane and at about 21 ∼ 20 s TWT (65.7 - 62.6 km) beneath the Qiangtang terrane. We speculate that Moho gets 9.4 km-12.5 km shallower from the Lhasa to the Qiangtang terrane rather than a 20 km offset. There is no obvious change of Moho depth across the Shung Hu suture.

Large heterogeneous structure beneath the Atotsugawa Fault, central Japan, revealed by seismic refraction and reflection experiments

A high-strain-rate zone termed the Niigata Kobe Tectonic Zone is located in central Japan and contains an active right-lateral fault called the Atotsugawa Fault. We present the results of an explosive-source seismic experiment that focused on identifying the formation mechanisms of the fault. The experiment used seven explosive sources and 1108 seismic stations, and was undertaken during October 2007 in the area of the Atotsugawa Fault. The seismic stations were set up on a survey line with a length of ~170km, oriented orthogonal to the strike of the fault. The linear array was used to identify the seismic structure of the fault zone and the deeper parts of the crust in the study area, yielding fine details of the seismic structure of the crust along the profile line, including lateral variations in P-wave velocity and the configuration of layers in the crust. A relatively low P-wave velocity reflective zone was detected beneath the fault at depths of 15–25km. This zone also has extremely low S-wave velocities, directly underlies three active faults, and contains three low-resistivity zones that underlie these faults. These reflective and low-resistivity zones are thought to represent fluid within the crust, and the data obtained during this study are consistent with a weak zone model for the formation of the active faults in the study area. This suggests that the reflective layer represents a weak zone in the lower crust that is responsible for the concentration of deformation within the upper crust. The presence of abundant fluids beneath the Atotsugawa Fault could also reduce the strength of the lower crust, again focusing deformation related to the regional stress regime in areas that overlie these fluid-containing zones. The results suggest that the reflective layer was an important factor in the formation of the Atotsugawa Fault.

Lithospheric velocity model across the Southern Central Iberian Zone (Variscan Iberian Massif): The ALCUDIA wide-angle seismic reflection transect

A P wave seismic velocity model has been obtained for the Central Iberian Zone, the largest continental fragment of the Iberian Variscan Belt. The spatially dense, high-resolution, wide-angle seismic reflection experiment, ALCUDIA-WA, was acquired in 2012 across central Iberia, aiming to constrain the lithospheric structure and resolve the physical properties of the crust and upper mantle. The seismic transect, ~310 km long, crossed the Central Iberian Zone from its suture with the Ossa-Morena Zone to the southern limit of the Central System mountain range. The energy generated by five shots was recorded by ~900 seismic stations. High-amplitude phases were identified in every shot gather for the upper crust (Pg and PiP) and Moho (PmP and Pn). In the upper crust, the P wave velocities increase beneath the Cenozoic Tajo Basin. The base of the upper crust varies from ~13 km to ~20 km between the southernmost Central Iberian Zone and the Tajo Basin. Lower crustal velocities are more homogeneous. From SW-NE, the traveltime of PmP arrivals varies from ~10.5 s to ~11.8 s, indicating lateral variations in the P wave velocity and the crustal thickness, reflecting an increase toward the north related with alpine tectonics and the isostatic response of the crust to the orogenic load. The results suggest that the high velocities of the upper crust near the Central System might correspond to igneous rocks and/or high-grade metamorphic rocks. The contrasting lithologies and the increase in the Moho depth to the north evidence differences in the Variscan evolution.

Local earthquake tomography of Scotland

Scotland is a relatively aseismic region for the use of local earthquake tomography, but 40 yr of earthquakes recorded by a good and growing network make it possible. A careful selection is made from the earthquakes located by the British Geological Survey (BGS) over the last four decades to provide a data set maximising arrival time accuracy and ray path coverage of Scotland. A large number of 1-D velocity models with different layer geometries are considered and differentiated by employing quarry blasts as ground-truth events. Then, SIMULPS14 is used to produce a robust 3-D tomographic P-wave velocity model for Scotland. In areas of high resolution the model shows good agreement with previously published interpretations of seismic refraction and reflection experiments. However, the model shows relatively little lateral variation in seismic velocity except at shallow depths, where sedimentary basins such as the Midland Valley are apparent. At greater depths, higher velocities in the northwest parts of the model suggest that the thickness of crust increases towards the south and east. This observation is also in agreement with previous studies. Quarry blasts used as ground truth events and relocated with the preferred 3-D model are shown to be markedly more accurate than when located with the existing BGS 1-D velocity model.

Seismic explosion sources on an ice cap – Technical considerations

Controlled source seismic investigation of crustal structure below ice covers is an emerging technique. We have recently conducted an explosive refraction/wide-angle reflection seismic experiment on the ice cap in east-central Greenland. The data-quality is high for all shot points and a full crustal model can be modelled. A crucial challenge for applying the technique is to control the sources. Here, we present data that describe the efficiency of explosive sources in the ice cover. Analysis of the data shows, that the ice cap traps a significant amount of energy, which is observed as a strong ice wave. The ice cap leads to low transmission of energy into the crust such that charges need be larger than in conventional onshore experiments to obtain reliable seismic signals. The strong reflection coefficient at the base of the ice generates strong multiples which may mask for secondary phases. This effect may be crucial for acquisition of reflection seismic profiles on ice caps. Our experience shows that it is essential to use optimum depth for the charges and to seal the boreholes carefully.

Crustal structure beneath the Rif Cordillera, North Morocco, from the RIFSIS wide‐angle reflection seismic experiment

AbstractThe different geodynamic models proposed since the late 1990s to account for the complex evolution of the Gibraltar Arc System lack definite constraints on the crustal structure of the Rif orogen. Here we present the first well‐resolved P‐wave velocity crustal models of the Rif Cordillera and its southern continuation toward the Atlas made using controlled‐source seismic data. Two 300+ km‐long wide‐angle reflection profiles crossed the Rif along NS and EW trends. The profiles recorded simultaneously five land explosions of 1Tn each using ∼850 high frequency seismometers. The crustal structure revealed from 2‐D forward modeling delineates a complex, laterally varying crustal structure below the Rif domains. The most surprising feature, seen on both profiles, is a ∼50 km deep crustal root localized beneath the External Rif. To the east, the crust thins rapidly by 20 km across the Nekkor fault, indicating that the fault is a crustal scale feature. On the NS profile the crust thins more gradually to 40 km thickness beneath Middle Atlas and 42 km beneath the Betics. These new seismic results are in overall agreement with regional trends of Bouguer gravity and are consistent with recent receiver function estimates of crustal thickness. The complex crustal structure of the Rif orogen in the Gibraltar Arc is a consequence of the Miocene collision between the Iberian and African plates. Both the abrupt change in crustal thickness at the Nekkor fault and the unexpectedly deep Rif crustal root can be attributed to interaction of the subducting Alboran slab with the North African passive margin at late Oligocene‐early Miocene times.

The Argentine Precordillera: A foreland thrust belt proximal to the subducted plate

The Precordillera thrust belt of western Argentina is anomalously close, both horizontally and vertically, to the coeval subduction zone of the Nazca plate. The thin-skinned part of the belt has an unusually deep decollement that is well defined by industry seismic reflection and recent broadband experiments. New area and line-length balanced cross sections show that the central Precordillera has accrued ∼90 ± 21 km of shortening since 13 Ma; much of that shortening occurred between 12 and 9 Ma. Fault-slip data generally show shortening approximately west-northwest–east-southeast, orthogonal to the traces of the thrust and folds in the Precordillera and oblique to the mean vector of local global positioning system (GPS) data. The GPS strain rate is –63 ± 9 × 10 –9 /yr, whereas strain rate in the thrust belt, averaged over 13 m.y., is –56 ± 4 × 10 –9 /yr. Although the decollement of the Precordillera cannot cut into Paleozoic Cuyania(?) terrane basement east of the crest of the high Andes, broadband receiver function data show that significant crustal thickening must occur beneath and even east of the thrust belt. We suggest that top-to-the-west shear and thickening of the lower crust due to flat subduction explains the distribution of crustal thickening.

Three-dimensional velocity structure of the outer fore arc of the Colombia-Ecuador subduction zone and implications for the 1958 megathrust earthquake rupture zone

In 2005, an onshore, offshore 3-D refraction and wide-angle reflection seismic experiment was conducted along the convergent margin at the border between Colombia and Ecuador, over the rupture zone of the 1958, Mw 7.6 subduction earthquake. A well-defined Vp velocity model of the plate boundary and upper and lower plates was constructed, down to 25 km depth, using first arrival traveltimes inversion. The model reveals a several kilometers thick, low-velocity zone in the upper plate, located immediately above the interplate contact. This low-velocity zone might be related to alteration and fracturing of the mafic and ultramafic rocks, which composed the upper plate in this area by fluids released by the lower plate with possible contributions from sediment underplating. Near the toe of the margin, the model shows a low-velocity gradient in the outer wedge, which is interpreted as highly faulted and fractured rocks. This low-velocity/low-gradient region appears to limit the oceanward extension of the rupture zones of the 1958 and 1979 earthquakes, possibly because coseismic deformation and uplift of the outer margin wedge dissipates most of the seismic energy.