The effect of the large-scale mantle flow field on the Iceland hotspot track

Fluid dynamical simulations were carried out in order to investigate the effect of the large-scale mantle flow field and the depth of the plume source on the structure of the Iceland plume through time. The time-dependent location and shape of the plume in the Earth's mantle was calculated in a global model and it was refined in the upper mantle using a 3D Cartesian model box. Global flow was computed based on density heterogeneities derived from seismic tomography. Plate motion history served as a velocity boundary condition in both models. Hotspot tracks of the plume conduits and the plume head were calculated and compared to actual bathymetry of the North Atlantic. If a plume source in the lowermost mantle is assumed, the calculated surface position of the plume conduit has a southward component of motion due to southward flow in the lower mantle. Depending on tomography model, assumed plume age and buoyancy the southward component is more or less dominating. Plume models having a source at the 660 km discontinuity are only influenced by flow in the upper mantle and transition zone and hence rather yield westward hotspot motion. Many whole-mantle plume models result in a V-shaped track, which does not match the straight Greenland–Iceland–Faroe ridge. Models without strong southward motion, such as for a plume source at 660 km depth, match actual bathymetry better. Plume tracks were calculated from both plume conduits and plume heads. A plume head of 120 K anomalous temperature gives the best match between plume head track and bathymetry.     

Geoelectric investigation of the Hellenic subduction zone using long period magnetotelluric data

The strongest evidence up to date for a subduction zone in the Hellenic region is a clearly identified Wadati-Benioff zone below the central Aegean Sea, to a maximum depth of 180 km. Alternative seismic tomography models suggest that subduction process continues deeper than the Wadati-Benioff zone to a maximum depth of at least 600 km. So far the lack of deep electrical studies in the region impeded scientists from imposing other control factors than seismic to the proposed models for the Hellenic Subduction Zone (HSZ). A Long Period Magnetotelluric (LMT) study was carried out in the southern part of the Greek mainland to study the deep electrical characteristics of the HSZ and examine whether prominent modelled features correlate with structures identified by the seismic methods. The study comprised collection, processing and modelling of magnetotelluric (MT) data in the period range 100–10000 s from ten sites located along a 250 km NE–SW trending profile. The dimensionality of the data was examined at a pre-modelling stage and it was found that they do not exhibit three-dimensional (3-D) features. The latter enabled to construct both one-dimensional (1-D) and two-dimensional (2-D) models. The proposed geoelectric model for HSZ was based on 2-D modelling, since it had better maximum depth resolution of about 400 km, and revealed structures not detected by 1-D modelling attempts. The model structure which was related to the African and Euro-Asian lithosphere is relatively resistive (> 800 Ω-m) and has an average thickness of 150–170 km. Although the bottom of the lithosphere is adequately resolved, the Wadati-Benioff zone that delineates the top of the subducting lithospheric slab is not identified by any electrical feature. The modelled structure associated with the subducting part of the African lithosphere penetrates a relatively conductive (< 200 Ω-m) asthenosphere with a dip angle of 42°. Intermediate electrical resistivities (200–800 Ω-m) are attributed to the ascending melting part of the lithosphere below the region of the Hellenic Volcanic Arc (HVA) and to a dipping zone below the south-western part of the profile, at 170–220 km depths.     

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

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

Implications for the lithospheric structure of Cambay rift zone,western India:Inferences from a magnetotelluric study

