La gran erupción de hace 4.2 ka cal en Cerro Blanco, Zona Volcánica Central, Andes: nuevos datos sobre los depósitos eruptivos holocenos en la Puna sur y regiones adyacentes ; The large eruption 4.2 ka cal BP in Cerro Blanco, Central Volcanic Zone, Andes: Insights to the Holocene eruptive deposit

The eruption of the Cerro Blanco Volcanic Complex, in the southern Puna, NW Argentina dated at 4410–4150 a cal BP, was investigated to produce new information on stratigraphy, geomorphology, physical volcanology, radiocarbon dating, petrography, and geochemistry. Identification of pre–, syn–, and post–caldera products allowed us to estimate the distribution of the Plinian fallout during the paroxysmal syn–caldera phase of the eruption. The new results provide evidence for a major rhyolitic explosive eruption that spread volcanic deposits over an area of ~500,000 km, accumulating >100 km of tephra (bulk volume). This last value exceeds the lower threshold of Volcanic Explosive Index (VEI) of 7. Ash-fall deposits mantled the region at distances >400 km from source and thick pyroclastic-flow deposits filled neighbouring valleys up to 35 km away. This eruption is the largest documented during the past five millennia in the Central Volcanic Zone of the Andes, and is probably one of the largest Holocene explosive eruptions in the world. We have also identified two additional rhyolitic eruptions in the region from two other eruptive sources: one during the Early–Holocene and another in the Late–Holocene. The identification and characterisation of these significant volcanic events provide new constraints into regional Holocene 2 J.L. Fernandez-Turiel et al. Estudios Geológicos, 75(1), enero-junio 2019, e088, ISSN-L: 0367-0449. https://doi.org/10.3989/egeol.43438.515 Introduction Tephra is fragmental material produced by an explosive volcanic eruption. Larger tephra (‘blocks’ or ‘bombs’) follow ballistic trajectories and pose a local threat (usually < 5 km). However, lapilli (2 to 64 mm in diameter) and ash (< 2 mm in diameter) are convected upwards within the eruption column and dispersed by wind or buoyancy forces. Ash particles can be carried hundreds or even thousands of kilometres away from source, being the most frequent, and often widespread, volcanic hazard (Jenkins, S. F. et al., 2015). In this way, tephra deposits form an isochron that directly links various sedimentary successions and enable us to synchronize and date geological, palaeoenvironmental or archaeological records (Lane et al., 2017; Lowe, 2011; Ponomareva et al., 2015). This is the rationale of the development of tephra studies and the establishment of tephrostratigraphy and tephrochronology. In addition, preserved tephra deposits are frequently used in the reconstruction of past volcanic eruptions and the evaluation of volcanic hazards to prevent the impact of similar eruptions possibly occurring in the future (Houghton, B. & Carey, R. J, 2015; Lowe, 2011). In South America, the explosive Holocene volcanism is mainly produced by the Andean Volcanic Arc, along the Andean Cordillera, and is segmented into four main volcanic zones separated from each other by volcanic gaps (Stern, 2004): the Northern (NVZ), the Central (CVZ), the Southern (SVZ), and the Austral Volcanic Zone (AVZ) (Fig. 1). In this setting, although mostly concentrated in SVZ and AVZ, tephrochronology studies are of critical importance because would help in reconstructing dispersal of ash during past explosive eruptions providing potential stratigraphic markers that contribute toward establishing a chronological framework for a range of applications (e.g., palaeoenvironment, palaeoclimate, archaeology, and rates of geomorphic change) (Fontijn et al., 2014; Stern, 2008). Several previous studies have focused on the widespread ash-fall deposits found in the uppermost Quaternary sequences of NW Argentina (Hermanns & Schellenberger, 2008; Malamud et al., 1996; Montero-López et al., 2009; Ruggieri et al., 2010; Trauth et al., 2003).). These studies attributed these ash-fall deposits to eruptions from volcanoes located in the Central Volcanic Zone of the Andes but no reliable correlation to a specific source has been proposed (Hermanns & Schellenberger, 2008; Malamud et al., 1996). This is mostly due to the lack of adequate geochemical data (e.g., single shards glass major and trace element composition, mineral chemistry) and reliable dating (often in contrast with each other) on proximal deposits of the potential geological and archaeological records, and offer extensive regional chronostratigraphic markers over a wide geographical area of South America.


Introduction
Tephra is fragmental material produced by an explosive volcanic eruption. Larger tephra ('blocks' or 'bombs') follow ballistic trajectories and pose a local threat (usually < 5 km). However, lapilli (2 to 64 mm in diameter) and ash (< 2 mm in diameter) are convected upwards within the eruption column and dispersed by wind or buoyancy forces. Ash particles can be carried hundreds or even thousands of kilometres away from source, being the most frequent, and often widespread, volcanic hazard (Jenkins, S. F. et al., 2015). In this way, tephra deposits form an isochron that directly links various sedimentary successions and enable us to synchronize and date geological, palaeoenvironmental or archaeological records (Lane et al., 2017;Lowe, 2011;Ponomareva et al., 2015). This is the rationale of the development of tephra studies and the establishment of tephrostratigraphy and tephrochronology. In addition, preserved tephra deposits are frequently used in the reconstruction of past volcanic eruptions and the evaluation of volcanic hazards to prevent the impact of similar eruptions possibly occurring in the future (Houghton, B. & Carey, R. J, 2015;Lowe, 2011).
In South America, the explosive Holocene volcanism is mainly produced by the Andean Volcanic Arc, along the Andean Cordillera, and is segmented into four main volcanic zones separated from each other by volcanic gaps (Stern, 2004): the Northern (NVZ), the Central (CVZ), the Southern (SVZ), and the Austral Volcanic Zone (AVZ) (Fig. 1). In this setting, although mostly concentrated in SVZ and AVZ, tephrochronology studies are of critical importance because would help in reconstructing dispersal of ash during past explosive eruptions providing potential stratigraphic markers that contribute toward establishing a chronological framework for a range of applications (e.g., palaeoenvironment, palaeoclimate, archaeology, and rates of geomorphic change) (Fontijn et al., 2014;Stern, 2008).
Several previous studies have focused on the widespread ash-fall deposits found in the uppermost Quaternary sequences of NW Argentina (Hermanns & Schellenberger, 2008;Malamud et al., 1996;Montero-López et al., 2009;Ruggieri et al., 2010;Trauth et al., 2003).). These studies attributed these ash-fall deposits to eruptions from volcanoes located in the Central Volcanic Zone of the Andes but no reliable correlation to a specific source has been proposed (Hermanns & Schellenberger, 2008;Malamud et al., 1996). This is mostly due to the lack of adequate geochemical data (e.g., single shards glass major and trace element composition, mineral chemistry) and reliable dating (often in contrast with each other) on proximal deposits of the potential geological and archaeological records, and offer extensive regional chronostratigraphic markers over a wide geographical area of South America.

