Orbital evolution and accretion of protoplanets tidally interacting with a gas disk: I. Effects of interaction with planetesimals and other protoplanets

We have performed N-body simulations on the stage of protoplanet formation from planetesimals, taking into account so-called “type-I migration,” and damping of orbital eccentricities and inclinations, as a result of tidal interaction with a gas disk without gap formation. One of the most serious problems in formation of terrestrial planets and jovian planet cores is that the migration time scale predicted by the linear theory is shorter than the disk lifetime (10⁶-10⁷ years). In this paper, we investigate retardation of type-I migration of a protoplanet due to a torque from a planetesimal disk in which a gap is opened up by the protoplanet, and torques from other protoplanets which are formed in inner and outer regions. In the first series of runs, we carried out N-body simulations of the planetesimal disk, which ranges from 0.9 to 1.1 AU, with a protoplanet seed in order to clarify how much retardation can be induced by the planetesimal disk and how long such retardation can last. We simulated six cases with different migration speeds. We found that in all of our simulations, a clear gap is not maintained for more than 10⁵ years in the planetesimal disk. For very fast migration, a gap cannot be created in the planetesimal disk. For migration slower than some critical speed, a gap does form. However, because of the growth of the surrounding planetesimals, gravitational perturbation of the planetesimals eventually becomes so strong that the planetesimals diffuse into the vicinity of the protoplanets, resulting in destruction of the gap. After the gap is destroyed, close encounters with the planetesimals rather accelerate the protoplanet migration. In this way, the migration cannot be retarded by the torque from the planetesimal disk, regardless of the migration speed. In the second series of runs, we simulated accretion of planetesimals in wide range of semimajor axis, 0.5 to 2-5 AU, starting with equal mass planetesimals without a protoplanet seed. Since formation of comparable-mass multiple protoplanets (“oligarchic growth”) is expected, the interactions with other protoplanets have a potential to alter the migration speed. However, inner protoplanets migrate before outer ones are formed, so that the migration and the accretion process of a runaway protoplanet are not affected by the other protoplanets placed inner and outer regions of its orbit. From the results of these two series of simulations, we conclude that the existence of planetesimals and multiple protoplanets do not affect type-I migration and therefore the migration shall proceed as the linear theory has suggested.     

Bidirectional shaking table tests of a low-cost friction sliding system with flat-inclined surfaces

A novel low-cost friction sliding system for bidirectional excitation is developed to improve the seismic performance of reinforced concrete (RC) bridge piers. The sliding system is a spherical prototype developed by combining a central flat surface with an inclined spherical segment, characterized by stable oscillation and a large reduction in response accelerations on the flat surface. The inclined part provides a restoring force that limits the residual displacements of the system. Conventional steel and concrete are employed to construct a flat-inclined spherical surface atop an RC pier. The seismic forces are dissipated through the frictions generated during the sliding movements; hence, the seismic resilience of bridges can be ensured with a low-cost design solution. The proposed system is fabricated utilizing a mold created by a three-dimensional printer, which facilitates the use of conventional concrete to construct spherical shapes. The concrete surface is lubricated with a resin material to prevent abrasion from multiple input ground motions. To demonstrate the effectiveness of the system, bidirectional shaking table tests are conducted in the longitudinal and transverse directions of a scaled bridge model. The effect of the inclination angle and the flat surface size is investigated. The results demonstrate a large decrease in response acceleration when the system exhibits circular sliding displacement. Furthermore, the inclination angle that generates the smallest residual displacement is identified experimentally.     

Cathodoluminescence of high-pressure feldspar minerals as a shock barometer

Cathodoluminescence (CL) analyses were carried out on maskelynite and lingunite in L6 chondrites of Tenham and Yamato-790729. Under CL microscopy, bright blue emission was observed in Na-lingunite in the shock veins. Dull blue-emitting maskelynite is adjacent to the shock veins, and aqua blue luminescent plagioclase lies farther away. CL spectroscopy of the Na-lingunite showed emission bands centered at ~330, 360–380, and ~590 nm. CL spectra of maskelynite consisted of emission bands at ~330 and ~380 nm. Only an emission band at 420 nm was recognized in crystalline plagioclase. Deconvolution of CL spectra from maskelynite successfully separated the UV–blue emission bands into Gaussian components at 3.88, 3.26, and 2.95 eV. For comparison, we prepared K-lingunite and experimentally shock-recovered feldspars at the known shock pressures of 11.1–41.2 GPa to measure CL spectra. Synthetic K-lingunite has similar UV–blue and characteristic yellow bands at ~550, ~660, ~720, ~750, and ~770 nm. The UV–blue emissions of shock-recovered feldspars and the diaplectic feldspar glasses show a good correlation between intensity and shock pressure after deconvolution. They may be assigned to pressure-induced defects in Si and Al octahedra and tetrahedra. The components at 3.88 and 3.26 eV were detectable in the lingunite, both of which may be caused by the defects in Si and Al octahedra, the same as maskelynite. CL of maskelynite and lingunite may be applicable to estimate shock pressure for feldspar-bearing meteorites, impactites, and samples returned by spacecraft mission, although we need to develop more as a reliable shock barometer.     

