Accretion rates of planetesimals by protoplanets embedded in nebular gas

When protoplanets growing by accretion of planetesimals have atmospheres, small planetesimals approaching the protoplanets lose their energy by gas drag from the atmospheres, which leads them to be captured within the Hill sphere of the protoplanets. As a result, growth rates of the protoplanets are enhanced. In order to study the effect of an atmosphere on planetary growth rates, we performed numerical integration of orbits of planetesimals for a wide range of orbital elements and obtained the effective accretion rates of planetesimals onto planets that have atmospheres. Numerical results are obtained as a function of planetesimals’ eccentricity, inclination, planet’s radius, and non-dimensional gas-drag parameters which can be expressed by several physical quantities such as the radius of planetesimals and the mass of the protoplanet. Assuming that the radial distribution of the gas density near the surface can be approximated by a power-law, we performed analytic calculation for the loss of planetesimals’ kinetic energy due to gas drag, and confirmed agreement with numerical results. We confirmed that the above approximation of the power-law density distribution is reasonable for accretion rate of protoplanets with 1-10 Earth masses, unless the size of planetesimals is too small. We also calculated the accretion rates of planetesimals averaged over a Rayleigh distribution of eccentricities and inclinations, and derived a semi-analytical formula of accretion rates, which reproduces the numerical results very well. Using the obtained expression of the accretion rate, we examined the growth of protoplanets in nebular gas. We found that the effect of atmospheric gas drag can enhance the growth rate significantly, depending on the size of planetesimals.     

Metamorphic Environments around the Cretaceous Kurihashi Pluton in the Kitakami Mountains, NE Japan: in Comparison to the Ganidake Pluton - Kamaishi Skarn Metasomatism

Abstract. Many granitic plutons of Early Cretaceous age are intruded on various scales in the Kitakami Mountains. The stock‐type Ganidake pluton accompanies enormous Fe‐Cu mineralization of the Kamaishi deposits, whereas the Kurihashi pluton accompanies less mineralization. To elucidate the cause of these differences, the metamorphic conditions and redox state of the contact metamorphic aureole around the Kurihashi pluton have been examined by the petrochemical study and gas analysis of the metamorphic rocks. A typical mineral assemblage in the pelitic rocks in the lowest‐grade part is biotite‐muscovite‐chlorite‐quartz‐plagio‐clase‐graphite, which occur more than 2 km away from the contact point with the Kurihashi pluton. Graphite disappears at the 1550 m point, and cordierite and garnet appear in the middle and highest‐grade parts, respectively. A typical mineral assemblage in the tuffaceous rocks in the lower‐grade part is chlorite‐actinolite‐biotite‐quartz‐plagioclase. Actinolite changes into hornblende near to the pluton. The CO₂/CH₄ ratios obtained in measurements by gas chromatography exceeds 100 in the pelitic rocks at the contact point with the pluton. The ratios decrease and become less than 0.1 with distance from the pluton. Equilibrium temperatures calculated from a garnet‐biotite pair in the pelitic rock and a hornblende‐plagioclase pair in the tuffaceous rock are 640d? and 681 d?C at the contact point, respectively. The log?o₂ values among these metamorphic aureoles estimated from the CO₂/CH₄ ratios are slightly lower than the FMQ‐buffer. Redox states of the contact metamorphic aureole are kept in an intermediate condition between oxidized magma of the Kurihashi pluton and graphite‐bearing pelitic country rocks. Judging from these metamorphic conditions around the Kurihashi pluton and from the re‐evaluation of the previous knowledge about contact metasomatism around the Ganidake pluton, the Kurihashi metamorphism has occurred at higher temperatures and dry conditions than the Ganidake metasomatism. These differences in the metamorphic conditions and presence or absence of a large limestone mass around the pluton might be the principal reasons why the Kurihashi pluton accompanies less mineralization and the Ganidake pluton accompanies gigantic Kamaishi skarn mineralization.     

Accuracy assessment of QuickBird stereo imagery

The Geographical Survey Institute of Japan has recently carried out an evaluation of the metric performance of QuickBird stereo satellite imagery. This paper describes the accuracy assessment of the sensor orientation and geopositioning phases of the study, the aim of which was twofold. First, it was desired to confirm the metric potential of QuickBird imagery for 1:25 000 scale topographic mapping. Second, a determination was to be made of the accuracy attainable from the Basic image product. The techniques of rational functions and affine bundle adjustment were employed, the former with bias compensation. The results obtained both reassert the high precision of the rational functions approach and cast doubt upon the applicability of the 3D affine model for accurate geopositioning from QuickBird imagery.     