Broad-band and long-period magnetotelluric(MT) data were acquired along an east-west trending traverse of nearly 200 km across the Kachchh,Cambay rift basins,and Aravalli-Delhi fold belt(ADFB),western India.The regional strike analysis of MT data indicated an approximate N59°E geoelectric strike direction under the traverse and it is in fair agreement with the predominant geological strike in the study area.The decomposed transverse electric(TE)-and transverse magnetic(TM)-data modes were inverted using a nonlinear conjugate gradient algorithm to image the electrical lithospheric structure across the Cambay rift basin and its surrounding regions.These studies show a thick(～1-5 km) layer of conductive Tertiary-Mesozoic sediments beneath the Kachchh and Cambay rift basins.The resistive blocks indicate presence of basic/ultrabasic volcanic intrusives,depleted mantle lithosphere,and different Precambrian structural units.The crustal conductor delineated within the ADFB indicates the presence of fluids within the fault zones,sulfide mineralization within polyphase metamorphic rocks,and/or Aravalli-Delhi sediments/metasediments.The observed conductive anomalies beneath the Cambay rift basin indicate the presence of basaltic underplating,volatile(CO_2,H_2 O) enriched melts and channelization of melt fractions/fluids into crustal depths that occurred due to plume-lithosphere interactions.The variations in electrical resistivity observed across the profile indicate that the impact of Reunion plume on lithospheric structures of the Cambay rift basin is more dominant at western continental margin of India(WCMI) and thus support the hypothesis proposed by Campbell Griffiths about the plume-lithosphere interactions.     

Hawaiian volcanic propagation and Hawaiian swell asymmetry: evidence of northwestward flow of the deep upper mantle

Bathymetry and the geoid anomaly of the northern flank of the Hawaiian swell is broader and higher than the southern flank, and it is characterized by higher heat flow than the axis or southern flank. It is here proposed that the northern flank of the Hawaiian swell has been augmented by heat conducted from the hotspot conduit into the upper mantle then transported northward of the volcanic axis by flow in the upper mantle (325°) that is more northerly than Pacific plate motion (292°). By assuming that the deep upper mantle is decoupled from the Pacific plate and is flowing at 325° to the northwest, changes in direction and rate of volcanic propagation and in geochemistry along individual volcanic segments of the Hawaiian volcanic chain can be interpreted in terms of tank experiment results showing that a volcanic hotspot conduit breaks into diapirs when tilted by mantle flow. Hawaiian volcanoes are aligned in en-echelon segments, and the Hawaiian Islands are the two most recent segments. For an individual segment, older northwestern volcanoes are aligned nearly parallel to the 292° plate motion direction, and they propagated to the southeast at approximately the same rate as the 92 km/m.y. speed of northwestward plate motion. In contrast, the alignment of the younger southeastern volcanoes is close to 325°, and they show a conspicuous acceleration in propagation of volcanism marked by out-of-sequence eruptions. Within the model proposed here, diapirs rise from instability nodes that develop along the tilted conduit of a mantle hotspot plume as it is sheared in the direction of deep upper-mantle flow and each diapir gives rise to a single volcanic center. As tilting progresses, diapirs form at lower levels along the conduit in more upstream positions of the mantle flow zone, rise sequentially into the decoupled lithosphere, erupt sequentially, and are translated in the direction of plate motion (older, northwestern Hawaiian Islands). Eventually, flow in the highly tilted conduit is impeded to the degree that the remaining upstream conduit breaks into a number of diapirs that rise together into the lithosphere. These late diapirs, translated as a group aligned in the direction of horizontal mantle flow, erupt over a relatively short time span and show out-of-sequence volcanism (younger, southeastern Hawaiian Islands). At this stage, a new cycle of rising and tilting will initiate the next en-echelon segment.     

