Diffuse Emission of CO2 from Showa-Shinzan, Hokkaido, Japan: A Sign of Volcanic Dome Degassing

Two soil CO₂ efflux surveys were carried out in September 1999 and June 2002 to study the spatial distribution of diffuse CO₂ degassing and estimate the total CO₂ output from Showa-Shinzan volcanic dome, Japan. Seventy-six and 81 measurements of CO₂ efflux were performed in 1999 and 2002, respectively, covering most of Showa-Shinzan volcano. Soil CO₂ efflux data showed a wide range of values up to 552 g m^-2 d^-1. Carbon isotope signatures of the soil CO₂ ranged from -0.9‰ to -30.9‰, suggesting a mixing between different carbon reservoirs. Most of the study area showed CO₂ efflux background values during the 1999 and 2002 surveys (B = 8.2 and 4.4 g m^-2 d^-1, respectively). The spatial distribution of CO₂ efflux anomalies for both surveys showed a good correlation with the soil temperature, indicating a similar origin for the extensive soil degassing generated by condensation processes and fluids discharged by the fumarolic system of Showa-Shinzan. The total diffuse CO₂ output of Showa-Shinzan was estimated to be about 14.0–15.6 t d^-1 of CO₂ for an area of 0.53 km².     

The Role of History-Dependent Rheology in Plate Boundary Lubrication for Generating One-Sided Subduction

We have developed a two-dimensional dynamical model of asymmetric subduction integrated into the mantle convection without imposed plate velocities. In this model we consider that weak oceanic crust behaves as a lubricator on the thrust fault at the plate boundary. We introduce a rheological layer that depends on the history of the past fracture to simulate the effect of the oceanic crust. The thickness of this layer is set to be as thin as the Earth's oceanic crust. To treat 1-kilometer scale structure at the plate boundary in the 1000-kilometer scale mantle convection calculation, we introduce a new numerical method to solve the hydrodynamic equations using a couple of uniform and nonuniform grids of control volumes. Using our developed models, we have systematically investigated effects of basic rheological parameters that determine the deformation strength of the lithosphere and the oceanic crust on the development of the subducted slab, with a focus on the plate motion controlling mechanism. In our model the plate subduction is produced when the friction coefficient (0.004–0.008) of the modeled oceanic crust and the maximum strength (400 MPa) of the lithosphere are in plausible range inferred from the observations on the plate driving forces and the plate deformation, and the rheology experiments. In this range of the plate strength, yielding induces the plate bending. In this case the speed of plate motion is controlled more by viscosity layering of the underlying mantle than by the plate strength. To examine the setting of the overriding plate, we also consider the two end-member cases in which the overriding plate is fixed or freely-movable. In the case of the freely-movable overriding plate, the trench motion considerably changes the dip angle of the deep slab. Especially in the case with a shallow-angle plate boundary, retrograde slab motion occurs to generate a shallow-angle deep slab.     

Orbital motion of a massive object in a dense stellar system

In this paper we show the positional oscillation of a massive object in a dense stellar system by numerical N -body simulations. We found that the central massive object, which at first is placed at rest at the centre of the surrounding spherical stellar system, promptly departs from the centre and rotates in accordance with the rotation of the stellar system, if the stellar system has an appreciable rotation. This oscillatory motion continues for a long time because of the absence of dynamical friction. Such a long-lasting oscillation may explain the asymmetric structure observed in the centres of M31 and NGC 4486B, may cause the secular flow of gaseous elements distributed in the central regions of galaxies on to the massive object, and may ignite activity in the centres of galaxies.     

Global shearing modes of galactic disks

A new scheme has been devised to calculate discrete unstable global shearing spiral modes of gaseous disk models of galaxies. The scheme makes use of the Legendre expansion of eigen functions and the problem of stability analysis is reduced to an eigenvalue problem of an infinite matrix. The spiral patterns of these shearing wave solutions of linearized equations change their form in the course of time due to the differential rotation of the equilibrium disk. These shearing wave solutions are presumed to have intermediate characteristics between so-called density-waves and material arms. Comparison between these shearing modes and the non-shearing normal modes for a series of disk models is presented.     

A new, automatic hydrothermal fluid sampler using a shape-memory alloy

A new hydrothermal fluid sampler has been developed to provide more maneuverability in underwater operation. The sampler characteristically employs a shape-memory alloy, which senses high temperature and actuates the suction mechanism. A shape-memory alloy is also used to switch the intake valve of the sampler, the intention being to avoid missampling when the inlet is in low temperature water. Prototype samplers were designed to collect the fluids hotter than 80°C. Test sampling was performed at hydrothermal vents (1372–1374 m deep) in the submarine volcano, Suiyo Seamount, Izu-Bonin Arc, northwestern Pacific. Observed fluid temperature was between 138 and 298°C, while the ambient seawater temperature was 3.1°C. Each prototype collected about 100 ml fluid as designed. The magnesium concentration in the samples indicated a seawater content of 47.5–90.8%, which indicates the entrainment of ambient seawater. Microscopic observatinn revealed the occurrence of microorganisms in the sample fluids at a population density of 10⁵ to 10⁶ cells ml⁻¹, which is 2–3 orders of magnitude higher than those in seawater at that depth. The use of the newly developed fluid sampler will greatly facilitate the collection of vent-associated microorganisms, which are of potential biological and biotechnological interest.     

Volume transports proceeding to the Kuroshio Extension region and recirculating in the Shikoku Basin

