Magnetometer array studies in India and the lithosphere

Two magnetometer array experiments were conducted in India during 1978–1980, under an Indo-Australian collaboration project, using 21 Australian three-component magnetometers of the Gough-Reitzel type. The first array study was made in the northwestern region covering the Aravallis, the Punjab, and the lesser Himalaya, while the second experiment was carried out in the southern peninsular shield area. Both these sets of geomagnetic deep sounding (GDS) observations yielded valuable results on the crustal and upper mantle structure in the two geologically and geophysically important regions of India.Geomagnetic induction patterns observed in northwest India have revealed a variety of electrical conductivity structures. The primary conductivity structure providing paths for induced currents is found to be striking at right-angles to the Himalayan Mountains. The conductivity structure is indicated to be a northward continuation of the Aravalli belt and, thus, suggesting the continuation of the Indian shield at depth into the base of the Himalayan foothills under the Ganga basin.The induction effects observed in the southern tip of peninsular India are by far the most complex geophysical phenomenon due to the simultaneous occurrence of the sea coast, the crustal and upper mantle conductivity anomalies between India and Sri Lanka under the sea, and the day-time equatorial electrojet as part of the external heterogeneous inducing field. It is further complicated by the existence of a conductive step, structure along the coastline at the Moho boundary and a “graben” structure in the Palk Strait, as revealed by the array observations.     

Structure of the crust and uppermost mantle beneath southern Finland revealed by analysis of local events registered by the SVEKALAPKO seismic array

In 1998–1999, a large-scale seismic array was deployed in Finland as a part of the EUROPROBE/SVEKALAPKO subproject, involving 14 European universities and research institutes. The objective of the project was to map the deep lithosphere structure and thickness beneath the Fennoscandian Shield by means of teleseismic events. In addition, about 580 local seismic events were registered during the data acquisition period. Among them, only eight local earthquakes were recorded, the rest being quarry blasts from mining sites in Russia, Finland, Estonia and Sweden. In this study, we present the analysis of the seismic wave field from the strongest local events registered by the majority of the stations of the SVEcofennian–KArelian–LAPland–KOla Transect (SVEKALAPKO) array with the aim of mapping the structure of the upper mantle beneath the array. For this purpose, we selected the events corresponding to a single source type and compared these recordings with those from wide-angle reflection and refraction experiments in the area to identify the regional phases. The record sections of selected events demonstrate strong reflections (PmP) from the Moho boundary. The refracted Pn phases can be seen as first arrivals at distances of about 200–400 km from the source. At offsets of about 400–800 km, phases reflected from inhomogeneities in the uppermost mantle (P1) and double reflections from the Moho boundary (PmPPmP) were recorded.Results from 2D forward ray trace modeling of reflected and refracted P-waves along four profile swathes from SVEKALAPKO stations demonstrate that the mantle reflections originate from two different groups of boundaries beneath the array: one group of phases arrive from subhorizontal and gently dipping reflectors below the Moho boundary at a depth of 70–90 km, while the other group are phases originating from a depth of 100 to 130 km. Based on the irregular character of the first group of reflections, their different spatial orientation and correlation with the Moho offsets, we interpret the boundaries of this group as relicts of ancient subduction and collision processes. The second group of reflections can be explained by a transition from mechanically strong to mechanically weak lithosphere.     

Integrated “plume winter” scenario for the double-phased extinction during the Paleozoic–Mesozoic transition: The G-LB and P-TB events from a Panthalassan perspective

