Tectonic inheritance of the Indian Shield: New insights from its elastic thickness structure

A new evaluation of the elastic thickness (Te) structure of the Indian Shield, derived from isotropic fan wavelet methodology, documents spatial variations of lithospheric deformation in different tectonic provinces correlated with episodic tectono-thermal events. The Te variations corroborated by shear velocity, crustal thickness, and seismogenic thickness reveal the heterogeneous rheology of the Indian lithosphere. The thinned, attenuated lithosphere beneath Peninsular India is considered to be the reason for its mechanically weak strength (<30 km), where a decoupled crust-mantle rheology under different surface/subsurface loading structures may explain the prominent low Te patterns. The arcuate Te structure of the Western Dharwar province and a NNE-trending band of low Te anomaly in the Southern Granulite Terrane are intriguing patterns. The average Te values (40-50 km) of the Central Indian Tectonic Zone, the Bastar Craton, and the northern Eastern Ghats Mobile Belt are suggestive of old, stable, Indian lithosphere, which was not affected by any major tectono-thermal events after cratonic stabilization. We propose that the anomalously high Te (60-85 km) and high S-wave velocity zone to the north of the Narmada-Son Lineament, mainly in NW Himalaya, and the northern Aravalli and Bundelkhand Cratons, suggest that Archean lithosphere characterized by a high velocity mantle keel supports the orogenic topographic loads in/near the Himalaya. The Te map clearly segments the volcanic provinces of the Indian Shield, where the signatures of the Reunion, Marion, and Kerguelen hotspots are indicated by significantly low Te patterns that correlate with plume- and rift-related thermal and mechanical rejuvenation, magmatic underplating, and crustal necking. The correlations between Te variations and the occurrence of seismicity over seismically active zones reveal different causal relationships, which led to the current seismogenic zonation of the Indian Shield. (C) 2013 Elsevier B.V. All rights reserved.

Archean tectonics and crustal evolution of the Biligiri Rangan Block, southern India

The Southern Granulite Terrain in India is a collage of crustal blocks ranging in age from Archean to Neoproterozoic. This study investigate the tectonic evolution of one of the northernmost block- the Biligiri Block (BRB) through a multidisciplinary approach involving field investigation, petrographic studies, LA-ICPMS zircon U-Pb geochronology, Hf isotopic analyses, metamorphic P-T phase diagram computations, and crustal thickness modeling. The garnet bearing quartzofeldspathic gneiss from the central BRB preserve Mesoarchean magmatic zircons with ages between 3207 and 2806 Ma and positive epsilon Hf value (+2.7) which possibly indicates vestiges of a Mesoarchean primitive continental crust. The occurrence of quartzite-iron formation intercalation as well as ultramafic lenses along the western boundary of the BRB is interpreted to indicate that the Kollegal structural lineament is a possible paleo-suture. Phase diagram computation of a metagabbro from the southwestern periphery of the Kollegal suture zone reveals high-pressure (similar to 18.5 kbar) and medium-temperature (similar to 840 degrees C) metamorphism, likely during eastward subduction of the Western Dharwar oceanic crust beneath the Mesoarchean BRB. In the model presented here, slab subduction, melting and underplating processes generated arc magmatism and subsequent charnockitization within the BRB between ca. 2650 Ma and ca. 2498 Ma. These results thus reveal Meso- to Neoarchean tectonic evolution of the BRB. The spatial variation of crustal thickness, derived from flexure inversion technique, provides additional constraints on the tectonic linkage of the BRB with its surrounding terrains. In conjunction with published data, the Moyar and the Kollegal suture zones are considered to mark the trace of ocean closure along which the Nilgiri and Biligiri Rangan Blocks accreted on to the Western Dharwar Craton. (C) 2016 Elsevier B.V. All rights reserved.

Genetic algorithm inversion of the average 1D crustal structure using local and regional earthquakes

Knowing the best 1D model of the crustal and upper mantle structure is useful not only for routine hypocenter determination, but also for linearized joint inversions of hypocenters and 3D crustal structure, where a good choice of the initial model can be very important. Here, we tested the combination of a simple GA inversion with the widely used HYPO71 program to find the best three-layer model (upper crust, lower crust, and upper mantle) by minimizing the overall P- and S-arrival residuals, using local and regional earthquakes in two areas of the Brazilian shield. Results from the Tocantins Province (Central Brazil) and the southern border of the Sao Francisco craton (SE Brazil) indicated an average crustal thickness of 38 and 43 km, respectively, consistent with previous estimates from receiver functions and seismic refraction lines. The GA + HYPO71 inversion produced correct Vp/Vs ratios (1.73 and 1.71, respectively), as expected from Wadati diagrams. Tests with synthetic data showed that the method is robust for the crustal thickness, Pn velocity, and Vp/Vs ratio when using events with distance up to about 400 km, despite the small number of events available (7 and 22, respectively). The velocities of the upper and lower crusts, however, are less well constrained. Interestingly, in the Tocantins Province, the GA + HYPO71 inversion showed a secondary solution (local minimum) for the average crustal thickness, besides the global minimum solution, which was caused by the existence of two distinct domains in the Central Brazil with very different crustal thicknesses. (C) 2010 Elsevier Ltd. All rights reserved.

