Lithosphere structure of European sedimentary basins from regional three-dimensional gravity modelling

A three-dimensional (3D) density model, approximated by two regional layers—the sedimentary cover and the crystalline crust (offshore, a sea-water layer was added), has been constructed in 1° averaging for the whole European continent. The crustal model is based on simplified velocity model represented by structure maps for main seismic horizons—the “seismic” basement and the Moho boundary. Laterally varying average density is assumed inside the model layers. Residual gravity anomalies, obtained by subtraction of the crustal gravity effect from the observed field, characterize the density heterogeneities in the upper mantle. Mantle anomalies are shown to correlate with the upper mantle velocity inhomogeneities revealed from seismic tomography data and geothermal data. Considering the type of mantle anomaly, specific features of the evolution and type of isostatic compensation, the sedimentary basins in Europe may be related into some groups: deep sedimentary basins located in the East European Platform and its northern and eastern margins (Peri-Caspian, Dnieper–Donets, Barents Sea Basins, Fore–Ural Trough) with no significant mantle anomalies; basins located on the activated thin crust of Variscan Western Europe and Mediterranean area with negative mantle anomalies of −150 to −200×10⁻⁵ ms⁻² amplitude and the basins associated with suture zones at the western and southern margins of the East European Platform (Polish Trough, South Caspian Basin) characterized by positive mantle anomalies of 50–150×10⁻⁵ ms⁻² magnitude. An analysis of the main features of the lithosphere structure of the basins in Europe and type of the compensation has been carried out.     

Crustal and mantle strengths in continental lithosphere: is the jelly sandwich model obsolete?

The relative importance of the contribution of the lower crust and of the lithospheric mantle to the total strength of the continental lithosphere is assessed systematically for realistic ranges of layer thickness, composition, and temperature. Results are presented as relative strength maps, giving the ratio of the lower crust to upper mantle contribution in terms of crustal thickness and surface heat flow. The lithosphere shows a “jelly sandwich” rheological layering for low surface heat flow, thin to average crustal thickness, and felsic or wet mafic lower crustal compositions. On the other hand, most of the total strength resides in the seismogenic crust in regions of high surface heat flow, crust of any thickness, and dry mafic lower crustal composition.     

Modelling the extension of heterogeneous hot lithosphere

The consequences of weak heterogeneities in the extension of soft and hot lithosphere without significant previous crustal thickening has been analysed in a series of centrifuge models. The experiments examined the effects of i) the location of heterogeneities in the ductile crust and/or in the lithospheric mantle, and ii) their orientation, perpendicular or oblique to the direction of bulk extension. The observed deformation patterns are all relevant to the so-called “wide rifting” mode of extension. Weak zones located in the ductile crust exert a more pronounced influence on localisation of deformation in the brittle layer than those located in the lithospheric mantle: the former localise faulting in the brittle crust whereas the latter tend to distribute faulting over a wider area. This latter behaviour depends in turn upon the decoupling provided by the ductile crust. Localised thinning in the brittle crust is accompanied by ductile doming of both crust and mantle. Domains of maximum thinning in the brittle crust and ductile crust and mantle are in opposition. Lateral differences in brittle crust thinning are accommodated by lateral flow in the ductile crust and mantle. This contrasts with “cold and strong” lithospheres whose high strength sub-Moho mantle triggers a necking instability at the lithosphere-scale. This also differs from the extension of thickened hot and soft lithospheres whose ductile crust is thick enough to give birth to metamorphic core complexes. Thus, for the given lithospheric rheology, the models have relevance to backarc type extensional systems, such as the Aegean and the Tyrrhenian domains.     

Effects of buoyancy and mechanical layering on collisional deformation of continental lithosphere: Results from physical modeling

