Triggered earthquakes as stress gauge: Implication for the uppercrust behavior in the Grenoble area,France

A sequence of moderate shallow earthquakes (3.5M _L5.3) was located within the Vercors massif (France) in the period 1961–1984. This subalpine massif has been a low seismic area for at least 5 centuries. During the period 1962–1963, 12 shallow earthquakes occurred in the neighborhood (10 km) of the Monteynard reservoir, 30 km south of the city of Grenoble. The latest fourM _L4.0 earthquakes occurred in 1979–1984 either at larger distance (35 km) or greater depth (10 km) from the reservoir. Two triggering mechanisms are suggested for this sequence: (i) the direct effect of elastic loading through either increased shear stress or strength reducing by increased pore pressure at depth; (ii) the pore pressure diffusion induced by poroelastic stress change due to the reservoir filling.The weekly water levels, local balanced geological cross sections, and focal mechanisms argue for two types of mechanical connection between the earthquake sequence and the filling cycles of the Monteynard reservoir. The seismic sequence started with the 1962–1963 shallow earthquakes that occurred during the first filling of the reservoir and are typical of the direct effect of elastic loading. The 1979 deeper earthquake is located at a 10 km depth below the reservoir. This event occurred 16 years after the initial reservoir impoundment, but one month after the previous 1963 maximum water level was exceeded. Moreover the yearly reservoir level increased gradually in the period 1962–1979 and has decreased since 1980. Accordingly we suggest that the gradual diffusion of water from reservoir to hypocentral depths decreases the strength of the rock matrices through increased pore pressure. The transition between the two types of seismic response is supported by the analysis ofM _L3.5 earthquakes which all occurred in the period 1964–1971, ranging between 10 and 30 km distance from the reservoir. The three other delayed earthquakes of the 1961–1984 seismic sequence (M _L4 during the 1979–1984 period) are all located 35 km away from the reservoir. Based on the seismic activity, the estimates for the hydraulic diffusivities range between 0.2–10 m²/s, except for the first event that occurred 30 km north of the reservoir, the filling just started. The lack ofin situ measurements of crustal hydrological properties in the area, shared by most of the Reservoir-Induced-Seismicity cases, prevents us from obtaining absolute evidence for the triggering processes. These observations and conceptual models attest that previous recurrence times for moderate natural shocks (4.5M _L5.5) estimated within this area using historical data, could be modified by 0.1–1 MPa stress changes. These small changes in deviatoric stress suggest that the upper crust is in this area nearly everywhere at a state of stress near failure. Although the paucity of both number and size of earthquakes in the French subalpine massif shows that aseismic displacements prevail, our study demonstrates that triggered earthquakes are important tools for assessing local seismic risk through mapping fault zones and identifying their possible seismic behavior.     

利用原地应力实测资料分析芦山地震震前浅部地壳应力积累

龙门山断裂带南段的芦山MS7.0级地震给雅安地区人民群众和当地社会带来巨大的损失,如何从地壳应力积累的角度来研究该次地震对于深入认识孕震发震机理将有很大的帮助.项目组在过去21年中在该区域积累了22个钻孔的应力实测资料,这些应力资料显示该区域的应力状态为逆冲断层应力状态,最大水平主压应力方向为N44°—64°W,而且在过去20多年里,该区域的应力状态是稳定的.文中提出利用平均差应力和平均有效应力之比μm来表征地壳应力积累的能力和水平.通过理论分析可知μm的物理意义与μ是近似的,完全可以用μm代替μ来开展分析.利用实测应力数据计算得到的μm范围为0.39~0.56,接近于Byerlee定律所定义的μ=0.6的下限值,该范围与龙门山断裂带震后断层泥稳态摩擦系数的范围完全吻合,芦山地震的发生说明研究区的地壳应处于摩擦极限平衡或者亚平衡状态.芦山地震的震源机制解反应的应力状态为逆冲断层应力状态,主压应力轴方向为122°,与该区域的应力测试数据所反应的应力状态一致.通过与他人的研究成果对比分析可知,芦山地震可看作汶川地震对龙门山断裂带南段区域的断层加载而诱发产生的一次地震.     

