In situ strength measurements on natural upper-mantle minerals

Using in situ strength measurements at pressures up to 10 GPa and at room temperature, 400, 600, and 700°C, we examined rheological properties of olivine, orthopyroxene, and chromian-spinel contained in a mantle-derived xenolith. Mineral strengths were estimated using widths of X-ray diffraction peaks as a function of pressure, temperature, and time. Differential stresses of all minerals increase with increasing pressure, but they decrease with increasing temperature because of elastic strain on compression and stress relaxation during heating. During compression at room temperature, all minerals deform plastically at differential stress of 4–6 GPa. During subsequent heating, thermally induced yielding is observed in olivine at 600°C. Neither orthopyroxene nor spinel shows complete stress relaxation, but both retain some stress even at 700°C. The strength of the minerals decreases in the order of chromian-spinel ≈ orthopyroxene > olivine for these conditions. This order of strength is consistent with the residual pressure of fluid inclusions in mantle xenoliths.     

The instability of viscous two-layer oscillatory flows

The instability of oscillatory flows in a two-layer fluid where the two layers differ in density and viscosity has been analysed using a perturbation method for long waves with special interest on effects of viscosity, time scale, density and depth of the fluid. The flow of a fluid with homogeneous density can be unstable, when the kinematic viscosity of the upper fluid layer is different from that of the lower one. Viscosity stratification results in unstable oscillatory flows. Two limiting cases of single-layer flow are also considered.     

Detection of systematic errors in XBT data and their correction

Comparison experiment between XBT of T-7 probe and CTD was conducted at 15 stations in the sea area centered on 29°N, 135°E in December 1985. There were systematic errors in XBT temperature profiles in comparison with CTD temperature profiles. The main cause of errors was attributed to an error in the free-fall speed of the XBT probes which was provided by the XBT maker. A previous equation for depth correction proposed by Heinmilleret al. (1983) could not give effective correction for our data. A new equation between the probe depth and the elapsed time from landing of the probe on the water was obtained by the method of adjusting temperature gradients of XBT profiles to those of CTD profiles. This equation agreed with the theoretical result given by Seaver and Kuleshov (1982) much better than that of Heinmilleret al. (1983). Systematic errors due to a scatter of values of the reference resistance and variation of B-constant of thermistors used in XBT also seemed to exist. After an adjustment using the temperature difference between XBT and CTD in the mixed layer with depths of about 100 m, the standard deviation of temperature difference between XBT and CTD from the surface to the depth of 750 m was 0.14°C.     

North Pacific Subtropical Mode Waters Observed in Long XBT Cross Sections along 32.5°N Line

The S/V Shoyo, of the Hydrographic Department, Japan Coast Guard, has conducted high-density expendable bathythermograph (XBT) measurements along the 32.5°N line in the North Pacific every year from 1990 to 1993 as a part of the Japanese-World Ocean Circulation Experiment (WOCE). These XBT data are analyzed here, focusing on year-to-year variations of the inventory and core layer temperature (CLT) of the North Pacific subtropical mode water (NPSTMW). Large year-to-year changes are found in the NPSTMW CLTs estimated in longitudes between 140°E and 160°E. CLT values were found of 17.4°C in 1990, 17.1°C in 1991, 17.3°C in 1992 and 17.6°C in 1993. Inspection of the wintertime westerlies over the formation area and sea surface temperature distribution revealed that this change in CLT can be qualitatively attributed to the strength of atmospheric cooling in the formation area in the previous winter. Although a large year-to-year variation of NPSTMW inventory was also found, it is hard to state any relationship between CLT and atmospheric forcing. There is a possibility that different observational seasons may affect the inventory. It has also been found that the thermocline depth in 1991 was shallower in the sea area east of 180° than in 1992 and 1993. Associated with this change, the North Pacific central mode water (NPCMW), characterized by thermostad with temperatures ranging from 14°C to 11°C, appears in the sea area east of 180° in the 1992 and 1993 cross sections. The 1993 cross section, which ranged from the Japanese coast to the west coast of North America, possessed another thermostad in the surface layer, with a temperature of about 17°C in the eastern part of the cross section, off California. This revised version was published online in August 2006 with corrections to the Cover Date.     

Introduction of accumulation-dispersal coefficient of marine organisms