We use two suites of lithospheric-scale physical experiments to investigate the manner in which deformation of the continental lithosphere is affected by both (1) variations of lithospheric density (quantified by the net buoyant mass per area in the lithospheric mantle layer, M_B), and (2) the degree of coupling between the crust and lithospheric mantle (characterized by a modified Ampferer ratio, Am). The dynamics of the experiments can be characterized with a Rayleigh–Taylor type ratio, _CLM. Models with a positively buoyant lithospheric mantle layer (M_B > 0 and _CLM > 0) result in distributed root formation and a wide deformation belt. In contrast, models with a negatively buoyant lithospheric mantle layer strongly coupled to the crust (M_B < 0, 0 > _CLM > ≈ − 0.2, and Am > ≈ 10^{− 3}) exhibit localized roots and narrow deformation belts. Syncollisional delamination of the model lithospheric mantle layer and a wide deformation belt is exhibited in models with negatively buoyant lithospheric mantle layers weakly coupled to the crust (M_B < 0, _CLM < 0, and Am < ≈ 10^{− 3}). Syncollisional delamination of the continental lithosphere may initiate due to buoyancy contrasts within the continental plate, instead of resulting from wedging by the opposing plate. Rayleigh–Taylor instabilities dominate the style of deformation in models with a negatively buoyant lithospheric mantle layer strongly coupled to the crust and a slow convergence rate (M_B < 0 and _CLM > ≈ − 0.2). The degree of coupling (Am) between the model crust and lithospheric mantle plays a lesser role in both the style of lower-lithospheric deformation and the width of the crustal deformed zone with increasing density of the lithospheric mantle layer.     

Lithospheric transition from the Variscan Iberian Massif to the Jurassic oceanic crust of the Central Atlantic

A 1000-km-long lithospheric transect running from the Variscan Iberian Massif (VIM) to the oceanic domain of the Northwest African margin is investigated. The main goal of the study is to image the lateral changes in crustal and lithospheric structure from a complete section of an old and stable orogenic belt—the Variscan Iberian Massif—to the adjacent Jurassic passive margin of SW Iberia, and across the transpressive and seismically active Africa–Eurasia plate boundary. The modelling approach incorporates available seismic data and integrates elevation, gravity, geoid and heat flow data under the assumptions of thermal steady state and local isostasy. The results show that the Variscan Iberian crust has a roughly constant thickness of 30 km, in opposition to previous works that propose a prominent thickening beneath the South Portuguese Zone (SPZ). The three layers forming the Variscan crust show noticeable thickness variations along the profile. The upper crust thins from central Iberia (about 20 km thick) to the Ossa Morena Zone (OMZ) and the NE region of the South Portuguese Zone where locally the thickness of the upper crust is <8 km. Conversely, there is a clear thickening of the middle crust (up to 17 km thick) under the Ossa Morena Zone, whereas the thickness of the lower crust remains quite constant (6 km). Under the margin, the thinning of the continental crust is quite gentle and occurs over distances of 200 km, resembling the crustal attitude observed further north along the West Iberian margins. In the oceanic domain, there is a 160-km-wide Ocean Transition Zone located between the thinned continental crust of the continental shelf and slope and the true oceanic crust of the Seine Abyssal Plain. The total lithospheric thickness varies from about 120 km at the ends of the model profile to less than 100 km below the Ossa Morena and the South Portuguese zones. An outstanding result is the mass deficit at deep lithospheric mantle levels required to fit the observed geoid, gravity and elevation over the Ossa Morena and South Portuguese zones. Such mass deficit can be interpreted either as a lithospheric thinning of 20–25 km or as an anomalous density reduction of 25 kg m⁻³ affecting the lower lithospheric levels. Whereas the first hypothesis is consistent with a possible thermal anomaly related to recent geodynamics affecting the nearby Betic–Rif arc, the second is consistent with mantle depletion related to ancient magmatic episodes that occurred during the Hercynian orogeny.     

Mantle structure and rifting processes in the Baikal–Mongolia region: geophysical data and evidence from xenoliths in volcanic rocks

