Comparison of Observations at ACE and <Emphasis Type="Italic">Ulysses</Emphasis> with Enlil Model Results: Stream Interaction Regions During Carrington Rotations 2016 – 2018

During the latitudinal alignment in 2004, ACE and Ulysses encountered two stream interaction regions (SIRs) each Carrington rotation from 2016 to 2018, at 1 and 5.4 AU, respectively. More SIR-driven shocks were observed at 5.4 AU than at 1 AU. Three small SIRs at 1 AU merged to form a strong SIR at 5.4 AU. We compare the Enlil model results with spacecraft observations from four aspects: i) the accuracy of the latest versions of models (WSA v2.2 and Enlil v2.7) vs. old versions (WSA v1.6 and Enlil v2.6), ii) the sensitivity to different solar magnetograms (MWO vs. NSO), iii) the sensitivity to different coronal models (WSA vs. MAS), iv) the predictive capability at 1 AU vs. 5.4 AU. We find the models can capture field sector boundaries with some time offset. Although the new versions have improved the SIR timing prediction, the time offset can be up to two days at 1 AU and four days at 5.4 AU. The models cannot capture some small-scale structures, including shocks and small SIRs at 1 AU. For SIRs, the temperature and total pressure are often underestimated, while the density compression is overestimated. For slow wind, the density is usually overestimated, while the temperature, magnetic field, and total pressure are often underestimated. The new versions have improved the prediction of the speed and density, but they need more robust scaling factors for magnetic field. The Enlil model results are very sensitive to different solar magnetograms and coronal models. It is hard to determine which magnetogram-coronal model combination is superior to others. Higher-resolution solar and coronal observations, a mission closer to the Sun, together with simulations of greater resolution and added physics, are ways to make progress for the solar wind modeling.     

Long term dynamical evolution and classification of classical TNOs

Classical trans-Neptunian objects (TNOs) are believed to represent the most dynamically pristine population in the trans-Neptunian belt (TNB) offering unprecedented clues about the formation of our Solar System. The long term dynamical evolution of classical TNOs was investigated using extensive simulations. We followed the evolution of more than 17000 particles with a wide range of initial conditions taking into account the perturbations from the four giant planets for 4 Gyr. The evolution of objects in the classical region is dependent on both their inclination and semimajor axes, with the inner (a<45 AU) and outer regions (a>45 AU) evolving differently. The reason is the influence of overlapping secular resonances with Uranus and Neptune (40–42 AU) and the 5:3 (a∼ ∼42.3 AU), 7:4 (a∼ ∼43.7 AU), 9:5 (a∼ ∼44.5 AU) and 11:6 (a∼ ∼ 45.0 AU) mean motion resonances strongly sculpting the inner region, while in the outer region only the 2:1 mean motion resonance (a∼ ∼47.7 AU) causes important perturbations. In particular, we found: (a) A substantial erosion of low-i bodies (i<10°) in the inner region caused by the secular resonances, except those objects that remained protected inside mean motion resonances which survived for billion of years; (b) An optimal stable region located at 45 AU<a<47 AU, q>40 AU and i>5° free of major perturbations; (c) Better defined boundaries for the classical region: 42–47.5 AU (q>38 AU) for cold classical TNOs and 40–47.5 AU (q>35 AU) for hot ones, with i=4.5° as the best threshold to distinguish between both populations; (d) The high inclination TNOs seen in the 40–42 AU region reflect their initial conditions. Therefore they should be classified as hot classical TNOs. Lastly, we report a good match between our results and observations, indicating that the former can provide explanations and predictions for the orbital structure in the classical region.     

The Radial Variation of Interplanetary Shocks in the Inner Heliosphere: Observations by <Emphasis Type="Italic">Helios</Emphasis>, MESSENGER,and STEREO

