The Holocene development of Kushu Lake on Rebun Island in Hokkaido, Japan

Diatom assemblages and sulfur content in sediments were analyzed to clarify changes in the sedimentary environment of Kushu Lake, a coastal lake on Rebun Island in Hokkaido, Japan. Salinity variations were assessed by means of a diatom-based index of paleosalinity and the sedimentary sulfur content. This paper discusses the Holocene development of the lake, in relation to Holocene relative sea-level change. For paleoenvironmental interpretation of the lake development, the rationale of the threshold method (Anundsen et al., 1994) was applied.At ca. 8000 yr BP, a coastal embayment (paleo-Kushu Bay) resulted from marine ingression. The threshold elevation at the mouth of the paleo-Kushu Bay kept pace with the rising sea-level, resulting in its enclosure at the culmination of Holocene marine transgression (ca. 6500–5000 yr BP). From predicted relative sea-level at ca. 6000 yr BP for Rebun Island (Nakada et al., 1991), the threshold may have been at least above –3 to –5 m altitude. A freshwater lake environment with strongly anoxic bottom conditions may have occurred from ca. 5500 to 5100 yr BP. After an important episode of marine ingression, the lake was isolated completely from the open sea at ca. 4900 yr BP. The diatom record suggests that the maximum lacustrine extent occurred at ca. 4900–3100 yr BP. Thereafter, water depth decreased at the lake margins.In Kushu Lake, the threshold elevation, due to a build-up of a coastal barrier, prevents us from determining the amplitude of sea-level changes, even though the age of isolation contacts corresponds to periods of regression and climatic deterioration. In spite of isostatic subsidence, the effective protection provided by the well-developed barrier did not allow registration of any relative sea-level fluctuations since its isolation.     

Crustal movement along a collision boundary of plates (case of eastern Taiwan)

Recent crustal movements in eastern Taiwan were investigated based on geodetic and seismological data as a case study of the tectonics along a colliding plate boundary. Secular horizontal strain along the Longitudinal Valley, a probable arc-continent collision boundary between the Philippine Sea and Eurasion plates, is compressive with a large rate of 2 microstrains/year, and the maximum compression axis trends in NW-SE direction. This means that nearly two-thirds of the convergence rate between the two plates is consumed in the narrow zone along the valley and confirms the view that the valley is the plate boundary between the two. The amount of uplift in the Central Range, deduced from trigonometric leveling, reaches more than one meter during the past 60 years. Such active orogeny is caused by the high strain rate in this area. Horizontal displacement of the Coastal Range, a leading edge of the Philippine Sea plate, has a left-lateral component, suggesting that the boundary is a region of oblique collision and the range moves toward the north, producing subduction of the Philippine Sea plate off the Coast of northeastern Taiwan.An extensional strain pattern is found northeastern Taiwan near Hualien, in northern Taiwan near Ilan, and in southeastern Taiwan south Taitung. In northeastern Taiwan, the direction of the maximum extension axis is nearly parallel to the isodipth contours of deep-focus earthquakes that occurred under the southern part of the Ryukyu Trench. From a model calculation using a finite-element method, this pattern of strain is interpreted as an edge effect of the collision-subduction junction in northeastern Taiwan. The direction of the maximum extension in northern Taiwan forms high angles with the direction of the tectonic line near Ilan stretching from the Okinawa Trough. This fact is harmonious with the view that the Okinawa Trough is spreading. Extensional strain observed in southeastern Taiwan may be related to the eastward subduction since the area is situated on a transition zone, from collision to subduction.     

Temporal change in scattering and attenuation associated with the earthquake occurrence—A review of recent studies on coda waves

The study of coda waves has recently attracted increasing attention from seismologists. This is due to the fact that it is viewed as a new means by which the stress accumulation stage preceding a large earthquake can be measured, since the scattering paths nearly uniformly cover a fairly large region around the focus and observation stations, compared with the direct ray paths. To date, we have had many reports on the temporal variation of the relation between coda duration and amplitude magnitude, and that of the coda attenuationQ _c ^–1 which is estimated from coda amplitude decay. Some of these have shown a precursor-like behavior; however, others seem to have shown a coseismic change. We have critically reviewed these reports, and discussed what these observational facts tell us about the change in the heterogeneous crust. We found significant temporal variations, not only in the mean but also in the scatter ofQ _c ^–1 , associated with the mainshock occurrence. The formation of new cracks, the reopening and growing of existing cracks, the interaction of these cracks, and the pore water movement through these cracks might correspond to such variations. In addition, we may expect an inhomogeneous distribution of crack clusters in a fairly large region, compared with the aftershock region. The gradual appearance of such crack clusters seems to be the most plausible mechanism by which coda decay gradients are caused to largely scatter in the stress accumulation stage.     