The origin of the Baikal rift zone (BRZ) has been debated between the advocates of passive and active rifting since the 1970s. A re-assessment of the relevant geological and geophysical data from Russian and international literature questions the concept of broad asthenospheric upwelling beneath the rift zone that has been the cornerstone of many “active rifting” models. Results of a large number of early and recent studies favour the role of far-field forces in the opening and development of the BRZ. This study emphasises the data obtained through studies of peridotite and pyroxenite xenoliths brought to the surface by alkali basaltic magmas in southern Siberia and central Mongolia. These xenoliths are direct samples of the upper mantle in the vicinity of the BRZ. Of particular importance are suites of garnet-bearing xenoliths that have been used to construct P–T- composition lithospheric cross-sections in the region for the depth range of 35–80 km.Xenolith studies have shown fundamental differences in the composition and thermal regime between the lithospheric mantle beneath the ancient Siberian platform (sampled by kimberlites) and beneath younger mobile belts south of the platform. The uppermost mantle in southern Siberia and central Mongolia is much hotter at similar levels than the mantle in the Siberian craton and also has significantly higher contents of ‘basaltic’ major elements (Ca, Al, Na) and iron, higher Fe/Si and Fe/Mg. The combination of the moderately high geothermal gradient and the fertile compositions in the off-cratonic mantle appears to be a determining factor controlling differences in sub-Moho seismic velocities relative to the Siberian craton. Chemical and isotopic compositions of the off-cratonic xenoliths indicate small-scale and regional mantle heterogeneities attributed to various partial melting and enrichment events, consistent with long-term evolution in the lithospheric mantle. Age estimates of mantle events based on Os–Sr–Nd isotopic data can be correlated with major regional stages of crustal formation and may indicate long-term crust–mantle coupling. The ratios of ^143/144Nd in many LREE-depleted xenoliths are higher than those in MORB or OIB source regions and are not consistent with a recent origin from asthenospheric mantle.Mantle xenoliths nearest to the rift basins (30–50 km south of southern Lake Baikal) show no unequivocal evidence for strong heating, unusual stress and deformation, solid state flow, magmatic activity or partial melting that could be indicative of an asthenospheric intrusion right below the Moho. Comparisons between xenoliths from older and younger volcanic rocks east of Lake Baikal, together with observations on phase transformations and mineral zoning in individual xenoliths, have indicated recent heating in portions of the lithospheric mantle that may be related to localised magmatic activity or small-scale ascent of deep mantle material. Overall, the petrographic, P–T, chemical and isotopic constraints from mantle xenoliths appear to be consistent with recent geophysical studies, which found no evidence for a large-scale asthenospheric upwarp beneath the rift, and lend support to passive rifting mechanism for the BRZ.     

The 1.80–1.74-Ga gabbro–anorthosite–rapakivi Korosten Pluton in the Ukrainian Shield: a 3-D geophysical reconstruction of deep structure

The deep structure of the gabbro–anorthosite–rapakivi granite (“AMCG-type”) Korosten Pluton (KP) in the northwestern Ukrainian Shield was studied by 3-D modelling of the gravity and magnetic fields together with previous seismic data. The KP occupies an area of ca. 12,500 km² and comprises several layered gabbro-anorthositic intrusions enveloped by large volumes of rapakivi-type granitoids. Between 1.80 and 1.74 Ga, the emplacement of mafic and associated granitoid melts took place in several pulses. The 3-D geophysical reconstruction included: (a) modelling of the density distribution in the crust using the observed Bouguer anomaly field constrained by seismic data on Moho depth, and (b) modelling of the magnetic anomaly field in order to outline rock domains of various magnetisation, size and shape in the upper and lower crust. The density modelling was referred to three depth levels of 0 to 5, 5 to 18, and 18 km to Moho, respectively. The 3-D reconstruction demonstrates close links between the subsurface geology of the KP and the structure of the lower crust. The existence of a non-magnetic body with anomalously high seismic velocity and density is documented. Most plausibly, it represents a gabbroic stock (a parent magma chamber) with a vertical extent of ca. 20 km, penetrating the entire lower crust. This stock has a half-cylindrical shape and a diameter of ca. 90 km. It appears to be connected with a crust–mantle transitional lens previously discovered by EUROBRIDGE seismic profiling. The position of the stock relative to the subsurface outlines of the KP is somewhat asymmetric. This may be due to a connection between the magmatism and sets of opposite-dipping faults initially developed during late Palaeoproterozoic collisional deformation in the Sarmatian crustal segment. Continuing movements and disturbances of the upper mantle and the lower crust during post-collisional tectonic events between 1.80 and 1.74 Ga may account for the long-lived, recurrent AMCG magmatism.     

Double-sided onshore–offshore seismic imaging of a plate boundary: “super-gathers” across South Island, New Zealand