RESUMEN
La erupción del Complejo Volcánico Cerro Blanco en el sur de la Puna, noroeste de Argentina (4410-4150 a BP) se investigó para obtener nueva información sobre estratigrafía, geomorfología, volcanología física, dataciones por radiocarbono, petrografía y geoquímica. La caracterización de los productos en relación a la evolución de la caldera de Cerro Blanco permitió estimar la distribución de los depósitos de ceniza de la fase paroxísmica Plineana de la erupción. Estos novedosos resultados evidencian una gran erupción explosiva riolítica que generó depósitos cineríticos en un área de aproximadamente 500.000 km 2 , acumulando >100 km 3 de tefra (volumen total). Este último valor supera el umbral inferior del Índice de Explosividad Volcánica (IEV) de 7. Los depósitos de caída de ceniza cubrieron la región, llegando a más de 400 km desde el Complejo Volcánico de Cerro Blanco, y los potentes depósitos de flujos piroclásticos rellenaron los valles vecinos alcanzando una distancia de 35 km. Esta erupción es la más grande documentada durante los últimos cinco milenios en la Zona Volcánica Central de los Andes y es probablemente una de las mayores erupciones explosivas holocenas del mundo. Además, se han identificado otras dos erupciones riolíticas en la región procedentes de otros dos centros eruptivos: una durante el Holoceno temprano y otra en el Holoceno tardío. La identificación y caracterización de estos grandes eventos volcánicos proporcionan nuevas guías para los registros geológicos y arqueológicos regionales del Holoceno, siendo marcadores cronostratigráficos de aplicación a una extensa área geográfica de América del Sur.   source volcanoes. This situation changed with the discovery of several upper Pleistocene deposits attributed to the Cerro Blanco Volcanic Complex (CBVC) (Arnosio et al., 2008) and, particularly, with identification of Holocene deposits associated with the CBVC (Montero-López et al., 2010a). In particular, the presence of Holocene volcanic deposits in the southern Puna led us to review all available data and to conduct new studies of recent pyroclastic fall deposits and ignimbrites in a large region in NW Argentina in order to better identify their sources (Fernandez-Turiel et al., 2013. We present here the results of geomorphology, stratigraphy and geographical distribution of tephra deposits, the detailed characterization of tephra (including particle size, mineralogy, and geochemical composition), and the radiocarbon dating of these deposits to test their possible association with large Holocene eruptions in the CBVC. This work mainly focuses on details of the Cerro Blanco pyroclastic fall deposits and ignimbrites, but the stratigraphic study of regional distribution requires us to include the analysis of bracketing Cueros de Purulla (older) and Bolsón de Fiambalá (younger) ash deposits.

Geological setting
The 3500-4700 m high Cerro Blanco Volcanic Complex (CBVC) is part of the Cordillera de San Buenaventura in NW Argentina (26 o 47' 24.46" S, 67 o 45' 40.68" W) (Global Volcanism Program, 2013b). It overlies the Miocene to Pliocene La Hoyada Volcanic Complex (~7 to ~2 Ma), which consists of ignimbrites, lava domes, and lava flows (Montero-López et al., 2010b). It is situated on the southern border of the Altiplano-Puna Plateau, in the Central Volcanic Zone of the Andes ( Fig. 1a and b), where voluminous mafic and silicic Miocene-Quaternary calcalkaline volcanism occurs within the plateau and along its margins (Bianchi et al., 2013;Kay et al., 2006Kay et al., , 2010. The underlying continental crust is composed of Palaeozoic-Mesozoic igneous and metasedimentary rocks (Montero-López et al., 2010a).

Methods
Identification, logging, and sampling of 62 outcrops of tephra sections and lava domes were performed during several field campaigns in a wide region of northwest Argentina (Fig. 1b), providing a set of more than 230 samples. The lithofacies descriptions are based on published methodologies (Branney & Kokelaar, 2002;. Petrography was studied by optical and scanning electron microscopy on polished thin sections and grain mounts on carbon stubs. Groundmass glass, phenocrysts, microphenocrysts, and microlites of glass shards were analysed on polished thin sections by electron probe microanalyser (EPMA) (Fernandez-Turiel et al., 2018a). Particle size distributions were determined in ash samples by laser diffraction (Fernandez-Turiel et al., 2018b).
Radiocarbon dating was done by Accelerator Mass Spectrometry (AMS) in Beta Analytic Inc. and NSF-Arizona AMS Facility, University of Arizona. All samples were taken below the surface, visually checked under a microscope, and handpicked cleaned. The 14 C ages were calibrated to calendar years BP using Bayesian statistics with OxCal v4.3.2 [104] (round to the nearest 10) (Bronk Ramsey, 2009) with the . CB 1 , precaldera block-and-ash flow deposits; CB 2 2, Plinian ash-fall deposits; CB 2 3, ignimbrites, CB 3 1, lava-domes; CB 3 2, lapilli fall deposits; CB 3 3, post-caldera block-and-ash flow deposits; CB 3 4, sinter deposits.  Table S1). The widely-used, advection-diffusion type, tephra fall simulation model, Tephra2 Connor & Connor, 2006), was used to model the paroxysmal phase of 4.2 ka cal BP eruption of CBVC. This deterministic modelling code allowed simulating affected area caused by CB 2 ash-fall. Tephra2 required three input datasets: a wind dataset, a grid dataset and a configuration dataset. The inverse problem method was used to set the input variables. Wind data (1980-2012; n=12,054 days) were downloaded from NCEP/NCAR Reanalysis site for the coordinates of CBVC (26 o 47' 24.46" S,67 o 45' 40.68" W). Different wind profiles were tested to be compatible with observations in the field (geographical distribution and ash layer thickness). The wind velocity profile retained was a prevailing wind direction towards 100 o E (Fernandez-Turiel et al., 2019; Table S2). The grid data were introduced in the model, taking as reference the position of the volcano. Thus, the total area was a rectangle centred in the volcano (coordinates 0, 0) and defined by 100 km to W, 1000 km to E, and 300 km to N and S, respectively. The elevation should be equal for all points on the grid. An elevation of 1000 m was used as reference due to the large differences in elevation in the affected region. Grids of 5×5 km and 10×10 km were checked and the results were very similar. To reduce computing time the 10×10 km grid was retained for modelling. The configuration data for Tephra2 simulation are in Fernandez-Turiel et al., 2019; Table S3.

Stratigraphy and description of deposits
The Cueros de Purulla sequence The Cueros de Purulla (CdP) sequence consists of the lower CdP 1 and upper CdP 2 units ( Table 1). The white CdP 1 unit mainly consists of fallout deposits, preserved in alluvial terrace and fan deposits that are elevated and isolated from active erosion channels (Figs. 3a and 3b). The CdP 1 unit is exposed up to several kilometres eastward from Cueros de Purulla volcano, with remnants recognised as far as 180 km from the probable vent. Proximal deposits consist of a series of 1-10 cm thick beds of lapilli and ash which together make up a total thickness between 20 and 70 cm. Some layers display normal grading of subangular pumice clasts, 1 to 60 mm in size. The content of lithic fragments is very low and includes exotic clasts of regional rocks that are generally smaller than 1 cm. Distal deposits of CdP 1 unit were recognized in the Tolombón-Cafayate area (Figs. 3a,3b,and 4). The grain-size distribution of the ash is bimodal (Fig. 5a). The coarser and finer modes are 180 and 20 µm in the thin-bedded layers and 320 and 20 µm in the unstratified layers (Fig. 5a).
The unit CdP 2 includes two depositional sub-units (CdP 2 1 and CdP 2 2) that together filled 10 km of the valley located east of Cueros de Purulla volcano (Fig. 3c). The lower CdP 2 1 sub-unit is in direct contact with the CdP 1 unit. This contact is sharp and sometimes erosional. Internal stratification is lacking in both sub-units but reverse grading of pumice and normal grading of lithics occur occasionally. The pumice clasts are sub-rounded, ranging from 2 to 6 cm. The lithic fragments generally consist of angular to sub-angular country rocks ranging from 1 to 3 cm. These deposits show no evidence of welding or vesicle collapse.

The Cerro Blanco sequence
The Cerro Blanco sequence (CB) consists of at least seven sub-units, which we have grouped into three packages defined largely by observable unconformities (Fernandez-Turiel et al., 2013 (Table 1).
The lowest unit CB 1 is exposed around the outer northern margin of the Cerro Blanco Caldera (Fig. 2), and is made up of angular decimetric rhyolitic blocks immersed in a rhyolitic lapilli and ash matrix with diffuse metrescale stratification showing reverse grading. The deposits are clast-supported and some blocks exceed 1 m in size. There are not pumice clasts.
The CB 2 unit comprises three sub-units (Table 1). The lower CB 2 1 sub-unit mantled palaeotopography and is preserved as patches on hill slopes and above alluvial terrace deposits incised by active erosional channels (Fig. 3d). The thickness of these deposits ranges from 26 to 5 cm at distances of 11 to 200 km from the CBVC (Fig. 6). It consists in alternating, parallel, very thin beds of lapilli and ash (Fig. 7). The CB 2 2 sub-unit directly overlies CB 2 1. The contact is sharp and not erosional. It consists of an ash deposit without grading, with poorly defined stratification. CB 2 2 is locally overlain by alluvium or colluvium and elsewhere the surface has incipient soil developed on top (Figs. 7b, 7c, 7e, 7f and 7i).
The deposit thins from several meters near the vent to ~30 cm thickness at 170 km, but then thickens to > 300 cm thickness at ~200 km (Tafí del Valle area) before thinning once more at greater distances, e.g., ~20 cm thickness at 370 km from the CBVC near Santiago del Estero (Figs. 6 and 7). As a first approach, an isopach map was constructed taking into account the field measurement of the fallout deposit thickness of CB 2 2 (Fig. 8). To validate these isopachs, the calculated areas were examined in a semilog plot of the square root of the contour area versus isopach thickness (Fierstein & Hildreth, 1992;Fierstein & Nathenson, 1992;Hildreth & Fierstein, 2012a, 2012bPyle, 1989), testing three different models for tephra deposition (exponential, power-law, and Weibull) using the AshCalc code (Daggitt et al., 2014). The lower mean relative squared error was obtained with the Weibull model ( Fig. 8) (Fernandez-Turiel et al., 2019; Table S4).
These numerical approaches allow us to have a rough estimation of the volume of the Plinian deposits, despite the uncertainties owing to different features, e.g., thickness reduction by erosion of the upper part of deposits or error associated with techniques of estimating volume of fallout deposits by extrapolating thickness versus distance plots and assuming continuous decrease in thickness regardless of secondary thickening. The Weibull model provides a volume estimate of 172 km 3 , higher than the power-law model (164 km 3 ) and lower than the exponential model (194 km 3 ) (Fernandez-Turiel et al., 2019; Table S2).
The particle size distributions of CB 2 1 and CB 2 2 are unimodal in proximal deposits and bimodal in distal deposits (>100 km from source) (Fig. 5b). In bimodal distributions, the lower size mode (~20 µm) does not change with distance (Fig. 5c). Instead, the upper size mode  shows a thinning from 600-700 to 100 µm particle size with distance up to ~200 km (Fig. 5c). Figure 5d shows the abundance of coarse and fine populations against distance for CB 2 1 and CB 2 2 sub-units. The coarse mode is less abundant with distance, but we cannot properly distinguish changes in volume of particles of finest mode (~20 µm) with distance.
The CB 2 3 sub-unit filled the valleys north, northwest, and south of the CBVC and can be followed for more than 35 km. It directly overlies sub-unit CB 2 1 in outcrops of The particle size bimodality observed in primary distal CB and CdP deposits is interpreted as evidence of particle aggregation in the eruptive plume, forming cored ash clusters that typically break on impact with the ground. The resulting fallout deposit is made up of core particles of hundreds of microns in size and shell particles of tens of microns in diameter that covered the former. (a) Typical particle size distributions of volcanic ash of CdP 1 unit of Cueros de Purulla sequence, and (b) CB 2 1 and CB 2 2 subunits of Cerro Blanco sequence. The particle size bimodality of Plinian ash-fall deposits is observed at >100 km from the CBVC: (c) Modes of unimodal and bimodal (coarse and fine) particle size distributions against distance for CB 2 1 and CB 2 2 sub-units of Cerro Blanco sequence; and (d) abundance of coarse and fine populations against distance for CB 2 1 and CB 2 2 sub-units of Cerro Blanco sequence.   metamorphic and volcanic rocks, 1-20 cm in diameter, are abundant at the basis of the unit in many places. Parallel and cross bedded, 0.2 to 2 m thick beds occurs particularly on palaeotopographic highs. The unit CB 3 consists of four sub-units (Table 1), all of them exposed close or within the caldera of Cerro Blanco. The exceptional preservation of these sub-units is mostly due to the arid environment of the region (Aulinas et al., 2015). The CB 3 1 sub-unit comprises crystal poor, very vesicular, rhyolite lava domes. These domes were described in detail by Báez et al. (2016). The largest dome appears in a rectangular shape of 2.7 × 1.4 km, and cut the SW caldera margin. In addition, three other smaller, pinkish, extra-caldera domes are aligned along a trend N100°E that includes the above-mentioned larger dome. The largest dome cut the caldera margin, indicating their post-caldera emplacement and a change in the eruptive style from highly explosive to effusive-moderately explosive.
The CB 3 2 sub-unit is made up of 3-30 cm thick, white lapilli and ash deposits organized in poorly to welldefined layers (Fig. 2a), forming the highest peaks of the CBVC, surpassing 4700 m a.s.l.; lava blocks are common within the deposit and, in some cases, show  a bread crust structure. The thickness reaches 150 m. These deposits are around the main lava dome, except on the side oriented towards the caldera centre where the lava dome is directly in contact with CB 3 3. These deposits are interpreted as fallout from eruptions associated with these structures (Báez et al., 2016). Within the unit also phreatomagmatic deposits with lapilli and ash can be observed, characterized by low-angle cross stratification. The deposits of CB 3 3 form several overlapping lobes reaching ~3.5 km from the main lava dome to the opposite wall of the caldera (Fig. 2), almost completely covering the bottom of the Cerro Blanco Caldera. They consist of 10-100 cm thick beds, diffusely stratified, and made up of decimetric angular blocks of pumice immersed in a coarse lapilli matrix. The total thickness reaches up to 50 m.
The Cerro Blanco sequence ends with 1-5 m thick deposits of siliceous sinter in the middle of the caldera (CB 3 4).

The Bolsón de Fiambalá sequence
A third sequence is defined on basis of the geomorphological, stratigraphical, mineralogical, and chemical characteristics of numerous outcrops occurring in the northwest of Bolsón de Fiambalá (BdF), 40-60 km to the south of CBVC (Fig. 1). We combine these deposits in only one lithostratigraphic unit called BdF 1 (Table 1). It mantles the palaeotopography with a series of unstratified and diffusely parallel medium bedded layers which together make a total thickness between 0.8 and 1.6 m. Some layers are pumice lapilli and others are largely ash. Occasionally, layering is highlighted by alignments of pumice of ~1 cm in size. BdF 1 is locally overlain by alluvium or hillslope colluvium, but also by aeolian deposits of lapilli and ash.

Compositional characteristics
The Cueros de Purulla sequence The ash of the CdP 1 unit in the Cueros de Purulla sequence is made up of glass shards and micropumices, with angular to sub-angular blocky shapes containing tubular vesicles (Fernandez-Turiel et al., 2019; Fig. S1). The dominant phenocrysts are biotite, which is the distinguishing feature of these deposits in the field, feldspars (plagioclase >> K-feldspar), quartz, titanomagnetite, and ilmenite. Amphiboles and clinopyroxenes are minor contents. Apatite, allanite-epidote, muscovite, titanite and zircon can be found in trace amounts.
The ash-dominant matrix in the CdP 2 unit is fine to coarse in size and consists of pumiceous glass shards, and very abundant feldspar (plagioclase >> K-feldspar), biotite, and quartz crystal fragments (Fernandez-Turiel et al.,  2019; Fig. S1).
The rare clinopyroxene found in the CdP 1 unit (n = 5) are diopsides (Fig. 14). Amphibole crystals are also rare in the Cueros de Purulla sequence and were only observed in CdP 1 unit (Fig. 14). Compositions include magnesiohornblende and pargasite, according the IMA 2012 recommendations (Locock, 2014). These samples showed values of Al# exceeding 0.21, the threshold at which the amphiboles are inferred to represent xenocrysts of crustal or mantle materials (Ridolfi et al., 2010).
The lava domes and the deposits of CB 3 2 and CB 3 3 are porphyritic and very vesicular; their glassy matrices include quartz and alkali feldspars as the main phenocrysts, and biotite and Fe-Ti oxides as minor phases. There are also rare apatite, allanite-epidote, and zircon phenocrysts. A detailed description of lava domes was perfomed by Báez et al., 2016. The glass composition of all CB units is rhyolitic (74.6-80.2 wt% SiO 2 , normalized to anhydrous contents) with alkali concentrations of 3.1-5.4 wt% K 2 O (Fig. 10), 2.1-6.2 wt% Na 2 O, and CaO < 0.9 wt%.
In the CB sequence, K-feldspar is more abundant than plagioclase and, while two groups of plagioclases are observed, only one group of K-feldspars can be distinguished (Fig. 11). Both groups of plagioclases coexist in the same sample. In one group of plagioclases, the An content in different crystals decreases from 44.3 to 28.5 mol %, with an increase in Or content from 2.3 to 8.1 mol % (n = 16). On the other hand, the second group shows An content in different crystals decreasing from 19.5 to 10.2 mol %, with an increase in Or content from 4.5 to 7.2 mol % (n = 138). The Or content in different crystals of K-feldspar ranges from 61.6 to 72.5 mol %, with an An content of < 1.3 mol % (n = 259).
The TiO 2 -MnO diagram of biotite provides a good chemical discrimination between two groups of biotites of the CB sequence: a high-MnO (1.11 ±0.06 wt%) and a low-TiO 2 (3.5 ±0.2 wt%) group and a low-MnO (0.21 ±0.07 wt%) and high-TiO 2 (4.7 ±0.7 wt%) group (Fig.  12). Both groups can be observed in the same sample.  The rest of the major elements show similar concentrations in both groups. The high-MnO group is relatively compact and includes samples from all sub-units of the CB sequence. However, the low-MnO group of biotites occurs only in CB 2 . The CB sequence contains two distinct titanomagnetite and ilmenite populations. Examples of both groups coexist in the same sample and exsolution features are not observed within the grains. One group is compositionally compact, showing minimal variability among CB sub-units with a TiO 2 content 2.6-3.7 wt% and MnO 2.4-3.2 wt% in titanomagnetite, and 41.3-45.4 wt% TiO 2 , and 4.7-7.9%, m/m MnO in ilmenite (Fig. 13). The other group was found only in the three sub-units of the CB 2 unit, i.e., those related to the paroxysmal phase, and it contains lower concentrations of MnO in both titanomagnetite (0.3-1.9 wt% MnO) and ilmenite (0.3-0.7 wt% MnO), and wide ranges of TiO 2 in titanomagnetite (1.5-8.2 wt%) and in ilmenite (0.3-0.7 wt%). The abundance ratio of titanomagnetite/ilmenite is ~50/50 in the high-Mn group, while the titanomagnetite dominates ilmenite in the low-Mn group.

The Cueros de Purulla sequence
Feldspar data were used to estimate the temperature. Three pairs of alkali feldspar-plagioclase crystals from CdP 1 , using the two-feldspar thermometer (eqn. 27b, Putirka 2008), yield temperatures between 780 and 800 °C for 5 kb. TiO 2 (wt%) Coexisting titanomagnetite-ilmenite pairs were evidenced in three samples, one each from CdP 1 , CdP 2 1, and CdP 2 2. Fe-Ti exchange temperatures, calculated using the Fe-Ti two oxide geothermometer and oxygen-barometer (Ghiorso & Evans, 2008), fall in the range of 750 to 800 °C. Results are also indicative of oxidizing conditions (log fO 2 relative to NNO between 1.2 and 1.4) and elevated TiO 2 activity (a TiO2 ~0.8) in the magma just prior to eruptive quenching.  The diopsides of CdP 1 unit allowed the use of the equilibrium pairs of clinopyroxenes and melt to determine both temperature and pressure conditions (eqn. 35 for T and eqn. 30 for P, Putirka, 2008). For melt composition, we used average glass composition of samples with clinopyroxene. Mineral-melt pairs were checked to ensure they satisfy the equilibrium condition K D (Mg-Fe) cpx-liq = 0.28 ±0.08 (eqn. 35, Putirka, 2008) (Gilbert et al., 2014;Schindlbeck et al., 2014). Clinopyroxene-liquid pairs satisfying the before-mentioned equilibrium condition (n = 5) yielded a temperature of 795-828 °C and a pressure of 3.7-8.3 kb with 0.21 < K D (Mg-Fe) cpx-liq < 0.23.
Fe-Ti oxides in CB 1 (n=1 pair) yielded an equilibration temperature of ~590 °C and a log fO 2 value of 1.2 (relative to NNO). CB 2 pairs (n = 10) provided from 580 to 725 °C and log fO 2 of 0.4-0.9. Only one pair of low-MnO titanomagnetite-ilmenite was found and it showed 794 °C and log fO 2 of 1.3. Finally, CB 3 3 pairs (n = 3) indicated ~610-630 °C and log fO 2 of 0.7-0.8. Physical-chemical conditions during amphibole formation were estimated by applying the thermobarometric method based on Al contents resulting in temperatures from 801 ±22 to 1009 ±22 °C, pressures from 108 ±12 to 554 ±61 MPa, water contents of the melt from 3.8 ±0.4 to 7.6 ±0.4 wt%, and log f (O 2 ) from -8.8 ±0.4 to -12.8 ±0.4. The reliability of this thermobarometer was discussed (Erdmann et al., 2014;Ridolfi & Renzulli, 2012), showing incorrect pressure values but reasonable estimations of temperature and oxygen fugacity (Pignatelli et al., 2016). These higher temperatures in comparison with feldspars and Fe-Ti oxides possibly indicate that these crystals are xenocrysts.

Eruption modelling of CB2 deposits
Tephra2 simulation results for CB 2 deposits indicates that tephra fall deposit was mainly distributed towards ESE (Fernandez-Turiel et al., 2019; Fig. S3). These results are compatible with the geographical ellipticity of the CB 2 deposits (Fig. 8), caused by the paroxysmal phase of 4.2 ka cal BP eruption of CBVC, strongly suggesting windy atmospheric conditions at the time of the eruption with a general dispersion towards the ESE. The simulation showed 0.5 m isopach contour reached as far as 400 km from the source (near Santiago del Estero). The preserved thickness in this area reaches 0.2-0.4 m (Figs. 6 and 7 i). The estimated plume height of the eruption was 27 km. The plume ratio of 0.9 indicated that pyroclastics were released from the upper 10% of the eruption column.

Volcanic meaning of studied sequences
The Cueros de Purulla sequence The CdP 1 unit is interpreted as air-fallout deposits from the buoyant eruption plume of a highly explosive eruption (Cas & Wright, 1988). The mineralogical and chemical correlation of proximal and distal deposits showed that ash-fall reached more than 180 km eastward Cueros de Purulla volcano (CdP). Reworking is discarded as the origin of grain-size bimodality, because this feature is observed in primary deposits. This bimodality is interpreted as evidence of aggregation of particles in the eruption cloud (Bagheri et al., 2016;Brazier et al., 1983;Brown et al., 2012;Sorem, 1982). The CdP 2 deposits are interpreted as ignimbrites emplaced in successive eruptive pulses.

The Cerro Blanco sequence
The Cerro Blanco sequence (CB) was interpreted largely by observable unconformities that developed in response to caldera collapse during the evolution of the Cerro Blanco Volcanic Complex (CBVC) (Báez et al., 2015(Báez et al., , 2016Fernandez-Turiel et al., 2013: the pre-caldera CB 1 unit, the syn-caldera CB 2 unit, and the post-caldera CB 3 unit ( Table 1).
Deposits of the lowest unit CB 1 were interpreted as produced by block-and-ash flows from pre-caldera dome structures (Báez et al., 2015(Báez et al., , 2016. The mantled topography, the poorly defined stratification or parallel thin-bedding as well as the particle size (lapilli and mostly ash) of CB 2 1 and CB 2 2 deposits are consistent with air-fallout deposits from the buoyant eruption plume of a highly explosive eruption (Bagheri et al., 2016;Cas & Wright, 1988). The lack of unconformities between CB 2 1 and CB 2 2 indicates no time gap between the emplacements of the deposits. The angular to sub-angular morphology of the pumice clasts and glass shards indicates a low degree of abrasion during transportation in the eruptive plume. The grain-size bimodality observed in primary distal CB 2 1 and CB 2 2 deposits (>100 km from source) is interpreted as evidence of particle aggregation and sedimentation of cored ash clusters. They are mostly sub-spherical fragile aggregates, observed with high-speed, high-resolution videos, that have never been observed in the deposits as they typically break at impact with the ground. They consist of a core particle (200-500 µm) fully covered by a thick shell of particles < 90 µm (Bagheri et al., 2016;Brown et al., 2012). The breaking of cored ash clusters and sedimentation of the integrating particles explain the bimodality in the distal CB 2 1 and CB 2 2 deposits, where the core particles are 100-400 µm in diameter and the shell particles are ~20 µm in diameter (Fig. 5c).
The CB 2 3 sub-unit is interpreted as an ignimbrite massive lapilli-ash (mLT) lithofacies with inner transitions to diffuse-stratified lapilli-ash (dsLT) lithofacies, and changes to parallel and cross bedded lapilli-ash (//sLT and xsLT) lithofacies towards the top and on topographic highs (Branney & Kokelaar, 2002) as well as in distal regions. Lateral and vertical variations are related to changes in the pyroclastic flow dynamics and their interactions with topography. The poorly defined stratification indicates pyroclastic flows under high depositional rates, transporting high concentrations of pumice and ash (Branney & Kokelaar, 2002;Cas et al., 2011). The pumice imbrication shows the shear sense within the flow at that location. The interaction of these flows with the valley bottoms produced the ingestion of metamorphic and volcanic clasts from alluvial beds and terraces and incorporation into the bottom of ignimbrite. The change from a poorly defined stratification to parallel and cross-bedded structures evidences the action of traction-dominated flows, especially onto palaeohills. The last deposits were interpreted as stratified topographic veneers (Báez, 2014;Báez et al., 2015).
Although there are no geochronological dates, the stratigraphy indicates that the paroxysmal phase of the eruption (CB 2 unit deposits) was followed by the emplacement of various lava domes (CB 3 1 sub-unit) with associated fallout and phreatomagmatic deposits (CB 3 2 sub-unit) (Báez et al., 2016). The dispersion direction of the lobes, the structure of deposits, the absence of pumice and exotic lithic clasts, and the same rhyolitic composition shared by all components, all indicate an origin of CB 3 3 sub-unit related to successive block-and-ash flows generated by gravitational or explosive lava-dome collapses (Charbonnier & Gertisser, 2008).
Finally, the siliceous sinters in the middle of the caldera (CB 3 4) are interpreted as result of amorphous silica deposited around hot springs.

The Bolsón de Fiambalá sequence
The BdF 1 unit is interpreted as fallout deposits from the buoyant eruption plume of a highly explosive eruption (Cas & Wright, 1988). The overlying lapilli and ash deposits with heterogeneous clasts of regional rocks indicate reworking by water or wind (Figs. 3e and 4). Aeolian pumice-rich dunes are very common in the central and western areas of Bolsón de Fiambalá. Between BdF 1 unit and these overlying sediments there is evidence of a time gap, recorded by an erosional contact and occasional alluvial deposits.
The Nevado Tres Cruces Volcanic Complex, 90 km W of the Bolsón de Fiambalá basin, is the strongest candidate for the source of the BdF 1 deposits. Although proximal samples are not available to test this hypothesis, the geographical distribution, the geomorphological features observed in satellite images, and the petrological information from previous works (Baker et al., 1987;Gardeweg et al., 2000) describing the main trends of the stratigraphy and mineralogy (abundance of hornblende and biotite, and absence of K-feldspars) in the younger proximal dacitic ash-fall deposits, ignimbrites, and lava-domes of the Nevado Tres Cruces complex, favours this edifice as the possible source.

Implications of tephrochronology
Our results on tephra deposits of Southern Puna and neighbouring areas in NW Argentina are all consistent with three large eruptive events during the Holocene in the southern Puna and adjacent regions. Compositional data of glass and minerals  and the unique occurrence of abundant large biotite crystals visible to the naked eye assist the association of the fallout deposits of Tolombón and Cafayate (CdP 1 unit) with the published Cerro Paranilla Ash of the Calchaquíes Valleys (Hermanns et al., 2000), whose age was established close to 7820 ±830 a BP. This age was obtained by 10 Be dating of a breakaway scarp of a rock-avalanche landslide directly overlying the ash (Hermanns & Schellenberger, 2008).
We found the fallout ash deposits of Cerro Blanco Volcanic Complex (CB 2 1 and CB 2 2) overlying a peat with an age of 4880 to 4780 a cal BP in Tafí del Valle). This age is consistent with the age of 4510 to 4240 a cal BP of a charred coal in a palaeosol underlying the ignimbrite (CB 2 3) in Las Papas, and with the age of 4410 to 4100 a cal BP of charred material within this ignimbrite in Laguna Aguada Alumbrera (Figs. 7h, 9a, 9c). The fallout ash deposits of CBVC are correlated by stratigraphic relationship, chemical composition, mineralogy, and geochronology with previously described ash layers in the region: the mid-Holocene ash (Frenguelli, 1936), Ash C (Malamud et al., 1996), volcanic ash layer (Wayne, 1999), Buey Muerto Ash (Hermanns & Schellenberger, 2008), and V1 ash layer (Sampietro-Vattuone & Peña-Monné, 2016) in the Calchaquíes Valleys and Tafí del Valle area, and with the volcanic ash layer found in the archaeological site of Cueva Salamanca I near Antofagasta de la Sierra (Pintar, 2014). Furthermore, the CB 2 3 ignimbrite dating is coherent with an age of an ignimbrite associated with CBVC (Montero-López et al., 2010a). Combining the published ages with our results, the age of the eruptive event is likely closer to 4.2 ka cal BP.
The most recent fallout ash deposits of BdF 1 were identified as the V2 ash layer in the Calchaquíes Valleys and Tafí del Valle area (Sampietro Vattuone et al., 2018;Sampietro-Vattuone & Peña-Monné, 2016). One fragment of a ceramic pot underlying the ash layer was TL (thermoluminiscence) dated 1772 ±112 a cal BP in Tafí del Valle (Sampietro-Vattuone & Peña-Monné, 2016). This age is coherent with those obtained from a charcoal (1520-1320 a cal BP) and a burnt basket (1400-1270 a cal BP), which were found beneath a reworked, 40 cm thick, pumiceous deposit of BdF 1 in the archaeological site of Palo Blanco (27°20'15.6"S, 67°44'36.4"W) (Fig. 3f) (Ratto, 2013). This deposit was identified by the amphibole phenocrysts. Finally, this ash layer was temporally located ca. 1500 AD in the Calchaquíes Valleys (Sampietro Vattuone et al., 2018), but dating was done in a soil profile different to the profile where the V2 ash was observed. On the other hand, the oldest chronicle known in the area, related to the expedition of Diego de Almagro to Chile in 1536, does not relate any eruption during expedition or in previous years in the region (Ramón Folch, 1954). In our opinion, more work is needed to precise the age of this eruption.

Eruptive style of studied sequences
The three eruptions studied in this work were highly explosive and generated voluminous ash-fall deposits, non-welded ignimbrites, and lava dome extrusions. The extensive area of the voluminous fallout deposits of the CdP, CB and BdF sequences points towards a Plinian eruption style related to a great column height. We are in a preliminary phase of the study of the CdP and BdF sequences, which does not allow going further into describing the involved eruptions. Instead, the results from the CB sequence have enabled the identification of the most relevant features of the 4.2 ka cal BP Cerro Blanco eruption. Our findings corroborate previous results on the last eruptive event of the CBVC (Báez et al., 2015;Fernandez-Turiel et al., 2013, distinguishing pre-, syn-, and post-caldera products, and allow us to estimate the large extent of the Plinian fallout during the paroxysmal phase of the eruption.
The pre-caldera phase of the Cerro Blanco eruption The post-depositional dissection of the blockand-ash flow deposits of CB 1 by the caldera edge allows the interpretation that it is a pre-caldera unit from an extrusive dome that is compositionally linked to CB. Therefore, this unit is included in the CB sequence assisted by its similar petrography and geochemistry to the other CB units that differentiate them from the older products of the CBVC (Báez et al., 2015;Montero-López et al., 2010b).

The syn-caldera phase of the Cerro Blanco eruption
The lack of geochronological data does not allow estimation of the time gap between CB 1 and CB 2 , i.e., the caldera collapse. The absence of lithic breccias or pyroclastic products around the caldera edges, as well as the low lithic content in pyroclastic products of CB 2 , suggest that a single vent in the SW of the caldera rim acted as emission conduit during the eruption, favouring optimal conditions for a sustained Plinian column and fallout deposition.
The fall deposit CB 2 1 of alternating layers, 1-3 cm thick, some of lapilli and some of ash are the first record of the Plinian column. The lower thickness of this deposit among CB 2 sub-units seems to be consistent with a relatively brief episode in comparison with those generating the overlying sub-units. Thin bedding suggests instability because of changing column heights and therefore variable vertical wind patterns (Cioni et al., 2000). The change to the thicker, unstratified rhyolitic ash (CB 2 2) indicates that the feeding system reached and maintained steady conditions during a single climactic column phase (Cioni et al., 2000;Pardo et al., 2012). The eruption column height for this phase reached 27 km, according the results of Tephra2 modelling, with an ESE dispersal on a large surface of the Chacopampean Plain (Fernandez-Turiel et al., 2019; Fig. S3).
Bimodal particle size distributions, as observed in CB 2 1 and CB 2 2, are a common feature in many distal ash deposits, such as in the eruptions of Campanian Ignimbrite (CI) 39 ka BP (Engwell et al., 2014), Pinatubo in 1991 (Wiesner et al., 2004), and Chaitén in 2008 Ruggieri et al., 2012;Watt et al., 2009). The occurrence of these particle size subpopulations in fallout deposits was explained by the propensity of aggregation in historic and ancient eruptions, independently of their style (Brown et al., 2012;Durant et al., 2009Durant et al., , 2012, but also by other processes without the need for aggregation, e.g., particle differences in texture and componentry, convective instabilities in eruptive plume, and lower atmospheric circulation patterns (Durant et al., 2009;Watt et al., 2015).
An interesting new approach to explain the particle size bimodality is through particulate cored clusters. They are sub-spherical fragile aggregates consisting of a core particle of hundreds of microns in size covered by a thick shell of particles of tens of microns in diameter (Bagheri et al., 2016). These aggregates are formed in the buoyant column and plume of the eruption. These ash clusters tend to disaggregate on impact with the ground and, as consequence of their very low preservation potential, they are difficult to find and document in the historic and ancient ash-fall deposits, although their presence can be inferred by polymodal particle size populations. In this way, the bimodality of the Plinian CB 2 1 and CB 2 2 sub-units from 100 km away from the Cerro Blanco through more than 350 km (Fig. 5c) could be due to the clustering phenomena of the ashes described in the literature.
The distal mass accumulation maximum of the Plinian CB 2 2 sub-unit reached several meters thick ~210 km from the volcano, in the eastern ranges that surrounds the southern Puna (Figs. 6,7e,7f, 2019; Fig. S4). Semi-arid conditions prevail west of this area, while warm, moist, and rainy conditions prevail over the easternmost slopes of these ranges and adjacent lowlands. The climate of Central Andes is characterized by a complex interplay of large-scale atmospheric circulation with local orographic effects (Garreaud et al., 2003;Garreaud, 2009;Neukom et al., 2015), and the main moisture source for the region is the easterly influx from the Amazon Basin. In this setting, the eastern upslope moist warm flow could accentuate the ice melting in aggregates, collapsing their structure and triggering ash-fall en masse.
The secondary thickening related to upslope moisture transport by easterly atmospheric flow, should also elucidate the presumable secondary thickenings of the Cueros de Purulla eruption near Cafayate, where there are several metre-thick fallout ash deposits of CdP 1 (Fig. 3b), and the BdF 1 ashes in the northwest of Bolsón de Fiambalá (Fernandez-Turiel et al., 2019; Fig. S4). However, the warm and moist easterlies hardly reach these areas implying that these deposits have not necessarily affected by wet disaggregation and that other process may be highly significant in driving the fallout of particles. If we observe the location of the secondary thickenings for CdP, CB and BdF fallout units in topographic profiles (Fernandez-Turiel et al., 2019; Fig. S4), the affected areas are in the eastern slopes of important ranges. This fact allows us to consider in the conceptual model the effect of the topographically induced turbulences in the disaggregation, e.g., the breaking of lee waves, generated by winds passing over elevated topography beneath the eruption plume (Watt et al., 2015). As a consequence, this spatial association of the secondary thickening and topography could help to locate affected areas in similar cases. It may also help to prevent hazards associated with the increased loading of ash on infrastructures and buildings (Wilson et al., 2014), because it allows us to constrain more precisely the most hazardous areas during large explosive eruptions around the world and, particularly, in the Central Volcanic Zone of Andes, which are the eastern slopes of the mountain ranges bordering the Altiplano-Puna plateau.
The CB 2 3 PDC's maybe were partially synchronous with the CB 2 2 ash-fall. Many Plinian columns generate regional fall deposits and large-volume ignimbrites at the same time. Marginal collapse forms PDC's while central column sends up regional fall (e.g., Bishop Tuff, 1912 Novarupta) (Hildreth & Fierstein, 2012a). The transition to non-optimal conditions for a buoyant column, possibly with a significant increase in discharge rates, may also produce large PDCs preceding caldera collapse (Legros, 2000).
Thus, on the whole, the CB 2 unit follows the idealised sequence of a caldera-forming eruption, displaying a Plinian fall deposit overlain by an ignimbrite (Druitt & Sparks, 1984;Legros et al., 2000). However, in the Cerro Blanco caldera this model is more complex, as is revealed by the two different compositions of plagioclases, biotite, and Fe-Ti oxides observed in the CB 2 syn-caldera sub-units which seem to indicate mixing of two magma batches during the eruption. The compositional homogeneity of CB 1 and CB 3 deposits suggests a single magma body. Within the paroxysm of the eruption, however, it is likely that another significant magma input occurred, recorded in the compositional bimodality of CB 2 deposits. The most frequent compositions in CB 2 are also the same as the most common observed in CB 1 and CB 3 (Figs. 10-14). A possible explanation for the minority compositions could be the mobilization of hotter magma, evidenced by geothermobarometry data, from a different part of the magmatic chamber. Similar mixing-derived bimodality was observed in titanomagnetite compositions of magmatic products from, for example, Mono Basin in California (Marcaida et al., 2014) and Usu volcano in Japan (Tomiya & Takahashi, 2005). There are no geochronological data to determine the time gap between the syn-caldera CB 2 and the post-caldera CB 3 units. The CB 3 unit represents a significant change in the eruption style of the CBVC, from an explosive phase (CB 2 ) to an essentially effusive magmatic phase characterised by the extrusion of the CB 3 1 lava domes in the SW edge of the Cerro Blanco caldera, and the occurrence of different fallout episodes of pyroclasts and block-and-ash flows related to them. Block-and-ash flows were topographically confined towards the caldera's interior, filling the bottom of this 5 km wide circular structure almost completely (Fig. 2a). These flows involved both single and multiple-collapse events such as those documented in the eruptions of Unzen in Japan, Colima in Mexico, Merapi in Indonesia, Arenal in Costa Rica, and Montserrat in Lesser Antilles (Charbonnier & Gertisser, 2008;Freundt et al., 2000). A post-eruption geothermal field produced white sinter deposits within the Cerro Blanco caldera. The associated hot springs, located in the central part of the caldera, are active in the present day . Finally, the Cerro Blanco caldera area showed a deformation pattern subsiding about 1-3 cm/year at least since 1992 (Pritchard & Simons, 2004, Brunori et al., 2013, Henderson & Pritchard, 2013, Yazdanparast et al., 2017.

Volcanic explosivity index of 4.2 cal BP CBVC eruption
Only proximal and medial fallout deposit thickness data, reaching 400 km eastward, are available for the Cerro Blanco eruption, which is common in many prehistoric deposits . Consequently, large uncertainties in volume estimates are expected. Despite this, the order of magnitude obtained for bulk volume estimates always exceeds 100 km 3 , i.e., the threshold value for a volcanic explosivity index or VEI of 7 (Fernandez-Turiel et al., 2013. Taking as a point of reference the density of the 2008 Chaitén rhyolitic ash deposits, that ranged from 1,000 to 1,250 kg·m -3 (Alfano et al., 2011;Watt et al., 2009), we have assumed a deposit density of 1,000 kg·m -3 excluding lithics. With this value and a denserock density of 2,300 kg·m -3 , the 172 km 3 of bulk volume amounts to a dense-rock equivalent (DRE) volume of 75 km 3 . Based on a model of the pre-eruptive topography, the volume of ignimbrite associated with Cerro Blanco accumulated 15 km 3 without lithics (Báez et al., 2015), corresponding to 8.5 km 3 DRE, assuming a bulk density of 1,300 kg·m -3 , and a denserock density of 2,300 kg·m -3 . Based on these results, the Cerro Blanco eruption is among the largest volcanic eruptions of the Holocene globally, similar to Kuril Lake (Kamchatka, Russia), Crater Lake (Oregon, USA), Tambora (Indonesia), Samalas (Indonesia), Aniakchak (Alaska, USA), Santorini (Greece), or Changbaishan (China-North Corea) (Brown et al., 2014;Crosweller et al., 2012;Lavigne et al., 2013;Vidal et al., 2015). It clearly exceeds the magnitude of the 1600 Huaynaputina eruption, the largest historical eruption known in the Central Volcanic Zone of the Andes (Global Volcanism Program, 2013b; Stern, 2004;Tilling, 2009).

Conclusions
Results confirm the existence of three major rhyolitic eruptive events during Holocene in the southern Puna and neighbouring areas, and provide a more complete picture of the style, frequency, distribution, and size of past explosive eruptions in the region. The Cerro Blanco Volcanic Complex (the focus of this work) is demonstrated as source of the 4.2 ka cal BP pyroclastic deposits, whereas Cueros de Purulla volcano is suggested as the origin of the Lower Holocene pyroclastic deposits, and the Nevado Tres Cruces Volcanic Complex is the strongest candidate for the source of the Upper Holocene pyroclastic deposits found in the Bolsón de Fiambalá basin.
Based on 1) field observations (stratigraphic and geomorphological relationships, structure, and thickness of deposits), 2) petrography and geochemistry of juvenile eruptive products, 3) 14 C ages, and 4) erupted volume, we conclude that the CBVC generated the largest documented eruption during the past five millennia in the Central Volcanic Zone of Andes and one of the largest around the world for this period. Cerro Blanco ash-fall deposits reached >400 km from the vent indicating that ~170 km 3 of ash were spread around 500,000 km 2 in Argentina, while pyroclastic flow deposits extended north, northwest, and south several tens of kilometres from the vent.
The implications of the findings of the present work reach far beyond having some chronostratigraphic markers. Further interdisciplinary research should be performed in order to draw general conclusions on these impacts in local environments and the disruptive consequences for local communities. This is invaluable not just for understanding how the system may have been affected over time, but also for evaluating volcanic hazards and risk mitigation measures related to potential future large explosive eruptions.