The influence of thermochemical convection on the fixity of mantle plumes

A general feature of both isochemical and thermochemical studies of mantle convection is that horizontal plume velocities tend to be smaller than typical convective velocities, however, it is not clear which system leads to a greater fixity of mantle plumes. We perform two- and three-dimensional numerical calculations and compare both thermochemical and isochemical cases with similar convective vigor to determine whether presence of a dense component in the mantle can lead to smaller ratios of horizontal plume velocity to surface velocity. We investigate different viscosity and density contrasts between chemical components in the thermochemical calculations, and we perform isochemical calculations with both free-slip and no-slip bottom boundary conditions. We then compare both visually and quantitatively the results of the thermochemical and isochemical calculations to determine which leads to greater plume fixity. We find that horizontal plume velocities for thermochemical calculations are similar to those from isochemical calculations with no-slip bottom boundary conditions. In addition, we find that plumes tend to be more fixed for isochemical cases with free-slip bottom boundary conditions for two-dimensional calculations, however, in three dimensions, we find that plume fixity is similar to that observed in thermochemical calculations.     

Helium isotope composition of the early Iceland mantle plume inferred from the Tertiary picrites of West Greenland

Picrites from the 61 million year old Vaigat Formation of the Nuussuaq Peninsula in West Greenland have ³He/⁴He ratios trapped in olivine phenocrysts which range up to 30 times the atmospheric ratio. These high values, measured during gas extraction by crushing in vacuum, are similar to the highest magmatic ³He/⁴He ratios found in young terrestrial volcanic rocks. By analogy with young basalts, in which crushing selectively extracts magmatic helium, any significant cosmogenic ³He appears to be absent in these picrites. Additional evidence for the absence of cosmogenic helium is provided by fusion results on the crushed olivine powders and by a single stepwise crushing experiment, in which only magmatic and radiogenic helium components are resolvable. The West Greenland picrites have Pb, Nd and Sr isotope compositions which overlap those found in picrites from Iceland and in basalts from Loihi Seamount, localities which today also have high ³He/⁴He ratios. Isotopic variations in He, Pb, Nd and Sr for the West Greenland picrites are interpreted to largely result from interaction of the early Iceland mantle plume with the upper mantle during plume ascent and dispersion beneath the continental lithosphere. The presence of high ³He/⁴He ratios in West Greenland, and the onset of magmatism across the North Atlantic Volcanic Province near 62 Ma, supports the hypothesis for very rapid dispersion (>1 m/year) of mantle plume head material during the earliest stages of plume impact, as predicted in recent numerical simulations of plume behavior during thermal mantle convection with non-Newtonian rheology.     

The thermal state and composition of the lithospheric mantle beneath the Leizhou Peninsula, South China

Two localities on the Leizhou Peninsula, southern China (Yingfengling and Tianyang basaltic volcanoes) yield a wide variety of mantle-derived xenoliths including Cr-diopside series mantle wall rocks and two distinct types of Al-augite series pyroxenites. Metapyroxenites have re-equilibrated granoblastic microstructures whereas pyroxenites with igneous microstructures have not thermally equilibrated to the mantle conditions. An abundant suite of megacrysts and megacrystic aggregates (including garnet, plagioclase, clinopyroxene, ilmenite and apatite) is interpreted as the pegmatitic equivalents of the igneous pyroxenite suite. Layered spinel lherzolite/spinel websterite xenoliths were formed by metamorphic differentiation caused by mantle deformation, inferred to be related to lithospheric thinning. Some metapyroxenites have garnet websterite assemblages that allow calculation of their mantle equilibration temperatures and pressures and the construction of the first xenolith geotherm for the southernmost China lithosphere. Heat flow data measured at the surface in this region yield model conductive geotherms (using average crustal conductivity values) that are consistent with the xenolith geotherm for the mantle. The calculated mean surface heat flux is 110 mW/m². This high heat flux and the high geotherm are consistent with young lithospheric thinning in southern China, and with recent tomography results showing shallow low-velocity zones in this region. The xenolith geotherm allows the construction of a lithospheric rock type section for the Leizhou region; it shows that the crust–mantle boundary lies at about 30 km, consistent with seismic data, and that the lithosphere–asthenosphere boundary lies at about 100 km.     