Active tectonics of the outer northern Apennines: Adriatic vs. Po Plain seismicity and stress fields

The recent earthquake sequences of 2012 (northern Italy) and 2013 (Marche offshore) provided new, fundamental constraints to the active tectonic setting of the outer northern Apennines. In contrast to the Po Plain area, where the 2012 northern Italy earthquakes confirmed active frontal thrusting, the new focal mechanisms obtained in this study for the 2013 Marche offshore earthquakes indicate that only minor thrust fault reactivation occurs in the Adriatic domain, even for a theoretically favourably oriented maximum horizontal compression. Recent seismicity in this domain appears to be mainly controlled by transcurrent crustal faults dissecting the Apennine thrust belt. The along-strike stress field variation from the Po Plain to the Adriatic area has been quantitatively investigated by applying the multiple inverse method (MIM) to the analysis of the entire seismicity recorded from January 1976 to August 2014, from the top 12 km of the crust (fault plane solutions from 127 earthquakes with M_W ≥ 4), allowing us to obtain a comprehensive picture of the state of stress over the outer zone of the fold and thrust belt. The present-day stress field has been defined for 39 cells of 1.5° × 1.5° surface area and 12 km depth. The obtained stress field maps point out that, although the entire outer northern Apennines belt is characterized by a sub-horizontal maximum compressive axis (σ₁), the minimum compression (σ₃) is sub-vertical only in the Po Plain area, becoming sub-horizontal in the Adriatic sector, thus confirming that the latter region is dominated by an active tectonic regime of strike-slip type.     

Comparison of the results obtained from different techniques such as overcoring, hydrofracturing, breakouts and focal mechanism solutions

In China, a lot of data derived from focal mechanism solutions, and especially, from overcoring and hydrofracturing stress measurements, as well as borehole breakouts carried out for various aims have contributed to understanding of the state of the crustal stress. We conducted everin situ stress measurements in a borehole with different methods, for the sake of comparison of various measuring techniques. Furthermore, we carried out measurements in different borehole with one or two kinds of techniques. All the observations taken out in these ways showed little variations. So we believe that the data obtained by the techniques can be used to study the tectonic stress field of certain region and thein situ state of the crustal stress. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,14, 149–155, 1992.     

Investigation of coupling between surface processes and induced flow in the lower continental crust as a cause of intraplate seismicity

Many studies have highlighted the role of coupling between surface processes and flow in the lower continental crust in deforming the crust and creating topographic relief over Quaternary timescales. On the basis of the rheological knowledge gained, it is suggested that intraplate seismicity can also be caused by coupling between surface processes and flow in the lower continental crust. This view is shown to be a natural consequence of the modern idea that isostatic equilibrium is maintained by flow in the weak lower crust in response to erosion and sedimentation. It is supported by a general correlation between the vigour of surface processes and rates of intraplate seismicity, and by instances of seasonal seismicity that correlates with seasonal climate. Human interference in the environment can affect surface loading: for instance, deforestation for agriculture or urban development can cause increased erosion rates; global warming is expected to cause increased storminess (and thus increased erosion rates) and/or global sea‐level rise. The possibility of increased rates of seismicity resulting from these processes should thus be considered in future hazard assessment. Copyright © 2006 John Wiley & Sons, Ltd.     

Recycling of the continental crust

