Olivine-spinel geospeedometry: Analysis of diffusion-controlled Mg-Fe2+ exchange

Subsolidus Mg-Fe²⁺ exchange between olivine and spinel is governed by Mg-Fe²⁺ interdiffusion. Incomplete exchange results in Mg-Fe²⁺ heterogeneity in both olivine and spinel, which provides information on the thermal histories of the host rocks. A composite sphere model has been developed to obtain quantitative cooling rates or heating duration from the Mg-Fe²⁺ heterogeneity. The model assumes that a spherical core of spinel and a surrounding semi-infinite spherical shell of olivine interact by diffusion-controlled exchange of Mg and Fe²⁺. The differential equations describing the model are solved numerically by finite difference approximations. The numerical solution reveals that cooling rates or heating duration can be estimated from the relationship between the grain size of spinel and temperature calculated from the chemical compositions of the core of a spinel grain and of olivine far away from it. The calculated temperature is employed in place of $\text{Mg}\text{(Mg + Fe}{}^{\mathrm{2+}}\text{)}$ at the center of spinel to obtain the absolute temperature of thermal events.This olivine-spinel geospeedometer has been applied to peridotites. gabbro, and picrites from some ophiolite complexes in Japan to estimate their cooling rates. The estimated cooling rates for the peridotites range from 10^?4 to 10^?1 °C/yr, and those for the picrites from 10³ to 10⁴ °C/yr. The geospeedometer has been extended to estimate the heating duration of lherzolite xenoliths in basalt from the Ichinomegata crater, northeast Japan. The estimated heating duration of the xenoliths is less than one day.     

Long-term variability of winter mixed layer depth and temperature along the Kuroshio jet in a high-resolution ocean general circulation model

To explore the causes of the winter shallow mixed layer and high sea surface temperature (SST) along the strong Kuroshio jet from the East China Sea to the upstream Kuroshio extension (25.5°N–150°E) during 1988–1994 when the Japanese sardine stocks collapsed, high-resolution ocean general circulation model (OGCM) hindcast data are analyzed with a bulk mixed layer model which traces particles at the mixed layer base. The shallow mixed layer and high SST along the Kuroshio jet are mainly caused by the acceleration of the Kuroshio current velocity and the reduction of the surface cooling. Because the acceleration reduces the time during which the mixed layer is exposed to wintertime cooling, deepening and cooling of the winter mixed layer are restricted. The weaker surface cooling due to less severe meteorological forcing also causes the shallow mixed layer and the high SST. The impact of the strong heat transport along the Kuroshio extends to the southern recirculation gyre of the Kuroshio/Kuroshio extension regions; previous indications that the Japanese sardine recruitment is correlated with the winter SST and the mixed layer depth (MLD) in the Kuroshio extension recirculation region could be related to the velocity, SST, and MLD near the Kuroshio axis which also could affect the variability of North Pacific subtropical water.     

Sea-ice retreat in the Sea of Okhotsk and the ice–ocean albedo feedback effect on it

Sea-ice retreat processes are examined in the Sea of Okhotsk. A heat budget analysis in the sea-ice zone shows that net heat flux from the atmosphere at the water surface is about 77 W m⁻² on average in the active ice melt season (April) due to large solar heating, while that at the ice surface is about 12 W m⁻² because of the difference in surface albedo. The temporal variation of the heat input into the upper ocean through the open water fraction corresponds well to that of the latent heat required for ice retreat. These results suggest that heat input into the ice–upper ocean system from the atmosphere mainly occurs at the open water fraction, and this heat input into the upper ocean is an important heat source for ice melting. The decrease in ice area in the active melt season (April) and the geostrophic wind just before the melt season (March) show a correlation: the decrease is large when the offshoreward wind is strong. This relationship can be explained by the following process. Once ice concentration is decreased (increased) by the offshoreward (onshoreward) wind just before the melt season, solar heating of the upper ocean through the increased (decreased) open water fraction is enhanced (reduced), leading to (suppressing) a further decrease in ice concentration. This positive feedback is regarded as the ice–ocean albedo feedback, and explains in part the large interannual variability of the ice cover in the ice melt season.     

Anisotropy of Mg isotopic fractionation during evaporation and Mg self-diffusion of forsterite in vacuum

