# Geologic setting of the Peña de Bernal Natural Monument, Querétaro, México: An endogenous volcanic dome

Gerardo J. Aguirre-Diaz, Alfredo Aguillón-Robles, Margarito

Tristán-González, Guillermo Labarthe-Hernández, Margarita Lopez-Martinez,

Hervé Bellon, Jorge Nieto-Obregon

► To cite this version:

Gerardo J. Aguirre-Diaz, Alfredo Aguillón-Robles, Margarito Tristán-González, Guillermo Labarthe-Hernández, Margarita Lopez-Martinez, et al. Geologic setting of the Peña de Bernal Natural Monument, Querétaro, México: An endogenous volcanic dome. Geosphere, 2013, 9 (3), pp.557-571. 10.1130/GES00843.1. insu-00841356

HAL Id: insu-00841356

https://insu.hal.science/insu-00841356v1

Submitted on 3 Sep 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Geologic setting of the Peña de Bernal Natural Monument,
Querétaro, México: An endogenous volcanic dome

Gerardo J. Aguirre-Díaz¹*, Alfredo Aguillón-Robles²*, Margarito Tristán-González², Guillermo Labarthe-Hernández², Margarita López-Martínez³*, Hervé Bellon⁴*, and Jorge Nieto-Obregón⁵*

¹Centro de Geociencias, Universidad Nacional Autónoma de México, campus Juriquilla, Boulevard Juriquilla 3001, Juriquilla, Querétaro, C.P. 76230, México

²Instituto de Geología, Universidad Autónoma de San Luis Potosí, San Luis Potosí, México

³División de Ciencias de la Tierra, Centro de Investigación Científica y de Educación Superior de Ensenada, Carretera Ensenada-Tijuana No 3918, Zona Playitas, C.P. 22860, Ensenada, México

⁴European University of Brittany, 6538, Domaines Océaniques, Institut Universitaire Européen de la Mer (IUEM), Université de Bretagne Occidentale, 29280 Plouzané, France

⁵Facultad de Ingeniería, Universidad Nacional Autónoma de México, Ciudad Universitaria, Delegación Coyoacán, 04510, México, D.F., México

ABSTRACT

Peña de Bernal is a natural monument located near the town of Bernal, in Querétaro State, central Mexico. It is one of the tallest monoliths of the world, with a maximum height of 433 m. Peña de Bernal was recently declared Intangible Cultural Heritage of Humanity Patrimony by United Nations Educational, Scientific, and Cultural Organization (UNESCO). In spite of being both a natural and cultural monument, little is known about its origin, physical characteristics, and chemical composition. It is a leucocratic-igneous rock intruding marine Mesozoic sedimentary rocks and has been misinterpreted as a pluton of Eocene or older age. However, this study shows that Peña de Bernal is a dacitic dome with $SiO_2$ = 67 wt% and an age of 8.7 ± 0.2 Ma. The complete Peña de Bernal body includes three plugs that crop out in an ~3.5 × 1.5 km area elongated N40°E. Texture of the rock is porphyritic, nearly holocrystalline (80 vol% crystals and 20 vol% glass), with a mineral assemblage of pyroxene, hornblende, biotite, plagioclase, and quartz, plus accessory apatite and zircon. Peña de Bernal dacite is a spine-type endogenous dome that was forcefully intruded through the Mesozoic sequence practically as a solid plug.

INTRODUCTION

Peña de Bernal, at 20°45'N, 99°57'W (Fig. 1), is a prominent peak that was declared a natural reserve area by the State Government of Querétaro in 2007 (Fig. 2). Since 30 September 2009, it has been included in UNESCO's Intangible Cultural Heritage sites (UNESCO, 2009) in order to preserve the traditions of Otomí and Chichimeca indigenous peoples as well as the Peña de Bernal itself. Furthermore, Peña de Bernal is known as the main entrance to the Sierra Gorda mountains area, a range that was declared in 2001 as a Biosphere Reserve by UNESCO (2001). In a few publications and in some web sites, including the official municipal (county) web page of the Bernal area and that of the National Institute of History and Anthropology (INAH) of Mexico, Peña de Bernal is identified as one of the highest monoliths in the world, with monolith defined as a prominent and large monolithologic rock (Fig. 2). The exact height of Peña de Bernal has not been well constrained, with reported heights ranging between 288 and 360 m (e.g., Coenraads and Koivula, 2007; INAH, 2009; Municipio Ezequiel Montes, 2009). From the topographic data of the official Mexican agency for National Geographic Studies and Information (Instituto Nacional de Geografía Estadística e Informática [INEGI], 2002) and from altimeter measurements carried out by us at the base and at the top of Peña de Bernal, we have calculated a height of 433 m on the northern side (Fig. 3). Height is 405 m on the southern side of Peña de Bernal with respect to the plain where Bernal town lies (2049 meters

above sea level [masl]). Either elevation chosen from these new values shows that the Peña de Bernal is apparently the highest monolith of the world.

Previous regional geologic mapping has considered Peña de Bernal as a Paleocene monzogranodioritic pluton (Segerstrom, 1961) or a subvolcanic Miocene felsic body (Carrillo-Martínez, 1998). However, the mechanism of emplacement has remained unknown. Before this work, there were no scientific studies specifically on Peña de Bernal's geologic nature and origin. The purpose of this study is therefore to provide detailed information on the geologic setting, including the stratigraphy, structure, geochronology, and geochemistry, to answer fundamental questions of when and how Peña de Bernal was formed and why it is topographically prominent with respect to its surroundings.

GEOLOGY

Peña de Bernal is located at the intersection of three major geologic provinces: the Miocene-Quaternary Mexican Volcanic Belt (Demant, 1978; Aguirre-Díaz et al., 1998; Siebe et al., 2006), the mid-Tertiary Sierra Madre Occidental ignimbrite-dominated sequence (McDowell and Clabaugh, 1979; Aguirre-Díaz et al., 2008), and the Mesozoic Sierra Madre Oriental fold-thrust belt (Suter, 1987; Bartolini et al., 1999), which in this area corresponds to the Sierra Gorda sector (Fig. 1). A simplified geologic map of the Peña de Bernal area is shown in Figure 4. The stratigraphy is briefly described below, from the oldest to the youngest unit (Fig. 5).

<footnote>*Emails: Aguirre-Díaz: ger@geociencias.unam.mx; Aguillón-Robles: aaguillonr@uaslp.mx; López-Martínez: marlopez@cices.mx; Bellon: herve.bellon@univ-brest.fr; Nieto-Obregón: nieto@servidor.unam.mx.</footnote>
<figure><img src="image_2.png" /><figcaption>Figure 1. Index map showing the position of Peña de Bernal with respect to some of the geological provinces of Mexico and the geography of the local area; the study area is highlighted in the inset.</figcaption></figure>

## Las Trancas Formation (Jt)

Las Trancas Formation was defined by Segerstrom (1961) as a Jurassic-Early Cretaceous (Kimmeridgian-Barremian) marine sedimentary sequence. It consists of dark-gray

fissile shales, calcareous limonites, and dirty limestones with pyrite and minor amounts of graywackes and flint (Segerstrom, 1961, 1962; Chauve et al., 1985; Carrillo-Martínez, 1990). In the Bernal area, this unit occurs as a thinly bedded sequence of fissile shales and limestones

with minor amounts of graywackes. Bedding is generally undulated, with thicknesses in the order of tens of meters. The shales and limestones of Las Trancas Formation were intensively deformed during the Cretaceous-Tertiary (K-T) Laramide orogeny (Chauve et al., 1985;
<figure><img src="image_3.png" /></figure>

<figure><img src="image_4.png" /></figure>

**Figure 2.** (A) Panoramic view toward the north of Peña de Bernal monolith and corresponding regional index map (inset). (B) Panoramic view toward the east of Peña de Bernal from summit of San Martín volcano. Note the elongated shape of the main plug of the Peña de Bernal dacitic dome and the rounded plug of Tortuga. The Northern plug is not shown in the image.

Padilla-y-Sánchez, 1985), producing incompetent folding, eastward thrust transportation, and a low-grade regional metamorphism (Suter, 1987; Carrillo-Martínez and Suter, 1990; Carrillo-Martínez, 1998; Bartolini et al., 1999) that gave a fissile and slate character to the shales and limestones, but without reaching the point of the schist facies.

### El Doctor Formation (Kd)

El Doctor Formation is a middle Cretaceous (Albian–Cenomanian) unit that was defined by Wilson et al. (1955) from a locality near the town of El Doctor, in the Sierra Gorda portion of the Sierra Madre Oriental fold-thrust belt, and ~35 km to the NE of Bernal (Fig. 1). This formation is characterized by its origin as a reef and, thus, is composed mainly of platform calcareous small fossils and microfossils. In the mapped area, El Doctor Formation crops out to the NW of San Antonio de la Cal (Fig. 4), as an elongated outcrop between Las Trancas and Soyatal formations. The outcrop

is ~100 m thick and consists mostly of gray fossiliferous limestones in thick to medium thick banks.

### Soyatal Formation (Ks)

Soyatal Formation is a Late Cretaceous (Turonian) marine sedimentary sequence composed mainly of limestones and flint lenses (White, 1948). Index fossils near Bernal and found by us include ammonites of *Nowakites* sp.

at the locality of Dedhó, and *Texamites* sp. in the Tolimán canyon, confirming the Late Cretaceous epoch of this unit. Soyatal Formation is thinly to medium bedded, showing boudinage primary bedding. It is dark gray when fresh, and light gray to white when weathered. As in Las Trancas, this unit was folded and northeastwardly tectonically transported during the K-T Laramide orogeny (Chauve et al., 1985). Las Trancas Formation was thrusted on Soyatal Formation during this orogeny (Suter, 1987; Carrillo-Martínez, 1998), forming klippes at top of younger Soyatal. In the Bernal area, the top of Soyatal Formation is eroded, and thus, a complete thickness estimate was not possible, but a minimum of 300 m is considered our best approach.

### Lower Tertiary Continental Clastic Sediments (Ts)

Continental detrital deposits crop out to the NW of the mapped area and nearby the town of San Martín (Fig. 1). Within the mapped area these deposits occur as small patches in the SW lower flank of Peña de Bernal too small to be shown in Figure 4, and to the southwest of Cerro Azul, in the northeastern portion of the area. Thickness ranges from 0 to 20 m. Near Bernal this unit consists mainly of conglomerates with clast sizes of 2–30 cm in diameter, mainly from shales of Las Trancas Formation and black flint and limestones of Soyatal Formation. Clasts are rounded to semi-rounded and are supported by a red clay matrix, in parts weathered to a calcareous soil. The deposits at Cerro Azul are finer and consist of cross-bedded sandstones, green and red, with no base exposed, and covered by younger pyroclastic rocks. These clastic deposits rest unconformably either on Soyatal or Las Trancas formations. Tentatively, we position this unit in the lower Tertiary, after the K-T orogeny, in correlation with similar continental deposits of central Mexico, such as the Eocene red conglomerates of Guanajuato and Zacatecas (Edwards, 1955; Tristán-González et al., 2009), and the Cenicera formation of San Luis Potosí (Labarthe-Hernández et al., 1982).

