Nonlinear inversion of local seismic travel times for the simultaneous determination of the 3D-velocity structure and hypocentres - application to the seismic zone Vrancea
Article Sidebar
Vols. 1-18 (1924-1944), ISSN 0044-2801
Main Article Content
Abstract
A new method for solving the full nonlinear problem of simultaneous inversion for 3D structure and hypocentres (SSH method) has been developed. In comparison with the simplified linear inversion procedure of Aki and Lee (1976), the method proposed here exhibits several extensions and improvements:
a) The seismic forward problem is solved exactly by a newly developed 3D ray tracing method. This procedure has been conceived as a shooting method. The 3D heterogeneous medium is parametrized by rectangular blocks.
b) The full nonlinear inverse problem is solved. The method works iteratively and was conceived in terms of the Levenberg - Marquardt (LM) algorithm. In each iteration step the seismic forward problem is solved with the ray tracing procedure and a linear inversion is performed.
c) In order to improve the stability and uniqueness of the inverse solution, a priori information on the model space may be used in the inversion process. For example, a known crustal structure may be exploited for the retrieval of the deeper lithosphere.
The SSH method has been tested on various theoretical 1D and 3D models and its capabilities and limitations are discussed. The results demonstrate the existence of a tradeoff between hypocentral depths, origin times and seismic velocity, particularly when the horizontal dimensions of the recording network are small relative to the depths of the earthquakes. The SSH method has been applied to real earthquake data of the Vrancea region, Romania. About 50 crustal and intermediate depth events were used in the inversion process. First, the crustal events were inverted to infer the crustal structure of the region. This crustal information was then used as an a priori constraint in the inversion of the intermediate depth events. The 3D-velocity model found for the Vrancea region reduces the RMS value of the travel-time residuals from 0.5 to 0.3 s. The most essential feature of the model is a high-velocity anomaly of about 4%-6% beneath the eastern Carpathian foredeeps, extending from about 80 to 160 km depth. This appears to be in agreement with results of teleseismic studies of other authors and may be understood in terms of the plate tectonic concept proposed for the Vrancea region. Due to lack of spatial resolution and high standard errors, the results obtained above should, nevertheless, be taken with some care and have to be substantiated by further investigations with higher quality data.
ARK: https://n2t.net/ark:/88439/y008846
Permalink: https://geophysicsjournal.com/article/240
Article Details
References
Aki, K., Richards, P.G. (1980) Quantitative seismology. W.H. Freeman, San Francisco
Aki, K., Christoffersson, A., Huseybye, E.S. (1977) Determination of the three-dimensional seismic structure of the lithosphere. J. Geophys. Res. 82:277-296
Backus, G.E., Gilbert, F. (1967) Numerical applications of a formalism for geophysical inverse problems. Geophys. J. R. Astron. Soc. 13:247-276
Backus, G.E., Gilbert, F. (1968) The resolving power of gross Earth data. Geophys. J. R. Astron. Soc. 16:169-205
Backus, G.E., Gilbert, F. (1970) Uniqueness in the inversion of inaccurate gross Earth data. Phil. Trans. R. Soc. Lond. A 266:123-192
Beck, J.V., Arnold, K.J. (1977) Parameter estimation in science and engineering. John Wiley & Sons, New York-London
Buland, R. (1976) The mechanics of locating earthquakes, Bull. Seismol. Soc. Am. 66:173-187
Cerveny, V., Molotkov, I.A., Psencik, I. (1977) Ray method in seismology. Univ. Karlova Press, Prague
Crosson, R.S. (1976a) Crustal structure modelling of earthquake data, 1) simultaneous least squares estimation of hypocenter and velocity parameters. J. Geophys. Res. 81:3036-3046
Crosson, R.S. (1976b) Crustal structure modelling of earthquake data, 2) velocity structure of the Pudget Sound region, Washington. J. Geophys. Res. 81:3047-3054
Draper, N.R., Smith, H. (1966) Applied regression analysis. John Wiley & Sons, Inc., New York-London-Sydney
Engdahl, E.R. (1973) Relocation of intermediate depth earthquakes in the central Aleutians by seismic ray tracing. Nature Phys. Sci. 245:23-25
Engdahl, E.R., Lee, W.H.K. (1976) Relocation of local earthquakes by ray tracing. J. Geophys. Res. 81:4300-4306
Fuchs, K., Bonjer, K.P., Bock, G., Cornea, I., Radu, C., Enescu, D., Jianu, D., Nourescu, A., Merkler, G., Moldoveanu, T., Tudorache, G. (1979) The Romanian earthquake of March 4, 1977. II. Aftershocks and migration of seismic activity. Tectonophysics 53:235-247
Gill, P.E., Murray, W., Wright, M.H. (1981) Practical optimization. Academic Press, London-New York
Hawley, P.W., Zandt, G., Smith, R.B. (1981) Simultaneous inversion for hypocenters and lateral velocity variations: an iterative solution with a layered model. J. Geophys. Res. 86:7073-7083
Hoerl, H.A., Kennard, R.W. (1970) Ridge regression: biased estimation for nonorthogonal problems, Technometrics 12:55-82
Horie, A., Aki, K. (1982) Three-dimensional velocity structure beneath the Kanto district, Japan. Journ. Phys. Earth 30:255-282
Hovland, J., Husebye, E.S. (1982) Upper mantle heterogeneities beneath eastern Europe. Tectonophysics 90:137-151
Huppert, L.N., Froehlich, C. (1981) The P-velocity within the Tonga Benioff zone determined from traced rays and observations. J. Geophys. Res. 86:3771-3782
Jackson, D.D. (1972) Interpretation of inaccurate, insufficient and inconsistent data. Geophys. J. R. Astron. Soc. 28:97-109
Jackson, D.D. (1979) The use of a priori data to resolve non-uniqueness in linear inversion. Geophys. J. R. Astron. Soc. 57:137-157
Jung, P. (1983) Herdparameter und Ausbreitungseffekte krustaler und mitteltiefer Erdbeben in den rumanischen Karpaten unter besonderer Berucksichtigung der Vrancea-Region. Diplomathesis, Univ. Karlsruhe
Koch, M. (1982) Seismicity and structural investigations of the Romanian Vrancea region: evidence for azimuthal variations of P-wave velocity and Poisson's ratio. Tectonophysics 90:91-115
Koch, M. (1983a) The determination of lateral velocity inhomogeneities by linear and nonlinear inversion of teleseismic and local seismic travel times - application to the seismic zone Vrancea, Romania (in German). Ph.D. Thesis, University of Karlsruhe
Koch, M. (1983b) A theoretical and numerical study on the determination of the JD-structure of the lithosphere by linear and nonlinear inversion of teleseismic travel times. Geophys. J. R. Astron. Soc.
Lanczos, C. (1961) Linear differential operators. D. Van Nostrand Co., London
Lawson, C.L., Hanson, R.J. (1974) Solving least squares problems. Prentice Hall, Englewood Cliffs, New Jersey
Lee, W.H.K., Lahr, J.C. (1972) HYPO 71, a computer program for determining hypocenter, magnitude and first motion pattern of local earthquakes. U.S. Geological Survey, Open File Report
Levenberg, K. (1944) A method for the solution of certain nonlinear problems in least squares. Quart. Appl. Math. 2:164-168
Marquardt, D.W. (1963) An algorithm for least squares of nonlinear parameters. J. Soc. Industr. Appl. Math. 11:431-441
Marquardt, D.W. (1970) Generalized inverses, ridge regression, biased linear estimation and nonlinear estimation, Technometrics 12:591-612
McKenzie, D. (1972) Active tectonics of the Mediterranian region. Geophys. J. R. Astron. Soc. 30:109-185
Oncescu, M.C., Burlacu, V., Anghel, M., Smalbergher, V. (1984) Threedimensional P-wave velocity image under the Carpathian arc. Tectonophysics 106:305-320
Pavlis, G.L., Booker, J.R. (1980) The mixed discrete-continuous inverse problem: application to the simultaneous determination of earthquake hypocenters and velocity structure. J. Geophys. Res. 85:4801-4810
Pavlis, G.L., Booker, J.R. (1983) A study of nonlinearity in the inversion of earthquake arrival time data for velocity structure. J. Geophys. Res. 88:5047-5055
Roecker, S. (1982) The velocity structure of the Pamir-Hindu Kush region: possible evidence of subducted crust. J. Geophys. Res. 87:945-959
Sollugub, V.B. (1969) Seismic crustal studies in southeastern Europe. In: The Earth's crust and upper mantle, pp 189-195, Washington
Spencer, C., Gubbins, D. (1980) Travel-time inversion for simultaneous earthquake location and velocity structure determination in laterally varying media. Geophys. J. R. Astron. Soc. 63:95-116
Takanami, T. (1982) Three-dimensional seismic structure of the crust and upper mantle beneath the Orogenic belts in Southern Hokkaido, Japan. Journ. Phys. Earth 30:87-104
Thurber, C.H., Ellsworth, W.L. (1980) Rapid solution of ray tracing problems in heterogeneous media. Bull. Seismol. Soc. Am. 70:1137-1148
Tykhonov, A.N., Arsenine, V. (1976) Methodes de resolution de problemes mal poses. Editions MIR, Moscow
Vinnik, L.P., Lenartovich, E. (1975) Horizontal inhomogeneities in the upper mantle of the Carpathians and Caucasus. Tectonophysics 28:275-291
Wiggins, R.A. (1972) The general linear inverse problem: Implications of surface waves and free oscillations for Earth structure. Rev. Geophys. Space Phys. 10:251-285