In order to understand the evolution of the crust-mantle system, it is important to recognize the role played by the recycling of continental crust. Crustal recycling can be considered as two fundamentally distinct processes: 1) intracrustal recycling and 2) crust-mantle recycling. Intracrustal recycling is the turnover of crustal material by processes taking place wholly within the crust and includes most sedimentary recycling, isotopic resetting (metamorphism), intracrustal melting and assimilation. Crust-mantle recycling is the transfer of crustal material to the mantle with possible subsequent return to the crust. Intracrustal recycling is important in interpreting secular changes in sediment composition through time. It also explains differences found in crustal area-age patterns measured by different isotopic systems and may also play a role in modeling crustal growth curves based on Nd-model ages. Crustal-mantle recycling, for the most part, is a subduction process and may be considered on three levels. The first is recycling with only short periods of time in the mantle (<10 m.y.). This may be important in explaining the origin of island-arc and related igneous rocks; there is growing agreement that 1–3% recycled sediment is involved in their origin. Components of recycled crustal material, with long-term storage (up to 2.5 b.y.) in the mantle as distinct entities, has been suggested for the origin of ocean island and ultrapotassic volcanics but there is considerably less agreement on this interpretation. A third proposal calls for the return of crustal material to the mantle with efficient remixing in order to swamp the geochemical and isotopic signature of the recycled component by the mantle. This type of recycling is required for steady-state models of crustal evolution where the mass of the continents remains constant over geological time. It is unlikely if crust-mantle recycling has exceeded 0.75 km³/yr over the past 1–2 Ga.Good evidence exists that selective recycling is an important process. Sedimentary rocks preserved in different tectonic settings are apparently recycled at different rates, resulting in a bias in the sediment types preserved in the geologic record. Selective recycling has important implications for the interpretation of Nd model ages of old sedimentary rocks and in the analysis of accreted terranes. Although there is evidence that continental crust was formed prior to 3.8 Ga, the oldest preserved rocks do not exceed this age. It is likely that the intense meteorite bombardment, which affected the earth during the period 4.56–3.8 Ga, coupled with rapid mantle convection, which resulted from greater heat production, caused the destruction and probable recycling into the mantle of any early formed crust.Although crust-mantle recycling is seen as a viable process, it is concluded that crustal growth has exceeded crust-mantle recycling since at least 3.8 Ga. Intracrustal recycling has not been given adequate consideration in models of crustal growth based on isotopic data (particularly Nd model ages). It is concluded that crustal growth curves based on Nd model ages, while vastly superior to those based on K/Ar or Rb/Sr, tend to underestimate the volume of old crust, due to crust-mantle and/or intracrustal recycling.     

Source parameters of mining-induced seismic events: An evaluation of homogeneous and inhomogeneous faulting models for assessing damage potential

Source parameter estimates based on the homogeneous and inhomogeneous source models have been examined for an anomalous sequence of seven mine-induced events located between 640 and 825 m depth at Strathcona mine, Ontario, and having magnitudes ranging betweenm _N 0.8 and 2.7. The derived Brune static stress drops were found to be similar to those observed for natural earthquakes (30 bars), whereas dynamic stress drops were found to range up to 250–300 bars. Source radii derived from Madariaga's model better fit documented evidence of underground damage. These values of source radii were similar to those observed for the inhomogeneous model. The displacement at the source, based on the observed attenuation relationship, was about 60 mm for three magnitude 2.7 events. This is in agreement with slip values calculated using peak velocities and assuming the asperity as a Brune source within itself (72 mm). By using Madariaga's model for the asperity, the slip was over 3 times larger than observed. Peak velocity and acceleration scaling relations with magnitude were investigated by incorporating available South African data, appropriately reduced to Canadian geophysical conditions. The dynamic stress drop scaled as the square root of the seismic moment, similar to reported results in the literature for crustal earthquakes. This behavior suggests that the size of the asperities responsible for the peak ground motion, with respect to the overall source size, follow distributions that may be similar over a wide range of magnitudes. Measurements of source rupture complexity (ranging from 2 to 4) were found to agree with estimates of overall source to asperity radii, suggesting, together with the observed low rupture velocities (0.3 to 0.6 ), that the sources were somewhat complex. Validation of source model appropriateness was achieved by direct comparison of the predicted ground motion level to observed underground damage in Creighton mine, located within the same regional stress and geological regime as Strathcona mine. Close to the source (<100 m), corresponding to relatively higher damage levels, a good agreement was found between the predicted peak particle velocities for the inhomogeneous model and velocities derived based on established geomechanical relationships. The similarity between asperity radii and the regions of the highest observed damage provided additional support for the use of the inhomogeneous source model in the assessment of damage potential.     

Zircon xenocrysts from the Kambalda volcanics: age constraints and direct evidence for older continental crust below the Kambalda-Norseman greenstones