<figure><img src="image_5.png" /></figure>

**Figure 3.** Topographic cross section outlining the height of Peña de Bernal. Map position of A–B section line is shown in Figure 4. masl—meters above sea level.
<figure><img src="image_6.png" /><figcaption>Figure 4. Geologic map of Peña de Bernal dacitic dome and surrounding area. Also indicated is the A–B section line shown in Figure 3.</figcaption></figure>

## Cerro Azul Epiclastic Deposits (Tcae)

Overlying the continental red bed deposits mentioned above, there is a sequence of epiclastic deposits consisting of green to yellow tuffaceous sands, in beds 20–100 cm thick, with cross bedding in places. Best exposures occur west of Cerro Azul, where the unit reaches up to 50 m in thickness (Fig. 4). The lower part of the sequence consists of beds of rounded to subangular clasts of quartz, plagioclase, and volcanic glass, supported by a clayey, chloritized matrix. The upper part is mainly a clay deposit, in beds of 2–20 cm. It is possible that these deposits, and particularly the upper clay beds, represent reworked basal ash deposits precursor to the emplacement of Cerro Azul ignimbrite, which conformably overlies this unit.

## Cerro Azul Ignimbrite (Tcai)

Cerro Azul ignimbrite is a widespread, densely welded, orange-brown unit, with up to 30 vol% of phenocrysts (2–3 mm) of quartz, sanidine, and biotite, supported in a matrix of devitrified glass shards. The thickness at Cerro Azul is up to 45 m, but more commonly it is around 30 m thick. The base of the ignimbrite includes a pink unwelded to poorly welded ash layer, up to 1 m thick, that lacks bedding, phenocrysts, and/or lithics. Conformably and continuously overlying this basal layer is the massive part of the welded ignimbrite. Above the ash layer, the ignimbrite contains sparse lithics of limestone and chert. K-Ar age dating performed on this unit yielded an age of 31.5 ± 0.7 Ma (Table 1).

Therefore, this ignimbrite corresponds to mid-Tertiary volcanism related to the Sierra Madre Occidental.

## Cerro Azul Andesite (Tcaa)

Andesitic lavas overlie Cerro Azul ignimbrite (Fig. 5). These lavas have been informally named as Cerro Azul andesite because the best exposures are found at Cerro Azul and at hills to the west of it. The andesite lavas usually form plateaus and apparently were fed from fissures. Some andesitic dikes with similar aspect and texture were observed nearby cutting the Cerro Azul epiclastic deposits and the Cerro Azul ignimbrite. Whole-rock K-Ar dating of this unit resulted in 29.4 ± 0.8 Ma (Table 1). It is generally a propylitized gray to green rock or
dark gray, when rare fresh andesite is found. In some sites, it has vesicles or amygdules of agate or chalcedony. Andesite is fine grained, with sparse phenocrysts with sizes of less than 1 mm of plagioclase, Fe-Ti oxides, and clinopyroxene (augite), which sometimes is bordered by a rim of amphibole.

## El Matón Ignimbrite (Tmi)

El Matón ignimbrite is one of the largest in distribution and volume of the region, and is commonly observed between Tolimán and San Miguel (outside the mapped area; Fig. 1). It is a partly welded, pink, rhyolitic ignimbrite, with eutaxitic texture. This is a crystal-rich ignimbrite, with up to 40 vol% of phenocrysts of quartz, sanidine, and oxidized biotite. Most crystals are subhedral and broken, with sizes between 1 and 3 mm. A matrix of devitrified glass shards is preserved. The degree of welding increases upward, from poorly to densely welded. It overlies Cerro Azul andesite, making the top of the hill with this name (Fig. 5). It also overlies Soyatal Formation and the lower Tertiary clastic sediments. K-Ar dating yields an age of 29.4 ± 0.7 Ma (Table 1). The ignimbrite is stratigraphically younger than 29.4 ± 0.8 Ma Cerro Azul andesite, although both resulted with the same isotopic age.

## San Martín Andesite (Tsma)

San Martín andesite includes a series of lavas erupted from San Martín volcano that are between the towns of San Martín and Bernal (Fig. 1). A couple of small outcrops can be observed at the southwestern part of the mapped area (Fig. 4). San Martín andesite overlies Las Trancas and Soyatal formations, and El Matón ignimbrite outside the mapped area. The lavas are dark gray to light gray and generally form thick lobes rather than plateaus. Lavas are porphyritic, with up to 20 vol% of phenocrysts of plagioclase and hornblende and a ground-

mass made of glass and microlites of the same phases plus sparse olivine. The structure of the lavas includes flow banding and flow foliation. K-Ar dating of two different lava units of San Martín volcano yielded ages of 11.0 ± 0.3 Ma and 10.2 ± 0.2 Ma (Table 1), suggesting a long-lived polygenetic history of this volcano. San Martín volcano is similar in age and composition with other middle Miocene stratovolcanoes in the region, such as La Joya and El Zamorano, which are related to the Mexican Volcanic Belt (Aguirre-Díaz, 2008).

## THE PEÑA DE BERNAL DACITIC DOME

### Description

Peña de Bernal is a highly crystalline and very resistant rock, gray to light gray when fresh and brown when weathered. It is composed of up to 80 vol% crystals and 20 vol% matrix glass; phenocrysts make up to 30 vol% of the total. Mineralogy includes plagioclase + orthopyroxene + hornblende + biotite + sanidine + quartz + Fe-Ti oxides + apatite + zircon (Fig. 6). It has a porphyritic texture with phenocryst content varying within 15–30 vol%. Amongst phenocrysts, plagioclase and hornblende are the most common and largest up to <2.5 mm; biotite and pyroxene are smaller and relatively less abundant. Mafic phases are generally oxidized either completely or along their rims (Fig. 6). The matrix is made of a microcrystalline mineral assemblage with all the phases mentioned above, plus interstitial glass. Plagioclase phenocrysts are generally zoned with complex oscillatory zoning and fritted zones (Fig. 6). Sanidine is rare, but it has been observed as phenocrysts. Glass is pale brown and occurs between the crystals. Textures suggest a long-term, complex cooling history of the Peña de Bernal magma, with resorption of phases and disequilibrium and reequilibration conditions, particularly recorded in the plagioclase. The large mineral assemblage

<figure><img src="image_7.png" /><figcaption>Figure 5. Composite stratigraphic column for the area shown in Figure 4. u—unconformity.</figcaption></figure>

includes most phases of the Bowen's series (except olivine), suggesting enough time for the nucleation and growth of all these phases with the liquid line of descent (Wilcox, 1979).

### Age

Two samples of Peña de Bernal were analyzed to determine its age. Samples were processed in separate laboratories and using different techniques in order to cross-check the results and to

<table><caption>TABLE 1. K-Ar and $^{40}$Ar-$^{39}$Ar ISOTOPIC AGES OF THE BERNAL AREA, QUERÉTARO</caption><thead><tr><th rowspan="2">Unit*</th><th rowspan="2">Sample number</th><th colspan="2">Coordinates</th><th rowspan="2">Age (±1σ)</th><th rowspan="2">$^{40}$Ar* (e-7 cc)</th><th rowspan="2">$^{40}$Ar* (%)</th><th rowspan="2">K<sub>2</sub>O (wt%)</th><th rowspan="2">Fraction†</th><th rowspan="2">Analyzed weight (g)</th></tr><tr><th>Latitude N</th><th>Longitude W</th></tr></thead><tbody><tr><td>Tpb</td><td>SM10-05</td><td>99°56.67'</td><td>20°45.58'</td><td>8.7 ± 0.2</td><td>7.24</td><td>71.4</td><td>2.13</td><td>GM</td><td>0.4759</td></tr><tr><td>Tpb</td><td>Qro-8Ar/Ar</td><td>99°56.73'</td><td>20°44.71'</td><td>8.7 ± 0.2<sup>Tpb</sup></td><td>4.52</td><td>60.5</td><td>—</td><td>PLG</td><td>0.6122</td></tr><tr><td>Tsma</td><td>SM09-05</td><td>99°59.04'</td><td>20°44.15'</td><td>10.2 ± 0.2</td><td>4.69</td><td>73.1</td><td>1.43</td><td>WR</td><td>0.5101</td></tr><tr><td>Tsma</td><td>SM03-05</td><td>99°58.39'</td><td>20°44.70'</td><td>11.0 ± 0.3</td><td>6.31</td><td>39.6</td><td>1.78</td><td>WR</td><td>0.6035</td></tr><tr><td>Tmi</td><td>SM14-05</td><td>99°57.16'</td><td>20°46.95'</td><td>29.4 ± 0.7</td><td>4.83</td><td>77.8</td><td>5.06</td><td>WR</td><td>0.5002</td></tr><tr><td>Tcaa</td><td>SM13-05</td><td>99°57.02'</td><td>20°47.02'</td><td>29.4 ± 0.8</td><td>11.5</td><td>49.8</td><td>1.21</td><td>WR</td><td>0.6223</td></tr><tr><td>Tica</td><td>SM12-05</td><td>99°56.76'</td><td>20°47.07'</td><td>31.5 ± 0.7</td><td>5.71</td><td>96.5</td><td>5.58</td><td>WR</td><td>0.5201</td></tr></tbody></table>

<footnote>*Units dated: Tpb—Peña Bernal dacite; Tsma—San Martín andesite; Tmi—El Matón ignimbrite; Tcaa—Cerro Azul andesite; Tica—Cerro Azul ignimbrite. All except Qro-8 are K-Ar ages done at Université de Bretagne Occidental, Brest, France, by Hervé Bellon (technique described in Bellon et al., 1981). $^{40}$Ar-$^{39}$Ar analysis done at División de Ciencias de la Tierra, Centro de Investigación Científica y de Educación Superior de Ensenada (CICESE), Baja California, México, by Margarita López-Martínez (technique described in Cerca-Martínez et al., 2000); Tp plateau age calculated with more than 75% of $^{39}$Ar released on four fractions. $^{40}$Ar*—extracted radiogenic $^{40}$Ar in standard cubic cm (cc) and percent (%).</footnote>

