Thermal equation of state of magnesiowüstite (Mg0.6Fe0.4)O

Volume measurements for magnesiowüstite (Mg_0.6Fe_0.4)O, were carried out up to pressures of 10.1 GPa in the temperature range 300–1273 K, using energy-dispersive synchrotron X-ray diffraction. These data allow reliable determination of the temperature dependence of the bulk modulus and good constraint on the thermal expansitivity at ambient pressure which was previously not known for magnesiowüstite. From these data, thermal and elastic parameters were derived from various approaches based on the Birch–Murnaghan equation of state (EOS) and on the relevant thermodynamic relations. The results from three different equations of state are remarkably consistent. With (∂K_T/∂P)_T fixed at 4, we obtained K₀=158(2) GPa, (∂K_T/∂T)_P=−0.029(3) GPa K⁻¹, (∂K_T/∂T)_V=−3.9(±2.3)×10⁻³ GPa K⁻¹, and α_T=3.45(18)×10⁻⁵+1.14(28)×10⁻⁸T. The K₀, (∂K_T/∂T)_P, and (∂K_T/∂T)_V values are in agreement with those of Fei et al. (1992) and are similar to previously determined values for MgO. The zero pressure thermal expansitivity of (Mg_0.6Fe_0.4)O is found to be similar to that for MgO (Suzuki, 1975). These results indicate that, for the compositional range x=0–0.4 in (Mg_1−xFe_x)O, the thermal and elastic properties of magnesiowüstite exhibit a dependence on the iron content that is negligibly small, within uncertainties of the experiments. They are consequently insensitive to the Fe–Mg partitioning between (Mg, Fe)SiO₃ perovskite and magnesiowüstite when applied to compositional models of the lower mantle. With the assumption that (Mg_0.6Fe_0.4)O is a Debye-like solid, a modified equation of heat capacity at constant pressure is proposed and thermodynamic properties of geophysically importance are calculated and tabulated at high temperatures.     

Single-crystal elasticity of stishovite: New experimental data obtained using high-frequency resonant ultrasound spectroscopy and a Gingham check structure model

Single-crystal elasticity of stishovite was examined using a new experimental technique and an empirical macroscopic model. Employing high-frequency resonant ultrasound spectroscopy, single-crystal elastic constants of stishovite were determined: C₁₁ = 443(3), C₃₃ = 781(4), C₁₂ = 193(2), C₂₃ = 199(2), C₄₄ = 256(2), and C₃₃ = 316(2) GPa. The frequency range of the resonant ultrasound spectroscopy techniques was 6–20 MHz, which is much lower than the ～10 GHz range of the Brillouin scattering technique. Of the elastic constants, the shear elastic constants C₄₄ and C₆₆ are consistent with the average value of the previously mentioned Brillouin scattering. Conversely, the four elastic constants, C₁₁, C₃₃, C₁₂, and C₂₃, slightly deviate outside the range of previous Brillouin scattering results. The present results, except those for C₁₂, are consistent with recent lattice dynamic analysis of inelastic X-ray scattering data. The adiabatic bulk modulus was calculated as 298 GPa, which is smaller and more consistent with the result of compression experiments than any other Brillouin scattering results (301–312 GPa). The present result shows greater P-wave velocity anisotropy (24.7%) than any preceding work. To understand the unique elastic properties of stishovite, the Gingham check model was proposed and examined. The result shows that the octahedron of 6-coordinated Si in stishovite crystal has stiffness comparable to that of diamond.     

Mechanisms and kinetics of the post-spinel transformation in Mg2SiO4

