The reflection seismic survey of project TIPTEQ—the inventory of the Chilean subduction zone at 38.2° S

We describe results of an active-source seismology experiment across the Chilean subduction zone at 38.2°S. The seismic sections clearly show the subducted Nazca plate with varying reflectivity. Below the coast the plate interface occurs at 25 km depth as the sharp lower boundary of a 2–5 km thick, highly reflective region, which we interpret as the subduction channel, that is, a zone of subducted material with a velocity gradient with respect to the upper and lower plate. Further downdip along the seismogenic coupling zone the reflectivity decreases in the area of the presumed 1960 Valdivia hypocentre. The plate interface itself can be traced further down to depths of 50–60 km below the Central Valley. We observe strong reflectivity at the plate interface as well as in the continental mantle wedge. The sections also show a segmented forearc crust in the overriding South American plate. Major features in the accretionary wedge, such as the Lanalhue fault zone, can be identified. At the eastern end of the profile a bright west-dipping reflector lies perpendicular to the plate interface and may be linked to the volcanic arc.     

Detailed 1×1° gravimetric Indian Ocean geoid and comparison with GEOS-3 radar altimeter geoid profiles

Summary. A new set of 1×1° mean free-air anomalies in the Indian Ocean is determined on the basis of previously published free-air anomaly maps (Talwani & Kahle) and the most recent Lamont surface ship gravity measurements. The data are then used to compute a (total) 1×1° gravimetric Indian Ocean geoid. The computation is carried out by combining the Goddard Space Flight Center (GSFC) GEM-6 geoid and a difference geoid that corresponds to the differences between the set of 1×1° surface gravity values and the GEM-6 gravity anomalies. The difference geoid is highest over the Madagascar Ridge (+ 20 m) and lowest over the Timor Trough (-30 m). The total geoid is compared with GEOS-3 radar altimeter derived geoid profiles and geophysical implications are discussed.     

Ccte d'Ivoire-Ghana margin: seismic imaging of passive rifted crust adjacent to a transform continental margin

During May 1990 and January-February 1991, an extensive geophysical data set was collected over the Côte d'Ivoire-Ghana continental margin, located along the equatorial coast of West Africa. The Ghana margin is a transform continental margin running subparallel to the Romanche Fracture Zone and its associated marginal ridge—the Côte d'Ivoire-Ghana Ridge. From this data set, an explosive refraction line running ∼ 150 km, ENE-WSW between 3°55'N, 3°21'W and 4°23'N, 2°4'W, has been modelled together with wide-angle airgun profiles, and seismic reflection and gravity data. This study is centred on the Côte d'Ivoire Basin located just to the north of the Côte d'Ivoire-Ghana Ridge, where bathymetric data suggest that a component of normal rifting occurred, rather than the transform motion observed along the majority of the equatorial West African margin.
Traveltime and amplitude modelling of the ocean-bottom seismometer data shows that the continental Moho beneath the margin rises in an oceanward direction, from ∼ 24 km below sea level to ∼ 17 km. In the centre of the line where the crust thins most rapidly, there exists a region of anomalously high velocity at the base of the crust, reaching some 8 km in thickness. This higher-velocity region is thought to represent an area of localized underplating related to rifting. Modelling of marine gravity data, collected coincident with the seismic line, has been used to test the best-fitting seismic model. This modelling has shown that the observed free-air anomaly is dominated by the effects of crustal thickness, and that a region of higher density is required at the base of the crust to fit the observed data. This higher-density region is consistent in size and location with the high velocities required to fit the seismic data.     

Analysis of gravity fields in the northwestern Himalayas and Kohistan region using deep seismic sounding data

