Imbrication, flow direction and possible source areas of the pumice-flow tuffs near Bend, Oregon, U.S.A.

Imbrication, indicating flow and source direction, occurs in three Pleistocene or upper Pliocene pumice-flow tuffs exposed in a 700-km² area on the east flank of the Cascade Range near Bend, Oregon, and shows the location of previously unknown source vents of these tuffs. The imbrication is formed by inclined elongate and/or flat pumice or lithic fragments and locally by elongate plagioclase crystals. Imbrication is best developed within the lower zones of individual flow units; the pumiceous top zones also locally show imbrication directions parallel to that in the lower zones. Moreover, the areal pattern of size distribution of lithic and pumice fragments in the flows is concordant with the flow direction pattern indicated by imbrication.The upper pumice flow shows a fan-shaped pattern of flow directions indicated by imbrication which points to a western source. A possible vent, about 20 km west of Bend in the highland near Broken Top Volcano, is marked by many silicic domes and basaltic cinder cones where there is a 6–8 mgal negative Bouguer gravity anomaly. In contrast, imbrication in the middle and lower pumice flows indicates flow from a source southwest of Bend. Vents in this direction are not obvious. Possible buried vents are located about 30 km and 45 km southwest of Bend near Sitkum Butte and Lookout Mountain, respectively.     

Water-Soluble dicarboxylic acids, ketoacids and dicarbonyls in the atmospheric aerosols over the southern ocean and western pacific ocean

Water-soluble dicarboxylic acids (DCAs), ketoacids, and α-dicarbonyls in the marine aerosol samples collected over the Southern Ocean and western Pacific Ocean were determined. Oxalic acid was the most abundant species, followed by malonic acid and then succinic acid. It is suggested that aerosol concentrations of the organics over the Southern Ocean in this work represent their global background levels. Over the Southern Ocean, total concentrations of DCAs ranged from 2.9 to 7.2 ng m⁻³ (average: 4.5 ng m⁻³), ketoacids from 0.14 to 0.40 ng m⁻³ (av.: 0.28 ng m⁻³), and dicarbonyls from 0.06 to 0.29 ng m⁻³ (av.: 0.11 ng m⁻³). Over the western Pacific, total concentrations of DCAs ranged from 1.7 to 170 ng m⁻³ (av.: 60 ng m⁻³), ketoacids from 0.08 to 5.3 ng m⁻³ (av.: 1.8 ng m⁻³), and dicarbonyls from 0.03 to 4.6 ng m⁻³ (av.: 0.95 ng m⁻³). DCAs over the western Pacific have constituted a large fraction of organic aerosols with a mean DCAs-C/TC (total carbon) of 7.0% (range: 0.59–14%). Such a high value was in contrast to the low DCAs-C/TC (av.: 1.8%; range: 0.89–4.0%) for the Southern Ocean aerosols. Based on the relative abundances and latitudinal distributions of these organics, we propose that long-range atmospheric transport is more important over the western Pacific Ocean, in contrast, in situ photochemical production is more significant over the Southern Ocean although absolute concentrations of the organics are much lower.     

Determination of three-dimensional in situ stresses from anelastic strain recovery measurement of cores at great depth

In order to examine whether the anelastic strain recovery (ASR) method can be applied for determining the in situ stress in hard rocks at great depths, the anelastic strain recovery of oriented cores was measured in six independent directions. The core specimens were taken from four depths within the range of about 2400–4500 m MD at the METI Niitsu well in Japan; the rock materials were mudstone, dolerite, basalt and andesite. For all the rocks the expansional anelastic strains were obtained, the magnitude of the strains in various directions continuously measured for 1 or 2 weeks was of the order of 1000 × 10^− 6 in mudstone; in contrast, strains of the other cores did not exceed a few hundred microstrains. These strains were used for a three-dimensional analysis of the principal in situ stresses. At the third depth, the principal stress directions were considered to be affected by fractures pre-existing near the core, and showed the features of a very local stress state. With the exception of this data, the directions determined by the ASR method were in agreement with those determined using other in situ measurement methods. Based on two assumptions, i.e., (i) the rock stress in vertical direction is equal to density-related gravitational overburden stress, (ii) the ratio of anelastic strain recovery compliance of shear deformation mode and the compliance of volumetric deformation mode is equal to 2, the values of the three principal stresses were estimated. The values of the minimum principal stress in the plane perpendicular to the well axis determined in this study were in agreement with those determined based on extended leak-off tests (ELOT or XLOT) conducted at the same well. Therefore, it can be said that the ASR method is well suited for use in directly determining the directions of principal in situ stresses in three dimensions and in estimating the magnitude of the stresses in isotropic rocks at great depths, such as those encountered when drilling deep into a submarine seismogenic zone.     

