Quaternary volcanism in the Yura Monogenetic Field near Arequipa city, southern Peru

Arequipa (Peru) is an area where volcanic activity has been persistent during the Quaternary. Studies carried out in this area have highlighted the emplacement of ignimbrite deposits, large volcanic clusters and stratovolcanoes. Monogenetic volcanism is also present, although poorly explored and studied. Due to its location over an ignimbrite plain and poor state of preservation, the only identified monogenetic cone in the Arequipa basin was the Nicholson volcano, while other monogenetic centers remained unknown. This lack of information about the recent volcanism can lead to inadequate definition of scenarios in a hazard assessment in the region. The present study has investigated monogenetic volcanism in the northwestern edge of the Arequipa basin based on geological survey, geochronology and geochemical data. Here, we report for the first time five small volcanic centers such as Yura Viejo, Ccapua, Uyupampa, El Chiral and Patacocha, which together with the Nicholson volcano form the Yura Monogenetic Field. Stratigraphic considerations and new 40Ar/39Ar ages allow us to place the eruptive activity in the Middle–Upper Pleistocene (c. 195–54 ka). Phreatomagmatic, Strombolian and effusive eruptions characterize the monogenetic activity of the field. As a result of these eruptions, small scoria cones, maars, and lava flows/coulées were generated. The eruptive products show ubiquitous olivine phenocryst-rich (<10 vol%) set in a fine pilotaxitic groundmass, suggesting rapid ascent of basaltic magmas to the surface controlled by the tectonic setting. The analyzed rocks lie in a narrow range of basaltic-andesite composition (50.9–55.6 wt% SiO2) being the most mafic Pleistocene - Recent volcanic products identified in the Arequipa basin, along with the least differentiated magmas from the nearby Chachani volcanic cluster. This work shows how monogenetic volcanism can occur contemporaneous and closely spaced to larger volcanic clusters and active stratovolcanoes. We hope the information provided here will contribute to improve the risk management by highlighting the scenario of monogenetic eruptions that should be considered in the hazard assessment.

Arequipa (Peru) is an area where volcanic activity has been persistent during the Quaternary. Studies carried out in this area have highlighted the emplacement of ignimbrite deposits, large volcanic clusters and stratovolcanoes. Monogenetic volcanism is also present, although poorly explored and studied. Due to its location over an ignimbrite plain and poor state of preservation, the only identified monogenetic cone in the Arequipa basin was the Nicholson volcano, while other monogenetic centers remained unknown. This lack of information about the recent volcanism can lead to inadequate definition of scenarios in a hazard assessment in the region. The present study has investigated monogenetic volcanism in the northwestern edge of the Arequipa basin based on geological survey, geochronology and geochemical data. Here, we report for the first time five small volcanic centers such as Yura Viejo, Ccapua, Uyupampa, El Chiral and Patacocha, which together with the Nicholson volcano form the Yura Monogenetic Field. Stratigraphic considerations and new 40 Ar/ 39 Ar ages allow us to place the eruptive activity in the . Phreatomagmatic, Strombolian and effusive eruptions characterize the monogenetic activity of the field. As a result of these eruptions, small scoria cones, maars, and lava flows/coulées were generated. The eruptive products show ubiquitous olivine phenocrystrich (<10 vol%) set in a fine pilotaxitic groundmass, suggesting rapid ascent of basaltic magmas to the surface controlled by the tectonic setting. The analyzed rocks lie in a narrow range of basaltic-andesite composition (50.9-55.6 wt% SiO 2 ) being the most mafic Pleistocene -Recent volcanic products identified in the Arequipa basin, along with the least differentiated magmas from the nearby Chachani volcanic cluster. This work shows how monogenetic volcanism can occur contemporaneous and closely spaced to larger volcanic clusters and active stratovolcanoes. We hope the information provided here will contribute to improve the risk management by highlighting the scenario of monogenetic eruptions that should be considered in the hazard assessment.

