Instability of gradient dependent elasto-viscoplasticity for clay

A gradient-dependent viscoplastic constitutive model for water saturated clay is proposed to describe the strain localization phenomena and pattern formation during deformation. Second- and fourth-order gradients of volumetric viscoplastic strain are introduced into the constitutive equations to account for the non-local effects due to the motion of microstructures. A linear perturbation analysis is applied to this model. The instability of the government equations (i.e. the constitutive equations and the equations of motion for the clay skeleton and pore water) is discussed for both the one-dimensional and the two-dimensional situations. In addition, issues concerned with the formulation of boundary value problems by finite element analysis in relation to the formulation and the boundary conditions are presented. Copyright © 1999 John Wiley & Sons, Ltd.     

Stability Dependence and Diurnal Change of Large-Scale Turbulence Structures in the Near-Neutral Atmospheric Boundary Layer Observed from a Meteorological Tower

Based on the analysis of observations from a 213-m tall meteorological tower at Tsukuba, Japan, we have investigated the favourable conditions for the predominant existence of large-scale turbulence structures in the near-neutral atmospheric boundary layer (ABL). From the wavelet variance spectrum for the streamwise velocity component ( $u$ ) measured by a sonic anemometer-thermometer at the highest level (200 m), large-scale structures (time-scale range of 100–300 s) predominantly exist under slightly unstable and close to neutral conditions. The emergence of large-scale structures also can be related to the diurnal evolution of the ABL. The large-scale structures play an important role in the overall flow structure of the lower boundary layer. For example, $u$ velocity components at the 200-m and 50-m levels show relatively high correlation with the existence of large-scale structures. Under slightly unstable (near-neutral) conditions, a low-speed region in advance of the high-speed structure shows a positive deviation of temperature and appears as the plume structure that is forced by buoyancy in the heated lower layer. In spite of the difference in buoyancy effects between the near-neutral and unstable cases, large-scale structures are frequently observed in both cases and the same vertical correlation of $u$ components is indicated. However, the vertical wind shear is smaller in the unstable cases. On the other hand, in near-neutral cases, the transport efficiency of momentum at the higher level and the flux contribution of sweep motions are larger than those in the unstable cases.     

Ductile deformation and development of andalusite microstructures in the Hongusan area: Constraints on the metamorphism and tectonics of the Ryoke Belt

Abstract   The main mesoscopic features of the ductile deformation in the Hongusan region of the Ryoke Belt can be explained in terms of a single deformation phase, Dm . The associated stretching lineation is roughly east–west and has an associated top-to-the-west sense of shear. This information fills an important gap in our knowledge about the deformational history of the Ryoke metamorphic belt. Andalusite porphyroblasts are well-developed in the Hongusan area and their microstructure can be used to link the deformational history with the metamorphism. However, the andalusite porphyroblasts show both intertectonic and post-tectonic microstructures. Some porphyroblasts have an intertectonic core and post-tectonic rim. Analysis of associated metamorphic minerals and bulk-rock chemistry suggests that the distinct types of microstructure can be explained in terms of distinct formation reactions related to regional and contact metamorphism. The intertectonic microstructure shows the presence of a deformation pre-dating Dm , referred to as Dc . These two deformation stages are sufficient to account for all the main deformation features. The distribution of post-tectonic andalusite growth shows a very broad thermal effect around an intrusion in the area. This can help place constraints on the thermal history of the Ryoke metamorphic belt.     

Compilation method for 1 km grid data of monthly mean air temperature for quantitative assessments of climate change impacts

A new method is proposed to compile 1 km grid data of monthly mean air temperature by dynamically downscaling general circulation model (GCM) data with a regional climate model (RCM). The downscaling method used is a technique referred to as the pseudoglobal warming method to reduce GCM bias. For the grid data, RCM data were corrected with data from an existing meteorological network. The correction model for the RCM bias was developed by stepwise multiple regression analysis using the difference in the monthly mean air temperatures between the observation and RCM output as a dependent variable and the geographical factors as independent variables. Our method corrected the RCM bias from 1.69°C to 0.58°C for the month of August in the 1990s (1990–1999).     

Self‐similar solution of a plane‐strain fracture driven by a power‐law fluid

This paper analyses the problem of a hydraulically driven fracture, propagating in an impermeable, linear elastic medium. The fracture is driven by injection of an incompressible, viscous fluid with power‐law rheology and behaviour index n?0. The opening of the fracture and the internal fluid pressure are related through the elastic singular integral equation, and the flow of fluid inside the crack is modelled using the lubrication theory. Under the additional assumptions of negligible toughness and no lag between the fluid front and the crack tip, the problem is reduced to self‐similar form. A solution that describes the crack length evolution, the fracture opening, the net fluid pressure and the fluid flow rate inside the crack is presented. This self‐similar solution is obtained by expanding the fracture opening in a series of Gegenbauer polynomials, with the series coefficients calculated using a numerical minimization procedure. The influence of the fluid index n in the crack propagation is also analysed. Copyright © 2002 John Wiley & Sons, Ltd.     

Empirical green's function corrected for source effect

One of the severe problems in the semi-empirical method for the prediction of strong ground motions is that there is no objective criterion for choosing empirical Green's functions. It is undesirable that synthesized strong ground motions are affected by the source process of an earthquake whose record is adopted as an empirical Green's function. Through the analysis of strong motion accelerograms of two aftershocks of the 1983 Japan Sea earthquakes, it is found that characteristics of the accelerograms are dependent on the moment rate function derived from teleseismic observations. A procedure is presented for removing the effect of the source process from observed strong motion accelerograms. The thus obtained empirical Green's function expresses approximately the impulse response of the medium between the earthquake source and the observation site.     

