Satellite detection of red tide in Ariake Sound, 1998–2001

High resolution SeaWiFS data was used to detect red tide events that occurred in the Ariake Sound, Japan, a small embayment known as one of the most productive areas in Japan. SeaWiFS chlorophyll data clearly showed that a large red tide event, which damaged seaweed (Nori) cultures, started early in December 2000 in Isahaya Bay, expanded to the whole sound and persisted to the end of February 2001. The monthly average of SeaWiFS data from May 1998 to December 2001 indicated that the chlorophyll peaks appeared twice a year, in early summer and in fall, after the peaks of rain and river discharge. The SeaWiFS data showed that the red tide event during 2000–2001 winter was part of the fall bloom; however, it started later and continued significantly longer than other years. Satellite ocean color data is useful to detect the red tide; however the algorithms require improvement to accurately estimate chlorophyll in highly turbid water and in red tide areas.     

Temporal change of clustered distribution of planktonic ciliates in Toyama Bay in summers of 1989 and 1990

Temporal change of clustered distribution in vertical profiles of three nutritional groups of planktonic ciliates, e.g. heterotrophic naked ciliates, mixotrophic naked ciliates and heterotrophic loricated ciliates, was investigated by following a drifting buoy in Toyama Bay on the Japan Sea coast of central Japan in summers of 1989 and 1990. Clustered distribution, represented as the mode of population density in the vertical plane, occurred mainly in the oligotrophic upper layer (0–50 m depth) above the subsurface chlorophyll-a maximum layer. Its clustered degree was stronger when the mode of population density in the vertical plane was formed at shallower depth, while its longevity was shorter as mentioned above. Vertical distribution of ciliates during summer in Toyama Bay is characterized by ephemeral clustered distribution, or in other wards, by rapid alternations of appearance and disappearance of the clustered distribution.     

Detecting fluctuations of the Kuroshio axis south of Japan using TOPEX/POSEIDON altimeter data

Sea surface dynamic topography (SSDT) can be divided into temporal mean SSDT and fluctuation SSDT. The former is approximated with a climatological mean SSDT and the latter is derived from satellite altimetry data, to give an approximated total SSDT (called a composite SSDT). The method is applied to detecting fluctuations of the Kuroshio axis south of Japan using TOPEX/POSEIDON altimeter data from the first year mission in 1992–1993. The fluctuation SSDT averaged over a wide area south of Japan clearly shows an annual cycle with an amplitude of about 15 cm. Temporal changes of SSDT along a subsatellite track crossing the Kuroshio compare moderately well with those estimated from repeated hydrographic observations, although there is a discrepancy of unknown origin. The composite SSDT also compares well with SSDT estimated from the same hydrographic data. Horizontal distribution of the surface geostrophic velocity component normal to subsatellite tracks is derived every ten days from the composite SSDT. Most locations of estimated strong eastward geostrophic velocities coincide well with locations of the Kuroshio axis determined every 15 days fromin situ surface velocity measurements on various vessels; for example, a fairly large meander of the Kuroshio south of Honshu is clearly detected. It is concluded that the composite SSDT can be used reliably to detect fluctuations of the Kuroshio axis south of Japan.     

Wave-induced porewater pressure in a seabed with inhomogeneous permeability

We have developed an analytical treatment for the wave-induced response of a seabed with inhomogeneous permeability, by employing the boundary layer approximation. By assuming an exponential function for vertically distributed permeability, the inhomogeneous consolidation equation was formulated as Bessel's differential equation. Subsequently, we proposed some representative quantities for examining the effects of inhomogeneity of permeability. The main effects of inhomogeneity are that the boundary layer becomes thinner and the phase delay in the direction of depth becomes larger as the inhomogeneity becomes significant. In addition, it was shown that the boundary layer is properly evaluated by using the averaged permeability.     

A laboratory experiment to monitor the contact state of a fault by transmission waves

