Fluctuation in volume transport distribution accompanied by the Kuroshio front migration in the Tokara Strait

A relation between migration of the Kuroshio front and fluctuation of distribution of volume transport in the Tokara Strait was described, using sea level records at five tide gauge stations around the strait and data which were composed of sea surface temperature, XBT casts, sea surface salinity and velocities at 20 m, 75 m and 150 m depths taken en route a ferryboat. The Kuroshio front extends to about 150 m depth. The sea surface salinity and the horizontal velocities abruptly change at the front. There is a good correlation in a period range from half a month to two months between the migration of the front, which is not only at the surface but also in the subsurface, and the sea level fluctuation at Nakano-shima. A northward migration of the front with a period range from 17 to 50 days decreases the transport in the southern strait between Naze and Nakano-shima but increases in the northern strait between Nakano-shima and Sata-misaki. The northward migration intensifies inflow into Kagoshima Bay and the Ohsumi Branch Current. Correlation between the transport in the northern strait and the Ohsumi Branch Current is significant in the period range from 30 to 50 days. In this significant period range, the former leads the latter by about 3 days.     

Response of Eucalanus bungii to oceanographic conditions in the western subarctic Pacific Ocean: Retrospective analysis of the Odate Collections

To evaluate their response to oceanographic conditions, interannual variations in seasonal abundance of Eucalanus bungii were investigated in zooplankton samples collected from the Oyashio Current system from 1960 to 2002. Large decadal changes were observed in seasonal timing and population age-structure. During the early 1970s and 1990s, E. bungii were abundant until mid-summer, but during the late 1970s and early 1980s, the season of maximum abundance was limited to spring and early summer. From the late 1970s to early 1980s, spring–summer abundance of newly recruited young copepodites (C1–C2) declined significantly, and an even more pronounced decline was observed for the abundance of the late copepodite stages (C3–C5). Monthly population structure showed that young of the year stopped development at C3 during the late 1970s to early 1980s, but molted into late copepodite stages in the other decades. Seasonal weakening of the Aleutian Low Pressure System estimated from North Pacific Index (NPI) was rapid during the late 1970s to early 1980s, and the NPI was positively correlated with phosphate concentrations at sea surface, spring–summer abundance of the young copepodites stages, and the extended duration of the season of high abundance. These results suggest that the decadal decline of copepod abundance originated at the early life stages, and was associated with a shift of atmospheric and oceanographic conditions. As possible biological mechanisms, we propose reduced egg production, lower survival for the portion of the annual cohort with late birth date, and overwintering of the survivors at younger stages.     

Two dimensional responses of fixed offshore platform in waves

Responses of a four-legged offshore platform to wave forces are analysed by taking account of one torsional motion around the vertical axis which passes the centre of gravity of the structure, and two translational motions in the horizontal plane. Applying the spring-mass system with one degree-of-freedom to each mode of motion, a calculation method for predicting resultant dynamic displacements of the topdeck is developed. Both the inline force and the transverse force are considered as external forces. The calculated results are found to agree with experiments. Furthermore, the influence of seven important factors on the dynamic response is examined by the calculations. It is shown that the combination pattern of the direction factor of the transverse force is of importance in dynamic analysis.     

Numerical Simulation of Biological Impact Caused by Direct Injection of Carbon Dioxide in the Ocean

The direct injection of CO₂ in the deep ocean is a promising way to mitigate global warming. One of the uncertainties in this method, however, is its impact on marine organisms in the near field before CO₂ is diluted widely in the ocean. Since field experiments cost enormously, computational simulations are expected to show detailed information on the dilution process near injection points and its impact on marine organisms. In general, the LC50 concept is widely applied for testing the acute impact of a toxic agent on organisms. As a biological impact model we therefore consider mortality, which reflects recent laboratory experiments on zooplankton at various concentrations of CO₂. Here we regard the sigmoid-transformed mortality as a linear function of time in the logarithmic scale, and not just of the concentration of CO₂ in the logarithmic scale. This model was installed in a computational simulation code for the reconstruction of small-scale ocean turbulence. The results suggest that the biological effect is not significant when the ship speed is 4 knots and CO₂ is injected at 0.1 ton/sec in the form of a spray through 100 nozzles provided vertically on a pipe at 10 m intervals. It is therefore considered that the moving-ship method is effective for direct CO₂ injection. This revised version was published online in July 2006 with corrections to the Cover Date.     

