Eruption Style, Emplacement Dynamics and Geometry of Peralkaline Ignimbrites: Insights From the Lajes-Angra Ignimbrite Formation, Terceira Island, Azores

Ignimbrites are relatively uncommon on ocean island volcanoes and yet they constitute a significant portion of the stratigraphy of Terceira Island (Azores). The Lajes-Angra Ignimbrite Formation (ca. 25 cal ka BP) contains the youngest ignimbrites on Terceira and records two ignimbrite-forming eruptions of Pico Alto volcano that occurred closely spaced in time. Here, we present the first detailed lithofacies analysis and architecture of the Angra and Lajes ignimbrites, complemented by petrographic, mineral chemical, whole rock and groundmass glass geochemical data. The two ignimbrites have the same comenditic trachyte composition, but show considerable variability in trace element and groundmass glass compositions, revealing complex petrogenetic processes in the Pico Alto magma reservoir prior to eruption. The Angra Ignimbrite has a high-aspect ratio and is massive throughout its thickness. It was formed by a small-volume but sustained pyroclastic density current (PDC) fed by a short-lived, low pyroclastic fountain. Overall, the PDC had high particle concentration, granular fluid-based flow conditions and was mostly channelled into a valley on the south part of Terceira. By contrast, the Lajes Ignimbrite has a low-aspect ratio and shows vertical and lateral lithofacies variations. It was formed by a sustained quasi-steady PDC generated from vigorous and prolonged pyroclastic fountaining. The ignimbrite architecture reveals that depositional conditions of the parent PDC evolved as the eruption waxed. The dilute front of the current rapidly changed to a high particle concentration, granular fluid-based PDC that extended to the north and south coasts, with limited capacity to surmount topographic highs. Contrary to what is commonly assumed, the low-aspect ratio of the Lajes Ignimbrite is interpreted to result from deposition of a relatively low velocity PDC over a generally flat topography. This work highlights that the geometry (aspect ratio) of ignimbrites does not necessarily reflect the kinetic energy of PDCs and thus should not be used as a proxy for PDC emplacement dynamics. Although the probability of an ignimbrite-forming eruption on Terceira is relatively low, such a scenario should not be underestimated, as a future event would have devastating consequences for the island’s 55,000 inhabitants.

The reconstruction of ignimbrite-forming eruptions on islands is challenging, as the deposits are often incompletely preserved on account of the small area of deposition, the typical rough topography, and the high erosion rates on volcanic islands. Small oceanic islands are particularly vulnerable to explosive eruptions and the risk is potentially aggravated by their geographical isolation and the weak socioeconomic conditions of local communities. Therefore, detailed studies are needed to widen our understanding of ignimbriteforming eruptions on oceanic islands and their potential impacts. Here, we focus on the youngest ignimbrite formation from Terceira, the Lajes-Angra Ignimbrite Formation (ca. 25 cal ka BP), hereafter abbreviated as LAI following Gertisser et al. (2010). The LAI is composed of two closely related ignimbrites (Lajes and Angra members) originating from Pico Alto volcano that together extend over nearly two-thirds of the island (Self, 1974;Self, 1976;Gertisser et al., 2010). Although closely related, the two ignimbrites contrast strongly in their spatial distribution, thickness and architecture.
In this paper, we use lithofacies analysis and ignimbrite architecture coupled with petrography, mineral chemistry, and whole rock and glass compositions to provide insights into the pre-eruptive magmatic system, and to infer the eruption style, transport and depositional processes of the parent PDCs of Angra and Lajes ignimbrites. We also reassess the use of aspect ratio as a proxy for emplacement dynamics of PDCs. The insights gained from this study are fundamental in improving our knowledge of the future eruptive behaviour of Pico Alto volcano and have significant implications for the assessment of volcanic hazard on Terceira and on other small islands prone to this type of volcanic activity.

Ignimbrite Geometry and Pyroclastic Density Current Emplacement Dynamics
Ignimbrite geometry, or aspect ratio (Walker et al., 1980;Walker, 1983), is one of the most common criteria used to classify ignimbrites. The aspect ratio is determined as the ratio of the average thickness of the deposit (H) to the diameter of a circle with an equal area to that covered by the ignimbrite (L). Accordingly, ignimbrites are subdivided into two endmembers: high-aspect ratio ignimbrites (HARIs) and low-aspect ratio ignimbrites (LARIs). The aspect ratio of ignimbrites has also been widely used as a proxy for PDC emplacement dynamics (i.e., as an estimate of kinetic energy, and hence velocity and competence).
HARIs are characteristically thick when compared to their horizontal extent (H/L 10 -2 -10 -3 ) and are usually confined to valleys or topographic depressions. Well known examples include the 0.76 Ma Bishop Tuff Ignimbrite, California (Wilson and Hildreth, 1997), the AD 1912 Valley of Ten Thousand Smokes Ignimbrite, Alaska (Hildreth and Fierstein, 2012), and the ignimbrite of the AD 1991 eruption of Mt. Pinatubo, Philippines (Scott et al., 1996). HARIs are thought to result from deposition of low energy PDCs (e.g., Walker, 1983;Druitt, 1998;Cas et al., 2011).
However, the aspect ratio as a simple geometric relationship does not correlate with lithofacies variations observed in ignimbrites (Giordano and Doronzo, 2017), which mainly reflect flow-boundary zone conditions in PDCs (Branney and Kokelaar, 2002). Therefore, it can be misleading to directly derive fundamental differences in emplacement dynamics of PDCs from the geometry of the deposits. Study of ignimbrite lithofacies and architecture is thus of paramount importance for better understanding of the temporal and spatial evolution of PDCs.