The event across the Paleozoic–Mesozoic transition involved the greatest mass extinction in history together with other unique geologic phenomena of global context, such as the onset of Pangean rifting and the development of superanoxia. The detailed stratigraphic analyses on the Permo-Triassic sedimentary rocks documented a two-stepped nature both of the extinction and relevant global environmental changes at the Guadalupian–Lopingian (Middle and Upper Permian) boundary (G-LB, ca. 260 Ma) and at the Permo-Triassic boundary (P-TB, ca. 252 Ma), suggesting two independent triggers for the global catastrophe. Despite the entire loss of the Permian–Triassic ocean floors by successive subduction, some fragments of mid-oceanic rocks were accreted to and preserved along active continental margins. These provide particularly important dataset for deciphering the Permo-Triassic paleo-environments of the extensive superocean Panthalassa that occupied nearly two thirds of the Earth’s surface. The accreted deep-sea pelagic cherts recorded the double-phased remarkable faunal reorganization in radiolarians (major marine plankton in the Paleozoic) both across the G-LB and the P-TB, and the prolonged deep-sea anoxia (superanoxia) from the Late Permian to early Middle Triassic with a peak around the P-TB. In contrast, the accreted mid-oceanic paleo-atoll carbonates deposited on seamounts recorded clear double-phased changes of fusuline (representative Late Paleozoic shallow marine benthos) diversity and of negative shift of stable carbon isotope ratio at the G-LB and the P-TB, in addition to the Paleozoic minimum in ⁸⁷Sr/⁸⁶Sr isotope ratio in the Capitanian (Late Guadalupian) and the paleomagnetic Illawarra Reversal in the late Guadalupian. These bio-, chemo-, and magneto-stratigraphical signatures are concordant with those reported from the coeval shallow marine shelf sequences around Pangea. The mid-oceanic, deep- and shallow-water Permian records indicate that significant changes have appeared twice in the second half of the Permian in a global extent. It is emphasized here that everything geologically unusual started in the Late Guadalupian; i.e., (1) the first mass extinction, (2) onset of the superanoxia, (3) sea-level drop down to the Phanerozoic minimum, (4) onset of volatile fluctuation in carbon isotope ratio, 5) ⁸⁷Sr/⁸⁶Sr ratio of the Paleozoic minimum, (6) extensive felsic alkaline volcanism, and (7) Illawarra Reversal.The felsic alkaline volcanism and the concurrent formation of several large igneous provinces (LIPs) in the eastern Pangea suggest that the Permian biosphere was involved in severe volcanic hazards twice at the G-LB and the P-TB. This episodic magmatism was likely related to the activity of a mantle superplume that initially rifted Pangea. The supercontinent-dividing superplume branched into several secondary plumes in the mantle transition zone (410–660 km deep) beneath Pangea. These secondary plumes induced the decompressional melting of mantle peridotite and pre-existing Pangean crust to form several LIPs that likely caused a “plume winter” with global cooling by dust/aerosol screens in the stratosphere, gas poisoning, acid rain damage to surface vegetation etc. After the main eruption of plume-derived flood basalt, global warming (plume summer) took over cooling, delayed the recovery of biodiversity, and intensified the ocean stratification. It was repeated twice at the G-LB and P-TB.A unique geomagnetic episode called the Illawarra Reversal around the Wordian–Capitanian boundary (ca. 265 Ma) recorded the appearance of a large instability in the geomagnetic dipole in the Earth’s outer core. This rapid change was triggered likely by the episodic fall-down of a cold megalith (subducted oceanic slabs) from the upper mantle to the D″ layer above the 2900 km-deep core-mantle boundary, in tight association with the launching of a mantle superplume. The initial changes in the surface environment in the Capitanian, i.e., the Kamura cooling event and the first biodiversity decline, were probably led by the weakened geomagnetic intensity due to unstable dipole of geodynamo. Under the low geomagnetic intensity, the flux of galactic cosmic radiation increased to cause extensive cloud coverage over the planet. The resultant high albedo likely drove the Kamura cooling event that also triggered the unusually high productivity in the superocean and also the expansion of O₂ minimum zone to start the superanoxia.The “plume winter” scenario is integrated here to explain the “triple-double” during the Paleozoic–Mesozoic transition interval, i.e., double-phased cause, process, and consequence of the greatest global catastrophe in the Phanerozoic, in terms of mantle superplume activity that involved the whole Earth from the core to the surface biosphere.     

Upper mantle structure of the Northern Eurasia from peaceful nuclear explosion data

Several long-range seismic profiles were carried out in Russia with Peaceful Nuclear Explosions (PNE). The data from 25 PNEs recorded along these profiles were used to compile a 3-D upper mantle velocity model for the central part of the Northern Eurasia. 2-D crust and upper mantle models were also constructed for all profiles using a common methodology for wavefield interpretation. Five basic boundaries were traced over the study area: N1 boundary (velocity level, V = 8.35 km/s; depth interval, D = 60–130 km), N2 (V = 8.4 km/s; D = 100–140 km), L (V = 8.5 km/s; D = 180–240 km) and H (V = 8.6 km/s; D = 300–330 km) and structural maps were compiled for each boundary. Together these boundaries describe a 3-D upper mantle model for northern Eurasia. A map characterised the velocity distribution in the uppermost mantle down to a depth of 60 km is also presented. Mostly horizontal inhomogeneity is observed in the uppermost mantle, and the velocities range from the average 8.0–8.1 km/s to 8.3–8.4 km/s in some blocks of the Siberian Craton. At a depth of 100–200 km, the local high velocity blocks disappear and only three large anomalies are observed: lower velocities in West Siberia and higher velocities in the East-European platform and in the central part of the Siberian Craton. In contrast, the depths to the H boundary are greater beneath the craton and lower beneath in the West Siberian Platform. A correlation between tectonics, geophysical fields and crustal structure is observed. In general, the old and cold cratons have higher velocities in the mantle than the young platforms with higher heat flows.Structural peculiarities of the upper mantle are difficult to describe in form of classical lithosphere–asthenosphere system. The asthenosphere cannot be traced from the seismic data; in contrary the lithosphere is suggested to be rheologically stratified. All the lithospheric boundaries are not simple discontinuities, they are heterogeneous (thin layering) zones which generate multiphase reflections. Many of them may be a result of fluids concentrated at some critical P–T conditions which produce rheologically weak zones. The most visible rheological variations are observed at depths of around 100 and 250 km.     

