Application of a Numerical Inverse Laplace Integration Method to Surface Loading on a Viscoelastic Compressible Earth Model

Normal mode approaches for calculating viscoelastic responses of self-gravitating and compressible spherical earth models have an intrinsic problem of determining the roots of the secular equation and the associated residues in the Laplace domain. To bypass this problem, a method based on numerical inverse Laplace integration was developed by Tanaka et al. (2006, 2007) for computations of viscoelastic deformation caused by an internal dislocation. The advantage of this approach is that the root-finding problem is avoided without imposing additional constraints on the governing equations and earth models. In this study, we apply the same algorithm to computations of viscoelastic responses to a surface load and show that the results obtained by this approach agree well with those obtained by a time-domain approach that does not need determinations of the normal modes in the Laplace domain. Using the elastic earth model PREM and a convex viscosity profile, we calculate viscoelastic load Love numbers (h, l, k) for compressible and incompressible models. Comparisons between the results show that effects due to compressibility are consistent with results obtained by previous studies and that the rate differences between the two models total 10–40%. This will serve as an independent method to confirm results obtained by time-domain approaches and will usefully increase the reliability when modeling postglacial rebound.     

An Application of Hydraulic Tomography to a Large‐Scale Fractured Granite Site,Mizunami, Japan

While hydraulic tomography (HT) is a mature aquifer characterization technology, its applications to characterize hydrogeology of kilometer‐scale fault and fracture zones are rare. This paper sequentially analyzes datasets from two new pumping tests as well as those from two previous pumping tests analyzed by Illman et al. (2009) at a fractured granite site in Mizunami, Japan. Results of this analysis show that datasets from two previous pumping tests at one side of a fault zone as used in the previous study led to inaccurate mapping of fracture and fault zones. Inclusion of the datasets from the two new pumping tests (one of which was conducted on the other side of the fault) yields locations of the fault zone consistent with those based on geological mapping. The new datasets also produce a detailed image of the irregular fault zone, which is not available from geological investigation alone and the previous study. As a result, we conclude that if prior knowledge about geological structures at a field site is considered during the design of HT surveys, valuable non‐redundant datasets about the fracture and fault zones can be collected. Only with these non‐redundant data sets, can HT then be a viable and robust tool for delineating fracture and fault distributions over kilometer scales, even when only a limited number of boreholes are available. In essence, this paper proves that HT is a new tool for geologists, geophysicists, and engineers for mapping large‐scale fracture and fault zone distributions.     

Frequency–Magnitude Distribution of −3.7 ≤ M W ≤ 1 Mining-Induced Earthquakes Around a Mining Front and b Value Invariance with Post-Blast Time

We investigated frequency-magnitude distribution (FMD) of acoustic emissions (AE) occurring near an active mining front in a South African gold mine, using a catalog developed from an AE network, which is capable of detecting AEs down to M _W  ?5. When records of blasts were removed, FMDs of AEs obeyed a Gutenberg?Richter law with similar b values, not depending on post-blasting time from the initial 1-min interval through more than 30 h. This result denies a suggestion in a previous study (Richardson and Jordan Bull Seismol Soc Am, 92:1766–1782, 2002) that new fractures generated by blasting disturb the size distribution of background events, which they interpreted as slip events on existing weak planes. Our AE catalog showed that the GR law with b ～ 1.2 was valid between M _W  ?3.7 and 0 for AEs around the mining front. Further, using the mine’s seismic catalog, which covers a longer time period of the same area, we could extend the validity range of the GR law with the same b value up to M _W 1.     

Influence of the Oyashio Current and Tsugaru Warm Current on the circulation and water properties of Otsuchi Bay,Japan

