From planetesimals to terrestrial planets: N-body simulations including the effects of nebular gas and giant planets

We present results from a suite of N-body simulations that follow the formation and accretion history of the terrestrial planets using a new parallel treecode that we have developed. We initially place 2000 equal size planetesimals between 0.5 and 4.0 AU and the collisional growth is followed until the completion of planetary accretion (>100 Myr). A total of 64 simulations were carried out to explore sensitivity to the key parameters and initial conditions. All the important effect of gas in laminar disks are taken into account: the aerodynamic gas drag, the disk-planet interaction including Type I migration, and the global disk potential which causes inward migration of secular resonances as the gas dissipates. We vary the initial total mass and spatial distribution of the planetesimals, the time scale of dissipation of nebular gas (which dissipates uniformly in space and exponentially in time), and orbits of Jupiter and Saturn. We end up with 1-5 planets in the terrestrial region. In order to maintain sufficient mass in this region in the presence of Type I migration, the time scale of gas dissipation needs to be 1-2 Myr. The final configurations and collisional histories strongly depend on the orbital eccentricity of Jupiter. If today’s eccentricity of Jupiter is used, then most of bodies in the asteroidal region are swept up within the terrestrial region owing to the inward migration of the secular resonance, and giant impacts between protoplanets occur most commonly around 10 Myr. If the orbital eccentricity of Jupiter is close to zero, as suggested in the Nice model, the effect of the secular resonance is negligible and a large amount of mass stays for a long period of time in the asteroidal region. With a circular orbit for Jupiter, giant impacts usually occur around 100 Myr, consistent with the accretion time scale indicated from isotope records. However, we inevitably have an Earth size planet at around 2 AU in this case. It is very difficult to obtain spatially concentrated terrestrial planets together with very late giant impacts, as long as we include all the above effects of gas and assume initial disks similar to the minimum mass solar nebular.     

Primary melt from Sannome-gata volcano,NE Japan arc: constraints on generation conditions of rear-arc magmas

The conditions under which rear-arc magmas are generated were estimated using primary basalts from the Sannome-gata volcano, located in the rear of the NE Japan arc. Scoriae from the volcano occur with abundant crustal and mantle xenoliths, suggesting that the magma ascended rapidly from the upper mantle. The scoriae show significant variations in their whole-rock compositions (7.9–11.1 wt% MgO). High-MgO scoriae (MgO > ~9.5 wt%) have mostly homogeneous ⁸⁷Sr/⁸⁶Sr ratios (0.70318–0.70320), whereas low-MgO scoriae (MgO < ~9 wt%) have higher ⁸⁷Sr/⁸⁶Sr ratios (>0.70327); ratios tend to increase with decreasing MgO content. The high-MgO scoriae are aphyric, containing ~5 vol% olivine microphenocrysts with Mg# [100 × Mg/(Mg + Fe²⁺)] of up to 90. In contrast, the low-MgO scoriae have crustal xenocrysts of plagioclase, alkali feldspar, and quartz, and the mineralogic modes correlate negatively with whole-rock MgO content. On the basis of these observations, it is inferred that the high-MgO scoriae represent primary or near-primary melts, while the low-MgO scoriae underwent considerable interaction with the crust. Using thermodynamic analysis of the observed petrological features of the high-MgO scoriae, the eruption temperature of the magmas was constrained to 1,160–1,220 °C. Given that the source mantle was depleted MORB-source mantle, the primary magma was plausibly generated by ~7 % melting of a garnet-bearing spinel peridotite; taking this into consideration, and considering the constraints of multi-component thermodynamics, we estimated that the primary Sannome-gata magma was generated in the source mantle with 0.5–0.6 wt% H₂O at 1,220–1,230 °C and at ~1.8 GPa, and that the H₂O content of the primary magma was 6–7 wt%. The rear-arc Sannome-gata magma was generated by a lower degree of melting of the mantle at greater depths and lower temperatures than the frontal-arc magma from the Iwate volcano, which was also estimated to be generated by ~15 % melting of the source mantle with 0.6–0.7 wt% H₂O at ~1,250 °C and at ~1.3 GPa.     

Jurassic plume-origin ophiolites in Japan: accreted fragments of oceanic plateaus