Concentration mechanisms of iron oxides and alumina in deep weathering crusts (Goshikidai, Kagawa, Western Japan)

Concentration of Fe-oxides and alumina in weathering processes are main geological reactions for lateritization and bauxitization, respectively. In western Japan, red-coloured soil formed by weathering processes developed in many places. This soil is composed of hydrous Fe-oxide minerals, hydrous alumina minerals and other minerals. It was formed in the upper part of deep weathering crust by weathering processes under some kind of sub-tropical climate, probably in the Pliocene. One of these occurrences is observed in the upper part of Goshikidai and Konodai, west part of Takamatsu city, northeast Shikoku Island, west Japan. A deep weathering crust is distributed on wide hilly plains ranging from 250 to 400 m a.m.s.l. in the northwestern region of Takamatsu city. Original rock of the weathering crust is bronzite andesite and glassy bronzite andesite, so-called ‘sanukite’. The andesites had been weathered under some special climate, and the geological age of the weathering is the same as above. The mineral assemblage and formation mechanism are similar to those of laterite and bauxite. The weathering crust developing in this region are subdivided into the three following zones: (1) A zone, composed of hydrous Fe-oxides and metahalloysite with small amounts of gibbsite and it is associated with white veins of metahalloysite; (2) B zone, composed of hydrous Fe-oxides and metahalloysite (some material is associated with -cristobalite); and C zone, composed of metahalloysite or halloysite and -cristobalite with relict crystals of feldspar and quartz, and some material is associated with montmorillonite. Chemical analyses of the materials of the three zones show the formation mechanisms of the weathering crust.     

Fractal dimension of fault systems in Japan: Fractal structure in rock fracture geometry at various scales

Based on fault maps, whether or not the fracture geometry of rocks is self-similar, was examined by using a box-counting algorithm. The statistical self-similarity (fractal structure) of the fault fracture systems holds well at the scale of about 2 to 20 km. The fractal dimension in Japan varied from 1.05 to 1.60. The fractal dimension is about 1.5–1.6 at the central part of the Japan Arc, and decreases with distance from the center. At a smaller scale, the fractal structure also holds well in the rock fracture geometry. The fractal dimension of the North Izu Peninsula fault system (branching faults) is 1.49 at the scale of 0.625 to 10 km, the fractal dimension of rock fracture geometry at the scale order of 10^–1 to 10^–2 meters is about 1.49–1.61. The upper limit of the fractal dimension of rock fracture geometry is about 1.6, judging from the estimation of fractal dimension on actual fracture geometry of rocks. This value may impose a restraint on modeling of faulting and the fracture process of rocks.     

Hydraulic and numerical study on the generation of a subaqueous landslide-induced tsunami along the coast

By carrying out the hydraulic experiments in a one-dimensional open channel and two-dimensional basin, we clarified the process of how a landslide on a uniform slope causes the generation of a tsunami. The effect of the interactive force that occurs between the debris flow layer and the tsunami is significant in the generation of a tsunami. The continuous flow of the debris into the water makes the wave period of the tsunami short. The present experiments apply numerical simulation using the two-layer model with shear stress models on the bottom and interface, and the results are compared. The simulated debris flow shows good agreement with the measured results and ensures the rushing process into the water. We propose that the model use a Manning coefficient of 0.01 for the smooth slope and 0.015 for the rough slope, and a horizontal viscosity of 0.01 m²/s for the landslide; an interactive force of 0.2 for each layer is recommended. The dispersion effect should be included in the numerical model for the propagation from the shore.     

Statistical Features of Vein Systems in the Hishikari Epithermal Gold Deposit, Japan