Estimativa de espessura crustal na Província Borborema (NE/Brasil) através de função do recpetor

Crustal thickness and VP/VS estimates are essential to the studies of subsurface geological structures and also to the understanding of the regional tectonic evolution of a given area. In this dissertation, we use the Langston´s (1979) Receiver Function Method using teleseismic events reaching the seismographic station with angles close to the vertical. In this method, the information of the geologic structures close to the station is isolated so that effects related to the instrument response and source mechanics are not present. The resulting time series obtained after the deconvolution between horizontal components contains the larger amplitude referring to the P arrival, followed by smaller arrival caused by the reverberation and conversion of the P-wave at the base of the crust. We also used the HK-Stacking after Zhu & Kanamori (2000) to obtain crustal thickness and Vp/VS estimates. This method works stacking receiver functions so that the best estimates of crustal thickness and Vp/VS are found when the direct P, the Ps wave and the first multiple are coherently stacked. We used five broadband seismographic stations distributed over the Borborema Province, NE Brazil. Crustal thickness and Vp/VS estimates are consistent with the crust-mantle interface obtained using gravity data. We also identified crutal thickening in the NW portion of the province, close to Sobral/CE. Towards the center-north portion of the province, there is an evident crustal thinning which coincides with a geological feature consisting of an alignment of sedimentary basins known as the Cariris-Potiguar trend. Towards the NE portion of the province, in Solânea/PB and Agrestina/PE regions, occurs a crustal thickening and a systematic increase in the VP/VS values which suggest the presence of mafic rocks in the lower crust also consistent with the hypothesis of underplating in the region

Estrutura crustal e mantélica da província Borborema através de funções do receptor e dispersão de ondas superficiais

The Borborema Province, located in northeastern Brazil, has a basement of Precambrian age and a tectonic framework structured at the Neoproterozoic (740-560 Ma). After separation between South America and Africa during the Mesozoic, a rift system was formed, giving rise to a number of marginal and inland basins in the Province. After continental breakup, episodes of volcanism and uplift characterized the evolution of the Province. Plateau uplift was initially related to magmatic underplating of mafic material at the base of the crust, perhaps related to the generation of young continental plugs (45-7 Ma) along the Macau-Queimadas Alignment (MQA), due to a small-scale convection at the continental edge. The goal of this study is to investigate the causes of intra-plate uplift and its relationship to MQA volcanism, by using broadband seismology and integrating our results with independent geophysical and geological studies in the Borborema Province. The investigation of the deep structure of the Province with broadband seismic data includes receiver functions and surface-wave dispersion tomography. Both the receiver functions and surface-wave dispersion tomography are methods that use teleseismic events and allow to develop estimates of crustal parameters such as crustal thickness, Vp/Vs ratio, and S-velocity structure. The seismograms used for the receiver function work were obtained from 52 stations in Northeast Brazil: 16 broadband stations from the RSISNE network (Rede Sismográfica do Nordeste do Brasil), and 21 short-period and 6 broadband stations from the INCT-ET network (Instituto Nacional de Ciência e Tecnologia – Estudos Tectônicos). These results add signifi- cantly to previous datasets collected at individual stations in the Province, which include station RCBR (GSN - Global Seismic Network), stations CAUB and AGBL (Brazilian Lithosphere Seismic Project IAG/USP), and 6 other broadband stations that were part of the Projeto Milênio - Estudos geofísicos e tectônicos na Província Borborema/CNPq. For the surface-wave vii tomography, seismograms recorde at 22 broadband stations were utilized: 16 broadband stations from the RSISNE network and 6 broadband stations from the Milênio project. The new constraints developed in this work include: (i) estimates of crustal thickness and bulk Vp/Vs ratio for each station using receiver functions; (ii) new measurements of surfassewave group velocity, which were integrated to existing measurementes from a continental-scale tomography for South America, and (iii) S-wave velocity models (1D) at various locations in the Borborema Province, developed through the simultaneous inversion of receiver functions and surface-wave dispersion velocities. The results display S-wave velocity structure down to the base of the crust that are consistent with the presence of a 5-7.5 km thick mafic layer. The mafic layer was observed only in the southern portion of the Plateau and absent in its northern portion. Another important observation is that our models divide the plateau into a region of thin crust (northern Plateau) and a region of thick crust (southern Plateau), confirming results from independent refraction surveys and receiver function analyses. Existing models of plateau uplift, nonetheless, cannot explain all the new observations. It is proposed that during the Brazilian orogeny a layer of preexisting mafic material was delaminated, as a whole or in part, from the original Brasiliano crust. Partial delamination would have happened in the southern portion of the plateau, where independent studies found evidence of a more resistant rheology. During Mesozoic rifting, thinning of the crust around the southern Plateau would have formed the marginal basins and the Sertaneja depression, which would have included the northern part of the Plateau. In the Cenozoic, uplift of the northern Plateau would have occurred, resulting in a northern Plateau without mafic material at the base of the crust and a southern Plateau with partially delaminated mafic layer.