Mechanisms and kinetics of the post-spinel transformation in Mg₂SiO₄ were examined at 22.7–28.2 GPa and 860–1200 °C by in situ X-ray diffraction experiments using synchrotron radiation combined with microstructural observations of the recovered samples. The post-spinel phases nucleated on spinel grain boundaries and grew with a lamellar texture. Under large overpressure conditions, reaction rims were formed along spinel grain boundaries at the initial stage of the transformation, whereas under small overpressure conditions, the transformation proceeded without formation of reaction rims. Mg₂SiO₄ spinel metastably dissociated into MgSiO₃ ilmenite and periclase, and stishovite and periclase as intermediate steps in the transformation into the stable assemblage of MgSiO₃ perovskite and periclase. Topotactic relationships were found in the transformation from spinel into ilmenite and periclase. Kinetic parameters in the Avrami rate equation, time taken to 10% completion, and the growth rate were estimated by analysis of the kinetic data obtained by in situ X-ray observations. The empirical activation energy for 10% transformation decreases with increasing pressure because the activation energy for nucleation becomes smaller at larger overpressure conditions. Extrapolations of the 10% transformation to ∼700 °C, which is the lowest temperature expected for the cold slabs at ∼700 km depth, suggest that overpressure of more than ∼1 GPa is needed for the transformation. Because the growth rate is estimated to be large even at low-temperatures of ∼700 °C and overpressures of 1 GPa, the depth of the post-spinel transformation in the cold slabs is possibly controlled by nucleation kinetics.     

High-pressure phase relations in the system CaAl4Si2O11–NaAl3Si3O11 with implication for Na-rich CAS phase in shocked Martian meteorites

High-pressure phase relations in the system NaAl₃Si₃O₁₁–CaAl₄Si₂O₁₁ were examined at 13–23 GPa and 1600–1900 °C, using a multianvil apparatus. A Ca-aluminosilicate with CaAl₄Si₂O₁₁ composition, designated CAS phase, is stable above about 13 GPa at 1600 °C. In the system NaAl₃Si₃O₁₁–CaAl₄Si₂O₁₁, the CAS phase dissolving NaAl₃Si₃O₁₁ component coexists with jadeite, corundum and stishovite below 22 GPa, above which the CAS phase coexists with Na-rich calcium ferrite, corundum and stishovite. At 1600 °C, the solubility of NaAl₃Si₃O₁₁ component in the CAS solid solution increases with increasing pressure up to about 50 mol% at about 22 GPa, above which the solubility decreases with pressure. The maximum solubility of NaAl₃Si₃O₁₁ component in the CAS phase increases with temperature up to around 70 mol% at 1900 °C at 22 GPa. The dissociation of NaAlSi₂O₆ jadeite to NaAlSiO₄ calcium ferrite plus stishovite occurs at about 22 GPa. Lattice parameters of the CAS phase with the hexagonal Ba-ferrite structure change with increase of the NaAl₃Si₃O₁₁ component: a-axis decreases and c-axis slightly increases, resulting in decrease of molar volume. Enthalpies of the CAS solid solutions were measured by high-temperature drop-solution calorimetry techniques. The results show that enthalpy of hypothetical NaAl₃Si₃O₁₁ CAS phase is much higher than the mixture of NaAlSi₂O₆ jadeite, corundum and stishovite and is close to that of the mixture of NaAlSiO₄ calcium ferrite, corundum and stishovite. When we adopt the Na:Ca ratio of 75:25 of the natural Na-rich CAS phase in a shocked Martian meteorite, Zagami, the phase relations determined above suggest that the natural CAS phase crystallized from melt at pressure around 22 GPa and temperature close to or higher than 2000–2200 °C. The inferred P, T conditions are consistent with those estimated using other high-pressure minerals in the shocked meteorite.     

Temperature dependence of the elastic moduli of ringwoodite

The elastic moduli of polycrystalline ringwoodite, (Mg_0.91Fe_0.09)₂SiO₄, were measured up to 470 K by means of the resonant sphere technique. The adiabatic bulk (K_S) and shear (μ) moduli were found to be 185.1(2) and 118.22(6) GPa at room temperature, and the average slopes of dK_S/dT and dμ/dT in the temperature range of the study were determined to be −0.0193(9) and −0.0148(3) GPa/K, respectively. Using these results, we estimate seismic wave velocity jumps for a pure olivine mantle model at 520 km depth. We find that the jump for the S-wave velocity is about 1.5 times larger than that for the P-wave velocity at this depth. This suggests that velocity jumps at the 520 km discontinuity are easier to detect using S-waves than P-waves.     

Spin state of ferric iron in MgSiO3 perovskite and its effect on elastic properties

