Author: ggan076 (page 2 of 2)

Resolution test for tomography of the AVF using regional events

Following a resolution test of tomography of the AVF using global teleseismic sources, we perform a separate test considering the local events with Nick Rawlinson’s FMTOMO package.

The events considered have paths that involve direct and up-going waves (little p and s, as opposed to P and S in the teleseismic tomography).

681 hypocentres and seismic stations in and around the AVF.

A depth slice at 90km resolves the checker-board pattern under the stations. Smearing is apparent in the latitudinal and longitudinal slices. There isn’t much difference in terms of the results between p- and s-wave inversion.

Using a finer checker-board

In order to resolve smaller scale structure in the subsurface,  a finer scale checker-board was tested.


The first 60 km in depth seem resolved in this test. However, deeper inversion results mostly add up to the effect of smearing. Unlike the S-wave results for teleseismic tomography, s-waves seem to benefit our regional tomography results.


  • Features on the order of  ~50km, are possibly recoverable up to 150 km deep, 174.2 to 176.4 longitudinally, and -36 to -37.8 latitudinally.
  • Features half this size may be retreived up to 60 km down, 174.2 to 177.2 longitudinally, and -36.6 to -37.4 latitudinally.
  • S wave tomography shows a better lateral resolution for the periphery of the stations. Most prominently with the finer scale checker-board.

Resolution test for tomography of the AVF using global events

The Auckland Volcanic Field (AVF) is intraplate and active. GEONET seismic stations are used to monitor the field.

Stations distribution. Distinguished between those enclosing the AVF and those in the outer ring.

We propose to perform P- (and possibly S-) wave tomography to infer the structure beneath the AVF, with Nick Rawlinson’s FMTT. Our first test is to establish the resolution that we can expect from this exercise. Below, we plot the hypocentres of the earthquakes used in this test:

We compute synthetic seismic P-wave arrival times for the AVF seismic stations, as if the earth structure under the AVF has a checker-board pattern (see images on the left side of the gallery below). For this test, this will be our “true” model of the Earth. And so the objective of FMTT is to invert these synthetic travel times to recover the “true” model checker-board pattern. Areas where the inverted results show a checker-board pattern suggest a good resolution.

Some of the sources need to be removed (17 out of 649 events) in order to run S wave tomography as the code could not determine the traveltimes of these 17 events from their respective sources to the base of the 3D model. Even still, S wave tomography shows a more prominent recovery overall. A major difference can be seen in the EW slice, where a band of positive anomaly cannot be recovered in the P wave equivalent.

Lateral resolution is better in NS rather than in EW orientation in particular at shallower depth. This is expected because the stations are more widely spread in this direction. In both though, resolution of the recovered checker-board is smeared, especially prominent in the lateral extent of the plot.

Directly below the AVF, we can see a pair of checker-board pattern being reflected. The EW slice at -37.0 shows this the best with minimal smearing. The pattern translates horizontally as well as shown in the depth slice. Also, some checker-board pattern persist even to the depth of >200 km.


  • 50-km features seem to be resolved for the very shallow (20km?) and very deep (200km+) parts under the seismic array.
  • S wave tomography shows a slighlty higher resolution, especially deeper  and at greater lateral extent.

If we consider closer events we can obtain better resolution at shallower depth. In order to do this we use Nick Rawlinson’s FMTOMO, which results can be seen in this post.