Estimation of the Crustal Bulk Properties Beneath Mainland Portugal from P-Wave Teleseismic Receiver Functions

In this work, we present results from teleseismic P-wave receiver functions (PRFs) obtained in Portugal, Western Iberia. A dense seismic station deployment conducted between 2010 and 2012, in the scope of the WILAS project and covering the entire country, allowed the most spatially extensive probing on the bulk crustal seismic properties of Portugal up to date. The application of the H-kappa stacking algorithm to the PRFs enabled us to estimate the crustal thickness (H) and the average crustal ratio of the P- and S-waves velocities V (p)/V (s) (kappa) for the region. Observations of Moho conversions indicate that this interface is relatively smooth with the crustal thickness ranging between 24 and 34 km, with an average of 30 km. The highest V (p)/V (s) values are found on the Mesozoic-Cenozoic crust beneath the western and southern coastal domain of Portugal, whereas the lowest values correspond to Palaeozoic crust underlying the remaining part of the subject area. An average V (p)/V (s) is found to be 1.72, ranging 1.63-1.86 across the study area, indicating a predominantly felsic composition. Overall, we systematically observe a decrease of V (p)/V (s) with increasing crustal thickness. Taken as a whole, our results indicate a clear distinction between the geological zones of the Variscan Iberian Massif in Portugal, the overall shape of the anomalies conditioned by the shape of the Ibero-Armorican Arc, and associated Late Paleozoic suture zones, and the Meso-Cenozoic basin associated with Atlantic rifting stages. Thickened crust (30-34 km) across the studied region may be inherited from continental collision during the Paleozoic Variscan orogeny. An anomalous crustal thinning to around 28 km is observed beneath the central part of the Central Iberian Zone and the eastern part of South Portuguese Zone.

Deep Crustal Structure of the North-West African Margin from Combined Wide-angle and Reflection Seismic Data (MIRROR Seismic Survey)

The structure of the Moroccan and Nova Scotia conjugate rifted margins is of key importance for understanding the Mesozoic break-up and evolution of the northern central Atlantic Ocean basin. Seven combined multichannel reflection (MCS) and wide-angle seismic (OBS) data profiles were acquired along the Atlantic Moroccan margin between the latitudes of 31.5° and 33° N during the MIRROR seismic survey in 2011, in order to image the transition from continental to oceanic crust, to study the variation in crustal structure and to characterize the crust under the West African Coast Magnetic Anomaly (WACMA). The data were modeled using a forward modeling approach. The final models image crustal thinning from 36 km thickness below the continent to approximately 8 km in the oceanic domain. A 100 km wide zone characterized by rough basement topography and high seismic velocities up to 7.4 km/s in the lower crust is observed westward of the West African Coast Magnetic Anomaly. No basin underlain by continental crust has been imaged in this region, as has been identified north of our study area. Comparison to the conjugate Nova Scotian margin shows a similar continental crustal thickness and layer geometry, and the existence of exhumed and serpentinized upper mantle material on the Canadian side only. The oceanic crustal thickness is lower on the Canadian margin.

The Southeast Indian Ridge between 88°E and 118°E: Gravity anomalies and crustal accretion at intermediate spreading rates