Recent studies have indicated that a significant amount of iron in MgSiO₃ perovskite (Pv) is Fe³⁺ (Fe³⁺/ΣFe = 10–60%) due to crystal chemistry effects at high pressure (P) and that Fe³⁺ is more likely than Fe²⁺ to undergo a high-spin (HS) to low-spin (LS) transition in Pv in the mantle. We have measured synchrotron Mössbauer spectroscopy (SMS), X-ray emission spectroscopy (XES), and X-ray diffraction (XRD) of Pv with all iron in Fe³⁺ in the laser-heated diamond-anvil cell to over 100 GPa. Fe³⁺ increases the anisotropy of the Pv unit cell, whereas Fe²⁺ decreases it. In Pv synthesized above 50 GPa, Fe³⁺ enters into both the dodecahedral (A) and octahedral (B) sites approximately equally, suggesting charge coupled substitution. Combining SMS and XES, we found that the LS population in the B site gradually increases with pressure up to 50–60 GPa where all Fe³⁺ in the B site becomes LS, while Fe³⁺ in the A site remains HS to at least 136 GPa. Fe³⁺ makes Pv more compressible than Mg-endmember below 50 GPa because of the gradual spin transition in the B site together with lattice compression. The completion of the spin transition at 50–60 GPa increases bulk modulus with no associated change in density. This elasticity change can be a useful seismic probe for investigating compositional heterogeneities associated with Fe³⁺.     

Importance of reversible attachment in predicting E. coli transport in saturated aquifers from column experiments

Drinking water wells indiscriminatingly placed adjacent to fecal contaminated surface water represents a significant but difficult to quantify health risk. Here we seek to understand mechanisms that limit the contamination extent by scaling up bacterial transport results from the laboratory to the field in a well constrained setting. Three pulses of Escherichia coli originating during the early monsoon from a freshly excavated pond receiving latrine effluent in Bangladesh were monitored in 6 wells and modeled with a two-dimensional (2-D) flow and transport model conditioned with measured hydraulic heads. The modeling was performed assuming three different modes of interaction of E. coli with aquifer sands: (1) irreversible attachment only (best-fit k_i = 7.6 day⁻¹); (2) reversible attachment only (k_a = 10.5 and k_d = 0.2 day⁻¹); and (3) a combination of reversible and irreversible modes of attachment (k_a = 60, k_d = 7.6, k_i = 5.2 day⁻¹). Only the third approach adequately reproduced the observed temporal and spatial distribution of E. coli, including a 4-log₁₀ lateral removal distance of ∼9 m. In saturated column experiments, carried out using aquifer sand from the field site, a combination of reversible and irreversible attachment was also required to reproduce the observed breakthrough curves and E. coli retention profiles within the laboratory columns. Applying the laboratory-measured kinetic parameters to the 2-D calibrated flow model of the field site underestimates the observed 4-log₁₀ lateral removal distance by less than a factor of two. This is promising for predicting field scale transport from laboratory experiments.     

Elasticity of MgO to 130 GPa: Implications for lower mantle mineralogy

The aggregate shear wave velocities of MgO (periclase) have been determined throughout Earth's lower mantle pressure regime approaching 130 GPa using Brillouin spectroscopy in conjunction with synchrotron X-ray diffraction technique in a diamond anvil cell apparatus. We found that the extrapolations of the high-pressure shear wave velocities and shear moduli to ambient pressure are highly consistent with earlier studies. However, the measurements over a wide pressure range revealed that the pressure derivative of the shear modulus (dG/dP = G₀′) of MgO is 1.92(2), which is distinctly lower than that of previous lower-pressure experiments. Compared with the previous results on (Mg,Fe)O ferropericlase, there is no clear correlation between iron content and G₀′. We calculate that the shear wave velocity profile of lower mantle along the adiabatic geotherm applied by the lower G₀′ value of periclase can remarkably well reproduce the global seismological 1-D velocity profile model with uniform composition model. The best-fitting result indicates the possibility of a lower mantle mineralogy with ~ 92 vol.% silicate perovskite phase, implying that the bulk composition of lower mantle is likely not to be pyrolitic but more chondritic. The present acoustic measurements performed over the large pressure range have thus led us to a better understanding of compositional model of the Earth's lower mantle.     

