Preparation of a New Geological Survey of Japan Geochemical Reference Material: Coral JCp-1

A new geochemical reference material, coral Porites sp. JCp-1 has been prepared by the Geological Survey of Japan (GSJ). Provisional values for twenty one major, minor and trace elements are presented. The homogeneity tests showed that all elements studied are considered to be homogeneously distributed.     

Reaction-induced grain boundary cracking and anisotropic fluid flow during prograde devolatilization reactions within subduction zones

Devolatilization reactions during prograde metamorphism are a key control on the fluid distribution within subduction zones. Garnets in Mn-rich quartz schist within the Sanbagawa metamorphic belt of Japan are characterized by skeletal structures containing abundant quartz inclusions. Each quartz inclusion was angular-shaped, and showed random crystallographic orientations, suggesting that these quartz inclusions were trapped via grain boundary cracking during garnet growth. Such skeletal garnet within the quartz schist formed related to decarbonation reactions with a positive total volume change (?V _t > 0), whereas the euhedral garnet within the pelitic schists formed as a result of dehydration reaction with negative ?V _t values. Coupled hydrological–chemical–mechanical processes during metamorphic devolatilization reactions were investigated by a distinct element method (DEM) numerical simulation on a foliated rock that contained reactive minerals and non-reactive matrix minerals. Negative ?V _t reactions cause a decrease in fluid pressure and do not produce fractures within the matrix. In contrast, a fluid pressure increase by positive ?V _t reactions results in hydrofracturing of the matrix. This fracturing preferentially occurs along grain boundaries and causes episodic fluid pulses associated with the development of the fracture network. The precipitation of garnet within grain boundary fractures could explain the formation of the skeletal garnet. Our DEM model also suggests a strong influence of reaction-induced fracturing on anisotropic fluid flow, meaning that dominant fluid flow directions could easily change in response to changes in stress configuration and the magnitude of differential stress during prograde metamorphism within a subduction zone.     

Periodic Intrusion of Cold Water from the Pacific Ocean into the Bottom Layer of the Bungo Channel in Japan

A cold-water intrusion, called a “bottom intrusion”, occurs in the lower layer of the Bungo Channel in Japan. It is an intrusion from the shelf slope region of the Pacific Ocean margin in the south of the channel. In order to reveal the fundamental characteristics of the bottom intrusion, we conducted long-term observations of water temperature at the surface and bottom layers of the channel and 15-day current observations at the bottom of the shelf-break region. The long-term water temperature data indicated that the bottom intrusion occurs repeatedly between early summer and late autumn, and its reiteration between early and mid-summer causes a local minimum of water temperature in the lower layer in mid-summer. Moreover, the data revealed that most of the bottom intrusions occurred in neap tidal periods. The current meter recorded a bottom intrusion with a speed of approximately 15 cm⋅s⁻¹. The current meter also revealed that the intruded cold water slowly retreated back to the shelf slope region after the intrusion. This revised version was published online in July 2006 with corrections to the Cover Date.     

A theory of semigeostrophic gravity waves and its application to the intrusion of a density current along a coast

Semigeostrophic gravity waves associated with a coastal boundary current, which has finite and uniform potential vorticity and is bounded away from the coastline by a density front on the ocean surface, are investigated. It is shown that the semigeostrophic coastal current has two waves which are named here the Semigeostrophic Coastal Wave (SCW) and the Semigeostrophic Frontal Wave (SFW). The SCW becomes an elementary Kelvin wave at some limit while the SFW is caused by the existence of the surface density front. The SCW appears mainly as variations in the upper layer depth at the coast and as alongshore velocity at the density front. On the other hand, the SFW appears mainly as variations in the width of the current. When the weak nonlinearity and ageostrophic effect are included, these semigeostrophic gravity waves satisfy the Kortweg- de Vries equation, which suggests that the local changes in the width and/or velocity of the semigeostrophic coastal current propagate as wave-like disturbances.     

K–Ar ages of the Akan-Shiretoko volcanic chain lying oblique to the Kurile trench: Implications for tectonic control of volcanism

The Akan‐Shiretoko volcanic chain, situated in the Southwestern Kurile arc, consists mainly of nine subaerial andesitic stratovolcanoes and three calderas. The chain extends in a SW–NE direction for 200 km, situated oblique to the Kurile trench at an angle of 25 degrees. Thirty‐seven new K–Ar ages, plus previous data, suggest that volcanic activity along the Akan‐Shiretoko volcanic chain began at ca 4 Ma at Akan, at the southwestern end of the chain, and systematically progressed northeastward, resulting in the southwest‐northeast‐trending volcanic chain. This spatial and temporal distribution of volcanoes can be explained by anticline development advancing northeastward from the Akan area, accompanied by magma rising through northeast‐trending fractures that developed along the anticlinal axis. The northeastward development of the anticline caused uplifting of the Akan‐Shiretoko area and changed the area from submarine to subaerial conditions. Anticline formation was likely due to deformation of the southwestern Kurile arc, with southwestward migration of the Kurile forearc sliver caused by oblique subduction of the Pacific plate. The echelon topographic arrangement of the Shiretoko, Kunashiri, Etorofu and Urup was formed at ca 1 Ma.     

Submarine slope–fan sedimentation in an ancient forearc related to contemporaneous magmatism: The Upper Cretaceous Izumi Group,southwestern Japan