The Hangingwall Basalt at Kambalda, Western Australia, contains zircons that have been shown by ion microprobe analyses to have very high U and Th contents and a wide variety of crystallization ages. Nearly all of these zircons certainly are xenocrysts; a few might relate to intrusive veinlets. The age of the youngest xenocrysts, 2693 ± 50Ma(2 σ), shows that the eruptive age of the basalt cannot exceed 2743 Ma. This confirms that the apparent SmNd isochron giving 3200 Ma [1,2] for Kambalda mafic and ultramafic rocks is a mixing-line [2] between unrelated components enriched and depleted in light rare earth elements. Mixing probably occurred at depth by erosion of 3200–3500 Ma old felsic crust from the walls of the HWB conduits. The zircon xenocryst ages are the first direct evidence for the presence of very old felsic crust in the eastern Yilgarn Block. The latter implies that the Kalgoorlie-Norseman greenstone sequences were formed in a continental rather than an oceanic environment.     

Seismogenic potential and earthquake hazard assessment in the Tell Atlas of Algeria

Seismotectonic zonation studies in the Tell Atlas of Algeria, a branch of the Africa-Eurasia plate boundary, provide a valuable input for deterministic seismic hazard calculations. We delineate a number of seismogenic zones from causal relationships established between geological structures and earthquakes and compile a working seismic catalogue mainly from readily available sources. To this catalogue, for a most rational and best-justified hazard analysis, we add estimates of earthquake size translated from active faulting characteristics. We assess the regional seismic hazard using a deterministic procedure based on the computation of complete synthetic seismograms (up to 1 Hz) by the modal summation technique. As a result, we generate seismic hazard maps of maximum velocity, maximum displacement, and design ground acceleration that blend information from geology, historical seismicity and observational seismology, leading to better estimates of the earthquake hazard throughout northern Algeria. Our analysis and the resulting maps illustrate how different the estimate of seismic hazard is based primarily on combined geologic and seismological data with respect to the one for which only information from earthquake catalogues has been used.     

Interaction between subducted continental crust and the mantle

The syncollisional mafic-ultramafic rocks with Nb, Zr, Ti negative anomalies in the North Dabie terrane have Sr and Nd isotopic compositions with EMI features. Their low and variable initial e_Nd values ranging from ?2 to ?18 are similar to those of their gneissic country rocks and the ultrahigh-pressure metamorphic rocks in the South Dabie terrane. These Sr and Nd isotopic features are difficult to be interpreted by mantle metasomatism related to oceanic subduction or crust assimilation, but is best explained by the mantle metasomatism related to continental subduction.     

Inversion of inline and broadside marine controlled‐source electromagnetic data with constraints derived from seismic data

We present a structural smoothing regularization scheme in the context of inversion of marine controlled‐source electromagnetic data. The regularizing hypothesis is that the electrical parameters have a structure similar to that of the elastic parameters observed from seismic data. The regularization is split into three steps. First, we ensure that our inversion grid conforms with the geometry derived from seismic. Second, we use a seismic stratigraphic attribute to define a spatially varying regularization strength. Third, we use an indexing strategy on the inversion grid to define smoothing along the seismic geometry. Enforcing such regularization in the inversion will encourage an inversion result that is more intuitive for the interpreter to deal with. However, the interpreter should also be aware of the bias introduced by using seismic data for regularization. We illustrate the method using one synthetic example and one field data example. The results show how the regularization works and that it clearly enforces the structure derived from seismic data. From the field data example we find that the inversion result improves when the structural smoothing regularization is employed. Including the broadside data improves the inversion results even more, due to a better balancing between the sensitivities for the horizontal and vertical resistivities.     

Interaction between subducted continental crust and the mantle

The syncollisional mafic-ultramafic intrusions in the North Dabie terrane are characterized by enriched LREE, Rb, Ba, and depletion of high field strength elements, such as Nb, Zr, Ti. The high Ti/& and Ti/Y ratios of their parent magma suggest that crust contamination of magma is not important to the above trace elements characters. They mainly reflect the features of their mantle source which has been modified by metasomatism in subduction environment.     

