Behavior of Cu, Ni and Cd during nutrient depletion in a spring bloom in Funka Bay

The concentrations of Cu, Ni and Cd were determined in Funka Bay during a spring phytoplankton bloom, consisting of diatoms. Just after the bloom, both dissolved Cd and nutrients were removed in the euphotic zone. However, the removal ratio of Cd to phosphate was very different from that in seawater. The removal of Cd took place at a Cd/phosphate ratio of 0.07×10⁻³, which was lower than in seawater before the bloom (0.25×10⁻³), leading to an increase in this ratio in seawater exceeding 0.7×10⁻³ at the end of the bloom. Elevated concentrations of Cd and phosphate were observed in the deeper layer after the bloom due to the decomposition of detrital materials produced in the bloom. The ratio of Cd/phosphate in the regeneration step was 0.24×10⁻³ which was different from the removal ratio of 0.07×10⁻³. These observations suggest that the high Cd/phosphate ratio in the regeneration would reflect a relatively high regeneration rate of Cd than that of phosphate. No significant decrease in Cu and Ni concentrations was observed during the development of the bloom, suggesting that biological removal of these metals was not so significant during the spring bloom. The concentrations of Cd, Cu and silicate in surface waters increased after the bloom with decreasing salinity due to the influence of a spring thaw.     

Determination of plastic regions around underground openings by a coupled boundary-element-characteristics method

A two-dimensional hybrid method for solving elastoplastic problems in engineering is presented by coupling two existing methods, namely, the boundary element method and the characteristics method. The formulation of this method is presented, as well as an excellent procedure for the determination of the boundary between elastic and plastic regions. It is shown not only that this method is a powerful and accurate method for evaluating the shape and extent of the plastic region around rock caverns, which is of prime importance for the construction of rock caverns, but also applicable to a given range of the initial stress field ratio where only compressive failure occurs. Then, some typical examples are solved in order to check the accuracy of the solution by this method. Furthermore, its successful applications are presented and discussed to determine the shape and the extent of the plastic regions around parallel, circular and rectangular openings.     

The phase transitions of superdense matter and supernova explosions

Effects of the phase transitions of superdense matter on supernova explosions are investigated with the aid of an idealized equation of state on the assumptions of adiabatic collapse. It is found that in the case of strong phase transitions explosions become weaker, while in the case of weak phase transitions explosions become stronger. However, the increment of the ejected energy is not so large as suggested by Migdalet al. (1979).Paper presented at the IAU Third Asian-Pacific Regional Meeting, held in Kyoto, Japan, between 30 September–6 October, 1984.     

Downstream transition of the Tsushima Current west of Kyushu in summer

In order to clarify detailed current structures west of Kyushu, ADCP measurements were carried out in July and September 1990 by the quadrireciprocal method (Katoh, 1988) for removing diurnal and semidiurnal tidal flows from observed flows. On the basis of these results, together with data of routine oceanographic observations, we study the downstream transition of the Tsushima Current west of Kyushu in summer. In the southwest of the Goto Islands, a northward current identified as the Tsushima Current was clearly found. In the south of Cheju, a westward current bifurcated from the Tsushima Current. In the northwest of the Goto Islands, the Tsushima Current narrowed and its velocity became strengthened. Salinity of the Tsushima Current water was much diluted by a current from the Cheju Strait. Near the west coast of the Goto Islands, a countercurrent bifurcating from the Tsushima Current often occurred. The volume transport of the Tsushima Current was 2.3 Sv (1 Sv=10⁶ m³s^–1) on the northern side of latitude 31°N. The substantial bifurcation of the Tsushima Current toward the Eastern and Western Channels of the Tsushima Strait occurred in the vicinity of Tsushima. The volume transport through the Western Channel was two to three times larger than that through the Eastern Channel. The baroclinic component in volume transport of the Tsushima Current west of Kyushu was much smaller than that in the Japan Sea.     

A Generation Mechanism for Mesoscale Eddies in the Kuril Basin of the Okhotsk Sea: Baroclinic Instability Caused by Enhanced Tidal Mixing

Mesoscale eddies, particularly anticyclonic ones, are dominant features in the Kuril Basin of the Okhotsk Sea. In 1999, both surface drifter and hydrographic observations caught the same anticyclonic eddy northwest of Bussol’ Strait, which has a diameter of ∼100 km, typical surface velocity of 0.2–0.3 m s⁻¹, and less dense core extending to a depth of ∼1200 m. Based on an idea that the generation of mesoscale eddies is caused by strong tidal mixing in and around Kuril Straits, we have conducted a series of three-dimensional numerical model experiments, in which strong tidal mixing is simply parameterized by increasing coefficients of vertical eddy viscosity and diffusivity along the eastern boundary. Initially, a regular series of disturbances with a wavelength of ∼70 km starts to develop. The disturbances can be clearly explained by a linear instability theory and regarded as the baroclinic instability associated with the near-surface front formed in the region between the enhanced mixing and offshore regions. In the mature phase, the disturbances grow large enough that some eddies pinch off and advect offshore (westward), with the scale of disturbances increasing gradually. Typical eddy scale and its westward propagation speed are ∼100 km and ∼0.6 km day⁻¹, respectively, which are consistent with the observations by satellites. The westward propagation can be explained partly due to nonlinear effect of self-offshore advection and partly due to the β-effect. With the inclusion of the upper ocean restoring, the dominance of anticyclonic eddy, extending from surface to a depth of ∼1200 m, can be reproduced.     