Coda-Q and its lapse time dependence analysis in the interaction zone of the Dinarides,the Alps and the Pannonian basin

The single backscattering model was used to estimate total attenuation of coda waves (Q_c) of local earthquakes recorded on eight seismological stations in the complex area of the western continental Croatia. We estimated Q₀ and n, parameters of the frequency dependent coda-Q using the relation Q_c = Q₀fⁿ. Lapse time dependence of these parameters was studied using a constant 30 s long time window that was slid along the coda of seismograms. Obtained Q_c were distributed into classes according to their lapse time, t_L. For t_L = 20–50 s we estimated Q₀ = 45–184 and n = 0.49–0.94, and for t_L = 60–100 s we obtained Q₀ = 119–316 and n = 0.37–0.82. There is a tendency of decrease of parameter n with increasing Q₀, and vice versa. The rates of change of both Q₀ and n seem to decrease for lapse times larger than 50–80 s, indicating an alteration in rock properties controlling coda attenuation at depths of about 100–160 km. A very good correlation was found between the frequency dependence parameter n and the Moho depths for lapse times of 50, 60 and 70 s.     

Elastic anisotropy of ferromagnesian post-perovskite in Earth's D” layer

We report experimental observation of a sizable elastic anisotropy in a polycrystalline sample of ferromagnesian silicate in post-perovskite (ppv) structure. Using a novel composite X-ray transparent gasket to contain and synthesize ppv in a panoramic diamond-anvil cell along with oblique X-ray diffraction geometry, we observed the anisotropic lattice strain and {1 0 0} or {1 1 0} slip-plane texture in the sample at 140 GPa. We deduced the elasticity tensor (c_ij), orientation-dependent compressional wave velocities, polarization-dependent shear-wave velocities, and the velocity anisotropy of the silicate ppv. Our results are consistent with calculations and indicate that with sufficient preferred orientation, the elastic anisotropy of this phase can produce large shear-wave splitting.     

Age estimation of the deposition of the cretaceous Upper Sandstones and a later remagnetization event,from Hirna area,SE Ethiopia

Outcrops of the Cretaceous Upper sandstone formation some 375 km to the East of Addis Ababa on the motor Highway to Harar was paleomagnetically investigated. About seventy core samples were collected at various stratigraphic levels from 250–300 meters thick sedimentary formation. After standard sample preparations in the laboratory the resulting specimens were subjected to routine paleomagnetic demagnetization protocol. In the first steps of demagnetizations process the recent and viscous magnetizations were removed by heating until a temperature of level of 300 °C. Further demagnetization of the samples resulted in the isolation of the final magnetization with stable line segments that is directed towards the origin, which is interpreted as Characteristic Remanent Magnetization (ChRM). Rock – magnetic experiments have identified goethite (αFeOOH), hematite (αFe₂O₃), detritial hematite, and magnetite as the magnetic mineral phases carrying the remanence. The ChRM identified resulted in an average value of (D_s = 0.5°, I_s = ?0.7°, α₉₅ = 4.3°, N = 34) for the red sandstones while an average value of (D_s = 335.8°, I_s = ?31.8°, α₉₅ = 4.7°, N = 14) for the limestone intercalations. The former ChRM in the red sandstone is determined to be secondary while the latter ChRM is known to be primary. Comparison of these directional results and their pole equivalents with the African plate Apparent Polar Wander Path curve established by Besse and Courtillot (2003) give ages of between 115–130 Million years for the limestone intercalation and ages of 30 million years for red sandstone unit. These are interpreted respectively as estimates of the age of deposition and a later remagnetization respectively.     

The oxidation state of sulfur in magmatic fluids