Onshore–offshore seismic refraction profiling allows for the determination of crustal and mantle structures in the transition between continental and oceanic environments. Islands and narrow landmasses have the unique geometry of allowing for double-sided onshore–offshore experiments that favor the construction of composite “super-gathers” using the acquisition of onshore–offshore and ocean-bottom seismometer receiver gathers, land explosion shot gathers, and near-vertical incidence multichannel seismic (MCS) profiling. A number of sites at plate boundaries are amenable to the application of double-sided onshore–offshore imaging, including the Indo-Australian/Pacific transform boundary on South Island, New Zealand. By comparing the ratio of island width to mantle refraction (Pn) “maximum” crossover distance, using nondimensional distances, we provide an indicator of raypath “coverage” for crustal illumination. Islands or narrow land masses whose widths are less than twice their maximum crossover distance are candidates for double-sided onshore–offshore experiments. The SIGHT (South Island GeopHysical invesTigation) experiment in New Zealand is located where the width of South Island is sufficiently narrow with respect to its crustal thickness that a double-sided onshore–offshore experiment allows for complete crustal imaging of the associated plate boundary.     

The isostatic state of the southern Urals crust

The Uralide orogen, in Central Russia, is the focus of intense geoscientific investigations during recent years. The international research is motivated by some unusual lithospheric features compared with other collisional belts including the preservation of (a) a collisional architecture with an orogenic root and a crustal thickness of 55–58 km, and (b) large volumes of very low-grade and non-metamorphic oceanic crust and island arc rocks in the upper crust of a low–relief mountain belt. The latter cause anomalous gravity highs along the thickened crust and the isostatic equilibrium inside the Uralides lithosphere as well as the overthrust high-metamorphic rocks. The integrated URSEIS '95 seismic experiment provides fundamentally new data revealing the lithospheric architecture of an intact Paleozoic collisional orogen that allows the construction of density models. In the Urals' lithosphere different velocity structures resolved by wide-angle seismic experiments along both the URSEIS '95- and the Troitsk profile. They can be used to constrain lithospheric density models: a first model consists of a deep subducted continental lower crust which has been highly eclogitized at depths of 60–90 km to a density of 3550 kg/m³. The second model shows a slightly eclogitized lower crust underlying the Uralide orogen with a crustal thickness of 60 km. The eclogitized lower crust causes a too-small impedance contrast to the lithospheric mantle resulting in a lack of reflectors in the area of the largest crustal thickness. Both models fit the measured gravity field. Analyzing the isostatic state of the southern Urals' lithosphere, both density models are in isostatic equilibrium.     

The structure, isostasy and gravity field of the Levant continental margin and the southeast Mediterranean area

A 3-D layered structure of the Levant and the southeastern Mediterranean lithospheric plates was constructed using interpretations of seismic measurements and borehole data. Structural maps of three principal interfaces, elevation, top basement and the Moho, were constructed for the area studied. This area includes the African, Sinai and Arabian plates, the Herodotus and the Levant marine basins and the Nile sedimentary cone. In addition, an isopach map of the Pliocene sediments, as well as the contemporaneous amount of denuded rock units, was prepared to enable setting up the structural map of the base Pliocene sediment. Variable density distributions are suggested for the sedimentary succession in accord with its composition and compaction. The spatial density distribution in the crystalline crust was calculated by weighting the thicknesses of the lower mafic and the upper felsic crustal layers, with densities of 2.9 g/cm³ and 2.77 g/cm³, respectively. Results of the local (Airy) isostatic modeling with compensation on the Moho interface show significant deviations from the local isostasy and require variable density distribution in the upper mantle. Moving the compensation level to the base of the lithosphere ( 100 km depth) and adopting density variations in the mantle lithosphere yielded isostatic compensation (± 200 m) over most of the area studied. The spatial pattern obtained of a density distribution with a range of ± 0.05 g/cm³ is supported by a regional heat flux. Simulations of the flexure (Vening Meinesz) isostasy related to the Pliocene to Recent sedimentary loading and unloading revealed concentric oscillatory negative and positive anomalies mostly related to the Nile sedimentary cone. Such anomalies may explain the rapid subsidence in the Levant Basin and the arching in central Israel, northern Sinai and Egypt during Pliocene–Recent times. Comparison between the observed (Bouguer) gravity and the calculated gravity for the constructed 3-D lithospheric structure, which has variable density distributions, provided a good match and an independent constraint for the large-scale structure suggested and confirmed an oceanic nature for the Levant Basin lithosphere.     

Teleseismic p-wave traveltime residuals and deep structure of the Aegean region