基于方位地震数据的地应力反演方法

在页岩油气藏的开发和勘探阶段,需要对储层进行水力压裂改造,形成有利于油气聚集和运移的裂缝.地应力是进行水力压裂改造的重要参数,能够决定裂缝的大小、方向以及分布形态,影响着压裂的增产效果,且最大和最小水平应力差异比（ODHSR,Orthorhombic Differential Horizontal Stress Ratio）是评价储层是否可压裂成网的重要因子.本文探讨了基于地震数据估算地应力的方法,以指导页岩气的水力压裂开发.首先,利用叠前方位地震数据反演得到地层的弹性参数和各向异性参数;其次,基于正交各向异性水平应力差异比近似公式,利用反演得到的弹性参数和各向异性参数估算地层的ODHSR;最后,选取某工区的裂缝型页岩储层的叠前方位地震数据对该方法进行实际应用.实际工区地震数据应用表明,基于叠前方位地震数据反演得到的ODHSR能够有效的识别储层中易于压裂成网的区域.

     

Gabbro and related rock emplacement beneath rifting continental crust: UPb geochronological and geochemical constraints for the Galicia passive margin (Spain)

The thinned continental crust of the west Galicia margin is bound by a belt of serpentinized peridotites (‘peridotite ridge’) lying about 300 km off the coast in the North Atlantic ocean. From this ridge, a gabbro and a chlorite rock were studied in an attempt to substantiate rift-related subcontinental magmatism, occurring prior to sea-floor spreading. U-Pb dating of 13 different zircon fractions yields a precise age of 122.1 ± 0.3 Ma (2σ) for the emplacement of the chlorite rock protolith, from which more than 50% of Si and alkali-calc-alkali elements were lost during greenschist facies tectonometamorphism. Sr and Nd isotope signatures suggest that the gabbro and chlorite rock protoliths were derived from mantle sources that were moderately depleted in LILE, relative to a chondritic reservoir. No evidence for the presence of continental material in the magma source regions can be observed. From the new zircon age of 122.1 ± 0.3 Ma, and earlier determined³⁹Ar⁴⁰Ar age of 122.0 ± 0.6 Ma for amphibole from the same locality, it can be documented that magma formation, solidification and unroofing of the mantle rocks occurred during a short period of time of about 3.4 Ma, which means that the peridotite ridge detached from the continent and rose to the surface immediately after, or even coevally with mantle melting.     

Thermal models for the magmatic accretion and subsequent metamorphism of continental crust

I present one-dimensional conductive relaxation models for the thermal evolution of continental crust which has undergone magmatic thickening in a predominantly recumbent tectonic regime. The sustained addition of magma, and associated heat, profoundly influences the conductive relaxation of the affected sialic crust. The mechanism of accretion, successively beneath (under-accretion), or successively above (over-accretion), previously added material, strongly affects the type of P-T-time history sustained by the rocks in the middle and deep crust. The pressure and temperature conditions which might be recorded by mineral assemblages equilibrating in such tectonothermal regimes may bear very little relation to conductive heat supplies from the upper mantle. P-T values obtained by geobarometry and geothermometry in such complexes must be used with caution to constrain steady state heat flows and geothermal gradients. In a contemporary plate tectonic context I consider the tectonothermal models developed here to apply particularly to island arcs and sites of continental collision associated with extensive plutonic igneous activity. Precambrian orthogneiss complexes, such as the Archaean of southern West Greenland, probably represent exhumed deep crustal elements involved in the intense magmatic thickening processes considered here.     

关于应力触发和应力影概念在地震预报中应用的一些思考

应力触发和应力影概念及其在地震预报中的应用,是学术界关注的热点,尚有一系列需要研究的问题。例如初始应力场的确定,合理流变模型的选取和计算,孔隙压力的变化和影响,断层强度随时间的变化等。特别在估计应力场变化时,人们现今只注意到大地震对应力场造成的变化。一些新的观测事实表明,大量的小地震效应,慢地震的效应和无震滑动的效应,也不容忽视。定量考虑它们的效应,是今后应力触发和应力影研究中应该注意的课题。     

