Basic theory of three-dimensional landslide modeling based on the logical model of geologic structure

Dynamic visualization of landslide cross-sections is important for understanding the structure and mechanism of landslide formation. Moreover, the modeling of geologic information plays an effective role in geo-hazard assessment and their mitigation. In this study, we developed the basic theory of a three-dimensional landslide modeling and applied it to the Nigawa landslide of the Hyogo Prefecture in central Japan. The construction of this model is based on the boundary surfaces of slump blocks and geologic units, and the hierarchical relationships between these surfaces. An application algorithm was validated and the model proved efficient in depicting the nature of landslides in the Nigawa area.     

Photographic observations of deep-sea infaunal ophiuroids in Suruga Bay,central Japan: an application of a free-fall system of time-lapse cameras and current meters

Underwater observations of infaunal amphiurid ophiuroids were made at a depth of about 480 m in Suruga Bay, central Japan, using a free-fall system which consists of time-lapse stereo-photography units and current meters. The megabenthos fauna was characterized by the dominance of infaunal echinoderms; in particular, amphiurid ophiuroids were numerically dominant. The density and biomass of the amphiurids were 170 m^–2 and 37 g m^–2, respectively. They buried their discs in the sediment and extended their arm tips out of the sediment surface. They protruded 2.2 arms per individual on the average. Strong bottom currents were observed, and the average velocity was 12 cm sec^–1 at 4 m above the sea floor. No arm tip was observed to be raised vertically into the water column for suspension feeding utilizing the bottom current, and amphiurids were considered to be primarily a surface deposit feeder at the present site.     

Horizontal displacements in Japan

Recent horizontal displacements in Japan are discussed by using the results of the first-order triangulation surveys. Large horizontal displacements are found in East Hokkaido, Tohoku, South Kanto, Tokai and Nankai districts (Fig.1). Crustal activity in these districts is also briefly discussed. The original triangulation survey in Japan neglected Laplace (azimuth) observations and consequently rotation and divergence around the assumed fixed point often appear as horizontal displacements. This is especially true in the central Japan where horizontal displacements result from cancelling the apparent rotation and divergence around the assumed fixed point (Fig.2).On the other hand, strain measurements avoid such shortcoming. Considering with the results of the first triangulation, the horizontal earth-strain can be calculated for every subsequent triangulation net. It is interesting that the velocity of the maximum shear strain is almost 2–3 · 10^/t7/year throughout Japan (Fig.3), even though the seismic activities show large regional discrepancies.     

Latitudinal changes in larger benthic foraminiferal assemblages in shallow-water reef sediments along the Ryukyu Islands, Japan

Abstract The taxonomic diversity of hermatypic corals decreases with increasing latitude, which correlates with sea‐surface temperatures. However, little is known about latitudinal changes in the taxonomic diversity and biogeographic patterns of larger benthic foraminifera, although their physiological requirements are similar to those of hermatypic corals because of their symbiotic relationships with microalgae. The present study examined how the abundance and taxonomic composition of larger foraminiferal assemblages in shallow‐water reef sediments change with latitude along the Ryukyu Islands (Ryukyus), which are located near the northern limit of coral‐reef distributions in the western Pacific Ocean. Three islands from different latitudes in the Ryukyus were selected to investigate latitudinal changes in larger foraminiferal assemblages: Ishigaki Island (24°20′N, 124°10′E), Kudaka Island (26°09′N, 127°54′E) and Tane‐ga‐shima Island (30°20′N, 131°E). Four sediment samples were taken at each of three topographic sites (beach, shallow lagoon and reef crest) on the reef flat of each island. Foraminiferal tests of a 2.0‐ to 0.5‐mm size fraction were selected, identified and counted. The variations in foraminiferal abundance in reef sediments from three latitudinally different islands exhibit two contrasting trends along reef flats: a shoreward decrease on Ishigaki and Tane‐ga‐shima Islands and a shoreward increase on Kudaka Island. A total of 25, 24 and 13 foraminiferal taxa were identified in Ishigaki, Kudaka and Tane‐ga‐shima Islands, respectively. Baculogypsina sphaerulata, Neorotalia calcar and Amphistegina spp. were dominant (i.e. >3% of foraminiferal assemblages) in the three islands. Calcarina gaudichaudii and Calcarina hispida were common on Ishigaki and Kudaka Islands but were absent on Tane‐ga‐shima Island. Larger foraminiferal assemblages from three different reef‐flat environments on Ishigaki Island can be distinguished, whereas those from the three environments on Kudaka and Tane‐ga‐shima Islands are similar in composition. These latitudinal changes in larger foraminiferal assemblages in reef sediments may possibly be caused by variations in the topography of reef flats, distributions and standing crops of living foraminifers on reef flats, and the northern limit of some calcarinid species in the northern Ryukyus.     

Seismicity and tectonics of the far Western Aleutian Islands

Study of focal mechanisms of earthquakes in the Near and Komandorsky Islands indicate that there are several distinct zones of tectonic activity. South of the Near Islands, normal faulting occurs in the trench east of 172°E and low-angle thrusting dominates the Aleutian ridge. Mechanisms indicate underthrusting as far west as Mednyy Island with strike-slip faulting restricted to the south and west of Beringa Island. A zone of northeast striking left-lateral faulting near 1645.°E is proposed to separate the Aleutian Ridge from Kamchatka Peninsula. This motion, as well as faulting north of the Komandorsky Islands, may be related to the existance of a buffer plate comprising the Aleutian Ridge in the Komandorsky Islands. Active subduction terminates near 173°E and the faulting north of the Komandorsky Islands may, in part, be due to the bouyancy of a remnant slab. Depth phase modelling indicates bulletin-reported depths are overestimated due to a misidentification of depth phases.     

