Local earthquake tomography of shallow subduction in north Chile: A combined onshore and offshore study

Selected travel time data from the aftershock series of the great Antofagasta earthquake (M w =8.0) have been inverted simultaneously for both hypocenter locations and three-dimensional V P and V P / V S struc- ture. The data were collected with a dense 44-station seismic network including ocean bottom hydrophones. We performed a series of inversions with increasing complexity: 1-D, 2-D, and 3-D. This approach was cho- sen to account for the heterogeneous seismicity distribution and to obtain a smooth regional model in areas of low resolution. Special efforts were made to assess the solution quality including standard resolution esti- mates and tests with synthetic travel times. The subducted plate is imaged between 20 and 50 km in depth as an eastward dipping high- V P feature. High V P / V S ratios within the oceanic crust possibly indicate elevated fluid content. Underplating of material eroded close to the trench is found beneath the Mejillones Peninsula. The lower crust of the overlying plate is characterized by an average V P of 6.8-6.9 km/s and an average to low V P / V S ratio. Large areas of anomalously high V P are found in the lower crust south of the city of Antofagasta; they are interpreted as remants of magmatic intrusions. A zone of high V P / V S ratios is found within the rup- ture area of the Antofagasta main shock, just above the subducted slab. Its location within the region of high- est stress release from the main shock suggests that the main shock rupture causes the high V P / V S ratio. The high V P / V S ratio could indicate postseismic fluid migration from the subducted oceanic crust into the overly- ing lower crust.

[1]  G. Asch,et al.  Three‐dimensional models of P wave velocity and P‐to‐S velocity ratio in the southern central Andes by simultaneous inversion of local earthquake data , 1999 .

[2]  P. Giese,et al.  Structures and crustal development of the central andes between 21° and 25°S , 1988 .

[3]  S. Barrientos,et al.  Coseismic coastal uplift and coralline algae record in Northern Chile: The 1995 Antofagasta earthquake case , 1996 .

[4]  Alberto Michelini,et al.  Seismological studies at Parkfield. I. Simultaneous inversion for velocity structure and hypocenters using cubic B-splines parameterization , 1991, Bulletin of the Seismological Society of America.

[5]  P. Podvin,et al.  Finite difference computation of traveltimes in very contrasted velocity models: a massively parallel approach and its associated tools , 1991 .

[6]  H. Negishi,et al.  The 1995 Kobe earthquake: Seismic image of the source zone and its implications for the rupture nucleation , 1998 .

[7]  Jean Virieux,et al.  Ray tracing in 3-D complex isotropic media: An analysis of the problem , 1991 .

[8]  R. L. Fisher,et al.  Topography and structure of the Peru-Chile trench , 1962 .

[9]  S. Pichowiak.  Early Jurassic to Early Cretaceous Magmatism in the Coastal Cordillera and the Central Depression of North Chile , 1994 .

[10]  E. Flueh,et al.  Crustal structure along the EDGE transect beneath the Kodiak shelf off Alaska derived from OBH seismic refraction data , 1997 .

[11]  Clifford H. Thurber,et al.  A fast algorithm for two-point seismic ray tracing , 1987 .

[12]  W. Weinrebe,et al.  Subduction erosion along the North Chile margin , 1999 .

[13]  C. Thurber.  Earth structure and earthquake locations in the Coyote Lake area, central California , 1981 .

[14]  N. Christensen.  Poisson's ratio and crustal seismology , 1996 .

[15]  Walter H. F. Smith,et al.  New version of the generic mapping tools , 1995 .

[16]  R. Dietmar Müller,et al.  Digital isochrons of the world's ocean floor , 1997 .

[17]  K. Reutter,et al.  Tectonic Development of the North Chilean Andes in Relation to Plate Convergence and Magmatism Since the Jurassic , 1994 .

[18]  J. Viramonte,et al.  Variation in the Crustal Structure of the Southern Central Andes Deduced from Seismic Refraction Investigations , 1994 .

[19]  M. Sobiesiak.  Fault plane structure of the Antofagasta, Chile Earthquake of 1995 , 2000 .

[20]  Douglas R. Toomey,et al.  Tomographic inversion of local earthquake data from the Hengill‐Grensdalur Central Volcano Complex, Iceland , 1989 .