<footnote>†Material used: WR—whole rock, GM—groundmass, PLG—plagioclase.</footnote>
<figure><img src="image_8.png" /></figure>

**Figure 6.** Representative photomicrographs of Peña de Bernal dacite, all with crossed polarizers (except F). (A) Phenocrysts of plagioclase (PLA), hornblende (HBL), and biotite (BIO) in a microcrystalline + glass matrix. (B) Large crystals of acicular hornblende (HBL) partly oxidized. (C) Phenocryst of rounded quartz (center). (D) Large phenocryst of plagioclase showing complex oscillatory zoning and fritted zones. (E) Euhedral crystals of orthopyroxene with oxidized rims. (F) Accessory apatite (APA) and glass patches between matrix crystals. Graphic scales in mm and microns indicated in each photograph.

have a good accuracy of the timing of formation. The K-Ar age obtained at the University of Brest, France, yielded $8.7 \pm 0.2$ Ma (Table 1). The $^{40}$Ar-$^{39}$Ar plateau age obtained at the Centro de Investigación Científica y de Educación Superior de Ensenada (CICESE) research center in Ensenada, Mexico, yielded also $8.7 \pm 0.2$ Ma on plagioclase separate (Fig. 7 and Table 2). Therefore, $8.7 \pm 0.2$ Ma should be the age of Peña de Bernal dacite. The middle Miocene age of Peña de Bernal is consistent with the general age of other silicic domes in the central-northern sector of the Mexican Volcanic Belt (Aguirre-Díaz et al., 1998; Aguirre-Díaz, 2008).

## Structure

Besides the main igneous mass that forms the peak of Peña de Bernal, there are two other plugs with lower elevations, and all three plugs crop out in a $4.8 \text{ km}^2$ area (Figs. 2 and 4). The main plug extends for ~3.5 km from the highest peak that forms Peña de Bernal to the southwesternmost part at San Antonio de la Cal (Fig. 4). The other two plugs are La Tortuga, $0.85 \times 0.45$ km, at the western portion of the main outcrop, and the Northern plug, at the northeasternmost part of the outcrop (Fig. 4). The rock is similar in all plugs with some textural or dome foliation

variations; textures vary mainly in crystal size with crystals bigger in the main plug and smaller in the other two. Apparently, plugs were emplaced simultaneously, as there are not sharp contacts between them but rather transitional contacts. Peña de Bernal peak shows dome foliation that varies from ~40° inward dips at the margins to nearly vertical dips at the center of the dome (Fig. 8). This foliation follows the general elongated trend of the main plug to the NE. Peña de Bernal dacite is in intrusive contact with the country rock, which is primarily the Las Trancas Formation, although it is also in contact with the Soyatal limestone in some places, in particular along the southwestern margin of the outcrop area (Fig. 4). It is not a thermal metamorphic aureole caused by the dome intrusion, and contact metamorphism of the country rock is also not observed. The metamorphism of Las Trancas Formation was caused by a regional tectonic event (Suter, 1987) and not by the intrusion of Peña de Bernal dome. The contact between the Peña de Bernal dacite and the country rock is sharp (Fig. 9A). Most notable is the fact that it looks like a normal fault contact, with vertical slickensides marked in the country rock, in particular in Las Trancas Formation (Fig. 10). The Peña de Bernal dacite near the contact is fractured but not brecciated, forming fractures with a peculiar orthogonal arrangement (Fig. 10A). Also observed along the margin of Peña de Bernal dome are sheared marginal zones that form wavy trends in map view between Peña de Bernal and the country rock (Fig. 10B). These observations indicate a rigid (fragile) behavior of the Peña de Bernal dacite along its margins during intrusion, which is discussed below.

## Chemistry

Two whole-rock analyses of the Peña de Bernal rock yielded a dacitic composition following the total-alkalis (TAS) classification of Le Bas et al. (1986) (Fig. 11A). Whole-rock analyses resulted in $\text{SiO}_2 = 67$ wt% and $\text{Na}_2\text{O} + \text{K}_2\text{O} = 6$ wt% (normalized volatile-free, Table 3), which are plotted in the corresponding TAS diagram, together with other volcanic units of the mapped area, such as the Oligocene rhyolitic ignimbrites and the Miocene San Martín andesite (Fig. 11A). For comparison, analyses for other units from the region are plotted too, such as the Miocene andesite of El Zamorano volcano and a Pliocene basaltic andesite from the Cenizas locality, both reported by Aguirre-Díaz and López-Martínez (2001).

Cerro Azul ignimbrite chemistry resulted in $\text{SiO}_2 = 72.5\%$, $\text{Al}_2\text{O}_3 = 11.7\%$, $\text{Na}_2\text{O} = 1.3\%$, $\text{K}_2\text{O} = 5.35\%$, with a high value for loss of ignition (LOI) of 4% (Table 3), probably because
<figure><img src="image_9.png" /><figcaption>Figure 7. $^{40}$Ar/$^{39}$Ar age spectra of Peña de Bernal dacite (sample Qro-8, Tables 1 and 2).</figcaption></figure>

<table><caption>TABLE 2. $^{40}$Ar-$^{39}$Ar DATA OF SAMPLE QRO-8 PLAGIOCLASE</caption><thead><tr><th>Temp.<br>(°C)</th><th>F <sup>39</sup>Ar</th><th><sup>40</sup>Ar*/<sup>39</sup>Ar<sub>K</sub></th><th>1σ</th><th>t<br>(Ma)</th><th>1σ</th><th></th><th><sup>40</sup>Ar*<br>(%)</th><th><sup>40</sup>Ar/<sup>36</sup>Ar</th><th><sup>37</sup>Ar<sub>Ca</sub>/<sup>39</sup>Ar<sub>K</sub></th><th>t<sub>i</sub><br>(Ma)</th><th>t<sub>p</sub><br>(Ma)</th><th>t<sub>c</sub><br>(Ma)</th><th><sup>40</sup>Ar/<sup>36</sup>Ar</th><th>MSWD/n</th></tr></thead><tbody><tr><td>800</td><td>0.1092</td><td>0.51</td><td>0.64</td><td>2.16</td><td>2.71</td><td>†</td><td>3.04</td><td>304.78</td><td>2.21</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>1000</td><td>0.4304</td><td>2.00</td><td>0.07</td><td>8.50</td><td>0.30</td><td>§</td><td>63.77</td><td>815.71</td><td>3.22</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>1300</td><td>0.2955</td><td>2.07</td><td>0.10</td><td>8.79</td><td>0.42</td><td>§</td><td>69.52</td><td>969.39</td><td>2.42</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>1500</td><td>0.1650</td><td>2.25</td><td>0.18</td><td>9.52</td><td>0.76</td><td>§</td><td>45.39</td><td>541.08</td><td>3.02</td><td>8.06 ± 0.41</td><td><b>8.69 ± 0.24</b></td><td>8.05 ± 0.67</td><td>331 ± 36</td><td>0.5/3</td></tr></tbody></table>

Note: Experiments were conducted with a MS-10 mass spectrometer on-line with a Modifications Ltd. Ta-furnace; 0.6122 g of plagioclase concentrate were used for the step-heating experiment. Temp. °C is the temperature used to release argon; t is the age of the individual fraction (it does not include the uncertainty in J); t<sub>c</sub>—integrated age; t<sub>p</sub>—plateau age calculated with the weighted mean of the fractions identified with the symbol §; t<sub>c</sub>—isochron age calculated ignoring the fraction identified with the symbol †; MSWD/n is the goodness of fit and number of points used to calculate the isochron age; all errors are given to 1σ level. J—0.002356 ± 0.000016. Preferred age is highlighted in bold typeface. MSWD—mean square of weighted deviates. To calculate the $^{40}$Ar-$^{39}$Ar ages, the constants recommended by Steiger and Jäeger (1977) were used. The equations presented in York et al. (2004) were used in all the straight line calculations. The internal standard used as irradiation monitor was calibrated with international standards, using the ages given in Renne et al. (1998).
<figure><img src="image_10.png" /><figcaption>Figure 8. (A) Dome foliation of Peña de Bernal in a view to the N, foliation dips from ~40° at borders to vertical in the central zone. (B) Detail of inward-dipping foliation on the western flank in a close-up of the area shown within white frame.</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>Figure 9. (A) Contact between Peña de Bernal dacite and Las Trancas Formation. Cooling joints of dacite shown in photograph are perpendicular to the dome foliation (not evident in photograph), which dips inward. (B) Vertical slickensides at sheared contact surface on Las Trancas Formation (slightly metamorphosed shales due to regional metamorphism). Location of site in photographs is showing in Figure 4.</figcaption></figure>
<figure><img src="image_12.png" /><figcaption>Figure 10. (A) Orthogonal fracturing observed along the inner margin of Peña de Bernal dome, near the intrusive contact with the country rock. Peña de Bernal dacite at its margin is fractured but not brecciated. (B) Plane view of the intensively sheared border zone of Peña de Bernal dacite next to the contact with the country rock caused by forceful intrusion (ruler for scale is 20 × 20 cm).</figcaption></figure>