Mechanics of seismic instabilities induced by the recovery of hydrocarbons

We review earthquake distributions associated with hydrocarbon fields in the context of pore pressure diffusion models, poroelastic stress transfer and isostasy theory. These three mechanisms trigger or induce seismic instabilities at both local scale (D5 km) and at regional scale (D20 km). The modeled changes in stress are small (1 MPa), whatever the tectonic setting. Each mechanism corresponds to different production processes. (1) Local hydraulic fracturing due to fluid injection induces seismic-slip on cracks (M _L3) within the injected reservoir through decreasing the effective stress. (2) Pure fluid withdrawal causes pore pressure to decrease within the reservoir. It triggers adjustments of the geological structure to perturbations related to the reservoir response to depletion. Poroelastic mechanisms transfer this stress change from the reservoir to the surrounding levels whereM _L5 seismic instabilities occur either above or below the reservoir. (3) Massive hydrocarbon recovery induces crustal readjustments due to the removal of load from the upper crust. It can induce larger earthquakes (M _L6) at greater distance from the hydrocarbon fields than the two other mechanisms.Due to the mechanical properties of the shallow rock matrices involved, seismic slip triggered either by mechanism (1) or (2), is a second-order process of the main elastoplastic deformation. for a minimum of 80% of commercially productive basins, most of the local deformation is reported as aseismic, i.e., there is no evidence forM _L3 earthquakes. Nevertheless, the induced stresses vary as a function of time in a manner that depends on the hydraulic diffusivity (i.e., permeability) of the reservoir and surrounding rocks. Because small earthquakes (M _L3) indicate changes in stress and pore pressure, monitoring of seismicity is a means of assessingin situ reservoir behavior.The less constrained seismic response to hydrocarbon recovery is the possible connection between local fluid manipulations, triggered earthquakes and major regional earthquakes. Positive feedback mechanisms suggest that the region of seismic hazard changes is much larger than the area where hydrocarbons are extracted. These observations and models testify that fluid movement and pore pressure changes (increase or decrease) play important roles in the mechanics of earthquakes and in the triggering of natural earthquakes.     

Regional stress field in the area of Stromboli (Italy): insights into structural data and crustal tectonic earthquakes

Crustal tectonic seismicity on Stromboli is here discussed along with the regional tectonic and geodynamic context of the Aeolian archipelago, Italy. The aim of this paper is to contribute to the reduction of the crucial gap in information concerning the tectonic activity originating at the volcano Stromboli in the last decade. It is indeed well known that Strombolian eruptions, related to the persistent state of activity on this volcano, allow a fast collection of seismic data which are mainly represented by explosion quakes and tremor. However, only some sporadic information concerns tectonic earthquakes due to their low occurrence rate. A catalog of 150 crustal tectonic earthquakes recorded at Stromboli from 1985 to 1996 is presented here. The analyses of the seismograms were carried out on the basis of the records of the station STR, located in the southwestern part of the island. The magnitude considered ranged from 1.5 to 3.7. Hypocentral location indicates foci at less than 12 km, although for the majority of the earthquakes, which are not recorded at the seismic stations deployed on the other islands of the Aeolian archipelago, a shallower (within a few kilometers) depth can be hypothesized. The strain release associated with the 150 tectonic earthquakes of our catalog was also evaluated. The low levels of occurrence and seismic energy release reached throughout the ca. twelve years analyzed here confirm the minor role played by processes of brittle failure of rocks. This aspect of seismic activity is nevertheless important as it can provide a contribution to the knowledge of the dynamics of this volcano. The focal mechanisms available, integrated with data from the literature, show the T-axes oriented about NE–SW and NW–SE. These findings are consistent with the deformative and tectonic characteristics of this volcanic edifice and of the Aeolian archipelago, which shows the same two directions of extension. This is in agreement with subduction, which is very old or which passively evolves towards the final stages, and is controlled by the transform activity of the Aeolian–Tindari–Giardini fault system.     

Structure and composition of the continental crust in East China

