Halfspace: a Python package for elastic halfspace calculations

(edited 4 October 2014)

Part 1: Introduction

I've been working the past year and a half on calculating topographic and tectonic stresses. This work has relied heavily on using elastic halfspace approximations, which are quite commonly used in the crustal deformation communities, particularly by geodesists. To this end, I have …

AGU and SSA have new open access policies

Open access to scientific publications and datasets is very important for both scientific progress and for the incorporation of scientific knowledge into the public sphere. The reasons for this are manifold, and have been enumerated so vociferously on various blogs, email lists, and internet fora that one fears the spittle …

Low-angle Normal Fault Likelihood IPython Notebook

This is an IPython Notebook version of a paper I have just submitted to Geophysical Research Letters. This is the complete paper, but with all the code embedded, and is executable. However, it does need a Python module called eq_stats.py as well as some of the typical Python scientific …

Evolution of fault stress through the earthquake cycle: a Mohr circle perspective

Topographic Stresses and the 2008 Wenchuan, China earthquake -- AGU 2013 Presentation

For the AGU Fall Meeting 2013, we (myself, Eric Hetland, Guohong Zhang) presented a talk on stresses resulting from high topography and how those stresses can influence faulting. The presentation focused on the 2008 Wenchuan, China earthquake. Here is the slideshow:

The influence of topographic stresses on faulting, emphasizing the …

Some tips for new scientific programmers

This post is slightly modified from an email that I wrote to some friends who were looking for tips on getting started with scientific computing with MATLAB. I'm not the world's expert on the subject, and I traded MATLAB for Python a couple years back, but the advice should be …

A new open-source GIS file for active(?) low-angle normal faults

TL;DR: I made a GIS file (SHP, KML, GeoJSON) of potentially-active low-angle normal faults. You can get it at github here, and contribute if you'd like.

A MathJax sample

This post is a quick mathjax test.

$$ P(A | \sim O) = \frac{P(\sim O|A) P(A)} {P(\sim O|A)P(A) + P(\sim O|\sim A)P(\sim A)} $$

That was easy.

Active Tectonics of the Andes (ATA) active fault database now on GitHub

This post originally appeared on my old blog.

In order to facilitate collaboration in development of the Active Tectonics of the Andes, a database of active faults in the (northern, for now) Andes, we have hosted it on GitHub. The development version is here, although the ‘release’ versions will continue …