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.     

Chemistry of atmospheres formed during accretion of the Earth and other terrestrial planets

We used chemical equilibrium and chemical kinetic calculations to model chemistry of the volatiles released by heating different types of carbonaceous, ordinary and enstatite chondritic material as a function of temperature and pressure. Our results predict the composition of atmospheres formed by outgassing during accretion of the Earth and other terrestrial planets. Outgassing of CI and CM carbonaceous chondritic material produces H₂O-rich (steam) atmospheres in agreement with the results of impact experiments. However, outgassing of other types of chondritic material produces atmospheres dominated by other gases. Outgassing of ordinary (H, L, LL) and high iron enstatite (EH) chondritic material yields H₂-rich atmospheres with CO and H₂O being the second and third most abundant gases. Outgassing of low iron enstatite (EL) chondritic material gives a CO-rich atmosphere with H₂, CO₂, and H₂O being the next most abundant gases. Outgassing of CV carbonaceous chondritic material gives a CO₂-rich atmosphere with H₂O being the second most abundant gas. Our results predict that the atmospheres formed during accretion of the Earth and Mars were probably H₂-rich unless the accreted material was dominantly CI and CM carbonaceous chondritic material. We also predict significant amounts of S, P, Cl, F, Na, and K in accretionary atmospheres at high temperatures (1500-2500 K). Finally, our results may be useful for interpreting spectroscopic observations of accreting extrasolar terrestrial planets.     

Lead angles and emitting electron energies of Io-controlled decameter radio arcs

The Io-controlled radio arcs are emissions in the decametric radio range which appear arc shaped in the time-frequency plane. Their occurrence is controlled by Io's position, so it has been for long inferred that they are powered by the Io-Jupiter electrodynamic interaction. Their frequency ranges correspond to the electron cyclotron frequencies along the Io Flux tube, so they are expected to be generated by cyclotron maser instability (CMI). The arc shape was proposed to be a consequence of the strong anisotropy of the decametric radio emissions beaming, combined with the topology of the magnetic field in the source and the observation geometry. Recent papers succeeded at reproducing the morphologies of a few typical radio arcs by modeling in three dimensions the observation geometry, using the best available magnetic field model and a beaming angle variation consistent with a loss-cone driven CMI. In the continuation of these studies, we present here the systematic modeling of a larger number of observations of the radio arcs emitted in Jupiter's southern hemisphere (including multiple arcs or arcs exhibiting abrupt changes of shape), which permits to obtain a statistical determination of the emitting field line localization (lead angle) relative to the instantaneous Io field line, and of the emitting particle velocities or energies. Variations of these parameters relative to Io's longitude are also measured and compared to the location of the UV footprints of the Io-Jupiter interaction. It is shown that the data are better organized in a reference frame attached to the UV spot resulting from the main Alfvén wing resulting from the Io-Jupiter interaction. It is proposed that the radio arcs are related to the first reflected Alfvén wing rather than to the main one.     

Abundances of Jupiter's trace hydrocarbons from Voyager and Cassini

The flybys of Jupiter by the Voyager spacecraft in 1979, and over two decades later by Cassini in 2000, have provided us with unique datasets from two different epochs, allowing the investigation of seasonal change in the atmosphere. In this paper we model zonal averages of thermal infrared spectra from the two instruments, Voyager 1 IRIS and Cassini CIRS, to retrieve the vertical and meridional profiles of temperature, and the abundances of the two minor hydrocarbons, acetylene (C₂H₂) and ethane (C₂H₆). The spatial variation of these gases is controlled by both chemistry and dynamics, and therefore their observed distribution gives us an insight into both processes. We find that the two gases paint quite different pictures of seasonal change. Whilst the 2-D cross-section of C₂H₆ abundance is slightly increased and more symmetric in 2000 (northern summer solstice) compared to 1979 (northern fall equinox), the major trend of equator to pole increase remains. For C₂H₂ on the other hand, the Voyager epoch exhibits almost no latitudinal variation, whilst the Cassini era shows a marked decrease polewards in both hemispheres. At the present time, these experimental findings are in advance of interpretation, as there are no published models of 2-D Jovian seasonal chemical variation available for comparison.     