[21]  R. Rutland,et al.  Andean Orogeny and Ocean Floor Spreading , 1971, Nature.

[22]  E. Papadimitriou,et al.  3D crustal structure from local earthquake tomography around the Gulf of Arta (Ionian region, NW Greece) , 1999 .

[23]  E. Kissling,et al.  Accurate hypocentre determination in the seismogenic zone of the subducting Nazca Plate in northern Chile using a combined on-/offshore network , 1999 .

[24]  L. Rivera,et al.  The Mw = 8.0 Antofagasta (northern Chile) earthquake of 30 July 1995: A precursor to the end of the large 1877 gap , 1997, Bulletin of the Seismological Society of America.

[25]  J. Mechie,et al.  Two‐dimensional velocity models of the Nazca Plate subduction zone between 19.5°S and 25°S from wide‐angle seismic measurements during the CINCA95 project , 1999 .

[26]  Akira Hasegawa,et al.  Tomographic imaging of P and S wave velocity structure beneath northeastern Japan , 1992 .

[27]  K. Reutter,et al.  Magmatic arc tectonics in the Central Andes between 21° and 25°S , 1992 .

[28]  D. Eberhart‐Phillips,et al.  A three-dimensional image of shallow subduction: crustal structure of the Raukumara Peninsula, New Zealand , 1999 .

[29]  H. Philip,et al.  Recent crustal deformation in the Antofagasta region (northern Chile) and the subduction process , 1998 .

[30]  H. Niemeyer,et al.  Evolución tectónica cenozoica del margen continental activo de Antofagasta, norte de Chile , 1996 .

[31]  Wim Spakman,et al.  Imaging algorithms, accuracy and resolution in delay time tomography , 1988 .

[32]  C. Zelt.  Lateral velocity resolution from three-dimensional seismic refraction data , 1998 .

[33]  Ricardo Thiele,et al.  Active faulting in northern Chile: ramp stacking and lateral decoupling along a subduction plate boundary? , 1990 .

[34]  D. Eberhart‐Phillips,et al.  Three-dimensional velocity structure in northern California Coast Ranges from inversion of local earthquake arrival times , 1986 .

[35]  Harmen Bijwaard,et al.  Closing the gap between regional and global travel time tomography , 1998 .

[36]  L. Rivera,et al.  The Andean subduction zone between 22 and 25°S (northern Chile): precise geometry and state of stress , 1996 .

[37]  D. Eberhart‐Phillips,et al.  Three‐dimensional P and S velocity structure in the Coalinga Region, California , 1990 .

[38]  J. Byerlee,et al.  Magnetic and electric fields associated with changes in high pore pressure in fault zones : Application to the Loma Prieta ULF emissions , 1995 .

[39]  Florian Haslinger,et al.  Velocity structure, seismicity and seismotectonics of northwestern Greece between the Gulf of Arta and Zakynthos , 1998 .

[40]  R. Huene,et al.  OBSERVATIONS AT CONVERGENT MARGINS CONCERNING SEDIMENT SUBDUCTION, SUBDUCTION EROSION, AND THE GROWTH , 1991 .

[41]  John C. Lahr,et al.  Three‐dimensional P and S wave velocity structure of Redoubt Volcano, Alaska , 1996 .

[42]  D. Eberhart‐Phillips,et al.  Continental subduction and three‐dimensional crustal structure: The northern South Island, New Zealand , 1997 .

[43]  Frank Scherbaum,et al.  The GIANT Analysis System (Graphical Interactive Aftershock Network Toolbox) , 1998 .

[44]  Edi Kissling,et al.  Geotomography with local earthquake data , 1988 .

[45]  J. Malavieille,et al.  Coulomb theory applied to accretionary and nonaccretionary wedges: Possible causes for tectonic erosion and/or frontal accretion , 1994 .

[46]  S. Stein,et al.  Decelerating Nazca‐South America and Nazca‐Pacific Plate motions , 1999 .

[47]  L. Rivera,et al.  Determination of seismogenic interplate contact zone and crustal seismicity around Antofagasta, northern Chile using local data , 1994 .

[48]  Clifford H. Thurber,et al.  Earthquake locations and three‐dimensional crustal structure in the Coyote Lake Area, central California , 1983 .