My research focuses on the interaction between earthquakes and groundwater. I use data from water wells, GPS stations and seismology along with computer models to better understand how stress changes, deformation, and pore fluid pressure affect each other.

Open Source Codes (updated 01/2017)

The codes I developed for my PhD are free to use and available here. The models are written in Python and use the FENICS computing platform. They solve the coupled set of poroelastic equations on both 2D and 3D subduction zone meshes. The 2D models can simulate synthetic earthquakes, slow slip events, subduction zone cycles and topographic driven flow. The 3D models can be used to simulate synthetic or data-driven earthquakes and slow slip events. Below is a short manual to help you download and run the fenics/docker programs as well as the models. Also available here