Risk Assessment of Potential Meteor Impacts

Sam O'Mullane
Posted on Apr 30, 2017

What are the chances 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.

Two parameters are generally used to classify these objects:

  1. Semi-major axis - the longer of the two radii that define an elliptical orbit
  2. Perihelion distance - the closer of the two distances from the object to the sun along the major axis

Orbital Parameter Definitions

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.

NEO orbit schematics

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.

kepler-law

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:

a-TJN-classq-TJN-classd-TJN-class

 

 

 

 

 

 

 

 

 

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.

 

orbit-TJN

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:

a-ATE-classq-ATE-classd-ATE-class

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!

orbit-ATE

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:

a-NEOq-NEOd-NEO

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.

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.

Screen Shot 2017-04-30 at 12.22.59 PM

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.

For a much less likely less scenario, let's see what happens if the largest PHA hits earth: 2010 AU118.

Screen Shot 2017-04-30 at 12.24.07 PM

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.

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.

impactProb

While this is not comforting, the following statistics should be:

  1. The probability of one small meteor (2010 RF12, diameter = 0.007 km) hitting Earth is 5%
  2. The probability of at least one meteor hitting Earth is rough 8.5%
  3. The most likely estimated year of impact is between 50 and 100 years from now

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

About Author

Sam O'Mullane

Sam O'Mullane

Sam O'Mullane is a multidisciplinary scientist who is passionate about making data science accessible to a wide audience. Sam received his MSc in Quantum Fields and Fundamental Forces from Imperial College London (Th. Physics) in 2013 and PhD...
View all posts by Sam O'Mullane >

Related Articles

Leave a Comment

Avatar
Google October 17, 2019
Google Here are some links to sites that we link to for the reason that we feel they are really worth visiting.
Avatar
Google September 26, 2019
Google The data talked about in the post are a few of the top out there.
Avatar
Louvenia July 14, 2017
This is a matter close to my heart cheers. Thanks
Avatar
Chance July 11, 2017
Great post. I am facing a couple of these problems.

View Posts by Categories


Our Recent Popular Posts


View Posts by Tags

2019 airbnb alumni Alumni Interview Alumni Spotlight alumni story Alumnus API Application artist aws beautiful soup Best Bootcamp Best Data Science 2019 Best Data Science Bootcamp Big Data Book Launch Book-Signing bootcamp Bootcamp Prep Bundles California Cancer Research capstone Career Career Day citibike clustering Coding Course Demo Course Report D3.js data Data Analyst data science Data Science Academy Data Science Bootcamp Data Scientist Data Scientist Jobs data visualization Deep Learning Demo Day Discount dplyr employer networking feature engineering Finance Financial Data Science Flask gbm Get Hired ggplot2 googleVis Hadoop higgs boson Hiring hiring partner events Industry Experts Job Jon Krohn JP Morgan Chase Kaggle Kickstarter lasso regression Lead Data Scienctist Lead Data Scientist leaflet linear regression Logistic Regression machine learning Maps matplotlib Medical Research meetup Networking neural network Neural networks New Courses nlp NYC NYC Data Science nyc data science academy NYC Open Data NYCDSA NYCDSA Alumni Open Data painter pandas Part-time Portfolio Development prediction Prework Programming PwC python python machine learning python scrapy python web scraping python webscraping Python Workshop R R language R Programming R Shiny r studio R Visualization R Workshop R-bloggers random forest recommendation recommendation system regression Scrapy scrapy visualization seaborn Selenium sentiment analysis Shiny Shiny Dashboard Spark Special Special Summer Sports statistics streaming Student Interview Student Showcase SVM Tableau TensorFlow Testimonial tf-idf Top Data Science Bootcamp twitter visualization web scraping Weekend Course What to expect word cloud word2vec XGBoost yelp