A ballistics analysis of the Deep Impact ejecta plume: Determining Comet Tempel 1's gravity, mass, and density

In July of 2005, the Deep Impact mission collided a 366 kg impactor with the nucleus of Comet 9P/Tempel 1, at a closing speed of 10.2 km s⁻¹. In this work, we develop a first-order, three-dimensional, forward model of the ejecta plume behavior resulting from this cratering event, and then adjust the model parameters to match the flyby-spacecraft observations of the actual ejecta plume, image by image. This modeling exercise indicates Deep Impact to have been a reasonably “well-behaved” oblique impact, in which the impactor-spacecraft apparently struck a small, westward-facing slope of roughly 1/3-1/2 the size of the final crater produced (determined from initial ejecta plume geometry), and possessing an effective strength of not more than . The resulting ejecta plume followed well-established scaling relationships for cratering in a medium-to-high porosity target, consistent with a transient crater of not more than 85-140 m diameter, formed in not more than 250-550 s, for the case of (gravity-dominated cratering); and not less than 22-26 m diameter, formed in not less than 1-3 s, for the case of (strength-dominated cratering). At , an upper limit to the total ejected mass of 1.8×¹⁰⁷ kg (1.5-2.2×¹⁰⁷ kg) is consistent with measurements made via long-range remote sensing, after taking into account that 90% of this mass would have stayed close to the surface and then landed within 45 min of the impact. However, at , a lower limit to the total ejected mass of 2.3×¹⁰⁵ kg (1.5-2.9×¹⁰⁵ kg) is also consistent with these measurements. The expansion rate of the ejecta plume imaged during the look-back phase of observations leads to an estimate of the comet's mean surface gravity of (0.17-0.90 mm s⁻²), which corresponds to a comet mass of mt=4.5×¹⁰¹³ kg (2.3-12.0×¹⁰¹³ kg) and a bulk density of (200-1000 kg m⁻³), where the large high-end error is due to uncertainties in the magnitude of coma gas pressure effects on the ejecta particles in flight.     

Near-Earth Interplanetary Coronal Mass Ejections During Solar Cycle 23 (1996 – 2009): Catalog and Summary of Properties

In a previous study (Cane and Richardson, J. Geophys. Res. 108(A4), SSH6-1, 2003), we investigated the occurrence of interplanetary coronal mass ejections in the near-Earth solar wind during 1996 – 2002, corresponding to the increasing and maximum phases of solar cycle 23, and provided a “comprehensive” catalog of these events. In this paper, we present a revised and updated catalog of the ≈300 near-Earth ICMEs in 1996 – 2009, encompassing the complete cycle 23, and summarize their basic properties and geomagnetic effects. In particular, solar wind composition and charge state observations are now considered when identifying the ICMEs. In general, these additional data confirm the earlier identifications based predominantly on other solar wind plasma and magnetic field parameters. However, the boundaries of ICME-like plasma based on charge state/composition data may deviate significantly from those based on conventional plasma/magnetic field parameters. Furthermore, the much studied “magnetic clouds”, with flux-rope-like magnetic field configurations, may form just a substructure of the total ICME interval.     

Propagation of the Bastille Day 2000 CME Shock in the Outer Heliosphere

The Bastille Day (14 July) 2000 CME is a fast, halo coronal mass ejection event headed earthward. The ejection reached Earth on 15 July 2000 and produced a very significant magnetic storm and widespread aurora. At 1 AU the Wind spacecraft recorded a strong forward shock with a speed jump from ∼ 600 to over 1000 km s⁻¹. About 6 months later, this CME-driven shock arrived at Voyager 2 (∼ 63 AU) on 12 January 2001 with a speed jump of ∼ 60 km s⁻¹. This provides a good opportunity to study the shock propagation in the outer heliosphere. In this study, we employ a 2.5-D MHD numerical model, which takes the interaction of solar wind protons and interstellar neutrals into account, to investigate the shock propagation in detail and compare the model predictions with the Voyager 2 observations. The Bastille Day CME shock undergoes a dramatic change in character from the inner to outer heliosphere. Its strength and propagation speed decay significantly with distance. The model results at the location of Voyager 2 are in good agreement with in-situ observations. Supplementary material to this paper is available in electronic form at http://dx.doi.org/10.1023/A:1014293527951     

