Merging active and passive data sets in traveltime tomography: The case study of Campi Flegrei caldera (Southern Italy)

We propose a strategy for merging both active and passive data sets in linearized tomographic inversion. We illustrate this in the reconstruction of 3D images of a complex volcanic structure, the Campi Flegrei caldera, located in the vicinity of the city of Naples, southern Italy. The caldera is occasionally the site of significant unrests characterized by large ground uplifts and seismicity. The P and S velocity models of the caldera structure are obtained by a tomographic inversion based on travel times recorded during two distinct experiments. The first data set is composed of 606 earthquakes recorded in 1984 and the second set is composed of recordings for 1528 shots produced during the SERAPIS experiment in 2001. The tomographic inversion is performed using an improved method based on an accurate finite‐difference traveltime computation and a simultaneous inversion of both velocity models and earthquake locations. In order to determine the adequate inversion parameters and relative data weighting factors, we perform massive synthetic simulations allowing one to merge the two types of data optimally. The proper merging provides high resolution velocity models, which allow one to reliably retrieve velocity anomalies over a large part of the tomography area. The obtained images confirm the presence of a high P velocity ring in the southern part of the bay of Pozzuoli and extends its trace inland as compared to previous results. This annular anomaly represents the buried trace of the rim of the Campi Flegrei caldera. Its shape at 1.5 km depth is in good agreement with the location of hydrothermalized lava inferred by gravimetric data modelling. The Vp/Vs model confirms the presence of two characteristic features. At about 1 km depth a very high Vp/Vs anomaly is observed below the town of Pozzuoli and is interpreted as due to the presence of rocks that contain fluids in the liquid phase. A low Vp/Vs body extending at about 3–4 km depth below a large part of the caldera is interpreted as the top of formations that are enriched in gas under supercritical conditions.     

Pyroclastic deposits as a guide for reconstructing the multi-stage evolution of the Somma-Vesuvius Caldera

 The evolution of the Somma-Vesuvius caldera has been reconstructed based on geomorphic observations, detailed stratigraphic studies, and the distribution and facies variations of pyroclastic and epiclastic deposits produced by the past 20,000 years of volcanic activity. The present caldera is a multicyclic, nested structure related to the emptying of large, shallow reservoirs during Plinian eruptions. The caldera cuts a stratovolcano whose original summit was at 1600–1900 m elevation, approximately 500 m north of the present crater. Four caldera-forming events have been recognized, each occurring during major Plinian eruptions (18,300 BP "Pomici di Base", 8000 BP "Mercato Pumice", 3400 BP "Avellino Pumice" and AD 79 "Pompeii Pumice"). The timing of each caldera collapse is defined by peculiar "collapse-marking" deposits, characterized by large amounts of lithic clasts from the outer margins of the magma chamber and its apophysis as well as from the shallow volcanic and sedimentary units. In proximal sites the deposits consist of coarse breccias resulting from emplacement of either dense pyroclastic flows (Pomici di Base and Pompeii eruptions) or fall layers (Avellino eruption). During each caldera collapse, the destabilization of the shallow magmatic system induced decompression of hydrothermal–magmatic and hydrothermal fluids hosted in the wall rocks. This process, and the magma–ground water interaction triggered by the fracturing of the thick Mesozoic carbonate basement hosting the aquifer system, strongly enhanced the explosivity of the eruptions. Received: 24 November 1997 / Accepted: 23 March 1999     

Seismic microzonation of Marsa Alam, Egypt using inversion HVSR of microtremor observations