Teleseismic P-wave traveltime residuals have been measured at the Greek seismic stations with respect to the Herrin 68 tables. In spite of the large scatter, some insight into crustal and upper-mantle structure of the Aegean region can be gained. The average absolute residuals (observed minus Herrin traveltimes) are of the order of + 2 s. The most plausible interpretation is an efficient low-velocity zone in the upper mantle. Simple estimates of densities and subsequently gravity with the aid of Birch's law suggest that the Aegean region is underlain by hot expanded upper mantle, perhaps involving partial melting. The relative P residuals (observed minus Herrin traveltime differences between a station and Athens) are generally positive and can be interpreted with lateral variations of the LVZ or of the crust. The latter interpretation is supported in some cases by seismic refraction data. The azimuthal variation of the relative residuals at stations on the non-volcanic arc bears a distinct relation with the arc orientation. At Archangelos (Rhodes) where we “see” through the Benioff zone, the residuals from N to W are between -1 and -2 s and indicate a high-velocity slab sinking below the Aegean sea. At Vamos (Crete), Valsamata (Kephallenia), and Joanina (Pindus Range) the largest (smallest) residuals are along directions parallel (perpendicular) to the arc. This can be interpreted by crustal thickening under the sedimentary arc and/or by velocity anisotropy with the maximum perpendicular to the arc. On the whole, our study supports the hypothesis that the Aegean region is a trench—island-arc—marginal-sea system.     

Deep structure under Yellowstone National Park U.S.A.: A continental “hot spot”

In order to understand the origin of long-lived loci of volcanism (sometimes called “hot spots”) and their possible role in global tectonic processes, it is essential to know their deep structure. Even though some work has been done on the crustal, upper-mantle, and deep-mantle structure under some of these “hot spots”, the picture is far from clear. In an attempt to study the structure under the Yellowstone National Park U.S.A., which is considered to be such a “hot spot”, we recorded teleseisms using 26 telemetered seismic stations and three groups of portable stations. The network was operated within a 150 km radius centered on the Yellowstone caldera, the major, Quaternary volcanic feature of the Yellowstone region. Teleseismic delays of about 1.5 sec are found inside the caldera, and the delays remain high over a 100 km wide area around the caldera. The spatial distribution and magnitude of the delays indicate the presence of a large body of low-velocity material with horizontal dimensions corresponding approximately to the caldera size (40 km × 80 km) near the surface and extending to a depth of 200–250 km under the caldera. Using ray-tracing and inversion techniques, it is estimated that the compressional velocity inside the anomalous body is lower than in the surrounding rock by about 15% in the upper crust and by 5% in the lower crust and upper mantle. It is postulated that the body is partly composed of molten rock with a high degree of partial melting at shallow depths and is responsible for the observed Yellowstone volcanism. The large size of the partially molten body, taken together with its location at the head of a 350 km zone of volcanic propagation along the axis of the Snake River Plain, indicates that the volcanism associated with Yellowstone has its origin below the lithosphere and is relatively stationary with respect to plate motion. Using our techniques, we are unable to detect any measurable velocity contrast in the mantle beneath the low-velocity body, and, hence, we are unable to determine whether the Yellowstone melting anomaly is associated with a deep heat source or with any deep phenomenon such as a convection plume, chemical plume, or gravitational anchor.     

Geochemistry of Cretaceous mafic rocks from the Lower Yangtze region, eastern China: Characteristics and evolution of the lithospheric mantle