Sulfur compounds in volcanic gases are responsible for the global cooling after explosive eruptions and they probably controlled the early evolution of the Earth's atmosphere. We have therefore studied the oxidation state of sulfur in aqueous fluids under the pressure and temperature conditions and oxygen fugacities typical for magma chambers (0.5–3 kbar, 650–950 °C, Ni–NiO to Re–ReO₂ buffer conditions). Sulfur speciation was determined by Raman spectroscopy of quenched fluids trapped as inclusions in quartz. Our results show that sulfur in hydrothermal fluids and volcanic gases is much more oxidized than previously thought and in particular, some explosive eruptions may release a significant fraction of sulfur as SO₃ or its hydrated forms. In the pressure range from 500 to 2000 bar, the equilibrium constant K₁ of the reaction 2H₂S + 3O₂ = 2SO₂ + 2H₂O in aqueous fluids can be described by lnK₁ = ?(57.1 ± 7.1) + (173,480 ± 7592)T^? 1, where T is temperature in Kelvin. The equilibrium constant K₂ for the reaction SO₂ + ½O₂ = SO₃ in aqueous fluids, where SO₃ may include hydrated forms, such as H₂SO₄, was found to be strongly pressure dependent, with lnK₂ = ?(5.2 ± 5.7) + (19,243 ± 5993)T^? 1 at 1500 bar; lnK₂ = ?(11.1 ± 1.3) + (25,383 ± 1371)T^? 1 at 2000 bar and lnK₂ = ?(22.1 ± 2.2) + (37,082 ± 2248)T^? 1 at 2500 bar. Our data imply that volcanoes may directly inject hexavalent sulfur in the form of H₂SO₄ into the atmosphere, not only on Earth, but possibly also on Venus and on Mars, when it was still tectonically active. Remote measurements from satellites may have underestimated the sulfur yield of some recent eruptions. Moreover, the mechanisms of the interaction of volcanic gases with the stratosphere need to be reconsidered.     

Balancing the global oceanic neodymium budget: Evaluating the role of groundwater

The distinctly different, ε_Nd(0) values of the Atlantic, Indian, and Pacific Oceans requires that the residence time of Nd in the ocean (i.e., τ_Nd) be on the order of, or less than, the ocean mixing time of ∼ 500–1500 yr. However, estimates of τ_Nd, based on river influxes, range from 4000 to 15,000 yr, thus exceeding the ocean mixing time. In order to reconcile the oceanic Nd budget and lower the residence time by roughly a factor of 10, an additional, as yet unidentified, and hence “missing Nd flux” to the ocean is necessary. Dissolution of materials deposited on continental margins has previously been proposed as a source of the missing flux. In this contribution, submarine groundwater discharge (SGD) is examined as a possible source of the missing Nd flux. Neodymium concentrations (n = 730) and ε_Nd(0) values (n = 58) for groundwaters were obtained from the literature in order to establish representative groundwater values. Mean groundwater Nd concentrations and ε_Nd(0) values were used along with recent estimates of the terrestrial (freshwater) component of SGD (6% of river discharge on a global basis) to test whether groundwater discharge to the coastal oceans could account for the missing flux. Employing mean Nd concentrations of the compiled data base (i.e., 31.8 nmol/kg for all 730 analyses and 11.3 nmol/kg for 141 groundwater samples from a coastal aquifer), the global, terrestrial-derived SGD flux of Nd is estimated to range between 2.9 × 10⁷ and 8.1 × 10⁷ mol/yr. These estimates are of the same order of magnitude, and within a factor of 2, of the missing Nd flux (i.e., 5.4 × 10⁷ mol/yr). Applying the SGD Nd flux estimates, the global average ε_Nd(0) of SGD is predicted to be − 9.1, which is similar to our estimate for the missing Nd flux (− 9.2), and in agreement with the mean (± S.D.) ε_Nd(0) measured in groundwaters (i.e., ε_Nd(0) = −8.9 ± 4.2). The similarities in the estimated SGD Nd flux and corresponding ε_Nd(0) values to the magnitude and isotope composition of the missing Nd flux are compelling, and suggest that discharge of groundwater to the oceans could account for the missing Nd flux. Future investigations should focus on quantifying the Nd concentrations and isotope compositions of groundwater from coastal aquifers from a variety of coastal settings, as well as the important geochemical reactions that effect Nd concentrations in subterranean estuaries in order to better constrain contributions of SGD to the oceanic Nd budget.     

An evaluation of methods for determining during-storm precipitation phase and the rain/snow transition elevation at the surface in a mountain basin