Significance of rare hydrous alkaline melts in Hawaiian xenoliths

 Two types of melt pockets occur in Hawaiian mantle xenoliths: amphibole-bearing (AMP) and spinel-bearing (SMP). AMPs contain amphibole (kaersutite), olivine (Fo₉₂), clinopyroxene (with 7–11% Al₂O₃), vesicles and glass. SMPs contain olivine, clinopyroxene, spinel, glass, and vesicles. The glasses in SMPs (SiO₂=44–45%, 11–12% alkalis, La=90–110 ppm) and AMPs (SiO₂=49–54%, 6–8.5% alkalis, La=8–14 ppm) are distinct in color and composition. Both glasses are generally characterized by LREE-enriched (chondrite-normalized) patterns. Amphibole and clinopyroxene have gently convex upward-to-moderately LREE-enriched patterns. Mineral/glass trace element abundance ratio plots show a strong negative Ti anomaly and a gentle negative Zr anomaly for clinopyroxene/glass; whereas amphibole/glass patterns show a distinctive positive Ti spike. The amphibole/glass trace element ratios are similar to published megacryst/lava values. An earlier study showed that the Hawaiian spinel lherzolites (lithosphere) have largely been metasomatized during post-erosional Honolulu magmatic activity. REE abundances of SMP glasses (melts) overlap the REE abundances calculated for such metasomes. The occurrence of hydrous, alkaline, mafic melt pockets in Hawaiian upper mantle xenoliths implies that (1) such hydrous liquids are generated in the upper mantle, and (2) water plays a role in magmatic activity associated with the Hawaiian plume. Although we are uncertain about the source (plume, lithosphere, or asthenosphere) of this water, we speculate that such melts and other alkalic lavas erupted on Oahu and on the sea-floor over the Hawaiian arch were generated from a broad „wet“ rim of a radially layered Hawaiian plume, whose hot and „dry“ core supplied the shield-forming magmas. Received: 6 February 1995 / Accepted: 28 August 1995     

The state of the upper mantle beneath southern Africa

We present a new upper mantle seismic model for southern Africa based on the fitting of a large (3622 waveforms) multi-mode surface wave data set with propagation paths significantly shorter (≤ 6000 km) than those in globally-derived surface wave models. The seismic lithosphere beneath the cratonic region of southern Africa in this model is about 175 ± 25 km thick, consistent with other recent surface wave models, but significantly thinner than indicated by teleseismic body-wave tomography. We determine the in situ geotherm from kimberlite nodules from beneath the same region and find the thermal lithosphere model that best fits the nodule data has a mechanical boundary layer thickness of 186 km and a thermal lithosphere thickness of 204 km, in very good agreement with the seismic measurement. The shear wave velocity determined from analyzes of the kimberlite nodule compositions agree with the seismic shear wave velocity to a depth of 150 km. However, the shear wave velocity decrease at the base of the lid seen in the seismic model does not correspond to a change in mineralogy. Recent experimental studies of the shear wave velocity in olivine as a function of temperature and period of oscillation demonstrate that this wave speed decrease can result from grain boundary relaxation at high temperatures at the period of seismic waves. This decrease in velocity occurs where the mantle temperature is close to the melting temperature (within 100 °C).     

Magnetotelluric images of the crust and mantle in the southwestern Taurides, Turkey

Inversion of the magnetotelluric data across the southwestern Taurides reveals two subzones of crust with varying thicknesses: conductive lower crust (<75 Ω m), overlain by resistive (>350 Ω m) upper crust, with four resistive cores (>2000 Ω m) separated by three relatively conductive vertical zones. The first and second vertical zones coincide with surface faults interpreted in Anatolia, such as Fethiye Burdur Fault Zone. The third one is the most conductive and lies in continuity with the Strabo Fault Zone in the Mediterranean Sea. A hypocentral cross section of earthquakes along the profile shows more dense seismic activity in the second resistive core where the conductive crust is not present beneath it. The depth of the crust/upper mantle boundary varies between 30 and 50 km and has an undulating character. The resistivity of the upper mantle reaches 500–1000 Ω m.     

North Arch volcanic fields near Hawaii are evidence favouring the restite‐root hypothesis for the origin of hotspot swells

A pronounced regional bathymetric swell is a common feature of oceanic hotspot volcanism. Recently, the hypothesis of a buoyant sublithospheric swell‐root has been favoured, the root being either a ‘refracted plume’ of hot, buoyant and hence low‐viscosity plume material embedded within surrounding higher viscosity asthenosphere, or a ‘restite‐root’ composed of hot and more viscous residues to hotspot melt‐extraction from the hottest central portion of the upwelling plume. In this article, we present numerical experiments that show that these two scenarios predict different flow and melting patterns whenever the plume passes beneath an obliquely oriented fracture zone where the base of the overriding lithosphere changes in thickness. We find that the restite‐root hypothesis predicts an asymmetric spatial pattern of persistent ‘arch volcanism’ strikingly similar to that found in the Hawaiian chain surrounding the Molokai Fracture Zone. In contrast, the ‘refracted plume’ hypothesis predicts more chaotic patterns quite different from the off‐chain pattern of Hawaiian volcanism.     

