Statistics – Computation
Scientific paper
Jul 2008
adsabs.harvard.edu/cgi-bin/nph-data_query?bibcode=2008geoji.174..295r&link_type=abstract
Geophysical Journal International, Volume 174, Issue 1, pp. 295-308.
Statistics
Computation
5
Inverse Theory, Seismic Tomography, Computational Seismology, Australia
Scientific paper
A new technique designed for generating multiple solutions to seismic tomography problems using gradient based inversion is presented. The basic principle is to exploit information gained from previous solutions to help drive the search for new models. This is achieved by adding a feedback or evolution term to the objective function that creates a local maximum at each point in parameter space occupied by the previously computed models. The advantage of this approach is that it only needs to produce a relatively small ensemble of solutions, since each model will substantially differ from all others to the extent permitted by the data. Common features present across the ensemble are, therefore, likely to be well constrained. A synthetic test using surface wave traveltimes and a highly irregular distribution of sources and receivers shows that a range of different velocity models are produced by the new technique. These models tend to be similar in regions of good path coverage, but can differ substantially elsewhere. A simple measure of the variation across the solution ensemble, given by one standard deviation of the velocity at each point, accurately reflects the robustness of the average solution model. Comparison with a standard bootstrap inversion method unequivocally shows that the new approach is superior in the presence of inhomogeneous data coverage that gives rise to under or mixed-determined inverse problems. Estimates of posterior covariance from linear theory correlate more closely with the dynamic objective function results, but require accurate knowledge of a priori model uncertainty. Application of the new method to traveltimes derived from long-term cross-correlations of ambient noise contained in passive seismic data recorded in the Australian region demonstrates its effectiveness in practice, with results well corroborated by prior information. The dynamic objective function scheme has several drawbacks, including a somewhat arbitrary choice for the shape of the evolution term, and no guarantee of a thorough exploration of parameter space. On the other hand, it is tolerant of non-linearity in the inverse problem, is relatively straightforward to implement, and appears to work well in practice. For many applications, it may be a useful addition to the suite of synthetic resolution tests that are commonly used.
Rawlinson Nicholas
Sambridge Malcolm
Saygin E.
No associations
LandOfFree
A dynamic objective function technique for generating multiple solution models in seismic tomography does not yet have a rating. At this time, there are no reviews or comments for this scientific paper.
If you have personal experience with A dynamic objective function technique for generating multiple solution models in seismic tomography, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and A dynamic objective function technique for generating multiple solution models in seismic tomography will most certainly appreciate the feedback.
Profile ID: LFWR-SCP-O-1317991