> 25 MeV Proton Events Observed by the High Energy Telescopes on the STEREO A and B Spacecraft and/or at Earth During the First ∼ Seven Years of the STEREO Mission

Using observations from the High Energy Telescopes (HETs) on the STEREO A and B spacecraft and similar observations from near-Earth spacecraft, we summarize the properties of more than 200 individual >?25 MeV solar proton events, some detected by multiple spacecraft, that occurred from the beginning of the STEREO mission in October 2006 to December 2013, and provide a catalog of these events and their solar sources and associations. Longitudinal dependencies of the electron and proton peak intensities and delays to onset and peak intensity relative to the solar event have been examined for 25 three-spacecraft particle events. Expressed as Gaussians, peak intensities fall off with longitude with σ=47±14^° for 0.7?–?4 MeV electrons, and σ=43±13^° for 14?–?24 MeV protons. Several particle events are discussed in more detail, including one on 3 November 2011, in which ～?25 MeV protons filled the inner heliosphere within 90 minutes of the solar event, and another on 7 March 2012, in which we demonstrate that the first of two coronal mass ejections that erupted from an active region within ～?1 hour was associated with particle acceleration. Comparing the current Solar Cycle 24 with the previous cycle, the first >?25 MeV proton event was detected at Earth in the current solar cycle around one year after smoothed sunspot minimum, compared with a delay of only two months in Cycle 23. Otherwise, solar energetic particle event occurrence rates were reasonably similar during the rising phases of Cycles 23 and 24. However, the rate declined in 2013, reflecting the decline in sunspot number since the peak in the northern-hemisphere sunspot number in November 2011. Observations in late 2013 suggest that the rate may be rising again in association with an increase in the southern sunspot number.     

Identification of Interplanetary Coronal Mass Ejections at Ulysses Using Multiple Solar Wind Signatures

Previous studies have discussed the identification of interplanetary coronal mass ejections (ICMEs) near the Earth based on various solar wind signatures. In particular, methods have been developed of identifying regions of anomalously low solar wind proton temperatures (T _p) and plasma compositional anomalies relative to the composition of the ambient solar wind that are frequently indicative of ICMEs. In this study, similar methods are applied to observations from the Ulysses spacecraft that was launched in 1990 and placed in a heliocentric orbit over the poles of the Sun. Some 279 probable ICMEs are identified during the spacecraft mission, which ended in 2009. The identifications complement those found independently in other studies of the Ulysses data, but a number of additional events are identified. The properties of the ICMEs detected at Ulysses and those observed near the Earth and in the inner heliosphere are compared.     

Improving Stable Boundary-Layer Height Estimation Using a Stability-Dependent Critical Bulk Richardson Number

For many decades, attempts have been made to find the universal value of the critical bulk Richardson number ( $Ri_{Bc}$ ; defined over the entire stable boundary layer). By analyzing an extensive large-eddy simulation database and various published wind-tunnel data, we show that $Ri_{Bc}$ is not a constant, rather it strongly depends on bulk atmospheric stability. A (qualitatively) similar dependency, based on the well-known resistance laws, was reported by Melgarejo and Deardorff (J Atmos Sci 31:1324–1333, 1974) about forty years ago. To the best of our knowledge, this result has largely been ignored. Based on data analysis, we find that the stability-dependent $Ri_{Bc}$ estimates boundary-layer height more accurately than the conventional constant $Ri_{Bc}$ approach. Furthermore, our results indicate that the common practice of setting $Ri_{Bc}$ as a constant in numerical modelling studies implicitly constrains the bulk stability of the simulated boundary layer. The proposed stability-dependent $Ri_{Bc}$ does not suffer from such an inappropriate constraint.     