The Mikabu and Sorachi–Yezo belts comprise Jurassic ophiolitic complexes in Japan, where abundant basaltic to picritic rocks occur as lavas and hyaloclastite blocks. In the studied northern Hamamatsu and Dodaira areas of the Mikabu belt, these rocks are divided into two geochemical types, namely depleted (D-) and enriched (E-) types. In addition, highly enriched (HE-) type has been reported from other areas in literature. The D-type picrites contain highly magnesian relic olivine phenocrysts up to Fo_93.5, and their Fo–NiO trend indicates fractional crystallization from a high-MgO primary magma. The MgO content is calculated as high as 25 wt%, indicating mantle melting at unusually high potential temperature (T _p) up to 1,650 °C. The E-type rocks represent the enrichment in Fe and LREE and the depletion in Mg, Al and HREE relative to the D-type rocks. These chemical characteristics are in good accordance with those of melts from garnet pyroxenite melting. Volcanics in the Sorachi–Yezo belts can be divided into the same types as the Mikabu belt, and the D-type picrites with magnesian olivines also show lines of evidence for production from high T _p mantle. Evidence for the high T _p mantle and geochemical similarities with high-Mg picrites and komatiites from oceanic and continental large igneous provinces (LIPs) indicate that the Mikabu and Sorachi–Yezo belts are accreted oceanic LIPs that were formed from hot large mantle plumes in the Late Jurassic Pacific Ocean. The E- and D-type rocks were formed as magmas generated by garnet pyroxenite melting at an early stage of LIP magmatism and by depleted peridotite melting at the later stage, respectively. The Mikabu belt characteristically bears abundant ultramafic cumulates, which could have been formed by crystal accumulation from a primary magma generated from Fe-rich peridotite mantle source, and the HE-type magma were produced by low degrees partial melting of garnet pyroxenite source. They should have been formed later and in lower temperatures than the E- and D-type rocks. The Mikabu and Sorachi Plateaus were formed in a low-latitude region of the Late Jurassic Pacific Ocean possibly near a subduction zone, partially experienced high P/T metamorphism during subduction, and then uplifted in association with (or without, in case of Mikabu) the supra-subduction zone ophiolite. The Mikabu and Sorachi Plateaus may be the Late Jurassic oceanic LIPs that could have been formed in brotherhood with the Shatsky Rise.     

Petrogenesis of the new lherzolitic shergottite Grove Mountains 99027: Constraints of petrography,mineral chemistry,and rare earth elements

Abstract— We report petrography, mineral chemistry, and microdistribution of rare earth elements (REE) in a new lherzolitic shergottite, Grove Mountains (GRV) 99027. The textural relationship and REE patterns of minerals suggest precipitation of cumulus olivine and chromite, followed by equilibrium crystallization of a closed system with a bulk composition of the inferred intercumulus melt. Subsolidus equilibrium temperatures of pyroxenes and olivine range from 1100 to 1210 °C, based on a two‐pyroxene thermometry and Ca partitioning between augite and olivine. Oxygen fugacity of the parent magma is 1.5–2.5 (av. 2.0 ± 0.4) log units below the quartz‐fayalite‐magnetite (QFM) buffer at 960–1360 °C, according to the olivine‐orthopyroxene‐chromite barometer. The ilmenite‐chromite barometer and thermometer show much wider ranges of oxygen fugacity (1.0–7.0 log unit below QFM) and temperature (1130–480 °C), suggesting subsolidus equilibration of the oxides at low temperatures, probably due to deep burial of GRV 99027 on Mars. The low oxygen fugacity and LREE depletion of the parent magma of GRV 99027 suggest low contamination by martian crust. Characteristics of GRV 99027 demonstrate similarity of lherzolitic shergottites, suggesting a high possibility of launch pairing or a homogeneous upper mantle of Mars if they were ejected by individual impact events. However, GRV 99027 probably experienced severe post‐shock thermal metamorphism in comparison with other lherzolitic shergottites, based on the re‐crystallization of maskelynite, the homogeneity of minerals, and the low subsolidus equilibrium temperatures between chromite and ilmenite.     

Spin rates of small moonlets embedded in planetary rings: I. Three-body calculations