Abstract. The Hishikari epithermal gold deposit consists of parallel vertical veins hosted by the Quaternary volcanic rocks (QVR) and Cretaceous sedimentary rocks (CSR) at the shallow and deep levels, respectively. In order to reveal the litholog‐ical influence on the evolution of vein systems in the Hishikari field, vein density, width, spacing and position were statistically examined using cores sampled from three boreholes. The statistical features of the veins identified were summarized as follows, (i) The QVR have lower mean vein densities than the CSR. (ii) The former has smaller deformation degrees (ratios of total vein width per length) than the latter, (iii) Cumulative frequencies of vein density obey negative‐exponential distributions in both the types of rocks, (iv) Cumulative frequencies of vein width obey power‐law distributions with exponents of about ‐1 in both the types of rocks, but consist of two segments in the CSR. (v) Cumulative frequencies of vein spacing obey lognormal distributions with coefficients of variation between 1.8 and 3.5 in both the types of rocks. And (vi) fractal dimensions of vein position are between 0.45 and 0.52 at measure lengths between 0.1 and 10 m without distinctive difference between the two types of rocks. The feature (i) indicates higher ductility of the QVR against fracturing than the CSR. The feature (ii) implies that large deformation was accumulated in rocks themselves (e.g., alteration zones) in addition to the veins examined in the QVR compared to only the veins in the CSR. The feature (iv) suggests that the vein systems in the CSR reached a stage at which the several selected large veins grew as main conduit, that potentially provided pathways for ore‐forming fluids to ascend from deep levels. The coefficients of variation of vein spacing and fractal dimensions of vein position (features v and vi) indicate clustered distribution of the veins in the field. This study might be helpful for a better understanding of lithological influence on the evolution of vein systems not only in the Hishikari field but also in other fields.     

Geochronological and petrological investigations of Miocene felsic igneous rocks in the Amakusa Islands,southwest Japan: Possible extension of the Setouchi Volcanic Belt

The opening of the Japan Sea led to the separation of southwest Japan from the Eurasian continent. Subsequent to this event, a diverse range of igneous activities occurred in southwest Japan. On the back-arc side of the region, igneous activity commenced at approximately 22 Ma and persisted for an extended period. In the trench-proximal region of southwest Japan, magmatism initiated around 15.6 Ma, immediately following the cessation of the Japan Sea opening, in correlation with the subduction of the Philippine Sea plate beneath southwest Japan. The Amakusa Islands in western Kyushu host felsic to intermediate igneous rocks with Miocene radiometric ages. There has been a debate regarding the attribution of the igneous rocks in Amakusa Island among the Miocene igneous rocks in southwest Japan. To address this issue, we conducted zircon U–Pb dating and analyzed the major- and trace-element compositions of felsic igneous rocks in the Amakusa Islands to elucidate their characteristics. The obtained U–Pb ages range from 14.5 to 14.8 Ma, suggesting contemporaneity between magmatism in the Amakusa Islands and the Setouchi Volcanic Rocks in the trench-proximal region of southwest Japan. The major and trace element compositions of the felsic igneous rocks exhibit similarities to the dacites of the Setouchi Volcanic Rocks. These findings support previous suggestions that the magmatism in the Amakusa Islands can be correlated with the Setouchi Volcanic Rocks, based on the discovery of a high-Mg andesite dike and paleo-stress analysis utilizing the direction of dikes and sills. Therefore, the Setouchi Volcanic Belt is proposed to extend further west than the previously identified Ohno volcanic rocks in eastern Kyushu. The subduction of the Shikoku Basin of the Philippine Sea plate toward western Kyushu supports the hypothesis that the Kyushu-Palau Ridge was positioned west of Kyushu at ~15 Ma.     

Dynamic loading test and simulation analysis of full‐scale semi‐active hydraulic damper for structural control

A semi‐active hydraulic damper (SHD) for a semi‐active damper system, which is useful for practical structural control especially for large earthquakes, has been developed. Its maximum damping force is set to 1 or 2 MN, and it is controlled by only 70 W of electric power. An SHD with a maximum damping force of 1 MN was applied to an actual building in 1998. This paper first presents the results of a dynamic loading test to confirm the control performance of the SHD. Next, an analytical model of SHDs (SHD model) is constructed with the same concept for two kinds of SHDs based on the test results. Through simulation analyses of the test results using the proposed SHD model, the dynamic characteristics of the SHD can be well represented within practical conditions. Simulation analyses are also carried out using a simple structure model with the SHD model. It is shown that this SHD model can be used to precisely evaluate the control effect of the semi‐active damper system and is useful in practical SHD design under the applied conditions. Copyright © 2000 John Wiley & Sons, Ltd.