Search USGSSearch

USGS · 70033010

The role of ridge subduction in determining the geochemistry and Nd–Sr–Pb isotopic evolution of the Kodiak batholith in southern Alaska

Abstract

The Paleocene Kodiak batholith, part of the Sanak&ndash;Baranof belt of Tertiary near-trench intrusive rocks, forms an elongate body (~ 150 km long) that transects Kodiak Island from SW to NE. The batholith consists of three zones (Southern, Central, and Northern) of kyanite-, muscovite-, and garnet-bearing biotite tonalite and granodiorite and less abundant granite that intruded an accretionary prism (Kodiak Formation, and Ghost Rocks Formation). Small and likely coeval bodies (Northern, Western, and Eastern satellite groups) of quartz gabbro, diorite, tonalite, granodiorite, and leucogranite flank the batholith. The batholith is calc-alkalic, has an aluminum saturation index of > 1.1, FeO t /(FeO t + MgO) ~ 0.65 (at SiO 2 = 65 wt.%), and increases in SiO 2 (~ 61 wt.%&ndash;73 wt.%) and decreases in TiO 2 (~ 0.9 wt.%&ndash;0.3 wt.%) from SW to NE. As a group, the granitic rocks have light REE-enriched chondrite-normalized patterns with small or no negative Eu anomalies, primitive mantle-normalized negative anomalies for Nb and Ti, and positive anomalies for Pb. Small to large negative anomalies for Th are also distinctive. The quartz gabbros and diorites are generally characterized by generally flat to light REE chondrite-normalized patterns (no Eu anomalies), and mantle-normalized negative anomalies for Nb, Ti, and P. Pb isotopic compositions ( 206 Pb/ 204 Pb = 18.850&ndash;18.960; 207 Pb/ 204 Pb = 15.575&ndash;15.694; 208 Pb/ 204 Pb = 38.350&ndash;39.039) are intermediate between depleted mantle and average continental crust. The Southern zone and a portion of the Central zone are characterized by negative &epsilon; Nd values of &minus; 3.7 to &minus; 0.3 and T DM ages ranging from ~ 838 Ma to 1011 Ma. Other granitic rocks from the Central and Northern zones have higher &epsilon; Nd values of &minus; 0.4 to + 4.7 and younger T DM ages of ~ 450 to 797 Ma. Granitic and mafic plutons from the Eastern satellites show a wide range of &epsilon; Nd values of &minus; 2.7 to + 6.4, and T DM ages from 204 Ma to 2124 Ma. 87 Sr/ 86 Sr values of the Southern and Central zones overlap and tend to be slightly more radiogenic ( 87 Sr/ 86 Sr > 0.70426) than the Northern zone ( 87 Sr/ 86 Sr < 0.70472). 206 Pb/ 204 Pb values increase slightly from the Southern and Central zones toward the Northern zone. There is no clear correlation of the major or trace elements with &epsilon; Nd , Pb or Sr isotopic values. Kodiak Formation and the Ghost Rocks Formation overlap the isotopic compositions (e.g., 206 Pb/ 204 Pb = 18.978 to 19.165, 87 Sr/ 86 Sr of 0.705715 to 0.707118, and &epsilon; Nd of &minus; 6.7 to &minus; 1.5 at 59 Ma) and T DM values (959 to 1489 Ma) of the batholith. Production of large volumes of granitic rocks in the Sanak&ndash;Baranof belt, and particularly on Kodiak Island, reflects a sequence of processes that includes underplating of mantle-derived mafic (possibly from the mantle wedge) and intermediate rocks under the accretionary flysch, interlayering of mantle-derived and flyschoid rocks, and partial melting of the mixed lithologic assemblages. Limited degrees of fractional crystallization or assimilation and fractional crystallization influenced compositions of the granitic rocks. The contribution of mantle-derived rocks that resided in the accretionary prism for only a short period of time prior to partial melting likely exceeds 40% (up to 80%). The balance (60 to 20%) is from a recently recycled crustal component represented by the Kodiak Formation. This type of progressive intracrustal melting from mixed sources controlled the geochemical character of the batholith and is most consistent with the hypothesis that the granitic rocks are associated with a slab-window produced by collision of a spreading oceanic center and a subduction zone and migration beneath the accretionary prism.

Explore related subjects

90° N90° S · 180° W ← longitude → 180° E
Source-reported bounding extent: 55.85064987433714° to 62.01121819833755° latitude; -155.36865234375° to -134.38476562499997° longitude. This indicates report coverage, not an exact sampling location. View area on OpenStreetMap.

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Robert A. Ayuso, Peter J. Haeussler, Dwight Bradley, David W. Farris, Nora K. Foley, Gregory A. Wandless. 2009. The role of ridge subduction in determining the geochemistry and Nd–Sr–Pb isotopic evolution of the Kodiak batholith in southern Alaska. https://doi.org/10.1016/j.tecto.2008.09.029

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

