Long nonlinear topographic planetary waves in a rotating stratified ocean

Long nonlinear topographic waves in a continuously stratified ocean with a linear bottom slope are investigated. It is shown that odd cross-channel modes are governed by the Korteweg-de Vries (K-dV) equation. The solitary waves are those of a low pressure type. The long waves are shown to be modulationally stable because of the nonlinear effect due to irrotational motion. All these results are missed if the conventional quasi-geostrophic approximation is adopted.     

Spatial Variability of Living Coccolithophore Distribution in the Western Subarctic Pacific and Western Bering Sea

Vertical distributions of coccolithophores were observed in the depth range 0–50 m in the western subarctic Pacific and western Bering Sea in summer, 1997. Thirty-five species of coccolithophores were collected. Overall, Emiliania huxleyi var. huxleyi was the most abundant taxon, accounting for 82.8% of all coccolithophores, although it was less abundant in the western Bering Sea. Maximum abundance of this species was found in an area south of 41°N and east of 175°E (Transition Zone) reaching >10,000 cells L⁻¹ in the water column. In addition to this species, Coccolithus pelagicus f. pelagicus, which accounted for 4.2% of the assemblage, was representative of the coccolithophore standing crop in the western part of the subarctic Pacific. Coccolithus pelagicus f. hyalinus was relatively abundant in the Bering Sea, accounting for 2.6% of the assemblage. Coccolithophore standing crops in the top 50 m were high south of 41°N (>241 × 10⁶ cells m⁻²) and east of 170°E (542 × 10⁶ cells m⁻²) where temperatures were higher than 12°C and salinities were greater than 34.2. The lowest standing crop was observed in the Bering Sea and Oyashio areas where temperatures were lower than 6–10°C and salinities were less than 33.0. From the coccolithophore volumes, the calcite stocks in the Transition, Subarctic, and the Bering Sea regions were estimated to be 73.0, 9.7, and 6.9 mg m⁻², respectively, corresponding to calcite fluxes of 3.6, 0.5, and 0.3 mg m⁻²d⁻¹ using Stoke's Law. This revised version was published online in July 2006 with corrections to the Cover Date.     

Four Decadal Ocean-Atmosphere Modes in the North Pacific Revealed by Various Analysis Methods

Various statistical methods (empirical orthogonal function (EOF), rotated EOF, singular value decomposition (SVD), principal oscillation pattern (POP), complex EOF (CEOF) and joint CEOF) were applied to low-pass filtered (>7 years) sea surface temperature (SST), subsurface temperature and 500 hPa geopotential height in order to reveal standing and propagating features of decadal variations in the North Pacific. Four decadal ocean-atmosphere covariant modes were found in this study. The first mode is the well-known ENSO-like mode associated with the “Pacific-North American” atmospheric pattern, showing SST variations reversed between the tropics and the extratropics. In the western tropical Pacific, subsurface temperature variations were found to be out of phase with the SST variations. The other three modes are related to the oceanic general circulation composed of the subtropical gyre, the Alaskan gyre and the subpolar gyre, respectively. The 1988/89 event in the northern North Pacific was found to be closely associated with the subtropical gyre mode, and the atmospheric pattern associated with this mode is the Arctic Oscillation. An upper ocean heat budget analysis suggests that the surface net heat flux and mean gyre advection are important to the Alaskan gyre mode. For the subpolar gyre mode, the mean gyre advection, local Ekman pumping and surface net heat flux play important roles. Possible air-sea interactions in the North Pacific are also discussed. The oceanic signals for these decadal modes occupy a thick layer in the North Pacific, so that accumulated heat content may in turn support long-term climate variations. This revised version was published online in July 2006 with corrections to the Cover Date.     

Effects of clear-cutting,meteorological, and physiological factors on evapotranspiration in the Kamabuchi experimental watershed in northern Japan