Plume-lithosphere interactions in the ocean basins: constraints from the source mineralogy

Trace element relationships of near-primary alkalic lavas from La Grille volcano, Grande Comore, in the Indian Ocean, as well as those of the Honolulu volcanic series, Oahu, Hawaii, show that their sources contain amphibole and/or phlogopite. Small amounts of each mineral (2% amphibole in the source of La Grille and 0.5% phlogopite plus some amphibole in the source of the Honolulu volcanics) and a range of absolute degrees of partial melting from 1 to 5% for both series are consistent with the observed trace element variation. Amphibole and phlogopite are not stable at the temperatures of convecting upper mantle or upwelling thermal plumes from the deep mantle; however, they are stable at pressure-temperature conditions of the oceanic lithospheric mantle. Therefore, the presence of amphibole and/or phlogopite in the magma source region of volcanics is strong evidence for lithospheric melting, and we conclude that the La Grille and the Honolulu series formed by melting of the oceanic lithospheric mantle.

The identification of amphibole ± phlogopite in the source region of both series implies that the metasomatism by fluids or volatile-rich melts occurred prior to melting. The presence of hydrous phases results in a lower solidus temperature of the lithospheric mantle, which can be reached by conductive heating by the thermal plumes. Isotope ratios of the La Grille and the Honolulu series display a restricted range in composition and represent compositional end-members for each island. Larger isotopic variations in shield lavas, represented by the contemporaneous Karthala volcano on Grande Comore and the older Koolau series on Oahu, reflect interaction of the upwelling thermal plumes with the lithospheric mantle rather than the heterogeneity of deep-seated mantle plume sources or entrainment of mantle material in the rising plume. Literature OsSr isotope ratio covariations constrain the process of plume-lithosphere interaction as occurring through mixing of plume melts and low-degree melts from the metasomatized oceanic lithospheric mantle.

The characterization of the lithospheric mantle signature allows the isotopic composition of the deep mantle plume components to be identified, and mixing relationships show that the Karthala and Koolau plume end-members have nearly uniform isotopic compositions. Based on independent arguments, isotopic variations on Heard and Easter islands have been shown to be a result of mixing between deep plume sources having distinct isotopic compositions with lithosphere or shallow asthenospheric mantle. To the extent that these case studies are representative of oceanic island volcanism, they indicate that interaction with oceanic lithospheric mantle plays an important role in the compositions of lavas erupted during the shield-building stage of plume magmatism, and that isotopic compositions of deep mantle plume sources are nearly uniform on the scale that they are sampled by melting.     

Transitions in the style of mantle convection at high Rayleigh numbers

The pattern and style of mantle convection govern the thermal evolution, internal dynamics, and large-scale surface deformation of the terrestrial planets. In order to characterize the nature of heat transport and convective behaviour at Rayleigh numbers, Ra, appropriate for planetary mantles (between 10⁴ and 10⁸), we perform a set of laboratory experiments. Convection is driven by a temperature gradient imposed between two rigid surfaces, and there is no internal heating. As the Rayleigh number is increased, two transitions in convective behaviour occur. First we observe a change from steady to time-dependent convection at Ra≈10⁵. A second transition occurs at higher Rayleigh numbers, Ra≈5×10⁶, with large-scale time-dependent flow being replaced by isolated rising and sinking plumes. Corresponding to the latter transition, the exponent β in the power law relating the Nusselt number Nu to the Rayleigh number (NuRa^β) is reduced. Both rising and sinking plumes always consist of plume heads followed by tails. There is no characteristic frequency for the formation of plumes.     