A Bouguer gravity anomaly map of the NW Himalayas and parts of the Kohistan/Hindukush region has been prepared using all available gravity data. Analysis of the gravity field has been carried out along a profile extending from Gujranwala (located near the edge of the Indian shield) to the Haramosh massif in a NNE–SSW direction. The gravity profile is located close to the DSS profile shot under the USSR–India scientific collaborative programme. Velocity information available along different parts of the profile has been used to infer values of crustal and upper mantle density.
The observed gravity field (Bouguer) has been interpreted in terms of Moho depth and density contrast between the crust and the mantle. The Moho depth is interpreted as increasing from nearly 35 km near the edge of the Indian shield to 75 km (below sea-level) underneath the Haramosh massif. A similar model is applicable to a profile passing to the west of Nanga Parbat massif, from Gujranwala to Ghizar, through the Kohistan region. However, along this profile high-density lower-crustal rocks appear to have been emplaced in the upper part along the main mantle thrust. The nature of isostatic compensation prevailing underneath the Himalayas has been discussed, as has the theory of lithospheric flexure proposed by Karner & Watts and Lyon-Caen & Molnar. It is felt that although these ideas explain the broad features of the Moho configuration as observed in the NW Himalayas, there are significant departures. The role of tectonic forces in shaping the Moho and causing changes in the density of the crust cannot be denied.     

Characterization of the depocenters and the basement structure,below the central Chile Andean Forearc: A 3D geophysical modelling in Santiago Basin area

Since the last century, several geological and geophysical studies have been developed in the Santiago Basin to understand its morphology and tectonic evolution. However, some uncertainties regarding sedimentary fill properties and possible density anomalies below the sediments/basement boundary remain. Considering that this is an area densely populated with more than 6 million inhabitants in a highly active seismotectonic environment, the physical properties of the Santiago Basin are important to study the geological and structural evolution of the Andean forearc and to characterize its seismic response and related seismic hazard. Two and three‐dimensional gravimetric models were developed, based on a database of 797 compiled and 883 newly acquired gravity stations. To produce a well‐constrained basement elevation model, a review of 499 wells and 30 transient electromagnetic soundings were used, which contribute with basement depth or minimum sedimentary thickness information. For the 2‐D modelling, a total of 49 gravimetric profiles were processed considering a homogeneous density contrast and independent regional trends. A strong positive gravity anomaly was observed in the centre of the basin, which complicated the modelling process but was carefully addressed with the available constrains. The resulting basement elevation models show complex basement geometry with, at least, eight recognizable depocenters with maximum sedimentary infill of ~ 500 m. The 3‐D density models show alignments in the basement that correlates well with important intrusive units of the Cenozoic and Mesozoic. Along with interpreted fault zones westwards and eastwards of the basin, the observations suggest a structural control of Santiago basin geometry, where recent deformation associated with the Andean contractional deformation front and old structures developed during the Cenozoic extension are superimposed to the variability of river erosion/deposition processes.     

Gravity constraints on the dynamics of the crust–mantle system during Calabrian subduction

The thermomechanic evolution of the lithosphere–upper mantle system during Calabrian subduction is analysed using a 2-D finite element approach, in which the lithosphere is compositionally stratified into crust and mantle. Gravity and topography predictions are cross-checked with observed gravity and topography patterns of the Calabrian region. Modelling results indicate that the gravity pattern in the arc-trench region is shaped by the sinking of light material, belonging to both the overriding and subduction plates. The sinking of light crustal material, up to depths of the order of 100–150 km is the ultimate responsible for the peculiar gravity signature of subduction, characterized by a minimum of gravity anomaly located at the trench, bounded by two highs located on the overriding and subducting plates, with a variation in magnitude of the order of 200 mGal along a wavelength of 200 km, in agreement with the isostatically compensated component of gravity anomaly observed along a transect crossing the Calabrian Arc, from the Tyrrhenian to the Ionian Seas. The striking agreement between the geodetic retrieved profiles and the modelled ones in the trench region confirms the crucial role of compositional stratification of the lithosphere in the subduction process and the correctness of the kinematic hypotheses considered in our modelling, that the present-day configuration of crust–mantle system below the Calabrian arc results from trench's retreat at a rate of about 3 cm yr⁻¹, followed by gravitational sinking of the subducted slab in the last 5 Myr.     

