SSBW - Python Tutorial 5: ObsPy and the 2018 Kilauea Eruption

Mike Brudzinski, Miami University-Oxford

Summary

In this interactive tutorial built for a learning management system, students 1) practice using the basic syntax of Python commands, 2) use the functionality of the ObsPy library to load, plot, and analyze data or metadata 3) use the ObsPy library to create stream objects for seismograms, and inventory objects for metadata, 4) evaluate seismograms to estimate seismicity rates, and 5) identify key timing and locations of events associated with the 2018 Kilauea eruption.

Seismology | College Upper (15-16)

Expand for more detail and links to related resources

Used this activity? Share your experiences and modifications

Context

Audience

The IRIS Seismology Skill Building Workshop (SSBW), which is a free, online, open-access, large-enrollment, 12-week summer workshop for upper level undergraduates.

Skills and concepts that students must have mastered

This assignment builds on prior experience using basic Python command syntax and variables, the seismological syntax of network, station, channel, and location codes for specifying a given recording, basic components of the ObsPy library (stream and inventory objects, Client functions), and the ability to retrieve relevant information about programming and seismology from internet research or help pages.

How the activity is situated in the course

This is assignment number 29 of 35 in the Seismology Skill Building Workshop (SSBW).

Goals

Content/concepts goals for this activity

Proficiency with Python, ObsPy, object structures, seismicity rates, and seismicity associated with eruptions.

Higher order thinking skills goals for this activity

Retrieving relevant information from digital sources to accomplish tasks, using correct programming syntax, evaluating and describing computing output in the context of science concepts.

Other skills goals for this activity

Description and Teaching Materials

Students will:
1. Identify the location of Puu Oo and its relationship to the Kilauea eruption in 2018
2. Estimate the timing of the Puu Oo collapse in multiple time formats using UTCDateTime
3. Retrieve an inventory of stations from the Hawaii Volcano Observatory network
4. Plot a map of stations and identify the stations closest to the lower Puna eruption
5. Retrieve waveforms from the closest station and make a traditional "drum plot"
6. Estimate the change in seismicity rate from the "drum plot" seismograms
7. Measure the time between Puu Oo collapse and the first detectable seismicity
8. Calculate the apparent speed of lava flowing from Puu Oo to the lower Puna eruption site

Show Student Learning Outcomes

Student Learning Outcomes (with corresponding Question #s that address them):

Retrieve relevant information about programming and seismology from internet research or help pages (Q1, 2, 3, 4, 5, 6, 7, 8, 9, 11, 13, 16, 18, 19, 26, 28)
Identify the locations of Puu Oo volcanic cone, the lower Puna eruption in 2018, and nearby seismic stations (Q1, 2, 3, 4, 14, 16, 17, 26)
Identify key timing of events associated with the Puu Oo collapse and lower Puna eruption (Q5, 6, 8, 20, 21, 22, 23, 24, 25)
Convert between different date and time formats including recognizing time zone differences, performing math on epoch times, and retrieving time from plots (Q6, 7, 10, 24, 25)
Use correct syntax associated with Python (Q7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 26)
Load appropriate modules for extending the functionality of the Python environment (Q7, 9, 26)
Use the ObsPy library and UTCDateTime functionality (Q7, 8, 10)
Identify differences between Python versions 2 and 3 (Q7)
Use variables and their attributes (Q8, 9, 10, 11, 13, 14, 15, 16, 17, 18, 19, 26)
Use the print command, including both variables and static text (Q8, 12, 15, 26)
Use a FDSN Client library and functions (get_events, get_stations, get_waveforms) to pull data or metadata from specific providers (Q9, 11, 12, 14, 15, 18)
Perform mathematical calculations using Python (Q10, 26, 27)
Specify network, station, channel, and location codes for a given recording (Q11, 14, 15, 18)
Use inventory objects and their functions (plot, get_coordinates) in Obspy (Q11, 12, 13, 14, 15, 16, 17)
Evaluate and describe the output of a Python command or program (Q12, 14, 15, 16, 17, 20, 21, 22, 23, 24)
Use stream objects and plot function in Obspy (Q18, 19, 20)
Estimate the seismicity rates associated with the 2018 Kilauea eruption (Q20, 21, 22, 23)
Evaluate seismograms (Q20, 21, 22, 23, 24)
Convert units (hours and minutes to seconds) (Q25)
Use a IRIS Client library and distaz function to perform distance and azimuth calculations (Q26)
Estimate the speed of a lava flow using the distance over time (Q27)

GIFT format file for SSBW Python Tutorial 5 -- educator-only file

Teaching Notes and Tips

This assignment was constructed in the Moodle learning management system, and has been exported in the GIFT format. More information about the syntax of this format can be found here: https://docs.moodle.org/en/GIFT_format

Assessment

This assignment is automatically graded by the learning management system. The number of questions of each type used are:
9 multiple choice questions
2 multiple answer questions
3 short answer (auto-graded) questions
14 numeric questions (geographic locations, output values, quantifying output, plot measurement, unit conversion, calculation)

References and Resources

During the Seismology Skill Building Workshop, students are provided with a virtual Linux machine that was tailored to include the software specifically needed to complete the assignments. This software is all freely available on the internet, but interested parties are encouraged to contact the instructor for access to this tailored virtual machine.

See more Search the Portal »