Professor Poirier | D760 | Spring 2024

OpenLab assignment: linear algebra for earthquakes!

Comment due Sunday, April 14

After last Friday’s magnitude 4.8 earthquake centered in New Jersey, I became curious about how linear algebra is applied in earthquake science. Most of what I found online was not quite relevant for us: the algebra often referred to is not linear (for example, when using triangulation to locate an earthquake’s epicenter, a system of quadratic equations is used… even ChatGPT got confused about this when I asked it for help).

However, I did find one topic that appears to use linear algebra: Geiger’s method for locating the hypocenter of an earthquake. We commonly refer to the epicenter of an earthquake as its location, but the epicenter is a point on the surface of the earth; an earthquake actually originates at some depth below the epicenter at its hypocenter. So, while the epicenter requires only two spacial coordinates, the hypocenter requires three spacial coordinates.

I am not an expert, but I have found a few resources that should help us understand Geiger’s method for locating the hypocenter of an earthquake:

  • This assignment from an earth sciences class at the Saint Louis University Earthquake Center (this looks like a draft, watch out for typos!),
  • This article (if you have seen partial derivatives in your differential equations or Calculus II class, this might make more sense to you),
  • This article (the section on Geiger’s method looks similar to the section on Geiger’s method in the previous article).

For this assignment, you do not have to read any of these resources thoroughly (unless you want to). Read just the section on Geiger’s method in one of these resources (or another one if you find one) and try to understand:

  1. what the variables represent,
  2. where the equations in the linear system and/or matrices come from,
  3. what the solutions of the linear system and/or matrix equations represent,
  4. the steps in the procedure (this is an iterative method, which means it’s probably applied more than once).

By the way, some of these resources mention finding the “least squares solution” of a matrix equation; you can read more about that procedure here.

Don’t worry about understanding everything completely or getting everything right. Comment on this post with your ideas about items 1-4 above. If you’re not confident about one or more of them, no problem, you can just say so. Read your classmates’ comments and comment if you agree or disagree with any of their ideas. This should be a conversation among you and your classmates to try to understand this application of linear algebra, even if that understanding is just at a surface level at first.

7 Comments

  1. Haris Rao

    1. Variables representation: In Geiger’s method, the variables represent the spatial coordinates (X, Y, Z) and the origin time (T) of an earthquake. Additionally, they include the coordinates of the seismic stations and the arrival times of seismic waves at these stations.
    2. Origin of equations: The equations in Geiger’s method come from the travel times of seismic waves from the earthquake to each station. The method uses the differences between observed and predicted arrival times to update the estimated earthquake location.
    3. Solutions representation: The solutions to the system of linear equations represent the adjustments to the initial estimates of the earthquake’s location and time, aiming to minimize the discrepancies between the observed and predicted seismic data.
    4. Steps in the procedure: The steps include an initial guess of the earthquake parameters, calculation of predicted arrival times, comparison with observed times to get residuals, forming a matrix equation involving derivatives of travel times, solving for adjustments, updating the parameters, and iterating this process until convergence is achieved​.
  2. Tsz Yan, Lam

    Geiger’s method works by solving a system of equations based on data from seismometers. It estimates the three-dimensional coordinates of the earthquake’s origin point and the origin time. The linear system comes from the equations modeling the relationship between the observed arrival times of seismic waves and the distances from seismometers to the hypocenter. Solving this system helps determine the location and time of the earthquake. It uses iteration to refine the estimation with each cycle.

  3. Roza

    Geiner’s method is used to find the source of the earthquake. As earthquakes has a starting point and it spreads linearly around the radius, we can use linear algebra to find the hypocenter and arrival time. Using Taylor polynomials we can write the arrival time function from hypocenter. To guess the best arrival time we write iteration of the function that self corrects after each iteration. We can put the iteration in a matrix and solve the equation.

  4. Wellington D Verduga

    1.- Geiger’s Method variables represent the coordinates of the earthquake’s hypocenter by spatial coorditantes (X, Y, Z) and the origin time (T). On top of that you have coordinates for stations(x_j, y_j, o), arrival time(t_j), travel time (TT).

    2.- The equations come from the study of the physics of wave propagation, They represent the travel time of the seismic waves from the hypocenter to seismic stations.

    3.- The solutions represent estimates of the earthquake’s hypocenter and origin time by minimizing the difference between observed travel times and the ones being calculated.

    4.- You take an initial guess for hypocenter and the origin time. You have to then calculate the travel time, which would adjust the hypocenter coordinates and origin time with the goal of minimizing the difference between the observed and calculated travel times. This is performed iteratively until the solution converges.

  5. Joseph Valentin

    The Geiger’s method can be used to find the center of an earthquake. The equations and matrices come from the seismic waves and are used to determine the location until it can converge to the epicenter. It’s a lot of calculations and requires the concepts that linear algebra teaches. The variables likely represent the coordinates of the earthquake’s epicenter

  6. Gurpreet

    1. Variables Representation: In the context of Geiger’s method, variables typically represent the spatial coordinates of the hypocenter (x, y, z) and the origin time of the earthquake
    • Gurpreet

      1. Equations and Matrices Origin: The core of Geiger’s method involves setting up a system of equations based on the travel times of seismic waves from the hypocenter to various seismograph stations. The differences between the observed travel times and the calculated travel times (based on an initial estimate of the hypocenter and origin time) are used to set up these equations. The goal is to minimize these differences, which is often framed as a least squares problem—a fundamental problem in linear algebra.
      2. Solutions Interpretation: Solving the linear system or matrix equations results in adjustments to the initial estimates of the hypocenter’s coordinates and the origin time. These solutions tell us how much we need to adjust our initial guesses to more accurately pinpoint the hypocenter.
      3. Iterative Procedure: Geiger’s method is iterative, meaning it repeats its process to refine the hypocenter’s coordinates. Initially, an estimate is made, and then adjustments are calculated using linear algebra techniques. The new estimates (initial guesses plus adjustments) are used in the next iteration, and this process continues until the changes become negligibly small, indicating that the solution has converged to an accurate location.

Leave a Reply

Your email address will not be published. Required fields are marked *