Development of regional simulation of seismic ground‐motion and induced liquefaction enhanced by GPU computing

Much research has been conducted for physics‐based ground‐motion simulation to reproduce seismic response of soil and structures precisely and to mitigate damages caused by earthquakes. We aimed at enabling physics‐based ground‐motion simulations of complex three‐dimensional (3D) models with multiple materials, such as a digital twin (high‐fidelity 3D model of the physical world that is constructed in cyberspace). To perform one case of such simulation requires high computational cost and it is necessary to perform a number of simulations for the estimation of parameters or consideration of the uncertainty of underground soil structure data. To overcome this problem, we proposed a fast simulation method using graphics processing unit computing that enables a simulation with small computational resources. We developed a finite‐element‐based method for large‐scale 3D seismic response analysis with small programming effort and high maintainability by using OpenACC, a directive‐based parallel programming model. A lower precision variable format was introduced to achieve further speeding up of the simulation. For an example usage of the developed method, we applied the developed method to soil liquefaction analysis and conducted two sets of simulations that compared the effect of countermeasures against soil liquefaction: grid‐form ground improvement to strengthen the earthquake resistance of existing houses and replacement of liquefiable backfill soil of river wharves for seismic reinforcement of the wharf structure. The developed method accelerates the simulation and enables us to quantitatively estimate the effect of countermeasures using the high‐fidelity 3D soil‐structure models on a small cluster of computers.     

Geophysical monitoring and reactive transport modeling of ureolytically-driven calcium carbonate precipitation

Ureolytically-driven calcium carbonate precipitation is the basis for a promising in-situ remediation method for sequestration of divalent radionuclide and trace metal ions. It has also been proposed for use in geotechnical engineering for soil strengthening applications. Monitoring the occurrence, spatial distribution, and temporal evolution of calcium carbonate precipitation in the subsurface is critical for evaluating the performance of this technology and for developing the predictive models needed for engineering application. In this study, we conducted laboratory column experiments using natural sediment and groundwater to evaluate the utility of geophysical (complex resistivity and seismic) sensing methods, dynamic synchrotron x-ray computed tomography (micro-CT), and reactive transport modeling for tracking ureolytically-driven calcium carbonate precipitation processes under site relevant conditions. Reactive transport modeling with TOUGHREACT successfully simulated the changes of the major chemical components during urea hydrolysis. Even at the relatively low level of urea hydrolysis observed in the experiments, the simulations predicted an enhanced calcium carbonate precipitation rate that was 3-4 times greater than the baseline level. Reactive transport modeling results, geophysical monitoring data and micro-CT imaging correlated well with reaction processes validated by geochemical data. In particular, increases in ionic strength of the pore fluid during urea hydrolysis predicted by geochemical modeling were successfully captured by electrical conductivity measurements and confirmed by geochemical data. The low level of urea hydrolysis and calcium carbonate precipitation suggested by the model and geochemical data was corroborated by minor changes in seismic P-wave velocity measurements and micro-CT imaging; the latter provided direct evidence of sparsely distributed calcium carbonate precipitation. Ion exchange processes promoted through NH4+ production during urea hydrolysis were incorporated in the model and captured critical changes in the major metal species. The electrical phase increases were potentially due to ion exchange processes that modified charge structure at mineral/water interfaces. Our study revealed the potential of geophysical monitoring for geochemical changes during urea hydrolysis and the advantages of combining multiple approaches to understand complex biogeochemical processes in the subsurface.     

Double seismic zones and stresses of intermediate depth earthquakes

Summary. Data from Japanese local seismograph networks suggest that the stresses in double seismic zones are in-plate compression for the upper zone and in-plate tension for the lower zone; the stresses do not necessarily appear to be down-dip. It may therefore be possible to identify other double seismic zones on the basis of data which indicate that events with differing orientations of in-plate stresses occur in a given segment of slab.
A global survey of published focal mechanisms for intermediate depth earthquakes suggests that the stress in the slab is controlled, at least in part, by the age of the slab and the rate of convergence. Old and slow slabs are under in-plate tensile stresses and the amount of in-plate compression in the slab increases with increasing convergence rate or decreasing slab age. Young and fast slabs are an exception to this trend; all such slabs are down-dip tensile. Since these slabs all subduct under continents, they may be bent by continental loading. Double seismic zones are not a feature common to all subduction zones and are only observed in slabs which are not dominated by tensile or compressive stresses.
Unbending of the lithosphere and upper mantle phase changes are unlikely to be the causes of the major features of double zones, although they may contribute to producing some of their characteristics. Sagging or thermal effects, possibly aided by asthenospheric relative motion, may produce the local deviatoric stresses that cause double zones.     

Mathematical modeling of biologically mediated redox processes of iron and arsenic release in groundwater

A one dimensional reactive transport model was developed in order to illustrate the biogeochemical behavior of arsenic and iron reduction and release to groundwater that accounts for the reaction coupling the major redox elements under reducing environment. Mass transport equation and the method of characteristics were used considering fundamental geochemical processes to simulate transport processes of different pollutants in mobile phase. The kinetic sub-model describes the heterotrophic metabolisms of several microorganisms. To model a complete redox sequence (aerobic or denitrifiers, Fe(III)-reduction, respiration bacteria of iron and arsenic compounds, and As(V) reduction) four functional bacterial groups (X ₁, X ₂, X ₃, and X ₄) were defined. Microbial growth was assumed to follow Monod type kinetics. The exchange between the different phases (mobile, bio, and matrix) was also considered in this approach. Results from a soil column experiment were used to verify the simulation results of the model. The model depicts the utilization of oxygen, nitrate, iron oxide and arsenic as electron acceptors for oxidation of organic carbon (OC) in a column. The OC as electron donor is one of the most important factors that affect the iron and arsenic reduction bacterial activity.