Geochemical and isotopic data for Cretaceous mafic rocks (basalt, gabbro, and diorite) from the Lower Yangtze region, northern Yangtze block, constrain the evolution of the lithospheric mantle. The mafic rocks, separated into the northeast and southwest groups, are alkaline and evolved, with low Mg# values (44–58) and variable SiO₂ contents (47.6–57.4 wt%). Enriched LREEs, LILEs, and Pb, together with depleted Nb, Zr, and Ti, suggest that the mantle sources were metasomatized by slab-derived fluid/melt. All samples show high radiogenic ²⁰⁷Pb/²⁰⁴Pb(t) (15.41–15.65) and ²⁰⁸Pb/²⁰⁴Pb(t) (37.66–38.51) ratios at given ²⁰⁶Pb/²⁰⁴Pb(t) (17.65–19.00) ratios, consistent with the mantle sources having been metasomatized by ancient slab-derived material. Mafic rocks of the southwest group show enriched Sr–Nd isotopic characteristics, with ⁸⁷Sr/⁸⁶Sr(t) ranging from 0.7056 to 0.7071 and ε_Nd(t) ranging from −5.3 to −8.3, indicating an origin from enriched lithospheric mantle. Mafic rocks of the northeast group, which record ⁸⁷Sr/⁸⁶Sr(t) ratios of between 0.7044 and 0.7050 and ε_Nd(t) of −2.8 to −0.7, possibly formed by the mixing of melts from isotopically enriched lithospheric mantle and isotopically depleted asthenospheric mantle. Taking into consideration the geochemical and isotopic characteristics of Cretaceous mafic rocks, Cenozoic basalts, and basalt-hosted peridotite xenoliths from the Lower Yangtze region, we propose that an isotopically enriched, subduction-modified lithospheric mantle was replaced by or transformed into an isotopically depleted “oceanic-type” mantle. Such a process appears to have occurred in the eastern North China Craton as well as the eastern Yangtze block, probably in response to subduction of the paleo-Pacific plate beneath East Asia.     

Gravity anomalies and associated tectonic features over the Indian Peninsular Shield and adjoining ocean basins

A combined gravity map over the Indian Peninsular Shield (IPS) and adjoining oceans brings out well the inter-relationships between the older tectonic features of the continent and the adjoining younger oceanic features. The NW–SE, NE–SW and N–S Precambrian trends of the IPS are reflected in the structural trends of the Arabian Sea and the Bay of Bengal suggesting their probable reactivation. The Simple Bouguer anomaly map shows consistent increase in gravity value from the continent to the deep ocean basins, which is attributed to isostatic compensation due to variations in the crustal thickness. A crustal density model computed along a profile across this region suggests a thick crust of 35–40 km under the continent, which reduces to 22/20–24 km under the Bay of Bengal with thick sediments of 8–10 km underlain by crustal layers of density 2720 and 2900/2840 kg/m³. Large crustal thickness and trends of the gravity anomalies may suggest a transitional crust in the Bay of Bengal up to 150–200 km from the east coast. The crustal thickness under the Laxmi ridge and east of it in the Arabian Sea is 20 and 14 km, respectively, with 5–6 km thick Tertiary and Mesozoic sediments separated by a thin layer of Deccan Trap. Crustal layers of densities 2750 and 2950 kg/m³ underlie sediments. The crustal density model in this part of the Arabian Sea (east of Laxmi ridge) and the structural trends similar to the Indian Peninsular Shield suggest a continent–ocean transitional crust (COTC). The COTC may represent down dropped and submerged parts of the Indian crust evolved at the time of break-up along the west coast of India and passage of Reunion hotspot over India during late Cretaceous. The crustal model under this part also shows an underplated lower crust and a low density upper mantle, extending over the continent across the west coast of India, which appears to be related to the Deccan volcanism. The crustal thickness under the western Arabian Sea (west of the Laxmi ridge) reduces to 8–9 km with crustal layers of densities 2650 and 2870 kg/m³ representing an oceanic crust.     

Geophysical constraints on understanding the origin of the Illinois basin and its underlying crust

Interpretation of reprocessed seismic reflection profiles reveals three highly coherent, layered, unconformity-bounded sequences that overlie (or are incorporated within) the Proterozoic “granite–rhyolite province” beneath the Paleozoic Illinois basin and extend down into middle crustal depths. The sequences, which are situated in east–central Illinois and west–central Indiana, are bounded by strong, laterally continuous reflectors that are mappable over distances in excess of 200 km and are expressed as broad “basinal” packages that become areally more restricted with depth. Normal-fault reflector offsets progressively disrupt the sequences with depth along their outer margins. We interpret these sequences as being remnants of a Proterozoic rhyolitic caldera complex and/or rift episode related to the original thermal event that produced the granite–rhyolite province. The overall thickness and distribution of the sequences mimic closely those of the overlying Mt. Simon (Late Cambrian) clastic sediments and indicate that an episode of localized subsidence was underway before deposition of the post-Cambrian Illinois basin stratigraphic succession, which is centered farther south over the “New Madrid rift system” (i.e., Reelfoot rift and Rough Creek graben). The present configuration of the Illinois basin was therefore shaped by the cumulative effects of subsidence in two separate regions, the Proterozoic caldera complex and/or rift in east–central Illinois and west–central Indiana and the New Madrid rift system to the south. Filtered isostatic gravity and magnetic intensity data preclude a large mafic igneous component to the crust so that any Proterozoic volcanic or rift episode must not have tapped deeply or significantly into the lower crust or upper mantle during the heating event responsible for the granite–rhyolite.     