Introduction
Monogenetic volcanic fields are important manifestations of volcanic activity that usually consist of a set of separate vents (volcanic emission centers) including several volcanic geoforms such as scoria cones, maars, tuff cones, tuff rings, lava domes, and lava flows (e.g., Németh, 2010;Németh and Kereszturi, 2015). The term 'monogenetic volcano' is used to designate a volcano that erupts during a single period of time (weeks, months or years) in which there is no evidence of a significant or longlasting temporary interruption in eruptive activity (Macdonald, 1972;Williams and McBirney, 1979;. It is important to note that even within a day the accumulating tephra may suffer some syn-eruptive erosion and remobilization of pyroclasts. This is especially valid in areas where rainfall or heavy wind action may coincide with the eruption. The generally small volume of magma emitted is produced by one short, sustained period of activity, then the volcano becomes extinct, and any further activity occurs through a different magma intrusion and new locations. Usually, in monogenetic volcanoes, the magma body rises rapidly (within hours to days) to the surface, so there is little interaction with crustal rocks; it is for this reason that most of these volcanoes emit basaltic magmas (Valentine and Gregg, 2008).
Monogenetic volcanism in the Arequipa region has been important in terms of volume and extension, as is the case of the Andahua (composed of at least 22 cones), Huambo and Orcopampa monogenetic fields (e.g., Delacour et al., 2007, their Figure 1; Galas, 2011) compared to other monogenetic fields in southern Peru. In the Arequipa tectonic depression, only the Nicholson monogenetic volcano was previously known to exist (Delacour et al., 2007;Global Volcanism Program, 2013). Emplaced over the Arequipa Airport Ignimbrite (Paquereau-Lebti et al., 2006), this scoria cone is easily recognizable in aerial photographs or satellite images (Figures 1, 5). Tephra-fall deposits cropping out near the Uyupampa town were mapped as part of older deposits of the Chachani volcanic cluster (Aguilar et al., 2022) since there was not an identified source and because these tephra layers were covered by lava flows from one of the Chachani's edifices.
In this work, we present an overview of the Quaternary monogenetic volcanism near the Arequipa basin. We herein document for the first time the Yura Viejo, Ccapua, Uyupampa, El Chiral and Patacocha volcanoes, which together with the Nicholson cone may form what we name the "Yura Monogenetic Field". These small volcanoes composing the monogenetic field, the Chachani volcanic cluster and Misti stratocone border the Arequipa basin following a NW-SE trend ( Figure 1B). Based on detailed fieldwork, satellite images and photogrammetry, whole-rock geochemistry and petrographic descriptions we characterize Nicholson, Ccapua, Yura Viejo, and Uyupampa volcanoes. Due to inaccessibility, we carried out a remotely-based analysis of El Chiral and Patacocha volcanic centers. Furthermore, combining field data and morphometric parameters we propose the eruption style and classification of each volcano in this monogenetic field, contributing to the understanding of monogenetic volcanism in southern Peru.
The Yura Monogenetic Field, located along the border of a tectonic basin, where the small monogenetic volcanoes were emplaced over ignimbrite deposits and near the Chachani large volcanic cluster and Misti stratovolcano, constitutes a very interesting example of monogenetic activity occurring in a temporal and spatial proximity to polygenetic volcanism with larger explosive and effusive eruptions. Reconstructing the history of volcanic activity around urban areas is fundamental for constraining the hazard scenarios and risk assessment.
The geodynamic context includes an active continental margin and an archetypal Andean subduction (Thorpe et al., 1982;Wilson, 1986). However, in southern Peru, the continental magmatic arc is also obliquely convergent (i.e., the plate margin is at about N74 ± 4°E angle to the subduction direction; Norabuena et al., 1998). This oblique convergence is accommodated since at least the Neogene by sinistral faults (Mering et al., 1996;Sempere et al., 2014). Neogene to Recent volcanism is often associated with trans-tensional regimes because these areas offer less resistance through normal and strike-slip faults, through which the magma can rise more easily to the surface (e.g., van Wyk de Vries and Merle 1998;Acocella et al., 1999).
The WNW-ESE oriented Arequipa depression ( Figure 1B) has been interpreted as a pull-apart basin associated with regional N130°transtensive faults (Mering et al., 1996;Thouret et al., 2001;Benavente et al., 2017). The basin is filled by Middle Miocene to Early Pleistocene ignimbrites (e.g., Paquereau-Lebti et al., 2006), and more recent volcanoclastic deposits and alluvial deposits. The basin boundary is formed by the Cincha-Lluta-Incapuquio Fault System ( Figure 1B) which has experienced different tectonic regimes (Vicente et al., 1982), delineating the western limit of the Western Cordillera, and has influenced the sedimentation of the Arequipa basin (e.g., Carlotto et al., 2009). Studies conducted around the Arequipa basin have shown that this fault system has operated during the Mesozoic in successive extensional and reverse regimes (e.g., Acosta et al., 2010). The system was interpreted as transpresive faults (Vicente et al., 1982), which now act in a trans-tensional regime.
The N130°oriented Ayo-Lluta-Arequipa Fault (also known as Huanca fault) is a regional structure that was formerly a thrust fault. This structure brings into contact Mesozoic sedimentary units (including the Yura Group) with the Cenozoic volcanosedimentary units (Figure 1). During the Cenozoic, the area was affected by major compressive phases. Since the Middle Miocene, extensional tectonics has had a major influence on volcanism in the Western Cordillera (Huaman Rodrigo et al., 1985in Mering et al., 1996. Since the Pliocene, the Ayo-Lluta-Arequipa system has acted as a normal and sinistral fault (Mering et al., 1996;Sempere and Jacay, 2006).

The bedrock stratigraphy of Yura Monogenetic Field
Several of the pyroclastic deposits in the Yura Monogenetic Field contain rock fragments extracted from the substrate on which the volcanoes were emplaced. In the area surrounding the Yura town, the bedrock is formed of Mesozoic sedimentary and Neogene volcanic deposits (Figure 1). A summary of stratigraphic units of the bedrock is presented in Table 1.

Geological mapping and sample collection
Fieldwork was carried out during short field campaigns between June 2018 and November 2021. We visited 114 field sites distributed in the monogenetic field along the Yura Valley ( Figure 1). The source of individual lava flows and pyroclastic deposits was distinguished by using stratigraphic and volcanological relationships between deposits and geoforms (scoria cones, maars, fissures that fed flows or topographic highs). Mapping of eruptive units (1:25000 scale) included lithologic criteria (type, size, and abundance of mineral phases) and morphologic criteria (e.g., Condit and Connor, 1996). Outcrops in natural surfaces together with road-cuts and quarries served to carry out the facies observations and descriptions. Four stratigraphic logs were measured in outcrops with greater exposure of stratigraphic units. The geological maps of the Yura Monogenetic Field were digitized using ArcMap ® 10.1 (https://desktop.arcgis.com) at 1:25000 scale, the same scale used by INGEMMET (Instituto Geológico Minero y Metalúrgico) for geological maps of active volcanoes in Peru.
The coordinate system of all maps presented herein are UTM WGS84, Zone 19S. Polygons delimiting the volcanoes were used to calculate their surface areas in ArcMap. In some areas, highresolution digital elevation models (UTM-WGS84 DEM, 60-cm spatial resolution) were obtained by oblique photogrammetry with a Phantom-4PRO drone. The methodology of photogrammetry was similar to that described by Diefenbach et al. (2012) and James et al. (2020). We used AgiSoft Photoscan ® software (https://www.agisoft.com) for photogrammetric analysis. Morphometric parameters described by Wood (1979) such as cone basal diameter, crater diameter and cone height were obtained from fieldwork data and high-resolution DEMs Frontiers in Earth Science frontiersin.org 04 using ArcMap. The average flank slopes were calculated trigonometrically using the morphometric parameters (e.g., Hasenaka and Carmichael, 1985).

Petrography, geochemistry and geochronology
Twelve lava and scoria samples were used for petrographic descriptions. Thin sections of those samples were studied in order to describe representative rocks from the volcanic centers. Ten of them have been taken in volcanic deposits (bombs and lavas), while two are xenoliths included in tephra-fall deposits. Given that the mineral assemblage is quite homogeneous, the number of samples is deemed representative. All analyzed samples are fresh, or at least show no conspicuous weathering. Modal analyses of the sample set are given in Table 2. The abundance of each main mineral phase in thin section was determined by using images taken with a polarized microscope in the most representative zones. Using Adobe ® Illustrator (https://www.adobe.com/la/products/illustrator.html), the images were examined to determine the percentage of each mineral. For determining the frequency (modal analysis), we used the point-counting method (e.g., Rodríguez Rey et al., 2004) following five steps: 1) we consider a rectangular frame inside the image to restrict the crystal counting to this area. Two sides of the frame have dashed lines and the other have continue lines. 2) All crystals that were completely included within the frame and the crystals that have some portion inside the dashed lines were counted (point counting). 3) Into the rectangular frame, we added a grid of 600 points. 4) Each point intersecting a crystal was manually counted to calculate the number of crystals contained in the frame. 5) The % abundance is obtained by: (# of crystals/Area of rectangular frame) × 100.
Major and trace element analyses of whole-rock samples were performed by ICP-AES and ICP-MS respectively, with an analytical precision (2σ) less than 1% for major elements and around 5% for trace elements at Instituto Geológico Minero y Metalúrgico (INGEMMET) in Peru. For these analyses, agate-crushed powders of seventeen samples were mixed with LiBO 2 , placed in a graphite crucible, and melted in an induction oven at 1050°C for 4.5 min, resulting in a homogeneous glass bead. The glass was then

Unit
Description Age Reference

Chachani Volcano Cluster
Consist of at least twelve volcanic edifices forming a group of stratovolcanoes, dome coulées and blocky-lava flow fields. Their deposits overly the Arequipa Airport Ignimbrite and the Yura Tuff, and a lava flow in its western flank partially covers the Ccapua monogenetic volcano.
1012 ± 53-56.5 ± 31.6 ka Aguilar et al. (2022) Yura Tuff Non-welded pumice and ash pyroclastic flows the volume of which was estimated at~1.5 km 3 . Yura Tuff outcrops are restricted to the north and west sides of Chachani complex. They fill a north-south elongated depression between sediments of the Yura Group and Pre-Chachani lava bedrock 1.278 ± 46 Ma Paquereau- Lebti et al. (2006Lebti et al. ( , 2008. Aguilar et al. (2022) Arequipa Airport Ignimbrite Consists of two units (18 km 3 ), a white indurated massive ash-andlapilli tuff overlain by a pink non-welded, massive and lithic-rich lapilli tuff. The source of this ignimbrite is buried under the Chachani Volcano Cluster, as indicated by Anisotropy of Magnetic Susceptibility measurements and the size of lithic fragments, which increases towards the north 1.65 ± 0.04 Ma Paquereau-Lebti et al. (2006 Tacaza Group Composed of volcanic breccia and andesitic lavas. These volcanic rocks crop out in the northwestern part of the study area Oligocene Jenks (1946) Arcurquina Fm Composed of marine, regularly-bedded, thickening-upward, grey to black, organic rich micritic limestones. These rocks are found to the northwestern side of the study area Middle-Upper Jurassic Benavides (1962) Chocolate Fm Composed of purplish volcanic and sedimentary (sandstone and limestone) rocks, some of which have been affected by metamorphism. The outcrops of Chocolate Fm. are found in the southern part of the area Lower Jurassic Vargas (1970) Frontiers in Earth Science frontiersin.org 05 dissolved in a solution of deionized water and nitric acid (HNO 3 ) and finally diluted by a factor of 2000. For comparison with surrounding volcanic products, these new data (Table 3) were combined with previously analysed rock samples from Chachani volcano cluster Aguilar et al., 2022), Arequipa Airport Ignimbrite and Yura Tuff (Paquereau-Lebti et al., 2006).
The chronology of volcanic activity in the Yura Monogenetic Field was constrained using the 40 Ar/ 39 Ar dating method at the Laboratoire des Sciences du Climat et de l'Environnement (LSCE/IPSL, Gif-sur-Yvette, France). Groundmass in the samples from Yura Viejo, Uyupampa, Ccapua and Nicholson volcanoes were used for age determinations. 40 Ar/ 39 Ar ages are summarized in Table 4, whereas, the corresponding analytical procedures and results are detailed in the Supplementary Material S1, S2.
4 Volcanism and crustal structure in the yura monogenetic field

Volcanic activity
The Yura Monogenetic Field is composed of six volcanic centers; three of them (Patacocha, Yura Viejo and Ccapua) are  Ages calculated relative Acs-2 at 1.1891 Ma (Niespolo et al., 2017) and the total 40 K decay constant of Renne et al. (2011). (1): Full external error. n/M: number of crystals retained in the age calculation over the number of analyzed crystals. F: Spreading factor of Jourdan et al. (2009). Values in bold indicate the ages considered for each volcano.

Frontiers in Earth Science
frontiersin.org distributed along~10 km in a N130°direction, and two volcanic centers (El Chiral and Uyupampa) follow a N140°trend separated by 6 km, aligned with the Ayo-Lluta-Arequipa Fault. The Nicholson volcano is located near the northwestern limit of the Arequipa basin. Eruptive activity has been characterized by eruptions that generated scoria flows, pyroclastic fallout and lava flow emissions.

Yura Viejo
This volcano was built above Yura Tuff and Mesozoic sedimentary rocks of the Yura Group (Figure 2). The crater (~520 m diameter) is situated in the current location of the town of Yura Viejo, which is marked by a semi-circular structure opened to the southwest. The preserved deposits cover an area~0 .74 km 2 surrounding the emission center. Three stratigraphic columns in the crater ring allowed us to describe its eruptive history, which is characterized by a succession of pyroclastic (PDC) scoria flows, scoria fallouts and reworked ash layers, ending with a lava flow deposit (Figures 2, 3). The pyroclastic components are mainly the same in all of the deposits, distributed in different sizes and proportions.
In the eastern flank (Figure 3, profiles 1 and 2), the lower successions, up to 12.50 m thick (thickness from the surveyed profiles), directly overlie the Yura Tuff along an erosional surface, and form the basal part of the ejecta ring on the eastern flank. These successions are characterized by scoria flows (10-150 cm thick), interbedded with bombs and lapilli fallout (20-80 cm thick), and ash layers (2-18 cm thick).

FIGURE 2
Geological map of the Yura Viejo volcano, embedded on a 60-cm resolution digital Elevation Model (DEM) map. Gray numbered dots indicate the location of stratigraphic sections shown in Figure 5. In the lower-right corner, a photograph of tephra fallout beds covered by lacustrine sediments.
Block and lapilli fallout beds are moderately to well sorted, clastsupported, non-agglutinated, and have variable thicknesses (20-80 cm). The juvenile scoria of these layers are dark gray; the size of scoria fragments varies from 2 to 50 cm, and rare dense lava fragments (1-10 cm) are also found. The lithics are composed of quartz sandstone and shale from Yura Group, with variable maximum size (10-40 cm) and scarce ignimbrite fragments depending on the layer's thickness. These rocks are also incorporated as xenoliths in juvenile blocks. The ash layers are pinkish-gray in color, and have variable thickness in the log, from 2 to 18 cm in the lower layers, and up to 18 cm in the upper layers.
The upper succession is 11.60 m thick (from surveyed profiles), and consists predominantly of scoria flows, interbedded with ash-fall (mainly remobilized) and few blockand-lapilli fallout layers. The scoria flows are massive, unsorted, matrix-supported and dark gray deposits composed of blocks (30-40%) embedded in a fine matrix (60-70%). Blocks are scoria (22-30%, 15-60 cm maximum size), and quartz-sandstone and shale lithics (3-10%, 20-75 cm maximum size). Two blockfallout deposits are 48 cm and 90 cm thick respectively. Scoria blocks have maximum sizes of 7 and 50 cm, while sedimentary lithics are 20 cm in diameter. Ash fall layers are pinkish gray in color, moderately welded, and their thickness vary from 5 to 62 cm. Bread-crust scoria bombs reaching diameters up to 50 cm are found at~1 km from the vent.

Frontiers in Earth Science
frontiersin.org 11 Described (successions) pyroclastic deposits in the eastern flank are not well-preserved in the southwestern flank (Figure 3, profile 3). However, scoria block fallouts interbedded with lahar flow deposits crop out over dark shale and sandstones strata from the Yura Group (Figures 2, 3) in a topographic low in this zone. Scoria-fall deposits are similar to those described above, whereas lahar flow deposits (1-2 m thick) are matrix-supported, compacted, pink and contain scoria fragments embedded in a muddy fine matrix. Superimposed on these deposits, a massive, dark-gray, a'a lava flow (~15 m thick) dated by 40 Ar/ 39 Ar at 167.8 ± 27.0 ka (Table 4) was emplaced filling a small valley oriented to the southwest, reaching a distance of~1 km from the crater (Figure 2). This lava flow shows a brown basal breccia of <70 cm thick. This volcano was probably emplaced in an ancient riverbed, forming a dam where lacustrine sediments were accumulated. These sediments overlap the scoria-fall deposits in the northeastern side of the volcano (Figure 2).
Due to its resistance to erosion, the lava emitted by the Yura Viejo volcano is well preserved, and the ravine that previously existed in the area opened up on its sides, forming an elevated plateau (relief inversion process). On the other hand, the tephra generated by the explosions are readily eroded, and are currently found only in the upper part of the surrounding hills. Deposits of Yura Viejo are cut by NW-SE active faults, in areas close to the vent (Figure 3).

Ccapua
This eruptive center (4.39 km 2 ) overlies the Arequipa Airport Ignimbrite. It is formed by two overlapping structures of low height, with two craters that lie below the surrounding ground level (Figure 4). The first emission center (crater I in Figure 4) has a cone basal diameter (Wco) of 640 m, the crater diameter (Wcr) reaches 360 m, and the average cone height (Hco) is 15 m, while the average flank slope is 6°. A basal deposit covering the ignimbrite is composed of scoria with mainly lapilli and scarce bombs (<1-m diameter), followed by a massive, unsorted, dark gray, and clast-supported (~70% blocks) deposit. The thickness of the ejecta ring is variable, with a major concentration in the southeastern flank.
The second emission center covers the northwestern half of the first crater, and is also partially covered on its northwestern edge by a blocky lava flow deposit (Uyu2) from the Chachani volcanic cluster (Uyupampa lava field; Aguilar et al., 2022). The Wco reaches 480 m, the Wcr is 320 m, the average Hco is 10 m, and the average flank slope is 10°. Lava flows and agglutinated scoria blocks form the rim of the crater (Figure 4). The base of the ejecta ring has moderately laminated layers of lapilli embedded in a whitish ash size matrix, overlain by a non-stratified deposit of dark red and black scoria blocks immersed in a medium to coarse brown ash matrix. The thickness of the uppermost scoria blocks deposit is irregular, and has its maximum value on the southwestern flank (~40 m). This last scoria deposit yielded a 40 Ar/ 39 Ar age of 97.7 ± 10.4 ka ( Table 4). The floors of both craters are partially filled with fine sediments (sand and silt) suggesting the presence of water in the past. One km west of the craters, a~20-m thick and dark-gray a'a lava flow crops out beneath tephra-fallout layers. The front and the upper surface of this lava flow is visible only along 120 m, making it difficult to determine the source vent. However, based on its proximity to the craters, the flow direction and the similarity in mineral content (mostly olivine phenocrysts) we assume that it was emitted from Ccapua volcano.
In addition, a massive deposit of dark brown and black scoria blocks and bombs embedded in a fine matrix was emplaced following a southwestern direction, and reaches distances as far as 2 km. This deposit has a maximum thickness of 60 m at 1.2 km from the vent. Ballistic bread-crust bombs attributed to Ccapua, mantle the Arequipa Airport Ignimbrite in the surrounding areas at distances as far as 2 km.

Nicholson
Nicholson volcano (1.42 km 2 ) was built at the southwestern edge of the Yura Monogenetic Field on a plain formed by thẽ 1.65 Ma Arequipa Airport Ignimbrite (Paquereau-Lebti et al., 2008; Figure 5), which is cut by the Yura River to the east. The volcanic edifice consists of two superimposed scoria cones. The older cone has a mean Wco of 540 m, Wcr of 200 m, Hco of 90 m, and average flank slope of 28°. The younger cone is located in the northern part, and its cone morphology is characterized by semicircular outlines, with a very well defined circular shape crater. This cone has a 520 m Wco, 210 m Wcr, 75-m Hco, and 26°a verage flank slope. Both cones were built by massive accumulations of black to dark red, poorly welded, unsorted, and clast-supported deposits, lacking any internal bedding. The constituents of these deposits are reddish and dark scoria clasts  (Table 4). In the lower part of the southeastern flank, 40-60 cm thick muddy and matrixsupported (50% of matrix) lahar deposits are observed, which incorporated angular scoria fragments (<10 cm in diameter). In most of its extension, Nicholson volcano is covered by a layer of reworked ash (~7-10 cm thick) from the 1600 B.P. eruption of Huaynaputina volcano, located 104 km to the southeast (e.g., Japura, 2018;Prival et al., 2020).

Uyupampa
Uyupampa volcano overlies the Yura Tuff, and its deposits cover 0.30 km 2 . This volcanic center seems to be strongly eroded and its deposits were used as quarries for housing construction. However, a semicircular structure of~300 m diameter is found in the present location of the town of Uyupampa, which may indicate the crater of the volcanic center. It is important to note that this monogenetic volcano is not related to the younger "Uyupampa lava-field" (blocky lavas) belonging to the Chachani Volcano Cluster (e.g., Aguilar et al., 2022) cropping out to the left side of La Paccha ravine (Figure 1).
The basal deposits of Uyupampa monogenetic volcano were emplaced over the Yura Tuff (Figure 6), and consist of pinkish ash layers (10-20 cm thick) interbedded with dark gray, moderately sorted scoria blocks/bombs -lapilli (35-130 cm thick) fallout deposits. At~120 m east of the crater, pinkish ash-fall layers with scarce scoria fragments (<10 cm diameter) grade to unsorted, clast-supported, brown-black scoria block/bomb fallout deposits (up to 70 cm size and 95% blocks). These layers are overlain by diffusely laminated and dark gray scoria and lithic ash-lapilli layers. The upper units consist of an accumulation of unsorted, clastsupported, brown and black scoria block/bombs layers. Most of the pyroclastic deposits contain bread-crust and cauliflower texture bombs that reach up to 100 cm in diameter (at 300 m from the crater). Scoria bombs and other tephra fall deposits have quartzsandstone and shale xenoliths and loose clasts (<30 cm).
The eruptive activity of this volcanic center culminated with the emplacement of a massive, olivine-rich a'a lava flow dated by 40 Ar/ 39 Ar method at 63.2 ± 9.0 ka (Table 4) which reached up to 500 m from the vent, following a southwest topographic low ( Figure 6).

El Chiral
This volcanic unit was emplaced near the Ayo-Lluta-Arequipa Fault, covering 1.54 km 2 . It is formed by a lobularshaped lava dome (Figure 7), which is cut by a~50-m thick and 1.60 km long lava flow deposit. This lava flow fills a small NE-SW depression where volcanic deposits of the Tacaza Group and quartz sandstones (Hualhuani Formation) of the Yura Group crop out. The geoforms of the vent and lava bodies of El Chiral suggest a cylindrical shape conduit, and emplacement of the lava in a gently inclined topography (e.g., Murcia and Németh, 2020).

Patacocha
This eruptive center (~0.45 km 2 ) shows an irregularly shaped and complex edifice formed by two coalescing cones and three Frontiers in Earth Science frontiersin.org 13 craters. The northwestern cone is capped by two coalescing craters (craters II and III), whereas the southeastern cone (crater I) has only one crater (Figure 8). Using Google Earth images and GIS software, we estimated the morphological parameters of the two cones: 1) Wco of 420 m, Wcr of 220 m, Hco of 24 m, and average flank slope of 14°for the northwestern cone. The smaller crater of 80-m diameter cuts the first larger crater. 2) Wco of 380 m, Wcr of 135 m, Hco of 40 m, and the average flank slope is 18°for the southeastern cone. The material expelled by this volcano covers deposits of the Tacaza Group, and apparently closed the channel of a stream, causing the damming of water and the formation of a lagoon called Patacocha (highland lagoon in Quechua). Due to the difficulty in accessing this area, it was not possible to take samples for laboratory analysis.

Volcano-structural setting
The local tectonic setting in which the Yura Monogenetic Field was emplaced consists of four groups of faults: N130°, N160°, N80°and N20°. This configuration is similar to the nearby Chachani volcano cluster described by Aguilar et al. (2022).

1) On satellite images, a N130°oriented lineament cuts the
Potrero-Airport domes and coincides with the collapse scar of the middle unit of this group of domes (Aguilar et al., 2022). The N120°-130°faults extend to the east of the city of Arequipa and offset pyroclastic deposits at the base of El Misti volcano (Upper Pleistocene) in the Chiguata zone (Thouret et al., 2001). Another remarkable fault is the dextral and normal fault that cuts across the western flank of El Misti volcano (Thouret et al., 2001) and extends towards the east flank of Estribo stratovolcano with a N130°orientation. 2) At a regional scale, N120°-130°and N160°faults and lineaments display an en-echelon fault system, which extends to El Misti stratocone and Pichu-Pichu volcanic complex. In the eastern central part of the Chachani volcanic cluster, this system cuts these edifices. The orientation of these lineaments is probably controlled by perpendicular N80°faults. 3) A group of long faults and lineaments (~20 km), with a prevailing N80°direction, allow us to define a NE-SW structure which delimits the northwestern edge of the Arequipa basin (Mering et al., 1996;Thouret et al., 2001;Benavente et al., 2017). This fault intersects the Ayo-Lluta-Arequipa Fault (~N130°strike-slip and normal fault; Figure 1B). 4) Finally, the last group of short lineaments follows a N20-40°d irection observed in the southern sector of the Chachani volcanic cluster. A remarkable feature in this area is the N20-40°predominant orientation of the majority of small eruptive vents, which have built the voluminous Airport-Potrero dome coulées. Uyupampa and Yura Viejo volcanoes in the Yura Monogenetic Field follow a similar trend.
The extensional zone along N80°and N135°-trending faults is associated with fissures that allow the ascent of magmas to the surface (e.g., Huaman-Rodrigo et al., 1993). Magnetic measurement data for Nicholson volcano identified two E-W trending tabular bodies beneath the edifice, suggesting dykeshaped feeding conduits (Rivera, 2021). E-W lineaments are also observed between Ccapua and Yura Viejo Volcanoes, and El Chiral and Patacocha volcanoes. Uyupampa and Yura Viejo volcanoes are aligned following a N40°trend. In the southwestern flank of Chachani, several small vents are also aligned along N10°and N40°eruptive fissures (Aguilar et al., 2022).
The N130°and N80°structures that we identified in the Yura Monogenetic Field and Chachani Volcano Cluster resemble those described by Huaman-Rodrigo et al. (1993) in the region of the Río Colca -Nevado Ampato-Sabancaya volcano, located about 50 km northwest of Yura (Maca and Huambo faults). In the area near Huambo NW of Yura, and in the northern flank of Chachani Volcano Cluster, N80°normal faults (Huaman-Rodrigo et al., 1993;Sempere and Jacay, 2006;Aguilar et al., 2022) are considered active because they cut recent Late Pleistocene lava flows and Holocene colluvial deposits. A destructive earthquake took place on one of them at Maca, on the left bank of the Río Colca canyon in 1992 (e.g., Antayhua et al., 2001). Other similar events occurred in August 2016 and March 2022, damaging houses and infrastructure.

Modal mineralogy
The mineral content of all analyzed samples from the volcanic centers of Yura Monogenetic Field is similar (Figure 9 and Table 2). All analyzed samples are porphyritic (7-10 vol%) and vesicular. In all samples, the dominant mineral is olivine, while the groundmass contains ortho-and clinopyroxene, plagioclase and glass. Olivine (7-10 vol%; ≤2 mm) occurs as subhedral phenocrysts with subrounded and embayed edges. In some cases, it has skeletal and breakdown textures (Figure 9). Pyroxene is present only in the groundmass (20-42 vol%; <400 μm; clinopyroxene > orthopyroxene) together with plagioclase microlites (26-53 vol %; <300 µm), giving the matrix a pilotaxitic texture by almost all samples (Table 2). Volcanic glass (<5 vol%) is found encompassing the microlites. The opaque minerals (<5 vol%; <800 µm) present anhedral and subhedral forms, and occur as accessory minerals as inclusions in olivine, or as replacement in crystal molds and disseminated in the groundmass ( Table 2). The vesicles (5-29 vol %) have irregular shapes.

Whole-rock major and trace element chemistry
Analyzed samples are given in Table 3. The scoria and lavas consist of high-K (1.60-2.01 wt% K 2 O) calc-alkaline rocks, which is a typical characteristic of the modern volcanic arcs in the Central Volcanic Zone (e.g., Wilson, 1986). Basaltic-andesite is an ubiquitous composition throughout the entire volcanic field, ranging from 50.9 to 55.6 wt% SiO 2 (Peccerillo and Taylor, 1976; Frontiers in Earth Science frontiersin.org 16 Figure 10). The least differentiated olivine-rich samples are found in the Ccapua and Nicholson volcanic centers. In Harker diagrams, some of which are shown in Figure 10, MgO, CaO, Fe 2 O 3 contents show negative correlation with SiO 2 indicating a common differentiation trend. Similarly, transition elements (e.g., Ni, Cr, V) and Dy are negatively correlated to SiO 2 suggesting that fractional crystallization is the major differentiation process in these magmas. Most incompatible elements such as Rb, U and Th together with Sm are positively correlated with increasing SiO 2 .

FIGURE 10
Chemical composition of Yura monogenetic volcanoes. Total Alkali-Silica versus SiO 2 diagram (Le Bas et al., 1986) showing the classification of Nicholson, Ccapua, Yura Viejo and Uyupampa volcanoes. Major oxides are shown as recalculated volatile-free values. Ca, MgO, Sr, and Rb versus SiO 2 diagrams show the compositional variation. Multi-element diagrams normalized to primitive mantle (Sun and McDonough, 1989) showing trace element patterns of rock samples.

Frontiers in Earth Science
frontiersin.org 6 Discussion 6.1 Eruptive style and classification of volcanoes in the yura monogenetic field At monogenetic volcanoes, the internal magma dynamics that occurs at shallow depths and the characteristics (form and dimension) of the conduit are linked to the way that magma is released at the surface, and the volcanic edifice construction (Murcia and Németh, 2020, and references therein). The topography where the magma is deposited also plays a role in the geoform of the volcano. The eruptive activity forming the monogenetic volcanoes usually consist of magmatic or phreatomagmatic explosions and in many cases accompanied by effusive emissions (Németh and Kereszturi, 2015). Field observations and description of the deposits in the Yura Monogenetic Field allow us to suggest: 1) Nicholson (and probably Patacocha) volcanoes were built by Strombolian (magmatic) eruptions. 2) Ccapua, Yura Viejo and Uyupampa volcanoes show ash to block particle sizes embedded in a muddy matrix, which suggests the interaction between magma and water, generating phreatomagmatic eruptions. However, accumulation of mainly scoria blocks/bombs and lava flows suggests Strombolian eruption phases. 3) El Chiral is composed of lava flows and coulées. Due to the lack of field data, it is not possible to characterize the type of eruption, although it was probably mainly effusive.
Basic measurements of morphometric parameters such as the cone basal diameter (Wco), crater diameter (Wcr) and cone height (Hco) allow obtaining morphometric signatures for classification of monogenetic volcanoes (e.g., Wood, 1979;Németh and Kereszturi, 2015). This classification includes spatter cones (Wco = 0.08 km, Wcr/Wco = 0.36 km and Hco/ Wco = 0.22 km), scoria/cinder cones (Wco = 0.8 km, Wcr/Wco = 0.4 km and Hco/Wco = 0.18 km), as well as maars and tuff rings (Wco = 1.38 km, Wcr/Wco = 0.6 km and Hco/Wco = 0.02 km) from average dimensions (Wood, 1979). Based on characterization parameters described by Németh and Kereszturi (2015), and Wood (1979), we suggest that the eruptive activity in the study area produced scoria cones (Nicholson and Patacocha), maars (Ccapua, Yura Viejo and Uyupampa), and a dome coulée and lava flow (El Chiral). A summary of the main characteristics and morphologic parameters is presented in Table 5. Figure 11 shows the comparison of Yura Monogenetic Field with surrounding Pleistocene volcanic rocks. Major element oxides show similar patterns among Yura Field, Chachani cluster, Yura Tuff and Arequipa Airport Ignimbrite. The Yura Field samples are the least differentiated magmas in the Arequipa basin; only some lava and scoria samples from the Chachani cluster exhibit similar composition. The Chachani volcanic cluster shows an intermediate and wider compositional range from basaltic andesite to rhyolite (54-71 SiO 2 wt%). The most differentiated samples belong to the Arequipa Airport Ignimbrite, which has a homogeneous silica content (74-76 SiO 2 wt%; Paquereau-Lebti et al., 2006). The Yura Tuff also shows highly differentiated compositions, but with a less homogeneous silica content (70-76 SiO 2 wt%).