The two major sources of collisionless shocks in the solar wind are interplanetary coronal mass ejections (ICMEs) and stream interaction regions (SIRs). Previous studies show that some SIR-associated shocks form between Venus and Earth while most form beyond 1 AU. Here we examine the high-resolution magnetometer records from Helios 1 and 2 obtained between 0.28 and 1 AU and from MESSENGER obtained between 0.3 and 0.7 AU to further refine our understanding as to where, and in what context, shocks are formed in the inner solar system. From Helios data (Helios 1 from 1974 to 1981 and Helios 2 from 1976 to 1980), we find there were only a few shocks observed inside the orbit of Venus with the closest shock to the Sun at 0.29 AU. We find that there is a strong correlation between shock occurrence and solar activity as measured by the sunspot number. Most of the shocks inside of the orbit of Venus appear to be associated with ICMEs. Even the ICME-associated shocks are quite weak inside the orbit of Venus. By comparing MESSENGER and STEREO results, from 2007 to 2009, we find that in the deep solar minimum, SIR-driven shocks began to form at about 0.4 AU and increased in number with heliocentric distance.     

Solar wind heating beyond 1 AU

The effect of an interplanetary atomic hydrogen gas on solar wind proton, electron and α-particle temperatures beyond 1 AU is considered. It is shown that the proton temperature (and probably also the α-particle temperature) reaches a minimum between 2 AU and 4 AU, depending on values chosen for solar wind and interstellar gas parameters. Heating of the electron gas depends primarily on the thermal coupling of the protons and electrons. For strong coupling (whenT _p ≳T _e ), the electron temperature reaches a minimum between 4 AU and 8 AU, but for weak coupling (Coulomb collisions only), the electron temperature continues to decrease throughout the inner solar system. A spacecraft travelling to Jupiter should be able to observe the heating effect of the solar wind-interplanetary hydrogen interaction, and from such observations it may be possible of infer some properties of the interstellar neutral gas. Currently a National Research Council Resident Research Associate.     

Stream Interactions and Interplanetary Coronal Mass Ejections at 0.72 AU

We present comprehensive surveys of 203 stream interaction regions (SIRs) and 124 interplanetary CMEs (ICMEs) during 1979 – 1988 using Pioneer Venus Orbiter (PVO) in situ solar-wind observations at 0.72 AU and examine the solar-cycle variations of the occurrence rate, shock association rate, duration, width, maximum total perpendicular pressure (P _t), maximum dynamic pressure, maximum magnetic field intensity, and maximum velocity change of these two large-scale solar-wind structures. The medians, averages, and histogram distributions of these parameters are also reported. Furthermore, we sort ICMEs into three groups based on the temporal profiles of P _t, and we investigate the variations of the fractional occurrence rate of three groups of ICMEs with solar activity. We find that the fractional occurrence rate of magnetic-cloud-like ICMEs declined with solar activity, consistent with our former 1-AU results. This study at 0.72 AU provides a point of comparison in the inner heliosphere for examining the radial evolution of SIRs and ICMEs. The width of SIRs and ICMEs increases by 0.04 and 0.1 AU, respectively, and the maximum P _t decreases to about 1/3 from Venus to Earth orbit. In addition, our work establishes the statistical properties of the solar-wind conditions at 0.72 AU that control the solar-wind interaction with Venus and its atmosphere loss by related processes. Electronic Supplementary Material The online version of this article () contains supplementary material, which is available to authorized users.     

The effects of nebula surface density profile and giant‐planet eccentricities on planetary accretion in the inner solar system

Abstract— We describe results of 32 N‐body planetary accretion simulations that investigate the dependence of terrestrial‐planet formation on nebula surface density profile σ and evolution of the eccentricities of Jupiter and Saturn e_j,s. Two surface density profiles are examined: a decaying profile with σ ∝ 1/a, where a is orbital semi‐major axis, and a peaked profile in which σ increases for a < 2 AU and decreases for a > 2 AU. The peaked profiles are generated by models of coagulation in an initially hot nebula. Models with initial e_j,s = 0.05 (the current value) and 0.1 are considered. Simulations using the decaying profile with e_j,s = 0.1 produce systems most like the observed planets in terms of mass‐weighted mean a and the absence of a planet in the asteroid belt. Simulations with doubled σ produce planets roughly twice as massive as the nominal case. Most initial embryos are removed in each simulation via ejection from the solar system or collision with the Sun. The asteroid belt is almost entirely cleared on a timescale of 10–100 Ma that depends sensitively on e_j,s. Most initial mass with a < 2 AU survives, with the degree of mass loss increasing with a. Mass loss from the terrestrial region occurs on a timescale that is long compared to the mass loss time for the asteroid belt. Substantial radial mixing of material occurs in all simulations, but is greater in simulations with initital e_j,s = 0.05. The degree of mixing is equivalent to a feeding zone of half width 1.5 and 0.9 AU for an Earth mass planet at 1 AU for the cases e_j,s = 0.05 and 0.1, respectively. In simulations with e_j,s = 0.05, roughly one‐third and 5–10% of the mass contained in final terrestrial planets originated in the region a > 2.5 AU for the decaying and peaked profiles, respectively. In the case e_j,s = 0.1, the median mass accreted from a > 2.5 AU is zero for both profiles.     