<figure><img src="image_13.png" /><figcaption>Figure 11. (A) Total alkalis-silica (TAS) plot of Le Bas et al. (1986) of analyzed samples of the Peña de Bernal area and neighboring volcanic rocks reported by Aguirre-Díaz and López-Martínez (2001). (B) SiO₂ versus K₂O plot of same samples shown in (A) with classification of Peccerillo and Taylor (1976).</figcaption></figure>
<table><caption>TABLE 3. MAJOR-ELEMENT CHEMISTRY ON WHOLE-ROCK SAMPLES OF THE BERNAL AREA, QUERETARO</caption><thead><tr><th>Sample</th><th>SM12-05</th><th>SM11-05</th><th>SM13-05</th><th>SM02-05</th><th>SM14-05</th><th>SM03-05</th><th>SM04-05</th><th>SM05-05</th><th>SM06-05</th><th>SM08-05</th><th>Qro-32</th><th>SM10-05</th><th>Qro-8</th><th>Qro-4</th></tr><tr><th>Unit</th><th>Tcai</th><th>Tcaa</th><th>Tcaa</th><th>Tmi</th><th>Tmi</th><th>Tsma</th><th>Tsma</th><th>Tsma</th><th>Tsma</th><th>Tsma</th><th>Tza</th><th>Tpb</th><th>Tpb</th><th>Tca</th></tr><tr><th>Rock type</th><th>R</th><th>BA</th><th>BA</th><th>R</th><th>R</th><th>A</th><th>D</th><th>A</th><th>D</th><th>A</th><th>A</th><th>D</th><th>D</th><th>BA</th></tr></thead><tbody><tr><td colspan="15"><strong>Major elements</strong></td></tr><tr><td>SiO<sub>2</sub></td><td>72.500</td><td>51.800</td><td>51.300</td><td>73.400</td><td>73.900</td><td>60.600</td><td>62.500</td><td>60.400</td><td>62.500</td><td>61.500</td><td>62.00</td><td>67.600</td><td>67.03</td><td>54.92</td></tr><tr><td>TiO<sub>2</sub></td><td>0.130</td><td>1.960</td><td>1.850</td><td>0.140</td><td>0.080</td><td>0.990</td><td>0.840</td><td>1.000</td><td>0.710</td><td>1.000</td><td>0.97</td><td>0.400</td><td>0.46</td><td>2.07</td></tr><tr><td>Al<sub>2</sub>O<sub>3</sub></td><td>11.760</td><td>15.220</td><td>15.200</td><td>13.850</td><td>11.800</td><td>16.150</td><td>15.600</td><td>17.000</td><td>17.350</td><td>16.100</td><td>17.61</td><td>15.400</td><td>16.39</td><td>16.85</td></tr><tr><td>Fe<sub>2</sub>O<sub>3</sub>T</td><td>2.410</td><td>10.420</td><td>10.390</td><td>2.050</td><td>1.350</td><td>5.320</td><td>4.750</td><td>5.380</td><td>4.640</td><td>5.160</td><td>5.36</td><td>2.850</td><td>3.02</td><td>9.25</td></tr><tr><td>MnO</td><td>0.030</td><td>0.180</td><td>0.160</td><td>0.020</td><td>0.020</td><td>0.080</td><td>0.080</td><td>0.080</td><td>0.040</td><td>0.070</td><td>0.09</td><td>0.050</td><td>0.05</td><td>0.14</td></tr><tr><td>MgO</td><td>0.740</td><td>5.070</td><td>5.750</td><td>0.150</td><td>0.690</td><td>2.610</td><td>2.230</td><td>2.740</td><td>0.990</td><td>2.390</td><td>2.54</td><td>1.390</td><td>1.32</td><td>3.71</td></tr><tr><td>CaO</td><td>0.700</td><td>8.120</td><td>7.680</td><td>0.820</td><td>0.550</td><td>6.300</td><td>5.890</td><td>6.720</td><td>5.280</td><td>6.320</td><td>5.66</td><td>3.420</td><td>4.07</td><td>7.48</td></tr><tr><td>Na<sub>2</sub>O</td><td>1.300</td><td>2.810</td><td>2.720</td><td>2.700</td><td>1.320</td><td>3.870</td><td>3.860</td><td>3.650</td><td>3.780</td><td>3.720</td><td>3.70</td><td>3.300</td><td>3.20</td><td>2.87</td></tr><tr><td>K<sub>2</sub>O</td><td>5.350</td><td>1.370</td><td>1.150</td><td>4.640</td><td>4.470</td><td>2.240</td><td>1.490</td><td>2.360</td><td>1.390</td><td>2.220</td><td>1.71</td><td>2.520</td><td>2.27</td><td>1.23</td></tr><tr><td>P<sub>2</sub>O<sub>5</sub></td><td>0.060</td><td>0.600</td><td>0.580</td><td>0.020</td><td>0.020</td><td>0.470</td><td>0.310</td><td>0.450</td><td>0.280</td><td>0.500</td><td>0.28</td><td>0.120</td><td>0.13</td><td>0.33</td></tr><tr><td>LOI</td><td>4.220</td><td>2.580</td><td>3.240</td><td>1.980</td><td>5.040</td><td>0.690</td><td>0.900</td><td>1.290</td><td>2.600</td><td>1.100</td><td>0.58</td><td>2.870</td><td>1.75</td><td>0.32</td></tr><tr><td>Total</td><td>99.200</td><td>100.130</td><td>100.010</td><td>99.770</td><td>99.240</td><td>99.320</td><td>99.850</td><td>99.670</td><td>99.640</td><td>100.080</td><td>100.50</td><td>99.920</td><td>99.72</td><td>99.16</td></tr><tr><td colspan="15"><strong>Major elements (free volatiles calculated with the SINCLAS program; Verma et al., 2002)</strong></td></tr><tr><td>SiO<sub>2</sub></td><td>76.430</td><td>53.470</td><td>53.382</td><td>75.136</td><td>78.508</td><td>61.631</td><td>63.339</td><td>61.587</td><td>64.639</td><td>62.318</td><td>62.243</td><td>69.766</td><td>68.545</td><td>55.888</td></tr><tr><td>TiO<sub>2</sub></td><td>0.137</td><td>2.023</td><td>1.925</td><td>0.143</td><td>0.085</td><td>1.007</td><td>0.851</td><td>1.020</td><td>0.734</td><td>1.013</td><td>0.974</td><td>0.413</td><td>0.478</td><td>2.106</td></tr><tr><td>Al<sub>2</sub>O<sub>3</sub></td><td>12.397</td><td>15.711</td><td>15.817</td><td>14.177</td><td>12.536</td><td>16.425</td><td>17.228</td><td>15.906</td><td>17.944</td><td>16.314</td><td>17.679</td><td>15.894</td><td>16.76</td><td>17.147</td></tr><tr><td>Fe<sub>2</sub>O<sub>3</sub></td><td>1.256</td><td>3.818</td><td>3.744</td><td>1.077</td><td>0.69</td><td>2.342</td><td>2.042</td><td>2.358</td><td>2.017</td><td>2.263</td><td>2.279</td><td>1.340</td><td>1.376</td><td>3.446</td></tr><tr><td>FeO</td><td>1.156</td><td>6.242</td><td>6.339</td><td>0.919</td><td>0.669</td><td>2.761</td><td>2.494</td><td>2.814</td><td>2.503</td><td>2.668</td><td>2.792</td><td>1.441</td><td>1.540</td><td>5.370</td></tr><tr><td>MnO</td><td>0.032</td><td>0.186</td><td>0.166</td><td>0.020</td><td>0.021</td><td>0.081</td><td>0.081</td><td>0.082</td><td>0.041</td><td>0.071</td><td>0.091</td><td>0.052</td><td>0.055</td><td>0.146</td></tr><tr><td>MgO</td><td>0.780</td><td>5.233</td><td>5.983</td><td>0.154</td><td>0.733</td><td>2.654</td><td>2.260</td><td>2.794</td><td>1.024</td><td>2.422</td><td>2.550</td><td>1.434</td><td>1.354</td><td>3.775</td></tr><tr><td>CaO</td><td>0.738</td><td>8.382</td><td>7.992</td><td>0.839</td><td>0.584</td><td>6.407</td><td>5.969</td><td>6.852</td><td>5.461</td><td>6.404</td><td>5.682</td><td>3.530</td><td>4.162</td><td>7.612</td></tr><tr><td>Na<sub>2</sub>O</td><td>1.370</td><td>2.900</td><td>2.830</td><td>2.764</td><td>1.402</td><td>3.936</td><td>3.912</td><td>3.722</td><td>3.909</td><td>3.770</td><td>3.714</td><td>3.406</td><td>3.272</td><td>2.921</td></tr><tr><td>K<sub>2</sub>O</td><td>5.640</td><td>1.414</td><td>1.197</td><td>4.750</td><td>4.749</td><td>2.278</td><td>1.510</td><td>2.406</td><td>1.438</td><td>2.250</td><td>1.717</td><td>2.601</td><td>2.321</td><td>1.252</td></tr><tr><td>P<sub>2</sub>O<sub>5</sub></td><td>0.063</td><td>0.619</td><td>0.604</td><td>0.020</td><td>0.021</td><td>0.478</td><td>0.314</td><td>0.459</td><td>0.290</td><td>0.507</td><td>0.279</td><td>0.124</td><td>0.134</td><td>0.337</td></tr><tr><td colspan="15"><strong>CIPW normative minerals (calculated with the SINCLAS program; Verma et al., 2002)</strong></td></tr><tr><td>Q</td><td>43.891</td><td>7.317</td><td>7.294</td><td>38.636</td><td>49.639</td><td>14.722</td><td>19.341</td><td>14.872</td><td>23.171</td><td>16.820</td><td>18.057</td><td>30.214</td><td>29.614</td><td>12.590</td></tr><tr><td>Or</td><td>33.330</td><td>8.356</td><td>7.074</td><td>28.071</td><td>28.065</td><td>13.462</td><td>8.924</td><td>14.218</td><td>8.498</td><td>13.297</td><td>10.147</td><td>15.371</td><td>13.716</td><td>7.399</td></tr><tr><td>Ab</td><td>11.592</td><td>24.539</td><td>23.946</td><td>23.388</td><td>11.863</td><td>33.305</td><td>33.102</td><td>31.494</td><td>33.077</td><td>31.900</td><td>31.427</td><td>28.82</td><td>27.687</td><td>24.716</td></tr><tr><td>An</td><td>3.250</td><td>25.675</td><td>26.919</td><td>4.032</td><td>2.760</td><td>20.421</td><td>24.988</td><td>19.587</td><td>25.198</td><td>20.946</td><td>26.366</td><td>16.703</td><td>19.773</td><td>29.977</td></tr><tr><td>C</td><td>2.847</td><td>—</td><td>5.328</td><td>3.011</td><td>4.078</td><td>—</td><td>—</td><td>—</td><td>0.722</td><td>—</td><td>0.047</td><td>1.354</td><td>—</td><td>—</td></tr><tr><td>Di</td><td>—</td><td>9.554</td><td>7.062</td><td>—</td><td>—</td><td>6.562</td><td>2.056</td><td>9.081</td><td>—</td><td>5.983</td><td>—</td><td>—</td><td>1.618</td><td>—</td></tr><tr><td>Hy</td><td>2.861</td><td>13.745</td><td>17.219</td><td>0.983</td><td>2.383</td><td>5.111</td><td>6.285</td><td>4.328</td><td>4.343</td><td>4.675</td><td>8.155</td><td>4.525</td><td>4.376</td><td>11.044</td></tr><tr><td>OI</td><td>—</td><td>5.535</td><td>—</td><td>1.561</td><td>1.000</td><td>—</td><td>2.96</td><td>3.418</td><td>2.924</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>Mt</td><td>1.821</td><td>—</td><td>5.428</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>3.281</td><td>3.304</td><td>1.943</td><td>1.995</td><td>4.996</td></tr><tr><td>Hm</td><td>—</td><td>3.842</td><td>—</td><td>0.272</td><td>—</td><td>1.912</td><td>1.616</td><td>1.937</td><td>1.394</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>II</td><td>0.260</td><td>1.434</td><td>3.666</td><td>0.046</td><td>0.161</td><td>1.107</td><td>0.728</td><td>1.063</td><td>0.672</td><td>1.324</td><td>1.850</td><td>0.784</td><td>0.908</td><td>4.000</td></tr><tr><td>Ap</td><td>0.146</td><td>0.003</td><td>1.399</td><td>0.023</td><td>0.049</td><td>0.007</td><td>0.003</td><td>0.005</td><td>0.005</td><td>0.175</td><td>0.640</td><td>0.287</td><td>0.310</td><td>0.781</td></tr></tbody></table>