An attempt was made to estimate seismic microzonation in Marsa Alam city, Egypt based on the analysis of seismic microtremor observations. Observations were carried out at 140 sites in the study region. Analysis and processing of microtremor were divided into two steps; the first one is to measure the horizontal-to-vertical spectral ratios (HVSR) for each site and picked predominant frequency and its amplitude for each site. The second step is to measure the average shear wave velocity in the upper 30 m (Vs30) of subsoil using inversion of HVSR technique. The results show well matching of theoretical HVSR and observed HVSR for body waves in all sites. The Vs30 parameter was used to classify subsoil into classes of soil properties converted to the slandered European soil code (Eurocode-8 (2002)) as follows; Vs < 180(class D), 180 ≤ Vs < 360 (class C), 360 ≤ Vs < 800 (class B) and Vs ≥ 800 m/s (class A). Our study exhibits that most of Marsa Alam city were covered by B and C classes with small portions of class D close to the shoreline and class A at the northeastern part of the region. The developed classification soil map of the study area was correlated with the distribution of the predominant frequency in view of the surface geology and given a good matching. The results of this study will be useful for planning the Marsa Alam region to be the future tourist dream for Egypt. The method used in assessment of seismic microzonation in Marsa Alam city could be the fast and inexpensive technique to measure the Vs30 based on the HVSR of microtremor and would be applied in many other areas in Egypt.     

Structure and evolution of the Bay of Pozzuoli (Italy) using marine seismic reflection data: implications for collapse of the Campi Flegrei caldera

Digital marine seismic reflection data acquired in 1973 in the Bay of Pozzuoli, and recently reprocessed, were used to study the volcanological evolution of the marine sector of Campi Flegrei Caldera during the last 37 ka. In order to gain more information, interpretation also involved estimation of the "pseudo-velocity" and the "pseudo-density" from the resistivity logs of two onshore deep exploration wells. The main results are: (1) discovery of ancient pre-18 ka and post-37 ka submarine and mainly effusive volcanic activity, along coeval emission centers located at the edges of Campi Flegrei Caldera; (2) confirmation that the caldera collapse in the marine sector of Campi Flegrei seems strongly controlled by regional NE–SW and NW–SE structural discontinuities; (3) the finding of at least two episodes of collapse in the bay; and (4) identification of a post-18 ka volcanic deflation phase that has caused about 150–200 m of subsidence in the central sector of the Bay of Pozzuoli in the last 18 ka.Editorial responsibilty: T. Druitt     

Gravity and aeromagnetic constraints on the structure of the Woods Mountains volcanic center, southeastern California

 The Woods Mountain volcanic center is a well-exposed, mildly alkaline volcanic center that formed during the Miocene in southeastern California. Detailed geologic mapping and geochemical studies have distinguished three major volcanic phases: precaldera, caldera forming, and postcaldera. Geologic mapping indicates that caldera formation occurred incrementally during eruptions of three large ignimbrites and continued into a period of voluminous intracaldera lava-flow eruptions. Rhyolitic ignimbrites and lava flows within the caldera are associated with large amplitude, circular gravity, and magnetic minima that are among the most prominent gravity and magnetic anomalies in southeastern California. Analysis of a Bouguer gravity anomaly map, reduced-to-the-pole magnetic intensity map, and three-dimensional gravity and magnetic models indicates that there is a single, funnel- to bowl-shaped caldera approximately 4 km thick and approximately 10 km wide at the surface. This model is consistent with other siliceous, pyroclastic-filled calderas on continental crust, except that most siliceous volcanic centers associated with more than one eruption are characterized by more than one caldera. Received: 20 December 1997 / Accepted: 15 October 1998     

Vp和Vp/Vs联合反演:在日本东北地区深部结构成像(英文)

本文提出了一种新的反演方法:通过采用纵、横波走时数据对(从相同的震源产生的P和S波被同一台站记录)来联合反演纵波速度(Vp)和纵、横波速度比(Vp/Vs),然后单独反演横波速度Vs,在反演过程中同时对地震参数进行定位。该方法不需要假设P和S波的射线路径一致,它是沿着P和S波射线路径计算相对慢度扰动值。该方法直接把Vp/Vs作为一个模型参数,由此能获得比采用从独立反演获得的Vp和Vs计算出Vp/Vs的方法更精确的速度比值。该新方法被应用到反演日本东北地区的壳幔速度及波速比结构的研究中,获得了较好的效果。反演结果表明,在日本东北地区,太平洋俯冲板块为一高Vp,高Vs和低Vp/Vs异常区,而在活火山下方的浅部地幔楔以及背弧深部地区为低Vp,低VS和高Vp/VS异常。虽然这些特征在前人的研究中已经报道过,但与前人的研究结果相比,本次研究所获得的Vp/Vs的空间分布具有较小的分散性,同时,它的分布特征能较好的与地震波速度结构相吻合。     

