Brian R. Jicha

Magmatic history and evolution of the Central American Land Bridge in Panama since Cretaceous times

Chemical compositions for 310 igneous rocks from the Cordillera de Panama and the Sona and Azuero peninsulas were supplemented by 40 Ar/ 39 Ar dating and Sr-, Nd-, Pb-, and O-isotope analysis to determine the magmatic evolution and oceanic plate interactions over the past 100 Ma in western Panama. An initial phase of intraplate magmatism, having geochemical characteristics of the Galapagos hotspot, formed the oceanic basement of the Caribbean large igneous province from 139 to 69 Ma. Younger accreted terranes with enriched trace element patterns (accreted ocean island basalt [OIB]) were amalgamated between 70 and 20 Ma. A second magmatic phase in the Azuero and Sona peninsulas has trace element patterns (Sona-Azuero arc) suggesting the initiation of subduction at 71–69 Ma. Arc magmatism continued in the Chagres basin region (Chagres-Bayano arc) from 68 to 40 Ma. A third phase formed discrete volcanic centers across the Cordillera de Panama (Cordilleran arc) from 19 to 5 Ma. The youngest phase consists of isolated volcanic centers of adakitic composition (Adakite suite) in the Cordillera de Panama that developed over the past 2 million years. Initiation of arc magmatism at 71 Ma coincides with the cessation of Galapagos plateau formation, suggesting a causal link. The transition from intraplate to arc magmatism occurred relatively quickly and introduced a new enriched mantle source. The arc magmatism involved progressive transition to more homogeneous intermediate mantle wedge compositions through mixing and homogenization of subarc magma sources through time and/or the replacement of the mantle wedge by a homogeneous, relatively undeleted asthenospheric mantle. Adakite volcanism started after a magmatic gap, enabled by the formation of a slab window.

Intercalibration of radioisotopic and astrochronologic time scales for the Cenomanian-Turonian boundary interval, Western Interior Basin, USA

We develop an intercalibrated astrochronologic and radioisotopic time scale for the Cenomanian-Turonian boundary (CTB) interval near the Global Stratotype Section and Point in Colorado, USA, where orbitally influenced rhythmic strata host bentonites that contain sanidine and zircon suitable for 40Ar/39Ar and U-Pb dating. Paired 40Ar/39Ar and U-Pb ages are determined from four bentonites that span the Vascoceras diartianum to Pseudaspidoceras flexuosum ammonite biozones, utilizing both newly collected material and legacy sanidine samples of J. Obradovich. Comparison of the 40Ar/39Ar and U-Pb results underscores the strengths and limitations of each system, and supports an astronomically calibrated Fish Canyon sanidine standard age of 28.201 Ma. The radioisotopic data and published astrochronology are employed to develop a new CTB time scale, using two statistical approaches: (1) a simple integration that yields a CTB age of 93.89 ± 0.14 Ma (2σ; total radioisotopic uncertainty), and (2) a Bayesian intercalibration that explicitly accounts for orbital time scale uncertainty, and yields a CTB age of 93.90 ± 0.15 Ma (95% credible interval; total radioisotopic and orbital time scale uncertainty). Both approaches firmly anchor the floating orbital time scale, and the Bayesian technique yields astronomically recalibrated radioisotopic ages for individual bentonites, with analytical uncertainties at the permil level of resolution, and total uncertainties below 2‰. Using our new results, the duration between the Cenomanian-Turonian and the Cretaceous-Paleogene boundaries is 27.94 ± 0.16 Ma, with an uncertainty of less than one-half of a long eccentricity cycle.

Eruptive history, geochronology, and magmatic evolution of the Puyehue-Cordón Caulle volcanic complex, Chile