The volume transport of the Kuroshio, the western boundary current of the North Pacific subtropical gyre, varies vigorously due to merging of disturbances propagating from the entire North Pacific. Taking into account the recirculation in the Shikoku Basin by the zonal observation line at 30°N to the west of the Izu–Ogasawara Ridge, we estimated the volume transport in the top 1,000 m layer toward the Kuroshio Extension region. The volume transport of the local recirculation gyre in the Shikoku Basin increases associated with the westward extension of the gyre, particularly in the period of the large meandering path of the Kuroshio south of Japan. Meanwhile, most of the transport variations toward the Kuroshio Extension region correspond to those of the Kuroshio transport on the continental slope south of Japan, which vary independently of those of the recirculation gyre.     

Analysis of formation pressure test results in the Mount Elbert methane hydrate reservoir through numerical simulation

Targeting the methane hydrate (MH) bearing units C and D at the Mount Elbert prospect on the Alaska North Slope, four MDT (Modular Dynamic Formation Tester) tests were conducted in February 2007. The C2 MDT test was selected for history matching simulation in the MH Simulator Code Comparison Study. Through history matching simulation, the physical and chemical properties of the unit C were adjusted, which suggested the most likely reservoir properties of this unit. Based on these properties thus tuned, the numerical models replicating “Mount Elbert C2 zone like reservoir”, “PBU L-Pad like reservoir” and “PBU L-Pad down dip like reservoir” were constructed. The long term production performances of wells in these reservoirs were then forecasted assuming the MH dissociation and production by the methods of depressurization, combination of depressurization and wellbore heating, and hot water huff and puff. The predicted cumulative gas production ranges from 2.16 × 10⁶ m³/well to 8.22 × 10⁸ m³/well depending mainly on the initial temperature of the reservoir and on the production method.This paper describes the details of modeling and history matching simulation. This paper also presents the results of the examinations on the effects of reservoir properties on MH dissociation and production performances under the application of the depressurization and thermal methods.     

Radiolaria‐dated Lower Permian clastic‐rock sequence in the Fukuji area of the Hida‐gaien terrane,central Japan,and its inter‐terrane correlation across Southwest Japan

The stratigraphy and radiolarian age of the Mizuyagadani Formation in the Fukuji area of the Hida‐gaien terrane, central Japan, represent those of Lower Permian clastic‐rock sequences of the Paleozoic non‐accretionary‐wedge terranes of Southwest Japan that formed in island arc–forearc/back‐arc basin settings. The Mizuyagadani Formation consists of calcareous clastic rocks, felsic tuff, tuffaceous sandstone, tuffaceous mudstone, sandstone, mudstone, conglomerate, and lenticular limestone. Two distinctive radiolarian faunas that are newly reported from the Lower Member correspond to the zonal faunas of the Pseudoalbaillella u‐forma morphotype I assemblage zone to the Pseudoalbaillella lomentaria range zone (Asselian to Sakmarian) and the Albaillella sinuata range zone (Kungurian). In spite of a previous interpretation that the Mizuyagadani Formation is of late Middle Permian age, it consists of Asselian to Kungurian tuffaceous clastic strata in its lower part and is conformably overlain by the Middle Permian Sorayama Formation. An inter‐terrane correlation of the Mizuyagadani Formation with Lower Permian tuffaceous clastic strata in the Kurosegawa terrane and the Nagato tectonic zone of Southwest Japan indicates the presence of an extensive Early Permian magmatic arc(s) that involved almost all of the Paleozoic non‐accretionary‐wedge terranes in Japan. These new biostratigraphic data provide the key to understanding the original relationships among highly disrupted Paleozoic terranes in Japan and northeast Asia.     

Factors affecting dominant peak-flow runoff-generation mechanisms among five neighbouring granitic headwater catchments

We examined how and why dominant peak-flow runoff-generation mechanisms differ among neighbouring headwater catchments. We monitored runoff and groundwater levels and performed terrain analyses in a granitic second-order catchment and its four neighbouring subcatchments in the Kiryu Experimental Watershed in Japan. Our analysis of lag times from peak rainfall to peak runoff suggests differences in the dominant peak-flow runoff-generation mechanisms among the five catchments. For two of the three zero-order catchments, with few perennial groundwater bodies, subsurface flow from hillslopes was the dominant mechanism at some events. However, the dominant mechanisms were channel precipitation and riparian runoff at almost all events in first- and second-order catchments and in the third zero-order catchment, which has a large perennial groundwater body over a bedrock depression in the riparian zone. In this zero-order catchment, the quick-flow ratio was the smallest of the five catchments because subsurface flow from the hillslope was buffered at the riparian zone. These facts suggest that the channel length, riparian buffering, and hillslope connectivity were the factors governing the different dominant peak-flow runoff-generation mechanisms among the catchments. Riparian buffering was affected, not only by surface topography, but also by bedrock topography and bedrock groundwater (BGW) dynamics. Our findings indicate that both of BGW dynamics and topography are important for catchment classification, and the relative importance of topography increases with the change from baseflow to stormflow. Furthermore, mismatching between a geographic source and a flow path resulted in different catchment classifications depending on the approach. Therefore, multiple approaches during both baseflow and stormflow periods are necessary for catchment classification to apply information obtained from one headwater catchment to other headwater catchments within the same region.     

Interactive visualization of 3D mantle convection

The current availability of thousands of processors at many high performance computing centers has made it feasible to carry out, in near real time, interactive visualization of 3D mantle convection temperature fields, using grid configurations having 10–100 million unknowns. We will describe the technical details involved in carrying out this endeavor, using the facilities available at the Laboratory of Computational Science and Engineering (LCSE) at the University of Minnesota. These technical details involve the modification of a parallel mantle convection program, ACuTEMan; the usage of client–server socket based programs to transfer upwards of a terabyte of time series scientific model data using a local network; a rendering system containing multiple nodes; a high resolution PowerWall display, and the interactive visualization software, DSCVR. We have found that working in an interactive visualizastion mode allows for fast and efficient analysis of mantle convection results. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.