The PT-phase relations of an MgO-rich Hawaiian tholeiite: the compositions of primary Hawaiian tholeiites

The PT-phase relations of a Hawaiian tholeiite with 18.2% MgO has olivine–orthopyroxene multiple saturation at 20.5 kbar and 1,550°C. This pressure is less than the pressure at the lithosphere/asthenosphere transition, and it is suggested that tholeiites with this and lesser MgO contents are fractionated. Assuming a harzburgitic residuum it is shown that Hawaiian primary tholeiites contain about 23% MgO, and are generated at 36±5 kbar and 1,680±50°C. This pressure is equivalent to a depth of 112 km, which is consistent with the thickness of the lithosphere and thermal plume modeling. The minimal MgO content of primary Hawaiian tholeiites is suggested as 19% MgO.     

A dynamic model of the Hellenic Arc deduced from geophysical data

Combined gravity and seismic data from Greece and the adjacent areas have been used to explain the high seismicity and tectonic activity of this area. Computed 2-D gravity models revealed that below the Aegean region a large “plume” of hot upper-mantle material is rising, causing strong attenuation of the crust. The hot “plume” extends to the base of the lithosphere and has very probably been mobilized through compressional processes that forced the lithosphere to sink into the asthenosphere. The above model is supported by: high heat flow in the Aegean region; low velocity of the compressional waves of 7.7 km/sec for the upper mantle; lower density than normal extending to the base of the lithosphere; teleseismic P-wave travel-time residuals of the order of +2 sec for seismic events recorded at the Greek seismic stations; volcanics in the Aegean area with a chemical composition which can be explained by assuming an assimilation of oceanic crust by the upper mantle; deep seismicity (200 km) which has been interpreted by various authors as a Benioff zone.     

Deep structure of the Longling–Ruili fault underneath Ruili basin near the eastern Himalayan syntaxis: insights from magnetotelluric imaging

Magnetotelluric (MT) geophysical profiling has been applied to the determination of the deep structure of the Longling–Ruili fault (LRF), part of a convergent strike-slip fault system, underneath thick Caenozoic cover in Ruili basin in southwestern Yunnan, China. The recorded MT data have been inverted using a two-dimensional (2-D) nonlinear conjugate gradients scheme with a variety of smooth starting models, and the resulting models show common subsurface conductivity structures that are deemed geological significant. The models show the presence of a conductive (5–60 Ω m) cover sequence that is thickest (1–1.5 km) in the centre of the basin and rapidly pinches out towards the margins. A half-graben structure is interpreted for the Ruili basin. This is underlain by about 7–10 km thick upper crustal layer of high resistivity (>200–4000 Ω m) that is dissected by steep faults, which we interpret to flatten at depth and root into an underlying mid-crustal conductive layer at about 10 km depth. The mid-crustal layer does not appear to have been severely affected by faulting; we interpret it as a zone of partial melt or intracrustal detachment. The MT models suggest SE directed thrusting of basement rocks in the area. The Longling–Ruili fault is interpreted as a NW-dipping feature bounding one of the identified upper crustal fragments underneath Ruili city. We suggest that MT imaging is a potent tool for deep subsurface mapping in this terrain.     

中国北方大陆下的地幔热柱与岩石圈运动

本文首次提出中国北方大陆下存在一个地幔热柱的论证,并提出亚热柱(sub—plume)的新概念。热柱的中心与边缘部分的隙间熔浆分别为苦橄质玄武岩与碱性玄武岩。在渐新世到中新世约18.4 Ma内,北方大陆以3.26cm/a的速率向东南飘移了约600km,使日本海、渤海—华北平原等脱离热柱。导致晚第三纪日本海扩张的停止,渤海—华北平原等早第三纪火山喷发的突然中止。火山喷发期间,在热柱头部若干个亚热柱的形成,好似若干个“铆钉”穿入岩石圈,有效地阻止了岩石圈的飘移(这时的飘移速率只有0.05cm/a),我们把火山喷发称为固定岩石圈的“铆钉效应”。     