火山岩地震屏蔽层的转换波叠前时间偏移成像

在反射地震转换波资料处理中,准确求取共转换点一直是一个难题,采用叠前时间偏移技术能避免共转换点道集的抽取,而且能够使转换波归位到真正的反射点上,实现准确成像.本文针对火山岩地震屏蔽层的转换波成像问题,通过对转换波共近似转换点道集进行速度分析,建立了转换波叠前时间偏移的初始速度场,通过速度扫描和纵、横波速度比值扫描确定最佳的偏移速度场和纵、横波速度比值,实现了在火山岩高速层覆盖区域的转换波偏移成像.实际资料的成像结果表明,本文采用的近似转换点计算以及转换波叠前时间偏移方法是有效的.     

Evaluating fracture patterns within a resurgent caldera: Campi Flegrei, Italy

Understanding deformation of active calderas allows their dynamics to be defined and their hazard mitigated. The Campi Flegrei resurgent caldera (Italy) is one of the most active and hazardous volcanoes in the world, characterized by post-collapse resurgence, eruptions, ground deformation, and seismicity. An original structural analysis provides an overview of the main fracture zones. NW-SE and NE-SW fractures (normal or transtensive faults and extensional fractures) predominate along the rim and within the caldera, suggesting a regional control, both during and after the collapses. While the NE-SW fractures are ubiquitous in the deposits of the last ∼37 ka, NW-SE fractures predominate in the last 4.5 ka, during resurgence. The most recently (<4.5 ka) strained area lies in the caldera center (Solfatara area), where the faults, with an overall ∼ENE-WSW extension direction, appear to be associated with the bending due to resurgence. Solfatara lies immediately to the east of the most uplifted part of the caldera (Pozzuoli area), where domes form and culminate both on the long-term (resurgence, accompanied by volcanic activity) and short-term deformation (1982–1984 bradyseism, accompanied by seismic and hydrothermal activity). Similar volcano-tectonic behavior characterizes the short- and long-term uplifts, and only the intensity of the tectonic and volcanic activity varies, being related to varying amounts of uplift. Seismicity and hydrothermal manifestations occur during the bradyseisms, with moderate uplift, while surface faulting and eruptions occur during resurgence, with higher uplift. The features observed at Campi Flegrei are found at other major calderas, suggesting consistent behavior of large magmatic systems.     

Modeling of Rayleigh waves dispersion in the Sannio region (Southern Italy) from seismic active refraction data

Short period surface waves, recorded during a seismic refractionsurvey in the Sannio region (Southern Italy), have been modeled to infera shallow velocity model for the area. Based on the decrease of resolutionwith depth, due to the bias on group velocity estimates arising frominterference of the Rayleigh waves with higher modes, we carried out aprocedure of fitting, with synthetic seismograms, of selected filtered traceswith a gaussian filter, having a width at half height equal to 1 Hz and acentral frequency lying in the range [1,4] Hz. We estimated the likelihoodbetween synthetic and observed seismograms by measuring their semblance.In this way we were able to infer a more refined local velocity modelcharacterized by a high Vp and Vs vertical gradient in the sedimentarycover. Two ad hoc resolution studies, based on group velocity andamplitude data respectively, indicate that the local velocity model is a goodvelocity model also for the entire studied area. The increase in the numberof available data when using amplitude information allows us to make amore selective choice in the model parameter space (Vp and Vs of eachlayer) and to solve for the Vp/Vs ratio. The inferred Vp velocity in thehalf-space is equal to 2.8 km/s. This value is in excellent agreement withthat inferred by other authors (3 km/s) by modeling P-wave travel timevs. distance. The best-fit model furnish low Vp/Vs for the sedimentarycover so indicating a high degree of the sediment's compaction in thestudied area. The inferred shallow high-velocity gradient indicates thatthe shallow sedimentary layer in the area could trap and focus the energytraveling into it.     

Iterative tomographic analysis based on automatic refined picking