Results of magnetotelluric and gravimetric measurements in western Nicaragua, Central America

Magnetotelluric and gravity data have been collected within a ca. 170 km long traverse running from the Pacific coast of Nicaragua in the west to the Nicaraguan Highland in the east. This part of Nicaragua is characterized by sedimentary rocks of the Pacific Coastal Plain, separated from the Tertiary volcanic rocks of the Highland by the NW-SE-trending Nicaraguan Depression. 2-D interpretation of the magnetotelluric (MT) data, collected at 13 stations, indicates four regions of high electrical conductivity in addition to the conductive coastal region. Two of these are associated with conducting sediments and pyroclastics in the upper part of the crust. Two other conductive structures have been defined at depth around 20 km and the one best defined is located below the depression. From the distribution of seismic events, volcanic activity in the depression and the similarity in geophysical characteristics with areas such as the Rio Grande Rift, this conductor is interpreted as a melt layer or a complex of magma chambers. Models of the upper lithosphere, constrained by the MT model, vertical electrical sounding (VES) data, seismic data and densities, have been tested using gravity data. A model that passes this test shows a gradual thickening of the crust eastwards from the Pacific coast. An anomaly centred over the depression is interpreted to have its origin in a thinning of the crust. In this model the melt layer is situated on top of the bulge of the lower lithosphere. A change in the composition of the crust, from the Pacific Coastal Plain to the Highland, is indicated from the change in character of the MT response and from the density distribution in the gravity model. This may support the hypothesis that the Pacific region is an accreted terrane. MT and gravity data indicate a depth to a resistive and high-density basement in the depression of ca. 2 km. On the basis of this, the vertical setting in the depression is estimated to be of the order of 2.5 km.     

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.     

A seismic study along the East Greenland margin from 72°N to 77°N

Around 4370 km of new seismic reflection data, collected along the East Greenland margin between 71°30'N and 77°N in 2003, provide a first detailed view of the sediment distribution and tectonic features along the East Greenland margin. After processing and converting the data to depth, we correlated ODP-Site 913 stratigraphy into the new seismic network. Unit GB-2 shows the greatest glacial sediment deposits beneath the East Greenland continental shelf. This unit is characterized by the beginning of prograding sequences and has, according to our stratigraphic correlation, a Middle Miocene age. It might have been caused by rapid changes in sea level and/or glacial erosion by an early ice sheet or glaciers along the coast. A basement high, presumably a 360 km long basement structure at 77°N–74°54'N, prevents continuous sediment transport from the shelf into the deep sea area in times before 15 Myr. The origin of this prominent structure remains speculative since no rock sample from this structure is available. Seaward dipping reflectors at the eastern flank of this structure strongly support that it is a volcanic construction and is most likely emplaced on continental or transitional crust. The compilation of sediment thickness provide an insight into the regional sediment distribution in the Greenland Basin. An average sediment thickness of 1 km is observed. The north bordering Boreas Basin has a sediment thickness of 1.8 km close to the Greenland fracture zone (GFZ).     

南极布兰斯菲尔德盆地重力异常研究

１９９１年，“海洋四号”科学考察船在南极布兰斯菲尔德海峡采集了大量的地球物理资料。本文在利用地震调查资料的基础上，根据现有的重力测量数据，研究引起重力异常的各种地质因素，对实测剖面重力异常进行正演拟合，建立了布兰斯菲尔德盆地的密度模型，并推导出该地区莫霍面深度－布格异常关系表达式，绘制出盆地的莫霍面深度图。     

The Mid-Atlantic Ridge: structure at 45° N

Summary. A structural model of the Mid-Atlantic Ridge at 45° N is proposed on the basis of travel-time data, amplitudes and synthetic seismograms. The crustal structure seems to be similar to that in the FAMOUS area (Fowler). At the ridge axis there is an absorptive zone in the upper mantle, the depth below the seabed to the top of this zone being about 6 km. Away from the ridge axis there is a positive velocity gradient of about 0.04 to 0.05 km/(skm) in the top 5 to 8 km of the upper mantle. Shear waves propagate across the ridge axis, suggesting that there is no sizeable crustal magma chamber. The shear-wave velocity of the uppermost mantle is 4.35 km/s.     