Population dynamics of earthquakes and mathematical modeling

We present a mathematical model that describes temporal variations of earthquakes. This model is represented as $$dn(t)/dt = n(t)\left[ {\alpha - \beta n(t) - \int_{ - \infty }^t {n(s)h(t - s)ds} } \right].$$ Heren(t) shows the numberof earthquakes per unit time in a certain region. α and β are constants. The functionh(t) denotes the hysteresis effect of the earthquake occurrences and can take the following forms depending on the physical conditions of the crusts; (A)h(t)=0: the equation represents a logistic type increase or decrease and approaches a stationary state asymptotically. This describes aftershock series of large earthquakes and earthquake swarms of large scale such as the Wakayama and Matsushiro swarms in Japan; (B)h(t)=constant (β=0): frequencyn(t) increases initially and then decreases gradually and shows some kind of volcanic swarms; (C)h(t) = κ · {exp(?γ₁ t) ? exp(γ₂ t)}, (γ₂ > γ₁): this denotes time delay effects and the model shows periodic patterns of bursts or “rhythms” of earthquakes, which are observed in earthquake swarms. When external effects are taken into consideration, the model is further generalized and can describe various seismic patterns. These effects represent various influences of the circumstances like the earth tide and fluctuations of plate motions, etc. Whenh(t) takes type (A) and the external effect is random, the equation displays repetitive random patterns with bursts. Particularly interesting cases may be those whenh(t) is type (C) and the external force is periodic like the earth tide. Various nonperiodic as well as periodic patterns of earthquakes appear. These are the phenomena of “chaos” and “entrainment”, etc. and can be commonly observed. Varieties of actual earthquake patterns seem to be, at least partly, explained by the nonlinear coupling between the tidal forces and autonomous rhythms of earthquakes.     

Estimation of lead sources in a Japanese cedar ecosystem using stable isotope analysis

Anthropogenic Pb affects the environment worldwide. To understand its effect on forest ecosystem, Pb isotope ratios were determined in precipitation, various components of vegetation, the forest floor, soil and parent material in a Japanese cedar (Cryptomeria japonica D. Don) forest stand. The average ²⁰⁶Pb/²⁰⁷Pb ratio in bulk precipitation was 1.14 ± 0.01 (mean ± SD), whereas that in the subsoil (20–130 cm) was 1.18 ± 0.01. Intermediate ratios ranging from 1.15 to 1.16 were observed in the vegetation, the forest floor, and the surface soil (0–10 cm). Using the ²⁰⁶Pb/²⁰⁷Pb ratios, the contribution of anthropogenic sources to Pb accumulated in the forest were estimated by the simple binary mixing model. Sixty-two percent of the Pb in the forest floor, 71% in the vegetation, and 55% in the surface soil (0–10 cm) originated from anthropogenic sources, but only 16% in the sub-surface soil (10–20 cm) was anthropogenic. These results suggest that internal Pb cycling occurs mainly between surface soil and vegetation in a Japanese cedar ecosystem, and that anthropogenic Pb strongly influences Pb cycling. Although the Japanese cedar ecosystem has a shallow forest floor, very little atmospherically derived Pb migrated downward over 10 cm in depth.     

Magnetotelluric resistivity modeling for 3D characterization of geothermal reservoirs in the Western side of Mt. Aso, SW Japan

Geothermal reservoirs are usually located at a depth range of 2 to 5 km, so to efficiently utilize such resources an advanced prospecting method is needed to detect these deep geologic structures. This study aimed to three-dimensionally characterize geothermal reservoirs by a combination of Magnetotelluric (MT) survey, inversion analysis of apparent resistivity, and interpolation of the resistivity data obtained. The western side of Mt. Aso crater, southwest Japan, was chosen as the case study area. Three hot springs exist there and a fault is assumed to go in the direction connecting them. A MT survey was carried out at 26 sites and the data processed by a remote reference method to reduce artificial noises. Based on skewness and Mohr circle analyses of the impedance tensor, the local geologic structure at each site could be approximated as horizontally layered and therefore, a one-dimensional inversion analysis was applied to the MT raw data. The resultant resistivity column data were then interpolated by the three-dimensional optimization principle method. The resistivity distributions obtained clarified continuous conductors with especially low resistivity (less than 10 Ω·m) at the hot springs along the fault. Because the resistivity decreases largely with an abundance of clay minerals, the conductors could be considered to correspond with the cap rocks. Thus, two geothermal reservoirs, whose shapes were estimated to be pillar, were detected under the cap rocks in an elevation range of − 1000 to − 3000 m. By comparing the resistivity distributions with the temperature distributions based on fluid-flow calculations at a steady state using FEM, the validity of the location and dimension of the estimated reservoirs were confirmed.     

