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.     

Radial and rotational velocities of young brown dwarfs and very low-mass stars in the Upper Scorpius OB association and the ρ Ophiuchi cloud core

We present the results of a radial velocity (RV) survey of 14 brown dwarfs (BDs) and very low-mass (VLM) stars in the Upper Scorpius OB association (UScoOB) and three BD candidates in the ρ Ophiuchi dark cloud core. We obtained high-resolution echelle spectra at the Very Large Telescope using the Ultraviolet and Visual Echelle Spectrograph (UVES) at two different epochs for each object, and measured the shifts in their RVs to identify candidates for binary/multiple systems in the sample. The average time separation of the RV measurements is 21.6 d, and our survey is sensitive to the binaries with separation <0.1 au. We found that four out of 17 objects (or 24⁺¹⁶₋₁₃ per cent by fraction) show a significant RV change in 4–33 d time-scale, and are considered as binary/multiple 'candidates'. We found no double-lined spectroscopic binaries in our sample, based on the shape of cross-correlation curves. The RV dispersion of the objects in UScoOB is found to be very similar to that of the BD and VLM stars in Chamaeleon I (Cha I). We also found the distribution of the mean rotational velocities ( v   sin   i ) of the UScoOB objects is similar to that of the Cha I, but the dispersion of v   sin   i is much larger than that of the Cha I objects.     

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.     

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.