The ever-growing size of data sets for active and passive seismic imaging makes the availability of automatic procedures for rapid analysis more and more valuable. Such procedures are especially important for time-critical applications like emergency decisions or re-orienting of ongoing seismic surveys. In this paper a new, iterative scheme for 3D traveltime tomography is presented. The technique, based on a tool originally developed for earthquake data, uses cross-correlation to examine waveform similarity and to adjust arrival times on seismic sections. A preliminary set of reference arrival times is first corrected by the cross-correlation lag and then used to build an initial 3D tomographic velocity model through a standard inversion code; traveltimes calculated from this model are then taken as new reference arrivals and the process of pick adjustment is repeated. The result is a tomographic image, upgraded and refined at each iteration of the procedure. The test performed on the waveform data set recorded during the 2001 SERAPIS active seismic survey in the gulfs of Naples and Pozzuoli (Southern Italy) shows that the 3D iterative tomography scheme produces a velocity image of the structure of the Campi Flegrei caldera which is consistent with the results from previous studies, employing just a fraction of the time needed by a human analyst to identify first breaks. We believe that this technique can be effectively employed for rapid analysis of large data-sets within time-critical or time-dependent tasks and for automatic 4D tomographic investigations.     

The role of regional-scale faults in controlling a trapdoor caldera,Coromandel Peninsula,New Zealand

Regional gravity data from an eroded Miocene to Pliocene volcanic arc exposed in the Coromandel Peninsula, New Zealand, delineate a circular − 26-mGal, 15-km-diameter gravity anomaly. This anomaly, which has steep gradients on its northern and western margins but shallow gradients elsewhere, correlates with relatively young volcanic and volcaniclastic rocks within a broad topographic depression. Gravity modelling, using an exponentially decreasing density contrast with depth profile, requires very low-density rocks (ca. 2280 kg m^− 3) in the near-surface to account for the observed anomaly, giving a total depth of ca. 2.8 km for these rocks. The northern and western margins of this body dip steeply inward at 70°, whereas the southern and eastern margins have shallow inward dips (20–30°). The western margin coincides with the regional-scale Mangakino Fault, but the northern margin, recognizable only in the geophysical data (and named here the Ohinemuri Fault), is partially buried under younger volcanic rocks. We interpret these deep and steeply bounded, low-density volcanics in terms of a trapdoor caldera, faulted on its northern and western margins, with its hinge on the southern and eastern margins. Epithermal deposits are spatially associated with the Mangakino and Ohinemuri Faults, suggesting that both structures may have influenced hydrothermal fluid flow. These deposits pre-date caldera fill, indicating that caldera development followed pre-existing regional faults. These results delineate the subsurface geometry of a trapdoor caldera and highlight the role of pre-existing, regional-scale faults in controlling such caldera location and collapse.     

3-D Modelling of Campi Flegrei Ground Deformations: Role of Caldera Boundary Discontinuities

Campi Flegrei is a caldera complex located west of Naples, Italy. The last eruption occurred in 1538, although the volcano has produced unrest episodes since then, involving rapid and large ground movements (up to 2 m vertical in two years), accompanied by intense seismic activity. Surface ground displacements detected by various techniques (mainly InSAR and levelling) for the 1970 to 1996 period can be modelled by a shallow point source in an elastic half-space, however the source depth is not compatible with seismic and drill hole observations, which suggest a magma chamber just below 4 km depth. This apparent paradox has been explained by the presence of boundary fractures marking the caldera collapse. We present here the first full 3-D modelling for the unrest of 1982–1985 including the effect of caldera bordering fractures and the topography. To model the presence of topography and of the complex caldera rim discontinuities, we used a mixed boundary elements method. The a priori caldera geometry is determined initially from gravimetric modelling results and refined by inversion. The presence of the caldera discontinuities allows a fit to the 1982–1985 levelling data as good as, or better than, in the continuous half-space case, with quite a different source depth which fits the actual magma chamber position as seen from seismic waves. These results show the importance of volcanic structures, and mainly of caldera collapses, in ground deformation episodes.     

Erratum to: Lithospheric structure of NW Iran from P and S receiver functions