Although slow spreading ridges characterized by a deep axial valley and fast spreading ridges characterized by an axial bathymetric high have been extensively studied, the transition between these two modes of axial morphology is not well understood. We conducted a geophysical-survey of the intermediate spreading rate Southeast Indian Ridge between 88 degrees E and 118 degrees E, a 2300-km-long section of the ridge located between the Amsterdam hot spot and the Australian-Antarctic Discordance where satellite gravity data suggest that the Southeast Indian Ridge (SEIR) undergoes a change from an axial high in the west to an axial valley in the east. A basic change in axial morphology is found near 103 degrees 30'E in the shipboard data; the axis to the west is marked by an axial high, while a valley is found to the east. Although a well-developed axial high, characteristic of the East Pacific Rise (EPR), is occasionally present, the more common observation is a rifted high that is lower and pervasively faulted, sometimes with significant (> 50 m throw) faults within a kilometer of the axis. A shallow axial valley (< 700 m deep) is observed from 104 degrees E to 114 degrees E with a sudden change to a deep (>1200 m deep) valley across a transform at 114 degrees E. The changes in axial morphology along the SEIR are accompanied by a 500 m increase in near-axis ridge flank depth from 2800 m near 88 degrees E to 3300 m near 114 degrees E and by a 50 mGal increase in the regional level of mantle Bouguer gravity anomalies over the same distance, The regional changes in depth and mantle Bouguer anomaly (MBA) gravity can be both explained by a 1.7-2.4 km change in crustal thickness or by a mantle temperature change of 50 degrees C-90 degrees C. In reality, melt supply (crustal thickness) and mantle temperature are linked, so that changes in both may occur simultaneously and these estimates serve as upper bounds. The along-axis MBA gradient is not uniform. Pronounced steps in the regional level of the MBA gravity occur at 103 degrees 30'E-104 degrees E and at 114 degrees E-116 degrees E and correspond to the changes in the nature of the axial morphology and in the amplitude of abyssal hill morphology suggesting that the different forms of morphology do not grade into each other but rather represent distinctly different forms of axial (s)tructure and tectonics with a sharp transition between them. The change from an axial high to an axial valley requires a threshold effect in which the strength of the lithosphere changes quickly. The presence or absence of a quasi-steady state magma chamber may provide such a mechanism. The different forms of axial morphology are also associated with different intrasegment MBA gravity patterns. Segments with an axial high have an MBA low located at a depth minimum near the center of the segment, At EPR-like segments, the MBA low is about 10 mGal with along-axis gradients of 0.15-0.25 mGal/km, similar to those observed at the EPR, Rifted highs have a shallower low and lower gradients suggesting an attenuated composite magma chamber and a reduced and perhaps episodic melt supply. Segments with a shallow axial valley have very flat along-axis MBA profiles with little correspondence between axial depth and axial MBA gravity.

25 million years to break a continent : early to middle Miocene rifting and syn-extensional magmatism in the southern Gulf of California

Cenozoic extension in western Mexico has been divided into two episodes separated by the change from convergence to oblique divergence at the plate boundary. The Gulf Extensional Province is thought to have started once subduction ended at ~12.5 Ma whereas early extension is classified as Basin and Range. Mid-Miocene volcanism of the Comondú group has been considered as a subduction-related arc, whereas post ~12.5 Ma volcanism would be extension-related. Our new integration of the continental onshore and offshore geology of the south-east Gulf region, backed by tens of Ar-Ar and U-Pb ages and geochemical studies, document an early-mid Miocene rifting and extension-related bimodal to andesitic magmatism prior to subduction termination. Between ~21 and 11 Ma a system of NNW-SSE high-angle extensional faults rifted the western side of the Sierra Madre Occidental (SMO) ignimbrite plateau. In Nayarit, rhyolitic domes and some basalts were emplaced along this extensional belt at 18-17 Ma. These rocks show strong antecrystic inheritance but an absence of Mesozoic and older xenocrysts, suggesting a genesis in the mid-upper crust triggered by extension-induced basaltic influx. In Sinaloa, large grabens were floored by huge dome complexes at ~21-17 Ma and filled by continental sediments with interlayered basalts dated at 15 Ma. Mid-Miocene volcanism, including the largely volcaniclastic Comondú strata in Baja California, was thus emplaced in rift basins and appears associated to decompression melting rather than subduction. Along the coast, flat-lying basaltic lava flows dated at 11-10 Ma are exposed just above the present sea level. Here crustal thickness is 25-20 Km, almost half that in the core of the SMO, implying significant lithosphere stretching before ~11 Ma. This mafic pulse, with relatively high Ti but still clear Nb-Ta negative spikes, may be related to the detachment of the lower part of the subducted slab, allowing asthenosphere to flow into parts of the mantle previously fluxed by subduction fluids. Very uniform OIB-like lavas appear in late Pliocene and Pleistocene, only 18 m.y. after the onset of rifting and ~9 m.y. after the end of subduction. Our study shows that rifting began much earlier than Late Miocene and progressively overwhelmed subduction in generating magmatism.

Late Oligocene to middle Miocene rifting and syn-extensional magmatism in the southwestern Sierra Madre occidental, Mexico: the beginning of the Gulf of California rift