Forty-three 40 Ar/ 39 Ar age determinations of lava flows, domes, ignimbrites, and dikes, plus 14 C dates from seven distal tephra layers, combined with stratigraphy, geochemistry, and Sr and Th isotope data, establish an eruptive chronology for the Puyehue-Cordon Caulle volcanic complex at 40.5° S in the Andean southern volcanic zone (SVZ). The complex preserves ~131 km 3 of lava and tephra that erupted from numerous vents widely separated in time and space. Approximately 80% of the total volume consists of basaltic to andesitic lava that formed two broad shield volcanoes between 314 and 70 ka. The modern Puyehue stratovolcano was built on the southerly shield during the past 69 k.y. following a hiatus of 25 k.y. Puyehue has erupted ~15 km 3 of basaltic to rhyolitic magma that spans the entire compositional range found in the southern SVZ and evolved via at least six phases including: (1) basaltic andesitic to dacitic lavas between 69 and 32 ka, (2) a shift to bimodal magma compositions that is first expressed by a rhyodacite mingled with inclusions of MgO-rich basaltic andesite at 34 ka, (3) dacitic to rhyolitic flows and domes from 19 to 12 ka, (4) basaltic to basaltic andesitic flows between 15 and 12 ka, (5) subsequent rhyolitic dome growth in several effusive and explosive stages between 7 and 5 ka, followed by (6) a powerful series of phreatomagmatic and sub-Plinian eruptions at ca. 1.1 ka that obliterated the preceding rhyolite domes and formed the present 2.5-km-diameter, 280-m-deep summit crater. Along the Cordon Caulle fissure zone, ~5 km 3 of rhyodacitic to rhyolitic lavas, domes, and cones have formed during the past ~16.5 k.y., including explosive and effusive eruptions in 1921–1922 and 1960. Eruptive rates were nonuniform over time, with background growth at 0.04 km 3 /k.y. or less, punctuated by spurts at up to 0.90 km 3 /k.y. The time-averaged rate, 0.42 km 3 /k.y., is nearly double that at the Tatara-San Pedro complex 500 km to the north during the past 300 k.y. These findings indicate that within a single arc the magmatic and eruptive fluxes at individual frontal volcanoes can be highly variable. The last three stratocone-building events on Puyehue began during periods of deglaciation, suggesting a relationship between unloading of ice and ease of magma ascent. Puyehue basalt exhibits subtle changes in 238 U- 230 Th, 87 Sr/ 86 Sr, and trace element composition over time that signal shifts in the composition and degree of melting of the mantle wedge, or the extent to which basalt was modified by assimilation of heterogeneous crustal melts. The complex has become exceptionally bimodal and more explosive over time with recent rhyolites evolving by extreme crystal fractionation of mafic magma and lesser volumes of andesite and dacite created via mixing of rhyolite and basalt. Despite the high flux of basalt during the past 300 k.y., no large silicic magma reservoir formed in the upper crust. Instead, 238 U- 230 Th data favor rapid ascent of several small bodies of basaltic and silicic magma from the lower crust, promoted perhaps by conduits that reflect strike-slip faulting beneath the complex.

Revised age of Aleutian Island Arc formation implies high rate of magma production

Radioisotopic dating of subaerial and submarine volcanic and plutonic rocks from the Aleutian Island Arc provides insight into the timing of arc formation in the middle Eocene. Twenty-eight 40 Ar/ 39 Ar ages constrain the duration of arc magmatism to the last 46 m.y. Basaltic lavas from the Finger Bay volcanics, the oldest exposed rocks in the arc, gave an isochron age of 37.4 ± 0.6 Ma, which is 12-17 m.y. younger than a widely cited age of 55-50 Ma. Three main pulses of arc-wide magmatism occurred at 38-29, 16-11, and 6-0 Ma, which coincide with periods of intense magmatism in other western Pacific island arcs. Using the geochronology and volumetric estimates of crust generated and eroded over the last 46 m.y., we calculate a time-averaged magma production rate for the entire arc that exceeds previous estimates by almost an order of magnitude.

Structural and temporal requirements for geomagnetic field reversal deduced from lava flows

Reversals of the Earths magnetic field reflect changes in the geodynamo—flow within the outer core—that generates the field. Constraining core processes or mantle properties that induce or modulate reversals requires knowing the timing and morphology of field changes that precede and accompany these reversals. But the short duration of transitional field states and fragmentary nature of even the best palaeomagnetic records make it difficult to provide a timeline for the reversal process. 40Ar/39Ar dating of lavas on Tahiti, long thought to record the primary part of the most recent ‘Matuyama–Brunhes’ reversal, gives an age of 795 ± 7 kyr, indistinguishable from that of lavas in Chile and La Palma that record a transition in the Earths magnetic field, but older than the accepted age for the reversal. Only the ‘transitional’ lavas on Maui and one from La Palma (dated at 776 ± 2 kyr), agree with the astronomical age for the reversal. Here we propose that the older lavas record the onset of a geodynamo process, which only on occasion would result in polarity change. This initial instability, associated with the first of two decreases in field intensity, began ∼18 kyr before the actual polarity switch. These data support the claim that complete reversals require a significant period for magnetic flux to escape from the solid inner core and sufficiently weaken its stabilizing effect.

40Ar/39Ar chronostratigraphy of Altiplano-Puna volcanic complex ignimbrites reveals the development of a major magmatic province