GEOLOGICAL SETTING AND EVOLUTION OF TERCEIRA ISLAND
The Azores archipelago consists of nine volcanic islands spread over 600 km in the Atlantic Ocean that sit astride the triple junction of the North American, Eurasian and Nubian lithospheric plates (Figure 1 inset). Terceira Island comprises four overlapping central volcanoes: Cinco Picos, Guilherme Moniz, Pico Alto, and Santa Bárbara; as well as a fissure zone ( Figure 1A; Self, 1974;Self, 1976;Madeira, 2005).
Cinco Picos volcano forms the eastern third of Terceira ( Figure 1A) and is dominated by a 7 km-wide caldera that is strongly eroded and buried by products of neighbouring volcanoes. The caldera floor is covered by basaltic lavas and scoria cones of the fissure zone. The products of the Cinco Picos Volcanic Complex are mostly hawaiitic to mugearitic lavas, with minor exposures of comenditic trachyte lavas and pyroclastic deposits (Self, 1974;Self, 1976;Self and Gunn, 1976;Madeira, 2005;Gertisser et al., 2010). The history of this extinct volcano is as poorly known as it is exposed; Cinco Picos is the oldest edifice of the island (dated at 401 ± 6 ka; Hildenbrand et al., 2014), its eruptive activity started in the Pleistocene and the caldera was formed less than 280 ka ago (Féraud et al., 1980). The NE flank of Cinco Picos is deeply dissected by the NW-SE-trending fault scarps of the Lajes Graben ( Figure 1A).
Guilherme Moniz volcano, in the centre of the island, is characterized by an elliptical caldera with a general NW-SE elongation ( Figure 1A), marked by near-vertical walls in the south, while the northern sector is almost entirely buried by recent lava domes and coulées from Pico Alto (Self, 1974;Self, 1976;Madeira, 2005;Gertisser et al., 2010). The floor of the caldera is filled by basaltic lavas of the fissure zone. Guilherme Moniz Volcanic Complex shows a wide range of products from hawaiitic lavas to comenditic coulées and ignimbrites (Self, 1974;Self, 1976;Self and Gunn, 1976;Gertisser et al., 2010). The eruptive activity of this extinct volcano started more than 270 ka and ceased sometime after 111 ka (Calvert et al., 2006).
Pico Alto volcano, located north of Guilherme Moniz, is characterized by a cluster of lava domes and coulées that completely fill and overflow a caldera ( Figure 1A). The Pico Alto Volcanic Complex comprises thick pyroclastic formations, including comenditic trachyte ignimbrites, comenditicpantelleritic trachyte lavas and pumice fall deposits (Self, 1974;Self, 1976;Self and Gunn, 1976;Mungall and Martin, 1995;Pimentel, 2006;Gertisser et al., 2010). The eruptive history of Pico Alto started more than 141 ka ago (Gertisser et al., 2010) and shows evidence of several violent eruptions. The last major ignimbrite-forming period is dated at ca. 25 cal ka BP and is recorded by the LAI. Since then, the volcano has produced many domes and coulées, often preceded by pumice fall deposits (Self, 1974;Self, 1976;Calvert et al., 2006;Pimentel, 2006). Santa Bárbara volcano occupies the western third of the island ( Figure 1A) and is a conical-shape edifice truncated by two small overlapping calderas. The 2 km-wide inner caldera is filled by lava domes and is nested within an older caldera. The slopes of the volcano show several alignments of lava domes and coulées (Self, 1974;Self, 1976;Madeira, 2005;Pimentel, 2006). The Santa Bárbara Volcanic Complex includes a broad range of products, from basaltic-hawaiitic lavas and pyroclastic deposits to comenditic lava domes and coulees with associated pumice fall deposits (Self, 1974;Self, 1976;Self and Gunn, 1976;Mungall and Martin, 1995;Pimentel, 2006). Santa Bárbara is the youngest central volcano on Terceira (dated at 65 ± 13 ka; Hildenbrand et al., 2014); its eruptive activity started during the Late Pleistocene and has continued until the present (Self, 1974;Self, 1976;Calvert et al., 2006).
The fissure zone intersects Terceira along a general WNW-ESE to NW-SE orientation ( Figure 1A). This volcanic system is characterized by a 2 km-wide diffuse alignment of vents (scoria cones, spatter cones and eruptive fissures) best seen in the central part of the island, which has produced basaltic-hawaiitic lava flows (Self, 1974;Self, 1976;Self and Gunn, 1976;Mungall and Martin, 1995;Zanon and Pimentel, 2015;Pimentel et al., 2016). Fissure volcanism extends offshore for several kilometres forming submarine ridges along the flanks of the island (Chiocci et al., 2013;Casalbore et al., 2015). The age of the fissure zone is poorly constrained to more than 43 ka (Calvert et al., 2006) in the SE segment but it is probably as old as the island itself, spanning from the Pleistocene to present. Three historical eruptions are known to have occurred: one paired subaerial eruption in the centre of the island in 1761  and two submarine eruptions along the Serreta Ridge in 1867-2001(Weston, 1964Gaspar et al., 2003;Casas et al., 2018).
The stratigraphy of Terceira is divided into two main groups ( Figure 1B), according to the scheme of Self (1974): the Lower Terceira Group (LTG) and the Upper Terceira Group (UTG). These are separated by the LAI which forms a widespread lithostratigraphic horizon across much of the island. The UTG is formed by the products of at least 116 intercalated eruptions of Pico Alto, Santa Bárbara, and the fissure zone (Self, 1974;Self, 1976).

OVERVIEW OF THE LAJES-ANGRA IGNIMBRITE FORMATION
The LAI is the youngest of at least seven ignimbrite-bearing formations that fall into a relatively narrow period of activity between ca. 86 ka and ca. 25 cal ka BP ago ( Figure 2). It is thought that most of the ignimbrite-forming events were associated with caldera-forming (or caldera-enlarging) phases of Pico Alto and, possibly, Guilherme Moniz volcanoes (Gertisser et al., 2010).
All ignimbrite-bearing formations pre-dating the LAI are part of the LTG. Most are named according to the type locations of the members found on the north and south coasts of the island ( Figure 2). These formations are bounded by major unconformities and have been interpreted to result from single eruptions or, more often, from a few events closely spaced in time (Gertisser et al., 2010). Most of the older ignimbrites on Terceira share similarities with the LAI, including mildly peralkaline compositions, lithofacies, and internal architecture, suggesting that the magmatic and eruptive processes involved in ignimbriteforming eruptions have been recurrent over the last ca. 100 kyr (Gertisser et al., 2010;Jeffery et al., 2017).
The LAI presently crops out along the north and south coasts, and in a few localities in the interior of the island, but originally may have covered nearly two-thirds of Terceira (Self, 1974;Self, 1976;Gertisser et al., 2010). The LAI consists of two members: the Angra Ignimbrite is named after the city of Angra do Heroísmo (in the south of the island, see Figure 1 for location) where thick ignimbrite exposures are found; the Lajes Ignimbrite is named after the village of Lajes (in the NE sector of the island, see Figure 1 for location) where it is best represented and was initially described. The Lajes Ignimbrite was one of the first examples of low-aspect ratio ignimbrites described in the literature (Self, 1971, followed by Walker et al., 1980). The two ignimbrites were initially interpreted by Self (1974), Self (1976) as resulting from two eruptions of Pico Alto volcano that occurred approximately 4 kyr apart. This interpretation was based on three radiocarbon ages (Shotton and Williams, 1973;Shotton et al., 1974) and a quarry outcrop in the centre of the island, where the Lajes Ignimbrite was interpreted possibly overlying the Angra Ignimbrite (Self, 1971). However, in the light of present-day knowledge of Terceira ignimbrite stratigraphy we now think that the lower unit of this outcrop was not the Angra Ignimbrite but an older ignimbrite (possibly the Linhares-Matela Ignimbrite; Gertisser et al., 2010). Unfortunately, this key exposure has been removed by subsequent quarrying.

FIELD AND LABORATORY METHODS
The LAI was studied in detail during field surveys, which combined geological mapping and lithofacies analysis. In total 83 stratigraphic sections were reconstructed throughout the island; each stratigraphic section was thoroughly described, measured and photographed. Maximum clast size (juveniles and lithics) was determined by averaging the length of the longest axis of the three largest clasts. The presence of major unconformities, including palaeosols and erosion surfaces, was used to delimit the deposits.
Lithofacies were described by non-genetic terms, based upon the lithological characteristics of the deposits (e.g., internal structures, grain size, sorting, component variations and intensity of welding). Nomenclature and abbreviations of lithofacies (Table 1) followed the schemes of Branney and   Branney and Kokelaar 2002;Sulpizio et al., 2007).

Name Abbreviation
Massive lapilli-ash mLA Eutaxitic massive lapilli-ash emLA Diffuse-stratified lapilli-ash dsLA Diffuse-bedded lapilli-ash dbLA Massive lithic breccia mlBr Massive ash mA Crystal-rich ash crA Lithic-rich ash lA Pumice-rich lapilli-blocks lens pLBlens Lithic-rich lapilli-blocks lens lLBlens Fines-poor lithic-rich pipe lpip Frontiers in Earth Science | www.frontiersin.org June 2021 | Volume 9 | Article 673686 Kokelaar (2002) and Sulpizio et al. (2007). Terminology used for the description of bed thickness, grain size, and sorting was adopted from Sohn and Chough (1989). Fifteen samples of fresh juvenile clasts were collected for petrographic, whole rock, glass and mineral chemistry analyses. Petrographic descriptions were performed by thinsection observation under an Olympus BX51 binocular microscope. Modal counting was carried out with a CPX-Solutions MicroStepper ID818 point-counter stage, counting at least 1,000 points on each thin-section.
Whole rock major and trace element analyses were undertaken at Actlabs (Activation Laboratories Ltd, Canada), using inductively coupled plasma optical emission spectroscopy (ICP-OES) and inductively coupled plasma mass spectrometry (ICP-MS). Major element compositions of groundmass glass and mineral phases were determined at the Dipartimento di Scienze della Terra "Ardito Desio", Università degli Studi di Milano (Italy) using a JEOL JXA 8200 Superprobe equipped with five wavelength-dispersive spectrometers, an energy-dispersive detector and cathodoluminescence system. Further details on analytical methods are presented in Supplementary Material S2.
Full whole rock and groundmass glass data are given in Supplementary Material S3.