Structure of the crust and uppermost mantle of the Yanshan Belt and adjacent regions at the northeastern boundary of the North China Craton from Rayleigh Wave Dispersion Analysis

We have used Rayleigh wave dispersion analysis and inversion to produce a high resolution S-wave velocity imaging profile of the crust and uppermost mantle structure beneath the northeastern boundary regions of the North China Craton (NCC). Using waveform data from 45 broadband NCISP stations, Rayleigh wave phase velocities were measured at periods from 10 to 48 s and utilized in subsequent inversions to solve for the S-wave velocity structure from 15 km down to 120 km depth. The inverted lower crust and uppermost mantle velocities, about 3.75 km/s and 4.3 km/s on average, are low compared with the global average. The Moho was constrained in the depth range of 30–40 km, indicating a typical crustal thickness along the profile. However, a thin lithosphere of no more than 100 km was imaged under a large part of the profile, decreasing to only ~ 60 km under the Inner Mongolian Axis (IMA) where an abnormally slow anomaly was observed below 60 km depth. The overall structural features of the study region resemble those of typical continental rift zones and are probably associated with the lithospheric reactivation and tectonic extension widespread in the eastern NCC during Mesozoic–Cenozoic time. Distinctly high velocities, up to ~ 4.6 km/s, were found immediately to the south of the IMA beneath the northern Yanshan Belt (YSB), extending down to > 100-km depth. The anomalous velocities are interpreted as the cratonic lithospheric lid of the region, which may have not been affected by the Mesozoic–Cenozoic deformation process as strongly as other regions in the eastern NCC. Based on our S-wave velocity structural image and other geophysical observations, we propose a possible lithosphere–asthenosphere interaction scenario at the northeastern boundary of the NCC. We speculate that significant undulations of the base of the lithosphere, which might have resulted from the uneven Mesozoic–Cenozoic lithospheric thinning, may induce mantle flows concentrating beneath the weak IMA zone. The relatively thick lithospheric lid in the northern YSB may serve as a tectonic barrier separating the on-craton and off-craton regions into different upper mantle convection systems at the present time.     

Possible contribution of lower mantle magnesiowüstite metallization into geomagnetic data

Phase change of dielectric magnesiowüstite in the lower mantle may leave signatures in geomagnetic records of the globally distributed array of observatories. We investigate theoretically which may be the contribution of magnesiowüstite metallization to geomagnetic data and how the variations of magnetic susceptibility associated with this phase change may influence the Earth's field. The modeling is performed using spherical harmonic analysis (SHA) of mantle electromagnetic (EM) responses in observatory geomagnetic data at periods of decades, 11 years, 1 year, and 27 days. The existence of a lower mantle conductor is checked against monthly means of real observatory records from 1920 through 2009 obtained by preliminary processing.     

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.     

On the nature of mantle heterogeneities and discontinuities: evidence from a very dense wide-angle shot record

A seismic survey with a receiver spacing of 50 m provided one of the most densely sampled wide-angle seismic reflection images of the lithosphere. This unique data set, recorded by an 18-km-long spread, reveals that at wide-angles the shallow subcrustal mantle features high amplitude reflectivity which contrasts with a lack of reflectivity at latter travel times. This change in the seismic signature is located at approximately 120–150 km depth, which correlates with the depth estimates of the lithosphere–asthenosphere boundary (LAB) of previous DSS studies. This seismic signature can be simulated by two-layer mantle model. Both layers with similar average velocities differ in their degree of heterogeneity. The shallow heterogeneous layer and the deeper and more homogeneous one correlate with the lithosphere and the asthenosphere, respectively. Studies involving surface outcrops of ultramafic massifs and mantle xenoliths infer that the upper mantle is a heterogeneous mixture of ultramafic rocks (lherzolites, harzburgites, pyroxenites, peridotites, dunites, and small amounts of eclogites). Laboratory measurements of physical properties of these mantle rocks indicate that compositional variations alone can account for the wide-angle reflectivity. A temperature increase would homogenize the mixture, decreasing the seismic reflection properties due to melting processes. It is proposed that this would take place below 120–150 km (1200 °C, the LAB).     