Anomalous high-frequency wave propagation from the Tonga–Kermadec seismic zone to New Zealand

Summary. Bulletins of the International Seismological Centre (ISC) show very large residuals, up to 15 s early, for arrivals from events in the Tonga–Kermadec subduction zone to the New Zealand network of seismometers. The very early arrivals are confined to events south of about 22°S, and shallower than about 350 km. The waveforms show two distinct phases: an early, emergent, first phase with energy in the high-frequency band 2–10 Hz, and a distinct second phase, containing lower frequency energy, arriving at about the time predicted by JB tables.
The residuals are attributed to propagation through the cold, subducted lithosphere, which has a seismic velocity 5 per cent faster, on average, than normal. Ray tracing shows that the ray paths lie very close to the slab for events south of 22°S, but pass well beneath the slab for events further north, corresponding to the change in residual pattern. This characteristic of the ray paths is due to the curved shape of the seismic zone, and in particular to the bend in the zone where the Louisville ridge intersects the trench at 25°S.
The residuals can only be explained if the high velocity anomaly extends to a depth of 450 km in the region of the gap in deep seismicity from 32 to 36°S. The very high-frequency character of the first phase requires the path from the bottom of the slab to the stations to be of high Q , and to transmit 2–10 Hz energy with little attenuation.
The absence of low-frequency energy in the first phase is due to the narrowness of the high-velocity slab, which transmits only short-wavelength waves. The second phase, which contains low frequencies, is identified as a P -wave travelling beneath the subducted slab in normal mantle. There is no need to invoke any special structures, such as low-velocity waveguides or reflectors, to explain any of the observations. The S -wave arrivals show similar effects.     

Modelling of sedimentary basins from gravity anomalies with variable density contrast

Summary. The decrease of density contrast with depth in sedimentary basins could be approximated by a quadratic function. The sedimentary basin is viewed as a number of outcropping vertical prisms put in juxtaposition at each observation point. Equations are derived for the gravity anomaly of a vertical prism for a quadratic density function. The initial model is obtained by calculating the thicknesses of the prisms at each station by an infinite slab formula. The differences between the observed anomalies and the anomalies of the initial model at each station are used to modify the initial model. A computer program is developed for initial estimates and refinement of the model. The method is applied to interpret the gravity anomaly of San Jacinto Graben, California, where sediments filling the graben have a density contrast that can be approximated to a quadratic function. Finally the convergence of the method is discussed.     

Why is there no Coastal Cordillera at the Arica Bend (Western Central Andes)?

The architecture of the Western Andes is remarkably constant between southern Peru and northern Chile. An exception, however, is present near Arica at 18°S, where the Andes change their strike direction by ca. 50° and the Coastal Cordillera is absent over an along‐strike width of 50 km. Although this feature has been mentioned in several previous studies, no effort has been made yet to describe and explain this peculiar morphology of the Western Central Andean forearc. Here, we propose a large‐scale model to explain the Myr‐long low uplift rate of the Arica Bend concerning seismic coupling and continental wedge‐top basin evolution. New geomorphic and sedimentologic data are integrated with seismicity and structural data from the literature to interpret the post‐Oligocene pattern of uplift, erosion and sediment transport to the trench. Results show that the Arica Bend has been marked by exceptionally low coastal uplift rates over post‐Oligocene timescales. In addition, this uplift anomaly at the Arica Bend correlates with relatively high sediment discharge to the corresponding trench segment since late Oligocene time. We interpret that before 25 Ma, the forming seaward concavity of the subduction zone induced trench‐parallel extension at the curvature apex of the overriding forearc. The subsequent low uplift rate would have then triggered a feedback mechanism, where the interplay between relatively low interplate friction, low coastal uplift and relatively high sediment discharge favoured Myr‐long relative subsidence at the Arica Bend, in contrast to Myr‐long uplift of the Coastal Cordillera north and south of it.     