Note: Rock type refers to the chemical classification according to total-alkalis versus silica diagram (Le Bas et al., 1986); Fe₂O₃-T-Fe₂O₃ as total iron; BA—basaltic-andesite; A—andesite; R—rhyolite; D—dacite. Unit refers to lithostratigraphic unit according to the geological map and stratigraphic column shown in Figures 4 and 5. Tcai—Cerro Azul ignimbrite, Tcaa—Cerro Azul andesite, Tmi—El Matón ignimbrite, Tsma—San Martín andesite, Tza—Zamorano andesite, Tca—Cenizas andesite, Tpb—Pena de Bernal dacite. Qro-4 and Qro-32 are samples reported by Aguirre-Díaz and López-Martinez (2001); these two and sample Qro-8 were analyzed at the laboratory of Instituto de Geologia at Universidad Nacional Autónoma de México (UNAM) using X-ray diffraction technique (Lozano-Santacruz and Bernal, 2005). The series SM samples were analyzed at the University of Brest, France, using inductively coupled plasma-atomic emission spectroscopy (Cotten et al., 1995).
we analyzed partly devitrified pumice fragments, indicating a rhyolitic composition in the TAS diagram (Fig. 11A). Major-element chemistry of Cerro Azul andesite yielded SiO₂ = 58.0%, Al₂O₃ = 15%, TiO₂ = 0.82%, Fe₂O₃ = 4.6%, Na₂O = 2%, and K₂O = 1.37% (Table 3), corresponding to a basaltic andesite composition in the TAS classification. Chemical analyses of El Matón ignimbrite yielded SiO₂ ~73%, Al₂O₃ = 11%–13%, Na₂O = 1.3%–2.7%, and K₂O ~4.4%. It classifies as rhyolite in the TAS diagram (Fig. 11A). However, it should be noted that the analyzed pumice was devitrified. Chemical data of two analyzed samples from San Martín andesite show SiO₂ = 60.4%–64.9%, Al₂O₃ = 15.4%–17.0%, TiO₂ = 0.4%–1.0%, Fe₂O₃ = 2.8%–5.3%, Na₂O ~3.6%, and K₂O = 1.3%–2.5% (Table 3). These values result in a TAS classification of andesite to dacite (Fig. 11A).

Peña de Bernal dacite as well as other volcanic units sampled from the Bernal area all are calc-alkaline, low-K silicic rocks, as shown in the SiO₂-K₂O plot with the classification of Peccerillo and Taylor (1976; Fig. 11B). According to these results, the suite of Oligocene, Miocene, and Pliocene rocks, and in particular the Peña de Bernal dacite, corresponds to a continental margin subduction zone setting (Wilson, 1989). This tectonic setting agrees well with the fact that the Peña de Bernal lies within the Mexican Volcanic Belt, which has been described as a volcanic province formed during a continental margin subduction zone tectonic regime (Aguirre-Díaz et al., 1998; Siebe et al., 2006).

Peña de Bernal dacite is relatively poorer in silica and alkali content with respect to other volcanic rocks in the area, such as the rhyolitic ignimbrites of El Matón and Cerro Azul. On the other hand, it is relatively more evolved than the neighboring andesites, such as those erupted from the San Martín volcano and the Cenizas andesite in the Amazcala caldera area (Fig. 12). This may be explained by the crystal-rich content of Bernal dacite (up to 80 vol%) with respect to these rhyolites, because crystals and in particular plagioclase tend to lower the silica and Na₂O + K₂O contents in the whole-rock composition. The glass composition described below confirms this point.

Microprobe analyses of glass show a higher Na₂O + K₂O content than the whole-rock composition, yielding a trachytic composition of SiO₂ = 65 wt% and Na₂O + K₂O = 13 wt% (Table 4; Fig. 13). Mineral chemistry (Table 5) indicates a wide range in plagioclase composition from An₂₀ to An₆₀ (Fig. 14). This range is confirmed in textures observed in large plagioclases with oscillatory zoning and fritted cores

(Fig. 6D). They generally have Ca-rich cores and Na-rich margins, suggesting a crystallization that followed the liquid line of descent of the Bowen's series (Wilcox, 1979). Two microprobe analyses confirm that orthopyroxene compositions correspond to enstatite (Fig. 14). Energy dispersive spectroscopy (EDS) analy-

ses (not shown in Table 5) confirm the presence of biotite and hornblende and of the accessory phases, apatite and zircon.

Oxygen fugacity and corresponding equilibration temperatures were calculated based upon microprobe analyses on Fe-Ti oxides and pyroxene (Table 5). The temperature and

<figure><img src="image_14.png" /><figcaption>Figure 12. Harker plots (silica versus X component) of volcanic rocks of the Peña de Bernal area and vicinity (Cenizas and El Zamorano andesites).</figcaption></figure>
<table><caption>TABLE 4. GLASS CHEMISTRY OF PEÑA DE BERNAL DACITE</caption><thead><tr><th>Sample</th><th>Qro-8A</th><th>Qro-8B</th><th>Qro-8C</th></tr></thead><tbody><tr><td>SiO<sub>2</sub></td><td>65.41</td><td>65.49</td><td>64.88</td></tr><tr><td>TiO<sub>2</sub></td><td>—</td><td>—</td><td>—</td></tr><tr><td>Al<sub>2</sub>O<sub>3</sub></td><td>18.82</td><td>19.41</td><td>19.07</td></tr><tr><td>FeO</td><td>0.02</td><td>0.06</td><td>0.20</td></tr><tr><td>MnO</td><td>—</td><td>—</td><td>—</td></tr><tr><td>MgO</td><td>—</td><td>—</td><td>—</td></tr><tr><td>CaO</td><td>0.38</td><td>0.45</td><td>0.43</td></tr><tr><td>Na<sub>2</sub>O</td><td>3.59</td><td>4.61</td><td>3.51</td></tr><tr><td>K<sub>2</sub>O</td><td>10.06</td><td>8.23</td><td>10.10</td></tr><tr><td>Cr<sub>2</sub>O<sub>3</sub></td><td>0.01</td><td>—</td><td>—</td></tr><tr><td>BaO</td><td>1.73</td><td>1.74</td><td>1.80</td></tr><tr><td>Total</td><td>100.01</td><td>100.00</td><td>100.00</td></tr></tbody></table>

Note: Glass analyses by scanning electron microprobe using a Jeol-100 at the Paris 6 University-Jussieu Place performed by G. Aguirre-Díaz.

<figure><img src="image_15.png" /><figcaption>Figure 13. Total alkalis-silica (TAS) plot of Le Bas et al. (1986) of glass analyses obtained with the microprobe from samples of Peña de Bernal dacite.</figcaption></figure>

oxygen fugacity (T-fO₂) parameters were calculated using the quartz-ilmenite-fayalite system proposed by Lindsley and Frost (1992) and the computer program QUILF (Andersen et al., 1993). Plotting the Peña de Bernal dacite data in a T-fO₂ diagram together with several buffers, a range of values is observed from T °C = 711.5, log fO₂ = -15.53 to T °C = 377, log fO₂ = -27.70, with several intermediate values (Fig. 15 and Table 6). The results plot between the nickel-nickel-oxide (NNO) and the hematite-

magnetite (HM) buffers, indicating relatively high fO₂ conditions, which are common in subduction zone settings.

DISCUSSION

The Peña de Bernal area includes rocks of three major geologic provinces of Mexico: the Mesozoic Sierra Madre Oriental fold-thrust belt, the mid-Tertiary Sierra Madre Occidental volcanic province, and the Neogene-Quaternary

Mexican Volcanic Belt. Focusing on the volcanic rocks, these include 31.5-29.4 Ma ignimbrites, 29.4 Ma andesites, 11-10 Ma andesites, and the 8.7 Ma Peña de Bernal dacite, which represents the youngest volcanic event in the mapped area (Fig. 4 and Table 1). The Oligocene ignimbrites and andesites crop out to the north of Peña de Bernal. At the base of this sequence is a layered succession of reworked pyroclastic deposits that in turn overlies continental red or green sandstones. The lower layered sequence resembles