Contrasting rifted margin styles south of Greenland: implications for mantle plume dynamics

We present new and reprocessed seismic reflection data from the area where the southeast and southwest Greenland margins intersected to form a triple junction south of Greenland in the early Tertiary. During breakup at 56 Ma, thick igneous crust was accreted along the entire 1300-km-long southeast Greenland margin from the Greenland Iceland Ridge to, and possibly 100 km beyond, the triple junction into the Labrador Sea. However, highly extended and thin crust 250 km to the west of the triple junction suggests that magmatically starved crustal formation occurred on the southwest Greenland margin at the same time. Thus, a transition from a volcanic to a non-volcanic margin over only 100–200 km is observed. Magmatism related to the impact of the Iceland plume below the North Atlantic around 61 Ma is known from central-west and southeast Greenland. The new seismic data also suggest the presence of a small volcanic plateau of similar age close to the triple junction. The extent of initial plume-related volcanism inferred from these observations is explained by a model of lateral flow of plume material that is guided by relief at the base of the lithosphere. Plume mantle is channelled to great distances provided that significant melting does not take place. Melting causes cooling and dehydration of the plume mantle. The associated viscosity increase acts against lateral flow and restricts plume material to its point of entry into an actively spreading rift. We further suggest that thick Archaean lithosphere blocked direct flow of plume material into the magma-starved southwest Greenland margin while the plume was free to flow into the central west and east Greenland margins. The model is consistent with a plume layer that is only moderately hotter, 100–200°C, than ambient mantle temperature, and has a thickness comparable to lithospheric thickness variations, 50–100 km. Lithospheric architecture, the timing of continental rifting and viscosity changes due to melting of the plume material are therefore critical parameters for understanding the distribution of magmatism.     

Lowermost mantle anisotropy beneath the Pacific: Imaging the source of the Hawaiian plume

We utilized recordings of seismic shear phases provided by several North American broadband seismometer arrays to provide unique constraints on shear wave anisotropy beneath the northern and central Pacific Ocean. Using a new analysis method that reduces measurement errors and enables the analysis of a larger number of available waveforms, we examined relative travel times of teleseismic S and Sdiff that sample a large area of lowermost mantle structure. The results of this study provide evidence for small-scale lateral and depth variations in shear wave anisotropy for a broad region of the lowermost mantle beneath the Pacific Ocean. In particular, we image a localized zone of anomalously strong anisotropy whose strength increases toward the top of D″ beneath Hawaii. Our results, combined with a previous study of V_P/V_SH ratios, indicate that ancient subducted slab material may be responsible for observations beneath the northern Pacific, while lenses or layers of core–mantle boundary reaction products or partial melt, oriented by horizontal inflow of mantle material to the Hawaiian plume source, can explain observations beneath the central Pacific.     

The effect of mantle composition on density in the extending lithosphere

The density distribution of the lithosphere is non-linear and discontinuous due to complex mineralogy and, most importantly, phase transitions. We evaluate the influence of changes in mantle composition on lithospheric density and its evolution during horizontal stretching, using thermodynamic calculations of the density as a function of pressure, temperature and composition. We also develop a simple parameterization based on end-member mineral reactions and geometric relationships between the geotherm and the phase boundary for comparison. The garnet–spinel peridotite transition leads to a moderate decrease in density of the mantle part of the lithospheric column at the initial stages of stretching. When the crust is sufficiently thinned and temperature is relatively high, plagioclase peridotite becomes stable in the upper part of the mantle. The density reduction due to the plagioclase-in reaction is controlled by bulk Al₂O₃ in the mantle and by the depth of the plagioclase-in reaction, which is mainly governed by the Na₂O/Al₂O₃ ratio. Since Na₂O and Al₂O₃ increase with the fertility of the mantle the phase transition effect is most pronounced for relatively fertile mantle (and strong extension) and can lead to 2.3% density reduction. This is equivalent to heating the entire lithosphere by 700 °C if only the effect of thermal expansion on density is taken into account. The formation of plagioclase peridotite can explain syn-rift uplift in sedimentary basins that experienced large mantle stretching without invoking an unrealistically strong increase in temperature. It might also be responsible for the break-up unconformity observed at continental margins.     