Determining surface precipitation phase is required to properly correct precipitation gage data for wind effects, to determine the hydrologic response to a precipitation event, and for hydrologic modeling when rain will be treated differently from snow. In this paper we present a comparison of several methods for determining precipitation phase using 12 years of hourly precipitation, weather and snow data from a long-term measurement site at Reynolds Mountain East (RME), a headwater catchment within the Reynolds Creek Experimental Watershed (RCEW), in the Owyhee Mountains of Idaho, USA. Methods are based on thresholds of (1) air temperature (T_a) at 0 °C, (2) dual T_a threshold, −1 to 3 °C, (3) dewpoint temperature (T_d) at 0 °C, and (4) wet bulb temperature (T_w) at 0 °C. The comparison shows that at the RME Grove site, the dual threshold approach predicts too much snow, while T_a, T_d and T_w are generally similar predicting equivalent snow volumes over the 12 year-period indicating that during storms the cloud level is at or close to the surface at this location. To scale up the evaluation of these methods we evaluate them across a 380 m elevation range in RCEW during a large mixed-phase storm event. The event began as snow at all elevations and over the course of 4 h transitioned to rain at the lowest through highest elevations. Using 15-minute measurements of precipitation, changes in snow depth (z_s), T_a, T_d and T_w, at seven sites through this elevation range, we found precipitation phase linked to the during-storm surface humidity. By measuring humidity along an elevation gradient during the storm we are able to track changes in T_d to reliably estimate precipitation phase and effectively track the elevation of the rain/snow transition during the event.     

Patterns of significant seismic quiescence on the Mexican Pacific coast

Many authors have proposed that the study of seismicity rates is an appropriate technique for evaluating how close a seismic gap may be to rupture. We designed an algorithm for identification of patterns of significant seismic quiescence by using the definition of seismic quiescence proposed by Schreider (1990). This algorithm shows the area of quiescence where an earthquake of great magnitude may probably occur. We have applied our algorithm to the earthquake catalog on the Mexican Pacific coast located between 14 and 21 degrees of North latitude and 94 and 106 degrees West longitude; with depths less than or equal to 60 km and magnitude greater than or equal to 4.3, which occurred from January, 1965 until December, 2014. We have found significant patterns of seismic quietude before the earthquakes of Oaxaca (November 1978, M_w = 7.8), Petatlán (March 1979, M_w = 7.6), Michoacán (September 1985, M_w = 8.0, and M_w = 7.6) and Colima (October 1995, M_w = 8.0). Fortunately, in this century earthquakes of great magnitude have not occurred in Mexico. However, we have identified well-defined seismic quiescences in the Guerrero seismic-gap, which are apparently correlated with the occurrence of silent earthquakes in 2002, 2006 and 2010 recently discovered by GPS technology.     

Recycling of harbor sediment as lightweight aggregate

Sediment sampled from Taichung Harbor was mixed with local reservoir sediment at different weight ratios to prepare lightweight aggregate at 1050, 1100, and 1150 °C. A pressure of 3000 or 5000 psi was used to shape the powder mixtures into pellets before the heating processes. The results indicate that the leaching levels of trace metals from the lightweight aggregate samples are considerably reduced to levels less than Taiwan Environmental Protection Administration regulatory limits. Increasing final process temperature tends to reduce the bulk density and crushing intensity of lightweight aggregate with a concomitant increase in water sorption capability. Lightweight aggregate with the lowest bulk density, 0.49 g cm⁻³ for the 5000 psi sample, was obtained with the heating process to 1150 °C. Based on the X-ray absorption near edge structure results, FeSO₄ decomposition with a concomitant release of SO_x (x = 2, 3) is suggested to play an important role for the bloating process in present study.     

Paleomagnetism of the Late Cretaceous ignimbrite from the Okhotsk-Chukotka Volcanic Belt,Kolyma-Omolon Composite Terrane: Tectonic implications

