Search

Located in the heart of the Lesser Caucasus mountains, where the Arabian and Eurasian tectonic plates collide, Armenia occupies an exceptional geological position shaped through millions of years of subduction and collision. It is a unique place on the Earth recording extensive intrusive and volcanic activity related to the long‐standing continental convergence. The volcanoes of Armenia provide a rare opportunity to study the sources and processes involved in this unusual type of magmatism. More than 500 Quaternary volcanoes have been mapped in Armenia, most of them formed from single eruptive episodes. Among several large composite volcanoes, the mighty Aragats stands out as the largest volcano in Armenia and the region altogether. Volcanic deposits testify to the range of eruptive styles—from the ignimbrites formed in eruptions as explosive and voluminous as any seen globally in the modern era to the enormous fissure‐fed lava flows that form the Southern Caucasus flood basalt province, the smallest and youngest Large Igneous Province in the world. Several pre‐historical and historical eruptions have been documented, highlighting the potential for future volcanic activity in the region. In recent years, research has focused on the volcanic hazards associated with the Armenian Nuclear Power Plant, located in the foothills of Aragats volcano. This article highlights some of the extraordinary volcanic and intrusive features observed in Armenia and summarizes aspects of recent volcanological and petrological research.

The Lesser Caucasus mountains sit on a transition within the Arabia-Eurasia collision zone between very thin lithosphere (< 100 km) to the west, under Eastern Anatolia, and a very thick lithospheric root (up to 200 km) in the east, under western Iran. A transect of volcanic highlands running from north-west to south-east in the Lesser Caucasus allows us to look at the effects of lithosphere thickness variations on the geochemistry of volcanic rocks in this continental collision zone. Volcanic rocks from across the region show a wide compositional range from basanites to rhyolites, and have arc-like geochemical characteristics, typified by ubiquitous negative Nb-Ta anomalies. Magmatic rocks from the south-east, where the lithosphere is thought to be thicker, are more enriched in incompatible trace elements, especially the light rare earth elements, Sr and P. They also have more radiogenic 87Sr/86Sr, and less radiogenic 143Nd/144Nd. Across the region, there is no correlation between SiO2 content and Sr-Nd isotope ratios, revealing a lack of crustal contamination. Instead, “spiky” mid-ocean-ridge basalt normalised trace element patterns are the result of derivation from a subduction-modified mantle source, which likely inherited its subduction component from subduction of the Tethys Ocean prior to the onset of continent-continent collision in the late Miocene. In addition to the more isotopically enriched mantle source, modelling of non-modal batch melting suggests lower degrees of melting and the involvement of garnet as a residual phase in the south-east. Melt thermobarometry calculations based on bulk-rock major elements confirm that melting in the south-east must occur at greater depths in the mantle. Temperatures of melting below 1200 °C, along with the subduction-modified source, suggest that melting occurred within the lithosphere. It is proposed that in the Northern Lesser Caucasus this melting occurs close to the base of the very thin lithosphere (at a depth of ∼45 km) as a result of small-scale delamination. A striking similarity between the conditions of melting in north-west Iran and the southern Lesser Caucasus (two regions between which the difference in lithosphere thickness is ∼ 100 km) suggests a common mechanism of melt generation in the mid-lithosphere (∼ 75 km). The southern Lesser Caucasus magmas result from mixing between partial melts of deep lithosphere (∼ 120 km in the south) and mid-lithosphere sources to give a composition intermediate between magmas from the northern Lesser Caucasus and north-west Iran. The mid-lithosphere magma source has a distinct composition compared to the base of the lithosphere, which is argued to be the result of the increased retention of metasomatic components in phases such as apatite and amphibole, which are stabilized by lower temperatures prior to magma generation.