IGNIMBRITE DISTRIBUTION
The spatial distribution of the Angra and Lajes ignimbrites was reconstructed combining direct field measurements with geological interpretation, accounting for the inferred palaeotopography and distance from source (following the method of Wilson, 1991). The Angra Ignimbrite crops out only in the southern part of the island, north of Angra do Heroísmo, where it reaches up to 14 m thickness, and along the coastal cliffs east of the city. Its distribution is constrained to one valley on the southern flank of Guilherme Moniz volcano ( Figure 3) and the reconstructed extent occupies an area of approximately 10 km 2 .
The Lajes Ignimbrite crops out on the northern part of the island, between Biscoitos and Lajes villages, and on the southern part at São Mateus and Porto Judeu villages ( Figure 4). It is best exposed on low-lying areas between Agualva and Lajes villages, where thick ignimbrite (up to 10 m in thickness) is ubiquitously found along the coastal cliffs. In the NE sector of the island, the floor of the Lajes Graben is largely filled by the Lajes Ignimbrite. North of Porto Judeu the outcrops extend into Cinco Picos caldera following a stream valley. Other minor patches of the ignimbrite crop out in the middle of the island in active or abandoned quarries. The reconstructed extent of the Lajes Ignimbrite occupies an area of at least 168 km 2 .
Remnants of a fine ash deposit, commonly weathered into orange-brown soil, crop out over large areas of Terceira ( Figure 4) at the same stratigraphic position of the Angra and Lajes ignimbrites. In places, the top of the ignimbrites can be traced to this ash deposit. The ash deposit occupies a much wider area than that of the ignimbrites and can be found covering topographic highs with a relatively uniform thickness (<50 cm). Where the ignimbrites are not present, this ash deposit forms a valuable lithostratigraphic marker horizon.
Overall, the LAI, including the associated ash deposit, extends radially outward from Pico Alto caldera over most of the island, except for the western third ( Figure 4). It is estimated that the deposits of this ignimbrite formation occupy a minimum area on land of 284 km 2 (i.e., >70% of the surface area of Terceira).

Massive Lapilli-Ash
Description: Massive lapilli-ash (mLA) is the most common lithofacies and typically constitutes the main part of the Angra and Lajes ignimbrites. It is matrix-supported, ranging from well sorted to very poorly sorted (Angra 1.5 < σϕ < 4.1 vs. Lajes 2.7 < σϕ < 5.0; Self, 1974), with variable proportions of juvenile and lithic clasts in a light grey to dark grey matrix of vitric ash and abundant feldspar crystals (up to 3 mm in length) ( Figures 5A,B). Juvenile clasts are rounded to subrounded, ranging from light grey pumice to dark grey/black scoria; black dense clasts are only found in the Lajes Ignimbrite. The juvenile clasts are usually coarsely porphyritic with feldspar crystals up to 5 mm in length (see detailed description of juvenile clasts in section "Petrography and mineral chemistry"). Lapilli-sized clasts are predominant but maximum clast sizes of 20-25 cm are common and locally blocks can reach up to 80 cm. Lithic clasts are subordinate, angular to subrounded and comprise basaltic and trachytic lavas, hydrothermally altered clasts, and syenite xenoliths, in decreasing order of abundance. The lithics are typically <10 cm in size but in places blocks reach up 65 cm. Massive lapilli-ash deposits lack internal stratification, although locally framework-supported pumice-rich or lithic-rich lapilli-blocks lenses (pLBlens, lLBlens) are recognized; fines-poor lithic-rich pipes (lpip) are rare. This lithofacies is non-welded, but is commonly lithified, and contains abundant carbonized wood fragments.
Frontiers in Earth Science | www.frontiersin.org June 2021 | Volume 9 | Article 673686 8 with a fluid escape-dominated flow-boundary zone, although the presence of clast lenses suggests the development of a granular flow-dominated flow-boundary zone (Druitt, 1998;Branney and Kokelaar, 2002;Scarpati et al., 2020). In a fluid escape-dominated flow-boundary zone clasts are supported by the upward flux of interstitial fluid escaping the matrix as a result of hindered settling, while in a granular flow-dominated flow-boundary zone clast interaction is the main support mechanism due to the high particle concentration and granular shear (Branney and Kokelaar, 2002;Sulpizio and Dellino, 2008).

Eutaxitic Massive Lapilli-Ash
Description: Eutaxitic massive lapilli-ash (emLA) is commonly seen in the lower part of the Lajes Ignimbrite but in places can represent more than half of its thickness. In the Angra Ignimbrite, this lithofacies is only found at one location on the costal cliff. The eutaxitic massive lapilli-ash lithofacies is welded and matrixsupported with juvenile and lithic clasts in varying proportions. The matrix is typically dark grey and crystal-rich, with feldspar crystals up to 3 mm in length. The welding intensity varies from poorly to densely welded with well-developed eutaxitic texture ( Figure 5C). Where densely welded, juvenile clasts are strongly flattened dense clasts, which form fiamme up to 30 cm in length. In poorly welded examples, juvenile clasts range from slightly flattened (oblate shapes) to subrounded clasts with scoriaceous texture. Lithic clasts have identical lithologies to those described in the non-welded massive lapilli-ash but are finer-grained (up to 3 cm in size) and less abundant. This lithofacies typically lacks internal stratification but diffuse juvenile-rich or lithic-rich horizons are recognized locally.
Interpretation: The massive nature of this lithofacies suggests that particles were predominately deposited from a current with fluid escape-dominated to granular flow-dominated flowboundary zones (Branney and Kokelaar, 2002;Scarpati et al., 2020), similarly to the non-welded lithofacies. The welding of the deposits, together with the coexistence of fiamme and oblate scoriaceous clasts, suggest rapid progressive aggradation from a high temperature PDC, where syn-depositional agglutination of pyroclasts occurred at temperatures exceeding the glass transition, rather than by load compaction (Russell and Quane, 2005;Dávila-Harris et al., 2013;Scarpati et al., 2020).

Diffuse-Stratified Lapilli-Ash
Description: Diffuse-stratified lapilli-ash (dsLA) is found locally in the Angra and Lajes ignimbrites. It consists of matrix-supported, poorly sorted deposits, with variable proportions of juvenile and lithic clasts in a light grey to dark grey ash matrix rich in feldspar crystals (up to 3 mm in length). This lithofacies is commonly lithified. The juvenile clasts are rounded to subrounded and comprise the same lithologies of those described in the massive lapilli-ash lithofacies. Lapilli-sized clasts are predominant, with maximum clast sizes up to 6 cm. Lithic clasts are angular to subrounded and have the same lithologies of those described in massive lithofacies but are finer-grained (up to 5 cm in size). The internal structure is marked by very thin layers (<1 cm thick) with diffuse parallel to sub-parallel stratification ( Figure 5D). Individual layers are composed of juvenile and/or lithic fine lapilli clasts that alternate with ash layers. In places, the fine lapilli layers show normal or inverse grading, defined by subtle grain size variations. The layers often show small thickness variations and are lateral discontinuous over a few meters.
Interpretation: Diffuse stratification and laterally discontinuous layering are indicative of subtle unsteadiness during deposition from a sustained high particle concentration PDC with transitional flow-boundary zones between fluid escape-dominated and granular flow-dominated, with short periods of tractiondominated conditions (Branney and Kokelaar, 2002;Brown and Branney, 2004;Brand et al., 2014). The unsteadiness within the flow-boundary zone likely arises from one or more mechanisms: 1) successive surges or waves in a fluctuating sustained PDC, 2) periodic impingement into the flow-boundary zone of turbulent eddies, 3) intrinsic frictional effects within the thickening granular flow-dominated flow-boundary zone (Branney and Kokelaar, 2002), and 4) flow-boundary zone perturbations travelling upcurrent in response to interaction with topography (Brand et al., 2014;Báez et al., 2020a).