Examining Periodic Solar-Wind Density Structures Observed in the SECCHI <Emphasis Type="Italic">Heliospheric Imagers</Emphasis>

We present an analysis of small-scale, periodic, solar-wind density enhancements (length scales as small as ≈ 1000 Mm) observed in images from the Heliospheric Imager (HI) aboard STEREO-A. We discuss their possible relationship to periodic fluctuations of the proton density that have been identified at 1 AU using in-situ plasma measurements. Specifically, Viall, Kepko, and Spence (J. Geophys. Res. 113, A07101, 2008) examined 11 years of in-situ solar-wind density measurements at 1 AU and demonstrated that not only turbulent structures, but also nonturbulent, periodic density structures exist in the solar wind with scale sizes of hundreds to one thousand Mm. In a subsequent paper, Viall, Spence, and Kasper (Geophys. Res. Lett. 36, L23102, 2009) analyzed the α-to-proton solar-wind abundance ratio measured during one such event of periodic density structures, demonstrating that the plasma behavior was highly suggestive that either temporally or spatially varying coronal source plasma created those density structures. Large periodic density structures observed at 1 AU, which were generated in the corona, can be observable in coronal and heliospheric white-light images if they possess sufficiently high density contrast. Indeed, we identify such periodic density structures as they enter the HI field of view and follow them as they advect with the solar wind through the images. The smaller, periodic density structures that we identify in the images are comparable in size to the larger structures analyzed in-situ at 1 AU, yielding further evidence that periodic density enhancements are a consequence of coronal activity as the solar wind is formed.     

Solar Wind Electron Parameters from Ulysses/URAP Quasi-Thermal Noise Measurements at Solar Maximum

We present the solar wind plasma parameters obtained from the Ulysses spacecraft during its second pole-to-pole fast latitude scan near the 2001 solar maximum. We study the solar wind properties from the electron density and core temperature measurements made by the radio receiver on Ulysses using the method of quasi-thermal noise spectroscopy. We analyze these parameters as functions of heliographic latitude and distance. We present their histograms normalized to 1 AU and find a bimodal distribution for the electron core temperature. The cooler population can be associated with the fast wind flow emanating from coronal holes present at various latitudes. We discuss a slight north/south asymmetry found for the electron density. Finally, we compare the present results to those obtained during the 1996 solar minimum and 1991 solar maximum.     

Stochastic wave growth, power balance, and beam evolution in type III solar radio sources

Energy-balance arguments are combined with the stochastic-growth theory of type III radio sources to determine the properties of the source in average dynamical equilibrium with the beam, and the beam's long-term evolution. Purely linear stochastic-growth theory has previously emphasized that the beam evolves to a state close to marginal stability. Small mean residual deviations from marginal stability are present at dynamical equilibrium and these lead to residual energy flows that feed the waves observed in situ and by remote receivers; consequently the beam energy is depleted. Here, dynamical equilibrium beam and wave levels are estimated for the first time and it is found that the main sink of beam-driven Langmuir waves is either via electrostatic decay into product Langmuir and ion-sound waves or via scattering by short-wavelength density fluctuations, depending on the conditions. Improved estimates of energy branching ratios imply that, at 1 AU from the Sun, typically 20% of the beam energy is converted to Langmuir waves that are scattered off low-frequency density fluctuations and then dissipated, with almost all the remaining waves undergoing electrostatic decay, although as little as one-third of the Langmuir waves may decay in atypical circumstances. Of order 10^–3 of the beam energy is converted into sound waves, which are mostly dissipated, and of order 10^–5 is converted into potentially observable electromagnetic waves. The mean lifetime of the Langmuir waves at 1 AU is 1–40 s, while that of the beam is of order 1000 s. The beam density decreases relative to that of the background as the beam propagates. For most parameters, analysis of energy losses from the beam to the waves shows that the beam velocity decreases at roughly the same rate as the thermal velocity of the background plasma. It is argued from these considerations, and from in situ observations at 1 AU, that these trends imply that only the densest and fastest type III beams will be able to penetrate much past 1 AU from the Sun. This implies a low-frequency cutoff to type III emission at roughly 10 kHz, in good agreement with recent Ulysses remote observations, showing their consistency with in situ measurements.     