Dynamic tensile cracking in slopes possibly induced by Rayleigh surface waves

The 1978 Miyagi-ken-oki earthquake generated crack openings in fill slopes in Sendai, Japan. Although the slopes were reinforced with steel pipe piles after this event, tensile cracking was induced again in the same slopes by the 2011 off the Pacific coast of Tohoku earthquake. This clearly indicates conventional countermeasures against dynamic slope failures, usually based on body wave interaction with slopes, are not always effective. Considering that similar open cracks in the top surfaces of slopes were also found in New Zealand and California, we may need to clarify the generation mechanism of such cracks with some new insights. Recently, it has been indicated that the effect of Rayleigh surface wave propagation on dynamic slope stability may become significant. In this contribution, therefore, by performing two-dimensional elastodynamic analyses of Rayleigh/body wave interaction with a geometrically simple model slope, we show Rayleigh waves may play a more dominant role than body waves in generating the open cracks, and emphasise the importance of taking into account the dynamic slope failures induced by Rayleigh waves. Surface waves may be produced also by nearby blasting, etc., and therefore, the results obtained here may be of crucial importance in comprehending the dynamic stability of slopes in general.     

A quadratic element method for evaluating groundwater flow by the inversion of surface tilt with application to the Tono Area, Japan

By modifying a previous method with constant elements, we developed a quadratic element method for more accurately estimating groundwater flow by the inversion of tilt data. In this method: (1) a region of groundwater flow is divided into quadratic elements in which the change in groundwater volume per unit volume of rock (Δv) and the Skempton coefficient (B) vary in a quadratic manner with the coordinates, (2) the values of Δv are set to zero at the boundaries of the region of groundwater flow and (3) the sum of the squared second derivatives of Δv is adopted as a constraining condition that is weighted and added to the sum of the squared errors in tilt. First, analyses were performed for a flow model to determine the accuracy of this method for estimating groundwater flow and also to clarify the effect of the assumed size of a region of groundwater flow. These analyses showed that the quadratic element method proposed in this study gives a much better estimation of Δv than the constant element method and that a large region of groundwater flow should be assumed, rather than a small region, since the values of Δv at points outside of the actual region of groundwater flow are estimated to be nearly zero when a large region is assumed while these values are greatly overestimated when an excessively small region is assumed. Finally, the quadratic element method was applied to the site of the Mizunami Underground Research Laboratory in the Tono area, Japan. Inverse analyses were performed for tilt data measured by four tiltmeters with a resolution of 10⁻⁹ radians during the excavation of two shafts under the assumption that the rock mass is an isotropic and homogeneous half- space. The results showed that the method proposed in this study reproduced the tilt data very accurately. Thus, the distribution of Δv was estimated without sacrificing the reproducibility of the tilt data. The contour maps of B(1 + ν)Δv (ν: Poisson’s ratio) showed that the heterogeneous flow of groundwater occurred at the site and that groundwater volume decreased mainly in the area surrounded by two faults. The latter result is consistent with the finding obtained by previous investigations that these faults have low permeability in the direction perpendicular to the strike and may act as a flow barrier.     

Non-monotonic optimal damper placement via steepest direction search

An efficient and systematic procedure is proposed for finding the optimal damper positioning to minimize the dynamic compliance of a 3-D shear building model. The dynamic compliance is expressed in terms of the transfer function amplitudes of the local interstorey drifts evaluated at the undamped fundamental natural frequency. The dynamic compliance is minimized subject to a constraint on the sum of the damping coefficients of added dampers. Optimality criteria are derived and the optimal damper positioning is determined via an original steepest direction search algorithm. This algorithm enables one to find an optimal damper positioning sequentially for gradually increasing damper capacity levels. A non-monotonic design path with respect to the total damper capacity level often appears in the application of this algorithm. A new augmented algorithm via parameter switching is devised to find this non-monotonic design path. Copyright © 1999 John Wiley & Sons, Ltd.     