Diffuse-Bedded Lapilli-Ash
Description: Diffuse-bedded lapilli-ash (dbLA) is only found in the upper part of the Lajes Ignimbrite, along the south coast. It is matrix-supported, poorly sorted, with juvenile and lithic clasts in varying proportions in a dark grey ash matrix with feldspar crystals (up to 3 mm in length). This lithofacies is commonly lithified. The juvenile clasts (rounded to subrounded) and lithic clasts (angular to subrounded) have the same lithologies of those described in the massive lapilli-ash lithofacies. Lapilli-sized clasts are predominant but maximum clast sizes of juveniles and lithics can reach up to 20 and 7 cm, respectively. The internal structure is defined by parallel to sub-parallel beds of juvenile clasts and few lithic clasts ( Figure 5E). Individual beds are centimetre to decimetre thick (typically <25 cm) and laterally discontinuous over a few metres; local thickness variations are common. Most beds are internally massive but some show inverse grading of juvenile clasts and/or normal grading of lithic clasts. The tops and bases of individual beds show diffuse boundaries.
Interpretation: Diffuse bedding suggests local current unsteadiness during deposition from a sustained high particle concentration PDC with granular flow-dominated to fluid escape-dominated flow-boundary zones (Branney and Kokelaar, 2002;Brown and Branney, 2004). The mechanisms that generate the unsteadiness within the flow-boundary zone are the same as described for the diffuse-stratified lithofacies. The occurrence of normal or inverse grading patterns within beds record marked current unsteadiness, possibly with brief periods of tractional deposition. The lateral discontinuous nature and thickness variations of the beds may result from current nonuniformity (Brown et al., 2007;Kokelaar et al., 2007).

Massive Lithic Breccia
Description: Massive lithic breccia (mlBr) is only found in the Lajes Ignimbrite, on the northern part of the island. It consists of matrix-supported, poorly sorted, lithic-rich deposits with a dark grey lapilli-ash matrix rich in feldspar crystals (up to 3 mm in length). Lithic clasts are angular to subrounded, with variable lithologies including basaltic and trachytic lavas, hydrothermally altered clasts (orange to reddish colour), syenite xenoliths, obsidian chips and basaltic scoria clasts, in decreasing order of abundance ( Figure 5F). Lapilli-sized lithic clasts are predominant but larger clasts with maximum clast sizes of 30-35 cm are common. Juvenile clasts are subordinate, rounded to subrounded and comprise scoria and dense clasts. The juvenile clasts reach up to 6 cm in size. The lithofacies is massive, without recognizable internal structures.
Interpretation: The massive nature and poor sorting suggest progressive aggradation from a sustained granular fluid-based lithic-rich PDC with transitional flow-boundary zones between fluid escape-dominated and granular flow-dominated (Branney and Kokelaar, 2002;Brown and Branney, 2004;Dávila-Harris et al., 2013;Báez et al., 2020a), where the high lithic clast concentration was enough to suppress basal turbulence and traction. It reflects a greater competence of the current with respect to the massive lapilli-ash lithofacies. The abundance of lithic clasts in these deposits may indicate significant conduit wall erosion or partial edifice (caldera) collapse (Branney and Kokelaar, 2002;Báez et al., 2020a).

Massive Ash
Description: Massive ash (mA) is only found in the lower part of the Lajes Ignimbrite. It consists of well sorted to poorly sorted (1.5 < σϕ < 3.0; Self, 1974) light grey vitric ash, that locally grades upward to dark grey, abundant feldspar crystals (up to 3 mm in length) and rare juvenile and lithic fine lapilli clasts ( Figure 5G). Juvenile clasts are rounded to subrounded and comprise pumice and scoria (up to 1 cm in size). Lithic clasts are angular to subrounded and have the same lithologies of other lithofacies (up to 0.6 cm in size). The massive ash deposits are structureless, locally lithified, and contain abundant carbonized wood fragments.
Interpretation: The massive nature, absence of internal structures, clast lenses and alignments, and fine grain size suggest steady progressive aggradation from a sustained high particle concentration, fine-grained granular fluid-based PDC with a fluid escape-dominated flow-boundary zone and a minimal granular flow-dominated component (Branney and Kokelaar, 2002;Brown et al., 2007;Kokelaar et al., 2007).

Crystal-Rich (Lithic-Rich) Ash
Description: Crystal-rich ash (crA) is commonly found at the base of Lajes Ignimbrite and locally at the base of the Angra Ignimbrite. It is clast-supported, well sorted, although rare poorly sorted examples exist (Angra 1.0 < σϕ < 2.2 vs. Lajes 1.1 < σϕ < 2.8; Self, 1974), and fines-depleted, mostly composed of feldspar crystals up to 3 mm in length ( Figure 5H). In places, it is rich in ash-sized to fine lapilli lithics (lA). Lithic clasts are angular to subangular and comprise basaltic and trachytic lavas, hydrothermally altered clasts, and syenite xenoliths, in decreasing order of abundance. Juvenile clasts are subordinate, rounded to subrounded, ash-sized to fine lapilli pumice clasts; rare dense clasts are only found in the Lajes Ignimbrite. Small, carbonized wood fragments are common. This lithofacies is typically massive (structureless) but plane-parallel bedding and imbrication of the denser clasts (crystals or lithics) are observed locally.
Interpretation: The typically good sorting, crystal-rich (lithicrich) nature and depletion of fine ash suggest rapid deposition of material from a dilute PDC or dilute part of a PDC (Branney and Kokelaar, 2002;Cas et al., 2011;Scarpati et al., 2015). The depletion of fine vitric ash and enrichment in crystals and lithic clasts implies density segregation during turbulent flow transport, where fine juvenile ash is elutriated. The ingestion and heating of cold air at the front of a PDC generates vigorous turbulence and consequently segregation of the denser clasts (crystals or lithics). The present of carbonized wood fragments reveals the current was hot enough to burn vegetation, which may have also enhanced elutriation of fine ash (Branney and Kokelaar, 2002;Scarpati et al., 2015). The crystal-rich (lithic-rich) ash bed at the base of the ignimbrites is interpreted as a ground layer (Sparks, 1976;Walker et al., 1981;Wilson, 1985;Scarpati et al., 2015).

