Heat-flow studies in the Peru Trench subduction zone

New heat-flow values were obtained in the central Peru Trench area during site surveys and drilling of Ocean Drilling Program (ODP) Leg 112 by measuring temperatures with ordinary surface heat-flow probes and in the drill holes and by estimating from bottom-simulating reflectors resulting from gas hydrates. The values determined by these methods are consistent with each other within the limits of error. When combined with existing data, heat-flow distribution from the trench to the coast was delineated. Heat flow is lower than 40 mW/m**2 at the bottom of the trench and 40 to 50 mW/m**2 on the landward slope. The low heat flow at the trench bottom can be explained partly by a high sedimentation rate. Heat flow is variable about where the Mendana Fracture Zone meets the trench. This low heat flow might result from hydrothermal circulation in the fracture zone, which some scientists believe is a new propagating rift. On the landward slope, no significant difference in heat flow is recognized between the northern side and the southern side of the fracture zone, in spite of differences in the age of the subducting plate and the tectonic history. Heat flow on the landward slope may be slightly higher than that in most other subduction zones.

Paleoseismological evidence of surface faulting along the northeastern Himalayan front, India: Timing, size, and spatial extent of great earthquakes

The similar to 2500 km long Himalayan arc has experienced three large to great earthquakes of M-w 7.8 to 8.4 during the past century, but none produced surface rupture. Paleoseismic studies have been conducted during the last decade to begin understanding the timing, size, rupture extent, return period, and mechanics of the faulting associated with the occurrence of large surface rupturing earthquakes along the similar to 2500 km long Himalayan Frontal Thrust (HFT) system of India and Nepal. The previous studies have been limited to about nine sites along the western two-thirds of the HFT extending through northwest India and along the southern border of Nepal. We present here the results of paleoseismic investigations at three additional sites further to the northeast along the HFT within the Indian states of West Bengal and Assam. The three sites reside between the meizoseismal areas of the 1934 Bihar-Nepal and 1950 Assam earthquakes. The two westernmost of the sites, near the village of Chalsa and near the Nameri Tiger Preserve, show that offsets during the last surface rupture event were at minimum of about 14 m and 12 m, respectively. Limits on the ages of surface rupture at Chalsa (site A) and Nameri (site B), though broad, allow the possibility that the two sites record the same great historical rupture reported in Nepal around A.D. 1100. The correlation between the two sites is supported by the observation that the large displacements as recorded at Chalsa and Nameri would most likely be associated with rupture lengths of hundreds of kilometers or more and are on the same order as reported for a surface rupture earthquake reported in Nepal around A.D. 1100. Assuming the offsets observed at Chalsa and Nameri occurred synchronously with reported offsets in Nepal, the rupture length of the event would approach 700 to 800 km. The easternmost site is located within Harmutty Tea Estate (site C) at the edges of the 1950 Assam earthquake meizoseismal area. Here the most recent event offset is relatively much smaller (<2.5 m), and radiocarbon dating shows it to have occurred after A.D. 1100 (after about A.D. 1270). The location of the site near the edge of the meizoseismal region of the 1950 Assam earthquake and the relatively lesser offset allows speculation that the displacement records the 1950 M-w 8.4 Assam earthquake. Scatter in radiocarbon ages on detrital charcoal has not resulted in a firm bracket on the timing of events observed in the trenches. Nonetheless, the observations collected here, when taken together, suggest that the largest of thrust earthquakes along the Himalayan arc have rupture lengths and displacements of similar scale to the largest that have occurred historically along the world's subduction zones.

Determination of source finiteness and depth of large earthquakes

In this thesis, a method to retrieve the source finiteness, depth of faulting, and the mechanisms of large earthquakes from long-period surface waves is developed and applied to several recent large events.