Crustal structures of nine broad tectonic units in China, except the Tarim craton, are derived from 18 seismic refraction profiles including 12 geoscience transects. Abundances of 63 major, trace and rare earth elements in the upper crust in East China are estimated. The estimates are based on sampling of 11 451 individual rock samples over an area of 950 000 km², from which 905 large composite samples are prepared and analyzed by 13 methods. The middle, lower and total crust compositions of East China are also estimated from studies of exposed crustal cross sections and granulite xenoliths and by correlation of seismic data with lithologies. All the tectonic units except the Tarim craton and the Qinling orogen show a four-layered crustal structure, consisting of the upper, middle, upper lower, and lowermost crusts. P-wave velocities of the bulk lower crust and total crust are 6.8–7.0 and 6_:4–6.5 km/s, respectively. They are slower by 0.2–0.4 km/s than the global averages. The bulk lower crust is suggested to be intermediate with 58% SiO₂ in East China. The results contrast with generally accepted global models of mafic lower crusi. The proposed total crust composition in East China is also more evolved than previous estimates and characterized by SiO₂=64%, a significant negative Eu anomaly (Eu/Eu* = 0.80), deficits in Sr and transition metals, a near-arc magma La/Nd ratio (3.0), and a calculatedμ(²³⁸U/²⁰⁴Pb) value of 5. In addition, it has the following ratios of element pairs exhibiting similar compatibility, which are identical or close to the primitive mantle values: Zr/Hf=37, Nb/Ta=17.5, Ba/Th=87, K/Pb=0.12x10⁴, Rb/Cs=25, Ba/Rb=8.94, Sn/Sm=0.31, Se/Cd=1.64, La/ As=10.3, Ce/Sb=271, Pb/Bi=57, Rb/TI=177, Er/Ag=52, Cu/Au=3.2×10⁴, Sm/Mo=7.5, Nd/W=40, CI/Li=10.8, F/Nd=21.9, and La/B=1.8. Project supported by the National Natural Science Foundation of China (Grant Nos. 49625305, 49573183, 49673184, 49794043), the State Comission of Education, the Ministry of Geology and Mineral Resources of China (Grant No. 850514), the Open Laboratory of Constitution, Interaction and Dynamics of the Crust-Mantle System, and the Alexander-von-Humboldt Foundation of Germany.     

Helium, volatile fluxes and the development of continental crust

Mantle-derived helium has a substantial primordial component and is readily distinguished from radiogenic “crustal” He by its isotopic composition. For some years it has been known to be escaping at mid-ocean ridges and more recently it has been shown to be escaping through the continental lithosphere in tectonically active areas, particularly those undergoing extension or volcanism. The C/³He value observed in ocean ridge basalts and continental gases that contain only mantle He, is close to 10⁹. This is believed to be a typical value for the upper mantle. Other continental gases have ratios that vary widely and are diluted with crustal carbon. The ratio C/⁴He decreases with time through the production of radiogenic⁴He, and depends on the C/(U + Th) value. Departures from the average may result from exceptional concentrations of U and Th or from C/He fractionation.There is circumstantial evidence for a steady-state flux of He through the continents that may be estimated from He accumulations in lakes and aquifers. The mantle component of such fluxes is calculated from their³He content. If the mantle component is accompanied by C in the proportion indicated above, and extensional areas make up as little as 10% of the crust at any one time, then about 10% of the present inventory of crustal C would have been added to the crust every Ga by this means. C/K values for the crust and mantle are today very similar, and K may therefore scale as C. K/U and K/Th vary within narrow limits and they may scale with C also.The most plausible means of scavenging He from the mantle is by partial melting: He is expected to enter the first few percent of liquid formed, and the loss of mantle He and C at the surface is associated with the emplacement of basaltic bodies in the lower crust carrying K, U and Th. Some limits are placed on the thickness of basalt added in extensional areas.Mantle-derived CO₂ has often been invoked as a means of dehydrating continental crust to produce granulites. However, the amounts of CO₂, estimated from mantle He fluxes, entering the crust in those active tectonic areas studied so far appears too small to produce dehydration on a regional scale. The addition of mantle-derived material to the crust in extensional zones is a first-order crustal growth process the importance of which has previously been underestimated.