Pb,Sr, Nd,and Hf isotopic constraints on the origin of Hawaiian basalts and evidence for a unique mantle source

Pb, Sr, Nd, and Hf isotopic relationships among basalts from the Hawaiian Islands suggest that these basalts were derived from three sources; the oceanic lithosphere (Kea end member), the depleted asthenosphere (posterosional end member) and a deep-mantle plume (Koolau end member).Hawaiian tholeiites are derived within the lithosphere and the isotopic trends collectively defined by the tholeiite data are interpreted as a plume-lithosphere mixing trend. The isotopic characteristics of late-stage basalts are derived from the tholeiite source (lithosphere + plume) with additional input from the lithosphere, asthenosphere, or both. These basalts probably originate from near the asthenosphere-lithosphere boundary. Posterosional basalts are derived from the depleted asthenosphere, but their isotopic characteristics have been slightly modified by either the plume or the source of previously erupted volcanics. The isotopic data require little or no mixing of asthenospheric material into the plume during tholeiite production and thus are consistent with the concept of a rapidly ascending, fluid-rich plume. In addition to providing a source of heat, the plume may supply volatiles to both the sources of tholeiites and posterosional basalts.The isotopic characteristics of the Koolau (plume) component are unique among OIB sources. If undifferentiated or “primitive” mantle material still exists, then the radiogenic-isotope data for Koolau in combination with rare gas data for Hawaiian basalts in general suggest that the Hawaiian plume may be derived from such material. In any case, the Hawaiian Islands data, when compared to those of other OIB, serve to illustrate the isotopically diverse nature of mantle sources.     

Influence of the precollisional stage on subduction dynamics and the buried crust thermal state: Insights from numerical simulations

At continental subduction initiation, the continental crust buoyancy may induce, first, a convergence slowdown, and second, a compressive stress increase that could lead to the forearc lithosphere rupture. Both processes could influence the slab surface P–T conditions, favoring on one side crust partial melting or on the opposite the formation of ultra-high pressure/low temperature (UHP-LT) mineral. We quantify these two effects by performing numerical simulations of subduction. Water transfers are computed as a function of slab dehydration/overlying mantle hydration reactions, and a strength decrease is imposed for hydrated mantle rocks. The model starts with an old oceanic plate ( 100 Ma) subducting for 145.5 Myr with a 5 cm/yr convergence rate. The arc lithosphere is thermally thinned between 100 km and 310 km away from the trench, due to small-scale convection occuring in the water-saturated mantle wedge. We test the influence of convergence slowdown by carrying on subduction with a decreased convergence rate (≤ 2 cm/yr). Surprisingly, the subduction slowdown yields not only a strong slab warming at great depth (> 80 km), but also a significant cooling of the forearc lithosphere at shallower depth. The convergence slowdown increases the subducted crust temperature at 90 km depth to 705 ± 62 °C, depending on the convergence rate reduction, and might thus favor the oceanic crust partial melting in presence of water. For subduction velocities ≤ 1 cm/yr, slab breakoff is triggered 20–32 Myr after slowdown onset, due to a drastic slab thermal weakening in the vicinity of the interplate plane base. At last, the rupture of the weakened forearc is simulated by imposing in the thinnest part of the overlying lithosphere a dipping weakness plane. For convergence with rates ≥ 1 cm/yr, the thinned forearc first shortens, then starts subducting along the slab surface. The forearc lithosphere subduction stops the slab surface warming by hot asthenosphere corner flow, and decreases in a first stage the slab surface temperature to 630 ± 20 °C at 80 km depth, in agreement with P–T range inferred from natural records of UHP-LT metamorphism. The subducted crust temperature is further reduced to 405 ± 10 °C for the crust directly buried below the subducting forearc. Such a cold thermal state at great depth has never been sampled in collision zones, suggesting that forearc subduction might not be always required to explain UHP-LT metamorphsim.     

Global 1° × 1° thermal model TC1 for the continental lithosphere: Implications for lithosphere secular evolution