Vertical Variations
The Angra Ignimbrite typically starts with a light grey, finegrained massive lapilli-ash base, although locally diffusestratified lapilli-ash is observed ( Figure 6A). In a few places a thin fines-poor, crystal-rich bed is present at the base. At one outcrop on the costal cliff (TERL72, see Figure 7 for location), the lowermost part of the ignimbrite is incipiently welded, showing weak eutaxitic texture. Upwards it grades diffusely into non-welded light grey, coarse-grained massive lapilli-ash. Pumice-rich lenses are present within the massive lapilli-ash. The topmost portion of the ignimbrite is typically finer grained.
The Lajes Ignimbrite typically starts with a thin basal finespoor, crystal-rich, locally lithic-rich bed (<4 cm thick on average). This basal crystal-rich bed is locally absent, for example along the south coast. The upper contact of this thin bed is usually sharp but in places can be diffuse. It is conformably overlain by light grey massive ash. It diffusely grades up into dark grey, finegrained massive lapilli-ash, which is commonly welded showing well-developed eutaxitic texture. Upwards it gradually passes through incipiently welded to non-welded dark grey, coarsegrained massive lapilli-ash ( Figure 6B). Pumice-rich or lithicrich lenses are locally present in the massive lapilli-ash. Locally it can grade from diffuse-stratified to diffuse-bedded lapilli-ash. In a few places in the north part of the island it includes a massive lithic breccia. The topmost portion of the ignimbrite is typically finer grained. The vertical contacts between the different lapilliash lithofacies are diffuse.
There is no evidence of unconformities (e.g., major erosion surfaces, reworked horizons, pumice/ash fall layers) in the ignimbrites. Neither the Lajes nor Angra ignimbrites have pumice fall deposits at the base or the top.

Lateral Variations
No significant lateral variations of lithofacies were observed in the Angra Ignimbrite as it has a very restricted outcropping area, which is topographically constrained to one valley (Figure 7).
The Lajes Ignimbrite shows lateral lithofacies variations along the north and south coastal areas (Figures 7, 8). On the north coast, where the Lajes Ignimbrite is best exposed, it is possible to observe lateral variations closely related to the underlying topography. Along the coastal cliffs, the deposit has a rather uniform thickness (3.5 m on average) and shows the same vertical lithofacies sequence (described previously) over several hundreds of metres ( Figure 9). However, thicker deposits (up to 10 m thick) are found in palaeovalleys and tend to be coarser-grained. Where thick, coarse-grained ignimbrite gradually thins against a topographic feature, it can be traced uphill to thin (typically <1 m) fine-grained massive lapilli-ash or massive ash deposit with crystal-rich base. This lateral transition is gradational and generally marked by the thinning and fining of the sequence.

Petrography and Mineral Chemistry
The juvenile clasts of the LAI include light grey pumices ( Figure 10A), which may be porphyritic, dark grey/black coarsely porphyritic scoriae ( Figure 10C), and black coarsely porphyritic dense clasts ( Figure 10E) that are exclusive to the Lajes Ignimbrite. Pumices are highly vesicular (53-81 vol. %) with low phenocryst contents (typically <10 vol. % but up to 20 vol. % in some instances). Phenocrysts are set in a whitish glassy groundmass with few to no microlites. Vesicles range from several millimetres to a few micrometres in size and have subcircular to slightly elongated shapes ( Figure 10B). Some examples show extreme vesicle elongation (i.e., tube pumice). Glomerophyric aggregates of phenocrysts are found sporadically.
Scoriae have lower vesicularity (30-47 vol. %) and higher phenocryst content (14-23 vol. %) than pumices. Coarsely porphyritic scoriae show large phenocrysts (4-5 mm) set in brown glassy, or more rarely microcrystalline, groundmass. Glomerophyric aggregates are more common than in pumices. Scoriae are typically microvesicular, with vesicles reaching up to a few millimetres in size ( Figure 10D), and subcircular to irregularly shaped (separated by thicker walls). Vesicles are frequently coalescent and in some cases are strongly stretched. Cryptocrystalline bands of devitrified glass are occasionally found alternating with microlite-free brown glass.
Dense clasts are characterized by low vesicularity (16-27 vol. %) and high phenocryst content (15-22 vol. %). They are coarsely porphyritic and characterized by brown glassy microlite-free or microlite-poor groundmasses. Glomerophyric aggregates of phenocrysts are common. The vesicles are small (not more than a few tens of micrometres in size), ranging from irregularly shaped, elongated to strongly stretched. In the most extreme cases, vesicles are completely collapsed ( Figure 10F). Alternating bands of dark brown and light brown glass are distinguishable and, in many cases, deformed, clearly draping around phenocrysts (i.e., flow structures). Isolated or clustered spherulites are present occasionally. The mineral contents, proportions, and compositions of the Angra and Lajes ignimbrites are similar. Anorthoclase (An 0-4 , Ab 64-71 , Or 25-35 ) phenocrysts represent the predominant mineral phase in the juvenile clasts of both ignimbrites (≤22 vol. %), reaching sizes of up to 5 mm in length, and exhibiting variable morphologies, including euhedral crystals ( Figure 11A), subhedral crystals with rounded edges, crystals with embayments, and large crystals with resorbed cores filled by glassy pockets (which are frequently vesicular; Figure 11B). Small microphenocrysts, and microlites when present, are euhedral and do not show disequilibrium features. Glomerocrysts comprising predominantly anorthoclase are common. Euhedral to subhedral augite (En 24-38 , Fs 20-34 , Wo 39-43 ) phenocrysts are less common (≤1.2 vol. %), reaching up to 3 mm in length. Large phenocrysts are frequently fractured and show disequilibrium features with the host rock (e.g., reaction rims and rounded edges). Microphenocrysts of clinopyroxene can also be found partially or fully included within alkali feldspars. Clinopyroxene phenocrysts are commonly associated with smaller Fe-Ti oxides ( Figure 11C) and acicular apatite. Olivine is rare (<1 vol. %) and occurs as phenocrysts (up to 800 μm), but more frequently as microphenocrysts. Phenocrysts of olivine commonly show disequilibrium features, such as anhedral morphologies, or partially resorbed crystals with reaction rims ( Figure 11D). Microphenocrysts tend to be subhedral. Olivine occurs as either isolated crystals or in glomerophyric aggregates with other mineral phases. Compositionally, they cover a narrow range of chemical compositions of Fo 30-41 . Fe-Ti oxides (≤0.6 vol. %) are euhedral to anhedral and range from phenocryst (up to 400 μm) to microlites. They are frequently found in association with, or as inclusions in, clinopyroxene and olivine (Figures 11C,D). Apatite is an accessory mineral phase and is found as isolated microphenocrysts (<300 μm) or as acicular inclusions within clinopyroxene and olivine ( Figure 11D). In compositional terms it can be classified as fluorapatite (3.2 wt. % F). In addition to mineral phases, LAI juvenile material contains rare, medium-grained syenitic xenoliths comprised predominantly of alkali feldspar (Gertisser et al., 2010;Jeffery et al., 2017).

