Formation of terrestrial planets in a dissipating gas disk with Jupiter and Saturn

We have performed N-body simulations on final accretion stage of terrestrial planets, including the eccentricity and inclination damping effect due to tidal interaction with a gas disk. We investigated the dependence on a depletion time scale of the disk, and the effect of secular perturbations by Jupiter and Saturn. In the final stage, terrestrial planets are formed through coagulation of protoplanets of about the size of Mars. They would collide and grow in a decaying gas disk. Kominami and Ida [Icarus 157 (2002) 43-56] showed that it is plausible that Earth-sized, low-eccentricity planets are formed in a mostly depleted gas disk. In this paper, we investigate the formation of planets in a decaying gas disk with various depletion time scales, assuming disk surface density of gas component decays exponentially with time scale of τ_gas. Fifteen protoplanets with are initially distributed in the terrestrial planet regions. We found that Earth-sized planets with low eccentricities are formed, independent of initial gas surface density, when the condition (τ_cross+τ_growth)/2?τ_gas?τ_cross is satisfied, where τ_cross is the time scale for initial protoplanets to start orbit crossing in a gas-free case and τ_growth is the time scale for Earth-sized planets to accrete during the orbit crossing stage. In the cases satisfying the above condition, the final masses and eccentricities of the largest planets are consistent with those of Earth and Venus. However, four or five protoplanets with the initial mass remain. In the final stage of terrestrial planetary formation, it is likely that Jupiter and Saturn have already been formed. When Jupiter and Saturn are included, their secular perturbations pump up eccentricities of protoplanets and tend to reduce the number of final planets in the terrestrial planet regions. However, we found that the reduction is not significant. The perturbations also shorten τ_cross. If the eccentricities of Jupiter and Saturn are comparable to or larger than present values (∼0.05), τ_cross become too short to satisfy the above condition. As a result, eccentricities of the planets cannot be damped to the observed value of Earth and Venus. Hence, for the formation of terrestrial planets, it is preferable that the secular perturbations from Jupiter and Saturn do not have significant effect upon the evolution. Such situation may be reproduced by Jupiter and Saturn not being fully grown, or their eccentricities being smaller than the present values during the terrestrial planets' formation. However, in such cases, we need some other mechanism to eliminate the problem that numerous Mars-sized planets remain uncollided.     

A survey of possible impact structures on 25143 Itokawa

We determined the morphologies and dimensions of possible impact craters on the surface of Asteroid 25143 Itokawa from images taken by the Hayabusa spacecraft. Circular depressions, circular features with flat floors or convex floors, and circular features with smooth surfaces were identified as possible craters. The survey identified 38 candidates with widely varying morphologies including rough, smooth and saddle-shaped floors, a lack of raised rims and fresh material exposures. The average depth/diameter ratio was 0.08±0.03: these craters are very shallow relative to craters observed on other asteroids. These shallow craters are a result of (1) target curvature influencing the cratering process, (2) raised rim not being generated by this process, and (3) fines infilling the craters. As many of the crater candidates have an unusual appearance, we used a classification scheme that reflects the likelihood of an observed candidate's formation by a hypervelocity impact. We considered a variety of alternative interpretations while developing this scheme, including inherited features from a proto-Itokawa, spall scars created by the disruption of the proto-Itokawa, spall scars following the formation of a large crater on Itokawa itself, and apparent depressions due to random arrangements of boulders. The size-frequency distribution of the crater candidates was close to the empirical saturation line at the largest diameter, and then decline with decreasing diameter.     

Application of a 3-D chemical fate prediction model (FATE3D) to predict dioxin concentrations in the Tokyo Bay

A 3-D chemical fate prediction model (FATE3D) was applied to predict the dioxin concentrations in the seawater of Tokyo Bay, Japan. The simulations were carried out for a period of one year (from September 2002 to August 2003). Parameters such as meteorological data, flow field conditions, concentrations and sinking rates of organic particulate matter, initial and boundary conditions, and loading fluxes and physico-chemical properties of dioxins were used as the model inputs.The simulation results compared favorably with the field measurements of dioxin concentrations in the bay for both the particulate and dissolved phases, indicating the validity and predictive capability of the model. Furthermore, the differences in the seasonal cycles and distributions between the particulate- and dissolved-phase dioxins in the bay were estimated from the simulation results.However, the particulate-phase dioxin concentrations in the bottom layers (+1 m from the bottom) were underestimated, probably because the resuspension process was not taken into account in the model. The improvement of the model's predictive capability, including the resuspension process, shall be the focus of our next study.     

The galactic warp and rotations of the fundamental system

The proper motion in galactic latitude of O-B stars enables us to detect the kinematic behaviour of an optical counterpart of the large-scale warp of the HI gas layer in our Galaxy. A selected set of the proper motions of about 350 O-B stars within 3kpc from the sun (R₀=8.5kpc) is analyzed on the proper motion systems of N30, FK4, and FK5. A remarkable differece in the kinematic behaviour of the warp appears between the old systems (N30 and FK4) and FK5-system. On the old systems, the O-B stars in the belt 8.5kpcR<9.5kpc exhibit a systematic z-motion upward from the galactic plane forl180° and downward forl>180° with the mean proper motions of about ±0".4/century, respectively. On the other hand, the results on the FK5-system show no meaningful systematic z-motion, even though the O-B star layer exterior to the solar circle is inclined (3°) with respect to the galactic plane. These findings can neither be inferred from the model of the oblique material flow nor from the concepts of the precessional stellar rings and of the bending oscillation of a stellar disk. The remarkable difference in the kinematic behaviour of the warp, appearing between the old and new systems, is caused mainly by the conversion of the proper motions on the old systems into those on the J2000.0 frame. The conversion near the galactic plane is given by µ_b(FK4(J2000.0))–µ_b(FK4)–0.50 sinl/century. The implication of this relation is discussed in connection with the warping motion of stars detected here.     