Chemical comparison with surrounding volcanoes
Incompatible-element ratios such as Ba/Th vs. Dy/Yb and Dy/Y vs. Sm/Yb (Figure 11) exhibit overlapping fields with higher values in Ba/Th and Dy/Yb in the Yura Monogenetic Field and Chachani cluster. By contrast, the more evolved magmas of Arequipa Airport Ignimbrite and Yura Tuff show lower Dy/Yb and Sm/Yb values.
Due to the proximity of their sources, magmas generating the ignimbrites, the Chachani cluster and the Yura Monogenetic Field, could have a similar source and have passed through a similar substrate in their way to the surface, so they would have similar contaminants (e.g., sedimentary rocks of the Yura

Frontiers in Earth Science
frontiersin.org 19 Group). However, the abundance of xenoliths and subangular bedrock fragments found in the Yura Monogenetic Field deposits indicate that the assimilation process was much less effective in them. These monogenetic rocks and less differentiated lavas from Chachani cluster show similar basaltic-andesite compositional patterns ( Figure 11). In general, rocks from the Yura Monogenetic Field may represent the mafic end-member of the Pleistocene-Holocene volcanic products in the Arequipa basin and surrounding areas (e.g., Aguilar et al., 2022). At a regional scale, these magmas show one of the least differentiated magmas in the Central Volcanic Zone (e.g., Delacour et al., 2007;Mamani et al., 2010;Torres et al., 2020).