Concentrations of 222 Rn, Its Short-Lived Daughters And 212 Pb And Their Ratios Under Complex Atmospheric Conditions And Topography

Atmospheric activity concentrations of ²¹²Pb and short-lived ²²²Rndaughters, together with meteorological elements, have been observed continuously atthree sites at Kamisaibara Village in Japan. In addition, atmospheric activity concentrationof ²²²Rn, equilibrium-equivalent concentration of ²²²Rn and conditionsof the lower atmosphere were observed for three intensive observation periods at Akawase,one of the three sites in Kamisaibara Village. The equilibrium-equivalent concentration of²²²Rn is almost the same as the atmospheric activity concentration of short-lived²²²Rn daughters.The activity concentrations of ²¹²Pb and the short-lived ²²²Rn daughtersand their ratio were low in the daytime owing to convective mixing, and high at nightowing to the surface-based inversion during periods of no precipitation. Their variationshave several patterns corresponding to the scale of the drainage wind or weak mixing.Mechanical mixing due to strong winds through both day and night during the first andsecond observation periods made the atmospheric activity concentrations of ²¹²Pb and the short-lived ²²²Rn daughters continuously low. However, their ratios werecontinuously high during the first period yet continuously low during the second period.This difference can be explained by the effect ofextraction of ²²⁰Rn and ²²²Rndue to strong winds and snow cover. There were also cases in which the ratio of theatmospheric activity concentration of ²¹²Pb to that of the short-lived ²²²Rndaughters at night was equal to or less than the ratio in the daytime. Thisinverse trend, asin the periods of no precipitation mentioned above, is considered to be due to near-neutralconditions on these nights.We find a difference in the ratio of the equilibrium-equivalent concentration of²²²Rn (the activity concentration of short-lived ²²²Rn daughters) tothe activity concentration of ²²²Rn during the first observation period and thatduring the second. The difference can be explained by snow cover on the ground. Wealso find differences among the ratios of the activity concentration of the short-lived²²²Rn daughters to that of ²²²Rn during the three observation periods.These differences can be explained by the submergence of paddy fields.     

Formation of the Yakuno ophiolite; accretionary subduction under medium-pressure-type metamorphic conditions

The notion that the Yakuno ophiolite and overlying Maizuru Group represents an accretionary prism formed during the Permian evolution of Japan on the Yakuno eruptive sequence, association of hemipelagic mudstone with silicic tuff, exotic fossiliferous limestones derived from previously accreted sea-mounts, upward coarsening of sequences terrigenous sandstone and conglomerate, and mildly deformed Permian and Triassic forearc basin formations. The most important indicator, however, is the seaward imbrication and repetition observed in both the Maizuru Group and the ophiolite itself. D1 deformation structures include axial–planar foliations (pressure-solution cleavage for the Maizuru Group and granulite–amphibolite metamorphic layering in the ophiolite), flattening type strain, symmetric pressure shadows and fringes, and isoclinal folds showing axial–planar foliations and thrust faulting at their overturned limb. The exceptional asymmetry observed indicates seaward-directed shearing near the thrust, while D1 structures in the Maizuru zone are explained by off-scraping, above the basal decollement. The later Jurassic D2 kink fold structure includes a first-order asymmetric kink with a brittle thrust at its overturned limb, more-or-less coeval with M2 retrograde metamorphism. Medium-pressure M1 prograde metamorphism in the Yakuno ophiolite produced layering of granulite and amphibolite, and in the Maizuru Group, formation of illite along pressure-solution cleavage of mudstones. The metamorphic grade is controlled by the stratigraphic relationships and appears typical of that in ocean floor regions. However, there was only one episode of M1 prograde metamorphism which occurred contemporaneously with D1 off-scraping. Given that subduction zones are normally characterized by high P/T metamorphic regimes, the observed P/T history appears to reflect relatively unusual conditions. Such high thermal gradients may plausibly reflect the approach of a young, hot oceanic plate which continued subducting beneath the Japanese arc. Accordingly, the Yakuno ophiolite was probably formed at the trench–trench–ridge triple junction.