| Title | Time-Dependent Seismic Tomography of Geothermal Systems |
|---|---|
| Authors | Julian, B.R.; Foulger, G.R. |
| Year | 2009 |
| Conference | Stanford Geothermal Workshop |
| Keywords | microearthquake; moment tensor; Enhanced Geothermal Systems; Coso Geothermal Field |
| Abstract | Temporal changes in seismic wave speeds in the Earth’s crust have been measured at several geothermal areas, notably The Geysers in California, in studies that used three-dimensional seismic tomography. These studies used conventional tomography methods to invert seismic-wave arrival time data sets for different epochs independently and assumed that any differences in the derived structures reflect real temporal variations. Such an assumption is dangerous because the results of repeated tomography experiments would differ even if the structure did not change, simply because of variation in the seismic ray distribution caused by the natural variation in earthquake locations. This problem can be severe when changes in the seismicity distribution are systematic, as, for example, when many data come from earthquake swarms. The sudden change in the ray distribution can produce artifacts that mimic changes in the seismic wave speeds at the time of a swarm. Even if the source locations did not change (if only explosion data were used, for example), derived structures would inevitably differ because of observational errors. In order to determine what changes are real and what are artifacts, we have developed a new tomography method, tomo4d, that inverts multiple data sets simultaneously, imposing constraints to minimize the differences between the models for different epochs. This problem is similar to that of seeking models similar to some a priori initial assumption, and a method similar to “damped least squares” can solve it. We present an algorithm for performing this computation efficiently. Using our program, inverting multiple epochs simultaneously is comparable in difficulty to inverting them independently. We illustrate the program’s performance using synthetic arrival times. |