Escape from a Crisis in Fokker-Planck Models

Recent N-body simulations have shown that there is a serious discrepancy between the results of N-body simulations and the results of Fokker-Planck simulations for the evolution of globular and rich open clusters under the influence of the galactic tidal field. In some cases, the lifetime obtained from Fokker-Planck calculations is more than an order of magnitude smaller than those from N-body simulations. In this paper we show that the principal cause for this discrepancy is an over-simplified treatment of the tidal field used in previous Fokker-Planck simulations. We performed new Fokker-Planck calculations using a more appropriate implementation for the boundary condition of the tidal field. The implementation is only possible with anisotropic Fokker-Planck models, while all previous Fokker-Planck calculations rely on the assumption of isotropy. Our new Fokker-Planck results agree well with N-body results. Comparison of the two types of simulations gives a better understanding of the evolution of such clusters. This revised version was published online in July 2006 with corrections to the Cover Date.     

Comparison of community structures between particle-associated and free-living prokaryotes in tropical and subtropical Pacific Ocean surface waters

The subtropical and tropical regions of the Pacific Ocean are less productive than other oceanic regions. Although particle association should be an important strategy for heterotrophic prokaryotes to survive in such environments, we have little information on particle-associated (PA) prokaryotes in these regions. The specific aim of this study was to determine bacterial and archaeal community structures in the PA assemblage in comparison to the free-living (FL) assemblage in the North Pacific Subtropical Gyre, the South Pacific Subtropical Gyre, and an eastern equatorial region of the Pacific Ocean. Community profiles and phylogenetic identities were obtained by denaturing gradient gel electrophoresis, 454-pyrosequencing, and cloning followed by Sanger sequencing of 16Sr RNA gene amplicons. The distribution patterns of some abundant groups in three regions and two lifestyles (PA and FL) are shown in this study. Also, the PA community structures of bacteria differed from the FL ones and exhibited higher diversity than the FL ones, while the archaeal community structures did not show significant differences between PA and FL assemblages. We found that specific phylotypes of Gammaproteobacteria and Flavobacteria were abundant in PA bacterial assemblages, suggesting that they prefer to attach and consume particulate organic matter. In summary, the surface seawater PA assemblages represent very different bacterial and archaeal community structures between three different oceanic regions, each of which had distinct PA and FL community structures. These results imply that environmental factors determine microbial community structures.     

The characteristic displacement in rate and state-dependent friction from a micromechanical point of view

The physical meaning of the characteristic displacement that has been observed in velocity-stepping friction experiments was investigated based on the micromechanics of asperity contact. It has been empirically found for bare rock surfaces that the magnitude of the characteristic displacement is dependent only on surface roughness and insensitive to both slip velocity and normal stress. Thus the characteristic displacement has been interpreted as the displacement required to change the population of contact points completely. Here arises a question about the physical mechanism by which the contact population changes. Because individual asperity contacts form, grow and are eliminated with displacement, there are at least two possible interpretations for the characteristic displacement: (1) it is the distance over which the contacts existing at the moment of the velocity change all fade away, being replaced by new asperity contacts, or (2) it is the distance required for a complete replacement in the real contact area that existed at the moment of the velocity change. In order to test these possibilities, theoretical models were developed based on the statistics of distributed asperity summits. A computer simulation was also performed to check the validity of the theoretical models using three-dimensional surface topography data with various surface roughnesses. The deformation was assumed to be elastic at each asperity contact. The results of both the simulation and the theoretical models show that the characteristic displacement in (1) is about three times longer than that in (2). Comparison of the results with the experimental observations obtained by others indicates that the possibility (2) is the correct interpretation. This means that the state in the rate and state variable friction law is memorized in a very confined area of real contact. Further, our results explain why the characteristic displacement is insensitive to normal stress: this comes from the fact that the microscopic properties such as the mean contact diameter are insensitive to normal stress. The approach based on the micromechanics of asperity contact is useful to investigate the underlying mechanism of various phenomena in rock friction.