Data Visualization on Risk Assessment of Meteor Impacts
The skills the author demoed here can be learned through taking Data Science with Machine Learning bootcamp with NYC Data Science Academy.
What are the chances, based on data,Β Earth is impacted by a meteor?
If a meteor does hit Earth, how much damage would it cause?
The meteor applicationΒ I created attempts to answer these questions as well as provideΒ background informationΒ on solar system objects.
Introduction
There are over 700,000 asteroids and comets that have been observedΒ in our solar system. All key information about these objects has been stored by the Jet Propulsion Laboratory (JPL) in their small-body database.
Parameters
Two parameters are generally used to classify these objects:
- Semi-major axis - the longer of the two radii that define an elliptical orbit
- Perihelion distance - theΒ closer of the two distances from the object to the sun along the major axis
Since the small-body database contains the class information about each object, the potential threat of each class can be assessed by charting the mean and median semi-major axis and perihelion distance.Β More specifically, an object makes a close approach when the fly-by distance from Earth is less than 0.05 AU (4.6 million miles) and is considered hazardous when the diameter is larger than 330 ft. Out of twentyΒ classes, only four are potentially dangerous, as pictured below (colors) crossing the Earth (inner) and Mars (outer) orbits.
More specific information about hazardousΒ objects is contained in another source: the Sentry - Earth impact monitoring database. This database contains information on object velocity and diameter which will be used to assess potential damage.
For the curious, the small-body database serves as a great demonstration of Kepler's third law. This relation between semi-major axis and orbital period served in part as a basis for Newton's work on gravity.
Asteroid and Comet Classification
As mentioned in the introduction, only a small subset of objects are potentially dangerous to Earth. Specifically, 95% of objects are contained within the asteroid belt between Mars and Jupiter while roughly 0.1% of objects have been determined to be potentially hazardous. My applicationΒ can be used to browse the various classes and explore why they are or are not potentially hazardous.
For example, theΒ Jupiter Trojan objects are described as hovering around Jupiter's Lagrangian points with semi-major axes around 5.2 AU (according to Wikipedia). By selecting the TJN class in the application, the followingΒ semi-major axis (left), perihelion distance (center) and diameter (right) plots are produced:
Observations
For reference, the dashed lines indicate mean of the distribution while the dotted lines indicate the median of the distribution. Using the mean values, we can plot the average expected orbit for each class to understand approximately where these objects travel on a macro level. For the TJN class (black), the following image is produced with Jupiter (purple), Mars (red) and Earth (blue) orbits are plotted for reference. We see that these objects are outside of the asteroid belt and are extremely far away from Earth, hence not dangerous in the least.
In contrast, the followingΒ semi-major axis (left), perihelion distance (center) and diameter (right) plots are produced for Aten class objects which are potentially dangerous, but much smaller in diameter:
Charting these orbital parameters in terms of mean orbits produces the image below. No wonder that these are hazardous, their average orbit crosses Earth's orbit!
Potentially Hazardous Objects
With an understanding of which objects are potentially hazardous asteroids (PHA), we can turn to objects which have a non-zero probability of impacting Earth. The distributions of these objects are plotted below:
Of main interest is the diameter distribution since the size of the meteor is proportional to the damage it will cause (assumingΒ that the object densities are similar). The median of this distribution is roughly 60 feet while the mean is roughly 130 feet.
Potential Impacts Page
In astronomical terms, these objects are tiny. How much damage would one of these median-size objects cause? Using the search feature in the datatable tab of the potential impacts page, we can locate an example asteroid: 2009 BE.
The crater application takes in a meteor, target and meteor material information, impact city (visualization only) and angle of impact to determine key figures such as crater radius, depth and time for formation while displaying key statistics about the meteor itself. For 2009 BE, we see that there is a relatively high impact probability of 43 thousandths of a percent though the resulting crater would only be 300 feet wide and 200 feet deep.
2010 AU118
For a much less likely less scenario, let's see what happens if the largest PHA hits earth: 2010 AU118.
We see that the crater will be two orders of magnitude larger and thus the destruction would likely be devastating. Not included in this calculation are the other effects that a meteor would cause such as pressure wave generation andΒ large scale dust generation and corresponding climate change. Given more time, this application could be updated to include more disaster-related calculations in addition to the crater generationΒ figures.
Impact Probability vs Diameter
So how likely is it that some massive meteor impacts Earth? The hope is that there is a negative correlation between impact probability and diameter. This plot is shown below (both on log scales) and while there is a general trend, there is a large spread about the least squares regression line with some notable outliers in high probability and large diameter.
While this is not comforting, the following statistics should be:
- The probability of one small meteor (2010 RF12, diameter = 0.007 km) hitting Earth is 5%
- TheΒ probability of at least one meteor hitting Earth is rough 8.5%
- The most likely estimated year of impact is between 50 and 100 years from now
Conclusion
In conclusion,Β there is very little to worry about regarding meteor impact especially considering the small fraction of the world that has significant population density. In the future, an estimation of the probability of a meteor impacting a populated area will be included.
The relevant equations and coefficients for the crater formation problem can be found in the following paper:
- Holsapple, K.A., 1993. The Scaling of Impact Processes in Planetary Sciences. Annu Rev Earth Pl Sc 21, 333β373