Saturn eddy momentum fluxes and convection: First estimates from Cassini images

We apply an automated cloud feature tracking algorithm to estimate eddy momentum fluxes in Saturn's southern hemisphere from Cassini Imaging Science Subsystem near-infrared continuum image sequences. Voyager Saturn manually tracked images had suggested no conversion of eddy to mean flow kinetic energy, but this was based on a small sample of <1000 wind vectors. The automated procedure we use for the Cassini data produces an order of magnitude more usable wind vectors with relatively unbiased sampling. Automated tracking is successful in and around the westward jet latitudes on Saturn but not in the vicinity of most eastward jets, where the linearity and non-discrete nature of cloud features produces ambiguous results. For the regions we are able to track, we find peak eddy fluxes and a clear positive correlation between eddy momentum fluxes and meridional shear of the mean zonal wind, implying that eddies supply momentum to eastward jets and remove momentum from westward jets at a rate . The behavior we observe is similar to that seen on Jupiter, though with smaller eddy-mean kinetic energy conversion rates per unit mass of atmosphere (). We also use the appearance and rapid evolution of small bright features at continuum wavelengths, in combination with evidence from weak methane band images where possible, to diagnose the occurrence of moist convective storms on Saturn. Areal expansion rates imply updraft speeds of over the convective anvil cloud area. As on Jupiter, convection preferentially occurs in cyclonic shear regions on Saturn, but unlike Jupiter, convection is also observed in eastward jet regions. With one possible exception, the large eddy fluxes seen in the cyclonic shear latitudes do not seem to be associated with convective events.     

Cassini图像中受散射光影响的7颗土星近环内卫星的天体测量

卡西尼号探测器携带了光学成像科学子系统(ISS: Imaging Science Subsystem), 在2004--2017年间拍摄了一些土星内卫星的图像. 部分图像中土星内卫星非常靠近土星环, 观测对象受土星环散射光影响导致测量精度差, 甚至无法测量. 由此提出一种适用于这类图像的背景消除算法, 可对受土星光环散射光影响的土星内卫星进行测量. 归算了7颗近环内卫星(土卫十(Janus)、土卫十一(Epimetheus)、土卫十五(Atlas)、土卫十六(Prometheus)、土卫十七(Pandora)、土卫三十二(Methone)和土卫四十九(Anthe))的70张ISS图像, 并与不消除散射光的方法进行了比较. 结果显示, 该方法至少可以提高43%的精度. 和喷气推进实验室的土星卫星历表\lk SAT415相比, 测量得到的赤经和赤纬方向的残差均值分别为0.72km和2.26km, 标准差分别为10.99km和11.36km.     

Tidal pattern instabilities on multi-moon planets

The equilibrium tide-generating forces in the lunar orbital plane of a planet of radius R are calculated for the case of N moons of mass Mi orbiting the planet at instantaneous polar coordinates (Di, αi). For the case of a single moon, there are only two high tides. For the case of two moons, it is found that there can exist a critical lunar orbital distance at which two high tides become unstable with respect to formation of three high tides. Bifurcation diagrams are presented which depict how the angular positions of the high and low tides on the planet vary with the lunar distances and lunar separation angle. Tidal stability diagrams, which illustrate the stability regions for various tidal patterns as a function of lunar distances and lunar separation angle, are presented for various values of D₂/D₁ and M₂/M₁. Generally speaking, the aforementioned tidal instability, and hence the propensity for formation of three high tides on a two-moon planet, exists over a significant range of lunar distances and separation angles provided that M₂/M₁∼³(D₂/D₁). For the case of N>2 moons, the tidal stability diagram becomes more complex, revealing a diversity of potential tidal patterns.     

大行星运动对太阳黑子活动的可能影响

本文计算了由太阳系大行星轨道运动引起的日心相对于太阳系质心的轨道运动角动量变化率j_⊙,在理论上对j_⊙作展开,表明它存在多项短周期变化,与太阳黑子资料的分析结果相比较,两者结果是符合的,它们具有一致的谱结构。因此,行星的轨道运动对太阳黑子活动存在动力作用的可能性又进一步得到了验证。