<table><caption>TABLE 5. MINERAL CHEMISTRY FOR PEÑA DE BERNAL SAMPLES</caption><thead><tr><th>Sample</th><th>1</th><th>2</th><th>3</th><th>4</th><th>5</th><th>6</th><th>7</th><th>8</th><th>9</th><th>10</th><th>11</th><th>12</th><th>13</th><th>14</th><th>15</th><th>16</th><th>17</th></tr></thead><tbody><tr><td>SiO<sub>2</sub></td><td>50.90</td><td>51.83</td><td>52.97</td><td>53.54</td><td>56.27</td><td>56.16</td><td>53.88</td><td>57.63</td><td>58.22</td><td>61.94</td><td>61.89</td><td>61.77</td><td>56.51</td><td>60.31</td><td>54.44</td><td>54.08</td><td>54.17</td></tr><tr><td>TiO<sub>2</sub></td><td>0.09</td><td>0.03</td><td>0.09</td><td>0.07</td><td>—</td><td>—</td><td>0.16</td><td>0.05</td><td>—</td><td>—</td><td>—</td><td>—</td><td>0.12</td><td>0.03</td><td>0.04</td><td>0.10</td><td>0.06</td></tr><tr><td>Al<sub>2</sub>O<sub>3</sub></td><td>30.29</td><td>29.98</td><td>28.75</td><td>28.82</td><td>27.03</td><td>27.27</td><td>28.20</td><td>26.10</td><td>25.31</td><td>23.16</td><td>22.54</td><td>22.90</td><td>26.60</td><td>24.78</td><td>28.22</td><td>28.00</td><td>28.00</td></tr><tr><td>FeO</td><td>0.35</td><td>0.23</td><td>0.36</td><td>0.18</td><td>0.29</td><td>0.38</td><td>0.41</td><td>0.29</td><td>0.32</td><td>0.35</td><td>0.39</td><td>0.32</td><td>0.02</td><td>0.04</td><td>0.41</td><td>0.20</td><td>0.33</td></tr><tr><td>MnO</td><td>—</td><td>0.15</td><td>—</td><td>0.09</td><td>0.05</td><td>0.02</td><td>0.02</td><td>0.06</td><td>—</td><td>—</td><td>0.01</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>0.01</td></tr><tr><td>MgO</td><td>0.06</td><td>0.06</td><td>—</td><td>0.04</td><td>—</td><td>—</td><td>0.03</td><td>0.04</td><td>0.03</td><td>—</td><td>—</td><td>—</td><td>—</td><td>0.00</td><td>0.03</td><td>0.01</td><td>—</td></tr><tr><td>CaO</td><td>13.40</td><td>12.37</td><td>11.62</td><td>11.15</td><td>9.50</td><td>9.16</td><td>10.04</td><td>8.01</td><td>5.90</td><td>4.14</td><td>4.64</td><td>4.31</td><td>8.71</td><td>6.39</td><td>10.74</td><td>10.76</td><td>10.76</td></tr><tr><td>Na<sub>2</sub>O</td><td>3.91</td><td>4.61</td><td>4.91</td><td>5.24</td><td>5.83</td><td>6.04</td><td>5.69</td><td>6.81</td><td>7.54</td><td>8.90</td><td>8.41</td><td>8.18</td><td>6.70</td><td>7.39</td><td>5.50</td><td>5.62</td><td>5.65</td></tr><tr><td>K<sub>2</sub>O</td><td>0.09</td><td>0.17</td><td>0.21</td><td>0.23</td><td>0.34</td><td>0.33</td><td>0.29</td><td>0.39</td><td>0.46</td><td>0.88</td><td>0.78</td><td>0.98</td><td>0.18</td><td>0.63</td><td>0.19</td><td>0.22</td><td>0.20</td></tr><tr><td>Cr<sub>2</sub>O<sub>3</sub></td><td>—</td><td>0.02</td><td>0.10</td><td>—</td><td>0.03</td><td>0.04</td><td>0.03</td><td>0.04</td><td>—</td><td>—</td><td>0.12</td><td>—</td><td>0.09</td><td>—</td><td>—</td><td>0.11</td><td>0.02</td></tr><tr><td>Total</td><td>99.09</td><td>99.45</td><td>99.01</td><td>99.36</td><td>99.34</td><td>99.40</td><td>98.75</td><td>99.42</td><td>97.78</td><td>99.37</td><td>98.78</td><td>98.46</td><td>98.93</td><td>99.57</td><td>99.57</td><td>99.10</td><td>99.20</td></tr><tr><td>Sample</td><td>18</td><td>19</td><td>20</td><td>21</td><td>22</td><td>23</td><td>24</td><td>25</td><td>26</td><td>27</td><td>28</td><td>29</td><td>30</td><td>31</td><td>32</td><td>33</td><td>34</td></tr><tr><td>SiO<sub>2</sub></td><td>53.36</td><td>53.74</td><td>54.42</td><td>53.53</td><td>55.99</td><td>56.08</td><td>55.86</td><td>53.90</td><td>54.56</td><td>53.88</td><td>55.80</td><td>55.97</td><td>55.49</td><td>53.93</td><td>55.53</td><td>51.96</td><td>51.36</td></tr><tr><td>TiO<sub>2</sub></td><td>0.07</td><td>—</td><td>—</td><td>0.05</td><td>0.06</td><td>0.11</td><td>—</td><td>0.02</td><td>—</td><td>0.05</td><td>0.05</td><td>—</td><td>—</td><td>0.08</td><td>0.09</td><td>0.26</td><td>0.21</td></tr><tr><td>Al<sub>2</sub>O<sub>3</sub></td><td>28.64</td><td>28.48</td><td>27.98</td><td>28.14</td><td>26.97</td><td>27.31</td><td>27.49</td><td>28.01</td><td>27.89</td><td>28.36</td><td>26.97</td><td>27.21</td><td>27.64</td><td>27.93</td><td>27.10</td><td>1.06</td><td>0.75</td></tr><tr><td>FeO</td><td>0.28</td><td>0.34</td><td>0.10</td><td>0.05</td><td>0.31</td><td>0.18</td><td>0.22</td><td>0.21</td><td>0.30</td><td>0.32</td><td>0.33</td><td>0.18</td><td>0.35</td><td>0.27</td><td>0.16</td><td>23.05</td><td>23.57</td></tr><tr><td>MnO</td><td>—</td><td>—</td><td>—</td><td>0.08</td><td>—</td><td>0.03</td><td>0.01</td><td>—</td><td>0.04</td><td>—</td><td>0.05</td><td>0.05</td><td>—</td><td>—</td><td>0.05</td><td>0.81</td><td>0.90</td></tr><tr><td>MgO</td><td>—</td><td>0.01</td><td>0.05</td><td>—</td><td>0.02</td><td>—</td><td>0.01</td><td>—</td><td>—</td><td>—</td><td>0.04</td><td>—</td><td>—</td><td>0.01</td><td>0.02</td><td>21.66</td><td>21.49</td></tr><tr><td>CaO</td><td>10.54</td><td>10.80</td><td>10.31</td><td>10.24</td><td>8.99</td><td>9.19</td><td>9.64</td><td>10.08</td><td>10.14</td><td>9.94</td><td>9.27</td><td>9.96</td><td>10.00</td><td>10.07</td><td>9.52</td><td>0.27</td><td>0.33</td></tr><tr><td>Na<sub>2</sub>O</td><td>5.49</td><td>5.67</td><td>6.02</td><td>5.48</td><td>6.42</td><td>6.10</td><td>6.02</td><td>5.55</td><td>5.58</td><td>5.48</td><td>6.52</td><td>5.76</td><td>5.58</td><td>5.62</td><td>6.32</td><td>—</td><td>—</td></tr><tr><td>K<sub>2</sub>O</td><td>0.17</td><td>0.20</td><td>0.20</td><td>0.29</td><td>0.41</td><td>0.26</td><td>0.33</td><td>0.22</td><td>0.29</td><td>0.28</td><td>0.40</td><td>0.28</td><td>0.33</td><td>0.30</td><td>0.42</td><td>—</td><td>—</td></tr><tr><td>Cr<sub>2</sub>O<sub>3</sub></td><td>—</td><td>—</td><td>—</td><td>—</td><td>0.02</td><td>—</td><td>0.05</td><td>—</td><td>—</td><td>0.10</td><td>0.05</td><td>—</td><td>—</td><td>0.10</td><td>—</td><td>0.04</td><td>0.12</td></tr><tr><td>Total</td><td>98.55</td><td>99.24</td><td>99.08</td><td>97.86</td><td>99.19</td><td>99.26</td><td>99.63</td><td>97.99</td><td>98.80</td><td>98.41</td><td>99.48</td><td>99.41</td><td>99.39</td><td>98.31</td><td>99.21</td><td>99.11</td><td>98.73</td></tr><tr><td>Sample</td><td>35</td><td>36</td><td>37</td><td>38</td><td>39</td><td>40</td><td>41</td><td>42</td><td>43</td><td>44</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>SiO<sub>2</sub></td><td>—</td><td>0.01</td><td>0.01</td><td>0.07</td><td>0.05</td><td>0.46</td><td>—</td><td>0.30</td><td>0.16</td><td>0.05</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>TiO<sub>2</sub></td><td>47.35</td><td>45.80</td><td>45.58</td><td>6.37</td><td>0.93</td><td>1.64</td><td>0.88</td><td>7.25</td><td>0.43</td><td>0.91</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Al<sub>2</sub>O<sub>3</sub></td><td>0.04</td><td>0.07</td><td>0.07</td><td>4.24</td><td>0.78</td><td>1.12</td><td>1.13</td><td>2.12</td><td>0.82</td><td>0.75</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>FeO</td><td>48.44</td><td>46.71</td><td>42.73</td><td>79.33</td><td>91.30</td><td>86.72</td><td>89.00</td><td>79.63</td><td>90.86</td><td>90.09</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>MnO</td><td>0.96</td><td>1.22</td><td>2.35</td><td>0.27</td><td>—</td><td>0.07</td><td>0.09</td><td>0.64</td><td>0.11</td><td>0.11</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>MgO</td><td>1.32</td><td>1.55</td><td>2.00</td><td>0.98</td><td>0.23</td><td>0.24</td><td>0.30</td><td>0.63</td><td>0.30</td><td>0.16</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>CaO</td><td>0.02</td><td>0.65</td><td>0.09</td><td>0.03</td><td>0.01</td><td>0.09</td><td>0.05</td><td>0.05</td><td>0.04</td><td>0.05</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Na<sub>2</sub>O</td><td>0.01</td><td>—</td><td>0.02</td><td>0.05</td><td>—</td><td>—</td><td>—</td><td>0.08</td><td>0.03</td><td>—</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>K<sub>2</sub>O</td><td>—</td><td>—</td><td>—</td><td>0.04</td><td>0.05</td><td>—</td><td>—</td><td>—</td><td>—</td><td>0.01</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Cr<sub>2</sub>O<sub>3</sub></td><td>0.06</td><td>0.02</td><td>0.01</td><td>0.12</td><td>0.03</td><td>—</td><td>0.03</td><td>—</td><td>—</td><td>—</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Total</td><td>98.20</td><td>96.03</td><td>92.86</td><td>91.50</td><td>93.38</td><td>90.34</td><td>91.48</td><td>90.70</td><td>92.75</td><td>92.13</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

Note: Samples 1 to 32 are feldspar analyses, 33 and 34 are pyroxene analyses, and 35 to 44 are Fe-Ti oxides analyses. Microprobe analyses performed at Paris 6 University-Jussieu Place by Gerardo Aguirre-Díaz.
Figure 14. (A) Ternary plot showing feldspar composition and classification from micro-probe analyses of plagioclases of Peña de Bernal dacite. (B) Quaternary plot of two pyroxene analyses and classification.

<figure><img src="image_16.png" /></figure>

<figure><img src="image_17.png" /></figure>

<figure><img src="image_18.png" /></figure>

Figure 15. Temperature versus oxygen fugacity plot of Peña de Bernal dacite, obtained from Fe-Ti oxide pairs microprobe geochemistry. See data in Table 5. NNO—nickel-nickel-oxide.

mid-Tertiary red continental deposits described in other places of central Mexico, such as the upper part of the Conglomerado Rojo of Guanajuato Formation, apparently of Eocene age (Edwards, 1955; Echegoyén et al., 1970). Furthermore, the contact between Guanajuato's red conglomerate with overlying green Loseros Formation (Echegoyén et al., 1970) is similar to the succession of continental red sandstones and Cerro Azul epiclastic deposits of the Bernal area. The ignimbrites include sanidine, quartz, and minor biotite, and they can be correlated with similar 30–31 Ma ignimbrites reported nearby and to the WNW of the Bernal area and next to, but unrelated to, the Amazcala caldera (Aguirre-Díaz and López-Martínez, 2001).