The Bastille day Shocks and Merged Interaction Region

Near 1 AU the solar wind structure associated with the solar flare of 14 July 2000 (Bastille Day) consisted of a large high-speed stream of 15 July and five nearby small streams during a 10-day period. At the leading edge of the large high-speed stream, in less than 6 hours, the flow speed increased from 600 km s⁻¹ to 1100 km s⁻¹, the magnetic field intensity increased from 10 nT to 60 nT, and an interaction region was identified. The interaction region was bounded between the pair of a forward shock F and a reverse shock R. Additional forward shocks were also identified at the leading edge of each of the five smaller streams. This paper presents a magnetohydrodynamics (MHD) simulation using ACE plasma and magnetic field data near 1 AU as input to study the radial evolution of the Bastille Day solar wind event. The two shocks, F and R, propagated in opposite directions away from each other in the solar wind frame and interacted with neighboring shocks and streams; the spatial and temporal extent of the interaction region continued to increase with the heliocentric distance. The solar wind was restructured from a series of streams at 1 AU to a huge merged interaction region (MIR) extending over a period of 12 days at 5.5 AU. Throughout the interior of the MIR bounded by the shock pair F and R the magnetic field intensity was a few times stronger than that outside the MIR. The simulation shows how merging of shocks, collision of shocks, and formation of new shocks contributed to the evolution process.     

An infrared imaging search for low‐mass companions to members of the young nearby β Pic and Tucana/Horologium associations