Seismic structure of the upper mantle beneath the western Philippine Sea

We investigate the seismic structure of the western Philippine Sea using two sets of seismological observations: ScS reverberations, which provide the layering framework for a regional upper mantle model, and observations of frequency-dependent phase delays for direct S waves, surface-reflected phases (sS, SS, sSS), and surface waves (R₁, G₁), which constrain the velocity and anisotropy structure within the layers. The combined data set, comprising 17 discontinuity amplitudes and layer travel times from the ScS-reverberation stack and more than 1000 frequency-dependent phase delays, was inverted for a path-averaged, radially anisotropic model. Mineralogical estimates of the bulk sound velocity and density are incorporated as complementary constraints. The final model, PHB3, is characterized by a 11.5-km thick crust, an anisotropic lid bounded by a sharp negative G discontinuity at 89 km, an anisotropic low-velocity layer extending to 166 km, a subjacent high-gradient region, and transition-zone discontinuities at depths of 408 km, 520 km and 664 km. The lid is slower than in a comparable model for the Tonga–Hawaii corridor (PA5), but also significantly thicker, requiring a compositional variation between the two regions. We explore the hypothesis that the thickness of the oceanic lid is controlled by the melting depth at the spreading centers during crust formation, and that the thicker crust and lid in the Philippine Sea results from deeper melting owing to a higher potential temperature and perhaps a higher water content in the upper mantle.     

On the dynamics of 3-D single thermal plumes at various Prandtl numbers and Rayleigh numbers

Three-dimensional (3-D) numerical simulations of single turbulent thermal plumes in the Boussinesq approximation are used to understand more deeply the interaction of a plume with itself and its environment. In order to do so, we varied the Rayleigh and Prandtl numbers from Ra?～?10⁵ to Ra?～?10⁸ and from Pr?～?0.025 to Pr?～?70. We found that thermal dissipation takes place mostly on the border of the plume. Moreover, the rate of energy dissipation per unit mass ε _T has a critical point around Pr?～?0.7. The reason is that at Pr greater than ～0.7, buoyancy dominates inertia and thermal advection dominates wave formation whereas this trend is reversed at Pr less than ～0.7. We also found that for large enough Prandtl number (Pr?～?70), the velocity field is mostly poloidal although this result was known for Rayleigh–Bénard convection (see Schmalzl et al. [On the validity of two-dimensional numerical approaches to time-dependent thermal convection. Europhys. Lett. 2004, 67, 390--396]). On the other hand, at small Prandtl numbers, the plume has a large helicity at large scale and a non-negligible toroidal part. Finally, as observed recently in details in weakly compressible turbulent thermal plume at Pr?=?0.7 (see Plourde et al. [Direct numerical simulations of a rapidly expanding thermal plume: structure and entrainment interaction. J. Fluid Mech. 2008, 604, 99--123]), we also noticed a two-time cycle in which there is entrainment of some of the external fluid to the plume, this process being most pronounced at the base of the plume. We explain this as a consequence of calculated Richardson number being unity at Pr?=?0.7 when buoyancy balance inertia.     

Reconciling strong slab pull and weak plate bending: The plate motion constraint on the strength of mantle slabs