The Izumi Group in southwestern Japan is considered to represent deposits in a forearc basin along an active volcanic arc during the late Late Cretaceous. The group consists mainly of felsic volcanic and plutonic detritus, and overlies a Lower to Upper Cretaceous plutono‐metamorphic complex (the Ryoke complex). In order to reconstruct the depositional environments and constrain the age of deposition, sedimentary facies and U–Pb dating of zircon grains in tuff were studied for a drilled core obtained from the basal part of the Izumi Group. On the basis of the lithofacies associations, the core was subdivided into six units from base to top, as follows: mudstone‐dominated unit nonconformably deposited on the Ryoke granodiorite; tuffaceous mudstone‐dominated unit; tuff unit; tuffaceous sandstone–mudstone unit; sandstone–mudstone unit; and sandstone‐dominated unit. This succession suggests that the depositional system changed from non‐volcanic muddy slope or basin floor, to volcaniclastic sandy submarine fan. Based on a review of published radiometric age data of the surrounding region of the Ryoke complex and the Sanyo Belt which was an active volcanic front during deposition of the Izumi Group, the U–Pb age (82.7 ±0.5 Ma) of zircon grains in the tuff unit corresponds to those of felsic volcanic and pyroclastic rocks in the Sanyo Belt.     

Preserved paleo-oceanic plateaus in accretionary complexes: Implications for the contributions of the Pacific superplume to global environmental change

We have reinvestigated the mid-Cretaceous plume pulse in relation to paleo-oceanic plateaus from accretionary prisms in the circum-Pacific region, and we have correlated the Pacific superplume activity with catastrophic environmental changes since the Neoproterozoic. The Paleo-oceanic plateaus are dated at 75–150 Ma; they were generated in the Pacific superplume region and are preserved in accretionary prisms. The volcanic edifice composed of both modern and paleo-oceanic plateaus is up to 10.7 × 10⁶ km² in area and 19.1 × 10⁷ km³ in volume. The degassing rate of CO₂ (0.82 − 1.1 × 10¹⁸ mol/m.y.) suggests a significant impact on Cretaceous global warming. The synchronous occurrence of paleo-oceanic plateaus in accretionary complexes indicates that Pacific superplume pulse activities roughly coincided at the Permo-Triassic boundary and the Vendian–Cambrian boundary interval. The CO₂ expelled by the Pacific superplume probably contributed to environmental catastrophes. The initiation of the Pacific superplume contributed to the snowball Earth event near the Vendian–Cambrian boundary; this was one of the most dramatic events in Earth's history. The scale of the Pacific superplume activity roughly corresponds to the scale of drastic environmental change.     

Methods for detecting apparent differences between spatial tessellations at different time points

This paper proposes methods for detecting apparent differences between spatial tessellations at two different points in time, with the objective of conflation of spatial tessellations at multiple time points. The methods comprise three steps. First, we eliminate systematic differences between tessellations using the affine transformation. Second, we match subregions between tessellations at two points in time and match boundaries based on matching relationships between the subregions. Third, we propose a distance metric for measuring differences between the matched boundaries and a method for determining whether the measured differences are apparent or not. We apply the proposed methods to a part of the US Census data for 1990 and 2000 and empirically demonstrate the effectiveness of these methods.     

Spatial Patterns in Nutrient and In Vivo Fluorescence Distributions in the Marginal Ice Zone and the Seasonally Open Oceanic Zone in the Indian Sector of the Antarctic Ocean, in Austral Summer

Surface temperature, salinity, concentrations of silicate (Si) and nitrate + nitrite (N), and in vivo fluorescence (Fluor) were investigated in the marginal ice zone (MIZ) and the seasonally open oceanic zone (SOOZ) (32–40°E, 64–69°S) from February 23 to 28 1992. In the MIZ the mean Si and N were 67.8 ± 2.2 M and 32.5 ± 1.7 M, respectively. There was a trend that low N values coincided with high Fluor values. Observation conducted at one point (64°S, 38°E) revealed a diel variation pattern in Fluor. Applying this pattern of deviation from noon value, all Fluor data were normalized to value at local noon. In the MIZ a significant negative correlation was observed between the normalized Fluor and N but not Si. On the other hand, Si decreased continuously from south to north in the SOOZ and was negatively correlated with the normalized Fluor. Difference in Si concentration was about 30 M between the sea around 64°S and the MIZ, while the difference in N concentration was estimated as less than 10 M. If diatoms take up silicate and nitrogen at an approximate ratio of 1:1, additional nitrogenous nutrients other than nitrate and nitrite (e.g. ammonia, urea etc.) would be required. In this case, an f-ratio of lower than 33% is obtained. It is suggested that in the MIZ abundance of phytoplankton community dominated by non-diatom increases utilizing nitrate while in the SOOZ abundance of phytoplankton community dominated by diatoms increases consuming Si and regenerated nitrogen.     

Genetic algorithm for fitting a mixed Bingham distribution to 3D orientations: a tool for the statistical and paleostress analyses of fracture orientations

The clustering of fracture orientations is important for tectonic studies and for geotechnical engineering. In this study, a real‐coded genetic algorithm was adopted to fitting a mixed Bingham distribution to orientation data by maximizing the log‐likelihood function of the distribution. The maximization is a difficult problem, because the function has multimodality and singularity. It was found that the algorithm was effective for this problem. Given the orientations of dilational fractures, the present method determines not only the stress axes and stress ratio of each of the fracture groups but also the maximum non‐dimensionalized fluid pressure at the time of their formation. In addition, the software calculates the 95 % error ellipses of the concentration axes. The present method found that the orientations of ore veins of the Akenobe Mine, SW Japan, should be partitioned into three clusters. It is shown that two of the groups had distinctive Zn and Sn contents, and that the ore fluids had overpressures only slightly greater than the minimum principal stress at the time of the deposition of Zn‐ and Sn‐rich veins.