GST-1: A high-resolution global sediment thickness model

Global Sediment Thickness 1 (GST-1) is a high-resolution sedimentary thickness model calculated on a 0.125° x 0.125° grid. It modifies the sediment thickness of the 1° x 1° Earth Crustal Model 1 (ECM1) by means of 3D inversions of free air gravity anomalies. GST-1 is calculated by performing structural inversions on high-density contrasts across two crustal boundaries: the sediment – basement interface and the crystalline crust – upper mantle interface. The inversions are calculated in each of ten overlapping 3D models that span the globe. These ten models are merged to obtain the GST-1 global model, providing an eight-fold increase in lateral spatial resolution in comparison with ECM1 and CRUST 1.0. Our sediment thickness model exploits the nearly continuous sampling of gravity data when compared to the irregular, sparse sampling of seismic refraction data. Sediment thickness values in GST-1 are in excellent agreement with independently derived cross sections from well-studied sedimentary basins, and within expected resolution limits of seismic refraction data. GST-1 offers a robust, high resolution global model of sedimentary thickness to support studies of sedimentary basins.

Tectonophysics

Fairweather transform boundary Oligocene to present orogenesis: Fairweather Range vertical extrusion and rotation of the Yakutat microplate at ca. 3 Ma

Oblique-slip along transform fault boundaries is often partitioned between a strike-slip system and thrust faults that accommodate contraction. However, topography along the Yakutat-North American transform (Fairweather fault), is asymmetric with low-terrain above active thrusts on the western, Yakutat side of the transform and high topography on the continental side with peaks >4500 m (Mount Fairweather: 4671 m) to the west of the Border Ranges fault, limited recorded earthquakes >M4, and no apparent reverse faults to generate the highest terrain. In this study we compile, for the first time, published U-Pb zircon, 40 Ar/ 39 Ar and K-Ar (hornblende, muscovite, and biotite) and U-Th/He and fission-track (zircon and apatite) bedrock ages (109) from 75 samples to investigate the exhumation history of the Fairweather Range region, complemented by a published detrital sample (ZFT and AFT) and 13 new 40 Ar/ 39 Ar (hornblende, biotite, and K-feldspar) ages on 9 bedrock samples from both sides of the Fairweather fault. Additionally, we examined published seismicity and geodetic data of the Fairweather region and assessed if plate paleo-vectors correlate with the cooling history of the Fairweather Range. Cooling age, seismic, and block-motion patterns indicate the Fairweather Range has been vertically extruded between the Fairweather and the Border Ranges faults as a coherent block since ca. 25 Ma. The pre-6 Ma Pacific plate motion (N30°W) aligns with the N33°W strike of the Fairweather Fault whereas a hypothetical pre-6 Ma Yakutat microplate paleo-vector of (N39°W) does not: indicating a post-6 Ma timing for Yakutat microplate counter-clockwise rotation (9°). We infer that rotation and impingement of the Yakutat microplate along the Fairweather fault at ca. 3 Ma led to the development of the Fairweather restraining bend and increased cooling rates. The resultant thickened Fairweather welt and the ∼30 km thick southeast end of the Yakutat microplate compounded double-indenter tectonics into Alaska's southeast convergent corner

Tectonophysics

A scaling relationship for the width of secondary deformation around strike-slip faults

Simple mechanical arguments suggest that slip along interlocked, rough faults, damages surrounding rocks. The same arguments require that the scale of secondary damage is proportional to the size of geometric irregularities along the main fault. This relationship could apply at all scales, but has, so far, been difficult to observe at the 10s to 100 s of km scales of large, natural faults, often because large-scale deformation is distributed across wide, complex plate-boundary fault systems, like the San Andreas Fault. The geometry and geology of another large-scale plate-boundary strike slip fault—the Queen Charlotte Fault (QCF)—is, in contrast, especially simple. Here, we show that observations of secondary deformation are well-aligned with predictions of stress variations caused by geometric irregularities along the QCF, suggesting a geometric relationship between primary fault geometry and secondary deformation. The analytic stress solution reveals that the highest stresses and highest likelihood of failure are confined to a zone of influence (ZOI) with a width quantified by ZOI = &#x3BB; / 2 &#x3C0; "> ZOI=λ/2π , where λ is the wavelength of geometric variations along the main fault. This simple model is consistent with ∼100-km-scale observations along the QCF and can theoretically be used to predict the width of secondary deformation at all scales.

Queen Charlotte Fault