Although subducting slabs undergo a bending deformation that resists tectonic plate motions, the magnitude of this resistance is not known because of poor constraints on slab strength. However, because slab bending slows the relatively rapid motions of oceanic plates, observed plate motions constrain the importance of bending. We estimated the slab pull force and the bending resistance globally for 207 subduction zone transects using new measurements of the bending curvature determined from slab seismicity. Predicting plate motions using a global mantle flow model, we constrain the viscosity of the bending slab to be at most ~ 300 times more viscous than the upper mantle; stronger slabs are intolerably slowed by the bending deformation. Weaker slabs, however, cannot transmit a pull force sufficient to explain rapid trenchward plate motions unless slabs stretch faster than seismically observed rates of ~ 10^− 15 s^− 1. The constrained bending viscosity (~ 2 × 10²³ Pa s) is larger than previous estimates that yielded similar or larger bending resistance (here ~ 25% of forces). This apparent discrepancy occurs because slabs bend more gently than previously thought, with an average radius of curvature of 390 km that permits subduction of strong slabs. This gentle bending may ultimately permit plate tectonics on Earth.     

Coupled Os–Os enrichments in the Earth’s mantle – core–mantle interaction or recycling of ferromanganese crusts and nodules?

¹⁸⁶Os enrichments in volcanic rocks and peridotite-derived iridosmine grains have been attributed to contributions from Earth’s outer core to the mantle, and apparently constrain the scale of mantle convection and an early timing for inner–outer core segregation more than 3.5 Gyr ago. Here, we highlight that marine ferromanganese crusts and nodules are characterised by high Pt/Os ratios and Pt–Os contents that develop much larger ¹⁸⁶Os excesses over geological time (≥0.2%/Gyr) than those hypothesised for Earth’s outer core (<0.005–0.01%/Gyr). ¹⁸⁷Os/¹⁸⁸Os ratios in ferromanganese crusts are radiogenic due to sequestering of continental Os from seawater. Similarly, ancient ferromanganese materials may have had ¹⁸⁶Os excesses (>0.1%) as a result of high Pt/Os ratios in continental crust, even prior to in-growth of ¹⁸⁶Os after formation due to their high Pt/Os ratios. Past recycling of small amounts of these materials into the Earth’s mantle will produce coupled ¹⁸⁷Os–¹⁸⁶Os excesses and little change in Re and platinum-group-element concentrations, as observed in Hawaiian picrites, and in contrast to the predicted result of outer core addition to the mantle. ¹⁸⁷Os and ¹⁸⁶Os enrichments in the Hawaiian mantle source are potentially consistent with it comprising recycled oceanic lithosphere, pelagic sediments and ferromanganese materials, and questions the notion that Os isotopes can be used to uniquely identify core–mantle interactions and the depth at which mantle sources for volcanism originate.     

Upper mantle flow beneath the Northwest of China and its lithospheric dynamics

ＵｐｐｅｒｍａｎｔｌｅｆｌｏｗｂｅｎｅａｔｈｔｈｅＮｏｒｔｈｗｅｓｔｏｆＣｈｉｎａａｎｄｉｔｓｌｉｔｈｏｓｐｈｅｒｉｃｄｙｎａｍｉｃｓＪＩＡＮ－ＨＵＡＨＵＡＮＧＩ（黄建华）；ＸＩＡ－ＨＵＡＣＨＡＮＧＩ（常筱华）ａｎｄＲＯＮＧ－ＳＨＡＮＦＵＩ傅容珊）（...     

Structure of hydrothermal plumes at the Logatchev vent field, 14°45′N, Mid-Atlantic Ridge: evidence from geochemical and geophysical data

In the Seventh cruise of R/V “Professor Logatchev” anomalies of natural electric field (EF), Eh and pS were discovered using a towed instrument package (RIFT) at 14°45′N on the MAR (Logatchev hydrothermal field). The anomalous zone (AZ) is situated close (10–35 m) to two low-temperature venting areas of degrading sulphides and a black smoker (Irina-Microsmoke) forming a distinct buoyant plume. Over or close to the main area of high-temperature venting situated to the south-east from the AZ, no EF or Eh anomalies were observed. According to the results of Mir dives the highly mineralised solutions from smoking craters at the main mound mostly form non-buoyant plumes (reverse-plumes). The buoyant plume structure shows the differentiation of the electrical and Eh fields within the plume. Maxima of the EF, Eh and E_H2S anomalies were revealed in the lower part (15 m) of the plume. The negative redox potential plume coupled with a sulphide anomaly is more localized in comparison with the EF. This observation indicates a distinct change in the composition of buoyant plume water, which may be due to the formation and fallout of early formed Fe sulphide particles soon after venting.     