A new tree-ring date for the “floating island” lava flow,Mount St. Helens,Washington

Anomalously narrow and missing rings in trees 12 m from Mount St. Helens' floating island lava flow, and synchronous growth increases in trees farther from the flow margin, are evidence that this andesitic flow was extruded between late summer 1799 and spring 1800 a.d., within a few months after the eruption of Mount St. Helens' dacitic layer T tephra. For ease of reference, we assign here an 1800 a.d. date to this flow. The new date shows that the start of Mount St. Helens' Goat Rocks eruptive period (1800–1857 a.d.) resembled the recent (1980–1986) activity in both petrochemical trends and timing. In both cases, an initial explosive eruption of dacite was quickly succeeded by the eruption of more mafic lavas; dacite lavas then reappeared during an extended concluding phase of activity. This behavior is consistent with a recently proposed fluid-dynamic model of magma withdrawal from a compositionally zoned magma chamber.     

Extracellular release of photosynthetic products by a pelagic blue-green alga,Trichodesmium thiebautii

Extracellular release of photosynthetic products by a pelagic blue-green alga,Trichodesmium thiebautii, was scrutinized in relation to pre-treatments in¹⁴C-method for measuring the primary productivity.The extracellular release increased markedly by the pre-treatments with such strong fixative solutions as formaldehyde and mercuric chloride. In such a case, the amount of extracellular fraction reached even almost 78 % of the total photosynthates. The extracellular release of photosynthates was dependent both on light intensity and on physiological state of algal cells. When photosynthesis was terminated by a fixative, the extracellular release was noticed even just after the fixation, and it increased with time. Thus, it appears to be difficult to deduce the real photosynthetic production capacity from the amount of cellular fraction only.The amount of extracellular fraction for unfixed samples filtered immediately after the incubation was less than 8 % of the total photosynthetic products inT. thiebautii. This value is nearly comparable to the amount of excretion reported in many cultured algae.     

Negligible effect of the benthic fauna on measuring the nutrient upward fluxes from coastal sediments

Although the influence of benthic fauna on estimating the nutrient upward fluxes from the coastal sediments are a concern, there are few comprehensive studies. We conducted monthly observations in a coastal embayment to examine the effect of benthic fauna on measuring the nutrient flux. Nutrient fluxes were measured by the core incubation method and the densities of organisms in benthic communities were also measured. Numerical composition of benthic communities, according to taxonomic group, was dominated by Polychaeta composed of a few Capitellida and Spionidae species. A negligible effect of the benthic communities on measuring nutrient fluxes was demonstrated.     

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.     

Bio-Optical Relationship of Case I Waters: The Difference between the Low- and Mid-Latitude Waters and the Southern Ocean

Both historic and currently operational chlorophyll algorithms of the satellite-borne ocean color sensors, such as SeaWiFS, were evaluated for in situ spectral radiation and chlorophyll data in some Case I waters, including the waters in the Indian Ocean sector of the Southern Ocean. Chlorophyll a concentration of the data set (n = 73) ranged from 0.04 to 1.01 mg m^–3. The algorithms had higher accuracy for the low- and mid-latitude waters (RMSE: 0.163–0.253), specifically the most recently developed algorithms of OCTS and Sea WiFS showed 0.163 and 0.170 of Root Mean Square Errors, respectively. However, these algorithms had large errors (0.422–0.621) for the Southern Ocean data set and underestimated the surface chlorophyll by more than a factor of 2.6. The absorption coefficients in the blue spectral region retrieved from remote sensing reflectance varied in a nonlinear manner with chlorophyll a concentration, and the value in the Southern Ocean was significantly lower than that in the low- and mid-latitude waters for each chlorophyll a concentration. The underestimation of chlorophyll a concentration in the Southern Ocean with these algorithms was caused by the lower specific absorption coefficient in the region compared with the low- and mid-latitude waters under the same chlorophyll a concentration.     

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.     

Audio frequency magneto-telluric survey of Norikura Volcano in central Japan

Norikura Volcano has not been active during the last 10,000 years in spite of the activity of the surrounding volcanic mountains. To study past volcanic activities, geological studies were carried out extensively. However, quite a few geophysical investigations were conducted to contribute to volcanology. Our objective is to detect the present subsurface structure of Norikura Volcano and to define volcanic stratifications. In the vicinity of Norikura Volcano, geothermal fields are still active. Subsurface volcanic rocks in this area have been exposed to geothermal activity and altered. To comprehend volcanic stratifications of Norikura and geothermal activity, we conducted audio frequency magneto-telluric (AMT) surveys around Norikura Volcano. AMT survey is useful in clearly defining the resistivity structure related to volcanic regions. The AMT data were acquired over a frequency range 10 Hz–10 kHz. Decomposition analysis was applied to the tensor impedance data. Subsequently, apparent resistivity and phase data were inverted using a two-dimensional magneto-telluric (MT) inversion and a model of Norikura was derived. The final model manifests that the surface resistors are in agreement with andesite lava or dacite lava. As for the deeper structure, a horizontal conductor is situated above resistive basements. The alteration of the conductor was weak, while basement rocks were strongly altered and/or heated through the thermal activity. The existence of these layers seems to indicate the degree of thermal activity of Norikura Volcano.