Evaporation of solid materials under low-pressure conditions could play important roles in chemical and isotopic fractionations in the early solar system. We have studied anisotropy of isotopic fractionation of ²⁶Mg and ²⁵Mg during kinetic evaporation of forsterite (Mg₂SiO₄), which is potentially a powerful tool to understand thermal histories of crystals in the early solar system. Ion-microprobe depth profiling revealed that the Mg isotopic zoning profiles of forsterite evaporated at 1500-1700 °C are notably differing along the a-, b-, and c-axes, which can be attributed to anisotropy in self-diffusion coefficient of Mg (D) and an isotopic fractionation factor for evaporation of Mg (α). The D and α were obtained from zoning profiles by applying the diffusion-controlled isotopic fractionation model of Wang et al. [1999. Evaporation of single crystal forsterite: Evaporation kinetics, magnesium isotope fractionation, and implications of mass-dependent isotopic fractionation of a diffusion-controlled reservoir. Geochim. Cosmochim. Acta 63(6), 953-966.].The D is largest and smallest along the a- and c-axes, respectively. The activation energy of 560-670 kJ/mol indicates that Mg diffusion at 1500-1700 °C occurred in the intrinsic diffusion regime.The α seems to be larger along the a- or c-axes than along the b-axis. The α along the a- or c-axes show weak temperature dependence. The α along all the crystallographic orientations is closer to unity than that expected from the kinetic theory of gases. These lines of evidence suggest that surface processes such as breaking of bonds and surface diffusion are responsible for the isotopic fractionation.     

Felsic dykes in the Neoproterozoic Nagar Parkar Igneous Complex,SE Sindh,Pakistan: geochemistry and tectonic settings

The Nagar Parkar Igneous Complex consists of Neoproterozoic igneous and metamorphic rocks dissected by mafic, felsic, and rhyolitic dykes. The latter can be classified broadly into porphyritic felsic dykes intruding gray and pink granites at Nagar Parkar and the surrounding areas, and the orthophyric felsic dykes intruding amphibolites, deformed pink granites, and the alkaline mafic dykes in the Dhedvero area, north of Nagar Parkar. The porphyritic felsic dykes are composed of perthites, quartz, and albitic plagioclase whereas the orthopheric felsic dykes contain K-feldspar (dominant), plagioclase, and minor quartz. Geochemically, the porphyritic and orthophyric felsic dykes are subalkaline and alkaline demonstrating post-orogenic A2- and OIB-A1-type characteristic on Nb–Y–Ce and Nb–Y–3Ga ternary plots, respectively. One orthophyric felsic dyke contains normative acmite and sodium metasilicate. This study suggests two distinct tectonic regimes for the origin of the felsic dykes of the area. The porphyritic felsic dykes show similarities with the ~800–700 Ma granites of the area, the rhyolite dykes of the Mount Abu, western Rajasthan in India, and the granites of the Seychelles microcontinent. The orthophyric felsic dykes show chemical resemblance with the Tavidar volcanic suite of western Rajasthan and the Silhouette and North islands of the Seychelles microcontinent. This study confirms spatial and temporal links among the Rodinian fragments exposed in the Nagar Parkar area of Pakistan, western Rajasthan of India, and the Seychelles microcontinent.     

Mineral classification from quantitative X-ray maps using neural network: Application to volcanic rocks

We developed a mineral classification technique of electron probe microanalyzer (EPMA) maps in order to reveal the mineral textures and compositions of volcanic rocks. In the case of lithologies such as basalt that include several kinds of minerals, X-ray intensities of several elements derived from EPMA must be considered simultaneously to determine the mineral map. In this research, we used a Kohonen self-organizing map (SOM) to classify minerals in the thin-sections from several X-ray intensity maps. The SOM is a type of artificial neural network that is trained using unsupervised training to produce a two-dimensional representation of multi-dimensional input data. The classified mineral maps of in situ oceanic basalts of the Juan de Fuca Plate allowed us to quantify mineralogical and textural differences among the marginal and central parts of the pillow basalts and the massive flow basalt. One advantage of mineral classification using a SOM is that relatively many minerals can be estimated from limited input elements. By applying our method to altered basalt which contains multiple minerals, we successfully classify eight minerals in thin-section.     

Feeding chronology of the vertically migrating gonostomatid fish,Vinciguerria nimbaria (Jordan and Williams), off southern Japan

