S. A. Silantyev1, O. G. Bogdanovskii2, P. I. Fedorov3, S. F. Karpenko1, and Yu. A. Kostitsyn1
1Vernadsky Institute of Geochemistry and Analytical Chemistry, Russian Academy of Sciences, Moscow, Russia
2Max-Plank-Institut für Chemie, Postfach 3060, D-55020 Mainz, Germany
3Geological Institute, Russian Academy of Sciences, Moscow, Russia
Figure 1 |
The De Long Islands were named after George Washington De Long, a US Navy lieutenant, who headed in 1879 an American expedition aimed to achieve the North Pole. De Long took his yacht "Jeannette'' through the Bering Strait and headed northwest. The "Jeannette'' was enclosed by the ice, froze in, and drifted in this state for 500 km northwest during 20 months. In the course of this drift, which ended with the crushing of the vessel and its sinking, the islands of Jeannette, Henrietta, and Bennett were discovered. Most participants of the expedition, including De Long himself, died of exposure and starvation when trying to reach Russian settlements in the Siberian Arctic shore near the Lena mouth. In 1902-1915, the islands and surrounding waters were visited by a series of Russian expeditions. During one of them, headed by E. V. Tol', the first rock samples were collected in Bennett Island, but Tol' himself and all his men vanished without a trace, and their rock collection was found later in Bennett Island. The expedition of two Russian icebreakers "Taymyr'' and "Vaigach'' (the chief of the 1913-1915 expeditions was B. A. Vil'kitskii, assistant to the chief A. N. Zhokhov), which was launched with the aim of exploring the prospects of using the Northern Sea Route, discovered the islands of Vil'kitskii (1913) and Zhokhov (1915). Later, fieldwork was sporadically conducted on the islands of the archipelago by AANII.
The De Long Islands lie practically exactly on the continuation of the trend axis of the Gakkel mid-oceanic ridge into the passive continental margin of the Laptev Sea shelf, a fact that led Ya. Ya. Gakkel to hypothesize in 1957 that the superstructure of the Mid-Atlantic Ridge (MAR), including the Gackcel Ridge, extends across the North Pole and De Long Islands to the basin of the Indigirka River [Gakkel, 1957]. Later geophysical data on this area in the Arctic basin were interpreted [Naryshkin, 1987] as indicating that the islands of Zhokhov and Vil'kitskii belong to an extensive submarine plateau, which reaches the continental slope of the Eurasian basin. Mafic magmatic rocks at De Long Islands were dated [Vol'nov and Sorokov, 1961; Vol'nov et al., 1970] as Cretaceous and younger. Newly obtained data on the composition and age of volcanics from Bennettt Island were reported in [Fedorov et al., 2002] and used in this paper.
The most exhaustive and systematic information on the composition and age of magmatic rocks and related mantle and crustal xenoliths from the central part of the De Long Islands were presented in a series of publications that synthesized the results of fieldwork on Zhokhov and Vil'kitskii islands conducted in 1986 and 1988 within the scope of the joint high-latitude expedition of Shirshov Institute of Oceanology and Vernadsky Institute of Geochemistry and Analytical Chemistry of the Russian Academy of Sciences [Bogdanovskii et al., 1992, 1993; Savostin et al., 1988; Silantyev et al., 1991, 2002]. The newly obtained data reported in these papers were supplemented with the latest information on magmatic rocks of Bennett Island and used as the basis of this research.
The main cause why we returned to the problem of the geodynamic environment of magmatism at De Long Islands was the obvious breakthrough in the currently adopted concepts of the geology of the Eurasian sector of the Arctic Ocean related to exploring the crest zone of the Gakkel Ridge between 5o W and 85o E conducted by the AMORE international expedition in 2001 aboard the well-equipped icebreaker-class vessels "Polarstern'' (Alfred-Wegener-Institut, Bremenhaven, Germany) and "Healy'' (USCGS, Seattle, United States). The results obtained by this expedition provided vast volumes of new geophysical and bathymetric information on the structure of the Gakkel Ridge and, what is particularly important in the context of our research, a large collection of mantle peridotites and basalts from throughout the axial zone of the ridge.
This paper also presents the very first data on the Sr and Nd isotopic composition of basalts from the rock collection sampled at lava flows on Bennett Island, one of the largest islands of the archipelago. When compared with preexisting data on Zhokhov and Vil'kitskii islands, this information makes it possible to reconstruct the geochemical nature of the mantle sources of basaltic magmatism in the continental shelf of the Laptev Sea immediately southeast of its intersection with the southern termination of the Gakkel Ridge. This comparative analysis of the isotopic characteristics of magmatism on the islands of the archipelago is instrumental in reproducing the major stages of intraplate volcanism in the passive continental margin during the gradual propagation of the ultraslow-spreading ridge. It should be mentioned that De Long Islands are among the world's few structures suitable for examining the geodynamic environment of interaction between a so-called "propagating ridge'' and a passive continental margin (such as the Red Sea Rift).
We considered it necessary to amend our earlier isotopic data published elsewhere. In our previous research conducted in the 1970s in compliance with the method developed by G. Wasserburg et al. at the California Institute of Technology, the 150Nd/ 142Nd isotopic ratio was normalized to 0.209647, and, thus, the 143Nd/ 144Nd ratio for CHUR was assumed equal to 0.511847. The data obtained in the course of this research by modern methods, currently employed at most isotopic laboratories around the world, were based on normalizing the 146Nd/144Nd ratio to 0.7219, which implies the value of the 143Nd/144Nd ratio for CHUR equal to 0.512638. Hence, to compare our earlier data [Bogdanovskii et al., 1992, 1993] with those obtained recently on volcanic rocks from Bennett Island, all of our earlier 143Nd/144Nd ratios were corrected using a coefficient equal to 1.00154. It is pertinent to mention that none of our 147Sm/144Nd ratios has changed.
Information available on the geology of the De Long Islands is scarce, almost all of it presented in the classic multivolume publication [Geology..., 1970].
According to the aforementioned publication, Bennett Island is composed of Cretaceous [Vol'nov and Sorokin, 1961] or Miocene [Geology..., 1984] basalt flows, which rest on compositionally variegated sedimentary rocks. The sedimentary complex consists of Cambrian mudstones with trilobite fossils and Ordovician mudstones and sandstones. The thickness of the olivine basalt flows amounts to approximately 500 m. It has been demonstrated [Fedorov et al., 2002] that the alkaline basalts composing stratified lava sheets on Bennett Island have an age of 106 Ma, while the magnesian basalts composing volcanic cones and lava flows were dated (K-Ar) at 109-124 Ma.
Earlier researchers suggested that both islands are made up of olivine and nepheline basalts [Backlund, 1920; Geology..., 1970]. The results of our fieldwork in 1986 and 1988 suggest that Zhokhov Island completely consists of basaltoid lava flows (mostly of olivine-phyric varieties). As can be seen in the highest (50 m) exposure in the coastal cliffs in the southeastern portion of Zhokhov Island, there are at least five lava sheets, whose tops are marked with cinderlike rocks of dark gray and reddish brown colors, often with traces of flow. Each of the flows is usually clearly stratified from top to bottom in the following manner: (1) highly porous olivine-phyric basalt, (2) basalt of intermediate porosity, (3) massive basalt with rare olivine crystals, and (4) variolitic basalt. The central, most elevated part of Zhokov Island, which is crowned with a small rocklet (called Devil's Finger), and the area north of it consist of black, reddish, and bright orange volcanic rocks with evidence of explosive genesis (abundant volcanic cinder and fragments of volcanic bombs). The low dome-shaped hillocks in the central part of the island seem to have been eruption centers from which lava flows were outpoured. All varieties of the Zhokov basaltoids contain small (no more than 2 3 5 cm) apple-green xenoliths, which are clearly discernible in the dark gray and black host rocks. These xenoliths resemble mantle xenoliths of lherzolites widespread in areas of intraplate magmatism. Lava flows on the southwestern shore of the island also occasionally contain xenoliths of mafic magmatic rocks and crustal quartzites.
When disembarking at Vil'kitskii Island in May of 1988, the authors of this paper managed to collect samples at an extensive exposure in the southern cliffy shore of this small (10 3 km) island. The central part of the exposure (that extended for about 4.5 km from west to east) was composed of massive basaltoids with very small xenoliths identical to those in the Zhokhov rocks. From east to west, the body of massive basaltoids was surrounded by highly porous lava flows resembling the volcanic rocks in the central part of Zhokhov Island. These porous lavas also contained rare very small xenoliths of apple-green color.
The data reported in [Geology..., 1970] indicate that the stratigraphic sequence of Henrietta Island comprises three major units: (1) lower, consisting of quartzites and clayey shales (150-200 m); (2) middle, made up of tuffites and diabases (700-900 m); and (3) upper, composed of conglomerates (200-300 m). A similar stratigraphy is also typical of Henrietta Island. The sedimentary complexes of Jeannette and Henrietta islands are of Cretaceous age. A similar age was determined for basalts with pyroxene (augite) phenocrysts that intrude the sedimentary sequence of Henrietta Island.
The contents of major and trace elements were determined on a PW-1600 (Philips) XRF spectrometer (analyst T. V. Romashova) at the Central Analytical Laboratory of Vernadsky Institute of Geochemistry and Analytical Chemistry, Russian Academy of Sciences. REE were analyzed in the rocks by neutron activation at the same analytical center. The Nd, Sr, S, C, and O isotopic ratios and the K-Ar age of the alkaline olivine basalts, limburgites, and crustal and mantle xenoliths from Zhokhov and Vil'kitskii islands were discussed in [Bogdanovskii et al., 1992, 1993]. We used data on the contents of major incompatible elements in basalts from Bennett Island and information on the mineralogy of these rocks compiled from [Fedorov et al., 2002].
Minerals in rocks from Zhokhov and Vil'kitskii islands were analyzed by N. N. Kononkova on a CAMEBAX electron microprobe at Vernadsky Institute. The standards were natural and synthetic minerals. The spot analyses were conducted at an accelerating voltage of 15-20 kV and a beam current of 35 mA.
The Sr and Nd isotopic compositions of volcanic rocks from Bennett Island were measured on a TRITON TI multicollector mass spectrometer at the Laboratory of Isotopic Geochemistry of Vernadsky Institute.
Figure 2 |
Figure 3 |
Figure 4 |
Figure 5 |
Figure 6 |
Figure 7 |
Figure 8 |
Figure 9 |
Figure 10 |
The central, steepest part of the cliffy exposure in the southern shore of Vil'kitskii Island is made up of massive dark gray volcanics cut by thin shattering zones, in which rocks acquire thin-platy parting. The volcanic pile ubiquitously contains small olivine phenocrysts and spinel lherzolite xenoliths. The volcanic rocks of Vil'kitskii Island are both holocrystalline and predominantly vitreous, with a devitrified matrix. The textures of these rocks and their mineralogy are very close to those of the volcanics described at Devil's Finger Hill on Zhokov Island (Figure 6b). In turn, the mineralogy of xenoliths from Vil'kitskii Island is completely identical to that of mantle xenoliths from Zhokhov Island.
Olivine is the most common mineral in magmatic rocks of the De Long Islands. According to its composition, this olivine was classified into the following groups: (1) Olivine in spinel lherzolite xenoliths, which has a high Mg#, low CaO content, and high NiO. Analogous characteristics are typical of olivine xenocrysts from volcanic rocks of Zhokhov and Vil'kitskii islands. As can be seen in Figure 3, the alkaline basaltoids of Bennett Island contain no xenocryst olivine. (2) Olivine phenocrysts contained in subordinate amounts in lava flows on Zhokhov and Vil'kitskii islands and occurring as major minerals in the basaltoids of Bennett Island. They are high in CaO, have a moderate Mg#, and are low in NiO. This group also includes olivine from the dolerite xenolith found on Zhokhov Island. The maximum CaO and minimum MgO contents are typical of olivine from basaltoids in Bennett Island. Within this compositional field, olivine from the magnesian silica-undersaturated alkaline basalts is more magnesian than this mineral from the alkaline (according to the systematics in [Fedorov et al., 2002]) basalts. (3) Some olivine phenocrysts in the lava complex of Zhokhov Island have compositions intermediate between those of the olivines of the two aforementioned groups (Figure 3). Obviously, these intermediate compositions were caused by the recrystallization of the mantle material during its interaction with the magmatic melt.
Figure 11 |
As the clinopyroxene, orthopyroxene phenocrysts in the lava complex and mantle xenoliths belong to two compositional groups with different Al2O3 contents: (1) 3.17-3.54 wt % and (2) 5.19-6.46 wt %. The high-Al compositions are also characterized by higher Na2O contents (Table 1).
Inasmuch as clinopyroxene phenocrysts in alkaline basalts typically have high Ti contents, there are good grounds to believe that the lavas of Zhokhov and Vil'kitskii islands contain only xenogenic clinopyroxene. Conversely, this mineral occurs in the basaltoids of Bennett Island exclusively as phenocrysts. Since orthopyroxene is atypical of the phenocryst assemblages of alkaline basalts and strongly undersaturated volcanic rocks [Coombs and Wilkinson, 1969], this phase in volcanic rocks from the central part of the De Long Islands can be attributed, along with the low-Ti clinopyroxene, to the xenocryst assemblage. It is worth mentioning that the alkaline basaltoids of Bennett Island contain no orthopyroxene at all. The fact that the spinel lherzolite xenoliths contain pyroxene of two populations having significantly different compositions suggests that there could be high-temperature interaction between the xenoliths and the melt entraining them.
Figure 12 |
The data presented above on the composition of rock-forming minerals in volcanic rocks and xenoliths from Zhokhov and Vil'kitskii islands provide evidence that olivine phenocrysts in the lava complex of these islands is of predominantly xenogenic (mantle) nature, whereas all pyroxene phenocrysts are xenogenic. At the same time, the mantle xenoliths contain both phases of undoubtedly mantle provenance and the assemblage of olivine and pyroxenes produced during the partial recrystallization of these xenoliths during their interaction with the melt. It is also reasonable to hypothesize that the basaltoids of Bennett Island provide no mineralogical indications that mantle material was entrapped by the parental melts.
Available data make it possible to evaluate the temperatures and pressure of the stability of the spinel lherzolites and xenocryst assemblage in volcanic rocks from the central part of the De Long Islands. For this purpose we used the compositions of coexisting clino- and orthopyroxene, olivine and clinopyroxene, and olivine and spinel. The pyroxene geothermometer making use of covariations in the Ca concentration and Mg# in pyroxene [Lindsley and Anderson, 1983] enables, according to [Vaganov and Sokolov, 1988], reliable temperature evaluations only for low and mildly aluminous compositions. More realistic estimates can be obtained by the two-pyroxene geothermobarometer [Mercier, 1980], which was proposed for pyroxenes with broadly varying Al# from lherzolites of the spinel depth facies. This method also allows for pressure evaluating, although, strictly speaking, these estimates are not completely independent, inasmuch as all the calculations are conducted with the same parameters from which the temperature is calculated. The P-T parameters calculated, according to [Mercier, 1980], for the crystallization of mantle xenoliths from Zhokhov Island are T = 1120o-970o C, P = 25-12.5 kbar. The olivine-clinopyroxene geothermometer [Mori and Green, 1978] enabled us to estimate the maximum crystallization temperature of the olivine-clinopyroxene assemblage in the mantle xenoliths, which was 1211o C. In addition to the geothermometers mentioned above, we also evaluated the temperatures by the orthopyroxene geothermometer [Sachtleben and Seck, 1981], olivine-spinel geothermometer [O'Neil and Wall, 1987], and the two-pyroxene geothermometer [Wells, 1977].
The P-T parameters calculated by the geothermometers and geobarometers listed above demonstrate that there were three major temperature ranges in which the equilibrium phase assemblages of the mantle xenoliths were formed: ~1100o-1200oC, ~1000o C, and ~760-810o C. The lowest temperatures were obtained by the olivine-spinel geothermometer, because the Ca-Mg exchange between pyroxenes in spinel peridotites is blocked at higher temperatures than those at which Fe-Mg exchange between olivine and spinel stops. Thus, for spinel lherzolite xenoliths from the De Long Islands, the olivine-spinel geothermometer yielded temperature values corresponding to the latest exhumation stages of the mantle material.
Figure 13 |
Figure 14 |
Figure 15 |
Figure 16 |
Figure 17 |
Figure 18 |
Figure 19 |
It can be concluded that the geochemical features of volcanic rocks from Bennett, Zhokhov, and Vil'kitskii islands make them similar to volcanics in within-plate magmatic provinces, including the basalts of oceanic islands (OIB) and continental rifts. It is pertinent to emphasize that within-plate volcanic products in the central part of the De Long Islands also show geochemical similarities with the young alkaline basalts of volcanic centers in the continental fringing of the Laptev Sea shelf southeast of the sea: Balagan-Tas Volcano in the Chersky Range, Moma Rift, and the Quaternary (?) dike complex at the Viliga River in the Sea of Okhotsk area, which contain spinel lherzolite xenoliths [Grachev, 1999]. The data presented in Figures 12-19 do not rule out the possibility that the alkaline basalts of Bennett Island were produced by a parental melt whose composition was similar to that of the parental melt of the olivine alkaline basalts in Zhokhov Island. The dolerite of the crustal xenolith form Zhokhov Island has many compositional parameters (the Sr and K2O concentrations and the (La/Sm)cn ratio] close to the analogous values for enriched MORB basaltoids, although the former rock is poorer in TiO2 and Na2O. Earlier, analysis of the distribution of major and incompatible elements in basalts from Bennett Island led to the conclusion that these rocks could be close to continental flood basalts or rift-related alkaline basalts in continental margins, and the magnesian alkaline basalts are close to the rocks of mildly alkaline series in continental rifts or OIB [Fedorov et al., 2002].
The geochemistry of spinel lherzolites from the De Long Islands make these rocks similar to peridotites, i.e., the material of the primitive (undepleted) mantle. This follows from the MgO/SiO2 and Al2O3/SiO2 ratios of these rocks: 1.19-1.25 and 0.07-0.08, respectively. Spinel lherzolites similar to those found on Zhokhov and Vil'kitskii islands are the most ubiquitous rocks of mantle xenoliths in areas of within-plate magmatism in both continents and oceanic islands, such as Victoria in Australia and Hawaii [Frey and Green, 1974; Irving, 1980; Jackson and Wright, 1970].
The K-Ar ages of the olivine alkaline basalts, limburgites, and dolerite xenolith from Zhokhov and Vil'kitskii islands were published in [Bogdanovskii et al., 1992] and are presented in Table 5. Judging from available data, limburgite volcanism on Vil'kitskii Island (0.89-0.4 Ma) can be correlated with the late evolutionary stages of the olivine alkaline basalt association on Zhokhov Island (6.1-0.4 Ma) but was significantly briefer. The K-Ar ages of the rocks composing the volcanic complex of the Bennett Island were reported in [Fedorov et al., 2002] and are 106 Ma for the alkaline basalts composing stratified lava sheets and 109-124 Ma for the magnesian alkaline basalts of the lava flows and cones.
The leucocratic and melanocratic fractions of dolerite from the crustal xenolith suite were dated by the K-Ar method at 152 5 and 100 3 Ma, respectively [Bogdanovskii et al., 1992]. Considering that the sample shows traces of secondary alterations and that plagioclase (which is abundant in the leucocratic fraction) is prone to entrap excess radiogenic Ar (whose content is higher in the leucocratic fraction), preference should be given to the younger age value, which was obtained for melanocratic minerals. It is, however, impossible to unambiguously correlate this value with any geologic event without employing additional information. Provisionally we assumed that the age of this dolerite is 100 3 Ma, but it cannot be ruled out that part of the radiogenic Ar was lost during reheating, so that the actual age of the rock should be much older.
There are, thus, good reasons to believe that volcanic activity at the De Long Islands was prone to become systematically younger from the offshore margin of the Laptev Sea shelf (Bennett Island--Early-Late Cretaceous boundary) toward its inner parts (Zhokhov Island) and farther southeastward (Vil'kitskii Island--Quaternary, Pleistocene). Proceeding from the assumption that the dolerite xenolith has isotopic and geochemical characteristics classifying it as a fragment of the crustal material of the Laptev Sea shelf, it should be assumed that this material has an early Cretaceous age.
Figure 20 |
Figure 21 |
The isotopic characteristics of spinel lherzolites from the De Long Islands and their host within-plate volcanic rocks allowed us to conclude that the young (6.1-0.4 Ma) volcanic rocks of this sector of the Arctic basin carry xenoliths of ancient (1.11 Ga) mantle material. Available data provide grounds to believe that mantle xenoliths from Zhokhov and Vil'kitskii islands are fragments of mantle rocks whose Nd and Sr isotopic composition corresponded to the depleted mantle. This material was not involved in the melting processes over at least the past 1.11 Ga.
The close spatial association of alkaline basalts and limburgites in the central part of the De Long Islands (on Zhokhov Island) suggests their genetic similarity and a common magmatic history. Data on the Sr and Nd isotopic composition of the volcanic rocks (Table 6) make it possible to reconstruct more confidently both the nature of the mantle sources of magmatism that gave rise to these rocks and the evolution of the parental melts. These rocks contain less radiogenic Nd and more radiogenic Sr than in the source of MORB, with the most radiogenic Sr detected in the dolerite of the crustal xenolith from Zhokhov Island.
The data presented above indicate that the alkaline basalts and limburgites of Bennett, Zhokhov, and Vil'kitskii islands have similar Sr and Nd isotopic composition and, thus, can be regarded as the derivatives of a single magmatic melt, which was derived by the melting of a mantle source (or sources) having similar Sr and Nd isotopic signatures, with the 87Sr/86Sr ratio varying within 0.0006 and e Nd within two units. However, when the Sr isotopic compositions of these rocks are examined more closely, they show small but systematic differences. These differences can be accounted for by (1) the insignificant isotopic heterogeneity of the source, (2) the contamination of the melt with the material of the oceanic crust, or (3) the interaction of the magmatic source with seawater. Below these three possibilities will be considered in the reverse order.
All of the rocks used in our research were very fresh, without any mineralogical or geochemical traces of interaction with seawater. An independent argument for the absence of this interaction is furnished by the identity of the Nd and Sr isotopic composition of the acid leachates and residues after the preliminary leaching of the samples.
Figure 22 |
As is seen in Figure 22a, the Sr and Nd isotopic ratios in the limburgites and alkaline olivine basalts are close to those in the enriched (plume-related) sources of MORB or OIB. The genesis of these enriched (with respect to the DMM source) rocks was a matter of heated discussions over the past decades and is believed to be somehow related to the recycling of the crustal material: via either the direct recycling of the isotopically anomalous material or the recycling of chemically enriched material that eventually gives rise to isotopic anomalies in the mantle. A comprehensive review on this problem was published in [Hofmann, 2003].
Figure 22a demonstrates mixing lines for a MORB melt from the hypothetical DMM source [Zindler and Hart, 1986] and a potentially possible crustal contaminant (represented by sample DLK-7). The Nd and Sr isotopic ratios in the limburgites and alkaline olivine basalts are in apparent conflict with the hypothesis of mixing of melts from the DMM source and oceanic crustal material.
The genesis of the isotopic ratios observed in the limburgites and alkaline basalts of the De Long Islands can be understood from Figure 21. The present variations in the Nd isotopic ratios in spinel lherzolite nodules are quite broad: from +1.7 to +14.7 e Nd units. This isotopic heterogeneity was caused by the in-situ enrichment of 143Nd in originally (at 1.11 Ga) isotopically homogeneous material as a natural consequence of the chemical heterogeneity of this material, including an uneven distribution of the Sm/Nd ratio. If the source of the melts that produced the De Long volcanics was characterized by a chemical heterogeneity similar to that in xenoliths of the analyzed samples (DLK-1 and DLK-3), then the occurrence of local domains enriched in lithophile elements, including LREE, for a long time (during which the Sm/Nd ratio decreased) should have inevitably brought about the development of isotopically enriched characteristics of these rocks. The melting of these isotopically and chemically anomalous zones could have simultaneously produce chemically enriched melts with lower Nd isotopic ratios.
Figure 23 |
Analyzing the spatial distribution of the geochemical types of volcanic complexes of different ages within the sector of the Arctic basin discussed in this paper enabled us to reproduce the geochemical nature of the upper mantle beneath the central part of the Laptev Sea continental shelf and the nature of the deep-seated mantle sources that produced the within-plate magmatic melt over the past 120 m.y.
Isotopic data testify that the spinel lherzolites entrained as xenoliths by magmas in the De Long Islands could not be captured in the mantle sources whose melting gave rise to the Mesozoic and Cenozoic magmatism on Bennett, Zhokhov, and Vil'kitskii islands. Available data indicate that the material represented by these xenoliths was not involved in melting processes over at least the past 1.11 Ga. At 1.11 Ga, the mantle material whose fragments are the spinel lherzolites was compositionally, but not isotopically, heterogeneous. However, this initial compositional heterogeneity triggered the development of isotopic heterogeneity in the same mantle material (see above: samples DLK-1 and DLK-3).
Mantle xenoliths with isotopic characteristics similar to those of the spinel lherzolites from the De Long Islands are known in the eastern and central parts of the Asian continent: these are garnet peridotites in the Mir kimberlite pipe in Yakutia and garnet and spinel peridotites from the Vitim Highlands and Mongolia [Ionov, 1988; Ionov and Jagoutz, 1988; Kovalenko et al., 1990]. It is possible that the De Long Islands mark the northernmost limit of the Asian mantle anomaly [Ionov and Jagoutz, 1988; Zhuravlev et al., 1991], which is located at depth of 45-75 km and has not been involved in melting processes since 1 Ga.
As was demonstrated above, petrological evidence indicates that the weakly depleted mantle source material represented by the spinel lherzolites is localized beneath the central part of the De Long Islands at a depth of approximately 60 km, and, thus, the sources of the parental melts of the alkaline basalts and limburgites should have been situated at even greater depths. A. E. Ringwood believed [Ringwood, 1975] that the melts parental for nephelinite (in our case, for the limburgite) series are derived at greater depths than the melts to which alkaline olivine basalts are related. At the same time, the two within-plate magmatic series were, perhaps, generated by the partial melting of a single mantle source, with limburgite-nephelinite and olivine alkaline basaltic melts produced at low and high degrees of melting, respectively. This scenario for the origin of the two major types of magmatic series in the De Long Islands seems to be the most plausible, because the data presented above demonstrate similarities between the isotopic-geochemical characteristics of the limburgites and alkaline olivine basalts. The enrichment of the limburgites in incompatible elements (for example, Sr and K) relative to the alkaline olivine basalts could hardly be caused by the fractionation of a common parental melt, because this scenario is at variance with the much more stronger LREE enrichment in the limburgites than in the alkaline olivine basalts (more than fourfold in some samples, see above).
It is quite difficult to interpret the genesis of the dolerite from the crustal xenolith suite on Zhokhov Island because of the ambiguity of the age estimates for this rock. If its age is 100 3 Ma, it is reasonable to assume that the dolerite was produced synchronously with the alkaline basalts of Bennett Island and, thus, composed, together with them, the Mesozoic mafic basement of the Laptev Sea shelf. This interpretation implies that the magmatic evolution of the continental shelf in the Laptev Sea at the Early and Late Cretaceous boundary was participated by at least two geochemically distinct mantle sources, one of which was responsible for the formation of plume-type melts (basalts of Bennett Island), and the other produced the parental melts of enriched MORB (dolerite from the crustal xenolith suite of Zhokhov Island). The geochemical similarities between the dolerite and enriched MORB is confirmed by similarities of some compositional parameters of the dolerite (see Figures 18a 18b) and basalts from the Gakkel Ridge. According to [Michael et al., 2003], the latter rocks are generally enriched more significantly than MORB in incompatible elements.
An alternative interpretation of the origin of the dolerite from xenoliths on Zhokhiv Island is based on the assumption that the radiogenic Ar contained in this rock was partly lost due to reheating, and the actual age of the dolerite is much older than 100 Ma. Then it should be assumed that the dolerite xenolith represents the ancient continental crust of the Laptev Sea shelf. A rough idea about its composition is provided by the Lower Paleozoic rocks of Bennett Island: Cambrian mudstones with trilobites and Ordovician mudstones and sandstones. The association of Paleozoic sedimentary and metasedimentary rocks can also include the quartzites found among xenoliths on Zhokhov Island.
The isotopic data presented above suggest that the alkaline basalts and limburgites of Bennett, Zhokhov, and Vil'kitskii islands could be generated by magmatic melts derived from a common mantle source or isotopically similar mantle sources. Within-plate melts with isotopic-geochemical signatures of the limburgites and alkaline olivine basalts from the De Long Islands could be produced by a mantle source whose chemical heterogeneity was comparable with the heterogeneity of the shallower sitting mantle material represented by spinel lherzolite mantle xenoliths on Zhokov and Vil'kistkii islands. The melting of such a mantle material, which included local zones enriched in lithophile elements, gave rise to chemically enriched melts with lower Nd isotopic ratios.
The results obtained by simulating the mixing process presented above are in conflict with the earlier hypothesis that the melts parental for the alkaline olivine basalts were contaminated with the material of the melanocratic basement of the Laptev Sea shelf [Bogdanovskii et al., 1992]. It is, thus, reasonable to believe that the magmatic system in which the within-plate volcanic associations of Bennett, Zhokhov, and Vil'kitskii islands were formed evolved without any notable chemical interaction with older crustal material.
As follows form the same model simulations, the Nd and Sr isotopic ratios observed in the limburgites and alkaline olivine basalts rule out the possibility that their parental melts mixed with the material of the DMM source. This conclusion implies that Mesozoic-Cenozoic magmatism at the De Long Islands was controlled exclusively by the melting of deep-seated enriched mantle sources, without involvement of shallower sources that produced N-MORB.
Judging from the isotopic-geochemical characteristics of volcanic rocks of the De Long Islands, these rocks are compositionally close to the most strongly enriched basaltoids in anomalous (plume) MAR segments or in oceanic islands (OIB). Hence, the alkaline basalts and limburgites of Bennett, Zhokhov, and Vil'kitskii islands can be classed with the products of plume magmatism. It is worth noting that the degree of enrichment of these rocks systematically increases from the alkaline basalts of Bennett Island to the alkaline olivine basalts and limburgites of Zhokhov Island and further to the limburgites of Vil'kitskii Island. This fact suggests that magmatism in the central part of the De Long Islands was contributed more significantly by the enriched mantle component as compare with North-West periphery of archipelago.
Figure 24 |
1. The Mesozoic mafic basement of the Laptev Sea shelf was produced at the boundary between the Early and Late Cretaceous. The origin of the magmatic complex during this stage of the geological history of the Laptev Sea shelf was participated by at least two geochemically distinct mantle sources. One of them gave rise to plume-type melts (the basalts of Bennett Island, which are the oldest magmatic rocks in this area), and the other was responsible for the origin of melts parental for the enriched MORB (dolerite from the crustal xenolith suite of Zhokhov Island). In the alternative variants, the dolerites should be ascribed to the material of the ancient continental crust of the Laptev Sea shelf, with this crustal material also including the Early Paleozoic sedimentary rocks of Bennett Island and quartzites of the crustal xenolith suite of Zhokhov Island.
2. After a fairly long (approximately 100 m.y.) interlude in the magmatic activity, the late stage magmatic of magmatic activity on the De Long Islands began. This stage was related to Miocene-Pliocene plume magmatism in the central part of the archipelago. The bulk of the volcanic pile of alkaline olivine basalts on Zhokhov Island was produced during the Late Miocene and Pliocene. The parental melts of these rocks had the same compositional characteristics as the melts that produced the within-plate basalts on Bennett Island. They were derived from a mantle source geochemically close to the reservoir that gave rise to the melts responsible for the Mesozoic within-plate magmatism in the northwestern termination of the De Long Islands. The alkaline basalt complex on Zhokhov Island was formed at relatively high degrees of melting, and the compositions of the rocks define fractionation trends typical of the picrite-alkaline basalt series of within-plate magmatism. When ascending to the surface, the magmatic melt from which the Zhokhov alkaline basalts crystallized entrained xenoliths of the ancient mantle material (spinel lherzolites and xenocrysts) and crustal rocks (quartzites, dolerites, and metavolcanic rocks). The "sampling'' of different depth levels in the crust and mantle by the magmatic melt was, perhaps, made possible by the high degree of melting of the mantle source and the low velocity of magma ascent to the surface.
The same time span (Pliocene) was marked by the evelopment of a small volcanic edifice on Zhokhov Island. This edifice consisted of limburgites, i.e., within-plate volcanic rocks atypical of Mesozoic plume magmatism at the De Long Islands. Judging from their isotopic signatures, these alkaline volcanics were derived from the same mantle source as the alkaline olivine basalts accompanying these rocks but at lower (possibly, much lower) degrees of melting.
3. The youngest manifestations of plume magmatism on the De Long Islands were discovered only on the southernmost of the examined islands, on Vil'kitskii Island. The limburgites that possibly make up the whole island are of Pliocene age, their composition is close to that of the limburgites of Zhokhov Island, and they were derived from the same mantle source as that for the Zhokhov limburgite. At the same time, the limburgites from Vil'kitskii Island show some geochemical characteristics suggesting that they are more enriched within-plate rocks than the analogous volcanics of Zhokhov Island. This could be caused by a lower degrees of melting at which the parental melts of the Vil'kitskii limburgites were derived and a higher degree of enrichment of their mantle source. Neither the limburgites in Vil'kitskii Island nor these rocks in the central part of Zhokhov Island show evidence of significant fractionation and are, perhaps, the most primitive volcanic rocks of the De Long Islands. The limburgites in the central part of the archipelago and related alkaline olivine basalts contain numerous mantle xenoliths of spinel lherzolites. No crustal xenoliths were found in the limburgites on either Zhokhov or Vil'kitskii islands.
The scenario proposed above for the magmatic evolution in the De Long Islands implies that it was closely related to the activity of a plume mantle source situated beneath the continental shelf in the Laptev Sea and responsible for the pulses of magmatic activity in this part of the eastern sector of the Arctic basin over the past 124 m.y. In this context, it is expedient to compare the spatial distribution of the geochemical types of volcanic rocks of different ages over the De Long Islands with data on Late Quaternary volcanism in Northeast Asia (see the review in [Grachev, 1999]). As is pointed out in this publication, the most striking feature of Late Quaternary volcanism in this vast area is widespread within-plate magmatism, which occurred at all types of the crust (oceanic or continental). Grachev [1999] suggested that Quaternary magmatism in Northeast Asia (including the De Long Islands) was controlled by the melting of mantle plume material with its insignificant lithospheric contamination. The geodynamic environment of magmatism in the shelf zone of the Laptev Sea, which was influenced by the propagation of the Gakkel mid-oceanic ridge into the shelf zone from the northwest, was referred to as pre-Afar [Grachev, 1999] and is thought to have corresponded to the early stages of the development of a mantle plume [Grachev, 1999].
The geodynamic style of magmatism of the De Long Islands makes it similar to magmatism in the Red Sea area, where the activity of two plumes (Afar at 45-0 Ma in the south and Levant-Nubian at 138-82 Ma in the north) predetermined the character of magmatism in the Red Sea Rift and north of it, along the Dead Sea and Jordan River valleys [Segev, 2000]. According to [Segev, 2000], these magmatic events were coupled with the extension of the continental lithosphere and the ascent of thermal anomalies corresponding to mantle plumes in the lower mantle. Thus, the development of rift zones in passive continental margins can be controlled by the ascent of one or several mantle plumes. It is worth noting that, according to [Sharkov, 2002], the northwestern part of the Arabian Peninsula (i.e., the area of the Levant-Nubian plume [Segev, 2000]) contains an association of Late Cenozoic within-plate volcanic rocks identical to the products of plume magmatism on the De Long Islands: alkaline basalts and basanites. As their within-plate analogues on Bennett and Vil'kitskii islands, these volcanic rocks, produced by magmatism triggered by the activity of a mantle plume and the opening of a young oceanic basin in the Red Sea area, bear mantle xenoliths of spinel lherzolites.
Compared with the magmatic complexes of the De Long Islands, the basaltoids in the Red Sea area are more heterogeneous chemically and comprise both oceanic tholeiites (including N-MORB), along with the acid derivatives of basaltoid magmatism (rhyolites and trachyandesites), and alkaline basalts [Barrat et al., 1990, 1993; Eissen et al., 1989; Schilling, 1969]. At the one hand, these differences could be caused by the more mature character of the spreading center in the Red Sea area, and, on the other hand, they could be controlled by kinematic factors. Indeed, magmatic complexes of the De Long Islands are compositionally similar to the plume-related rocks association in the northernmost part of the influence area of the Red Sea Rift, where there are no obvious traces of any oceanic spreading center (Levant-Nubian plume). According to [Segev, 2000], plume magmatism in this area becomes systematically younger from south to north starting from 244 Ma (Sinai) to 0.23 Ma (Golan Heights) or the Pleistocene (Shin and Syrian-Jordan plateaus [Sharkov, 2002]). Similar systematic variations in age are also typical of the volcanic products on the De Long Islands, where these variations are observed from northwest to southeast, from Bennett Island (124 Ma) to Vil'kitskii Island (0.4 Ma). It is pertinent to recall that, according to the data recently obtained in the course of the joint Russian-German seismic research in the Laptev Sea, its central part (in which we identified the youngest manifestations of plume magmatism) is marked by a Moho uplift [Piskarev, 2001].
Available data testify that the Red Sea spreading center, which controls the kinematics of lithospheric plates in the Middle East and West Asia, is characterized by extension velocities of about 1.6-2.2 cm/yr [Dubinin and Ushakov, 2001]. At the same time, the extension velocity in the greatest spreading zone of the Arctic Ocean, the Gakkel Ridge, near its intersection with the Laptev Sea shelf northwest of Bennett island is 0.7 cm/yr [Michael et al., 2003]. Evidently, the geodynamics of the Gakkel Ridge (which belongs, according to [Dick et al., 2003], to a new class of mid-oceanic ridges with ultraslow spreading) can significantly affect the style of the within-plate magmatism (first and foremost, the degree of melting and fractionation of the melts), which is related to its interaction with the lithosphere of the passive continental margin. Another distinctive feature of magmatism in this geodynamic environment is, perhaps, its very strong links with deep-seated enriched mantle sources and the absence of any indications that shallower DMM sources could have been involved in this process. Hence, the geochemical characteristics of magmatic products of the De Long Islands, which differentiate these rocks from magmatic associations in young fast-spreading rifts in passive continental margins, can be regarded as indirect arguments in support of the idea that the geological evolution of the continental shelf in the Laptev Sea was affected by the Gakkel mid-oceanic ridge propagating into this area.
Summarizing data on the magmatic complexes of the De Long Islands and related xenolith suite, it is logical to draw the conclusion that plume magmatism initiated rifting within the passive continental margin of the Laptev Sea shelf and served as a "conduit'' for the spreading center propagating into this area. Another important outcome of our research is the conclusion that there is a certain dependence between the spreading velocity in a mid-oceanic ridge interacting with a passive continental margin and the geochemistry of magmatic rocks produced in this margin. It is possible that the geographic position of the mantle plume center determines the trajectory along which the mid-oceanic ridge penetrates into the passive continental margin and, thus, harbingers the development of a would-be divergent boundary of lithospheric plates in the distant future.
Aoki, K.(1984), Petrology of materials derived from the upper mantle, in Materials Science of the Earth's Interior, edited by I. Sunagawa, 415 pp., Terrapub, Tokyo.
Backlund, H. (1920), On the eastern part of the Arctic basalt plateau, Med. Acad. Abo, Geol.-Mineral. Inst., 1, Abo.
Barrat, J. A., B. M. Jahn, S. Fourcade, and J. L. Joron (1993), Magma genesis in an ongoing rifting zone: the Tadjoura Gulf (Afar area), Geochim. Cosmochim. Acta, 57, 2291-2302.
Barrat, J.-A., B. M. Jahn, J. L. Joron, et al. (1990), Mantle heterogeneity in northeastern Africa: evidence from Nd isotopic compositions and hygromagmaphile element geochemistry of basaltic rocks from the Gulf of Tadjoura and southern Red Sea regions, Earth Planet. Sci. Lett., 101, 233-247.
Bogdanovskii, O. G., S. D. Mineev, S. S. Assonov, S. A. Silantyev, S. F. Karpenko, Yu. A. Shukolyukov, and L. A. Savostin (1992), Magmatism on the De Long Islands, eastern Arctic: isotopic geochemistry and geochronology, Geokhimiya, 1, 47-57.
Bogdanovskii, O. G., S. A. Silantyev, S. F. Karpenko, S. D. Mineev, and L. A. Savostin (1993), Ancient mantle xenoliths in young volcanic rocks on Zhokhov Island, De Long Islands, Dokl. Russ. Akad. Nauk., 330, 750-753.
Coombs, D. C., and J. F. G. Wilkinson (1969), Lineages and Fractionation trends in undersaturated rocks from the East Otago Volcanic Province, New Zealand, and related rocks, J. Petrol., 10, 440-501.
Dick, H. J. B., J. Lin, and H. Shouten (2003), An ultraslow-spreading class of ocean ridge, Nature,
Dosso, L., H. Bougault, C. Langmuir, et al. (1999), The age and distribution of mantle heterogeneity along the Mid-Atlantic Ridge (31o-41o N), Earth Planet. Sci. Lett., 170, 269-286.
Dubinin, E. P., and S. A. Ushakov (2001), Oceanic riftogenesis, 292 pp., GEOS, Moscow.
Eissen, J.-P., T. Juteau, J.-L. Joron, et al. (1989), Petrology and geochemistry of basalts from the Red Sea axial rift at 18o north, J. Petrol., 30, 791-839.
Fedorov, P. I., G. B. Flerov, and D. I. Golovin (2002), Volcanism on Bennett Island, East Arctic: new data on the composition and age of the rocks, in School on the Geochemistry of Magmatic Rocks, pp. 93-94, Vernadsky Institute of Geochemistry and Analytical Chemistry, Russ. Acad. Sci., Moscow.
Frey, F. A., and D. H. Green (1974), The mineralogy, geochemistry and origin of lherzolite inclusions in Victorian basanites, Geochim. Cosmochim. Acta, 38, 1023-1059.
Gakkel, Ya. Ya. (1957), Science and Exploring of the Arctic, 133 pp., Morskoi Transport, Leningrad.
Geology of the USSR (1984), vol. IX, 550 pp.
Geology of the USSR, Islands of the Soviet Arctic (1970), vol. XXVI, pp. 324-374.
Grachev, A. F. (1999), Quaternary volcanism and problems of geodynamics in Northeast Asia, Fiz. Zemli, 9, 19-37.
Hellebrand, E., J. E. Snow and R. Muhe (2002), Mantle melting beneath the Gakkel Ridge (Arctic Ocean): abyssal peridotite spinel composition, Chem. Geol., 182, 227-235.
Hofmann, A. W. (2003), Sampling mantle heterogeneity through oceanic basalts: isotopes and trace elements, in Treatise on Geochemistry, vol. 2, 61-101.
Ionov, D. A. (1988), Xenoliths in continental basalts, Ultramafic rocks of the earth's deep zones, in Magmatic Rocks, Ultramafic Rocks, pp. 311-332, Nauka, Moscow.
Ionov, D. A., and E. Jagoutz (1988), Sm and Nd isotopic composition of minerals in garnet and spinel peridotite xenoliths from the Vitim Highlands: the first data on mantle xenoliths in the USSR, Doklady AN SSSR, 301, 1195-1199.
Irving, A. J. (1980), Petrology and geochemistry of composite ultramafic xenoliths in alkalic basalts and implications for magmatic processes within the mantle, Am. J. Sci., 280-A, 389-426.
Jackson, E. D., and T. L. Wright (1970), Xenoliths in the Honolulu volcanic series, Hawaii, J. Petrol., 11, 405-430.
Kovalenko, V. I., D. A. Ionov, V. V. Yarmolyuk, E. Jagoutz, G. Lugmair, and H. G. Stosh (1990), Mantle evolution and its correlation with the crustal evolution in some areas in Central Asia: isotopic data, Geokhimiya, 9, 1308-1319.
Le Maitre, R. W. et al. (Eds.) (1989), Classification of Igneous Rocks and Glossary of Terms, Part I, Recommendations of the International Union of Geological Sciences Subcommission on the Systematics of Igneous Rocks, 150 pp. Blackwell, Oxford-London-Edinburgh-Boston-Melbourne.
Lindsley, D. H., and D. J. Anderson (1983), A two-pyroxene thermometer. Proc. Lunar Planet. Sci. Conf., Part 2: J. Geophys. Res., 88A, 887-906.
Mercier, J. C. (1980), Single-pyroxene thermobarometry, Tectonophysics, 70, 1-37.
Michael, P. J., C. H. Langmuir, H. J. B. Dick, J. E. Snow, S. L. Goldstein, D. W. Graham, K. Lehnert, G. Kurras, W. Jokat, R. Muhe, and H. N. Edmonds (2003), Magmatic and amagmatic seafloor generation at the ultraslow-spreading Gakkel Ridge, Arctic Ocean, Nature, 423, 956-961.
Mori, T., and D. H. Green (1978), Laboratory duplication of phase equilibria observed in natural garnet lherzolites, J. Geol., 86, 83-9.
Naryshkin, G. D. (1987), The mid-oceanic ridge of the Eurasian basin in the Arctic Ocean, in The Results of Investigations Within the Scope of International Geophysical Projects, 70 pp., Moscow.
O'Neill, H. St. C., and V. J. Wall (1987), The olivine-orthopyroxene-spinel oxygen geobarometer: the nickel precipitation curve and oxygen fugacity of the Earth's upper mantle, J. Petrol., 28, 1169-1193.
Piskarev, A. L. (2001), Sedimentation dynamics in the Laptev Sea as an indicator of changes in the spreading kinematics in the Eurasian basin of the Arctic Ocean, in Geology and Geophysics of Mid-Oceanic Ridges, Trans. of S. P. Mashchenkov workshop of the Russian branch of the InterRidge international project, p. 44, St. Petersburg.
Ringwood, A. E. (1975), Composition and Petrology of the Earth's Mantle, 570 pp., McGraw-Hill, New York-London-Toronto-Paris-Tokyo.
Sachtleben, Th., and H. A. Seck (1981), Chemical control of Al solubility in orthopyroxene and its implications for pyroxene geothermometry, Contrib. Mineral. Petrol, 78, 157-165.
Saggerson, E. P., and L. A. S. Williams (1964), Ngurumanite from southern Kenya and its bearing on the origin of rocks in the northern Tanganyika alkaline district, J. Petrol., 5, 40-81.
Savostin, L. A., S. A. Silantyev, and O. G. Bogdanovskii (1988), New data on volcanism at Zhokhov Island, De Long Islands, Arctic basin, Doklady AN SSSR, 302, 1443-1447.
Schilling, J.-G. (1969), Red Sea floor origin: rare-earth evidence, Science, 165, 1357-1360.
Segev, A. (2000), Synchronous magmatic cycles during the fragmentation of Gondwana: radiometric ages from the Levant and other provinces, Tectonophysics, 325, 257-277.
Sharkov, E. V. (2002), Interaction between the roof of a mantle plume and the earth crust during within-plate deformations: an example of the Late Cenozoic evolution of the Syrian territory, in Int. symposium on mantle plumes and metallogeny, pp. 285-287, Petrozavodsk-Moscow.
Silantyev, S. A., O. G. Bogdanovskii, and L. A. Savostin (2002), Geochemistry of within-plate oceanic rocks and related xenoliths in the central part of the De Long Islands as a clue to reconstructing the magmatic history of the Laptev Sea shelf, East Arctic basin, in Int. symposium on mantle plumes and metallogeny, pp. 223-225, Petrozavodsk-Moscow.
Silantyev, S. A., O. G. Bogdanovskii, L. A. Savostin, and N. N. Kononkova (1991), Magmatism at the De Long Islands, East Arctic: petrology and geochemistry of volcanic rocks and related xenoliths (Zhokhov and Vil'kitskii islands), Geokhimiya, 2, 267-277.
Sun, S.-S., and W. F. McDonough (1989), Chemical and isotopic systematics of oceanic basalts: implications for mantle composition and processes, in Magmatism in Ocean Basins. Geol. Soc. Spec. Publ. London, vol. 42, edited by A. D. Saunders and M. J. Norry, pp. 313-345.
Vaganov, V. I., and S. V. Sokolov (1988), Thermobarometry of Mafic Mineral Assemblages, 354 pp., Nedra, Moscow.
Vol'nov, D. A., and D. S. Sorokov (1961), Geology of Bennett Island, Trans. Inst. Geol. Arktiki, Antarktiki, Akad. Nauk USSR, 123, 25 pp.
Vol'nov, D. A., D. A. Voitsekhovskii, V. N. Voitsekhovskii, O. A. Ivanov, D. S. Sorokov, and D. S. Yashin (1970), New Siberian Islands, in Geology of the USSR, Islands of the Soviet Arctic, vol. XXVI, pp. 324-374.
Wass, S. Y. (1980), Geochemistry and origin of xenolith-bearing and related alkali basaltic rocks from the Southern Highlands, New South Wales, Australia, Am. J. Sci., 280-A, 639-666.
Wells, P. R. A. (1977), Pyroxene thermometry in simple and complex systems, Contrib. Mineral. Petrol., 62, 129-139.
Zhuravlev, A. Z., E. E. Laz'ko, and A. I. Ponomarenko (1991), Radiogenic isotopes and REE in minerals from garnet peridotite xenoliths in kimberlites of the Mir pipe, Yakutia, Geokhimiya, 7, 982-994.
Zindler, A., and S. Hart (1986), Chemical geodynamics, Ann. Rev. Earth Planet. Sci., 14, 439-571.