Did the 1999 Earthquake Swarm on Gakkel Ridge Open a Volcanic Conduit? a Detailed Teleseismic Data Analysis Carsten Riedel, Vera Schlindwein
Total Page:16
File Type:pdf, Size:1020Kb
Did the 1999 earthquake swarm on Gakkel Ridge open a volcanic conduit? A detailed teleseismic data analysis Carsten Riedel, Vera Schlindwein To cite this version: Carsten Riedel, Vera Schlindwein. Did the 1999 earthquake swarm on Gakkel Ridge open a volcanic conduit? A detailed teleseismic data analysis. Journal of Seismology, Springer Verlag, 2009, 14 (3), pp.505-522. 10.1007/s10950-009-9179-6. hal-00537995 HAL Id: hal-00537995 https://hal.archives-ouvertes.fr/hal-00537995 Submitted on 20 Nov 2010 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. J Seismol (2010) 14:505–522 DOI 10.1007/s10950-009-9179-6 ORIGINAL ARTICLE Did the 1999 earthquake swarm on Gakkel Ridge open a volcanic conduit? A detailed teleseismic data analysis Carsten Riedel · Vera Schlindwein Received: 19 February 2009 / Accepted: 2 November 2009 / Published online: 20 November 2009 © Springer Science + Business Media B.V. 2009 Abstract In 1999, a seismic swarm of 237 tele- ascending toward the potential conduit during the seismically recorded events marked a submarine beginning of April 1999, indicating an opening of eruption along the Arctic Gakkel Ridge, later on the vent. also analyzed by sonar, bathymetric, hydrother- mal, and local seismic studies. We relocated the Keywords Gakkel Ridge · HYPOSAT · swarm with the global location algorithm HY- Relocation · Station residuals · Spatiotemporal POSAT and analyzed the waveforms of the sta- organization · 1999 swarm · Polar research · tions closest to the events by cross-correlation. Double difference · Cross-correlation · Magma We find event locations scattered around 85◦35 ascent · Oden volcano · Ultraslow-spreading ridge Nand85◦ E at the southern rift wall and inside the rift valley of the Gakkel Ridge. Waveforms of three highly correlating events indicate a vol- 1 Introduction umetric moment tensor component and highly precise referenced double-difference arrival times The Arctic Gakkel Ridge (Fig. 1) is one of the lead us to believe that they occur at the same slowest-spreading zones of the Earth’s lithosphere geographical position and mark the conduit lo- (Dick et al. 2003) with a full spreading rate lower cated further southeast close to a chain of recently than 1.2 cm/year. Along the rift axis, volcanic and imaged volcanic cones. This result is supported amagmatic units alternate (Michael et al. 2003). by station residual anomalies in the direction of Details of the active rifting processes are still the potential conduit. Seismicity is focused at the unknown due to the remote location of Gakkel crust–mantle boundary at 16–20 km depth, but Ridge and its cover with Arctic sea ice which only few research vessels can tackle. However, in 1999, a sequence of earthquakes of body wave C. Riedel (B) · V. Schlindwein magnitude 3.3 (International Seismological Cen- Alfred-Wegener-Institut für Polar- und tre [ISC]) and higher was recorded worldwide Meeresforschung, Am alten Hafen 26, 27568, Bremerhaven, Germany (Müller and Jokat 2000; Tolstoy et al. 2001), which e-mail: [email protected] is so far the largest teleseismically recorded earth- quake swarm at any mid-ocean ridge (MOR) vol- Present Address: canic system and occurred in the area surrounding C. Riedel ◦ ◦ KGS Rastede, WIlhelmstr. 5, 85 40 N, 85 E (Fig. 1). Global centroid moment 26180 Rastede, Germany tensor solutions with large isotropic components 506 J Seismol (2010) 14:505–522 270˚ 240˚ 300˚ 210˚ 330˚ 80˚ 85 m below sea surface ˚48' 82˚ 84˚ 86˚ 0˚ 88˚ 90˚ 180˚ 85˚48' 30˚ 150˚ 85˚36' 60 ˚ ˚ 120 90˚ 85˚36' 70 60 50 85˚24' 40 85˚24' 30 km 20 0 5 10 10 85˚12' 80˚ 0 82˚ 3 4 5 6 84˚ 86˚ 85˚12' 88˚ 90˚ Fig. 1 The ISC locations of the 237 events mentioned size of the circles in the location graph is proportional to in the text have been plotted as dark gray circles on a the square values of magnitude (scale see below magnitude bathymetric map of Gakkel Ridge in the survey area on distribution) to ease the view of high-magnitude events. the left side of the figure. A scale for the bathymetry where A small white triangle (volcano) marks the position of the black is related to topographic lows and white is related to crater of Oden volcano which was discovered by Sohn topographic highs is displayed at the top. On the right side et al. (2008). All figures including this one have been the position of the swarm area is marked as a white box plotted using Generic Mapping Tools (Wessel and Smith on a plot of the high latitudinal areas of the northern hemi- 1998) sphere. Below, a histogram of magnitudes is presented. The and high b values indicated a typical volcanic (Schlindwein et al. 2005; Schlindwein and Riedel character (Müller and Jokat 2000) and sonar im- 2009) even 2 years after the major earthquake ages from the USS Hawkbill only 4 months after swarm. In 2007, the AGAVE cruise revealed thick the main activity showed areas of high backscatter blankets of volcaniclastic sediments in the area interpreted as fresh lava flows (Edwards et al. which contain Limu o Pele, i.e., very thin glassy 2001). crystallized rims of gas bubbles in magma and, During the joint AMORE 2001 cruise on the thus, clear evidence for considerable explosive USCGC Healy and RV Polarstern rock sampling, volcanism at water depths of about 4,000 m and bathymetric surveys and hydrothermal studies dis- a chain of small volcanic peaks with Oden volcano covered that the seismically active area repre- (white triangle in Fig. 1) at its easternmost loca- sents a major volcanic complex (Michael et al. tion (Sohn et al. 2008). 2003). A giant hydrothermal plume in the wa- While MOR volcanism is generally considered ter column indicated sustained volcanic activity to be effusive, the relevance of explosive erup- (Edmonds et al. 2003; Baker et al. 2004). A local tions at MORs has only recently been recognized seismic experiment with seismic arrays on ice floes (Clague et al. 2008). Understanding the seismo- recorded a swarm of sea bottom explosive signals volcanic episode at the 85◦ E volcanic complex J Seismol (2010) 14:505–522 507 at Gakkel Ridge may, therefore, yield valuable ability of data for some of the decisive events from new information on deep submarine eruption stations farther east is scarce. Thus, we have rela- processes. In this light, we intensively analyze tively huge azimuthal gaps for waveform analysis the teleseismically recorded earthquake swarm and did not try a reanalysis of moment tensor of 1999, relocating it with an advanced classical solutions. location tool, HYPOSAT, using various veloc- ity models and exploring new ways of describ- ing the dataset in terms of its spatiotemporal 3 Methods and results organization. 3.1 HYPOSAT locations of the swarm events 2 Dataset After intercomparison with several other loca- tion tools such as HYPO71 (Lee and Lahr 1975), We used all seismic events between January and HYPOCENTER (Lienert and Havskov 1995), September 1999 registered in the ISC online bul- and ISCloc (Luckett and Storchak 2001), we used letins (we downloaded them on the 6 September HYPOSAT (Schweitzer 2001) as one of the ad- 2007 from http://www.isc.ac.uk)inanareafrom vanced classical global location tools for reloca- 85◦ Nto87◦ Nand70◦ Eto95◦ E. This dataset tion of the swarm: is not completely identical to previous studies – It is classical in the way it uses the Geiger (Müller and Jokat 2000; Tolstoy et al. 2001), which method for locating earthquakes. used a preliminary dataset available at that time. – It allows the selection of confidence levels for Our dataset contains 237 events in contrast to 252 error ellipses that we chose as 95%, the error events of Tolstoy et al. (2001) with body wave is handled as rms values as well as distance magnitudes (m ) between 3.3 and 5 (Fig. 1). b errors in all spatial directions. The original dataset was located using ISCloc – It documents the search for location solu- in the ISC standard procedures as outlined in tions far away from solutions already found Storchak (2006) using the Jeffreys–Bullen travel in previous iterations and thus documents the time tables and a complicated weighting scheme robustness of the hypocenter solution (better that delivers unrealistically small errors. The root than HYPOCENTER). mean square (rms) travel time residual for the – It uses an arbitrary regional velocity model original dataset is ∼1.14 s with only four of these and can also use a global velocity model events not fixed in depth. in farther distances from the source location For further analysis, we also inspected pub- (better than ISCloc). In our case, we used licly available waveforms of the nearest stations the IASPEI-91-model (IASP91; Kennett and that recorded the events (see Fig. 2), taking Engdahl 1991; Kennett et al. 1995)orthe care to realize an equiazimuthal coverage. This preliminary reference earth model (PREM; set of waveforms includes data from the North Dziewonski and Anderson 1981) which has Atlantic sector (Kings Bay [KBS], Spitsbergen ar- been shown to work better in marine en- ray [SPITS], Daneborg–Greenland [DAG], North vironments (Reves-Sohn, personal communi- Norway [ARCES], and Lovozero–European Rus- cation) as the global model.