Chronology of monogenetic volcanism in the yura valley
Magmatic activity in the Arequipa region occurred continuously during the Pleistocene. The Arequipa Airport Ignimbrite marks a moderate to large magnitude explosive event (~18 km 3 ) that occurred c. 1.65 Ma (Paquereau-Lebti et al., 2006). This event emplaced pyroclastic density currents that filled the Arequipa depression, and caused the formation of a possible caldera (García et al., 1997). The geographic location of the ignimbrite source is likely covered by the edifices composing the Chachani volcano cluster (Paquereau-Lebti et al., 2006Aguilar et al., 2022). After caldera formation, andesitic lava flows were emplaced in the western edge of Chachani, with outcrops now located in the Quebrada La Paccha (Figure 1). Yura Tuff (c. 1.28 Ma; Paquereau-Lebti et al., 2006) overlies those lavas and fills the Yura valley (formed between the current northwestern side of Chachani and the wall formed by Yura Group units; Figure 1). The Chachani volcanic cluster started its activity at 1.01 Ma on the northeastern side, having recorded continuous activity until c. 641.8 ka with a possible hiatus that lasted until c. 463 ka. After c. 463 ka, the eruptive activity migrated to the central and southern part of the whole cluster (Aguilar et al., 2022). Towards the lower western flank of Chachani volcanic cluster and covering ignimbrite deposits, the Uyupampa lava field is composed of porphyritic lava flows, one of which yield an age of c. 232 ka. The youngest dated rock in the Chachani Volcanic Cluster belong to a block-and-ash deposit from Cabreria dome (c. 56.5 ka; Aguilar et al., 2022).
The oldest lithological units, on which the volcanoes of the Yura Monogenetic Field (e.g., El Chiral and Patacocha) were emplaced, belong to the Tacaza Group (c. 30-24 Ma; e.g., Mamani et al., 2010) in the northwestern sector of the field. However, the degree of preservation of volcanic structures such as the craters (cones), lobes and ridges (lava flows) suggests that Patacocha and El Chiral are much younger than c. 24 Ma. The ages obtained for the Yura Viejo, Ccapua, Nicholson and Uyupampa volcanoes match the stratigraphy of the area since they cover older deposits of the Arequipa Airport Ignimbrite and the Yura Tuff. The c. 168 ka age of Yura Viejo is in accordance with the degree of preservation of its deposits. Nicholson and Ccapua volcanoes were built over deposits of the Arequipa Airport Ignimbrite, while the Uyupampa and Yura Viejo volcanoes directly overlie the Yura Tuff. Based on the degree of preservation of the Patacocha scoria cones and El Chiral lava flow-dome coulée, we suggest that these volcanoes are relatively contemporaneous to the dated ones. Thus, the Yura Monogenetic Field would have been active in a time window of 141 kyr, between c. 195 ka and c. 54 ka. We note that this range of ages is younger than expected since the second crater of Ccapua volcano was partially covered by a lava flow from the Uyupampa field (Chachani volcanic cluster, Figure 1). Detailed fieldwork however, allowed us to determine that Ccapua was formed after the emplacement of the lava flow (Uyu1) dated by 40 Ar/ 39 Ar at 231.7 ± 36.2 (on groundmass; Aguilar et al., 2022), and was later covered by another lava flow (Uyu2) which should be younger than 97.7 ± 10.4 ka. The new ages obtained allow us to suggest that the eruptive activity of the Yura Monogenetic Field temporally coincided with the youngest edifices of Chachani volcanic cluster (Aguilar et al., 2022) and the basal units of the Misti composite cone (≤112 ka; Thouret et al., 2001). Future radiometric dating of the volcanic products from Patacocha and El Chiral are necessary to better constrain the age and duration of the magmatism in this monogenetic field.
6.4 Yura monogenetic field as potential geosite for education and awareness of volcanic hazards Recent volcanic activity in southern Peru has caused significant impacts to the population and economic activity for the communities settled in surrounding areas. Recent volcanic eruptions of the Ubinas volcano (2006-2009, 2014-2017 and 2019) have affected agriculture, critical infrastructure, and forced authorities to evacuate the people to shelters in safe areas (e.g., Rivera et al., 2010). The Sabancaya volcano began a new long-lasting vulcanian eruption process in 2016, and continues to this day, with no signs of ending in the near future. The ash emitted by frequent explosions of Sabancaya have caused damage to agriculture, livestock, and has even affected air traffic in the city of Arequipa (e.g., Aguilar, 2019;INGEMMET, 2021). In addition, during rainy seasons, the occurrence of lahars in the ravines that descend from the Misti and Chachani volcanoes has caused the destruction of homes, damage to basic infrastructure, and casualties (e.g., Rivera et al., 2018). These events remind us that volcanic eruptions and associated processes can negatively impact the development of communities; therefore, it is necessary that the population Frontiers in Earth Science frontiersin.org 20 and decision-makers develop measures to reduce the impact of these hazards. According to the record of eruptive activity, the hazards related to monogenetic volcanoes in the Yura Monogenetic Field are tephra-fall, ballistic fall, scoria flow, lava flow and lahar. Thus, the study of structures and deposits of the Yura Monogenetic Field can allow the identification of potential geosites for scientific, educational and geoturism uses and related to volcanic hazards. Nicholson is an easily accessible volcano to explain Strombolian eruptions and the formation of scoria cones. Ccapua volcano clearly shows the superposition of two craters and the formation of a maar which is partially covered by a younger lava flow. The Yura Viejo volcano shows explosive Strombolian-type activity combined with lahar flows, culminating with the emplacement of a lava flow. In addition, outcrops of ignimbrites that represent violent explosive eruptions can be seen in the area, as well as very representative sequences of sedimentary rocks of the Yura Group. The small monogenetic volcanoes also provide superb views of the surrounding larger volcanoes, and this, combined with their accessible outcrops and clear geology, provides an especially rich educational and geotouristic resource. This is especially important for communicating and understanding hazards and risk, as the urban centers and human infrastructure are visible in the landscape. In recently published press releases, it was indicated that the town of Yura Viejo is located in the crater of an extinct monogenetic volcano. This information was wellreceived by the local inhabitants and has attracted the attention of visitors, whose flow is increasing, not only to Yura Viejo, but also to the nearby Ccapua waterfalls (La Paccha ravine near Ccapua maar; Figure 1).
Due to its proximity to the urban area of Arequipa, and the good preservation of its deposits, the Yura Monogenetic Field has a great potential to promote geotourism for recreational, scientific and educational purposes in the face of the volcanic hazards (e.g., Galas et al., 2018;Arias et al., 2021). The valorization of these geosites could help in the popularization of these places and above all in their preservation by using an information simple to understand. There are examples of monogenetic fields that are growing tourist attractions, such as the Andagua monogenetic volcanoes (e.g. Galas et al., 2018), which are located within the Peruvian "Colca and Andagua Volcanoes" UNESCO Global Geopark. Other geosites are La Garrotxa monogenetic field in Spain (e.g., Planagumà-Guàrdia et al., 2022), and the Chaine de Puys in Clermont-Ferrand, France (e.g. Vereb et al., 2020).