Although Basin and Range style extension affected several areas of western Mexico since the Late Eocene, extension in the Gulf of California region (the Gulf Extensional Province GEP) is thought to have started as subduction waned and ended at ~14 12.5 Ma. A general consensus also exists in considering the mid Miocene Comondú group as a suprasubduction volcanic arc. Our new integration of the geology of the south east Gulf region, backed by 43 new Ar Ar and U Pb mineral ages and geochemical studies, document a widespread phase of extension in the southern GEP between latest Oligocene and Early Miocene that subsequently focused in the region of the future Gulf in the Middle Miocene. Upper Oligocene to Lower Miocene rocks across the southern Sierra Madre Occidental (SMO)(northern Nayarit and southern Sinaloa) were affected by major ~N S to NNW striking normal faults prior to ~21 Ma. Then, between ~21 and 11 Ma, a system of NNW-SSE high angle extensional faults continued extending the southwestern side of the SMO. Rhyolitic domes, shallow intrusive bodies, and lesser basalts were emplaced along this extensional belt at 20-17 Ma. In northern Sinaloa, large grabens were floored by huge dome complexes at ~21-17 Ma and filled by continental sediments with interlayered basalts dated at 15-14 Ma, a setting and timing very similar to Sonora. Early to Middle Miocene volcanism, including the largely volcaniclastic Comondú strata in Baja California Sur, was thus emplaced in rift basins and was likely associated to decompression melting of upper mantle (inducing crustal partial melting) rather than to fluxing by fluids from the young subducting plate. Along the Nayarit and Sinaloa coast, flatlying basaltic lava flows dated at 11-10 Ma are exposed just above the present sea level. Here, crustal thickness is almost half that in the unextended core of the SMO, implying significant lithosphere stretching before ~11 Ma. Our study shows that rifting began much earlier than Late Miocene and provided a fundamental control on the style and composition of volcanism from at least 30 Ma. We envision a sustained period of lithospheric stretching and magmatism during which the pace and breadth of extension changed at ~20-18 Ma to be narrower and likely more rapid, and again at ~12.5 Ma, when the kinematics of rifting became more oblique.

New insights into magmatism and changing rift styles in the evolution of the Gulf of California

The Gulf of California (GoC) has been an important focus site for understanding the spatial and temporal evolution of rifts, with recent studies concluding: 1) rapid crustal rupturing within 10 Myrs; 2) surprisingly abrupt variations in rifting style and magmatism with apparently wide magma-poor and narrow, magmatic rift segments; and 3) that high sedimentation rates may promote switching from wide to narrow rift modes or thermally blanket the crust to enhance rift magmatism. Critical to these conclusions is the onset of rifting at~12 Ma following the cessation of subduction. New field-based volcanostratigraphic and geochronologic studies along the southeastern GoC margin reveal Early Miocene (~25-18 Ma) bimodal volcanism in wide rifting mode (~400 km width), followed by a mid-Miocene (~18-12 Ma) phase of dominantly intermediate composition magmatism in and around the nascent GoC with lavas/domes often emplaced into actively subsiding basins, but contemporaneous with bimodal volcanism regionally. Flat-lying intraplate basaltic lava fields emplaced ~12-10 Ma along the GoC east coast abut tilted blocks of ~20 Ma ignimbrites onshore, and also occur offshore. The reduction in crustal thickness from ~55 to 20 km along the eastern GoC edge must have been largely achieved by 12 Ma. Extension has demonstrably began earlier than previously thought, downplaying rapid rifting and any thermal effects from <6 Ma sedimentation. New age data from onshore indicate significant structurally controlled corridors of magmatism during 18-12 Ma extension in apparently magma-poor rift segments, and this magmatism temporally coincides with the switch from wide to narrow rifting.

Late Oligocene to Middle Miocene rifting and synextensional magmatism in the southwestern Sierra Madre Occidental, Mexico : the beginning of the Gulf of California rift

Although Basin and Range–style extension affected large areas of western Mexico after the Late Eocene, most consider that extension in the Gulf of California region began as subduction waned and ended ca. 14–12.5 Ma. A general consensus also exists in considering Early and Middle Miocene volcanism of the Sierra Madre Occidental and Comondú Group as subduction related, whereas volcanism after ca. 12.5 Ma is extension related. Here we present a new regional geologic study of the eastern Gulf of California margin in the states of Nayarit and Sinaloa, Mexico, backed by 43 new Ar-Ar and U-Pb mineral ages, and geochemical data that document an earlier widespread phase of extension. This extension across the southern and central Gulf Extensional Province began between Late Oligocene and Early Miocene time, but was focused in the region of the future Gulf of California in the Middle Miocene. Late Oligocene to Early Miocene rocks across northern Nayarit and southern Sinaloa were affected by major approximately north-south– to north-northwest– striking normal faults prior to ca. 21 Ma. Between ca. 21 and 11 Ma, a system of north-northwest–south-southeast high angle extensional faults continued extending the southwestern side of the Sierra Madre Occidental. Rhyolitic domes, shallow intrusive bodies, and lesser basalts were emplaced along this extensional belt at 20–17 Ma. Rhyolitic rocks, in particular the domes and lavas, often show strong antecrystic inheritance but only a few Mesozoic or older xenocrysts, suggesting silicic magma generation in the mid-upper crust triggered by an extension induced basaltic infl ux. In northern Sinaloa, large grabens were occupied by huge volcanic dome complexes ca. 21–17 Ma and filled by continental sediments with interlayered basalts dated as 15–14 Ma, a stratigraphy and timing very similar to those found in central Sonora (northeastern Gulf of California margin). Early to Middle Miocene volcanism occurred thus in rift basins, and was likely associated with decompression melting of upper mantle (inducing crustal partial melting) rather than with fluxing by fluids from the young and slow subducting microplates. Along the eastern side of the Gulf of California coast, from Farallón de San Ignacio island offshore Los Mochis, Sinaloa, to San Blas, Nayarit, a strike distance of >700 km, flat lying basaltic lavas dated as ca. 11.5–10 Ma are exposed just above the present sea level. Here crustal thickness is almost half that in the unextended core of the adjacent Sierra Madre Occidental, implying signifi cant lithosphere stretching before ca. 11 Ma. This mafic pulse, with subdued Nb-Ta negative spikes, may be related to the detachment of the lower part of the subducted slab, allowing an upward asthenospheric flow into an upper mantle previously modified by fluid fluxes related to past subduction. Widespread eruption of very uniform oceanic island basalt–like lavas occurred by the late Pliocene and Pleistocene, only 20 m.y. after the onset of rifting and ~9 m.y. after the end of subduction, implying that preexisting subduction-modified mantle had now become isolated from melt source regions. Our study shows that rifting across the southern-central Gulf Extensional Province began much earlier than the Late Miocene and provided a fundamental control on the style and composition of volcanism from at least 30 Ma. We envision a sustained period of lithospheric stretching and magmatism during which the pace and breadth of extension changed ca. 20–18 Ma to be narrower, and again after ca. 12.5 Ma, when the kinematics of rifting became more oblique.

