Geochemistry and isotopes of fluids from sulphur springs, Valles Caldera, New Mexico

Detailed geochemistry supported by geologic mapping has been used to investigate Sulphur Springs, an acid-sulfate hot spring system that issues from the western flank of the resurgent dome inside Valles Caldera. The most intense activity occurs at the intersection of faults offsetting caldera-fill deposits and post-caldera rhyolites. Three geothermal wells in the area have encountered pressures <1 MPa and temperatures of 200°C at depths of 600 to 1000 m. Hot spring and fumarole fluids may discharge at boiling temperatures with pH 1.0 and SO₄ 8000 mg/l. These conditions cause argillic alterations throughout a large area.Non-condensible gases consist of roughly 99% CO₂ with minor amounts of H₂S, H₂, and CH₄. Empirical gas geothermometry suggests a deep reservoir temperature of 215 to 280°C. Comparison of ¹³C and ¹⁸O between CaCO₃ from well cuttings and CO₂ from fumarole steam indicates a fractionation temperature between 200 and 300°C by decarbonation of hydrothermally altered Paleozoic limestone and vein calcite in the reservoir rocks. Tritium concentrations obtained from steam condensed in a mudpot and deep reservoir fluids (Baca #13, 278°C) are 2.1 and 1.0 T.U. respectively, suggesting the steam originates from a reservoir whose water is mostly >50 yrs old. Deuterium contents of fumarole steam, deep reservoir fluid, and local meteoric water are practically identical even though ¹⁸O contents range through 4‰, thus, precipitation on the resurgent dome of the caldera could recharge the hydrothermal system by slow percolation. From analysis of D and ¹⁸O values between fumarol steam and deep reservoir fluid, steam reaches the surface either (1) by vaporizing relatively shallow groundwater at 200°C or (2) by means of a two-stage boiling process through an intermediate level reservoir at roughly 200°C.Although many characteristics of known vapor-dominated geothermal systems are found at Sulphur Springs, fundamental differences exist in temperature and pressure of our postulated vapor-zone. We propose that the reservoir beneath Sulphur Springs is too small or too poorly confined to sustain a “true” vapor-dominated system and that the Sulphur Springs system may be a “dying” vapor-dominated system that has practically boiled itself dry.     

A shallow attenuating anomaly inside the ring fracture of the Valles Caldera, New Mexico

Spectral ratios of teleseismic direct and scattered P waves observed in the Valles Caldera, New Mexico, show a systematic pattern of low amplitudes at sites inside the caldera relative to sites on or outside the ring fracture. Waveforms recorded at caldera stations are considerably more complex than those recorded outside the caldera. The data used in this study were collected during a passive seismic monitoring experiment conducted in 1987. Twenty-four teleseismic events were recorded on two linear arrays spanning the caldera. To first order, the pattern of low amplitudes did not vary with source incidence angle or azimuth of approach, and could not be explained by anomalous amplification at the ring fracture. This observation suggests the presence of a shallow, attenuating zone associated with the caldera fill material inside of the ring fracture. We estimated the general features of the caldera's near-surface structure for the two-dimensional vertical cross section beneath the array, using a modification of the Aki-Larner discrete-wavenumber method to forward model the observed amplitude variations. Our results indicate that the caldera fill material must be subdivided into at least two distinct zones: a strongly attenuating lower zone, extending to depths in excess of 4 km, and a mildly attenuating surface layer. To fit the data we had to assign an unrealistically low value to seismic Q in the deeper attenuating anomaly. We attribute this to the inability of the Aki-Larner method to account for strong re-direction of energy away from the caldera due to local heterogeneity that we could not include within the low-Q anomaly. This interpretation is consistent with the pervasive, fractured hydrothermal system that is known to exist in the caldera fill material.     

Geochemistry of the volcanics of central Mexico

The Tertiary and Recent volcanics of Mexico occur in two provinces. The Cordillera Province is made up of about 1700 m of ignimbrite sheets overlain and intercalated in the upper part by olivine basalt and basaltic andesite. The Rio Lerma Province extends transversely across Mexico and in the Valley of Mexico the lavas consist mainly of andesite and dacite, 68 % of those analysed having 62 | 4.7 % SiO₂. A total of 108 chemical analyses were made for the major elements, 90 of these including determinations of Cr, Ni, Cu, Zn, Rb, Sr, Ba, Th and Pb for two areas, the Valley of Mexico in the Rio Lerma Province and the Guadalajara region which lies at the intersection of the two provinces. Computer constructions of normative components in the basalt tetrahedron and other projections support an origin of partial melting of tholeiitic to pyrolitic material for the production of andesite. The Guadalajara lavas have consistently higher K/SiO₂ and K/Rb ratios and lower Mg/SiO₂ ratio than the Valley of Mexico rocks suggesting generation at greater depth.     

Caldera formation at Volcán Colima, Mexico, by a large large holocene volcanic debris avalanche

About 4,300 years ago, 10 km³ of the upper cone of ancestral Volcán Colima collapsed to the southwest leaving a horseshoe-shaped caldera 4 km in diameter. The collapse produced a massive volcanic debris avalanche deposit covering over 1550 km² on the southern flanks of the volcano and extending at least 70 km from the former summit. The avalanche followed a steep topographic gradient unobstructed by barriers, resulting in an unusually high area/volume ratio for the Colima deposit. The apparent coefficient of friction (fall height/distance traveled) for the Colima avalanche is 0.06, a low value similar to those of other large-volume deposits. The debris avalanche deposit contains 40–75% angular volcanic clasts from the ancestral cone, a small proportion of vesicular blocks that may be juvenile, and in distal exposures, rare carbonate clasts plucked from the underlying surface by the moving avalanche. Clasts range in size to over 20 m in diameter and are brecciated to different degrees, pulverized, and surrounded by a rock-flour matrix. The upper surface of the deposit shows prominent hummocky topography with closed depressions and surface boulders. A thick, coarse-grained, compositionally zoned scoria-fall layer on the upper northeastern slope of the volcano may have erupted at the time of collapse. A fine-grained surge layer is present beneath the avalanche deposit at one locality, apparently representing an initial blast event. Most of the missing volume of the ancestral volcano has since been restored at an average rate of 0.002 km³/yr through repeated eruptions from the post-caldera cone. As a result, the southern slope of Volcán Colima may again be susceptible to collapse. Over 200,000 people are now living on primary or secondary deposits of the debris avalanche, and a repetition of this event would constitute a volcanic disaster of great magnitude.Ancestral Volcán Colima grew on the southern, trenchward flank of the earlier and larger volcano Nevado de Colima. Trenchward collapse was favored by the buttressing effect of Nevado, the rapid elevation drop to the south, and the intrusion of magma into the southern flank of the ancestral volcano. Other such trenchward-younging, paired volcanoes are known from Mexico, Guatemala, El Salvador, Chile, and Japan. The trenchward slopes of the younger cones are common sites for cone collapse to form avalanche deposits, as occurred at Colima and Popocatepetl in Mexico and at San Pedro Volcano in Chile.     

ESR dating of quartz phenocrysts in the El Cajete and Battleship Rock Members of Valles Rhyolite, Valles Caldera, New Mexico

The electron spin resonance (ESR) dating method was employed on quartz phenocrysts separated from pumice of the El Cajete and Battleship Rock Members of the Valles Rhyolite in the Valles caldera, New Mexico. The results of heating experiments indicate that Ti impurity centers have two components; a thermally stable one and a less stable, temperature sensitive one. ESR dates using the stable Ti center yield eruption ages of 59 ± 6 ka for the Battleship Rock Member and 53 ± 6 ka for the El Cajete Member while recent ¹⁴C dates (S. Reneau and J. Gardner, unpub. data) from carbonized logs in the El Cajete pumice indicate that its age is older than 50 ka. Our results indicate that volcanism in the Valles caldera is much younger than previously thought (≥ 130 ka) and that recent revisions to the post-0.5 Ma stratigraphy of Valles caldera are probably in error. The results suggest that ESR dating of quartz may be a useful method for obtaining ages of units in other Quaternary volcanic areas.     

Shallow Miocene basaltic magma reservoirs in the Bahia de Los Angeles basin, Baja California, Mexico

The basement in the Bahía de Los Angeles basin consists of Paleozoic metamorphic rocks and Cretaceous granitoids. The Neogene stratigraphy overlying the basement is formed, from the base to the top, by andesitic lava flows and plugs, sandstone and conglomeratic horizons, and Miocene pyroclastic flow units and basaltic flows. Basaltic dikes also intrude the whole section. To further define its structure, a detailed gravimetric survey was conducted across the basin about 1 km north of the Sierra Las Flores. In spite of the rough and lineal topography along the foothills of the Sierra La Libertad, we found no evidence for large-scale faulting. Gravity data indicates that the basin has a maximum depth of 120 m in the Valle Las Tinajas and averages 75 m along the gravimetric profile. High density bodies below the northern part of the Sierra Las Flores and Valle Las Tinajas are interpreted to be part of basaltic dikes. The intrusive body located north of the Sierra Las Flores is 2.5 km wide and its top is about 500 m deep. The lava flows of the top of the Sierra Las Flores, together with the distribution of basaltic activity north of this sierra, suggests that this intrusive body continues for 20 km along a NNW-trending strike. Between the sierras Las Flores and Las Animas, a 0.5-km-wide, 300-m-thick intrusive body is interpreted at a depth of about 100 m. This dike could be part of the basaltic activity of the Cerro Las Tinajas and the small mounds along the foothills of western Sierra Las Animas. The observed local normal faulting in the basin is inferred to be mostly associated with the emplacement of the shallow magma reservoirs below Las Flores and Las Tinajas.     

Eruptive history of Fisher Caldera,Alaska, USA

New field, compositional, and geochronologic data from Fisher Caldera, the largest of 12 Holocene calderas in Alaska, provide insights into the eruptive history and formation of this volcanic system. Prior to the caldera-forming eruption (CFE) 9400 years ago, the volcanic system consisted of a cluster of several small (∼3 km³) stratocones, which were independently active between 66±144 and 9.4±0.2 ka. Fisher Caldera formed through a single eruption, which produced a thick dacitic fall deposit and two pyroclastic-flow deposits, a small dacitic flow and a compositionally mixed basaltic-dacitic flow. Thickness and grain-size data indicate that the fall deposit was dispersed primarily to the northeast, whereas the two flows were oppositely directed to the south and north. After the cataclysmic eruption, a lake filled much of the caldera during what may have been a significant quiescent period. Volcanic activity from intracaldera vents gradually resumed, producing thick successions of scoria fall interbedded with lake sediments. Several Holocene stratocones have developed; one of which has had a major collapse event. The caldera lake catastrophically drained when a phreatomagmatic eruption generated a large wave that overtopped and incised the southwestern caldera wall. Multiple accretionary-lapilli-bearing deposits inside and outside the caldera suggest significant Holocene phreatomagmatic activity. The most recent eruptive activity from the Fisher volcanic system was a small explosive eruption in 1826, and current activity is hydrothermal. Late Pleistocene to Holocene magma eruption rates range from 0.03 to 0.09 km³ ky⁻¹ km⁻¹, respectively. The Fisher volcanic system is chemically diverse, ∼48–72 wt.% SiO₂, with at least seven dacitic eruptions over the last 82±14 ka that may have become more frequent over time. Least squares calculations suggest that prior to the CFE, Fisher Volcano products were not derived from a single, large magma reservoir, and were likely erupted from multiple, compositionally independent magma reservoirs. After the CFE, the majority of products appear to have derived from a single reservoir in which magma mixing has occurred.     

Recent crustal subsidence at Yellowstone Caldera,Wyoming

Following a period of net uplift at an average rate of 15±1 mm/year from 1923 to 1984, the east-central floor of Yellowstone Caldera stopped rising during 1984–1985 and then subsided 25±7 mm during 1985–1986 and an additional 35±7 mm during 1986–1987. The average horizontal strain rates in the northeast part of the caldera for the period from 1984 to 1987 were: ₁ = 0.10 ± 0.09 strain/year oriented N33° E±9° and ₂ = 0.20 ± 0.09 strain/year oriented N57° W±9° (extension reckoned positive). A best-fit elastic model of the 1985–1987 vertical and horizontal displacements in the eastern part of the caldera suggests deflation of a horizontal tabular body located 10±5 km beneath Le Hardys Rapids, i.e., within a deep hydrothermal system or within an underlying body of partly molten rhyolite. Two end-member models each explain most aspects of historical unrest at Yellowstone, including the recent reversal from uplift to subsidence. Both involve crystallization of an amount of rhyolitic magma that is compatible with the thermal energy requirements of Yellowstone's vigorous hydrothermal system. In the first model, injection of basalt near the base of the rhyolitic system is the primary cause of uplift. Higher in the magmatic system, rhyolite crystallizes and releases all of its magmatic volatiles into the shallow hydrothermal system. Uplift stops and subsidence starts whenever the supply rate of basalt is less than the subsidence rate produced by crystallization of rhyolite and associated fluid loss. In the second model, uplift is caused primarily by pressurization of the deep hydrothermal system by magmatic gas and brine that are released during crystallization of rhyolite and them trapped at lithostatic pressure beneath an impermeable self-sealed zone. Subsidence occurs during episodic hydrofracturing and injection of pore fluid from the deep lithostatic-pressure zone into a shallow hydrostatic-pressure zone. Heat input from basaltic intrusions is required to maintain Yellowstone's silicic magmatic system and shallow hydrothermal system over time scales longer than about 10⁵ years, but for the historical time period crystallization of rhyolite can account for most aspects of unrest at Yellowstone, including seismicity, uplift, subsidence, and hydrothermal activity.     

Evolution of the Olympus Mons Caldera,Mars

Synoptic images of the Martian volcano Olympus Mons are of a quality and quantity that are unique for mars and, somewhat surprisingly, are appreciably better than image data that exist for many volcanoes on Earth. Useful information about the evolution of shield volcanoes on Earth can thus be derived from the investigation of this extraterrestrial example. We have used shadow-length measurements and photoclinometrically derived profiles to supplement and refine the topographic map of the Olympus Mons caldera. As much as 2.5 km of collapse took place within the 80×65 km diameter caldera and the elevation of the caldera rim varies by almost 2.0 km (low around the oldest collapse events, high around the youngest). An eight-stage evolutionary sequence for the caldera of Olympus Mons is identified which shows that caldera subsidence was a longterm process rather than the near-instantaneous event that has been interpreted from comparable terrestrial examples. Tectonic features on the caldera floor indicate a transition from an extensional environment (graben formation) around the perimeter of the caldera to compression (ridge formation) towards the caldera center. This transition from a compressional to extensional environment is surprisingly sudden, occurs at a radial distance of 17 km from the caldera center, and is import because it can be used to infer that the magma chamber was relatively shallow (thought to be at a depth of <16 km beneath the caldera floor; Zuber and Mouginis-Mark 1990). Ample evidence is also found within the Olympus Mons caldera for solidified lava lakes more than 30 km in width, and for the localzed overturning and/or withdrawal of lava within these lakes.     

Geochemistry of Cenezoic volcanic rocks, Baja California, Mexico: Implications for the petrogenesis of post-subduction magmas

Late Cenozoic volcanism in Baja California records the effects of cessation of subduction at a previously convergent, plate margin. Prior to 12.5 m.y., when subduction along the margin of Baja ceased, the predominant volcanic activity had a calc-alkaline signature, ranging in composition from basalt to rhyolite. Acidic pyroclastic activity was common, and possibly represented the westermost, distal edge of the Sierra Madre Occidental province. After 12.5 m.y., however, the style and composition of the magmatic products changed dramatically. The dominant rock type within the Jaraguay and San Borja volcanic fields is a magnesian andesite, with up to 8% MgO at 57% SiO₂, low Fe/Mg ratios, and high Na/K ratios. These rocks have unusual trace-element characteristics, with high abundances of Sr (up to 3000 ppm), low contents of Rb; K/Rb ratios are very high (usually over 1000, and up to 2500), and Rb/Sr ratios are low (less than 0.01). Furthermore, La_n/Yb_n ratios are high, consistent with derivation from a mantle source with fractionated REE patterns. ⁸⁷Sr/⁸⁶Sr ratios are less than 0.7048, and usually less than 0.7040, whereas the pre-12.5 m.y. lavas have ⁸⁷Sr/⁸⁶Sr ratios between 0.7038 and 0.7063. We have previously termed these rocks bajaites, in order to distinguish them from other magnesian andesites. Bajaites also occur in southernmost Chile and the Aleutian Islands, areas which also have histories of attempted or successful ridge subduction.It is proposed that the bajaite series is produced during the unusual physico-chemical conditions operating during the subduction of young oceanic lithosphere, or subduction of a spreading centre. During normal subduction, the oceanic crust dehydrates, releasing volatiles (water, Rb and other large-ion lithophile elements) into the overlying wedge. Subduction of younger crust will result in a progressive decrease, and eventual cessation of the transfer of volatiles when subduction stops. Thermal rebound of the mantle may cause the slab to melt, perhaps under eclogitestable conditions. The resulting melt will be heavy-REE-depleted, perhaps dacitic, but will otherwise inherit MORB-like Rb/Sr and K/Rb ratios. The ascending melt will react with the mantle to form the source of the bajaitic rocks. Furthermore, any amphibole in the mantle, stabilised during the higher PH₂O conditions of earlier subduction, will break down and contribute a high-K/Rb ratio component.The implications of this study are that firstly, the subducted slab does not contribute a highly fractionated REE component in most modern arcs (i.e. the slab does not melt); secondly, Rb has a very short residence time in the mantle, and its abundance in arc rocks is a direct reflection of the input from the dehydrating slab; and thirdly, bajaitelike rocks may provide recognition of attempted or successful ridge subduction in the geologic past.     

Spatial Distribution,Scaling and Self-similar Behavior of Fracture Arrays in the Los Planes Fault,Baja California Sur,Mexico

We present a case of detailed analysis of fracture arrays spanning four orders of magnitude in length; all of them measured at a single natural site by acquiring images at progressively larger scales. There is a high dispersion of cumulative-length exponents, box dimensions and fracture densities. However, the fractal analysis supports the fractal nature of fracture arrays. Our data indicate the existence of an upper limit for the density parameters, as similarly reported by other authors. We prove that box dimension is in inverse relation with fracture concentration and in direct relation with fracture density. These relations are also observed in our data and additionally there is an upper limit for the box dimensions. We interpret the dispersion in our results as more fundamental than methodological problems. It represents a truncation in the complete evolution of the fracture systems because in natural cases strain initiates overprinting of previous fracture arrays. Considering that larger fractures accommodate strain more efficiently than small fractures, the generation of small fractures is inhibited in the presence of pre-existing larger fractures. Maximum values of fracture density prevent accommodating an excess of strain in a single or restricted range of scales; we claim this condition produces migration of fracturing to larger scales originating fracture scaling.     

Evidence from fluid inclusions for a paleogeothermal gradient at the geothermal test well sites, Los Alamos, New Mexico

Homogenization temperatures of individual fluid inclusions from the geothermal test well sites near Los Alamos, New Mexico, systematically change as a function of depth in the cores. Inclusions in samples from depths between 1.5 and 3.0 km have re-equilibrated to thermal gradients higher than the present gradient of 50–60°C/km. The loci of maximum temperatures attained has a slope of about 70°C/km; the deepest sample has cooled to 200°C from a maximum of 230°C. The wide range of salinities (0.0 wt.% equivalent NaCl to more than 25 wt.% equivalent NaCl) observed in each sample indicates a large amount of pervasive fluid circulation had not occurred at the time of re-equilibration of these inclusions. The results are relevant to calculations for the thermal history of the test site.     

Los Azufres silicic center (Mexico): inference of caldera structural elements from gravity, aeromagnetic, and geoelectric data

Los Azufres geothermal field is located within a silicic volcanic complex in central Mexico. The complex is one of the major silicic centers in the Trans-Mexican Volcanic Belt (TMVB). Pradal and Robin (1985) first suggested the existence of the Los Azufres caldera, and Ferrari et al. (1991) recognized the existence of a collapse structure. According to Pradal and Robin this is a caldera of resurgent type. This geophysical study aims to contribute to the knowledge of the structure of the Los Azufres area. Gravity, aeromagnetic, magnetotelluric (MT) and d.c. vertical electric-resistivity soundings were analyzed. Results show that Los Azufres is a very structurally complex setting with relatively thin crust caused by the extensional tectonics characterizing this central sector of the TMVB. Faults belonging to the E-W to NE-SW (extensional neotectonics) and NW-SE (Basin and Range province) systems are observed to affect the geologic units of Los Azufres. According to our study, the Los Azufres geothermal field is located in a structural high located in the middle of a sub-circular depression delimited to the north-northeast by the Santa Ines Range, and to the southwest by the Mil Cumbres formation. The larger depression consists of two narrow, deep depressions that correspond to La Venta and to the Valley of Juarez. They are separated by the above mentioned structural high. These sub-depressions are believed to be the sites of a maximum caldera collapse, and the structural high is interpreted to be at least in part the caldera's resurgent dome. Geoelectric structure of the caldera derived from d.c. resistivity indicates that the brines of the Los Azufres geothermal system ascend along faults, both bounding and internally disrupting the structural high/resurgent dome. A reasonable correlation is observed between gravity and aeromagnetic data.     

Topographic correlations with soil and regolith thickness from shallow‐seismic refraction constraints across upland hillslopes in the Valles Caldera,New Mexico

How rock is weathered physically and chemically into transportable material is a fundamental question in critical‐zone science. In addition, the distribution of weathered material (soil and intact regolith) across upland landscapes exerts a first‐order control on the hydrology of watersheds. In this paper we present the results of six shallow seismic‐refraction surveys in the Redondo Mountain region of the Valles Caldera, New Mexico. The P‐wave velocities corresponding to soil (≤ 0.6 km s^?1) were inferred from a seventh seismic survey where soil‐thickness data were determined by pit excavation. Using multivariable regression, we quantified the relationships among slope gradient, aspect, and topographic wetness index (TWI) on soil and regolith (soil plus intact regolith) thicknesses. Our results show that both soil and regolith thicknesses vary inversely with TWI in all six survey areas while varying directly with slope aspect (i.e. thicker beneath north‐facing slopes) and inversely with slope gradient (i.e. thinner beneath steep slopes) in the majority of the survey areas. An empirical model based on power‐law relationships between regolith thickness and its correlative variables can fit our inferred thicknesses with R² ‐values up to 0.880 for soil and 0.831 for regolith in areas with significant topographic variations. These results further demonstrate the efficacy of shallow seismic refraction for mapping and determining how soil and regolith variations correlate with topography across upland landscapes. Copyright © 2016 John Wiley & Sons, Ltd.     

Asymmetric,multiple-block collapse at Rotorua Caldera,Taupo Volcanic Zone,New Zealand

Calderas worldwide have been classified according to their dominant collapse styles, although there is a good deal of speculation about the processes involved. Recent laboratory experiments have tried to constrain these processes by modelling magma withdrawal and observing the effects on overlying materials. However, many other factors also contribute to final caldera morphology. Rotorua Caldera formed during the eruption of the Mamaku Ignimbrite. Collapse structure and evolution of Rotorua Caldera is interpreted based its geophysical response, geology and geomorphology, and the stratigraphy of the Mamaku Ignimbrite. Rotorua Caldera is situated at the edge of the extensional Taupo Volcanic Zone, in which major faults strike NE-SW. A second, less dominant fault set strikes NW-SE. These two fault sets have a strong influence on the morphology of Rotorua Caldera. No one style of collapse can be applied to Rotorua Caldera; it was formed during a single eruption, but subsided as many blocks and shows features of trapdoor, piecemeal and downsag types of collapse. Here Rotorua Caldera is described, according to its composition, activity and geometry, as a rhyolitic, single event, asymmetric, multiple-block, single locus collapse structure. The Mamaku Ignimbrite is the only ignimbrite to have erupted from Rotorua Caldera. Extracaldera thickness of the Mamaku Ignimbrite is up to 145 m, whereas inside the caldera it may be greater than 1 km thick. The Mamaku Ignimbrite can be separated into a basal tephra sequence and main ignimbrite sequence. The main ignimbrite sequence contains no observable flow unit boundaries but can be split into lower, middle and upper parts (LMI, mMI, uMI respectively) based on crystal content, welding, jointing, devitrification and vapour phase alteration. Juvenile clasts within the ignimbrite comprise three consanguineous silicic pumice types and andesitic fragments. Only the most evolved pumice type occurs in the basal tephra sequence. All three pumice types occur together throughout the main ignimbrite sequence, whereas the andesitic fragments are only present in uMI. Lithic lag breccias in uMI indicate a late stage of caldera collapse. Concentration of lithic fragments increases towards the middle of the ignimbrite, and may also reflect increased subsidence rate during an earlier stage. Collapse of Rotorua Caldera is thought to have occurred throughout the eruption of the main ignimbrite sequence of the Mamaku Ignimbrite, allowing simultaneous eruption of all the different pumice types and causing the abrupt transition from deposition of the basal tephra sequence to the main ignimbrite sequence.