Survey
* Your assessment is very important for improving the workof artificial intelligence, which forms the content of this project
* Your assessment is very important for improving the workof artificial intelligence, which forms the content of this project
The Estimated Population of Small NEOs Alan Harris MoreData! Inc. Target NEO 2 Washington, DC, July 9, 2013 Estimating completion from re-detection ratio Current survey completion Completion or Re-detection ratio 1.0 Survey sim. re-detections Survey sim. completion 0.8 0.6 0.4 0.2 0 -8 -7 -6 -5 -4 -3 -2 -1 0 dm 1 2 3 4 5 6 7 8 For a computer modeled survey with a known number of synthetic NEAs, we can tabulate the re-detection ratio vs. size (dm) as well as the real completion. For a real survey, we can measure redetection ratio, but not completion, since we don’t know the total population. Estimating completion from re-detection ratio Current survey completion Completion or Re-detection ratio 1.0 Survey sim. re-detections Survey sim. completion Re-detections, 2010-2012 0.8 0.6 0.4 0.2 dm 0 -8 -7 -6 -5 -4 -3 -2 -1 0 1 2 3 4 5 6 7 8 29 28 27 26 25 24 23 22 21 20 19 18 17 16 15 14 13 H magnitude The actual redetection ratios for the combination of LINEAR, Catalina, and Siding Spring can be adjusted horizontally to match the model curve to within the uncertainties in the survey data. We thus adopt the model completion curve as adjusted to match the dm scale as representing current survey completion. Completion Extrapolation to very small size The observed re-detection ratio becomes uncertain below about 0.1 (that is, H greater than about 22) due to the low number of re-detections. However, having “calibrated” the completion curve in the range of good re-detection statistics, we can extend to still smaller sizes by assuming that the computer completion curve accurately models actual completion. This works until the number of “detections” in the computer model falls below a statistically useful number, say about 100 “detections” out of the 100,000 model asteroids, or a completion of about 10-3. This corresponds to about dm = -4.0, or on the scaled curve to about H = 25.0. Fortunately, below dm of ~3.0, -1 detections are close to the Earth and 10 can be modeled with rectilinear motion rather than accounting for orbital -3 10 motion. An analytical completion function can be matched to the computer completion curve and -5 10 Computer model completion extrapolated to arbitrarily small size. Theoretical log(C) -0.8dm -7 10 -11 -10 -9 -8 -7 -6 -5 -4 -3 -2 -1 0 dm = Vlim - H 1 2 3 4 5 6 7 With these extensions, we now have an estimate of completion over the entire size range of observed objects. Differential Population Diameter, km 0.01 0.1 1 10 8 10 n io e at tim es n(H) t la pu Po 6 10 4 10 2 10 Total discovered 0 10 31 30 29 28 27 26 25 24 23 22 21 20 19 18 17 16 15 14 13 12 11 10 9 H Plotted here are the numbers in each halfmagnitude interval, in red the total number discovered as of August 2012, and in blue the estimated total population in that size range, based on the completion curves of the previous graphs. Cumulative Population Impact Energy, MT -1 10 2 10 10 5 10 8 10 10 8 10 0 10 2 10 6 K-T Impactor N(<H) 10 4 10 2 10 8 0 31 30 29 28 27 26 25 24 23 22 21 20 19 18 17 16 15 14 13 12 11 10 9 0.01 0.1 1 Diameter, Km 6 10 10 Absolute Magnitude, H 10 4 10 10 Impact Interval, years Tunguska Estimated, 2012 Annual bolide event Constant power law Discovered to August, 2012 The cumulative population is the running sum of the differential population, from the previous plot. The number N is the total number of NEAs larger than the specified size (H or Diameter). Population of very small missionsuitable NEAs In the material presented so far, it can be seen that there are approximately 108 NEAs in the size range 5-10 m diameter, of which only about 100 have been found. This would imply that there are abundant numbers of suitable mission targets that could be discovered. But not all of the 108 NEAs are suitable mission targets, in fact only a very tiny fraction are: • Over half do not cross the Earth’s orbit or even approach very closely. Only about 1/5 are PHA-class (MOID < 0.05 AU). • The median encounter velocity with the Earth (v) is ~20 km/sec; the fraction with suitably low v (< 2.6 km/sec) is at most a few tenths of a percent. Velocity (v) distribution of small Earthcrossing asteroids (ECAs) Velocity distribution of ECAs 1.0 Fraction < v 0.8 0.6 > 1 km intrinsic distribution (model) > 1 km discovered fraction (actual) < 40 m discovered fraction (model) < 40 m discovered fraction (actual) Maxwellian distribution 0.4 0.2 0 0 10 20 30 v 40 We expect the velocity distribution of the entire population to be roughly homologous over size. The dramatic difference among discovered small objects is mostly or entirely due to the relative ease of discovering low v objects. Very low v NEAs Most ECAs arrive from the main asteroid belt through perturbations leading to multiple planet orbit crossings, hence most have relative encounter velocities with the Earth that are sufficient, if properly directed, to reach to the neighboring planets, Mars or Venus. That turns out to be, for both neighbors, about v 2.5 km/sec. Multiple encounters with only the Earth-moon system can evolve v to lower values, but require a “Maxwell’s Demon” (e.g., spacecraft navigation) to do so effectively. Thus, there is a “hole” in the velocity distribution below ~2.5 km/sec. There are about 30 discovered objects with v 2.5 km/sec, so the phase space is not entirely empty. However, one must be mindful that there are three plausible sources of such objects: • “Thermalized” tail of the v distribution of regular NEAs from the main asteroid belt • Lunar ejecta (Gladman et al.1996) • Space debris v Close-up of fitted Maxwellian distribution Close-up of low range of fitted Maxwellian distribution Fraction < v 0.010 Minimum v to get to/from Venus or Mars 0.005 0 0 1 2 3 4 v This suggests that the “thermalized” tail of the velocity distribution of real NEAs with v < 2.5 km/sec is only ~0.002 of the total population. EarthEarthcollision time scale vs. v collision time scale for low v ECAs Earth collision time scale, years 10 8 10 7 10 6 10 5 10 4 10 3 0.5 Minimum v to get to/from Venus or Mars 1 2 5 v But the Earth collision time scale of v < 2.5 km/sec objects is shorter than the diffusion time scale to populate that phase space. So the Earth is gobbling up much of what chances to make it there, and the steady state population is likely even less 10 than “Maxwellian”. Evolution of orbits of Lunar ejecta Gladman et al., Science 271, 1387-1392, 1996. Orbits of low-v NEAs a,e distribution of low-v NEAs 1.3 On the left is the (a,e) distribution of discovered NEAs. On the right is a plot of Gladman et al. orbits of Lunar ejecta after 105 years, matched to the same scale. The lack of discovered NEAs interior to the Earth’s orbit is likely a selection effect against discovering such objects. 1.2 a 1.1 1.0 0.9 0.8 0 0.05 0.10 0.15 e Evolution of v 2.5 km/sec orbits Gladman et al. consider the transfer of Martian ejecta to the Earth, which yields some insight into the diffusion of orbits from the main belt into Earthcrossing orbits of low v. Unfortunately a coordinate is missing, inclination, so we can only infer minimum values of v, corresponding to zero inclination. Only one object in the above panels, the one at a = 0.95, e = 0.05 in the 106 year panel, would have a v < 3 km/sec even at zero inclination, and that one object would have v > 2.5 km/sec if i > 5. Thus, out of the 200 particles simulated, at most one, and probably none, evolved into an Earth-crossing orbit with v < 3 km/sec. The prospects for objects of low v coming from the main asteroid belt are perhaps even less likely. Space Debris We know at least one discovered object of low v has been confirmed to be space debris; several others are likely old rocket bodies. Sources of ultra-low v NEAs 1. 2. 3. 4. Lunar ejecta (most) Space debris (some) Main-belt asteroids (almost none) Mars ejecta (almost none) The take-home message is that in choosing a very low v target, you need to have very good physical characterization of the object if you want to be sure you aren’t bringing a piece of the moon back to its home, or even an old rocket body. Binary and spin properties of small NEAs 4,775 reliable asteroid rotation periods (November, 2012) 0.001 • There are no binaries below ~300 m diameter. • All objects smaller than ~20 m are super-fast rotators. • Many small objects are tumblers. Rotation period, hours 4596 single, PA rotators 33 Tumblers 146 Binaries 0.1 10 1000 0.001 0.1 10 Diameter, km 1000 Better know before you go.