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.     

Bayesian and Neural Network Approaches to Estimate Deep Temperature Distribution for Assessing a Supercritical Geothermal System: Evaluation Using a Numerical Model

The temperature distribution at depth is a key variable when assessing the potential of a supercritical geothermal resource as well as a conventional geothermal resource. Data-driven estimation by a machine-learning approach is a promising way to estimate temperature distributions at depth in geothermal fields. In this study, we developed two methodologies—one based on Bayesian estimation and the other on neural networks—to estimate temperature distributions in geothermal fields. These methodologies can be used to supplement existing temperature logs, by estimating temperature distributions in unexplored regions of the subsurface, based on electrical resistivity data, observed geological/mineralogical boundaries, and microseismic observations. We evaluated the accuracy and characteristics of these methodologies using a numerical model of the Kakkonda geothermal field, Japan, where a temperature above 500 °C was observed below a depth of about 3.7 km. When using geological and geophysical knowledge as prior information for the machine learning methods, the results demonstrate that the approaches can provide subsurface temperature estimates that are consistent with the temperature distribution given by the numerical model. Using a numerical model as a benchmark helps to understand the characteristics of the machine learning approaches and may help to identify ways of improving these methods.

     

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.     

Trace Fossils and Sulfur Isotopes in Mudstones around the Kuroko Deposits in the Hokuroku Basin, Northeast Japan: An Attempt to Delineate the Depositional Environment

Abstract. A comprehensive investigation was carried out on the distribution of both trace fossils and sulfur isotopes in mud-stones in the Hokuroku district, northeast Japan, in the hope of delineating the depositional environment of the mudstones in which the Kuroko deposits are embedded. The mudstones are generally massive in structure and usually contain large trace fossils, being indicative of an aerobic biofacies. On the other hand, some mudstones in and above the Kuroko ore horizon are partly laminated and usually contain smaller trace fossils, being assignable to an anaerobic or dysaerobic biofacies. The δ³⁴S values of sulfides in the mudstones above and below the ore horizon range from -40 to -12 %o, indicating mostly oxic depositional conditions in equilibrium with the inferred aerobic biofacies. In the mudstones in the ore horizon, the δ³⁴S values exhibit regionally discriminated variations: -44 to -12 %o in areas far (>1 km) from the known Kuroko deposits and -24 to +6 %o in areas closer to them. The latter high δ³⁴S group implies the temporal occurrence of local anoxic basins in the vicinity of the known Kuroko deposits. At the time of late Nishikurosawa Stage (i.e. the currently assumed Kuroko metallogenic epoch), an intense oceanic stagnation is suggested to have taken place to form the local anoxic basins responsible for the formation and preservation of Kuroko deposits. This oceanic environmental event is considered to be most likely due to increasing biological productivity primarily triggered and enhanced by upwelling of NADW in the paleo-Sea of Japan at that time.     

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.     

Paleoceanography and Heavy Metal Enrichment of Mudstones in the Green Tuff Region, Northeastern Japan

Abstract. A series of mudstones corresponding to N.8 to N.17 of Blow from the Uyashinai Mudstone Member, the Onnagawa Formation and the Funakawa Formation in the Taiheizan area of the central Green Tuff region, northeast Japan, was examined by chemical analysis. According to conventional chemical systematics, the mudstones are analyzed in terms of detrital, biogenic-A (siliceous), biogenic-B (calcareous) and hydrogenous components. The relative contribution of the biogenic-A (siliceous) component increases in upward succession in the Uyashinai and Onnagawa mudstones, whereas the contribution of the biogenic-B (calcareous) component is restricted to the lower Uyashinai mudstones. The contribution of hydrogenous components Zn, Pb, Cu and Ba tends to increase during the Nishikurosawa stage, and decreases near the Nishikurosawa/Onnagawa boundary (Pb and Ba) or in the Onnagawa stage (Zn and Cu). The observed enrichment of heavy metals in the Uyashinai mudstones is attributed to the onset of vigorous upwelling of deep water with higher nutrient and metal content associated with a contemporaneous change in global deep-water circulation.     

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.