Partitioning of elements between majorite garnet and melt and implications for petrogenesis of komatiite

Partitioning of elements between majorite garnet and ultrabasic melt has been studied at 16 GPa and 1950° C. Ca, Ti, La, Sm, Gd, Zr, Hf, Fe, Ni, Mn, K, and Na are enriched in the melt, whereas Al, Cr, V, Sc and Yb are concentrated in majorite garnet. Thus, majorite garnet fractionation by partial melting could produce chemical heterogeneities in these elements deviating from chondritic abundance. Using the partitioning behaviour of elements between majorite garnet and ultrabasic melt, the petrogenesis of komatiite is discussed. A simple model to explain the chemical varieties of komatiites is as follows. Aluminadepleted komatiite was generated by partial melting of the primitive mantle at 200–650 km depth, and alumina-enriched komatiite is the product of remelting of the residual solid at the same depths, whereas alumina-undepleted komatiite was formed by partial melting of the primitive upper mantle at depths shallower than 200 km. We suggest the possibility of large-scale chemical layering or heterogeneity in the early Archean upper mantle as an alternative model for komatiite genesis; shallower mantle depleted in majorite garnet and the underlying mantle enriched in majorite garnet. Alumina-depleted and alumina-enriched komatiites in the early Archean might be generated by a high degree of partial melting of the layered mantle. Such chemical layering could have been homogenized by the late Archean. This explains the observations that alumina-depleted and alumina-enriched komatiites were generally formed in the early Archean but alumina-undepleted komatiite was erupted in the late Archean.     

Recent and historical tsunami deposits from Lake Tokotan,eastern Hokkaido,Japan, inferred from nondestructive,grain size,and radioactive cesium analyses

Geological evidence of recent tsunamis from sediment samples collected from Lake Tokotan, a coastal lagoon in eastern Hokkaido, northern Japan, was detected using computed tomography (CT) and soft X-ray images, grain size, and radionuclide profiles. Initial field observations revealed that sediments had no discernable sedimentary structures at the top of the core. However, results of CT imaging, soft X-ray, and grain size analyses show evidence for three invisible sand layers that are intercalated with mud layers. These sand layers exhibit trends of landward fining and thinning. Furthermore, the distribution of sand layers was limited to the center and seaward parts of the lake. Vertical profiles of cesium and lead concentrations in combination with recent eyewitness accounts indicated that these sand layers are correlated with the 1973 Nemuro-oki, 1960 Chilean, and 1952 Tokachi-oki tsunami events. The deeper part of the sediment cores includes three volcanic ash layers and three prehistoric coarse sand layers. The prehistoric layers are correlated with unusually large tsunamis that were geologically identified in previous studies from eastern Hokkaido. These findings suggest that nondestructive techniques, in combination with radionuclide analysis, allow for detection of frequent but faint tsunami deposits. This technique allows for an improved understanding of the history of tsunami inundation in Lake Tokotan and of other locations for which stratigraphic evidence for faint tsunamis layers is not readily apparent from field assessments.

     

Automated classifications of topography from DEMs by an unsupervised nested-means algorithm and a three-part geometric signature

An iterative procedure that implements the classification of continuous topography as a problem in digital image-processing automatically divides an area into categories of surface form; three taxonomic criteria–slope gradient, local convexity, and surface texture–are calculated from a square-grid digital elevation model (DEM). The sequence of programmed operations combines twofold-partitioned maps of the three variables converted to greyscale images, using the mean of each variable as the dividing threshold. To subdivide increasingly subtle topography, grid cells sloping at less than mean gradient of the input DEM are classified by designating mean values of successively lower-sloping subsets of the study area (nested means) as taxonomic thresholds, thereby increasing the number of output categories from the minimum 8 to 12 or 16. Program output is exemplified by 16 topographic types for the world at 1-km spatial resolution (SRTM30 data), the Japanese Islands at 270 m, and part of Hokkaido at 55 m. Because the procedure is unsupervised and reflects frequency distributions of the input variables rather than pre-set criteria, the resulting classes are undefined and must be calibrated empirically by subsequent analysis. Maps of the example classifications reflect physiographic regions, geological structure, and landform as well as slope materials and processes; fine-textured terrain categories tend to correlate with erosional topography or older surfaces, coarse-textured classes with areas of little dissection. In Japan the resulting classes approximate landform types mapped from airphoto analysis, while in the Americas they create map patterns resembling Hammond's terrain types or surface-form classes; SRTM30 output for the United States compares favorably with Fenneman's physical divisions. Experiments are suggested for further developing the method; the Arc/Info AML and the map of terrain classes for the world are available as online downloads.     

Time variation of the X-ray spectrum and optical luminosity of SCO X-1

Results of rocket observations of SCO X-1 over the spectral range of 220 keV are presented. The observations have been performed partly in India and partly in Japan under the collaboration of the three groups. The present results are compared with results of similar observations carried out by the LRL (Lawrence Radiation Laboratory) group. Some of these X-ray observations were accompanied by simultaneous optical observations. Relationships between the hardness of the X-ray spectrum and the X-ray intensity and between the hardness and the optical luminosity are compiled. The relationships among the parameters (temperature, density and size) which characterize the postulated isothermal cloud model of SCO X-1 are given. They indicate that SCO X-1 is characterized by a temperature of about 10⁷–10⁸K, a density of about 10¹⁶–10¹⁷ cm^–3 and a radius of about 10⁸–10⁹ cm respectively. We further show that the temperature is inversely correlated with the size of the source; an increase in temperature corresponds to a decrease in the radius and an increase in density.