In Chapter 1, source finiteness parameters of eleven large earthquakes were determined from long-period Rayleigh waves recorded at IDA and GDSN stations. The basic data set is the seismic spectra of periods from 150 to 300 sec. Two simple models of source finiteness are studied. The first model is a point source with finite duration. In the determination of the duration or source-process times, we used Furumoto's phase method and a linear inversion method, in which we simultaneously inverted the spectra and determined the source-process time that minimizes the error in the inversion. These two methods yielded consistent results. The second model is the finite fault model. Source finiteness of large shallow earthquakes with rupture on a fault plane with a large aspect ratio was modeled with the source-finiteness function introduced by Ben-Menahem. The spectra were inverted to find the extent and direction of the rupture of the earthquake that minimize the error in the inversion. This method is applied to the 1977 Sumbawa, Indonesia, 1979 Colombia-Ecuador, 1983 Akita-Oki, Japan, 1985 Valparaiso, Chile, and 1985 Michoacan, Mexico earthquakes. The method yielded results consistent with the rupture extent inferred from the aftershock area of these earthquakes.

In Chapter 2, the depths and source mechanisms of nine large shallow earthquakes were determined. We inverted the data set of complex source spectra for a moment tensor (linear) or a double couple (nonlinear). By solving a least-squares problem, we obtained the centroid depth or the extent of the distributed source for each earthquake. The depths and source mechanisms of large shallow earthquakes determined from long-period Rayleigh waves depend on the models of source finiteness, wave propagation, and the excitation. We tested various models of the source finiteness, Q, the group velocity, and the excitation in the determination of earthquake depths.

The depth estimates obtained using the Q model of Dziewonski and Steim (1982) and the excitation functions computed for the average ocean model of Regan and Anderson (1984) are considered most reasonable. Dziewonski and Steim's Q model represents a good global average of Q determined over a period range of the Rayleigh waves used in this study. Since most of the earthquakes studied here occurred in subduction zones Regan and Anderson's average ocean model is considered most appropriate.

Our depth estimates are in general consistent with the Harvard CMT solutions. The centroid depths and their 90 % confidence intervals (numbers in the parentheses) determined by the Student's t test are: Colombia-Ecuador earthquake (12 December 1979), d = 11 km, (9, 24) km; Santa Cruz Is. earthquake (17 July 1980), d = 36 km, (18, 46) km; Samoa earthquake (1 September 1981), d = 15 km, (9, 26) km; Playa Azul, Mexico earthquake (25 October 1981), d = 41 km, (28, 49) km; El Salvador earthquake (19 June 1982), d = 49 km, (41, 55) km; New Ireland earthquake (18 March 1983), d = 75 km, (72, 79) km; Chagos Bank earthquake (30 November 1983), d = 31 km, (16, 41) km; Valparaiso, Chile earthquake (3 March 1985), d = 44 km, (15, 54) km; Michoacan, Mexico earthquake (19 September 1985), d = 24 km, (12, 34) km.

In Chapter 3, the vertical extent of faulting of the 1983 Akita-Oki, and 1977 Sumbawa, Indonesia earthquakes are determined from fundamental and overtone Rayleigh waves. Using fundamental Rayleigh waves, the depths are determined from the moment tensor inversion and fault inversion. The observed overtone Rayleigh waves are compared to the synthetic overtone seismograms to estimate the depth of faulting of these earthquakes. The depths obtained from overtone Rayleigh waves are consistent with the depths determined from fundamental Rayleigh waves for the two earthquakes. Appendix B gives the observed seismograms of fundamental and overtone Rayleigh waves for eleven large earthquakes.

Earthquakes of the Nepal Himalaya : towards a physical model of the seismic cycle