Lithospheric evolution of the Andean fold–thrust belt, Bolivia, and the origin of the central Andean plateau

We combine geological and geophysical data to develop a generalized model for the lithospheric evolution of the central Andean plateau between 18° and 20° S from Late Cretaceous to present. By integrating geophysical results of upper mantle structure, crustal thickness, and composition with recently published structural, stratigraphic, and thermochronologic data, we emphasize the importance of both the crust and upper mantle in the evolution of the central Andean plateau. Four key steps in the evolution of the Andean plateau are as follows. 1) Initiation of mountain building by 70 Ma suggested by the associated foreland basin depositional history. 2) Eastward jump of a narrow, early fold–thrust belt at 40 Ma through the eastward propagation of a 200–400-km-long basement thrust sheet. 3) Continued shortening within the Eastern Cordillera from 40 to 15 Ma, which thickened the crust and mantle and established the eastern boundary of the modern central Andean plateau. Removal of excess mantle through lithospheric delamination at the Eastern Cordillera–Altiplano boundary during the early Miocene appears necessary to accommodate underthrusting of the Brazilian shield. Replacement of mantle lithosphere by hot asthenosphere may have provided the heat source for a pulse of mafic volcanism in the Eastern Cordillera and Altiplano at 24–23 Ma, and further volcanism recorded by 12–7 Ma crustal ignimbrites. 4) After 20 Ma, deformation waned in the Eastern Cordillera and Interandean zone and began to be transferred into the Subandean zone. Long-term rates of shortening in the fold–thrust belt indicate that the average shortening rate has remained fairly constant (8–10 mm/year) through time with possible slowing (5–7 mm/year) in the last 15–20 myr. We suggest that Cenozoic deformation within the mantle lithosphere has been focused at the Eastern Cordillera–Altiplano boundary where the mantle most likely continues to be removed through piecemeal delamination.     

Seismic studies of the Brasília fold belt at the western border of the São Francisco Craton, Central Brazil, using receiver function, surface-wave dispersion and teleseismic tomography

The Tocantins Province in Central Brazil is composed of a series of SSW–NNE trending terranes of mainly Proterozoic ages, which stabilized in the Neoproterozoic in the final collision between the Amazon and São Francisco cratons. No previous information on crustal seismic properties was available for this region. Several broadband stations were used to study the regional patterns of crustal and upper mantle structure, extending the results of a recent E–W seismic refraction profile. Receiver functions and surface wave dispersion showed a thin crust (33–37 km) in the Neoproterozoic Magmatic Arc terrane. High average crustal Vp/Vs ratios (1.74–1.76) were consistently observed in this unit. The foreland domain of the Brasília foldbelt, on the other hand, is characterized by thicker crust (42–43 km). Low Vp/Vs ratios (1.70–1.72) were observed in the low-grade foreland fold and thrust zone of the Brasília belt adjacent to the São Francisco craton. Teleseismic P-wave tomography shows that the lithospheric upper mantle has lower velocities beneath the Magmatic Arc and Goiás Massif compared with the foreland zone of the belt and São Francisco craton. The variations in crustal thickness and upper mantle velocities observed with the broadband stations correlate well with the measurements along the seismic refraction profile. The integration of all seismic observations and gravity data indicates a strong lithospheric contrast between the Goiás Massif and the foreland domain of the Brasília belt, whereas little variation was found across the foldbelt/craton surface boundary. These results support the hypothesis that the Brasília foreland domain and the São Francisco craton were part of a larger São Francisco-Congo continental plate in the final collision with the Amazon plate.     

Rheologic properties of the lithosphere and applications to passive continental margins