This study investigated the effects of clear-cutting and the meteorological and physiological factors on forest evapotranspiration (ET), by using the water-budget method in the Kamabuchi experimental watershed (KMB; 38° 56′ 21″ N, 140° 15′ 58″ E) in northern Japan. Meteorological and discharge data collected during no-snow periods (from June to October) from 1939 were used to compare ET in three sub-watersheds: No. 1, where the forest had been left undisturbed, and No. 2 and No. 3, where Cryptomeria japonica was planted after clear-cutting. Paired watershed experiments revealed that clear-cutting caused ET to decrease by approximately 100 mm yr⁻¹, and this reduction continued for more than 20 years, even after C. japonica was planted. ET fluctuated similarly across all watersheds, regardless of clear-cutting or planting. This fluctuation is mainly caused by solar radiation and temperature. Intrinsic water-use efficiency (iWUE) calculated using δ¹³C of tree-ring cellulose in C. japonica increased due to elevated atmospheric CO₂ concentration. We estimated annual carbon fixation in a single tree as the annual net photosynthesis (A). Subsequently, transpiration (E) was calculated from the relationship between iWUE and A. The results showed that A and E per tree increased as the tree grew older; however, the trees' responses to increasing c_a suppress the increase in ET. Moreover, the fluctuation of ET from the watershed was small compared to the fluctuation of P during the observation periods because the increase and decrease in E and interception loss complemented each other.     

Fast computation of Jacobian elliptic functions and incomplete elliptic integrals for constant values of elliptic parameter and elliptic characteristic

In order to accelerate the numerical evaluation of torque-free rotation of triaxial rigid bodies, we present a fast method to compute various kinds of elliptic functions for a series of the elliptic argument when the elliptic parameter and the elliptic characteristic are fixed. The functions we evaluate are the Jacobian elliptic functions and the incomplete elliptic integral of the second and third kinds regarded as a function of that of the first kind. The key technique is the utilization of the Maclaurin series expansion and the addition theorems with respect to the elliptic argument. The new method is around 25 times faster than the method using the incomplete elliptic integral of general kind and around 70 times faster than the method using mathematical libraries given in the latest version of Numerical Recipes.     

An earthquake swarm in the northern part of the Aso caldera and migration of their foci

A seismicity gap was observed in the early stage of the activity of the earthquake swarm which has been taken place in the northern part of the Aso caldera during January 1975. In the following stage, small earthquakes began to occur, until at last the largest earthquake occured within the seismicity gap. The b-value appreciably increased after the largest earthquake. The above-mentioned characteristics of this earthquake swarm may be a typical example of the fracturing phenomena in relatively heterogeneous rocks. From the estimation of the ultimate strain energy that can be stored up in the seismicity gap, it may be concluded that the above-mentioned largest earthquake took place by the release of ultimate strain energy stored within the seismicity gap.     

New Indices for Characterizing Spatial Models of Ore Deposits by the Use of a Sensitivity Vector and an Influence Factor

Precise spatial estimation of ore grades and impurity contents from sample data limited in amount and location is indispensable to metallic and nonmetallic resource exploration. One of the advantages of using geostatistics for this purpose is that it can incorporate multivariate data into spatial estimation of one variable. However, there are two weak points concerning technical and post-processing problems. First is the difficulty in application to geologic data in which spatial correlations are not clear because of intrinsic nonlinear behavior. Second is the absence of indices to interpret the mechanisms and factors which govern the spatial distribution. To address these problems, a spatial method of modeling based on a feedforward neural network, SLANS, which recognizes the relationship between the data value and location by considering supplementary attributes such as lithology and biostratigraphy, and a sensitivity analysis using this network were developed. These methods were applied to two case studies, genetic mechanisms of kuroko deposits and quality assessment of a limestone mine. The first case study is a spatial analysis of principal metals of kuroko deposits (volcanogenic massive sulfide deposits) in the Hokuroku district, northern Japan. It was clarified that upward and downward sensitivity vectors were distinguished near the deposits inside and outside the tectonic basin, respectively. Sensitivity analysis for the second case study showed a strong effect of crystalline limestone on the important impurity, P₂O₅ contents. Hydrothermal alteration, which could cause leaching and secondary concentration of phosphorus, is considered to have produced this effect.     

Transformation from Cartesian to Geodetic Coordinates Accelerated by Halley’s Method

By using Halley’s third-order formula to find the root of a non-linear equation, we develop a new iterative procedure to solve an irrational form of the “latitude equation”, the equation to determine the geodetic latitude for given Cartesian coordinates. With a limit to one iteration, starting from zero height, and minimizing the number of divisions by means of the rational form representation of Halley’s formula, we obtain a new non-iterative method to transform Cartesian coordinates to geodetic ones. The new method is sufficiently precise in the sense that the maximum error of the latitude and the relative height is less than 6 micro-arcseconds for the range of height, −10 km ≤ h ≤ 30,000 km. The new method is around 50% faster than our previous method, roughly twice as fast as the well-known Bowring’s method, and much faster than the recently developed methods of Borkowski, Laskowski, Lin and Wang, Jones, Pollard, and Vermeille.