The Lipez region of southwest Bolivia is the locus of a major Neogene ignimbrite flare-up, and yet it is the least studied portion of the Altiplano-Puna volcanic complex of the Central Andes. Recent mapping and laser-fusion 40 Ar/ 39 Ar dating of sanidine and biotite from 56 locations, coupled with paleomagnetic data, refine the timing and volumes of ignimbrite emplacement in Bolivia and northern Chile to reveal that monotonous intermediate volcanism was prodigious and episodic throughout the complex. The new results unravel the eruptive history of the Pastos Grandes and Guacha calderas, two large multicyclic caldera complexes located in Bolivia. These two calderas, together with the Vilama and La Pacana caldera complexes and smaller ignimbrite shields, were the dominant sources of the ignimbrite-producing eruptions during the ∼10 m.y. history of the Altiplano-Puna volcanic complex. The oldest ignimbrites erupted between 11 and 10 Ma represent relatively small volumes (approximately hundreds of km 3 ) of magma from sources distributed throughout the volcanic complex. The first major pulse was manifest at 8.41 Ma and 8.33 Ma as the Vilama and Sifon ignimbrites, respectively. During pulse 1, at least 2400 km 3 of dacitic magma was erupted over 0.08 m.y. Pulse 2 involved near-coincident eruptions from three of the major calderas resulting in the 5.60 Ma Pujsa, 5.65 Ma Guacha, and 5.45 Ma Chuhuilla ignimbrites, for a total minimum volume of 3000 km 3 of magma. Pulse 3, the largest, produced at least 3100 km 3 of magma during a 0.1 m.y. period centered at 4 Ma, with the eruption of the 4.09 Ma Puripicar, 4.00 Ma Chaxas, and 3.96 Ma Atana ignimbrites. This third pulse was followed by two more volcanic explosivity index (VEI) 8 eruptions, producing the 3.49 Ma Tara (800 km 3 DRE) and 2.89 Ma Pastos Grandes (1500 km 3 DRE) ignimbrites. In addition to these large caldera-related eruptions, new age determinations refine the timing of two little-known ignimbrite shields, the 5.23 Ma Alota and 1.98 Ma Laguna Colorada centers. Moreover, 40 Ar/ 39 Ar age determinations of 13 ignimbrites from northern Chile previously dated by the K-Ar method improve the overall temporal resolution of Altiplano-Puna volcanic complex development. Together with the updated volume estimates, the new age determinations demonstrate a distinct onset of Altiplano-Puna volcanic complex ignimbrite volcanism with modest output rates, an episodic middle phase with the highest eruption rates, followed by a decline in volcanic output. The cyclic nature of individual caldera complexes and the spatiotemporal pattern of the volcanic field as a whole are consistent with both incremental construction of plutons as well as a composite Cordilleran batholith.

Dynamics of a large, restless, rhyolitic magma system at Laguna del Maule, southern Andes, Chile

