Skip to main content
Geosciences LibreTexts

4.6: Earthquake Catalog Data Analysis Problem Set

  • Page ID
    7009
  • \( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

    \( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

    \( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)

    ( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)

    \( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

    \( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)

    \( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

    \( \newcommand{\Span}{\mathrm{span}}\)

    \( \newcommand{\id}{\mathrm{id}}\)

    \( \newcommand{\Span}{\mathrm{span}}\)

    \( \newcommand{\kernel}{\mathrm{null}\,}\)

    \( \newcommand{\range}{\mathrm{range}\,}\)

    \( \newcommand{\RealPart}{\mathrm{Re}}\)

    \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

    \( \newcommand{\Argument}{\mathrm{Arg}}\)

    \( \newcommand{\norm}[1]{\| #1 \|}\)

    \( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

    \( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)

    \( \newcommand{\vectorA}[1]{\vec{#1}}      % arrow\)

    \( \newcommand{\vectorAt}[1]{\vec{\text{#1}}}      % arrow\)

    \( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

    \( \newcommand{\vectorC}[1]{\textbf{#1}} \)

    \( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)

    \( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)

    \( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)

    \( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

    \( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

    \(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)

    For the following problem set, you will work with the seismicity catalog maintained by the University of Memphis for the New Madrid region in order to make your own frequency-magnitude diagrams and calculate a recurrence interval for a large earthquake at the NMSZ. You will also compare NMSZ data to seismicity catalogs for southern California and the world. The point of this comparison is that you will see that the overall shape of a frequency-magnitude diagram is scale-independent. It doesn't matter how big your regional area is, or how many years your catalog covers, the same basic -1 slope coupled to two sections on either end that deviate from -1 will always be there. What changes is the place on the diagram where the deviation occurs? Pay attention to this when you make your different plots.

    Part 1: Get familiar with the NMSZ and its seismicity

    Go to the New Madrid Earthquake Catalog Search.

    Once you get there, follow my directions to make a 1-year catalog of NMSZ seismicity:

    Video: How to Make a 1-year Catalog of NMSZ Seismicity (2:18)

    How to Make a 1-year Catalog of NMSZ Seismicity

    Click here for the transcript of How to Make a 1-year Catalog of NMSZ Seismicity.

    PRESENTER: So you've landed on the New Madrid Earthquake Catalog search page. And here's how we're going to get our first catalog of earthquakes out of here. First, I'm going to change the year to 2011, and then in this box, I'm going to change it to 2012. And you'll notice that this gives me a full year of earthquakes. It's actually not crucial that you pick the exact same year I did, but in this example, I'm just going to pick 2011.

    I do want you to get a full year. So now in the magnitude range, we'll just leave this alone. The lat long box search, we are going to change these numbers around slightly from the defaults. The minimum latitude we're going to use is 34 degrees, and the maximum latitude we're going to use is 39 degrees. The minimum longitude we are going to use is negative 92 degrees. That's 92 degrees west. And the maximum longitude is going to be negative 87 degrees. So that's 87 degrees west.

    Now, leave the radial search numbers alone, and we are ready to just click Begin Search. So now we pop up the catalog search results page. You can click here, and you'll see your catalog. You can go to the File menu and save the page as some plain text suitable for importing into your favorite plotting program.

    You go back to the catalog search results page by clicking the Back button on your browser, you can also read a description of the catalog information which tells you what's in every column of your text file. And I also want you to generate a map. You can play around with all of these options on your own, but clicking generate map will generate a map of the catalog you've just made. Looks like this. I want you to save this map somewhere to your computer and paste it into your problem set for turning in later on. That's all there is to it.

    Credit: Dutton

    Here is video above as plain text

    Create a word processing document (Microsoft Word, Macintosh Pages, Google Docs, or PDF) to record your work for this problem set.

    On the worksheet, paste in your map of the 1-year catalog you made. Then answer the following questions:

    1.1 How many earthquakes are in your catalog?

    1.2 What is the largest magnitude earthquake in your catalog? How many earthquakes are there of this magnitude in your catalog?

    1.3 What's the smallest magnitude earthquake in your catalog? How many earthquakes are there of this magnitude in your catalog?

    1.4 Describe your map in a few sentences. (What part of the country is it? Are the earthquakes sprinkled randomly about or do they cluster in patterns? If the latter, describe what the patterns look like.)

    Double-Check! Your worksheet should now have a map and answers to the Part 1 questions. If it does, you are ready to take on Part 2.

    Part 2: Making frequency-magnitude plots of NMSZ seismicity

    Make three different frequency-magnitude plots using the New Madrid Earthquake Catalog data: The first plot will use the one-year catalog you made in Part 1 of the problem set. In the second plot, you will add curves that correspond to a 10-year, 20-year, and 30-year catalog. The third plot will depict the same data as in the second plot, except that you will normalize all the catalogs to one year. Specific directions follow:

    Plot 1: Make a frequency-magnitude plot using the one-year-long catalog that you just made in Part 1 of this problem set.

    For this plot, all you need are the magnitudes, which are in column 7 of your plain text catalog file. Post to Questions if you need help isolating that column. My recommendation is to create a new file in whatever plotting/spreadsheet program you like and then you will type your counting statistics into this new file. In your new file, you want to create the values that will be on the x-axis of your plot. They will be magnitudes from your catalog's lowest to your catalog's highest in 0.1-unit bins. To create the values that will be on the y-axis of your plot, count how many earthquakes in your catalog are equal to or greater than each value of magnitude. These are your y-values: cumulative frequency. It's easier to count if you sort your magnitudes first, which is fine to do because we do not care what order they are in for this plot.

    Then plot cumulative frequency vs. magnitude.

    Use a logarithmic (base 10) scale for the y-axis or take the log of your cumulative frequency data and plot that on a linear axis. You can use a linear x-axis because magnitude is already a power of 10. [Is this confusing? Post to the Questions forum if you need help.]

    If you are having trouble, see my example plot for a one-year catalog of earthquakes in the New Madrid Seismic Zone (NMSZ). I used the year 1975, so this plot will not be precisely the same as yours, but it should look pretty close.

    Plot 2: Make a plot that has four different frequency-magnitude curves on it.

    Start with the plot you just made.

    Go back to the New Madrid Earthquake Catalog and make a catalog for a ten-year time period (you can choose any ten-year period). Overlay the frequency-magnitude curve for this ten-year catalog onto the curve you made in the first plot.

    Repeat for a twenty-year period.

    Repeat for a thirty-year period.

    You should now have one plot with four curves on it.

    Make sure each curve is distinguishable (by color or linestyle) and labeled.

    Plot 3: Make a plot that has four different normalized frequency-magnitude curves on it.

    Start with the original one-year catalog plot you made.

    Overlay the curve for the ten-year catalog, but normalize the curve to one year. This is accomplished by dividing each of your y-values by 10.

    Repeat for the twenty-year catalog (divide y-values by 20)

    Repeat for the thirty-year catalog (divide y-values by 30)

    You should now have one plot with four curves on it.

    Use the same distinguishing color or linestyle for each curve as you did in your second plot.

    Dust off your worksheet from Part 1. On the worksheet, first paste your three plots into the worksheet, then answer the following questions:

    2.1 Look at the first frequency-magnitude plot you made of the one-year catalog. Approximately what is the lower magnitude threshold for this catalog? (Follow your data from right to left and tell me at about what magnitude the line flattens out to having zero slope?)

    2.2 Now look at the other two plots you made. Do the other curves show a significantly different lower magnitude threshold? From this observation, what do you conclude about the relationship between catalog timespan and lower magnitude sensitivity?

    2.3 Look at the second plot you made. Describe the differences and similarities among the four curves in a few sentences. For example, are the curves of the same shape? Where are the x- and y-intercepts relative to each other? What makes the y-intercepts different? What causes the x-intercepts to be different?

    2.4 Look at the second plot you made. Imagine having a catalog that spans 100 years. Using the four curves you made as a guide, extrapolate where the x and y intercepts would each be for a 100-year catalog. What is the largest earthquake you would expect for a 100-year catalog?

    2.5 Look at the third plot you made. Extrapolate your curves and predict how often a magnitude 7 earthquake will occur in this region and how often a magnitude 8 will occur in this region. I want you to make a reasonable eyeball-fit. I am not asking you to calculate a best fit line. **If the answer is a fraction less than one, then you can take the reciprocal and predict how many years go by in between magnitude 7's and in between magnitude 8's.

    2.6 You have just used frequency-magnitude relationships to predict a recurrence interval for a large New Madrid earthquake. Cool! What are the sources of uncertainty in the prediction you made in problem 2.5? (One way to realize just how much uncertainty there is in an extrapolation like this one is to try making several slightly different fits to the data that all look "pretty good" to you and see how different your final answers end up being.)

    Double-Check! Now your worksheet should have the map and answers from Part 1 as well as three plots and answers for Part 2. Hang on to the worksheet and use it for Part 3.

    Part 3: Comparing the NMSZ to Southern Cal and to the world

    Plot: Southern California

    Use the Southern California Earthquake Center Web site to make a seismicity catalog, map, and frequency-magnitude plot.

    Go to the Southern California Earthquake Data Center's Earthquake Catalog Search page.

    Video: How to Make a One-Year Catalog of Seismicity (2:10)

    How to Make a One-Year Catalog of Seismicity

    Click here for the transcript of How to Make a One-Year Catalog of Seismicity.

    PRESENTER: Here we are at the Southern California Earthquake Data Center catalog search page. What we're about to do is create a catalog that is the same length in time as our one year New Madrid catalog, and also covers the same latitude and longitude box area. So that way, we can really compare how active this area is in terms of number of earthquakes produced compared to New Madrid, because we'll be looking at the same size and length of time.

    All right, so the catalog to search is this one. The output format default is fine. Let's make a one year catalog for the same year as we did before. So let's make it 2011, day one, month one. Like that.

    We'll leave the magnitude and depth defaults alone. For latitude, let's change it to 33 and 38. And for longitude, we're going to change it to West 124 and West 119. I'm going to Shift Click to select local, regional, and teleseisms.

    And I want to send this output to a file. So you're going to save this as plain text somewhere on your computer so that you can manipulate it to make a frequency magnitude plot later on. Let's get rid of that right now.

    Now, I also want you to make a map. So in that case, what you do is go back up here to the output format, leaving everything else the same. And change the output to a Google map. Now you're going to want to output this to a web page. Submit request. Stare off into space while this thing happens.

    And you have a map. You might want to zoom out or center it differently. And then take a screenshot of this map and paste it into your problems set.

    Credit: Dutton

    Or follow along as a plain text to make a one-year catalog of seismicity.

    Take a screenshot of the google map you made and paste it into your problem set.

    Make a frequency-magnitude plot for the events in this catalog:

    • Start with your frequency-magnitude plot for the one-year NMSZ catalog.
    • Overlay the frequency-magnitude curve for your Southern California catalog on the same axes. The magnitudes are in column 5 of the SCEDC catalog file, and they are all you need.
    • Make sure each of your curves are distinguishable and labeled.
    • Save this plot because you are going to add another curve to it.

    Add to your plot: The World

    Use the United States Geologic Survey catalog to make a one-year global seismicity catalog, and add this data to your frequency-magnitude plot that has Southern Cal and the NMSZ on it.

    Go to the USGS Earthquake Search page.

    Once you are there, follow my plain text directions for making a one-year global catalog.

    Add to your plot!

    • Start with the frequency-magnitude plot you made that has a curve for the one-year NMSZ catalog and a curve for the SCEDC catalog.
    • Overlay the frequency-magnitude curve for the global catalog on the same axes.
    • Make sure each of your curves are distinguishable and labeled.
    • Paste this plot into your problem set.

    A few more questions

    3.1 How many earthquakes are in your one-year catalog for Southern California? What is the largest magnitude earthquake in the catalog? How many earthquakes are there of this magnitude?

    3.2 How many earthquakes are in your one-year catalog for the world? Are you surprised by this number? What is the largest magnitude earthquake in the world catalog? How many earthquakes are there of this magnitude? Remember that we cut off our global catalog at a minimum magnitude of 4.5. Look at your frequency-magnitude curve for the global catalog and estimate how many earthquakes there would be in your catalog if we had gone all the way down to zero for the minimum magnitude. Translate that into an approximate number of earthquakes per day in the world. (wow, huh!)

    3.3 Look at your map of Southern Californian earthquakes. Describe it in a few sentences (i.e., Where are the earthquakes? Do they cluster in space? Beware of artificial clustering that we induced by where we set our search parameters.).

    3.4 Look at the map of one year's worth of earthquakes made by the Advanced National Seismic System and describe it in a few sentences. How do the earthquakes cluster?

    3.5 Compare the frequency-magnitude curves for New Madrid and Southern California. Which one of the two catalogs is has its lower magnitude threshold at a smaller magnitude? Which region is more seismically active in terms of the number of earthquakes? Which region is more seismically active in terms of earthquake magnitude?

    3.6 Compare all three frequency-magnitude curves. How often does a big earthquake (M > 7 or so) happen in the global catalog vs. in the two regional catalogs? Why is this?

    Submitting your work

    Now your worksheet should contain the map and answers from Part 1, the three plots and answers from Part 2, and the map, plot, and answers from Part 3. And the green grass grows all around and the green grass grows all around. Haha, but seriously, save your file in the following format:
    L4_catalog_AccessAccountID_LastName.doc (or .yourExtension)
    For example, Cardinal pitcher Michael Wacha's file would be named "L4_catalog_mjw52_wacha.doc"

    Create one document that contains:

    1. Part 1 NMSZ map
    2. Part 1 answers to follow-up questions
    3. Part 2 three frequency-magnitude diagrams
    4. Part 2 answers to the follow-up questions
    5. Part 3 a map for southern Cal
    6. Part 3 one frequency-magnitude diagram with the NMSZ, SCEC, and USGS catalog data on it
    7. Part 3 answers to the follow-up questions.

    Once you've finished this whole problem set, submit it to the "Earthquake catalog problem set" assignment in Canvas by the due date listed on the table on the first page of this lesson.

    Grading rubric

    I will use my general rubric for grading problem sets to grade this activity.


    This page titled 4.6: Earthquake Catalog Data Analysis Problem Set is shared under a CC BY-NC-SA 4.0 license and was authored, remixed, and/or curated by Eliza Richardson (John A. Dutton: e-Education Institute) via source content that was edited to the style and standards of the LibreTexts platform.