We computed P and S receiver functions to investigate the lithospheric structure beneath the northwest Iran and compute the Vp/Vs ratio within the crust of this seismologically active area. Our results enabled us to map the lateral variations of the Moho as well as those of the lithosphere–asthenosphere boundary (LAB) beneath this region. We selected data from teleseismic events (M_b > 5.5, epicentral distance between 30° and 95° for P receiver functions and M_b > 5.7, epicentral distance between 60° and 85° for S receiver functions) recorded from 1995 to 2008 at 8 three-component short-period stations of Tabriz Telemetry Seismic Network. Our results obtained from P receiver functions indicate clear conversions at the Moho boundary. The Moho depth was firstly estimated from the delay time of the Moho converted phase relative to the direct P wave. Then we used the H-Vp/Vs stacking algorithm of Zhu and Kanamori to estimate the crustal thickness and Vp/Vs ratio underneath the stations with clear Moho multiples. We found an average Moho depth of 48 km, which varies between 38.5 and 53 km. The Moho boundary showed a significant deepening towards east and north. This may reveal a crustal thickening towards northeast possibly due to the collision between the Central Iran and South Caspian plates. The obtained average Vp/Vs ratio was estimated to be 1.76, which varies between 1.73 and 1.82. The crustal structure was also determined by modeling of P receiver functions. We obtained a three-layered model for the crust beneath this area. The thickness of the layers is estimated to be 6–11, 18–35, and 38–53 km, respectively. The average of the shear wave velocity was calculated to be 3.4 km/s in the crust and reaches 4.3 km/s below the Moho discontinuity. The crustal thickness values obtained from P receiver functions are in good agreement with those derived by S receiver functions. In addition, clear conversions with negative polarity were observed at ~8.7 s in S receiver functions, which could be related to the conversion at the LAB. This may show a relatively thin continental lithosphere of about 85 km implying that the lithosphere was influenced by various geodynamical reworking processes in the past.     

A re-appraisal of the stratigraphy and volcanology of the Cerro Galán volcanic system,NW Argentina

From detailed fieldwork and biotite ⁴⁰Ar/³⁹Ar dating correlated with paleomagnetic analyses of lithic clasts, we present a revision of the stratigraphy, areal extent and volume estimates of ignimbrites in the Cerro Galán volcanic complex. We find evidence for nine distinct outflow ignimbrites, including two newly identified ignimbrites in the Toconquis Group (the Pitas and Vega Ignimbrites). Toconquis Group Ignimbrites (~5.60–4.51 Ma biotite ages) have been discovered to the southwest and north of the caldera, increasing their spatial extents from previous estimates. Previously thought to be contemporaneous, we distinguish the Real Grande Ignimbrite (4.68 ± 0.07 Ma biotite age) from the Cueva Negra Ignimbrite (3.77 ± 0.08 Ma biotite age). The form and collapse processes of the Cerro Galán caldera are also reassessed. Based on re-interpretation of the margins of the caldera, we find evidence for a fault-bounded trapdoor collapse hinged along a regional N-S fault on the eastern side of the caldera and accommodated on a N-S fault on the western caldera margin. The collapsed area defines a roughly isosceles trapezoid shape elongated E-W and with maximum dimensions 27 × 16 km. The Cerro Galán Ignimbrite (CGI; 2.08 ± 0.02 Ma sanidine age) outflow sheet extends to 40 km in all directions from the inferred structural margins, with a maximum runout distance of ~80 km to the north of the caldera. New deposit volume estimates confirm an increase in eruptive volume through time, wherein the Toconquis Group Ignimbrites increase in volume from the ~10 km³ Lower Merihuaca Ignimbrite to a maximum of ~390 km³ (Dense Rock Equivalent; DRE) with the Real Grande Ignimbrite. The climactic CGI has a revised volume of ~630 km³ (DRE), approximately two thirds of the commonly quoted value.     

P and S reflection and P refraction: An integration for characterising shallow subsurface

Characterization of shallow structures was performed by using different approaches analysing both P- and S-wave seismic data with different resolution. The refraction tomography provided P and S velocity models of the first 80 m, while the reflection seismic processing gives a reasonable stacking velocity field until 300 m depth for both P- and S-wave data. So, we estimated the V_p/V_s ratio and an empirical relationship between the two velocities. We characterised the shallow layers using tomographic velocity models and the deeper layers using seismic images with different resolution. The seismic images were obtained by conventional CMP reflection seismic processing and by a novel multi-refractor imaging technique.     