The evolution of He Isotopes in the convecting mantle and the preservation of high He/He ratios

A key requirement for any model of mantle evolution is accounting for the high ³He/⁴He ratios of many ocean island basalts compared to those of mid-ocean ridge basalts. The early, popular paradigm of primitive, undegassed mantle stored in a convectively isolated lower mantle is incompatible with geophysical constraints that imply whole mantle convection. Thus it has been suggested more recently that domains with high ³He/U ratios have been created continuously from the bulk mantle throughout Earth history. Such models require that the ³He/⁴He ratio of the convecting mantle was at least as high as the highest values seen in OIB at the time the OIB source was generated. These domains must also be created with sufficient He to impart distinctive He isotopic signatures to ocean island basalts. However, the He isotope evolution of the mantle has not been consistently quantified to determine if such scenarios are plausible.

Here a simple model of the He evolution of the whole mantle is examined. Using a wide range of possible histories of continental extraction and He degassing, the bulk convecting mantle was found to have had ³He/⁴He ratios as high as those seen in the Iceland hotspot only prior to 3 Ga. Such high ³He/⁴He ratios can only be preserved if located in domains that are not modified by convective mixing or diffusive homogenisation since that time. Further, there are difficulties in producing, with commonly invoked magmatic processes, domains with sufficiently high ³He/U ratios and enough ³He to be able to impart this signature to ocean island basalts. The results are consistent with models that store such He signatures in the core or a deep layer in the mantle, but are hard to reconcile with models that continuously generate high ³He/⁴He domains within the mantle.     

Intracratonic asthenosphere upwelling and lithosphere rejuvenation beneath the Hoggar swell (Algeria): Evidence from HIMU metasomatised lherzolite mantle xenoliths

The mantle xenoliths included in Quaternary alkaline volcanics from the Manzaz-district (Central Hoggar) are proto-granular, anhydrous spinel lherzolites. Major and trace element analyses on bulk rocks and constituent mineral phases show that the primary compositions are widely overprinted by metasomatic processes. Trace element modelling of the metasomatised clinopyroxenes allows the inference that the metasomatic agents that enriched the lithospheric mantle were highly alkaline carbonate-rich melts such as nephelinites/melilitites (or as extreme silico-carbonatites). These metasomatic agents were characterized by a clear HIMU Sr–Nd–Pb isotopic signature, whereas there is no evidence of EM1 components recorded by the Hoggar Oligocene tholeiitic basalts. This can be interpreted as being due to replacement of the older cratonic lithospheric mantle, from which tholeiites generated, by asthenospheric upwelling dominated by the presence of an HIMU signature. Accordingly, this rejuvenated lithosphere (accreted asthenosphere without any EM influence), may represent an appropriate mantle section from which deep alkaline basic melts could have been generated and shallower mantle xenoliths sampled, respectively. The available data on lherzolite xenoliths and alkaline lavas (including He isotopes, Ra < 9) indicate that there is no requirement for a deep plume anchored in the lower mantle, and that sources in the upper mantle may satisfactorily account for all the geochemical/petrological/geophysical evidence that characterizes the Hoggar swell. Therefore the Hoggar volcanism, as well as other volcanic occurrences in the Saharan belt, are likely to be related to passive asthenospheric mantle uprising and decompression melting linked to tensional stresses in the lithosphere during Cenozoic reactivation and rifting of the Pan–African basement. This can be considered a far-field foreland reaction of the Africa–Europe collisional system since the Eocene.     