Geomagnetic induction responses of anisotropic conducting mantle

Phase change of dielectric magnesiowustite in the lower mantle may leave signatures in geomagnetic records of the globally distributed array of observatories. The related features appear in EM induction responses of lower mantle, which are studied theoretically. The surface EM field corresponding to a response of the earth with conductivity anisotropy in a mantle spherical layer is presented as the sum of the magnetic and electric modes. Equations for the fields of both modes and their relationship in a weakly anisotropic earth are obtained by the perturbation method. The two field modes are analyzed jointly and separately to characterize the conductivity tensor of the anisotropic lower mantle. The tensor elements corresponding to the tangential components of the field can be estimated from the magnetic mode alone recorded currently by the global network of geomagnetic observatories. For the tensor data to be complete, observatory data on lateral variations of the electric field are required in addition to three-component geomagnetic records.     

Nonlocal response functions: application to processing BEAR data

It has been suggested that regional array MTS data be processed with an algorithm based on the nonlocal response of laterally heterogeneous subsurface to an arbitrary electromagnetic excitation. The algorithm is tested with synthetic models which show that the inversion quality can be improved by accumulation of information through repeated measurements, at either changing or fixed array configurations. Besides the array configuration, the inversion quality depends on the selected size of the survey area. The choice of the area size defines the inversion conditions in terms of choice between the magnetic and electric mode for correlation of electromagnetic potentials. The algorithm has been applied to data acquired during the BEAR Project in the Fennoscandian shield. Preliminary results indicate a spatial correlation between conductivity minimums and Moho depth maximums. The apparent resistivity maps obtained with the suggested algorithm highlight the boundaries between the Lapland–Kola, Karelia, and Svecofennia tectonic provinces.     

The lithospheric density structure of the Eastern Alps

The three-dimensional (3D) lithospheric density structure of the Eastern Alps was investigated by integrating results from reflection seismics, receiver function analyses and tomography. The modelling was carried out with respect to the Bouguer gravity and the geoid undulations and emphasis were laid on the investigations of the importance of deep lithospheric features. Although the influence of inhomogeneities at the lithosphere–asthenosphere boundary on the potential field is not neglectable, they are overprinted by the response of the density contrast at the crust–mantle boundary and intra-crustal density anomalies. The uncertainties in the interpretations are in the same order of magnitude as the gravity field generated by the deep lithosphere.After including the deep lithospheric geometry from the tomographic model it is shown that full isostatic equilibrium is not achieved below the Eastern Alps. However, calculation of the isostatic lithospheric thickness shows two areas of lithospheric thickening along the central axis of the Eastern Alps with a transition zone below the area of the TRANSALP profile. This is in agreement with the tomographic model, which features a change in lithospheric subduction direction.     

Thin-sheet modelling of lithospheric deformation and surface mass transport

We study the effects of incorporating surface mass transport and the gravitational potential energy of both crust and lithospheric mantle to the viscous thin sheet approach. Recent 2D (cross-section) numerical models show that surface erosion and sediment transport can play a major role in shaping the large-scale deformation of the crust. In order to study these effects in 3D (planform view), we develop a numerical model in which both the dynamics of lithospheric deformation and surface processes are fully coupled. Deformation is calculated as a thin viscous layer with a vertically-averaged rheology and subjected to plane stresses. The coupled system of equations for momentum and energy conservation is solved numerically. This model accounts for the isostatic and potential-energy effects due to crustal and lithospheric thickness variations. The results show that the variations of gravitational potential energy due to the lateral changes of the lithosphere–asthenosphere boundary can modify the mode of deformation of the lithosphere. Surface processes, incorporated to the model via a diffusive transport equation, rather than just passively reacting to changes in topography, play an active role in controlling the lateral variations of the effective viscosity and hence of the deformation of the lithosphere.     

Double-planed deep seismic zone and anomalous structure in the upper mantle beneath northeastern Honshu (Japan)

