Hydrogen and carbon isotope studies on the graphite-bearing metapelites in the northern Kiso district of central Japan

Measurements were made of the hydrogen isotope ratios of hydrous silicates (mica and amphibole) and whole rocks, and the carbon isotope ratios of graphite and carbonaceous matter in the metamorphic rocks from the northern Kiso district in central Japan.D values of hydrous silicates in the graphite-bearing metapelites are always higher than those in graphite-free schists, even though the sample localities of the two rock-types are very close. Hydrogen isotopic equilibrium has been attained between the coexisting minerals.D/H ratios of water in the metamorphic fluids seem to depend strongly on the presence or absence of graphite and seem to be not constant throughout the district. The district is divided into three areas of low (metamorphic zones I, II), medium (zones IIIa–V) and high ¹³C_gr value (zones VIa–VII) areas. In the high ¹³C_gr values area, the carbon contents of the graphite-bearing rocks decrease slightly from zones VIa to VII, whereas the ¹³C_gr values increase sharply from the upper part of zone VIa to VIb. TheD values of biotite in these graphite-bearing rocks are higher than those in the medium ¹³C_gr area. This suggests that methane enriched inH and¹²C is produced and liberated by the devolatilization reactions between muscovite, graphite and water. The fluid produced is composed of water, methane and a subordinate amount of carbon dioxide, and its logfO₂ value is deduced to be about 1.2 lower than that defined by the FMQ buffer. In the medium ¹³C_gr area, the ¹³C values of graphite are nearly constant (–20.8), while the Fe₂O₃/(Fe₂O₃ + FeO) ratio of the graphite-bearing rocks apparently decreases with increasing metamorphic grade.D differences in hydrous silicates between graphite-bearing and graphite-free rocks are observed. These facts are interpreted to mean that methane was produced in addition to water and carbon dioxide, and that its generation ( ratio of the fluid was about 2) had practically no isotope effect on the graphite. In the low ¹³C_gr area, the carbon contents of the rocks decrease clearly from zones I to IIIa. TheD and ¹³C_gr values of the non-metamorphosed shales are much lower than those of the low grade graphite-bearing metapelites. This suggests that methane is produced and liberated from the rocks even at the incipient stage of metamorphism.     

Structure and lattice vibrations of Mg-Al spinel solid solution

X-ray structure refinements have been made for nonstoichiometric (MgO · 3Al₂O₃) and stoichiometric Mg-Al spinels. Several structure variations with chemical composition have been observed and are discussed in relation to Al substitution in tetrahedral sites. Infrared reflection and Raman spectra of the single crystal of the nonstoichiometric spinel (MgO · 3Al₂O₃) have been measured and analyzed. The results obtained are compared with those reported for the stoichiometric sample. From the infrared and Raman frequencies reported for the stoichiometric Mg-Al spinel, which are partly complemented with our results, the effective ionic charges of the ions in MgAl₂O₄ have been estimated on the basis of the rigid ion model.     

Bending of spherical lithosphere–axisymmetric case

Effects of sphericity are commonly ignored in the lithospheric bending problem. In order to examine its effects, I solve a simple axisymmetric spherical-shell model. The full solution and the asymptotic solution are derived from the basic equations, and their relationship to the flat-plate solution is examined. For displacement, effects of sphericity are small, and use of the flat-plate solution produces results that are numerically indistinguishable from those of the spherical solution. The most significant effect of sphericity appears in the stress, in particular the normal stress along the strike direction of the trench. This stress is approximately given by Eu_r/R , where E is Young's modulus, u_r is the vertical deformation of the shell and R is its radius of curvature. If the shell (lithosphere) is bent downwards and reaches 30 km, this stress can become about 5 kbar in the Earth. While plastic behaviour may set in under such high pressure conditions and analysis beyond elasticity theory may be required, sphericity may be a cause of large compressive stress in the trench strike direction. This stress may play an important role in forming the overall shape of the Earth's subduction zones.     

Phase equilibria for the Fe2SiO4–Fe3O4 system under high pressure

High pressure phase relation of the system Fe₂SiO₄–Fe₃O₄ was investigated by synthesis experiments using multi-anvil high pressure apparatus. A complete solid solution with spinel structure along Fe₂SiO₄–Fe₃O₄ join occurs above 9 GPa at 1200 °C. Lattice constants of the solid solution show almost linear variation with composition. A spinelloid phase is stable for intermediate compositions in the pressure range from 3 to 9 GPa. the synthesized spinelloid phase is successfully indexed assuming nickel aluminosilicate V type structure. Received: October 16, 1995 / Revised, accepted: March 19, 1997     

An ecological-physical coupled model applied to Station Papa

A vertical one-dimensional ecosystem model was constructed and applied to Station Papa. The model has seven compartments (phytoplankton, nitrate, ammonium, zooplankton, particulate organic matters, dissolved organic matters, dissolved oxygen) and was coupled with a mixed layer model for calculating diffusion coefficient which appears in the governing equations. The mixed layer model was driven by SST, SSS data observed at Station Papa in 1980 and ECMWF wind data for 1980, and the ecosystem model was driven by fixing nitrate concentration in deep layer to an observational value. The phytoplankton maximum in March was reproduced by the model although the maximum in fall-winter could not be reproduced. The model also suggests the importance of studying nitrification. As a whole, the model could reproduce characteristic features at Station Papa such as the summer ammonium maximum at 50 m depth, the summer dissolved oxygen maximum at 70 m depth and the absence of remarkable phytoplankton bloom.     