中国地幔结构及物性研究的进展

对2004——2007年中国地球物理学家在地幔内部结构和物质性质方面的研究工作进行了总结. 对地幔结构的地震波速度成像,地幔介质的各向异性,地幔间断面,地幔对流,以及地幔介质物性进行了综述,指出了各方面的主要内容,使用的主要方法和主要结果. 从这4年的研究可以发现,一些原有的研究领域工作更加深入,方法更加先进,而且进行了广泛的国际合作,合作的范围也逐渐扩大,方式多样,并出现了一些新的研究方向.      

Geochemistry and Sr-Nd-Pb isotopic characteristics of the Mugouriwang Cenozoic volcanic rocks from Tibetan Plateau:Constraints on mantle source of the underplated basic magma

The Mugouriwang Cenozoic volcanic rocks exposed in the north Qiangtang Block of Tibetan Plateau are mainly composed of basalt and andesitic-basalt,both characterized by the lower SiO2 (51%―54%),high refractory elements (i.e. Mg,Cr,Ni) as well as the moderate enrichment in light rare earth elements (LREE) relative to a slight depleted in Eu and high strength field elements (HFSE,i.e. Nb,Ta,Ti). Be-sides,the fairly low Sm/Yb value (3.07―4.35) could signify that the rocks should be derived directly from partial melting of the spinel lherzolite at the upper part of the asthenosphere. These rocks have radiogenic Sr and Pb (87Sr/86Sr = 0.705339 to 0.705667; 208Pb/204Pb = 38.8192 to 38.8937; 207Pb/204Pb = 15.6093 to 15.6245; 206Pb/204Pb = 18.6246 to 18.6383),and non-radiogenic Nd (143Nd/144Nd = 0.512604 to 0.512639; εNd = 0.02 to -0.66) in agreement with those values of the BSE mantle reservoir. The DUPAL anomaly of the rocks can be evidently attested by the △8/4Pb = 66.82 to 74.53 ,△7/4Pb = 9.88 to 11.42,△Sr>50,implying that the Mugouriwang volcanic rock is likely to be generated by partial melting of a Gondwana-bearing asthenospheric mantle ever matasomatised by the fluid from subduction zone. Depending on the previous study on the high-K calc-alkaline intermediate-felsic volcanics in the study area,this paper proposed that the fluids derived from the subducted Lhasa Block metasomatised the asthenosphere beneath the Qiangtang Block,and induced its partial melting,and then the melt under-plated the thickened Qiangtang lithosphere and caused the generation of the Cenozoic adakite-like felsic magmas in the Qiangtang region.     

Viscoelastic relaxation of the upper mantle and afterslip following the 2014 MW8.1 Iquique earthquake

An improved understanding of postseismic crustal deformation following large subduction earthquakes may help to better understand the rheological properties of upper mantle and the slip behavior of subduction interface. Here we construct a three-dimensional viscoelastic finite element model to study the postseismic deformation of the 2014 M_W8.1 Iquique, Chile earthquake. Elastic units in the model include the subducting slab, continental and oceanic lithospheres. Rheological units include the mantle wedge, the oceanic asthenosphere and upper mantle. We use a 2 ​km thick weak shear zone attached to the subduction fault to simulate the time-dependent stress-driven afterslip. The viscoelastic relaxation in the rheological units is represented by the Burgers rheology. We carry out grid-searches on the shear zone viscosity, thickness and viscosity of the asthenosphere, and they are determined to be 10¹⁷ ​Pa ​s, 110 ​km and 2 ​× ​10¹⁸ ​Pa ​s, respectively. The stress-driven afterlsip within the first two years is up to ~47 ​cm and becomes negligible after two years (no more than 5 ​cm/yr). Our results suggest that a thin, low-viscosity oceanic asthenosphere together with a weak shear zone attached to the fault are required to better reproduce the observed postseismic deformation.