Gravitational deformation and inherited structural control on slope morphology in the subduction zone of north‐central Chile (ca. 29–33°S)

Subduction zones provide direct insight into plate boundary deformation and by studying these areas we better understand tectonic processes and variability over time. We studied the structure of the offshore subduction zone system of the Pampean flat‐slab segment (ca. 29–33°S) of the Chilean margin using seismic and bathymetric constraints. Here, we related and analysed the structural styles of the offshore and onshore western fore‐arc. Overlying the acoustic top of the continental basement, two syn‐extensional seismic sequences were recognised and correlated with onshore geological units and the Valparaíso Forearc Basin seismic sequences: (SII) Pliocene‐Pleistocene and (SI) Miocene‐Pliocene (Late Cretaceous (?) to Miocene‐Pliocene) syn‐extensional sequences. These sequences are separated by an unconformity (i.e. Valparaíso Unconformity). Seismic reflection data reveal that the eastward dipping extensional system (EI) recognised at the upper slope can be extended to the middle slope and controlled the accumulation of the older seismic package (SI). The westward dipping extensional system (EII) is essentially restricted to the middle slope. Here, EII cuts the eastward dipping extensional system (EI), preferentially parallel to the inclination of the older sequences (SI), and controlled a series of middle slope basins which are filled by the Pliocene‐Pleistocene seismic sequence (SII). At the upper slope and in the western Coastal Cordillera, the SII sequence is controlled by eastward dipping faults (EII) which are the local reactivation of older extensional faults (EI). The tectonic boundary between the middle (eastern outermost forearc block) and upper continental slope (western coastal block) is a prominent system of trenchward dipping normal fault scarps (ca. 1 km offset) that resemble a major trenchward dipping extensional fault system. This prominent structural feature can be readily detected along the Chilean erosive margin as well as the two extensional sets (EI and EII). Evidence of slumping, thrusting, reactivated faults and mass transport deposits, were recognised in the slope domain and locally restricted to some eastern dipping faults. These features could be related to gravitational effects or slope deformation due to coseismic deformation. The regional inclination of the pre‐Pliocene sequences favoured the gravitational collapse of the outermost forearc block. We propose that the structural configuration of the study area is dominantly controlled by tectonic erosion as well as the uplift of the Coastal Cordillera, which is partially controlled by pre‐Pliocene architecture.     

Regional three-dimensional gravity modelling of the NE Atlantic margin

A series of three‐dimensional models has been constructed for the structure of the crust and upper mantle over a large region spanning the NE Atlantic passive margin. These incorporate isostatic and flexural principles, together with gravity modelling and integration with seismic interpretations. An initial isostatic model was based on known bathymetric/topographic variations, an estimate of the thickness and density of the sedimentary cover, and upper mantle densities based on thermal modelling. The thickness of the crystalline crust in this model was adjusted to equalise the load at a compensation depth lying below the zone of lateral mantle density variations. Flexural backstripping was used to derive alternative models which tested the effect of varying the strength of the lithosphere during sediment loading. The models were analysed by comparing calculated and observed gravity fields and by calibrating the predicted geometries against independent (primarily seismic) evidence. Further models were generated in which the thickness of the sedimentary layer and the crystalline crust were modified in order to improve the fit to observed gravity anomalies. The potential effects of igneous underplating and variable upper mantle depletion were explored by a series of sensitivity trials. The results provide a new regional lithospheric framework for the margin and a means of setting more detailed, local investigations in their regional context. The flexural modelling suggests lateral variations in the strength of the lithosphere, with much of the margin being relatively weak but areas such as the Porcupine Basin and parts of the Rockall Basin having greater strength. Observed differences between the model Moho and seismic Moho along the continental margin can be interpreted in terms of underplating. A Moho discrepancy to the northwest of Scotland is ascribed to uplift caused by a region of upper mantle with anomalously low density, which may be associated with depletion or with a temperature anomaly.     