Whole Rock and Glass Composition
Whole rock major element geochemical data indicate that both the Angra and Lajes ignimbrites are classified as peralkaline trachyte, with SiO 2 contents that range from 65 to 66 wt. %, total alkali contents of 11-12 wt. % ( Figure 12A), and peralkalinity indices (PI molar Na 2 O + K 2 O/Al 2 O 3 ) between 1.00 and 1.12. Groundmass glass analyses for both ignimbrites show slightly more chemical diversity and some minor differences between the two (e.g., Angra total alkalis 9.6 to 11.4 vs. Lajes total alkalis 10.1 to 13.1 wt. %; Angra PI 0.88 to 0.99 vs. Lajes PI 0.91 to 1.29; Figure 12B). In the Al 2 O 3 vs. FeO T classification scheme of Macdonald (1974), all the peralkaline samples of this study (whole rock and groundmass glass) are comenditic trachyte, excluding a small number of glass analyses from the Lajes Ignimbrite which plot as comendite ( Figure 12C). On the basis of major element chemistry alone, the two ignimbrites are difficult to distinguish from each other. However, trace element compositions yield further information. Zirconium concentrations, which are used here as a differentiation index, show a considerable range between 435 and 963 ppm. Plotted against Zr, FeO T concentrations increase with differentiation, whilst MgO and CaO contents decrease ( Figure 13). When plotted against other incompatible elements such as Nb and Rb, the LAI is defined by a positive, linear trend ( Figures 13D,E). Groundmass glass compositions from the Lajes Ignimbrite extend across the full length of this range, whereas whole rock analyses are limited to <760 ppm. The whole rock data for the Angra Ignimbrite are similar to those of the Lajes Ignimbrite; however, the groundmass glass compositions are bimodal, with the majority of data clustering around 790 ppm, and a smaller cluster around 520 ppm. Strontium concentrations are characteristically low for peralkaline magmas (<20 ppm) and exhibit a scattered but broadly negative correlation with Zr ( Figure 13F). Barium exhibits a similar relationship with Zr but also shows considerable scatter and range (55-737 ppm; Figure 13G). Chrondrite-normalised REE patterns are typically LREE enriched relative to the HREE, with both ignimbrites having average La N /Yb N ratios of 9.7 ( Figure 14). Variably negative Eu anomalies occur in both the Angra and Lajes ignimbrites, with average Eu/Eu* of 0.61 and 0.51, respectively. Furthermore, Eu/Eu* values show a negative correlation with Zr concentrations ( Figure 13H).

Ignimbrite Aspect Ratio and Volume
The aspect ratio of the Angra Ignimbrite was estimated based on an average thickness of 8 m and the diameter of a circle with an equal area to that reconstructed for the ignimbrite (10 km 2 ; Figure 3). This yields an aspect ratio of 1:500 (H/L 2 × 10 -3 ) which is in the range of typical values considered for HARIs (Walker et al., 1980;Giordano and Doronzo, 2017). The aspect ratio of the Lajes Ignimbrite was estimated assuming an average thickness of 3.5 m and the diameter of a circle with the same area as that reconstructed for the ignimbrite (168 km 2 ; Figure 4). This corresponds to an aspect ratio of 1:4200 (H/L 2.4 × 10 -4 ). Alternatively, if the maximum run-out distance of the Lajes Ignimbrite is considered, from Pico Alto caldera to Porto Judeu coastal cliffs (15 km; Figure 4), then the aspect ratio of the Lajes Ignimbrite is in the order of 1:8600 (H/L 1.2 × 10 -4 ). Both these values are within the typical range considered for LARIs (Walker et al., 1980;Dade, 2003;Giordano and Doronzo, 2017).
Volumes of the Angra and Lajes ignimbrites were conservatively estimated by multiplying the average thickness of the deposits by the area of the reconstructed extents (Figures 3,  4). The bulk on-land volume estimated for the Angra Ignimbrite is 0.08 km 3 . This converts to a dense rock equivalent (DRE) onland volume of 0.04-0.06 km 3 , based on ignimbrite densities of 1,300-1800 kg/m 3 and a trachytic magma density of 2,500 kg/m 3 (Self, 1974). The Lajes Ignimbrite yields a bulk on-land volume of 0.59 km 3 . This corresponds to a DRE on-land volume of 0.31-0.42 km 3 using the same range of ignimbrite densities and trachytic magma density as for the Angra Ignimbrite.
Estimating physical parameters of old pyroclastic deposits on islands such as Terceira is challenging due to various factors, including the small size of the island, the thick cover of younger UTG products, the luxuriant vegetation and the absence of deeply dissecting streams, which limit the number of exposures. Taking into account the thickness of the outcrops along the north and south coastal cliffs it is reasonable to assume that a significant (unknown) amount of erupted material was deposited at sea, as commonly reported for other islands such as Montserrat (Lesser Antilles), Stromboli (Aeolian Islands), Faial and São Miguel (Azores archipelago) (e.g., Trofimovs et al., 2006;Andronico and Pistolesi, 2010;Pimentel et al., 2015;Porreca et al., 2018). Therefore, the values obtained for the Lajes and Angra ignimbrites are minimal.
Nevertheless, the volume of erupted magma estimated for the Lajes Ignimbrite (0.31-0.42 km 3 DRE) is within the same order of magnitude of other major explosive eruptions recorded in the Azores, such as the 4.6 ka Fogo A Plinian eruption (Fogo volcano, São Miguel) with 0.64 km 3 DRE (Booth et al., 1978) and the 16 ka Santa Bárbara caldera-forming eruption (Sete Cidades volcano, São Miguel) with 0.25-0.54 km 3 DRE (Kueppers et al., 2019). However, these differ from the Lajes Ignimbrite-forming eruption as they produced thick pumice fall deposits in addition to ignimbrites (Pensa et al., 2015;Queiroz et al., 2015).

Magmatic System and Eruption Style
The Angra and Lajes ignimbrites share the same comenditic trachyte whole rock major element composition ( Figure 12) and mineral assemblage dominated by anorthoclase. However, their whole rock trace element and groundmass glass compositions vary considerably (Figure 13). Such variability is also observed in trace element composition of groundmass glasses and in textures and compositions of anorthoclase crystals of the Lajes Ignimbrite (D'Oriano et al., 2017). The data presented in this study are consistent with the presence of a zoned magma reservoir which fed the Angra and Lajes ignimbrite-forming eruptions. This is a recurrent theme for the magma system of Pico Alto volcano (and potentially Guilherme Moniz), where the interplay of complex petrogenetic processes acting in the magma reservoir has been identified (D'Oriano et al., 2017;Jeffery et al., 2017). These include the derivation of the ignimbrite-forming comenditic trachyte from a zoned magma reservoir, with thermometric and hygrometric gradients, dominated by extended fractional crystallization and important contributions from in-situ crystallization, remobilization of cumulate material (crystal mush) and magma mixing processes. All these processes are, to some extent, masked by the seemingly homogeneous trachytic major element compositions of the eruptive products.
Chemically zoned magma reservoirs are common in other peralkaline volcanic systems, such as Furnas volcano on São Miguel (Jeffery et al., 2016), and are recorded by peralkaline ignimbrites such as the 13.9 Ma Ignimbrite "TL" (Sumner and Branney, 2002), the 13.6 Ma Ignimbrite "A" (Troll and Schmincke, 2002), and the 11.8 Ma Ayagaures ignimbrite (Jutzeler et al., 2010), on Gran Canaria, and the 45 ka Green Tuff Ignimbrite, Pantelleria (Williams, 2010), as well as a range of other ocean island systems (Jeffery and Gertisser, 2018 and references therein). Despite this, the data presented in this study cannot be accounted for adequately by a traditional zoned magma reservoir model featuring a single magma body with an upward gradational transition from less to more evolved magma, and may instead highlight a more complex magmatic plumbing system, potentially comprising multiple magma bodies.
Although the eruption triggering mechanism is uncertain, D' Oriano et al. (2017) proposed that the magmatic system of Pico Alto was disturbed by a heat source, most probably induced by the arrival of a hotter magma from depth, but without mixing with the resident comenditic trachyte magma reservoir. The thermal destabilization at the base of the reservoir and the onset of eruption caused syn-eruptive mixing between different portions of the zoned magma reservoir, which contributed to the homogenization of the resident magma. We speculate that the first eruptive event tapped the uppermost part of the reservoir, which is recorded by the Angra Ignimbrite. Some (unknown) time after the initial event, a second eruptive event took place involving the withdrawal of a larger volume of magma and probably reaching the lower parts of the reservoir. The latter is recorded by the Lajes Ignimbrite.
The Angra Ignimbrite-forming eruption was characterized by short-lived low pyroclastic fountaining, most probably located in the  (Macdonald, 1974).
Frontiers in Earth Science | www.frontiersin.org June 2021 | Volume 9 | Article 673686 southern part of Pico Alto caldera. The event did not form a convective eruptive column, as no pumice fall deposit was found in the ignimbrite sequence. The pyroclastic fountain generated a small but sustained PDC that seems to have been mostly channelled along a valley on the southern part of Terceira ( Figure 3). The eruption of the Lajes Ignimbrite was dominated by vigorous and prolonged low pyroclastic fountaining (or boiling-over) from Pico Alto caldera generating a sustained PDC, which spread radially towards both the north and south coasts of the island (Figure 4).
Other ignimbrites worldwide have also been interpreted to result from boiling-over eruptions (e.g., Cas et al., 2011;Willcock et al., 2013;Pacheco-Hoyos et al., 2018). Like the Angra Ignimbrite, the Lajes Ignimbrite-forming eruption was neither preceded nor followed by the formation of a convective eruptive column, as there are no pumice fall deposits associated with this ignimbrite sequence.

