Numerical analysis on seepage failures of dike due to water level-up and rainfall using a water–soil-coupled smoothed particle hydrodynamics model

For seepage failures of dike due to water level-up and rainfall, surface infiltration and strength change induced by suction reduction are important factors; thus, numerical analysis should consider the coupling of water and soil, as well as the effect of saturation to obtain more precise failure mechanism. Based on the advanced smoothed particle hydrodynamics (SPH) method, this work proposed a two-phase-coupled SPH model in coordination with a novel constitutive model for unsaturated soils. Then, a triaxial compression test is simulated to check the applicability of the SPH method on the soil phase. After that, the failure test of a dike due to water level-up is discretized and simulated, from which the seepage process, the distribution of maximum shear strain, the slip surface, and pore water pressure are obtained. The two-phase-coupled SPH model is also applied to a slope failure test of heavy rainfall, and the results are compared to the model test. Finally, a dike failure test due to rainfall is analyzed using the proposed SPH model to reproduce the surface infiltration and suction reduction. The proposed SPH model provides several insights of seepage failures and can be a helpful tool for the analysis of dike failures induced by water level-up and rainfall.     

Use of Statistical Models to Analyze Periodic Seismicity Observed for Clusters in the Kanto Region, Central Japan

—A periodic pattern of seismicity has been reported for the Kinugawa cluster in the Kanto region, where several earthquake clusters are observed at depths between 40 and 90 km. To analyze this periodicity, statistical studies are performed for the Kinugawa cluster together with eight other clusters. Hypocentral parameters of the earthquakes with magnitudes 4.5 and larger for the period between 1950 and 1995 are taken from the JMA catalogue. The simple sinusoidal function, the exponential of sinusoidal function and the stress release model are applied as the intensity function. Model parameters are determined by the maximum likelihood method and the best model for each cluster is selected by using the Akaike Information Criterion (AIC). In six cases the sinusoidal model or the exponential of the sinusoidal model is selected as the best option and achieves AIC reductions of values between 2.4 and 13.2 units from the simple Poisson model. The stress release model is selected for two clusters. The three clusters, the Kinugawa, Kasumigaura, and Choshi clusters, have a similar optimal period of about 10 years, and align in the northwest–southeast direction at a similar depth range of 40 to 70 km. A model modified from the stress release model is applied to the three clusters so to analyze the relationship among them. In the modified model, an earthquake occurrence in one zone increases the stress in the other zone, which is different from the original stress release model which assumes a linear increase with time. Applying the modified model to the Kinugawa cluster, an AIC reduction from the Poisson model is significantly larger than the value obtained with the sinusoidal model. This suggests that the periodic seismicity observed for the Kinugawa cluster can be explained with the more comprehensive model than the sinusoidal model.     

Time Distribution of Immediate Foreshocks Obtained by a Stacking Method

—We apply a stacking method to investigate the time distribution of foreshock activity immediately before a mainshock. The foreshocks are searched for events with M≥ 3.0 within a distance of 50 km and two days from each mainshock with M≥ 5.0, in the JMA catalog from 1977 through 1997/9/30. About 33% of M≥ 5.0 earthquakes are preceded by foreshocks, and 50–70% in some areas. The relative location and time of three types of representative foreshocks, that is, the largest one, the nearest one to the mainshock in distance, and the nearest one in time, are stacked in reference to each mainshock. The statistical test for stacked time distribution of foreshocks within 30km from and two days before mainshocks shows that the inverse power-law type of a probability density time function is a significantly better fit than the exponential one for all three types of representative foreshocks. Two explanations possibly interpret the results. One is that foreshocks occur as a result of a stress change in the region, and the other one is that a foreshock is the cause of a stress change in the region and it triggers a mainshock. The second explanation is compatible with the relationship between a mainshock and aftershocks, when an aftershock happens to become larger than the mainshock. However the values of exponent of the power law obtained for stacked foreshocks are significantly smaller than those for similarly stacked aftershocks. Therefore the foreshock–mainshock relation should not be explained as a normal aftershock activity. Probably an increase of stress during foreshock activity results in apparently smaller values of the exponent, if the second explanation is the case.     

Horizontal distribution and population dynamics of the dominant mysid Hyperacanthomysis longirostris along a temperate macrotidal estuary (Chikugo River estuary, Japan)

The estuarine turbidity maximum (ETM) that develops in the lower salinity areas of macrotidal estuaries has been considered as an important nursery for many fish species. Mysids are one of the dominant organisms in the ETM, serving as a key food source for juvenile fish. To investigate the horizontal distribution and population dynamics of dominant mysids in relation to the fluctuation of physical conditions (temperature, salinity, turbidity, and freshwater discharge), we conducted monthly sampling (hauls of a ring net in the surface water) along the macrotidal Chikugo River estuary in Japan from May 2005 to December 2006. Hyperacanthomysis longirostris was the dominant mysid in the estuary, usually showing peaks of density and biomass in or close to the ETM (salinity 1–10). In addition, intra-specific differences (life-cycle stage, sex, and size) in horizontal distribution were found along the estuary. Larger males and females, particularly gravid females, were distributed upstream from the center of distribution where juveniles were overwhelmingly dominant. Juveniles increased in size toward the sea in marked contrast with males and females. The findings suggest a possible system of population maintenance within the estuary; gravid females release juveniles in the upper estuary, juveniles grow during downstream transport, young males and females mature during the upstream migration. Density and biomass were primarily controlled by seasonal changes of temperature, being high at intermediate temperatures (ca. 15–25 °C in late spring and fall) and being low at the extreme temperatures (ca. 10 °C in midwinter and 30 °C in midsummer). High density (up to 666 ind. m⁻³) and biomass (up to 168 mg dry weight m⁻³) of H. longirostris were considered to be comparable with those of copepods in the estuary.     