This paper reports a new 1° × 1° global thermal model for the continental lithosphere (TC1). Geotherms for continental terranes of different ages (> 3.6 Ga to present) constrained by reliable data on borehole heat flow measurements (Artemieva, I.M., Mooney, W.D. 2001. Thermal structure and evolution of Precambrian lithosphere: a global study. J. Geophys. Res 106, 16387–16414.), are statistically analyzed as a function of age and are used to estimate lithospheric temperatures in continental regions with no or low-quality heat flow data (ca. 60% of the continents). These data are supplemented by cratonic geotherms based on electromagnetic and xenolith data; the latter indicate the existence of Archean cratons with two characteristic thicknesses, ca. 200 and > 250 km. A map of tectono-thermal ages of lithospheric terranes complied for the continents on a 1° × 1° grid and combined with the statistical age relationship of continental geotherms (z = 0.04  t + 93.6, where z is lithospheric thermal thickness in km and t is age in Ma) formed the basis for a new global thermal model of the continental lithosphere (TC1). The TC1 model is presented by a set of maps, which show significant thermal heterogeneity within continental upper mantle, with the strongest lateral temperature variations (as large as 800 °C) in the shallow mantle. A map of the depth to a 550 °C isotherm (Curie isotherm for magnetite) in continental upper mantle is presented as a proxy to the thickness of the magnetic crust; the same map provides a rough estimate of elastic thickness of old (> 200 Ma) continental lithosphere, in which flexural rigidity is dominated by olivine rheology of the mantle.Statistical analysis of continental geotherms reveals that thick (> 250 km) lithosphere is restricted solely to young Archean terranes (3.0–2.6 Ga), while in old Archean cratons (3.6–3.0 Ga) lithospheric roots do not extend deeper than 200–220 km. It is proposed that the former were formed by tectonic stacking and underplating during paleocollision of continental nuclei; it is likely that such exceptionally thick lithospheric roots have a limited lateral extent and are restricted to paleoterrane boundaries. This conclusion is supported by an analysis of the growth rate of the lithosphere since the Archean, which does not reveal a peak in lithospheric volume at 2.7–2.6 Ga as expected from growth curves for juvenile crust.A pronounced peak in the rate of lithospheric growth (10–18 km³/year) at 2.1–1.7 Ga (as compared to 5–8 km³/year in the Archean) well correlates with a peak in the growth of juvenile crust and with a consequent global extraction of massif-type anorthosites. It is proposed that large-scale variations in lithospheric thickness at cratonic margins and at paleoterrane boundaries controlled anorogenic magmatism. In particular, mid-Proterozoic anorogenic magmatism at the cratonic margins was caused by edge-driven convection triggered by a fast growth of the lithospheric mantle at 2.1–1.7 Ga. Belts of anorogenic magmatism within cratonic interiors can be caused by a deflection of mantle heat by a locally thickened lithosphere at paleosutures and, thus, can be surface manifestations of exceptionally thick lithospheric roots. The present volume of continental lithosphere as estimated from the new global map of lithospheric thermal thickness is 27.8 (± 7.0) × 10⁹ km³ (excluding submerged terranes with continental crust); preserved continental crust comprises ca. 7.7 × 10⁹ km³. About 50% of the present continental lithosphere existed by 1.8 Ga.     

The lithospheric transition between the Delamerian and Lachlan orogens in western Victoria: new insights from 3D magnetotelluric imaging