Explosive eruptions of large-volume rhyolitic magma systems are common in the geologic record and pose a major potential threat to society. Unlike other natural hazards, such as earthquakes and tsunamis, a large rhyolitic volcano may provide warning signs long before a caldera-forming eruption occurs. Yet, these signs—and what they imply about magma-crust dynamics—are not well known. This is because we have learned how these systems form, grow, and erupt mainly from the study of ash flow tuffs deposited tens to hundreds of thousands of years ago or more, or from the geophysical imaging of the unerupted portions of the reservoirs beneath the associated calderas. The Laguna del Maule Volcanic Field, Chile, includes an unusually large and recent concentration of silicic eruptions. Since 2007, the crust there has been inflating at an astonishing rate of at least 25 cm/yr. This unique opportunity to investigate the dynamics of a large rhyolitic system while magma migration, reservoir growth, and crustal deformation are actively under way is stimulating a new international collaboration. Findings thus far lead to the hypothesis that the silicic vents have tapped an extensive layer of crystal-poor, rhyolitic melt that began to form atop a magmatic mush zone that was established by ca. 20 ka with a renewed phase of rhyolite eruptions during the Holocene. Modeling of surface deformation, magnetotelluric data, and gravity changes suggest that magma is currently intruding at a depth of ~5 km. The next phase of this investigation seeks to enlarge the sets of geophysical and geochemical data and to use these observations in numerical models of system dynamics. INTRODUCTION Caldera-scale rhyolitic volcanoes can rapidly deposit hundreds of cubic kilometers of ash over several million square kilometers, threatening people and agriculture at the scale of an entire continent (Sparks et al., 2005; Lowenstern et al., 2006; Self, 2006). Sooner or later, Earth will experience another eruption of this magnitude (Lowenstern et al., 2006; Self and Blake, 2008); consequently, there is a need to gather comprehensive information and create multi-scale models that realistically capture the dynamics leading to these destructive events. Most of our current understanding of this type of volcanic system has been gleaned from the study of eruptive products long after the catastrophic eruption, including voluminous ash flow deposits, such as the Bishop, Bandelier, Huckleberry Ridge, and Oruanui Tuffs (Lowenstern et al., 2006; Hildreth and Wilson, 2007; Bachmann and Bergantz, 2008; Wilson, 2008). The most recent rhyolitic “super-eruption” produced the Oruanui Tuff 26,500 years ago in New Zealand. Even in this relatively recent case, the geologic evidence has been partly obliterated by caldera-collapse, erosion, and burial (Wilson et al., 2005). Moreover, probing the present-day structures beneath a number of calderas using seismic tomography (e.g., Romero et al., 1993; Steck et al., 1998; Farrell et al., 2014) or other geophysical measures (e.g., Lowenstern et al., 2006; Battaglia et al., 2003; Tizzani et al., 2009) has not detected eruptible domains of crystal-poor melt in the shallow crust, nor has it captured the dynamics that preceded these large eruptions. This paper focuses on the Laguna del Maule Volcanic Field, Chile, a large, potentially hazardous, rhyolitic magmatic system, where an alarming rate of surface uplift for the past seven years and concentrated swarms of shallow earthquakes prompted Observatorio Volcanologico de los Andes del Sur (OVDAS) to declare in March 2013 a yellow alert, signaling a potential eruption within months or years. Straddling the Andean range crest at 36° S (Fig. 1A), this volcanic field features: (1) 13 km of rhyolite that erupted both explosively and effusively during the past 20 k.y.; (2) a zone of low electrical resistivity in the shallow crust below the deforming area; (3) widespread elevated CO 2 concentrations; and (4) a negative (~10 mGal) Bouguer anomaly and preliminary evidence for a positive dynamic gravity signal indicating mass addition. The underlying magma system has been sampled by eruptions numerous times since its apparent inception in the late Pleistocene, including a dozen crystal-poor, glassy rhyolitic lavas during the Holocene. Linking the assembly and evolution of this

Integrating 40Ar/39Ar, U-Pb, and astronomical clocks in the Cretaceous Niobrara Formation, Western Interior Basin, USA

This study revises and improves the chronostratigraphic framework for late Turonian through early Campanian time based on work in the Western Interior U.S. and introduces new methods to better quantify uncertainties associated with the development of such time scales. Building on the unique attributes of the Western Interior Basin, which contains abundant volcanic ash beds and rhythmic strata interpreted to record orbital cycles, we integrate new radioisotopic data of improved accuracy with a recently published astrochronologic framework for the Niobrara Formation. New 40Ar/39Ar laser fusion ages corresponding to eight different ammonite biozones are determined by analysis of legacy samples, as well as newly collected material. These results are complemented by new U-Pb (zircon) chemical abrasion–isotope dilution–thermal ionization mass spectrometry ages from four biozones in the study interval. When combined with published radioisotopic data from the Cenomanian-Turonian boundary, paired 206Pb/238U and 40Ar/39Ar ages spanning Cenomanian to Campanian time support an astronomically calibrated Fish Canyon sanidine standard age of 28.201 Ma. Stage boundary ages are estimated via integration of new radioisotopic data with the floating astrochronology for the Niobrara Formation. The ages are determined by anchoring the long eccentricity bandpass from spectral analysis of the Niobrara Formation to radioisotopic ages with the lowest uncertainty proximal to the boundary, and adding or subtracting time by parsing the 405 k.y. cycles. The new stage boundary age determinations are: 89.75 ± 0.38 Ma for the Turonian-Coniacian, 86.49 ± 0.44 Ma for the Coniacian-Santonian, and 84.19 ± 0.38 Ma for the Santonian-Campanian boundary. The 2σ uncertainties for these estimates include systematic contributions from the radioisotopic measurements, astrochronologic methods, and geologic uncertainties (related to stratigraphic correlation and the presence of hiatuses). The latter geologic uncertainties have not been directly addressed in prior time scale studies and their determination was made possible by critical biostratigraphic observations. Each methodological approach employed in this study—new radioisotopic analysis, stratigraphic correlation, astrochronology, and ammonite and inoceramid biostratigraphy—was critical for achieving the final result.