Petrological and Geochemical Characteristics of Intrusive Rocks Related to Porphyry Copper Mineralization and the Implications for the Genesis of Deposits in the Namosi area, Viti Levu, Republic of the Fiji Islands

In the Namosi district, of Viti Levu Island, Republic of the Fiji Islands, we conducted chemical analysis on intrusive rocks and attempted to clarify the relationship between petrological characteristics and mineralization. Samples were collected from the Waisoi copper deposit, Waivaka corridor, Waisomo, Wailutelevu and Wainabama Southwest localities. The intrusive rocks in the district show the overall characteristics of calc-alkaline magnetite series. Using the normative composition diagram for granite classification, most of the intrusive rocks plot in the tonalite, granodiorite, quartz monzodiorite and quartz diorite fields. Rare earth element (REE) analysis revealed that at Waisoi the intrusive rocks tend towards an Eu positive anomaly, whereas most of those at Waisomo tend to have a negative anomaly. The (Ce_cn/Yb_cn)/Yb_cn ratio shows that the intrusive rocks in the Namosi district may be derived from a mantle origin magma possibly contaminated by crustal materials. Waisoi rocks have less contamination with continental crust and show higher gold content compared to those of the Waivaka corridor. The observed differences between Waisoi and Waivaka corridor in the Cu/Au ratio might also be caused by varying amounts of contamination by crustal materials. There is no clear correlation between the La/Yb ratio and K–Ar ages, however, the La/Yb ratio of the Waivaka corridor rocks has a wide range, suggesting a variable amount of crustal contamination of the magma from 7 to 5 Ma. The less crustal contamination and generation of relatively high silica contents under relatively oxidized conditions in the Waisoi magmas is correlated with Au-rich porphyry copper mineralization and in contrast more crustal contamination produced Au-poor porphyry copper mineralization at the Waivaka Corridor in the Namosi district.     

Palaeotemperature curve for the Late Cretaceous of the northwestern circum-Pacific

In the northwestern circum-Pacific, two main trends in Late Cretaceous temperatures can be recognized. (1) In general, a recurrent warming trend is thought to have begun in the Turonian–Campanian, reaching temperature maxima in the early Late Santonian and early Late Campanian, and temperature minima in the earliest Santonian and perhaps early Campanian. (2) During the Maastrichtian, temperatures dropped sharply, with only a slight warming in the early Late Maastrichtian. The existence of a thermal maximum at the Coniacian–Santonian transition has previously been expected, but is not confirmed by new isotopic results.     

Experiment and Theoretical Study of Earthquake Detection Capability by Means of Microwave Passive Sensors on a Satellite

This letter presents the possibility of detecting earthquakes (EQs) from microwaves emitted when rock fractures. The method is based on an experiment in which microwave emission was detected from rock fracturing in a laboratory for the first time in the world. First, the method of calibrating emitted microwave power from experimental data is presented. A model of microwave emission and propagation to a satellite is then proposed. An advantage of microwaves is that they penetrate the Earth's ionosphere, unlike radiowaves of frequencies lower than several tens of megahertz. The power received by a satelliteborne receiver is estimated by assuming parameters of a radiometer currently operating in orbit. The result indicates that a satelliteborne receiver can detect microwave signals generated by an EQ. Based on this result, we attempted to detect some features associated with an actual EQ from the data of the Advanced Microwave Scanning Radiometer for Earth Observation System aboard the remote sensing satellite Aqua.     

FDM Simulation of Seismic Waves, Ocean Acoustic Waves, and Tsunamis Based on Tsunami-Coupled Equations of Motion

We have developed a new, unified modeling technique for the total simulation of seismic waves, ocean acoustic waves, and tsunamis resulting from earthquakes, based on a finite difference method simulation of the 3D equations of motion. Using the equilibrium between the pressure gradient and gravity in these equations, tsunami propagation is naturally incorporated in the simulation based on the equations of motion. The performance of the parallel computation for the newly developed tsunami-coupled equations using a domain partitioning procedure shows a high efficiency coefficient with a large number of CPU cores. The simulation results show how the near-field term associated with seismic waves produced by shallow earthquakes leads to a permanent coseismic deformation of the ground surface, which gives rise to the initial tsunami on the sea surface. Propagation of the tsunami along the sea surface as a gravity wave, and ocean acoustic waves in seawater with high-frequency multiple P-wave reflections between the free surface and sea bottom, are also clearly demonstrated by the present simulations. We find a good agreement in the tsunami waveform between our results and those obtained by other simulations based on an analytical model and the Navier–Stokes equations, demonstrating the effectiveness of the tsunami-coupling simulation model. Based on this simulation, we show that the ratio of the amplitude of ocean acoustic waves to the height of the tsunami, both of which are produced by the earthquake, strongly depends on the rise time of the earthquake rupture. This ratio can be used to obtain a more detailed understanding of the source rupture processes of subduction zone earthquakes, and for implementing an improved tsunami alert system for slow tsunami earthquakes.