Seismicity located by using the most recent data obtained from the high-gain seismograph network of Tohoku University shows that the deep seismic zone beneath northeastern Honshu, Japan, is composed of two thin planes which are parallel to each other and are 30–40 km apart. Focal mechanisms derived from the earthquakes in the upper plane are reverse-faulting, or, some of them, down-dip compression. As a contrast, those in the lower plane are down-dip extension. The location of the upper boundary of the descending lithospheric slab, inferred from the arrival-time difference between ScS and ScSp waves and from the travel-time anomaly of intermediate-depth earthquakes observed at the small-scale seismic array, coincides exactly with the upper plane of the double-planed deep seismic zone. Anelasticity (1/Q) structure of the upper mantle consists of three distinct zones: a high-Q (Qs− 1500) inclined lithospheric slab, an intennediate-Q (Qs−350) land-side mantle between the Pacific coast and the volcanic front, and a low-Q (Qs − 100) land-side mantle between the volcanic front and the coast of the Japan Sea.The evidence obtained here provides valuable information as to the definition of the type of mechanism producing the plate motion beneath island arcs.     

Permo–Triassic intraplate magmatism and rifting in Eurasia: implications for mantle plumes and mantle dynamics

At the transition from the Permian to the Triassic, Eurasia was the site of voluminous flood-basalt extrusion and rifting. Major flood-basalt provinces occur in the Tunguska, Taymyr, Kuznetsk, Verkhoyansk–Vilyuy and Pechora areas, as well as in the South Chinese Emeishen area. Contemporaneous rift systems developed in the West Siberian, South Kara Sea and Pyasina–Khatanga areas, on the Scythian platform and in the West European and Arctic–North Atlantic domain. At the Permo–Triassic transition, major extensional stresses affected apparently Eurasia, and possibly also Pangea, as evidenced by the development of new rift systems. Contemporaneous flood-basalt activity, inducing a global environmental crisis, is interpreted as related to the impingement of major mantle plumes on the base of the Eurasian lithosphere. Moreover, the Permo–Triassic transition coincided with a period of regional uplift and erosion and a low-stand in sea level. Permo–Triassic rifting and mantle plume activity occurred together with a major reorganization of plate boundaries and plate kinematics that marked the transition from the assembly of Pangea to its break-up. This plate reorganization was possibly associated with a reorganization of the global mantle convection system. On the base of the geological record, we recognize short-lived and long-lived plumes with a duration of magmatic activity of some 10–20 million years and 100–150 million years, respectively. The Permo–Triassic Siberian and Emeishan flood-basalt provinces are good examples of “short-lived” plumes, which contrast with such “long lived” plumes as those of Iceland and Hawaii. The global record indicates that mantle plume activity occurred episodically. Purely empirical considerations indicate that times of major mantle plume activity are associated with periods of global mantle convection reorganization during which thermally driven mantle convection is not fully able to facilitate the necessary heat transfer from the core of the Earth to its surface. In this respect, we distinguish between two geodynamically different scenarios for major plume activity. The major Permo–Triassic plume event followed the assembly Pangea and the detachment of deep-seated subduction slabs from the lithosphere. The Early–Middle Cretaceous major plume event, as well as the terminal–Cretaceous–Paleocene plume event, followed a sharp acceleration of global sea-floor spreading rates and the insertion of new subduction zone slabs deep into the mantle. We conclude that global plate kinematics, driven by mantle convection, have a bearing on the development of major mantle plumes and, to a degree, also on the pattern of related flood-basalt magmatism.     

High field strength element ratios in Archean basalts: a window to evolving sources of mantle plumes?

In terms of high field strength element ratios Nb/Th, Zr/Nb, Nb/Y and Zr/Y, most basalts from non-arc type Archean greenstones are similar to oceanic plateau basalts, suggestive of mantle plume sources. A large number of these basalts have ratios similar to primitive mantle composition. Perhaps the Archean mantle was less fractionated than at present and “primitive mantle” comprised much of the deep mantle and made a significant contribution to mantle plumes. The near absence of Archean greenstone basalts similar to NMORB in composition is also consistent with a relatively unfractionated mantle in which a shallow depleted source (DM) was volumetrically insignificant. The element ratios in basalts also indicate the existence of recycled components (HIMU, EM1, EM2) in the mantle by the Late Archean. This suggests that oceanic lithosphere was recycled into the deep mantle and became incorporated in some mantle plumes by the Late Archean. High field strength element ratios also indicate an important contribution of continental crust or/and subcontinental lithosphere to some non-arc Archean greenstone basalts. This implies that at least thin continental lithosphere was relatively widespread in the Archean.     