Gravimetric determination of an intrusive complex under the Island of Faial (Azores): some methodological improvements

We present some improvements of a gravity inversion method to determine the geometry of the anomalous bodies for priori density contrasts. The 3-D method is based on an exploratory process applied, not for the global model, but for the steps of a growth approach. The (positive and/or negative) anomalous structure is described by successive aggregation of cells, while its corresponding gravity field remains nearly proportional to the observed one. Moreover, a simple (e.g. linear) regional trend can be simultaneously adjusted. The corresponding program is applied to new gravity data on the volcanic island of Faial (Azores archipelago). The inversion approach shows a subsurface anomalous structure for the island, the main feature being an elongated high-density body. The body is interpreted as a compact sheeted dyke swarm, emplaced along Faial-Pico Fracture Zone, a leaky transform structure that forms the current boundary between Eurasian and African plates in the Azores area. The new results in this paper are (1) a Bouguer gravity anomaly map, (2) several improvements in the inversion process (robust process, optimal balance fitness/model magnitude), (3) a new gravimetric method for estimating the mean terrain density, (4) a 3-D model for subsurface mass anomalies in Faial and (5) some interpretative conclusions about a main intrusive complex detected under the island as a wall-like structure extending from a depth of 0.5 to 6 km b.s.l., with a N100°E trend and corresponding to an early fissural volcanic episode controlled by the regional tectonics.     

Density inhomogeneities of the upper mantle of the Central Balkans

A 3-D density model was created for the Central Balkans area down to a depth of 670  km on the basis of seismic (both artificial sources and earthquakes) and gravity data. This model is based on density columns constructed for the main geological units of the study region. The densities for these columns were obtained using the density variation method. This method makes it possible to extrapolate the density distribution from the well-studied uppermost layers to the deeper levels of the Earth. The constructed 3-D density model was interpreted in relation to the available data on the heat flow and the seismicity of the region. The subdivision of the region by the Maritza fault into two parts—the southern part including the Rhodope massif and the northern part including the structures of Alpine activation of Srednogorie and the Balkanides—was confirmed. The upraised position of the 400  km boundary within the upper mantle, which was established from the density modelling, is assumed to be a sign of development of recent geodynamical processes over the Srednogorie block. From the viewpoint of seismicity prediction, a finding of mantle inhomogeneities orthogonal to the Maritza suture is of great importance.     

Lithoprobe seismic reflection profiling across Vancouver Island: results from reprocessing

Summary. Multichannel seismic reflection sections recorded across Vancouver Island have revealed two extensive zones of deep seismic reflections that dip gently to the northeast, and a number of moderate northeasterly dipping reflections that can be traced to the surface where major faults are exposed. Based on an integrated interpretation of these data with information from gravity, heat flow, seismicity, seismic refraction, magnetotelluric and geological studies it is concluded that the lower zone of gently dipping reflections is due to underplated oceanic sediments and igneous rocks associated with the current subduction of the Juan de Fuca plate, and that the upper zone represents a similar sequence of accreted rocks associated with an earlier episode of subduction. The high density/high velocity material between the two reflection zones is either an underplated slab of oceanic lithosphere or an imbricated package of mafic rocks. Reprocessing of data from two of the seismic lines has produced a remarkable image of the terrane bounding Leech River fault, with its dip undulating from >60° near the surface to 20° at 3 km depth and ∼38° at 6 km depth.     

南极半岛和东南极间南威德尔海陆架的地壳结构

重力资料结合磁性体埋深计算和地震折射资料的分析 ,表明在南极半岛和东南极克拉通之间的菲尔希内尔和罗纳冰架之下存在一减薄的地壳。推测该区地壳的减薄主要由南极半岛和东南极克拉通之间的近东西向拉张作用所造成 ,该作用很可能发生在陆架盆地发展的早期。