We present deep high dynamic range infrared images of young nearby stars in the Tucana/Horologium and β Pic associations, all ∼10 to 35 Myrs young and at ∼10 to 60 pc distance. Such young nearby stars are well‐suited for direct imaging searches for brown dwarf and even planetary companions, because young sub‐stellar objects are still self‐luminous due to contraction and accretion. We performed our observations at the ESO 3.5m NTT with the normal infrared imaging detector SofI and the MPE speckle camera Sharp‐I. Three arc sec north of GSC 8047‐0232 in Horologium a promising brown dwarf companion candidate is detected, which needs to be confirmed by proper motion and/or spectroscopy. Several other faint companion candidates are already rejected by second epoch imaging. Among 21 stars observed in Tucana/Horologium, there are not more than one to five brown dwarf companions outside of 75 AU (1.5″ at 50 pc); most certainly only ≤5% of the Tuc/HorA stars have brown dwarf companions (13 to 78 Jupiter masses) outside of 75 AU. For the first time, we can report an upper limit for the frequency of massive planets (∼10 M_jup) at wide separations (∼100 AU) using a meaningfull and homogeneous sample: Of 11 stars observed sufficiently deep in β Pic (12 Myrs), not more than one has a massive planet outside of ∼100 AU, i.e. massive planets at large separations are rare (≤9%). (© 2003 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Formation rate of triple stars relative to double stars

A series of our and McAlisteret al.'s speckle observations of spectroscopic binaries show that certain fractions of them are triple systems. A typical configuration of those systems is that of two stars separated by about 0.1 AU with a third component in an orbit at a distance about 10 AU from the central stars. These results suggest that the third star was formed from the outer circumstellar envelope after forming the central double.     

Near-parabolic cometary flux in the outer solar system

The flux of near-parabolic comets in the outer-planetary region is estimated on the presumption that the major planets and the galactic tide control the dynamics of comets. It is found that the flux of the Oort cloud comets (semi-major axis > 20000 AU) is similar to the case of a strong comet shower derived on the presumption that the galactic tidal force were not operative. On the other hand, the flux of comets with semi-major axes <- 20000 AU is found to be an increasing function of q (perihelion distance) until q reaches 20 AU, while for a 45000 AU it is a rapidly increasing function for q 12 AU. In other words, for comets of the inner extension of the Oort cloud the planetary perturbation acts as a strong barrier for them to penetrate into the inner planetary region.     

Some requirements of a colliding comet source of gamma ray bursts

Colliding comets in the Solar System may be an important source of gamma ray bursts. The spherical gamma ray comet cloud required by the results of the Venera Satellites (Mazets and Golenetskii, 1987) and the BATSE detector on the Compton Satellite (Meeganet al., 1992a, b) is neither the Oort Cloud nor the Kuiper Belt. To satisfy observations ofN(>P _max) vsP _max for the maximum gamma ray fluxes,P _max > 10^–5 erg cm^–2 s^–1 (about 30 bursts yr^–1), the comet density,n, should increase asn a ¹ from about 40 to 100 AU wherea is the comet heliocentric distance. The turnover above 100 AU requiresn a ^–1/2 to 200 AU to fit the Venera results andn a ^1/4 to 400 AU to fit the BATSE data. Then the masses of comets in the 3 regions are from: 40–100 AU, about 9 earth masses,m _E; 100–200 AU about 25m _E; and 100–400 AU, about 900m _E. The flux of 10^–5 erg cm^–2 s^–1 corresponds to a luminosity at 100 AU of 3 × 10²⁶ erg s^–1. Two colliding spherical comets at a distance of 100 AU, each with nucleus of radiusR of 5 km, density of 0.5 g cm^–3 and Keplerian velocity 3 km s^–1 have a combined kinetic energy of 3 × 10²⁸ erg, a factor of about 100 greater than required by the burst maximum fluxes that last for one second. Betatron acceleration in the compressed magnetic fields between the colliding comets could accelerate electrons to energies sufficient to produce the observed high energy gamma rays. Many of the additional observed features of gamma ray bursts can be explained by the solar comet collision source.     

Electromagnetic lower hybrid instability driven by solar wind heat flux

Under a fully electromagnetic treatment, the threshold for excitation of the lower hybrid instability driven by solar wind electron heat flux is found to be much higher than that predicted by electrostatic approximation. For average solar wind conditions at 1 AU, the fully electromagnetic lower hybrid instability is excited when the core electron drift speed is about 8V _A, whereV _A is the Alfvén speed. The region between the Sun and 1 AU is expected to be more favourable than 1 AU for this instability.     

The multiple modes of interaction of the solar wind with a comet as it approaches the Sun

A quasi-steady 1-D hydrodynamic model, with mass addition, has been used to study the various modes of interaction of the solar wind with a medium-bright, H₂O-dominated comet (such as P/Halley) approaching the Sun.At large heliocentric distances (d 5 AU) the solar wind penetrates unimpeded on to the surface. As the comet moves further in, mass loading of the solar wind by heavy ions from the fledgling cometary atmosphere causes it to slow down, thereby causing a significant enhancement of the interplanetary field. Still further in at d 3.14 AU, as the mass loading reaches a critical value, a collision-less standing shock is formed in the solar wind upstream of the nucleus. As d decreases further, the distance of this shock from the nucleus increases. The cometary atmosphere becomes dense enough to stand off the solar wind ahead of the nuclear surface and form a well defined tangential discontinuity surface (or ionopause) only when d reaches the value 2.65 AU. When d 2.65 AU an inner shock could, in principle, also form within the cometary ionosphere, although its existence would depend on the detailed thermodynamics of the cometary ionosphere. Resolution of this question is beyond the scope of the present analysis.The conclusions of the present study would be qualitatively valid for other comets having sizes, surface optical properties and chemical compositions, different from those adopted here. The helio-centric distances at which the various transitions take place from the one mode of solar wind interaction to another, would, of course, be different, with all these distances being smaller for less active comets.     

Interplanetary baseline observations of type III solar radio bursts

Simultaneous observations of type III radio bursts from spacecraft separated by 0.43 AU have been made using the solar orbiters HELIOS-A and HELIOS-B. The burst beginning at 19:22 UT on March 28, 1976 has been located from the intersection of the source directions measured at each spacecraft, and from burst arrival time differences. The source positions range from 0.03 AU from the Sun at 3000 kHz to 0.08 AU at 585 kHz. The electron density along the burst trajectory, and the exciter velocity (=0.13c) were determined directly, without the need to assume a density model as has been done with single-spacecraft observations. The separation of HELIOS-A and -B has also provided the first measurements of burst directivity at low frequencies. For the March 28 burst the intensity observed from near the source longitude (HELIOS-B) was 3–10 dB greater than that from 60° west of the source (HELIOS-A).     

Solar radio burst and in situ determination of interplanetary electron density

We review and discuss a few interplanetary electron density scales which have been derived from the analysis of interplanetary solar radio bursts, and we compare them to a model derived from 1974–1980 Helios 1 and 2 in situ density observations made in the 0.3–1.0 AU range. The Helios densities were normalized to 1976 with the aid of IMP and ISEE data at 1 AU, and were then sorted into 0.1 AU bins and logarithmically averaged within each bin. The best fit to these 1976-normalized, bin averages is N(R _AU) = 6.1R ^-2.10 cm^-3. This model is in rather good agreement with the solar burst determination if the radiation is assumed to be on the second harmonic of the plasma frequency. This analysis also suggests that the radio emissions tend to be produced in regions denser than the average where the density gradient decreases faster with distance than the observed R ^-2.10.NAS/NRC Postdoctoral Research Associate on leave from Laboratory Associated with CNRS No. 264, Paris Observatory, France.     

A model of the galactic tidal interaction with the oort comet cloud

We consider a model of the in situ Oort cloud which is isotropic with a random distrihution of perihelia directions and angular momenta. The energy distribution adopted has a continuous range of values appropriate for long-period (>200 yr) comets. Only the tidal torque of the Galaxy is included as a perturbation of comet orbits and it is approximated to be that due to a quasi-steady state distribution of matter with disk-like symmetry. The time evolution of all orbital elements can be analytically obtained for this case. In particular, the change in the perihelion distance per orbit and its dependence on other orbital elements is readily found. We further make the assumption that a comet whose perihelion distance was beyond 15 AU during its last passage through the Solar System would have orbit parameters that are essentially unchanged by planetary perturbations. Conversely, if the prior passage was inside 15 AU we assume that planetary perturbations would have removed the comet from the in situ energy distribution accessible by the galactic tide. Comets which had their perihelia changed from beyond 15 AU to within 5 AU in a single orbit are taken to be observable. We are able to track the evolution of 10⁶ comets as they are made observable by the galactic tidal touque. Detailed results are obtained for the predicted distribution of new (0 < 1/ < 10^–4 AU^–1) comets. Further, correlations between orbital elements can be studied. We present predictions of observed distributions and compare them with the random in situ results as well as with the actual observed distributions of class I comets. The predictions are in reasonable agreement with actual observations and, in many cases, are significantly different from random when perihelia directions are separated into galactic northern and southern hemispheres. However the well-known asymmetry in the north-south populations of perihelia remains to be explained. Such an asymmetry is consistent with the dominance of tidal torques today if a major stochastic event produced it in the past since tidal torques are unable to cause the migration of perihelia across the latitude barriers ±26°.6 in the disk model.     

On the dynamic connection of comets with Uranus

We consider the connection with Uranus for: (1) 945 near-parabolic comets (the period P > 200 years, the perihelion distance q > 0.1 AU), (2) 1277 Kreutz comets (P > 200 years, q < 0.01 AU), and (3) 414 short-period comets (P < 200 years). It turns out that none of near-parabolic comets passed through Uranus’s activity sphere, none of the Kreutz comets approach Uranus closer than 11 AU, and only two short-period comets, C/2006 U7 and C/2006 F2, could have a close approach to Uranus during 5000 years.