Home to hundreds of millions of souls and land of excessiveness, the Himalaya is also the locus of a unique seismicity whose scope and peculiarities still remain to this day somewhat mysterious. Having claimed the lives of kings, or turned ancient timeworn cities into heaps of rubbles and ruins, earthquakes eerily inhabit Nepalese folk tales with the fatalistic message that nothing lasts forever. From a scientific point of view as much as from a human perspective, solving the mysteries of Himalayan seismicity thus represents a challenge of prime importance. Documenting geodetic strain across the Nepal Himalaya with various GPS and leveling data, we show that unlike other subduction zones that exhibit a heterogeneous and patchy coupling pattern along strike, the last hundred kilometers of the Main Himalayan Thrust fault, or MHT, appear to be uniformly locked, devoid of any of the “creeping barriers” that traditionally ward off the propagation of large events. The approximately 20 mm/yr of reckoned convergence across the Himalaya matching previously established estimates of the secular deformation at the front of the arc, the slip accumulated at depth has to somehow elastically propagate all the way to the surface at some point. And yet, neither large events from the past nor currently recorded microseismicity nearly compensate for the massive moment deficit that quietly builds up under the giant mountains. Along with this large unbalanced moment deficit, the uncommonly homogeneous coupling pattern on the MHT raises the question of whether or not the locked portion of the MHT can rupture all at once in a giant earthquake. Univocally answering this question appears contingent on the still elusive estimate of the magnitude of the largest possible earthquake in the Himalaya, and requires tight constraints on local fault properties. What makes the Himalaya enigmatic also makes it the potential source of an incredible wealth of information, and we exploit some of the oddities of Himalayan seismicity in an effort to improve the understanding of earthquake physics and cipher out the properties of the MHT. Thanks to the Himalaya, the Indo-Gangetic plain is deluged each year under a tremendous amount of water during the annual summer monsoon that collects and bears down on the Indian plate enough to pull it away from the Eurasian plate slightly, temporarily relieving a small portion of the stress mounting on the MHT. As the rainwater evaporates in the dry winter season, the plate rebounds and tension is increased back on the fault. Interestingly, the mild waggle of stress induced by the monsoon rains is about the same size as that from solid-Earth tides which gently tug at the planets solid layers, but whereas changes in earthquake frequency correspond with the annually occurring monsoon, there is no such correlation with Earth tides, which oscillate back-and-forth twice a day. We therefore investigate the general response of the creeping and seismogenic parts of MHT to periodic stresses in order to link these observations to physical parameters. First, the response of the creeping part of the MHT is analyzed with a simple spring-and-slider system bearing rate-strengthening rheology, and we show that at the transition with the locked zone, where the friction becomes near velocity neutral, the response of the slip rate may be amplified at some periods, which values are analytically related to the physical parameters of the problem. Such predictions therefore hold the potential of constraining fault properties on the MHT, but still await observational counterparts to be applied, as nothing indicates that the variations of seismicity rate on the locked part of the MHT are the direct expressions of variations of the slip rate on its creeping part, and no variations of the slip rate have been singled out from the GPS measurements to this day. When shifting to the locked seismogenic part of the MHT, spring-and-slider models with rate-weakening rheology are insufficient to explain the contrasted responses of the seismicity to the periodic loads that tides and monsoon both place on the MHT. Instead, we resort to numerical simulations using the Boundary Integral CYCLes of Earthquakes algorithm and examine the response of a 2D finite fault embedded with a rate-weakening patch to harmonic stress perturbations of various periods. We show that such simulations are able to reproduce results consistent with a gradual amplification of sensitivity as the perturbing period get larger, up to a critical period corresponding to the characteristic time of evolution of the seismicity in response to a step-like perturbation of stress. This increase of sensitivity was not reproduced by simple 1D-spring-slider systems, probably because of the complexity of the nucleation process, reproduced only by 2D-fault models. When the nucleation zone is close to its critical unstable size, its growth becomes highly sensitive to any external perturbations and the timings of produced events may therefore find themselves highly affected. A fully analytical framework has yet to be developed and further work is needed to fully describe the behavior of the fault in terms of physical parameters, which will likely provide the keys to deduce constitutive properties of the MHT from seismological observations.

Exploiting seismic waveforms of ambient noise and earthquakes

In this thesis, I apply detailed waveform modeling to study noise correlations in different environments, and earthquake waveforms for source parameters and velocity structure.