Avoiding the subject for fish accumulation, the traditional view in fish population dynamics has ascribed immigration and emigration of fish to dispersal of fish. The main purpose of this paper is to find a quantity that represents the time rate of accumulation-dispersal of marine organisms, and also has some relation to the horizontal convergence of current velocity of the surrounding water. For this, the accumulation-dispersal coefficient is introduced not in the form of diffusion, but in the same form as the convergence. Since the accumulation-dispersal of organisms is a factor that changes its distribution density, all factors causing the change are first classified to locate the position occupied by the accumulation dispersal. Each factor corresponds to each coefficient appearing in a linearized equation describing the rate of change in the density, averaged over a region or a group. The immigration-emigration coefficient is divided into three coefficients of passage, accumulation dispersal and diffusion velocity. For the organisms ranging in a nearly horizontal layer, the accumulation-dispersal coefficient is shown to equal the area-averaged horizontal convergence of organismal velocity relative to land, which is the sum of the area-averaged horizontal convergences of swimming velocity relative to water and of current velocity. However, the area-averaged convergence of current velocity associated with the accumulation-dispersal coefficient for a region is shown to be somewhat different from the usual one.     

A method of simultaneous analysis of wind and temperature profile data over the water surface by use of the method of least squares

Priestley’s 1959 discussion concerning the estimations of momentum flux and of heat flux from profile observations of wind and temperature, stating that the closeness of the curve fitting is deceptive, is first critically re-examined. Then a method is proposed to estimate vertical flux of momentum and heat over the water surface, from the combined data of of wind speed and temperature are assumed to have similar log-linear profiles, and the most probable profiles are determined by applying the method of least squares simultaneously to wind and temperature data. Consequently, the most probable values of vertical flux of momentum and heat may be estimated to satisfy as much as possible observed data of both wind and temperature simultaneously. The coefficient of the linear term of the log-linear profile, which is treated as an indeterminate coefficient in this method, may be determined from each observed data as a function of stability length. By tentatively applying the method to Rider’s 1954 data it is found that the coefficient shows a characteristic behavior with the stability length.     

A New Fall-Rate Equation for T-5 Expendable Bathythermograph (XBT) by TSK

The accuracy of the manufacturer’s fall-rate equation for the T-5 Model of expendable bathythermograph (XBT) has been investigated based on about 300 collocated pairs of XBT-CTD (Conductivity-Temperature-Depth profiler) measurements in various climatological regions. We found that the equation systematically overestimates depth by about 5% for the T-5 produced by Tsurumi Seiki, Co. Ltd. (TSK), but almost no bias is associated with the T-5 produced by Sippican, Inc., in USA. The cause of this difference is not clear, because the two manufacturers’ T-5 probes are reported to have identical shape and weight in water. We propose a new fall-rate equation for the TSK T-5: z(t) = 6.54071t - 0.0018691t ², where z(t) is depth in meters at time, t, in seconds.     

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.     

Baroclinic Buoyancy-Inertia Joint Stability Parameter

At present, the barotropic buoyant stability parameter has been derived from a vertical virtual displacement of a water parcel. The barotropic inertial stability parameter in the eccentrically cyclogeostrophic, basic current field was derived in 2003 from a horizontal cross-stream virtual displacement of a parcel. By expressing acceleration of a parcel due to a virtual displacement, which is arbitrarily sloping within a vertical section across the basic current, in terms of natural coordinates, we derived the vertical component of baroclinic buoyant stability parameter B ₂ ², the horizontal component of baroclinic inertial stability parameter I ₂ ², the baroclinic joint stability parameter J ², its buoyant component B ² and its inertial component I ². B ² is far greater than I ₂ ², and when neglecting relative vorticity except for vertical shear, a downward convex curve of J ² plotted against the slope of a virtual displacement follows a trend of B ² curve. If a parcel displaces along a horizontal surface or an isopycnal surface, however, B ² vanishes, and J ² becomes equal to I ². Actual parcel is apt to displace not only along the bottom slope, but also along the sea surface and an isopycnal interfacial surface, which is approximately equivalent to an isentropic surface, preferred by lateral mixing and exchange of momentum. Such actual displacement makes B ² vanishing, and grants I ² an important role. The present analysis of I ² examining effects due to curvature and horizontal and vertical shear vorticities are useful in deepening our understanding of baroclinic instability in actual oceanic streams.     

Prediction of ocean waves based on the single-parameter growth equation of wind waves

A method for the prediction of ocean waves was developed on the basis of the single-parameter growth equation of wind waves, proposed byToba (1978) on the basis of similarity in growing wind waves. The applicability of the method to actual problems was tested by hindcasting the wave characteristics with the method, for two cases with differing time and space scales, one in Kii Channel Approach, Japan, and the other in the North Atlantic Ocean. The results showed that the present method can predict waves within an error of 1.3 m in wave heights, which ranged from 3 to 12 m.