San Martín andesite was derived from the neighboring San Martín volcano 3 km to the west of Bernal (Fig. 4). It was dated at 11–10 Ma (Table 1) and is one more of several middle Miocene andesitic stratovolcanoes at the north-central border of the Mexican Volcanic Belt, including San Pedro, Palo Huérfano, Chichimequillas, La Joya, and El Zamorano (Valdez-Moreno et al., 1998; Verma and Carrasco-Núñez, 2003; Aguirre-Díaz, 2008).

El Zamorano volcano includes a spine-type lava dome in its crater that is similar to the Peña de Bernal. It is a light-gray, porphyritic rock, almost completely crystalline, and with phenocrysts of plagioclase, hornblende, and quartz. Dome foliation jointing varies from vertical at the central part of the dome to subhorizontal and inward dipping at the borders, similar to the foliation observed in Peña de Bernal. The dacitic dome intrudes and is in high angular unconformity with surrounding outward-dipping lavas of the crater of the El Zamorano volcano. With an altitude of 3360 m above sea level, this spine forms the highest point of the El Zamorano volcano, which is also the highest peak of the State of Querétaro.

Peña de Bernal dacite is interpreted in this work as a lava dome with intrusive contacts with the adjacent rocks. This interpretation is based upon the porphyritic texture and fabric of the rock and the dome-type foliation geometry (e.g., Duffield et al., 1995). A thermal aureole is not present in country rock surrounding the Peña de Bernal dome. The slight metamorphism of Las Trancas Formation was caused by the earlier regional tectonic event of the Laramide orogeny (Suter, 1987) and not by the intrusion of Peña de Bernal dome.

Peña de Bernal dacite dome is intensively fractured along its margins, but it is not brecciated. Shearing of Peña de Bernal dacite along the contact surface formed an intensively fractured rock along some margins, and vertical slickensides were also observed in some places due
<table><caption>TABLE 6. CALCULATED TEMPERATURE (T) AND OXYGEN FUGACITY ($fO_2$) FOR PEÑA DE BERNAL SAMPLES</caption><thead><tr><th>Pair</th><th>Log $fO_2$</th><th>T (°C)</th></tr></thead><tbody><tr><td>01</td><td>-15.53</td><td>711</td></tr><tr><td>02</td><td>-27.70</td><td>377</td></tr><tr><td>03</td><td>-22.75</td><td>444</td></tr><tr><td>04</td><td>-25.56</td><td>430</td></tr><tr><td>05</td><td>-22.67</td><td>377</td></tr><tr><td>06</td><td>-22.65</td><td>445</td></tr><tr><td>07</td><td>-20.38</td><td>512</td></tr><tr><td>08</td><td>-14.60</td><td>747</td></tr><tr><td>09</td><td>-24.76</td><td>391</td></tr><tr><td>10</td><td>-22.76</td><td>443</td></tr><tr><td>11</td><td>-23.54</td><td>431</td></tr><tr><td>12</td><td>-25.46</td><td>381</td></tr><tr><td>13</td><td>-15.51</td><td>718</td></tr><tr><td>14</td><td>-23.26</td><td>437</td></tr><tr><td>15</td><td>-14.76</td><td>766</td></tr></tbody></table>

Note: Mineral microprobe analyses of Fe-Ti oxides analyzed at Paris 6 University-Jussieu Place, by G. Aguirre-Díaz.

to this shearing. These observations, together with the highly crystallized nature of the rock, suggest a forceful intrusion of the dome when it was nearly solid (80% crystalline) and relatively cold (Fig. 16). This is confirmed by geothermometry obtained by microprobe analyses on minerals, which yielded temperatures as high as 711 °C and as low as 377 °C, following the hematite-magnetite buffer (Fig. 15). Therefore, the crystal-rich dacitic magma was almost totally crystallized and relatively cold when it was emplaced at or near surface.

Peña de Bernal dacite was formed through three simultaneous pulses, each constructing a plug; these are the Main plug forming Peña de Bernal peak, the Tortuga plug, and the Northern plug. The Tortuga and Northern plugs are much smaller than the main one. There is no evidence for magmas reaching a vent and flowing over

the surface as lava coulees. Instead, Peña de Bernal dacite was emplaced as an endogenous lava dome; that is, the dome grew as a dacitic magma batch that was emplaced and cooled below but near surface.

In summary, Peña de Bernal dacite was apparently intruded as a spine dome, following Blake's (1990) classification, and it probably has retained its peculiar steep peak shape since its formation. However, in the case of Peña de Bernal, the spine dome did not form inside a crater of a volcano, as normally occurs in similar cases, for instance, El Zamorano, Sangangüey, or Tequila volcanoes (Nelson and Carmichael, 1984; Wallace and Carmichael, 1994). Instead, Peña de Bernal was emplaced through a marine Mesozoic sedimentary sequence, and without any evident relationship with a preexisting volcano. There is no evidence that Peña de Bernal vented to the surface, either in volcanic deposits of appropriate age or in dikes that might have fed surface eruptions. The presence of Miocene, relatively intact volcanoes nearby, such as El Zamorano and San Martin, suggests that erosion would not have removed a volcanic edifice above Peña de Bernal dacite, if such had existed.

## CONCLUSION

Peña de Bernal is a prominent igneous peak, and it emphasizes the particular geographic and cultural setting where it lies, at the entrance of UNESCO's *Sierra Gorda Biosphere Reserve*, as part of UNESCO's *Intangible Cultural Heritage Patrimony List*, and next to the *Magic Town* of Bernal. With a measured height of 433 m from its base, it is apparently the highest monolith in the world. Peña de Bernal includes three plugs

with a general N40°E trend and covers an area of ~4.8 km². We conclude that Peña de Bernal was emplaced as a spine-type, endogenous dome, at ca. 8.7 ± 0.2 Ma. Compositionally it corresponds to a dacitic lava with $SiO_2$ = 67 wt% and $Na_2O + K_2O$ = 6 wt%. The rock is 80 vol% crystalline and only 20 vol% glass. It is porphyritic with plagioclase, pyroxene, hornblende, biotite, quartz, magnetite, ilmenite, apatite, and zircon. On the basis of temperatures as low as 377 °C, the nearly crystalline texture, and the fragile intrusion observations, such as slickensides at the contact with the country rock, it is concluded that Peña de Bernal rock was emplaced by forceful intrusion when nearly solid.

## ACKNOWLEDGMENTS

We thank the reviewers Don F. Parker, Fred W. McDowell, and Geosphere's associate editor Megan Elwood Madden and science editor Carol Frost for their comments and recommendations that substantially improved the final version. We are grateful to Isaac A. Farraz-Montes and Azalea J. Ortiz-Rodríguez for their help with figures and digital elevation data analysis. We thank Compañía Minera Peña de Bernal, S.A. de C.V. for allowing the use of some data from the unpublished internal report Estudio geológico-estratigráfico y estructural, en el área del proyecto San Martín performed in 2006 by researchers of Universidad Autónoma de San Luis Potosí (UASLP). We thank Andrea Di Muro for arranging access to the microprobe of the Laboratoire de Pétrologie Modélisation des Matériaux et Processus, of Universite Pierre et Marie Curie, Paris 6, France. This work was partially supported by Consejo Nacional de Ciencia y Tecnología (CONACYT) and Programa de Apoyo a Proyectos de Investigación e Innovación Tecnológica—Universidad Nacional Autónoma de México (PAPIIT-UNAM) grants to the first author, numbers IN-114606, IN-106109, IN-112312, P-46005F and 33084T.

## REFERENCES CITED

* Aguirre-Díaz, G.J., 2008, Historia volcánica del entorno del valle de Querétaro, in *Lozano-Guzmán, A.*, and Arzate, J., eds., El Valle de Querétaro y su Geoentorno, Querétaro: Publicación Especial Consejo Nacional de Ciencia y Tecnología del Estado de Querétaro [CONCYTEQ], Querétaro, México y Centro de Ciencias, UNAM-Juriquilla, p. 27-43.

* Aguirre-Díaz, G.J., and López-Martínez, M., 2001, The Amazcala caldera, Querétaro, central Mexican Volcanic Belt, México: Geology and geochronology: Journal of Volcanology and Geothermal Research, v. 111, p. 203-218, doi:10.1016/S0377-0273(01)00227-X.

* Aguirre-Díaz, G.J., Ferrari, L., Nelson, S.A., Carrasco-Núñez, G., López-Martínez, M., and Urrutia-Fucugauchi, J., 1998, El Cinturón Volcánico Mexicano: Un Nuevo Proyecto Multidisciplinario: Unión Geofísica Mexicana: Geos, v. 18, p. 131-138.

* Aguirre-Díaz, G.J., Labarthe-Hernández, G., Tristán-González, M., Nieto-Obregón, J., and Gutiérrez-Palomares, I., 2008, in *Gottsmann, J.*, and Martí, J., eds., Ignimbrite Flare-Up and Graben-Calderas of the Sierra Madre Occidental, Mexico, Caldera Volcanism: Analysis, Modelling and Response: Amsterdam, Elsevier, Developments in Volcanology, v. 10, p. 143-180.

* Andersen, D.J., Lindsley, L.D.H., and Davidson, P.M., 1993, QUILF, A Pascal program to assess equilibria among Fe-Mg-Mn-Ti oxides, pyroxenes, olivine, and quartz:

<figure><img src="image_19.png" /><figcaption>Figure 16. Schematic model showing the emplacement evolution of Peña de Bernal dome, from (A) intrusion of a dacitic magma body at depth, (B) ascent and crystallization of magma body to form an endogenous lava dome near the surface, and (C) forceful intrusion of a nearly solid plug as a spine-type lava dome. From A to B, the magma is richer in crystal content, reaching 80 vol% in the last stage when it completely cools down. Concentric dome foliation typical of endogenous lava domes is indicated by the vertical to subhorizontal white light lines.</figcaption></figure>
Computers & Geosciences, v. 19, p. 1333-1350, doi: 10.1016/0098-3004(93)90033-2.

Bartolini, C., Wilson, J.L., and Lawton, T.F., eds., 1999, Mesozoic Sedimentary and Tectonic History of North Central Mexico: Boulder, Colorado, Geological Society of America Special Paper 340, 380 p.