Volcanic history and magmatic evolution of Seguam Island, Aleutian Island arc, Alaska

New 40 Ar/ 39 Ar dating coupled with detailed field mapping, stratigraphy, and chemical analyses have established an eruptive chronology that reveals and constrains the compositional and volumetric evolution of Seguam Island in the Aleutian Island arc, Alaska. Sixty new 40 Ar/ 39 Ar ages from lavas, domes, and pyroclastic deposits were obtained using furnace incremental-heating techniques on replicate samples of whole-rock and ground-mass separates, and they constrain the duration of Pleistocene to Holocene subaerial volcanism to 318 k.y. The 40 Ar/ 39 Ar plateau ages indicate that over 85% of the complex, ∼68 km 3 of material, was erupted almost continuously between 318 ka and 9 ka. At ca. 9 ka, a stratocone on the eastern half of the island partially collapsed producing a 4-km-wide crater. Rhyolitic dome-forming eruptions followed from vents in the newly created crater, and were likely contemporaneous with 8.0 km 3 of basaltic and basaltic andesitic effusions from Pyre Peak, and a 1.4 km 3 basaltic eruption from a monogenetic cone on the far eastern end of the island. Geochemical changes over the last 318 k.y. are subtle. Most notably, the earliest eruptions from 318 to 142 ka produced no andesite, and basalt from this period has larger ranges in Zr/Rb and La/Yb than younger basalts. Small volumes of dacitic to rhyolitic magma were produced from basalt by a monotonic crystal-liquid fractionation process that varied only slightly in successive eruptive phases over 318 k.y. We identified minor geochemical changes in magma composition during each of the three main stages of volcanism, but overall the monotonic variations in major- and trace-element compositions of basaltic andesitic to rhyolitic lavas are consistent with an origin via closed-system fractional crystallization of basalt. Using the 40 Ar/ 39 Ar geochronology and estimates of individual flow volumes, we calculated a time-averaged eruptive rate at Seguam that is similar to growth rates of other well-dated arc volcanoes in the Cascades and Chilean Southern volcanic zone but less than that of Mount Katmai and Mount Mageik, which are located on the Alaska Peninsula. The eruptive flux at Seguam has been highly variable, fluctuating more than an order of magnitude, from 0.07 km 3 /k.y. during the early history of bimodal volcanism to 1.18 km 3 /k.y. over the past 9 k.y.

New age constraints for the Salamanca Formation and lower Río Chico Group in the western San Jorge Basin, Patagonia, Argentina: Implications for cretaceous-paleogene extinction recovery and land mammal age correlations

The Salamanca Formation of the San Jorge Basin (Patagonia, Argentina) preserves critical records of Southern Hemisphere Paleocene biotas, but its age remains poorly resolved, with estimates ranging from Late Cretaceous to middle Paleocene. We report a multi-disciplinary geochronologic study of the Salamanca Formation and overlying Rio Chico Group in the western part of the basin. New constraints include (1) an 40Ar/39Ar age determination of 67.31 ± 0.55 Ma from a basalt flow underlying the Salamanca Formation, (2) micropaleontological results indicating an early Danian age for the base of the Salamanca Formation, (3) laser ablation HR-MC-ICP-MS (high resolution-multi collector-inductively coupled plasma-mass spectrometry) U-Pb ages and a high-resolution TIMS (thermal ionization mass spectrometry) age of 61.984 ± 0.041(0.074)[0.100] Ma for zircons from volcanic ash beds in the Penas Coloradas Formation (Rio Chico Group), and (4) paleomagnetic results indicating that the Salamanca Formation in this area is entirely of normal polarity, with reversals occurring in the Rio Chico Group. Placing these new age constraints in the context of a sequence stratigraphic model for the basin, we correlate the Salamanca Formation in the study area to Chrons C29n and C28n, with the Banco Negro Inferior (BNI), a mature widespread fossiliferous paleosol unit at the top of the Salamanca Formation, corresponding to the top of Chron C28n. The diverse paleobotanical assemblages from this area are here assigned to C28n (64.67–63.49 Ma), ∼2–3 million years older than previously thought, adding to growing evidence for rapid Southern Hemisphere floral recovery after the Cretaceous-Paleogene extinction. Important Peligran and “Carodnia” zone vertebrate fossil assemblages from coastal BNI and Penas Coloradas exposures are likely older than previously thought and correlate to the early Torrejonian and early Tiffanian North American Land Mammal Ages, respectively.