Wind fluctuations near a cold vortex-tropopause funnel system observed by the MU radar

Vertical and temporal variations of three-dimensional wind velocity associated with an upper-tropospheric cold vortex-tropopause funnel system were observed by an MST radar in Japan (the MU radar). Marked changes of vertical velocity and horizontal wind direction were found between the inside and outside of the cold vortex. The vertical velocity activity outside the vortex was asymmetric; it was most active in a sector before the vortex. Unsaturated internal gravity waves in their generation stage contribute predominantly to the vertical velocity activity, suggesting that tropospheric occluded cyclones may be a possible source of middle-atmospheric gravity waves through the geostrophic adjustment process.     

Internal gravity wave selection in the upper troposphere and lower stratosphere observed by the MU radar: Preliminary results

Marked wavelike variations of the lower stratospheric wind observed on 7–10 May, 1985 by an MST radar in Japan (by the MU radar) are analyzed assuming that they are induced by monochromatic internal inertio-gravity waves. These variations are mainly composed of two modes (periods: 22 and 24 hours), both of which have zonal phase velocities (C _X ) slower than the mean westerly wind (). A statistical analysis of the zonal phase velocity shows thatC _X above andC _X below the tropopause jet stream, which is considered to be a vivid proof of wave selection due to the tropospheric mean flow and upward wave emission from the tropopause jet. A comparison between the MU radar results and routine meteorological observations leads to the conclusion that the marked waves appear when the jet stream takes a maximum wind speed.     

Volcanic and structural history of the Hohi volcanic zone,central Kyushu,Japan

More than 5000 km³ of magmatic material was erupted in Pliocene-Pleistocene times in a volcano-tectonic depression, i. e., the Hohi volcanic zone (HVZ) in central Kyushu, Japan. The eruptive deposits consist mainly of andesite lava flows and large-scale pyroclastic-flow deposits. Their eruptions were accompanied by the formation of an EW-oriented graben (70 km × 45 km) under regional NS extensional stress. Pre-Tertiary basement rocks are absent on the surface of the graben but occur at depth, having subsided up to 3 km. Radiometric ages of volcanic rocks on the surface show zoned isochrons from 5 Ma at the margin to 0.3 Ma in the center of the HVZ. The youngest center of age zonation coincides with a 30 mgal negative Bouguer gravity anomaly. Radiometric ages of rocks from drill cores are older toward the bottom of the graben, reaching a maximum of at least 4 Ma. Volcanic activity concentrated over time toward the center of the graben and buried successively erupted material. Areas of active volcanism in the HVZ became smaller and changed in style during the 5-Ma history of activity. Volcanism of the early stage (5-2 Ma) was characterized by voluminous eruptions of andesitic lava flows that formed lava plateaus and were intruded by EW-oriented feeder dikes, perhaps related to fissure eruptions. In contrast, late-stage volcanism (2-0 Ma) resulted primarily in andesitic to dacitic lava domes with features of monogenetic volcanoes produced at low eruption rates. The HVZ shows unimodal volcanism dominated by andesitic and dacitic lavas with a small amount of rhyolite and only traces of basalt; these characteristics differ from those that typify volcanism in most other extensional areas. Erupted material in the HVZ is of the calc-alkali and high-alkali tholeiite series and shows no significant chemical changes over 5 Ma, except for an increase in K₂O after 1.6 Ma. The net horizontal displacement along normal faults indicates that the HVZ widened by about 10%–20% across the graben at an average rate of 0.1 cm/yr. I interpret the HVZ to be neither a pull-apart structure of the pre-Tertiary basement nor the result of propagation of the Okinawa Trough, but rather the earliest stage of rifting when vertical subsidence caused by normal faulting is compensated by filling with volcanic material.     

Low frequency electrical impedance of partially molten gabbro: The effect of melt geometry on electrical properties

