My main source of data was the real time USGS GeoJSON earthquake feed. This feed provides a feature collection of points indicating the epicenters of recent earthquakes. From this data, I was able to retrieve intensity readings for each survey center that detected the earthquake. Loading this data into the dataset editor, I then created polygons representing the different intensity zones.
At this point the data was still too hard to read, so I made a dataset of hexagons covering the affected area. Each hexagon was then correlated to the highest intensity zone the hexagon covered from my first pass at the data. The result is what you see below.