Suomen maankuoren kivilajikoostumuksen tulkinta seismisistä aineistoista

Aim of this study is to investigate composition of the crust in Finland using seismic wide-angle velocity models and laboratory measurements on P- and S-wave velocities of different rock types. The velocities adopted from wide-angle velocity models were compared with laboratory velocities of different rock types corrected for the crustal PT conditions in the study area. The wide-angle velocity models indicate that the P-wave velocity does not only increase step-wise at boundaries of major crustal layers, but there is also gradual increase of velocity within the layers. On the other hand, the laboratory measurements of velocities indicate that no single rock type is able to provide the gradual downward increasing trends. Thus, there must be gradual vertical changes in rock composition. The downward increase of velocities indicates that the composition of the crust becomes gradually more mafic with increasing depth. Even though single rock types cannot simulate the wide-angle model velocities, it can be done with a mixture of rock types. There are a large number of rock type mixtures giving the correct P-wave velocities. Therefore, the inverse solution of rock types and their proportions from velocities is a non-unique problem if only P-wave velocities is available. Amount of the possible rock type mixtures can be limitted using S-wave velocities, reflection seismic results and other geological and geophysical results of the study area. Crustal model FINMIX-2 is presented in this study and it suggest that the crustal velocity profiles can be simulated with rock type mixtures, where the upper crust consists of felsic gneisses and granitic-granodioritic rocks with a minor contribution of quartzite, amphibolite and diabase. In the middle crust the amphibolite proportion increases. The lower crust consists of tonalitic gneiss, mafic garnet granulite, hornblendite, pyroxenite and minor mafic eclogite. This composition model is in agreement with deep crustal kimberlite-hosted xenolith data in eastern Finland and reflectivity of the FIRE (Finnish Reflection Experiment). According to FINMIX-2 model the Moho is deeper and the crustal composition is a more mafic than an average global continental model would suggest. Composition models of southern Finland are quite similar than FINMIX-2 model. However, there are minor differencies between the models, which indicates areal differences of composition. Models of northern Finland shows that the crustal thickness is smaller than southern Finland and composition of the upper crust is different. Density profiles calculated from the lithological models suggest that there is practically no density contrast at Moho in areas of the high-velocity lower crust. This implies that crustal thickness in the central Fennoscandian Shield may have been controlled by the densities of the lower crustal and upper mantle rocks.

青藏高原及邻区面波频散反演--岩石圈-软流层体系结构研究