We performed a series of laboratory experiments in which elastic waves were transmitted across a simulated fault. Two types of experiments were carried out: (1) Normal Stress Holding Test (NSHT): normal stress was kept constant for about 3 h without shear stress and transmission waves were observed. (2) Shear Stress Increasing Test (SSIT): shear stress was gradually increased until a stick-slip event occurred. Transmission waves were continuously observed throughout the process of stress accumulation. We focused on the change in transmission waves during the application of shear stress and especially during precursory slips.It was found in NSHT that the amplitude of transmission waves linearly increased with the logarithm of stationary contact time. The increase amounted to a few percent after about 3 h. Creep at asperity contacts is responsible for this phenomenon. From a theoretical consideration, it was concluded that the real contact area increased with the logarithm of stationary contact time.We observed in SSIT a significant increase in wave amplitude with shear stress application. This phenomenon cannot be attributed to the time effect observed in NSHT. Instead, it can be explained by the mechanism of “junction growth” proposed by Tabor. Junction growth yields an increase in real contact area. It is required for junction growth to occur that the material in contact is already plastic under a purely normal loading condition. A computer simulation confirmed that this requirement was satisfied in our experiments. We also found that the rate at which the amplitude increased was slightly reduced prior to a stick-slip event. The onset time of the reduction well coincides with the onset of precursory slip. The cause of the reduction is attributed to the reset of stationary contact time due to displacement. This interpretation is supported by the result of NSHT. Taking the time of stationary contact in SSIT into account, we may expect the change in wave amplitude to be, at most, only a few percent. The observed slight reduction in increasing rate is, in this sense, reasonable. The static stiffness of the fault also decreases with precursory slip. It was also found that low frequency waves are a better indicator of precursory slip than high frequency waves. This might suggest that low frequency waves with longer wavelength are a better indicator of average behavior of faults. The problem, however, merits a further investigation. The shifts in phase were also found to be a good indicator of the change in contact state of the fault. The changes in both amplitude and phase of transmission waves are unifyingly understood through the theory of transmission coefficient presented by Pyrak-Nolte et al. Rough surfaces have a tendency to give larger stick-slips than smooth surfaces. The amount of precursory slip is larger for rough surfaces than for smooth surfaces. Although it was confirmed by a computer simulation that rough surfaces have larger contact diameters than smooth surfaces, the rigorous relationship between the surface roughness (contact diameter) and the amount of precursory slips was not established.     

Carbon isotopes and petrography of kerogens in ∼3.5-Ga hydrothermal silica dikes in the North Pole area, Western Australia

More than 600 specimens of ∼3.5 Ga-old hydrothermal silica dikes from the North Pole area, Pilbara craton, Western Australia, have been studied petrographically. The kerogens in 44 samples have been analyzed isotopically (C and N) and chemically (C, N, and H). The silica dikes are composed mainly of fine-grained silica (modal abundance: >97%) and are classified into two types by minor mineral assemblages: B(black)-type and G(gray)-type. The B-type silica dikes contain kerogen (0.37 to 6.72 mgC/g; average 2.44 mgC/g, n = 21) and disseminated sulfides, dominantly pyrite and Fe-poor sphalerite. In some cases, carbonate and apatite are also present. Their silica-dominated and sulfide-poor mineral assemblages suggest precipitation from low-temperature reducing hydrothermal fluid (likely 100-200°C). On the other hand, the G-type silica dikes are sulfide-free and concentrations of kerogen are relatively low (0.05 to 0.41 mgC/g; average 0.17 mgC/g, n = 13). They typically contain Fe-oxide (mainly hematite) which commonly replaces cubic pyrite and rhombic carbonate. Some G-types occur along secondary quartz veins. These textures indicate that the G-type silica dikes were formed by postdepositional metasomatism (oxidation) of the B-types, and that the B-types probably possess premetasomatic signatures. The δ¹³C values of kerogen in the B-types are −38.1 to −33.1‰ (average −35.9‰, n = 21), which are ∼4‰ lower than those of the G-types (−34.5 to −30.0‰; average −32.2‰, n = 19), and ∼6‰ lower than bedded chert (−31.2 to −29.4‰; average −30.5‰, n = 4). This indicates the preferential loss of ¹²C during the metasomatism (estimated fractionation factor: 0.9985). Considering the metasomatic effect on carbon isotopes with probably minor diagenetic and metamorphic overprints, we conclude that the original δ¹³C values of the kerogen in the silica dikes would have been heterogeneous (∼5‰) and at least some material had initial δ¹³C values of ≤ −38‰. The inferred ¹³C-depletions of organic carbon could have been produced by anaerobic chemoautotrophs such as methanogen, but not by aerobic photoautotrophs. This is consistent with the estimated physical and chemical condition of the hydrothermal fluid, which was probably habitable for anaerobic and thermophilic/hyperthermophilic chemoautotrophs. Alternatively, the organic matter may have been possibly produced by abiological reaction such as Fischer-Tropsch Type (FTT) synthesis under the hydrothermal condition. However, the estimated condition is inconsistent with the presence of the effective catalysts for the FTT reaction (i.e., Fe-Ni alloy, magnetite, and hematite). These lines of evidence suggest the possible existence of biosphere in the ∼3.5 Ga sub-seafloor hydrothermal system.     

Effectiveness of Mg–Al-layered double hydroxide for heavy metal removal from mine wastewater and sludge volume reduction

Health hazards from heavy metal pollution in water systems are a global environmental problem. Of similar concern is sludge that results from wastewater treatment due to unsatisfactory sludge management technology. Therefore, the effectiveness of using Mg–Al-layered double hydroxide in the removal of heavy metals from mine wastewater was tested and compared with that of calcium hydroxide [Ca(OH)₂], which is a common treatment method for heavy metal removal. Initially, the mine wastewater contained cations of the heavy metals iron (Fe), zinc (Zn), copper (Cu), and lead (Pb). The Mg–Al-layered double hydroxides were able to remove 371, 7.2, 121, and 0.4 mg/L of these pollutants, respectively, using the co-precipitation method. The removal of these metals is most effective using 0.5 g Mg–Al-layered double hydroxide (Mg/Al molar ratio 4) and 20 min of shaking. Zn was removed by the formation of Zn(NO₃)(OH)·H₂O and Zn₅(NO₃)₂(OH)₈ when LDH, Mg/Al molar ratios of 4 and 2, respectively, were used. Similarly, Fe, Cu, and Pb were removed by the formation of Fe–Al-layered double hydroxide, Cu₂(OH)₃·NO₃ and Pb₄(OH)₄(NO₃)₄, respectively. While Ca(OH)₂ is also capable of reducing the heavy metal concentrations below the Japanese recommended values, this analysis shows that using 0.5 g Mg–Al-layered double hydroxide is a better treatment condition for mine wastewater, because it generates lower sludge volumes than 0.1 g of Ca(OH)₂. The measured sludge volume was 1.5 mL for Mg–Al-layered double hydroxide and 2.5 mL for Ca(OH)₂, a nearly twofold further reduction.     