Thermal and petrologic models of the crust and upper mantle are used for calculating effective viscosities on the basis of constant creep rates. Viscosity—depth models together with pressure—depth models are calculated for continental and oceanic blocks facing each other at continental margins. It is found from these “static models” that the overburden pressure in the lower crust and uppermost mantle causes a stress which is directed from the ocean to the continent. The generally low viscosity of 10²⁰–10²³ poise in this region should permit a creep process which could finally lead to a “silent” subduction. In the upper crust static stresses act in the opposite direction, i.e. from the continent to the ocean, favouring tension which could produce normal faulting in the continent. Differences between observations and the results obtained from the static models are attributed to dynamical forces.     

Late Pleistocene crustal uplift and gravity anomaly in the eastern part of Kyushu, Japan, and its geophysical implications

The Miyazaki Plain, eastern part of Kyushu, Japan, is characterized by both significant negative gravity anomalies and aseismic crustal uplifting (1 mm/year) in the Late Pleistocene and Holocene. We examine the relationship between these two phenomena, which may provide important constraints on the interaction between the collision and/or subduction of the Kyushu-Palau Ridge and the forearc. We estimate the mass deficiency below 11-km depth by using the gravity anomalies and P-wave velocity structure of the upper crust. The onset of the load accumulation, 0.5–0.4 Ma, is inferred from the movement of the fluvial terraces considering the tephrochronology. The loading history is assumed to be a linear function of time. We evaluate the crustal rebound by assuming a viscoelastic plate deformation with an underplating load existing at 20- or 30-km depth. The predicted crustal movement for models with a lithospheric (crustal) viscosity of 10²³–10²⁴ Pa s can explain the observed altitudes of the shoreline of the marine terraces formed at the Last Interglacial of about 125 kyr BP and the middle Holocene of 5–6 kyr BP. Although we cannot restrict the origin of the buoyant body, the subduction of the Kyushu-Palau Ridge, remnant arc associated with back-arc opening of the Shikoku Basin, may be related to the buoyancy for the uplifting region examined here. On the other hand, the buoyant body off the Miyazaki Plain probably plays an important role in the interaction between the subducting oceanic slab and the overriding forearc crust. Thus, the observed lateral variation of the interplate coupling on the convergent boundary along the Nankai Trough may be attributed to the existence of the buoyant body.     

Collisional plateaus

In the fifteen years since the importance of collisional plateaus with thickened continental crust began to be recognized as one of the inevitable consequences of the processes of plate tectonics, rapid progress in their understanding has come from studies of the world's only active terminal collision zones in the Himalayan-Tibetan and Turkish-Iranian plateaus.Ancient collisional plateaus are being recognized throughout the geological record (back to 3.8 Ga) from the occurrence of extensive areas (typically > 500,000 km²) of 8 kbar metamorphism in granulite facies or from the occurrence of extensive areas of higher level minimum-melt composition granite rocks whose isotopic signatures indicate reactivation of existing continental crust rather than addition of new crust from the mantle at the time of collision. Recognition of strike-slip faulting in the ancient collisional plateau areas indicates that “tectonic escape” may have been as important in the past as it is today.Earth may not be the only planet on which collisional plateaus are important. The highlands of Venus (approximately 7% of the surface with elevations over 1.5 km above mean planetary radius) can only exist as a result of crustal thickening, and not as a product of lithospheric thinning. Most of these highlands can be explained by models involving volcanic construction. However, the highest peaks, including Maxwell Montes, the highest mapped area of Venus rising over 10 km above mean planetary radius, require much greater crustal thickening to support them than can reasonably be explained by a volcanic mechanism. Geological features of Maxwell Montes inferred from radar images suggest some analogy between Maxwell Montes and the Tibetan plateau.It is somewhat paradoxical that extensional tectonics are commonly associated with continental collision, and that collision-related rifts may be the only sites where the uppermost layers of a collision-thickened crust are preserved from erosion. Extensional stress fields are generated during continental collision, primarily in areas associated with strike-slip faulting and “tectonic escape”. Additional extensional stresses are gravitationally generated associated with the topography and thickened crust in a collision zone. Tectonically thickened crust is particularly susceptible to rifting as its lithosphere is weak as a result of heating associated with magmatism. This lithosphere is also compositionally weak because of the relatively thick crust, dominated by a weak quartz rheology, and thin mantle lithosphere, dominated by a strong olivine rheology, in comparison with a lithosphere with a more normal crustal thickness. Thus, the common association of rifts and collision zones may be a consequence of both stresses generated during collision and modification of the lithosphere by collision.