The Qinghai-Tibet Plateau lies in the place of the continent-continent collision between Indian and Eurasian plates. Because of their interaction the shallow and deep structures are very complicated. The force system forming the tectonic patterns and driving tectonic movements is effected together by the deep part of the lithosphere and the asthenosphere. It is important to study the 3-D velocity structures, the spheres and layers structures, material properties and states of the lithosphere and the asthenosphere for getting knowledge of their formation and evolution, dynamic process, layers coupling and exchange of material and energy. Based on the Rayleigh wave dispersion theory, we study the 3-D velocity structures, the depths of interfaces and thicknesses of different layers, including the crust, the lithosphere and the asthenosphere, the lithosphere-asthenosphere system in the Qinghai-Tibet Plateau and its adjacent areas. The following tasks include: (1)The digital seismic records of 221 seismic events have been collected, whose magnitudes are larger than 5.0 over the Qinghai-Tibet Plateau and its adjacent areas. These records come from 31 digital seismic stations of GSN , CDSN、NCDSN and part of Indian stations. After making instrument response calibration and filtering, group velocities of fundamental mode of Rayleigh waves are measured using the frequency-time analysis (FTAN) to get the observed dispersions. Furthermore, we strike cluster average for those similar ray paths. Finally, 819 dispersion curves (8-150s) are ready for dispersion inversion. (2)From these dispersion curves, pure dispersion data in 2°×2° cells of the areas (18°N-42°N, 70°E-106°E) are calculated by using function expansion method, proposed by Yanovskaya. The average initial model has been constructed by taking account of global AK135 model along with geodetic, geological, geophysical, receiving function and wide-angle reflection data. Then, initial S-wave velocity structures of the crust and upper mantle in the research areas have been obtained by using linear inversion (SVD) method. (3)Taking the results of the linear inversion as the initial model, we simultaneously invert the S wave velocities and thicknesses by using non-linear inversion (improved Simulated Annealing algorithm). Moreover, during the temperature dropping the variable-scale models are used. Comparing with the linear results, the spheres and layers by the non-linear inversion can be recognized better from the velocity value and offset. (4)The Moho discontinuity and top interface of the asthenosphere are recognized from the velocity value and offset of the layers. The thicknesses of the crust, lithosphere and asthenosphere are gained. These thicknesses are helpful to studying the structural differentia between the Qinghai-Tibet Plateau and its adjacent areas and among geologic units of the plateau. The results of the inversion will provide deep geophysical evidences for studying deep dynamical mechanism and exploring metal mineral resource and oil and gas resources. The following conclusions are reached by the distributions of the S wave velocities and thicknesses of the crust, lithosphere and asthenosphere, combining with previous researches. (1)The crust is very thick in the Qinghai-Tibet Plateau, varying from 60 km to 80 km. The lithospheric thickness in the Qinghai-Tibet Plateau is thinner (130-160 km) than its adjacent areas. Its asthenosphere is relatively thicker, varies from 150 km to 230 km, and the thickest area lies in the western Qiangtang. India located in south of Main Boundary thrust has a thinner crust (32-38 km), a thicker lithosphere of about 190 km and a rather thin asthenosphere of only 60 km. Sichuan and Tarim basins have the crust thickness less than 50km. Their lithospheres are thicker than the Qinghai-Tibet Plateau, and their asthenospheres are thinner. (2)The S-wave velocity variation pattern in the lithosphere-asthenosphere system has band-belted distribution along east-westward. These variations correlate with geology structures sketched by sutures and major faults. These sutures include Main Boundary thrust (MBT), Yarlung-Zangbo River suture (YZS), Bangong Lake-Nujiang suture (BNS), Jinshajiang suture (JSJS), Kunlun edge suture (KL). In the velocity maps of the upper and middle crust, these sutures can be sketched. In velocity maps of 250-300 km depth, MBT, BNS and JSJS can be sketched. In maps of the crustal thickness, the lithospheric thickness and the asthenospheric thickness, these sutures can be still sketched. In particular, MBT can be obviously resolved in these velocity maps and thickness maps. (3)Since the collision between India and Eurasian plate, the “loss” of surface material arising from crustal shortening is caused not only by crustal thickening but also by lateral extrusion material. The source of lateral extrusion lies in the Qiangtang block. These materials extrude along the JSJS and BNS with both rotation and dispersion in Daguaiwan. Finally, it extends toward southeast direction. (4)There is the crust-mantle transition zone of no distinct velocity jump in the lithosphere beneath the Qiangtang Terrane. It has thinner lithosphere and developed thicker asthenosphere. It implies that the crust-mantle transition zone of partial melting is connected with the developed asthenosphere. The underplating of asthenosphere may thin the lithosphere. This buoyancy might be the main mechanism and deep dynamics of the uplift of the Qinghai-Tibet hinterland. At the same time, the transport of hot material with low velocity intrudes into the upper mantle and the lower crust along cracks and faults forming the crust-mantle transition zone.

大别山岩石圈有效弹性厚度初步研究