CHIME age dating of monazites from metamorphic rocks and granitic rocks of the Ryoke belt in the Iwakuni area, Southwest Japan

Abstract CHIME (chemical Th-U-total Pb isochron method) ages were determined for monazite from gneisses and granitoids of the Ryoke belt in the Iwakuni area. The CHIME monazite ages are 99.6 ± 2.4, 98.9 ± 2.1 and 98.2 ± 5.7 Ma for the Ryoke gneiss, 90.7 ± 2.2, 89.7 ± 2.0 and 89.3 ± 2.2 Ma for the Tajiri Granite, 91.0 ± 3.2, 90.6 ± 3.2 and 89.9 ± 3.2 Ma for the Namera Granite, 89.3 ±3.3 and 88.6 ± 5.6 Ma for a small stock at Shimizu, and 87.3 ± 1.6 and 86.6 ± 2.1 Ma for the post-tectonic Shimokuhara Granite. The CHIME monazite ages, interpreted as the time of the first attainment at the amphibolite facies conditions for the gneisses and as the time of emplacement for the granites, respectively, agree with the field intrusive sequence. The present dating documented that the Ryoke metamorphism in the Iwakuni area reached the amphibolite facies conditions at ∼98 Ma, was complete at -87 Ma, the time of emplacement of the post-tectonic Shimokuhara Granite.     

Effect of ocean acidification on coastal phytoplankton composition and accompanying organic nitrogen production

The effect of ocean acidification, caused by the increase in pCO₂ in seawater, on phytoplankton population and on related organic nitrogen production was experimentally examined by use of a natural coastal microbial population. pCO₂ and pH were controlled by aeration with air in which pCO₂ was at the current level (control), for which ambient air was used, and with air in which pCO₂ was ??800?? and ??1200?? ppm, in 500-L culture vessels. The experiment was continued for 15?days after addition of the inorganic nutrients such as nitrate, phosphate, and silicate. During most of the experimental period, a minor increase in phytoplankton biomass was noted, probably because of low irradiance, an increase in phytoplankton biomass was observed at the end of the experiment. Flow cytometric and microscopic observations revealed that this increase was because of Chrysochromulina sp. (Haptophyceae). The growth of Chrysochromulina sp. was most obvious in the control vessel, and tended to be obscured by increasing pCO₂ (decrease in pH), indicating the possibility that ocean acidification inhibits the growth of specific phytoplankton groups, for example Chrysochromulina sp. Production of particulate organic nitrogen (PON), determined by the ¹⁵N tracer method, also diminished under acidified conditions compared with that at the current level.     

Strain localization analysis by elasto-viscoplastic softening model for frozen sand

The aim of the present paper is to numerically analyse the behaviour of frozen sand by using a viscoplastic constitutive model with strain softening. A constitutive model has been developed introducing the stress history tensor which is a functional of the stress history, with respect to a generalized time measure. It is shown that Adachi and Oka's model is applicable to the results of triaxial tests on a frozen Toyoura sand at different strain rates. First, the instability of the model is discussed within the framework of bifurcation theory. The model is then implemented into a FEM code to numerically simulate the behaviour under plane strain conditions. From the numerical results, it is revealed that the formation of shear bands is possible and the characteristics of strain localization, such as shear banding, depend on the strain rates.     

Multiple growth of garnet,sillimanite/kyanite and monazite during amphibolite facies metamorphism: implications for the P–T–t and tectonic evolution of the western Altai Range,Mongolia

Four amphibolite facies pelitic gneisses from the western Mongolian Altai Range exhibit multistage aluminosilicate formation and various chemical‐zoning patterns in garnet. Two of them contain kyanite in the matrix and sillimanite inclusions in garnet, and the others have kyanite inclusions in garnet with sillimanite or kyanite in the matrix. The Ca‐zoning patterns of the garnet are different in each rock type. U–Th–Pb monazite geochronology revealed that all rock units experienced a c. 360 Ma event, and three of them were also affected by a c. 260 Ma event. The variations in the microstructures and garnet‐zoning profiles are caused by the differences in the (i) whole‐rock chemistry, (ii) pressure conditions during garnet growth at c. 360 Ma and (iii) equilibrium temperatures at c. 260 Ma. The garnet with sillimanite inclusions records an increase in pressure at low‐P (~5.2–7.2 kbar) and moderate temperature conditions (~620–660 °C) at c. 360 Ma. The garnet with kyanite inclusions in the other rock types was also formed during an increase in pressure but at higher pressure conditions (~7.0–8.9 kbar at ~600–640 °C). The detrital zircon provenance of all the rock types is similar and is consistent with that from the sedimentary rocks in the Altai Range, suggesting that the provenance of all the rock types was a surrounding accretionary wedge. One possible scenario for the different thermal gradient is Devonian ridge subduction beneath the Altai Range, as proposed by several researchers. The subducting ridge could have supplied heat to the accretionary wedge and elevated the geotherm at c. 360 Ma. The differences in the thermal gradients that resulted in varying prograde P–T paths might be due to variations in the thermal regimes in the upper plate that were generated by the subducting ridge. The c. 260 Ma event is characterized by a relatively high‐T/P gradient (~25 °C km^?1) and may be due to collision‐related granitic activity and re‐equilibrium at middle crustal depths, which caused the variations in the aluminosilicates in the matrix between the rock units.