Anomalous distribution of suspended matter and some chemical compositions in seawater near the seabed: Transport processes

Vertical fine distributions of suspended matter and some chemical compositions in seawater were measured in the layer near the seabed. Distributions of organic suspended matter were almost uniform throughout the layer, but others were anomalous. However, with respect to their peak heights and distribution patterns, the tendency of profiles was the same on the whole. We suggested that these anomalous distributions were primarily caused by the destruction of the large flocculent particle on the seabed and by upward diffusion of disintegrated particles. The similarity of profiles between chemical compositions of dissolved substances and suspended matter of non-dissolved substances, was explained by assuming that pore water in the large particle was simultaneously transported with disintegrated particles after the destruction. An equation for the vertical distribution of suspended matter near the seabed was derived, provided that the rate of destruction of large flocculent particle was in proportion to the current velocity on the seabed. The equation represented the existence of the anomalous distribution, which was continuously changing its pattern. Measurements of vertical profiles of suspended matter showed almost the same tendency with the theory. From the characteristics of the theoretical equation, it was expected that the eddy diffusivity near the seabed was 1050 cm²/sec.     

Moho depth variation beneath southwestern Japan revealed from the velocity structure based on receiver function inversion

The Philippine Sea plate is subducting under the Eurasian plate beneath the Chugoku-Shikoku region, southwestern Japan. We have constructed depth contours for the continental and oceanic Mohos derived from the velocity structure based on receiver function inversion. Receiver functions were calculated using teleseismic waveforms recorded by the high-density seismograph network in southwestern Japan. In order to determine crustal velocity structure, we first improved the linearized time-domain receiver function inversion method. The continental Moho is relatively shallow ( 30 km) at the coastline of the Sea of Japan and at the Seto Inland Sea, and becomes deeper–greater than 40 km–around 35°N and 133.8°E. Near the Seto Inland Sea, a low-velocity layer of thickness 10 km lies under the continental Moho. This low-velocity layer corresponds to the subducting oceanic crust of the Philippine Sea plate. The oceanic Moho continues to descend from south to northwest and exhibits complicated ridge and valley features. The oceanic Moho runs around 25 km beneath the Pacific coast and 45 km beneath the Seto Inland Sea, and it extends to at least to 34.5°N. The depth variation of the Moho discontinuities is in good qualitative agreement with the concept of isostasy. From the configurations of both the continental and oceanic Mohos, we demonstrate that the continental lower crust and the subducting oceanic crust overlap beneath the southern and central part of Shikoku and that a mantle wedge may exist beneath the western and eastern part of Shikoku. The southern edge of the overlapping region coincides with the downdip limit of the slip area of a megathrust earthquake.     

Estimation prospects of the source number density of ultra-high-energy cosmic rays

We discuss the possibility of accurately estimating the source number density of ultra-high-energy cosmic rays (UHECRs) using small-scale anisotropy in their arrival distribution. The arrival distribution has information on their source and source distribution. We calculate the propagation of UHE protons in a structured extragalactic magnetic field (EGMF) and simulate their arrival distribution at the Earth using our previously developed method. The source number density that can best reproduce observational results by Akeno Giant Air Shower Array is estimated at about 10⁻⁵ Mpc⁻³ in a simple source model. Despite having large uncertainties of about one order of magnitude, due to small number of observed events in current status, we find that more detection of UHECRs in the Auger era can sufficiently decrease this so that the source number density can be more robustly estimated. Two hundred event observation above 4 × 10¹⁹ eV in a hemisphere can discriminate between 10⁻⁵ and 10⁻⁶ Mpc⁻³. Number of events to discriminate between 10⁻⁴ and 10⁻⁵ Mpc⁻³ is dependent on EGMF strength. We also discuss the same in another source model in this paper.     

Neutron star models

Neutron star models are calculated using an equation of state discussed in an earlier paper. A maximum mass for a neutron star of 1.74 solar masses is found. The central density of thie star is 3.3×10¹⁵ g/cm³. The lightest stars have masses of 0.02 (resp. 0.03) solar masses with central densities 2.2×10¹⁴ g/cm³ (resp. 1.9×10¹⁴).