As powerful tools to study the lithosphere dynamics, the effective elastic thickness (Te) as well as the envelope of yielding stress of lithosphere have been attracted great attention of geoscientists in the past thirty years. The oceanic lithosphere, contrary to the continental lithosphere, has more fruits for its simple structures and evolution process. In continent, the lithosphere commonly is complex and variable in the rheological, thermal structures, and has a complicated history. Therefore, the application of the effective elastic thickness in continent is still a subject to learn in a long time. Te, with the definition of the thickness of an elastic plate in theory flexured by the equal benging of the real stress in the lithosphere plate (Turcotte, 1982), marks the depth of transition between elastic and fluid behaviors of rocks subjected to stress exceeding 100 MPa over the geological timescales (McNutt, 1990). There are three methods often adapted: admittance or isostatic response function, coherence and forwarding. In principle, the models of Te consist of thermal-rheological, non-linear Maxwell, non-linear work hardening and rheological layered models. There is a tentative knowledge of Te that it is affected by the following factors: crustal thickness, crust-mantle decoupling, plate bending, boundary conditions of plate (end forces and bending moments), stress state, sedimentary layer, faulting effect, variation in the mountain belts' strike, foreland basin, inheritance of tectonic evolution, convection of mantle, seismic depth and lithosphere strength. In this thesis, the author introduces the geological sketch of the Dabie collisional orogenic belt and the Hefei Basin. The Dabie Mts. is famous for the ultra-high pressure metamorphism. The crustal materials subducted down to the depth of at least 100 km and exhumed. So that the front subjects arise such as the deeply subduction of continent, and the post-collisional crust-mantle interaction. In a geological journey at June of 1999, the author found the rarely variolitic basaltic andesite in the Dabie Mts. It occurs in Susong Group, near Zhifenghe Countryside, Susong County, Anhui Province. It is just to the south of the boundary between the high-grade Susong melange and the ultra-high grade South Dabie melange. It has a noticeable knobby or pitted appearance in the surface. The size of the varioles is about 1-4 mm. In hand-specimen and under microscope, there are distinct contacts between the varioles and the matrice. The mineralogy of the varioles is primarily radiate plagioclase, with little pyroxene, hornblende and quartz. The pyroxene, hornblende and quartz are in the interstices between plagioclase. The matrix is consisted of glass, and micro-crystals of chlorite, epidote and zoisite. It is clearly subjected and extensive alteration. The andesite has an uncommon chemical composition. The SiO_2 content is about 56.8%, TiO_2 = 0.9%, MgO = 6.4%, (Fe_2O_3)_(Total) = 6.7% ～ 7.6%, 100 Mg/(Mg+Fe) = 64.1 ～ 66.2. Mg# is significantly high. The andesite has higher abundances of large-lithophile trace elements (e.g. K, Ba, Sr, LREE), e.g. La/Nd = 5.56-6.07, low abundances of high-strength-field elements (HFSE, e.g. Ta, Nb, P, Ti), particularly Ta and Nb strongly depleted. These are consistent with the characteristics of subducted-related magmas (Pearcce, 1982; Sun and McDonaugh, 1989). In the spider diagram of trace elements, from Ce to right hand, the abundances of elements decrease quickly, showing a characteristic of the continental margins (Pearce, 1982). There has a strongly enrichment of light-rare-earth elements, with a significant diffraction of REEs (the mean value of (La/Yb)_N is 32.84). No Eu anomaly, but there are anomaly high (La/Yb)_N = 28.63-36.74, (La/Y)_N = 70.33 - 82.84. The elements Y and Yb depleted greatly, Y < 20 ppm, Y_N = 2.74-2.84, Yb_N = 2.18 - 2.35. From the La-(La/Sm) diagram, the andesite is derived from partial melting. But the epsilone value of Nd is -18.7 ～ -19.2, so that the material source may be the mantle materials affected by the crustal materials. The Nd model age is 1.9 Ga indicating that the basaltic andesite was resulted from the post-collisional crust-mantle interaction between the subducted Yangze carton and the mantle of Sino-Korea carton. To obtain the Te of the lithosphere beneath the Dabie Mts. and the Hefei Basin, the author applies the coherence method in this thesis. The author makes two topography-gravity profiles (profiles 7~(th) and 9~(th)) across the Dabie Mts. and the Hefei Basin, and calculates the auto-coherence, across coherence, power spectrum, across power spectrum of the topography and gravity of the two profiles. From the relationships between the coherence and the wave-number of profiles. From the relationships between the coherence and the wave-number of profiles 7~(th) and 9~(th), it is obtained that the characteristic wavelengths respectively are 157 km and 126 km. Consequently the values of effective elastic thickness are 6.5 km and 4.8 km, respectively. However, the Te values merely are the minimum value of the lithosphere because the coherencemethod in a relative small region will generate a systemic underestimation. Why there is a so low Te value? In order to check the strength of the lithosphere beneath the Dabie Mts., the authore tries to outline the yielding-stress envelope of the lithosphere. It is suggested that the elastic layers in the crust and upper mantle are 18 km and 35 km, respectively. Since there exist a low viscosity layer about 3-5 km thickness, so it is reasonable that the decoupling between the crust and mantle occurred. So the effective thickness of the lithosphere can be estimated from the two elastic layers. Te is about 34 km. This is the maximum strength of the lithosphere. We can make an approximately estimation about the strength of the lithosphere beneath the Dabie Mts.: Te is about 20-30 km. The author believes that the following factors should be responsible for the low Te value: (1) the Dabie Mts. has elevated strongly since K_3-J_1. The north part of the Dabie Mts. elevates faster than the south part today; (2) there occur large active striking faults in this area. And in the east, the huge Tan-Lu striking fault anyway tends to decrease the lithosphere strength; (3) the lithosphere beneath the Dabie Mts. is heter-homogeneous in spatio-temporal; (4) the study area just locates in the adjacent region between the eastern China where the lithosphere thickness is significantly reduced and the normal western China. These factors will decrease the lithosphere strength.