We investigate the spin rates of moonlets embedded in planetary rings, subject to collisions with surrounding small particles, using three-body integrations including friction and spins. All successive impacts of the particle with the moonlet are followed, including a possible sliding phase after the initial inelastic rebounds. Two methods for treating impacts, (1) as instantaneous velocity changes and (2) using an impact force model, are applied after Salo (1995, Icarus 117, 287). Conducting a series of integrations with various initial summed spin velocity of the moonlet and the particle, we determine the equilibrium spin rate for which the averaged torque vanishes. This equilibrium spin rate corresponds to the final spin rate of the moonlet if the moonlet is much larger than the surrounding particles; it also corresponds to the mean spin rate for a ring composed of identical particles. We find that the equilibrium spin rate is enhanced by sliding orbits as compared with the spin rate determined by considering only the first impacts of the particles with the moonlet. If the random velocities of incident particles are small enough, the resulting equilibrium spin rate of the moonlet can be larger than the synchronous rotation rate, for r_p∼1, where r_p denotes the sum of radii of the colliding pair normalized by their mutual Hill radius. In this special case aggregates without internal strength may become rotationally unstable. However, the equilibrium spin rate decreases with increasing random velocity, and aggregates are always rotationally stable in the more likely case where the relative velocities are comparable to the mutual escape velocity. We also compare our results with the mean spin rates found in previous N-body simulations, and find a good agreement for optically thin rings; however the spin rates for optically thick rings are significantly larger than those predicted by our three-body calculations.     

Effects of experimental aqueous alteration on the abundances of argon‐rich noble gases in the Ningqiang carbonaceous chondrite

Abstract— Ar‐rich noble gases, the so‐called “subsolar” noble gases, are a major component of heavy primordial noble gases in unequilibrated ordinary chondrites and some classes of anhydrous carbonaceous chondrites, whereas they are almost absent in hydrous carbonaceous chondrites that suffered extensive aqueous alteration. To understand the effects of aqueous alteration on the abundance of Ar‐rich noble gases, we performed an aqueous alteration experiments on the Ningqiang type 3 carbonaceous chondrite that consists entirely of anhydrous minerals and contains Ar‐rich noble gases. Powdered samples and deionized neutral water were kept at 200 °C for 10 and 20 days, respectively. Mineralogical analyses show that, during the 10‐day alteration, serpentine and hematite formed at the expense of olivine, low‐Ca pyroxene, and sulfide. Noble gas analyses show that the 10‐day alteration of natural Ningqiang removed 79% of the primordial ³⁶Ar, 68% of the ⁸⁴Kr, and 60% of the ¹³²Xe, but only 45% of the ⁴He and 53% of the primordial ²⁰Ne. Calculated elemental ratios of the noble gases removed during the 10‐day alteration are in the range of those of Ar‐rich noble gases. These results indicate that Ar‐rich noble gases are located in materials that are very susceptible to aqueous alteration. In contrast, heavy primordial noble gases remaining in the altered samples are close to Q gas in elemental and isotope compositions. This indicates that phase Q is much more resistant to aqueous alteration than the host phases of Ar‐rich noble gases. In the 20‐day sample, the mineralogical and noble gas signatures are basically similar to those of the 10‐day sample, indicating that the loss of Ar‐rich noble gases was completed within the 10‐day alteration. Our results suggest that almost all of the Ar‐rich noble gases were lost from primitive asteroids during early, low‐temperature aqueous alteration.     

Compilation method for 1 km grid data of monthly mean air temperature for quantitative assessments of climate change impacts

A new method is proposed to compile 1 km grid data of monthly mean air temperature by dynamically downscaling general circulation model (GCM) data with a regional climate model (RCM). The downscaling method used is a technique referred to as the pseudoglobal warming method to reduce GCM bias. For the grid data, RCM data were corrected with data from an existing meteorological network. The correction model for the RCM bias was developed by stepwise multiple regression analysis using the difference in the monthly mean air temperatures between the observation and RCM output as a dependent variable and the geographical factors as independent variables. Our method corrected the RCM bias from 1.69°C to 0.58°C for the month of August in the 1990s (1990–1999).     

The impact of the Tokyo Metropolitan Emissions Trading Scheme on reducing greenhouse gas emissions: findings from a facility-based study