Binary near-Earth asteroid formation: Rubble pile model of tidal disruptions

We present numerical simulations of near-Earth asteroid (NEA) tidal disruption resulting in bound, mutually orbiting systems. Using a rubble pile model we have constrained the relative likelihoods for possible physical and dynamical properties of the binaries created. Overall 110,500 simulations were run, with each body consisting of ∼1000 particles. The encounter parameters of close approach distance and velocity were varied, as were the bodies' spin, elongation and spin axis direction. The binary production rate increases for closer encounters, at lower speeds, for more elongated bodies, and for bodies with greater spin. The semimajor axes for resultant binaries are peaked between 5 to 20 primary radii, and there is an overall trend for high eccentricity, with 97% of binaries having e > 0.1. The secondary-to-primary size ratios of the simulated binaries are peaked between 0.1 and 0.2, similar to trends among observed asteroid binaries. The spin rates of the primary bodies are narrowly distributed between 3.5- and 6-h periods, whereas the secondaries' periods are more evenly distributed and can exceed 15-h periods. The spin axes of the primary bodies are very closely aligned with the angular momenta of the binary orbits, whereas the secondary spin axes are nearly random. The shapes of the primaries show a large distribution of axis ratios, where those with low elongation (ratio of long and short axis) are both oblate and prolate, and nearly all with large elongation are prolate. This work presents results that suggest tidal disruption of gravitational aggregates can make binaries physically similar to those currently observed in the NEA population. As well, tidal disruption may create an equal number of binaries with qualities different from those observed, mostly binaries with large separation and with elongated primaries.     

A fast method for finding bound systems in numerical simulations: Results from the formation of asteroid binaries

We present a new code (companion) that identifies bound systems of particles in O(NlogN) time. Simple binaries consisting of pairs of mutually bound particles and complex hierarchies consisting of collections of mutually bound particles are identifiable with this code. In comparison, brute force binary search methods scale as O(N²) while full hierarchy searches can be as expensive as O(N³), making analysis highly inefficient for multiple data sets with N?¹⁰³. A simple test case is provided to illustrate the method. Timing tests demonstrating O(NlogN) scaling with the new code on real data are presented. We apply our method to data from asteroid satellite simulations [Durda et al., 2004. Icarus 167, 382-396; Erratum: Icarus 170, 242; reprinted article: Icarus 170, 243-257] and note interesting multi-particle configurations. The code is available at http://www.astro.umd.edu/zoe/companion/ and is distributed under the terms and conditions of the GNU Public License.     

Snow density variations: consequences for ground‐penetrating radar

Reliable hydrological forecasts of snowmelt runoff are of major importance for many areas. Ground‐penetrating radar (GPR) measurements are used to assess snowpack water equivalent for planning of hydropower production in northern Sweden. The travel time of the radar pulse through the snow cover is recorded and converted to snow water equivalent (SWE) using a constant snowpack mean density from the drainage basin studied. In this paper we improve the method to estimate SWE by introducing a depth‐dependent snowpack density. We used 6 years measurements of peak snow depth and snowpack mean density at 11 locations in the Swedish mountains. The original method systematically overestimates the SWE at shallow depths (+25% for 0·5 m) and underestimates the SWE at large depths (?35% for 2·0 m). A large improvement was obtained by introducing a depth–density relation based on average conditions for several years, whereas refining this by using separate relations for individual years yielded a smaller improvement. The SWE estimates were substantially improved for thick snow covers, reducing the average error from 162 ± 23 mm to 53 ± 10 mm for depth range 1·2–2·0 m. Consequently, the introduction of a depth‐dependent snow density yields substantial improvements of the accuracy in SWE values calculated from GPR data. Copyright © 2005 John Wiley & Sons, Ltd.