Green's functions from ambient noise correlations have primarily been used for travel-time measurement. In Part I of this thesis, by detailed waveform modeling of noise correlation functions, I retrieve both surface waves and crustal body waves from noise, and use them in improving earthquake centroid locations and regional crustal structures. I also present examples in which the noise correlations do not yield Green's functions, yet the results are still interesting and useful after case-by-case analyses, including non-uniform distribution of noise sources, spurious velocity changes, and noise correlations on the Amery Ice Shelf.

In Part II of this thesis, I study teleseismic body waves of earthquakes for source parameters or near-source structure. With the dense modern global network and improved methodologies, I obtain high-resolution earthquake locations, focal mechanisms and rupture processes, which provide critical insights to earthquake faulting processes in shallow and deep parts of subduction zones. Waveform modeling of relatively simple subduction zone events also displays new constraints on the structure of subducted slabs.

In summary, behind my approaches to the relatively independent problems, the philosophy is to bring observational insights from seismic waveforms in critical and simple ways.

Travel time tomographic studies of seismic structures around subducted lithospheric slabs

The nature of the subducted lithospheric slab is investigated seismologically by tomographic inversions of ISC residual travel times. The slab, in which nearly all deep earthquakes occur, is fast in the seismic images because it is much cooler than the ambient mantle. High resolution three-dimensional P and S wave models in the NW Pacific are obtained using regional data, while inversion for the SW Pacific slabs includes teleseismic arrivals. Resolution and noise estimations show the models are generally well-resolved.

The slab anomalies in these models, as inferred from the seismicity, are generally coherent in the upper mantle and become contorted and decrease in amplitude with depth. Fast slabs are surrounded by slow regions shallower than 350 km depth. Slab fingering, including segmentation and spreading, is indicated near the bottom of the upper mantle. The fast anomalies associated with the Japan, Izu-Bonin, Mariana and Kermadec subduction zones tend to flatten to sub-horizontal at depth, while downward spreading may occur under parts of the Mariana and Kuril arcs. The Tonga slab appears to end around 550 km depth, but is underlain by a fast band at 750-1000 km depths.

The NW Pacific model combined with the Clayton-Comer mantle model predicts many observed residual sphere patterns. The predictions indicate that the near-source anomalies affect the residual spheres less than the teleseismic contributions. The teleseismic contributions may be removed either by using a mantle model, or using teleseismic station averages of residuals from only regional events. The slab-like fast bands in the corrected residual spheres are are consistent with seismicity trends under the Mariana Tzu-Bonin and Japan trenches, but are inconsistent for the Kuril events.

The comparison of the tomographic models with earthquake focal mechanisms shows that deep compression axes and fast velocity slab anomalies are in consistent alignment, even when the slab is contorted or flattened. Abnormal stress patterns are seen at major junctions of the arcs. The depth boundary between tension and compression in the central parts of these arcs appears to depend on the dip and topology of the slab.

Dynamical models for the Earth's geoid

The Earth's largest geoid anomalies occur at the lowest spherical harmonic degrees, or longest wavelengths, and are primarily the result of mantle convection. Thermal density contrasts due to convection are partially compensated by boundary deformations due to viscous flow whose effects must be included in order to obtain a dynamically consistent model for the geoid. These deformations occur rapidly with respect to the timescale for convection, and we have analytically calculated geoid response kernels for steady-state, viscous, incompressible, self-gravitating, layered Earth models which include the deformation of boundaries due to internal loads. Both the sign and magnitude of geoid anomalies depend strongly upon the viscosity structure of the mantle as well as the possible presence of chemical layering.

Correlations of various global geophysical data sets with the observed geoid can be used to construct theoretical geoid models which constrain the dynamics of mantle convection. Surface features such as topography and plate velocities are not obviously related to the low-degree geoid, with the exception of subduction zones which are characterized by geoid highs (degrees 4-9). Recent models for seismic heterogeneity in the mantle provide additional constraints, and much of the low-degree (2-3) geoid can be attributed to seismically inferred density anomalies in the lower mantle. The Earth's largest geoid highs are underlain by low density material in the lower mantle, thus requiring compensating deformations of the Earth's surface. A dynamical model for whole mantle convection with a low viscosity upper mantle can explain these observations and successfully predicts more than 80% of the observed geoid variance.

