Finite-element Simulation of Seismic Ground Motion with a Voxel Mesh

— Accurate simulation of seismic ground motion for three-dimensionally complex topography and structures is one of the most important goals of strong motion seismology. The finite-element method (FEM) is well suited for this kind of simulation, since traction-free conditions are already included in the formulation, and the Courant condition is less strict than for the finite-difference method (FDM). However, the FEM usually requires both large memory and computation time. These limitations can be overcome by using a mesh consisting of voxels (rectangular prisms) with isotropy built into the explicit formulation of the dynamic matrix equation. Since operators in the voxel FEM are the combinations of ordinary FDM operators and additional terms, the method keeps accuracy of the same order as FDM and the terms relax the Courant condition. The voxel FEM requires a similar amount of memory and only takes 1.2～1.4 times longer computation time. The voxel mesh can be generated considerably faster than the popular tetrahedral mesh. Both ground motions and static displacements due to a point or line source can be calculated using the voxel FEM approach. Comparisons with the reflectivity method and theoretical solutions demonstrate the successful implementation of the method, which is then applied to more complex problems.     

Relationship of topography to surface water chemistry with particular focus on nitrogen and organic carbon solutes within a forested watershed in Hokkaido,Japan

We studied the relationships between streamwater chemistry and the topography of subcatchments in the Dorokawa watershed in Hokkaido Island, northern Japan, to examine the use of topography as a predictor of streamwater chemistry in a watershed with relatively moderate terrain compared with other regions of Japan. Topographic characteristics of the Dorokawa watershed and its subcatchments were expressed as topographic index (TI) values, which ranged from 4·5 to 20·4 for individual grid cells (50 × 50 m²), but averaged from 6·4 to 7·4 for the 20 subcatchments. Streamwater samples for chemical analyses were collected four times between June and October 2002 from 20 locations in the watershed. The pH of water that passed through the watershed increased from ～5·0 to 7·0, with major increases in Na⁺ and Ca²⁺ and marked decreases in NO₃^? and SO . Distinctive spatial patterns were observed for dissolved organic carbon (DOC), dissolved organic nitrogen (DON), and NO₃^? concentrations of streamwater across the watershed. Statistical analyses indicated significant linear relationships between the average TI values of subcatchments and DOC, DON, and NO₃^? concentrations. Furthermore, the proportion of DOC in streamwaters in the wet season increased with TI values relative to other nitrogen species, whereas NO₃^? concentrations decreased with TI. The gradients of soil wetness and the presence of wetlands explained many of the observed spatial and temporal patterns of DOC, DON, and NO₃^? concentrations in the surface waters of the Dorokawa watershed. Our results suggest that the TI is especially useful for predicting the spatial distribution of DOC, DON and NO₃^? in the surface waters of Hokkaido, where topographical relief is moderate and wetlands more common than in other regions of Japan. Copyright © 2006 John Wiley & Sons, Ltd.     

Computer-processed potentiometric titration for the determination of calcium and magnesium in sea water

An improved potentiometric titration method for the determination of calcium and magnesium in sea water has been newly devised. In this method, a mini-computer is used for the automation of titrations, and ion-selective electrodes are used as an end-point detector. Calcium is determined by titration with EGTA, and total alkaline earth metals (magnesium + calcium + strontium) by titration with EDTA. Magnesium can be determined by the difference, strontium having been determined by a suitable method. In the present method, calcium and magnesium in sea water can be determined with a precision of 0.1 %.     

The characterization and distribution of aerosol and gaseous species in the winter monsoon over the western pacific ocean

In order to obtain a better understanding of the behavior of aerosols and SO₂ in the longrange transport through a marine boundary layer, a simple box-model is applied to the evaluation of the residence times of the species from the concentrations of gases and aerosols measured simulataneously on two islands in the West Pacific Ocean in the north-west monsoon. For Aitken and large particles, the residence time is varied from 3.7 to 7.4 days depending on the particle size, and their flux to the sea is equal to or slightly smaller than that of the free atmosphere. The residence time of giant particles is about one day and their flux to the sea is three times larger than that of the free atmosphere. The residence time of SO₂ is 15 hr, and the relative SO₂ mass flows of the deposition to the sea, of the diffusion to the free atmosphere and of the transformation to SO₄ ^2- are approximately 4, 1 and 1, respectively.     

Isotopic Evaluation on the Genesis of the Kanamaru Pegmatite Deposit, Niigata Prefecture, Japan

Abstract: The Kanamaru pegmatite is one of the two largest pegmatite deposits in Japan and is rich in K–feldspar, which is maximum microcline. The mica minerals and microcline from this ore deposit were studied for K-Ar ages and δ¹⁸O ratios. Biotite gives K-Ar age of 78. 5 Ma, but microcline is much younger in age of 66. 0 and 58. 3 Ma. Oxygen isotopic equilibrium temperature of the biotite-quartz pair is 625°C. Using closure temperature of the K-Ar system for biotite and microcline, cooling rate of the pegmatite is calculated to be 88°C/m. y. at the beginning but 12°C/m. y. in the later period of the cooling history. The microcline is slightly depleted in ¹⁸O, yielding apparent equilibrium temperatures of 230°C and 155°C on quartz-microcline pairs. It is suggested that the fluid phase, magmatic and partly meteoric in origin, contained in the solidifying pegmatite, was interacted with the original orthoclase and converted it to microcline during the subsolidus stage.     

Calcium-alkalinity relationship in the North Pacific

The dissolution of calcium carbonate in deep ocean water causes variation in calcium concentration (Ca) and alkalinity (TA) in the ratio of one to two. The decomposition of organic matter generates nitric acid, phosphoric acid and sulfuric acid. A proton flux which is derived from this process also changes alkalinity. Using the variation in nitrate concentration (NO₃) as an index of the proton flux, the relationship betweenCa,TA andNO₃ is expressed asCa=0.5TA+0.63NO₃ The values of Ca obtained from direct measurements in the North Pacific are in good agreement with the values estimated from this equation.     

Calcium-alkalinity-nitrate relationship in the north pacific and the Japan Sea

When biogenic particles of calcium carbonate and organic matter are formed or decomposed in sea water, calcium concentration, alkalinity and nitrate concentration are changed according to the following equation: ΔCa=0.5 ΔTA+0.63 ΔNO ₃, where ΔCa, ΔTA and ΔNO ₃ are the differences in calcium, alkalinity and nitrate concentration, respectively, between two waters. This relationship was pointed out in a previous paper of ours, and it holds well for data obtained from the North Pacific and the Japan Sea. Furthermore, in sea water containing little or no nutrients, the variations in calcium concentration and alkalinity are consistent with the variations expected from the formation and dissolution of calcium carbonate alone.     

Finite-element simulations of long-period ground motions: Japanese subduction-zone earthquakes and the 1906 San Francisco earthquake

Large earthquakes at shallow depths commonly excite long-period ground motions in distant sedimentary basins, thereby inflicting damage upon large-scale structures. For example, the 2003 Tokachi-oki earthquake, Japan, damaged oil tanks in the Yufutsu Basin, located 250 km from the epicenter. Similar long-range effects were also observed during the 2004 earthquake off Kii Peninsula, Japan. In this study, we present the results of simulations of these earthquakes undertaken using the finite element method (FEM) with a voxel mesh. In addition, to examine whether the 1906 San Francisco earthquake excited long-period ground motions in the Los Angeles-area basins, we performed long-period ground motion simulations of most of the California region. The FEM simulations confirmed the importance of path effects for the development of long-period ground motions.