Mooring and hydrographic observations were conducted from September 2012 to May 2014 at the mouth of Otsuchi Bay, a ria along the Pacific coast of Japan. Our observations quantitatively demonstrated that the circulation and the water properties of Otsuchi Bay are strongly influenced by the Tsugaru Warm Current (TWC) and Oyashio Current (OY) at seasonal and subseasonal time scales. Two bottom-mounted velocity profilers and temperature and salinity measurements beneath the near-surface halocline showed a counterclockwise lateral circulation pattern related to the TWC, which was enhanced from summer to autumn. From winter to early spring, the lateral circulation patterns related to the TWC weakened and the influence of the OY occasionally increased. When the OY was weak, surface flows became an overturning structure, with outflows in the upper layer and inflows in the lower layer. When the OY was strong and passed close to the Sanriku coast, the circulation became highly variable and intermittent. Intrusions of the markedly low-salinity OY water were observed on two occasions and persisted for periods of several weeks to several months. Salinity was sometimes less than 33.7, the lower limit of the typical TWC from late summer to autumn even when the TWC dominates. We suggested that this is the seasonal fluctuations of the TWC itself, as the upstream current of the Tsushima Warm Current is freshened in summer as a result of the influence of the Changjiang River. The surface water was generally fresher in the south of the bay than in the north, suggesting the Coriolis deflection of the river plume.     

Geotechnical sealing material for coastal disposal facility for soils and wastes contaminated by radioactive cesium

Safe disposal of wastes produced due to the process of decontamination in and around Fukushima No. 1 nuclear power plant is an urgent requirement. The purpose of this study was to develop a sealing material which can be used as an engineered barrier for a final disposal facility for the soils and wastes contaminated by radioactive cesium. The analyses conducted based on 4-m-thick sealing layer revealed that the hydraulic conductivity of the sealing material needs to be kept below 5.0?×?10^?10?m/s to avoid the seepage of contaminants below the environmentally safe limits. Sealing material was developed using marine clay–bentonite mixture and the engineering characteristics were examined. The results of laboratory experiments showed that, with the addition of bentonite, the hydraulic conductivity equal to or less than 5.0?×?10^?10?m/s was achieved when the effective consolidation stress is equal to or more than 27?kPa. From the tests for adsorption properties for cesium, it was found that the sealing material showed the significant capacity of adsorption for cesium in seawater. It was concluded that the construction of waste disposal facility on the sandy seafloor is feasible using the sealing layer proposed in this study.     

Assessing the role of forests in mitigating eutrophication downstream of pasture during spring snowmelt

Little research has examined whether forests reduce stream water eutrophication in agricultural areas during spring snowmelt periods. This study evaluated the role of forests in ameliorating deteriorated stream water quality in agricultural areas, including pasture, during snowmelt periods. Temporal variation in stream water quality at a mixed land‐use basin (565 ha: pasture 13%, forestry 87%), northern Japan, was monitored for 7 years. Synoptic stream water sampling was also conducted at 16 sites across a wide range of forest and agricultural areas in a basin (18.3 km²) in spring, summer and fall. Atmospheric nitrogen (N) and phosphorus (P) deposition were measured for 4 years. The results showed that concentration pulses of nitrate, organic N and total P in stream water were observed when discharge increased during spring snowmelt. Their concentrations were high when silicate concentrations were low, suggesting surface water exported from pasture largely contributed to stream water pollution during snowmelt. Atmospheric N and P deposition (4.1 kg N ha^?1 y^?1; 0.09 kg P ha^?1 y^?1, respectively) was too low to affect the background concentrations of N and P in streams from forested areas. Reduction of eutrophication caused by nutrients from pasture was mainly due to dilution by water containing low concentrations of N and P exported from forested areas, whereas in‐stream reduction was not a dominant process. Results indicate that forests have a limited capacity to reduce the concentration pulses of N and P in stream water during snowmelt in this study basin. Copyright © 2014 John Wiley & Sons, Ltd.     

Identification of 1771 Meiwa Tsunami deposits using a combination of radiocarbon dating and oxygen isotope microprofiling of emerged massive Porites boulders

