A mechanism of shape-transformation of quartz inclusions in albite of regional metamorphic rocks

The process of shape-transformation of quartz inclusions from polyhedral to spherical grains in albite single crystals during metamorphism is mainly controlled by the grain boundary diffusion of oxygen along the quartz/albite interface to reduce the interfacial free energy. The rate of the process, which is represented by the growth rate of the curvature of the edge surface of the grain, depends significantly on temperature and on the grain size of the quartz inclusion. The relations between temperature, T, the time, t_r, and the critical radius, R_c, which is equal to the radius of maximum spherical grains, are given by log R_c = −0.11E_b/RT + 0.25log t_r + C, in which E_b is the activation energy of the grain boundary diffusion of oxygen along the quartz/albite interface and C is a material constant.

The mean critical radius of spherical quartz inclusions in albite is 5 μm for the upper chlorite zone and garnet zone, 10 μm for the lower biotite zone, and 20 μm for the upper biotite zone in the Sambagawa metamorphic terrain. The mean values of the critical radii of spherical quartz inclusions in oligoclase of the Ryoke metamorphic rocks is about 5 μm for the chlorite zone and about 10–20 μm for the sillimanite zone.

Assuming temperatures of about 350°C for the upper chlorite and garnet zones, 400°C for the lower biotite zone, 550°C for the upper biotite zone, and 700°C for the sillimanite zone, the activation energy for the grain boundary diffusion of oxygen along the quartz/plagioclase interfase is estimated to be about 30 kcal/mol.     

Evolution of Mount Fuji,Japan: Inference from drilling into the subaerial oldest volcano,pre‐Komitake

A buried, old volcanic body (pre‐Komitake Volcano) was discovered during drilling into the northeastern flank of Mount Fuji. The pre‐Komitake Volcano is characterized by hornblende‐bearing andesite and dacite, in contrast to the porphyritic basaltic rocks of Komitake Volcano and to the olivine‐bearing basaltic rocks of Fuji Volcano. K‐Ar age determinations and geological analysis of drilling cores suggest that the pre‐Komitake Volcano began with effusion of basaltic lava flows around 260 ka and ended with explosive eruptions of basaltic andesite and dacite magma around 160 ka. After deposition of a thin soil layer on the pre‐Komitake volcanic rocks, successive effusions of lava flows occurred at Komitake Volcano until 100 ka. Explosive eruptions of Fuji Volcano followed shortly after the activity of Komitake. The long‐term eruption rate of about 3 km³/ka or more for Fuji Volcano is much higher than that estimated for pre‐Komitake and Komitake. The chemical variation within Fuji Volcano, represented by an increase in incompatible elements at nearly constant SiO₂, differs from that within pre‐Komitake and other volcanoes in the northern Izu‐Bonin arc, where incompatible elements increase with increasing SiO₂. These changes in the volcanism in Mount Fuji may have occurred due to a change in regional tectonics around 150 ka, although this remains unproven.     

CMT data inversion using a Bayesian information criterion to estimate seismogenic stress fields

We developed an inversion method to estimate the stress fields related to earthquake generation (seismogenic stress fields) from the centroid moment tensors (CMT) of seismic events by using Akaike's Bayesian information criterion (ABIC). On the idea that the occurrence of an earthquake releases some part of the seismogenic stress field around its hypocentre, we define the CMT of a seismic event by a weighted volume integral of the true but unknown seismogenic stress field. Representing each component of the seismogenic stress field by the superposition of a finite number of 3-D basis functions (tri-cubic B-splines), we obtain a set of linear observation equations to be solved for the expansion coefficients (model parameters). We introduce prior constraint on the roughness of the seismogenic stress field and combine it with observed data to construct a Bayesian model with hierarchic, highly flexible structure controlled by hyper-parameters. The optimum values of the hyper-parameters are objectively determined form observed data by using ABIC. Given the optimum values of the hyper-parameters, we can obtain the best estimates of model parameters by using a maximum likelihood algorithm. We tested the validity of the inversion method through numerical experiments on two synthetic CMT data sets, assuming the distribution of fault orientations to be aligned with the maximum shear stress plane in one case and to be random in the other case. Then we applied the inversion method to actual CMT data in northeast Japan, and obtained the pattern of the seismogenic stress field consistent with geophysical and geological observations.     

Metaproteomic characterization of dissolved organic matter in coastal seawater

We performed a comprehensive metaproteomic analysis of the dissolved organic matter (DOM) in Japanese coastal waters using liquid chromatography–tandem mass spectrometry and demonstrated that these proteomes were characterized by proteins with various functions, including metabolic enzymes, membranes, and photosynthetic proteins. The protein sources included cyanobacteria, heterotrophic bacteria, and eukaryotic phytoplankton. Most of the components were similar among samples and also similar to pelagic components. We also observed differences in the compositions of the microbial communities of origin among the different dissolved protein samples and differences in the relative abundance of specific dissolved protein types (e.g., cytoskeletal proteins), possibly indicating potential dynamics in the coastal DOM pool.     

Correction of atmospheric effect on ADEOS/OCTS ocean color data: Algorithm description and evaluation of its performance

This paper first describes the atmospheric correction algorithm for OCTS visible band data used at NASDA/EOC. Sharing a basic structure with Gordon and Wang’s Sea WiFS algorithm, it uses 10 candidate aerosol models including the “Asian dust model” introduced in consideration of the unique feature of aerosols over the east Asian waters. Based on the observations at 670 and 865 nm bands, the algorithm selects a pair of aerosol models that account best for the observed spectral reflectances, and synthesizes the aerosol reflectance used for the atmospheric correction. Two different schemes for determining the value of the parameter for the aerosol model selection are presented and their anticipated estimation error is analyzed in terms of retrieved water reflectance at 443 nm. The results of our numerical simulation show that the standard deviation of the estimation error of the “weighted average” scheme is mostly within the permissible level of ±0.002, reducing the error by 18% on average compared to the “simple average” scheme. The paper further discusses the expected error under the old CZCS-type atmospheric correction, which assumes constant aerosol optical properties throughout the given image. Although our algorithm has a better performance than the CZCS algorithm, further analysis shows that the error induced by the assumption taken in the algorithm that the water-leaving radiance at 670 nm band is negligibly small may be large in high pigment concentration waters, indicating the necessity for future improvements.