Experimental constraints on the effects of pressure and H2O on the fractional crystallization of high-Mg island arc basalt

H₂O-undersaturated melting experiments of synthesized basalt (SiO₂ = 50.7 wt.%, MgO = 8.3 wt.%, Mg# = 60) were conducted at fO₂ corresponding to NNO+1 and NNO−1 to clarify the effects of pressure (2–7 kbar) and H₂O on fractional crystallization in island arcs. H₂O content was ranged from nominally anhydrous to 4.4 wt.%. Differentiation trends, namely the liquid lines of descent, change sensitively according to pressure-H₂O relations. Tholeiitic differentiation trends are reproduced with H₂O ≤ ∼2 wt.% in primary magma. With such quantities of H₂O, fractional crystallization is controlled by olivine + plagioclase at 2 kbar. Increasing the pressure from 2 to ≥4 kbar induces early crystallization of orthopyroxene instead of olivine and therefore SiO₂ enrichment in the residual melts is suppressed. Increasing H₂O (≥ ∼2 wt.% in primary magma) stabilizes clinopyroxene relative to orthopyroxene and/or magnetite. Although the phase relations and proportions strongly depend on fO₂ and H₂O content, differentiation trends are always calc-alkaline.     

An Ecosystem Model Coupled with Nitrogen-Silicon-Carbon Cycles Applied to Station A7 in the Northwestern Pacific

A model based on that of Kishi et al. (2001) has been extended to 15 compartments including silicon and carbon cycles. This model was applied to Station A7 off Hokkaido, Japan, in the Northwestern Pacific. The model successfully simulated the observations of: 1. a spring bloom of diatoms; 2. large seasonal variations of nitrate and silicate concentrations in the surface water; and 3. large inter-annual variations in chlorophyll-a. It also reproduced the observed features of the seasonal variations of carbon dioxide partial pressure (pCO₂)—a peak in pCO₂ in winter resulting from deep winter convection, a rapid decrease in pCO₂ as a result of the spring bloom, and an almost constant pCO₂ from summer through fall (when the effect of increasing temperature cancels the effect of biological production). A comparison of cases with and without silicate limitation shows that including silicate limitation in the model results in: 1. decreased production by diatoms during summer; and 2. a transition in the dominant phytoplankton species, from diatoms to other species that do not take up silicate. Both of these phenomena are observed at Station A7, and our results support the hypothesis that they are caused by silicate limitation of diatom growth. This revised version was published online in July 2006 with corrections to the Cover Date.     

An energy conversion method in the ocean using the density difference of water

A new method which produces energy from the ocean by utilizing the density difference of water, by means of a “chimney effect”, is proposed. Density difference of water in the ocean occurs in two ways, namely differences of consistency and water temperature. For instance, fresh river and melting floes and icebergs are pointed out as some origins of the former, while thermal effects of volcanoes and hot springs may account for the latter. Furthermore, sunlight power transmitted into the sea through glass fibre power lines can be an artificial source in the latter case.This paper concentrates on regions near estuaries where low density river water can be piped into the high density ocean. An analysis estimates that the proposed method extracts energy of several thousands or more kW in the above case.     

Full-range equation for wave boundary layer thickness

Full-range equation covering all the flow regimes in a wave boundary layer is proposed for the boundary layer thickness. The results are compared with the available experimental data and good agreement has been found. In case of wave boundary layers, there are three definitions of boundary layer thickness in use. Therefore, the full-range equation is derived for three of the definitions. The findings of this study may be useful in calculating suspended sediment transport in coastal environments and studying wave–current combined motion.     

Air-Sea Gas Transfer Velocity in Stormy Winter Estimated from Radon Deficiency

Thirteen vertical profiles of ²²⁶Ra and ²²²Rn in the near-surface water were obtained in the western North Pacific in winter, and the gas transfer velocities across the air-sea interface were estimated. The transfer velocities found by applying a steady state model varied widely from 2.1 to 30.2 m day⁻¹ with a mean of 9.4 m day⁻¹. The mean value is almost 5 times higher than that in summer in other oceans, and the maximum value is a record high for world oceans. This is partly due to the inadequacy of the steady state model, which overestimates when stronger winds blow in more recent days than the ²²²Rn half-life of about 4 days. In fact, a strong low pressure zone passed through the station about 2 days earlier, which was one of the low pressure zones that with a period of develop once a week or so in the northwestern North Pacific in winter. Instead of steady-state removal, if half of the radon removal occurred sporadically every 7 days, and the last removal took place two days before the observation, the transfer velocity would be 26 m day⁻¹. Our mean transfer velocity, which is less than 20% different from the steady state value including both overestimated and underestimated values, 9.4 ± 4.8 m day⁻¹, seems to represent the mean state of this region in winter. This suggests that the gas exchange fluxes under extremely rough conditions in the open ocean are larger than those estimated by using a transfer velocity equation with a linear or quadratic relationship with wind speed. This revised version was published online in July 2006 with corrections to the Cover Date.     

A Reconstruction of Observed Profiles in the Sea East of Japan Using Vertical Coupled Temperature-Salinity EOF Modes

It is important to estimate hard-to-observe parameters in the ocean interior from easy-to-observe parameters. This study therefore demostrates a reconstruction of observed temperature and salinity profiles of the sea east of Japan (30°≈40°N, 140°≈150°E). The reconstruction was done by estimating suboptimal state from several values of the observed profiles and/or sea surface dynamic height (SDH) calculated from the profiles. The estimation used a variational method with vertical coupled temperature-salinity empirical orthogonal function (EOF) modes. Profiles of temperature and salinity in the subtropical region are effectively reconstructed from in situ temperature profile data, or sea surface temperature (SST) and SDH. For example, the analyzed temperature field from SST and SDH has an accuracy to within 1°C in the subtropical region. Salinity in the sea north of Kuroshio, however, is difficult to estimate because of its complex variability which is less correlated with temperature than in the subtropical region. Sea surface salinity is useful to estimate the subsurface structure. We also show the possibility that the estimation is improved by considering nonlinearity in the equation calculating SDH from temperature and salinity analysis values in order to examine the misfit between analysis and observation. Analysis using TOPEX/POSEIDON altimetry data instead of SDH was also performed. This revised version was published online in July 2006 with corrections to the Cover Date.     

Preconditioned Optimizing Utility for Large-dimensional analyses (POpULar)

I present the derivation of the Preconditioned Optimizing Utility for Large-dimensional analyses (POpULar), which is developed for adopting a non-diagonal background error covariance matrix in nonlinear variational analyses (i.e., analyses employing a non-quadratic cost function). POpULar is based on the idea of a linear preconditioned conjugate gradient method widely adopted in ocean data assimilation systems. POpULar uses the background error covariance matrix as a preconditioner without any decomposition of the matrix. This preconditioning accelerates the convergence. Moreover, the inverse of the matrix is not required. POpULar therefore allows us easily to handle the correlations among deviations of control variables (i.e., the variables which will be analyzed) from their background in nonlinear problems. In order to demonstrate the usefulness of POpULar, we illustrate two effects which are often neglected in studies of ocean data assimilation before. One is the effect of correlations among the deviations of control variables in an adjoint analysis. The other is the nonlinear effect of sea surface dynamic height calculation required when sea surface height observation is employed in a three-dimensional ocean analysis. As the results, these effects are not so small to neglect.     

Estimating the geostrophic velocity obtained by HF radar observations in the upstream area of the Kuroshio

A method to extract geostrophic current in the daily mean HF radar data in the Kuroshio upstream region is established by comparison with geostrophic velocity determined from the along-track altimetry data. The estimated Ekman current in the HF velocity is 1.2% (1.5%) and 48° (38°)-clockwise rotated with respect to the daily mean wind in (outside) the Kuroshio. Furthermore, additional temporal smoothing is found necessary to remove residual ageostrophic currents such as the inertial oscillation. After removal of the ageostrophic components, the HF geostrophic velocity agrees well with that from the altimetry data with rms difference 0.14 (0.12) m/s in (outside) the Kuroshio.