Temperature variations associated with density anomalies in the man tie cause lateral viscosity variations whose effects are not included in the analytical models. However, perturbation theory and numerical tests show that broad-scale lateral viscosity variations are much less important than radial variations; in this respect, geoid models, which depend upon steady-state surface deformations, may provide more reliable constraints on mantle structure than inferences from transient phenomena such as postglacial rebound. Stronger, smaller-scale viscosity variations associated with mantle plumes and subducting slabs may be more important. On the basis of numerical modelling of low viscosity plumes, we conclude that the global association of geoid highs (after slab effects are removed) with hotspots and, perhaps, mantle plumes, is the result of hot, upwelling material in the lower mantle; this conclusion does not depend strongly upon plume rheology. The global distribution of hotspots and the dominant, low-degree geoid highs may correspond to a dominant mode of convection stabilized by the ancient Pangean continental assemblage.

Sri Lanka: is it a mid-plate platelet?

Two observations suggest the possibility that Sri Lanka is acting as a small-mid-plate platelet moving very slowly within and relative to the larger Indian plate. First, sediments of the Bengal Deep-Sea Fan off the SSE continental margin are folded and uplifted in a manner similar to the deformation from front of accretionary prisms where thick sediment columns are passing into subduction zones. And second, subsidence rates in the area of presumed spreading or continuing stretching of continental crust, the Cauvery-Palk Strait-Gulf of Mannar Basin, have not decreased during the Cenozoic as would be predicted by an aborted rift or aulacogen model, but instead appear to have accelerated during the Neogene. Information available on other phenomena which re predicted by the model is at the present time inadequate for evaluation.

Mechanical and hydraulic properties of Nankai accretionary prism sediments: Effect of stress path

We have conducted triaxial deformation experiments along different loading paths on prism sediments from the Nankai Trough. Different load paths of isotropic loading, uniaxial strain loading, triaxial compression (at constant confining pressure, Pc), undrained Pc reduction, drained Pc reduction, and triaxial unloading at constant Pc, were used to understand the evolution of mechanical and hydraulic properties under complicated stress states and loading histories in accretionary subduction zones. Five deformation experiments were conducted on three sediment core samples for the Nankai prism, specifically from older accreted sediments at the forearc basin, underthrust slope sediments beneath the megasplay fault, and overthrust Upper Shikoku Basin sediments along the frontal thrust. Yield envelopes for each sample were constructed based on the stress paths of Pc-reduction using the modified Cam-clay model, and in situ stress states of the prism were constrained using the results from the other load paths and accounting for horizontal stress. Results suggest that the sediments in the vicinity of the megasplay fault and frontal thrust are highly overconsolidated, and thus likely to deform brittle rather than ductile. The porosity of sediments decreases as the yield envelope expands, while the reduction in permeability mainly depends on the effective mean stress before yield, and the differential stress after yield. An improved understanding of sediment yield strength and hydromechanical properties along different load paths is necessary to treat accurately the coupling of deformation and fluid flow in accretionary subduction zones. © 2012 American Geophysical Union All Rights Reserved.

俯冲带深部流体特征——以西南天山变质俯冲杂岩为例