Argon isotopic analyses were funded by NERC grant to I. Savov, IP-1690-1116. PJS was funded through a NERC studentship as part of the Leeds York Spheres Doctoral Training Partnership (DTP) (grant number NE/L002574/1). Part of the field work and research were supported by base funding of the Institute of Geological Sciences (IGS) and a thematic project by the Science Committee of the Armenian Ministry of Education and Science (project #18 T-1E368). The post-collisional Syunik and Vardenis volcanic highlands, located in the southern Lesser Caucasus mountains (part of the Arabia-Eurasia collision zone) are host to over 200 monogenetic volcanoes, as well as 2 large Quaternary polygenetic volcanoes in the Syunik highland. The latter are overlain by lavas from the monogenetic volcanoes, suggesting there was a transition in the style of volcanic activity from large-volume central vent eruptions to dispersed small-volume eruptions. 12 new high quality 40Ar/39Ar ages are presented here, with 11 ages calculated by step-heating experiments on groundmass separates, and the final age obtained from total fusions of a population of sanidines. All the ages were younger than 1.5 Ma, except for one ignimbrite deposit whose sanidines gave an age of 6 Ma. While the bulk of the exposed products of post-collisional volcanism relate to Pleistocene activity, it is clear there has been active volcanism in the region since at least the late Miocene. All ages for monogenetic volcanoes in the Syunik highland are younger than 1 Ma, but to the north in Vardenis there is geochronological evidence of monogenetic volcanism at 1.4 and 1.3 Ma. An age of 1.3 Ma is determined for a lava flow from one of the polygenetic volcanoes- Tskhouk, and when combined with other ages helps constrain the timing of the polygenetic to monogenetic transition to around 1 Ma. The new ages illustrate a degree of spatio-temporal coupling in the formation of new vents, which could be related to pull-apart basins focussing ascending magmas. This coupling means that future eruptions are particularly likely to occur close to the sites of the most recent Holocene activity. The polygenetic to monogenetic transition is argued to be the result of a decreasing magma supply based on: (i) volume estimates for Holocene eruptions and for all monogenetic volcanoes and their lava flows in Syunik; and (ii) the volcanic stratigraphy of the Lesser Caucasus region which shows late Pliocene- early Pleistocene continental flood basalts being succeeded by a few large andesite-dacite volcanoes and then the most recent deposits consisting of small-volume scoria cones. The Syunik highland has the highest density of monogenetic centres in the Lesser Caucasus, which is taken to indicate this region has the highest magma flux, and was therefore the last location to transition to monogenetic volcanism, which is why the transition is most clearly seen there. There is no evidence from Sr-Nd-B isotope measurements for the exhaustion of fusible slab components in the mantle source, showing that an inherited slab signature can survive for millions of years after the end of subduction. Although volcanism in the Lesser Caucasus is currently waning, a future pulse of activity is possible. Publisher PDF

This work was supported as part of Patrick Sugden's PhD studentship funded through the Leeds-York SPHERES NERC doctoral training partnership (NE/L002754/1). Part of the fieldwork and research was funded by the University of Leeds, the Carnegie Institution of Washington, the ERASMUS exchange programme (for IS) and the Science Committee of the Armenian Ministry of Education and Science (project #18T-1E368). The majority of the B isotope analytical work was supported by IGG-CNR funds P1600514. We present the first boron abundance and δ11B data for young (1.5-0 Ma) volcanic rocks formed in an active continent-continent collision zone. The δ11B of post-collisional volcanic rocks (−5 to +2‰) from the Armenian sector of the Arabia-Eurasia collision zone are heavier than mid-ocean ridge basalts (MORB), confirming trace element and isotope evidence for their derivation from a subduction-modified mantle source. Based on the low B/Nb (0.03-0.25 vs 0.2-90 in arc magmas), as well as low Ba/Th and Pb/Ce, this source records a subduction signature which is presently fluid-mobile element depleted relative to most arc settings. The heavier than MORB δ11B of post-collision volcanic rocks argues against derivation of their subduction signature from a stalled slab, which would be expected to produce a component with a lighter than MORB B, due to previous fluid depletion. Instead, the similarity of δ11B in Plio-Pleistocene post-collision to 41 Ma alkaline igneous rocks also from Armenia (and also presented in this study), suggests that the subduction signature is inherited from Mesozoic-Paleogene subduction of Neotethys oceanic slabs. The slab component is then stored in the mantle lithosphere in amphibole, which is consistent with the low [B] in both Armenian volcanic rocks and metasomatic amphibole in mantle xenoliths. Based on trace element and radiogenic isotope systematics, this slab component is thought to be dominated by sediment melts (or supercritical fluids). Previously published δ11B of metasediments suggests a sediment-derived metasomatic agent could produce the B isotope composition observed in Armenian volcanic rocks. The lack of evidence for aqueous fluids preserved over the 40 Myr since initial collision supports observations that this latter component is transitory, while the lifetime of sediment melts/supercritical fluids can be extended to >40 Myr. Publisher PDF