New Late Cretaceous paleomagnetic results from the Okhotsk-Chukotka Volcanic Belt in the Kolyma-Omolon Composite Terrane yield stable and consistent remanent directions. The Late Cretaceous (86–81 Ma) ignimbrites from the Kholchan and Ola suites were sampled at 19 sites in the Magadan area (60.4° N, 151.0° E). We isolated the characteristic paleomagnetic directions from 16 sampled sites using an alternating field demagnetization procedure. The primary nature of these directions is ascertained by dual polarities and positive fold tests. A tilt-corrected mean direction (D = 42.8°, I = 84.7°, k = 46.0, α₉₅ = 10.0°) yields a paleomagnetic pole of 66.7° N, 168.5° E (A₉₅ = 18.8°) which appears almost identical to the 90–67 Ma pole reported from the Lake El’gygytgyn area of the Okhotsk-Chukotka Volcanic Belt (Chukotka Terrane). This consistency suggests that the Kolyma-Omolon Composite Terrane and Chukotka Terrane has acted as a single tectonic unit since 80 Ma without any significant internal deformation. Accordingly, we calculate a combined 80 Ma characteristic paleomagnetic pole (Long. = 164.7° E, Lat. = 68.0°, A₉₅ = 10.9°, N = 12) for the Kolyma-Omolon-Chukotka Block which falls 16.5–17.5° south of the same age poles from Europe and East Asia. We ascribe this discrepancy in pole positions to tectonic activity in the area and infer a southward displacement of 1640 ± 1380 km for the Kolyma-Omolon-Chukotka Block with respect to the North American and Eurasian blocks since 80 Ma; more than 260 km of it is attributed to tectonic displacement in the Arctic Ocean due to the opening of the Canadian Basin.     

Formation of the Late Permian Panzhihua plutonic-hypabyssal-volcanic igneous complex: Implications for the genesis of Fe–Ti oxide deposits and A-type granites of SW China

The Late Permian (260 Ma) Emeishan large igneous province of SW China contains numerous magmatic Fe–Ti oxide deposits. The Fe–Ti oxide deposits occur in the lower parts of evolved layered gabbroic intrusions which are spatially and temporally associated with A-type granitic rocks. The 260 Ma Panzhihua layered gabbroic intrusion hosts one of the largest magmatic Fe–Ti oxide deposits in China and is coeval with a peralkaline A-type granitic pluton. The granite has intruded the overlying Emeishan flood basalts and fed at least one dyke which erupted onto the surface producing columnar jointed trachytes. The presence of syenodiorite between the layered gabbro and granite is evidence for compositional evolution from mafic to intermediate to felsic rocks. The syenodiorites have intermediate to felsic composition with SiO₂ = 61 to 65 wt.%, MgO = 0.27 to 0.6 wt.% and CaO = 1.0 to 2.5 wt.% as compared to the granite SiO₂ = 65 to 72 wt.%, MgO = 0.1 to 0.4 wt.%, CaO = < 1.0 wt.%. Primitive-mantle-normalized incompatible element plots show corresponding reciprocal patterns between the mafic and felsic rocks. The chondrite-normalized REE patterns show Eu anomalies changing from > 1(Eu/Eu? = 1.1 to 2.6) in the gabbroic intrusion, to < 1 in the syenodiorite (Eu/Eu? = 0.75 to 0.83), granites and trachytes (Eu/Eu? = 0.55–0.87). Previously published εNd_(T) values from clinopyroxenes (εNd_(T) = + 1.1 to + 3.2) of the gabbroic intrusion match the whole-rock values of the syenodiorite (εNd_(T) = + 2.1 to + 2.5), granite and trachyte (εNd_(T) = + 2.2 to + 2.9), suggesting that all rock types originated from the same mantle source. MELTS and trace element modeling confirm that all rock types can be generated by fractional crystallization of high-Ti Emeishan basalt. The jump in SiO₂ from the gabbro to the syenodiorite is attributed to the en masse crystallization of the Fe–Ti oxides. The geological and geochemical data indicate that fractional crystallization of a common parental magma produced the layered gabbroic intrusion and Fe–Ti oxide deposit, the syenodiorite, granites and trachyte of the Panzhihua region, which thus form a genetically related plutonic-hypabyssal-volcanic complex. Other granite–gabbro complexes in the region likely formed in a similar manner.