High-reliability zircon separation for hunting the oldest material on Earth: An automatic zircon separator with image-processing/microtweezers-manipulating system and double-step dating

Despite the recent development in radiometric dating of numerous zircons by LA-ICPMS, mineral separation still remains a major obstacle, particularly in the search for the oldest material on Earth. To improve the efficiency in zircon separation by an order of magnitude, we have designed/developed a new machine-an automatic zircon separator(AZS). This is designed particularly for automatic pick-up of100 μm-sized zircon grains out of a heavy mineral fraction after conventional separation procedures. The AZS operates in three modes:(1) image processing to choose targeted individual zircon grains out of all heavy minerals spread on a tray,(2) automatic capturing of the individual zircon grains with microtweezers, and(3) placing them one-by-one in a coordinated alignment on a receiving tray. The automatic capturing was designed/created for continuous mineral selecting without human presence for many hours. This software also enables the registration of each separated zircon grain for dating, by recording digital photo-image, optical(color) indices, and coordinates on a receiving tray. We developed two new approaches for the dating; i.e.(1) direct dating of zircons selected by LA-ICPMS without conventional resin-mounting/polishing,(2) high speed U-Pb dating, combined with conventional sample preparation procedures using the new equipment with multiple-ion counting detectors(LA-MIC-ICPMS).With the first approach, Pb-Pb ages obtained from the surface of a mineral were crosschecked with the interior of the same grain after resin-mounting/polishing. With the second approach, the amount of time required for dating one zircon grain is ca. 20 s, and a sample throughput of 150 grains per hour can be achieved with sufficient precision(ca. 0.5%).We tested the practical efficiency of the AZS, by analyzing an Archean Jack Hills conglomerate in Western Australia with the known oldest(4.3 Ga) zircon on Earth. Preliminary results are positive; we were able to obtain more than 194 zircons that are over 4.0 Ga out of ca. 3800 checked grains, and 9 grains were over 4300 Ma with the oldest at 4371 ± 7 Ma. This separation system by AZS, combined with the new approaches, guarantees much higher yield in the hunt for old zircons.     

Preliminary results of a study on the responses of sedimentary rocks to shaft excavation

A shaft with a diameter of 6 m and a depth of 150 m was excavated in sedimentary rocks. Anin situ study was carried out to identify the size and rock properties of the zone disturbed by the excavation and to estimate the applicability of existing methods for measuring and modelling such a zone. The study detected displacements and changes in the properties of the rock within about 1 m of the shaft wall. Numerical analysis using a model of continuous rock mass can simulate some features of disturbance due to excavation, although slight discrepancies exist between the measured and analysed results. The excavation disturbance detected in this study is probably due to fracturing, the redistribution of stress and undersaturation. However, the relationships between the disturbance due to excavation and these processes have not been quantified. The existing methods for measurements and numerical analysis can provide important information on the disturbance due to excavation but need to be improved to understand disturbance due to excavation further.     

Slab weakening: Mechanical and thermal-mechanical consequences for slab detachment

Abstract Slab detachment is a geophysical instability whose manifestation can be revealed by seismic tomography. Evidence of this phenomenon is in the Dinarides/Hellenic and the New Hebrides subduction zones. Subducted slabs in these regions are torn horizontally at depths ranging from 100 to 300 km. We constructed a viscoelastic three-dimensional finite element model and investigated the state of stress. We found that an area with high stress concentration of the order of several hundred MPa is formed near the tip of the tear inside the slab, which can cause lateral migration of the tear. Favorable conditions for slab detachment are characterized by large interplate frictional force at a subduction zone and small slab resistance force deeper down. Stress concentration increases with the down-dip tension inside the slab. The phenomenon of slab weakening has also been studied from a thermal-mechanical standpoint, using a two-dimensional convection model with non-Newtonian, temperature-dependent rheology. The stress-dependent rheology plays an important role in causing local weakening of the descending slab. In strongly time-dependent situations the fast descending slab is not strong everywhere but has a weak region in the middle, making it vulnerable to slab detachment. The presence of viscous heating will enhance slab detachment tendency by further weakening the interior by the frictional heating. Besides these effects, there are other mechanisms which can also weaken the slab interior and help to make slabs more pliable and susceptible to detachment.