The P-wave velocity structure of Deception Island,Antarctica, from two-dimensional seismic tomography

Deception Island is a volcanic island with a flooded caldera that has a complex geological setting in Bransfield Strait, Antarctica. We use P-wave arrivals recorded on land and seafloor seismometers from airgun shots within the caldera and around the island to invert for the P-wave velocity structure along two orthogonal profiles. The results show that there is a sharp increase in velocity to the north of the caldera which coincides with a regional normal fault that defines the northwestern boundary of the Bransfield Strait backarc basin. There is a low-velocity region beneath the caldera extending from the seafloor to > 4 km depth with a maximum negative anomaly of 1 km/s. Refracted arrivals are consistent with a 1.2-km-thick layer of low-velocity sediments and pyroclastites infilling the caldera. Synthetic inversions show that this layer accounts for only a small portion of the velocity anomaly, implying that there is a significant region of low velocities at greater depths. Further synthetic inversions and melt fraction calculations are consistent with, but do not require, the presence of an extensive magma chamber beneath the caldera that extends downwards from ≤ 2 km depth.     

Structural deformation and sedimentation in an active Caldera, Rabaul, Papua New Guinea

Recent seismic and tectonic activity in Rabaul Caldera, Papua New Guinea, suggests that magma is accumulating at a shallow depth beneath this partially submerged structure and that a new volcano may be developing. Changes in onshore elevation since 1971 (as much as 2 m on south Matupit Island) indicate that rapid and large-scale uplifts have occurred on the seafloor near the center of the caldera. The frequency of seismic events within the caldera has also increased during this period. Earthquake locations define an elliptical ring surrounding the center of this uplift within the caldera.A marine geophysical survey in 1982 by the U.S. Geological Survey's R/V “S.P. Lee” in Rabaul Caldera shows the development of a bulge in the seafloor near the center of the caldera. High-resolution seismic reflection profiles show that this bulge consists of two domal uplifts bounded and separated by two major north-south-trending fault zones. Deformed sediments overlie these zones; a prominent slump flanks the area of the bulge.Five major acoustic units were identified in the seismic reflection profiles: an acoustic basement and four sedimentary units consisting of irregularly layered, cross-layered, contorted, and well-layered sequences. The acoustic basement is probably composed of crystalline volcanic rocks, and the layered acoustic units are probably sediments, primarily ash deposited in different environments. The cross-layered, irregularly layered, and contorted units appear to have been deposited in a dynamic environment subjected to strong currents, seismicity, and/or mass wasting, while the well-layered units were deposited in a low-energy environment. Locally, well-layered sequences interfinger with the other sedimentary units, indicating a transitional environment that alternated between high-energy and low-energy depositional processes.A submarine channel cuts most of the acoustic units and appears to be the conduit for sediment transport out of the caldera; it occupies an older buried channel north of the caldera that is presently being exhumed. In the south, active erosion of well-layered sediments is taking place. What are believed to be several young volcanic cones also disrupt the depositional layers.We conclude that the bulge in the seafloor and the associated fault zones are a result of emplacement of magma at a shallow depth. Contorted sediment and slumps adjacent to the bulge are probably the result of uplift and seismic activity. The pattern of seismicity appears to reflect increased magma pressure at depth beneath the caldera floor. This activity may eventually lead to an eruption.     

Three-dimensional shallow velocity structure beneath Taal Volcano,Philippines

Based on its numerous historical explosive eruptions and high potential hazards to nearby population of millions, Taal Volcano is one of the most dangerous “Decade Volcanoes” in the world. To provide better investigation on local seismicity and seismic structure beneath Taal Volcano, we deployed a temporary seismic network consisting of eight stations from March 2008 to March 2010. In the preliminary data processing stage, three periods showing linear time-drifting of internal clock were clearly identified from noise-derived empirical Green’s functions. The time-drifting errors were corrected prior to further data analyses. By using VELEST, 2274 local earthquakes were manually picked and located. Two major earthquake groups are noticed, with one lying beneath the western shore of Taal Lake showing a linear feature, and the other spreading around the eastern flank of Taal Volcano Island at shallower depths. We performed seismic tomography to image the 3D structure beneath Taal Volcano using the LOTOS algorithm. Some interesting features are revealed from the tomographic results, including a solidified magma conduit below the northwestern corner of Taal Volcano Island, indicated by high Vp, Vs, and low Vp/Vs ratio, and a large potential hydrothermal reservoir beneath the center of Taal Volcano Island, suggested by low Vs and high Vp/Vs ratio. Furthermore, combining earthquake distributions and tomographic images, we suggest potential existence of a hydrothermal reservoir beneath the southwestern corner of Taal Lake, and a fluid conduit extending to the northwest. These seismic features have never been proposed in previous studies, implying that new hydrothermal activity might be formed in places away from the historical craters on Taal Volcano Island.     