The Meiwa Tsunami, one of the largest tsunamis recorded in historical documents in Japan struck Ishigaki Island and neighboring islands of the Ryukyu Arc in April 1771 AD, killing more than 12 000 people. An enormous number of massive Porites coral boulders are scattered on the shore and in the reef moat of eastern Ishigaki Island. Although these boulders likely were cast ashore by the Meiwa Tsunami, a detailed examination has not yet been conducted. When the marine reservoir effect is taken into account, one of mode values of calibrated radiocarbon dates possibly corresponds to the time of the 1771 event. However, the range of calibrated radiocarbon ages indicates that the transport of the boulders cannot be ascribed solely to the tsunami. Oxygen isotope microprofiling, which indicates sea-surface temperature variation, was therefore conducted to further investigate the mechanism of transport. The results suggest that the skeletal growth of most coral colonies was interrupted in summer or autumn; hence, tropical storms and typhoons are also very likely to be agents of transport. Thus, by combining radiocarbon dating with oxygen isotope microprofiling to investigate Porites coral boulders, it is possible to separate paleotsunami boulders from those transported by storm events as far as tsunamis occurred during the non-storm season.     

Hawaiian double volcanic chain triggered by an episodic involvement of recycled material: Constraints from temporal Sr–Nd–Hf–Pb isotopic trend of the Loa-type volcanoes

The two parallel loci of recent Hawaiian volcanoes, Kea and Loa, have been regarded as the best targets to interpret the chemical structure of an upwelling mantle plume derived from the lower mantle. Here we show that the Sr–Nd–Hf–Pb isotopic data of the shield-building lavas along the Loa locus form a systematic trend from the main shield stage of Koolau (> 2.9 Ma) to the active Loihi volcanoes. During the growth of the Koolau volcano, the dominant material in the melting region successively changed from the proposed KEA, DMK (depleted Makapuu), to EMK (enriched Makapuu) components. The proportion of EMK, dominated by a recycled mafic component, is typified by some Koolau Makapuu-stage and some Lanai lavas. Subsequently, the EMK component decreased and LOIHI component increased toward the Loihi lavas. The temporal coincidence between the episodically elevated magma production rate and the abrupt appearance of the typical Loa-type lavas that is restricted to the last 3 Myr should be linked to magma genesis. We suggest that the abrupt appearance of Loa-type magmatism should be attributed to the transient incorporation of the relatively dense recycled material and surrounding less degassed lower mantle material that accumulated near the core–mantle boundary into the upwelling plume. This episodic involvement could have been trigged by episodic thermal pulses and buoyancy increases in the plume. The continuous appearance of Kea-type lavas during the long history of Hawaiian-chain magmatism and the larger magma volume of Kea-type lavas relative to that of the Loa-type lavas in the last 3 Myr indicate that the Kea locus is closer to the thermal centre of the Hawaiian plume relative to that of the Loa locus.     

Eocene extensional tectonics in the Amakusa region,northern Ryukyu arc

Paleogene surface tectonics in Japan is not well understood because of the paucity of onshore Paleogene stratigraphic records except for those from accretionary complexes. Paralic Paleogene formations remaining in SW Japan are usually so thin that it is difficult to decipher the tectonics from them. However, the Eocene paralic sedimentary package with a thickness of kilometers indicates syn-depositional tectonic subsidence by a few kilometers in the Amakusa archipelago, west of Kyushu Island. Thus, we made a detailed geological map of the Eocene formations in an area of ~50 square kilometers in the northwestern part of the archipelago. We identified NE-SW and NW-SE trending normal faults, most of which were recognized by previous researchers, and also discovered low-angle faults. NW-SE trending ones are known to be of the Miocene. NE-SW trending and low-angle normal faults are the oldest map-scale structures in the Eocene ones. It is not obvious within the above-mentioned area whether those normal faults are accompanied by growth strata. However, the significant southeastward thickening of the Eocene formations across the Amakusa archipelago suggests that they filled a large half graben with the basin margin fault along the eastern side of the archipelago. This basin model is consistent with the N-S to NW-SE transport directions of the low-angle and NE-SW trending normal faults. Since many NE-SW to EW trending Eocene grabens were formed in the offshore regions west of Kyushu Island and in the East China Sea, the Amakusa region was probably a northeastern branch of the rift system. The geologic structures and depositional ages of the Eocene formations indicate that the Eocene extensional tectonics removed the overlying strata to some extent for the high-P/T Takahama Metamorphic Rocks which crops out to the south of our study area.