The subduction zone is an important site of the fluid activity and recycling of chemical elements. The fluid characteristic of deep subduction zones is a top scientific problem attracting the petrologists, geochemists and tectonists. In this dissertation, the characteristics of fluid activity within a deep subduction zone have been explored on the basis of the studies on the petrography, mineral chemistry, fluid inclusions, geochemistry and metamorphic P–T conditions of the omphacite-bearing high-pressure veins and related hosts from the low-temperature/high-pressure metamorphic belt in southwestern Tianshan, China. Multiple high-pressure veins are exposed in host eclogites and blueschists. The veins are composed predominantly of omphacite, garnet, quartz, and other minerals. Some veins contain cm-sized rutiles. In general, the vein can be divided into three types, the ‘in situ dehydration’ vein, the ‘external transport’ vein and the ‘composite’ vein. The omphacites within the veins and related host rocks contain lots of two-phase or three-phase primary fluid inclusions. The final melting temperature (Tfm) of fluid inclusions varies mainly from -0.6 to -4.3 °C, the homogeneous temperature (Th) varies from 185 to 251 °C, the salinity varies from 1.1 to 6.9 wt.% NaCl equivalent and the density varies from 0.81 to 0.9 g/cm3. The fluids were released under the conditions of T = 520–580°C and P = 15–19 kbar at blueschist facies to eclogite facies transition. The fluids include not only Li, Be, LILE, La, Pb-enriched and HFSE- and HREE-depleted aqueous fluids but also HFSE (Ti-Nb-Ta)-rich aqueous fluids. The complex composed of aluminosilicate polymers and F was the catalyst which had caused the Ti-Nb-Ta to be dissolved into the fluids. During the transport of the LILE-rich and HFSE- and HREE-poor fluids, they can exchange some chemical elements with country rocks and leach some trace elements in some extent. The rutile could be precipitated from the HFSE (Ti-Nb-Ta)-rich aqueous fluids when CO2 was added into the fluids. The host rocks could obtain some elements, such as Ca, Cs, Rb, Ba and Th, from the external fluids. The fluids with complex composition had been released within the deep subduction zone (>50 km) in Early Carboniferous during the subduction of the South Tianshan Ocean under the Yili–Central Tianshan Plate. The results obtained in this dissertation have made new progress compared with the published data (e.g. Tatsumi, 1989; Becker et al., 1999; Scambelluri and Philippot, 2001; Manning, 2004; Hermann et al., 2006; Spandler and Hermann, 2006).

基于三重震相波形模拟研究亚洲东北部上地幔过渡带速度结构

Both the global and regional P wave tomographic studies have revealed significant deep structural heterogeneities in subduction zone regions. In particular, low-velocity anomalies have been observed beneath the descending high-velocity slabs in a number of subduction zones. The limited resolution at large depths and possible trade-off between the high and low velocities, however, make it difficult to substantiate this feature and evaluate the vertical extent of the low-velocity structure. From broadband waveform modeling of triplicated phases near the 660-km discontinuity for three deep events, we constrained both the P and SH wave velocity structures around the base of the upper mantle in northeast Asia. For the two events beneath the southern Kurile, the rays traveled through the lowermost transition zone and uppermost lower mantle under the descending Pacific slab. Our preferred models consistently suggest normal-to-lower P and significantly low SH velocities above and below the 660-km discontinuity extending to about 760-km depth compared with the global IASP91 model, corroborating previous observations for a slow structure underneath the slab. In contrast, both high P and SH velocity anomalies are shown in our preferred models for the Japan subduction zone region, likely reflecting the structural feature of a slab stagnant above the 660-km discontinuity. The velocity jumps across the 660-km discontinuity were found to be on average 4.5% and 7% for P and S waves under the south Kurile, and 3% and 6% under the Japan subduction zone. The respective velocity contrasts in the two regions are consistent with mineralogical models for colder slab interior and hotter under-slab areas. Based on mineral physics data, the depth-averaged ~1.5% P and ~2.5% SH velocity differences in the depth range of 560-760 km between the two regions could be primarily explained by a 350~450K temperature variation, although the presence of about 0.5wt%~1wt% water might also contribute to the subtle velocity variations near the base of the transition zone in the southern Kurile. From our modeling results, we speculate that the slow structure in the southern Kurile may be correlated to the low velocity zone observed previously around the 410-km discontinuity under Northern Honshu. Both are probably associated with a thermal anomaly rooted in the lower mantle beneath the subduction zone in northeast Asia.