Non-explosive,constructional evolution of the ice-filled caldera at Volcán Sollipulli,Chile

 A radar and gravity survey of the ice-filled caldera at Volcán Sollipulli, Chile, indicates that the intra-caldera ice has a thickness of up to 650 m in its central part and that the caldera harbours a minimum of 6 km³ of ice. Reconnaissance geological observations show that the volcano has erupted compositions ranging from olivine basalt to dacite and have identified five distinct volcanic units in the caldera walls. Pre- or syn-caldera collapse deposits (the Sharkfin pyroclastic unit) comprise a sequence which evolved from subglacial to subaerial facies. Post-caldera collapse products, which crop out along 17 of the 20 km length of the caldera wall, were erupted almost exclusively along the caldera margins in the presence of a large body of intra-caldera ice. The Alpehué crater, formed by an explosive eruption between 2960 and 2780 a. BP, in the southwest part of the caldera is shown to post date formation of the caldera. Sollipulli lacks voluminous silicic pyroclastic rocks associated with caldera formation and the collapse structure does not appear to be a consequence of a large-magnitude explosive eruption. Instead, lateral magma movement at depth resulting in emptying of the magma chamber may have generated the caldera. The radar and gravity data show that the central part of the caldera floor is flat but, within a few hundred metres of the caldera walls, the floor has a stepped topography with relatively low-density rock bodies beneath the ice in this region. This, coupled with the fact that most of the post-caldera eruptions have taken place along the caldera walls, implies that the caldera has been substantially modified by subglacial marginal eruptions. Sollipulli caldera has evolved from a collapse to a constructional feature with intra-caldera ice playing a major role. The post-caldera eruptions have resulted in an increase in height of the walls and concomitant deepening of the caldera with time. Received: 12 June 1995 / Accepted: 7 December 1995     

Correlation of welded ignimbrites on Pantelleria (Strait of Sicily) using paleomagnetism

Although the oldest volcanic rocks exposed at Pantelleria (Strait of Sicily) are older than 300 ka, most of the island is covered by the 45–50 ka Green Tuff ignimbrite, thought to be related to the Cinque Denti caldera, and younger lavas and scoria cones. Pre-50 ka rocks (predominantly rheomorphic ignimbrites) are exposed at isolated sea cliffs, and their stratigraphy and chronology are not completely resolved. Based on volcanic stratigraphy and K/Ar dating, it has been proposed that the older La Vecchia caldera is related to ignimbrite Q (114 ka), and that ignimbrites F, D, and Z (106, 94, and 79 ka, respectively) were erupted after caldera formation. We report here the paleomagnetic directions obtained from 23 sites in ignimbrite P (133 ka) and four younger ignimbrites, and from an uncorrelated (and loosely dated) welded lithic breccia thought to record a caldera-forming eruption. The paleosecular variation of the geomagnetic field recorded by ignimbrites is used as correlative tool, with an estimated time resolution in the order of 100 years. We find that ignimbrites D and Z correspond, in good agreement with recent Ar/Ar ages constraining the D/Z eruption to 87 ka. The welded lithic breccia correlates with a thinner breccia lying just below ignimbrite P at another locality, implying that collapse of the La Vecchia caldera took place at ~130–160 ka. This caldera was subsequently buried by ignimbrites P, Q, F, and D/Z. Paleomagnetic data also show that the northern caldera margin underwent a ~10° west–northwest (outwards) tilting after emplacement of ignimbrite P, possibly recording magma resurgence in the crust.