Lower lithospheric structure beneath the Trans-European Suture Zone from POLONAISE''97 seismic profiles

The large-scale POLONAISE'97 seismic experiment investigated the velocity structure of the lithosphere in the Trans-European Suture Zone (TESZ) region between the Precambrian East European Craton (EEC) and Palaeozoic Platform (PP). In the area of the Polish Basin, the P-wave velocity is very low (Vp <6.1 km/s) down to depths of 15–20 km, and the consolidated basement (Vp5.7–5.8 km/s) is 5–12 km deep. The thickness of the crust is 30 km beneath the Palaeozoic Platform, 40–45 km beneath the TESZ, and 40–50 km beneath the EEC. The compressional wave velocity of the sub-Moho mantle is >8.25 km/s in the Palaeozoic Platform and 8.1 km/s in the Precambrian Platform. Good quality record sections were obtained to the longest offsets of about 600 km from the shot points, with clear first arrivals and later phases of waves reflected/refracted in the lower lithosphere. Two-dimensional interpretation of the reversed system of travel times constrains a series of reflectors in the depth range of 50–90 km. A seismic reflector appears as a general feature at around 10 km depth below Moho in the area, independent of the actual depth to the Moho and sub-Moho seismic velocity. “Ringing reflections” are explained by relatively small-scale heterogeneities beneath the depth interval from 90 to 110 km. Qualitative interpretation of the observed wave field shows a differentiation of the reflectivity in the lower lithosphere. The seismic reflectivity of the uppermost mantle is stronger beneath the Palaeozoic Platform and TESZ than the East European Platform. The deepest interpreted seismic reflector with zone of high reflectivity may mark a change in upper mantle structure from an upper zone characterised by seismic scatterers of small vertical dimension to a lower zone with vertically larger seismic scatterers, possible caused by inclusions of partial melt.     

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.     

Andean subduction-related mantle xenoliths: Isotopic evidence of Sr–Nd decoupling during metasomatism

Sr–Nd isotopic analyses on some mantle xenolith samples from the Northern, Southern and Austral Andean volcanic zones exhibit radiogenic Sr enrichment without dramatic changing of the Nd isotopic composition. This anomalous effect (Sr–Nd decoupling) makes these samples plot displaced to the right side of the “mantle array” trend (here called the “MORB–OIB–BSE trend”) in the ⁸⁷Sr/⁸⁶Sr vs. ¹⁴³Nd/¹⁴⁴Nd isotopic diagram. Such behavior reflects processes that took place in the mantle and can be related to: i) the mixture of a depleted mantle and an enriched source (enriched mantle II—EMII); ii) the mixture of a depleted mantle and a mixture of mantle-derived and slab-derived melts; and iii) a chromatographic process that occurs during the percolation of a metasomatic agent through the mantle.     

Petrogenesis and tectonic significance of Jurassic IAT magma types in the Hellenide ophiolites as deduced from the Rhodiani ophiolites (Pelagonian zone, Greece)

The Rhodiani ophiolites are represented by two tectonically superimposed ophiolitic units: the “lower” Ultramafic unit and the “upper” Volcanic unit, both bearing calcareous sedimentary covers. The Ultramafic unit consists of mantle harzburgites with dunite pods and chromitite ores, and represents the typical mantle section of supra-subduction zone (SSZ) settings. The Volcanic unit is represented by a sheeted dyke complex overlain by a pillow and massive lava sequence, both including basalts, basaltic andesites, andesites, and dacites. Chemically, the Volcanic unit displays low-Ti affinity typical of island arc tholeiite (IAT) ophiolitic series from SSZ settings, having, as most distinctive chemical features, low Ti/V ratios (< 20) and depletion in high field strength elements and light rare earth elements.The rare earth element and incompatible element composition of the more primitive basaltic andesites from the Rhodiani ophiolites can be successfully reproduced with about 15% non-modal fractional melting of depleted lherzolites, which are very common in the Hellenide ophiolites. The calculated residua correspond to the depleted harzburgites found in the Rhodiani and Othrys ophiolites. Both field and chemical evidence suggest that the whole sequence of the Rhodiani Volcanic unit (from basalt to dacite) originated by low-pressure fractional crystallization under partially open-system conditions. The modelling of mantle source, melt generation, and mantle residua carried out in this paper provides new constraints for the tectono-magmatic evolution of the Mirdita–Pindos oceanic basin.     

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.