This paper provides a detailed analysis of the Tokyo Metropolitan Emissions Trading Scheme (Tokyo ETS), Japan’s first emissions trading scheme with mandatory cap initiated by the government of Tokyo. Unlike trading schemes in other countries, the Tokyo ETS covers indirect emissions from the commercial sector. It is well known that a variety of market barriers impede full realization of energy efficiency opportunities, especially in the commercial sector. Experiences with the Tokyo ETS should therefore provide important lessons for the design of climate change mitigation policies, especially when targeting the commercial sector. The emissions from covered entities have been drastically reduced from those at the scheme’s outset, with an average 14% reduction as of the end of the first commitment period of five years (2010–2014) compared with 2009 levels. This paper shows that the Tokyo ETS alone did not cause these reductions; there were other drivers. Among them, the energy savings triggered by the Great East Japan Earthquake in 2011 were crucial. The contribution of credit trading, in contrast, was limited since most of the covered entities reduced emissions by themselves. Through an investigation of official reports, an assessment of the emissions data from the covered entities compared to those of uncovered entities and in-depth interviews with firms covered by the scheme, this paper confirms that the main drivers of emissions reductions by covered entities were separate from the ETS. In fact, the advisory aspect of the scheme seems to be much more important in encouraging energy-saving actions.

Key policy insights

• Most of the observed emission reductions were not caused by the Tokyo ETS alone.

• An advisory instrument was crucial to the effectiveness of the Tokyo ETS.

• The experience of the Tokyo ETS suggests that making full use of the advantages of emissions trading is difficult in the case of the commercial sector.

• Price signals have not provided a stimulus to climate change mitigation actions, which implies that establishing a cap to yield effective carbon prices poses a challenge.

     

A pair of Langmuir cells in two laboratory tanks (II) on generation mechanism

Vertical and cross-wind profiles of mean currents were measured systematically in vertical cross-sections of two wind-wave tanks with aspect ratios of order one to study the secondary flow in the tanks. A pair of Langmuir cells turned out to be driven by a close combination of the pressure gradient along the tank and the side-wall effects. That is, part of the adverse pressure gradient produced a parabolic cross-wind profile with the smallest downwind current at the centerline and the largest current along the two sidewalls. As a result, upwelling occurred in the center zone where the return flow was strongest, probably because of the entrainment action of the wind-driven current. In order to compensate for this upwelling, downwelling occurred along the two side-walls from the flow continuity. The resulting vertical circulation formed a pair of Langmuir circulations across the span and served to maintain the parabolic profile formed by the pressure gradient. A positive feedback mechanism is thus found between the primary and secondary circulations through upwelling of the return flow in the center zone. Vertical shears of the span-averaged downwind current measured in two tanks were found to be systematically different from each other. This difference seems to depend on the magnitude of the advective Reynolds stresses in the two tanks.     

Strength and deformation behavior of the Shimanto accretionary complex across the Nobeoka thrust

A rapid reduction in sediment porosity from 60 to 70 % at seafloor to less than 10 % at several kilometers depth can play an important role in deformation and seismicity in the shallow portion of subduction zones. We conducted deformation experiments on rocks from an ancient accretionary complex, the Shimanto Belt, across the Nobeoka Thrust to understand the deformation behaviors of rocks along plate boundary faults at seismogenic depth. Our experimental results for phyllites in the hanging wall and shale‐tuff mélanges in the footwall of the Nobeoka Thrust indicate that the Shimanto Belt rocks fail brittlely accompanied by a stress drop at effective pressures < 80 MPa, whereas they exhibit strain hardening at higher effective pressures. The transition from brittle to ductile behavior in the shale–tuff mélanges lies on the same trend in effective stress–porosity space as that for clay‐rich and tuffaceous sediments subducting into the modern Nankai subduction zone. Both the absolute yield strength and the effective pressure at the brittle–ductile transition for the phyllosilicate‐rich materials are much lower than for sandstones. These results suggest that as the clay‐rich or tuffaceous sediments subduct and their porosities are reduced, their deformation behavior gradually transitions from ductile to brittle and their yield strength increases. Our results also suggest that samples of the ancient Shimanto accretionary prism can serve as an analog for underthrust rocks at seismogenic depth in the modern Nankai Trough.