The feeding habits of the gonostomatid fish,Vinciguerria nimbaria (Jordan andWilliams), from off southern Japan were studied in relation to its food organisms, process of digestion, and diel vertical migration. Food organisms were composed mainly of small- to moderate-sized copepods which live almost entirely in the 0–200 m epipelagic zone.V. nimbaria appears to feed twice a day; after having migrated to the epipelagic zone in the evening, all individuals fed intensively from 6 to 8 p.m. They then spent the remainder of the night in digestion. Many, but not all, of the individuals foraged again in the morning.     

Cr and Al diffusion in chromite spinel: experimental determination and its implication for diffusion creep

Diffusion coefficients of Cr and Al in chromite spinel have been determined at pressures ranging from 3 to 7 GPa and temperatures ranging from 1,400 to 1,700°C by using the diffusion couple of natural single crystals of MgAl₂O₄ spinel and chromite. The interdiffusion coefficient of Cr–Al as a function of Cr# (=Cr/(Cr + Al)) was determined as D _Cr–Al = D ₀ exp {−(Q′ + PV*)/RT}, where D ₀ = exp{(10.3 ± 0.08) × Cr#^0.54±0.02} + (1170 ± 31.2) cm²/s, Q′ = 520 ± 81 kJ/mol at 3 GPa, and V* = 1.36 ± 0.25 cm³/mol at 1,600°C, which is applicable up to Cr# = 0.8. The estimation of the self-diffusion coefficients of Cr and Al from Cr–Al interdiffusion shows that the diffusivity of Cr is more than one order of magnitude smaller than that of Al. These results are in agreement with patterns of multipolar Cr–Al zoning observed in natural chromite spinel samples deformed by diffusion creep.     

Sound velocities of Na0.4Mg0.6Al1.6Si0.4O4 NAL and CF phases to 73 GPa determined by Brillouin scattering method

The sound velocities of two aluminum-rich phases in the lower mantle, hexagonal new Al-rich phase (NAL) and its corresponding high-pressure polymorph orthorhombic Ca-ferrite-type phase (CF), were determined with the Brillouin scattering method in a pressure range from 9 to 73 GPa at room temperature. Both NAL and CF samples have identical chemical composition of Na_0.4Mg_0.6Al_1.6Si_0.4O₄ (40 % NaAlSiO₄–60 % MgAl₂O₄). Infrared laser annealing in the diamond anvil cell was performed to minimize the stress state of the sample and obtain the high-quality Brillouin spectra. The results show shear modulus at zero pressure G ₀ = 121.960 ± 0.087 GPa and its pressure derivative G’ = 1.961 ± 0.009 for the NAL phase, and G ₀ = 129.653 ± 0.059 GPa and G’ = 2.340 ± 0.004 for the CF phase. The zero-pressure shear velocities of the NAL and CF phases are obtained to be 5.601 ± 0.005 km/sec and 5.741 ± 0.001 km/sec, respectively. We also found that shear velocity increases by 2.5 % upon phase transition from NAL to CF at around 40 GPa.     

Isochemical breakdown of garnet in orogenic garnet peridotite and its implication to reaction kinetics

An isochemical kelyphite (orthopyroxene+spinel+plagioclase) that has nearly the same bulk chemical composition as the precursor garnet was found within a matrix of ordinary kelyphites (orthopyroxene+clinopyroxene+spinel±amphibole) in garnet peridotites from the Czech part of the Moldanubian Zone. It was shown that the kelyphitization of garnet took place in three stages: (1) the garnet-olivine reaction, accompanied by a long-range material transfer across the reaction zone, and (2) the isochemical breakdown of garnet, essentially in a chemically-closed system, and finally, (3) an open-system hydration reaction producing a thin hydrous zone (amphibole+spinel+plagioclase), which is located between the isochemical kelyphite and relict garnet. The presence of relict garnet suggests that this breakdown reaction of the second stage did not proceed to a completion probably being hindered by the formation of the hydrous zone at the reaction front. It was found by electron back-scattered diffraction method that orthopyroxene and spinel do not show any topotaxic relationship in the first type of kelyphite; whereas they show locally topotaxic relationship in the isochemical kelyphite. The transition from the first type to the second type of kelyphite is discussed on the basis of the detailed observations in the transition zone between the two kelyphites. More widespread occurrence of isochemical kelyphite is expected to occur in orogenic peridotites as well as from xenoliths brought by volcanics.