Mapping Earthquakes with ArcGIS Pro

Jeff Oppong


When an earthquake strikes, it can lead to loss of lives, properties, and resources. Geospatial tools like ArcGIS Pro can be used by map out both fault lines and seismic activity.

What to expect with this tutorial on earthquake mapping

This GIS tutorial teaches you how to use USGS data and ArcGIS Pro to map earthquakes and geological faults using the area around Dallas, Texas as an example. This tutorial uses Landsat 8 images and earthquake data from the USGS website. PCI Geomatica software is used to narrow the geographic area within the image and pull out features like geological fault lines.

This fault line GIS data, along with the earthquake data from the USGS, can then be imported into ArcGIS Pro software. Next, using ArcGIS Pro, the earthquake data will be interpolated using the “spline,” tool in order to visually represent earthquake magnitude as a raster layer. Finally, we will create a map that shows both the earthquake strength and the fault lines. Figure 1 shows the flow map of the tools and data you will use to make this map in ArcGIS Pro using Landsat 8 images and earthquake data.

A black and white flow chart for the tools and data needed to map earthquakes.
Figure 1: Geospatial tools and data used to create an earthquake map showing geological faults.

Landsat 8 and preparing geological faults

Landsat 8 imagery can be accessed freely from the United States Geological Surveys (USGS)’s website after a free sign up. For this tutorial Landsat 8 (collection 2 level 1) imagery is used. For more information on how to access the imagery, watch this video. Geological faults is extracted from Landsat 8 imagery using the PCI Geomatica software. For information on how a free trial and how to used to download the PCI software, visit Catalyst Professional – CATALYST.Earth.

Free weekly newsletter

Fill out your e-mail address to receive our newsletter!

Extracting geological faults from PCI Geomatica

  1. Launch and create a new project. 
A screenshot showing an empty project in PCI geomatica.
Figure 2: PCI Geomatica interface.
  1. From the menu bar, click tools.
  2. Click clipping/subsetting from the set of options displayed. 
  3. From the clipping/subsetting window, set input file to the location of Landsat 8 band 10.
  4. Set output file to name and location of your clipped dataset.
  5. Be sure to enable the input data, by checking the box of the raster (Landsat 8 band 10) under the Available layers section. 
  6. Under the output raster size, click square to define the boundary of the clipped dataset.
  7. Click clip.  
Screenshot from PCI Geomatic showing a raster image that is a grayscale gradient.
Figure 3: Clipping/Subsetting in PCI Geomatica.

After clipping, geological faults are extracted from the clipped data by following the steps outlined below:

  1. From tools on the menu bar, select “algorithm librarian” from the drop-down menu.
  2. Select Analysis.
  3. Select Geological/geophysical analysis.
  4. Select Line/lineament extraction.
  5. Select Open.
  6. Set input file to the saved clipped file. 
  7. Click Run.
Screenshot showing the tools interface for extracting faults lines in PCI geomatica.
Figure 4: Line extraction in PCI Geomatica.
  1. From the results, right-click on “Extracted lineaments” and select save as. 
  2. From the output location section, set output name and location.
  3. Be sure to set file type to “ArcView shapefile.” Click Save.  
Screenshot showing the tool to use to export a shape filed from PCI Geomatica.
Figure 5: Saving geological faults as a shapefile in PCI Geomatica.

Finding earthquake GIS data for a region of interest

For this tutorial, geological fault data helps to validate the earthquake data for a geographic area. Earthquake data can be retrieved from the USGS’s Earthquake Catalog using the steps outlined below:

  1. Set magnitude to 2.5+.
  2. Use the Draw rectangle to specify the geographic region.
Screenshot showing the USGS Earthquake Catalog search interface.
Figure 5: The USGS Earthquake Catalog.
  1. From the output options, set format to CSV.
  2. Set order by to Time-Oldest First.
  3. Click Search.

Creating an earthquake map showing geological faults in ArcGIS Pro

After retrieving the earthquake data from USGS and extracting the geological fault lines from Landsat 8 (band 10), use the following steps to create a map in ArcGIS Pro:

  1. Open and create a new project in ArcGIS Pro. 
  2. Select Add data. 
  3. Toggle to the location of the USGS earthquake data that you retrieved.
  4. Click Ok. 
Screenshot showing the browser window for adding data in ArcGIS Pro.
Figure 6: Add data in ArcGIS Pro.
  1. Right-click on the query table.
  2. From the Display XY data, set X field to longitude. Set Y field to latitude. 
  3. Click Ok. 
  4. Type and search Spline from the search bar of geoprocessing toolbox. 
  5. From the spline window, set input point features to downloaded earthquake data.
  6. Set Z value to “mag.”
  7. Click Run to create a raster.

Now add the geological fault data processed from Landsat 8 (band 10) as a overlay to your raster dataset of earthquake magnitude:

  1. Click the Add data tool from the Map tab of the menu bar.
  2. Toggle to the location of the fault data.
  3. Click Ok. 
Screenshot showing a raster layer interpolated from the magnitude value of an earthquake data set with fault lines.
Figure 9: Map showing geological faults and a raster layer of the magnitude of earthquakes in the Dallas, Texas area.
Photo of author
About the author
Jeff Oppong
Jeff Oppong holds a BSc in Geomatic Engineering and currently a graduate student at Hohai University in China, where he's studying MSc. Harbor, Coastal, and Offshore Engineering. Jeff is a prolific researcher and a GIS/Remote sensing expert who aspires to be a change-agent and a renowned Engineer. Jeff Oppong can be contacted via email