One Dimensional Ecosystem Model Simulation of the Effects of Vertical Dilution by the Winter Mixing on the Spring Diatom Bloom

A one-dimensional ecosystem model has been used to investigate the processes relevant to the spring diatom bloom which play important roles in the biogeochemical cycle in the western subarctic Pacific. The model represents the plankton dynamics and the nutrient cycles in the spring diatom bloom; its results show the importance of dilution by deep mixing in winter. It is supposed that the vertically integrated biomass of phytoplankton decreases in the winter due to the decrease of photosynthesis, because the deep mixing transports phytoplankton to a layer with a low light level. However, the observed integrated diatom biomass increases as the mixed layer deepens. This is because the decrease of concentration due to dilution by mixing causes the diatom grazed pressure to be less significant than diatom photosynthesis. In other words, the effect of dilution on the grazed rate is more significant than the effect on the photosynthesis rate because the grazed rate depends on the concentrations of both diatom and grazer, whereas the photosynthesis rate depends only diatom concentration. The average specific diatom grazed rate, defined as grazed rate divided by diatom biomass, decreases by 35% associated with the deepening, while the average specific photosynthesis rate of diatom decreases by 11%. As a result, the average specific net diatom growth rate during the deep mixing is about 70% of its maximum during the spring diatom bloom. The deep mixing significantly affects the amplitude of the spring diatom bloom not only by the supply of nutrients but also by the dilution which drastically decreases the grazed pressure. This revised version was published online in July 2006 with corrections to the Cover Date.     

Attempting Consistent Simulations of Stn. ALOHA with a Multi-Element Ecosystem Model

We used a one dimensional, multi-element model to simulate the primary production (PP), recycling and export of organic matter at Stn. ALOHA, near Hawaii. We compared versions of the model with and without the cycling of dissolved organic matter (DOM) via the Microbial Food Web (MFW). We incorporated recently published measurements of high C:N ratios for uptake by diazotrophs. For other phytoplankton we included a formulation for overflow production of dissolved organic carbon (DOC), which occurs under nutrient-limited, light-replete conditions. We were able to match the observed mean DOC profile near the surface with both models, by tuning only the fraction of overflow DOC that is labile. The simulated bulk C:N remineralization ratio from the MFW model agreed well with a data-based estimate for the North Pacific subtropical gyre, but that from the Base model was too low. This is because the MFW model includes bacteria, with their low-C:N biomass. Simulated mean PP was lower than observed by 10% (Base) and 27% (MFW). This is consistent with the expectation that the ¹⁴C-method measures something greater than net production. DOC accounted for approximately half of simulated PP, most of this being overflow DOC. We find that overflow production and the MFW are key processes for reconciling the various data and PP measurements at this oligotrophic site. The impact of bacteria on the C:N remineralization ratio is an important link between ecosystem structure and the cycling of carbon.     

Sedimentary Responses to an Abrupt Change of Biogenic Silica Flux by a Sediment Model for Long Timescale Simulations

A one-dimensional sediment model was developed by introducing the CIP-CSL2 scheme in the advection term. This enables us to represent fluctuations of the vertical profiles of 100 cm depth, which needs integrations on a longer timescale than 10 kyr, because this scheme avoids smoothing of the vertical profiles due to no numerical diffusion, instead of a commonly used forward scheme. Using the models with 10 and 100 cm depths (hereafter, called the 10 cm model and the 100 cm model, respectively) to predict the contents of biogenic silica and clay, we have evaluated the sedimentary responses to an abrupt change in the biogenic silica rain flux. When the rain flux of biogenic silica abruptly decreases, the 10 cm model underestimates the dissolution flux of biogenic silica compared with the 100 cm model. This is because only clay is supplied from the bottom boundary associated with excess dissolution and also because dissolution of biogenic silica below 10 cm depth is neglected by the 10 cm model. When we consider dissolved matter inputs from sediments to marine biogiochemical cycles, simulations by the 100 cm model would lead to more realistic responses in sediments.     

Interfacial energies for quartz and albite in pelitic schist

Interfacial energies of quartz/quartz (qz/qz), albite/albite (ab/ab), and quartz/albite (qz/ab) boundaries in low-grade pelitic schist were determined based on measured values of dihedral angles. Three kinds of microstructures were investigated, and the interfacial energies were obtained in two independent ways. (1) Relative values of interfacial energy were calculated from dihedral angles formed at quartz and albite triple junctions. (2) Subgrain boundary energy was calculated using the Read-Shockley theory for a boundary connected to an intergranular pore. Dihedral angles formed at the corners of intergranular pores were measured. From the interfacial tension balance equation, the value of the qz/qz grain boundary energy was then obtained. (3) Dihedral angles formed at intersections of either pericline or albite twin boundaries with either ab/ab or qz/ab boundaries were measured. The twin boundary energy was calculated based on a previously derived equation using Landau potential, twin wall thickness, and critical temperature for a phase transition in albite. With a modified interfacial tension balance equation for a twin boundary fixed to a facet orientation, the interfacial energies of ab/ab and qz/ab boundaries were obtained. Energies obtained by methods of (2) and (3) are in good agreement. The interfacial energies for qz/qz, ab/ab, and qz/ab boundaries obtained in this study are 270뀶, 300끞, and 250끀 mJ/m², respectively.