Concluding remarks
The study of the eruptive history of monogenetic volcanic centers is useful to define the sequence of events and the dynamics of the eruptions that built these volcanic edifices.
Such volcanoes are key for understanding the geological evolution of an area, as they reflect the magma dynamics, the tectonics and the subsequent surface processes. Therefore, the combination of volcanic stratigraphy, geochemical and petrographic studies of volcanic activity in the Yura Monogenetic Field is important to better understand the evolution of volcanism in southern Peru. The study of the deposits has allowed us to reconstruct the history of this monogenetic field emplaced after medium-size ignimbrites formation in the Arequipa tectonic depression. Our new data suggest that volcanism in the Yura Monogenetic Field was active between c. 195 ka and c. 54 ka, which is more recent than previously expected (c. 1278 ka-c. 232 ka). The stratigraphic position and the degree of preservation of cone craters, lobes and ridges in lava flow/coulées indicate that El Chiral and Patacocha are much younger than Tacaza Group (Oligocene) as they cover volcanic units assigned to this group ( Figure 1). However, due to their geographical proximity, similarity in their petrographic and chemical characteristics these volcanoes are contemporaneous to the dated ones. Hence, the volcanism in the area lasted for~141 ka, beginning with the Yura Viejo volcano and probably finishing with Uyupampa volcano.
Bulk-rock major and trace elements highlight the fact that compositions have not varied with the locations of vents along the monogenetic field. These volcanic centers have relatively homogeneous basaltic-andesite compositions (50.9-55.6 wt% SiO 2 ), being among the least differentiated Quaternary rocks in the surrounding areas of the Arequipa basin. The dynamic of eruptions is marked by Strombolian and phreatomagmatic explosive activity, combined with scattered lava flows and small pyroclastic density current deposits (scoria flows). In most of the volcanoes, scoria flows including glassy breadcrust and cauliflower bombs suggest phreatomagmatic activity that crop out together with lahar deposits.
The Yura Viejo volcano shows one of the most complex histories in the Yura Monogenetic Field. The observed deposits indicate that its activity began with phreatomagmatic evolving to Strombolian eruptions, which expelled rock fragments (tephra) up to 1 km away from the crater. Subsequently, it generated a lava flow that was emplaced in a small gorge reaching a distance of 1 km towards the current Yura River, ending its eruptive activity.
This contribution has revealed the presence of a monogenetic field in the northwestern edge of the Arequipa basin by locating, studying and dating unknown volcanic centers. The results obtained provide important information to improve the knowledge of recent volcanic activity in the Arequipa region, which is relevant for volcanic hazard assessment and related risk management. This information could be also useful for the knowledge of the territory and identification of geosites, promote their conservation and propose new areas for education about geohazards and risk.
Frontiers in Earth Science frontiersin.org 21