Electrical impedance of partially molten gabbro was measured as a function of the frequency of the applied electric field. The impedance of the partially molten gabbro was mainly controlled by the volume and geometry of partial melt, which were observed with an optical microscope study of the quenched sample. The experimental results are interpreted by a theory that formulates the electrical property of rocks containing partial melt, taking ionic diffusion in the melt into account. The theory evaluates the contributions from two types of melt-geometry distributed in the rock matrix, namely:

1. (a) nearly isolated melt pockets, and

2. (b) connected melt along grain boundaries.

According to the theory, connected melt mainly contributes to the increase of net conductivity, while melt pockets increase displacement current and cause the frequency dependence of impedance. The characteristic ionic diffusivity in melt pockets was determined from the peak frequency of the imaginary part of impedance. The magnitude of diffusion coefficient is estimated to be of the order of 10⁻⁷ to 10⁻⁶ cm²/sec. This suggests that sodium, iron, magnesium and/or calcium ions mainly carry the electric charge in partial melt. The observed frequency dependence of impedance is considered to reflect ionic polarization in melt pockets due to the migration of these alkali ions toward solid-melt interfaces.The present experimental and theoretical work suggests that electrical properties of the earth's upper mantle may involve an intrinsic frequency dependence associated with ionic diffusion in the frequency range of geomagnetic observation, if a partial melt zone exists. The analysis of the earth's conductivity structure may be inadequate if such intrinsic frequency dependence of mantle media is neglected.     

Early Paleozoic Structural and Metamorphic Evolution of Western Sierra de San Luis (Argentina), in Relation to Cuyania Accretion

Structural, metamorphic and isotopic data obtained from the Nogoli Metamorphic Complex of western Sierra de San Luis indicate that the Early Paleozoic Famatinian Orogeny overprinted an already structured and metamorphosed older basement. The older geological features are relict NW trending fabric associated with high-grade (amphibolite facies) regional metamorphism preserved within thin strips of schists and paragneisses and in the core of mafic to ultramafic lenses. Arc magmatism, medium P (Barrovian type)/high T (amphibolite to granulite facies) regional metamorphism and penetrative NNE to NE trending foliation are related to the building of the Famatinian orogenic belt. The P-T conditions of the Famatinian prograde metamorphism reached a pressure peak of ca. 8 kb, with a thermal peak from -750°C up to -820°C. U-Pb conventional and chemical dating and Ar-Ar plateau ages constrain the peak of the main orogenic phase related to the Famatinian belt to 470–457 Ma (Early to Mid-Ordovician). Greenschist facies retrograde metamorphism closely associated with shear zones and secondary Ar-Ar plateau and Sm-Nd ages suggest that a late to post-orogenic phase of the Famatinian belt was active at least since -445 Ma. This phase continued during the Silurian to Late Devonian times through multiple reactivation of early shear zones. The Famatinian Orogeny reset a previous thermal history and therefore, the timing of the relict fabric could not be constrained conclusively with radiometric dates. Despite this difficulty, a range of 520 to 490 Ma suggests some inheritance from Pampean events registered by the older NW-SE fabric. The Early to Mid-Ordovician regional metamorphism and ductile deformation of the western Sierra de San Luis is interpreted as the orogenic effects of the collision of the allochthonous Cuyania terrane with the autochthonous proto-Pacific margin of Gondwana during the Famatinian Orogeny.     

A practical numerical method for large strain liquefaction analysis of saturated soils

Accurate prediction of the liquefaction of saturated soils is based on strong coupling between the pore fluid phase and soil skeleton. A practical numerical method for large strain dynamic analysis of saturated soils is presented. The u–p formulation is used for the governing equations that describe the coupled problem in terms of soil skeleton displacement and excess pore pressure. A mixed finite element and finite difference scheme related to large strain analysis of saturated soils based on the updated Lagrangian method is given. The equilibrium equation of fluid-saturated soils is spatially discretized by the finite element method, whereas terms associated with excess pore pressure in the continuity equation are spatially discretized by the finite difference method. An effective cyclic elasto-plastic constitutive model is adopted to simulate the non-linear behavior of saturated soils under dynamic loading. Several numerical examples that include a saturated soil column and caisson-type quay wall are presented to verify the accuracy of the method and its usefulness and applicability to solutions of large strain liquefaction analysis of saturated soils in practical problems.