地幔对流理论及其在解释地球物理观测的应用

With the improving of mantle convection theory, the developing of computing method and increasing of the measurement data, we can numerically simulate more clearly about the effects on some geophysical observed phenomenons such as the global heat flow and global lithospheric stress field in the Earth's surface caused by mantle convection, which is the primary mechanism for the transport of heat from the Earth's deep interior to its surface and the underlying force mechanism of dynamics in the Earth.Chapter 1 reviews the historical background and present research state of mantle convection theory.In Chapter 2, the basic conception of thermal convection and the basic theory about mantle flow.The effects on generation and distribution of global lithospheric stres s field induced by mantle flow are the subject of Chapter 3. Mantle convection causes normal stress and tangential stresses at the bottom of the lithosphere, and then the sublithospheric stress field induces the lithospheric deformation as sixrface force and results in the stress field within the lithosphere. The simulation shows that the agreement between predictions and observations is good in most regions. Most of subduction zones and continental collisions are under compressive. While ocean ridges, such as the east Pacific ridge, the Atlantic ridge and the east African rift valley, are under tensile. And most of the hotspots preferentially occur in regions where calculated stress is tensile. The calculated directions of the most compressive principal horizontal stress are largely in accord with that of the observation except for some regions such as the NW-Pacifie subduction zone and Qinghai-Tibet Plateau, in which the directions of the most compressive principal horizontal stress are different. It shows that the mantel flow plays an important role in causing or affecting the large-scale stress field within the lithosphere.The global heat flow simulation based on a kinematic model of mantle convection is given in Chapter 4. Mantle convection velocities are calculated based on the internal loading theory at first, the velocity field is used as the input to solve the thermal problem. Results show that calculated depth derivatives of the near surface temperature are closely correlated to the observed surface heat flow pattern. Higher heat flow values around midocean ridge systems can be reproduced very well. The predicted average temperature as a function of function of depth reveals that there are two thermal boundary layers, one is close to the surface and another is close to the core-mantle boundary, the rest of the mantle is nearly isothermal. Although, in most of the mantle, advection dominates the heat transfer, the conductive heat transfer is still locally important in the boundary layers and plays an important role for the surface heat flow pattern. The existence of surface plates is responsible for the long wavelength surface heat flow pattern.In Chapter 5, the effects on present-day crustal movement in the China Mainland resulted from the mantle convection are introduced. Using a dynamic method, we present a quantitative model for the present-day crustal movement in China. We consider not only the effect of the India-Eurasia collision, the gravitational potential energy difference of the Tibet Plateau, but also the contribution of the shear traction on the bottom of the lithosphere induced by the global mantle convection. The comparison between our results and the velocity field obtained from the GPS observation shows that our model satisfactorily reproduces the general picture of crustal deformation in China. Numerical modeling results reveal that the stress field on the base of the lithosphere induced by the mantle flow is probably a considerable factor that causes the movement and deformation of the lithosphere in continental China with its eflfcet focuing on the Eastern China A numerical research on the small-scale convection with variable viscosity in the upper mantle is introduced in Chapter 6. Based on a two-dimensional model, small-scale convection in the mantle-lithosphere system with variable viscosity is researched by using of finite element method. Variation of viscosity in exponential form with temperature is considered in this paper The results show that if viscosity is strongly temperature-dependent, the upper part of the system does not take a share in the convection and a stagnant lid, which is identified as lithosphere, is formed on the top of system because of low temperature and high viscosity. The calculated surface heat flow, topography and gravity anomaly are associated well with the convection pattern, namely, the regions with high heat flow and uplift correspond to the upwelling flow, and vice versa.In Chapter 7, we give a brief of future research subject: The inversion of lateral density heterogeneity in the mantle by minimizing the viscous dissipation.

Source depletion and extent of melting in the Tongan sub-arc mantle