Pyroclastic Density Current Emplacement Dynamics
The vertical and lateral lithofacies variations observed in the Angra and Lajes ignimbrites reflect different conditions of the parent PDCs during deposition. Vertical lithofacies variations ( Figure 6) record changes in PDC deposition through time (i.e., current unsteadiness), while lateral lithofacies variations (Figures 7, 8) record spatial variability of the PDC during deposition and its interaction with the topography (i.e., current non-uniformity) (Branney and Kokelaar, 2002). The Angra Ignimbrite is composed of non-welded massive lapilli-ash throughout most of its thickness ( Figure 6A). The massive nature and small volume of this ignimbrite indicate that it was deposited from a small but sustained high particle concentration granular fluid-based PDC. The dilute front of the current, which left in places a thin, fines-poor and crystalrich bed, rapidly changed to high particle concentration, granular fluid-based flow conditions, overall dominated by steady deposition from fluid escape-dominated to granular flowdominated flow-boundary zones (Branney and Kokelaar, 2002). Locally, subtle unsteadiness developed within the flowboundary zone, leading to deposition of diffuse-stratified lapilliash. No unconformities were observed within the ignimbrite, which implies that the emplacement was continuous without significant time breaks, representing a single depositional unit. The Angra Ignimbrite has a high-aspect ratio geometry and limited areal extent, cropping out in only one valley south of Guilherme Moniz caldera, suggesting that the parent PDC may have been mostly channelled along this valley ( Figure 3). However, it cannot be excluded that the Angra Ignimbrite is buried under younger deposits in the middle of the island, particularly on the floor of Guilherme Moniz and Cinco Picos calderas. The features of the Angra Ignimbrite are typical of HARIs deposited from relatively low velocity PDCs confined to topographic depressions (e.g., Walker, 1983;Cas and Wright, 1987;Cas et al., 2011;Martí et al., 2019).
The vertical lithofacies variations observed in the Lajes Ignimbrite reveal an overall inverse grading pattern, starting as a finer-grained base and grading to a coarser-grained upper part. However, the topmost portion of the ignimbrite is typically finer grained ( Figure 6B). From its characteristics, we can infer that the ignimbrite was mostly deposited from a sustained, quasi-steady, high particle concentration granular fluid-based PDC that progressively waxed in competence during emplacement and then rapidly waned at the end (Kneller and Branney, 1995;Branney and Kokelaar, 2002). This reveals that the PDC dynamics changed as the eruption intensified. The initial dilute front of the PDC deposited a basal, fines-poor and crystal-rich bed (similar to the ground layer of other ignimbrites; Sparks, 1976;Scarpati et al., 2015), but it rapidly evolved to a high particleconcentration, granular fluid-based current dominated by steady deposition from fluid escape-dominated to granular flowdominated flow-boundary zones. The high temperature of the PDC promoted the syn-depositional agglutination of pyroclasts and thus welding of the deposit, which is recorded by the eutaxitic massive lapilli-ash. As the PDC progressively waxed in competence, it deposited non-welded, coarse-grained massive lapilli-ash. The increase in the concentration and size of lithic clasts had a cooling effect on the current (Eichelberger and Koch, 1979;Martí et al., 1991), sufficient to inhibit welding of the upper part of the deposit. Local current unsteadiness, with transitional flow-boundary zones between fluid escape-dominated and granular flow-dominated, with brief periods of traction-  et al. (2017). Chondrite values taken from Sun and McDonough (1989). dominated conditions (Branney and Kokelaar, 2002), led to deposition of diffuse-stratified and diffuse-bedded lapilli-ash.
Overall, the vertical sequence of the Lajes Ignimbrite indicates that this thin ignimbrite sheet (3.5 m average thickness) corresponds to a single depositional unit from a sustained quasi-steady PDC. This contrasts with previous studies (Self, 1974;Self, 1976;Gertisser et al., 2010) that interpreted the presence of two depositional layers based on strong pumice concentration zones and sharp boundaries between welded and non-welded lithofacies. However, these pumice concentrations are not ubiquitous; they are found in a few places in the northern part of the island, and the sharp welded/non-welded boundaries are only seen locally at São Mateus (south coast). Several lines of evidence support the single depositional unit, including the gradual transitions of lithofacies (i.e., diffuse boundaries), the gradual variation of welding intensity, and the absence of unconformities such as intra-ignimbrite pumice or ash fall layers and erosion surfaces.
The lateral variations of lithofacies observed in the Lajes Ignimbrite (Figures 7, 8) cannot be interpreted as typical valley-pond ignimbrite to veneer deposit transitions commonly seen in LARIs. The thin fine-grained massive lapilli-ash or massive ash deposit with crystal-rich base observed where the ignimbrite thins against topography differs from the stratified veneer lithofacies described for intermediate to small volume ignimbrites (e.g., Walker et al., 1980;Walker et al., 1981;Wilson and Walker, 1982;Wilson, 1985;Giordano et al., 2002). In the present case, the lateral transition corresponds to a condensed sequence of the Lajes Ignimbrite, including the lateral extent of the ground layer. The absence of stratified lithofacies suggests that the upper part of the current was not as turbulent as in other more dynamic PDCs. This sequence clearly differs from the fine ash deposit (commonly weathered into soil) that displays mantle bedding and is widespread across the island, including on topographic highs (Figure 4), which we interpret to correspond to the co-ignimbrite ash fall deposit.
Thus, the lateral lithofacies variations can be interpreted as recording non-uniformity within the flow-boundary zone (Branney and Kokelaar, 2002). The simultaneous deposition of coarse-grained massive lapilli-ash in a certain location and massive ash in another location may result from current instabilities and relatively non-uniform deposition. This is evident where the Lajes Ignimbrite parent PDC interacted with more irregular topography, such as stream valleys, fault scarps or caldera walls, which led to marked lateral variations of lithofacies in the space of just a few metres distance.
Although the Lajes Ignimbrite has a low-aspect ratio geometry, the parent PDC generally flowed along areas of low-lying ground with limited capacity to overcome topographic highs, such as the faults scarps of the Lajes Graben (10-70 m high) and the caldera walls of Cinco Picos and Guilherme Moniz volcanoes (up to 200 m high). In the western part of the island, the PDC was limited by the slopes of Santa Bárbara volcano (Figure 4). This contrasts with the typical flow behaviour of PDCs that generate LARIs, which are interpreted as highly energetic currents capable of surmounting significant topographic barriers (Walker et al., 1980;Cas and Wright, 1987;Wilson et al., 1995;Cas et al., 2011). Therefore, the low-aspect ratio of the Lajes Ignimbrite may simply result from deposition over a generally flat topography such as the coastal areas and the floor of the Lajes Graben, rather than from a highly energetic PDC. We propose that the ignimbrite was emplaced from a sustained high particle concentration, granular-fluid-based PDC with relatively low velocity and low heat loss. The low energy, placid style of emplacement helped conserve heat and promote welding in the Lajes Ignimbrite. Similar emplacement dynamics have been proposed for large volume ignimbrites such as the 4.3 Ma Kizilkaya ignimbrite, Turkey (Schumacher and Mues-Schumacher, 1996), the 2.08 Ma Cerro Galán Ignimbrite, Argentina (Cas et al., 2011), and the <73 ka Campo de la Piedra Pómez ignimbrite, Argentina (Báez et al., 2020b), but this mechanism apparently also holds for smaller volume ignimbrites.
The different flow paths from Pico Alto caldera towards the north and south coasts of the island (Figure 4) mostly reflect the interaction of the PDC with palaeotopography, in particular with the caldera walls of Guilherme Moniz and Cinco Picos volcanoes, and the fault scarps of the Lajes Graben. The parent PDC seems to have been deflected by the fault scarps of the graben to reach the eastern end of the island at Praia da Vitória. The detailed geochemical reconstruction of the Lajes Ignimbrite indicates that the PDC was not simultaneously emplaced radially but rather different flow paths were emplaced at different times (D'Oriano et al., 2017). Similar flow behaviour has been recognized for the 45 ka Green Tuff Ignimbrite, Pantelleria (Williams et al., 2014) and for the 39 ka Campanian Ignimbrite, Italy (Fedele et al., 2016) through detailed chemical stratigraphy of the deposit.

Implications for Hazard Assessment
The geological record of Terceira shows that several ignimbriteforming eruptions coming from Pico Alto and, possibly, Guilherme Moniz volcanoes have occurred in the last ca. 100 kyr (Gertisser et al., 2010). Although the frequency of ignimbrite-forming eruptions is relatively low and the last events occurred ca. 25 cal ka BP ago (recorded by the LAI), Pico Alto is an active central volcano and the possibility of other similar eruptions in the future should not be underestimated. The probability of an ignimbrite-forming eruption occurring from Santa Bárbara volcano is low, as its stratigraphy does not include ignimbrites, although a lithic-rich PDC deposit has been reported (Calvert et al., 2006). However, Santa Bárbara volcano has the potential to produce explosive and, therefore, ignimbrite-forming eruptions, thus a future event cannot be ruled out.
A future ignimbrite-forming eruption of Pico Alto volcano could strongly impact the island, with devastating and long-term socio-economic consequences for its 55,000 inhabitants (PORDATA, 2020) and for the entire Azores archipelago (Medeiros et al., 2021). An eruption with similar characteristics to the Lajes Ignimbrite-forming event would be a worst-case scenario in terms of hazard assessment for a small island such as Terceira (400 km 2 ). If a future event behaved similarly to the Lajes Ignimbrite parent PDC, it would affect all of the central part of island and flow towards the sea on the north and south coasts (Figure 4). The relatively low velocity of the PDC might mean it would have limited ability to surmount topographic barriers, such as the fault scarps of the Lajes Graben or the caldera walls of Cinco Picos and Guilherme Moniz volcanoes, but it could be channelled to more distal areas maintaining high temperatures.
The well-populated villages along the north coast of the island between Biscoitos and Lajes would be strongly impacted due to their proximity to Pico Alto caldera ( Figure 4). However, other densely populated areas on the eastern part of the island between Lajes and Praia da Vitória, and several villages on the south coast between São Mateus and Porto Judeu, including the city of Angra do Heroísmo could also be affected by channelled PDCs. In this scenario the western part of the island would not be affected and could possibly serve as a retreat area in case of evacuation of the affected villages. Moreover, the hazard would be especially compounded by the lack of an opening Plinian fallout phase that could give warning of impending PDCs and during which evacuations or other mitigation operations could be prepared.
Apart from dwellings, all critical infrastructure of the island might be expected to be damaged or destroyed. The geothermal power plant located on the southwest flank of Pico Alto caldera would most likely also be destroyed. All of the island's communication hubs, including the only airport and NATO Air Base on the floor of Lajes Graben and the main harbour located on Praia da Vitória could suffer significant damage (see Figure 4 for location). This would make Terceira difficult to be reached by air or maritime transport, with major implications for evacuations and emergency response. The only hospital and the headquarters of the Azores Regional Civil Protection and Fire Services, both located in the city of Angra do Heroísmo (see Figure 4 for location), could also be severely affected.
In the light of this and despite the relatively low probability of such an event, this worst-case scenario should be considered by regional and national authorities in risk management plans. The potential risk is aggravated by the small size, remote location, high population density and weak economy of Terceira, as is commonly the case of volcanic islands (Pelling and Uitto, 2001;Wilkinson et al., 2016;Pimentel et al., 2020), such as demonstrated recently by the 2021 explosive eruption of La Soufrière volcano, St. Vincent (Lesser Antilles).

CONCLUSIONS
The Lajes-Angra Ignimbrite Formation (LAI) is the youngest (ca. 25 cal ka BP) of at least seven major formations containing peralkaline ignimbrites on Terceira Island (Azores). Most of the older ignimbrites share physical and geochemical similarities with the LAI, suggesting that the processes involved in ignimbrite-forming eruptions have been recurrent in the last ca. 100 kyr, with significant implications for volcanic hazard assessment. In this paper we combined lithofacies analysis and ignimbrite architecture with petrographic and geochemical data to provide insights into the pre-eruptive magmatic system of Pico Alto volcano and infer the eruption style and emplacement dynamics of the parent PDCs of the Angra and Lajes ignimbrites.
The two ignimbrites are difficult to distinguish based on whole rock major element chemistry alone, as both have comenditic trachyte compositions. However, whole rock trace element and groundmass glass compositions show considerable variability, which together with disequilibrium textures and compositions of anorthoclase crystals, reveal the occurrence of complex petrogenetic processes in the Pico Alto magma reservoir prior to eruption. Although the two ignimbrite-forming events were closely spaced in time (unresolvable by radiocarbon dating), their PDCs had dissimilar eruptive and depositional processes. The Angra Ignimbrite was formed by a small-volume sustained PDC that was fed by a short-lived, low pyroclastic fountain, which generated a high particle concentration, granular fluid-based PDC. This current seems to have been mostly channelled along a single valley on the south part of the island, forming a high-aspect ratio ignimbrite. On the other hand, the Lajes Ignimbrite resulted from a sustained quasi-steady PDC that was fed by vigorous and prolonged pyroclastic fountaining. Ignimbrite architecture suggests evolving depositional conditions of the parent PDC through time and space. The dilute front of the current rapidly evolved to a high particle concentration, granular fluid-based PDC as the eruption waxed. The PDC spread radially from Pico Alto caldera to both the north and south coasts with limited capacity to overcome topographic highs. The low-aspect ratio of the Lajes Ignimbrite most probably resulted from deposition of a relatively low velocity PDC over a generally flat topography, rather than from a highly energetic current. This placid style of emplacement allowed the PDC to conserve heat, which in turn, prompted welding in this thin ignimbrite sheet. We highlight that ignimbrite geometry (aspect ratio) does not directly reflect the kinetic energy of a PDC and therefore should not be used to infer on PDC emplacement dynamics.
A future ignimbrite-forming eruption of Pico Alto volcano would have devastating consequences for Terceira. In this worstcase scenario, the central part of the island and the villages located along the north and south coasts would be affected by PDCs, including all critical infrastructure and communication hubs. Despite the relatively low probability, such an event presents a considerable challenge for civil protection authorities. The findings of this study can provide a framework to better understand ignimbrite-forming eruptions on other active volcanic islands, as well as their potential impacts in case of future volcanic unrest.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.