The magnetotelluric (MT) method was used to image the crust and upper mantle beneath the Delamerian and Lachlan orogens in western Victoria, Australia. During the Cambrian time period, this region changed from being the extended passive margin of Proterozoic Australia into an Andean-style convergent margin that progressively began to accrete younger oceanic terranes. Several broadband MT transects, which were collected in stages along coincident deep (full crust imaging) seismic reflection lines, have now been combined to create a continuous 500 km east–west transect over the Delamerian–Lachlan transition region in the Stawell Zone. We present the electrical resistivity structure of the lithosphere using both 3D and 2D inversion methods. Additionally, 1D inversions of long-period AusLAMP (Australian Lithospheric Architecture Magnetotelluric Project) MT data on a 55 km regionally spaced grid were used to provide starting constraints for the 3D inversion of the 2D profile. The Delamerian to Lachlan Orogen transition region coincides with the Mortlake Discontinuity, which marks an isotopic discontinuity in Cenozoic basalts, with higher strontium isotope enrichment ratios in the Lachlan Orogen relative to the Delamerian Orogen. Phase tensor ellipses of the MT data reveal a distinct change in electrical resistivity structure near the location of the Mortlake Discontinuity, and results of 3D and 2D inversions along the MT profile image a more conductive lower crust and upper mantle beneath the Lachlan Orogen than the Delamerian Orogen. Increased conductivity is commonly ascribed to mantle enrichment and thus supports the notion that the isotope enrichment of the Cenozoic basalts at least partially reflects an enriched mantle source rather than crustal contamination. Fault slivers of the lower crust from the more conductive Lachlan region expose Cambrian boninites and island arc andesites indicative of subduction, a process that can enrich the mantle isotopically, and also electrically, by introducing carbon (graphite) and water (hydrogen).     

Paleomagnetic modeling of seamounts near the Hawaiian–Emperor bend

The Hawaiian–Emperor Seamount chain records the motion of the Pacific Plate relative to the Hawaiian mantle hotspot for 80 m.y. A notable feature of the chain is the pronounced bend at its middle. This bend had been widely credited to a change in plate motion, but recent research suggests a change in hotspot motion as an alternative. Existing paleomagnetic data from the Emperor Chain suggest that the hotspot moved south during the Late Cretaceous and Early Tertiary, but reached its current latitude by the age of the bend. Thus, data from area of the bend are important for understanding changes in plume latitude. In this study, we analyze the magnetic anomalies of five seamounts (Annei, Daikakuji-W, Daikakuji- E, Abbott, and Colahan) in the region of the bend. These particular seamounts were chosen because they have been recently surveyed to collect multibeam bathymetry and magnetic data positioned with GPS navigation. Inversions of the magnetic and bathymetric data were performed to determine the mean magnetization of each seamount and from these results, paleomagnetic poles and paleolatitudes were calculated. Three of the five seamounts have reversed magnetic polarities (two are normal) and four contain a small volume of magnetic polarity opposite to the main body, consistent with formation during the Early Cenozoic, a time of geomagnetic field reversals. Although magnetization inhomogene ties can degrade the accuracy of paleomagnetic poles calculated from such models, the seamounts give results consistent with one another and with other Pacific paleomagnetic data of approximately the same age. Seamount paleolatitudes range from 13.7 to 23.7, with an average of 19.4 ± 7.4 (2σ). These values are indistinguishable from the present-day paleolatitude of the Hawaiian hotspot. Together with other paleomagnetic and geologic evidence, these data imply that the Hawaiian hotspot has moved little in latitude during the past 45 m.y.     

Crustal accretion and dynamic feedback on mantle melting of a ridge centred plume: The Iceland case