The fluid immobile High Field Strength Elements (HFSE) Nb and Ta can be used to distinguish between the effects of variable extents of melting and prior source depletion of the Tongan sub-arc mantle. Melting of spinel Iherzolite beneath the Lau Basin back-arc spreading centres has the ability to fractionate Nb from Ta due to the greater compatibility of the latter in clinopyroxene. The identified spatial variation in plate velocities and separation of melt extraction zones, combined with extremely depleted lavas make Tonga an ideal setting in which to test models for arc melt generation and the role of back-arc magmatism. We present new data acquired by laser ablation-ICPMS of fused sample glasses produced without the use of a melt fluxing agent. The results show an arc trend towards strongly sub-chondritic Nb/Ta (

Magma Evolution in the Primitive, Intra-oceanic Tonga Arc: Petrogenesis of Basaltic Andesites at Tofua Volcano

Tofua volcano is situated midway along the Tonga oceanic arc and has undergone two phases of ignimbrite-forming activity. The eruptive products are almost entirely basaltic andesites (52 center dot 5-57 wt % SiO2) with the exception of a volumetrically minor pre-caldera dacite. The suite displays a strong tholeiitic trend with K2O <1 wt %. Phenocryst assemblages typically comprise plagioclase + clinopyroxene +/- orthopyroxene with microlites of Ti-magnetite. Olivine (Fo(83-88)) is rare and believed to be dominantly antecrystic. An increase in the extent and frequency of reverse zoning in phenocrysts, sieve-textured plagioclase and the occurrence of antecrystic phases in post-caldera lavas record a shift to dynamic conditions, allowing the interaction of magma batches that were previously distinct. Pyroxene thermobarometry suggests crystallization at 950-1200 degrees C and 0 center dot 8-1 center dot 8 kbar. Volatile measurements of glassy melt inclusions indicate a maximum H2O content of 4 center dot 16 wt % H2O, and CO2-H2O saturation curves indicate that crystallization occurred at two levels, at depths of 4-5 center dot 5 km and 1 center dot 5-2 center dot 5 km. Major and trace element models suggest that the compositions of the majority of the samples represent a differentiation trend whereby the dacite was produced by 65% fractional crystallization of the most primitive basaltic andesite. Trace element models suggest that the sub-arc mantle source is the residuum of depleted Indian mid-ocean ridge basalt mantle (IDMM-1% melt), whereas radiogenic isotope data imply addition of 0 center dot 2% average Tongan sediment melt and a fluid component derived from the subducted altered Pacific oceanic crust. A horizontal array on the U-Th equiline diagram and Ra excesses of up to 500% suggest fluid addition to the mantle wedge within the last few thousand years. Time-integrated (Ra-226/Th-230) vs Sr/Th and Ba/Th fractionation models imply differentiation timescales of up to 4500 years for the dacitic magma compositions at Tofua.

Recent contribution of sediments and fluids to the mantle's volatile budget

Subduction modifies the cycling of Earth's volatile elements. Fluid-rich sediments and hydrated oceanic lithosphere enter the convecting mantle at subduction zones. Some of the sediments and volatile components are released from the subducting slab, promote mantle melting and are returned to the surface by volcanism. The remainder continue into the deeper mantle. Quantification of the fate of these volatiles requires an understanding of both the nature and timing of fluid release and mantle melting(1). Here we analyse the trace element and isotopic geochemistry of fragments of upper mantle rocks that were transported to the surface by volcanic eruptions above the Batan Island subduction zone, Philippines. We find that the mantle fragments exhibit extreme disequilibrium between their U-Th-Ra isotopic ratios, which we interpret to result from the interaction of wet sediment melts and slab-derived fluids with rocks in the overlying mantle wedge. We infer that wet sediments were delivered from the slab to the mantle wedge between 8,000 and 10,000 years ago, whereas aqueous fluids were delivered separately much later. We estimate that about 625 ppm of water is retained in the wedge. A significant volume of water could therefore be delivered to the mantle transition zone at the base of the upper mantle, or even to the deeper mantle.