Bellon, H., Quoc-Buï, N., Chaumont, J., and Philpott, J.C., 1981, Implantation ionique d'argon dans une cible support: Application au traçage isotopique de l'argon contenu dans les minéraux et les roches: Les Comptes Rendus de l'Académie de Paris, v. 292, p. 977-980.

Blake, S., 1990, Viscoplastic models of lava domes, in Fink, J. H., ed., Lava Flow and Domes: Emplacement, Mechanisms and Hazard Implications: Berlin, Springer-Verlag, p. 88-126.

Carrillo-Martínez, M., 1990, Geometría estructural de la Sierra Madre Oriental entre Peñamiller y Jalpan, Estado de Querétaro: Sociedad Geológica Mexicana, Convenición Geológica Nacional, 6, México, D.F., Libro-guía de la excursión geológica a la región de Zimapán y áreas circundantes, Estados de Hidalgo y Querétaro, p. 1-20.

Carrillo-Martínez, M., 1998, Hoja Zimapán, Carta Geológica de México Serie 1:100,000, Hoja 14Q-e(7): Instituto de Geología, Universidad Nacional Autónoma de México, 1 sheet.

Carrillo-Martínez, M., and Suter, M., 1990, Tectónica de los alrededores de Zimapán, Hidalgo y Querétaro: Sociedad Geológica Mexicana, Convenición Geológica Nacional, 6, México, D.F., Libro-guía de la excursión geológica a la región Zimapán y áreas circundantes, estados Hidalgo y Querétaro, p. 1-20.

Cerca-Martínez, L.M., Aguirre-Díaz, G., and López-Martínez, M., 2000, The geologic evolution of southern Sierra de Guanajuato, México: A documented example of the transition from the Sierra Madre Occidental to the Mexican Volcanic Belt: International Geology Review, v. 42, p. 131-151, doi:10.1080/00206810009465073.

Chauve, P., Fourçade, E., and Carrillo-Martínez, M., 1985, Le raports structuraux entre les domaines cordillérain et mésogéen dans la partie centrale du Mexique: Les Comptes Rendus de l'Académie de Paris, v. 301, p. 335-340.

Coenraads, R.R., and Koivula, J.I., 2007, Earth's dynamic forces: Helenora Heights, Australia, Geologica, Millennium House Pty. Ltd., 576 p.

Cotten, J., Le Dez, A., Bau, M., Caroff, M., Maury, R., Dulski, P., Fourçade, S., Bohn, M., and Brousse, R., 1995, Origin of anomalous rare-earth element and yttrium enrichments in subaerially exposed basalts: Evidence from French Polynesia: Chemical Geology, v. 119, p. 115-138, doi:10.1016/0009-2541(94)00102-E.

Demant, A., 1978, Interpretación geodinámica del volcánismo del Eje Neovolcánico Transmexicano: Revista Mexicana de Ciencias Geológicas, v. 5, p. 217-222.

Duffield, W.A., Richter, D.H., and Priest, S.S., 1995, Physical volcanology of silicic lava domes as exemplified by the Taylor Creek Rhyolite, Catron and Sierra Counties, New Mexico: U.S. Geological Survey Map I-2399, scale 1:50,000.

Echegoyén, J., Romero-Martínez, S., and Velázquez-Silva, S., 1970, Geología y yacimientos minerales de la parte

central del distrito minero de Guanajuato: Boletín Consejo de Recursos Minera no Renovables, v. 75, 35 p.

Edwards, D.J., 1955, Studies of some early Tertiary red conglomerates of central Mexico: U.S. Geological Survey Professional Paper 264-H, p. 153-185.

INAH, 2009, Instituto Nacional de Antropología e Historia, México: Official web page www.inah.gob.mx, note retrieved on September 9, 2009.

INEGI, 2002, Carta topográfica escala 1:50,000, CD con datos vectoriales, Datum: ITRF92, Elipsoide: GRS80, Cartas F14C57 y F14C67.

Labarthe-Hernández, G., Tristán-González, M., and Aranda-Gómez, J.J., 1982, Revisión estratigráfica del Cenozoico de la parte central del Estado de San Luis Potosí: Universidad Autónoma de San Luis Potosí, Instituto de Geología y Metalurgia, Folletó Técnico 85, 208 p.

Le Bas, M.J., Le Maitre, R.W., Streckeisen, A., and Zanettin, B., 1986, A chemical classification of volcanic rocks based on the total alkali-silica diagram: Journal of Petrology, v. 27, p. 745-750, doi:10.1093/petrology/27.3.745.

Lindsley, D.H., and Frost, B.R., 1992, Equilibria among Fe-Ti oxides, pyroxenes, olivine, and quartz: Part I. Theory: The American Mineralogist, v. 77, p. 987-1003.

Lozano-Santacruz, R., and Bernal, J.P., 2005, Characterization of a new set of eight geochemical reference materials for XRF major and trace element analysis: Revista Mexicana de Ciencias Geológicas, v. 22, p. 329-344.

McDowell, F.W., and Clabaugh, S.E., 1979, Ignimbrites of the Sierra Madre Occidental and their relation to the tectonic history of western Mexico: Geological Society of America Special Paper 118, p. 113-124.

Municipio Ezequiel Montes, 2009, Official web page of Municipio de Ezequiel Montes, Querétaro: http://www.qqueretaro-mexico.com.mx/montes/, retrieved on September 9, 2009.

Nelson, S.A., and Carmichael, I.S.E., 1984, Pleistocene to Recent alkali volcanism in the region of Sanganguey volcano, Nayarit, Mexico: Contributions to Mineralogy and Petrology, v. 85, p. 321-335, doi:10.1007/BF01150290.

Padilla y Sánchez, R.J., 1985, Las estructuras de la Curvatura de Monterrey, Estados de Coahuila, Nuevo León, Zacatecas y San Luis Potosí: Revista Mexicana de Ciencias Geológicas, v. 6, p. 1-20.

Peccerillo, A., and Taylor, S.R., 1976, Geochemistry of Eocene calc-alkaline volcanic rocks from the Camentou area, northern Turkey: Contributions to Mineralogy and Petrology, v. 58, p. 63-81, doi:10.1007/BF00384745.

Renne, P.R., Swisher, C.C., Deino, A.L., Karner, D.B., Owens, T.L., and DePaolo, D.J., 1998, Intercalibration of standards, absolute ages and uncertainties in $^{40}$Ar/$^{39}$Ar dating: Chemical Geology, v. 145, p. 117-152, doi:10.1016/S0009-2541(97)00159-9.

Segerstrom, K., 1961, Geology of the Bernal-Jalpan Area Estado de Querétaro Mexico: U.S. Geological Survey Bulletin 1104-B, p. 19-82.

Segerstrom, K., 1962, Geology of south-central Hidalgo and northeastern Mexico, Mexico: U.S. Geological Survey Bulletin 1104-C, p. 87-162.

Siebe, C., Macías, J.L., and Aguirre-Díaz, G.J., eds., 2006, Neogene-Quaternary Continental Margin Volcanism: A Perspective from Mexico: Boulder Colorado, Geological Society of America Special Paper 402, 329 p.

Steiger, R.H., and Jäger, C., 1977, Subcommission on Geochronology: Convention on the use of decay constants in Geo and Cosmochronology: Earth and Planetary Science Letters, v. 36, p. 359-362.

Suter, M., 1987, Structural traverse across the Sierra Madre Oriental fold-thrust belt in east-central Mexico: Geological Society of America Bulletin, v. 98, p. 249-264, doi:10.1130/0016-7606(1987)98:249:STATSM-2.0:CO2.

Tristán-González, M., Aguirre-Díaz, G.J., Labarthe-Hernández, G., Torres-Hernández, J.R., and Bellon, H., 2009, Post-Laramide and pre-Basin and Range deformation and implications for Paleogene (55-25 Ma) volcanism in central Mexico: A geological basis for a volcano-technical stress model: Tectonophysics, v. 471, p. 136-152.

UNESCO, 2001, Sierra Gorda (Mexico) and Waterberg (South Africa) join world network of UNESCO biosphere reserves, The MAB United Nations Educational, Scientific and Cultural Organization, UNESCO MAT Program, UNESCO Press, March 2001 issue: http://www.unesco.org/bpi/eng/unescopress/2001/01-46e.shtml (March, 2001).

UNESCO, 2009, Places of memory and living traditions of the Otomí-Chichimecas people of Tolimán: The Peña de Bernal, guardian of a sacred territory, UNESCO Culture Sector, Intangible Heritage: http://www.unesco.org/culture/ich/RL/00174 (September 2009).

Valdez-Moreno, G., Aguirre-Díaz, G.J., and López-Martínez, M., 1998, El Volcán La Joya, Edos. de Querétaro y Guanajuato: Un Estratovolcán antiguo del Cinturón Volcánico Mexicano: Revista Mexicana de Ciencias Geológicas, v. 15, p. 181-197.

Verma, S.P., and Carrasco-Núñez, G., 2003, Reappraisal of the geology and geochemistry of volcan Zamorano, Central Mexico: Implications for discriminating the Sierra Madre Occidental and Mexican Volcanic Belt provinces: International Geology Review, v. 45, p. 724-752, doi:10.2747/0020-6814.45.8.724.

Verma, S.P., Sotelo-Rodríguez, Z.T., and Torres-Alvarado, I.S., 2002, SINCLAS: Standard Igneous Norm and Volcanic Rock Classification System: Computers and Geosciences, v. 28, p. 711-715.

Wallace, P.J., and Carmichael, I.S.E., 1994, Petrology of Volcan Tequila, Jalisco, Mexico: Disequilibrium phenocryst assemblages and evolution of the subvolcanic magma system: Contributions to Mineralogy and Petrology, v. 117, p. 345-361, doi:10.1007/BF00307270.

White, D.E., 1948, Antimony deposits of the Soyatal District, State of Querétaro, Mexico: U.S. Geological Survey Bulletin 960-B, p. 35-88.

Wilcox, R.E., 1979, The liquid line of descent and variation diagrams, in Yoder, H.S., Jr., ed., The Evolution of the Igneous Rocks: Princeton University Press, p. 206-232.

Wilson, M., 1989, Igneous Petrogenesis: A Global Tectonic Approach: London, Unwin Hyman, 466 p.

Wilson, B.W., Hernández-M., J.P., and Meave, T.E., 1955, Un banco calizó del Cretácico en la parte oriental del estado de Querétaro: Boletín de la Sociedad Geológica Mexicana, v. 18, p. 1-10.

York, D., Evensen, N.M., López-Martínez, M., and De Basabe-Delgado, J., 2004, Unified equations for the slope, intercept, and standard errors of the best straight line: American Journal of Physics, v. 72, no. 3, p. 367-375.