A major consequence of the interaction of a plume with an oceanic ridge is the enhanced melt production and associated crust generation. In the case of Iceland crustal thickness as large as 20 to 40 km has been reported. Crustal seismic velocities are high, and have to be explained by thermal or chemical effects. In the first part of the paper we address the question whether extraction of melt out of the plume beneath a slowly spreading ridge and deposition of extracted basalt volumes at the surface produces a dynamic feedback mechanism on mantle melting. To study this question we solve the convection equations for a ridge centred plume with non-Newtonian rheology including melting, melt extraction associated with deposition of cold crust at the surface of the model, and using a simplified approach for compaction. The assumption of cold crust is justified if the thickness of each deposited basaltic layer is less than roughly 1 km. Depending on the buoyancy flux of the plume, crustal thicknesses between 10 and 40 km are modelled, showing characteristic dipping structures resembling the rift-ward dipping basaltic layers of East- and Western Iceland. Comparing the resulting crustal thickness and magma generation rate with models in which the dynamic effect of crust deposition has been suppressed indicates, that melt generation beneath a slowly spreading ridge is considerably damped by the dynamic feedback mechanism if the plume buoyancy flux exceeds 400 to 600 kg/s. Based on the observed crustal thickness of Iceland our models predict a plume buoyancy flux of 1140 kg/s.In the second part we study the accretion of the Icelandic crust by a thermo-mechanical model in more detail based on the Navier–Stokes-, the heat transport and the mass conservation equations including volumetric sources. Hot (1200 °C) molten crustal material is injected into the newly forming crust with a constant rate at different crustal source regions: a) deep, widespread emplacement of dykes and sills including crustal underplating, b) magma chambers at shallow to mid-crustal level, and c) surface extrusions and intrusions in fissure swarms at shallow depth connected to volcanic centres. We identify the material from the different source regions by a marker approach. Varying the relative dominance of these source regions, characteristic crustal structures evolve, showing shallow dipping upper crustal layers with dip angles between 10 and 15°. The thermal structure of the crust varies between cold crust (shallow-source region dominating) and hot crust (deep-source region dominating). We use observations of maximum depth of seismicity to constrain the depth of the 650 °C isotherm and seismological inferences on the lower crust to constrain temperatures in that region. The best agreement with our models is achieved for crust formation dominated by deep dykes and underplating with a considerable influence of magma chamber accretion.     

Petrogenesis of Hawaiian tholeiites: 2, aspects of dynamic melt segregation

To Hawaiian magma genesis, dynamic melt segregation offers a potential resolution of conflict arising between trace-element evidence and phase-equilibria evidence, for deep garnet-present melting versus shallow garnet-absent melting. In this study comprehensive dynamic melting models, which incorporate phase-equilibria constraints and variable partition coefficients, have been applied in efforts to simulate decompression melting of a mantle plume. These models specifically endeavour to reproduce Hawaiian REE (rare-earth-element) patterns from a peridotitic upper mantle source with chondritic relative abundances of middle and HREE (heavy REE). If the flow of both melt and solid mantle is vertical through the partially molten source region, and melting proceeds beyond the stability limit of garnet in peridotite, dynamic melting processes are unable to produce the fractionated REE patterns of Hawaiian tholeiites. Instead, three-dimensional dynamic melting modles need to be invoked, in which lateral migration of the melt relative to the residual matrix also takes place. This enables the derivation of small garnet-equilibrated melt fractions from a larger source volume than that supplying more extensive melt fractions from shallower garnet-absent levels of melting (i.e melting shapes with a mean degree of melting smaller than the maximum extent of melting). This can be achieved by either drawing small-degree melt fractions, formed in the presence of garnet at the plume peripheries, toward the plume centre, or by advecting the mantle residue away from the plume centre as it ascends. Fluid dynamic theory supports a plume model incorporating the latter, with melt flow occurring vertically through a matrix flow which is deflected by the lithosphere and diverges away from the plume centre. In this framework, the generation of melting shapes dominated by small-degree garnetpresent melt fractions, requires a decrease in the rate of melting with progressive melting and height along melt-flow paths within the plume. This is consistent with a decrease in vertical velocity of the matrix (and thus decompression melting rate) upwards through the plume and, with diminishing melting rates upon exhaustion of garnet and clinopyroxene as melting progresses. Providing melt segregation occurs by percolation, equilibrium between the segregating melt and residual peridotite matrix may be maintained throughout the plume. In this way, primary melts extracted from the Hawaiian plume have their bulk compositions determined by phase equilibrium with the extensively melted matrix residue (harzburgite) at the plume top and shallowest level of melting (2.0 GPa), and their incompatible-trace-element characteristics determined by smaller-degree melt fractions derived from deeper, garnet-present levels of melting (3.0 GPa). Simple unidimensional models for melt segregation by percolation or via channels are shown to produce incompatible-trace-element abundances and ratios which are similar to those generated by equivalent degrees of batch melting. Moreover, contrary to a common belief held for dynamic melting, the enrichment of more-incompatible elements over less-incompatible elements is not always greater than that produced by an equivalent amount of batch melting.